Epigenetic and transcriptomic characterization reveals progression markers and essential pathways in clear cell renal cell carcinoma.

Identifying tumor-cell-specific markers and elucidating their epigenetic regulation and spatial heterogeneity provides mechanistic insights into cancer etiology. Here, we perform snRNA-seq and snATAC-seq in 34 and 28 human clear cell renal cell carcinoma (ccRCC) specimens, respectively, with matched bulk proteogenomics data. By identifying 20 tumor-specific markers through a multi-omics tiered approach, we reveal an association between higher ceruloplasmin (CP) expression and reduced survival. CP knockdown, combined with spatial transcriptomics, suggests a role for CP in regulating hyalinized stroma and tumor-stroma interactions in ccRCC. Intratumoral heterogeneity analysis portrays tumor cell-intrinsic inflammation and epithelial-mesenchymal transition (EMT) as two distinguishing features of tumor subpopulations. Finally, BAP1 mutations are associated with widespread reduction of chromatin accessibility, while PBRM1 mutations generally increase accessibility, with the former affecting five times more accessible peaks than the latter. These integrated analyses reveal the cellular architecture of ccRCC, providing insights into key markers and pathways in ccRCC tumorigenesis.

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.

Nuclear pore protein NUP210 depletion suppresses metastasis through heterochromatin-mediated disruption of tumor cell mechanical response.

Mechanical signals from the extracellular microenvironment have been implicated in tumor and metastatic progression. Here, we identify nucleoporin NUP210 as a metastasis susceptibility gene for human estrogen receptor positive (ER+) breast cancer and a cellular mechanosensor. Nup210 depletion suppresses lung metastasis in mouse models of breast cancer. Mechanistically, NUP210 interacts with LINC complex protein SUN2 which connects the nucleus to the cytoskeleton. In addition, the NUP210/SUN2 complex interacts with chromatin via the short isoform of BRD4 and histone H3.1/H3.2 at the nuclear periphery. In Nup210 knockout cells, mechanosensitive genes accumulate H3K27me3 heterochromatin modification, mediated by the polycomb repressive complex 2 and differentially reposition within the nucleus. Transcriptional repression in Nup210 knockout cells results in defective mechanotransduction and focal adhesion necessary for their metastatic capacity. Our study provides an important role of nuclear pore protein in cellular mechanosensation and metastasis.

RUNX1-mediated alphaherpesvirus-host trans-species chromatin interaction promotes viral transcription.

Like most DNA viruses, herpesviruses precisely deliver their genomes into the sophisticatedly organized nuclei of the infected host cells to initiate subsequent transcription and replication. However, it remains elusive how the viral genome specifically interacts with the host genome and hijacks host transcription machinery. Using pseudorabies virus (PRV) as model virus, we performed chromosome conformation capture assays to demonstrate a genome-wide specific trans-species chromatin interaction between the virus and host. Our data show that the PRV genome is delivered by the host DNA binding protein RUNX1 into the open chromatin and active transcription zone. This facilitates virus hijacking host RNAPII to efficiently transcribe viral genes, which is significantly inhibited by either a RUNX1 inhibitor or RNA interference. Together, these findings provide insights into the chromatin interaction between viral and host genomes and identify new areas of research to advance the understanding of herpesvirus genome transcription.

Macrophage mediated recognition and clearance of Borrelia burgdorferi elicits MyD88-dependent and -independent phagosomal signals that contribute to phagocytosis and inflammation.

BACKGROUND: Macrophages play prominent roles in bacteria recognition and clearance, including Borrelia burgdorferi (Bb), the Lyme disease spirochete. To elucidate mechanisms by which MyD88/TLR signaling enhances clearance of Bb by macrophages, we studied wildtype (WT) and MyD88-/- Bb-stimulated bone marrow-derived macrophages (BMDMs). RESULTS: MyD88-/- BMDMs exhibit impaired uptake of spirochetes but comparable maturation of phagosomes following internalization of spirochetes. RNA-sequencing of infected WT and MyD88-/- BMDMs identified a large cohort of differentially expressed MyD88-dependent genes associated with re-organization of actin and cytoskeleton during phagocytosis along with several MyD88-independent chemokines involved in inflammatory cell recruitment. We computationally generated networks which identified several MyD88-dependent intermediate proteins (Rhoq and Cyfip1) that are known to mediate inflammation and phagocytosis respectively. CONCLUSION: Our findings show that MyD88 signaling enhances, but is not required, for bacterial uptake or phagosomal maturation and provide mechanistic insights into how MyD88-mediated phagosomal signaling enhances Bb uptake and clearance.

