Mapping of transcription factor binding regions in mammalian cells by ChIP: comparison of array- and sequencing-based technologies

Genomic transcription factor binding site selection is edited by the chromatin remodeling factor CHD4

Biologically precise enhancer licensing by lineage-determining transcription factors enables activation of transcripts appropriate to biological demand and prevents deleterious gene activation. This essential process is challenged by the millions of matches to most transcription factor binding motifs present in many eukaryotic genomes, leading to questions about how transcription factors achieve the exquisite specificity required. The importance of chromatin remodeling factors to enhancer activation is highlighted by their frequent mutation in developmental disorders and in cancer. Here, we determine the roles of CHD4 in enhancer licensing and maintenance in breast cancer cells and during cellular reprogramming. In unchallenged basal breast cancer cells, CHD4 modulates chromatin accessibility. Its depletion leads to redistribution of transcription factors to previously unoccupied sites. During cellular reprogramming induced by the pioneer factor GATA3, CHD4 activity is necessary to prevent inappropriate chromatin opening. Mechanistically, CHD4 promotes nucleosome positioning over GATA3 binding motifs to compete with transcription factor-DNA interaction. We propose that CHD4 acts as a chromatin proof-reading enzyme that prevents unnecessary gene expression by editing chromatin binding activities of transcription factors.

Research Topics of the Bioinformatics of Gene Regulation

The study of gene expression regulation raises the challenge of developing bioinformatics tools and algorithms, demanding data integration [...].

Systematic identification and characterization of lncRNAs and lncRNA-miRNA-mRNA networks in the liver of turbot (Scophthalmus maximus L.) induced with Vibrio anguillarum

Long noncoding RNAs (lncRNAs), can regulate mRNA by targeting miRNA in a competing endogenous RNA network, have become a hot topic in the research of fish immune mechanism recent years. While in turbot (Scophthalmus maximus L.), an economically important marine fish, there are limited researches about the role of lncRNAs in its immune response to bacterial infection. In this study, a total of 184 differentially expressed lncRNAs (DElncRNAs) were systematically identified and characterized using whole-transcriptome sequencing of the liver of turbot challenged with Vibrioanguillarum at 0 h (control) and three different time points post infection (2 h, 12 h and 24 h, respectively). Subsequently, GO and KEGG signaling pathways of differentially expressed lncRNAs were analyzed to predict their function. We found that lncRNAs in our results were significantly enriched in several immune-related signaling pathways, including the NOD-like receptor signaling pathway, Toll-like receptor signaling pathway, Cytokine-cytokine receptor, MAPK signaling pathway, phagosome, PPAR signaling pathway and the regulation of autophagy. In addition, a total of 492 DE lncRNA - DE miRNA -DE mRNA networks were identified at three different time points post infection, which were consisted of 102 networks at 2 h, 122 networks at 12 h and 81 networks at 24 h post infection, respectively. Noticeably, 92 of these regulated networks were immune-related. These observations suggested that lncRNAs can regulate the expression of immune-related genes in the response to bacterial infection in turbot. Moreover, our findings would provide a new insight into the immune response of turbot to pathogen infection and lay a foundation for future study.

Posterior cingulate cortex reveals an expression profile of resilience in cognitively intact elders

The posterior cingulate cortex, a key hub of the default mode network, underlies autobiographical memory retrieval and displays hypometabolic changes early in Alzheimer disease. To obtain an unbiased understanding of the molecular pathobiology of the aged posterior cingulate cortex, we performed RNA sequencing (RNA-seq) on tissue obtained from 26 participants of the Rush Religious Orders Study (11 males/15 females; aged 76-96 years) with a pre-mortem clinical diagnosis of no cognitive impairment and post-mortem neurofibrillary tangle Braak Stages I/II, III, and IV. Transcriptomic data were gathered using next-generation sequencing of RNA extracted from posterior cingulate cortex generating an average of 60 million paired reads per subject. Normalized expression of RNA-seq data was calculated using a global gene annotation and a microRNA profile. Differential expression (DESeq2, edgeR) using Braak staging as the comparison structure isolated genes for dimensional scaling, associative network building and functional clustering. Curated genes were correlated with the Mini-Mental State Examination and semantic, working and episodic memory, visuospatial ability, and a composite Global Cognitive Score. Regulatory mechanisms were determined by co-expression networks with microRNAs and an overlap of transcription factor binding sites. Analysis revealed 750 genes and 12 microRNAs significantly differentially expressed between Braak Stages I/II and III/IV and an associated six groups of transcription factor binding sites. Inputting significantly different gene/network data into a functional annotation clustering model revealed elevated presynaptic, postsynaptic and ATP-related expression in Braak Stages III and IV compared with Stages I/II, suggesting these pathways are integral for cognitive resilience seen in unimpaired elderly subjects. Principal component analysis and Kruskal-Wallis testing did not associate Braak stage with cognitive function. However, Spearman correlations between genes and cognitive test scores followed by network analysis revealed upregulation of classes of synaptic genes positively associated with performance on the visuospatial perceptual orientation domain. Upregulation of key synaptic genes suggests a role for these transcripts and associated synaptic pathways in cognitive resilience seen in elders despite Alzheimer disease pathology and dementia.

Identification of phenotype-specific networks from paired gene expression-cell shape imaging data

The morphology of breast cancer cells is often used as an indicator of tumor severity and prognosis. Additionally, morphology can be used to identify more fine-grained, molecular developments within a cancer cell, such as transcriptomic changes and signaling pathway activity. Delineating the interface between morphology and signaling is important to understand the mechanical cues that a cell processes in order to undergo epithelial-to-mesenchymal transition and consequently metastasize. However, the exact regulatory systems that define these changes remain poorly characterized. In this study, we used a network-systems approach to integrate imaging data and RNA-seq expression data. Our workflow allowed the discovery of unbiased and context-specific gene expression signatures and cell signaling subnetworks relevant to the regulation of cell shape, rather than focusing on the identification of previously known, but not always representative, pathways. By constructing a cell-shape signaling network from shape-correlated gene expression modules and their upstream regulators, we found central roles for developmental pathways such as WNT and Notch, as well as evidence for the fine control of NF-kB signaling by numerous kinase and transcriptional regulators. Further analysis of our network implicates a gene expression module enriched in the RAP1 signaling pathway as a mediator between the sensing of mechanical stimuli and regulation of NF-kB activity, with specific relevance to cell shape in breast cancer.

Chromatin Profiling of Human Naïve Pluripotent Stem Cells

Chromatin immunoprecipitation combined with high-throughput sequencing (ChIP-sequencing) facilitates the genome-wide mapping of DNA sequences that are enriched for specific chromatin-binding proteins or histone post-translational modifications. More recently developed chromatin profiling methods called Cleavage Under Targets and Release Using Nuclease (CUT&RUN) and Cleavage Under Targets and Tagmentation (CUT&Tag) have adapted the ChIP-sequencing approach to produce similar data from a smaller amount of starting material, and while overcoming many of the conventional drawbacks of ChIP-sequencing. Here, we present detailed protocols for ChIP-seq, CUT&RUN, and CUT&Tag to profile genome-wide protein-DNA interactions in naïve human pluripotent stem cells.

Disease-related blood-based differential methylation in cystic fibrosis and its representation in lung cancer revealed a regulatory locus in PKP3 in lung epithelial cells

