High-resolution Hi-C maps highlight multiscale 3D epigenome reprogramming during pancreatic cancer metastasis

A review on trends in development and translation of omics signatures in cancer

The field of cancer genomics and transcriptomics has evolved from targeted profiling to swift sequencing of individual tumor genome and transcriptome. The steady growth in genome, epigenome, and transcriptome datasets on a genome-wide scale has significantly increased our capability in capturing signatures that represent both the intrinsic and extrinsic biological features of tumors. These biological differences can help in precise molecular subtyping of cancer, predicting tumor progression, metastatic potential, and resistance to therapeutic agents. In this review, we summarized the current development of genomic, methylomic, transcriptomic, proteomic and metabolic signatures in the field of cancer research and highlighted their potentials in clinical applications to improve diagnosis, prognosis, and treatment decision in cancer patients.

Histone Lactylation Dynamics: Unlocking the Triad of Metabolism, Epigenetics, and Immune Regulation in Metastatic Cascade of Pancreatic Cancer

Cancer cells rewire metabolism to sculpt the immune tumor microenvironment (TME) and propel tumor advancement, which intricately tied to post-translational modifications. Histone lactylation has emerged as a novel player in modulating protein functions, whereas little is known about its pathological role in pancreatic ductal adenocarcinoma (PDAC) progression. Employing a multi-omics approach encompassing bulk and single-cell RNA sequencing, metabolomics, ATAC-seq, and CUT&Tag methodologies, we unveiled the potential of histone lactylation in prognostic prediction, patient stratification and TME characterization. Notably, "LDHA-H4K12la-immuno-genes" axis has introduced a novel node into the regulatory framework of "metabolism-epigenetics-immunity," shedding new light on the landscape of PDAC progression. Furthermore, the heightened interplay between cancer cells and immune counterparts via Nectin-2 in liver metastasis with elevated HLS unraveled a positive feedback loop in driving immune evasion. Simultaneously, immune cells exhibited altered HLS and autonomous functionality across the metastatic cascade. Consequently, the exploration of innovative combination strategies targeting the metabolism-epigenetics-immunity axis holds promise in curbing distant metastasis and improving survival prospects for individuals grappling with challenges of PDAC.

Cell type- and transcription-independent spatial proximity between enhancers and promoters

Cell type-specific enhancers are critically important for lineage specification. The mechanisms that determine cell-type specificity of enhancer activity, however, are not fully understood. Most current models for how enhancers function invoke physical proximity between enhancer elements and their target genes. Here, we use an imaging-based approach to examine the spatial relationship of cell type-specific enhancers and their target genes with single-cell resolution. Using high-throughput microscopy, we measure the spatial distance from target promoters to their cell type-specific active and inactive enhancers in individual pancreatic cells derived from distinct lineages. We find increased proximity of all promoter-enhancer pairs relative to non-enhancer pairs separated by similar genomic distances. Strikingly, spatial proximity between enhancers and target genes was unrelated to tissue-specific enhancer activity. Furthermore, promoter-enhancer proximity did not correlate with the expression status of target genes. Our results suggest that promoter-enhancer pairs exist in a distinctive chromatin environment but that genome folding is not a universal driver of cell-type specificity in enhancer function.

LTR retrotransposon-derived LncRNA LINC01446 promotes hepatocellular carcinoma progression and angiogenesis by regulating the SRPK2/SRSF1/VEGF axis

The causal link between long terminal repeat (LTR) retrotransposon-derived lncRNAs and hepatocellular carcinoma (HCC) remains elusive and whether these cancer-exclusive lncRNAs contribute to the effectiveness of current HCC therapies is yet to explore. Here, we investigated the activation of LTR retrotransposon-derived lncRNAs in a broad range of liver diseases. We found that LTR retrotransposon-derived lncRNAs are mainly activated in HCC and is correlated with the proliferation status of HCC. Furthermore, we discovered that an LTR retrotransposon-derived lncRNA, LINC01446, exhibits specific expression in HCC. HCC patients with higher LINC01446 expression had shorter overall survival times. In vitro and in vivo assays showed that LINC01446 promoted HCC growth and angiogenesis. Mechanistically, LINC01446 bound to serine/arginine protein kinase 2 (SRPK2) and activated its downstream target, serine/arginine splicing factor 1 (SRSF1). Furthermore, activation of the SRPK2-SRSF1 axis increased the splicing and expression of VEGF isoform A165 (VEGFA165). Notably, inhibiting LINC01446 expression dramatically impaired tumor growth in vivo and resulted in better therapeutic outcomes when combined with antiangiogenic agents. In addition, we found that the transcription factor MESI2 bound to the cryptic MLT2B3 LTR promoter and drove LINC01446 transcription in HCC cells. Taken together, our findings demonstrate that LTR retrotransposon-derived LINC01446 promotes the progression of HCC by activating the SRPK2/SRSF1/VEGFA165 axis and highlight targeting LINC01446 as a potential therapeutic strategy for HCC patients.

