Search

Only a handful of somatic alterations have been linked to endocrine therapy resistance in hormone-dependent breast cancer, potentially explaining ∼40% of relapses. If other mechanisms underlie the evolution of hormone-dependent breast cancer under adjuvant therapy is currently unknown. In this work, we employ functional genomics to dissect the contribution of cis-regulatory elements (CRE) to cancer evolution by focusing on 12 megabases of noncoding DNA, including clonal enhancers, gene promoters, and boundaries of topologically associating domains. Parallel epigenetic perturbation (CRISPRi) in vitro reveals context-dependent roles for many of these CREs, with a specific impact on dormancy entrance and endocrine therapy resistance. Profiling of CRE somatic alterations in a unique, longitudinal cohort of patients treated with endocrine therapies identifies a limited set of noncoding changes potentially involved in therapy resistance. Overall, our data uncover how endocrine therapies trigger the emergence of transient features which could ultimately be exploited to hinder the adaptive process. Significance: This study shows that cells adapting to endocrine therapies undergo changes in the usage or regulatory regions. Dormant cells are less vulnerable to regulatory perturbation but gain transient dependencies which can be exploited to decrease the formation of dormant persisters.

Activated Wnt/β-catenin pathway is a key genetic event in liver cancer development. Solute carrier (SLC) transporters are promising drug targets. Here, we identify SLC13A3 as a drug-targetable effector downstream of β-catenin in liver cancer. SLC13A3 expression is elevated in human liver cancer samples with gain of function (GOF) mutant CTNNB1, the gene encoding β-catenin. Activation of β-catenin up-regulates SLC13A3, leading to intracellular accumulation of endogenous SLC13A3 substrates. SLC13A3 is identified as a low-affinity transporter for glutathione (GSH). Silencing of SLC13A3 downregulates the leucine transporter SLC7A5 via c-MYC signaling, leading to leucine depletion and mTOR inactivation. Furthermore, silencing of SLC13A3 depletes GSH and induces autophagic ferroptosis in β-catenin-activated liver cancer cells. Importantly, both genetic inhibition of SLC13A3 and a small molecule SLC13A3 inhibitor suppress β-catenin-driven hepatocarcinogenesis in mice. Altogether, our study suggests that SLC13A3 could be a promising therapeutic target for treating human liver cancers with GOF CTNNB1 mutations.

EWS fusion oncoproteins underlie several human malignancies including Desmoplastic Small Round Cell Tumor (DSRCT), an aggressive cancer driven by EWS-WT1 fusion proteins. Here we combine chromatin occupancy and 3D profiles to identify EWS-WT1-dependent gene regulation networks and target genes. We show that EWS-WT1 is a powerful chromatin activator controlling an oncogenic gene expression program that characterizes primary tumors. Similar to wild type WT1, EWS-WT1 has two isoforms that differ in their DNA binding domain and we find that they have distinct DNA binding profiles and are both required to generate viable tumors that resemble primary DSRCT. Finally, we identify candidate EWS-WT1 target genes with potential therapeutic implications, including CCND1, whose inhibition by the clinically-approved drug Palbociclib leads to marked tumor burden decrease in DSRCT PDXs in vivo. Taken together, our studies identify gene regulation programs and therapeutic vulnerabilities in DSRCT and provide a mechanistic understanding of the complex oncogenic activity of EWS-WT1.

Resistance to the current Androgen Receptor Signaling Inhibitor (ARSI) therapies has led to higher incidences of therapy-induced neuroendocrine-like prostate cancer (t-NEPC). This highly aggressive subtype with predominant small-cell-like characteristics is resistant to taxane chemotherapies and has a dismal overall survival. t-NEPCs are mostly treated with platinum-based drugs with a combination of etoposide or taxane and have less selectivity and high systemic toxicity, which often limit their clinical potential. During t-NEPC transformation, adenocarcinomas lose their luminal features and adopt neuro-basal characteristics. Whether the adaptive neuronal characteristics of t-NEPC are responsible for such taxane resistance remains unknown. Pathway analysis from patient gene-expression databases indicates that t-NEPC upregulates various neuronal pathways associated with enhanced cellular networks. To identify transcription factor(s) (TF) that could be important for promoting the gene expression for neuronal characters in t-NEPC, we performed ATAC-Seq, acetylated-histone ChIP-seq, and RNA-seq in our NE-like cell line models and analyzed the promoters of transcriptionally active and significantly enriched neuroendocrine-like (NE-like) cancer-specific genes. Our results indicate that Pax5 could be an important transcription factor for neuronal gene expression and specific to t-NEPC. Pathway analysis revealed that Pax5 expression is involved in axonal guidance, neurotransmitter regulation, and neuronal adhesion, which are critical for strong cellular communications. Further results suggest that depletion of Pax5 disrupts neurite-mediated cellular communication in NE-like cells and reduces surface growth factor receptor activation, thereby, sensitizing them to docetaxel therapies. Moreover, t-NEPC-specific hydroxymethylation of Pax5 promoter CpG islands favors Pbx1 binding to induce Pax5 expression. Based on our study, we concluded that continuous exposure to ARSI therapies leads to epigenetic modifications and Pax5 activation in t-NEPC, which promotes the expression of genes necessary to adopt taxane-resistant NE-like cancer. Thus, targeting the Pax5 axis can be beneficial for reverting their taxane sensitivity.

PURPOSE: The purpose of this study was to evaluate RB1 expression and survival across ovarian carcinoma histotypes and how co-occurrence of BRCA1 or BRCA2 (BRCA) alterations and RB1 loss influences survival in tubo-ovarian high-grade serous carcinoma (HGSC). EXPERIMENTAL DESIGN: RB1 protein expression was classified by immunohistochemistry in ovarian carcinomas of 7,436 patients from the Ovarian Tumor Tissue Analysis consortium. We examined RB1 expression and germline BRCA status in a subset of 1,134 HGSC, and related genotype to overall survival (OS), tumor-infiltrating CD8+ lymphocytes, and transcriptomic subtypes. Using CRISPR-Cas9, we deleted RB1 in HGSC cells with and without BRCA1 alterations to model co-loss with treatment response. We performed whole-genome and transcriptome data analyses on 126 patients with primary HGSC to characterize tumors with concurrent BRCA deficiency and RB1 loss. RESULTS: RB1 loss was associated with longer OS in HGSC but with poorer prognosis in endometrioid ovarian carcinoma. Patients with HGSC harboring both RB1 loss and pathogenic germline BRCA variants had superior OS compared with patients with either alteration alone, and their median OS was three times longer than those without pathogenic BRCA variants and retained RB1 expression (9.3 vs. 3.1 years). Enhanced sensitivity to cisplatin and paclitaxel was seen in BRCA1-altered cells with RB1 knockout. Combined RB1 loss and BRCA deficiency correlated with transcriptional markers of enhanced IFN response, cell-cycle deregulation, and reduced epithelial-mesenchymal transition. CD8+ lymphocytes were most prevalent in BRCA-deficient HGSC with co-loss of RB1. CONCLUSIONS: Co-occurrence of RB1 loss and BRCA deficiency was associated with exceptionally long survival in patients with HGSC, potentially due to better treatment response and immune stimulation.

