Glioblastoma evolution and heterogeneity from a 3D whole-tumor perspective.

Treatment failure for the lethal brain tumor glioblastoma (GBM) is attributed to intratumoral heterogeneity and tumor evolution. We utilized 3D neuronavigation during surgical resection to acquire samples representing the whole tumor mapped by 3D spatial coordinates. Integrative tissue and single-cell analysis revealed sources of genomic, epigenomic, and microenvironmental intratumoral heterogeneity and their spatial patterning. By distinguishing tumor-wide molecular features from those with regional specificity, we inferred GBM evolutionary trajectories from neurodevelopmental lineage origins and initiating events such as chromothripsis to emergence of genetic subclones and spatially restricted activation of differential tumor and microenvironmental programs in the core, periphery, and contrast-enhancing regions. Our work depicts GBM evolution and heterogeneity from a 3D whole-tumor perspective, highlights potential therapeutic targets that might circumvent heterogeneity-related failures, and establishes an interactive platform enabling 360° visualization and analysis of 3D spatial patterns for user-selected genes, programs, and other features across whole GBM tumors.

Subtype-specific 3D genome alteration in acute myeloid leukaemia.

Acute myeloid leukaemia (AML) represents a set of heterogeneous myeloid malignancies, and hallmarks include mutations in epigenetic modifiers, transcription factors and kinases1-5. The extent to which mutations in AML drive alterations in chromatin 3D structure and contribute to myeloid transformation is unclear. Here we use Hi-C and whole-genome sequencing to analyse 25 samples from patients with AML and 7 samples from healthy donors. Recurrent and subtype-specific alterations in A/B compartments, topologically associating domains and chromatin loops were identified. RNA sequencing, ATAC with sequencing and CUT&Tag for CTCF, H3K27ac and H3K27me3 in the same AML samples also revealed extensive and recurrent AML-specific promoter-enhancer and promoter-silencer loops. We validated the role of repressive loops on their target genes by CRISPR deletion and interference. Structural variation-induced enhancer-hijacking and silencer-hijacking events were further identified in AML samples. Hijacked enhancers play a part in AML cell growth, as demonstrated by CRISPR screening, whereas hijacked silencers have a downregulating role, as evidenced by CRISPR-interference-mediated de-repression. Finally, whole-genome bisulfite sequencing of 20 AML and normal samples revealed the delicate relationship between DNA methylation, CTCF binding and 3D genome structure. Treatment of AML cells with a DNA hypomethylating agent and triple knockdown of DNMT1, DNMT3A and DNMT3B enabled the manipulation of DNA methylation to revert 3D genome organization and gene expression. Overall, this study provides a resource for leukaemia studies and highlights the role of repressive loops and hijacked cis elements in human diseases.

A map of cis-regulatory elements and 3D genome structures in zebrafish.

The zebrafish (Danio rerio) has been widely used in the study of human disease and development, and about 70% of the protein-coding genes are conserved between the two species1. However, studies in zebrafish remain constrained by the sparse annotation of functional control elements in the zebrafish genome. Here we performed RNA sequencing, assay for transposase-accessible chromatin using sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing, whole-genome bisulfite sequencing, and chromosome conformation capture (Hi-C) experiments in up to eleven adult and two embryonic tissues to generate a comprehensive map of transcriptomes, cis-regulatory elements, heterochromatin, methylomes and 3D genome organization in the zebrafish Tübingen reference strain. A comparison of zebrafish, human and mouse regulatory elements enabled the identification of both evolutionarily conserved and species-specific regulatory sequences and networks. We observed enrichment of evolutionary breakpoints at topologically associating domain boundaries, which were correlated with strong histone H3 lysine 4 trimethylation (H3K4me3) and CCCTC-binding factor (CTCF) signals. We performed single-cell ATAC-seq in zebrafish brain, which delineated 25 different clusters of cell types. By combining long-read DNA sequencing and Hi-C, we assembled the sex-determining chromosome 4 de novo. Overall, our work provides an additional epigenomic anchor for the functional annotation of vertebrate genomes and the study of evolutionarily conserved elements of 3D genome organization.

Genome-wide colocalization of RNA-DNA interactions and fusion RNA pairs.

