Decoding heterogeneous and coordinated tissue architecture in glioblastoma using spatial transcriptomics.

Glioblastoma multiforme (GBM) is one of the most lethal brain tumors, characterized by profound heterogeneity. While single-cell transcriptomic studies have revealed extensive intra-tumor heterogeneity, shed light on intra-tumor diversity, spatial intricacies remain largely unexplored. Leveraging clinical GBM specimens, this study employs spatial transcriptomics technology to delve into gene expression heterogeneity. Our investigation unveils a significant enrichment of tissue stem cell signature in regions bordering necrosis and the peritumoral area, positively correlated with the mesenchymal subtype signature. Moreover, upregulated genes in these regions are linked with extracellular matrix (ECM)-receptor interaction, proteoglycans, as well as vascular endothelial growth factor (VEGF) and angiopoietin-Tie (ANGPT) signaling pathways. In contrast, signatures related to glycogen metabolism and oxidative phosphorylation show no relevance to pathological zoning, whereas creatine metabolism signature is notably exclusive to vascular-enriched areas. These spatial profiles not only offer valuable references but also pave the way for future in-depth functional and mechanistic investigations into GBM progression.

Landscape of enhancer disruption and functional screen in melanoma cells.

BACKGROUND: The high mutation rate throughout the entire melanoma genome presents a major challenge in stratifying true driver events from the background mutations. Numerous recurrent non-coding alterations, such as those in enhancers, can shape tumor evolution, thereby emphasizing the importance in systematically deciphering enhancer disruptions in melanoma. RESULTS: Here, we leveraged 297 melanoma whole-genome sequencing samples to prioritize highly recurrent regions. By performing a genome-scale CRISPR interference (CRISPRi) screen on highly recurrent region-associated enhancers in melanoma cells, we identified 66 significant hits which could have tumor-suppressive roles. These functional enhancers show unique mutational patterns independent of classical significantly mutated genes in melanoma. Target gene analysis for the essential enhancers reveal many known and hidden mechanisms underlying melanoma growth. Utilizing extensive functional validation experiments, we demonstrate that a super enhancer element could modulate melanoma cell proliferation by targeting MEF2A, and another distal enhancer is able to sustain PTEN tumor-suppressive potential via long-range interactions. CONCLUSIONS: Our study establishes a catalogue of crucial enhancers and their target genes in melanoma growth and progression, and illuminates the identification of novel mechanisms of dysregulation for melanoma driver genes and new therapeutic targeting strategies.

Single-Cell DNA Sequencing Reveals Punctuated and Gradual Clonal Evolution in Hepatocellular Carcinoma.

BACKGROUND & AIMS: Copy number alterations (CNAs), elicited by genome instability, are a major source of intratumor heterogeneity. How CNAs evolve in hepatocellular carcinoma (HCC) remains unknown. METHODS: We performed single-cell DNA sequencing (scDNA-seq) on 1275 cells isolated from 10 patients with HCC, ploidy-resolved scDNA-seq on 356 cells from 1 additional patient, and single-cell RNA sequencing on 27,344 cells from 3 additional patients. Three statistical fitting models were compared to investigate the CNA accumulation pattern. RESULTS: Cells in the tumor were categorized into the following 3 subpopulations: euploid, pseudoeuploid, and aneuploid. Our scDNA-seq analysis revealed that CNA accumulation followed a dual-phase copy number evolution model, that is, a punctuated phase followed by a gradual phase. Patients who exhibited prolonged gradual phase showed higher intratumor heterogeneity and worse disease-free survival. Integrating bulk RNA sequencing of 17 patients with HCC, published datasets of 1196 liver tumors, and immunohistochemical staining of 202 HCC tumors, we found that high expression of CAD, a gene involved in pyrimidine synthesis, was correlated with rapid tumorigenesis and reduced survival. The dual-phase copy number evolution model was validated by our single-cell RNA sequencing data and published scDNA-seq datasets of other cancer types. Furthermore, ploidy-resolved scDNA-seq revealed the common clonal origin of diploid- and polyploid-aneuploid cells, suggesting that polyploid tumor cells were generated by whole genome doubling of diploid tumor cells. CONCLUSIONS: Our work revealed a novel dual-phase copy number evolution model, showed HCC with longer gradual phase was more severe, identified CAD as a promising biomarker for early recurrence of HCC, and supported the diploid origin of polyploid HCC.