The development of the retina is under tight temporal and spatial control. To gain insights into the molecular basis of this process, we generate a single-nuclei dual-omic atlas of the human developing retina with approximately 220,000 nuclei from 14 human embryos and fetuses aged between 8 and 23-weeks post-conception with matched macular and peripheral tissues. This atlas captures all major cell classes in the retina, along with a large proportion of progenitors and cell-type-specific precursors. Cell trajectory analysis reveals a transition from continuous progression in early progenitors to a hierarchical development during the later stages of cell type specification. Both known and unrecorded candidate transcription factors, along with gene regulatory networks that drive the transitions of various cell fates, are identified. Comparisons between the macular and peripheral retinae indicate a largely consistent yet distinct developmental pattern. This atlas offers unparalleled resolution into the transcriptional and chromatin accessibility landscapes during development, providing an invaluable resource for deeper insights into retinal development and associated diseases.

HER2, encoded by the ERBB2 gene, is an important druggable driver of human cancer gaining increasing importance as a therapeutic target in urothelial carcinoma (UC). The genomic underpinnings of HER2 overexpression in ERBB2 nonamplified UC are poorly defined. To address this knowledge gap, we investigated 172 UC tumors from patients treated at the University of California San Francisco, using immunohistochemistry and next-generation sequencing. We found that GATA3 and PPARG copy number gains individually predicted HER2 protein expression independently of ERBB2 amplification. To validate these findings, we interrogated the Memorial Sloan Kettering/The Cancer Genome Atlas (MSK/TCGA) dataset and found that GATA3 and PPARG copy number gains individually predicted ERBB2 mRNA expression independently of ERBB2 amplification. Our findings reveal a potential link between the luminal marker HER2 and the key transcription factors GATA3 and PPARG in UC and highlight the utility of examining GATA3 and PPARG copy number states to identify UC tumors that overexpress HER2 in the absence of ERBB2 amplification. In summary, we found that an increase in copy number of GATA3 and PPARG was independently associated with higher ERBB2 expression in patient samples of UC. This finding provides a potential explanation for HER2 overexpression in UC tumors without ERBB2 amplification and a way to identify these tumors for HER2-targeted therapies.

Phosphorus is essential in all cells structural, metabolic and regulatory functions. For fungal cells that import inorganic phosphate (Pi) up a steep concentration gradient, surface Pi transporters are critical capacitators of growth. Fungi must deploy Pi transporters that enable optimal Pi uptake in pH and Pi concentration ranges prevalent in their environments. Single, triple and quadruple mutants were used to characterize the four Pi transporters we identified for the human fungal pathogen Candida albicans, which must adapt to alkaline conditions during invasion of the host bloodstream and deep organs. A high-affinity Pi transporter, Pho84, was most efficient across the widest pH range while another, Pho89, showed high-affinity characteristics only within one pH unit of neutral. Two low-affinity Pi transporters, Pho87 and Fgr2, were active only in acidic conditions. Only Pho84 among the Pi transporters was clearly required in previously identified Pi-related functions including Target of Rapamycin Complex 1 signaling, oxidative stress resistance and hyphal growth. We used in vitro evolution and whole genome sequencing as an unbiased forward genetic approach to probe adaptation to prolonged Pi scarcity of two quadruple mutant lineages lacking all 4 Pi transporters. Lineage-specific genomic changes corresponded to divergent success of the two lineages in fitness recovery during Pi limitation. Initial, large-scale genomic alterations like aneuploidies and loss of heterozygosity eventually resolved, as populations gained small-scale mutations. Severity of some phenotypes linked to Pi starvation, like cell wall stress hypersensitivity, decreased in parallel to evolving populations fitness recovery in Pi scarcity, while severity of others like membrane stress responses diverged from Pi scarcity fitness. Among preliminary candidate genes for contributors to fitness recovery, those with links to TORC1 were overrepresented. Since Pi homeostasis differs substantially between fungi and humans, adaptive processes to Pi deprivation may harbor small-molecule targets that impact fungal growth, stress resistance and virulence.

UNLABELLED: Antagonism of the PD-1/PD-L1 axis is a critical therapeutic strategy for patients with advanced bladder cancer. IFNγ functions as a key regulator of PD-L1 in both immune as well as cancer cells. Forkhead box P3 (FOXP3) is a transcription factor synonymous in T regulatory cell function but with increasingly described functions in cancer cells. Here, we investigated the relationship between FOXP3 and PD-L1 in bladder cancer. We showed that FOXP3 is critical in the ability for IFNγ to activate PD-L1 in bladder cancer cells. FOXP3 can bind to the PD-L1 promoter and induces a gene program that leads to regulation of multiple immune-related genes and genes involved in epithelial-to-mesenchymal transition (EMT). Using in vitro and in vivo human and murine models, we showed that FOXP3 can influence bladder cancer EMT as well as promote cancer metastases. Furthermore, FOXP3 may be a convergent factor for multiple activators of PD-L1, including the chemotherapeutic drug cisplatin. SIGNIFICANCE: Historically a key transcription factor driving T regulatory cell function, FOXP3 has an increasingly recognized role in cancer cells. In bladder cancer, we defined a novel mechanism whereby FOXP3 mediates the activation of the immune checkpoint PD-L1 by the cytokine IFNγ. We also showed that FOXP3 induces other immune checkpoints as well as genes involved in EMT, promoting immune resistance and cancer metastases.

Macrophages sense pathogens and orchestrate specific immune responses. Stimulus specificity is thought to be achieved through combinatorial and dynamical coding by signaling pathways. While NFκB dynamics are known to encode stimulus information, dynamical coding in other signaling pathways and their combinatorial coordination remain unclear. Here, we established live-cell microscopy to investigate how NFκB and p38 dynamics interface in stimulated macrophages. Information theory and machine learning revealed that p38 dynamics distinguish cytokine TNF from pathogen-associated molecular patterns and high doses from low, but contributed little to information-rich NFκB dynamics when both pathways are considered. This suggests that immune response genes benefit from decoding immune signaling dynamics or combinatorics, but not both. We found that the heterogeneity of the two pathways is surprisingly uncorrelated. Mathematical modeling revealed potential sources of uncorrelated heterogeneity in the branched pathway network topology and predicted it to drive gene expression variability. Indeed, genes dependent on both p38 and NFκB showed high scRNAseq variability and bimodality. These results identify combinatorial signaling as a mechanism to restrict NFκB-AND-p38-responsive inflammatory cytokine expression to few cells.

The impact of variations in the three-dimensional structure of the genome has been recognized, but solid cancer tissue studies are limited. Here, we performed integrated deep Hi-C sequencing with matched whole-genome sequencing, whole-genome bisulfite sequencing, 5-hydroxymethylcytosine (5hmC) sequencing and RNA sequencing across a cohort of 80 biopsy samples from patients with metastatic castration-resistant prostate cancer. Dramatic differences were present in gene expression, 5-methylcytosine/5hmC methylation and in structural variation versus mutation rate between A and B (open and closed) chromatin compartments. A subset of tumors exhibited depleted regional chromatin contacts at the AR locus, linked to extrachromosomal circular DNA (ecDNA) and worse response to AR signaling inhibitors. We also identified topological subtypes associated with stark differences in methylation structure, gene expression and prognosis. Our data suggested that DNA interactions may predispose to structural variant formation, exemplified by the recurrent TMPRSS2-ERG fusion. This comprehensive integrated sequencing effort represents a unique clinical tumor resource.

