ENCODE data in the UCSC Genome Browser: year 5 update

Tracing human stem cell lineage during development using DNA methylation

Stem cell maturation is a fundamental, yet poorly understood aspect of human development. We devised a DNA methylation signature deeply reminiscent of embryonic stem cells (a fetal cell origin signature, FCO) to interrogate the evolving character of multiple human tissues. The cell fraction displaying this FCO signature was highly dependent upon developmental stage (fetal versus adult), and in leukocytes, it described a dynamic transition during the first 5 yr of life. Significant individual variation in the FCO signature of leukocytes was evident at birth, in childhood, and throughout adult life. The genes characterizing the signature included transcription factors and proteins intimately involved in embryonic development. We defined and applied a DNA methylation signature common among human fetal hematopoietic progenitor cells and have shown that this signature traces the lineage of cells and informs the study of stem cell heterogeneity in humans under homeostatic conditions.

Assessment of rosacea symptom severity by genome-wide association study and expression analysis highlights immuno-inflammatory and skin pigmentation genes

Rosacea is a common, chronic skin disease of variable severity with limited treatment options. The cause of rosacea is unknown, but it is believed to be due to a combination of hereditary and environmental factors. Little is known about the genetics of the disease. We performed a genome-wide association study (GWAS) of rosacea symptom severity with data from 73 265 research participants of European ancestry from the 23andMe customer base. Seven loci had variants associated with rosacea at the genome-wide significance level (P < 5 × 10-8). Further analyses highlighted likely gene regions or effector genes including IRF4 (P = 1.5 × 10-17), a human leukocyte antigen (HLA) region flanked by PSMB9 and HLA-DMB (P = 2.2 × 10-15), HERC2-OCA2 (P = 4.2 × 10-12), SLC45A2 (P = 1.7 × 10-10), IL13 (P = 2.8 × 10-9), a region flanked by NRXN3 and DIO2 (P = 4.1 × 10-9), and a region flanked by OVOL1and SNX32 (P = 1.2 × 10-8). All associations with rosacea were novel except for the HLA locus. Two of these loci (HERC-OCA2 and SLC45A2) and another precedented variant (rs1805007 in melanocortin 1 receptor) with an association P value just below the significance threshold (P = 1.3 × 10-7) have been previously associated with skin phenotypes and pigmentation, two of these loci are linked to immuno-inflammation phenotypes (IL13 and PSMB9-HLA-DMA) and one has been associated with both categories (IRF4). Genes within three loci (PSMB9-HLA-DMA, HERC-OCA2 and NRX3-DIO2) were differentially expressed in a previously published clinical rosacea transcriptomics study that compared lesional to non-lesional samples. The identified loci provide specificity of inflammatory mechanisms in rosacea, and identify potential pathways for therapeutic intervention.

Non-coding genetic variation shaping mental health

•

Gene expression determined by the genome mediating a response to cell environment.

•

Genetic variation results in distinct individual response in gene expression.

•

Non-coding DNA is an important site for such functional genetic variation.

•

Gene expression is a major modulator of brain chemistry and thus behavior.

Over 98% of our genome is non-coding and is now recognised to have a major role in orchestrating the tissue specific and stimulus inducible gene expression pattern which underpins our wellbeing and mental health. The non-coding genome responds functionally to our environment at all levels, encompassing the span from psychological to physiological challenge. The gene expression pattern, termed the transcriptome, ultimately gives us our neurochemistry. Therefore a major modulator of mental wellbeing is how our genes are regulated in response to life experiences. Superimposed on the aforementioned non-coding DNA framework is a vast body of genetic variation in the elements that control response to challenges. These differences, termed polymorphisms, allow for a differential response from a specific DNA element to the same challenge thus potentially allowing ‘individuality’ in the modulation of our transcriptome. This review will focus on a fundamental mechanism defining our psychological and psychiatric wellbeing, namely how genetic variation can be correlated with differential gene expression in response to specific challenges, thus resulting in altered neurochemistry which consequently may shape behaviour.

Genomic functions of developmental pluripotency associated factor 4 (Dppa4) in pluripotent stem cells and cancer

Developmental pluripotency associated factor 4 (Dppa4) is a highly specific marker of pluripotent cells, and is also overexpressed in certain cancers, but its function in either of these contexts is poorly understood. In this study, we use ChIP-Seq to identify Dppa4 binding genome-wide in three distinct cell types: mouse embryonic stem cells (mESC), embryonal carcinoma cells, and 3T3 fibroblasts ectopically expressing Dppa4. We find a core set of Dppa4 binding sites shared across cell types, and also a substantial number of sites unique to each cell type. Across cell types Dppa4 shows a preference for binding to regions with active chromatin signatures, and can influence chromatin modifications at target genes. In 3T3 fibroblasts with enforced Dppa4 expression, Dppa4 represses the cell cycle inhibitor Cdkn2c and activates Ets family transcription factor Etv4, leading to alterations in the cell cycle that likely contribute to the oncogenic phenotype. Dppa4 also directly regulates Etv4 in mESC but represses it in this context, and binds with Oct4 to a set of shared targets that are largely independent of Sox2 and Nanog, indicating that Dppa4 functions independently of the core pluripotency network in stem cells. Together these data provide novel insights into Dppa4 function in both pluripotent and oncogenic contexts.

Alternative DNA secondary structure formation affects RNA polymerase II promoter-proximal pausing in human

BACKGROUND

Alternative DNA secondary structures can arise from single-stranded DNA when duplex DNA is unwound during DNA processes such as transcription, resulting in the regulation or perturbation of these processes. We identify sites of high propensity to form stable DNA secondary structure across the human genome using Mfold and ViennaRNA programs with parameters for analyzing DNA.

RESULTS

The promoter-proximal regions of genes with paused transcription are significantly and energetically more favorable to form DNA secondary structure than non-paused genes or genes without RNA polymerase II (Pol II) binding. Using Pol II ChIP-seq, GRO-seq, NET-seq, and mNET-seq data, we arrive at a robust set of criteria for Pol II pausing, independent of annotation, and find that a highly stable secondary structure is likely to form about 10-50 nucleotides upstream of a Pol II pausing site. Structure probing data confirm the existence of DNA secondary structures enriched at the promoter-proximal regions of paused genes in human cells. Using an in vitro transcription assay, we demonstrate that Pol II pausing at HSPA1B, a human heat shock gene, is affected by manipulating DNA secondary structure upstream of the pausing site.

CONCLUSIONS

Our results indicate alternative DNA secondary structure formation as a mechanism for how GC-rich sequences regulate RNA Pol II promoter-proximal pausing genome-wide.

Nuclear respiratory factor 1 promotes spheroid survival and mesenchymal transition in mammary epithelial cells

Epithelial cells aggregate into spheroids when deprived of matrix, and the proclivity for spheroid formation and survival is a hallmark of normal and tumorigenic mammary stem cells. We show here that Nuclear Respiratory Factor 1 (NRF1) is a spheroid promoter by in silico identification of this transcription factor as highly connected to top shRNA-hits deduced from re-iterative selections for shRNAs enriched in MCF10A spheroids. NRF1-promoted spheroid survival is linked to its stimulation of mitochondrial OXPHOS, cell migration, invasion, and mesenchymal transition. Conversely, NRF1 knockdown in breast cancer MDA-MB-231 cells reduced spheroids, migration, invasion, and mesenchymal marker expression. NRF1 knockdown also reduced tumor burden in mammary fat pads and lungs of orthotopic- or tail vein-transplanted mice. With the Luminal A subtype of breast cancer, higher NRF1 expression is associated with lower survival. These results show that NRF1, an activator of mitochondrial metabolism, supports mammary spheroid survival and tumor development.

Abnormal islet sphingolipid metabolism in type 1 diabetes

AIMS/HYPOTHESIS

Sphingolipids play important roles in beta cell physiology, by regulating proinsulin folding and insulin secretion and in controlling apoptosis, as studied in animal models and cell cultures. Here we investigate whether sphingolipid metabolism may contribute to the pathogenesis of human type 1 diabetes and whether increasing the levels of the sphingolipid sulfatide would prevent models of diabetes in NOD mice.

