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

Integrative analysis of the molecular mechanisms, immunological features and immunotherapy response of ferroptosis regulators across 33 cancer types

Ferroptosis is a recently described mode of cell death caused by the accumulation of intracellular iron and lipid reactive oxygen species (ROS), which play critical roles in tumorigenesis and cancer progression. However, the underlying molecular mechanisms and promising biomarkers of ferroptosis among cancers remain to be elucidated. In this study, 30 ferroptosis regulators in ferroptosis-related signaling pathways were identified and analyzed in 33 cancer types. We found transcriptomic aberrations and evaluated the prognostic value of ferroptosis regulators across 33 cancer types. Then, we predicted and validated potential transcription factors (including E2F7, KLF5 and FOXM1) and therapeutic drugs (such as cyclophosphamide, vinblastine, and gefitinib) that target ferroptosis regulators in cancer. Moreover, we explored the molecular mechanisms of ferroptosis and found that signaling pathways such as the IL-1 and IL-2 pathways are closely associated with ferroptosis. Additionally, we found that ferroptosis regulators have a close relationship with immunity-related parameters, including the immune score, immune cell infiltration level, and immune checkpoint protein level. Finally, we determined a ferroptosis score using GSVA method. We found that the ferroptosis score effectively predicted ferroptotic cell death in tumor samples. And ferroptosis score is served as an independent prognostic indicator for the incidence and recurrence of cancers. More importantly, patients with high ferroptosis scores received greater benefit from immunotherapy. We aslo created an online webserver based on the nomogram prognostic model to predict the survival in immunotherapy cohort. The reason for this outcome is partially the result of patients with a high ferroptosis rate also having high immune scores, HLA-related gene expression and immune checkpoint protein expression, such as PDL2 and TIM3. Moreover, patients with high ferroptosis scores exhibited CD8 T cell and TIL infiltration and immune-related signaling pathway enrichment. In summary, we systematically summarize the molecular characteristics, clinical relevance and immune features of ferroptosis across cancers and show that the ferroptosis score can be used as a prognostic factor and for the evaluation of immunotherapy effects.

Clinicopathological significance and underlying molecular mechanism of downregulation of basonuclin 1 expression in ovarian carcinoma

In this study, we aim to identify the clinical significance of basonuclin 1 (BNC1) expression in ovarian carcinoma (OV) and to explore its latent mechanisms. Via integrating in-house tissue microarrays, gene chips, and RNA-sequencing data, we explored the expression and clinical value of BNC1 in OV. Immunohistochemical staining was utilized to confirm the protein expression status of BNC1. A combined SMD of -2.339 (95% CI: -3.649 to -1.028, P < 0.001) identified that BNC1 was downregulated based on 1346 samples, and the sROC (AUC = 0.93) showed a favorable discriminatory ability of BNC1 in OV patients. We used univariate and multivariate Cox regulation to evaluate the prognostic role of BNC1 for OV patients, and a combined hazard ratio of 0.717 (95% CI: 0.445-0.989, P < 0.001) revealed that BNC1 was a protective factor for OV. Furthermore, the fraction of infiltrating naive B cells, memory B cells, and other immune cells showed statistical differences between the high- and low-BNC1 expression groups through cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) algorithm. Enrichment analysis showed that BNC1 may have a relationship with immune-related items in OV. By predicting the potential regulatory transcription factors (TFs) of BNC1, friend leukemia virus integration 1 (FLI1) may be a potential upstream TF of BNC1. Corporately, a decreasing trend of BNC1 may serve as a tumor suppressor and prognostic biomarker in OV patients. Moreover, BNC1 may take part in immune-related pathways and influence the fraction of tumor-infiltrating immune cells.

CTCF Expression and Dynamic Motif Accessibility Modulates Epithelial-Mesenchymal Gene Expression

Epithelial-mesenchymal transition (EMT) and its reversal, mesenchymal-epithelial transition (MET) drive tissue reorganization critical for early development. In carcinomas, processing through EMT, MET, or partial states promotes migration, invasion, dormancy, and metastatic colonization. As a reversible process, EMT is inherently regulated at epigenetic and epigenomic levels. To understand the epigenomic nature of reversible EMT and its partial states, we characterized chromatin accessibility dynamics, transcriptomic output, protein expression, and cellular phenotypes during stepwise reversible EMT. We find that the chromatin insulating protein machinery, including CTCF, is suppressed and re-expressed, coincident with broad alterations in chromatin accessibility, during EMT/MET, and is lower in triple-negative breast cancer cell lines with EMT features. Through an analysis of chromatin accessibility using ATAC-seq, we identify that early phases of EMT are characterized by enrichment for AP-1 family member binding motifs, but also by a diminished enrichment for CTCF binding motifs. Through a loss-of-function analysis, we demonstrate that the suppression of CTCF alters cellular plasticity, strengthening the epithelial phenotype via the upregulation of epithelial markers E-cadherin/CDH1 and downregulation of N-cadherin/CDH2. Conversely, the upregulation of CTCF leads to the upregulation of EMT gene expression and an increase in mesenchymal traits. These findings are indicative of a role of CTCF in regulating epithelial-mesenchymal plasticity and gene expression.

Prognostic biomarker SMARCC1 and its association with immune infiltrates in hepatocellular carcinoma

Background

Epigenetic alterations contribute greatly to metastasis and dissemination in hepatocellular carcinoma (HCC). SMARCC1, as a SWI/SNF chromatin remodeling factor, has been reported to play important roles in many cancers. For the first time, with the bioinformatics analysis and wet-bench experiments, we explored the biological significance of SMARCC1 and its potential as putative therapeutic target in HCC.

Methods

The mRNA expression profiles and prognostic value of SMARCC1 were analyzed in the Oncomine, UALCAN and Kaplan–Meier Plotter databases. The expression of SMARCC1 and associated clinicopathological factors were further evaluated using a tissue microarray. Differentially expressed genes associated with SMARCC1 in HCC were obtained and analyzed via the LinkedOmics and GEPIA databases and Cytoscape software. To verify the important role of SMARCC1 in HCC, we knocked down and overexpressed SMARCC1 in different hepatic cell lines and conducted several functional experiments. Then, we evaluated the mutation profiles and transcriptional regulators of SMARCC1 using the cBioPortal, COSMIC, CistromeDB and TCGA databases. Finally, we addressed the relationship of SMARCC1 expression with immune cell infiltration via TIMER database analysis.

Results

Through data mining and tissue microarray verification, we found that the protein and mRNA levels of SMARCC1 are high in tumor tissues, which has remarkable diagnostic value in HCC patients. SMARCC1 and its hub genes showed prognostic value in HCC. Furthermore, we confirmed that SMARCC1 influenced the proliferation, migration, and invasion of HCC cells. Moreover, correlation analyses revealed that SMARCC1 expression was positively correlated with ZBTB40 transcription factors and negatively correlated with the DNA methylation level. Overall, we found that SMARCC1 affects immune infiltration and plays a tumor-promoting role in HCC.

Conclusions

SMARCC1 is overexpressed and is a putative prognostic predictor in HCC. Due to the tumor-promoting role of SMARCC1, treatments inhibiting DNA methyltransferases and transcription factors or weakening the role of SMARCC1 in immune infiltration might improve the survival of HCC patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02413-w.

MCM2-7 in Clear Cell Renal Cell Carcinoma: MCM7 Promotes Tumor Cell Proliferation

Background

Clear cell renal cell carcinoma (ccRCC) accounts for 60-70% of renal cell carcinoma (RCC) cases. Finding more therapeutic targets for advanced ccRCC is an urgent mission. The minichromosome maintenance proteins 2-7 (MCM2-7) protein forms a stable heterohexamer and plays an important role in DNA replication in eukaryotic cells. In the study, we provide a comprehensive study of MCM2-7 genes expression and their potential roles in ccRCC.

Methods

The expression and prognosis of the MCM2-7 genes in ccRCC were analyzed using data from TCGA, GEO and ArrayExpress. MCM2-7 related genes were identified by weighted co-expression network analysis (WGCNA) and Metascape. CancerSEA and GSEA were used to analyze the function of MCM2–7 genes in ccRCC. The gene effect scores (CERES) of MCM2-7, which reflects carcinogenic or tumor suppressor, were obtained from DepMap. We used clinical and expression data of MCM2-7 from the TCGA dataset and the LASSO Cox regression analysis to develop a risk score to predict survival of patients with ccRCC. The correlations between risk score and other clinical indicators such as gender, age and stage were also analyzed. Further validation of this risk score was engaged in another cohort, E-MTAB-1980 from the ArrayExpress dataset.

Results

The mRNA and protein expression of MCM2-7 were increased in ccRCC compared with normal tissues. High MCM2, MCM4, MCM6 and MCM7 expression were associated with a poor prognosis of ccRCC patients. Functional enrichment analysis revealed that MCM2-7 might influence the progress of ccRCC by regulating the cell cycle. Knockdown of MCM7 can inhibit the proliferation of ccRCC cells. A two-gene risk score including MCM4 and MCM6 can predict overall survival (OS) of ccRCC patients. The risk score was successfully verified by further using Arrayexpress cohort.

Conclusion

