3D genome organization in the epithelial-mesenchymal transition spectrum

Inter3D: Capture of TAD Reorganization Endows Variant Patterns of Gene Transcription

Topologically associating domain (TAD) reorganization commonly occurs in the cell nucleus and contributes to gene activation and inhibition through the separation or fusion of adjacent TADs. However, functional genes impacted by TAD alteration and the underlying mechanism of TAD reorganization regulating gene transcription remain to be fully elucidated. Here, we first developed a novel approach termed Inter3D to specifically identify genes regulated by TAD reorganization. Our study revealed that the segregation of TADs led to the disruption of intrachromosomal looping at the myosin light chain 12B (MYL12B) locus, via the meticulous reorganization of TADs mediating epigenomic landscapes within tumor cells, thereby exhibiting a significant correlation with the down-regulation of its transcriptional activity. Conversely, the fusion of TADs facilitated intrachromosomal interactions, suggesting a potential association with the activation of cytochrome P450 family 27 subfamily B member 1 (CYP27B1). Our study provides comprehensive insight into the capture of TAD rearrangement-mediated gene loci and moves toward understanding the functional role of TAD reorganization in gene transcription. The Inter3D pipeline developed in this study is freely available at https://github.com/bm2-lab/inter3D and https://ngdc.cncb.ac.cn/biocode/tool/BT7399.

Deciphering Pre-existing and Induced 3D Genome Architecture Changes involved in Constricted Melanoma Migration

Metastatic cancer cells traverse constricted spaces that exert forces on their nucleus and the genomic contents within. Cancerous tumors are highly heterogeneous and not all cells within them can achieve such a feat. Here, we investigated what initial genome architecture characteristics favor the constricted migratory ability of cancer cells and which arise only after passage through multiple constrictions. We identified a cell surface protein (ITGB4) whose expression correlates with increased initial constricted migration ability in human melanoma A375 cells. Sorting out this subpopulation allowed us to identify cellular and nuclear features that pre-exist and favor migration, as well as alterations that only appear after cells have passed through constrictions. We identified specific genomic regions that experienced altered genome spatial compartment profiles only after constricted migration. Our study reveals 3D genome structure contributions to both selection and induction mechanisms of cell fate change during cancer metastasis.

Regulating epithelial-mesenchymal plasticity from 3D genome organization

Epithelial-mesenchymal transition (EMT) is a dynamic process enabling polarized epithelial cells to acquire mesenchymal features implicated in development and carcinoma progression. As our understanding evolves, it is clear the reversible execution of EMT arises from complex epigenomic regulation involving histone modifications and 3-dimensional (3D) genome structural changes, leading to a cascade of transcriptional events. This review summarizes current knowledge on chromatin organization in EMT, with a focus on hierarchical structures of the 3D genome and chromatin accessibility changes.

A partial epithelial-mesenchymal transition signature for highly aggressive colorectal cancer cells that survive under nutrient restriction

Partial epithelial-mesenchymal transition (p-EMT) has recently been identified as a hybrid state consisting of cells with both epithelial and mesenchymal characteristics and is associated with the migration, metastasis, and chemoresistance of cancer cells. Here, we describe the induction of p-EMT in starved colorectal cancer (CRC) cells and identify a p-EMT gene signature that can predict prognosis. Functional characterisation of starvation-induced p-EMT in HCT116, DLD1, and HT29 cells showed changes in proliferation, morphology, and drug sensitivity, supported by in vivo studies using the chorioallantoic membrane model. An EMT-specific quantitative polymerase chain reaction (qPCR) array was used to screen for deregulated genes, leading to the establishment of an in silico gene signature that was correlated with poor disease-free survival in CRC patients along with the CRC consensus molecular subtype CMS4. Among the significantly deregulated p-EMT genes, a triple-gene signature consisting of SERPINE1, SOX10, and epidermal growth factor receptor (EGFR) was identified. Starvation-induced p-EMT was characterised by increased migratory potential and chemoresistance, as well as E-cadherin processing and internalisation. Both gene signature and E-cadherin alterations could be reversed by the proteasomal inhibitor MG132. Spatially resolving EGFR expression with high-resolution immunofluorescence imaging identified a proliferation stop in starved CRC cells caused by EGFR internalisation. In conclusion, we have gained insight into a previously undiscovered EMT mechanism that may become relevant when tumour cells are under nutrient stress, as seen in early stages of metastasis. Targeting this process of tumour cell dissemination might help to prevent EMT and overcome drug resistance. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

M6A reader YTHDF1 promotes malignant progression of laryngeal squamous carcinoma through activating the EMT pathway by EIF4A3