3D Chromatin Alteration by Disrupting β-Catenin/CBP Interaction Is Enriched with Insulin Signaling in Pancreatic Cancer

The therapeutic potential of targeting the β-catenin/CBP interaction has been demonstrated in a variety of preclinical tumor models with a small molecule inhibitor, ICG-001, characterized as a β-catenin/CBP antagonist. Despite the high binding specificity of ICG-001 for the N-terminus of CBP, this β-catenin/CBP antagonist exhibits pleiotropic effects. Our recent studies found global changes in three-dimensional (3D) chromatin architecture in response to disruption of the β-catenin/CBP interaction in pancreatic cancer cells. However, an understanding of how the functional crosstalk between the antagonist and the β-catenin/CBP interaction affects changes in 3D chromatin architecture and, thereby, gene expression and downstream effects remains to be elucidated. Here, we perform Hi-C analyses on canonical and patient-derived pancreatic cancer cells before and after treatment with ICG-001. In addition to global alteration of 3D chromatin domains, we unexpectedly identify insulin signaling genes enriched in the altered chromatin domains. We further demonstrate that the chromatin loops associated with insulin signaling genes are significantly weakened after ICG-001 treatment. We finally elicit the deletion of a looping of IRS1-a key insulin signaling gene-significantly impeding pancreatic cancer cell growth, indicating that looping-mediated insulin signaling might act as an oncogenic pathway to promote pancreatic cancer progression. Our work shows that targeting aberrant insulin chromatin looping in pancreatic cancer might provide a therapeutic benefit.

Immune-related signature identifies IL1R2 as an immunological and prognostic biomarker in pancreatic cancer

Objective

Pancreatic cancer is one of the most aggressive malignancies, a robust prognostic signature and novel biomarkers are urgently needed for accurate stratification of the patients and optimization of clinical decision-making.

Methods

A list of bioinformatic analysis were applied in public dataset to construct an immune-related signature. Furthermore, the most pivotal gene in the signature was identified. The potential mechanism of the core gene function was revealed through GSEA, CIBERSORT, ESTIMATE, immunophenoscore (IPS) algorithm, single-cell analysis, and functional experiment.

Results

An immune-related prognostic signature and associated nomogram were constructed and validated. Among the genes constituting the signature, interleukin 1 receptor type II (IL1R2) was identified as the gene occupying the most paramount position in the risk signature. Meanwhile, knockdown of IL1R2 significantly inhibited the proliferation, invasion, and migration ability of pancreatic cancer cells. Additionally, high IL1R2 expression was associated with reduced CD8+ T cell infiltration in pancreatic cancer microenvironment, which may be due to high programmed cell death-ligand-1 (PD-L1) expression in cancer cells. Finally, the IPS algorithm proved that patients with high IL1R2 expression possessed a higher tumor mutation burden and a higher probability of benefiting from immunotherapy.

Conclusion

In conclusion, our study constructed an efficient immune-related prognostic signature and identified the key role of IL1R2 in the development of pancreatic cancer, as well as its potential to serve as a biomarker for immunotherapy efficacy prediction for pancreatic cancer.

Augmentation of hepatocellular carcinoma malignancy by annexin A5 through modulation of invasion and angiogenesis

BACKGROUND

Hepatocellular carcinoma (HCC) continues to play a substantial role in cancer-related morbidity and mortality, largely owing to its pronounced tumor heterogeneity and propensity for recurrence. This underscores the pressing need for in-depth examination of its highly malignant mechanisms. Annexin A5 (ANXA5), recognized as a hallmark tumor protein, has emerged as a focal point of interest because of its ambiguous function and mechanism in HCC prognosis. This study aimed to provide a comprehensive understanding of the role of ANXA5 in the malignant progression of human HCC cells by employing an integrative approach that combines conventional experimental methods with RNA sequencing.

METHODS

Differences in ANXA5 expression between HCC tissues and corresponding nontumor tissues were evaluated using immunofluorescence (n = 25). Correlation analysis was subsequently performed to assess the association between ANXA5 expression and clinicopathological features (n = 65). The role of ANXA5 in human HCC cell lines with ANXA5 gene knockout and overexpression was explored in vitro using migration and invasion assays and Ki-67 indices and in vivo based on node mice xenograft model. A tube formation assay using human umbilical vein endothelial cells (HUVECs) was conducted to demonstrate the angiogenic effects of ANXA5 in HCC. Single-cell and bulk RNA sequencing was used to further investigate the underlying mechanisms involved.