We analyze MCM2-7 mRNA and protein levels in ccRCC. MCM7 is determined to promote tumor proliferation. Meanwhile, our study has determined a risk score model composed of MCM2-7 can predict the prognosis of ccRCC patients, which may help future treatment strategies.

Interrogating cell type-specific cooperation of transcriptional regulators in 3D chromatin

Context-specific activities of transcription regulators (TRs) in the nucleus modulate spatiotemporal gene expression precisely. Using the largest ChIP-seq data and chromatin loops in the human K562 cell line, we initially interrogated TR cooperation in 3D chromatin via a graphical model and revealed many known and novel TRs manipulating context-specific pathways. To explore TR cooperation across broad tissue/cell types, we systematically leveraged large-scale open chromatin profiles, computational footprinting, and high-resolution chromatin interactions to investigate tissue/cell type-specific TR cooperation. We first delineated a landscape of TR cooperation across 40 human tissue/cell types. Network modularity analyses uncovered the commonality and specificity of TR cooperation in different conditions. We also demonstrated that TR cooperation information can better interpret the disease-causal variants identified by genome-wide association studies and recapitulate cell states during neural development. Our study characterizes shared and unique patterns of TR cooperation associated with the cell type specificity of gene regulation in 3D chromatin.

•

Computational inference of transcriptional regulator (TR) cooperation in 3D chromatin

•

A landscape of 3D TR cooperation across 40 human tissue/cell types

•

TR cooperation can better interpret the disease-causal variants identified by GWAS

•

Cooperation of certain TRs shapes context-specific gene regulation in cell development

Epigenetic loss of heterogeneity from low to high grade localized prostate tumours

Identifying precise molecular subtypes attributable to specific stages of localized prostate cancer has proven difficult due to high levels of heterogeneity. Bulk assays represent a population-average, which mask the heterogeneity that exists at the single-cell level. In this work, we sequence the accessible chromatin regions of 14,424 single-cells from 18 flash-frozen prostate tumours. We observe shared chromatin features among low-grade prostate cancer cells are lost in high-grade tumours. Despite this loss, high-grade tumours exhibit an enrichment for FOXA1, HOXB13 and CDX2 transcription factor binding sites, indicating a shared trans-regulatory programme. We identify two unique genes encoding neuronal adhesion molecules that are highly accessible in high-grade prostate tumours. We show NRXN1 and NLGN1 expression in epithelial, endothelial, immune and neuronal cells in prostate cancer using cyclic immunofluorescence. Our results provide a deeper understanding of the active gene regulatory networks in primary prostate tumours, critical for molecular stratification of the disease.

Reduced NCOR2 expression accelerates androgen deprivation therapy failure in prostate cancer

This study addresses the roles of nuclear receptor corepressor 2 (NCOR2) in prostate cancer (PC) progression in response to androgen deprivation therapy (ADT). Reduced NCOR2 expression significantly associates with shorter disease-free survival in patients with PC receiving adjuvant ADT. Utilizing the CWR22 xenograft model, we demonstrate that stably reduced NCOR2 expression accelerates disease recurrence following ADT, associates with gene expression patterns that include neuroendocrine features, and induces DNA hypermethylation. Stably reduced NCOR2 expression in isogenic LNCaP (androgen-sensitive) and LNCaP-C4–2 (androgen-independent) cells revealed that NCOR2 reduction phenocopies the impact of androgen treatment and induces global DNA hypermethylation patterns. NCOR2 genomic binding is greatest in LNCaP-C4–2 cells and most clearly associates with forkhead box (FOX) transcription factor FOXA1 binding. NCOR2 binding significantly associates with transcriptional regulation most when in active enhancer regions. These studies reveal robust roles for NCOR2 in regulating the PC transcriptome and epigenome and underscore recent mutational studies linking NCOR2 loss of function to PC disease progression.

Long et al. show that reduced levels of NCOR2 lead to accelerated prostate cancer recurrence during androgen withdrawal in a patient-derived xenograft model. NCOR2 reduction is characterized by incomplete response to androgen withdrawal, and recurrent tumors show increased neuroendocrine traits. These phenotypic changes are associated with hypermethylated enhancers.

Insights into high-risk multiple myeloma from an analysis of the role of PHF19 in cancer

Despite  improvements in outcome, 15-25% of newly diagnosed multiple myeloma (MM) patients have treatment resistant high-risk (HR) disease with a poor survival. The lack of a genetic basis for HR has focused attention on the role played by epigenetic changes. Aberrant expression and somatic mutations affecting genes involved in the regulation of tri-methylation of the lysine (K) 27 on histone 3 H3 (H3K27me3) are common in cancer. H3K27me3 is catalyzed by EZH2, the catalytic subunit of the Polycomb Repressive Complex 2 (PRC2). The deregulation of H3K27me3 has been shown to be involved in oncogenic transformation and tumor progression in a variety of hematological malignancies including MM. Recently we have shown that aberrant overexpression of the PRC2 subunit PHD Finger Protein 19 (PHF19) is the most significant overall contributor to HR status further focusing attention on the role played by epigenetic change in MM. By modulating both the PRC2/EZH2 catalytic activity and recruitment, PHF19 regulates the expression of key genes involved in cell growth and differentiation. Here we review the expression, regulation and function of PHF19 both in normal and the pathological contexts of solid cancers and MM. We present evidence that strongly implicates PHF19 in the regulation of genes important in cell cycle and the genetic stability of MM cells making it highly relevant to HR MM behavior. A detailed understanding of the normal and pathological functions of PHF19 will allow us to design therapeutic strategies able to target aggressive subsets of MM.

Master lineage transcription factors anchor trans mega transcriptional complexes at highly accessible enhancer sites to promote long-range chromatin clustering and transcription of distal target genes

The term 'super enhancers' (SE) has been widely used to describe stretches of closely localized enhancers that are occupied collectively by large numbers of transcription factors (TFs) and co-factors, and control the transcription of highly-expressed genes. Through integrated analysis of >600 DNase-seq, ChIP-seq, GRO-seq, STARR-seq, RNA-seq, Hi-C and ChIA-PET data in five human cancer cell lines, we identified a new class of autonomous SEs (aSEs) that are excluded from classic SE calls by the widely used Rank Ordering of Super-Enhancers (ROSE) method. TF footprint analysis revealed that compared to classic SEs and regular enhancers, aSEs are tightly bound by a dense array of master lineage TFs, which serve as anchors to recruit additional TFs and co-factors in trans. In addition, aSEs are preferentially enriched for Cohesins, which likely involve in stabilizing long-distance interactions between aSEs and their distal target genes. Finally, we showed that aSEs can be reliably predicted using a single DNase-seq data or combined with Mediator and/or P300 ChIP-seq. Overall, our study demonstrates that aSEs represent a unique class of functionally important enhancer elements that distally regulate the transcription of highly expressed genes.

Identification of Key Immune-Related Genes in the Progression of Septic Shock

Objective: Septic shock is the severe complication of sepsis, with a high mortality. The inflammatory response regulates the immune status and mediates the progression of septic shock. In this study, we aim to identify the key immune-related genes (IRGs) of septic shock and explore their potential mechanism. Methods: Gene expression profiles of septic shock blood samples and normal whole blood samples were retrieved from the Gene Expression Omnibus (GEO) and Genotype-Tissue Expression Portal (GTEx). The differential expression genes (DEGs) and septic shock-specific immune-related genes (SSSIRGs) were evaluated and identified, along with the immune components by "cell type identification by estimating relative subsets of RNA transcripts (CIBERSORT, version x)" algorithm. Additionally, in order to explore the key regulatory network, the relationship among SSSIRGs, upstream transcription factors (TFs), and downstream signaling pathways were also identified by Gene Set Variation Analysis (GSVA) and co-expression analysis. Moreover, the Connectivity Map (CMap) analysis was applied to find bioactive small molecules against the members of regulation network while Chromatin Immunoprecipitation sequencing (ChIP-seq) and Assay for Targeting Accessible-Chromatin with high-throughput sequencing (ATAC-seq) data were used to validate the regulation mechanism of the network. Results: A total of 14,843 DEGs were found between 63 septic shock blood samples and 337 normal whole blood samples. Then, we identified septic shock-specific 839 IRGs as the intersection of DEGs and IRGs. Moreover, we uncovered the regulatory networks based on co-expression analysis and found 28 co-expression interaction pairs. In the regulation network, protein phosphatase 3, catalytic subunit, alpha isozyme (PPP3CA) may regulate late estrogen response, glycolysis and TNFα signaling via NFκB and HLA; Kirsten rat sarcoma viral oncogene homolog (KRAS) may be related to late estrogen response and HLA; and Toll-like receptor 8 (TLR8) may be associated with TNFα signaling via NFκB. And the regulation mechanisms between TFs and IRGs (TLR8, PPP3CA, and KRAS) were validated by ChIP-seq and ATAC-seq. Conclusion: Our data identify three SSSIRGs (TLR8, PPP3CA, and KRAS) as candidate therapeutic targets for septic shock and provide constructed regulatory networks in septic shock to explore its potential mechanism.

Genetic fusions favor tumorigenesis through degron loss in oncogenes