Cystic fibrosis (CF) is a monogenic disease, characterized by massive chronic lung inflammation. The observed variability in clinical phenotypes in monozygotic CF twins is likely associated with the extent of inflammation. This study sought to investigate inflammation-related aberrant DNA methylation in CF twins and to determine to what extent acquired methylation changes may be associated with lung cancer.

Blood-based genome-wide DNA methylation analysis was performed to compare the DNA methylomes of monozygotic twins, from the European CF Twin and Sibling Study with various degrees of disease severity. Putatively inflammation-related and differentially methylated positions were selected from a large lung cancer case-control study and investigated in blood by targeted bisulphite next-generation-sequencing. An inflammation-related locus located in the Plakophilin-3 (PKP3) gene was functionally analysed regarding promoter and enhancer activity in presence and absence of methylation using luciferase reporter assays.

We confirmed in a unique cohort that monozygotic twins, even if clinically discordant, have only minor differences in global DNA methylation patterns and blood cell composition. Further, we determined the most differentially methylated positions, a high proportion of which are blood cell-type-specific, whereas others may be acquired and thus have potential relevance in the context of inflammation as lung cancer risk factors. We identified a sequence in the gene body of PKP3 which is hypermethylated in blood from CF twins with severe phenotype and highly variably methylated in lung cancer patients and controls, independent of known clinical parameters, and showed that this region exhibits methylation-dependent promoter activity in lung epithelial cells.

Potential impact and mechanism of Long Non-coding RNAs on cancer and associated T cells

The discovery of many aberrant expressions of long non-coding RNAs (lncRNAs) in various cancers has focused attention on the effects of lncRNA on cancer cells themselves, including cell proliferation, growth inhibition, cell migration, cell immortality, vascular regeneration and cell viability. But with the increasing role of immunotherapy in cancer therapy, a large number of studies have revealed that the regulatory role of lncRNAs in immunity such as differentiation of immune cells can also influence the development and progression of cancer. In particular, recent publications have suggested that lncRNAs play critical roles in T-lymphocyte activation, proliferation, differentiation, function, apoptosis and metabolism. To elucidate the actual functions of lncRNAs at the molecular level of cancer pathogenesis, we summarize some of the current lncRNA regulatory mechanisms associated with T cell to discuss their effects in cancer in the hope of providing potential cancer therapeutic targets or cancer biomarkers. However, we all know that the differentiation and function of T cells is an extremely complex process that involves the expression and regulation of multiple lncRNAs. As a result, more regulatory mechanisms of lncRNAs need to be further studied.

CRED9: A differentially expressed elncRNA regulates expression of transcription factor CEBPA

Enhancer RNAs (eRNA) are non-coding transcripts produced from active enhancers and have potential gene regulatory function. CCAAT enhancer-binding protein alpha (CEBPA) is a transcription factor generally involved in metabolism, cell cycle inhibition, hematopoiesis, adipogenesis, hepatogenesis, and is associated with tumorigenesis. In this study, we demonstrate that an enhancer-associated long non-coding RNA (elncRNA), transcribed from an enhancer located 9kb downstream from the transcriptional start site (TSS) of CEBPA, positively regulates the expression of CEBPA. As a result, we named this elncRNA 'CEBPA regulatory elncRNA downstream 9kb' or 'CRED9'. CRED9 expression level positively correlates with CEBPA mRNA expression across multiple cell lines as detected by RT droplet digital PCR. Knockdown of CRED9 resulted in a reduction of CEBPA mRNA expression in Hep3B cells. Additionally, CRED9 knockdown in Hep3B and HepG2 cells resulted in lower CEBPA protein expression. We also found that knockdown of CRED9 in Hep3B cells caused a 57.8% reduction in H3K27ac levels at the +9kb CEBPA enhancer. H3K27ac has previously been described as a marker of active enhancers. Taken together, the evidence presented here supports a previously proposed model whereby, in some contexts, eRNA transcripts are necessary to amplify and maintain H3K27ac levels at a given enhancer. Ultimately, this study adds to the growing body of evidence that elncRNA transcripts have important roles in enhancer function and gene regulation.

Pan-cancer chromatin analysis of the human vtRNA genes uncovers their association with cancer biology

Background: The vault RNAs (vtRNAs) are a class of 84-141-nt eukaryotic non-coding RNAs transcribed by RNA polymerase III, associated to the ribonucleoprotein complex known as vault particle. Of the four human vtRNA genes, vtRNA1-1, vtRNA1-2 and vtRNA1-3, clustered at locus 1, are integral components of the vault particle, while vtRNA2-1 is a more divergent homologue located in a second locus. Gene expression studies of vtRNAs in large cohorts have been hindered by their unsuccessful sequencing using conventional transcriptomic approaches.

Methods: VtRNA expression in The Cancer Genome Atlas (TCGA) Pan-Cancer cohort was estimated using the genome-wide DNA methylation and chromatin accessibility data (ATAC-seq) of their genes as surrogate variables. The association between vtRNA expression and patient clinical outcome, immune subtypes and transcriptionally co-regulated gene programs was analyzed in the dataset.

Results: VtRNAs promoters are enriched in transcription factors related to viral infection. VtRNA2-1 is likely the most independently regulated homologue. VtRNA1-1 has the most accessible chromatin, followed by vtRNA1-2, vtRNA2-1 and vtRNA1-3. VtRNA1-1 and vtRNA1-3 chromatin status does not significantly change in cancer tissues. Meanwhile, vtRNA2-1 and vtRNA1-2 expression is widely deregulated in neoplastic tissues and its alteration is compatible with a broad oncogenic role for vtRNA1-2, and both tumor suppressor and oncogenic functions for vtRNA2-1. Yet, vtRNA1-1, vtRNA1-2 and vtRNA2-1 promoter DNA methylation predicts a shorter patient overall survival cancer-wide. In addition, gene ontology analyses of vtRNAs co-regulated genes identify a chromosome regulatory domain, epithelial differentiation, immune and thyroid cancer gene sets for specific vtRNAs. Furthermore, vtRNA expression patterns are associated with cancer immune subtypes and vtRNA1-2 expression is positively associated with cell proliferation and wound healing.

Conclusions: Our study presents the landscape of vtRNA chromatin status cancer-wide, identifying co-regulated gene networks and ontological pathways associated with the different vtRNA genes that may account for their diverse roles in cancer.

Genomic Space of MGMT in Human Glioma Revisited: Novel Motifs, Regulatory RNAs, NRF1, 2, and CTCF Involvement in Gene Expression

Background: The molecular regulation of increased MGMT expression in human brain tumors, the associated regulatory elements, and linkages of these to its epigenetic silencing are not understood. Because the heightened expression or non-expression of MGMT plays a pivotal role in glioma therapeutics, we applied bioinformatics and experimental tools to identify the regulatory elements in the MGMT and neighboring EBF3 gene loci. Results: Extensive genome database analyses showed that the MGMT genomic space was rich in and harbored many undescribed RNA regulatory sequences and recognition motifs. We extended the MGMT’s exon-1 promoter to 2019 bp to include five overlapping alternate promoters. Consensus sequences in the revised promoter for (a) the transcriptional factors CTCF, NRF1/NRF2, GAF, (b) the genetic switch MYC/MAX/MAD, and (c) two well-defined p53 response elements in MGMT intron-1, were identified. A putative protein-coding or non-coding RNA sequence was located in the extended 3′ UTR of the MGMT transcript. Eleven non-coding RNA loci coding for miRNAs, antisense RNA, and lncRNAs were identified in the MGMT-EBF3 region and six of these showed validated potential for curtailing the expression of both MGMT and EBF3 genes. ChIP analysis verified the binding site in MGMT promoter for CTCF which regulates the genomic methylation and chromatin looping. CTCF depletion by a pool of specific siRNA and shRNAs led to a significant attenuation of MGMT expression in human GBM cell lines. Computational analysis of the ChIP sequence data in ENCODE showed the presence of NRF1 in the MGMT promoter and this occurred only in MGMT-proficient cell lines. Further, an enforced NRF2 expression markedly augmented the MGMT mRNA and protein levels in glioma cells. Conclusions: We provide the first evidence for several new regulatory components in the MGMT gene locus which predict complex transcriptional and posttranscriptional controls with potential for new therapeutic avenues.

