Three-dimensional disorganization of the cancer genome occurs coincident with long-range genetic and epigenetic alterations

CDK9 inhibition inhibits multiple oncogenic transcriptional and epigenetic pathways in prostate cancer

BACKGROUND

Cyclin-dependent kinase 9 (CDK9) stimulates oncogenic transcriptional pathways in cancer and CDK9 inhibitors have emerged as promising therapeutic candidates.

METHODS

The activity of an orally bioavailable CDK9 inhibitor, CDKI-73, was evaluated in prostate cancer cell lines, a xenograft mouse model, and patient-derived tumor explants and organoids. Expression of CDK9 was evaluated in clinical specimens by mining public datasets and immunohistochemistry. Effects of CDKI-73 on prostate cancer cells were determined by cell-based assays, molecular profiling and transcriptomic/epigenomic approaches.

RESULTS

CDKI-73 inhibited proliferation and enhanced cell death in diverse in vitro and in vivo models of androgen receptor (AR)-driven and AR-independent models. Mechanistically, CDKI-73-mediated inhibition of RNA polymerase II serine 2 phosphorylation resulted in reduced expression of BCL-2 anti-apoptotic factors and transcriptional defects. Transcriptomic and epigenomic approaches revealed that CDKI-73 suppressed signaling pathways regulated by AR, MYC, and BRD4, key drivers of dysregulated transcription in prostate cancer, and reprogrammed cancer-associated super-enhancers. These latter findings prompted the evaluation of CDKI-73 with the BRD4 inhibitor AZD5153, a combination that was synergistic in patient-derived organoids and in vivo.

CONCLUSION

Our work demonstrates that CDK9 inhibition disrupts multiple oncogenic pathways and positions CDKI-73 as a promising therapeutic agent for prostate cancer, particularly aggressive, therapy-resistant subtypes.

The effect of trisomic chromosomes on spatial genome organization and global transcription in embryonic stem cells

Aneuploidy frequently occurs in cancer and developmental diseases such as Down syndrome, with its functional consequences implicated in dosage effects on gene expression and global perturbation of stress response and cell proliferation pathways. However, how aneuploidy affects spatial genome organization remains less understood. In this study, we addressed this question by utilizing the previously established isogenic wild-type (WT) and trisomic mouse embryonic stem cells (mESCs). We employed a combination of Hi-C, RNA-seq, chromosome painting and nascent RNA imaging technologies to compare the spatial genome structures and gene transcription among these cells. We found that trisomy has little effect on spatial genome organization at the level of A/B compartment or topologically associating domain (TAD). Inter-chromosomal interactions are associated with chromosome regions with high gene density, active histone modifications and high transcription levels, which are confirmed by imaging. Imaging also revealed contracted chromosome volume and weakened transcriptional activity for trisomic chromosomes, suggesting potential implications for the transcriptional output of these chromosomes. Our data resources and findings may contribute to a better understanding of the consequences of aneuploidy from the angle of spatial genome organization.

Evaluation of the gene fusion landscape in early onset sporadic rectal cancer reveals association with chromatin architecture and genome stability

Gene fusions represent a distinct class of structural variants identified frequently in cancer genomes across cancer types. Several gene fusions exhibit gain of oncogenic function and thus have been the focus of development of efficient targeted therapies. However, investigation of fusion landscape in early-onset sporadic rectal cancer, a poorly studied colorectal cancer subtype prevalent in developing countries, has not been performed. Here, we present a comprehensive landscape of gene fusions in EOSRC and CRC using patient derived tumor samples and data from The Cancer Genome Atlas, respectively. Gene Ontology analysis revealed enrichment of unique biological process terms associated with 5'- and 3'- fusion partner genes. Extensive network analysis highlighted genes exhibiting significant promiscuity in fusion formation and their association with chromosome fragile sites. Investigation of fusion formation in the context of global chromatin architecture unraveled a novel mode of gene activation that arose from fusion between genes located in orthogonal chromatin compartments. The study provides novel evidence linking fusions to genome stability and architecture and unearthed a hitherto unidentified mode of gene activation in cancer.

Estrogen-induced chromatin looping changes identify a subset of functional regulatory elements

Transcriptional enhancers can regulate individual or multiple genes through long-range three-dimensional (3D) genome interactions, and these interactions are commonly altered in cancer. Yet, the functional relationship between changes in 3D interactions associated with regulatory regions and differential gene expression appears context-dependent. In this study, we used HiChiP to capture changes in 3D genome interactions between active regulatory regions of endometrial cancer cells in response to estrogen treatment and uncovered significant differential long-range interactions that are strongly enriched for estrogen receptor α (ER) bound sites (ERBS). The ERBS anchoring differential loops with either a gene's promoter or distal regions were correlated with larger transcriptional responses to estrogen compared to ERBS not involved in differential interactions. To functionally test this observation, CRISPR-based Enhancer-i was used to deactivate specific ERBS, which revealed a wide range of effects on the transcriptional response to estrogen. However, these effects are only subtly and not significantly stronger for ERBS in differential loops. In addition, we observed an enrichment of 3D interactions between the promoters of estrogen up-regulated genes and found that looped promoters can work together cooperatively. Overall, our work suggests that changes in 3D genome structure upon estrogen treatment identify some functionally important regulatory regions; however, these changes aren't required for a transcriptional response to E2 in endometrial cancer cells.

Transcription factor shapes chromosomal conformation and regulates gene expression in bacterial adaptation

Genomic mutations allow bacteria to adapt rapidly to adverse stress environments. The three-dimensional conformation of the genome may also play an important role in transcriptional regulation and environmental adaptation. Here, using chromosome conformation capture, we investigate the high-order architecture of the Zymomonas mobilis chromosome in response to genomic mutation and ambient stimuli (acetic acid and furfural, derived from lignocellulosic hydrolysate). We find that genomic mutation only influences the local chromosome contacts, whereas stress of acetic acid and furfural restrict the long-range contacts and significantly change the chromosome organization at domain scales. Further deciphering the domain feature unveils the important transcription factors, Ferric uptake regulator (Fur) proteins, which act as nucleoid-associated proteins to promote long-range (>200 kb) chromosomal communications and regulate the expression of genes involved in stress response. Our work suggests that ubiquitous transcription factors in prokaryotes mediate chromosome organization and regulate stress-resistance genes in bacterial adaptation.

A comprehensive benchmarking with interpretation and operational guidance for the hierarchy of topologically associating domains

Topologically associating domains (TADs), megabase-scale features of chromatin spatial architecture, are organized in a domain-within-domain TAD hierarchy. Within TADs, the inner and smaller subTADs not only manifest cell-to-cell variability, but also precisely regulate transcription and differentiation. Although over 20 TAD callers are able to detect TAD, their usability in biomedicine is confined by a disagreement of outputs and a limit in understanding TAD hierarchy. We compare 13 computational tools across various conditions and develop a metric to evaluate the similarity of TAD hierarchy. Although outputs of TAD hierarchy at each level vary among callers, data resolutions, sequencing depths, and matrices normalization, they are more consistent when they have a higher similarity of larger TADs. We present comprehensive benchmarking of TAD hierarchy callers and operational guidance to researchers of life science researchers. Moreover, by simulating the mixing of different types of cells, we confirm that TAD hierarchy is generated not simply from stacking Hi-C heatmaps of heterogeneous cells. Finally, we propose an air conditioner model to decipher the role of TAD hierarchy in transcription.

Somatic structural variants drive distinct modes of oncogenesis in melanoma