RESULTS

This study revealed that ANXA5 is highly expressed in patients with HCC and correlates with poor prognosis. Assays for migration, invasion, and proliferation based on ANXA5 gene knockout and overexpression systems in human HCC cell lines have demonstrated that ANXA5 enhances HCC malignancy in vitro and in vivo. Tube formation assays of HUVECs indicated that ANXA5 facilitates angiogenesis and recruits endothelial cells to HCC cells. Single-cell and bulk RNA sequencing data analysis further confirmed that ANXA5 expression in HCC is associated with hepatocyte metabolism, immune response activation, and various oncogenic signaling pathways.

CONCLUSIONS

This study revealed a meaningful association between elevated ANXA5 expression in tumor tissues and an unfavorable prognosis in patients with HCC. In addition, ANXA5 promotes HCC malignancy by promoting invasion and angiogenesis. Thus, ANXA5 has emerged as a promising therapeutic target for HCC and has the potential to improve patient outcomes.

Crosstalk between metabolic remodeling and epigenetic reprogramming: A new perspective on pancreatic cancer

Pancreatic cancer is a highly malignant solid tumor with a poor prognosis and a high mortality rate. Thus, exploring the mechanisms underlying the development and progression of pancreatic cancer is critical for identifying targets for diagnosis and treatment. Two important hallmarks of cancer-metabolic remodeling and epigenetic reprogramming-are interconnected and closely linked to regulate one another, creating a complex interaction landscape that is implicated in tumorigenesis, invasive metastasis, and immune escape. For example, metabolites can be involved in the regulation of epigenetic enzymes as substrates or cofactors, and alterations in epigenetic modifications can in turn regulate the expression of metabolic enzymes. The crosstalk between metabolic remodeling and epigenetic reprogramming in pancreatic cancer has gained considerable attention. Here, we review the emerging data with a focus on the reciprocal regulation of metabolic remodeling and epigenetic reprogramming. We aim to highlight how these mechanisms could be applied to develop better therapeutic strategies.

Alteration of chromosome structure impacts gene expressions implicated in pancreatic ductal adenocarcinoma cells

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) is a lethal malignancy with a five-year survival rate of approximately 10%. Genetic mutations are pivotal drivers in PDAC pathogenesis, but recent investigations also revealed the involvement of non-genetic alterations in the disease development. In this study, we undertook a multi-omics approach, encompassing ATAC-seq, RNA-seq, ChIP-seq, and Hi-C methodologies, to dissect gene expression alterations arising from changes in chromosome accessibility and chromatin three-dimensional interactions in PDAC.

RESULTS

Our findings indicate that chromosomal structural alterations can lead to abnormal expressions on key genes during PDAC development. Notably, overexpression of oncogenes FGFR2, FOXA2, CYP2R1, and CPOX can be attributed to the augmentation of promoter accessibility, coupled with long-range interactions with distal elements. Additionally, our findings indicate that chromosomal structural alterations caused by genomic instability can lead to abnormal expressions in PDACs. As an example, by analyzing chromosomal changes, we identified a putative oncogenic gene, LPAR1, which shows upregulated expression in both PDAC cell lines and clinical samples. The overexpression is correlated with alterations in LPAR1-associated 3D genome structure and chromatin state. We further demonstrated that high LPAR1 activity is required for enhanced PDAC cell migration in vitro.

CONCLUSIONS

Collectively, our findings reveal that the chromosomal conformational alterations, in addition to the well-known genetic mutations, are critical for PDAC tumorigenesis.

Chromosome compartmentalization: causes, changes, consequences, and conundrums

The spatial segregation of the genome into compartments is a major feature of 3D genome organization. New data on mammalian chromosome organization across different conditions reveal important information about how and why these compartments form and change. A combination of epigenetic state, nuclear body tethering, physical forces, gene expression, and replication timing (RT) can all influence the establishment and alteration of chromosome compartments. We review the causes and implications of genomic regions undergoing a 'compartment switch' that changes their physical associations and spatial location in the nucleus. About 20-30% of genomic regions change compartment during cell differentiation or cancer progression, whereas alterations in response to a stimulus within a cell type are usually much more limited. However, even a change in 1-2% of genomic bins may have biologically relevant implications. Finally, we review the effects of compartment changes on gene regulation, DNA damage repair, replication, and the physical state of the cell.

Epigenetic regulation in cancer

