Cistrome Data Browser: expanded datasets and new tools for gene regulatory analysis

C/EBPβ promotes poly(ADP-ribose) polymerase inhibitor resistance by enhancing homologous recombination repair in high-grade serous ovarian cancer

PARP inhibitors (PARPi) are efficacious in treating high-grade serous ovarian cancer (HG-SOC) with homologous recombination (HR) deficiency. However, they exhibit suboptimal efficiency in HR-proficient cancers. Here, we found that the expression of CCAAT/enhancer-binding protein β (C/EBPβ), a transcription factor, was inversely correlated with PARPi sensitivity in vitro and in vivo, both in HR-proficient condition. High C/EBPβ expression enhanced PARPi tolerance; PARPi treatment in turn induced C/EBPβ expression. C/EBPβ directly targeted and upregulated multiple HR genes (BRCA1, BRIP1, BRIT1, and RAD51), thereby inducing restoration of HR capacity and mediating acquired PARPi resistance. C/EBPβ is a key regulator of the HR pathway and an indicator of PARPi responsiveness. Targeting C/EBPβ could induce HR deficiency and rescue PARPi sensitivity accordingly. Our findings indicate that HR-proficient patients may benefit from PARPi via targeting C/EBPβ, and C/EBPβ expression levels enable predicting and tracking PARPi responsiveness during treatment.

FKBP-CaN-NFAT pathway polymorphisms selected by in silico biological function prediction are associated with tacrolimus efficacy in renal transplant patients

AIM

The aim of the present study was to investigate the potential effects of genetic variations in the FKBP-CaN-NFAT pathway on clinical events associated with tacrolimus efficacy in Chinese renal transplant patients.

METHODS

One hundred and forty Chinese renal transplant patients of Han ethnicity with over five years of follow-up were enrolled in our study. A pool of single nucleotide polymorphisms (SNPs) (1284 SNPs) was extracted from the Ensembl database according to chromosomal regions of the candidate genes. Next, 109 SNPs were screened out from this pool using multiple bioinformatics tools for subsequent genotyping using the MALDI-TOF-MS method. The associations of these candidate SNPs with acute rejection, nephrotoxicity, pneumonia and post-transplant estimated glomerular filtration rate (eGFR) were explored.

RESULTS

Fourty-four SNPs were found to be associated with tacrolimus-related clinical drug response. Specifically, eight SNPs were associated with the incidence of biopsy-proven acute rejection, four SNPs were associated with the rate of nephrotoxicity, 16 SNPs were correlated with the onset of pneumonia, and 26 SNPs were found to significantly influence post-transplant eGFR trend. An elaborate scoring system was implemented to prioritize the validation of these potentially causal SNPs. In particular, NFATC2 rs150348438 (G>T) performed well during integrative scoring (P_total=23.8) and was significantly associated with the occurrence of pneumonia (P = 0.0035, HR=0.91, 95% CI=0.85-0.97) and post-transplant eGFR levels (P = 0.000003).

CONCLUSIONS

NFATC2 rs150348438, rs6013219, rs1052653, and NFATC1 rs754093, ranking high in scoring, significantly affected the post-transplant eGFR and the incidence of pneumonia, acute rejection, and nephrotoxicity in renal transplant patients taking tacrolimus. Those SNPs may alter the expression and regulation of FKBP-CaN-NFAT pathway by influencing transcription regulation, mature mRNA degradation and RNA splicing, or protein coding. Critical SNPs of high ranking may serve as PD-associated pharmacogenetic biomarkers indicating individual response variability of TAC, and thus aid the clinical management of renal transplant patients.

SPC25 overexpression promotes tumor proliferation and is prognostic of poor survival in hepatocellular carcinoma

Background: The nuclear division cycle 80 (NDC80) complex assures proper chromosome segregation during the cell cycle progression. SPC25 is a crucial component of NDC80, and its role in hepatocellular carcinoma (HCC) has been explored recently. This study characterized the differential expression of SPC25 in HCC patients of different races and HBV infection status.

Methods: Expression patterns of SPC25 were evaluated in TCGA and Chinese HCC patients. Kaplan-Meier analysis was applied to examine the predictive value of SPC25. In vitro and in vivo functional assays were conducted to explore the role of SPC25 in HCC. Bioinformatics methods were applied to investigate the regulatory mechanisms of SPC25.

Findings: The mRNA levels of SPC25 were up-regulated in HCC. SPC25 has a significantly higher transcriptional level in Asian patients than Caucasian patients. SPC25 promoted HCC cell proliferation in vitro and tumor growth in vivo by accelerating the cell cycle. We identified transcription factors, miRNAs, and immune cells that may interact with SPC25.

Interpretation: The findings suggest that increased expression of SPC25 is associated with poor prognosis of HCC and enhances the proliferative capacity of HCC cells. SPC25 could serve as a valuable prognostic marker and a novel treatment target for HCC.

MFAP2 Promotes the Proliferation of Cancer Cells and Is Associated With a Poor Prognosis in Hepatocellular Carcinoma

Backgrounds:

Microfibril-associated protein 2 (MFAP2) is an extracellular matrix protein that regulates the function of microfibrils by interacting with fibrillin. MFAP2 has been reported to play an important role in metabolic diseases and has been shown to be significantly overexpressed in head and neck squamous cell carcinoma and Hepatocellular carcinoma (HCC). However, the molecular function and prognostic value of MFAP2 have never been reported in HCC or other tumors.

Methods:

In the present study, expression characteristics of MFAP2 in HCC, its influence on the development of HCC, as well as its function and potential mechanism in HCC were verified by Quantitative reverse transcription-polymerase chain reaction, bioinformatics data mining and in vitro cell experiments.

Results:

MFAP2 was prominently high-expressed in HCC and associated with cancer stages. HCC patients with higher MFAP2 expression displayed lower overall survival (OS) and disease-specific survival(DSS), while there was no significant difference in recurrence-free survival (RFS). In vitro experiments showed that downregulation of MFAP2 inhibited proliferation, migration level of HCC cells. Transcription factors, DNA methyltransferases, immune factors may interact with MFAP2 mRNA to promote tumor progression in HCC.

Conclusion:

These findings suggest that MFAP2 may play a key role in the development of HCC. Therefore, MFAP2 may be a valuable prognostic marker and an effective anticancer target in HCC.

Dissecting the DNA binding landscape and gene regulatory network of p63 and p53

The transcription factor p53 is the best-known tumor suppressor, but its sibling p63 is a master regulator of epidermis development and a key oncogenic driver in squamous cell carcinomas (SCC). Despite multiple gene expression studies becoming available, the limited overlap of reported p63-dependent genes has made it difficult to decipher the p63 gene regulatory network. Particularly, analyses of p63 response elements differed substantially among the studies. To address this intricate data situation, we provide an integrated resource that enables assessing the p63-dependent regulation of any human gene of interest. We use a novel iterative de novo motif search approach in conjunction with extensive ChIP-seq data to achieve a precise global distinction between p53-and p63-binding sites, recognition motifs, and potential co-factors. We integrate these data with enhancer:gene associations to predict p63 target genes and identify those that are commonly de-regulated in SCC representing candidates for prognosis and therapeutic interventions.

How to Use Online Tools to Generate New Hypotheses for Mammary Gland Biology Research: A Case Study for Wnt7b