The androgen receptor (AR) is the main driver in the development of castration-resistant prostate cancer, where the emergence of AR splice variants leads to treatment-resistant disease. Through detailed molecular studies of the marine alkaloid manzamine A (MA), we identified transcription factor E2F8 as a previously unknown regulator of AR transcription that prevents AR synthesis in prostate cancer cells. MA significantly inhibited the growth of various prostate cancer cell lines and was highly effective in inhibiting xenograft tumor growth in mice without any pathophysiological perturbations in major organs. MA suppressed the full-length AR (AR-FL), its spliced variant AR-V7, and the AR-regulated prostate-specific antigen (PSA; also known as KLK3) and human kallikrein 2 (hK2; also known as KLK2) genes. RNA sequencing (RNA-seq) analysis and protein modeling studies revealed E2F8 interactions with DNA as a potential novel target of MA, suppressing AR transcription and its synthesis. This novel mechanism of blocking AR biogenesis via E2F8 may provide an opportunity to control therapy-resistant prostate cancer over the currently used AR antagonists designed to target different parts of the AR gene.

BACKGROUND: Deeper insights into ERBB2-driven cancers are essential to develop new treatment approaches for ERBB2+ breast cancers (BCs). We employed the Collaborative Cross (CC) mouse model to unearth genetic factors underpinning Erbb2-driven mammary tumour development and metastasis. METHODS: 732 F1 hybrid female mice between FVB/N MMTV-Erbb2 and 30 CC strains were monitored for mammary tumour phenotypes. GWAS pinpointed SNPs that influence various tumour phenotypes. Multivariate analyses and models were used to construct the polygenic score and to develop a mouse tumour susceptibility gene signature (mTSGS), where the corresponding human ortholog was identified and designated as hTSGS. The importance and clinical value of hTSGS in human BC was evaluated using public datasets, encompassing TCGA, METABRIC, GSE96058, and I-SPY2 cohorts. The predictive power of mTSGS for response to chemotherapy was validated in vivo using genetically diverse MMTV-Erbb2 mice. FINDINGS: Distinct variances in tumour onset, multiplicity, and metastatic patterns were observed in F1-hybrid female mice between FVB/N MMTV-Erbb2 and 30 CC strains. Besides lung metastasis, liver and kidney metastases emerged in specific CC strains. GWAS identified specific SNPs significantly associated with tumour onset, multiplicity, lung metastasis, and liver metastasis. Multivariate analyses flagged SNPs in 20 genes (Stx6, Ramp1, Traf3ip1, Nckap5, Pfkfb2, Trmt1l, Rprd1b, Rer1, Sepsecs, Rhobtb1, Tsen15, Abcc3, Arid5b, Tnr, Dock2, Tti1, Fam81a, Oxr1, Plxna2, and Tbc1d31) independently tied to various tumour characteristics, designated as a mTSGS. hTSGS scores (hTSGSS) based on their transcriptional level showed prognostic values, superseding clinical factors and PAM50 subtype across multiple human BC cohorts, and predicted pathological complete response independent of and superior to MammaPrint score in I-SPY2 study. The power of mTSGS score for predicting chemotherapy response was further validated in an in vivo mouse MMTV-Erbb2 model, showing that, like findings in human patients, mouse tumours with low mTSGS scores were most likely to respond to treatment. INTERPRETATION: Our investigation has unveiled many new genes predisposing individuals to ERBB2-driven cancer. Translational findings indicate that hTSGS holds promise as a biomarker for refining treatment strategies for patients with BC. FUNDING: The U.S. Department of Defense (DoD) Breast Cancer Research Program (BCRP) (BC190820), United States; MCIN/AEI/10.13039/501100011039 (PID2020-118527RB-I00, PDC2021-121735-I00), the European Union Next Generation EU/PRTR, the Regional Government of Castile and León (CSI144P20), European Union.

IGHMBP2 is a nonessential, superfamily 1 DNA/RNA helicase that is mutated in patients with rare neuromuscular diseases SMARD1 and CMT2S. IGHMBP2 is implicated in translational and transcriptional regulation via biochemical association with ribosomal proteins, pre-rRNA processing factors, and tRNA-related species. To uncover the cellular consequences of perturbing IGHMBP2, we generated full and partial IGHMBP2 deletion K562 cell lines. Using polysome profiling and a nascent protein synthesis assay, we found that IGHMBP2 deletion modestly reduces global translation. We performed Ribo-seq and RNA-seq and identified diverse gene expression changes due to IGHMBP2 deletion, including ATF4 up-regulation. With recent studies showing the integrated stress response (ISR) can contribute to tRNA metabolism-linked neuropathies, we asked whether perturbing IGHMBP2 promotes ISR activation. We generated ATF4 reporter cell lines and found IGHMBP2 knockout cells demonstrate basal, chronic ISR activation. Our work expands upon the impact of IGHMBP2 in translation and elucidates molecular mechanisms that may link mutant IGHMBP2 to severe clinical phenotypes.

UNLABELLED: Vibrio cholerae, the causative agent of the diarrheal disease cholera, poses an ongoing health threat due to its wide repertoire of horizontally acquired elements (HAEs) and virulence factors. New clinical isolates of the bacterium with improved fitness abilities, often associated with HAEs, frequently emerge. The appropriate control and expression of such genetic elements is critical for the bacteria to thrive in the different environmental niches they occupy. H-NS, the histone-like nucleoid structuring protein, is the best-studied xenogeneic silencer of HAEs in gamma-proteobacteria. Although H-NS and other highly abundant nucleoid-associated proteins (NAPs) have been shown to play important roles in regulating HAEs and virulence in model bacteria, we still lack a comprehensive understanding of how different NAPs modulate transcription in V. cholerae. By obtaining genome-wide measurements of protein occupancy and active transcription in a clinical isolate of V. cholerae, harboring recently discovered HAEs encoding for phage defense systems, we show that a lack of H-NS causes a robust increase in the expression of genes found in many HAEs. We further found that TsrA, a protein with partial homology to H-NS, regulates virulence genes primarily through modulation of H-NS activity. We also identified few sites that are affected by TsrA independently of H-NS, suggesting TsrA may act with diverse regulatory mechanisms. Our results demonstrate how the combinatorial activity of NAPs is employed by a clinical isolate of an important pathogen to regulate recently discovered HAEs. IMPORTANCE: New strains of the bacterial pathogen Vibrio cholerae, bearing novel horizontally acquired elements (HAEs), frequently emerge. HAEs provide beneficial traits to the bacterium, such as antibiotic resistance and defense against invading bacteriophages. Xenogeneic silencers are proteins that help bacteria harness new HAEs and silence those HAEs until they are needed. H-NS is the best-studied xenogeneic silencer; it is one of the nucleoid-associated proteins (NAPs) in gamma-proteobacteria and is responsible for the proper regulation of HAEs within the bacterial transcriptional network. We studied the effects of H-NS and other NAPs on the HAEs of a clinical isolate of V. cholerae. Importantly, we found that H-NS partners with a small and poorly characterized protein, TsrA, to help domesticate new HAEs involved in bacterial survival and in causing disease. A proper understanding of the regulatory state in emerging isolates of V. cholerae will provide improved therapies against new isolates of the pathogen.

Patterns of transcriptional activity are encoded in our genome through regulatory elements such as promoters or enhancers that, paradoxically, contain similar assortments of sequence-specific transcription factor (TF) binding sites1-3. Knowledge of how these sequence motifs encode multiple, often overlapping, gene expression programs is central to understanding gene regulation and how mutations in non-coding DNA manifest in disease4,5. Here, by studying gene regulation from the perspective of individual transcription start sites (TSSs), using natural genetic variation, perturbation of endogenous TF protein levels and massively parallel analysis of natural and synthetic regulatory elements, we show that the effect of TF binding on transcription initiation is position dependent. Analysing TF-binding-site occurrences relative to the TSS, we identified several motifs with highly preferential positioning. We show that these patterns are a combination of a TF's distinct functional profiles-many TFs, including canonical activators such as NRF1, NFY and Sp1, activate or repress transcription initiation depending on their precise position relative to the TSS. As such, TFs and their spacing collectively guide the site and frequency of transcription initiation. More broadly, these findings reveal how similar assortments of TF binding sites can generate distinct gene regulatory outcomes depending on their spatial configuration and how DNA sequence polymorphisms may contribute to transcription variation and disease and underscore a critical role for TSS data in decoding the regulatory information of our genome.