Chromosomal rearrangements can generate genetic fusions composed of two distinct gene sequences, many of which have been implicated in tumorigenesis and progression. Our study proposes a model whereby oncogenic gene fusions frequently alter the protein stability of the resulting fusion products, via exchanging protein degradation signal (degron) between gene sequences. Computational analyses of The Cancer Genome Atlas (TCGA) identify 2,406 cases of degron exchange events and reveal an enrichment of oncogene stabilization due to loss of degrons from fusion. Furthermore, we identify and experimentally validate that some recurrent fusions, such as BCR-ABL, CCDC6-RET and PML-RARA fusions, perturb protein stability by exchanging internal degrons. Likewise, we also validate that EGFR or RAF1 fusions can be stabilized by losing a computationally-predicted C-terminal degron. Thus, complementary to enhanced oncogene transcription via promoter swapping, our model of degron loss illustrates another general mechanism for recurrent fusion proteins in driving tumorigenesis.

Enhancer architecture-dependent multilayered transcriptional regulation orchestrates RA signaling-induced early lineage differentiation of ESCs

Signaling pathway-driven target gene transcription is critical for fate determination of embryonic stem cells (ESCs), but enhancer-dependent transcriptional regulation in these processes remains poorly understood. Here, we report enhancer architecture-dependent multilayered transcriptional regulation at the Halr1–Hoxa1 locus that orchestrates retinoic acid (RA) signaling-induced early lineage differentiation of ESCs. We show that both homeobox A1 (Hoxa1) and Hoxa adjacent long non-coding RNA 1 (Halr1) are identified as direct downstream targets of RA signaling and regulated by RARA/RXRA via RA response elements (RAREs). Chromosome conformation capture-based screens indicate that RA signaling promotes enhancer interactions essential for Hoxa1 and Halr1 expression and mesendoderm differentiation of ESCs. Furthermore, the results also show that HOXA1 promotes expression of Halr1 through binding to enhancer; conversely, loss of Halr1 enhances interaction between Hoxa1 chromatin and four distal enhancers but weakens interaction with chromatin inside the HoxA cluster, leading to RA signaling-induced Hoxa1 overactivation and enhanced endoderm differentiation. These findings reveal complex transcriptional regulation involving synergistic regulation by enhancers, transcription factors and lncRNA. This work provides new insight into intrinsic molecular mechanisms underlying ESC fate determination during RA signaling-induced early differentiation.

Ginsenoside Rh2 reduces m6A RNA methylation in cancer via the KIF26B-SRF positive feedback loop

Background

The underlying mechanisms of the potential tumor-suppressive effects of ginsenoside Rh2 are complex. N6-methyladenosine (m6A) RNA methylation is usually dysregulated in cancer. This study explored the regulatory effect of ginsenoside Rh2 on m6A RNA methylation in cancer.

Methods: m6A RNA quantification and gene-specific m6A RIP-qPCR assays were applied to assess total and gene-specific m6A RNA levels. Co-immunoprecipitation, fractionation western blotting, and immunofluorescence staining were performed to detect protein interactions and distribution. QRT-PCR, dual-luciferase, and ChIP-qPCR assays were conducted to check the transcriptional regulation.

Results

Ginsenoside Rh2 reduces m6A RNA methylation and KIF26B expression in a dose-dependent manner in some cancers. KIF26B interacts with ZC3H13 and CBLL1 in the cytoplasm of cancer cells and enhances their nuclear distribution. KIF26B inhibition reduces m6A RNA methylation level in cancer cells. SRF bound to the KIF26B promoter and activated its transcription. SRF mRNA m6A abundance significantly decreased upon KIF26B silencing. SRF knockdown suppressed cancer cell proliferation and growth both in vitro and in vivo, the effect of which was partly rescued by KIF26B overexpression.

Conclusion: ginsenoside Rh2 reduces m6A RNA methylation via downregulating KIF26B expression in some cancer cells. KIF26B elevates m6A RNA methylation via enhancing ZC3H13/CBLL1 nuclear localization. KIF26B-SRF forms a positive feedback loop facilitating tumor growth.

ReMap 2022: a database of Human, Mouse, Drosophila and Arabidopsis regulatory regions from an integrative analysis of DNA-binding sequencing experiments