An increasing number of '-omics' datasets, generated by labs all across the world, are becoming available. They contain a wealth of data that are largely unexplored. Not every scientist, however, will have access to the required resources and expertise to analyze such data from scratch. Fortunately, a growing number of investigators is dedicating their time and effort to the development of user friendly, online applications that allow researchers to use and investigate these datasets. Here, we will illustrate the usefulness of such an approach. Using regulation of Wnt7b expression as an example, we will highlight a selection of accessible tools and resources that are available to researchers in the area of mammary gland biology. We show how they can be used for in silico analyses of gene regulatory mechanisms, resulting in new hypotheses and providing leads for experimental follow up. We also call out to the mammary gland community to join forces in a coordinated effort to generate and share additional tissue-specific '-omics' datasets and thereby expand the in silico toolbox.

Comprehensive Mapping of Key Regulatory Networks that Drive Oncogene Expression

Gene expression is controlled by the collective binding of transcription factors to cis-regulatory regions. Deciphering gene-centered regulatory networks is vital to understanding and controlling gene misexpression in human disease; however, systematic approaches to uncovering regulatory networks have been lacking. Here we present high-throughput interrogation of gene-centered activation networks (HIGAN), a pipeline that employs a suite of multifaceted genomic approaches to connect upstream signaling inputs, trans-acting TFs, and cis-regulatory elements. We apply HIGAN to understand the aberrant activation of the cytidine deaminase APOBEC3B, an intrinsic source of cancer hypermutation. We reveal that nuclear factor κB (NF-κB) and AP-1 pathways are the most salient trans-acting inputs, with minor roles for other inflammatory pathways. We identify a cis-regulatory architecture dominated by a major intronic enhancer that requires coordinated NF-κB and AP-1 activity with secondary inputs from distal regulatory regions. Our data demonstrate how integration of cis and trans genomic screening platforms provides a paradigm for building gene-centered regulatory networks.

Genes Associated with Calcium Signaling are Involved in Alcohol-Induced Breast Cancer Growth

BACKGROUND

Alcohol consumption is a risk factor for breast cancer, contributing to up to nearly 23,000 new cases each year. Mechanistic studies show that alcohol increases tumor aggressiveness and metastatic potential, promotes angiogenesis, induces chronic inflammation, and dysregulates RNA polymerase III-related genes. Alcohol has also been shown to affect estrogen signaling in breast cancer, including in our study of the transcriptomic effects of alcohol in breast cancer cells.

METHODS

To elucidate mechanisms of action of alcohol in breast cancer, we carried out secondary analyses of our alcohol-responsive transcriptome data using gene ontology and pathway databases and analysis tools and cistromic data analysis of candidate transcription factors which may mediate the transcriptomic alterations. Predicted alcohol-responsive pathways and mechanisms were perturbed and examined experimentally in breast cancer cells. The clinical relevance of identified genes was determined by expression profiles in patient samples and correlation with disease outcomes and alcohol consumption in previously published study cohorts.

RESULTS

Gene ontology analysis showed that alcohol alters the expression of many metabolism-related genes, and cistromic data of differentially expressed genes revealed the potential involvement of nuclear factor of activated T cells 3 (NFATC3) in mediating the transcriptomic effects of alcohol. Pathway analysis also predicted regulation of calcium signaling by alcohol in breast cancer cells. Chemical perturbation of this pathway reversed the effect of alcohol on breast cancer cell growth and reduced the elevated cytosolic Ca²⁺ levels induced by alcohol. Expression levels of alcohol-responsive genes in tumor samples from breast cancer patients are associated with poor disease outcomes. Moreover, expression of some of these genes was altered in breast cancer patients who consumed alcohol previously as compared to those who did not drink.

CONCLUSION

Alcohol alters expression of genes that regulate intracellular calcium levels and downstream signaling pathways which drive breast cancer cell proliferation and disease progression.

GREG-studying transcriptional regulation using integrative graph databases