Treatment-induced neuroendocrine prostate cancer (t-NEPC) often arises from adenocarcinoma via lineage plasticity in response to androgen receptor signaling inhibitors, such as enzalutamide. However, the specific regulators and targets involved in the transition to NEPC are not well understood. Plexin D1 (PLXND1) is a cellular receptor of the semaphorin (SEMA) family that plays important roles in modulating the cytoskeleton and cell adhesion. Here, we found that PLXND1 was highly expressed and positively correlated with neuroendocrine markers in patients with NEPC. High PLXND1 expression was associated with poorer prognosis in prostate cancer patients. Additionally, PLXND1 was upregulated and negatively regulated by androgen receptor signaling in enzalutamide-resistant cells. Knockdown or knockout of PLXND1 inhibited neural lineage pathways, thereby suppressing NEPC cell proliferation, patient derived xenograft (PDX) tumor organoid viability, and xenograft tumor growth. Mechanistically, the heat shock protein 70 (HSP70) regulated PLXND1 protein stability through degradation, and inhibition of HSP70 decreased PLXND1 expression and NEPC organoid growth. In summary, our findings indicate that PLXND1 could serve as a promising therapeutic target and molecular marker for NEPC.

Melanoma clinical outcomes emerge from incompletely understood genetic mechanisms operating within the tumor and its microenvironment. Here, we used single-cell RNA-based spatial molecular imaging (RNA-SMI) in patient-derived archival tumors to reveal clinically relevant markers of malignancy progression and prognosis. We examined spatial gene expression of 203,472 cells inside benign and malignant melanocytic neoplasms, including melanocytic nevi and primary invasive and metastatic melanomas. Algorithmic cell clustering paired with intratumoral comparative two-dimensional analyses visualized synergistic, spatial gene signatures linking cellular proliferation, metabolism, and malignancy, validated by protein expression. Metastatic niches included up-regulation of CDK2 and FABP5, which independently predicted poor clinical outcome in 473 patients with melanoma via Cox regression analysis. More generally, our work demonstrates a framework for applying single-cell RNA-SMI technology toward identifying gene regulatory landscapes pertinent to cancer progression and patient survival.

The prevalence and progression of cancer differ in males and females, and thus, sexual dimorphism in tumor development directly impacts clinical research and medicine. Long non-coding RNAs (lncRNAs) are increasingly recognized as important players in gene expression and various cellular processes, including cancer development and progression. In recent years, lncRNAs have been implicated in the differences observed in cancer incidence, progression, and treatment responses between men and women. Here, we present a brief overview of the current knowledge regarding the role of lncRNAs in cancer sex dimorphism, focusing on how they affect epigenetic processes in male and female mammalian cells. We discuss the potential mechanisms by which lncRNAs may contribute to sex differences in cancer, including transcriptional control of sex chromosomes, hormonal signaling pathways, and immune responses. We also propose strategies for studying lncRNA functions in cancer sex dimorphism. Furthermore, we emphasize the importance of considering sex as a biological variable in cancer research and the need to investigate the role lncRNAs play in mediating these sex differences. In summary, we highlight the emerging link between lncRNAs and cancer sex dimorphism and their potential as therapeutic targets.

Glioblastoma (GBM) is the most common malignant primary brain tumor and remains incurable. Previous work has shown that systemic administration of Decitabine (DAC) induces sufficient expression of cancer-testis antigens (CTA) in GBM for targeting by adoptive T-cell therapy in vivo. However, the mechanisms by which DAC enhances immunogenicity in GBM remain to be elucidated. Using New York esophageal squamous cell carcinoma 1 (NY-ESO-1) as a representative inducible CTA, we demonstrate in patient tissue, immortalized glioma cells, and primary patient-derived gliomaspheres that basal CTA expression is restricted by promoter hypermethylation in gliomas. DAC treatment of glioma cells specifically inhibits DNA methylation silencing to render NY-ESO-1 and other CTA into inducible tumor antigens at single-cell resolution. Functionally, NY-ESO-1 T-cell receptor-engineered effector cell targeting of DAC-induced antigen in primary glioma cells promotes specific and polyfunctional T-cell cytokine profiles. In addition to induction of CTA, DAC concomitantly reactivates tumor-intrinsic human endogenous retroviruses, interferon response signatures, and MHC-I. Overall, we demonstrate that DAC induces targetable tumor antigen and enhances T-cell functionality against GBM, ultimately contributing to the improvement of targeted immune therapies in glioma.SignificanceThis study dissects the tumor-intrinsic epigenetic and transcriptional mechanisms underlying enhanced T-cell functionality targeting decitabine-induced cancer-testis antigens in glioma. Our findings demonstrate concomitant induction of tumor antigens, reactivation of human endogenous retroviruses, and stimulation of interferon signaling as a mechanistic rationale to epigenetically prime human gliomas to immunotherapeutic targeting.

The high heritability of amyotrophic lateral sclerosis (ALS) contrasts with its low molecular diagnosis rate post-genetic testing, pointing to potential undiscovered genetic factors. To aid the exploration of these factors, we introduced EpiOut, an algorithm to identify chromatin accessibility outliers that are regions exhibiting divergent accessibility from the population baseline in a single or few samples. Annotation of accessible regions with histone chromatin immunoprecipitation sequencing and Hi-C indicates that outliers are concentrated in functional loci, especially among promoters interacting with active enhancers. Across different omics levels, outliers are robustly replicated, and chromatin accessibility outliers are reliable predictors of gene expression outliers and aberrant protein levels. When promoter accessibility does not align with gene expression, our results indicate that molecular aberrations are more likely to be linked to post-transcriptional regulation rather than transcriptional regulation. Our findings demonstrate that the outlier detection paradigm can uncover dysregulated regions in rare diseases. EpiOut is available at github.com/uci-cbcl/EpiOut.

Long noncoding RNAs (LncRNAs) are key regulators of gene expression and can mediate their effects in both the nucleus and cytoplasm. Some of the best-characterized lncRNAs are localized within the nucleus, where they modulate the nuclear architecture and influence gene expression. In this review, we discuss the role of lncRNAs in nuclear architecture in the context of their gene regulatory functions in innate immunity. Here, we discuss various approaches to functionally characterize nuclear-localized lncRNAs and the challenges faced in the field.

Extrachromosomal DNA is a common cause of oncogene amplification in cancer. The non-chromosomal inheritance of ecDNA enables tumors to rapidly evolve, contributing to treatment resistance and poor outcome for patients. The transcriptional context in which ecDNAs arise and progress, including chromosomally-driven transcription, is incompletely understood. We examined gene expression patterns of 870 tumors of varied histological types, to identify transcriptional correlates of ecDNA. Here, we show that ecDNA-containing tumors impact four major biological processes. Specifically, ecDNA-containing tumors up-regulate DNA damage and repair, cell cycle control, and mitotic processes, but down-regulate global immune regulation pathways. Taken together, these results suggest profound alterations in gene regulation in ecDNA-containing tumors, shedding light on molecular processes that give rise to their development and progression.