ReMap (https://remap.univ-amu.fr) aims to provide manually curated, high-quality catalogs of regulatory regions resulting from a large-scale integrative analysis of DNA-binding experiments in Human, Mouse, Fly and Arabidopsis thaliana for hundreds of transcription factors and regulators. In this 2022 update, we have uniformly processed >11 000 DNA-binding sequencing datasets from public sources across four species. The updated Human regulatory atlas includes 8103 datasets covering a total of 1210 transcriptional regulators (TRs) with a catalog of 182 million (M) peaks, while the updated Arabidopsis atlas reaches 4.8M peaks, 423 TRs across 694 datasets. Also, this ReMap release is enriched by two new regulatory catalogs for Mus musculus and Drosophila melanogaster. First, the Mouse regulatory catalog consists of 123M peaks across 648 TRs as a result of the integration and validation of 5503 ChIP-seq datasets. Second, the Drosophila melanogaster catalog contains 16.6M peaks across 550 TRs from the integration of 1205 datasets. The four regulatory catalogs are browsable through track hubs at UCSC, Ensembl and NCBI genome browsers. Finally, ReMap 2022 comes with a new Cis Regulatory Module identification method, improved quality controls, faster search results, and better user experience with an interactive tour and video tutorials on browsing and filtering ReMap catalogs.

Multi-omics mapping of human papillomavirus integration sites illuminates novel cervical cancer target genes

BACKGROUND

Integration of human papillomavirus (HPV) into the host genome is a dominant feature of invasive cervical cancer (ICC), yet the tumorigenicity of cis genomic changes at integration sites remains largely understudied.

METHODS

Combining multi-omics data from The Cancer Genome Atlas with patient-matched long-read sequencing of HPV integration sites, we developed a strategy for using HPV integration events to identify and prioritise novel candidate ICC target genes (integration-detected genes (IDGs)). Four IDGs were then chosen for in vitro functional studies employing small interfering RNA-mediated knockdown in cell migration, proliferation and colony formation assays.

RESULTS

PacBio data revealed 267 unique human-HPV breakpoints comprising 87 total integration events in eight tumours. Candidate IDGs were filtered based on the following criteria: (1) proximity to integration site, (2) clonal representation of integration event, (3) tumour-specific expression (Z-score) and (4) association with ICC survival. Four candidates prioritised based on their unknown function in ICC (BNC1, RSBN1, USP36 and TAOK3) exhibited oncogenic properties in cervical cancer cell lines. Further, annotation of integration events provided clues regarding potential mechanisms underlying altered IDG expression in both integrated and non-integrated ICC tumours.

CONCLUSIONS

HPV integration events can guide the identification of novel IDGs for further study in cervical carcinogenesis and as putative therapeutic targets.

An emergent Wnt5a/YAP/TAZ regulatory circuit and its possible role in cancer

Wnt5a is a ligand that plays several roles in development, homeostasis, and disease. A growing body of evidence indicates that Wnt5a is involved in cancer progression. Despite extensive research in this field, our knowledge about how Wnt5a is precisely involved in cancer is still incomplete. It is usually thought that certain combinations of Frizzled receptors and co-receptors might explain the observed effects of Wnt5a either as a tumor suppressor or by promoting migration and invasion. While accepting this 'receptor context' model, this review proposes that Wnt5a is integrated within a larger regulatory circuit involving β-catenin, YAP/TAZ, and LATS1/2. Remarkably, WNT5A and YAP1 are transcriptionally regulated by the Hippo and Wnt pathways, respectively, and might form a regulatory circuit acting through LATS kinases and secreted Wnt/β-catenin inhibitors, including Wnt5a itself. Therefore, understanding the precise role of Wnt5a and YAP in cancer requires a systems biology perspective.

Comprehensive Analysis of microRNAs in Human Adult Erythropoiesis

MicroRNAs (miRNAs) are small non-coding RNAs, which play an important role in various cellular and developmental processes. The study of miRNAs in erythropoiesis is crucial to uncover the cellular pathways that are modulated during the different stages of erythroid differentiation. Using erythroid cells derived from human CD34+ hematopoietic stem and progenitor cells (HSPCs)and small RNA sequencing, our study unravels the various miRNAs involved in critical cellular pathways in erythroid maturation. We analyzed the occupancy of erythroid transcription factors and chromatin accessibility in the promoter and enhancer regions of the differentially expressed miRNAs to integrate miRNAs in the transcriptional circuitry of erythropoiesis. Analysis of the targets of the differentially expressed miRNAs revealed novel pathways in erythroid differentiation. Finally, we described the application of Clustered regularly interspaced short palindromic repeats-Cas9 (CRISPR-Cas9) based editing of miRNAs to study their function in human erythropoiesis.

Comprehensive understanding of Tn5 insertion preference improves transcription regulatory element identification

Tn5 transposase, which can efficiently tagment the genome, has been widely adopted as a molecular tool in next-generation sequencing, from short-read sequencing to more complex methods such as assay for transposase-accessible chromatin using sequencing (ATAC-seq). Here, we systematically map Tn5 insertion characteristics across several model organisms, finding critical parameters that affect its insertion. On naked genomic DNA, we found that Tn5 insertion is not uniformly distributed or random. To uncover drivers of these biases, we used a machine learning framework, which revealed that DNA shape cooperatively works with DNA motif to affect Tn5 insertion preference. These intrinsic insertion preferences can be modeled using nucleotide dependence information from DNA sequences, and we developed a computational pipeline to correct for these biases in ATAC-seq data. Using our pipeline, we show that bias correction improves the overall performance of ATAC-seq peak detection, recovering many potential false-negative peaks. Furthermore, we found that these peaks are bound by transcription factors, underscoring the biological relevance of capturing this additional information. These findings highlight the benefits of an improved understanding and precise correction of Tn5 insertion preference.

Opposing transcriptional programs of KLF5 and AR emerge during therapy for advanced prostate cancer

Endocrine therapies for prostate cancer inhibit the androgen receptor (AR) transcription factor. In most cases, AR activity resumes during therapy and drives progression to castration-resistant prostate cancer (CRPC). However, therapy can also promote lineage plasticity and select for AR-independent phenotypes that are uniformly lethal. Here, we demonstrate the stem cell transcription factor Krüppel-like factor 5 (KLF5) is low or absent in prostate cancers prior to endocrine therapy, but induced in a subset of CRPC, including CRPC displaying lineage plasticity. KLF5 and AR physically interact on chromatin and drive opposing transcriptional programs, with KLF5 promoting cellular migration, anchorage-independent growth, and basal epithelial cell phenotypes. We identify ERBB2 as a point of transcriptional convergence displaying activation by KLF5 and repression by AR. ERBB2 inhibitors preferentially block KLF5-driven oncogenic phenotypes. These findings implicate KLF5 as an oncogene that can be upregulated in CRPC to oppose AR activities and promote lineage plasticity.

Introducing ADNP and SIRT1 as new partners regulating microtubules and histone methylation

Activity-dependent neuroprotective protein (ADNP) is essential for brain formation and function. As such, de novo mutations in ADNP lead to the autistic ADNP syndrome and somatic ADNP mutations may drive Alzheimer's disease (AD) tauopathy. Sirtuin 1 (SIRT1) is positively associated with aging, the major risk for AD. Here, we revealed two key interaction sites for ADNP and SIRT1. One, at the microtubule end-binding protein (EB1 and EB3) Tau level, with EB1/EB3 serving as amplifiers for microtubule dynamics, synapse formation, axonal transport, and protection against tauopathy. Two, on the DNA/chromatin site, with yin yang 1, histone deacetylase 2, and ADNP, sharing a DNA binding motif and regulating SIRT1, ADNP, and EB1 (MAPRE1). This interaction was linked to sex- and age-dependent altered histone modification, associated with ADNP/SIRT1/WD repeat-containing protein 5, which mediates the assembly of histone modification complexes. Single-cell RNA and protein expression analyses as well as gene expression correlations placed SIRT1-ADNP and either MAPRE1 (EB1), MAPRE3 (EB3), or both in the same mouse and human cell; however, while MAPRE1 seemed to be similarly regulated to ADNP and SIRT1, MAPRE3 seemed to deviate. Finally, we demonstrated an extremely tight correlation for the gene transcripts described above, including related gene products. This correlation was specifically abolished in affected postmortem AD and Parkinson's disease brain select areas compared to matched controls, while being maintained in blood samples. Thus, we identified an ADNP-SIRT1 complex that may serve as a new target for the understanding of brain degeneration.

FGFR-inhibitor-mediated dismissal of SWI/SNF complexes from YAP-dependent enhancers induces adaptive therapeutic resistance

How cancer cells adapt to evade the therapeutic effects of drugs targeting oncogenic drivers is poorly understood. Here we report an epigenetic mechanism leading to the adaptive resistance of triple-negative breast cancer (TNBC) to fibroblast growth factor receptor (FGFR) inhibitors. Prolonged FGFR inhibition suppresses the function of BRG1-dependent chromatin remodelling, leading to an epigenetic state that derepresses YAP-associated enhancers. These chromatin changes induce the expression of several amino acid transporters, resulting in increased intracellular levels of specific amino acids that reactivate mTORC1. Consistent with this mechanism, addition of mTORC1 or YAP inhibitors to FGFR blockade synergistically attenuated the growth of TNBC patient-derived xenograft models. Collectively, these findings reveal a feedback loop involving an epigenetic state transition and metabolic reprogramming that leads to adaptive therapeutic resistance and provides potential therapeutic strategies to overcome this mechanism of resistance.

Identification of Key eRNAs for Spinal Cord Injury by Integrated Multinomial Bioinformatics Analysis

Background: Spinal cord injury (SCI) is a severe neurological deficit affecting both young and older people worldwide. The potential role of key enhancer RNAs (eRNAs) in SCI remains elusive, which is a prominent challenge in the trauma repair process. This study aims to investigate the roles of key eRNAs, transcription factors (TFs), signaling pathways, and small-molecule inhibitors in SCI using multi-omics bioinformatics analysis. Methods: Microarray data of peripheral blood mononuclear cell (PBMC) samples from 27 healthy volunteers and 25 chronic-phase SCI patients were retrieved from the Gene Expression Omnibus database. Differentially expressed transcription factors (DETFs), differentially expressed enhancer RNAs (DEeRNAs), and differentially expressed target genes (DETGs) were identified using the Linear Models for Microarray Data (limma) package. Fraction of immune cells was estimated using CIBERSORT algorithm. Gene Set Variation Analysis (GSVA) was applied to identify the downstream signaling pathways. The eRNA regulatory network was constructed based on the correlation results. Connectivity Map (CMap) database was used to find potential drugs for SCI patients. The cellular communication analysis was performed to explore the molecular regulation mechanism of SCI based on single-cell RNA sequencing (scRNA-seq) data. Chromatin immunoprecipitation sequencing (ChIP-seq) and Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) data were used to validate the key regulatory mechanisms. scRNA-seq dataset was used to validate the cell subtype localization of the key eRNAs. Results: In total, 21 DETFs, 24 DEeRNAs, and 829 DETGs were identified. A regulatory network of 13 DETFs, six DEeRNAs, seven DETGs, two hallmark pathways, two immune cells, and six immune pathways was constructed. The link of Splicing factor proline and glutamine rich (SFPQ) (TF) and vesicular overexpressed in cancer prosurvival protein 1 (VOPP1) (eRNA) (R = 0.990, p < 0.001, positive), VOPP1 (eRNA) and epidermal growth factor receptor (EGFR) (target gene) (R = 0.974, p < 0.001, positive), VOPP1, and T helper (Th) cells (R = -0.987, p < 0.001, negative), and VOPP1 and hallmark coagulation (R = 0.937, p < 0.001, positive) was selected. Trichostatin A was considered the best compound target to SCI-related eRNAs (specificity = 0.471, p < 0.001). Conclusion: VOPP1, upregulated by SFPQ, strengthened the transient expression of EGFR. Th cells and coagulation were the potential downstream pathways of VOPP1. This regulatory network and potential inhibitors provide novel diagnostic biomarkers and therapeutic targets for SCI.

Screening and cellular validation of prognostic genes regulated by super enhancers in oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) is the leading cause of death in patients with head and neck cancer. Reliable biomarkers to guide treatment decisions for OSCC remain scarce. The purpose of this study was to identify novel prognostic markers regulated by super enhancers in OSCC. Eight modules were obtained by weighted gene co-expression network analysis (WGCNA), among which MEblue module had the highest correlation with tumor stage, alcohol consumption and smoking. There were 41 genes regulated by super enhancers in MEblue module. Functional analysis showed that 41 super enhancer-regulated genes were involved in cancer progression. A total of twenty transcription factors of the 41 genes were predicted. Prognostic analysis of the 41 genes and the top 5 transcription factors showed that patients with high expression of AHCY, KCMF1, MANBAL and TFDP1 had a poor prognosis. Immunohistochemical analysis showed that AHCY, KCMF1 and MANBAL were highly expressed in OSCC tissue. Cellular experiment demonstrated that TFDP1 promoted AHCY, KCMF1 and MANBAL expression by binding to the super enhancers of these genes. Knockdown of TFDP1, AHCY, KCMF1 and MANBAL inhibited the proliferation of OSCC cells. In conclusion, AHCY, KCMF1 and MANBAL were recognized as super enhancer-regulated prognostic biomarkers regulated by TFDP1 in OSCC.

Laryngeal Squamous Cell Carcinoma: Clinical Significance and Potential Mechanism of Cell Division Cycle 45

Background: Cell division cycle 45 (CDC45) plays an important role in the occurrence and development of numerous carcinomas, but its effect in laryngeal squamous cell carcinoma (LSCC) remains unclear. Materials and Methods: The messenger RNA and protein expression levels of CDC45 in LSCC were evaluated with a t test and the standard mean difference (SMD). The ability of CDC45 expression to distinguish the LSCC was assessed through receiver operating characteristic (ROC) curves. Gene set enrichment analysis (GSEA), protein-protein interaction, public databases, and online tools were used to explore the potential molecular mechanism of CDC45 in LSCC. Results: A high expression of CDC45 was identified in LSCC (SMD = 2.61, 95% confidence interval [1.62-3.61]). Through ROC curves, the expression of CDC45 makes it feasible to distinguish the LSCC group from the non-LSCC counterpart. CDC45 was relevant to the progression-free interval of LSCC patients (log-rank p = 0.03). GSEAs show that CDC45 is related to the cell cycle. CDC45, CDC6, KIF2C, and AURKB were identified as hub genes of LSCC. E2F1 may be the regulatory transcription factor of CDC45. Conclusions: High expression of CDC45 likely demonstrates carcinogenic effects in LSCC, and CDC45 is a potential target in screening and treatment of LSCC.