The diversity of structural variants (SVs) in melanoma and how they impact oncogenesis are incompletely known. We performed harmonized analysis of SVs across melanoma histologic and genomic subtypes, and we identified distinct global properties between subtypes. These included the frequency and size of SVs and SV classes, their relation to chromothripsis events, and the impact on cancer-related genes of SVs that alter topologically associated domain (TAD) boundaries. Following our prior identification of double-stranded break repair deficiency in a subset of triple-wild-type cutaneous melanoma, we identified MRE11 and NBN loss-of-function SVs in melanomas with this mutational signature. Experimental knockouts of MRE11 and NBN, followed by olaparib cell viability assays in melanoma cells, indicated that dysregulation of each of these genes may cause sensitivity to PARP inhibitors in cutaneous melanomas. Broadly, harmonized analysis of melanoma SVs revealed distinct global genomic properties and molecular drivers, which may have biological and therapeutic impact.

Three-dimensional chromatin landscapes in somatotroph tumour

BACKGROUND

The three-dimensional (3D) genome architecture plays a critical role inregulating gene expression. However, the specific alterations in thisarchitecture within somatotroph tumors and their implications for gene expression remain largely unexplored.

METHODS

We employed Hi-C and RNA-seq analyses to compare the 3D genomic structures of somatotroph tumors with normal pituitary tissue. This comprehensive approachenabled the characterization of A/B compartments, topologically associateddomains (TADs), and chromatin loops, integrating these with gene expression patterns.

RESULTS

We observed a decrease in both the frequency of chromosomal interactions andthe size of TADs in tumor tissue compared to normal tissue. Conversely, the number of TADs and chromatin loops was found to be increased in tumors. Integrated analysis of Hi-C and RNA-seq data demonstrated that changes inhigher-order chromat in structure were associated with alterations in gene expression. Specifically, genes in A compartments showed higher density and increased expression relative to those in B compartments. Moreover, the weakand enhanced insulation boundaries were identified, and the associated genes were enriched in the Wnt/β-Catenin signaling pathway. We identified the gainedand lost loops in tumor and integrated these differences with transcriptional changes to examine the functional relevance of the identified loops. Notably, we observed an enhanced insulation boundary and a greater number of loops in the TCF7L2 gene region within tumors, which was accompanied by an upregulation of TCF7L2 expression. Subsequently, TCF7L2 expression was confirmed through qRT-PCR, and upregulated TCF7L2 prompted cell proliferation and growth hormone (GH) secretion in vitro.

CONCLUSION

Our results provide comprehensive 3D chromatin architecture maps of somatotroph tumors and offer a valuable resource for furthering the understanding of the underlying biology and mechanisms of gene expression regulation.

Chromatin activation with H3K36me2 and compartment shift in metastatic castration-resistant prostate cancer

Epigenetic modifiers are upregulated during the process of prostate cancer, acquiring resistance to castration therapy and becoming lethal metastatic castration-resistant prostate cancer (CRPC). However, the relationship between regulation of histone modifications and chromatin structure in CRPC has yet not fully been validated. Here, we reanalyzed publicly available clinical transcriptome and clinical outcome data and identified NSD2, a histone methyltransferase that catalyzes H3K36me2, as an epigenetic modifier that was upregulated in CRPC and whose increased expression in prostate cancer correlated with higher recurrence rate. We performed ChIP-seq, RNA-seq, and Hi-C to conduct comprehensive epigenomic and transcriptomic analyses to identify epigenetic reprogramming in CRPC. In regions where H3K36me2 was increased, H3K27me3 was decreased, and the compartment was shifted from inactive to active. In these regions, 68 aberrantly activated genes were identified as candidate downstream genes of NSD2 in CRPC. Among these genes, we identified KIF18A as critical for CRPC growth. Under NSD2 upregulation in CRPC, epigenetic alteration with H3K36me2-gain and H3K27me3-loss occurs accompanying with an inactive-to-active compartment shift, suggesting that histone modification and chromatin structure cooperatively change prostate carcinogenesis.

Emergence of enhancers at late DNA replicating regions

Enhancers are fast-evolving genomic sequences that control spatiotemporal gene expression patterns. By examining enhancer turnover across mammalian species and in multiple tissue types, we uncover a relationship between the emergence of enhancers and genome organization as a function of germline DNA replication time. While enhancers are most abundant in euchromatic regions, enhancers emerge almost twice as often in late compared to early germline replicating regions, independent of transposable elements. Using a deep learning sequence model, we demonstrate that new enhancers are enriched for mutations that alter transcription factor (TF) binding. Recently evolved enhancers appear to be mostly neutrally evolving and enriched in eQTLs. They also show more tissue specificity than conserved enhancers, and the TFs that bind to these elements, as inferred by binding sequences, also show increased tissue-specific gene expression. We find a similar relationship with DNA replication time in cancer, suggesting that these observations may be time-invariant principles of genome evolution. Our work underscores that genome organization has a profound impact in shaping mammalian gene regulation.

Machine and deep learning methods for predicting 3D genome organization

Three-Dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, Topologically Associating Domains (TADs), and A/B compartments play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers, Transcription Factor Binding Site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, TAD boundaries) and analyze their pros and cons. We also point out obstacles of computational prediction of 3D interactions and suggest future research directions.

Sterile inflammation via TRPM8 RNA-dependent TLR3-NF-kB/IRF3 activation promotes antitumor immunity in prostate cancer

Inflammation is a common condition of prostate tissue, whose impact on carcinogenesis is highly debated. Microbial colonization is a well-documented cause of a small percentage of prostatitis cases, but it remains unclear what underlies the majority of sterile inflammation reported. Here, androgen- independent fluctuations of PSA expression in prostate cells have lead us to identify a prominent function of the Transient Receptor Potential Cation Channel Subfamily M Member 8 (TRPM8) gene in sterile inflammation. Prostate cells secret TRPM8 RNA into extracellular vesicles (EVs), which primes TLR3/NF-kB-mediated inflammatory signaling after EV endocytosis by epithelial cancer cells. Furthermore, prostate cancer xenografts expressing a translation-defective form of TRPM8 RNA contain less collagen type I in the extracellular matrix, significantly more infiltrating NK cells, and larger necrotic areas as compared to control xenografts. These findings imply sustained, androgen-independent expression of TRPM8 constitutes as a promoter of anticancer innate immunity, which may constitute a clinically relevant condition affecting prostate cancer prognosis.

Epigenetic disruption of the RARγ complex impairs its function to bookmark AR enhancer interactions required for enzalutamide sensitivity in prostate cancer

The current study in prostate cancer (PCa) focused on the genomic mechanisms at the cross-roads of pro-differentiation signals and the emergence of lineage plasticity. We explored an understudied cistromic mechanism involving RARγ's ability to govern AR cistrome-transcriptome relationships, including those associated with more aggressive PCa features. The RARγ complex in PCa cell models was enriched for canonical cofactors, as well as proteins involved in RNA processing and bookmarking. Identifying the repertoire of miR-96 bound and regulated gene targets, including those recognition elements marked by m6A, revealed their significant enrichment in the RARγ complex. RARγ significantly enhanced the AR cistrome, particularly in active enhancers and super-enhancers, and overlapped with the binding of bookmarking factors. Furthermore, RARγ expression led to nucleosome-free chromatin enriched with H3K27ac, and significantly enhanced the AR cistrome in G2/M cells. RARγ functions also antagonized the transcriptional actions of the lineage master regulator ONECUT2. Similarly, gene programs regulated by either miR-96 or antagonized by RARγ were enriched in alternative lineages and more aggressive PCa phenotypes. Together these findings reveal an under-investigated role for RARγ, modulated by miR-96, to bookmark enhancer sites during mitosis. These sites are required by the AR to promote transcriptional competence, and emphasize luminal differentiation, while antagonizing ONECUT2.

Three-dimensional chromatin landscapes in hepatocellular carcinoma associated with hepatitis B virus

BACKGROUND

Three-dimensional (3D) chromatin architecture frequently altered in cancer. However, its changes during the pathogenesis of hepatocellular carcinoma (HCC) remained elusive.

METHODS