The lncRNA Toolkit: Databases and In Silico Tools for lncRNA Analysis

Long non-coding RNAs (lncRNAs) are a rapidly expanding field of research, with many new transcripts identified each year. However, only a small subset of lncRNAs has been characterized functionally thus far. To aid investigating the mechanisms of action by which new lncRNAs act, bioinformatic tools and databases are invaluable. Here, we review a selection of computational tools and databases for the in silico analysis of lncRNAs, including tissue-specific expression, protein coding potential, subcellular localization, structural conformation, and interaction partners. The assembled lncRNA toolkit is aimed primarily at experimental researchers as a useful starting point to guide wet-lab experiments, mainly containing multi-functional, user-friendly interfaces. With more and more new lncRNA analysis tools available, it will be essential to provide continuous updates and maintain the availability of key software in the future.

Comprehensive analysis of prognostic gene signatures based on immune infiltration of ovarian cancer

BACKGROUND

Ovarian cancer (OV) is one of the most common malignant tumors of gynecology oncology. The lack of effective early diagnosis methods and treatment strategies result in a low five-year survival rate. Also, immunotherapy plays an important auxiliary role in the treatment of advanced OV patient, so it is of great significance to find out effective immune-related tumor markers for the diagnosis and treatment of OV.

METHODS

Based on the consensus clustering analysis of single-sample gene set enrichment analysis (ssGSEA) score transformed via The Cancer Genome Atlas (TCGA) mRNA profile, we obtained two groups with high and low levels of immune infiltration. Multiple machine learning methods were conducted to explore prognostic genes associated with immune infiltration. Simultaneously, the correlation between the expression of mark genes and immune cells components was explored.

RESULTS

A prognostic classifier including 5 genes (CXCL11, S1PR4, TNFRSF17, FPR1 and DHRS95) was established and its robust efficacy for predicting overall survival was validated via 1129 OV samples. Some significant variations of copy number on gene loci were found between two risk groups and it showed that patients with fine chemosensitivity has lower risk score than patient with poor chemosensitivity (P = 0.013). The high and low-risk groups showed significantly different distribution (P < 0.001) of five immune cells (Monocytes, Macrophages M1, Macrophages M2, T cells CD4 menory and T cells CD8).

CONCLUSION

The present study identified five prognostic genes associated with immune infiltration of OV, which may provide some potential clinical implications for OV treatment.

Induced pluripotent stem cell reprogramming-associated methylation at the GABRA2 promoter and chr4p12 GABAA subunit gene expression in the context of alcohol use disorder

Twin studies indicate that there is a significant genetic contribution to the risk of developing alcohol use disorder (AUD). With the exception of coding variants in ADH1B and ALDH2, little is known about the molecular effects of AUD-associated loci. We previously reported that the AUD-associated synonymous polymorphism rs279858 within the GABAA α2 receptor subunit gene, GABRA2, was associated with gene expression of the chr4p12 GABAA subunit gene cluster in induced pluripotent stem cell (iPSC)-derived neural cultures. Based on this and other studies that showed changes in GABRA2 DNA methylation associated with schizophrenia and aging, we examined methylation in GABRA2. Specifically, using 69 iPSC lines and neural cultures derived from 47 of them, we examined whether GABRA2 rs279858 genotype predicted methylation levels and whether methylation was related to GABAA receptor subunit gene expression. We found that the GABRA2 CpG island undergoes random stochastic methylation during reprogramming and that methylation is associated with decreased GABRA2 gene expression, an effect that extends to the GABRB1 gene over 600 kb distal to GABRA2. Further, we identified additive effects of GABRA2 CpG methylation and GABRA2 rs279858 genotype on expression of the GABRB1 subunit gene in iPSC-derived neural cultures.

Enhancer occlusion transcripts regulate the activity of human enhancer domains via transcriptional interference: a computational perspective

Analysis of ENCODE long RNA-Seq and ChIP-seq (Chromatin Immunoprecipitation Sequencing) datasets for HepG2 and HeLa cell lines uncovered 1647 and 1958 transcripts that interfere with transcription factor binding to human enhancer domains. TFBSs (Transcription Factor Binding Sites) intersected by these 'Enhancer Occlusion Transcripts' (EOTrs) displayed significantly lower relative transcription factor (TF) binding affinities compared to TFBSs for the same TF devoid of EOTrs. Expression of most EOTrs was regulated in a cell line specific manner; analysis for the same TFBSs across cell lines, i.e. in the absence or presence of EOTrs, yielded consistently higher relative TF/DNA-binding affinities for TFBSs devoid of EOTrs. Lower activities of EOTr-associated enhancer domains coincided with reduced occupancy levels for histone tail modifications H3K27ac and H3K9ac. Similarly, the analysis of EOTrs with allele-specific expression identified lower activities for alleles associated with EOTrs. ChIA-PET (Chromatin Interaction Analysis by Paired-End Tag Sequencing) and 5C (Carbon Copy Chromosome Conformation Capture) uncovered that enhancer domains associated with EOTrs preferentially interacted with poised gene promoters. Analysis of EOTr regions with GRO-seq (Global run-on) data established the correlation of RNA polymerase pausing and occlusion of TF-binding. Our results implied that EOTr expression regulates human enhancer domains via transcriptional interference.

HSF1: Primary Factor in Molecular Chaperone Expression and a Major Contributor to Cancer Morbidity

Heat shock factor 1 (HSF1) is the primary component for initiation of the powerful heat shock response (HSR) in eukaryotes. The HSR is an evolutionarily conserved mechanism for responding to proteotoxic stress and involves the rapid expression of heat shock protein (HSP) molecular chaperones that promote cell viability by facilitating proteostasis. HSF1 activity is amplified in many tumor contexts in a manner that resembles a chronic state of stress, characterized by high levels of HSP gene expression as well as HSF1-mediated non-HSP gene regulation. HSF1 and its gene targets are essential for tumorigenesis across several experimental tumor models, and facilitate metastatic and resistant properties within cancer cells. Recent studies have suggested the significant potential of HSF1 as a therapeutic target and have motivated research efforts to understand the mechanisms of HSF1 regulation and develop methods for pharmacological intervention. We review what is currently known regarding the contribution of HSF1 activity to cancer pathology, its regulation and expression across human cancers, and strategies to target HSF1 for cancer therapy.

Biophysical Techniques for Target Validation and Drug Discovery in Transcription-Targeted Therapy