In vivo CRISPR screens identify the E3 ligase Cop1 as a modulator of macrophage infiltration and cancer immunotherapy target

Despite remarkable clinical efficacy of immune checkpoint blockade (ICB) in cancer treatment, ICB benefits for triple-negative breast cancer (TNBC) remain limited. Through pooled in vivo CRISPR knockout (KO) screens in syngeneic TNBC mouse models, we found that deletion of the E3 ubiquitin ligase Cop1 in cancer cells decreases secretion of macrophage-associated chemokines, reduces tumor macrophage infiltration, enhances anti-tumor immunity, and strengthens ICB response. Transcriptomics, epigenomics, and proteomics analyses revealed that Cop1 functions through proteasomal degradation of the C/ebpδ protein. The Cop1 substrate Trib2 functions as a scaffold linking Cop1 and C/ebpδ, which leads to polyubiquitination of C/ebpδ. In addition, deletion of the E3 ubiquitin ligase Cop1 in cancer cells stabilizes C/ebpδ to suppress expression of macrophage chemoattractant genes. Our integrated approach implicates Cop1 as a target for improving cancer immunotherapy efficacy in TNBC by regulating chemokine secretion and macrophage infiltration in the tumor microenvironment.

Regulatory association of long noncoding RNAs and chromatin accessibility facilitates erythroid differentiation

Erythroid differentiation is a dynamic process regulated by multiple factors, while the interaction between long non-coding RNAs and chromatin accessibility and its influence on erythroid differentiation remains unclear. To elucidate this interaction, we employed hematopoietic stem cells, multipotent progenitor cells, common myeloid progenitor cells, megakaryocyte-erythroid progenitor cells, and erythroblasts from human cord blood as an erythroid differentiation model to explore the coordinated regulatory functions of lncRNAs and chromatin accessibility by integrating RNA-Seq and ATAC-Seq data. We revealed that the integrated network of chromatin accessibility and lncRNAs exhibits stage-specific changes throughout the erythroid differentiation process, and that the changes at the EB stage of maturation are dramatic. We identified a subset of stage-specific lncRNAs and transcription factors (TFs) that associate with chromatin accessibility during erythroid differentiation, in which lncRNAs are key regulators of terminal erythroid differentiation via a lncRNA-TF-gene network. LncRNA PCED1B-AS1 was revealed to regulate terminal erythroid differentiation by coordinating GATA1 dynamically binding to the chromatin and interacting with cytoskeleton network during erythroid differentiation. DANCR, another lncRNA that is highly expressed at the MEP stage, was verified to promote erythroid differentiation by compromising megakaryocyte differentiation and coordinating with chromatin accessibility and TFs, such as RUNX1. Overall, our results identified the associated network of lncRNAs and chromatin accessibility in erythropoiesis and provide novel insights into erythroid differentiation and abundant resources for further study.

SIRT2-mediated deacetylation and deubiquitination of C/EBPβ prevents ethanol-induced liver injury

Protein acetylation has emerged to play pivotal roles in alcoholic liver disease (ALD). Sirutin 2 (SIRT2) is a nicotinamide adenine dinucleotide (NAD+)-dependent deacetylase involved in the regulation of aging, metabolism, and stress. However, the role of SIRT2 in ALD remains unclear. Here, we report that the SIRT2-mediated deacetylation-deubiquitination switch of CCAAT/enhancer-binding protein beta (C/EBPβ) prevents ALD. Our results showed that hepatic SIRT2 protein expression was negatively correlated with the severity of alcoholic liver injury in ALD patients. Liver-specific SIRT2 deficiency sensitized mice to ALD, whereas transgenic SIRT2 overexpression in hepatocytes significantly prevented ethanol-induced liver injury via normalization of hepatic steatosis, lipid peroxidation, and hepatocyte apoptosis. Mechanistically, we identified C/EBPβ as a critical substrate of SIRT2 implicated in ALD. SIRT2-mediated deacetylation at lysines 102 and 211 decreased C/EBPβ ubiquitination, resulting in enhanced protein stability and subsequently increased transcription of C/EBPβ-target gene LCN2. Importantly, hepatic deacetylated C/EBPβ and LCN2 compensation reversed SIRT2 deletion-induced ALD aggravation in mice. Furthermore, C/EBPβ protein expression was positively correlated with SIRT2 and LCN2 expression in the livers of ALD patients and was inversely correlated with ALD development. Therefore, activating SIRT2-C/EBPβ-LCN2 signaling pathway is a potential therapy for ALD.

The Expression and Prognostic Value of SUMO1-Activating Enzyme Subunit 1 and Its Potential Mechanism in Triple-Negative Breast Cancer

Background: Triple-negative breast cancer (TNBC) is the most invasive and metastatic subtype of breast cancer. SUMO1-activating enzyme subunit 1 (SAE1), an E1-activating enzyme, is indispensable for protein SUMOylation. SAE1 has been found to be a relevant biomarker for progression and prognosis in several tumor types. However, the role of SAE1 in TNBC remains to be elucidated.

Methods: In the research, the mRNA expression of SAE1 was analyzed via the cancer genome atlas (TCGA) and gene expression omnibus (GEO) database. Cistrome DB Toolkit was used to predict which transcription factors (TFs) are most likely to increase SAE1 expression in TNBC. The correlation between the expression of SAE1 and the methylation of SAE1 or quantity of tumor-infiltrating immune cells was further invested. Single-cell analysis, using CancerSEA, was performed to query which functional states are associated with SAE1 in different cancers in breast cancer at the single-cell level. Next, weighted gene coexpression network (WGCNA) was applied to reveal the highly correlated genes and coexpression networks of SAE1 in TNBC patients, and a prognostic model containing SAE1 and correlated genes was constructed. Finally, we also examined SAE1 protein expression of 207 TNBC tissues using immunohistochemical (IHC) staining.

Results: The mRNA and protein expression of SAE1 were increased in TNBC tissues compared with adjacent normal tissues, and the protein expression of SAE1 was significantly associated with overall survival (OS) and disease-free survival (DFS). Correlation analyses revealed that SAE1 expression was positively correlated with forkhead box M1 (FOXM1) TFs and negatively correlated with SAE1 methylation site (cg14042711) level. WGCNA indicated that the genes coexpressed with SAE1 belonged to the green module containing 1,176 genes. Through pathway enrichment analysis of the module, 1,176 genes were found enriched in cell cycle and DNA repair. Single-cell analysis indicated that SAE1 and its coexpression genes were associated with cell cycle, DNA damage, DNA repair, and cell proliferation. Using the LASSO COX regression, a prognostic model including SAE1 and polo-like kinase 1 (PLK1) was built to accurately predict the likelihood of DFS in TNBC patients.

Conclusion: In conclusion, we comprehensively analyzed the mRNA and protein expression, prognosis, and interaction genes of SAE1 in TNBC and constructed a prognostic model including SAE1 and PLK1. These results might be important for better understanding of the role of SAE1 in TNBC. In addition, DNA methyltransferase and TFs inhibitor treatments targeting SAE1 might improve the survival of TNBC patients.

Meta-analysis and Consolidation of Farnesoid X Receptor Chromatin Immunoprecipitation Sequencing Data Across Different Species and Conditions