Hi-C and RNA-seq were applied to study the 3D chromatin landscapes and gene expression of HCC and ANHT. Hi-C Pro was used to generate genome-wide raw interaction matrices, which were normalized via iterative correction (ICE). Moreover, the chromosomes were divided into different compartments according to the first principal component (E1). Furthermore, topologically associated domains (TADs) were visualized via WashU Epigenome Browser. Furthermore, differential expression analysis of ANHT and HCC was performed using the DESeq2 R package. Additionally, dysregulated genes associated with 3D genome architecture altered were confirmed using TCGA, qRT-PCR, immunohistochemistry (IHC), etc. RESULTS: First, the intrachromosomal interactions of chr1, chr2, chr5, and chr11 were significantly different, and the interchromosomal interactions of chr4-chr10, chr13-chr21, chr15-chr22, and chr16-chr19 are remarkably different between ANHT and HCC, which resulted in the up-regulation of TP53I3 and ZNF738 and the down-regulation of APOC3 and APOA5 in HCC. Second, 49 compartment regions on 18 chromosomes have significantly switched (A-B or B-A) during HCC tumorigenesis, contributing to up-regulation of RAP2A. Finally, a tumor-specific TAD boundary located on chr5: 6271000-6478000 and enhancer hijacking were identified in HCC tissues, potentially associated with the elevated expression of MED10, whose expression were associated with poor prognosis of HCC patients.

CONCLUSION

This study demonstrates the crucial role of chromosomal structure variation in HCC oncogenesis and potential novel biomarkers of HCC, laying a foundation for cancer precision medicine development.

3D genome structural variations play important roles in regulating seed oil content of Brassica napus

Dissecting the complex regulatory mechanism of seed oil content (SOC) is one of the main research goals in Brassica napus. Increasing evidence suggests that genome architecture is linked to multiple biological functions. However, the effect of genome architecture on SOC regulation remains unclear. Here, we used high-throughput chromatin conformation capture to characterize differences in the three-dimensional (3D) landscape of genome architecture of seeds from two B. napus lines, N53-2 (with high SOC) and Ken-C8 (with low SOC). Bioinformatics analysis demonstrated that differentially accessible regions and differentially expressed genes between N53-2 and Ken-C8 were preferentially enriched in regions with quantitative trait loci (QTLs)/associated genomic regions (AGRs) for SOC. A multi-omics analysis demonstrated that expression of SOC-related genes was tightly correlated with genome structural variations in QTLs/AGRs of B. napus. The candidate gene BnaA09g48250D, which showed structural variation in a QTL/AGR on chrA09, was identified by fine-mapping of a KN double-haploid population derived from hybridization of N53-2 and Ken-C8. Overexpression and knockout of BnaA09g48250D led to significant increases and decreases in SOC, respectively, in the transgenic lines. Taken together, our results reveal the 3D genome architecture of B. napus seeds and the roles of genome structural variations in SOC regulation, enriching our understanding of the molecular mechanisms of SOC regulation from the perspective of spatial chromatin structure.

Rethinking cancer initiation: The role of large-scale mutational events

Cancer initiation is revisited in light of recent discoveries in cancer pathogenesis. Of note is the detection of mutated cancer genes in benign conditions. More significantly, somatic clones, which harbor mutations in cancer genes, arise in normal tissues from early development through adulthood, but seldom do they transform into cancer. Further, clustered mutational events-kataegis, chromothripsis and chromoplexy-are widespread in cancer, generating point mutations and chromosomal rearrangements in a single cellular catastrophe. These observations are contrary to the prevailing somatic mutation theory, which states that a cancer is caused by the gradual accumulation of mutations over time. A different perspective is proposed within the framework of Waddington's epigenetic landscape wherein tumorigenesis is viewed primarily as a disruption of cell development. Cell types are defined by their specific gene-expression profiles, determined by the gene regulatory network, and can be regarded as attractor states of the network dynamics: they represent specific, self-stabilizing patterns of gene activities across the genome. However, large-scale mutational events reshape the landscape topology, creating abnormal "unphysiological" attractors. This is the crux of the process of initiation. Initiation primes the cell for conversion into a tumor phenotype by oncogenes and tumor suppressor genes, which drive cell proliferation and clonal diversification. This view of tumorigenesis calls for a different approach to therapy.

The MYC axis in advanced prostate cancer is impacted through concurrent targeting of ERβ and AR using a novel ERβ-selective ligand alongside Enzalutamide

We have dissected the role of Estrogen receptor beta (ERβ) in prostate cancer (PCa) with a novel ERβ ligand, OSU-ERb-12. Drug screens revealed additive interactions between OSU-ERB-12 and either epigenetic inhibitors or the androgen receptor antagonist, Enzalutamide (Enza). Clonogenic and cell biolody studies supported the potent additive effects of OSU-ERB-12 (100nM) and Enza (1μM). The cooperative behavior was in PCa cell lines treated with either OSU-ERB-12 plus Enza or combinations involving 17β-estradiol (E2). OSU-ERb-12 plus Enza uniquely impacted the transcriptiome, accessible chromatin, and the AR, MYC and H3K27ac cistromes. This included skewed transcriptional responses including suppression of the androgen and MYC transcriptomes, and repressed MYC protein. OSU-ERb-12 plus Enza uniquely impacted chromatin accessibility at approximately 3000 nucleosome-free sites, enriched at enhancers, enriched for basic Helix-Loop-Helix motifs. CUT&RUN experiments revealed combination treatment targeting of MYC, AR, and H3K27ac again shaping enhancer accessibility. Specifically, it repressed MYC interactions at enhancer regions enriched for bHLH motifs, and overlapped with publicly-available bHLH cistromes. Finally, cistrome-transcriptome analyses identified ~200 genes that distinguished advanced PCa tumors in the SU2C cohort with high androgen and low neuroendocrine scores.

The Role of Epigenetic Modifier Mutations in Peripheral T-Cell Lymphomas

Peripheral T-cell lymphomas (PTCLs) are a group of diseases with a low incidence, high degree of heterogeneity, and a dismal prognosis in most cases. Because of the low incidence of these diseases, there have been few therapeutic novelties developed over time. Nevertheless, this fact is changing presently as epigenetic modifiers have been shown to be recurrently mutated in some types of PTCLs, especially in the cases of PTCLs not otherwise specified (PTCL-NOS), T follicular helper (TFH), and angioimmunoblastic T-cell lymphoma (AITL). These have brought about more insight into PTCL biology, especially in the case of PTCLs arising from TFH lymphocytes. From a biological perspective, it has been observed that ten-eleven translocators (TET2) mutated T lymphocytes tend to polarize to TFH, while Tregs lose their inhibitory properties. IDH2 R172 was shown to have inhibitory effects on TET2, mimicking the effects of TET2 mutations, as well as having effects on histone methylation. DNA methyltransferase 3A (DNMT3A) loss-of-function, although it was shown to have opposite effects to TET2 from an inflammatory perspective, was also shown to increase the number of T lymphocyte progenitors. Aside from bringing about more knowledge of PTCL biology, these mutations were shown to increase the sensitivity of PTCLs to certain epigenetic therapies, like hypomethylating agents (HMAs) and histone deacetylase inhibitors (HDACis). Thus, to answer the question from the title of this review: We found the Achilles heel, but only for one of the Achilles.

From multi-omics approaches to personalized medicine in myocardial infarction

Myocardial infarction (MI) is a prevalent cardiovascular disease characterized by myocardial necrosis resulting from coronary artery ischemia and hypoxia, which can lead to severe complications such as arrhythmia, cardiac rupture, heart failure, and sudden death. Despite being a research hotspot, the etiological mechanism of MI remains unclear. The emergence and widespread use of omics technologies, including genomics, transcriptomics, proteomics, metabolomics, and other omics, have provided new opportunities for exploring the molecular mechanism of MI and identifying a large number of disease biomarkers. However, a single-omics approach has limitations in understanding the complex biological pathways of diseases. The multi-omics approach can reveal the interaction network among molecules at various levels and overcome the limitations of the single-omics approaches. This review focuses on the omics studies of MI, including genomics, epigenomics, transcriptomics, proteomics, metabolomics, and other omics. The exploration extended into the domain of multi-omics integrative analysis, accompanied by a compilation of diverse online resources, databases, and tools conducive to these investigations. Additionally, we discussed the role and prospects of multi-omics approaches in personalized medicine, highlighting the potential for improving diagnosis, treatment, and prognosis of MI.