Laryngeal squamous cell carcinoma (LSCC) is one of the common malignant tumors in the head and neck region, and its high migration and invasion seriously threaten the survival and health of patients. In cancer development, m6A RNA modification plays a crucial role in regulating gene expression and signaling. This study delved into the function and mechanism of the m6A reading protein YTHDF1 in LSCC. It was found that YTHDF1 was highly expressed in the GEO database and LSCC tissues. Cell function experiments confirmed that the downregulation of YTHDF1 significantly inhibited the proliferation, migration, and invasion ability of LSCC cells. Further studies revealed that EIF4A3 was a downstream target gene of YTHDF1, and knockdown of EIF4A3 similarly significantly inhibited the malignant progression of LSCC in both in vivo and in vitro experiments. The molecular mechanism studies suggested that YTHDF1-EIF4A3 may promote the malignant development of LSCC by activating the EMT signaling pathway. This study provides important clues for an in-depth understanding of the pathogenesis of LSCC and is a solid foundation for the discovery of new therapeutic targets and approaches.

DiffDomain enables identification of structurally reorganized topologically associating domains

Topologically associating domains (TADs) are critical structural units in three-dimensional genome organization of mammalian genome. Dynamic reorganizations of TADs between health and disease states are associated with essential genome functions. However, computational methods for identifying reorganized TADs are still in the early stages of development. Here, we present DiffDomain, an algorithm leveraging high-dimensional random matrix theory to identify structurally reorganized TADs using high-throughput chromosome conformation capture (Hi-C) contact maps. Method comparison using multiple real Hi-C datasets reveals that DiffDomain outperforms alternative methods for false positive rates, true positive rates, and identifying a new subtype of reorganized TADs. Applying DiffDomain to Hi-C data from different cell types and disease states demonstrates its biological relevance. Identified reorganized TADs are associated with structural variations and epigenomic changes such as changes in CTCF binding sites. By applying to a single-cell Hi-C data from mouse neuronal development, DiffDomain can identify reorganized TADs between cell types with reasonable reproducibility using pseudo-bulk Hi-C data from as few as 100 cells per condition. Moreover, DiffDomain reveals differential cell-to-population variability and heterogeneous cell-to-cell variability in TADs. Therefore, DiffDomain is a statistically sound method for better comparative analysis of TADs using both Hi-C and single-cell Hi-C data.

TGF-β, EMT, and resistance to anti-cancer treatment

Transforming growth factor-β (TGF-β) signaling regulates cell-specific programs involved in embryonic development, wound-healing, and immune homeostasis. Yet, during tumor progression, these TGF-β-mediated programs are altered, leading to epithelial cell plasticity and a reprogramming of epithelial cells into mesenchymal lineages through epithelial-to-mesenchymal transition (EMT), a critical developmental program in morphogenesis and organogenesis. These changes, in turn, lead to enhanced carcinoma cell invasion, metastasis, immune cell differentiation, immune evasion, and chemotherapy resistance. Here, we discuss EMT as one of the critical programs associated with carcinoma cell plasticity and the influence exerted by TGF-β on carcinoma status and function. We further explore the composition of carcinoma and other cell populations within the tumor microenvironment, and consider the relevant outcomes related to the programs associated with cancer treatment resistance.

Role of epithelial-mesenchymal transition factor SNAI1 and its targets in ovarian cancer aggressiveness

Ovarian cancer remains the most lethal gynecologic malignancy in the USA. For over twenty years, epithelial-mesenchymal transition (EMT) has been characterized extensively in development and disease. The dysregulation of this process in cancer has been identified as a mechanism by which epithelial tumors become more aggressive, allowing them to survive and invade distant tissues. This occurs in part due to the increased expression of the EMT transcription factor, SNAI1 (Snail). In the case of epithelial ovarian cancer, Snail has been shown to contribute to cancer invasion, stemness, chemoresistance, and metabolic changes. Thus, in this review, we focus on summarizing current findings on the role of EMT (specifically, factors downstream of Snail) in determining ovarian cancer aggressiveness.

3D genome alterations and editing in pathology

The human genome is folded into a multi-level 3D structure that controls many nuclear functions including gene expression. Recently, alterations in 3D genome organization were associated with several genetic diseases and cancer. As a consequence, experimental approaches are now being developed to modify the global 3D genome organization and that of specific loci. Here, we discuss emerging experimental approaches of 3D genome editing that may prove useful in biomedicine.

AP-2δ Is the Most Relevant Target of AP-2 Family-Focused Cancer Therapy and Affects Genome Organization