The cellular components of tumors and their microenvironment play pivotal roles in tumor progression, patient survival, and the response to cancer treatments. Unveiling a comprehensive cellular profile within bulk tumors via single-cell RNA sequencing (scRNA-seq) data is crucial, as it unveils intrinsic tumor cellular traits that elude identification through conventional cancer subtyping methods. Our contribution, scBeacon, is a tool that derives cell-type signatures by integrating and clustering multiple scRNA-seq datasets to extract signatures for deconvolving unrelated tumor datasets on bulk samples. Through the employment of scBeacon on the The Cancer Genome Atlas (TCGA) cohort, we find cellular and molecular attributes within specific tumor categories, many with patient outcome relevance. We developed a tumor cell-type map to visually depict the relationships among TCGA samples based on the cell-type inferences.

Meningiomas, although mostly benign, can be recurrent and fatal. World Health Organization (WHO) grading of the tumor does not always identify high-risk meningioma, and better characterizations of their aggressive biology are needed. To approach this problem, we combined 13 bulk RNA sequencing (RNA-seq) datasets to create a dimension-reduced reference landscape of 1,298 meningiomas. The clinical and genomic metadata effectively correlated with landscape regions, which led to the identification of meningioma subtypes with specific biological signatures. The time to recurrence also correlated with the map location. Further, we developed an algorithm that maps new patients onto this landscape, where the nearest neighbors predict outcome. This study highlights the utility of combining bulk transcriptomic datasets to visualize the complexity of tumor populations. Further, we provide an interactive tool for understanding the disease and predicting patient outcomes. This resource is accessible via the online tool Oncoscape, where the scientific community can explore the meningioma landscape.

Many autism spectrum disorder (ASD)-associated genes act as transcriptional regulators (TRs). Chromatin immunoprecipitation sequencing (ChIP-seq) was used to identify the regulatory targets of ARID1B, BCL11A, FOXP1, TBR1, and TCF7L2, ASD-associated TRs in the developing human and mouse cortex. These TRs shared substantial overlap in the binding sites, especially within open chromatin. The overlap within a promoter region, 1-2,000 bp upstream of the transcription start site, was highly predictive of brain-expressed genes. This signature was observed in 96 out of 102 ASD-associated genes. In vitro CRISPRi against ARID1B and TBR1 delineated downstream convergent biology in mouse cortical cultures. After 8 days, NeuN+ and CALB+ cells were decreased, GFAP+ cells were increased, and transcriptomic signatures correlated with the postmortem brain samples from individuals with ASD. We suggest that functional convergence across five ASD-associated TRs leads to shared neurodevelopmental outcomes of haploinsufficient disruption.

Intratumor heterogeneity underlies cancer evolution and treatment resistance, but targetable mechanisms driving intratumor heterogeneity are poorly understood. Meningiomas are the most common primary intracranial tumors and are resistant to all medical therapies, and high-grade meningiomas have significant intratumor heterogeneity. Here we use spatial approaches to identify genomic, biochemical and cellular mechanisms linking intratumor heterogeneity to the molecular, temporal and spatial evolution of high-grade meningiomas. We show that divergent intratumor gene and protein expression programs distinguish high-grade meningiomas that are otherwise grouped together by current classification systems. Analyses of matched pairs of primary and recurrent meningiomas reveal spatial expansion of subclonal copy number variants associated with treatment resistance. Multiplexed sequential immunofluorescence and deconvolution of meningioma spatial transcriptomes using cell types from single-cell RNA sequencing show decreased immune infiltration, decreased MAPK signaling, increased PI3K-AKT signaling and increased cell proliferation, which are associated with meningioma recurrence. To translate these findings to preclinical models, we use CRISPR interference and lineage tracing approaches to identify combination therapies that target intratumor heterogeneity in meningioma cell co-cultures.

A quarter of humanity is estimated to have been exposed to Mycobacterium tuberculosis (Mtb) with a 5-10% risk of developing tuberculosis (TB) disease. Variability in responses to Mtb infection could be due to host or pathogen heterogeneity. Here, we focused on host genetic variation in a Peruvian population and its associations with gene regulation in monocyte-derived macrophages and dendritic cells (DCs). We recruited former household contacts of TB patients who previously progressed to TB (cases, n = 63) or did not progress to TB (controls, n = 63). Transcriptomic profiling of monocyte-derived DCs and macrophages measured the impact of genetic variants on gene expression by identifying expression quantitative trait loci (eQTL). We identified 330 and 257 eQTL genes in DCs and macrophages (False Discovery Rate (FDR) < 0.05), respectively. Four genes in DCs showed interaction between eQTL variants and TB progression status. The top eQTL interaction for a protein-coding gene was with FAH, the gene encoding fumarylacetoacetate hydrolase, which mediates the last step in mammalian tyrosine catabolism. FAH expression was associated with genetic regulatory variation in cases but not controls. Using public transcriptomic and epigenomic data of Mtb-infected monocyte-derived dendritic cells, we found that Mtb infection results in FAH downregulation and DNA methylation changes in the locus. Overall, this study demonstrates effects of genetic variation on gene expression levels that are dependent on history of infectious disease and highlights a candidate pathogenic mechanism through pathogen-response genes. Furthermore, our results point to tyrosine metabolism and related candidate TB progression pathways for further investigation.

Chondrocyte differentiation controls skeleton development and stature. Here we provide a comprehensive map of chondrocyte-specific enhancers and show that they provide a mechanistic framework through which non-coding genetic variants can influence skeletal development and human stature. Working with fetal chondrocytes isolated from mice bearing a Col2a1 fluorescent regulatory sensor, we identify 780 genes and 2'704 putative enhancers specifically active in chondrocytes using a combination of RNA-seq, ATAC-seq and H3K27ac ChIP-seq. Most of these enhancers (74%) show pan-chondrogenic activity, with smaller populations being restricted to limb (18%) or trunk (8%) chondrocytes only. Notably, genetic variations overlapping these enhancers better explain height differences than those overlapping non-chondrogenic enhancers. Finally, targeted deletions of identified enhancers at the Fgfr3, Col2a1, Hhip and, Nkx3-2 loci confirm their role in regulating cognate genes. This enhancer map provides a framework for understanding how genes and non-coding variations influence bone development and diseases.

Cis-regulatory elements (CREs) interact with trans regulators to orchestrate gene expression, but how transcriptional regulation is coordinated in multi-gene loci has not been experimentally defined. We sought to characterize the CREs controlling dynamic expression of the adjacent costimulatory genes CD28, CTLA4 and ICOS, encoding regulators of T cell-mediated immunity. Tiling CRISPR interference (CRISPRi) screens in primary human T cells, both conventional and regulatory subsets, uncovered gene-, cell subset- and stimulation-specific CREs. Integration with CRISPR knockout screens and assay for transposase-accessible chromatin with sequencing (ATAC-seq) profiling identified trans regulators influencing chromatin states at specific CRISPRi-responsive elements to control costimulatory gene expression. We then discovered a critical CCCTC-binding factor (CTCF) boundary that reinforces CRE interaction with CTLA4 while also preventing promiscuous activation of CD28. By systematically mapping CREs and associated trans regulators directly in primary human T cell subsets, this work overcomes longstanding experimental limitations to decode context-dependent gene regulatory programs in a complex, multi-gene locus critical to immune homeostasis.