Farnesoid X receptor (FXR) is a nuclear receptor that controls gene regulation of different metabolic pathways and represents an upcoming drug target for various liver diseases. Several data sets on genome-wide FXR binding in different species and conditions exist. We have previously reported that these data sets are heterogeneous and do not cover the full spectrum of potential FXR binding sites. Here, we report the first meta-analysis of all publicly available FXR chromatin immunoprecipitation sequencing (ChIP-seq) data sets from mouse, rat, and human across different conditions using a newly generated analysis pipeline. All publicly available single data sets were biocurated in a standardized manner and compared on every relevant level from raw reads to affected functional pathways. Individual murine data sets were then virtually merged into a single unique "FXR binding atlas" spanning all potential binding sites across various conditions. Comparison of the single biocurated data sets showed that the overlap of FXR binding sites between different species is modest and ranges from 48% (mouse-human) to 55% (mouse-rat). Moreover, in vivo data among different species are more similar than human in vivo data compared to human in vitro data. The consolidated murine global FXR binding atlas virtually increases sequencing depth and allows recovering more and novel potential binding sites and signaling pathways that were missed in the individual data sets. The FXR binding atlas is publicly searchable (https://fxratlas.tugraz.at). Conclusion: Published single FXR ChIP-seq data sets and large-scale integrated omics data sets do not cover the full spectrum of FXR binding. Combining different individual data sets and creating an "FXR super-binding atlas" enhances understanding of FXR signaling capacities across different conditions. This is important when considering the potential wide spectrum for drugs targeting FXR in liver diseases.

Subtype heterogeneity and epigenetic convergence in neuroendocrine prostate cancer

Neuroendocrine carcinomas (NEC) are tumors expressing markers of neuronal differentiation that can arise at different anatomic sites but have strong histological and clinical similarities. Here we report the chromatin landscapes of a range of human NECs and show convergence to the activation of a common epigenetic program. With a particular focus on treatment emergent neuroendocrine prostate cancer (NEPC), we analyze cell lines, patient-derived xenograft (PDX) models and human clinical samples to show the existence of two distinct NEPC subtypes based on the expression of the neuronal transcription factors ASCL1 and NEUROD1. While in cell lines and PDX models these subtypes are mutually exclusive, single-cell analysis of human clinical samples exhibits a more complex tumor structure with subtypes coexisting as separate sub-populations within the same tumor. These tumor sub-populations differ genetically and epigenetically contributing to intra- and inter-tumoral heterogeneity in human metastases. Overall, our results provide a deeper understanding of the shared clinicopathological characteristics shown by NECs. Furthermore, the intratumoral heterogeneity of human NEPCs suggests the requirement of simultaneous targeting of coexisting tumor populations as a therapeutic strategy.

Neuroendocrine carcinomas (NECs) arise from different anatomic sites, but have similar histological and clinical features. Here, the authors show that the epigenetic landscape of a range of NECs converges towards a common epigenetic state, while distinct subtypes occur within neuroendocrine prostate cancer contributing to intratumor heterogeneity in clinical samples.

Productive visualization of high-throughput sequencing data using the SeqCode open portable platform

Large-scale sequencing techniques to chart genomes are entirely consolidated. Stable computational methods to perform primary tasks such as quality control, read mapping, peak calling, and counting are likewise available. However, there is a lack of uniform standards for graphical data mining, which is also of central importance. To fill this gap, we developed SeqCode, an open suite of applications that analyzes sequencing data in an elegant but efficient manner. Our software is a portable resource written in ANSI C that can be expected to work for almost all genomes in any computational configuration. Furthermore, we offer a user-friendly front-end web server that integrates SeqCode functions with other graphical analysis tools. Our analysis and visualization toolkit represents a significant improvement in terms of performance and usability as compare to other existing programs. Thus, SeqCode has the potential to become a key multipurpose instrument for high-throughput professional analysis; further, it provides an extremely useful open educational platform for the world-wide scientific community. SeqCode website is hosted at http://ldicrocelab.crg.eu, and the source code is freely distributed at https://github.com/eblancoga/seqcode.

PFKFB4 is overexpressed in clear-cell renal cell carcinoma promoting pentose phosphate pathway that mediates Sunitinib resistance

Background

Kinases play critical role in clear-cell renal cell carcinoma (ccRCC). We aim to exploit novel kinase that is both protumorigenic and drugable in ccRCC.

Methods

Reproduction of public datasets with validation using microarray was performed to identify candidate gene. Functionality was studied using multi-omics with validation in vitro and in vivo.

Results

6-Phosphofructo-2-Kinase/Fructose-2,6-Biphosphatase 4 (PFKFB4) was differentially expressed showing significantly higher expression in tumor than in normal kidney. PFKFB4 overexpression was associated with advanced tumor grade, stage and worsened prognosis. PFKFB4-knockdown significantly impaired fitness in cell proliferation, migration and wound healing. Despite being recurrently deleted on 3p, PFKFN4 mRNA remained actively transcribed by HIF1α. Metabolomics showed overexpressed PFKFB4 showed enriched metabolites in pentose phosphate pathway (PPP). Phosphoproteomics and immunoprecipitation showed PFKFB4 also phosphorylated NCOA3 which interacted with FBP1 to counteract overactive PPP flux, forming a regulatory loop. PFKFB4-knockdown overcame resistance to Sunitinib in vitro and in vivo both in xenograft and tail-vein injection murine models.

Conclusion

We concluded PFKFB4 was associated with PPP activity and the fine-tuning of which was mediated by its phosphorylation of NCOA3. Targeting PFKFB4 held promise to combat resistance to Sunitinib.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-021-02103-5.

Transcription factors Krüppel-like factor 4 and paired box 5 regulate the expression of the Grainyhead-like genes

Genes from the Grainyhead-like (GRHL) family code for transcription factors necessary for the development and maintenance of various epithelia. These genes are also very important in the development of many types of cancer. However, little is known about the regulation of expression of GRHL genes. Previously, there were no systematic analyses of the promoters of GRHL genes or transcription factors that bind to these promoters. Here we report that the Krüppel-like factor 4 (KLF4) and the paired box 5 factor (PAX5) bind to the regulatory regions of the GRHL genes and regulate their expression. Ectopic expression of KLF4 or PAX5 alters the expression of GRHL genes. In KLF4-overexpressing HEK293 cells, the expression of GRHL1 and GRHL3 genes was upregulated by 32% and 60%, respectively, whereas the mRNA level of GRHL2 gene was lowered by 28% when compared to the respective controls. The levels of GRHL1 and GRHL3 expression were decreased by 30% or 33% in PAX5-overexpressing HEK293 cells. The presence of minor frequency allele of single nucleotide polymorphism rs115898376 in the promoter of the GRHL1 gene affected the binding of KLF4 to this site. The evidence presented here suggests an important role of KLF4 and PAX5 in the regulation of expression of GRHL1-3 genes.

VannoPortal: multiscale functional annotation of human genetic variants for interrogating molecular mechanism of traits and diseases

Interpreting the molecular mechanism of genomic variations and their causal relationship with diseases/traits are important and challenging problems in the human genetic study. To provide comprehensive and context-specific variant annotations for biologists and clinicians, here, by systematically integrating over 4TB genomic/epigenomic profiles and frequently-used annotation databases from various biological domains, we develop a variant annotation database, called VannoPortal. In general, the database has following major features: (i) systematically integrates 40 genome-wide variant annotations and prediction scores regarding allele frequency, linkage disequilibrium, evolutionary signature, disease/trait association, tissue/cell type-specific epigenome, base-wise functional prediction, allelic imbalance and pathogenicity; (ii) equips with our recent novel index system and parallel random-sweep searching algorithms for efficient management of backend databases and information extraction; (iii) greatly expands context-dependent variant annotation to incorporate large-scale epigenomic maps and regulatory profiles (such as EpiMap) across over 33 tissue/cell types; (iv) compiles many genome-scale base-wise prediction scores for regulatory/pathogenic variant classification beyond protein-coding region; (v) enables fast retrieval and direct comparison of functional evidence among linked variants using highly interactive web panel in addition to plain table; (vi) introduces many visualization functions for more efficient identification and interpretation of functional variants in single web page. VannoPortal is freely available at http://mulinlab.org/vportal.

Histone H3 lysine 27 acetylation profile undergoes two global shifts in undernourished children and suggests altered one-carbon metabolism

BACKGROUND

Stunting is a condition in which a child does not reach their full growth potential due to chronic undernutrition. It arises during the first 2 years of a child's life and is associated with developmental deficiencies and life-long health problems. Current interventions provide some benefit, but new approaches to prevention and treatment grounded in a molecular understanding of stunting are needed. Epigenetic analyses are critical as they can provide insight into how signals from a poor environment lead to changes in cell function.

RESULTS

Here we profiled histone H3 acetylation on lysine 27 (H3K27ac) in peripheral blood mononuclear cells (PBMCs) of 18-week-old (n = 14) and 1-year-old children (n = 22) living in an urban slum in Dhaka, Bangladesh. We show that 18-week-old children destined to become stunted have elevated levels of H3K27ac overall, functional analysis of which indicates activation of the immune system and stress response pathways as a primary response to a poor environment with high pathogen load. Conversely, overt stunting at 1-year-of age is associated with globally reduced H3K27ac that is indicative of metabolic rewiring and downregulation of the immune system and DNA repair pathways that are likely secondary responses to chronic exposure to a poor environment with limited nutrients. Among processes altered in 1-year-old children, we identified one-carbon metabolism, the significance of which is supported by integrative analysis with results from histone H3 trimethylation on lysine 4 (H3K4me3). Together, these results suggest altered one-carbon metabolism in this population of stunted children.

CONCLUSIONS

The epigenomes of stunted children undergo two global changes in H3K27ac within their first year of life, which are associated with probable initial hyperactive immune responses followed by reduced metabolic capacity. Limitation of one-carbon metabolites may play a key role in the development of stunting. Trial registration ClinicalTrials.gov NCT01375647. Registered 17 June 2011, retrospectively registered, https://clinicaltrials.gov/ct2/show/NCT01375647 .

Polymorphic estrogen receptor binding site causes Cd2-dependent sex bias in the susceptibility to autoimmune diseases

Complex autoimmune diseases are sexually dimorphic. An interplay between predisposing genetics and sex-related factors probably controls the sex discrepancy in the immune response, but the underlying mechanisms are unclear. Here we positionally identify a polymorphic estrogen receptor binding site that regulates Cd2 expression, leading to female-specific differences in T cell-dependent mouse models of autoimmunity. Female mice with reduced Cd2 expression have impaired autoreactive T cell responses. T cells lacking Cd2 costimulation upregulate inhibitory Lag-3. These findings help explain sexual dimorphism in human autoimmunity, as we find that CD2 polymorphisms are associated with rheumatoid arthritis and 17-β-estradiol-regulation of CD2 is conserved in human T cells. Hormonal regulation of CD2 might have implications for CD2-targeted therapy, as anti-Cd2 treatment more potently affects T cells in female mice. These results demonstrate the relevance of sex-genotype interactions, providing strong evidence for CD2 as a sex-sensitive predisposing factor in autoimmunity.

Integrative Epigenome Map of the Normal Human Prostate Provides Insights Into Prostate Cancer Predisposition

Cells of all tissues in the human body share almost the exact same DNA sequence, but the epigenomic landscape can be drastically distinct. To improve our understanding of the epigenetic abnormalities in prostate-related diseases, it is important to use the epigenome of normal prostate as a reference. Although previous efforts have provided critical insights into the genetic and transcriptomic features of the normal prostate, a comprehensive epigenome map has been lacking. To address this need, we conducted a Roadmap Epigenomics legacy project integrating six histone marks (H3K4me1, H3K4me3, H3K9me3, H3K36me3, H3K27me3, and H3K27ac) with complete DNA methylome, transcriptome, and chromatin accessibility data to produce a comprehensive epigenome map of normal prostate tissue. Our epigenome map is composed of 18 chromatin states each with unique signatures of DNA methylation, chromatin accessibility, and gene expression. This map provides a high-resolution comprehensive annotation of regulatory regions of the prostate, including 105,593 enhancer and 70,481 promoter elements, which account for 5.3% of the genome. By comparing with other epigenomes, we identified 7,580 prostate-specific active enhancers associated with prostate development. Epigenomic annotation of GWAS SNPs associated with prostate cancers revealed that two out of nine SNPs within prostate enhancer regions destroyed putative androgen receptor (AR) binding motif. A notable SNP rs17694493, might decouple AR's repressive effect on CDKN2B-AS1 and cell cycle regulation, thereby playing a causal role in predisposing cancer risk. The comprehensive epigenome map of the prostate is valuable for investigating prostate-related diseases.

A prelude to the proximity interaction mapping of CXXC5

CXXC5 is a member of the zinc-finger CXXC family proteins that interact with unmodified CpG dinucleotides through a conserved ZF-CXXC domain. CXXC5 is involved in the modulation of gene expressions that lead to alterations in diverse cellular events. However, the underlying mechanism of CXXC5-modulated gene expressions remains unclear. Proteins perform their functions in a network of proteins whose identities and amounts change spatiotemporally in response to various stimuli in a lineage-specific manner. Since CXXC5 lacks an intrinsic transcription regulatory function or enzymatic activity but is a DNA binder, CXXC5 by interacting with proteins could act as a scaffold to establish a chromatin state restrictive or permissive for transcription. To initially address this, we utilized the proximity-dependent biotinylation approach. Proximity interaction partners of CXXC5 include DNA and chromatin modifiers, transcription factors/co-regulators, and RNA processors. Of these, CXXC5 through its CXXC domain interacted with EMD, MAZ, and MeCP2. Furthermore, an interplay between CXXC5 and MeCP2 was critical for a subset of CXXC5 target gene expressions. It appears that CXXC5 may act as a nucleation factor in modulating gene expressions. Providing a prelude for CXXC5 actions, our results could also contribute to a better understanding of CXXC5-mediated cellular processes in physiology and pathophysiology.

Transcriptional and epigenetic landscape of Ca2+-signaling genes in hepatocellular carcinoma

Calcium (Ca2+) signaling has a major role in regulating a wide range of cellular mechanisms, including gene expression, proliferation, metabolism, cell death, muscle contraction, among others. Recent evidence suggests that ~ 1600 genes are related to the Ca2+ signaling. Some of these genes' expression is altered in several pathological conditions, including different cancer types, and epigenetic mechanisms are involved. However, their expression and regulation in hepatocellular carcinoma (HCC) and the liver are barely known. Here, we aimed to explore the expression of genes involved in the Ca2+-signaling in HCC, liver regeneration, and hepatocyte differentiation, and whether their expression is regulated by epigenetic mechanisms such as DNA methylation and histone posttranslational modifications (HPM). Results show that several Ca2+-signaling genes' expression is altered in HCC samples; among these, a subset of twenty-two correlate with patients' survival. DNA methylation correlates with eight of these genes' expression, and Guadecitabine, a hypomethylating agent, regulates the expression of seven down-regulated and three up-regulated genes in HepG2 cells. The down-regulated genes displayed a marked decrease of euchromatin histone marks, whereas up-regulated genes displayed gain in these marks. Additionally, the expression of these genes is modulated during liver regeneration and showed similar profiles between in vitro differentiated hepatocytes and liver-derived hepatocytes. In conclusion, some components of the Ca2+-signaling are altered in HCC and displayed a correlation with patients' survival. DNA methylation and HMP are an attractive target for future investigations to regulate their expression. Ca2+-signaling could be an important regulator of cell proliferation and differentiation in the liver.

Systematic identification of non-canonical transcription factor motifs

Sequence-specific transcription factors (TFs) recognize motifs of related nucleotide sequences at their DNA binding sites. Upon binding at these sites, TFs regulate critical molecular processes such as gene expression. It is widely assumed that a TF recognizes a single “canonical” motif, although recent studies have identified additional “non-canonical” motifs for some TFs. A comprehensive approach to identify non-canonical DNA binding motifs and the functional importance of those motifs’ matches in the human genome is necessary for fully understanding the mechanisms of TF-regulated molecular processes in human cells. To address this need, we developed a statistical pipeline for in vitro HT-SELEX data that identifies and characterizes the distributions of non-canonical TF motifs in a stringent manner. Analyzing ~170 human TFs’ HT-SELEX data, we found non-canonical motifs for 19 TFs (11%). These non-canonical motifs occur independently of the TFs’ canonical motifs. Non-canonical motif occurrences in the human genome show similar evolutionary conservation to canonical motif occurrences, explain TF binding in locations without canonical motifs, and occur within gene promoters and epigenetically marked regulatory sequences in human cell lines and tissues. Our approach and collection of non-canonical motifs expand current understanding of functionally relevant DNA binding sites for human TFs.

Upregulation of BIRC5 plays essential role in esophageal squamous cell carcinoma

BACKGROUND

Esophageal squamous cell carcinoma (ESCC) is one of the most common cancers in the world, the detection and prognosis of which are still unsatisfactory. Thus, it is essential to explore the factors that may identify ESCC and evaluate the prognosis of ESCC patients.

RESULTS

Both protein and mRNA expression levels of BIRC5 are upregulated in ESCC group rather than non-ESCC group (standardized mean difference > 0). BIRC5 mRNA expression is related to the age, tumor location, lymph node stage and clinical stage of ESCC patients (p < 0.05). BIRC5 expression makes it feasible to distinguish ESCC from non-ESCC (area under the curve > 0.9), and its high expression is related to poor prognosis of ESCC patients (restrictive survival time difference = -0.036, p < 0.05). BIRC5 may play an important role in ESCC by influencing the cell cycle pathway, and CDK1, MAD2L and CDC20 may be the hub genes of this pathway. The transcription factors-MAZ and TFPD1 -are likely to regulate the transcription of BIRC5, which may be one of the factors for the high expression of BIRC5 in ESCC.

CONCLUSIONS

The current study shows that upregulation of BIRC5 may have essential clinical value in ESCC, and contributes to the understanding of the pathogenesis of ESCC.

ELK4 promotes the development of gastric cancer by inducing M2 polarization of macrophages through regulation of the KDM5A-PJA2-KSR1 axis

Background

We tried to elaborate the molecular mechanism of ETS-like transcription factor 4 (ELK4) affecting gastric cancer (GC) progression through M2 polarization of macrophages mediated by lysine-specific demethylase 5A (KDM5A)-Praja2 (PJA2)-kinase suppressor of ras 1 (KSR1) axis.

Methods

GC expression dataset was obtained from GEO database, and the downstream regulatory mechanism of ELK4 was predicted. Tumor-associated macrophages (TAMs) were isolated from GC tissues. The interaction among ELK4, KDM5A, PJA2 and KSR1 was analyzed by dual luciferase reporter gene, ChIP and Co-IP assays. The stability of KSR1 protein was detected by cycloheximide (CHX) treatment. After TAMs were co-cultured with HGC-27 cells, HGC-27 cell biological processes were assessed through gain- and loss-of function assays. Tumorigenicity was detected by tumorigenicity test in nude mice.

Results

In GC and TAMs, ELK4, KDM5A and KSR1 were highly expressed, while PJA2 was lowly expressed. M2 polarization of macrophages promoted the development of GC. ELK4 activated KDM5A by transcription and promoted macrophage M2 polarization. KDM5A inhibited the expression of PJA2 by removing H3K4me3 of PJA2 promoter, which promoted M2 polarization of macrophages. PJA2 reduced KSR1 by ubiquitination. ELK4 promoted the proliferative, migrative and invasive potentials of GC cells as well as the growth of GC xenografts by regulating KSR1.

Conclusion

ELK4 may reduce the PJA2-dependent inhibition of KSR1 by transcriptional activation of KDM5A to promote M2 polarization of macrophages, thus promoting the development of GC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-02915-1.

Transcription factor RFX7 governs a tumor suppressor network in response to p53 and stress

Despite its prominence, the mechanisms through which the tumor suppressor p53 regulates most genes remain unclear. Recently, the regulatory factor X 7 (RFX7) emerged as a suppressor of lymphoid neoplasms, but its regulation and target genes mediating tumor suppression remain unknown. Here, we identify a novel p53-RFX7 signaling axis. Integrative analysis of the RFX7 DNA binding landscape and the RFX7-regulated transcriptome in three distinct cell systems reveals that RFX7 directly controls multiple established tumor suppressors, including PDCD4, PIK3IP1, MXD4, and PNRC1, across cell types and is the missing link for their activation in response to p53 and stress. RFX7 target gene expression correlates with cell differentiation and better prognosis in numerous cancer types. Interestingly, we find that RFX7 sensitizes cells to Doxorubicin by promoting apoptosis. Together, our work establishes RFX7’s role as a ubiquitous regulator of cell growth and fate determination and a key node in the p53 transcriptional program.

A CpG island promoter drives the CXXC5 gene expression

CXXC5 is a member of the zinc-finger CXXC family that binds to unmethylated CpG dinucleotides. CXXC5 modulates gene expressions resulting in diverse cellular events mediated by distinct signaling pathways. However, the mechanism responsible for CXXC5 expression remains largely unknown. We found here that of the 14 annotated CXXC5 transcripts with distinct 5' untranslated regions encoding the same protein, transcript variant 2 with the highest expression level among variants represents the main transcript in cell models. The DNA segment in and at the immediate 5'-sequences of the first exon of variant 2 contains a core promoter within which multiple transcription start sites are present. Residing in a region with high G-C nucleotide content and CpG repeats, the core promoter is unmethylated, deficient in nucleosomes, and associated with active RNA polymerase-II. These findings suggest that a CpG island promoter drives CXXC5 expression. Promoter pull-down revealed the association of various transcription factors (TFs) and transcription co-regulatory proteins, as well as proteins involved in histone/chromatin, DNA, and RNA processing with the core promoter. Of the TFs, we verified that ELF1 and MAZ contribute to CXXC5 expression. Moreover, the first exon of variant 2 may contain a G-quadruplex forming region that could modulate CXXC5 expression.

Nuclear FGFR1 Regulates Gene Transcription and Promotes Antiestrogen Resistance in ER+ Breast Cancer

PURPOSE

FGFR1 overexpression has been associated with endocrine resistance in ER⁺ breast cancer. We found FGFR1 localized in the nucleus of breast cancer cells in primary tumors resistant to estrogen suppression. We investigated a role of nuclear FGFR1 on gene transcription and antiestrogen resistance.

EXPERIMENTAL DESIGN

Tumors from patients treated with letrozole were subjected to Ki67 and FGFR1 IHC. MCF7 cells were transduced with FGFR1(SP-)(NLS) to promote nuclear FGFR1 overexpression. FGFR1 genomic activity in ER⁺/FGFR1-amplified breast cancer cells ± FOXA1 siRNA or ± the FGFR tyrosine kinase inhibitor (TKI) erdafitinib was examined by chromatin immunoprecipitation sequencing (ChIP-seq) and RNA sequencing (RNA-seq). The nuclear and chromatin-bound FGFR1 interactome was investigated by mass spectrometry (MS).

RESULTS

High nuclear FGFR1 expression in ER⁺ primary tumors positively correlated with post-letrozole Ki67 values. Nuclear FGFR1 overexpression influenced gene transcription and promoted resistance to estrogen suppression and to fulvestrant in vivo. A gene expression signature induced by nuclear FGFR1 correlated with shorter survival in the METABRIC cohort of patients treated with antiestrogens. ChIP-Seq revealed FGFR1 occupancy at transcription start sites, overlapping with active transcription histone marks. MS analysis of the nuclear FGFR1 interactome identified phosphorylated RNA-Polymerase II and FOXA1, with FOXA1 RNAi impairing FGFR1 recruitment to chromatin. Treatment with erdafitinib did not impair nuclear FGFR1 translocation and genomic activity.

CONCLUSIONS

These data suggest nuclear FGFR1 contributes to endocrine resistance by modulating gene transcription in ER⁺ breast cancer. Nuclear FGFR1 activity was unaffected by FGFR TKIs, thus supporting the development of treatment strategies to inhibit nuclear FGFR1 in ER+/FGFR1 overexpressing breast cancer.

A novel RNA binding protein-associated prognostic model to predict overall survival in hepatocellular carcinoma patients

Hepatocellular carcinoma (HCC) is 1 of the deadliest malignancies worldwide. Despite significant advances in diagnosis and treatment, the mortality rate from HCC persists at a substantial level. Construction of a prognostic model that can reliably predict HCC patients' overall survival is urgently needed.Two RNA-seq dataset (the Cancer Genome Atlas and International Cancer Genome Consortium) and 1 microarray dataset (GSE14520) were included in our study. RNA-binding proteins (RBPs) in HCC patients was examined by differentially expressed genes analysis, functional enrichment analysis and protein-protein interaction network analysis. Subsequently, the Cancer Genome Atlas dataset was randomly divided into training and testing cohort with a prognostic model developed in the training cohort. In order to evaluate the prognostic value of the model, a comprehensive survival assessment was conducted.Five RBPs (ribosomal protein L10-like, enhancer of zeste homolog 2 (EZH2), peroxisome proliferator-activated receptor gamma coactivator 1 alpha (PPARGC1A), zinc finger protein 239, interferon-induced protein with tetratricopeptide repeats 1) were used to construct the model. The model accurately predicted the prognosis of liver cancer patients in both the training cohort and validation cohort. HCC patients could be assigned into a high-risk group and a low-risk group by this model, and the overall survival of these 2 groups was significantly different (P  < .05). Furthermore, the risk scores obtained by this model were highly correlated with immune cell infiltration.The prognostic model helps to identify HCC patients at high risk of mortality, which optimizes decision-making for individualized treatment.

Modest changes in Spi1 dosage reveal the potential for altered microglial function as seen in Alzheimer's disease

Genetic association studies have identified multiple variants at the SPI1 locus that modify risk and age of onset for Alzheimer's Disease (AD). Reports linking risk variants to gene expression suggest that variants denoting higher SPI1 expression are likely to have an earlier AD onset, and several other AD risk genes contain PU.1 binding sites in the promoter region. Overall, this suggests the level of SPI1 may alter microglial phenotype potentially impacting AD. This study determined how the microglial transcriptome was altered following modest changes to Spi1 expression in primary mouse microglia. RNA-sequencing was performed on microglia with reduced or increased Spi1/PU.1 expression to provide an unbiased approach to determine transcriptomic changes affected by Spi1. In summary, a reduction in microglial Spi1 resulted in the dysregulation of transcripts encoding proteins involved in DNA replication pathways while an increased Spi1 results in an upregulation of genes associated with immune response pathways. Additionally, a subset of 194 Spi1 dose-sensitive genes was identified and pathway analysis suggests that several innate immune and interferon response pathways are impacted by the concentration of Spi1. Together these results suggest Spi1 levels can alter the microglial transcriptome and suggests interferon pathways may be altered in individuals with AD related Spi1 risk SNPs.

Comprehensive analysis of BTN3A1 in cancers: mining of omics data and validation in patient samples and cellular models

Butyrophilin 3A1 (BTN3A1), a major histocompatibility complex‐associated gene that encodes a membrane protein with two extracellular immunoglobulin domains and an intracellular B30.2 domain, is critical in T‐cell activation and adaptive immune response. Here, the expression of BTN3A1 in cancers was analyzed in eight databases comprising 86 733 patients of 33 cancers, and the findings were validated in patient samples and cell models. We showed that BTN3A1 was expressed in most cancers, and its expression level was strongly correlated with clinical outcome of 13 cancers. Mutations of BTN3A1 were detected, and the mutations were distributed throughout the entire gene. Gene set enrichment analysis showed that BTN3A1 co‐expression genes and interacting proteins were enriched in immune regulation‐related pathways. BTN3A1 was associated with tumor‐infiltrating immune cells and was co‐expressed with multiple immune checkpoints in patients with breast cancer (BRCA) and non‐small cell lung cancer (NSCLC). We reported that BTN3A1 was downregulated in 46 of 65 (70.8%) NSCLCs, and its expression level was inversely associated with clinical outcome of the patients. BTN3A1 in tumor samples was lower than in counterpart normal tissues in 31 of 38 (81.6%) BRCAs. Bioinformatics analyses showed that BTN3A1 could be a target gene of transcription factor Spi‐1 proto‐oncogene (SPI1), and our ‘wet’ experiments showed that ectopic expression of SPI1 upregulated, whereas silencing of SPI1 downregulated, BTN3A1 expression in cells. These results suggest that BTN3A1 may function as a tumor suppressor and may serve as a potential prognostic biomarker in NSCLCs and BRCAs.