A gene regulatory process is the result of the concerted action of transcription factors, co-factors, regulatory non-coding RNAs (ncRNAs) and chromatin interactions. Therefore, the combination of protein-DNA, protein-protein, ncRNA-DNA, ncRNA-protein and DNA-DNA data in a single graph database offers new possibilities regarding generation of biological hypotheses. GREG (The Gene Regulation Graph Database) is an integrative database and web resource that allows the user to visualize and explore the network of all above-mentioned interactions for a query transcription factor, long non-coding RNA, genomic range or DNA annotation, as well as extracting node and interaction information, identifying connected nodes and performing advanced graphical queries directly on the regulatory network, in a simple and efficient way. In this article, we introduce GREG together with some application examples (including exploratory research of Nanog's regulatory landscape and the etiology of chronic obstructive pulmonary disease), which we use as a demonstration of the advantages of using graph databases in biomedical research. Database URL: https://mora-lab.github.io/projects/greg.html, www.moralab.science/GREG/.

An integrative multi-omics network-based approach identifies key regulators for breast cancer

Although genome-wide association studies (GWASs) have successfully identified thousands of risk variants for human complex diseases, understanding the biological function and molecular mechanisms of the associated SNPs involved in complex diseases is challenging. Here we developed a framework named integrative multi-omics network-based approach (IMNA), aiming to identify potential key genes in regulatory networks by integrating molecular interactions across multiple biological scales, including GWAS signals, gene expression-based signatures, chromatin interactions and protein interactions from the network topology. We applied this approach to breast cancer, and prioritized key genes involved in regulatory networks. We also developed an abnormal gene expression score (AGES) signature based on the gene expression deviation of the top 20 rank-ordered genes in breast cancer. The AGES values are associated with genetic variants, tumor properties and patient survival outcomes. Among the top 20 genes, RNASEH2A was identified as a new candidate gene for breast cancer. Thus, our integrative network-based approach provides a genetic-driven framework to unveil tissue-specific interactions from multiple biological scales and reveal potential key regulatory genes for breast cancer. This approach can also be applied in other complex diseases such as ovarian cancer to unravel underlying mechanisms and help for developing therapeutic targets.

Ultrafast and scalable variant annotation and prioritization with big functional genomics data

The advances of large-scale genomics studies have enabled compilation of cell type–specific, genome-wide DNA functional elements at high resolution. With the growing volume of functional annotation data and sequencing variants, existing variant annotation algorithms lack the efficiency and scalability to process big genomic data, particularly when annotating whole-genome sequencing variants against a huge database with billions of genomic features. Here, we develop VarNote to rapidly annotate genome-scale variants in large and complex functional annotation resources. Equipped with a novel index system and a parallel random-sweep searching algorithm, VarNote shows substantial performance improvements (two to three orders of magnitude) over existing algorithms at different scales. It supports both region-based and allele-specific annotations and introduces advanced functions for the flexible extraction of annotations. By integrating massive base-wise and context-dependent annotations in the VarNote framework, we introduce three efficient and accurate pipelines to prioritize the causal regulatory variants for common diseases, Mendelian disorders, and cancers.

Clonal tracing reveals diverse patterns of response to immune checkpoint blockade

BACKGROUND

Immune checkpoint blockade (ICB) therapy has improved patient survival in a variety of cancers, but only a minority of cancer patients respond. Multiple studies have sought to identify general biomarkers of ICB response, but elucidating the molecular and cellular drivers of resistance for individual tumors remains challenging. We sought to determine whether a tumor with defined genetic background exhibits a stereotypic or heterogeneous response to ICB treatment.

RESULTS

We establish a unique mouse system that utilizes clonal tracing and mathematical modeling to monitor the growth of each cancer clone, as well as the bulk tumor, in response to ICB. We find that tumors derived from the same clonal populations showed heterogeneous ICB response and diverse response patterns. Primary response is associated with higher immune infiltration and leads to enrichment of pre-existing ICB-resistant cancer clones. We further identify several cancer cell-intrinsic gene expression signatures associated with ICB resistance, including increased interferon response genes and glucocorticoid response genes. These findings are supported by clinical data from ICB treatment cohorts.

CONCLUSIONS

Our study demonstrates diverse response patterns from the same ancestor cancer cells in response to ICB. This suggests the value of monitoring clonal constitution and tumor microenvironment over time to optimize ICB response and to design new combination therapies. Furthermore, as ICB response may enrich for cancer cell-intrinsic resistance signatures, this can affect interpretations of tumor RNA-seq data for response-signature association studies.

EZH2 overexpression dampens tumor-suppressive signals via an EGR1 silencer to drive breast tumorigenesis

The mechanism underlying EZH2 overexpression in breast cancer and its involvement in tumorigenesis remain poorly understood. In this study, we developed an approach to systematically identify the trans-acting factors regulating the EZH2 expression, and identified more than 20 such factors. We revealed reciprocal regulation of early growth response 1 (EGR1) and EZH2: EGR1 activates the expression of EZH2, and EZH2 represses EGR1 expression. Using CRISPR-mediated genome/epigenome editing, we demonstrated that EHZ2 represses EGR1 expression through a silencer downstream of the EGR1 gene. Deletion of the EGR1 silencer resulted in reduced cell growth, invasion, tumorigenicity of breast cancer cells, and extensive changes in gene expression, such as upregulation of GADD45, DDIT3, and RND1; and downregulation of genes encoding cholesterol biosynthesis pathway enzymes. We hypothesize that EZH2/PRC2 acts as a "brake" for EGR1 expression by targeting the EGR1 silencer, and EZH2 overexpression dampens tumor-suppressive signals mediated by EGR1 to drive breast tumorigenesis.

Tissue-specific usage of transposable element-derived promoters in mouse development

BACKGROUND

Transposable elements (TEs) are a significant component of eukaryotic genomes and play essential roles in genome evolution. Mounting evidence indicates that TEs are highly transcribed in early embryo development and contribute to distinct biological functions and tissue morphology.

RESULTS

We examine the epigenetic dynamics of mouse TEs during the development of five tissues: intestine, liver, lung, stomach, and kidney. We found that TEs are associated with over 20% of open chromatin regions during development. Close to half of these accessible TEs are only activated in a single tissue and a specific developmental stage. Most accessible TEs are rodent-specific. Across these five tissues, 453 accessible TEs are found to create the transcription start sites of downstream genes in mouse, including 117 protein-coding genes and 144 lincRNA genes, 93.7% of which are mouse-specific. Species-specific TE-derived transcription start sites are found to drive the expression of tissue-specific genes and change their tissue-specific expression patterns during evolution.

CONCLUSION

Our results suggest that TE insertions increase the regulatory potential of the genome, and some TEs have been domesticated to become a crucial component of gene and regulate tissue-specific expression during mouse tissue development.

Systematic screening of CTCF binding partners identifies that BHLHE40 regulates CTCF genome-wide distribution and long-range chromatin interactions

CTCF plays a pivotal role in mediating chromatin interactions, but it does not do so alone. A number of factors have been reported to co-localize with CTCF and regulate CTCF loops, but no comprehensive analysis of binding partners has been performed. This prompted us to identify CTCF loop participants and regulators by co-localization analysis with CTCF. We screened all factors that had ChIP-seq data in humans by co-localization analysis with human super conserved CTCF (hscCTCF) binding sites, and identified many new factors that overlapped with hscCTCF binding sites. Combined with CTCF loop information, we observed that clustered factors could promote CTCF loops. After in-depth mining of each factor, we found that many factors might have the potential to promote CTCF loops. Our data further demonstrated that BHLHE40 affected CTCF loops by regulating CTCF binding. Together, this study revealed that many factors have the potential to participate in or regulate CTCF loops, and discovered a new role for BHLHE40 in modulating CTCF loop formation.

EAT-UpTF: Enrichment Analysis Tool for Upstream Transcription Factors of a Group of Plant Genes

EAT-UpTF (Enrichment Analysis Tool for Upstream Transcription Factors of a group of plant genes) is an open-source Python script that analyzes the enrichment of upstream transcription factors (TFs) in a group of genes-of-interest (GOIs). EAT-UpTF utilizes genome-wide lists of TF-target genes generated by DNA affinity purification followed by sequencing (DAP-seq) or chromatin immunoprecipitation followed by sequencing (ChIP-seq). Unlike previous methods based on the two-step prediction of cis-motifs and DNA-element-binding TFs, our EAT-UpTF analysis enabled a one-step identification of enriched upstream TFs in a set of GOIs using lists of empirically determined TF-target genes. The tool is designed particularly for plant researches, due to the lack of analytic tools for upstream TF enrichment, and available at https://github.com/sangreashim/EAT-UpTF and http://chromatindynamics.snu.ac.kr:8080/EatupTF.

STARR-seq identifies active, chromatin-masked, and dormant enhancers in pluripotent mouse embryonic stem cells

BACKGROUND

Enhancers are distal regulators of gene expression that shape cell identity and control cell fate transitions. In mouse embryonic stem cells (mESCs), the pluripotency network is maintained by the function of a complex network of enhancers, that are drastically altered upon differentiation. Genome-wide chromatin accessibility and histone modification assays are commonly used as a proxy for identifying putative enhancers and for describing their activity levels and dynamics.

RESULTS

Here, we applied STARR-seq, a genome-wide plasmid-based assay, as a read-out for the enhancer landscape in "ground-state" (2i+LIF; 2iL) and "metastable" (serum+LIF; SL) mESCs. This analysis reveals that active STARR-seq loci show modest overlap with enhancer locations derived from peak calling of ChIP-seq libraries for common enhancer marks. We unveil ZIC3-bound loci with significant STARR-seq activity in SL-ESCs. Knock-out of Zic3 removes STARR-seq activity only in SL-ESCs and increases their propensity to differentiate towards the endodermal fate. STARR-seq also reveals enhancers that are not accessible, masked by a repressive chromatin signature. We describe a class of dormant, p53 bound enhancers that gain H3K27ac under specific conditions, such as after treatment with Nocodazol, or transiently during reprogramming from fibroblasts to pluripotency.

CONCLUSIONS

In conclusion, loci identified as active by STARR-seq often overlap with those identified by chromatin accessibility and active epigenetic marking, yet a significant fraction is epigenetically repressed or display condition-specific enhancer activity.

Etv1 Controls the Establishment of Non-overlapping Motor Innervation of Neighboring Facial Muscles during Development

The somatotopic motor-neuron projections onto their cognate target muscles are essential for coordinated movement, but how that occurs for facial motor circuits, which have critical roles in respiratory and interactive behaviors, is poorly understood. We report extensive molecular heterogeneity in developing facial motor neurons in the mouse and identify markers of subnuclei and the motor pools innervating specific facial muscles. Facial subnuclei differentiate during migration to the ventral hindbrain, where neurons with progressively later birth dates—and evolutionarily more recent functions—settle in more-lateral positions. One subpopulation marker, ETV1, determines both positional and target muscle identity for neurons of the dorsolateral (DL) subnucleus. In Etv1 mutants, many markers of DL differentiation are lost, and individual motor pools project indifferently to their own and neighboring muscle targets. The resulting aberrant activation patterns are reminiscent of the facial synkinesis observed in humans after facial nerve injury.

Tenney et al. demonstrate that embryonic facial motor neurons are transcriptionally diverse as they establish somatotopic innervation of the facial muscles, a process that requires the transcription factor ETV1. Facial-motor axon-targeting errors in Etv1 mutants cause coordination of whisking and eyeblink evocative of human blepharospasm.

Calcium signaling and epigenetics: A key point to understand carcinogenesis

Calcium (Ca²⁺) signaling controls a wide range of cellular processes, including the hallmarks of cancer. The Ca²⁺ signaling system encompasses several types of proteins, such as receptors, channels, pumps, exchangers, buffers, and sensors, of which several are mutated or with altered expression in cancer cells. Since epigenetic mechanisms are disrupted in all stages of carcinogenesis, and reversibly regulate gene expression, they have been studied by different research groups to understand their role in Ca²⁺ signaling remodeling in cancer cells and the carcinogenic process. In this review, we link Ca²⁺ signaling, cancer, and epigenetics fields to generate a comprehensive landscape of this complex group of diseases.

LncRNA DSCAM-AS1 interacts with YBX1 to promote cancer progression by forming a positive feedback loop that activates FOXA1 transcription network

Rationale: The forkhead box A1 (FOXA1) is a crucial transcription factor in initiation and development of breast, lung and prostate cancer. Previous studies about the FOXA1 transcriptional network were mainly focused on protein-coding genes. Its regulatory network of long non-coding RNAs (lncRNAs) and their role in FOXA1 oncogenic activity remains unknown.

Methods: The Cancer Genome Atlas (TCGA) data, RNA-seq and ChIP-seq data were used to analyze FOXA1 regulated lncRNAs. RT-qPCR was used to detect the expression of DSCAM-AS1, RT-qPCR and Western blotting were used to determine the expression of FOXA1, estrogen receptor α (ERα) and Y box binding protein 1 (YBX1). RNA pull-down and RIP-qPCR were employed to investigate the interaction between DSCAM-AS1 and YBX1. The effect of DSCAM-AS1 on malignant phenotypes was examined through in vitro and in vivo assays.

Results: In this study, we conducted a global analysis of FOXA1 regulated lncRNAs. For detailed analysis, we chose lncRNA DSCAM-AS1, which is specifically expressed in lung adenocarcinoma, breast and prostate cancer. The expression level of DSCAM-AS1 is regulated by two super-enhancers (SEs) driven by FOXA1. High expression levels of DSCAM-AS1 was associated with poor prognosis. Knockout experiments showed DSCAM-AS1 was essential for the growth of xenograft tumors. Moreover, we demonstrated DSCAM-AS1 can regulate the expression of the master transcriptional factor FOXA1. In breast cancer, DSCAM-AS1 was also found to regulate ERα. Mechanistically, DSCAM-AS1 interacts with YBX1 and influences the recruitment of YBX1 in the promoter regions of FOXA1 and ERα.

Conclusion: Our study demonstrated that lncRNA DSCAM-AS1 was transcriptionally activated by super-enhancers driven by FOXA1 and exhibited lineage-specific expression pattern. DSCAM-AS1 can promote cancer progression by interacting with YBX1 and regulating expression of FOXA1 and ERα.

Multiomics data integration unveils core transcriptional regulatory networks governing cell-type identity

A plethora of computational approaches have been proposed for reconstructing gene regulatory networks (GRNs) from gene expression data. However, gene regulatory processes are often too complex to predict from the transcriptome alone. Here, we present a computational method, Moni, that systematically integrates epigenetics, transcriptomics, and protein–protein interactions to reconstruct GRNs among core transcription factors and their co-factors governing cell identity. We applied Moni to 57 datasets of human cell types and lines and demonstrate that it can accurately infer GRNs, thereby outperforming state-of-the-art methods.

Mapping the Effects of Genetic Variation on Chromatin State and Gene Expression Reveals Loci That Control Ground State Pluripotency

Mouse embryonic stem cells (mESCs) cultured in the presence of LIF occupy a ground state with highly active pluripotency-associated transcriptional and epigenetic circuitry. However, ground state pluripotency in some inbred strain backgrounds is unstable in the absence of ERK1/2 and GSK3 inhibition. Using an unbiased genetic approach, we dissect the basis of this divergent response to extracellular cues by profiling gene expression and chromatin accessibility in 170 genetically heterogeneous mESCs. We map thousands of loci affecting chromatin accessibility and/or transcript abundance, including 10 QTL hotspots where genetic variation at a single locus coordinates the regulation of genes throughout the genome. For one hotspot, we identify a single enhancer variant ∼10 kb upstream of Lifr associated with chromatin accessibility and mediating a cascade of molecular events affecting pluripotency. We validate causation through reciprocal allele swaps, demonstrating the functional consequences of noncoding variation in gene regulatory networks that stabilize pluripotent states in vitro.

Integrative analyses of single-cell transcriptome and regulome using MAESTRO

We present Model-based AnalysEs of Transcriptome and RegulOme (MAESTRO), a comprehensive open-source computational workflow ( http://github.com/liulab-dfci/MAESTRO ) for the integrative analyses of single-cell RNA-seq (scRNA-seq) and ATAC-seq (scATAC-seq) data from multiple platforms. MAESTRO provides functions for pre-processing, alignment, quality control, expression and chromatin accessibility quantification, clustering, differential analysis, and annotation. By modeling gene regulatory potential from chromatin accessibilities at the single-cell level, MAESTRO outperforms the existing methods for integrating the cell clusters between scRNA-seq and scATAC-seq. Furthermore, MAESTRO supports automatic cell-type annotation using predefined cell type marker genes and identifies driver regulators from differential scRNA-seq genes and scATAC-seq peaks.

A novel risk score based on a combined signature of 10 immune system genes to predict bladder cancer prognosis

Bladder cancer (BC) is a common internal malignant tumor with a poor prognosis worldwide. There is an urgent need to better understand the pathogenesis and progression of BC and to find useful biomarkers for diagnosis and prognosis. This study was aimed at developing a potential immunogenomic prognostic signature for BC patients. To identify possible immune-system-related genes (IRGs) whose parameters predict the survival of BC patients, we chose 371 BC patients and analyzed differentially expressed IRGs from The Cancer Genome Atlas (TCGA) datasets. We then derived a 10-IRG formula, including MMP9, RBP7, PDGFRA, AHNAK, OAS1, OLR1, RAC3, SLIT2, IGF1, and AGTR1, to estimate BC prognosis. To validate the mRNA levels of these IRGs, we performed quantitative PCR and found that the expression of these genes almost matched the corresponding mRNA expression levels in TCGA. Furthermore, we validated the prognostic value of the new risk model using two external datasets from Gene Expression Omnibus: GSE13507 (n = 165) and GSE32894 (n = 224). Our data pointed to a significant correlation between the risk model and patients' prognosis. Bioinformatic analysis revealed that products of the IRGs have possible effects on tumor immune processes such as an inflammatory response and cytokine-cytokine receptor interaction. Finally, assessment of the clinical value of the immune-system-based risk signature showed that several of these IRGs were differentially expressed between patients with different clinical characteristics: a high risk score positively correlated with female sex, advanced tumor stage, more advanced T stage, and lymph node metastasis. This immunogenomic signature may represents a reliable prognostic tool for BC and can help to design an individualized immunotherapy.

The hepatocyte clock and feeding control chronophysiology of multiple liver cell types

Most cells of the body contain molecular clocks, but the requirement of peripheral clocks for rhythmicity and their effects on physiology are not well understood. We show that deletion of core clock components REV-ERBα and REV-ERBβ in adult mouse hepatocytes disrupts diurnal rhythms of a subset of liver genes and alters the diurnal rhythm of de novo lipogenesis. Liver function is also influenced by nonhepatocytic cells, and the loss of hepatocyte REV-ERBs remodels the rhythmic transcriptomes and metabolomes of multiple cell types within the liver. Finally, alteration of food availability demonstrates the hierarchy of the cell-intrinsic hepatocyte clock mechanism and the feeding environment. Together, these studies reveal previously unsuspected roles of the hepatocyte clock in the physiological coordination of nutritional signals and cell-cell communication controlling rhythmic metabolism.

Inter-Strain Epigenomic Profiling Reveals a Candidate IAP Master Copy in C3H Mice

Insertions of endogenous retroviruses cause a significant fraction of mutations in inbred mice but not all strains are equally susceptible. Notably, most new Intracisternal A particle (IAP) ERV mutagenic insertions have occurred in C3H mice. We show here that strain-specific insertional polymorphic IAPs accumulate faster in C3H/HeJ mice, relative to other sequenced strains, and that IAP transcript levels are higher in C3H/HeJ embryonic stem (ES) cells compared to other ES cells. To investigate the mechanism for high IAP activity in C3H mice, we identified 61 IAP copies in C3H/HeJ ES cells enriched with H3K4me3 (a mark of active promoters) and, among those tested, all are unmethylated in C3H/HeJ ES cells. Notably, 13 of the 61 are specific to C3H/HeJ and are members of the non-autonomous 1Δ1 IAP subfamily that is responsible for nearly all new insertions in C3H. One copy is full length with intact open reading frames and hence potentially capable of providing proteins in trans to other 1Δ1 elements. This potential “master copy” is present in other strains, including 129, but its 5’ long terminal repeat (LTR) is methylated in 129 ES cells. Thus, the unusual IAP activity in C3H may be due to reduced epigenetic repression coupled with the presence of a master copy.

Cohesin-Dependent and -Independent Mechanisms Mediate Chromosomal Contacts between Promoters and Enhancers

It is currently assumed that 3D chromosomal organization plays a central role in transcriptional control. However, depletion of cohesin and CTCF affects the steady-state levels of only a minority of transcripts. Here, we use high-resolution Capture Hi-C to interrogate the dynamics of chromosomal contacts of all annotated human gene promoters upon degradation of cohesin and CTCF. We show that a majority of promoter-anchored contacts are lost in these conditions, but many contacts with distinct properties are maintained, and some new ones are gained. The rewiring of contacts between promoters and active enhancers upon cohesin degradation associates with rapid changes in target gene transcription as detected by SLAM sequencing (SLAM-seq). These results provide a mechanistic explanation for the limited, but consistent, effects of cohesin and CTCF depletion on steady-state transcription and suggest the existence of both cohesin-dependent and -independent mechanisms of enhancer-promoter pairing.

Epigenetic profiling demarcates molecular subtypes of muscle-invasive bladder cancer

Muscle-invasive bladder cancer (MIBC) is a heterogeneous disease that often recurs despite aggressive treatment with neoadjuvant chemotherapy and (radical) cystectomy. Basal and luminal molecular subtypes have been identified that are linked to clinical characteristics and have differential sensitivities to chemotherapy. While it has been suggested that epigenetic mechanisms play a role in defining these subtypes, a thorough understanding of the biological mechanisms is lacking. This report details the first genome-wide analysis of histone methylation patterns of human primary bladder tumours by chromatin immunoprecipitations and next-generation sequencing (ChIP-seq). We profiled multiple histone marks: H3K27me3, a marker for repressed genes, and H3K4me1 and H3K4me3, which are indicators of active enhancers and active promoters. Integrated analysis of ChIP-seq data and RNA sequencing revealed that H3K4 mono-methylation demarcates MIBC subtypes, while no association was found for the other two histone modifications in relation to basal and luminal subtypes. Additionally, we identified differentially methylated H3K4me1 peaks in basal and luminal tumour samples, suggesting that active enhancers play a role in defining subtypes. Our study is the first analysis of histone modifications in primary bladder cancer tissue and provides an important resource for the bladder cancer community.

Prediction and identification of immune genes related to the prognosis of patients with colon adenocarcinoma and its mechanisms

Background

Colon adenocarcinoma (COAD) is a gastrointestinal tumor with a high degree of malignancy. Its deterioration process is closely related to the tumor microenvironment, and transcription factors (TF) play a regulatory role in this process. Currently, there is a lack of exploration between the genes related to the COAD tumor microenvironment and the survival prognosis of patients. Models composed of multiple genes usually predict the survival prognosis of patients more accurately than single genes. We can analyze the multigene models that can predict the prognosis of COAD from the current database.

Methods

The limma package of the R programming language is used for gene differential expression analysis. Kaplan-Meier curve is used to analyze the relationship between the patient risk score model and survival data. The hazard model is used to analyze the relationship between the risk score and the clinical data of COAD patients. The information of immune genes and immune cells is obtained from IMMPORT database and TIMER database. Receiver operating characteristic (ROC) curve is used to judge the stability of the model.

Results

We found 7 immune genes, which can built a risk score model to predict the survival prognosis of COAD. According to univariate and multivariate analysis, the risk score can be used as an independent predictor. The content of some immune microenvironment cells will also increase as the risk score increases.

Conclusions

We found 7 immune genes, such as SLC10A2 (solute carrier family 10 member 2), CXCL3 (C-X-C motif chemokine ligand 3), IGHV5-51 (immunoglobulin heavy variable 5-51), INHBA (inhibin subunit beta A), STC1 (stanniocalcin 1), UCN (urocortin), and OXTR (oxytocin receptor), can constitute a model for predicting the prognosis of COAD. They may provide potential therapeutic targets for clinical treatment of COAD.

Intragenic recruitment of NF-κB drives splicing modifications upon activation by the oncogene Tax of HTLV-1

Chronic NF-κB activation in inflammation and cancer has long been linked to persistent activation of NF-κB–responsive gene promoters. However, NF-κB factors also massively bind to gene bodies. Here, we demonstrate that recruitment of the NF-κB factor RELA to intragenic regions regulates alternative splicing upon NF-κB activation by the viral oncogene Tax of HTLV-1. Integrative analyses of RNA splicing and chromatin occupancy, combined with chromatin tethering assays, demonstrate that DNA-bound RELA interacts with and recruits the splicing regulator DDX17, in an NF-κB activation-dependent manner. This leads to alternative splicing of target exons due to the RNA helicase activity of DDX17. Similar results were obtained upon Tax-independent NF-κB activation, indicating that Tax likely exacerbates a physiological process where RELA provides splice target specificity. Collectively, our results demonstrate a physical and direct involvement of NF-κB in alternative splicing regulation, which significantly revisits our knowledge of HTLV-1 pathogenesis and other NF-κB-related diseases.

The nuclear factors κB (NF-κB) is a transcription factor involved in immune functions, inflammation, and cancer. Here, the authors show that the NF-κB factor RELA regulates splicing of target genes by recruiting DDX17 on chromatin upon expression of the viral oncogene Tax.

Folate Receptor β (FRβ) Expression in Tissue-Resident and Tumor-Associated Macrophages Associates with and Depends on the Expression of PU.1

As macrophages exhibit a huge functional plasticity under homeostasis and pathological conditions, they have become a therapeutic target for chronic inflammatory diseases. Hence, the identification of macrophage subset-specific markers is a requisite for the development of macrophage-directed therapeutic interventions. In this regard, the macrophage-specific Folate Receptor β (FRβ, encoded by the FOLR2 gene) has been already validated as a target for molecular delivery in cancer as well as in macrophage-targeting therapeutic strategies for chronic inflammatory pathologies. We now show that the transcriptome of human macrophages from healthy and inflamed tissues (tumor; rheumatoid arthritis, RA) share a significant over-representation of the “anti-inflammatory gene set”, which defines the gene profile of M-CSF-dependent IL-10-producing human macrophages (M-MØ). More specifically, FOLR2 expression has been found to strongly correlate with the expression of M-MØ-specific genes in tissue-resident macrophages, tumor-associated macrophages (TAM) and macrophages from inflamed synovium, and also correlates with the presence of the PU.1 transcription factor. In fact, PU.1-binding elements are found upstream of the first exon of FOLR2 and most M-MØ-specific- and TAM-specific genes. The functional relevance of PU.1 binding was demonstrated through analysis of the proximal regulatory region of the FOLR2 gene, whose activity was dependent on a cluster of PU.1-binding sequences. Further, siRNA-mediated knockdown established the importance of PU.1 for FOLR2 gene expression in myeloid cells. Therefore, we provide evidence that FRβ marks tissue-resident macrophages as well as macrophages within inflamed tissues, and its expression is dependent on PU.1.

The road towards data integration in human genomics: players, steps and interactions

Thousands of new experimental datasets are becoming available every day; in many cases, they are produced within the scope of large cooperative efforts, involving a variety of laboratories spread all over the world, and typically open for public use. Although the potential collective amount of available information is huge, the effective combination of such public sources is hindered by data heterogeneity, as the datasets exhibit a wide variety of notations and formats, concerning both experimental values and metadata. Thus, data integration is becoming a fundamental activity, to be performed prior to data analysis and biological knowledge discovery, consisting of subsequent steps of data extraction, normalization, matching and enrichment; once applied to heterogeneous data sources, it builds multiple perspectives over the genome, leading to the identification of meaningful relationships that could not be perceived by using incompatible data formats. In this paper, we first describe a technological pipeline from data production to data integration; we then propose a taxonomy of genomic data players (based on the distinction between contributors, repository hosts, consortia, integrators and consumers) and apply the taxonomy to describe about 30 important players in genomic data management. We specifically focus on the integrator players and analyse the issues in solving the genomic data integration challenges, as well as evaluate the computational environments that they provide to follow up data integration by means of visualization and analysis tools.

Downexpression of HSD17B6 correlates with clinical prognosis and tumor immune infiltrates in hepatocellular carcinoma

Background

Hydroxysteroid 17-Beta Dehydrogenase 6 (HSD17B6), a key protein involved in synthetizing dihydrotestosterone, is abundant in the liver. Previous studies have suggested a role for dihydrotestosterone in modulating progress of various malignancies, and HSD17B6 dysfunction was associated with lung cancer and prostate cancer. However, little is known about the detailed role of HSD17B6 in hepatocellular carcinoma (HCC).

Methods

Clinical implication and survival data related to HSD17B6 expression in patients with HCC were obtained through TCGA, ICGC, ONCOMINE, GEO and HPA databases. Survival analysis plots were drawn with Kaplan–Meier Plotter. The ChIP-seq data were obtained from Cistrome DB. Protein–Protein Interaction and gene functional enrichment analyses were performed in STRING database. The correlations between HSD17B6 and tumor immune infiltrates was investigated via TIMER and xCell. The proliferation, migration and invasion of liver cancer cells transfected with HSD17B6 were evaluated by the CCK8 assay, wound healing test and transwell assay respectively. Expression of HSD17B6, TGFB1 and PD-L1 were assessed by quantitative RT-PCR.

Results

HSD17B6 expression was lower in HCC compared to normal liver and correlated with tumor stage and grade. Lower expression of HSD17B6 was associated with worse OS, PFS, RFS and DSS in HCC patients. HNF4A bound to enhancer and promoter regions of HSD17B6 gene, activating its transcription, and DNA methylation of HSD17B6 promoter negatively controlled the expression. HSD17B6 and its interaction partners were involved in androgen metabolism and biosynthesis in liver. HSD17B6 inhibited tumor cell proliferation, migration and invasion in liver cancer cells and low expression of HSD17B6 correlated with high immune cells infiltration, relative reduction of immune responses and multiple immune checkpoint genes expression in HCC, probably by regulating the expression of TGFB1.

Conclusions

This study indicate that HSD17B6 could be a new biomarker for the prognosis of HCC and an important negative regulator of immune responses in HCC.

CAUSALdb: a database for disease/trait causal variants identified using summary statistics of genome-wide association studies

Genome-wide association studies (GWASs) have revolutionized the field of complex trait genetics over the past decade, yet for most of the significant genotype-phenotype associations the true causal variants remain unknown. Identifying and interpreting how causal genetic variants confer disease susceptibility is still a big challenge. Herein we introduce a new database, CAUSALdb, to integrate the most comprehensive GWAS summary statistics to date and identify credible sets of potential causal variants using uniformly processed fine-mapping. The database has six major features: it (i) curates 3052 high-quality, fine-mappable GWAS summary statistics across five human super-populations and 2629 unique traits; (ii) estimates causal probabilities of all genetic variants in GWAS significant loci using three state-of-the-art fine-mapping tools; (iii) maps the reported traits to a powerful ontology MeSH, making it simple for users to browse studies on the trait tree; (iv) incorporates highly interactive Manhattan and LocusZoom-like plots to allow visualization of credible sets in a single web page more efficiently; (v) enables online comparison of causal relations on variant-, gene- and trait-levels among studies with different sample sizes or populations and (vi) offers comprehensive variant annotations by integrating massive base-wise and allele-specific functional annotations. CAUSALdb is freely available at http://mulinlab.org/causaldb.

ERα-related chromothripsis enhances concordant gene transcription on chromosome 17q11.1-q24.1 in luminal breast cancer

BACKGROUND

Chromothripsis is an event of genomic instability leading to complex chromosomal alterations in cancer. Frequent long-range chromatin interactions between transcription factors (TFs) and targets may promote extensive translocations and copy-number alterations in proximal contact regions through inappropriate DNA stitching. Although studies have proposed models to explain the initiation of chromothripsis, few discussed how TFs influence this process for tumor progression.

METHODS

This study focused on genomic alterations in amplification associated regions within chromosome 17. Inter-/intra-chromosomal rearrangements were analyzed using whole genome sequencing data of breast tumors in the Cancer Genome Atlas (TCGA) cohort. Common ERα binding sites were defined based on MCF-7, T47D, and MDA-MB-134 breast cancer cell lines using univariate K-means clustering methods. Nanopore sequencing technology was applied to validate frequent rearrangements detected between ATC loci on 17q23 and an ERα hub on 20q13. The efficacy of pharmacological inhibition of a potentially druggable target gene on 17q23 was evaluated using breast cancer cell lines and patient-derived circulating breast tumor cells.

RESULTS

There are five adjoining regions from 17q11.1 to 17q24.1 being hotspots of chromothripsis. Inter-/intra-chromosomal rearrangements of these regions occurred more frequently in ERα-positive tumors than in ERα-negative tumors. In addition, the locations of the rearrangements were often mapped within or close to dense ERα binding sites localized on these five 17q regions or other chromosomes. This chromothriptic event was linked to concordant upregulation of 96 loci that predominantly regulate cell-cycle machineries in advanced luminal tumors. Genome-editing analysis confirmed that an ERα hub localized on 20q13 coordinately regulates a subset of these loci localized on 17q23 through long-range chromosome interactions. One of these loci, Tousled Like Kinase 2 (TLK2) known to participate in DNA damage checkpoint control, is an actionable target using phenothiazine antipsychotics (PTZs). The antiproliferative effect of PTZs was prominent in high TLK2-expressing cells, compared to low expressing cells.

CONCLUSION

This study demonstrates a new approach for identifying tumorigenic drivers from genomic regions highly susceptible to ERα-related chromothripsis. We found a group of luminal breast tumors displaying 17q-related chromothripsis for which antipsychotics can be repurposed as treatment adjuncts.

A single ChIP-seq dataset is sufficient for comprehensive analysis of motifs co-occurrence with MCOT package

Recognition of composite elements consisting of two transcription factor binding sites gets behind the studies of tissue-, stage- and condition-specific transcription. Genome-wide data on transcription factor binding generated with ChIP-seq method facilitate an identification of composite elements, but the existing bioinformatics tools either require ChIP-seq datasets for both partner transcription factors, or omit composite elements with motifs overlapping. Here we present an universal Motifs Co-Occurrence Tool (MCOT) that retrieves maximum information about overrepresented composite elements from a single ChIP-seq dataset. This includes homo- and heterotypic composite elements of four mutual orientations of motifs, separated with a spacer or overlapping, even if recognition of motifs within composite element requires various stringencies. Analysis of 52 ChIP-seq datasets for 18 human transcription factors confirmed that for over 60% of analyzed datasets and transcription factors predicted co-occurrence of motifs implied experimentally proven protein-protein interaction of respecting transcription factors. Analysis of 164 ChIP-seq datasets for 57 mammalian transcription factors showed that abundance of predicted composite elements with an overlap of motifs compared to those with a spacer more than doubled; and they had 1.5-fold increase of asymmetrical pairs of motifs with one more conservative 'leading' motif and another one 'guided'.

Small-molecule activating SIRT6 elicits therapeutic effects and synergistically promotes anti-tumor activity of vitamin D3 in colorectal cancer

Colorectal cancer (CRC) is the leading cause of cancer death; however, targets with broad anti-CRC effects are limited. Sirtuin6 (SIRT6) is a conserved nicotinamide adenine dinucleotide (NAD⁺)-dependent deacetylase that is widely pathologically downregulated in CRC, but its pharmacological effect in CRC remains undefined due to the lack of small-molecule SIRT6 activators. We searched for a compound activating SIRT6 and investigated its anti-CRC effect in various models. Methods: We identified an allosteric SIRT6 activator, MDL-811. Its ability to enhance SIRT6 deacetylation at protein and cellular levels was evaluated by Fluor de Lys (FDL) and western blots. We assessed the proliferation of 26 CRC cell lines and patient-derived organoids (PDOs) treated with MDL-811. In vivo efficacy of MDL-811 was evaluated in HCT116 cell line- and patient-derived xenografts as well as a spontaneous CRC model. RNA sequencing and real-time quantitative PCR assays were performed to analyze gene expression changes in MDL-811-treated HCT116 cells. Along with controls in SIRT6-overexpressing HCT116 cells, the SIRT6-mediated histone H3 deacetylation at the Cytochrome P450 family 24 subfamily A member 1 (CYP24A1) gene locus was assessed by chromatin immunoprecipitation (ChIP) in MDL-811-treated HCT116 cells. A combination therapy against CRC based on the downstream gene of SIRT6 activation was evaluated in cells and mouse models. Results: MDL-811 significantly activated SIRT6 histone H3 deacetylation (H3K9Ac, H3K18Ac, and H3K56Ac) in vitro and had broad antiproliferative effects on diverse CRC cell lines and PDOs. More importantly, the in vivo anti-tumor efficacy of MDL-811 was demonstrated across cell line- and patient-derived xenografts and in the APC^min/+ spontaneous CRC model. Mechanically, we identified a new downstream target gene of SIRT6 in CRC, CYP24A1. Based on these findings, a combination drug strategy with MDL-811 to synergistically enhance the anti-CRC effect of vitamin D₃ was validated in vitro and in vivo. Conclusions: Our data provide proof of concept that targeting SIRT6 using a small-molecule activator is an attractive therapeutic strategy for CRC and that MDL-811 could be a promising lead compound for further preclinical and clinical studies of treatments for CRC.

ncHMR detector: a computational framework to systematically reveal non-classical functions of histone modification regulators

Recently, several non-classical functions of histone modification regulators (HMRs), independent of their known histone modification substrates and products, have been reported to be essential for specific cellular processes. However, there is no framework designed for identifying such functions systematically. Here, we develop ncHMR detector, the first computational framework to predict non-classical functions and cofactors of a given HMR, based on ChIP-seq data mining. We apply ncHMR detector in ChIP-seq data-rich cell types and predict non-classical functions of HMRs. Finally, we experimentally reveal that the predicted non-classical function of CBX7 is biologically significant for the maintenance of pluripotency.

Misregulation of ELK1, AP1, and E12 Transcription Factor Networks Is Associated with Melanoma Progression

Melanoma is among the most malignant cutaneous cancers and when metastasized results in dramatically high mortality. Despite advances in high-throughput gene expression profiling in cancer transcriptomic studies, our understanding of mechanisms driving melanoma progression is still limited. We present here an in-depth bioinformatic analysis of the melanoma RNAseq, chromatin immunoprecipitation (ChIP)seq, and single-cell (sc)RNA seq data to understand cancer progression. Specifically, we have performed a consensus network analysis of RNA-seq data from clinically re-grouped melanoma samples to identify gene co-expression networks that are conserved in early (stage 1) and late (stage 4/invasive) stage melanoma. Overlaying the fold-change information on co-expression networks revealed several coordinately up or down-regulated subnetworks that may play a critical role in melanoma progression. Furthermore, by incorporating histone lysine-27 acetylation information and highly expressed genes identified from the single-cell RNA data from melanoma patient samples, we present a comprehensive list of pathways, putative protein-protein interactions (PPIs) and transcription factor (TF) networks that are driving cancer progression. From this analysis, we have identified Elk1, AP1 and E12 TF networks that coordinately change expression in late melanoma when compared to early melanoma, implicating these TFs in melanoma progression. Additionally, the sumoylation-associated interactome is upregulated in invasive melanoma. Together, this bioinformatic analysis potentially implicates a combination of TF networks and PPIs in melanoma progression, which if confirmed in the experimental systems, could be used as targets for drug intervention in melanoma.

epiCOLOC: Integrating Large-Scale and Context-Dependent Epigenomics Features for Comprehensive Colocalization Analysis

High-throughput genome-wide epigenomic assays, such as ChIP-seq, DNase-seq and ATAC-seq, have profiled a huge number of functional elements across numerous human tissues/cell types, which provide an unprecedented opportunity to interpret human genome and disease in context-dependent manner. Colocalization analysis determines whether genomic features are functionally related to a given search and will facilitate identifying the underlying biological functions characterizing intricate relationships with queries for genomic regions. Existing colocalization methods leveraged diverse assumptions and background models to assess the significance of enrichment, however, they only provided limited and predefined sets of epigenomic features. Here, we comprehensively collected and integrated over 44,385 bulk or single-cell epigenomic assays across 53 human tissues/cell types, such as transcription factor binding, histone modification, open chromatin and transcriptional event. By classifying these profiles into hierarchy of tissue/cell type, we developed a web portal, epiCOLOC (http://mulinlab.org/epicoloc or http://mulinlab.tmu.edu.cn/epicoloc), for users to perform context-dependent colocalization analysis in a convenient way.

Single-cell RNA sequencing of immune cells in gastric cancer patients

Cancer immunotherapy has achieved positive clinical responses in the treatment of various cancers, including gastric cancer (GC). In this study, we characterized the heterogeneity of T cells isolated from GC patients at the single-cell level using single-cell RNA sequencing. We identified different immune cell subtypes and their heterogeneous transcription factors and depicted their developmental trajectories. In particular, we focused on exhausted CD8⁺ cells and Tregs and discovered that, as compared to control, the IRF8 transcription factor was downregulated in CD8⁺ tumour-infiltrating lymphocytes (TILs) from GC tissues, and that GC patients with lower IRF8 levels in blood CD8⁺ T cells tended to be a at a more advanced disease stage. These findings provide a theoretical basis for targeted immune therapy in GC.

BedSect: An Integrated Web Server Application to Perform Intersection, Visualization, and Functional Annotation of Genomic Regions From Multiple Datasets

A large number of genomic regions, such as transcription factor binding sites (TFBSs) captured from next generation sequencing (NGS) data analyses or those available from the public resource database ENCODE, are generally overlapped to answer a variety of biological questions. Though several command-line tools are available to perform such an analysis, there is a notable lack of an integrated webserver application with which to identify genomic region intersections, generate publication-ready plots depicting subsets of the overlapped regions, and perform functional annotation. Thus, there is an ardent need for a comprehensive and user-friendly webserver application that allows the users to either upload multiple datasets or select from the integrated Gene Transcription Regulation Database (GTRD). We thus introduce BedSect (http://imgsb.org/bedsect/.), which not only fulfils the above criteria but also performs intersection analysis along with visualization of the intersection regions as an UpSet and correlation plot using the integrated Shiny application. Moreover, analyses, including functional annotation, gene ontology, and biological pathways enrichment for the identified unique and intersected genomic regions, can also be performed using the integrated GREAT tool. To view the genomic regions in the genome browser, the inbuilt hyperlink for UCSC can redirect the user to visualize the results as custom tracks.

CRISPR Screens Identify Essential Cell Growth Mediators in BRAF Inhibitor-resistant Melanoma

BRAF is a serine/threonine kinase that harbors activating mutations in ∼7% of human malignancies and ∼60% of melanomas. Despite initial clinical responses to BRAF inhibitors, patients frequently develop drug resistance. To identify candidate therapeutic targets for BRAF inhibitor resistant melanoma, we conduct CRISPR screens in melanoma cells harboring an activating BRAF mutation that had also acquired resistance to BRAF inhibitors. To investigate the mechanisms and pathways enabling resistance to BRAF inhibitors in melanomas, we integrate expression, ATAC-seq, and CRISPR screen data. We identify the JUN family transcription factors and the ETS family transcription factor ETV5 as key regulators of CDK6, which together enable resistance to BRAF inhibitors in melanoma cells. Our findings reveal genes contributing to resistance to a selective BRAF inhibitor PLX4720, providing new insights into gene regulation in BRAF inhibitor resistant melanoma cells.

Senescence-associated genes and non-coding RNAs function in pancreatic cancer progression

Pancreatic cancer is a major cause of mortality with a poor diagnosis and prognosis that most often occurs in elderly patients. Few studies, however, focus on the interplay of age and pancreatic cancer at the transcriptional level. Here we evaluated the possible roles of age-dependent, differentially expressed genes (DEGs) in pancreatic cancer. These DEGs were used to construct a correlation network and clustered in six gene modules, among which two modules were highly correlated with patients' survival time. Integrating different datasets, including ATAC-Seq and ChIP-Seq, we performed multi-parallel analyses and identified eight age-dependent protein coding genes and two non-coding RNAs as potential candidates. These candidates, together with KLF5, a potent functional transcription factor in pancreatic cancer, are likely to be key elements linking cellular senescence and pancreatic cancer, providing insights on the balance between them, as well as on diagnosis and subsequent prognosis of pancreatic cancer.

TFregulomeR reveals transcription factors' context-specific features and functions

Transcription factors (TFs) are sequence-specific DNA binding proteins, fine-tuning spatiotemporal gene expression. Since genomic occupancy of a TF is highly dynamic, it is crucial to study TF binding sites (TFBSs) in a cell-specific context. To date, thousands of ChIP-seq datasets have portrayed the genomic binding landscapes of numerous TFs in different cell types. Although these datasets can be browsed via several platforms, tools that can operate on that data flow are still lacking. Here, we introduce TFregulomeR (https://github.com/benoukraflab/TFregulomeR), an R-library linked to an up-to-date compendium of cistrome and methylome datasets, implemented with functionalities that facilitate integrative analyses. In particular, TFregulomeR enables the characterization of TF binding partners and cell-specific TFBSs, along with the study of TF’s functions in the context of different partnerships and DNA methylation levels. We demonstrated that TFs’ target gene ontologies can differ notably depending on their partners and, by re-analyzing well characterized TFs, we brought to light that numerous leucine zipper TFBSs derived from ChIP-seq experiments documented in current databases were inadequately characterized, due to the fact that their position weight matrices were assembled using a mixture of homodimer and heterodimer binding sites. Altogether, analyses of context-specific transcription regulation with TFregulomeR foster our understanding of regulatory network-dependent TF functions.

Long noncoding RNA SOX2OT promotes the proliferation of pancreatic cancer by binding to FUS

Pancreatic ductal adenocarcinoma (PDAC) is one of the most malignant tumors has one of the worst prognoses, and the role of long noncoding RNAs (lncRNAs) in the biological and pathological processes of pancreatic cancer, including tumor cell proliferation, is a popular topic in tumor research. Our previous study revealed the correlation between high levels of the lncRNA-SOX2OT (SOX2OT) with poor survival outcomes. Cell Counting Kit-8, EdU, Flow cytometry and Colony formation assays as well as Xenograft growth of PDAC cells in mice were used for the detection of PDAC cells proliferation progression. Fluorescence in situ hybridization, RNA-binding protein pulldown and RNA immunoprecipitation assays were also used to identify the putative mechanisms of SOX2OT participating in the tumor progression. SOX2OT and its potential downstream targets were verified by Western blot and quantitative real-time polymerase chain reaction (qRT-PCR). SOX2OT was confirmed to promote the proliferation of PDAC cells. It was found to directly physically bind to FUS and we also demonstrated that FUS protein stability was affected by binding with SOX2OT and FUS could suppressed PDAC tumor by regulating cell cycle-associated factors CCND1 and p27. Our findings suggest that SOX2OT may act as a tumor promoter in PDAC through physically binding FUS and regulating its downstream cell cycle-associated factors CCND1 and p27. It may serve as an effective target for antitumor treatment for pancreatic cancer.