METHODS

We examined the amount and distribution of sulfatide in human pancreatic islets by immunohistochemistry, immunofluorescence and electron microscopy. Transcriptional analysis was used to evaluate expression of sphingolipid-related genes in isolated human islets. Genome-wide association studies (GWAS) and a T cell proliferation assay were used to identify type 1 diabetes related polymorphisms and test how these affect cellular islet autoimmunity. Finally, we treated NOD mice with fenofibrate, a known activator of sulfatide biosynthesis, to evaluate the effect on experimental autoimmune diabetes development.

RESULTS

We found reduced amounts of sulfatide, 23% of the levels in control participants, in pancreatic islets of individuals with newly diagnosed type 1 diabetes, which were associated with reduced expression of enzymes involved in sphingolipid metabolism. Next, we discovered eight gene polymorphisms (ORMDL3, SPHK2, B4GALNT1, SLC1A5, GALC, PPARD, PPARG and B4GALT1) involved in sphingolipid metabolism that contribute to the genetic predisposition to type 1 diabetes. These gene polymorphisms correlated with the degree of cellular islet autoimmunity in a cohort of individuals with type 1 diabetes. Finally, using fenofibrate, which activates sulfatide biosynthesis, we completely prevented diabetes in NOD mice and even reversed the disease in half of otherwise diabetic animals.

CONCLUSIONS/INTERPRETATION

These results indicate that islet sphingolipid metabolism is abnormal in type 1 diabetes and suggest that modulation may represent a novel therapeutic approach.

DATA AVAILABILITY

The RNA expression data is available online at https://www.dropbox.com/s/93mk5tzl5fdyo6b/Abnormal%20islet%20sphingolipid%20metabolism%20in%20type%201%20diabetes%2C%20RNA%20expression.xlsx?dl=0 . A list of SNPs identified is available at https://www.dropbox.com/s/yfojma9xanpp2ju/Abnormal%20islet%20sphingolipid%20metabolism%20in%20type%201%20diabetes%20SNP.xlsx?dl=0 .

Identification of an enhancer region within the TP63/LEPREL1 locus containing genetic variants associated with bladder cancer risk

Purpose

Genome-wide association studies (GWAS) have led to the identification of a bladder cancer susceptibility variant (rs710521) in a non-coding intergenic region between the TP63 and LEPREL1 genes on chromosome 3q28, suggesting a role in the transcriptional regulation of these genes. In this study, we aimed to functionally characterize the 3q28 bladder cancer risk locus.

Methods

Fine-mapping was performed by focusing on the region surrounding rs710521, and variants were prioritized for further experiments using ENCODE regulatory data. The enhancer activity of the identified region was evaluated using dual-luciferase assays. CRISPR/Cas9-mediated deletion of the enhancer region was performed and the effect of this deletion on cell proliferation and gene expression levels was evaluated using CellTiter-Glo and RT-qPCR, respectively.

Results

Fine-mapping of the GWAS signal region led to the identification of twenty SNPs that showed a stronger association with bladder cancer risk than rs710521. Using publicly available data on regulatory elements and sequences, an enhancer region containing the bladder cancer risk variants was identified. Through reporter assays, we found that the presence of the enhancer region significantly increased ΔNTP63 promoter activity in bladder cancer-derived cell lines. CRISPR/Cas9-mediated deletion of the enhancer region reduced the viability of bladder cancer cells by decreasing the expression of ΔNTP63 and p63 target genes.

Conclusions

Taken together, our data show that bladder cancer risk-associated variants on chromosome 3q28 are located in an active enhancer region. Further characterization of the allele-specific activity of the identified enhancer and its target genes may lead to the identification of novel signaling pathways involved in bladder carcinogenesis.

Electronic supplementary material

The online version of this article (10.1007/s13402-018-0393-5) contains supplementary material, which is available to authorized users.

Towards a map of cis-regulatory sequences in the human genome

Accumulating evidence indicates that transcription factor (TF) binding sites, or cis-regulatory elements (CREs), and their clusters termed cis-regulatory modules (CRMs) play a more important role than do gene-coding sequences in specifying complex traits in humans, including the susceptibility to common complex diseases. To fully characterize their roles in deriving the complex traits/diseases, it is necessary to annotate all CREs and CRMs encoded in the human genome. However, the current annotations of CREs and CRMs in the human genome are still very limited and mostly coarse-grained, as they often lack the detailed information of CREs in CRMs. Here, we integrated 620 TF ChIP-seq datasets produced by the ENCODE project for 168 TFs in 79 different cell/tissue types and predicted an unprecedentedly completely map of CREs in CRMs in the human genome at single nucleotide resolution. The map includes 305 912 CRMs containing a total of 1 178 913 CREs belonging to 736 unique TF binding motifs. The predicted CREs and CRMs tend to be subject to either purifying selection or positive selection, thus are likely to be functional. Based on the results, we also examined the status of available ChIP-seq datasets for predicting the entire regulatory genome of humans.

The in vitro functional analysis of single-nucleotide polymorphisms associated with growth hormone (GH) response in children with GH deficiency

Response to recombinant human growth hormone (r-hGH) in the first year of therapy has been associated with single-nucleotide polymorphisms (SNPs) in children with GH deficiency (GHD). Associated SNPs were screened for regulatory function using a combination of in silico techniques. Four SNPs in regulatory sequences were selected for the analysis of in vitro transcriptional activity (TA). There was an additive effect of the alleles in the four genes associated with good growth response. For rs3110697 within IGFBP3, rs1045992 in CYP19A1 and rs2888586 in SOS1, the variant associated with better growth response showed higher TA with r-hGH treatment. For rs1024531 in GRB10, a negative regulator of IGF-I signalling and growth, the variant associated with better growth response had a significantly lower TA on r-hGH stimulation. These results indicate that specific SNP variants have effects on TA that provide a rationale for their clinical impact on growth response to r-hGH therapy.

Computational approaches for predicting key transcription factors in targeted cell reprogramming (Review)

There is a need for specific cell types in regenerative medicine and biological research. Frequently, specific cell types may not be easily obtained or the quantity obtained is insufficient for study. Therefore, reprogramming by the direct conversion (transdifferentiation) or re‑induction of induced pluripotent stem cells has been used to obtain cells expressing similar profiles to those of the desired types. Therefore, a specific cocktail of transcription factors (TFs) is required for induction. Nevertheless, identifying the correct combination of TFs is difficult. Although certain computational approaches have been proposed for this task, their methods are complex, and corresponding implementations are difficult to use and generalize for specific source or target cell types. In the present review four computational approaches that have been proposed to obtain likely TFs were compared and discussed. A simplified view of the computational complexity of these methods is provided that consists of three basic ideas: i) The definition of target and non‑target cell types; ii) the estimation of candidate TFs; and iii) filtering candidates. This simplified view was validated by analyzing a well‑documented cardiomyocyte differentiation. Subsequently, these reviewed methods were compared when applied to an unknown differentiation of corneal endothelial cells. The generated results may provide important insights for laboratory assays. Data and computer scripts that may assist with direct conversions in other cell types are also provided.

A genome-wide association study of IgM antibody against phosphorylcholine: shared genetics and phenotypic relationship to chronic lymphocytic leukemia

Phosphorylcholine (PC) is an epitope on oxidized low-density lipoprotein (oxLDL), apoptotic cells and several pathogens like Streptococcus pneumoniae. Immunoglobulin M against PC (IgM anti-PC) has the ability to inhibit uptake of oxLDL by macrophages and increase clearance of apoptotic cells. From our genome-wide association studies (GWASs) in four European-ancestry cohorts, six single nucleotide polymorphisms (SNPs) in 11q24.1 were discovered (in 3002 individuals) and replicated (in 646 individuals) to be associated with serum level of IgM anti-PC (the leading SNP rs35923643-G, combined β = 0.19, 95% confidence interval 0.13-0.24, P = 4.3 × 10-11). The haplotype tagged by rs35923643-G (or its proxy SNP rs735665-A) is also known as the top risk allele for chronic lymphocytic leukemia (CLL), and a main increasing allele for general IgM. By using summary GWAS results of IgM anti-PC and CLL in the polygenic risk score (PRS) analysis, PRS on the basis of IgM anti-PC risk alleles positively associated with CLL risk (explained 0.6% of CLL variance, P = 1.2 × 10-15). Functional prediction suggested that rs35923643-G might impede the binding of Runt-related transcription factor 3, a tumor suppressor playing a central role in the immune regulation of cancers. Contrary to the expectations from the shared genetics between IgM anti-PC and CLL, an inverse relationship at the phenotypic level was found in a nested case-control study (30 CLL cases with 90 age- and sex-matched controls), potentially reflecting reverse causation. The suggested function of the top variant as well as the phenotypic association between IgM anti-PC and CLL risk needs replication and motivates further studies.