In the post-genome era, pathologies become associated with specific gene expression profiles and defined molecular lesions can be identified. The traditional therapeutic strategy is to block the identified aberrant biochemical activity. However, an attractive alternative could aim at antagonizing key transcriptional events underlying the pathogenesis, thereby blocking the consequences of a disorder, irrespective of the original biochemical nature. This approach, called transcription therapy, is now rendered possible by major advances in biophysical technologies. In the last two decades, techniques have evolved to become key components of drug discovery platforms, within pharmaceutical companies as well as academic laboratories. This review outlines the current biophysical strategies for transcription manipulation and provides examples of successful applications. It also provides insights into the future development of biophysical methods in drug discovery and personalized medicine.

Long non-coding RNAs in diseases related to inflammation and immunity

Long non-coding RNAs (lncRNAs) have been key regulators of gene expression in innate and adaptive immunity. Although lncRNAs have been reported to be associated with some diseases, its expression and function in diseases related to inflammation and immunity are still unknown. We reviewed how lncRNA regulated transcription and controlled the function and balance of the cells in the immune response. In addition, we discussed the impacts and challenges of lncRNAs on immunity in diseases.

Insights from resolving protein-DNA interactions at near base-pair resolution

One of the central goals in molecular biology is to understand how cell-type-specific expression patterns arise through selective recruitment of RNA polymerase II (Pol II) to a subset of gene promoters. Pol II needs to be recruited to a precise genomic position at the proper time to produce messenger RNA from a DNA template. Ostensibly, transcription is a relatively simple cellular process; yet, experimentally measuring and then understanding the combinatorial possibilities of transcriptional regulators remain a daunting task. Since its introduction in 1985, chromatin immunoprecipitation (ChIP) has remained a key tool for investigating protein-DNA contacts in vivo. Over 30 years of intensive research using ChIP have provided numerous insights into mechanisms of gene regulation. As functional genomic technologies improve, they present new opportunities to address key biological questions. ChIP-exo is a refined version of ChIP-seq that significantly reduces background signal, while providing near base-pair mapping resolution for protein-DNA interactions. This review discusses the evolution of the ChIP assay over the years; the methodological differences between ChIP-seq, ChIP-exo and ChIP-nexus; and highlight new insights into epigenetic and transcriptional mechanisms that were uniquely enabled with the near base-pair resolution of ChIP-exo.

PML/RARA inhibits expression of HSP90 and its target AKT

Essential for cell survival, the 90 kD Heat Shock Proteins (HSP90) are molecular chaperons required for conformational stabilization and trafficking of numerous client proteins. Functional HSP90 is required for the stability of AKT, a serine-threonine kinase phosphorylated in response to growth factor stimulation. AKT plays a crucial regulatory role in differentiation, cell cycle, transcription, translation, metabolism and apoptosis. Acute promyelocytic leukaemia (APL) is characterized by the presence of the promyelocytic leukaemia/retinoic acid receptor alpha (PML/RARA) fusion protein, which deregulates expression of several genes involved in differentiation and apoptosis. Here, we report inhibition of HSP90AA1 and HSP90AB1 isomer transcription in blasts isolated from patients with APL, associated with reduction of HSP90 protein expression and loss of control on AKT protein phosphorylation. We show that in vitro treatment of PML/RARA expressing cells with all-trans retinoic acid (ATRA) up-regulates HSP90 expression and stabilizes AKT. Addition of the HSP90-inhibitor 17-(allylamino)-17-demethoxygeldanamycin in combination with ATRA, blocks upregulation of AKT protein, indicating that HSP90 is necessary for ATRA action on AKT. This is the first report proving that expression of HSP90 isomers are directly and differentially repressed by PML/RARA, with critical results on cellular homeostasis of target proteins, such as AKT, in APL blasts.

LnChrom: a resource of experimentally validated lncRNA-chromatin interactions in human and mouse

Long non-coding RNAs (lncRNAs) constitute an important layer of chromatin regulation that contributes to various biological processes and diseases. By interacting with chromatin, many lncRNAs can regulate that state of chromatin by recruiting chromatin-modifying complexes and thus control large-scale gene expression programs. However, the available information on interactions between lncRNAs and chromatin is hidden in a large amount of dispersed literature and has not been extensively collected. We established the LnChrom database, a manually curated resource of experimentally validated lncRNA–chromatin interactions. The current release of LnChrom includes 382 743 interactions in human and mouse. We also manually collected detailed metadata for each interaction pair, including those of chromatin modifying factors, epigenetic marks and disease associations. LnChrom provides a user-friendly interface to facilitate browsing, searching and retrieving of lncRNA–chromatin interaction data. Additionally, a large amount of multi-omics data was integrated into LnChrom to aid in characterizing the effects of lncRNA–chromatin interactions on epigenetic modifications and transcriptional expression. We believe that LnChrom is a timely and valuable resource that can greatly motivate mechanistic research into lncRNAs.

Database URL: http://biocc.hrbmu.edu.cn/LnChrom/

Exopolysaccharides from Lactobacillus plantarum NCU116 Regulate Intestinal Barrier Function via STAT3 Signaling Pathway

Lactic acid bacteria (LAB) and their exopolysaccharides (EPS) are recognized to promote intestinal barrier function by mechanisms that remain incompletely understood. Herein, we sought to identify the roles of exopolysaccharides from Lactobacillus plantarum NCU116 (EPS116) in intestinal barrier function. Our data showed that EPS116 attenuated dextran sodium sulfate (DSS) induced colitis and promoted epithelial barrier function and the expression of tight junction (TJ) proteins in vivo and in vitro. Moreover, chromatin immunoprecipitation data showed that EPS116 facilitated STAT3 (signal transducer and activator of transcription 3) binding to the promoter of occludin and ZO-1. Furthermore, knockdown of STAT3 in Caco-2 cell with EPS116 treatment led to decreased expression of occludin and ZO-1 and increased intestinal permeability, suggesting that the regulation of epithelial barrier function by EPS116 should be STAT3 dependent. Thus, our data revealed a novel mechanism that EPS116 inhibited intestinal inflammation via regulating intestinal epithelial barrier function.

Analysis of sensitive information leakage in functional genomics signal profiles through genomic deletions

Functional genomics experiments, such as RNA-seq, provide non-individual specific information about gene expression under different conditions such as disease and normal. There is great desire to share these data. However, privacy concerns often preclude sharing of the raw reads. To enable safe sharing, aggregated summaries such as read-depth signal profiles and levels of gene expression are used. Projects such as GTEx and ENCODE share these because they ostensibly do not leak much identifying information. Here, we attempt to quantify the validity of this statement, measuring the leakage of genomic deletions from signal profiles. We present information theoretic measures for the degree to which one can genotype these deletions. We then develop practical genotyping approaches and demonstrate how to use these to identify an individual within a large cohort in the context of linking attacks. Finally, we present an anonymization method removing much of the leakage from signal profiles.

Transcriptome analysis reveals a long non-coding RNA signature to improve biochemical recurrence prediction in prostate cancer

Despite highly successful treatments for localized prostate cancer (PCa), prognostic biomarkers are needed to improve patient management and prognosis. Accumulating evidence suggests that long noncoding RNAs (lncRNAs) are key regulators with biological and clinical significance. By transcriptome analysis, we identified a set of consistently dysregulated lncRNAs in PCa across different datasets and revealed an eight-lncRNA signature that significantly associated with the biochemical recurrence (BCR)-free survival. Based on the signature, patients could be classified into high- and low-risk groups with significantly different survival (HR = 2.19; 95% CI = 1.67-2.88; P < 0.0001). Validations in the validation cohorts and another independent cohort confirmed its prognostic value for recurrence prediction. Multivariable analysis showed that the signature was independent of common clinicopathological features and stratified analysis further revealed its role in elevating risk stratification of current prognostic models. Additionally, the eight-lncRNA signature was able to improve on the CAPRA-S score for the prediction of BCR as well as to reflect the metastatic potential of PCa. Functional characterization suggested that these lncRNAs which showed PCa-specific expression patterns may involve in critical processes in tumorigenesis. Overall, our results demonstrated potential application of lncRNAs as novel independent biomarkers. The eight-lncRNA signature may have clinical potential for facilitating further stratification of more aggressive patients who would benefit from adjuvant therapy.