Bimodal regulation of the PRC2 complex by USP7 underlies tumorigenesis.

Although overexpression of EZH2, a catalytic subunit of the polycomb repressive complex 2 (PRC2), is an eminent feature of various cancers, the regulation of its abundance and function remains insufficiently understood. We report here that the PRC2 complex is physically associated with ubiquitin-specific protease USP7 in cancer cells where USP7 acts to deubiquitinate and stabilize EZH2. Interestingly, we found that USP7-catalyzed H2BK120ub1 deubiquitination is a prerequisite for chromatin loading of PRC2 thus H3K27 trimethylation, and this process is not affected by H2AK119 ubiquitination catalyzed by PRC1. Genome-wide analysis of the transcriptional targets of the USP7/PRC2 complex identified a cohort of genes including FOXO1 that are involved in cell growth and proliferation. We demonstrated that the USP7/PRC2 complex drives cancer cell proliferation and tumorigenesis in vitro and in vivo. We showed that the expression of both USP7 and EZH2 elevates during tumor progression, corresponding to a diminished FOXO1 expression, and the level of the expression of USP7 and EZH2 strongly correlates with histological grades and prognosis of tumor patients. These results reveal a dual role for USP7 in the regulation of the abundance and function of EZH2, supporting the pursuit of USP7 as a therapeutic target for cancer intervention.

Ultrafast and scalable variant annotation and prioritization with big functional genomics data.

The advances of large-scale genomics studies have enabled compilation of cell type-specific, genome-wide DNA functional elements at high resolution. With the growing volume of functional annotation data and sequencing variants, existing variant annotation algorithms lack the efficiency and scalability to process big genomic data, particularly when annotating whole-genome sequencing variants against a huge database with billions of genomic features. Here, we develop VarNote to rapidly annotate genome-scale variants in large and complex functional annotation resources. Equipped with a novel index system and a parallel random-sweep searching algorithm, VarNote shows substantial performance improvements (two to three orders of magnitude) over existing algorithms at different scales. It supports both region-based and allele-specific annotations and introduces advanced functions for the flexible extraction of annotations. By integrating massive base-wise and context-dependent annotations in the VarNote framework, we introduce three efficient and accurate pipelines to prioritize the causal regulatory variants for common diseases, Mendelian disorders, and cancers.

PHF20L1 as a H3K27me2 reader coordinates with transcriptional repressors to promote breast tumorigenesis.

TUDOR domain-containing proteins (TDRDs) are chiefly responsible for recognizing methyl-lysine/arginine residue. However, how TDRD dysregulation contributes to breast tumorigenesis is poorly understood. Here, we report that TUDOR domain-containing PHF20L1 as a H3K27me2 reader exerts transcriptional repression by recruiting polycomb repressive complex 2 (PRC2) and Mi-2/nucleosome remodeling and deacetylase (NuRD) complex, linking PRC2-mediated methylation and NuRD-mediated deacetylation of H3K27. Furthermore, PHF20L1 was found to serve as a potential MYC and hypoxia-driven oncogene, promoting glycolysis, proliferation, and metastasis of breast cancer cells by directly inhibiting tumor suppressors such as HIC1, KISS1, and BRCA1. PHF20L1 expression was also strongly correlated with higher histologic grades of breast cancer and markedly up-regulated in several cancers. Meanwhile, Phf20l1 deletion not only induces growth retardation and mammary ductal outgrowth delay but also inhibits tumorigenesis in vivo. Our data indicate that PHF20L1 promotes tumorigenesis, supporting the pursuit of PHF20L1 as a target for cancer therapy.

The novel miR-1269b-regulated protein SVEP1 induces hepatocellular carcinoma proliferation and metastasis likely through the PI3K/Akt pathway.