MAX Mutations in Endometrial Cancer: Clinicopathologic Associations and Recurrent MAX p.His28Arg Functional Characterization

Background

Genomic studies have revealed that multiple genes are mutated at varying frequency in endometrial cancer (EC); however, the relevance of many of these mutations is poorly understood. An EC-specific recurrent mutation in the MAX transcription factor p.His28Arg was recently discovered. We sought to assess the functional consequences of this hotspot mutation and determine its association with cancer-relevant phenotypes.

Methods

MAX was sequenced in 509 endometrioid ECs, and associations between mutation status and clinicopathologic features were assessed. EC cell lines stably expressing MAXH28R were established and used for functional experiments. DNA binding was examined using electrophoretic mobility shift assays and chromatin immunoprecipitation. Transcriptional profiling was performed with microarrays. Murine flank (six to 11 mice per group) and intraperitoneal tumor models were used for in vivo studies. Vascularity of xenografts was assessed by MECA-32 immunohistochemistry. The paracrine pro-angiogenic nature of MAXH28R-expressing EC cells was tested using microfluidic HUVEC sprouting assays and VEGFA enzyme-linked immunosorbent assays. All statistical tests were two-sided.

Results

Twenty-two of 509 tumors harbored mutations in MAX, including 12 tumors with the p.His28Arg mutation. Patients with a MAX mutation had statistically significantly reduced recurrence-free survival (hazard ratio = 4.00, 95% confidence interval = 1.15 to 13.91, P = .03). MAXH28R increased affinity for canonical E-box sequences, and MAXH28R-expressing EC cells dramatically altered transcriptional profiles. MAXH28R-derived xenografts statistically significantly increased vascular area compared with MAXWT and empty vector tumors (P = .003 and P = .008, respectively). MAXH28R-expressing EC cells secreted nearly double the levels of VEGFA compared with MAXWT cells (P = .03, .005, and .005 at 24, 48, and 72 hours, respectively), and conditioned media from MAXH28R cells increased sprouting when applied to HUVECs.

Conclusion

These data highlight the importance of MAX mutations in EC and point to increased vascularity as one mechanism contributing to clinical aggressiveness of EC.

UFD1 contributes to MYC-mediated leukemia aggressiveness through suppression of the proapoptotic unfolded protein response

Despite the pivotal role of MYC in tumorigenesis, the mechanisms by which it promotes cancer aggressiveness remain incompletely understood. Here we show that MYC transcriptionally upregulates the ubiquitin fusion degradation 1 (UFD1) gene in T-cell acute lymphoblastic leukemia (T-ALL). Allelic loss of ufd1 in zebrafish induces tumor-cell apoptosis and impairs MYC-driven T-ALL progression but does not affect general health. As the E2 component of an endoplasmic reticulum (ER)-associated degradation (ERAD) complex, UFD1 facilitates the elimination of misfolded/unfolded proteins from the ER. We found that UFD1 inactivation in human T-ALL cells impairs ERAD, exacerbates ER stress, and induces apoptosis. Moreover, we show that UFD1 inactivation promotes the proapoptotic unfolded protein response (UPR) mediated by protein kinase RNA-like ER kinase (PERK). This effect is demonstrated by an upregulation of PERK and its downstream effector C/EBP homologous protein (CHOP), as well as a downregulation of BCL2 and BCLxL. Indeed, CHOP inactivation or BCL2 overexpression is sufficient to rescue tumor-cell apoptosis induced by UFD1 knockdown. Together, our studies identify UFD1 as a critical regulator of the ER stress response and a novel contributor to MYC-mediated leukemia aggressiveness, with implications for targeted therapy in T-ALL and likely other MYC-driven cancers.

The multiple myeloma risk allele at 5q15 lowers ELL2 expression and increases ribosomal gene expression

Recently, we identified ELL2 as a susceptibility gene for multiple myeloma (MM). To understand its mechanism of action, we performed expression quantitative trait locus analysis in CD138⁺ plasma cells from 1630 MM patients from four populations. We show that the MM risk allele lowers ELL2 expression in these cells (P_combined = 2.5 × 10⁻²⁷; β_combined = −0.24 SD), but not in peripheral blood or other tissues. Consistent with this, several variants representing the MM risk allele map to regulatory genomic regions, and three yield reduced transcriptional activity in plasmocytoma cell lines. One of these (rs3777189-C) co-locates with the best-supported lead variants for ELL2 expression and MM risk, and reduces binding of MAFF/G/K family transcription factors. Moreover, further analysis reveals that the MM risk allele associates with upregulation of gene sets related to ribosome biogenesis, and knockout/knockdown and rescue experiments in plasmocytoma cell lines support a cause–effect relationship. Our results provide mechanistic insight into MM predisposition.

ELL2 was recently discovered as a susceptibility gene for multiple myeloma (MM). Here, they show that the MM risk allele lowers ELL2 expression in plasma cells, that it also upregulates gene sets related to ribosome biogenesis, and that one of the linked variants reduces binding of MAFF/G/K family transcription factors.

LDSplitDB: a database for studies of meiotic recombination hotspots in MHC using human genomic data

Background

Meiotic recombination happens during the process of meiosis when chromosomes inherited from two parents exchange genetic materials to generate chromosomes in the gamete cells. The recombination events tend to occur in narrow genomic regions called recombination hotspots. Its dysregulation could lead to serious human diseases such as birth defects. Although the regulatory mechanism of recombination events is still unclear, DNA sequence polymorphisms have been found to play crucial roles in the regulation of recombination hotspots.

Method

To facilitate the studies of the underlying mechanism, we developed a database named LDSplitDB which provides an integrative and interactive data mining and visualization platform for the genome-wide association studies of recombination hotspots. It contains the pre-computed association maps of the major histocompatibility complex (MHC) region in the 1000 Genomes Project and the HapMap Phase III datasets, and a genome-scale study of the European population from the HapMap Phase II dataset. Besides the recombination profiles, related data of genes, SNPs and different types of epigenetic modifications, which could be associated with meiotic recombination, are provided for comprehensive analysis. To meet the computational requirement of the rapidly increasing population genomics data, we prepared a lookup table of 400 haplotypes for recombination rate estimation using the well-known LDhat algorithm which includes all possible two-locus haplotype configurations.

Conclusion

To the best of our knowledge, LDSplitDB is the first large-scale database for the association analysis of human recombination hotspots with DNA sequence polymorphisms. It provides valuable resources for the discovery of the mechanism of meiotic recombination hotspots. The information about MHC in this database could help understand the roles of recombination in human immune system.

Database URL

http://histone.scse.ntu.edu.sg/LDSplitDB

ViFi: accurate detection of viral integration and mRNA fusion reveals indiscriminate and unregulated transcription in proximal genomic regions in cervical cancer

The integration of viral sequences into the host genome is an important driver of tumorigenesis in many viral mediated cancers, notably cervical cancer and hepatocellular carcinoma. We present ViFi, a computational method that combines phylogenetic methods with reference-based read mapping to detect viral integrations. In contrast with read-based reference mapping approaches, ViFi is faster, and shows high precision and sensitivity on both simulated and biological data, even when the integrated virus is a novel strain or highly mutated. We applied ViFi to matched genomic and mRNA data from 68 cervical cancer samples from TCGA and found high concordance between the two. Surprisingly, viral integration resulted in a dramatic transcriptional upregulation in all proximal elements, including LINEs and LTRs that are not normally transcribed. This upregulation is highly correlated with the presence of a viral gene fused with a downstream human element. Moreover, genomic rearrangements suggest the formation of apparent circular extrachromosomal (ecDNA) human-viral structures. Our results suggest the presence of apparent small circular fusion viral/human ecDNA, which correlates with indiscriminate and unregulated expression of proximal genomic elements, potentially contributing to the pathogenesis of HPV-associated cervical cancers. ViFi is available at https://github.com/namphuon/ViFi.