Multi-focal control of mitochondrial gene expression by oncogenic MYC provides potential therapeutic targets in cancer

Despite ubiquitous activation in human cancer, essential downstream effector pathways of the MYC transcription factor have been difficult to define and target. Using a structure/function-based approach, we identified the mitochondrial RNA polymerase (POLRMT) locus as a critical downstream target of MYC. The multifunctional POLRMT enzyme controls mitochondrial gene expression, a process required both for mitochondrial function and mitochondrial biogenesis. We further demonstrate that inhibition of this newly defined MYC effector pathway causes robust and selective tumor cell apoptosis, via an acute, checkpoint-like mechanism linked to aberrant electron transport chain complex assembly and mitochondrial reactive oxygen species (ROS) production. Fortuitously, MYC-dependent tumor cell death can be induced by inhibiting the mitochondrial gene expression pathway using a variety of strategies, including treatment with FDA-approved antibiotics. In vivo studies using a mouse model of Burkitt's Lymphoma provide pre-clinical evidence that these antibiotics can successfully block progression of MYC-dependent tumors.

Transcriptional interference by small transcripts in proximal promoter regions

Proximal promoter regions (PPR) are heavily transcribed yielding different types of small RNAs. The act of transcription within PPRs might regulate downstream gene expression via transcriptional interference (TI). For analysis, we investigated capped and polyadenylated small RNA transcripts within PPRs of human RefSeq genes in eight different cell lines. Transcripts of our datasets overlapped with experimentally determined transcription factor binding sites (TFBS). For TFBSs intersected by these small RNA transcripts, we established negative correlation of sRNA expression levels and transcription factor (TF) DNA binding affinities; suggesting that the transcripts acted via TI. Accordingly, datasets were designated as TFbiTrs (TF-binding interfering transcripts). Expression of most TFbiTrs was restricted to certain cell lines. This facilitated the analysis of effects related to TFbiTr expression for the same RefSeq genes across cell lines. We consistently uncovered higher relative TF/DNA binding affinities and concomitantly higher expression levels for RefSeq genes in the absence of TFbiTrs. Analysis of corresponding chromatin landscapes supported these results. ChIA-PET revealed the participation of distal enhancers in TFbiTr transcription. Enhancers regulating TFbiTrs, in effect, act as repressors for corresponding downstream RefSeq genes. We demonstrate the significant impact of TI on gene expression using selected small RNA datasets.

Functional variant in the promoter region of IL-27 alters gene transcription and confers a risk for ulcerative colitis in northern Chinese Han

Ulcerative colitis (UC) is a chronic inflammatory disorder of unknown etiology and a polygenic disease. IL-27 encodes p28, a subunit of IL-12 family cytokines, and has been implicated in the pathogenesis of UC. The aims of the present study were to evaluate the genetic association of a variant of the IL-27 gene with UC and to further characterize the functional variant in the IL-27 gene that influences the risk for UC. Our data demonstrated that the genetic variant rs153109 in the 5' upstream region of IL-27 is significantly associated with UC in Chinese Han individuals. Analysis of IL-27 transcripts demonstrated that individuals carrying the risk allele of rs153109 display reduced transcription of IL-27 in PBMCs. Luciferase activity assays demonstrated that the risk allele rs153109 results in decreased promoter activity compared to a non-risk allele in a tissue specific manner. Mechanistic characterization of histone modifications in the promoter region revealed that the risk haplotype tagged by the risk allele of rs153109 reduces the levels of H3K3me3 and H3K27ac.

Mathematical Modeling of Avidity Distribution and Estimating General Binding Properties of Transcription Factors from Genome-Wide Binding Profiles

The shape of the experimental frequency distributions (EFD) of diverse molecular interaction events quantifying genome-wide binding is often skewed to the rare but abundant quantities. Such distributions are systematically deviated from standard power-law functions proposed by scale-free network models suggesting that more explanatory and predictive probabilistic model(s) are needed. Identification of the mechanism-based data-driven statistical distributions that provide an estimation and prediction of binding properties of transcription factors from genome-wide binding profiles is the goal of this analytical survey. Here, we review and develop an analytical framework for modeling, analysis, and prediction of transcription factor (TF) DNA binding properties detected at the genome scale. We introduce a mixture probabilistic model of binding avidity function that includes nonspecific and specific binding events. A method for decomposition of specific and nonspecific TF-DNA binding events is proposed. We show that the Kolmogorov-Waring (KW) probability function (PF), modeling the steady state TF binding-dissociation stochastic process, fits well with the EFD for diverse TF-DNA binding datasets. Furthermore, this distribution predicts total number of TF-DNA binding sites (BSs), estimating specificity and sensitivity as well as other basic statistical features of DNA-TF binding when the experimental datasets are noise-rich and essentially incomplete. The KW distribution fits equally well to TF-DNA binding activity for different TFs including ERE, CREB, STAT1, Nanog, and Oct4. Our analysis reveals that the KW distribution and its generalized form provides the family of power-law-like distributions given in terms of hypergeometric series functions, including standard and generalized Pareto and Waring distributions, providing flexible and common skewed forms of the transcription factor binding site (TFBS) avidity distribution function. We suggest that the skewed binding events may be due to a wide range of evolutionary processes of creating weak avidity TFBS associated with random mutations, while the rare high-avidity binding sites (i.e., high-avidity evolutionarily conserved canonical e-boxes) rarely occurred. These, however, may be positively selected in microevolution.

Resolving the Genomic Localization of the Kollerin Cohesin-Loader Complex

The kollerin complex, consisting of Scc2/Scc4 in yeast and Nipbl/Mau2 in vertebrates, is crucial for the chromatin-association of the cohesin complex and therefore for the critical functions of cohesin in cell division, transcriptional regulation and chromatin organisation. Despite the recent efforts to determine the genomic localization of the kollerin complex in different cell lines, major questions still remain unresolved, for instance where cohesin is actually loaded onto chromatin. Further, Nipbl seems to have also additional roles, for instance as transcription factor.This chapter summarizes our current knowledge on kollerin function and the recent studies on the genomic localization of Scc2, highlighting and critically discussing controversial data.

Omics Profiling in Precision Oncology

Cancer causes significant morbidity and mortality worldwide, and is the area most targeted in precision medicine. Recent development of high-throughput methods enables detailed omics analysis of the molecular mechanisms underpinning tumor biology. These studies have identified clinically actionable mutations, gene and protein expression patterns associated with prognosis, and provided further insights into the molecular mechanisms indicative of cancer biology and new therapeutics strategies such as immunotherapy. In this review, we summarize the techniques used for tumor omics analysis, recapitulate the key findings in cancer omics studies, and point to areas requiring further research on precision oncology.

Role and Regulation of Myeloid Zinc Finger Protein 1 in Cancer