Decreased intercellular adhesion is a key step in the metastasis and recurrence of many cancers, including hepatocellular carcinoma (HCC). SVEP1 is an important cell adhesion molecule that plays a key role in regulating intercellular adhesion and embryonic lymphatic development. However, the expression patterns and roles of SVEP1 in HCC are still largely unknown. We identified SVEP1 expression by analyzing 220 HCC samples from our cancer center. TCGA and GEO online-databases were used for data calibration and validation. SVEP1 was differentially expressed in two groups of HCCs with different risks of recurrence and was deemed as an independent risk factor for the prognosis of HCC. The expression of SVEP1 is negatively related to the proliferation and metastasis of HCC. Downregulation of SVEP1 expression promoted in vitro HCC cell migration, chemotaxis, invasion and proliferation, as well as in vivo tumor growth, local invasion and metastasis in a mouse model. Bioinformatic analysis and RT-PCR results showed that miR-1269b expression is negatively correlated with the SVEP1 expression and the prognosis of HCC patients. Further experiments showed that miR-1269b directly targets and downregulates the expression of SVEP1, which further induces the phosphorylation of Akt at thr308. These regulatory effects ultimately mediate the proliferation and metastasis of HCC cells. SVEP1 could serve as a promising prognostic marker of HCC. MiR-1269b downregulates SVEP1 expression and promotes HCC proliferation and metastasis likely through the PI3k/Akt signaling pathway.

HRP2-DPF3a-BAF complex coordinates histone modification and chromatin remodeling to regulate myogenic gene transcription.

Functional crosstalk between histone modifications and chromatin remodeling has emerged as a key regulatory mode of transcriptional control during cell fate decisions, but the underlying mechanisms are not fully understood. Here we discover an HRP2-DPF3a-BAF epigenetic pathway that coordinates methylated histone H3 lysine 36 (H3K36me) and ATP-dependent chromatin remodeling to regulate chromatin dynamics and gene transcription during myogenic differentiation. Using siRNA screening targeting epigenetic modifiers, we identify hepatoma-derived growth factor-related protein 2 (HRP2) as a key regulator of myogenesis. Knockout of HRP2 in mice leads to impaired muscle regeneration. Mechanistically, through its HIV integrase binding domain (IBD), HRP2 associates with the BRG1/BRM-associated factor (BAF) chromatin remodeling complex by interacting directly with the BAF45c (DPF3a) subunit. Through its Pro-Trp-Trp-Pro (PWWP) domain, HRP2 preferentially binds to H3K36me2. Consistent with the biochemical studies, ChIP-seq analyses show that HRP2 colocalizes with DPF3a across the genome and that the recruitment of HRP2/DPF3a to chromatin is dependent on H3K36me2. Integrative transcriptomic and cistromic analyses, coupled with ATAC-seq, reveal that HRP2 and DPF3a activate myogenic genes by increasing chromatin accessibility through recruitment of BRG1, the ATPase subunit of the BAF complex. Taken together, these results illuminate a key role for the HRP2-DPF3a-BAF complex in the epigenetic coordination of gene transcription during myogenic differentiation.

The role of bone marrow-derived cells in the origin of liver cancer revealed by single-cell sequencing.

Objective: Epithelial cancers often originate from progenitor cells, while the origin of hepatocellular carcinoma (HCC) is still controversial. HCC, one of the deadliest cancers, is closely linked with liver injuries and chronic inflammation, which trigger massive infiltration of bone marrow-derived cells (BMDCs) during liver repair. Methods: To address the possible roles of BMDCs in HCC origination, we established a diethylnitrosamine (DEN)-induced HCC model in bone marrow transplanted mice. Immunohistochemistry and frozen tissue immunofluorescence were used to verify DEN-induced HCC in the pathology of the disease. The cellular origin of DEN-induced HCC was further studied by single cell sequencing, single-cell nested PCR, and immunofluorescence-fluorescence in situ hybridization. Results: Studies by using single cell sequencing and biochemical analysis revealed that HCC cells in these mice were coming from donor mice BMDCs, and not from recipient mice. Furthermore, the copy numbers of mouse orthologs of several HCC-related genes previously reported in human HCC were also altered in our mouse model. DEN-induced HCCs exhibited a similar histological phenotype and genomic profile as human HCCs. Conclusions: These results suggested that BMDCs are an important origin of HCC, which provide important clues to HCC prevention, detection, and treatments.

