ENCODE data in the UCSC Genome Browser: year 5 update

Evaluating Chromatin Accessibility Differences Across Multiple Primate Species Using a Joint Modeling Approach

Changes in transcriptional regulation are thought to be a major contributor to the evolution of phenotypic traits, but the contribution of changes in chromatin accessibility to the evolution of gene expression remains almost entirely unknown. To address this important gap in knowledge, we developed a new method to identify DNase I Hypersensitive (DHS) sites with differential chromatin accessibility between species using a joint modeling approach. Our method overcomes several limitations inherent to conventional threshold-based pairwise comparisons that become increasingly apparent as the number of species analyzed rises. Our approach employs a single quantitative test which is more sensitive than existing pairwise methods. To illustrate, we applied our joint approach to DHS sites in fibroblast cells from five primates (human, chimpanzee, gorilla, orangutan, and rhesus macaque). We identified 89,744 DHS sites, of which 41% are identified as differential between species using the joint model compared with 33% using the conventional pairwise approach. The joint model provides a principled approach to distinguishing single from multiple chromatin accessibility changes among species. We found that nondifferential DHS sites are enriched for nucleotide conservation. Differential DHS sites with decreased chromatin accessibility relative to rhesus macaque occur more commonly near transcription start sites (TSS), while those with increased chromatin accessibility occur more commonly distal to TSS. Further, differential DHS sites near TSS are less cell type-specific than more distal regulatory elements. Taken together, these results point to distinct classes of DHS sites, each with distinct characteristics of selection, genomic location, and cell type specificity.

Differentially methylated regions in bipolar disorder and suicide

The addition of a methyl group to, typically, a cytosine-guanine dinucleotide (CpG) creates distinct DNA methylation patterns across the genome that can regulate gene expression. Aberrant DNA methylation of CpG sites has been associated with many psychiatric disorders including bipolar disorder (BD) and suicide. Using the SureSelect^XT system, Methyl-Seq, we investigated the DNA methylation status of CpG sites throughout the genome in 50 BD individuals (23 subjects who died by suicide and 27 subjects who died from other causes) and 31 nonpsychiatric controls. We identified differentially methylated regions (DMRs) from three analyses: (a) BD subjects compared to nonpsychiatric controls (BD-NC), (b) BD subjects who died by suicide compared to nonpsychiatric controls (BDS-NC), and (c) BDS subjects compared to BD subjects who died from other causes (BDS-BDNS). One DMR from the BDS-NC analysis, located in ARHGEF38, was significantly hypomethylated (23.4%) in BDS subjects. This finding remained significant after multiple testing (P_Bootstrapped = 9.0 × 10^-3 ), was validated using pyrosequencing, and was more significant in males. A secondary analysis utilized Ingenuity Pathway Analysis to identify enrichment in nominally significant DMRs. This identified an association with several pathways including axonal guidance signaling, calcium signaling, β-adrenergic signaling, and opioid signaling. Our comprehensive study provides further support that DNA methylation alterations influence the risk for BD and suicide. However, further investigation is required to confirm these associations and identify their functional consequences.

Developmentally regulated Shh expression is robust to TAD perturbations

Topologically associating domains (TADs) have been proposed to both guide and constrain enhancer activity. Shh is located within a TAD known to contain all its enhancers. To investigate the importance of chromatin conformation and TAD integrity on developmental gene regulation, we have manipulated the Shh TAD - creating internal deletions, deleting CTCF sites, and deleting and inverting sequences at TAD boundaries. Chromosome conformation capture and fluorescence in situ hybridisation assays were used to investigate the changes in chromatin conformation that result from these manipulations. Our data suggest that these substantial alterations in TAD structure have no readily detectable effect on Shh expression patterns or levels of Shh expression during development - except where enhancers are deleted - and result in no detectable phenotypes. Only in the case of a larger deletion at one TAD boundary could ectopic influence of the Shh limb enhancer be detected on a gene (Mnx1) in the neighbouring TAD. Our data suggests that, contrary to expectations, the developmental regulation of Shh expression is remarkably robust to TAD perturbations.

A variant of the castor zinc finger 1 (CASZ1) gene is differentially associated with the clinical classification of chronic venous disease

Recent reports have suggested a reproducible association between the rs11121615 SNP, located within an intron of the castor zinc finger 1 (CASZ1) gene, and varicose veins. This study aimed to determine if this variant is also differentially associated with the various clinical classifications of chronic venous disease (CVD). The rs11121615 SNP was genotyped in two independent cohorts from New Zealand (n = 1876 controls /1606 CVD cases) and the Netherlands (n = 1626/2966). Participants were clinically assessed using well-established CVD criteria. The association between the rs11121615 C-allele and varicose veins was validated in both cohorts. This was strongest in those with higher clinical severity classes and was not significant in those with non-varicose vein CVD. Functional analysis of the rs11121615 variant demonstrated that the risk allele was associated with increased enhancer activity. This study demonstrates that the CASZ1 gene associated C-allele of rs11121615 has a significant, reproducible, association with CVD (CEAP C ≥ 2 meta-odds ratio 1.31, 95% CI 1.27–1.34, P = 1 × 10⁻⁹⁸, PHet = 0.25), but not with non-varicose vein (CEAP C1, telangiectasia or reticular veins) forms of venous disease. The effect size of this association therefore appears to be susceptible to influence by phenotypic heterogeneity, particularly if a cohort includes a large number of cases with lower severity CVD.

DNA methylation loss promotes immune evasion of tumours with high mutation and copy number load

Mitotic cell division increases tumour mutation burden and copy number load, predictive markers of the clinical benefit of immunotherapy. Cell division correlates also with genomic demethylation involving methylation loss in late-replicating partial methylation domains. Here we find that immunomodulatory pathway genes are concentrated in these domains and transcriptionally repressed in demethylated tumours with CpG island promoter hypermethylation. Global methylation loss correlated with immune evasion signatures independently of mutation burden and aneuploidy. Methylome data of our cohort (n = 60) and a published cohort (n = 81) in lung cancer and a melanoma cohort (n = 40) consistently demonstrated that genomic methylation alterations counteract the contribution of high mutation burden and increase immunotherapeutic resistance. Higher predictive power was observed for methylation loss than mutation burden. We also found that genomic hypomethylation correlates with the immune escape signatures of aneuploid tumours. Hence, DNA methylation alterations implicate epigenetic modulation in precision immunotherapy.

Demethylation of the genome is found in cancer. Here, the authors show that genomic demethylation entails changes in promoter methylation and gene expression associated with immune escape and suggest that the epigenetic alterations may be an important determinant of responses to immunotherapy.

Reconstruction of the lncRNA-miRNA-mRNA network based on competitive endogenous RNA reveal functional lncRNAs in Cerebral Infarction

Functioning as miRNA sponges, long non-coding RNA (lncRNA) exert its pharmacological action via regulating expression of protein-coding genes. However, the lncRNA-mediated ceRNA in cerebral Infarction (CI) remains unclear. In this study, the expression recordsets of mRNA, lncRNA and miRNA of CI samples were obtained from the NCBI GEO datasets separately. The differentially expressed lncRNAs (DELs), miRNAs (DEMis) and mRNAs (DEMs) were identified by limma package in R platform. A total of 267 DELs, 26 DEMis, and 760 DEMs were identified as differentially expressed profiles, with which we constructed the ceRNA network composed of DELs-DEMis-DEMs. Further, clusterProfiler package in R platform is employed for performing Gene Ontology (GO) and KEGG pathway analysis. An aberrant ceRNA network was constructed according to node degrees in CI, including 28 DELs, 19 DEMs and 12 DEMis, from which we extracted the core network, in which 9 nodes were recognized as kernel genes including Tspan3, Eif4a2, rno-miR-208a-3p, rno-miR-194-5p, Pdpn, H3f3b, Stat3, Cd63 and Sdc4. Finally, with the DELs-DEMis-DEMs ceRNA network provided above, we can improve our understanding of the pathogenesis of CI mediated by lncRNA.

CRISPR/Cas9-mediated knock-in of an optimized TetO repeat for live cell imaging of endogenous loci