Fusion transcripts are used as biomarkers in companion diagnoses. Although more than 15,000 fusion RNAs have been identified from diverse cancer types, few common features have been reported. Here, we compared 16,410 fusion transcripts detected in cancer (from a published cohort of 9,966 tumor samples of 33 cancer types) with genome-wide RNA-DNA interactions mapped in two normal, noncancerous cell types [using iMARGI, an enhanced version of the mapping of RNA-genome interactions (MARGI) assay]. Among the top 10 most significant RNA-DNA interactions in normal cells, 5 colocalized with the gene pairs that formed fusion RNAs in cancer. Furthermore, throughout the genome, the frequency of a gene pair to exhibit RNA-DNA interactions is positively correlated with the probability of this gene pair to present documented fusion transcripts in cancer. To test whether RNA-DNA interactions in normal cells are predictive of fusion RNAs, we analyzed these in a validation cohort of 96 lung cancer samples using RNA sequencing (RNA-seq). Thirty-seven of 42 fusion transcripts in the validation cohort were found to exhibit RNA-DNA interactions in normal cells. Finally, by combining RNA-seq, single-molecule RNA FISH, and DNA FISH, we detected a cancer sample with EML4-ALK fusion RNA without forming the EML4-ALK fusion gene. Collectively, these data suggest an RNA-poise model, where spatial proximity of RNA and DNA could poise for the creation of fusion transcripts.

Virological characteristics of SARS-CoV-2 Omicron BA.5.2.48.

With the prevalence of sequentially-emerged sublineages including BA.1, BA.2 and BA.5, SARS-CoV-2 Omicron infection has transformed into a regional epidemic disease. As a sublineage of BA.5, the BA.5.2.48 outbroke and evolved into multi-subvariants in China without clearly established virological characteristics. Here, we evaluated the virological characteristics of two isolates of the prevalent BA.5.2.48 subvariant, DY.2 and DY.1.1 (a subvariant of DY.1). Compared to the normal BA.5 spike, the double-mutated DY.1.1 spike demonstrates efficient cleavage, reduced fusogenicity and higher hACE2 binding affinity. BA.5.2.48 demonstrated enhanced airborne transmission capacity than BA.2 in hamsters. The pathogenicity of BA.5.2.48 is greater than BA.2, as revealed in Omicron-lethal H11-K18-hACE2 rodents. In both naïve and convalescent hamsters, DY.1.1 shows stronger fitness than DY.2 in hamster turbinates. Thus regional outbreaking of BA.5.2.48 promotes the multidirectional evolution of its subvariants, gaining either enhanced pathogenicity or a fitness in upper airways which is associated with higher transmission.

Molecular characterization and genome-wide mutations in porcine anal atresia candidate gene GLI2.

Anal atresia (i.e., anorectal malformations) is a severe disorder that occurs during the development of the distal hindgut in infants, swine, and many other mammals and has an unclear genetic background. Recently, the Shh-responsive transcription factor GLI2 has been shown as essential to the normal development of the hindgut, and QTL studies in pigs revealed that this gene may be an important candidate for anal atresia (AA). We used the pig as the model to study the contribution of GLI2 to AA. We revealed the genomic structure of the porcine GLI2 gene with 14 exons and obtained the porcine GLI2 mRNA sequence with a 4,656-bp ORF coding a 1,551-amino acid protein. We further scanned the genome-wide mutations in this gene by direct sequencing using three genomic DNA pools from the AA pigs, full-sibs of AA pigs, and unaffected pigs, respectively. Finally, 30 single nucleotide polymorphisms (SNPs) and one intronic 9-nucleotide (nt) deletion were identified. Of these SNPs, 23 are intronic, 6 are synonymous, and 1 (446 G>A) in exon 8 is nonsynonymous (365Met >Ile). NCOI-RFLP of the 446 G>A polymorphism suggested that the predominant genotypes were all GG and AG in the three pig groups. In addition, there was no significant difference among the three groups in allele frequencies, which demonstrated that this locus was not associated with AA in pigs. However, the 12 SNPs encompassing exon 4 to exon 8 showed strong linkage disequilibrium in the AA pigs, which indicated that the mutations somewhere in this region may contribute to AA in pigs. Therefore, further investigation in this region is needed to elucidate the underlying mutations involved in the porcine AA.

Fluorescence Polarization Assay Based on a New Recognition Motif QepA for the One-Step Detection of Fluoroquinolones in Eggs.