Long-range gene regulation in hormone-dependent cancer

The human genome is organized into multiple structural layers, ranging from chromosome territories to progressively smaller substructures, such as topologically associating domains (TADs) and chromatin loops. These substructures, collectively referred to as long-range chromatin interactions (LRIs), have a significant role in regulating gene expression. TADs are regions of the genome that harbour groups of genes and regulatory elements that frequently interact with each other and are insulated from other regions, thereby preventing widespread uncontrolled DNA contacts. Chromatin loops formed within TADs through enhancer and promoter interactions are elastic, allowing transcriptional heterogeneity and stochasticity. Over the past decade, it has become evident that the 3D genome structure, also referred to as the chromatin architecture, is central to many transcriptional cellular decisions. In this Review, we delve into the intricate relationship between steroid receptors and LRIs, discussing how steroid receptors interact with and modulate these chromatin interactions. Genetic alterations in the many processes involved in organizing the nuclear architecture are often associated with the development of hormone-dependent cancers. A better understanding of the interplay between architectural proteins and hormone regulatory networks can ultimately be exploited to develop improved approaches for cancer treatment.

Transposable elements in mammalian chromatin organization

Transposable elements (TEs) are mobile DNA elements that comprise almost 50% of mammalian genomic sequence. TEs are capable of making additional copies of themselves that integrate into new positions in host genomes. This unique property has had an important impact on mammalian genome evolution and on the regulation of gene expression because TE-derived sequences can function as cis-regulatory elements such as enhancers, promoters and silencers. Now, advances in our ability to identify and characterize TEs have revealed that TE-derived sequences also regulate gene expression by both maintaining and shaping 3D genome architecture. Studies are revealing how TEs contribute raw sequence that can give rise to the structures that shape chromatin organization, and thus gene expression, allowing for species-specific genome innovation and evolutionary novelty.

Small extrachromosomal circular DNAs as biomarkers for multi-cancer diagnosis and monitoring

BACKGROUND

Small extrachromosomal circular DNAs (eccDNAs) have the potential to be cancer biomarkers. However, the formation mechanisms and functions of small eccDNAs selected in carcinogenesis are not clear, and whether the small eccDNA profile in the plasma of cancer patients represents that in cancer tissues remains to be elucidated.

METHODS

A novel sequencing workflow based on the nanopore sequencing platform was used to sequence naturally existing full-length small eccDNAs in tissues and plasma collected from 25 cancer patients (including prostate cancer, hepatocellular carcinoma and colorectal cancer), and from an independent validation cohort (including 7 cancer plasma and 14 healthy plasma).

RESULTS

Compared with those in non-cancer tissues, small eccDNAs detected in cancer tissues had a significantly larger number and size (P = 0.040 and 2.2e-16, respectively), along with more even distribution and different formation mechanisms. Although small eccDNAs had different general characteristics and genomic annotation between cancer tissues and the paired plasma, they had similar formation mechanisms and cancer-related functions. Small eccDNAs originated from some specific genes had great multi-cancer diagnostic value in tissues (AUC ≥ 0.8) and plasma (AUC > 0.9), especially increasing the accuracy of multi-cancer prediction of CEA/CA19-9 levels. The high multi-cancer diagnostic value of small eccDNAs originated from ALK&ETV6 could be extrapolated from tissues (AUC = 0.804) to plasma and showed high positive predictive value (100%) and negative predictive value (82.35%) in a validation cohort.

CONCLUSIONS

As independent and stable circular DNA molecules, small eccDNAs in both tissues and plasma can be used as ideal biomarkers for cost-effective multi-cancer diagnosis and monitoring.

Identifying Cancer Type-Specific Transcriptional Programs through Network Analysis

Identifying cancer type-specific genes that define cell states is important to develop effective therapies for patients and methods for detection, early diagnosis, and prevention. While molecular mechanisms that drive malignancy have been identified for various cancers, the identification of cell-type defining transcription factors (TFs) that distinguish normal cells from cancer cells has not been fully elucidated. Here, we utilized a network biology framework, which assesses the fidelity of cell fate conversions, to identify cancer type-specific gene regulatory networks (GRN) for 17 types of cancer. Through an integrative analysis of a compendium of expression data, we elucidated core TFs and GRNs for multiple cancer types. Moreover, by comparing normal tissues and cells to cancer type-specific GRNs, we found that the expression of key network-influencing TFs can be utilized as a survival prognostic indicator for a diverse cohort of cancer patients. These findings offer a valuable resource for exploring cancer type-specific networks across a broad range of cancer types.

Integrative Identification by Hi-C Revealed Distinct Advanced Structural Variations in Lung Adenocarcinoma Tissue

Advanced three-dimensional structure variations of chromatin in large genome fragments, such as conversion of A/B compartment, topologically associated domains (TADs) and chromatin loops are related closely to occurrence of malignant tumors. However, the structural characteristics of lung cancer still remain uncovered. In this study, we used high-throughput chromosome (Hi-C) conformation capture technology to detect the advanced structural variations in chromatin of two non-smoking lung adenocarcinoma (LUAD) tumor and paired normal tissues. The results indicate that significant chromatin variations are detected in tumor tissues compared with normal tissues. At compartment scale, the main conversion type of compartment is A → B in tumor tissues, which are concentrated mainly on chromosome 3 (Chr3) (33.6%). A total of 216 tumor-specific TADs are identified in tumor tissues, which are distributed mainly in Chr1 (19), Chr2 (15) and Chr3 (17). Forty-one distinct enhancer-promoter loops are observed in tumor tissue, which are associated closely to tumor-related pathways including mitogen-activated protein kinase (MAPK), Phosphatidylinositol-3-kinase-Protein kinase B (PI3K-AKT), Ras, Wnt and Ras1. The most important observation in this study is that we identify five important genes (SYT16, NCEH1, NXPE3, MB21D2, and DZIP1L), which are detected in both A → B compartment, TADs and chromatin loops in tumor samples, and four of these genes (NCEH1, NXPE3, MB21D2, and DZIP1L) locate on q arm of Chr3. Further gene expression and invasion experiment analysis show that NCEH1, MB21D2 and SYT16 are involved in the tumor development. Thus, we provide a comprehensive overview of advanced structures in LUAD for the first time and provide a basis for further research on the genetic variation of this tumor.

Improving comparative analyses of Hi-C data via contrastive self-supervised learning

Hi-C is a widely applied chromosome conformation capture (3C)-based technique, which has produced a large number of genomic contact maps with high sequencing depths for a wide range of cell types, enabling comprehensive analyses of the relationships between biological functionalities (e.g. gene regulation and expression) and the three-dimensional genome structure. Comparative analyses play significant roles in Hi-C data studies, which are designed to make comparisons between Hi-C contact maps, thus evaluating the consistency of replicate Hi-C experiments (i.e. reproducibility measurement) and detecting statistically differential interacting regions with biological significance (i.e. differential chromatin interaction detection). However, due to the complex and hierarchical nature of Hi-C contact maps, it remains challenging to conduct systematic and reliable comparative analyses of Hi-C data. Here, we proposed sslHiC, a contrastive self-supervised representation learning framework, for precisely modeling the multi-level features of chromosome conformation and automatically producing informative feature embeddings for genomic loci and their interactions to facilitate comparative analyses of Hi-C contact maps. Comprehensive computational experiments on both simulated and real datasets demonstrated that our method consistently outperformed the state-of-the-art baseline methods in providing reliable measurements of reproducibility and detecting differential interactions with biological meanings.