Myeloid zinc finger 1 (MZF1) belongs to the SCAN-Zinc Finger (SCAN-ZF) transcription factor family that has recently been implicated in a number of types of cancer. Although the initial studies concentrated on the role of MZF1 in myeloid differentiation and leukemia, the factor now appears to be involved in the etiology of major solid tumors such as lung, cervical, breast, and colorectal cancer. Here we discuss the regulation of MZF1 that mediated its recruitment and activation in cancer, concentrating on posttranslational modification by phosphorylation, and sumoylation, formation of promyelocytic leukemia nuclear bodies and its association with co-activators and co-repressors.

B cell receptor induced Fc receptor-like 5 expression is mediated by multiple signaling pathways converging on NF-κB and NFAT

Fc receptor-like (FCRL) proteins are novel regulators of the B cell response to antigen. Human FCRL5 binds intact IgG and modifies the strength of antigen receptor (BCR) signaling. Altering FCRL5 expression could therefore regulate the B cell response to antigen. In this study, we found that FCRL5 expression is induced specifically upon BCR stimulation and dissected the molecular mechanism. FCRL5 mRNA and cell surface protein expression required prolonged BCR stimulation and de novo protein synthesis. Using chemical inhibitors and activators, we identified roles for several signaling pathways, indicating a complex mechanism. Specifically, the PI3K/AKT, JNK, PKC and IKK2-dependent classical NF-κB pathways were involved in induced FCRL5 expression. Furthermore, induced FCRL5 expression required elevation of intracellular Ca(++) and was partially blocked by cyclosporine A, a calcineurin inhibitor. The importance of the transcription factors NF-κB, NFAT and CREB-binding protein was revealed based on sensitivity to inhibitors. Using reporter gene assays, we showed that the core FCRL5 promoter was sufficient to drive induced gene expression. Mutations of two predicted NF-κB sites or an NFAT site in the core promoter abrogated induced gene expression, suggesting direct regulation of the FCRL5 gene by NF-κB and NFAT. In support, we detected binding of NF-κB and NFAT family proteins to oligonucleotides corresponding to the predicted sites. We propose that the identified intricate mechanism serves to ensure that FCRL5 is expressed on B cells at a precise time following antigen encounter, with potential implications regarding regulation of the B cell response.

Unlinking an lncRNA from Its Associated cis Element

Long non-coding (lnc) RNAs can regulate gene expression and protein functions. However, the proportion of lncRNAs with biological activities among the thousands expressed in mammalian cells is controversial. We studied Lockd (lncRNA downstream of Cdkn1b), a 434-nt polyadenylated lncRNA originating 4 kb 3' to the Cdkn1b gene. Deletion of the 25-kb Lockd locus reduced Cdkn1b transcription by approximately 70% in an erythroid cell line. In contrast, homozygous insertion of a polyadenylation cassette 80 bp downstream of the Lockd transcription start site reduced the entire lncRNA transcript level by >90% with no effect on Cdkn1b transcription. The Lockd promoter contains a DNase-hypersensitive site, binds numerous transcription factors, and physically associates with the Cdkn1b promoter in chromosomal conformation capture studies. Therefore, the Lockd gene positively regulates Cdkn1b transcription through an enhancer-like cis element, whereas the lncRNA itself is dispensable, which may be the case for other lncRNAs.

Genetic variants associated with autoimmunity drive NFκB signaling and responses to inflammatory stimuli

The transcription factor nuclear factor κB (NFκB) is a central regulator of inflammation, and genome-wide association studies in subjects with autoimmune disease have identified a number of variants within the NFκB signaling cascade. In addition, causal variant fine-mapping has demonstrated that autoimmune disease susceptibility variants for multiple sclerosis (MS) and ulcerative colitis are strongly enriched within binding sites for NFκB. We report that MS-associated variants proximal to NFκB1 and in an intron of TNFRSF1A (TNFR1) are associated with increased NFκB signaling after tumor necrosis factor-α (TNFα) stimulation. Both variants result in increased degradation of inhibitor of NFκB α (IκBα), a negative regulator of NFκB, and nuclear translocation of p65 NFκB. The variant proximal to NFκB1 controls signaling responses by altering the expression of NFκB itself, with the GG risk genotype expressing 20-fold more p50 NFκB and diminished expression of the negative regulators of the NFκB pathway: TNFα-induced protein 3 (TNFAIP3), B cell leukemia 3 (BCL3), and cellular inhibitor of apoptosis 1 (CIAP1). Finally, naïve CD4 T cells from patients with MS express enhanced activation of p65 NFκB. These results demonstrate that genetic variants associated with risk of developing MS alter NFκB signaling pathways, resulting in enhanced NFκB activation and greater responsiveness to inflammatory stimuli. As such, this suggests that rapid genetic screening for variants associated with NFκB signaling may identify individuals amenable to NFκB or cytokine blockade.

Long non-coding RNAs in innate and adaptive immunity

Long noncoding RNAs (lncRNAs) represent a newly discovered class of regulatory molecules that impact a variety of biological processes in cells and organ systems. In humans, it is estimated that there may be more than twice as many lncRNA genes than protein-coding genes. However, only a handful of lncRNAs have been analyzed in detail. In this review, we describe expression and functions of lncRNAs that have been demonstrated to impact innate and adaptive immunity. These emerging paradigms illustrate remarkably diverse mechanisms that lncRNAs utilize to impact the transcriptional programs of immune cells required to fight against pathogens and maintain normal health and homeostasis.

DNA methylation reactivates GAD1 expression in cancer by preventing CTCF-mediated polycomb repressive complex 2 recruitment

Levels of γ-aminobutyric acid (GABA) and glutamic acid decarboxylase 1 (GAD1), the enzyme that synthesizes GABA, are significantly increased in neoplastic tissues. However, the mechanism underlying this increase remains elusive. Instead of silencing gene transcription, we showed that the GAD1 promoter was hypermethylated in both colon and liver cancer cells, leading to the production of high levels of GAD1. GAD1 is a target gene that is silenced by H3K27me3. The key locus responsible for GAD1 reactivation was mapped to a DNA methylation-sensitive CTCF-binding site (CTCF-BS3) within the third intron of GAD1. Chromosome configuration capture (3C) analysis indicated that an intrachromosomal loop was formed by CTCF self-dimerisation in normal cells (CTCF binds to both unmethylated CTCF-BS3 and CTCF-BS2). The CTCF dimer then interacted with suppressor of zeste 12 homologue (SUZ12), which is a domain of Polycomb repressive complex 2 (PRC2), promoting the methylation of H3K27 and the silencing of GAD1 expression. This silencing was shown to be inhibited by DNA methylation in cancer cells. These findings strongly suggest that GAD1 is reactivated by DNA methylation, which provided a model for DNA methylation and the active orchestration of oncogenic gene expression by CTCF in cancer cells.

Long noncoding RNAs in T lymphocytes

Long noncoding RNAs are recently discovered regulatory RNA molecules that do not code for proteins but influence a vast array of biologic processes. In vertebrates, the number of long noncoding RNA genes is thought to greatly exceed the number of protein-coding genes. It is also thought that long noncoding RNAs drive the biologic complexity observed in vertebrates compared with that in invertebrates. Evidence of this complexity has been found in the T-lymphocyte compartment of the adaptive immune system. In the present review, we describe our current level of understanding of the expression of specific long or large intergenic or intervening long noncoding RNAs during T-lymphocyte development in the thymus and differentiation in the periphery and highlight the mechanisms of action that specific long noncoding RNAs employ to regulate T-lymphocyte function, both in vitro and in vivo.