Nuclear organization has an important role in determining genome function; however, it is not clear how spatiotemporal organization of the genome relates to functionality. To elucidate this relationship, a method for tracking any locus of interest is desirable. Recently clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) or transcription activator-like effectors were adapted for imaging endogenous loci; however, they are mostly limited to visualization of repetitive regions. Here, we report an efficient and scalable method named SHACKTeR (Short Homology and CRISPR/Cas9-mediated Knock-in of a TetO Repeat) for live cell imaging of specific chromosomal regions without the need for a pre-existing repetitive sequence. SHACKTeR requires only two modifications to the genome: CRISPR/Cas9-mediated knock-in of an optimized TetO repeat and its visualization by TetR-EGFP expression. Our simplified knock-in protocol, utilizing short homology arms integrated by polymerase chain reaction, was successful at labeling 10 different loci in HCT116 cells. We also showed the feasibility of knock-in into lamina-associated, heterochromatin regions, demonstrating that these regions prefer non-homologous end joining for knock-in. Using SHACKTeR, we were able to observe DNA replication at a specific locus by long-term live cell imaging. We anticipate the general applicability and scalability of our method will enhance causative analyses between gene function and compartmentalization in a high-throughput manner.

Mechanistic interpretation of non-coding variants for discovering transcriptional regulators of drug response

BACKGROUND

Identification of functional non-coding variants and their mechanistic interpretation is a major challenge of modern genomics, especially for precision medicine. Transcription factor (TF) binding profiles and epigenomic landscapes in reference samples allow functional annotation of the genome, but do not provide ready answers regarding the effects of non-coding variants on phenotypes. A promising computational approach is to build models that predict TF-DNA binding from sequence, and use such models to score a variant's impact on TF binding strength. Here, we asked if this mechanistic approach to variant interpretation can be combined with information on genotype-phenotype associations to discover transcription factors regulating phenotypic variation among individuals.

RESULTS

We developed a statistical approach that integrates phenotype, genotype, gene expression, TF ChIP-seq, and Hi-C chromatin interaction data to answer this question. Using drug sensitivity of lymphoblastoid cell lines as the phenotype of interest, we tested if non-coding variants statistically linked to the phenotype are enriched for strong predicted impact on DNA binding strength of a TF and thus identified TFs regulating individual differences in the phenotype. Our approach relies on a new method for predicting variant impact on TF-DNA binding that uses a combination of biophysical modeling and machine learning. We report statistical and literature-based support for many of the TFs discovered here as regulators of drug response variation. We show that the use of mechanistically driven variant impact predictors can identify TF-drug associations that would otherwise be missed. We examined in depth one reported association-that of the transcription factor ELF1 with the drug doxorubicin-and identified several genes that may mediate this regulatory relationship.

CONCLUSION

Our work represents initial steps in utilizing predictions of variant impact on TF binding sites for discovery of regulatory mechanisms underlying phenotypic variation. Future advances on this topic will be greatly beneficial to the reconstruction of phenotype-associated gene regulatory networks.

Stable enhancers are active in development, and fragile enhancers are associated with evolutionary adaptation

Background

Despite continual progress in the identification and characterization of trait- and disease-associated variants that disrupt transcription factor (TF)-DNA binding, little is known about the distribution of TF binding deactivating mutations (deMs) in enhancer sequences. Here, we focus on elucidating the mechanism underlying the different densities of deMs in human enhancers.

Results

We identify two classes of enhancers based on the density of nucleotides prone to deMs. Firstly, fragile enhancers with abundant deM nucleotides are associated with the immune system and regular cellular maintenance. Secondly, stable enhancers with only a few deM nucleotides are associated with the development and regulation of TFs and are evolutionarily conserved. These two classes of enhancers feature different regulatory programs: the binding sites of pioneer TFs of FOX family are specifically enriched in stable enhancers, while tissue-specific TFs are enriched in fragile enhancers. Moreover, stable enhancers are more tolerant of deMs due to their dominant employment of homotypic TF binding site (TFBS) clusters, as opposed to the larger-extent usage of heterotypic TFBS clusters in fragile enhancers. Notably, the sequence environment and chromatin context of the cognate motif, other than the motif itself, contribute more to the susceptibility to deMs of TF binding.

Conclusions

This dichotomy of enhancer activity is conserved across different tissues, has a specific footprint in epigenetic profiles, and argues for a bimodal evolution of gene regulatory programs in vertebrates. Specifically encoded stable enhancers are evolutionarily conserved and associated with development, while differently encoded fragile enhancers are associated with the adaptation of species.

Electronic supplementary material

The online version of this article (10.1186/s13059-019-1750-z) contains supplementary material, which is available to authorized users.

Chronic activation of FXR-induced liver growth with tissue-specific targeting Cyclin D1

The nuclear receptor (FXR) plays essential roles in maintaining bile acid and lipid homeostasis by regulating diverse target genes. And its agonists were promising agents for treating various liver diseases. Nevertheless, the potential side effect of chronic FXR activation by specific agonists is not fully understood. In this study, we investigated the mechanism of FXR agonist WAY-362450 induced liver enlargement during treating liver diseases. We demonstrated that chronic ingestion of WAY-362450 induced liver hypertrophy instead of hyperplasia in mouse. Global transcriptional pattern was also examined in mouse livers after treatment with WAY-362450 by RNA-seq assay. Through GO and KEGG enrichment analyses, we demonstrated that the expression of Cyclin D1 (Ccnd1) among the cell cycle-regulating genes was notably increased in WAY-362450-treated mouse liver. Activation of FXR-induced Ccnd1 expression in hepatocyte in a time-dependent manner in vivo and in vitro. Through bioinformatics analysis and ChIP assay, we identified FXR as a direct transcriptional activator of Ccnd1 through binding to a potential enhancer, which was specifically active in livers. We also found active histone acetylation was essential for Ccnd1 induction by FXR. Thus, our study indicated that activation of FXR-induced harmless liver hypertrophy with spatiotemporal modulation of Ccnd1. With a better understanding of the mechanism of tissue-specific gene regulation by FXR, it is beneficial for development and appropriate application of its specific agonist in preventing hepatic diseases.

Nanofiber-expanded human CD34+ cells heal cutaneous wounds in streptozotocin-induced diabetic mice

Despite advances in diabetic wound care, the significant number of amputations that occur every year demands more effective therapeutics. Herein, we offer an aminated polyethersulfone nanofiber-expanded human umbilical cord blood-derived CD34⁺ cells (henceforth CD34⁺ cells) effective therapy, tested in cutaneous wounds developed in streptozotocin-induced diabetic NOD/SCID mice. We show that systemic administration of CD34⁺ cells homed to the wound site and significantly accelerated wound closure. Wound closure was associated with improved re-epithelialization and increased neovascularization; and with decreased sustained pro-inflammatory activity of NF-κB and its downstream effector molecules TNF-α, IL-1β, and IL-6 at the wound bed. This finding was further supported by the observation of a decreased number of myeloperoxidase positive neutrophils, and concomitantly increased levels of IL-10. In addition, improved granulation tissue formation was observed along with higher collagen deposition and myofibroblasts and decreased expressions of MMP-1. Mechanistically, CD34⁺ cells reduced the level of MMP-1 expression by inhibiting recruitment of NF-κB to the MMP-1 promoter site in dermal fibroblasts. In summary, we provide evidence of a novel nanofiber-expanded CD34⁺ stem cell therapeutic development for treating diabetic wounds by defining their cellular and molecular mechanisms.

Haplotype-aware diplotyping from noisy long reads

Current genotyping approaches for single-nucleotide variations rely on short, accurate reads from second-generation sequencing devices. Presently, third-generation sequencing platforms are rapidly becoming more widespread, yet approaches for leveraging their long but error-prone reads for genotyping are lacking. Here, we introduce a novel statistical framework for the joint inference of haplotypes and genotypes from noisy long reads, which we term diplotyping. Our technique takes full advantage of linkage information provided by long reads. We validate hundreds of thousands of candidate variants that have not yet been included in the high-confidence reference set of the Genome-in-a-Bottle effort.

De Deyn P, C. van Swieten J, M. van Duijn C, van der Zee J, Sleegers K, Van Broeckhoven C. Loss of DPP6 in neurodegenerative dementia: a genetic player in the dysfunction of neuronal excitability