Epigenetic modifications are defined as heritable changes in gene activity that do not involve changes in the underlying DNA sequence. The oncogenic process is driven by the accumulation of alterations that impact genome's structure and function. Genetic mutations, which directly disrupt the DNA sequence, are complemented by epigenetic modifications that modulate gene expression, thereby facilitating the acquisition of malignant characteristics. Principals among these epigenetic changes are shifts in DNA methylation and histone mark patterns, which promote tumor development and metastasis. Notably, the reversible nature of epigenetic alterations, as opposed to the permanence of genetic changes, positions the epigenetic machinery as a prime target in the discovery of novel therapeutics. Our review delves into the complexities of epigenetic regulation, exploring its profound effects on tumor initiation, metastatic behavior, metabolic pathways, and the tumor microenvironment. We place a particular emphasis on the dysregulation at each level of epigenetic modulation, including but not limited to, the aberrations in enzymes responsible for DNA methylation and histone modification, subunit loss or fusions in chromatin remodeling complexes, and the disturbances in higher-order chromatin structure. Finally, we also evaluate therapeutic approaches that leverage the growing understanding of chromatin dysregulation, offering new avenues for cancer treatment.

Topological reorganization and functional alteration of distinct genomic components in gallbladder cancer

Altered three-dimensional architecture of chromatin influences various genomic regulators and subsequent gene expression in human cancer. However, knowledge of the topological rearrangement of genomic hierarchical layers in cancer is largely limited. Here, by taking advantage of in situ Hi-C, RNA-sequencing, and chromatin immunoprecipitation sequencing (ChIP-seq), we investigated structural reorganization and functional changes in chromosomal compartments, topologically associated domains (TADs), and CCCTC binding factor (CTCF)-mediated loops in gallbladder cancer (GBC) tissues and cell lines. We observed that the chromosomal compartment A/B switch was correlated with CTCF binding levels and gene expression changes. Increased inter-TAD interactions with weaker TAD boundaries were identified in cancer cell lines relative to normal controls. Furthermore, the chromatin short loops and cancer unique loops associated with chromatin remodeling and epithelial-mesenchymal transition activation were enriched in cancer compared with their control counterparts. Cancer-specific enhancer-promoter loops, which contain multiple transcription factor binding motifs, acted as a central element to regulate aberrant gene expression. Depletion of individual enhancers in each loop anchor that connects with promoters led to the inhibition of their corresponding gene expressions. Collectively, our data offer the landscape of hierarchical layers of cancer genome and functional alterations that contribute to the development of GBC.

Cancer epigenetics: from laboratory studies and clinical trials to precision medicine

Epigenetic dysregulation is a common feature of a myriad of human diseases, particularly cancer. Defining the epigenetic defects associated with malignant tumors has become a focus of cancer research resulting in the gradual elucidation of cancer cell epigenetic regulation. In fact, most stages of tumor progression, including tumorigenesis, promotion, progression, and recurrence are accompanied by epigenetic alterations, some of which can be reversed by epigenetic drugs. The main objective of epigenetic therapy in the era of personalized precision medicine is to detect cancer biomarkers to improve risk assessment, diagnosis, and targeted treatment interventions. Rapid technological advancements streamlining the characterization of molecular epigenetic changes associated with cancers have propelled epigenetic drug research and development. This review summarizes the main mechanisms of epigenetic dysregulation and discusses past and present examples of epigenetic inhibitors in cancer diagnosis and treatment, with an emphasis on the development of epigenetic enzyme inhibitors or drugs. In the final part, the prospect of precise diagnosis and treatment is considered based on a better understanding of epigenetic abnormalities in cancer.

Reprogramming of 3D chromatin domains by antagonizing the β-catenin/CBP interaction attenuates insulin signaling in pancreatic cancer

The therapeutic potential of targeting the β-catenin/CBP interaction has been demonstrated in a variety of preclinical tumor models with a small molecule inhibitor, ICG-001, characterized as a β-catenin/CBP antagonist. Despite the high binding specificity of ICG-001 for the N-terminus of CBP, this β-catenin/CBP antagonist exhibits pleiotropic effects. Our recent studies found global changes in three-dimensional (3D) chromatin architecture in response to disruption of the β-catenin/CBP interaction in pancreatic cancer cells. However, an understanding of the functional crosstalk between antagonizing the β-catenin/CBP interaction effect changes in 3D chromatin architecture and thereby gene expression and downstream effects remains to be elucidated. Here we perform Hi-C analyses on canonical and patient-derived pancreatic cancer cells before and after the treatment with ICG-001. In addition to global alteration of 3D chromatin domains, we unexpectedly identify insulin signaling genes enriched in the altered chromatin domains. We further demonstrate the chromatin loops associated with insulin signaling genes are significantly weakened after ICG-001 treatment. We finally elicit the deletion of a looping of IRS1, a key insulin signaling gene, significantly impede pancreatic cancer cell growth, indicating that looping-mediated insulin signaling might act as an oncogenic pathway to promote pancreatic cancer progression. Our work shows that targeting aberrant insulin chromatin looping in pancreatic cancer might provide a therapeutic benefit.