Silencing of transposable elements may not be a major driver of regulatory evolution in primate iPSCs

Transposable elements (TEs) comprise almost half of primate genomes and their aberrant regulation can result in deleterious effects. In pluripotent stem cells, rapidly evolving KRAB-ZNF genes target TEs for silencing by H3K9me3. To investigate the evolution of TE silencing, we performed H3K9me3 ChIP-seq experiments in induced pluripotent stem cells from 10 human and 7 chimpanzee individuals. We identified four million orthologous TEs and found the SVA and ERV families to be marked most frequently by H3K9me3. We found little evidence of inter-species differences in TE silencing, with as many as 82% of putatively silenced TEs marked at similar levels in humans and chimpanzees. TEs that are preferentially silenced in one species are a similar age to those silenced in both species and are not more likely to be associated with expression divergence of nearby orthologous genes. Our data suggest limited species-specificity of TE silencing across 6 million years of primate evolution.

Integrated analysis sheds light on evolutionary trajectories of young transcription start sites in the human genome

Understanding the molecular mechanisms and evolution of the gene regulatory system remains a major challenge in biology. Transcription start sites (TSSs) are especially interesting because they are central to initiating gene expression. Previous studies revealed widespread transcription initiation and fast turnover of TSSs in mammalian genomes. Yet, how new TSSs originate and how they evolve over time remain poorly understood. To address these questions, we analyzed ∼200,000 human TSSs by integrating evolutionary (inter- and intra-species) and functional genomic data, particularly focusing on evolutionarily young TSSs that emerged in the primate lineage. TSSs were grouped according to their evolutionary age using sequence alignment information as a proxy. Comparisons of young and old TSSs revealed that (1) new TSSs emerge through a combination of intrinsic factors, like the sequence properties of transposable elements and tandem repeats, and extrinsic factors such as their proximity to existing regulatory modules; (2) new TSSs undergo rapid evolution that reduces the inherent instability of repeat sequences associated with a high propensity of TSS emergence; and (3) once established, the transcriptional competence of surviving TSSs is gradually enhanced, with evolutionary changes subject to temporal (fewer regulatory changes in younger TSSs) and spatial constraints (fewer regulatory changes in more isolated TSSs). These findings advance our understanding of how regulatory innovations arise in the genome throughout evolution and highlight the genomic robustness and evolvability in these processes.

Tissue-specific impact of FADS cluster variants on FADS1 and FADS2 gene expression

Omega-6 (n-6) and omega-3 (n-3) long (≥ 20 carbon) chain polyunsaturated fatty acids (LC-PUFAs) play a critical role in human health and disease. Biosynthesis of LC-PUFAs from dietary 18 carbon PUFAs in tissues such as the liver is highly associated with genetic variation within the fatty acid desaturase (FADS) gene cluster, containing FADS1 and FADS2 that encode the rate-limiting desaturation enzymes in the LC-PUFA biosynthesis pathway. However, the molecular mechanisms by which FADS genetic variants affect LC-PUFA biosynthesis, and in which tissues, are unclear. The current study examined associations between common single nucleotide polymorphisms (SNPs) within the FADS gene cluster and FADS1 and FADS2 gene expression in 44 different human tissues (sample sizes ranging 70-361) from the Genotype-Tissue Expression (GTEx) Project. FADS1 and FADS2 expression were detected in all 44 tissues. Significant cis-eQTLs (within 1 megabase of each gene, False Discovery Rate, FDR<0.05, as defined by GTEx) were identified in 12 tissues for FADS1 gene expression and 23 tissues for FADS2 gene expression. Six tissues had significant (FDR< 0.05) eQTLs associated with both FADS1 and FADS2 (including artery, esophagus, heart, muscle, nerve, and thyroid). Interestingly, the identified eQTLs were consistently found to be associated in opposite directions for FADS1 and FADS2 expression. Taken together, findings from this study suggest common SNPs within the FADS gene cluster impact the transcription of FADS1 and FADS2 in numerous tissues and raise important questions about how the inverse expression of these two genes impact intermediate molecular (such a LC-PUFA and LC-PUFA-containing glycerolipid levels) and ultimately clinical phenotypes associated with inflammatory diseases and brain health.

Identification and characterization of cis-regulatory elements 'insulator and repressor' in PPARD gene

AIM

Identification and functional characterization of cis-regulatory elements in human PPARD gene.

METHODS

We used various bioinformatic tools on the publicly available human genome and Encyclopedia of DNA Elements databases to explore potential cis-regulatory elements in PPARD gene region.

RESULTS

We predicted an insulator and an enhancer element in intron 2 of PPARD gene. Functional characterization using transient transfection, reporter assay and CTCF binding confirmed the insulator status. However, the predicted enhancer element showed repressor/silencer activity. Finally, we observed a potential interaction between these two cis-regulatory elements which is in agreement with 5C-Encyclopedia of DNA Elements data.

CONCLUSION

We report two functionally validated cis-regulatory elements in PPARD gene which will aid in understanding its regulation and role in metabolic functions.

annoPeak: a web application to annotate and visualize peaks from ChIP-seq/ChIP-exo-seq

Summary

We developed annoPeak, a web application to annotate, visualize and compare predicted protein-binding regions derived from ChIP-seq/ChIP-exo-seq experiments using human and mouse cells. Users can upload peak regions from multiple experiments onto the annoPeak server to annotate them with biological context, identify associated target genes and categorize binding sites with respect to gene structure. Users can also compare multiple binding profiles intuitively with the help of visualization tools and tables provided by annoPeak. In general, annoPeak will help users identify patterns of genome wide transcription factor binding profiles, assess binding profiles in different biological contexts and generate new hypotheses.

Availability and Implementation

The web service is freely accessible through URL: http://ccc-annopeak.osumc.edu/annoPeak . Source code is available at https://github.com/XingTang2014/annoPeak .

Contact

gustavo.leone@osumc.edu or kun.huang@osumc.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

A point mutation in the pre-ZRS disrupts sonic hedgehog expression in the limb bud and results in triphalangeal thumb-polysyndactyly syndrome

PURPOSE

The zone of polarizing activity regulatory sequence (ZRS) is an enhancer that regulates sonic hedgehog during embryonic limb development. Recently, mutations in a noncoding evolutionary conserved sequence 500 bp upstream of the ZRS, termed the pre-ZRS (pZRS), have been associated with polydactyly in dogs and humans. Here, we report the first case of triphalangeal thumb-polysyndactyly syndrome (TPT-PS) to be associated with mutations in this region and show via mouse enhancer assays how this mutation leads to ectopic expression throughout the developing limb bud.

METHODS

We used linkage analysis, whole-exome sequencing, Sanger sequencing, fluorescence in situ hybridization, multiplex ligation-dependent probe amplification, single-nucleotide polymorphism array, and a mouse transgenic enhancer assay.

RESULTS

Ten members of a TPT-PS family were included in this study. The mutation was linked to chromosome 7q36 (LOD score 3.0). No aberrations in the ZRS could be identified. A point mutation in the pZRS (chr7:156585476G>C; GRCh37/hg19) was detected in all affected family members. Functional characterization using a mouse transgenic enhancer essay showed extended ectopic expression dispersed throughout the entire limb bud (E11.5).

CONCLUSION

Our work describes the first mutation in the pZRS to be associated with TPT-PS and provides functional evidence that this mutation leads to ectopic expression of this enhancer within the developing limb.

Induction of anti-VEGF therapy resistance by upregulated expression of microseminoprotein (MSMP)