A recognition motif is vital in determining the specificity and sensitivity of the fluorescence polarization assay (FPA) for detecting chemical contaminants in food. Four candidates (Gyrase, GyrBA, TopIV, and QepA) were prepared for this study. The applicability of QepA was confirmed through DNA cleavage assay, inhibition effects, and mechanism investigations using molecular docking, compared to other counterparts. Finally, a novel FPA based on QepA and a CIP-FITC tracer for the detection of fluoroquinolones (FQs) in eggs was developed. The limits of detection (LODs) for eight fluoroquinolones ranged from 2.2 to 5.1 ng g-1, with enrofloxacin, danofloxacin, and difloxacin meeting the maximum residue limits (MRLs). The spiked recoveries ranged from 65.8 to 103.6% with coefficients of variation (CVs) of 5.4-12.8%. Therefore, a new recognition motif for FQs that did not belong to conventional antibodies was identified, and QepA-based FPA could be a potential tool for rapid, homogeneous, and sensitive monitoring of the residue of FQs in eggs.

Engineered MED12 mutations drive leiomyoma-like transcriptional and metabolic programs by altering the 3D genome compartmentalization.

Nearly 70% of Uterine fibroid (UF) tumors are driven by recurrent MED12 hotspot mutations. Unfortunately, no cellular models could be generated because the mutant cells have lower fitness in 2D culture conditions. To address this, we employ CRISPR to precisely engineer MED12 Gly44 mutations in UF-relevant myometrial smooth muscle cells. The engineered mutant cells recapitulate several UF-like cellular, transcriptional, and metabolic alterations, including altered Tryptophan/kynurenine metabolism. The aberrant gene expression program in the mutant cells is, in part, driven by a substantial 3D genome compartmentalization switch. At the cellular level, the mutant cells gain enhanced proliferation rates in 3D spheres and form larger lesions in vivo with elevated production of collagen and extracellular matrix deposition. These findings indicate that the engineered cellular model faithfully models key features of UF tumors and provides a platform for the broader scientific community to characterize genomics of recurrent MED12 mutations.

Enhancer Coamplification and Hijacking Promote Oncogene Expression in Liposarcoma.

SIGNIFICANCE: Comprehensive profiling of the enhancer landscape and 3D genome structure in liposarcoma identifies extensive enhancer-oncogene coamplification and enhancer hijacking events, deepening the understanding of how oncogenes are regulated in cancer.

Engineered MED12 mutations drive uterine fibroid-like transcriptional and metabolic programs by altering the 3D genome compartmentalization.

Uterine fibroid (UF) tumors originate from a mutated smooth muscle cell (SMC). Nearly 70% of these tumors are driven by hotspot recurrent somatic mutations in the MED12 gene; however, there are no tractable genetic models to study the biology of UF tumors because, under culture conditions, the non-mutant fibroblasts outgrow the mutant SMC cells, resulting in the conversion of the population to WT phenotype. The lack of faithful cellular models hampered our ability to delineate the molecular pathways downstream of MED12 mutations and identify therapeutics that may selectively target the mutant cells. To overcome this challenge, we employed CRISPR knock-in with a sensitive PCR-based screening strategy to precisely engineer cells with mutant MED12 Gly44, which constitutes 50% of MED12 exon two mutations. Critically, the engineered myometrial SMC cells recapitulate several UF-like cellular, transcriptional and metabolic alterations, including enhanced proliferation rates in 3D spheres and altered Tryptophan/kynurenine metabolism. Our transcriptomic analysis supported by DNA synthesis tracking reveals that MED12 mutant cells, like UF tumors, have heightened expression of DNA repair genes but reduced DNA synthesis rates. Consequently, these cells accumulate significantly higher rates of DNA damage and are selectively more sensitive to common DNA-damaging chemotherapy, indicating mutation-specific and therapeutically relevant vulnerabilities. Our high-resolution 3D chromatin interaction analysis demonstrates that the engineered MED12 mutations drive aberrant genomic activity due to a genome-wide chromatin compartmentalization switch. These findings indicate that the engineered cellular model faithfully models key features of UF tumors and provides a novel platform for the broader scientific community to characterize genomics of recurrent MED12 mutations and discover potential therapeutic targets.

Breakage fusion bridge cycles drive high oncogene copy number, but not intratumoral genetic heterogeneity or rapid cancer genome change.