Emerging evidence suggested a converging mechanism in neurodegenerative brain diseases (NBD) involving early neuronal network dysfunctions and alterations in the homeostasis of neuronal firing as culprits of neurodegeneration. In this study, we used paired-end short-read and direct long-read whole genome sequencing to investigate an unresolved autosomal dominant dementia family significantly linked to 7q36. We identified and validated a chromosomal inversion of ca. 4 Mb, segregating on the disease haplotype and disrupting the coding sequence of dipeptidyl-peptidase 6 gene (DPP6). DPP6 resequencing identified significantly more rare variants-nonsense, frameshift, and missense-in early-onset Alzheimer's disease (EOAD, p value = 0.03, OR = 2.21 95% CI 1.05-4.82) and frontotemporal dementia (FTD, p = 0.006, OR = 2.59, 95% CI 1.28-5.49) patient cohorts. DPP6 is a type II transmembrane protein with a highly structured extracellular domain and is mainly expressed in brain, where it binds to the potassium channel Kv4.2 enhancing its expression, regulating its gating properties and controlling the dendritic excitability of hippocampal neurons. Using in vitro modeling, we showed that the missense variants found in patients destabilize DPP6 and reduce its membrane expression (p < 0.001 and p < 0.0001) leading to a loss of protein. Reduced DPP6 and/or Kv4.2 expression was also detected in brain tissue of missense variant carriers. Loss of DPP6 is known to cause neuronal hyperexcitability and behavioral alterations in Dpp6-KO mice. Taken together, the results of our genomic, genetic, expression and modeling analyses, provided direct evidence supporting the involvement of DPP6 loss in dementia. We propose that loss of function variants have a higher penetrance and disease impact, whereas the missense variants have a variable risk contribution to disease that can vary from high to low penetrance. Our findings of DPP6, as novel gene in dementia, strengthen the involvement of neuronal hyperexcitability and alteration in the homeostasis of neuronal firing as a disease mechanism to further investigate.

SP140 regulates the expression of immune-related genes associated with multiple sclerosis and other autoimmune diseases by NF-κB inhibition

SP140 locus has been associated with multiple sclerosis (MS) as well as other autoimmune diseases by genome-wide association studies (GWAS). The causal variant of these associations (rs28445040-T) alters the splicing of the SP140 gene transcripts reducing the protein expression. We aimed to understand why the reduction of SP140 expression produced by the risk variant can increase the susceptibility to MS. To this end, we determined by RNA sequencing (RNA-seq) analysis the differentially expressed genes after SP140 silencing in lymphoblastoid cell lines (LCLs). We analyzed these genes by gene ontology (GO), comparative transcriptome profiles, enrichment of transcription factors (TFs) in the promoters of these genes and colocalization with GWAS risk variants. We also monitored the activity of the nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) in SP140-silenced cells by luciferase reporter system. We identified 100 genes that were up-regulated and 22 genes down-regulated in SP140-silenced LCLs. GO analysis revealed that genes affected by SP140 were involved in regulation of cytokine production, inflammatory response and cell-cell adhesion. We observed enrichment of NF-κB TF in the promoter of up-regulated genes and NF-κB-increased activity in SP140-silenced cell lines. We showed enrichment of genes regulated by SP140 in GWAS-detected risk loci for MS (14.63 folds), Crohn's disease (4.82 folds) and inflammatory bowel disease (4.47 folds), not observed in other unrelated immune diseases. Our findings showed that SP140 is an important repressor of genes implicated in inflammation, suggesting that decreased expression of SP140, promoted by the rs28445040-T risk variant, may lead to up-regulation of these genes by means of NF-κB inhibition in B cells.

Hierarchical cooperation of transcription factors from integration analysis of DNA sequences, ChIP-Seq and ChIA-PET data

Background

Chromosomal architecture, which is constituted by chromatin loops, plays an important role in cellular functions. Gene expression and cell identity can be regulated by the chromatin loop, which is formed by proximal or distal enhancers and promoters in linear DNA (1D). Enhancers and promoters are fundamental non-coding elements enriched with transcription factors (TFs) to form chromatin loops. However, the specific cooperation of TFs involved in forming chromatin loops is not fully understood.

Results

Here, we proposed a method for investigating the cooperation of TFs in four cell lines by the integrative analysis of DNA sequences, ChIP-Seq and ChIA-PET data. Results demonstrate that the interaction of enhancers and promoters is a hierarchical and dynamic complex process with cooperative interactions of different TFs synergistically regulating gene expression and chromatin structure. The TF cooperation involved in maintaining and regulating the chromatin loop of cells can be regulated by epigenetic factors, such as other TFs and DNA methylation.

Conclusions

Such cooperation among TFs provides the potential features that can affect chromatin’s 3D architecture in cells. The regulation of chromatin 3D organization and gene expression is a complex process associated with the hierarchical and dynamic prosperities of TFs.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5535-2) contains supplementary material, which is available to authorized users.

Genetic variants in the circadian rhythm pathway as indicators of prostate cancer progression

Background

To determine the association between circadian pathway genetic variants and the risk of prostate cancer progression.

Methods

We systematically evaluated 79 germline variants in nine circadian pathway genes in a cohort of 458 patients with localized prostate cancer as the discovery phase. We then replicated the significant findings in another cohort of 324 men with more advanced disease. The association of each variant with prostate cancer progression was evaluated by a log-rank test and Cox regression.

Results

A single nucleotide polymorphism of the neuronal PAS domain protein 2 (NPAS2) gene (rs6542993 A>T) was found to be associated with a significantly higher risk of disease progression in both localized (P = 0.001) and advanced (P = 0.039) prostate cancer cases. In silico analysis revealed decreased expression levels of NPAS2 in carriers of the T allele of rs6542993 compared with those carrying the A allele. Consistently, downregulation of NPAS2 expression was associated with more aggressive prostate cancer and poor progression-free survival (log-rank P = 0.002).

Conclusions

The NPAS2 rs6542993 polymorphism may be a promising biomarker, and may shed light on the pathways that govern prostate cancer progression.

In vivo Application of the REMOTE-control System for the Manipulation of Endogenous Gene Expression

Here we describe a protocol for implementing the REMOTE-control system (Reversible Manipulation of Transcription at Endogenous loci), which allows for reversible and tunable expression control of an endogenous gene of interest in living model systems. The REMOTE-control system employs enhanced lac repression and tet activation systems to achieve down- or upregulation of a target gene within a single biological system. Tight repression can be achieved from repressor binding sites flexibly located far downstream of a transcription start site by inhibiting transcription elongation. Robust upregulation can be attained by enhancing the transcription of an endogenous gene by targeting tet transcriptional activators to the cognate promoter. This reversible and tunable expression control can be applied and withdrawn repeatedly in organisms. The potency and versatility of the system, as demonstrated for endogenous Dnmt1 here, will allow more precise in vivo functional analyses by enabling investigation of gene function at various expression levels and by testing the reversibility of a phenotype.

High Throughput Sequencing-Based Approaches for Gene Expression Analysis

Next-generation sequencing has emerged as the method of choice to answer fundamental questions in biology. The massively parallel sequencing technology for RNA-Seq analysis enables better understanding of gene expression patterns in model and nonmodel organisms. Sequencing per se has reached the stage of commodity level while analyzing and interpreting huge amount of data has been a significant challenge. This chapter is aimed at discussing the complexities involved in sequencing and analysis, and tries to simplify sequencing based gene expression analysis. Biologists and experimental scientists were kept in mind while discussing the methods and analysis workflow.

Role of Epigenomics in Bone and Cartilage Disease

Phenotypic variation in skeletal traits and diseases is the product of genetic and environmental factors. Epigenetic mechanisms include information-containing factors, other than DNA sequence, that cause stable changes in gene expression and are maintained during cell divisions. They represent a link between environmental influences, genome features, and the resulting phenotype. The main epigenetic factors are DNA methylation, posttranslational changes of histones, and higher-order chromatin structure. Sometimes non-coding RNAs, such as microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), are also included in the broad term of epigenetic factors. There is rapidly expanding experimental evidence for a role of epigenetic factors in the differentiation of bone cells and the pathogenesis of skeletal disorders, such as osteoporosis and osteoarthritis. However, different from genetic factors, epigenetic signatures are cell- and tissue-specific and can change with time. Thus, elucidating their role has particular difficulties, especially in human studies. Nevertheless, epigenomewide association studies are beginning to disclose some disease-specific patterns that help to understand skeletal cell biology and may lead to development of new epigenetic-based biomarkers, as well as new drug targets useful for treating diffuse and localized disorders. Here we provide an overview and update of recent advances on the role of epigenomics in bone and cartilage diseases. © 2019 American Society for Bone and Mineral Research.

Sleeping Beauty Screen Identifies RREB1 and Other Genetic Drivers in Human B-cell Lymphoma