Enhancers dysfunction in the 3D genome of cancer cells

Eukaryotic genomes are spatially organized inside the cell nucleus, forming a threedimensional (3D) architecture that allows for spatial separation of nuclear processes and for controlled expression of genes required for cell identity specification and tissue homeostasis. Hence, it is of no surprise that mis-regulation of genome architecture through rearrangements of the linear genome sequence or epigenetic perturbations are often linked to aberrant gene expression programs in tumor cells. Increasing research efforts have shed light into the causes and consequences of alterations of 3D genome organization. In this review, we summarize the current knowledge on how 3D genome architecture is dysregulated in cancer, with a focus on enhancer highjacking events and their contribution to tumorigenesis. Studying the functional effects of genome architecture perturbations on gene expression in cancer offers a unique opportunity for a deeper understanding of tumor biology and sets the basis for the discovery of novel therapeutic targets.

Three-dimensional genome structure and function

Linear DNA undergoes a series of compression and folding events, forming various three-dimensional (3D) structural units in mammalian cells, including chromosomal territory, compartment, topologically associating domain, and chromatin loop. These structures play crucial roles in regulating gene expression, cell differentiation, and disease progression. Deciphering the principles underlying 3D genome folding and the molecular mechanisms governing cell fate determination remains a challenge. With advancements in high-throughput sequencing and imaging techniques, the hierarchical organization and functional roles of higher-order chromatin structures have been gradually illuminated. This review systematically discussed the structural hierarchy of the 3D genome, the effects and mechanisms of cis-regulatory elements interaction in the 3D genome for regulating spatiotemporally specific gene expression, the roles and mechanisms of dynamic changes in 3D chromatin conformation during embryonic development, and the pathological mechanisms of diseases such as congenital developmental abnormalities and cancer, which are attributed to alterations in 3D genome organization and aberrations in key structural proteins. Finally, prospects were made for the research about 3D genome structure, function, and genetic intervention, and the roles in disease development, prevention, and treatment, which may offer some clues for precise diagnosis and treatment of related diseases.

Epigenetic reprogramming-induced guanidinoacetic acid synthesis promotes pancreatic cancer metastasis and transcription-activating histone modifications

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) tends to undergo distant metastasis, especially liver metastasis, leading to a poor prognosis. Metabolic remodelling and epigenetic reprogramming are two important hallmarks of malignant tumours and participate in regulating PDAC tumorigenesis and metastasis. However, the interaction between these two processes during PDAC metastasis has not been fully elucidated.

METHODS

We performed metabolomics analysis to identify the critical metabolites associated with PDAC liver metastasis and focused on guanidinoacetic acid (GAA). Intracellular GAA content was significantly increased in liver metastatic PDAC cells compared to primary cancer cells in mouse xenograft tumour models. The effects of GAA supplementation and glycine amidinotransferase (GATM) knockdown on PDAC metastasis were assessed by analysing cell migration, filopodia formation, epithelial-mesenchymal transition (EMT), and in vivo metastasis in different cell and animal models. Next, ChIP‒qPCR, 3C‒qPCR, and CRISPRi/dCas9-KRAB experiments were used to validate the "epigenome-metabolome" mechanism. Finally, the results of in vitro approaches, including RNA-seq, CUT&RUN, RT‒qPCR, and western blot analyses, as well as luciferase reporter gene assay and transwell assay, revealed the GAA-c-Myc-HMGA axis and transcription-activating histone modifications reprogramming.

RESULTS

A high level of intracellular GAA was associated with PDAC liver metastasis. GAA could promote the migration, EMT, and liver metastasis of pancreatic cancer cells in vitro and in vivo. Next, we explored the role of GATM-mediated de novo GAA synthesis in pancreatic cancer metastasis. High expression of GATM was positively correlated with advanced N stage in PDAC. Knockdown of GATM significantly reduced the intracellular level of GAA, suppressed EMT, and inhibited PDAC liver metastasis, and these effects were attenuated by GAA supplementation. Mechanistically, we identified the active enhancers looped to the Gatm gene locus that promoted GATM expression and PDAC liver metastasis. Furthermore, we found that GAA promoted cell migration and EMT by regulating c-Myc-mediated high mobility group AT-hook protein expression. Moreover, GAA increased the H3K4me3 modification level by upregulating histone methyltransferases, which induced the transcription of metastasis-related genes, including Myc.

CONCLUSIONS

These findings revealed the critical role of the epigenome-metabolome interaction in regulating PDAC liver metastasis and suggested potential therapeutic strategies targeting GAA metabolism and epigenetic regulatory mechanisms.