Oncogenic extrachromosomal DNA functions as mobile enhancers to globally amplify chromosomal transcription.

Extrachromosomal, circular DNA (ecDNA) is emerging as a prevalent yet less characterized oncogenic alteration in cancer genomes. We leverage ChIA-PET and ChIA-Drop chromatin interaction assays to characterize genome-wide ecDNA-mediated chromatin contacts that impact transcriptional programs in cancers. ecDNAs in glioblastoma patient-derived neurosphere and prostate cancer cell cultures are marked by widespread intra-ecDNA and genome-wide chromosomal interactions. ecDNA-chromatin contact foci are characterized by broad and high-level H3K27ac signals converging predominantly on chromosomal genes of increased expression levels. Prostate cancer cells harboring synthetic ecDNA circles composed of characterized enhancers result in the genome-wide activation of chromosomal gene transcription. Deciphering the chromosomal targets of ecDNAs at single-molecule resolution reveals an association with actively expressed oncogenes spatially clustered within ecDNA-directed interaction networks. Our results suggest that ecDNA can function as mobile transcriptional enhancers to promote tumor progression and manifest a potential synthetic aneuploidy mechanism of transcription control in cancer.

Chromatin topology reorganization and transcription repression by PML-RARα in acute promyeloid leukemia.

BACKGROUND: Acute promyeloid leukemia (APL) is characterized by the oncogenic fusion protein PML-RARα, a major etiological agent in APL. However, the molecular mechanisms underlying the role of PML-RARα in leukemogenesis remain largely unknown. RESULTS: Using an inducible system, we comprehensively analyze the 3D genome organization in myeloid cells and its reorganization after PML-RARα induction and perform additional analyses in patient-derived APL cells with native PML-RARα. We discover that PML-RARα mediates extensive chromatin interactions genome-wide. Globally, it redefines the chromatin topology of the myeloid genome toward a more condensed configuration in APL cells; locally, it intrudes RNAPII-associated interaction domains, interrupts myeloid-specific transcription factors binding at enhancers and super-enhancers, and leads to transcriptional repression of genes critical for myeloid differentiation and maturation. CONCLUSIONS: Our results not only provide novel topological insights for the roles of PML-RARα in transforming myeloid cells into leukemia cells, but further uncover a topological framework of a molecular mechanism for oncogenic fusion proteins in cancers.

Coordinated demethylation of H3K9 and H3K27 is required for rapid inflammatory responses of endothelial cells.

Histone H3 lysine-9 di-methylation (H3K9me2) and lysine-27 tri-methylation (H3K27me3) are linked to repression of gene expression, but the functions of repressive histone methylation dynamics during inflammatory responses remain enigmatic. Here, we report that lysine demethylases 7A (KDM7A) and 6A (UTX) play crucial roles in tumor necrosis factor (TNF)-α signaling in endothelial cells (ECs), where they are regulated by a novel TNF-α-responsive microRNA, miR-3679-5p. TNF-α rapidly induces co-occupancy of KDM7A and UTX at nuclear factor kappa-B (NF-κB)-associated elements in human ECs. KDM7A and UTX demethylate H3K9me2 and H3K27me3, respectively, and are both required for activation of NF-κB-dependent inflammatory genes. Chromosome conformation capture-based methods furthermore uncover increased interactions between TNF-α-induced super enhancers at NF-κB-relevant loci, coinciding with KDM7A and UTX recruitments. Simultaneous pharmacological inhibition of KDM7A and UTX significantly reduces leukocyte adhesion in mice, establishing the biological and potential translational relevance of this mechanism. Collectively, these findings suggest that rapid erasure of repressive histone marks by KDM7A and UTX is essential for NF-κB-dependent regulation of genes that control inflammatory responses of ECs.

An AR-ERG transcriptional signature defined by long-range chromatin interactomes in prostate cancer cells.