Follicular lymphoma and diffuse large B-cell lymphoma (DLBCL) are the most common non-Hodgkin lymphomas distinguishable by unique mutations, chromosomal rearrangements, and gene expression patterns. Here, it is demonstrated that early B-cell progenitors express 2',3'-cyclic-nucleotide 3' phosphodiesterase (CNP) and that when targeted with Sleeping Beauty (SB) mutagenesis, Trp53R270H mutation or Pten loss gave rise to highly penetrant lymphoid diseases, predominantly follicular lymphoma and DLBCL. In efforts to identify the genetic drivers and signaling pathways that are functionally important in lymphomagenesis, SB transposon insertions were analyzed from splenomegaly specimens of SB-mutagenized mice (n = 23) and SB-mutagenized mice on a Trp53R270H background (n = 7) and identified 48 and 12 sites with statistically recurrent transposon insertion events, respectively. Comparison with human data sets revealed novel and known driver genes for B-cell development, disease, and signaling pathways: PI3K-AKT-mTOR, MAPK, NFκB, and B-cell receptor (BCR). Finally, functional data indicate that modulating Ras-responsive element-binding protein 1 (RREB1) expression in human DLBCL cell lines in vitro alters KRAS expression, signaling, and proliferation; thus, suggesting that this proto-oncogene is a common mechanism of RAS/MAPK hyperactivation in human DLBCL. IMPLICATIONS: A forward genetic screen identified new genetic drivers of human B-cell lymphoma and uncovered a RAS/MAPK-activating mechanism not previously appreciated in human lymphoid disease. Overall, these data support targeting the RAS/MAPK pathway as a viable therapeutic target in a subset of human patients with DLBCL.

Epigenome-wide association study in whole blood on type 2 diabetes among sub-Saharan African individuals: findings from the RODAM study

BACKGROUND

Type 2 diabetes (T2D) results from a complex interplay between genetics and the environment. Several epigenome-wide association studies (EWAS) have found DNA methylation loci associated with T2D in European populations. However, data from African populations are lacking. We undertook the first EWAS for T2D among sub-Saharan Africans, aiming at identifying ubiquitous and novel DNA methylation loci associated with T2D.

METHODS

The Illumina 450k DNA-methylation array was used on whole blood samples of 713 Ghanaian participants (256 with T2D, 457 controls) from the cross-sectional Research on Obesity and Diabetes among African Migrants (RODAM) study. Differentially methylated positions (DMPs) for T2D and HbA1c were identified through linear regression analysis adjusted for age, sex, estimated cell counts, hybridization batch, array position and body mass index (BMI). We also did a candidate analysis of previously reported EWAS loci for T2D in non-African populations, identified through a systematic literature search.

RESULTS

Four DMPs [cg19693031 (TXNIP), cg04816311 (C7orf50), cg00574958 (CPT1A), cg07988171 (TPM4)] were associated with T2D after correction for inflation by possible systematic biases. The most strongly associated DMP-cg19693031, TXNIP (P = 2.6E-19) -showed hypomethylation in T2D cases compared with controls. Two out of the four DMPs [cg19693031 (TXNIP), cg04816311 (C7orf50)] remained associated with T2D after adjustment for BMI, and one locus [cg07988171 (TPM4)] that has not been reported previously.

CONCLUSIONS

In this first EWAS for T2D in sub-Saharan Africans, we have identified four DMPs at epigenome-wide level, one of which is novel. These findings provide insight into the epigenetic loci that underlie the burden of T2D in sub-Saharan Africans.

The Third Intron of IRF8 Is a Cell-Type-Specific Chromatin Priming Element during Mouse Embryonal Stem Cell Differentiation

Interferon regulatory factor 8 (IRF8) is a nuclear transcription factor that plays a key role in the hierarchical differentiation of hematopoietic stem cells toward monocyte/dendritic cell lineages. Therefore, its expression is mainly limited to bone marrow-derived cells. The molecular mechanisms governing this cell-type-restricted expression have been described. However, the molecular mechanisms that are responsible for its silencing in non-hematopoietic cells are elusive. Recently, we demonstrated a role for IRF8 third intron in restricting its expression in non-hematopoietic cells. Furthermore, we showed that this intron alone is sufficient to promote repressed chromatin a cell-type-specific manner. Here we demonstrate the effect of the IRF8 third intron on chromatin conformation during murine embryonal stem cell differentiation. Using genome editing, we provide data showing that the third intron plays a key role in priming the chromatin state of the IRF8 locus during cell differentiation. It mediates dual regulatory effects in a cell-type-specific mode. It acts as a repressor element governing chromatin state of the IRF8 locus during embryonal stem cell differentiation to cardiomyocytes that are expression-restrictive cells. Conversely, it functions as an activator element that is essential for open chromatin structure during the differentiation of these cells to dendritic cells that are expression-permissive cells. Together, these results point to the role of IRF8 third intron as a cell-type-specific chromatin priming element during embryonal stem cell differentiation. These data add another layer to our understanding of the molecular mechanisms governing misexpression of a cell-type-specific gene such as IRF8.

Decreased Nuclear Ascorbate Accumulation Accompanied with Altered Genomic Methylation Pattern in Fibroblasts from Arterial Tortuosity Syndrome Patients

Ascorbate requiring Fe²⁺/2-oxoglutarate-dependent dioxygenases located in the nucleoplasm have been shown to participate in epigenetic regulation of gene expression via histone and DNA demethylation. Transport of dehydroascorbic acid is impaired in the endomembranes of fibroblasts from arterial tortuosity syndrome (ATS) patients, due to the mutation in the gene coding for glucose transporter GLUT10. We hypothesized that altered nuclear ascorbate concentration might be present in ATS fibroblasts, affecting dioxygenase activity and DNA demethylation. Therefore, our aim was to characterize the subcellular distribution of vitamin C, the global and site-specific changes in 5-methylcytosine and 5-hydroxymethylcytosine levels, and the effect of ascorbate supplementation in control and ATS fibroblast cultures. Diminished nuclear accumulation of ascorbate was found in ATS fibroblasts upon ascorbate or dehydroascorbic acid addition. Analyzing DNA samples of cultured fibroblasts from controls and ATS patients, a lower global 5-hydroxymethylcytosine level was found in ATS fibroblasts, which could not be significantly modified by ascorbate addition. Investigation of the (hydroxy)methylation status of specific regions in six candidate genes related to ascorbate metabolism and function showed that ascorbate addition could stimulate hydroxymethylation and active DNA demethylation at the PPAR-γ gene region in control fibroblasts only. The altered DNA hydroxymethylation patterns in patient cells both at the global level and at specific gene regions accompanied with decreased nuclear accumulation of ascorbate suggests the epigenetic role of vitamin C in the pathomechanism of ATS. The present findings represent the first example for the role of vitamin C transport in epigenetic regulation suggesting that ATS is a compartmentalization disease.

A tale of two cities: The genetic mechanisms governing calvarial bone development

The skull bones must grow in a coordinated, three-dimensional manner to coalesce and form the head and face. Mammalian skull bones have a dual embryonic origin from cranial neural crest cells (CNCC) and paraxial mesoderm (PM) and ossify through intramembranous ossification. The calvarial bones, the bones of the cranium which cover the brain, are derived from the supraorbital arch (SOA) region mesenchyme. The SOA is the site of frontal and parietal bone morphogenesis and primary center of ossification. The objective of this review is to frame our current in vivo understanding of the morphogenesis of the calvarial bones and the gene networks regulating calvarial bone initiation in the SOA mesenchyme.

A novel form of JARID2 is required for differentiation in lineage-committed cells

Polycomb repressive complex‐2 (PRC2) is a group of proteins that play an important role during development and in cell differentiation. PRC2 is a histone‐modifying complex that catalyses methylation of lysine 27 of histone H3 (H3K27me3) at differentiation genes leading to their transcriptional repression. JARID2 is a co‐factor of PRC2 and is important for targeting PRC2 to chromatin. Here, we show that, unlike in embryonic stem cells, in lineage‐committed human cells, including human epidermal keratinocytes, JARID2 predominantly exists as a novel low molecular weight form, which lacks the N‐terminal PRC2‐interacting domain (ΔN‐JARID2). We show that ΔN‐JARID2 is a cleaved product of full‐length JARID2 spanning the C‐terminal conserved jumonji domains. JARID2 knockout in keratinocytes results in up‐regulation of cell cycle genes and repression of many epidermal differentiation genes. Surprisingly, repression of epidermal differentiation genes in JARID2‐null keratinocytes can be rescued by expression of ΔN‐JARID2 suggesting that, in contrast to PRC2, ΔN‐JARID2 promotes activation of differentiation genes. We propose that a switch from expression of full‐length JARID2 to ΔN‐JARID2 is important for the up‐regulation differentiation genes.