Comprehensive multiomics analyses reveal pervasive involvement of aberrant cohesin binding in transcriptional and chromosomal disorder of cancer cells

Chromatin organization, whose malfunction causes various diseases including cancer, is fundamentally controlled by cohesin. While cancer cells have been found with mutated or misexpressed cohesin genes, there is no comprehensive survey about the presence and role of abnormal cohesin binding in cancer cells. Here, we systematically identified ∼1% of cohesin-binding sites (701-2,633) as cancer-aberrant binding sites of cohesin (CASs). We integrated CASs with large-scale transcriptomics, epigenomics, 3D genomics, and clinical information. CASs represent tissue-specific epigenomic signatures enriched for cancer-dysregulated genes with functional and clinical significance. CASs exhibited alterations in chromatin compartments, loops within topologically associated domains, and cis-regulatory elements, indicating that CASs induce dysregulated genes through misguided chromatin structure. Cohesin depletion data suggested that cohesin binding at CASs actively regulates cancer-dysregulated genes. Overall, our comprehensive investigation suggests that aberrant cohesin binding is an essential epigenomic signature responsible for dysregulated chromatin structure and transcription in cancer cells.

High-resolution Hi-C maps highlight multiscale chromatin architecture reorganization during cold stress in Brachypodium distachyon

BACKGROUND

The adaptation of plants to cold stress involves changes in gene expression profiles that are associated with epigenetic regulation. Although the three-dimensional (3D) genome architecture is considered an important epigenetic regulator, the role of 3D genome organization in the cold stress response remains unclear.

RESULTS

In this study, we developed high-resolution 3D genomic maps using control and cold-treated leaf tissue of the model plant Brachypodium distachyon using Hi-C to determine how cold stress affects the 3D genome architecture. We generated ~ 1.5 kb resolution chromatin interaction maps and showed that cold stress disrupts different levels of chromosome organization, including A/B compartment transition, a reduction in chromatin compartmentalization and the size of topologically associating domains (TADs), and loss of long-range chromatin loops. Integrating RNA-seq information, we identified cold-response genes and revealed that transcription was largely unaffected by the A/B compartment transition. The cold-response genes were predominantly localized in compartment A. In contrast, transcriptional changes are required for TAD reorganization. We demonstrated that dynamic TAD events were associated with H3K27me3 and H3K27ac state alterations. Moreover, a loss of chromatin looping, rather than a gain of looping, coincides with alterations in gene expression, indicating that chromatin loop disruption may play a more important role than loop formation in the cold-stress response.

CONCLUSIONS

Our study highlights the multiscale 3D genome reprogramming that occurs during cold stress and expands our knowledge of the mechanisms underlying transcriptional regulation in response to cold stress in plants.

Dynamic chromatin architectures provide insights into the genetics of cattle myogenesis

BACKGROUND

Sharply increased beef consumption is propelling the genetic improvement projects of beef cattle in China. Three-dimensional genome structure is confirmed to be an important layer of transcription regulation. Although genome-wide interaction data of several livestock species have already been produced, the genome structure states and its regulatory rules in cattle muscle are still limited.

RESULTS

Here we present the first 3D genome data in Longissimus dorsi muscle of fetal and adult cattle (Bos taurus). We showed that compartments, topologically associating domains (TADs), and loop undergo re-organization and the structure dynamics were consistent with transcriptomic divergence during muscle development. Furthermore, we annotated cis-regulatory elements in cattle genome during myogenesis and demonstrated the enrichments of promoter and enhancer in selection sweeps. We further validated the regulatory function of one HMGA2 intronic enhancer near a strong sweep region on primary bovine myoblast proliferation.

CONCLUSIONS

Our data provide key insights of the regulatory function of high order chromatin structure and cattle myogenic biology, which will benefit the progress of genetic improvement of beef cattle.

HNF1B-driven three-dimensional chromatin structure for molecular classification in pancreatic cancers

The molecular subtypes of pancreatic cancer (PC), either classical/progenitor-like or basal/squamous-like, are currently a major topic of research because of their direct association with clinical outcomes. Some transcription factors (TFs) have been reported to be associated with these subtypes. However, the mechanisms by which these molecular signatures of PCs are established remain unknown. Epigenetic regulatory processes, supported by dynamic changes in the chromatin structure, are essential for transcriptional profiles. Previously, we reported the importance of open chromatin profiles in the biological features and transcriptional status of PCs. Here, we aimed to analyze the relationships between three-dimensional (3D) genome structures and the molecular subtypes of human PCs using Hi-C analysis. We observed a correlation of the specific elements of 3D genome modules, including compartments, topologically associating domains, and enhancer-promoter loops, with the expression of related genes. We focused on HNF1B, a TF that is implicated in the progenitor subtype. Forced expression of HNF1B in squamous-type PC organoids induced the upregulation and downregulation of genes associated with progenitor and squamous subtypes, respectively. Long-range genomic interactions induced by HNF1B were accompanied by compartment modulation and H3K27ac redistribution. We also found that these HNF1B-induced changes in subtype-related gene expression required an intrinsically disordered region, suggesting a possible involvement of phase separation in compartment modulation. Thus, mapping of 3D structural changes induced by TFs, such as HNF1B, may become a useful resource for further understanding the molecular features of PCs.