Anti-vascular endothelial growth factor (VEGF) therapy has demonstrated efficacy in treating human metastatic cancers, but therapeutic resistance is a practical limitation and most tumors eventually become unresponsive. To identify microenvironmental factors underlying the resistance of cancer to antiangiogenesis therapy, we conducted genomic analyses of intraperitoneal ovarian tumors in which adaptive resistance to anti-VEGF therapy (B20 antibody) developed. We found that expression of the microseminoprotein, prostate-associated (MSMP) gene was substantially upregulated in resistant compared with control tumors. MSMP secretion from cancer cells was induced by hypoxia, triggering MAPK signaling in endothelial cells to promote tube formation in vitro. Recruitment of the transcriptional repressor CCCTC-binding factor (CTCF) to the MSMP enhancer region was decreased by histone acetylation under hypoxic conditions in cancer cells. MSMP siRNA, delivered in vivo using the DOPC nanoliposomes, restored tumor sensitivity to anti-VEGF therapy. In ovarian cancer patients treated with bevacizumab, serum MSMP concentration increased significantly only in non-responders. These findings imply that MSMP inhibition combined with the use of antiangiogenesis drugs may be a new strategy to overcome resistance to antiangiogenesis therapy.

A null variant in the apolipoprotein L3 gene is associated with non-diabetic nephropathy

Background

Inheritance of apolipoprotein L1 gene (APOL1) renal-risk variants in a recessive pattern strongly associates with non-diabetic end-stage kidney disease (ESKD). Further evidence supports risk modifiers in APOL1-associated nephropathy; some studies demonstrate that heterozygotes possess excess risk for ESKD or show earlier age at ESKD, relative to those with zero risk alleles. Nearby loci are also associated with ESKD in non-African Americans.

Methods

We assessed the role of the APOL3 null allele rs11089781 on risk of non-diabetic ESKD. Four cohorts containing 2781 ESKD cases and 2474 controls were analyzed.

Results

Stratifying by APOL1 risk genotype (recessive) and adjusting for African ancestry identified a significant additive association between rs11089781 and ESKD in each stratum and in a meta-analysis [meta-analysis P  =  0.0070; odds ratio (OR) = 1.29]; ORs were consistent across APOL1 risk strata. The biological significance of this association is supported by the finding that the APOL3 gene is co-regulated with APOL1, and that APOL3 protein was able to bind to APOL1 protein.

Conclusions

Taken together, the genetic and biological data support the concept that other APOL proteins besides APOL1 may also influence the risk of non-diabetic ESKD.

SiNoPsis: Single Nucleotide Polymorphisms selection and promoter profiling

MOTIVATION

The selection of a single nucleotide polymorphism (SNP) using bibliographic methods can be a very time-consuming task. Moreover, a SNP selected in this way may not be easily visualized in its genomic context by a standard user hoping to correlate it with other valuable information. Here we propose a web form built on top of Circos that can assist SNP-centered screening, based on their location in the genome and the regulatory modules they can disrupt. Its use may allow researchers to prioritize SNPs in genotyping and disease studies.

RESULTS

SiNoPsis is bundled as a web portal. It focuses on the different structures involved in the genomic expression of a gene, especially those found in the core promoter upstream region. These structures include transcription factor binding sites (for promoter and enhancer signals), histones and promoter flanking regions. Additionally, the tool provides eQTL and linkage disequilibrium (LD) properties for a given SNP query, yielding further clues about other indirectly associated SNPs. Possible disruptions of the aforementioned structures affecting gene transcription are reported using multiple resource databases. SiNoPsis has a simple user-friendly interface, which allows single queries by gene symbol, genomic coordinates, Ensembl gene identifiers, RefSeq transcript identifiers and SNPs. It is the only portal providing useful SNP selection based on regulatory modules and LD with functional variants in both textual and graphic modes (by properly defining the arguments and parameters needed to run Circos).

AVAILABILITY AND IMPLEMENTATION

SiNoPsis is freely available at https://compgen.bio.ub.edu/SiNoPsis/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Ensemble genomic analysis in human lung tissue identifies novel genes for chronic obstructive pulmonary disease

BACKGROUND

Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) significantly associated with chronic obstructive pulmonary disease (COPD). However, many genetic variants show suggestive evidence for association but do not meet the strict threshold for genome-wide significance. Integrative analysis of multiple omics datasets has the potential to identify novel genes involved in disease pathogenesis by leveraging these variants in a functional, regulatory context.

RESULTS

We performed expression quantitative trait locus (eQTL) analysis using genome-wide SNP genotyping and gene expression profiling of lung tissue samples from 86 COPD cases and 31 controls, testing for SNPs associated with gene expression levels. These results were integrated with a prior COPD GWAS using an ensemble statistical and network methods approach to identify relevant genes and observe them in the context of overall genetic control of gene expression to highlight co-regulated genes and disease pathways. We identified 250,312 unique SNPs and 4997 genes in the cis(local)-eQTL analysis (5% false discovery rate). The top gene from the integrative analysis was MAPT, a gene recently identified in an independent GWAS of lung function. The genes HNRNPAB and PCBP2 with RNA binding activity and the gene ACVR1B were identified in network communities with validated disease relevance.

CONCLUSIONS

The integration of lung tissue gene expression with genome-wide SNP genotyping and subsequent intersection with prior GWAS and omics studies highlighted candidate genes within COPD loci and in communities harboring known COPD genes. This integration also identified novel disease genes in sub-threshold regions that would otherwise have been missed through GWAS.

Contribution of transposable elements and distal enhancers to evolution of human-specific features of interphase chromatin architecture in embryonic stem cells

Transposable elements have made major evolutionary impacts on creation of primate-specific and human-specific genomic regulatory loci and species-specific genomic regulatory networks (GRNs). Molecular and genetic definitions of human-specific changes to GRNs contributing to development of unique to human phenotypes remain a highly significant challenge. Genome-wide proximity placement analysis of diverse families of human-specific genomic regulatory loci (HSGRL) identified topologically associating domains (TADs) that are significantly enriched for HSGRL and designated rapidly evolving in human TADs. Here, the analysis of HSGRL, hESC-enriched enhancers, super-enhancers (SEs), and specific sub-TAD structures termed super-enhancer domains (SEDs) has been performed. In the hESC genome, 331 of 504 (66%) of SED-harboring TADs contain HSGRL and 68% of SEDs co-localize with HSGRL, suggesting that emergence of HSGRL may have rewired SED-associated GRNs within specific TADs by inserting novel and/or erasing existing non-coding regulatory sequences. Consequently, markedly distinct features of the principal regulatory structures of interphase chromatin evolved in the hESC genome compared to mouse: the SED quantity is 3-fold higher and the median SED size is significantly larger. Concomitantly, the overall TAD quantity is increased by 42% while the median TAD size is significantly decreased (p = 9.11E-37) in the hESC genome. Present analyses illustrate a putative global role for transposable elements and HSGRL in shaping the human-specific features of the interphase chromatin organization and functions, which are facilitated by accelerated creation of novel transcription factor binding sites and new enhancers driven by targeted placement of HSGRL at defined genomic coordinates. A trend toward the convergence of TAD and SED architectures of interphase chromatin in the hESC genome may reflect changes of 3D-folding patterns of linear chromatin fibers designed to enhance both regulatory complexity and functional precision of GRNs by creating predominantly a single gene (or a set of functionally linked genes) per regulatory domain structures. Collectively, present analyses reveal critical evolutionary contributions of transposable elements and distal enhancers to creation of thousands primate- and human-specific elements of a chromatin folding code, which defines the 3D context of interphase chromatin both restricting and facilitating biological functions of GRNs.

Genome-Wide Association Studies and Meta-Analyses for Congenital Heart Defects

BACKGROUND

Maternal and inherited (ie, case) genetic factors likely contribute to the pathogenesis of congenital heart defects, but it is unclear whether individual common variants confer a large risk.

METHODS AND RESULTS