Genome-wide association study (GWAS) of ovarian cancer in Japanese predicted regulatory variants in 22q13.1

Genome-wide association studies (GWAS) have identified greater than 30 variants associated with ovarian cancer, but most of these variants were investigated in European populations. Here, we integrated GWAS and subsequent functional analyses to identify the genetic variants with potential regulatory effects. We conducted GWAS for ovarian cancer using 681 Japanese cases and 17,492 controls and found that rs137672 on 22q13.1 exhibited a strong association with a P-value of 1.05 × 10⁻⁷ and an odds ratio of 0.573 with a 95% confidence interval of 0.466–0.703. In addition, three previously reported SNPs, i.e., rs10088218, rs9870207 and rs1400482, were validated in the Japanese population (P < 0.05) with the same risk allele as noted in previous studies. Functional studies including regulatory feature analysis and electrophoretic mobility shift assay (EMSA) revealed two regulatory SNPs in 22q13.1, rs2072872 and rs6509, that affect the binding affinity to some nuclear proteins in ovarian cancer cells. The plausible regulatory proteins whose motifs could be affected by the allele changes of these two SNPs were also proposed. Moreover, the protective G allele of rs6509 was associated with a decreased SYNGR1 expression level in normal ovarian tissues. Our findings elucidated the regulatory variants in 22q13.1 that are associated with ovarian cancer risk.

Significantly altered peripheral blood cell DNA methylation profile as a result of immediate effect of metformin use in healthy individuals

BACKGROUND

Metformin is a widely prescribed antihyperglycemic agent that has been also associated with multiple therapeutic effects in various diseases, including several types of malignancies. There is growing evidence regarding the contribution of the epigenetic mechanisms in reaching metformin's therapeutic goals; however, the effect of metformin on human cells in vivo is not comprehensively studied. The aim of our study was to examine metformin-induced alterations of DNA methylation profiles in white blood cells of healthy volunteers, employing a longitudinal study design.

RESULTS

Twelve healthy metformin-naïve individuals where enrolled in the study. Genome-wide DNA methylation pattern was estimated at baseline, 10 h and 7 days after the start of metformin administration. The whole-genome DNA methylation analysis in total revealed 125 differentially methylated CpGs, of which 11 CpGs and their associated genes with the most consistent changes in the DNA methylation profile were selected: POFUT2, CAMKK1, EML3, KIAA1614, UPF1, MUC4, LOC727982, SIX3, ADAM8, SNORD12B, VPS8, and several differentially methylated regions as novel potential epigenetic targets of metformin. The main functions of the majority of top-ranked differentially methylated loci and their representative cell signaling pathways were linked to the well-known metformin therapy targets: regulatory processes of energy homeostasis, inflammatory responses, tumorigenesis, and neurodegenerative diseases.

CONCLUSIONS

Here we demonstrate for the first time the immediate effect of short-term metformin administration at therapeutic doses on epigenetic regulation in human white blood cells. These findings suggest the DNA methylation process as one of the mechanisms involved in the action of metformin, thereby revealing novel targets and directions of the molecular mechanisms underlying the various beneficial effects of metformin.

TRIAL REGISTRATION

EU Clinical Trials Register, 2016-001092-74. Registered 23 March 2017, https://www.clinicaltrialsregister.eu/ctr-search/trial/2016-001092-74/LV .

YY1 negatively regulates the XAF1 gene expression in prostate cancer

XAF1 is a tumor suppressor gene with low or absent expression in cancer. Since transcriptional reactivation or ectopic-mediated expression of XAF1 inhibits tumor growth, it is of great interest to elucidate the molecular mechanisms leading to XAF1 silencing. YY1 is a transcription factor that acts as a repressor or an activator to modulate several cancer-associated cellular processes. Both YY1 and XAF1 have key roles in prostate cancer (PCa) progression and are associated with worse clinical outcomes. To assess whether YY1 regulates the transcriptional activation of the XAF1 gene, we performed gene-reporter assays coupled with site-directed mutagenesis, which showed that YY1 is able to mediate XAF1 silencing. Concordantly, ChIP-qPCR assays showed that YY1 interacts with the XAF1 promoter in PC3 cells that lacks XAF1 expression. This association was lost after exposure to epigenetic modulators that induce XAF1 expression. Further supporting the YY1's repressive role, we found transcriptional reactivation of the XAF1 gene by YY1 downregulation. As expected by previous reports showing that HDAC1 is needed for YY1-mediated repressive actions, we observed XAF1 re-expression after either inhibition or downregulation of the HDAC1 gene. Finally, expression data retrieved from the TCGA consortium showed that PCa samples presented lower XAF1 and higher HDAC expression levels than normal tissues. Thus, our results support a model in which YY1 is able to silence tumor suppressor genes such as XAF1 through HDAC1 in PCa.

HDAC is indispensable for IFN-γ-induced B7-H1 expression in gastric cancer

Background

B7 homolog 1 (B7-H1) overexpression on tumor cells is an important mechanism of immune evasion in gastric cancer (GC). Elucidation of the regulation of B7-H1 expression is urgently required to guide B7-H1-targeted cancer therapy. Interferon gamma (IFN-γ) is thought to be the main driving force behind B7-H1 expression, and epigenetic factors including histone acetylation are recently linked to the process. Here, we investigated the potential role of histone deacetylase (HDAC) in IFN-γ-induced B7-H1 expression in GC. The effect of Vorinostat (SAHA), a small molecular inhibitor of HDAC, on tumor growth and B7-H1 expression in a mouse GC model was also evaluated.

Results

RNA-seq data from The Cancer Genome Atlas revealed that expression of B7-H1, HDAC1–3, 6–8, and 10 and SIRT1, 3, 5, and 6 was higher, and expression of HDAC5 and SIRT4 was lower in GC compared to that in normal gastric tissues; that HDAC3 and HDAC1 expression level significantly correlated with B7-H1 in GC with a respective r value of 0.42 (p < 0.001) and 0.21 (p < 0.001). HDAC inhibitor (Trichostatin A, SAHA, and sodium butyrate) pretreatment suppressed IFN-γ-induced B7-H1 expression on HGC-27 cells. HDAC1 and HDAC3 gene knockdown had the same effect. SAHA pretreatment or HDAC knockdown resulted in impaired IFN-γ signaling, demonstrated by the reduction of JAK2, p-JAK1, p-JAK2, and p-STAT1 expression and inefficient STAT1 nuclear translocation. Furthermore, SAHA pretreatment compromised IFN-γ-induced upregulation of histone H3 lysine 9 acetylation level in B7-H1 gene promoter. In the grafted mouse GC model, SAHA treatment suppressed tumor growth, inhibited B7-H1 expression, and elevated the percentage of tumor-infiltrating CD8+ T cells.

Conclusion

HDAC is indispensable for IFN-γ-induced B7-H1 in GC. The study suggests the possibility of targeting B7-H1 using small molecular HDAC inhibitors for cancer treatment.

Electronic supplementary material

The online version of this article (10.1186/s13148-018-0589-6) contains supplementary material, which is available to authorized users.

Glucocorticoid Receptor Binding Inhibits an Intronic IL33 Enhancer and is Disrupted by rs4742170 (T) Allele Associated with Specific Wheezing Phenotype in Early Childhood

Interleukin 33 (IL-33) is a cytokine constitutively expressed by various cells of barrier tissues that contribute to the development of inflammatory immune responses. According to its function as an alarmin secreted by lung and airway epithelium, IL-33 plays a significant role in pathogenesis of allergic disorders. IL-33 is strongly involved in the pathogenesis of asthma, anaphylaxis, allergy and dermatitis, and genetic variations in IL33 locus are associated with increased susceptibility to asthma. Genome-wide association studies have identified risk "T" allele of the single-nucleotide polymorphism rs4742170 located in putative IL33 enhancer area as susceptible variant for development of specific wheezing phenotype in early childhood. Here, we demonstrate that risk "T" rs4742170 allele disrupts binding of glucocorticoid receptor (GR) transcription factor to IL33 putative enhancer. The IL33 promoter/enhancer constructs containing either 4742170 (T) allele or point mutations in the GR-binding site, were significantly more active and did not respond to cortisol in a pulmonary epithelial cell line. At the same time, the constructs containing rs4742170 (C) allele with a functional GR-binding site were less active and further inhibitable by cortisol. The latter effect was GR-dependent as it was completely abolished by GR-specific siRNA. This mechanism may explain the negative effect of the rs4742170 (T) risk allele on the development of wheezing phenotype that strongly correlates with allergic sensitization in childhood.