From Samples to Germline and Somatic Sequence Variation: A Focus on Next-Generation Sequencing in Melanoma Research

Next-generation sequencing (NGS) applications have flourished in the last decade, permitting the identification of cancer driver genes and profoundly expanding the possibilities of genomic studies of cancer, including melanoma. Here we aimed to present a technical review across many of the methodological approaches brought by the use of NGS applications with a focus on assessing germline and somatic sequence variation. We provide cautionary notes and discuss key technical details involved in library preparation, the most common problems with the samples, and guidance to circumvent them. We also provide an overview of the sequence-based methods for cancer genomics, exposing the pros and cons of targeted sequencing vs. exome or whole-genome sequencing (WGS), the fundamentals of the most common commercial platforms, and a comparison of throughputs and key applications. Details of the steps and the main software involved in the bioinformatics processing of the sequencing results, from preprocessing to variant prioritization and filtering, are also provided in the context of the full spectrum of genetic variation (SNVs, indels, CNVs, structural variation, and gene fusions). Finally, we put the emphasis on selected bioinformatic pipelines behind (a) short-read WGS identification of small germline and somatic variants, (b) detection of gene fusions from transcriptomes, and (c) de novo assembly of genomes from long-read WGS data. Overall, we provide comprehensive guidance across the main methodological procedures involved in obtaining sequencing results for the most common short- and long-read NGS platforms, highlighting key applications in melanoma research.

Multilevel view on chromatin architecture alterations in cancer

Chromosomes inside the nucleus are not located in the form of linear molecules. Instead, there is a complex multilevel genome folding that includes nucleosomes packaging, formation of chromatin loops, domains, compartments, and finally, chromosomal territories. Proper spatial organization play an essential role for the correct functioning of the genome, and is therefore dynamically changed during development or disease. Here we discuss how the organization of the cancer cell genome differs from the healthy genome at various levels. A better understanding of how malignization affects genome organization and long-range gene regulation will help to reveal the molecular mechanisms underlying cancer development and evolution.

A cis-regulatory lexicon of DNA motif combinations mediating cell-type-specific gene regulation

Gene expression is controlled by transcription factors (TFs) that bind cognate DNA motif sequences in cis-regulatory elements (CREs). The combinations of DNA motifs acting within homeostasis and disease, however, are unclear. Gene expression, chromatin accessibility, TF footprinting, and H3K27ac-dependent DNA looping data were generated and a random-forest-based model was applied to identify 7,531 cell-type-specific cis-regulatory modules (CRMs) across 15 diploid human cell types. A co-enrichment framework within CRMs nominated 838 cell-type-specific, recurrent heterotypic DNA motif combinations (DMCs), which were functionally validated using massively parallel reporter assays. Cancer cells engaged DMCs linked to neoplasia-enabling processes operative in normal cells while also activating new DMCs only seen in the neoplastic state. This integrative approach identifies cell-type-specific cis-regulatory combinatorial DNA motifs in diverse normal and diseased human cells and represents a general framework for deciphering cis-regulatory sequence logic in gene regulation.

Downregulation of TUSC3 promotes EMT and hepatocellular carcinoma progression through LIPC/AKT axis

Background

Hepatocellular carcinoma (HCC) is one of the most common and malignant tumors in the digestive tract. Tumor Suppressor Candidate 3 (TUSC3) is one subunit of the endoplasmic reticulum Oligosaccharyl transferase (OST) complex, which plays an important role in N-glycosylation during the protein folding process. However, the role of TUSC3 in the initiation and progression of HCC has not been mentioned yet. In the present study, we aim to investigate the effects of TUSC3 on the initiation and progression of HCC.

Methods

Immunohistochemical assay and qRT-PCR were used to detect the expression of TUSC3 and lipase C hepatic type (LIPC) in HCC tissue and cells. Loss-of-function and gain-of-function were applied to detect the function of TUSC3 and LIPC in vivo and in vitro. Immunofluorescence assay and co-immunoprecipitation were used to detect the relationship between TUSC3 and LPC. Western blot was applied to detect the expression of epithelial–mesenchymal transition (EMT) markers and the Akt signaling pathway.

Results