To evaluate the relationship between individual common maternal/inherited genotypes and risk for heart defects, we conducted genome-wide association studies in 5 cohorts. Three cohorts were recruited at the Children's Hospital of Philadelphia: 670 conotruncal heart defect (CTD) case-parent trios, 317 left ventricular obstructive tract defect (LVOTD) case-parent trios, and 406 CTD cases (n=406) and 2976 pediatric controls. Two cohorts were recruited through the Pediatric Cardiac Genomics Consortium: 355 CTD trios and 192 LVOTD trios. We also conducted meta-analyses using the genome-wide association study results from the CTD cohorts, the LVOTD cohorts, and from the combined CTD and LVOTD cohorts. In the individual genome-wide association studies, several genome-wide significant associations (P≤5×10^-8) were observed. In our meta-analyses, 1 genome-wide significant association was detected: the case genotype for rs72820264, an intragenetic single-nucleotide polymorphism associated with LVOTDs (P=2.1×10^-8).

CONCLUSIONS

We identified 1 novel candidate region associated with LVOTDs and report on several additional regions with suggestive evidence for association with CTD and LVOTD. These studies were constrained by the relatively small samples sizes and thus have limited power to detect small to moderate associations. Approaches that minimize the multiple testing burden (eg, gene or pathway based) may, therefore, be required to uncover common variants contributing to the risk of these relatively rare conditions.

The SIN3A histone deacetylase complex is required for a complete transcriptional response to hypoxia

Cells adapt to environmental changes, including fluctuations in oxygen levels, through the induction of specific gene expression programs. To identify genes regulated by hypoxia at the transcriptional level, we pulse-labeled HUVEC cells with 4-thiouridine and sequenced nascent transcripts. Then, we searched genome-wide binding profiles from the ENCODE project for factors that correlated with changes in transcription and identified binding of several components of the Sin3A co-repressor complex, including SIN3A, SAP30 and HDAC1/2, proximal to genes repressed by hypoxia. SIN3A interference revealed that it participates in the downregulation of 75% of the hypoxia-repressed genes in endothelial cells. Unexpectedly, it also blunted the induction of 47% of the upregulated genes, suggesting a role for this corepressor in gene induction. In agreement, ChIP-seq experiments showed that SIN3A preferentially localizes to the promoter region of actively transcribed genes and that SIN3A signal was enriched in hypoxia-repressed genes, prior exposure to the stimulus. Importantly, SINA3 occupancy was not altered by hypoxia in spite of changes in H3K27ac signal. In summary, our results reveal a prominent role for SIN3A in the transcriptional response to hypoxia and suggest a model where modulation of the associated histone deacetylase activity, rather than its recruitment, determines the transcriptional output.

Identification of specific feed-forward apoptosis mechanisms and associated higher survival rates for low grade glioma and lung squamous cell carcinoma

The mechanisms of cell proliferation due to the overexpression of certain transcription factors (TFs) have been well documented in the cancer setting. However, many of these same TFs have pro-apoptotic effects, particularly when expressed or activated at high levels, a process referred to as feed-forward apoptosis (FFA). To determine whether cancers could be stratified on the basis of specific FFA signatures, RNASeq data representing samples from the cancer genome atlas were analyzed, revealing that high expression of the pro-proliferative TFs, MYC and YY1, is associated with a favorable outcome in low-grade glioma (LGG) and lung squamous cell carcinoma (LUSC), respectively. Analysis of the RNASeq data also led to the identification of specific apoptosis-effector genes whose expression levels correlate with increased survival rates, for both LGG and LUSC. Although FFA has been demonstrated as a general effect in cancer, in this report, for the first time, results identify specific TFs and their responsive effector genes that distinguish subsets of cancer samples undergoing more or less of a FFA process in a way that is associated with distinct patient survival rates.

HEDD: Human Enhancer Disease Database

Enhancers, as specialized genomic cis-regulatory elements, activate transcription of their target genes and play an important role in pathogenesis of many human complex diseases. Despite recent systematic identification of them in the human genome, currently there is an urgent need for comprehensive annotation databases of human enhancers with a focus on their disease connections. In response, we built the Human Enhancer Disease Database (HEDD) to facilitate studies of enhancers and their potential roles in human complex diseases. HEDD currently provides comprehensive genomic information for ∼2.8 million human enhancers identified by ENCODE, FANTOM5 and RoadMap with disease association scores based on enhancer-gene and gene-disease connections. It also provides Web-based analytical tools to visualize enhancer networks and score enhancers given a set of selected genes in a specific gene network. HEDD is freely accessible at http://zdzlab.einstein.yu.edu/1/hedd.php.

TERT structural rearrangements in metastatic pheochromocytomas

Pheochromocytomas (PC) and paragangliomas (PGL) are endocrine tumors for which the genetic and clinicopathological features of metastatic progression remain incompletely understood. As a result, the risk of metastasis from a primary tumor cannot be predicted. Early diagnosis of individuals at high risk of developing metastases is clinically important and the identification of new biomarkers that are predictive of metastatic potential is of high value. Activation of TERT has been associated with a number of malignant tumors, including PC/PGL. However, the mechanism of TERT activation in the majority of PC/PGL remains unclear. As TERT promoter mutations occur rarely in PC/PGL, we hypothesized that other mechanisms - such as structural variations - may underlie TERT activation in these tumors. From 35 PC and four PGL, we identified three primary PCs that developed metastases with elevated TERT expression, each of which lacked TERT promoter mutations and promoter DNA methylation. Using whole genome sequencing, we identified somatic structural alterations proximal to the TERT locus in two of these tumors. In both tumors, the genomic rearrangements led to the positioning of super-enhancers proximal to the TERT promoter, that are likely responsible for the activation of the normally tightly repressed TERT expression in chromaffin cells.

The Sorghum bicolor reference genome: improved assembly, gene annotations, a transcriptome atlas, and signatures of genome organization

Sorghum bicolor is a drought tolerant C4 grass used for the production of grain, forage, sugar, and lignocellulosic biomass and a genetic model for C4 grasses due to its relatively small genome (approximately 800 Mbp), diploid genetics, diverse germplasm, and colinearity with other C4 grass genomes. In this study, deep sequencing, genetic linkage analysis, and transcriptome data were used to produce and annotate a high-quality reference genome sequence. Reference genome sequence order was improved, 29.6 Mbp of additional sequence was incorporated, the number of genes annotated increased 24% to 34 211, average gene length and N50 increased, and error frequency was reduced 10-fold to 1 per 100 kbp. Subtelomeric repeats with characteristics of Tandem Repeats in Miniature (TRIM) elements were identified at the termini of most chromosomes. Nucleosome occupancy predictions identified nucleosomes positioned immediately downstream of transcription start sites and at different densities across chromosomes. Alignment of more than 50 resequenced genomes from diverse sorghum genotypes to the reference genome identified approximately 7.4 M single nucleotide polymorphisms (SNPs) and 1.9 M indels. Large-scale variant features in euchromatin were identified with periodicities of approximately 25 kbp. A transcriptome atlas of gene expression was constructed from 47 RNA-seq profiles of growing and developed tissues of the major plant organs (roots, leaves, stems, panicles, and seed) collected during the juvenile, vegetative and reproductive phases. Analysis of the transcriptome data indicated that tissue type and protein kinase expression had large influences on transcriptional profile clustering. The updated assembly, annotation, and transcriptome data represent a resource for C4 grass research and crop improvement.

Pharmacogenetics of Adverse Drug Reactions

Adverse drug reactions (ADRs) are an important cause of morbidity and mortality. Genetic factors predispose to many ADRs, affecting susceptibility to both type A and type B reactions. The overall contribution of genetics will vary according to drug and ADR, and should be considered when attempting to predict and prevent ADRs. Genetic risk factors are considered in detail for a number of type A ADRs, especially those relating to warfarin and thiopurines, and type B ADRs affecting skin, the liver, and the heart. As the availability of whole genome sequencing increases, it is likely that prospective genotype for particular ADRs prior to drug prescription will become more common in the future. Current examples of genetic testing to prevent ADRs which have already been implemented and future prospects for developments in the field are discussed in detail.

MicroRNA-transcription factor interactions and their combined effect on target gene expression in colon cancer cases