Diversity spectrum analysis identifies mutation-specific effects of cancer driver genes.

Mutation-specific effects of cancer driver genes influence drug responses and the success of clinical trials. We reasoned that these effects could unbalance the distribution of each mutation across different cancer types, as a result, the cancer preference can be used to distinguish the effects of the causal mutation. Here, we developed a network-based framework to systematically measure cancer diversity for each driver mutation. We found that half of the driver genes harbor cancer type-specific and pancancer mutations simultaneously, suggesting that the pervasive functional heterogeneity of the mutations from even the same driver gene. We further demonstrated that the specificity of the mutations could influence patient drug responses. Moreover, we observed that diversity was generally increased in advanced tumors. Finally, we scanned potentially novel cancer driver genes based on the diversity spectrum. Diversity spectrum analysis provides a new approach to define driver mutations and optimize off-label clinical trials.

regBase: whole genome base-wise aggregation and functional prediction for human non-coding regulatory variants.

Predicting the functional or pathogenic regulatory variants in the human non-coding genome facilitates the interpretation of disease causation. While numerous prediction methods are available, their performance is inconsistent or restricted to specific tasks, which raises the demand of developing comprehensive integration for those methods. Here, we compile whole genome base-wise aggregations, regBase, that incorporate largest prediction scores. Building on different assumptions of causality, we train three composite models to score functional, pathogenic and cancer driver non-coding regulatory variants respectively. We demonstrate the superior and stable performance of our models using independent benchmarks and show great success to fine-map causal regulatory variants on specific locus or at base-wise resolution. We believe that regBase database together with three composite models will be useful in different areas of human genetic studies, such as annotation-based casual variant fine-mapping, pathogenic variant discovery as well as cancer driver mutation identification. regBase is freely available at https://github.com/mulinlab/regBase.

Enterohemorrhagic Escherichia coli Tir inhibits TAK1 activation and mediates immune evasion.

Many pathogens infect hosts through various immune evasion strategies. However, the molecular mechanisms by which pathogen proteins modulate and evade the host immune response remain unclear. Enterohemorrhagic Escherichia coli (EHEC) is a pathological strain that can induce mitogen-activated protein (MAP) kinase (Erk, Jnk and p38 MAPK) and NF-κB pathway activation and proinflammatory cytokine production, which then causes diarrheal diseases such as hemorrhagic colitis and hemolytic uremic syndrome. Transforming growth factor ß-activated kinase-1 (TAK1) is a key regulator involved in distinct innate immune signalling pathways. Here we report that EHEC translocated intimin receptor (Tir) protein inhibits the expression of EHEC-induced proinflammatory cytokines by interacting with the host tyrosine phosphatase SHP-1, which is dependent on the phosphorylation of immunoreceptor tyrosine-based inhibition motifs (ITIMs). Mechanistically, the association of EHEC Tir with SHP-1 facilitated the recruitment of SHP-1 to TAK1 and inhibited TAK1 phosphorylation, which then negatively regulated K63-linked polyubiquitination of TAK1 and downstream signal transduction. Taken together, these results suggest that EHEC Tir negatively regulates proinflammatory responses by inhibiting the activation of TAK1, which is essential for immune evasion and could be a potential target for the treatment of bacterial infection.

Targeting Super-Enhancer-Driven Oncogenic Transcription by CDK7 Inhibition in Anaplastic Thyroid Carcinoma.