Long noncoding RNAs (lncRNAs) account for the largest portion of RNA from the transcriptome, yet most of their functions remain unknown. Here, we performed two independent high-throughput CRISPRi screens to understand the role of lncRNAs in monocyte function and differentiation. The first was a reporter-based screen to identify lncRNAs that regulate TLR4-NFkB signaling in human monocytes and the second screen identified lncRNAs involved in monocyte to macrophage differentiation. We successfully identified numerous noncoding and protein-coding genes that can positively or negatively regulate inflammation and differentiation. To understand the functional roles of lncRNAs in both processes, we chose to further study the lncRNA LOUP [lncRNA originating from upstream regulatory element of SPI1 (also known as PU.1)], as it emerged as a top hit in both screens. Not only does LOUP regulate its neighboring gene, the myeloid fate-determining factor SPI1, thereby affecting monocyte to macrophage differentiation, but knockdown of LOUP leads to a broad upregulation of NFkB-targeted genes at baseline and upon TLR4-NFkB activation. LOUP also harbors three small open reading frames capable of being translated and are responsible for LOUPs ability to negatively regulate TLR4/NFkB signaling. This work emphasizes the value of high-throughput screening to rapidly identify functional lncRNAs in the innate immune system.

AbstractThe transcription factor RUNX1 is a master regulator of hematopoiesis and is frequently mutated in myeloid malignancies. Mutations in its runt homology domain (RHD) frequently disrupt DNA binding and result in loss of RUNX1 function. However, it is not clearly understood how other RUNX1 mutations contribute to disease development. Here, we characterized RUNX1 mutations outside of the RHD. Our analysis of the patient data sets revealed that mutations within the C-terminus frequently occur in hematopoietic disorders. Remarkably, most of these mutations were nonsense or frameshift mutations and were predicted to be exempt from nonsense-mediated messenger RNA decay. Therefore, this class of mutation is projected to produce DNA-binding proteins that contribute to the pathogenesis in a distinct manner. To model this, we introduced the RUNX1R320∗ mutation into the endogenous gene locus and demonstrated the production of RUNX1R320∗ protein. Expression of RUNX1R320∗ resulted in the disruption of RUNX1 regulated processes such as megakaryocytic differentiation, through a transcriptional signature different from RUNX1 depletion. To understand the underlying mechanisms, we used Global RNA Interactions with DNA by deep sequencing (GRID-seq) to examine enhancer-promoter connections. We identified widespread alterations in the enhancer-promoter networks within RUNX1 mutant cells. Additionally, we uncovered enrichment of RUNX1R320∗ and FOXK2 binding at the MYC super enhancer locus, significantly upregulating MYC transcription and signaling pathways. Together, our study demonstrated that most RUNX1 mutations outside the DNA-binding domain are not subject to nonsense-mediated decay, producing protein products that act in concert with additional cofactors to dysregulate hematopoiesis through mechanisms distinct from those induced by RUNX1 depletion.

Neuropsychiatric genome-wide association studies (GWASs), including those for autism spectrum disorder and schizophrenia, show strong enrichment for regulatory elements in the developing brain. However, prioritizing risk genes and mechanisms is challenging without a unified regulatory atlas. Across 672 diverse developing human brains, we identified 15,752 genes harboring gene, isoform, and/or splicing quantitative trait loci, mapping 3739 to cellular contexts. Gene expression heritability drops during development, likely reflecting both increasing cellular heterogeneity and the intrinsic properties of neuronal maturation. Isoform-level regulation, particularly in the second trimester, mediated the largest proportion of GWAS heritability. Through colocalization, we prioritized mechanisms for about 60% of GWAS loci across five disorders, exceeding adult brain findings. Finally, we contextualized results within gene and isoform coexpression networks, revealing the comprehensive landscape of transcriptome regulation in development and disease.

Nucleotide changes in gene regulatory elements are important determinants of neuronal development and diseases. Using massively parallel reporter assays in primary human cells from mid-gestation cortex and cerebral organoids, we interrogated the cis-regulatory activity of 102,767 open chromatin regions, including thousands of sequences with cell type-specific accessibility and variants associated with brain gene regulation. In primary cells, we identified 46,802 active enhancer sequences and 164 variants that alter enhancer activity. Activity was comparable in organoids and primary cells, suggesting that organoids provide an adequate model for the developing cortex. Using deep learning we decoded the sequence basis and upstream regulators of enhancer activity. This work establishes a comprehensive catalog of functional gene regulatory elements and variants in human neuronal development.

Osteosarcoma (OS) is the most common malignant bone tumor with a poor prognosis. Here, we show that the nuclear receptor RORγ may serve as a potential therapeutic target in OS. OS exhibits a hyperactivated oxidative phosphorylation (OXPHOS) program, which fuels the carbon source to promote tumor progression. We found that RORγ is overexpressed in OS tumors and is linked to hyperactivated OXPHOS. RORγ induces the expression of PGC-1β and physically interacts with it to activate the OXPHOS program by upregulating the expression of respiratory chain component genes. Inhibition of RORγ strongly inhibits OXPHOS activation, downregulates mitochondrial functions, and increases ROS production, which results in OS cell apoptosis and ferroptosis. RORγ inverse agonists strongly suppressed OS tumor growth and progression and sensitized OS tumors to chemotherapy. Taken together, our results indicate that RORγ is a critical regulator of the OXPHOS program in OS and provides an effective therapeutic strategy for this deadly disease.

Mediator kinases CDK19 and CDK8, pleiotropic regulators of transcriptional reprogramming, are differentially regulated by androgen signaling, but both kinases are upregulated in castration-resistant prostate cancer (CRPC). Genetic or pharmacological inhibition of CDK8 and CDK19 reverses the castration-resistant phenotype and restores the sensitivity of CRPC xenografts to androgen deprivation in vivo. Prolonged CDK8/19 inhibitor treatment combined with castration not only suppressed the growth of CRPC xenografts but also induced tumor regression and cures. Transcriptomic analysis revealed that Mediator kinase inhibition amplified and modulated the effects of castration on gene expression, disrupting CRPC adaptation to androgen deprivation. Mediator kinase inactivation in tumor cells also affected stromal gene expression, indicating that Mediator kinase activity in CRPC molded the tumor microenvironment. The combination of castration and Mediator kinase inhibition downregulated the MYC pathway, and Mediator kinase inhibition suppressed a MYC-driven CRPC tumor model even without castration. CDK8/19 inhibitors showed efficacy in patient-derived xenograft models of CRPC, and a gene signature of Mediator kinase activity correlated with tumor progression and overall survival in clinical samples of metastatic CRPC. These results indicate that Mediator kinases mediated androgen-independent in vivo growth of CRPC, supporting the development of CDK8/19 inhibitors for the treatment of this presently incurable disease.