Oncogene amplification is a major driver of cancer pathogenesis. Breakage fusion bridge (BFB) cycles, like extrachromosomal DNA (ecDNA), can lead to high copy numbers of oncogenes, but their impact on intratumoral heterogeneity, treatment response, and patient survival are not well understood due to difficulty in detecting them by DNA sequencing. We describe a novel algorithm that detects and reconstructs BFB amplifications using optical genome maps (OGMs), called OM2BFB. OM2BFB showed high precision (>93%) and recall (92%) in detecting BFB amplifications in cancer cell lines, PDX models and primary tumors. OM-based comparisons demonstrated that short-read BFB detection using our AmpliconSuite (AS) toolkit also achieved high precision, albeit with reduced sensitivity. We detected 371 BFB events using whole genome sequences from 2,557 primary tumors and cancer lines. BFB amplifications were preferentially found in cervical, head and neck, lung, and esophageal cancers, but rarely in brain cancers. BFB amplified genes show lower variance of gene expression, with fewer options for regulatory rewiring relative to ecDNA amplified genes. BFB positive (BFB (+)) tumors showed reduced heterogeneity of amplicon structures, and delayed onset of resistance, relative to ecDNA(+) tumors. EcDNA and BFB amplifications represent contrasting mechanisms to increase the copy numbers of oncogene with markedly different characteristics that suggest different routes for intervention.

Subtype-associated epigenomic landscape and 3D genome structure in bladder cancer.

Muscle-invasive bladder cancers are characterized by their distinct expression of luminal and basal genes, which could be used to predict key clinical features such as disease progression and overall survival. Transcriptionally, FOXA1, GATA3, and PPARG are shown to be essential for luminal subtype-specific gene regulation and subtype switching, while TP63, STAT3, and TFAP2 family members are critical for regulation of basal subtype-specific genes. Despite these advances, the underlying epigenetic mechanisms and 3D chromatin architecture responsible for subtype-specific regulation in bladder cancer remain unknown. RESULT: We determine the genome-wide transcriptome, enhancer landscape, and transcription factor binding profiles of FOXA1 and GATA3 in luminal and basal subtypes of bladder cancer. Furthermore, we report the first-ever mapping of genome-wide chromatin interactions by Hi-C in both bladder cancer cell lines and primary patient tumors. We show that subtype-specific transcription is accompanied by specific open chromatin and epigenomic marks, at least partially driven by distinct transcription factor binding at distal enhancers of luminal and basal bladder cancers. Finally, we identify a novel clinically relevant transcription factor, Neuronal PAS Domain Protein 2 (NPAS2), in luminal bladder cancers that regulates other subtype-specific genes and influences cancer cell proliferation and migration. CONCLUSION: In summary, our work identifies unique epigenomic signatures and 3D genome structures in luminal and basal urinary bladder cancers and suggests a novel link between the circadian transcription factor NPAS2 and a clinical bladder cancer subtype.

Function and regulation mechanism of Chk1 during meiotic maturation in porcine oocytes.

Checkpoint 1 (Chk1), as an important member of DNA replication checkpoint and DNA damage response, has an important role during the G2/M stage of mitosis. In this study, we used porcine oocyte as a model to investigate the function of Chk1 during porcine oocyte maturation. Chk1 was expressed from germinal vesicle (GV) to metaphase II (MII) stages, mainly localized in the cytoplasm at GV stage and moved to the spindle after germinal vesicle breakdown (GVBD). Chk1 depletion not only induced oocytes to be arrested at MI stage with abnormal chromosomes arrangement, but also inhibited the degradation of Cyclin B1 and decreased the expression of Mitotic Arrest Deficient 2-Like 1 (Mad2L1), one of spindle assembly checkpoint (SAC) proteins, and cadherin 1 (Cdh1), one of coactivation for anaphase-promoting complex/cyclosome (APC/C). Moreover, Chk1 overexpression delayed GVBD. These results demonstrated that Chk1 facilitated the timely degradation of Cyclin B1 at anaphase I (AI) and maintained the expression of Mad2L1 and Cdh1, which ensured that all chromosomes were accurately located in a line, and then oocytes passed metaphase I (MI) and AI and exited from the first meiotic division successfully. In addition, we proved that Chk1 had not function on GVBD of porcine oocytes, which suggested that maturation of porcine oocytes did not need the DNA damage checkpoint, which was different from the mouse oocyte maturation.