TUSC3 was aberrantly decreased in hepatocellular carcinoma tissues compared to the matched adjacent normal tissues, which resulted in bigger size of tumor (P = 0.001, Table 2), worse differentiation (P = 0.006, Table 2) and an advanced BCLC stage. Down-regulation of TUSC3 led to the enhanced proliferation and migration of hepatocellular carcinoma cells in vivo and vitro, whereas the opposite effect could be observed in the TUSC3-overexpression group. The analysis of TUSC3 microarray showed that LIPC, a glycoprotein primarily synthesized and secreted by hepatocytes, was a downstream target of TUSC3, and it negatively modulated the development of HCC. The morphological changes in HCC cells indicated that TUSC3 regulated the epithelial-mesenchymal transition (EMT). Mechanistically, TUSC3 inhibited EMT progression through the LIPC/AKT axis.

Conclusion

Down-regulation of TUSC3 promotes EMT progression by activating AKT signaling via targeting LIPC in HCC, which is probably the possible mechanism driving TUSC3-deficient hepatocellular carcinoma cells toward a malignant phenotype.

3D chromatin remodeling potentiates transcriptional programs driving cell invasion

The contribution of deregulated chromatin architecture, including topologically associated domains (TADs), to cancer progression remains ambiguous. CCCTC-binding factor (CTCF) is a central regulator of higher-order chromatin structure that undergoes copy number loss in over half of all breast cancers, but the impact of this defect on epigenetic programming and chromatin architecture remains unclear. We find that under physiological conditions, CTCF organizes subTADs to limit the expression of oncogenic pathways, including phosphatidylinositol 3-kinase (PI3K) and cell adhesion networks. Loss of a single CTCF allele potentiates cell invasion through compromised chromatin insulation and a reorganization of chromatin architecture and histone programming that facilitates de novo promoter-enhancer contacts. However, this change in the higher-order chromatin landscape leads to a vulnerability to inhibitors of mTOR. These data support a model whereby subTAD reorganization drives both modification of histones at de novo enhancer-promoter contacts and transcriptional up-regulation of oncogenic transcriptional networks.

Construction of immune-related signature and identification of S100A14 determining immune-suppressive microenvironment in pancreatic cancer

Pancreatic cancer (PC) is a highly lethal and aggressive disease with its incidence and mortality quite discouraging. A robust prognostic signature and novel biomarkers are urgently needed for accurate stratification of the patients and optimization of clinical decision-making. Since the critical role of immune microenvironment in the progression of PC, a prognostic signature based on seven immune-related genes was established, which was validated in The Cancer Genome Atlas (TCGA) training set, TCGA testing set, TCGA entire set and GSE71729 set. Furthermore, S100A14 (S100 Calcium Binding Protein A14) was identified as the gene occupying the most paramount position in risk signature. According to the GSEA, CIBERSORT and ESTIMATE algorithm, S100A14 was mainly associated with lower proportion of CD8 + T cells and higher proportion of M0 macrophages in PC tissue. Meanwhile, analysis of single-cell dataset CRA001160 revealed a significant negative correlation between S100A14 expression in PC cells and CD8 + T cell infiltration, which was further confirmed by tissue microenvironment landscape imaging and machine learning-based analysis in our own PUMCH cohort. Additionally, analysis of a pan-pancreatic cancer cell line illustrated that S100A14 might inhibit CD8 + T cell activation via the upregulation of PD-L1 expression in PC cells, which was also verified by the immunohistochemical results of PUMCH cohort. Finally, tumor mutation burden analysis and immunophenoscore algorithm revealed that patients with high S100A14 expression had a higher probability of responding to immunotherapy. In conclusion, our study established an efficient immune-related prediction model and identified the potential role of S100A14 in regulating the immune microenvironment and serving as a biomarker for immunotherapy efficacy prediction.

Multifaceted regulation of enhancers in cancer

Enhancer is one kind of cis-elements regulating gene transcription, whose activity is tightly controlled by epigenetic enzymes and histone modifications. Active enhancers are classified into typical enhancers, super-enhancers and over-active enhancers, according to the enrichment and location of histone modifications. Epigenetic factors control the level of histone modifications on enhancers to determine their activity, such as histone methyltransferases and acetylases. Transcription factors, cofactors and mediators co-operate together and are required for enhancer functions. In turn, abnormalities in these trans-acting factors affect enhancer activity. Recent studies have revealed enhancer dysregulation as one of the important features for cancer. Variations in enhancer regions and mutations of enhancer regulatory genes are frequently observed in cancer cells, and altering the activity of onco-enhancers is able to repress oncogene expression, and suppress tumorigenesis and metastasis. Here we summarize the recent discoveries about enhancer regulation in cancer and discuss their potential application in diagnosis and treatment.