Molecular analysis of NPAS3 functional domains and variants

Background

NPAS3 encodes a transcription factor which has been associated with multiple human psychiatric and neurodevelopmental disorders. In mice, deletion of Npas3 was found to cause alterations in neurodevelopment, as well as a marked reduction in neurogenesis in the adult mouse hippocampus. This neurogenic deficit, alongside the reduction in cortical interneuron number, likely contributes to the behavioral and cognitive alterations observed in Npas3 knockout mice. Although loss of Npas3 has been found to affect proliferation and apoptosis, the molecular function of NPAS3 is largely uncharacterized outside of predictions based on its high homology to bHLH–PAS transcription factors. Here we set out to characterize NPAS3 as a transcription factor, and to confirm whether NPAS3 acts as predicted for a Class 1 bHLH–PAS family member.

Results

Through these studies we have experimentally demonstrated that NPAS3 behaves as a true transcription factor, capable of gene regulation through direct association with DNA. NPAS3 and ARNT are confirmed to directly interact in human cells through both bHLH and PAS dimerization domains. The C-terminus of NPAS3 was found to contain a functional transactivation domain. Further, the NPAS3::ARNT heterodimer was shown to directly regulate the expression of VGF and TXNIP through binding of their proximal promoters. Finally, we assessed the effects of three human variants of NPAS3 on gene regulatory function and do not observe significant deficits.

Conclusions

NPAS3 is a true transcription factor capable of regulating expression of target genes through their promoters by directly cooperating with ARNT. The tested human variants of NPAS3 require further characterization to identify their effects on NPAS3 expression and function in the individuals that carry them. These data enhance our understanding of the molecular function of NPAS3 and the mechanism by which it contributes to normal and abnormal neurodevelopment and neural function.

Electronic supplementary material

The online version of this article (10.1186/s12867-018-0117-4) contains supplementary material, which is available to authorized users.

Antiviral activity of bone morphogenetic proteins and activins

Understanding the control of viral infections is of broad importance. Chronic hepatitis C virus (HCV) infection causes decreased expression of the iron hormone hepcidin, which is regulated by hepatic bone morphogenetic protein (BMP)/SMAD signalling. We found that HCV infection and the BMP/SMAD pathway are mutually antagonistic. HCV blunted induction of hepcidin expression by BMP6, probably via tumour necrosis factor (TNF)-mediated downregulation of the BMP co-receptor haemojuvelin. In HCV-infected patients, disruption of the BMP6/hepcidin axis and genetic variation associated with the BMP/SMAD pathway predicted the outcome of infection, suggesting that BMP/SMAD activity influences antiviral immunity. Correspondingly, BMP6 regulated a gene repertoire reminiscent of type I interferon (IFN) signalling, including upregulating interferon regulatory factors (IRFs) and downregulating an inhibitor of IFN signalling, USP18. Moreover, in BMP-stimulated cells, SMAD1 occupied loci across the genome, similar to those bound by IRF1 in IFN-stimulated cells. Functionally, BMP6 enhanced the transcriptional and antiviral response to IFN, but BMP6 and related activin proteins also potently blocked HCV replication independently of IFN. Furthermore, BMP6 and activin A suppressed growth of HBV in cell culture, and activin A inhibited Zika virus replication alone and in combination with IFN. The data establish an unappreciated important role for BMPs and activins in cellular antiviral immunity, which acts independently of, and modulates, IFN.

A novel enhancer near the Pitx1 gene influences development and evolution of pelvic appendages in vertebrates

Vertebrate pelvic reduction is a classic example of repeated evolution. Recurrent loss of pelvic appendages in sticklebacks has previously been linked to natural mutations in a pelvic enhancer that maps upstream of Pitx1. The sequence of this upstream PelA enhancer is not conserved to mammals, so we have surveyed a large region surrounding the mouse Pitx1 gene for other possible hind limb control sequences. Here we identify a new pelvic enhancer, PelB, that maps downstream rather than upstream of Pitx1. PelB drives expression in the posterior portion of the developing hind limb, and deleting the sequence from mice alters the size of several hind limb structures. PelB sequences are broadly conserved from fish to mammals. A wild stickleback population lacking the pelvis has an insertion/deletion mutation that disrupts the structure and function of PelB, suggesting that changes in this ancient enhancer contribute to evolutionary modification of pelvic appendages in nature.

ShapeGTB: the role of local DNA shape in prioritization of functional variants in human promoters with machine learning

Motivation

The identification of functional sequence variations in regulatory DNA regions is one of the major challenges of modern genetics. Here, we report results of a combined multifactor analysis of properties characterizing functional sequence variants located in promoter regions of genes.

Results

We demonstrate that GC-content of the local sequence fragments and local DNA shape features play significant role in prioritization of functional variants and outscore features related to histone modifications, transcription factors binding sites, or evolutionary conservation descriptors. Those observations allowed us to build specialized machine learning classifier identifying functional single nucleotide polymorphisms within promoter regions—ShapeGTB. We compared our method with more general tools predicting pathogenicity of all non-coding variants. ShapeGTB outperformed them by a wide margin (average precision 0.93 vs. 0.47–0.55). On the external validation set based on ClinVar database it displayed worse performance but was still competitive with other methods (average precision 0.47 vs. 0.23–0.42). Such results suggest unique characteristics of mutations located within promoter regions and are a promising signal for the development of more accurate variant prioritization tools in the future.

Runs of homozygosity associate with decreased risks of lung cancer in never-smoking East Asian females

Although genome-wide association studies (GWASs) have identified some risk single-nucleotide polymorphisms in East Asian never-smoking females, the unexplained missing heritability is still required to be investigated. Runs of homozygosity (ROHs) are thought to be a type of genetic variation acting on human complex traits and diseases. We detected ROHs in 8,881 East Asian never-smoking women. The summed ROHs were used to fit a logistic regression model which noteworthily revealed a significant association between ROHs and the decreased risk of lung cancer (P < 0.05). We identified 4 common ROHs regions located at 2p22.1, which were significantly associated with decreased risk of lung cancer (P = 2.00 × 10^-4 - 1.35 × 10^-4). Functional annotation was conducted to investigate the regulatory function of ROHs. The common ROHs were overlapped with potential regulatory elements, such as active epigenome elements and chromatin states in lung-derived cell lines. SOS1 and ARHGEF33 were significantly up-regulated as the putative target genes of the identified ROHs in lung cancer samples according to the analysis of differently expressed genes. Our results suggest that ROHs could act as recessive contributing factors and regulatory elements to influence the risk of lung cancer in never-smoking East Asian females.

ENHO, RXRA, and LXRA polymorphisms and dyslipidaemia, related comorbidities and survival in haemodialysis patients

BACKGROUND

The energy homeostasis-associated gene (ENHO), retinoid X receptor alpha gene (RXRA), and liver X receptor alpha gene (LXRA) are involved in adipogenic/lipogenic regulation. We investigated whether single-nucleotide polymorphisms in these genes (ENHO rs2281997, rs72735260; RXRA rs749759, rs10776909, rs10881578; LXRA rs2279238, rs7120118, rs11039155) are associated with dyslipidaemia, related comorbidities and survival of haemodialysis (HD) patients also tested for T-helper (Th) cell interleukin genes (IL).

METHODS

The study was carried out in 873 HD patients. Dyslipidaemia was diagnosed by the recommendations of the Kidney Disease Outcomes Quality Initiative (K/DOQI) guidelines (2003); atherogenic dyslipidaemia was referred to if the TG/HDL cholesterol ratio was equal to or higher than 3.8. Genotyping of ENHO SNPs, LXRA SNPs, and IL12A rs568408 was carried out using HRM analysis. RXRA SNPs, IL12B rs3212227, and IL18 rs360719 were genotyped using PCR-RFLP analysis. The circulating adropin concentration was determined in 126 patients by enzyme-linked immunosorbent assay. Survival probability was analysed using the Kaplan-Meier method in 440 patients followed through 7.5 years.

RESULTS