Formerly hailed as "undruggable" proteins, transcription factors (TFs) are now under investigation for targeted therapy. In cancer, this may alter, inter alia, immune evasion or replicative immortality, which are implicated in genome organization, a process that accompanies multi-step tumorigenesis and which frequently develops in a non-random manner. Still, targeting-related research on some TFs is scarce, e.g., among AP-2 proteins, which are known for their altered functionality in cancer and prognostic importance. Using public repositories, bioinformatics tools, and RNA-seq data, the present study examined the ligandability of all AP-2 members, selecting the best one, which was investigated in terms of mutations, targets, co-activators, correlated genes, and impact on genome organization. AP-2 proteins were found to have the conserved "TF_AP-2" domain, but manifested different binding characteristics and evolution. Among them, AP-2δ has not only the highest number of post-translational modifications and extended strands but also contains a specific histidine-rich region and cleft that can receive a ligand. Uterine, colon, lung, and stomach tumors are most susceptible to AP-2δ mutations, which also co-depend with cancer hallmark genes and drug targets. Considering AP-2δ targets, some of them were located proximally in the spatial genome or served as co-factors of the genes regulated by AP-2δ. Correlation and functional analyses suggested that AP-2δ affects various processes, including genome organization, via its targets; this has been eventually verified in lung adenocarcinoma using expression and immunohistochemistry data of chromosomal conformation-related genes. In conclusion, AP-2δ affects chromosomal conformation and is the most appropriate target for cancer therapy focused on the AP-2 family.

Novel biological insights revealed from the investigation of multiscale genome architecture

Gene expression and cell fate determination require precise and coordinated epigenetic regulation. The complex three-dimensional (3D) genome organization plays a critical role in transcription in myriad biological processes. A wide range of architectural features of the 3D genome, including chromatin loops, topologically associated domains (TADs), chromatin compartments, and phase separation, together regulate the chromatin state and transcriptional activity at multiple levels. With the help of 3D genome informatics, recent biochemistry and imaging approaches based on different strategies have revealed functional interactions among biomacromolecules, even at the single-cell level. Here, we review the occurrence, mechanistic basis, and functional implications of dynamic genome organization, and outline recent experimental and computational approaches for profiling multiscale genome architecture to provide robust tools for studying the 3D genome.

Every gene everywhere all at once: High-precision measurement of 3D chromosome architecture with single-cell Hi-C

The three-dimensional (3D) structure of chromosomes influences essential biological processes such as gene expression, genome replication, and DNA damage repair and has been implicated in many developmental and degenerative diseases. In the past two centuries, two complementary genres of technology-microscopy, such as fluorescence in situ hybridization (FISH), and biochemistry, such as chromosome conformation capture (3C or Hi-C)-have revealed general principles of chromosome folding in the cell nucleus. However, the extraordinary complexity and cell-to-cell variability of the chromosome structure necessitate new tools with genome-wide coverage and single-cell precision. In the past decade, single-cell Hi-C emerges as a new approach that builds upon yet conceptually differs from bulk Hi-C assays. Instead of measuring population-averaged statistical properties of chromosome folding, single-cell Hi-C works as a proximity-based "biochemical microscope" that measures actual 3D structures of individual genomes, revealing features hidden in bulk Hi-C such as radial organization, multi-way interactions, and chromosome intermingling. Single-cell Hi-C has been used to study highly dynamic processes such as the cell cycle, cell-type-specific chromosome architecture ("structure types"), and structure-expression interplay, deepening our understanding of DNA organization and function.

3D genome organization in the epithelial-mesenchymal transition spectrum

Background

The plasticity along the epithelial-mesenchymal transition (EMT) spectrum has been shown to be regulated by various epigenetic repertoires. Emerging evidence of local chromatin conformation changes suggests that regulation of EMT may occur at a higher order of three-dimensional genome level.

Results

We perform Hi-C analysis and combine ChIP-seq data across cancer cell lines representing different EMT states. We demonstrate that the epithelial and mesenchymal genes are regulated distinctively. We find that EMT genes are regulated within their topologically associated domains (TADs), with only a subset of mesenchymal genes being influenced by A/B compartment switches, indicating topological remodeling is required in the transcriptional regulation of these genes. At the TAD level, epithelial and mesenchymal genes are associated with different regulatory trajectories. The epithelial gene-residing TADs are enriched with H3K27me3 marks in the mesenchymal-like states. The mesenchymal gene-residing TADs, which do not show enrichment of H3K27me3 in epithelial-like states, exhibit increased interaction frequencies with regulatory elements in the mesenchymal-like states.

Conclusions

We propose a novel workflow coupling immunofluorescence and dielectrophoresis to unravel EMT heterogeneity at single-cell resolution. The predicted three-dimensional structures of chromosome 10, harboring Vimentin, identify cell clusters of different states. Our results pioneer a novel avenue to decipher the complexities underlying the regulation of EMT and may infer the barriers of plasticity in the 3D genome context.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02687-x.