Critical role of antioxidant programs in enzalutamide-resistant prostate cancer

Therapy resistance to second-generation androgen receptor (AR) antagonists, such as enzalutamide, is common in patients with advanced prostate cancer (PCa). To understand the metabolic alterations involved in enzalutamide resistance, we performed metabolomic, transcriptomic, and cistromic analyses of enzalutamide-sensitive and -resistant PCa cells, xenografts, patient-derived organoids, patient-derived explants, and tumors. We noted dramatically higher basal and inducible levels of reactive oxygen species (ROS) in enzalutamide-resistant PCa and castration-resistant PCa (CRPC), in comparison to enzalutamide-sensitive PCa cells or primary therapy-naive tumors respectively. Unbiased metabolomic evaluation identified that glutamine metabolism was consistently upregulated in enzalutamide-resistant PCa cells and CRPC tumors. Stable isotope tracing studies suggest that this enhanced glutamine metabolism drives an antioxidant program that allows these cells to tolerate higher basal levels of ROS. Inhibition of glutamine metabolism with either a small-molecule glutaminase inhibitor or genetic knockout of glutaminase enhanced ROS levels, and blocked the growth of enzalutamide-resistant PCa. The critical role of compensatory antioxidant pathways in maintaining enzalutamide-resistant PCa cells was validated by targeting another antioxidant program driver, ferredoxin 1. Taken together, our data identify a metabolic need to maintain antioxidant programs and a potentially targetable metabolic vulnerability in enzalutamide-resistant PCa.

3D genome organization and epigenetic regulation in autoimmune diseases

Three-dimensional (3D) genomics is an emerging field of research that investigates the relationship between gene regulatory function and the spatial structure of chromatin. Chromatin folding can be studied using chromosome conformation capture (3C) technology and 3C-based derivative sequencing technologies, including chromosome conformation capture-on-chip (4C), chromosome conformation capture carbon copy (5C), and high-throughput chromosome conformation capture (Hi-C), which allow scientists to capture 3D conformations from a single site to the entire genome. A comprehensive analysis of the relationships between various regulatory components and gene function also requires the integration of multi-omics data such as genomics, transcriptomics, and epigenomics. 3D genome folding is involved in immune cell differentiation, activation, and dysfunction and participates in a wide range of diseases, including autoimmune diseases. We describe hierarchical 3D chromatin organization in this review and conclude with characteristics of C-techniques and multi-omics applications of the 3D genome. In addition, we describe the relationship between 3D genome structure and the differentiation and maturation of immune cells and address how changes in chromosome folding contribute to autoimmune diseases.

Mapping Multi-factor-mediated Chromatin Interactions to Assess Dysregulation of Lung Cancer-related Genes

Studies on the lung cancer genome are indispensable for developing a cure for lung cancer. Whole-genome resequencing, genome-wide association studies, and transcriptome sequencing have greatly improved our understanding of the cancer genome. However, dysregulation of long-range chromatin interactions in lung cancer remains poorly described. To better understand the three-dimensional (3D) genomic interaction features of the lung cancer genome, we used the A549 cell line as a model system and generated high-resolution chromatin interactions associated with RNA polymerase II (RNAPII), CCCTC-binding factor (CTCF), enhancer of zeste homolog 2 (EZH2), and histone 3 lysine 27 trimethylation (H3K27me3) using long-read chromatin interaction analysis by paired-end tag sequencing (ChIA-PET). Analysis showed that EZH2/H3K27me3-mediated interactions further repressed target genes, either through loops or domains, and their distributions along the genome were distinct from and complementary to those associated with RNAPII. Cancer-related genes were highly enriched with chromatin interactions, and chromatin interactions specific to the A549 cell line were associated with oncogenes and tumor suppressor genes, such as additional repressive interactions on FOXO4 and promoter-promoter interactions between NF1 and RNF135. Knockout of an anchor associated with chromatin interactions reversed the dysregulation of cancer-related genes, suggesting that chromatin interactions are essential for proper expression of lung cancer-related genes. These findings demonstrate the 3D landscape and gene regulatory relationships of the lung cancer genome.

LPAD: using network construction and label propagation to detect topologically associating domains from Hi-C data

With the development of chromosome conformation capture technique, the study of spatial conformation of a genome based on Hi-C technique has made a quantum leap. Previous studies reveal that genomes are folded into hierarchy of three-dimensional (3D) structures associated with topologically associating domains (TADs), and detecting TAD boundaries is of great significance in the chromosome-level analysis of 3D genome architecture. In this paper, we propose a novel TAD identification method, LPAD, which first extracts node correlations from global interactions of chromosomes based on the random walk with restart and then builds an undirected graph from Hi-C contact matrix. Next, LPAD designs a label propagation-based approach to discover communities and generates TADs. Experimental results verify the effectiveness and quality of TAD detections compared with existing methods. Furthermore, experimental evaluation of chromatin immunoprecipitation sequencing data shows that LPAD performs high enrichment of histone modifications remarkably nearby the TAD boundaries, and these results demonstrate LPAD's advantages on TAD identification accuracy.

A graph neural network-based interpretable framework reveals a novel DNA fragility-associated chromatin structural unit

BACKGROUND

DNA double-strand breaks (DSBs) are among the most deleterious DNA lesions, and they can cause cancer if improperly repaired. Recent chromosome conformation capture techniques, such as Hi-C, have enabled the identification of relationships between the 3D chromatin structure and DSBs, but little is known about how to explain these relationships, especially from global contact maps, or their contributions to DSB formation.

RESULTS

Here, we propose a framework that integrates graph neural network (GNN) to unravel the relationship between 3D chromatin structure and DSBs using an advanced interpretable technique GNNExplainer. We identify a new chromatin structural unit named the DNA fragility-associated chromatin interaction network (FaCIN). FaCIN is a bottleneck-like structure, and it helps to reveal a universal form of how the fragility of a piece of DNA might be affected by the whole genome through chromatin interactions. Moreover, we demonstrate that neck interactions in FaCIN can serve as chromatin structural determinants of DSB formation.

CONCLUSIONS

Our study provides a more systematic and refined view enabling a better understanding of the mechanisms of DSB formation under the context of the 3D genome.

Epigenetic Modifications in Prostate Cancer Metastasis and Microenvironment

The gradual evolution of prostate tissue from benign tumor to malignant lesion or distant metastasis is driven by intracellular epigenetic changes and the tumor microenvironment remodeling. With the continuous study of epigenetic modifications, these tumor-driving forces are being discovered and are providing new treatments for cancer. Here we introduce the classification of epigenetic modification and highlight the role of epigenetic modification in tumor remodeling and communication of the tumor microenvironment.

Rewiring of the 3D genome during acquisition of carboplatin resistance in a triple-negative breast cancer patient-derived xenograft

Changes in the three-dimensional (3D) structure of the genome are an emerging hallmark of cancer. Cancer-associated copy number variants and single nucleotide polymorphisms promote rewiring of chromatin loops, disruption of topologically associating domains (TADs), active/inactive chromatin state switching, leading to oncogene expression and silencing of tumor suppressors. However, little is known about 3D changes during cancer progression to a chemotherapy-resistant state. We integrated chromatin conformation capture (Hi-C), RNA-seq, and whole-genome sequencing obtained from triple-negative breast cancer patient-derived xenograft primary tumors (UCD52) and carboplatin-resistant samples and found increased short-range (< 2 Mb) interactions, chromatin looping, formation of TAD, chromatin state switching into a more active state, and amplification of ATP-binding cassette transporters. Transcriptome changes suggested the role of long-noncoding RNAs in carboplatin resistance. Rewiring of the 3D genome was associated with TP53, TP63, BATF, FOS-JUN family of transcription factors and led to activation of aggressiveness-, metastasis- and other cancer-related pathways. Integrative analysis highlighted increased ribosome biogenesis and oxidative phosphorylation, suggesting the role of mitochondrial energy metabolism. Our results suggest that 3D genome remodeling may be a key mechanism underlying carboplatin resistance.