The aberrant activities of transcription factors such as the androgen receptor (AR) underpin prostate cancer development. While the AR cis-regulation has been extensively studied in prostate cancer, information pertaining to the spatial architecture of the AR transcriptional circuitry remains limited. In this paper, we propose a novel framework to profile long-range chromatin interactions associated with AR and its collaborative transcription factor, erythroblast transformation-specific related gene (ERG), using chromatin interaction analysis by paired-end tag (ChIA-PET). We identified ERG-associated long-range chromatin interactions as a cooperative component in the AR-associated chromatin interactome, acting in concert to achieve coordinated regulation of a subset of AR target genes. Through multifaceted functional data analysis, we found that AR-ERG interaction hub regions are characterized by distinct functional signatures, including bidirectional transcription and cotranscription factor binding. In addition, cancer-associated long noncoding RNAs were found to be connected near protein-coding genes through AR-ERG looping. Finally, we found strong enrichment of prostate cancer genome-wide association study (GWAS) single nucleotide polymorphisms (SNPs) at AR-ERG co-binding sites participating in chromatin interactions and gene regulation, suggesting GWAS target genes identified from chromatin looping data provide more biologically relevant findings than using the nearest gene approach. Taken together, our results revealed an AR-ERG-centric higher-order chromatin structure that drives coordinated gene expression in prostate cancer progression and the identification of potential target genes for therapeutic intervention.

The Tandem Duplicator Phenotype Is a Prevalent Genome-Wide Cancer Configuration Driven by Distinct Gene Mutations.

The tandem duplicator phenotype (TDP) is a genome-wide instability configuration primarily observed in breast, ovarian, and endometrial carcinomas. Here, we stratify TDP tumors by classifying their tandem duplications (TDs) into three span intervals, with modal values of 11 kb, 231 kb, and 1.7 Mb, respectively. TDPs with ∼11 kb TDs feature loss of TP53 and BRCA1. TDPs with ∼231 kb and ∼1.7 Mb TDs associate with CCNE1 pathway activation and CDK12 disruptions, respectively. We demonstrate that p53 and BRCA1 conjoint abrogation drives TDP induction by generating short-span TDP mammary tumors in genetically modified mice lacking them. Lastly, we show how TDs in TDP tumors disrupt heterogeneous combinations of tumor suppressors and chromatin topologically associating domains while duplicating oncogenes and super-enhancers.

The Energetics and Physiological Impact of Cohesin Extrusion.

Cohesin extrusion is thought to play a central role in establishing the architecture of mammalian genomes. However, extrusion has not been visualized in vivo, and thus, its functional impact and energetics are unknown. Using ultra-deep Hi-C, we show that loop domains form by a process that requires cohesin ATPases. Once formed, however, loops and compartments are maintained for hours without energy input. Strikingly, without ATP, we observe the emergence of hundreds of CTCF-independent loops that link regulatory DNA. We also identify architectural "stripes," where a loop anchor interacts with entire domains at high frequency. Stripes often tether super-enhancers to cognate promoters, and in B cells, they facilitate Igh transcription and recombination. Stripe anchors represent major hotspots for topoisomerase-mediated lesions, which promote chromosomal translocations and cancer. In plasmacytomas, stripes can deregulate Igh-translocated oncogenes. We propose that higher organisms have coopted cohesin extrusion to enhance transcription and recombination, with implications for tumor development.

Digestion-ligation-only Hi-C is an efficient and cost-effective method for chromosome conformation capture.

Chromosome conformation capture (3C) technologies can be used to investigate 3D genomic structures. However, high background noise, high costs, and a lack of straightforward noise evaluation in current methods impede the advancement of 3D genomic research. Here we developed a simple digestion-ligation-only Hi-C (DLO Hi-C) technology to explore the 3D landscape of the genome. This method requires only two rounds of digestion and ligation, without the need for biotin labeling and pulldown. Non-ligated DNA was efficiently removed in a cost-effective step by purifying specific linker-ligated DNA fragments. Notably, random ligation could be quickly evaluated in an early quality-control step before sequencing. Moreover, an in situ version of DLO Hi-C using a four-cutter restriction enzyme has been developed. We applied DLO Hi-C to delineate the genomic architecture of THP-1 and K562 cells and uncovered chromosomal translocations. This technology may facilitate investigation of genomic organization, gene regulation, and (meta)genome assembly.