Background: Anaplastic thyroid carcinoma (ATC) is one of the most aggressive malignancies, with no effective treatment currently available. The molecular mechanisms of ATC carcinogenesis remain poorly understood. The objective of this study was to investigate the mechanisms and functions of super-enhancer (SE)-driven oncogenic transcriptional addiction in the progression of ATC and identify new drug targets for ATC treatments. Methods: High-throughput chemical screening was performed to identify new drugs inhibiting ATC cell growth. Cell viability assay, colony formation analysis, cell-cycle analysis, and animal study were used to examine the effects of drug treatments on ATC progression. Chromatin immunoprecipitation sequencing was conducted to establish a SE landscape of ATC. Integrative analysis of RNA sequencing, chromatin immunoprecipitation sequencing, and CRISPR/Cas9-mediated gene editing was used to identify THZ1 target genes. Drug combination analysis was performed to assess drug synergy. Patient samples were analyzed to evaluate candidate biomarkers of prognosis in ATC. Results: THZ1, a covalent inhibitor of cyclin-dependent kinase 7 (CDK7), was identified as a potent anti-ATC compound by high-throughput chemical screening. ATC cells, but not papillary thyroid carcinoma cells, are exceptionally sensitive to CDK7 inhibition. An integrative analysis of both gene expression profiles and SE features revealed that the SE-mediated oncogenic transcriptional amplification mediates the vulnerability of ATC cells to THZ1 treatment. Combining this integrative analysis with functional assays led to the discovery of a number of novel cancer genes of ATC, including PPP1R15A, SMG9, and KLF2. Inhibition of PPP1R15A with Guanabenz or Sephin1 greatly suppresses ATC growth. Significantly, the expression level of PPP1R15A is correlated with CDK7 expression in ATC tissue samples. Elevated expression of PPP1R15A and CDK7 are both associated with poor clinical prognosis in ATC patients. Importantly, CDK7 or PPP1R15A inhibition sensitizes ATC cells to conventional chemotherapy. Conclusions: Taken together, these findings demonstrate transcriptional addiction in ATC pathobiology and identify CDK7 and PPP1R15A as potential biomarkers and therapeutic targets for ATC.

Downregulation of glycine decarboxylase enhanced cofilin-mediated migration in hepatocellular carcinoma cells.

Metabolic reprogramming is a hallmark of cancer. Glycine decarboxylase (GLDC), an oxidoreductase, plays an important role in amino acid metabolism. While GLDC promotes tumor initiation and proliferation in non-small cell lung cancer and glioma and it was reported as a putative tumor suppressor gene in gastric cancer, the role of GLDC in hepatocellular carcinoma (HCC) is unknown. In the current study, microarray-based analysis suggested that GLDC expression was low in highly malignant HCC cell lines, and clinicopathological analysis revealed a decrease in GLDC in HCC tumor samples. While the knockdown of GLDC enhanced cancer cell migration and invasion, GLDC overexpression inhibited them. Mechanistic studies revealed that GLDC knockdown increased the levels of reactive oxygen species (ROS) and decreased the ratio of glutathione/oxidized glutathione (GSH/GSSG), which in turn dampened the ubiquitination of cofilin, a key regulator of actin polymerization. Consequently, the protein level of cofilin was elevated, which accounted for the increase in cell migration. The overexpression of GLDC reversed the phenotype. Treatment with N-acetyl-L-cysteine decreased the protein level of cofilin while treatment with H2O2 increased it, further confirming the role of ROS in regulating cofilin degradation. In a tumor xenographic transplant nude mouse model, the knockdown of GLDC promoted intrahepatic metastasis of HCC while GLDC overexpression inhibited it. Our data indicate that GLDC downregulation decreases ROS-mediated ubiquitination of cofilin to enhance HCC progression and intrahepatic metastasis.

Multiomics Integration Reveals the Landscape of Prometastasis Metabolism in Hepatocellular Carcinoma.

The systematic investigation of gene mutation and expression is important to discover novel biomarkers and therapeutic targets in cancers. Here, we integrated genomics, transcriptomics, proteomics, and metabolomics to analyze three hepatocellular carcinoma (HCC) cell lines with differential metastatic potentials. The results revealed the profile of the prometastasis metabolism potentially associated with HCC metastasis. The multiomic analysis identified 12 genes with variations at multiple levels from three metabolic pathways, including glycolysis, starch, and sucrose metabolism, and glutathione metabolism. Furthermore, uridine diphosphate (UDP)-glucose pyrophosphorylase 2 (UGP2), was observed to be persistently up-regulated with increased metastatic potential. UGP2 overexpression promoted cell migration and invasion and enhanced glycogenesis in vitro The role of UGP2 in metastasis was further confirmed using a tumor xenograft mouse model. Taken together, the compendium of multiomic data provides valuable insights in understanding the roles of shifted cellular metabolism in HCC metastasis.