Variations in potassium channel genes are associated with distinct trajectories of persistent breast pain after breast cancer surgery

Persistent pain after breast cancer surgery is a common clinical problem. Given the role of potassium channels in modulating neuronal excitability, coupled with recently published genetic associations with preoperative breast pain, we hypothesized that variations in potassium channel genes will be associated with persistent postsurgical breast pain. In this study, associations between 10 potassium channel genes and persistent breast pain were evaluated. Using growth mixture modeling (GMM), 4 distinct latent classes of patients, who were assessed before and monthly for 6 months after breast cancer surgery, were identified previously (ie, No Pain, Mild Pain, Moderate Pain, Severe Pain). Genotyping was done using a custom array. Using logistic regression analyses, significant differences in a number of genotype or haplotype frequencies were found between: Mild Pain vs No Pain and Severe Pain vs No Pain classes. Seven single-nucleotide polymorphisms (SNPs) across 5 genes (ie, potassium voltage-gated channel, subfamily A, member 1 [KCNA1], potassium voltage-gated channel, subfamily D, member 2 [KCND2], potassium inwardly rectifying channel, subfamily J, members 3 and 6 (KCNJ3 and KCNJ6), potassium channel, subfamily K, member 9 [KCNK9]) were associated with membership in the Mild Pain class. In addition, 3 SNPs and 1 haplotype across 4 genes (ie, KCND2, KCNJ3, KCNJ6, KCNK9) were associated with membership in the Severe Pain class. These findings suggest that variations in potassium channel genes are associated with both mild and severe persistent breast pain after breast cancer surgery. Although findings from this study warrant replication, they provide intriguing preliminary information on potential therapeutic targets.

Functional diversity of CTCFs is encoded in their binding motifs

BACKGROUND

The CCCTC-binding factor (CTCF) has diverse regulatory functions. However, the definitive characteristics of the CTCF binding motif required for its functional diversity still remains elusive.

RESULTS

Here, we describe a new motif discovery workflow by which we have identified three CTCF binding motif variations with highly divergent functionalities. Supported by transcriptomic, epigenomic and chromatin-interactomic data, we show that the functional diversity of the CTCF binding motifs is strongly associated with their GC content, CpG dinucleotide coverage and relative DNA methylation level at the 12th position of the motifs. Further analysis suggested that the co-localization of cohesin, the key factor in cohesion of sister chromatids, is negatively correlated with the CpG coverage and the relative DNA methylation level at the 12th position. Finally, we present evidences for a hypothetical model in which chromatin interactions between promoters and distal regulatory regions are likely mediated by CTCFs binding to sequences with high CpG.

CONCLUSION

These results demonstrate the existence of definitive CTCF binding motifs corresponding to CTCF's diverse functions, and that the functional diversity of the motifs is strongly associated with genetic and epigenetic features at the 12th position of the motifs.

DIS3 shapes the RNA polymerase II transcriptome in humans by degrading a variety of unwanted transcripts

Human DIS3, the nuclear catalytic subunit of the exosome complex, contains exonucleolytic and endonucleolytic active domains. To identify DIS3 targets genome-wide, we combined comprehensive transcriptomic analyses of engineered HEK293 cells that expressed mutant DIS3, with Photoactivatable Ribonucleoside-Enhanced Cross-Linking and Immunoprecipitation (PAR-CLIP) experiments. In cells expressing DIS3 with both catalytic sites mutated, RNAs originating from unannotated genomic regions increased ∼2.5-fold, covering ∼70% of the genome and allowing for thousands of novel transcripts to be discovered. Previously described pervasive transcription products, such as Promoter Upstream Transcripts (PROMPTs), accumulated robustly upon DIS3 dysfunction, representing a significant fraction of PAR-CLIP reads. We have also detected relatively long putative premature RNA polymerase II termination products of protein-coding genes whose levels in DIS3 mutant cells can exceed the mature mRNAs, indicating that production of such truncated RNA is a common phenomenon. In addition, we found DIS3 to be involved in controlling the formation of paraspeckles, nuclear bodies that are organized around NEAT1 lncRNA, whose short form was overexpressed in cells with mutated DIS3. Moreover, the DIS3 mutations resulted in misregulation of expression of ∼50% of transcribed protein-coding genes, probably as a secondary effect of accumulation of various noncoding RNA species. Finally, cells expressing mutant DIS3 accumulated snoRNA precursors, which correlated with a strong PAR-CLIP signal, indicating that DIS3 is the main snoRNA-processing enzyme. EXOSC10 (RRP6) instead controls the levels of the mature snoRNAs. Overall, we show that DIS3 has a major nucleoplasmic function in shaping the human RNA polymerase II transcriptome.

Regulation and function of the human HSP90AA1 gene

Heat shock protein 90α (Hsp90α), encoded by the HSP90AA1 gene, is the stress inducible isoform of the molecular chaperone Hsp90. Hsp90α is regulated differently and has different functions when compared to the constitutively expressed Hsp90β isoform, despite high amino acid sequence identity between the two proteins. These differences are likely due to variations in nucleotide sequence within non-coding regions, which allows for specific regulation through interaction with particular transcription factors, and to subtle changes in amino acid sequence that allow for unique post-translational modifications. This article will specifically focus on the expression, function and regulation of Hsp90α.

Mining for single nucleotide variants (SNVs) at the kallikrein locus with predicted functional consequences

Kallikreins (KLKs) are a group of 15 serine proteases encoded by the KLK locus on chromosome 19. Certain single nucleotide variants (SNVs) within the KLK locus have been linked to human disease. Next-generation sequencing of large human cohorts enables reexamination of genomic variation at the KLK locus. We aimed to identify all KLK-related SNVs and examine their impact on gene regulation and function. To this end, we mined KLK SNVs across Ensembl and Exome Variant Server, with exome-sequencing data from 6503 individuals. PolyPhen-2-based prediction of damaging SNVs and population frequencies of these SNVs were examined. Damaging SNVs were plotted on protein sequence and structure. We identified 4866 SNVs, the largest number of KLK-related SNVs reported. Fourteen percent of noncoding SNVs overlapped with transcription factor binding sites. We identified 602 missense coding SNVs, among which 148 were predicted to be damaging. Nine missense SNVs were common (>1% frequency) and displayed significantly different frequencies between European-American and African-American populations. SNVs predicted to be damaging appeared to alter tertiary structure of KLK1 and KLK6. Similarly, these missense SNVs may affect KLK function, resulting in disease phenotypes. Our study represents a mine of information for those studying KLK-related SNVs and their associations with diseases.

Global analysis of transcription factor-binding sites in yeast using ChIP-Seq

Transcription factors influence gene expression through their ability to bind DNA at specific regulatory elements. Specific DNA-protein interactions can be isolated through the chromatin immunoprecipitation (ChIP) procedure, in which DNA fragments bound by the protein of interest are recovered. ChIP is followed by high-throughput DNA sequencing (Seq) to determine the genomic provenance of ChIP DNA fragments and their relative abundance in the sample. This chapter describes a ChIP-Seq strategy adapted for budding yeast to enable the genome-wide characterization of binding sites of transcription factors (TFs) and other DNA-binding proteins in an efficient and cost-effective way.Yeast strains with epitope-tagged TFs are most commonly used for ChIP-Seq, along with their matching untagged control strains. The initial step of ChIP involves the cross-linking of DNA and proteins. Next, yeast cells are lysed and sonicated to shear chromatin into smaller fragments. An antibody against an epitope-tagged TF is used to pull down chromatin complexes containing DNA and the TF of interest. DNA is then purified and proteins degraded. Specific barcoded adapters for multiplex DNA sequencing are ligated to ChIP DNA. Short DNA sequence reads (28-36 base pairs) are parsed according to the barcode and aligned against the yeast reference genome, thus generating a nucleotide-resolution map of transcription factor-binding sites and their occupancy.