The SUVmax is a measure of FDG uptake and is related with tumor aggressiveness in thyroid cancer, however, its association with molecular pathways is unclear. Here, we investigated the relationship between SUVmax and gene expression profiles in 80 papillary thyroid cancer (PTC) patients. We conducted an analysis of DEGs and enriched pathways in relation to SUVmax and tumor size. SUVmax showed a positive correlation with tumor size and correlated with glucose metabolic process. The genes that indicate thyroid differentiation, such as SLC5A5 and TPO, were negatively correlated with SUVmax. Unsupervised analysis revealed that SUVmax positively correlated with DNA replication(r = 0.29, p = 0.009), pyrimidine metabolism(r = 0.50, p

During development, neural stem cells in the cerebral cortex, also known as radial glial cells (RGCs), generate excitatory neurons, followed by production of cortical macroglia and inhibitory neurons that migrate to the olfactory bulb (OB). Understanding the mechanisms for this lineage switch is fundamental for unraveling how proper numbers of diverse neuronal and glial cell types are controlled. We and others recently showed that Sonic Hedgehog (Shh) signaling promotes the cortical RGC lineage switch to generate cortical oligodendrocytes and OB interneurons. During this process, cortical RGCs generate intermediate progenitor cells that express critical gliogenesis genes Ascl1, Egfr, and Olig2. The increased Ascl1 expression and appearance of Egfr+ and Olig2+ cortical progenitors are concurrent with the switch from excitatory neurogenesis to gliogenesis and OB interneuron neurogenesis in the cortex. While Shh signaling promotes Olig2 expression in the developing spinal cord, the exact mechanism for this transcriptional regulation is not known. Furthermore, the transcriptional regulation of Olig2 and Egfr has not been explored. Here, we show that in cortical progenitor cells, multiple regulatory programs, including Pax6 and Gli3, prevent precocious expression of Olig2, a gene essential for production of cortical oligodendrocytes and astrocytes. We identify multiple enhancers that control Olig2 expression in cortical progenitors and show that the mechanisms for regulating Olig2 expression are conserved between the mouse and human. Our study reveals evolutionarily conserved regulatory logic controlling the lineage switch of cortical neural stem cells.

We explored the dysregulation of G-protein-coupled receptor (GPCR) ligand systems in cancer transcriptomics datasets to uncover new therapeutics opportunities in oncology. We derived an interaction network of receptors with ligands and their biosynthetic enzymes. Multiple GPCRs are differentially regulated together with their upstream partners across cancer subtypes and are associated to specific transcriptional programs and to patient survival patterns. The expression of both receptor-ligand (or enzymes) partners improved patient stratification, suggesting a synergistic role for the activation of GPCR networks in modulating cancer phenotypes. Remarkably, we identified many such axes across several cancer molecular subtypes, including many involving receptor-biosynthetic enzymes for neurotransmitters. We found that GPCRs from these actionable axes, including, e.g., muscarinic, adenosine, 5-hydroxytryptamine, and chemokine receptors, are the targets of multiple drugs displaying anti-growth effects in large-scale, cancer cell drug screens, which we further validated. We have made the results generated in this study freely available through a webapp (gpcrcanceraxes.bioinfolab.sns.it).

Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation.

UNLABELLED: Patients with estrogen receptor-positive breast cancer receive adjuvant endocrine therapies (ET) that delay relapse by targeting clinically undetectable micrometastatic deposits. Yet, up to 50% of patients relapse even decades after surgery through unknown mechanisms likely involving dormancy. To investigate genetic and transcriptional changes underlying tumor awakening, we analyzed late relapse patients and longitudinally profiled a rare cohort treated with long-term neoadjuvant ETs until progression. Next, we developed an in vitro evolutionary study to record the adaptive strategies of individual lineages in unperturbed parallel experiments. Our data demonstrate that ETs induce nongenetic cell state transitions into dormancy in a stochastic subset of cells via epigenetic reprogramming. Single lineages with divergent phenotypes awaken unpredictably in the absence of recurrent genetic alterations. Targeting the dormant epigenome shows promising activity against adapting cancer cells. Overall, this study uncovers the contribution of epigenetic adaptation to the evolution of resistance to ETs. SIGNIFICANCE: This study advances the understanding of therapy-induced dormancy with potential clinical implications for breast cancer. Estrogen receptor-positive breast cancer cells adapt to endocrine treatment by entering a dormant state characterized by strong heterochromatinization with no recurrent genetic changes. Targeting the epigenetic rewiring impairs the adaptation of cancer cells to ETs. See related commentary by Llinas-Bertran et al., p. 704. This article is featured in Selected Articles from This Issue, p. 695.

Phenotypic variation among species is a product of evolutionary changes to developmental programs1,2. However, how these changes generate novel morphological traits remains largely unclear. Here we studied the genomic and developmental basis of the mammalian gliding membrane, or patagium-an adaptative trait that has repeatedly evolved in different lineages, including in closely related marsupial species. Through comparative genomic analysis of 15 marsupial genomes, both from gliding and non-gliding species, we find that the Emx2 locus experienced lineage-specific patterns of accelerated cis-regulatory evolution in gliding species. By combining epigenomics, transcriptomics and in-pouch marsupial transgenics, we show that Emx2 is a critical upstream regulator of patagium development. Moreover, we identify different cis-regulatory elements that may be responsible for driving increased Emx2 expression levels in gliding species. Lastly, using mouse functional experiments, we find evidence that Emx2 expression patterns in gliders may have been modified from a pre-existing program found in all mammals. Together, our results suggest that patagia repeatedly originated through a process of convergent genomic evolution, whereby regulation of Emx2 was altered by distinct cis-regulatory elements in independently evolved species. Thus, different regulatory elements targeting the same key developmental gene may constitute an effective strategy by which natural selection has harnessed regulatory evolution in marsupial genomes to generate phenotypic novelty.

UNLABELLED: Combinations of immune checkpoint inhibitors (ICI, including anti-PD-1/PD-L1) and chemotherapy have been FDA approved for metastatic and early-stage triple-negative breast cancer (TNBC), but most patients do not benefit. B7-H4 is a B7 family ligand with proposed immunosuppressive functions being explored as a cancer immunotherapy target and may be associated with anti-PD-L1 resistance. However, little is known about its regulation and effect on immune cell function in breast cancers. We assessed murine and human breast cancer cells to identify regulation mechanisms of B7-H4 in vitro. We used an immunocompetent anti-PD-L1-sensitive orthotopic mammary cancer model and induced ectopic expression of B7-H4. We assessed therapy response and transcriptional changes at baseline and under treatment with anti-PD-L1. We observed B7-H4 was highly associated with epithelial cell status and transcription factors and found to be regulated by PI3K activity. EMT6 tumors with cell-surface B7-H4 expression were more resistant to immunotherapy. In addition, tumor-infiltrating immune cells had reduced immune activation signaling based on transcriptomic analysis. Paradoxically, in human breast cancer, B7-H4 expression was associated with survival benefit for patients with metastatic TNBC treated with carboplatin plus anti-PD-L1 and was associated with no change in response or survival for patients with early breast cancer receiving chemotherapy plus anti-PD-1. While B7-H4 induces tumor resistance to anti-PD-L1 in murine models, there are alternative mechanisms of signaling and function in human cancers. In addition, the strong correlation of B7-H4 to epithelial cell markers suggests a potential regulatory mechanism of B7-H4 independent of PD-L1. SIGNIFICANCE: This translational study confirms the association of B7-H4 expression with a cold immune microenvironment in breast cancer and offers preclinical studies demonstrating a potential role for B7-H4 in suppressing response to checkpoint therapy. However, analysis of two clinical trials with checkpoint inhibitors in the early and metastatic settings argue against B7-H4 as being a mechanism of clinical resistance to checkpoints, with clear implications for its candidacy as a therapeutic target.

BACKGROUND: Loss of AZGP1 expression is a biomarker associated with progression to castration resistance, development of metastasis, and poor disease-specific survival in prostate cancer. However, high expression of AZGP1 cells in prostate cancer has been reported to increase proliferation and invasion. The exact role of AZGP1 in prostate cancer progression remains elusive. METHOD: AZGP1 knockout and overexpressing prostate cancer cells were generated using a lentiviral system. The effects of AZGP1 under- or over-expression in prostate cancer cells were evaluated by in vitro cell proliferation, migration, and invasion assays. Heterozygous AZGP1± mice were obtained from European Mouse Mutant Archive (EMMA), and prostate tissues from homozygous knockout male mice were collected at 2, 6 and 10 months for histological analysis. In vivo xenografts generated from AZGP1 under- or over-expressing prostate cancer cells were used to determine the role of AZGP1 in prostate cancer tumor growth, and subsequent proteomics analysis was conducted to elucidate the mechanisms of AZGP1 action in prostate cancer progression. AZGP1 expression and microvessel density were measured in human prostate cancer samples on a tissue microarray of 215 independent patient samples. RESULT: Neither the knockout nor overexpression of AZGP1 exhibited significant effects on prostate cancer cell proliferation, clonal growth, migration, or invasion in vitro. The prostates of AZGP1-/- mice initially appeared to have grossly normal morphology; however, we observed fibrosis in the periglandular stroma and higher blood vessel density in the mouse prostate by 6 months. In PC3 and DU145 mouse xenografts, over-expression of AZGP1 did not affect tumor growth. Instead, these tumors displayed decreased microvessel density compared to xenografts derived from PC3 and DU145 control cells, suggesting that AZGP1 functions to inhibit angiogenesis in prostate cancer. Proteomics profiling further indicated that, compared to control xenografts, AZGP1 overexpressing PC3 xenografts are enriched with angiogenesis pathway proteins, including YWHAZ, EPHA2, SERPINE1, and PDCD6, MMP9, GPX1, HSPB1, COL18A1, RNH1, and ANXA1. In vitro functional studies show that AZGP1 inhibits human umbilical vein endothelial cell proliferation, migration, tubular formation and branching. Additionally, tumor microarray analysis shows that AZGP1 expression is negatively correlated with blood vessel density in human prostate cancer tissues. CONCLUSION: AZGP1 is a negative regulator of angiogenesis, such that loss of AZGP1 promotes angiogenesis in prostate cancer. AZGP1 likely exerts heterotypical effects on cells in the tumor microenvironment, such as stromal and endothelial cells. This study sheds light on the anti-angiogenic characteristics of AZGP1 in the prostate and provides a rationale to target AZGP1 to inhibit prostate cancer progression.

Molecular profiling of clear cell renal cell carcinoma (ccRCC) tumors of patients in a clinical trial has identified distinct transcriptomic signatures with predictive value, yet data in non-clear cell variants (nccRCC) are lacking. We examined the transcriptional profiles of RCC tumors representing key molecular pathways, from a multi-institutional, real-world patient cohort, including ccRCC and centrally reviewed nccRCC samples. ccRCC had increased angiogenesis signature scores compared with the heterogeneous group of nccRCC tumors, while cell cycle, fatty acid oxidation/AMPK signaling, and fatty acid synthesis/pentose phosphate signature scores were increased in one or more nccRCC subtypes. Among both ccRCC and nccRCC tumors, T effector scores statistically correlated with increased immune cell infiltration and were more commonly associated with immunotherapy-related markers (PD-L1+/TMBhi/MSIhi). In conclusion, this study provides evidence of differential gene transcriptional profiles among ccRCC versus nccRCC tumors, providing insights for optimizing personalized and histology-specific therapeutic strategies for patients with advanced RCC.

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using assay for transposase-accessible chromatin coupled to self-transcribing active regulatory region (ATAC-STARR) sequencing. In addition to thousands of cis changes, we discover an unexpected number (∼10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie ∼37% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-in regulatory differences between species.

Although often located at a distance from their target gene promoters, enhancers are the primary genomic determinants of temporal and spatial transcriptional specificity in metazoans. Since the discovery of the first enhancer element in simian virus 40, there has been substantial interest in unraveling the mechanism(s) by which enhancers communicate with their partner promoters to ensure proper gene expression. These research efforts have benefited considerably from the application of increasingly sophisticated sequencing- and imaging-based approaches in conjunction with innovative (epi)genome-editing technologies; however, despite various proposed models, the principles of enhancer-promoter interaction have still not been fully elucidated. In this review, we provide an overview of recent progress in the eukaryotic gene transcription field pertaining to enhancer-promoter specificity. A better understanding of the mechanistic basis of lineage- and context-dependent enhancer-promoter engagement, along with the continued identification of functional enhancers, will provide key insights into the spatiotemporal control of gene expression that can reveal therapeutic opportunities for a range of enhancer-related diseases.

Neuroblastoma is the most common extracranial malignant tumor of childhood, accounting for 15% of all pediatric cancer deaths. Despite significant advances in our understanding of neuroblastoma biology, five-year survival rates for high-risk disease remain less than 50%, highlighting the importance of identifying novel therapeutic targets to combat the disease. MYCN amplification is the most frequent and predictive molecular aberration correlating with poor outcome in neuroblastoma. N-Myc is a short-lived protein primarily due to its rapid proteasomal degradation, a potentially exploitable vulnerability in neuroblastoma. AF1q is an oncoprotein with established roles in leukemia and solid tumor progression. It is normally expressed in brain and sympathetic neurons and has been postulated to play a part in neural differentiation. However, no role for AF1q in tumors of neural origin has been reported. In this study, we found AF1q to be a universal marker of neuroblastoma tumors. Silencing AF1q in neuroblastoma cells caused proteasomal degradation of N-Myc through Ras/ERK and AKT/GSK3β pathways, activated p53 and blocked cell cycle progression, culminating in cell death via the intrinsic apoptotic pathway. Moreover, silencing AF1q attenuated neuroblastoma tumorigenicity in vivo signifying AF1qs importance in neuroblastoma oncogenesis. Our findings reveal AF1q to be a novel regulator of N-Myc and potential therapeutic target in neuroblastoma.

Mechanisms underlying phenotypic divergence across species remain unresolved. In this issue of Cell Genomics, Hansen, Fong, et al.1 systematically dissect human and rhesus macaque gene expression divergence by screening tens of thousands of orthologous elements for enhancer activity in lymphoblastoid cell lines, revealing a much greater role for trans divergence at levels equal to those of cis effects, counter to the prevailing consensus in the field.

IL-17C is an epithelial cell-derived proinflammatory cytokine whose transcriptional regulation remains unclear. Analysis of the IL17C promoter region identified TCF4 as putative regulator, and siRNA knockdown of TCF4 in human keratinocytes (KCs) increased IL17C. IL-17C stimulation of KCs (along with IL-17A and TNF-α stimulation) decreased TCF4 and increased NFKBIZ and ZC3H12A expression in an IL-17RA/RE-dependent manner, thus creating a feedback loop. ZC3H12A (MCPIP1/Regnase-1), a transcriptional immune-response regulator, also increased following TCF4 siRNA knockdown, and siRNA knockdown of ZC3H12A decreased NFKBIZ, IL1B, IL36G, CCL20, and CXCL1, revealing a proinflammatory role for ZC3H12A. Examination of lesional skin from the KC-Tie2 inflammatory dermatitis mouse model identified decreases in TCF4 protein concomitant with increases in IL-17C and Zc3h12a that reversed following the genetic elimination of Il17c, Il17ra, and Il17re and improvement in the skin phenotype. Conversely, interference with Tcf4 in KC-Tie2 mouse skin increased Il17c and exacerbated the inflammatory skin phenotype. Together, these findings identify a role for TCF4 in the negative regulation of IL-17C, which, alone and with TNF-α and IL-17A, feed back to decrease TCF4 in an IL-17RA/RE-dependent manner. This loop is further amplified by IL-17C-TCF4 autocrine regulation of ZC3H12A and IL-17C regulation of NFKBIZ to promote self-sustaining skin inflammation.