African American Prostate Cancer Displays Quantitatively Distinct Vitamin D Receptor Cistrome-transcriptome Relationships Regulated by BAZ1A

African American (AA) prostate cancer associates with vitamin D3 deficiency, but vitamin D receptor (VDR) genomic actions have not been investigated in this context. We undertook VDR proteogenomic analyses in European American (EA) and AA prostate cell lines and four clinical cohorts. Rapid immunoprecipitation mass spectrometry of endogenous protein (RIME) analyses revealed that nonmalignant AA RC43N prostate cells displayed the greatest dynamic protein content in the VDR complex. Likewise, in AA cells, Assay for Transposase-Accessible Chromatin using sequencing established greater 1α,25(OH)2D3-regulated chromatin accessibility, chromatin immunoprecipitation sequencing revealed significant enhancer-enriched VDR cistrome, and RNA sequencing identified the largest 1α,25(OH)2D3-dependent transcriptome. These VDR functions were significantly corrupted in the isogenic AA RC43T prostate cancer cells, and significantly distinct from EA cell models. We identified reduced expression of the chromatin remodeler, BAZ1A, in three AA prostate cancer cohorts as well as RC43T compared with RC43N. Restored BAZ1A expression significantly increased 1α,25(OH)2D3-regulated VDR-dependent gene expression in RC43T, but not HPr1AR or LNCaP cells. The clinical impact of VDR cistrome-transcriptome relationships were tested in three different clinical prostate cancer cohorts. Strikingly, only in AA patients with prostate cancer, the genes bound by VDR and/or associated with 1α,25(OH)2D3-dependent open chromatin (i) predicted progression from high-grade prostatic intraepithelial neoplasia to prostate cancer; (ii) responded to vitamin D3 supplementation in prostate cancer tumors; (iii) differentially responded to 25(OH)D3 serum levels. Finally, partial correlation analyses established that BAZ1A and components of the VDR complex identified by RIME significantly strengthened the correlation between VDR and target genes in AA prostate cancer only. Therefore, VDR transcriptional control is most potent in AA prostate cells and distorted through a BAZ1A-dependent control of VDR function.

Significance

Our study identified that genomic ancestry drives the VDR complex composition, genomic distribution, and transcriptional function, and is disrupted by BAZ1A and illustrates a novel driver for AA prostate cancer.

Defining the separation landscape of topological domains for decoding consensus domain organization of the 3D genome

Topologically associating domains (TADs) have emerged as basic structural and functional units of genome organization and have been determined by many computational methods from Hi-C contact maps. However, the TADs obtained by different methods vary greatly, which makes the accurate determination of TADs a challenging issue and hinders subsequent biological analyses about their organization and functions. Obvious inconsistencies among the TADs identified by different methods indeed make the statistical and biological properties of TADs overly depend on the chosen method rather than on the data. To this end, we use the consensus structural information captured by these methods to define the TAD separation landscape for decoding the consensus domain organization of the 3D genome. We show that the TAD separation landscape could be used to compare domain boundaries across multiple cell types for discovering conserved and divergent topological structures, decipher three types of boundary regions with diverse biological features, and identify consensus TADs (ConsTADs). We illustrate that these analyses could deepen our understanding of the relationships between the topological domains and chromatin states, gene expression, and DNA replication timing.

Emerging regulatory mechanisms of noncoding RNAs in topologically associating domains

Topologically associating domains (TADs) are integral to spatial genome organization, instructing gene expression, and cell fate. Recently, several advances have uncovered roles for noncoding RNAs (ncRNAs) in the regulation of the form and function of mammalian TADs. Phase separation has also emerged as a potential arbiter of ncRNAs in the regulation of TADs. In this review we discuss the implications of these novel findings in relation to how ncRNAs might structurally and functionally regulate TADs from two perspectives: moderating loop extrusion through interactions with architectural proteins, and facilitating TAD phase separation. Additionally, we propose future studies and directions to investigate these phenomena.

Epigenetic regulation in hematopoiesis and its implications in the targeted therapy of hematologic malignancies

Hematologic malignancies are one of the most common cancers, and the incidence has been rising in recent decades. The clinical and molecular features of hematologic malignancies are highly heterogenous, and some hematologic malignancies are incurable, challenging the treatment, and prognosis of the patients. However, hematopoiesis and oncogenesis of hematologic malignancies are profoundly affected by epigenetic regulation. Studies have found that methylation-related mutations, abnormal methylation profiles of DNA, and abnormal histone deacetylase expression are recurrent in leukemia and lymphoma. Furthermore, the hypomethylating agents and histone deacetylase inhibitors are effective to treat acute myeloid leukemia and T-cell lymphomas, indicating that epigenetic regulation is indispensable to hematologic oncogenesis. Epigenetic regulation mainly includes DNA modifications, histone modifications, and noncoding RNA-mediated targeting, and regulates various DNA-based processes. This review presents the role of writers, readers, and erasers of DNA methylation and histone methylation, and acetylation in hematologic malignancies. In addition, this review provides the influence of microRNAs and long noncoding RNAs on hematologic malignancies. Furthermore, the implication of epigenetic regulation in targeted treatment is discussed. This review comprehensively presents the change and function of each epigenetic regulator in normal and oncogenic hematopoiesis and provides innovative epigenetic-targeted treatment in clinical practice.

A Comprehensive Investigation of Genomic Variants in Prostate Cancer Reveals 30 Putative Regulatory Variants

Prostate cancer (PC) is the most frequently diagnosed non-skin cancer in the world. Previous studies have shown that genomic alterations represent the most common mechanism for molecular alterations responsible for the development and progression of PC. This highlights the importance of identifying functional genomic variants for early detection in high-risk PC individuals. Great efforts have been made to identify common protein-coding genetic variations; however, the impact of non-coding variations, including regulatory genetic variants, is not well understood. Identification of these variants and the underlying target genes will be a key step in improving the detection and treatment of PC. To gain an understanding of the functional impact of genetic variants, and in particular, regulatory variants in PC, we developed an integrative pipeline (AGV) that uses whole genome/exome sequences, GWAS SNPs, chromosome conformation capture data, and ChIP-Seq signals to investigate the potential impact of genomic variants on the underlying target genes in PC. We identified 646 putative regulatory variants, of which 30 significantly altered the expression of at least one protein-coding gene. Our analysis of chromatin interactions data (Hi-C) revealed that the 30 putative regulatory variants could affect 131 coding and non-coding genes. Interestingly, our study identified the 131 protein-coding genes that are involved in disease-related pathways, including Reactome and MSigDB, for most of which targeted treatment options are currently available. Notably, our analysis revealed several non-coding RNAs, including RP11-136K7.2 and RAMP2-AS1, as potential enhancer elements of the protein-coding genes CDH12 and EZH1, respectively. Our results provide a comprehensive map of genomic variants in PC and reveal their potential contribution to prostate cancer progression and development.

Three-dimensional genome landscape comprehensively reveals patterns of spatial gene regulation in papillary and anaplastic thyroid cancers: a study using representative cell lines for each cancer type

BACKGROUND

Spatial chromatin structure is intricately linked with somatic aberrations, and somatic mutations of various cancer-related genes, termed co-mutations (CoMuts), occur in certain patterns during cancer initiation and progression. The functional mechanisms underlying these genetic events remain largely unclear in thyroid cancer (TC). With discrepant differentiation, papillary thyroid cancer (PTC) and anaplastic thyroid cancer (ATC) differ greatly in characteristics and prognosis. We aimed to reveal the spatial gene alterations and regulations between the two TC subtypes.