Emerging applications of read profiles towards the functional annotation of the genome

Functional annotation of the genome is important to understand the phenotypic complexity of various species. The road toward functional annotation involves several challenges ranging from experiments on individual molecules to large-scale analysis of high-throughput sequencing (HTS) data. HTS data is typically a result of the protocol designed to address specific research questions. The sequencing results in reads, which when mapped to a reference genome often leads to the formation of distinct patterns (read profiles). Interpretation of these read profiles is essential for their analysis in relation to the research question addressed. Several strategies have been employed at varying levels of abstraction ranging from a somewhat ad hoc to a more systematic analysis of read profiles. These include methods which can compare read profiles, e.g., from direct (non-sequence based) alignments to classification of patterns into functional groups. In this review, we highlight the emerging applications of read profiles for the annotation of non-coding RNA and cis-regulatory elements (CREs) such as enhancers and promoters. We also discuss the biological rationale behind their formation.

Quantitative analysis of the TNF-α-induced phosphoproteome reveals AEG-1/MTDH/LYRIC as an IKKβ substrate

The inhibitor of the nuclear factor-κB (IκB) kinase (IKK) complex is a key regulator of the canonical NF-κB signalling cascade and is crucial for fundamental cellular functions, including stress and immune responses. The majority of IKK complex functions are attributed to NF-κB activation; however, there is increasing evidence for NF-κB pathway-independent signalling. Here we combine quantitative mass spectrometry with random forest bioinformatics to dissect the TNF-α-IKKβ-induced phosphoproteome in MCF-7 breast cancer cells. In total, we identify over 20,000 phosphorylation sites, of which ∼1% are regulated up on TNF-α stimulation. We identify various potential novel IKKβ substrates including kinases and regulators of cellular trafficking. Moreover, we show that one of the candidates, AEG-1/MTDH/LYRIC, is directly phosphorylated by IKKβ on serine 298. We provide evidence that IKKβ-mediated AEG-1 phosphorylation is essential for IκBα degradation as well as NF-κB-dependent gene expression and cell proliferation, which correlate with cancer patient survival in vivo.

Inflammatory cytokines such as TNF-α influence inflammation, apoptosis and tumour development through regulation of the kinase IKKβ. Krishnan and Nolte et al. apply quantitative proteomics to identify potential IKKβ targets, and reveal phosphorylation of AEG-1 by IKKβ as a mechanism controlling NF-κB signalling.

MACE: model based analysis of ChIP-exo

Understanding the role of a given transcription factor (TF) in regulating gene expression requires precise mapping of its binding sites in the genome. Chromatin immunoprecipitation-exo, an emerging technique using λ exonuclease to digest TF unbound DNA after ChIP, is designed to reveal transcription factor binding site (TFBS) boundaries with near-single nucleotide resolution. Although ChIP-exo promises deeper insights into transcription regulation, no dedicated bioinformatics tool exists to leverage its advantages. Most ChIP-seq and ChIP-chip analytic methods are not tailored for ChIP-exo, and thus cannot take full advantage of high-resolution ChIP-exo data. Here we describe a novel analysis framework, termed MACE (model-based analysis of ChIP-exo) dedicated to ChIP-exo data analysis. The MACE workflow consists of four steps: (i) sequencing data normalization and bias correction; (ii) signal consolidation and noise reduction; (iii) single-nucleotide resolution border peak detection using the Chebyshev Inequality and (iv) border matching using the Gale-Shapley stable matching algorithm. When applied to published human CTCF, yeast Reb1 and our own mouse ONECUT1/HNF6 ChIP-exo data, MACE is able to define TFBSs with high sensitivity, specificity and spatial resolution, as evidenced by multiple criteria including motif enrichment, sequence conservation, direct sequence pileup, nucleosome positioning and open chromatin states. In addition, we show that the fundamental advance of MACE is the identification of two boundaries of a TFBS with high resolution, whereas other methods only report a single location of the same event. The two boundaries help elucidate the in vivo binding structure of a given TF, e.g. whether the TF may bind as dimers or in a complex with other co-factors.

SeAMotE: a method for high-throughput motif discovery in nucleic acid sequences

Background

The large amount of data produced by high-throughput sequencing poses new computational challenges. In the last decade, several tools have been developed for the identification of transcription and splicing factor binding sites.

Results

Here, we introduce the SeAMotE (Sequence Analysis of Motifs Enrichment) algorithm for discovery of regulatory regions in nucleic acid sequences. SeAMotE provides (i) a robust analysis of high-throughput sequence sets, (ii) a motif search based on pattern occurrences and (iii) an easy-to-use web-server interface. We applied our method to recently published data including 351 chromatin immunoprecipitation (ChIP) and 13 crosslinking immunoprecipitation (CLIP) experiments and compared our results with those of other well-established motif discovery tools. SeAMotE shows an average accuracy of 80% in finding discriminative motifs and outperforms other methods available in literature.

Conclusions

SeAMotE is a fast, accurate and flexible algorithm for the identification of sequence patterns involved in protein-DNA and protein-RNA recognition. The server can be freely accessed at http://s.tartaglialab.com/new_submission/seamote.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-925) contains supplementary material, which is available to authorized users.

Functional annotation signatures of disease susceptibility loci improve SNP association analysis

BACKGROUND

Genetic association studies are conducted to discover genetic loci that contribute to an inherited trait, identify the variants behind these associations and ascertain their functional role in determining the phenotype. To date, functional annotations of the genetic variants have rarely played more than an indirect role in assessing evidence for association. Here, we demonstrate how these data can be systematically integrated into an association study's analysis plan.

RESULTS

We developed a Bayesian statistical model for the prior probability of phenotype-genotype association that incorporates data from past association studies and publicly available functional annotation data regarding the susceptibility variants under study. The model takes the form of a binary regression of association status on a set of annotation variables whose coefficients were estimated through an analysis of associated SNPs in the GWAS Catalog (GC). The functional predictors examined included measures that have been demonstrated to correlate with the association status of SNPs in the GC and some whose utility in this regard is speculative: summaries of the UCSC Human Genome Browser ENCODE super-track data, dbSNP function class, sequence conservation summaries, proximity to genomic variants in the Database of Genomic Variants and known regulatory elements in the Open Regulatory Annotation database, PolyPhen-2 probabilities and RegulomeDB categories. Because we expected that only a fraction of the annotations would contribute to predicting association, we employed a penalized likelihood method to reduce the impact of non-informative predictors and evaluated the model's ability to predict GC SNPs not used to construct the model. We show that the functional data alone are predictive of a SNP's presence in the GC. Further, using data from a genome-wide study of ovarian cancer, we demonstrate that their use as prior data when testing for association is practical at the genome-wide scale and improves power to detect associations.

CONCLUSIONS

We show how diverse functional annotations can be efficiently combined to create 'functional signatures' that predict the a priori odds of a variant's association to a trait and how these signatures can be integrated into a standard genome-wide-scale association analysis, resulting in improved power to detect truly associated variants.