Transcription factors (TFs) and microRNAs (miRNAs) regulate gene expression: TFs by influencing messenger RNA (mRNA) transcription and miRNAs by influencing mRNA translation and transcript degradation. Additionally, miRNAs and TFs alter each other's expression, making it difficult to ascertain the effect either one has on target gene (TG) expression. In this investigation, we use a two‐way interaction model with the TF and miRNA as independent variables to investigate whether miRNAs and TFs work together to influence TG expression levels in colon cancer subjects. We used known TF binding sites and validated miRNA targets to determine potential miRNA‐TF‐TG interactions, restricting interactions to those with a TF previously associated with altered risk of colorectal cancer death. We analyzed interactions using normal colonic mucosa expression as well as differential expression, which is measured as colonic carcinoma expression minus normal colonic mucosa expression. We analyzed 3518 miRNA‐TF‐TG triplets using normal mucosa expression and 617 triplets using differential expression. Normal colonic RNA‐Seq data were available for 168 individuals; of these, 159 also had carcinoma RNA‐Seq data. Thirteen unique miRNA‐TF‐TG interactions, comprising six miRNAs, four TFs, and 11 TGs, were statistically significant after adjustment for multiple comparisons in normal colonic mucosa, and 14 unique miRNA‐TF‐TG interactions, comprising two miRNAs, two TFs, and 13 TGs, were found for carcinoma‐normal differential expression. Our results show that TG expression is influenced by both miRNAs as well as TFs, and the influence of one regulator impacts the effect of the other on the shared TG expression.

Allele-Specific DNA Methylation and Its Interplay with Repressive Histone Marks at Promoter-Mutant TERT Genes

A mutation in the promoter of the Telomerase Reverse Transcriptase (TERT) gene is the most frequent noncoding mutation in cancer. The mutation drives unusual monoallelic expression of TERT, allowing immortalization. Here, we find that DNA methylation of the TERT CpG island (CGI) is also allele-specific in multiple cancers. The expressed allele is hypomethylated, which is opposite to cancers without TERT promoter mutations. The continued presence of Polycomb repressive complex 2 (PRC2) on the inactive allele suggests that histone marks of repressed chromatin may be causally linked to high DNA methylation. Consistent with this hypothesis, TERT promoter DNA containing 5-methyl-CpG has much increased affinity for PRC2 in vitro. Thus, CpG methylation and histone marks appear to collaborate to maintain the two TERT alleles in different epigenetic states in TERT promoter mutant cancers. Finally, in several cancers, DNA methylation levels at the TERT CGI correlate with altered patient survival.

"DNA Methylation signatures in panic disorder"

Panic disorder (PD) affects about four million Europeans, with women affected twice as likely as men, causing substantial suffering and high economic costs. The etiopathogenesis of PD remains largely unknown, but both genetic and environmental factors contribute to risk. An epigenome-wide association study (EWAS) was conducted to compare medication-free PD patients (n = 89) with healthy controls (n = 76) stratified by gender. Replication was sought in an independent sample (131 cases, 169 controls) and functional analyses were conducted in a third sample (N = 71). DNA methylation was assessed in whole blood using the Infinium HumanMethylation450 BeadChip. One genome-wide association surviving FDR of 5% (cg07308824, P = 1.094 × 10-7, P-adj = 0.046) was identified in female PD patients (N = 49) compared to controls (N = 48). The same locus, located in an enhancer region of the HECA gene, was also hypermethylated in female PD patients in the replication sample (P = 0.035) and the significance of the association improved in the meta-analysis (P-adj = 0.004). Methylation at this CpG site was associated with HECA mRNA expression in another independent female sample (N = 71) both at baseline (P = 0.046) and after induction by dexamethasone (P = 0.029). Of 15 candidates, 5 previously reported as associated with PD or anxiety traits also showed differences in DNA methylation after gene-wise correction and included SGK1, FHIT, ADCYAP1, HTR1A, HTR2A. Our study examines epigenome-wide differences in peripheral blood for PD patients. Our results point to possible sex-specific methylation changes in the HECA gene for PD but overall highlight that this disorder is not associated with extensive changes in DNA methylation in peripheral blood.

The RNA processing factors THRAP3 and BCLAF1 promote the DNA damage response through selective mRNA splicing and nuclear export

mRNA splicing and export plays a key role in the regulation of gene expression, with recent evidence suggesting an additional layer of regulation of gene expression and cellular function through the selective splicing and export of genes within specific pathways. Here we describe a role for the RNA processing factors THRAP3 and BCLAF1 in the regulation of the cellular DNA damage response (DDR) pathway, a key pathway involved in the maintenance of genomic stability and the prevention of oncogenic transformation. We show that loss of THRAP3 and/or BCLAF1 leads to sensitivity to DNA damaging agents, defective DNA repair and genomic instability. Additionally, we demonstrate that this phenotype can be at least partially explained by the role of THRAP3 and BCLAF1 in the selective mRNA splicing and export of transcripts encoding key DDR proteins, including the ATM kinase. Moreover, we show that cancer associated mutations within THRAP3 result in deregulated processing of THRAP3/BCLAF1-regulated transcripts and consequently defective DNA repair. Taken together, these results suggest that THRAP3 and BCLAF1 mutant tumors may be promising targets for DNA damaging chemotherapy.

5-Azacytidine engages an IRE1α-EGFR-ERK1/2 signaling pathway that stabilizes the LDL receptor mRNA

Hepatic low-density lipoprotein receptor (LDLR) is the primary conduit for the clearance of plasma LDL-cholesterol and increasing its expression represents a central goal for treating cardiovascular disease. However, LDLR mRNA is unstable and undergoes rapid turnover mainly due to the three AU-rich elements (ARE) in its proximal 3'-untranslated region (3'-UTR). Herein, our data revealed that 5-azacytidine (5-AzaC), an antimetabolite used in the treatment of myelodysplastic syndrome, stabilizes the LDLR mRNA through a previously unrecognized signaling pathway resulting in a strong increase of its protein level in human hepatocytes in culture. 5-AzaC caused a sustained activation of the inositol-requiring enzyme 1α (IRE1α) kinase domain and c-Jun N-terminal kinase (JNK) independently of endoplasmic reticulum stress. This resulted in activation of the epidermal growth factor receptor (EGFR) and extracellular signal-regulated kinase1/2 (ERK1/2) that, in turn, stabilized LDLR mRNA. Systematic mutation of the AREs (ARE1-3) in the LDLR 3'UTR and expression of each mutant coupled to a luciferase reporter in Huh7 cells demonstrated that ARE1 is required for rapid LDLR mRNA decay and 5-AzaC-induced mRNA stabilization via the IRE1α-EGFR-ERK1/2 signaling cascade. The characterization of this pathway will help to reveal potential targets to enhance plasma LDL clearance and novel cholesterol-lowering therapeutic strategies.

Chromatin accessibility dynamics reveal novel functional enhancers in C. elegans

Chromatin accessibility, a crucial component of genome regulation, has primarily been studied in homogeneous and simple systems, such as isolated cell populations or early-development models. Whether chromatin accessibility can be assessed in complex, dynamic systems in vivo with high sensitivity remains largely unexplored. In this study, we use ATAC-seq to identify chromatin accessibility changes in a whole animal, the model organism Caenorhabditis elegans, from embryogenesis to adulthood. Chromatin accessibility changes between developmental stages are highly reproducible, recapitulate histone modification changes, and reveal key regulatory aspects of the epigenomic landscape throughout organismal development. We find that over 5000 distal noncoding regions exhibit dynamic changes in chromatin accessibility between developmental stages and could thereby represent putative enhancers. When tested in vivo, several of these putative enhancers indeed drive novel cell-type- and temporal-specific patterns of expression. Finally, by integrating transcription factor binding motifs in a machine learning framework, we identify EOR-1 as a unique transcription factor that may regulate chromatin dynamics during development. Our study provides a unique resource for C. elegans, a system in which the prevalence and importance of enhancers remains poorly characterized, and demonstrates the power of using whole organism chromatin accessibility to identify novel regulatory regions in complex systems.

Hypoxic Induction of Vasorin Regulates Notch1 Turnover to Maintain Glioma Stem-like Cells