ASXL1 interacts with the cohesin complex to maintain chromatid separation and gene expression for normal hematopoiesis.

ASXL1 is frequently mutated in a spectrum of myeloid malignancies with poor prognosis. Loss of Asxl1 leads to myelodysplastic syndrome-like disease in mice; however, the underlying molecular mechanisms remain unclear. We report that ASXL1 interacts with the cohesin complex, which has been shown to guide sister chromatid segregation and regulate gene expression. Loss of Asxl1 impairs the cohesin function, as reflected by an impaired telophase chromatid disjunction in hematopoietic cells. Chromatin immunoprecipitation followed by DNA sequencing data revealed that ASXL1, RAD21, and SMC1A share 93% of genomic binding sites at promoter regions in Lin-cKit+ (LK) cells. We have shown that loss of Asxl1 reduces the genome binding of RAD21 and SMC1A and alters the expression of ASXL1/cohesin target genes in LK cells. Our study underscores the ASXL1-cohesin interaction as a novel means to maintain normal sister chromatid separation and regulate gene expression in hematopoietic cells.

The genome sequence of Sea-Island cotton (Gossypium barbadense) provides insights into the allopolyploidization and development of superior spinnable fibres.

Gossypium hirsutum contributes the most production of cotton fibre, but G. barbadense is valued for its better comprehensive resistance and superior fibre properties. However, the allotetraploid genome of G. barbadense has not been comprehensively analysed. Here we present a high-quality assembly of the 2.57 gigabase genome of G. barbadense, including 80,876 protein-coding genes. The double-sized genome of the A (or At) (1.50 Gb) against D (or Dt) (853 Mb) primarily resulted from the expansion of Gypsy elements, including Peabody and Retrosat2 subclades in the Del clade, and the Athila subclade in the Athila/Tat clade. Substantial gene expansion and contraction were observed and rich homoeologous gene pairs with biased expression patterns were identified, suggesting abundant gene sub-functionalization occurred by allopolyploidization. More specifically, the CesA gene family has adapted differentially temporal expression patterns, suggesting an integrated regulatory mechanism of CesA genes from At and Dt subgenomes for the primary and secondary cellulose biosynthesis of cotton fibre in a "relay race"-like fashion. We anticipate that the G. barbadense genome sequence will advance our understanding the mechanism of genome polyploidization and underpin genome-wide comparison research in this genus.

Recurrent Fusion Genes in Gastric Cancer: CLDN18-ARHGAP26 Induces Loss of Epithelial Integrity.

Genome rearrangements, a hallmark of cancer, can result in gene fusions with oncogenic properties. Using DNA paired-end-tag (DNA-PET) whole-genome sequencing, we analyzed 15 gastric cancers (GCs) from Southeast Asians. Rearrangements were enriched in open chromatin and shaped by chromatin structure. We identified seven rearrangement hot spots and 136 gene fusions. In three out of 100 GC cases, we found recurrent fusions between CLDN18, a tight junction gene, and ARHGAP26, a gene encoding a RHOA inhibitor. Epithelial cell lines expressing CLDN18-ARHGAP26 displayed a dramatic loss of epithelial phenotype and long protrusions indicative of epithelial-mesenchymal transition (EMT). Fusion-positive cell lines showed impaired barrier properties, reduced cell-cell and cell-extracellular matrix adhesion, retarded wound healing, and inhibition of RHOA. Gain of invasion was seen in cancer cell lines expressing the fusion. Thus, CLDN18-ARHGAP26 mediates epithelial disintegration, possibly leading to stomach H(+) leakage, and the fusion might contribute to invasiveness once a cell is transformed.

B cell super-enhancers and regulatory clusters recruit AID tumorigenic activity.

The antibody gene mutator activation-induced cytidine deaminase (AID) promiscuously damages oncogenes, leading to chromosomal translocations and tumorigenesis. Why nonimmunoglobulin loci are susceptible to AID activity is unknown. Here, we study AID-mediated lesions in the context of nuclear architecture and the B cell regulome. We show that AID targets are not randomly distributed across the genome but are predominantly grouped within super-enhancers and regulatory clusters. Unexpectedly, in these domains, AID deaminates active promoters and eRNA(+) enhancers interconnected in some instances over megabases of linear chromatin. Using genome editing, we demonstrate that 3D-linked targets cooperate to recruit AID-mediated breaks. Furthermore, a comparison of hypermutation in mouse B cells, AID-induced kataegis in human lymphomas, and translocations in MEFs reveals that AID damages different genes in different cell types. Yet, in all cases, the targets are predominantly associated with topological complex, highly transcribed super-enhancers, demonstrating that these compartments are key mediators of AID recruitment.