METHODS

We systematically investigated and compared the spatial co-mutations between ATC (8305C), PTC (BCPAP and TPC-1), and normal thyroid cells (Nthy-ori-3-1). We constructed a framework integrating whole-genome sequencing (WGS), high-throughput chromosome conformation capture (Hi-C), and transcriptome sequencing, to systematically detect the associations between the somatic co-mutations of cancer-related genes, structural variations (SVs), copy number variations (CNVs), and high-order chromatin conformation.

RESULTS

Spatial co-mutation hotspots were enriched around topologically associating domains (TADs) in TC. A common set of 227 boundaries were identified in both ATC and PTC, with significant overlaps between them. The spatial proximities of the co-mutated gene pairs in the two TC types were significantly greater than in the gene-level and overall backgrounds, and ATC cells had higher TAD contact frequency with CoMuts > 10 compared with PTC cells. Compared with normal thyroid cells, in ATC the number of the created novel three-dimensional chromatin structural domains increased by 10%, and the number of shifted TADs decreased by 7%. We found five TAD blocks with CoMut genes/events specific to ATC with certain mutations in genes including MAST-NSUN4, AM129B/TRUB2, COL5A1/PPP1R26, PPP1R26/GPSM1/CCDC183, and PRAC2/DLX4. For the majority of ATC and PTC cells, the HOXA10 and HIF2α signals close to the transcription start sites of CoMut genes within TADs were significantly stronger than those at the background. CNV breakpoints significantly overlapped with TAD boundaries in both TC subtypes. ATCs had more CNV losses overlapping with TAD boundaries, and noncoding SVs involved in intrachromosomal SVs, amplified inversions, and tandem duplication differed between ATC and PTC. TADs with short range were more abundant in ATC than PTC. More switches of A/B compartment types existed in ATC cells compared with PTC. Gene expression was significantly synchronized, and orchestrated by complex epigenetics and regulatory elements.

CONCLUSION

Chromatin interactions and gene alterations and regulations are largely heterogeneous in TC. CNVs and complex SVs may function in the TC genome by interplaying with TADs, and are largely different between ATC and PTC. Complexity of TC genomes, which are highly organized by 3D genome-wide interactions mediating mutational and structural variations and gene activation, may have been largely underappreciated. Our comprehensive analysis may provide key evidence and targets for more customized diagnosis and treatment of TC.

The dynamics of three-dimensional chromatin organization and phase separation in cell fate transitions and diseases

Cell fate transition is a fascinating process involving complex dynamics of three-dimensional (3D) chromatin organization and phase separation, which play an essential role in cell fate decision by regulating gene expression. Phase separation is increasingly being considered a driving force of chromatin folding. In this review, we have summarized the dynamic features of 3D chromatin and phase separation during physiological and pathological cell fate transitions and systematically analyzed recent evidence of phase separation facilitating the chromatin structure. In addition, we discuss current advances in understanding how phase separation contributes to physical and functional enhancer-promoter contacts. We highlight the functional roles of 3D chromatin organization and phase separation in cell fate transitions, and more explorations are required to study the regulatory relationship between 3D chromatin organization and phase separation. 3D chromatin organization (shown by Hi-C contact map) and phase separation are highly dynamic and play functional roles during early embryonic development, cell differentiation, somatic reprogramming, cell transdifferentiation and pathogenetic process. Phase separation can regulate 3D chromatin organization directly, but whether 3D chromatin organization regulates phase separation remains unclear.

Three-Dimensional Organization of Chicken Genome Provides Insights into Genetic Adaptation to Extreme Environments

The high-throughput chromosome conformation capture (Hi-C) technique is widely used to study the functional roles of the three-dimensional (3D) architecture of genomes. However, the knowledge of the 3D genome structure and its dynamics during extreme environmental adaptations remains poor. Here, we characterized 3D genome architectures using the Hi-C technique for chicken liver cells. Upon comparing Lindian chicken (LDC) liver cells with Wenchang chicken (WCC) liver cells, we discovered that environmental adaptation contributed to the switching of A/B compartments, the reorganization of topologically associated domains (TADs), and TAD boundaries in both liver cells. In addition, the analysis of the switching of A/B compartments revealed that the switched compartmental genes (SCGs) were strongly associated with extreme environment adaption-related pathways, including tight junction, notch signaling pathway, vascular smooth muscle contraction, and the RIG-I-like receptor signaling pathway. The findings of this study advanced our understanding of the evolutionary role of chicken 3D genome architecture and its significance in genome activity and transcriptional regulation.

3D-Epigenomic Regulation of Gene Transcription in Hepatocellular Carcinoma

The fundamental cause of transcription dysregulation in hepatocellular carcinoma (HCC) remains elusive. To investigate the underlying mechanisms, comprehensive 3D-epigenomic analyses are performed in cellular models of THLE2 (a normal hepatocytes cell line) and HepG2 (a hepatocellular carcinoma cell line) using integrative approaches for chromatin topology, genomic and epigenomic variation, and transcriptional output. Comparing the 3D-epigenomes in THLE2 and HepG2 reveal that most HCC-associated genes are organized in complex chromatin interactions mediated by RNA polymerase II (RNAPII). Incorporation of genome-wide association studies (GWAS) data enables the identification of non-coding genetic variants that are enriched in distal enhancers connecting to the promoters of HCC-associated genes via long-range chromatin interactions, highlighting their functional roles. Interestingly, CTCF binding and looping proximal to HCC-associated genes appear to form chromatin architectures that overarch RNAPII-mediated chromatin interactions. It is further demonstrated that epigenetic variants by DNA hypomethylation at a subset of CTCF motifs proximal to HCC-associated genes can modify chromatin topological configuration, which in turn alter RNAPII-mediated chromatin interactions and lead to dysregulation of transcription. Together, the 3D-epigenomic analyses provide novel insights of multifaceted interplays involving genetics, epigenetics, and chromatin topology in HCC cells.

Trans-interaction of risk loci 6p24.1 and 10q11.21 is associated with endothelial damage in coronary artery disease

BACKGROUND AND AIMS

Single nucleotide polymorphism rs6903956 has been identified as one of the genetic risk factors for coronary artery disease (CAD). However, rs6903956 lies in a non-coding locus on chromosome 6p24.1. We aim to interrogate the molecular basis of 6p24.1 containing rs6903956 risk alleles in endothelial disease biology.

METHODS AND RESULTS

We generated induced pluripotent stem cells (iPSCs) from CAD patients (AA risk genotype at rs6903956) and non-CAD subjects (GG non-risk genotype at rs6903956). CRISPR-Cas9-based deletions (Δ63-89bp) on 6p24.1, including both rs6903956 and a short tandem repeat variant rs140361069 in linkage disequilibrium, were performed to generate isogenic iPSC-derived endothelial cells. Edited CAD endothelial cells, with removal of 'A' risk alleles, exhibited a global transcriptional downregulation of pathways relating to abnormal vascular physiology and activated endothelial processes. A CXC chemokine ligand on chromosome 10q11.21, CXCL12, was uncovered as a potential effector gene in CAD endothelial cells. Underlying this effect was the preferential inter-chromosomal interaction of 6p24.1 risk locus to a weak promoter of CXCL12, confirmed by chromatin conformation capture assays on our iPSC-derived endothelial cells. Functionally, risk genotypes AA/AG at rs6903956 were associated significantly with elevated levels of circulating damaged endothelial cells in CAD patients. Circulating endothelial cells isolated from patients with risk genotypes AA/AG were also found to have 10 folds higher CXCL12 transcript copies/cell than those with non-risk genotype GG.

CONCLUSIONS

Our study reveals the trans-acting impact of 6p24.1 with another CAD locus on 10q11.21 and is associated with intensified endothelial injury.

Collaborative Sense-Making in Genomic Research: The Role of Visualisation