Tumor hypoxia is associated with poor patient survival and is a characteristic of glioblastoma. Notch signaling is implicated in maintaining glioma stem-like cells (GSCs) within the hypoxic niche, although the molecular mechanisms linking hypoxia to Notch activation have not been clearly delineated. Here we show that Vasorin is a critical link between hypoxia and Notch signaling in GSCs. Vasorin is preferentially induced in GSCs by a HIF1α/STAT3 co-activator complex and stabilizes Notch1 protein at the cell membrane. This interaction prevents Numb from binding Notch1, rescuing it from Numb-mediated lysosomal degradation. Thus, Vasorin acts as a switch to augment Notch signaling under hypoxic conditions. Vasorin promotes tumor growth and reduces survival in mouse models of glioblastoma, and its expression correlates with increased aggression of human gliomas. These findings provide mechanistic insights into how hypoxia promotes Notch signaling in glioma and identify Vasorin as a potential therapeutic target.

Methylomic predictors demonstrate the role of NF-κB in old-age mortality and are unrelated to the aging-associated epigenetic drift

Changes in the DNA methylation (DNAm) landscape have been implicated in aging and cellular senescence. To unravel the role of specific DNAm patterns in late-life survival, we performed genome-wide methylation profiling in nonagenarians (n=111) and determined the performance of the methylomic predictors and conventional risk markers in a longitudinal setting. The survival model containing only the methylomic markers was superior in terms of predictive accuracy compared with the model containing only the conventional predictors or the model containing conventional predictors combined with the methylomic markers. At the 2.55-year follow-up, we identified 19 mortality-associated (false-discovery rate <0.5) CpG sites that mapped to genes functionally clustering around the nuclear factor kappa B (NF-κB) complex. Interestingly, none of the mortality-associated CpG sites overlapped with the established aging-associated DNAm sites. Our results are in line with previous findings on the role of NF-κB in controlling animal life spans and demonstrate the role of this complex in human longevity.

Chromatin accessibility dynamics reveal novel functional enhancers in C. elegans

Chromatin accessibility, a crucial component of genome regulation, has primarily been studied in homogeneous and simple systems, such as isolated cell populations or early-development models. Whether chromatin accessibility can be assessed in complex, dynamic systems in vivo with high sensitivity remains largely unexplored. In this study, we use ATAC-seq to identify chromatin accessibility changes in a whole animal, the model organism Caenorhabditis elegans, from embryogenesis to adulthood. Chromatin accessibility changes between developmental stages are highly reproducible, recapitulate histone modification changes, and reveal key regulatory aspects of the epigenomic landscape throughout organismal development. We find that over 5000 distal noncoding regions exhibit dynamic changes in chromatin accessibility between developmental stages and could thereby represent putative enhancers. When tested in vivo, several of these putative enhancers indeed drive novel cell-type- and temporal-specific patterns of expression. Finally, by integrating transcription factor binding motifs in a machine learning framework, we identify EOR-1 as a unique transcription factor that may regulate chromatin dynamics during development. Our study provides a unique resource for C. elegans, a system in which the prevalence and importance of enhancers remains poorly characterized, and demonstrates the power of using whole organism chromatin accessibility to identify novel regulatory regions in complex systems.

MYC drives overexpression of telomerase RNA (hTR/TERC) in prostate cancer

Telomerase consists of at least two essential elements, an RNA component hTR or TERC that contains the template for telomere DNA addition and a catalytic reverse transcriptase (TERT). While expression of TERT has been considered the key rate-limiting component for telomerase activity, increasing evidence suggests an important role for the regulation of TERC in telomere maintenance and perhaps other functions in human cancer. By using three orthogonal methods including RNAseq, RT-qPCR, and an analytically validated chromogenic RNA in situ hybridization assay, we report consistent overexpression of TERC in prostate cancer. This overexpression occurs at the precursor stage (e.g. high-grade prostatic intraepithelial neoplasia or PIN) and persists throughout all stages of disease progression. Levels of TERC correlate with levels of MYC (a known driver of prostate cancer) in clinical samples and we also show the following: forced reductions of MYC result in decreased TERC levels in eight cancer cell lines (prostate, lung, breast, and colorectal); forced overexpression of MYC in PCa cell lines, and in the mouse prostate, results in increased TERC levels; human TERC promoter activity is decreased after MYC silencing; and MYC occupies the TERC locus as assessed by chromatin immunoprecipitation (ChIP). Finally, we show that knockdown of TERC by siRNA results in reduced proliferation of prostate cancer cell lines. These studies indicate that TERC is consistently overexpressed in all stages of prostatic adenocarcinoma and that its expression is regulated by MYC. These findings nominate TERC as a novel prostate cancer biomarker and therapeutic target. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Chromatin-state discovery and genome annotation with ChromHMM

Noncoding DNA regions have central roles in human biology, evolution, and disease. ChromHMM helps to annotate the noncoding genome using epigenomic information across one or multiple cell types. It combines multiple genome-wide epigenomic maps, and uses combinatorial and spatial mark patterns to infer a complete annotation for each cell type. ChromHMM learns chromatin-state signatures using a multivariate hidden Markov model (HMM) that explicitly models the combinatorial presence or absence of each mark. ChromHMM uses these signatures to generate a genome-wide annotation for each cell type by calculating the most probable state for each genomic segment. ChromHMM provides an automated enrichment analysis of the resulting annotations to facilitate the functional interpretations of each chromatin state. ChromHMM is distinguished by its modeling emphasis on combinations of marks, its tight integration with downstream functional enrichment analyses, its speed, and its ease of use. Chromatin states are learned, annotations are produced, and enrichments are computed within 1 d.

Beyond genome-wide scan: Association of a cis-regulatory NCR3 variant with mild malaria in a population living in the Republic of Congo

Linkage studies have revealed a linkage of mild malaria to chromosome 6p21 that contains the NCR3 gene encoding a natural killer cell receptor, whereas NCR3-412G>C (rs2736191) located in its promoter region was found to be associated with malaria in Burkina Faso. Here we confirmed the association of rs2736191 with mild malaria in a Congolese cohort and investigated its potential cis-regulatory effect. Luciferase assay results indicated that rs2736191-G allele had a significantly increased promoter activity compared to rs2736191-C allele. Furthermore, EMSAs demonstrated an altered binding of two nuclear protein complexes to the rs2736191-C allele in comparison to rs2736191-G allele. Finally, after in silico identification of transcription factor candidates, pull-down western blot experiments confirmed that both STAT4 and RUNX3 bind the region encompassing rs2736191 with a higher affinity for the G allele. To our knowledge, this is the first report that explored the functional role of rs2736191. These results support the hypothesis that genetic variation within natural killer cell receptors alters malaria resistance in humans.

Coordinated increase of nuclear tension and lamin-A with matrix stiffness outcompetes lamin-B receptor that favors soft tissue phenotypes

Matrix stiffness that is sensed by a cell or measured by a purely physical probe reflects the intrinsic elasticity of the matrix and also how thick or thin the matrix is. Here, mesenchymal stem cells (MSCs) and their nuclei spread in response to thickness-corrected matrix microelasticity, with increases in nuclear tension and nuclear stiffness resulting from increases in myosin-II and lamin-A,C. Linearity between the widely varying projected area of a cell and its nucleus across many matrices, timescales, and myosin-II activity levels indicates a constant ratio of nucleus-to-cell volume, despite MSCs' lineage plasticity. Nuclear envelope fluctuations are suppressed on the stiffest matrices, and fluctuation spectra reveal a high nuclear tension that matches trends from traction force microscopy and from increased lamin-A,C. Transcriptomes of many diverse tissues and MSCs further show that lamin-A,C's increase with tissue or matrix stiffness anti-correlates with lamin-B receptor (LBR), which contributes to lipid/sterol biosynthesis. Adipogenesis (a soft lineage) indeed increases LBR:lamin-A,C protein stoichiometry in MSCs versus osteogenesis (stiff). The two factors compete for lamin-B in response to matrix elasticity, knockdown, myosin-II inhibition, and even constricted migration that disrupts and segregates lamins in situ. Matrix stiffness-driven contractility thus tenses the nucleus to favor lamin-A,C accumulation and suppress soft tissue phenotypes.