Dyslipidaemia by K/DOQI was diagnosed in 459 patients (91% revealed hyper-LDL- cholesterolaemia), atherogenic dyslipidaemia was diagnosed in 454 patients, and 231 patients were free of dyslipidaemia by both criteria. The variant allele (T) of ENHO rs2281997 was associated with the hyper-LDL cholesterolaemic pattern of dyslipidaemia by K/DOQI. The frequency of atherogenic dyslipidaemia was lower in T-allele bearers than in CC-genotype patients. The rs2281997 T allele was associated with lower cardiovascular mortality in HD patients showing atherogenic dyslipidaemia. ENHO, RXRA, and LXRA showed epistatic interactions in dyslipidaemia. Circulating adropin was lower in atherogenic dyslipidaemia than in non-atherogenic conditions. RXRA rs10776909 was associated with myocardial infarction. Bearers of LXRA rs2279238, rs7120118 or rs11039155 minor alleles showed higher mortality. ENHO SNP positions fell within the same DNase 1 hypersensitivity site expressed in the Th1 cell line. Epistatic interactions occurred between rs2281997 and Th1 IL SNPs (rs360719, rs568408).

CONCLUSIONS

Atherogenic dyslipidaemia occurs in HD patients in whom ENHO encodes less adropin. ENHO, RXRA, and LXRA SNPs, separately or jointly, are associated with dyslipidaemia, myocardial infarction, and survival in HD patients. Differences in the availability of transcription binding sites may contribute to these associations.

Meta-analysis of epigenome-wide association studies of cognitive abilities

Cognitive functions are important correlates of health outcomes across the life-course. Individual differences in cognitive functions are partly heritable. Epigenetic modifications, such as DNA methylation, are susceptible to both genetic and environmental factors and may provide insights into individual differences in cognitive functions. Epigenome-wide meta-analyses for blood-based DNA methylation levels at ~420,000 CpG sites were performed for seven measures of cognitive functioning using data from 11 cohorts. CpGs that passed a Bonferroni correction, adjusting for the number of CpGs and cognitive tests, were assessed for: longitudinal change; being under genetic control (methylation QTLs); and associations with brain health (structural MRI), brain methylation and Alzheimer's disease pathology. Across the seven measures of cognitive functioning (meta-analysis n range: 2557-6809), there were epigenome-wide significant (P < 1.7 × 10-8) associations for global cognitive function (cg21450381, P = 1.6 × 10-8), and phonemic verbal fluency (cg12507869, P = 2.5 × 10-9). The CpGs are located in an intergenic region on chromosome 12 and the INPP5A gene on chromosome 10, respectively. Both probes have moderate correlations (~0.4) with brain methylation in Brodmann area 20 (ventral temporal cortex). Neither probe showed evidence of longitudinal change in late-life or associations with white matter brain MRI measures in one cohort with these data. A methylation QTL analysis suggested that rs113565688 was a cis methylation QTL for cg12507869 (P = 5 × 10-5 and 4 × 10-13 in two lookup cohorts). We demonstrate a link between blood-based DNA methylation and measures of phonemic verbal fluency and global cognitive ability. Further research is warranted to understand the mechanisms linking genomic regulatory changes with cognitive function to health and disease.

MC1R CpG island regulates MC1R expression and is methylated in a subset of melanoma tumours

Melanocortin 1 receptor (MC1R) is a G protein-coupled receptor expressed in melanocytes where it plays an important role in skin pigmentation and in the UV response, and has implications in melanoma development. Here we show that methylation of a CpG island (CGI) within the MC1R gene can control expression of MC1R in melanoma. This CGI overlaps with a potential enhancer region, and is unmethylated in normal melanocytes but highly methylated in other skin cells, suggesting a melanocyte specific function. Analysis showed that MC1R was the only gene significantly differentially expressed by methylation of this region. Within several data sets, this region is methylated in a subset of melanoma tumours (55%-74% of tumours) and results in reduced MC1R expression and significantly longer overall survival.

Topological alternate centrality measure capturing drug targets in the network of MAPK pathways

A new centrality of the nodes in the network is proposed called alternate centrality, which can isolate effective drug targets in the complex signalling network. Alternate centrality metric defined over the network substructure (four nodes - motifs). The nodes involving in alternative activation in the motifs gain in metric values. Targeting high alternative centrality nodes hypothesised to be destructive free to the network due to their alternative activation mechanism. Overlapping and crosstalk among the gene products in the conserved network of MAPK pathways selected for the study. In silico knock-out of high alternate centrality nodes causing rewiring in the network is investigated using MCF-7 breast cancer cell line-based data. Degree of top alternate centrality nodes lies between the degree of bridging and pagerank nodes. Node deletion of high alternate centrality on the centralities such as eccentricity, closeness, betweenness, stress, centroid and radiality causes low perturbation. The authors identified the following alternate centrality nodes ERK1, ERK2, MEKK2, MKK5, MKK4, MLK3, MLK2, MLK1, MEKK4, MEKK1, TAK1, P38alpha, ZAK, DLK, LZK, MLTKa/b and P38beta as efficient drug targets for breast cancer. Alternate centrality identifies effective drug targets and is free from intertwined biological processes and lethality.

Replication of a rare risk haplotype on 1p36.33 for autism spectrum disorder

Hundreds of genes have been implicated in autism spectrum disorders (ASDs). In genetically heterogeneous conditions, large families with multiple affected individuals provide strong evidence implicating a rare variant, and replication of the same variant in multiple families is unusual. We previously published linkage analyses and follow-up exome sequencing in seven large families with ASDs, implicating 14 rare exome variants. These included rs200195897, which was transmitted to four affected individuals in one family. We attempted replication of those variants in the MSSNG database. MSSNG is a unique resource for replication of ASD risk loci, containing whole genome sequence (WGS) on thousands of individuals diagnosed with ASDs and family members. For each exome variant, we obtained all carriers and their relatives in MSSNG, using a TDT test to quantify evidence for transmission and association. We replicated the transmission of rs200195897 to four affected individuals in three additional families. rs200195897 was also present in three singleton affected individuals, and no unaffected individuals other than transmitting parents. We identified two additional rare variants (rs566472488 and rs185038034) transmitted with rs200195897 on 1p36.33. Sanger sequencing confirmed the presence of these variants in the original family segregating rs200195897. To our knowledge, this is the first example of a rare haplotype being transmitted with ASD in multiple families. The candidate risk variants include a missense mutation in SAMD11, an intronic variant in NOC2L, and a regulatory region variant close to both genes. NOC2L is a transcription repressor, and several genes involved in transcription regulation have been previously associated with ASDs.

Demethylator phenotypes in acute myeloid leukemia

Acute myeloid leukemia (AML) often harbors mutations in epigenetic regulators, and also has frequent DNA hypermethylation, including the presence of CpG island methylator phenotypes (CIMPs). Although global hypomethylation is well known in cancer, the question of whether distinct demethylator phenotypes (DMPs) exist remains unanswered. Using Illumina 450k arrays for 194 patients from The Cancer Genome Atlas, we identified two distinct DMPs by hierarchical clustering: DMP.1 and DMP.2. DMP.1 cases harbored mutations in NPM1 (94%), FLT3 (71%) and DNMT3A (61%). Surprisingly, only 40% of patients with DNMT3A mutations were DMP.1, which has implications for mechanisms of transformation by this mutation. In contrast, DMP.2 AML was comprised of patients with t(8;21), inv(16) or t(15;17), suggesting common methylation defects connect these disparate rearrangements. RNA-seq revealed upregulated genes functioning in immune response (DMP.1) and development (DMP.2). We confirmed these findings by integrating independent 450k data sets (236 additional cases), and found prognostic effects by DMP status, independent of age and cytogenetics. The existence of DMPs has implications for AML pathogenesis and may augment existing tools in risk stratification.

Genome-wide association analysis of common genetic variants of resistant hypertension