Genomic research emerges from collaborative work within and across different scientific disciplines. A diverse range of visualisation techniques has been employed to aid this research, yet relatively little is known as to how these techniques facilitate collaboration. We conducted a case study of collaborative research within a biomedical institute to learn more about the role visualisation plays in genomic mapping. Interviews were conducted with molecular biologists (N = 5) and bioinformaticians (N = 6). We found that genomic research comprises a variety of distinct disciplines engaged in complex analytic tasks that each resist simplification, and their complexity influences how visualisations were used. Visualisation use was impacted by group-specific interactions and temporal work patterns. Visualisations were also crucial to the scientific workflow, used for both question formation and confirmation of hypotheses, and acted as an anchor for the communication of ideas and discussion. In the latter case, two approaches were taken: providing collaborators with either interactive or static imagery representing a viewpoint. The use of generic software for simplified visualisations, and quick production and curation was also noted. We discuss these findings with reference to group-specific interactions and present recommendations for improving collaborative practices through visual analytics.

Structural variants drive context-dependent oncogene activation in cancer

Higher-order chromatin structure is important for the regulation of genes by distal regulatory sequences^1,2. Structural variants (SVs) that alter three-dimensional (3D) genome organization can lead to enhancer-promoter rewiring and human disease, particularly in the context of cancer³. However, only a small minority of SVs are associated with altered gene expression^4,5, and it remains unclear why certain SVs lead to changes in distal gene expression and others do not. To address these questions, we used a combination of genomic profiling and genome engineering to identify sites of recurrent changes in 3D genome structure in cancer and determine the effects of specific rearrangements on oncogene activation. By analysing Hi-C data from 92 cancer cell lines and patient samples, we identified loci affected by recurrent alterations to 3D genome structure, including oncogenes such as MYC, TERT and CCND1. By using CRISPR-Cas9 genome engineering to generate de novo SVs, we show that oncogene activity can be predicted by using 'activity-by-contact' models that consider partner region chromatin contacts and enhancer activity. However, activity-by-contact models are only predictive of specific subsets of genes in the genome, suggesting that different classes of genes engage in distinct modes of regulation by distal regulatory elements. These results indicate that SVs that alter 3D genome organization are widespread in cancer genomes and begin to illustrate predictive rules for the consequences of SVs on oncogene activation.

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.

Constricted migration is associated with stable 3D genome structure differences in cancer cells

To spread from a localized tumor, metastatic cancer cells must squeeze through constrictions that cause major nuclear deformations. Since chromosome structure affects nucleus stiffness, gene regulation, and DNA repair, here, we investigate the relationship between 3D genome structure and constricted migration in cancer cells. Using melanoma (A375) cells, we identify phenotypic differences in cells that have undergone multiple rounds of constricted migration. These cells display a stably higher migration efficiency, elongated morphology, and differences in the distribution of Lamin A/C and heterochromatin. Hi-C experiments reveal differences in chromosome spatial compartmentalization specific to cells that have passed through constrictions and related alterations in expression of genes associated with migration and metastasis. Certain features of the 3D genome structure changes, such as a loss of B compartment interaction strength, are consistently observed after constricted migration in clonal populations of A375 cells and in MDA-MB-231 breast cancer cells. Our observations suggest that consistent types of chromosome structure changes are induced or selected by passage through constrictions and that these may epigenetically encode stable differences in gene expression and cellular migration phenotype.

Reorganization of 3D chromatin architecture in doxorubicin-resistant breast cancer cells

Background: Doxorubicin resistance remains a major therapeutic challenge leading to poor survival prognosis and treatment failure in breast cancer. Although doxorubicin induces massive changes in the transcriptional landscape are well known, potential diagnostic or therapeutic targets associated with the reorganization of three-dimensional (3D) chromatin architecture have not yet been systematically investigated.

Methods: Here we performed in situ high-throughput chromosome conformation capture (Hi-C) on parental and doxorubicin-resistant MCF7 (MCF7-DR) human breast cancer cells, followed by integrative analysis of HiC, ATAC-seq, RNA-seq and TCGA data.

Results: It revealed that A/B compartment switching was positively correlated to genome-wide differential gene expression. The genome of MCF7-DR cells was spatially reorganized into smaller topologically associating domains (TADs) and chromatin loops. We also revealed the contribution of increased chromatin accessibility and potential transcription factor families, including CTCF, AP-1 and bHLH, to gained TADs or loops. Intriguingly, we observed two condensed genomic regions (∼20 kb) with decreased chromatin accessibility flanking TAD boundaries, which might play a critical role in the formation or maintenance of TADs. Finally, combining data from TCGA, we identified a number of gained and lost enhancer-promoter interactions and their corresponding differentially expressed genes involved in chromatin organization and breast cancer signaling pathways, including FA2H, FOXA1 and JRKL, which might serve as potential treatment targets for breast cancer.

Conclusion: These data uncovered a close connection between 3D genome reorganization, chromatin accessibility as well as gene transcription and provide novel insights into the epigenomic mechanisms involving doxorubicin resistance in breast cancer.

Neuropilin-2 promotes lineage plasticity and progression to neuroendocrine prostate cancer

Neuroendocrine prostate cancer (NEPC), a lethal subset of prostate cancer, is characterized by loss of AR signaling and resulting resistance to AR-targeted therapy during neuroendocrine transdifferentiation, for which the molecular mechanisms remain unclear. Here, we report that neuropilin 2 (NRP2) is upregulated in both de novo and therapy-induced NEPC, which induces neuroendocrine markers, neuroendocrine cell morphology, and NEPC cell aggressive behavior. NRP2 silencing restricted NEPC tumor xenograft growth. Mechanistically, NRP2 engages in reciprocal crosstalk with AR, where NRP2 is transcriptionally inhibited by AR, and in turn suppresses AR signaling by downregulating the AR transcriptional program and confers resistance to enzalutamide. Moreover, NRP2 physically interacts with VEGFR2 through the intracellular SEA domain to activate STAT3 phosphorylation and subsequently SOX2, thus driving NEPC differentiation and growth. Collectively, these results characterize NRP2 as a driver of NEPC and suggest NRP2 as a potential therapeutic target in NEPC.

Androgen receptor signaling and spatial chromatin organization in castration-resistant prostate cancer

Prostate cancer is one of the leading causes of cancer death and affects millions of men in the world. The American Cancer Society estimated about 34,500 deaths from prostate cancer in the United States in year 2022. The Androgen receptor (AR) signaling is a major pathway that sustains local and metastatic prostate tumor growth. Androgen-deprivation therapy (ADT) is the standard of care for metastatic prostate cancer patient and can suppress the tumor growth for a median of 2–3 years. Unfortunately, the malignancy inevitably progresses to castration-resistant prostate cancer (CRPC) which is more aggressive and no longer responsive to ADT. Surprisingly, for most of the CPRC patients, cancer growth still depends on androgen receptor signaling. Accumulating evidence suggests that CRPC cells have rewired their transcriptional program to retain AR signaling in the absence of androgens. Besides AR, other transcription factors also contribute to the resistance mechanism through multiple pathways including enhancing AR signaling pathway and activating other complementary signaling pathways for the favor of AR downstream genes expression. More recent studies have shown the role of transcription factors in reconfiguring chromatin 3D structure and regulating topologically associating domains (TADs). Pioneer factors, transcription factors and coactivators form liquid-liquid phase separation compartment that can modulate transcriptional events along with configuring TADs. The role of AR and other transcription factors on chromatin structure change and formation of condensate compartment in prostate cancer cells has only been recently investigated and appreciated. This review intends to provide an overview of transcription factors that contribute to AR signaling through activation of gene expression, governing 3D chromatin structure and establishing phase to phase separation. A more detailed understanding of the spatial role of transcription factors in CRPC might provide novel therapeutic targets for the treatment of CRPC.