Body mass index and risk of lung cancer: Systematic review and dose-response meta-analysis.

Questions remain about the significance of the dose-response relationship between body mass index (BMI) and lung cancer (LC) risk. Pertinent studies were identified through a search in EMBASE and PUBMED from July 2014 until March 2015. The summary relative risk (SRR) and confidence interval (CI) were estimated. The dose-response relationship was assessed using a restricted cubic spline. The overall meta-analysis showed evidence of a nonlinear association between BMI and LC risk (Pnonlinearity < 0.001). The SRR were 0.98 (95%CI: 0.95-1.01) for 25 kg/m(2), 0.91 (95%CI: 0.85-0.98) for 30 kg/m(2) and 0.81 (95% CI: 0.72-0.91) for 35 kg/m(2), with mild between-study heterogeneity (I(2) = 5%). The results of the stratified analysis by gender were comparable to those of the overall meta-analysis. When stratified by smoking status, linear dose-response associations were observed for current smokers, ex-smokers and non-smokers (Pnonlinearity > 0.05), whereas the effects were attenuated when restricting analysis to non-smokers, and at the point of 30 kg/m(2), the SRR was 0.96 (95%CI: 0.86-1.07) for males and 0.95 (95%CI: 0.89-1.02) for females. This meta-analysis provides quantitative evidence that increasing BMI is a protective factor against LC. Keeping normal-to-moderate BMI should be prescribed as an evidence-based lifestyle tip for LC prevention in smokers.

Exome sequencing identifies frequent mutation of MLL2 in non-small cell lung carcinoma from Chinese patients.

Lung cancer is the most common cause of cancer mortality worldwide, with an estimated 1.4 million deaths each year. Here we report whole-exome sequencing of nine tumor/normal tissue pairs from Chinese patients with non-small cell lung carcinoma (NSCLC). This allows us to identify a number of significantly mutated genes in NSCLC, which were highly enriched in DNA damage repair, NF-κB pathway, JAK/STAT signaling and chromatin modification. Notably, we identify a histone-lysine methyltransferase gene, namely, MLL2, as one of the most significantly mutated genes in our screen. In a following validation study, we identify deleterious mutations of MLL2 in 12 out of 105 (11.4%) NSCLC patients. Additionally, reduced or lost expression of MLL2 was commonly observed in tumor tissues as compared with paired adjacent non-tumor tissues regardless of mutation status. Together, our study defines the landscape of somatic mutations in Chinese NSCLC and supports the role of MLL2 mutation in the pathogenesis of the disease.

Combined analysis with copy number variation identifies risk loci in lung cancer.

BACKGROUND: Lung cancer is the most important cause of cancer mortality worldwide, but the underlying mechanisms of this disease are not fully understood. Copy number variations (CNVs) are promising genetic variations to study because of their potential effects on cancer. METHODOLOGY/PRINCIPAL FINDINGS: Here we conducted a pilot study in which we systematically analyzed the association of CNVs in two lung cancer datasets: the Environment And Genetics in Lung cancer Etiology (EAGLE) and the Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial datasets. We used a preestablished association method to test the datasets separately and conducted a combined analysis to test the association accordance between the two datasets. Finally, we identified 167 risk SNP loci and 22 CNVs associated with lung cancer and linked them with recombination hotspots. Functional annotation and biological relevance analyses implied that some of our predicted risk loci were supported by other studies and might be potential candidate loci for lung cancer studies. CONCLUSIONS/SIGNIFICANCE: Our results further emphasized the importance of copy number variations in cancer and might be a valuable complement to current genome-wide association studies on cancer.