Resistant hypertension (RHTN), defined as uncontrolled blood pressure (BP) ≥ 140/90 using three or more drugs or controlled BP (<140/90) using four or more drugs, is associated with adverse outcomes, including decline in kidney function. We conducted a genome-wide association analysis in 1194 White and Hispanic participants with hypertension and coronary artery disease from the INternational VErapamil-SR Trandolapril STudy-GENEtic Substudy (INVEST-GENES). Top variants associated with RHTN at p < 10^-4 were tested for replication in 585 White and Hispanic participants with hypertension and subcortical strokes from the Secondary Prevention of Subcortical Strokes GENEtic Substudy (SPS3-GENES). A genetic risk score for RHTN was created by summing the risk alleles of replicated RHTN signals. rs11749255 in MSX2 was associated with RHTN in INVEST (odds ratio (OR) (95% CI) = 1.50 (1.2-1.8), p = 7.3 × 10^-5) and replicated in SPS3 (OR = 2.0 (1.4-2.8), p = 4.3 × 10^-5), with genome-wide significance in meta-analysis (OR = 1.60 (1.3-1.9), p = 3.8 × 10^-8). Other replicated signals were in IFLTD1 and PTPRD. IFLTD1 rs6487504 was associated with RHTN in INVEST (OR = 1.90 (1.4-2.5), p = 1.1 × 10^-5) and SPS3 (OR = 1.70 (1.2-2.5), p = 4 × 10^-3). PTPRD rs324498, a previously reported RHTN signal, was among the top signals in INVEST (OR = 1.60 (1.3-2.0), p = 3.4 × 10^-5) and replicated in SPS3 (OR = 1.60 (1.1-2.4), one-sided p = 0.005). Participants with the highest number of risk alleles were at increased risk of RHTN compared to participants with a lower number (p-trend = 1.8 × 10^-15). Overall, we identified and replicated associations with RHTN in the MSX2, IFLTD1, and PTPRD regions, and combined these associations to create a genetic risk score.

Reduced monocyte and macrophage TNFSF15/TL1A expression is associated with susceptibility to inflammatory bowel disease

Chronic inflammation in inflammatory bowel disease (IBD) results from a breakdown of intestinal immune homeostasis and compromise of the intestinal barrier. Genome-wide association studies have identified over 200 genetic loci associated with risk for IBD, but the functional mechanisms of most of these genetic variants remain unknown. Polymorphisms at the TNFSF15 locus, which encodes the TNF superfamily cytokine commonly known as TL1A, are associated with susceptibility to IBD in multiple ethnic groups. In a wide variety of murine models of inflammation including models of IBD, TNFSF15 promotes immunopathology by signaling through its receptor DR3. Such evidence has led to the hypothesis that expression of this lymphocyte costimulatory cytokine increases risk for IBD. In contrast, here we show that the IBD-risk haplotype at TNFSF15 is associated with decreased expression of the gene by peripheral blood monocytes in both healthy volunteers and IBD patients. This association persists under various stimulation conditions at both the RNA and protein levels and is maintained after macrophage differentiation. Utilizing a "recall-by-genotype" bioresource for allele-specific expression measurements in a functional fine-mapping assay, we localize the polymorphism controlling TNFSF15 expression to the regulatory region upstream of the gene. Through a T cell costimulation assay, we demonstrate that genetically regulated TNFSF15 has functional relevance. These findings indicate that genetically enhanced expression of TNFSF15 in specific cell types may confer protection against the development of IBD.

Histone H2A.Z Suppression of Interferon-Stimulated Transcription and Antiviral Immunity Is Modulated by GCN5 and BRD2

Type I interferon (IFN)-stimulated gene (ISG) expression requires interaction between a transcription factor complex, ISGF3, and target gene promoters to initiate transcription and protection against infection. To uncover chromatin regulatory features of this antiviral immune response, IFN-induced nucleosome and histone dynamics of human ISG loci were examined. ISGF3 recruitment after IFN stimulation was accompanied by nucleosome reorganization at promoters and gene bodies. IFN stimulation induced loss of core histones H2B, H3, and H4, as well as H2A.Z at ISG promoters. A strong correlation was found between H2A.Z occupancy and ISGF3 target sites, and IFN-stimulated H2A.Z removal requires STAT1, STAT2, and IRF9. Neither INO80 nor SWI/SNF participate in IFN-driven H2A.Z eviction, but GCN5 and BRD2 are required. Interference with H2A.Z expression enhanced ISGF3 recruitment to ISG promoters, ISG mRNA expression, and IFN-stimulated antiviral immunity. This indicates that H2A.Z nucleosomes at ISG promoters restrict optimal ISGF3 engagement and modulate the biological response to IFN.

MiR-9, miR-153 and miR-124 are down-regulated by acute exposure to cocaine in a dopaminergic cell model and may contribute to cocaine dependence

Cocaine is one of the most used psychostimulant drugs worldwide. MicroRNAs are post-transcriptional regulators of gene expression that are highly expressed in brain, and several studies have shown that cocaine can alter their expression. In a previous study, we identified several protein-coding genes that are differentially expressed in a dopaminergic neuron-like model after an acute exposure to cocaine. Now, we used the prediction tool WebGestalt to identify miRNA molecules potentially involved in the regulation of these genes. Using the same cellular model, we found that seven of these miRNAs are down-regulated by cocaine: miR-124-3p, miR-124-5p, miR-137, miR-101-3p, miR-9-5p, miR-369-3p and miR-153-3p, the last three not previously related to cocaine. Furthermore, we found that three of the miRNA genes that are differentially expressed in our model (hsa-miR-9-1, hsa-miR-153-1 and hsa-miR-124-3) are nominally associated with cocaine dependence in a case-control study (2,085 cases and 4,293 controls). In summary, we highlighted novel miRNAs that may be involved in those cocaine-induced changes of gene expression that underlie addiction. Moreover, we identified genetic variants that contribute to cocaine dependence in three of these miRNA genes, supporting the idea that genes differentially expressed under cocaine may play an important role in the susceptibility to cocaine dependence.

Disease-Associated Short Tandem Repeats Co-localize with Chromatin Domain Boundaries

More than 25 inherited human disorders are caused by the unstable expansion of repetitive DNA sequences termed short tandem repeats (STRs). A fundamental unresolved question is why some STRs are susceptible to pathologic expansion, whereas thousands of repeat tracts across the human genome are relatively stable. Here, we discover that nearly all disease-associated STRs (daSTRs) are located at boundaries demarcating 3D chromatin domains. We identify a subset of boundaries with markedly higher CpG island density compared to the rest of the genome. daSTRs specifically localize to ultra-high-density CpG island boundaries, suggesting they might be hotspots for epigenetic misregulation or topological disruption linked to STR expansion. Fragile X syndrome patients exhibit severe boundary disruption in a manner that correlates with local loss of CTCF occupancy and the degree of FMR1 silencing. Our data uncover higher-order chromatin architecture as a new dimension in understanding repeat expansion disorders.

Continuous-Trait Probabilistic Model for Comparing Multi-species Functional Genomic Data

A large amount of multi-species functional genomic data from high-throughput assays are becoming available to help understand the molecular mechanisms for phenotypic diversity across species. However, continuous-trait probabilistic models, which are key to such comparative analysis, remain under-explored. Here we develop a new model, called phylogenetic hidden Markov Gaussian processes (Phylo-HMGP), to simultaneously infer heterogeneous evolutionary states of functional genomic features in a genome-wide manner. Both simulation studies and real data application demonstrate the effectiveness of Phylo-HMGP. Importantly, we applied Phylo-HMGP to analyze a new cross-species DNA replication timing (RT) dataset from the same cell type in five primate species (human, chimpanzee, orangutan, gibbon, and green monkey). We demonstrate that our Phylo-HMGP model enables discovery of genomic regions with distinct evolutionary patterns of RT. Our method provides a generic framework for comparative analysis of multi-species continuous functional genomic signals to help reveal regions with conserved or lineage-specific regulatory roles.

Curating Clinically Relevant Transcripts for the Interpretation of Sequence Variants

Variant interpretation depends on accurate annotations using biologically relevant transcripts. We have developed a systematic strategy for designating primary transcripts and have applied it to 109 hearing loss-associated genes that were divided into three categories. Category 1 genes (n = 38) had a single transcript; category 2 genes (n = 33) had multiple transcripts, but a single transcript was sufficient to represent all exons; and category 3 genes (n = 38) had multiple transcripts with unique exons. Transcripts were curated with respect to gene expression reported in the literature and the Genotype-Tissue Expression Project. In addition, high-frequency loss-of-function variants in the Genome Aggregation Database and disease-causing variants in ClinVar and the Human Gene Mutation Database across the 109 genes were queried. These data were used to classify exons as clinically significant, insignificant, or of uncertain significance. Interestingly, 6% of all exons, containing 124 reportedly disease-causing variants, were of uncertain significance. Finally, we used exon-level next-generation sequencing quality metrics generated at two clinical laboratories and identified a total of 43 technically challenging exons in 20 different genes that had inadequate coverage and/or homology issues that might lead to false-variant calls. We have demonstrated that transcript analysis plays a critical role in accurate clinical variant interpretation.