Detection of chromosomal breakpoints in patients with developmental delay and speech disorders.

Delineating candidate genes at the chromosomal breakpoint regions in the apparently balanced chromosome rearrangements (ABCR) has been shown to be more effective with the emergence of next-generation sequencing (NGS) technologies. We employed a large-insert (7-11 kb) paired-end tag sequencing technology (DNA-PET) to systematically analyze genome of four patients harbouring cytogenetically defined ABCR with neurodevelopmental symptoms, including developmental delay (DD) and speech disorders. We characterized structural variants (SVs) specific to each individual, including those matching the chromosomal breakpoints. Refinement of these regions by Sanger sequencing resulted in the identification of five disrupted genes in three individuals: guanine nucleotide binding protein, q polypeptide (GNAQ), RNA-binding protein, fox-1 homolog (RBFOX3), unc-5 homolog D (C.elegans) (UNC5D), transmembrane protein 47 (TMEM47), and X-linked inhibitor of apoptosis (XIAP). Among them, XIAP is the causative gene for the immunodeficiency phenotype seen in the patient. The remaining genes displayed specific expression in the fetal brain and have known biologically relevant functions in brain development, suggesting putative candidate genes for neurodevelopmental phenotypes. This study demonstrates the application of NGS technologies in mapping individual gene disruptions in ABCR as a resource for deciphering candidate genes in human neurodevelopmental disorders (NDDs).

Dissecting the chromatin interactome of microRNA genes.

Our knowledge of the role of higher-order chromatin structures in transcription of microRNA genes (MIRs) is evolving rapidly. Here we investigate the effect of 3D architecture of chromatin on the transcriptional regulation of MIRs. We demonstrate that MIRs have transcriptional features that are similar to protein-coding genes. RNA polymerase II-associated ChIA-PET data reveal that many groups of MIRs and protein-coding genes are organized into functionally compartmentalized chromatin communities and undergo coordinated expression when their genomic loci are spatially colocated. We observe that MIRs display widespread communication in those transcriptionally active communities. Moreover, miRNA-target interactions are significantly enriched among communities with functional homogeneity while depleted from the same community from which they originated, suggesting MIRs coordinating function-related pathways at posttranscriptional level. Further investigation demonstrates the existence of spatial MIR-MIR chromatin interacting networks. We show that groups of spatially coordinated MIRs are frequently from the same family and involved in the same disease category. The spatial interaction network possesses both common and cell-specific subnetwork modules that result from the spatial organization of chromatin within different cell types. Together, our study unveils an entirely unexplored layer of MIR regulation throughout the human genome that links the spatial coordination of MIRs to their co-expression and function.

Whole-genome reconstruction and mutational signatures in gastric cancer.

BACKGROUND: Gastric cancer is the second highest cause of global cancer mortality. To explore the complete repertoire of somatic alterations in gastric cancer, we combined massively parallel short read and DNA paired-end tag sequencing to present the first whole-genome analysis of two gastric adenocarcinomas, one with chromosomal instability and the other with microsatellite instability. RESULTS: Integrative analysis and de novo assemblies revealed the architecture of a wild-type KRAS amplification, a common driver event in gastric cancer. We discovered three distinct mutational signatures in gastric cancer--against a genome-wide backdrop of oxidative and microsatellite instability-related mutational signatures, we identified the first exome-specific mutational signature. Further characterization of the impact of these signatures by combining sequencing data from 40 complete gastric cancer exomes and targeted screening of an additional 94 independent gastric tumors uncovered ACVR2A, RPL22 and LMAN1 as recurrently mutated genes in microsatellite instability-positive gastric cancer and PAPPA as a recurrently mutated gene in TP53 wild-type gastric cancer. CONCLUSIONS: These results highlight how whole-genome cancer sequencing can uncover information relevant to tissue-specific carcinogenesis that would otherwise be missed from exome-sequencing data.