Aging limits stemness and tumorigenesis in the lung by reprogramming iron homeostasis.

Aging is associated with a decline in the number and fitness of adult stem cells 1-4 . Aging-associated loss of stemness is posited to suppress tumorigenesis 5,6 , but this hypothesis has not been tested in vivo . Here, using physiologically aged autochthonous genetically engineered mouse models and primary cells 7,8 , we demonstrate aging suppresses lung cancer initiation and progression by degrading stemness of the alveolar cell of origin. This phenotype is underpinned by aging-associated induction of the transcription factor NUPR1 and its downstream target lipocalin-2 in the cell of origin in mice and humans, leading to a functional iron insufficiency in the aged cells. Genetic inactivation of the NUPR1-lipocalin-2 axis or iron supplementation rescue stemness and promote tumorigenic potential of aged alveolar cells. Conversely, targeting the NUPR1- lipocalin-2 axis is detrimental to young alveolar cells via induction of ferroptosis. We find that aging-associated DNA hypomethylation at specific enhancer sites associates with elevated NUPR1 expression, which is recapitulated in young alveolar cells by inhibition of DNA methylation. We uncover that aging drives a functional iron insufficiency, which leads to loss of stemness and tumorigenesis, but promotes resistance to ferroptosis. These findings have significant implications for the therapeutic modulation of cellular iron homeostasis in regenerative medicine and in cancer prevention. Furthermore, our findings are consistent with a model whereby most human cancers initiate in young individuals, revealing a critical window for such cancer prevention efforts.

TAD hierarchy restricts poised LTR activation and loss of TAD hierarchy promotes LTR co-option in cancer.

Transposable elements (TEs) are abundant in the human genome, and they provide the sources for genetic and functional diversity. The regulation of TEs expression and their functional consequences in physiological conditions and cancer development remain to be fully elucidated. Previous studies suggested TEs are repressed by DNA methylation and chromatin modifications. The effect of 3D chromatin topology on TE regulation remains elusive. Here, by integrating transcriptome and 3D genome architecture studies, we showed that haploinsufficient loss of NIPBL selectively activates alternative promoters at the long terminal repeats (LTRs) of the TE subclasses. This activation occurs through the reorganization of topologically associating domain (TAD) hierarchical structures and recruitment of proximal enhancers. These observations indicate that TAD hierarchy restricts transcriptional activation of LTRs that already possess open chromatin features. In cancer, perturbation of the hierarchical chromatin topology can lead to co-option of LTRs as functional alternative promoters in a context-dependent manner and drive aberrant transcriptional activation of novel oncogenes and other divergent transcripts. These data uncovered a new layer of regulatory mechanism of TE expression beyond DNA and chromatin modification in human genome. They also posit the TAD hierarchy dysregulation as a novel mechanism for alternative promoter-mediated oncogene activation and transcriptional diversity in cancer, which may be exploited therapeutically.

An alternative NURF complex sustains acute myeloid leukemia by regulating the accessibility of insulator regions.

Efficient treatment of acute myeloid leukemia (AML) patients remains a challenge despite recent therapeutic advances. Here, using a CRISPRi screen targeting chromatin factors, we identified the nucleosome-remodeling factor (NURF) subunit BPTF as an essential regulator of AML cell survival. We demonstrate that BPTF forms an alternative NURF chromatin remodeling complex with SMARCA5 and BAP18, which regulates the accessibility of a large set of insulator regions in leukemic cells. This ensures efficient CTCF binding and boundary formation between topologically associated domains that is essential for maintaining the leukemic transcriptional programs. We also demonstrate that the well-studied PHD2-BROMO chromatin reader domains of BPTF, while contributing to complex recruitment to chromatin, are dispensable for leukemic cell growth. Taken together, our results uncover how the alternative NURF complex contributes to leukemia and provide a rationale for its targeting in AML.

Exportin 1 inhibition prevents neuroendocrine transformation through SOX2 down-regulation in lung and prostate cancers.

In lung and prostate adenocarcinomas, neuroendocrine (NE) transformation to an aggressive derivative resembling small cell lung cancer (SCLC) is associated with poor prognosis. We previously described dependency of SCLC on the nuclear transporter exportin 1. Here, we explored the role of exportin 1 in NE transformation. We observed up-regulated exportin 1 in lung and prostate pretransformation adenocarcinomas. Exportin 1 was up-regulated after genetic inactivation of TP53 and RB1 in lung and prostate adenocarcinoma cell lines, accompanied by increased sensitivity to the exportin 1 inhibitor selinexor in vitro. Exportin 1 inhibition prevented NE transformation in different TP53/RB1-inactivated prostate adenocarcinoma xenograft models that acquire NE features upon treatment with the aromatase inhibitor enzalutamide and extended response to the EGFR inhibitor osimertinib in a lung cancer transformation patient-derived xenograft (PDX) model exhibiting combined adenocarcinoma/SCLC histology. Ectopic SOX2 expression restored the enzalutamide-promoted NE phenotype on adenocarcinoma-to-NE transformation xenograft models despite selinexor treatment. Selinexor sensitized NE-transformed lung and prostate small cell carcinoma PDXs to standard cytotoxics. Together, these data nominate exportin 1 inhibition as a potential therapeutic target to constrain lineage plasticity and prevent or treat NE transformation in lung and prostate adenocarcinoma.

ChromaFold predicts the 3D contact map from single-cell chromatin accessibility.

The identification of cell-type-specific 3D chromatin interactions between regulatory elements can help to decipher gene regulation and to interpret the function of disease-associated non-coding variants. However, current chromosome conformation capture (3C) technologies are unable to resolve interactions at this resolution when only small numbers of cells are available as input. We therefore present ChromaFold, a deep learning model that predicts 3D contact maps and regulatory interactions from single-cell ATAC sequencing (scATAC-seq) data alone. ChromaFold uses pseudobulk chromatin accessibility, co-accessibility profiles across metacells, and predicted CTCF motif tracks as input features and employs a lightweight architecture to enable training on standard GPUs. Once trained on paired scATAC-seq and Hi-C data in human cell lines and tissues, ChromaFold can accurately predict both the 3D contact map and peak-level interactions across diverse human and mouse test cell types. In benchmarking against a recent deep learning method that uses bulk ATAC-seq, DNA sequence, and CTCF ChIP-seq to make cell-type-specific predictions, ChromaFold yields superior prediction performance when including CTCF ChIP-seq data as an input and comparable performance without. Finally, fine-tuning ChromaFold on paired scATAC-seq and Hi-C in a complex tissue enables deconvolution of chromatin interactions across cell subpopulations. ChromaFold thus achieves state-of-the-art prediction of 3D contact maps and regulatory interactions using scATAC-seq alone as input data, enabling accurate inference of cell-type-specific interactions in settings where 3C-based assays are infeasible.

TET2 guards against unchecked BATF3-induced CAR T cell expansion.

Further advances in cell engineering are needed to increase the efficacy of chimeric antigen receptor (CAR) and other T cell-based therapies1-5. As T cell differentiation and functional states are associated with distinct epigenetic profiles6,7, we hypothesized that epigenetic programming may provide a means to improve CAR T cell performance. Targeting the gene that encodes the epigenetic regulator ten-eleven translocation 2 (TET2)8 presents an interesting opportunity as its loss may enhance T cell memory9,10, albeit not cause malignancy9,11,12. Here we show that disruption of TET2 enhances T cell-mediated tumour rejection in leukaemia and prostate cancer models. However, loss of TET2 also enables antigen-independent CAR T cell clonal expansions that may eventually result in prominent systemic tissue infiltration. These clonal proliferations require biallelic TET2 disruption and sustained expression of the AP-1 factor BATF3 to drive a MYC-dependent proliferative program. This proliferative state is associated with reduced effector function that differs from both canonical T cell memory13,14 and exhaustion15,16 states, and is prone to the acquisition of secondary somatic mutations, establishing TET2 as a guardian against BATF3-induced CAR T cell proliferation and ensuing genomic instability. Our findings illustrate the potential of epigenetic programming to enhance T cell immunity but highlight the risk of unleashing unchecked proliferative responses.

Targeting Lysine-Specific Demethylase 1 Rescues Major Histocompatibility Complex Class I Antigen Presentation and Overcomes Programmed Death-Ligand 1 Blockade Resistance in SCLC.

INTRODUCTION: SCLC is a highly aggressive neuroendocrine tumor that is characterized by early acquired therapeutic resistance and modest benefit from immune checkpoint blockade (ICB). Repression of the major histocompatibility complex class I (MHC-I) represents a key mechanism driving resistance to T cell-based immunotherapies. METHODS: We evaluated the role of the lysine-specific demethylase 1 (LSD1) as a determinant of MHC-I expression, functional antigen presentation, and immune activation in SCLC in vitro and in vivo through evaluation of both human SCLC cell lines and immunocompetent mouse models. RESULTS: We found that targeted inhibition of LSD1 in SCLC restores MHC-I cell surface expression and transcriptionally activates genes encoding the antigen presentation pathway. LSD1 inhibition further activates interferon signaling, induces tumor-intrinsic immunogenicity, and sensitizes SCLC cells to MHC-I-restricted T cell cytolysis. Combination of LSD1 inhibitor with ICB augments the antitumor immune response in refractory SCLC models. Together, these data define a role for LSD1 as a potent regulator of MHC-I antigen presentation and provide rationale for combinatory use of LSD1 inhibitors with ICB to improve therapeutic response in SCLC. CONCLUSIONS: Epigenetic silencing of MHC-I in SCLC contributes to its poor response to ICB. Our study identifies a previously uncharacterized role for LSD1 as a regulator of MHC-I antigen presentation in SCLC. LSD1 inhibition enables MHC-I-restricted T cell cytolysis, induces immune activation, and augments the antitumor immune response to ICB in SCLC.

WEE1 inhibition enhances the antitumor immune response to PD-L1 blockade by the concomitant activation of STING and STAT1 pathways in SCLC.

Small cell lung cancers (SCLCs) have high mutational burden but are relatively unresponsive to immune checkpoint blockade (ICB). Using SCLC models, we demonstrate that inhibition of WEE1, a G2/M checkpoint regulator induced by DNA damage, activates the STING-TBK1-IRF3 pathway, which increases type I interferons (IFN-α and IFN-ß) and pro-inflammatory chemokines (CXCL10 and CCL5), facilitating an immune response via CD8+ cytotoxic T cell infiltration. We further show that WEE1 inhibition concomitantly activates the STAT1 pathway, increasing IFN-γ and PD-L1 expression. Consistent with these findings, combined WEE1 inhibition (AZD1775) and PD-L1 blockade causes remarkable tumor regression, activation of type I and II interferon pathways, and infiltration of cytotoxic T cells in multiple immunocompetent SCLC genetically engineered mouse models, including an aggressive model with stabilized MYC. Our study demonstrates cell-autonomous and immune-stimulating activity of WEE1 inhibition in SCLC models. Combined inhibition of WEE1 plus PD-L1 blockade represents a promising immunotherapeutic approach in SCLC.

Inhibition of XPO1 Sensitizes Small Cell Lung Cancer to First- and Second-Line Chemotherapy.

Small cell lung cancer (SCLC) is an aggressive malignancy characterized by early metastasis and extreme lethality. The backbone of SCLC treatment over the past several decades has been platinum-based doublet chemotherapy, with the recent addition of immunotherapy providing modest benefits in a subset of patients. However, nearly all patients treated with systemic therapy quickly develop resistant disease, and there is an absence of effective therapies for recurrent and progressive disease. Here we conducted CRISPR-Cas9 screens using a druggable genome library in multiple SCLC cell lines representing distinct molecular subtypes. This screen nominated exportin-1, encoded by XPO1, as a therapeutic target. XPO1 was highly and ubiquitously expressed in SCLC relative to other lung cancer histologies and other tumor types. XPO1 knockout enhanced chemosensitivity, and exportin-1 inhibition demonstrated synergy with both first- and second-line chemotherapy. The small molecule exportin-1 inhibitor selinexor in combination with cisplatin or irinotecan dramatically inhibited tumor growth in chemonaïve and chemorelapsed SCLC patient-derived xenografts, respectively. Together these data identify exportin-1 as a promising therapeutic target in SCLC, with the potential to markedly augment the efficacy of cytotoxic agents commonly used in treating this disease. SIGNIFICANCE: CRISPR-Cas9 screening nominates exportin-1 as a therapeutic target in SCLC, and exportin-1 inhibition enhances chemotherapy efficacy in patient-derived xenografts, providing a novel therapeutic opportunity in this disease.

Comprehensive molecular characterization of lung tumors implicates AKT and MYC signaling in adenocarcinoma to squamous cell transdifferentiation.

BACKGROUND: Lineage plasticity, the ability to transdifferentiate among distinct phenotypic identities, facilitates therapeutic resistance in cancer. In lung adenocarcinomas (LUADs), this phenomenon includes small cell and squamous cell (LUSC) histologic transformation in the context of acquired resistance to targeted inhibition of driver mutations. LUAD-to-LUSC transdifferentiation, occurring in up to 9% of EGFR-mutant patients relapsed on osimertinib, is associated with notably poor prognosis. We hypothesized that multi-parameter profiling of the components of mixed histology (LUAD/LUSC) tumors could provide insight into factors licensing lineage plasticity between these histologies. METHODS: We performed genomic, epigenomics, transcriptomics and protein analyses of microdissected LUAD and LUSC components from mixed histology tumors, pre-/post-transformation tumors and reference non-transformed LUAD and LUSC samples. We validated our findings through genetic manipulation of preclinical models in vitro and in vivo and performed patient-derived xenograft (PDX) treatments to validate potential therapeutic targets in a LUAD PDX model acquiring LUSC features after osimertinib treatment. RESULTS: Our data suggest that LUSC transdifferentiation is primarily driven by transcriptional reprogramming rather than mutational events. We observed consistent relative upregulation of PI3K/AKT, MYC and PRC2 pathway genes. Concurrent activation of PI3K/AKT and MYC induced squamous features in EGFR-mutant LUAD preclinical models. Pharmacologic inhibition of EZH1/2 in combination with osimertinib prevented relapse with squamous-features in an EGFR-mutant patient-derived xenograft model, and inhibition of EZH1/2 or PI3K/AKT signaling re-sensitized resistant squamous-like tumors to osimertinib. CONCLUSIONS: Our findings provide the first comprehensive molecular characterization of LUSC transdifferentiation, suggesting putative drivers and potential therapeutic targets to constrain or prevent lineage plasticity.

Multiomic Analysis of Lung Tumors Defines Pathways Activated in Neuroendocrine Transformation.

Lineage plasticity is implicated in treatment resistance in multiple cancers. In lung adenocarcinomas (LUAD) amenable to targeted therapy, transformation to small cell lung cancer (SCLC) is a recognized resistance mechanism. Defining molecular mechanisms of neuroendocrine (NE) transformation in lung cancer has been limited by a paucity of pre/posttransformation clinical samples. Detailed genomic, epigenomic, transcriptomic, and protein characterization of combined LUAD/SCLC tumors, as well as pre/posttransformation samples, supports that NE transformation is primarily driven by transcriptional reprogramming rather than mutational events. We identify genomic contexts in which NE transformation is favored, including frequent loss of the 3p chromosome arm. We observed enhanced expression of genes involved in the PRC2 complex and PI3K/AKT and NOTCH pathways. Pharmacologic inhibition of the PI3K/AKT pathway delayed tumor growth and NE transformation in an EGFR-mutant patient-derived xenograft model. Our findings define a novel landscape of potential drivers and therapeutic vulnerabilities of NE transformation in lung cancer. SIGNIFICANCE: The difficulty in collection of transformation samples has precluded the performance of molecular analyses, and thus little is known about the lineage plasticity mechanisms leading to LUAD-to-SCLC transformation. Here, we describe biological pathways dysregulated upon transformation and identify potential predictors and potential therapeutic vulnerabilities of NE transformation in the lung. See related commentary by Meador and Lovly, p. 2962. This article is highlighted in the In This Issue feature, p. 2945.

HiC-DC+ enables systematic 3D interaction calls and differential analysis for Hi-C and HiChIP.

Recent genome-wide chromosome conformation capture assays such as Hi-C and HiChIP have vastly expanded the resolution and throughput with which we can study 3D genomic architecture and function. Here, we present HiC-DC+, a software tool for Hi-C/HiChIP interaction calling and differential analysis using an efficient implementation of the HiC-DC statistical framework. HiC-DC+ integrates with popular preprocessing and visualization tools and includes topologically associating domain (TAD) and A/B compartment callers. We found that HiC-DC+ can more accurately identify enhancer-promoter interactions in H3K27ac HiChIP, as validated by CRISPRi-FlowFISH experiments, compared to existing methods. Differential HiC-DC+ analyses of published HiChIP and Hi-C data sets in settings of cellular differentiation and cohesin perturbation systematically and quantitatively recovers biological findings, including enhancer hubs, TAD aggregation, and the relationship between promoter-enhancer loop dynamics and gene expression changes. HiC-DC+ therefore provides a principled statistical analysis tool to empower genome-wide studies of 3D chromatin architecture and function.

Complex-dependent histone acetyltransferase activity of KAT8 determines its role in transcription and cellular homeostasis.

Acetylation of lysine 16 on histone H4 (H4K16ac) is catalyzed by histone acetyltransferase KAT8 and can prevent chromatin compaction in vitro. Although extensively studied in Drosophila, the functions of H4K16ac and two KAT8-containing protein complexes (NSL and MSL) are not well understood in mammals. Here, we demonstrate a surprising complex-dependent activity of KAT8: it catalyzes H4K5ac and H4K8ac as part of the NSL complex, whereas it catalyzes the bulk of H4K16ac as part of the MSL complex. Furthermore, we show that MSL complex proteins and H4K16ac are not required for cell proliferation and chromatin accessibility, whereas the NSL complex is essential for cell survival, as it stimulates transcription initiation at the promoters of housekeeping genes. In summary, we show that KAT8 switches catalytic activity and function depending on its associated proteins and that, when in the NSL complex, it catalyzes H4K5ac and H4K8ac required for the expression of essential genes.

Fostering bioinformatics education through skill development of professors: Big Genomic Data Skills Training for Professors.

Bioinformatics has become an indispensable part of life science over the past 2 decades. However, bioinformatics education is not well integrated at the undergraduate level, especially in liberal arts colleges and regional universities in the United States. One significant obstacle pointed out by the Network for Integrating Bioinformatics into Life Sciences Education is the lack of faculty in the bioinformatics area. Most current life science professors did not acquire bioinformatics analysis skills during their own training. Consequently, a great number of undergraduate and graduate students do not get the chance to learn bioinformatics or computational biology skills within a structured curriculum during their education. To address this gap, we developed a module-based, week-long short course to train small college and regional university professors with essential bioinformatics skills. The bioinformatics modules were built to be adapted by the professor-trainees afterward and used in their own classes. All the course materials can be accessed at https://github.com/TheJacksonLaboratory/JAXBD2K-ShortCourse.

The cis conformation of proline leads to weaker binding of a p53 peptide to MDM2 compared to trans.

The cis and trans conformations of the Xaa-Pro (Xaa: any amino acid) peptide bond are thermodynamically stable while other peptide bonds strongly prefer trans. The effect of proline cis-trans isomerization on protein binding has not been thoroughly investigated. In this study, computer simulations were used to calculate the absolute binding affinity for a p53 peptide (residues 17-29) to MDM2 for both cis and trans isomers of the p53 proline in position 27. Results show that the cis isomer of p53(17-29) binds more weakly to MDM2 than the trans isomer, and that this is primarily due to the difference in the free energy cost associated with the loss of conformational entropy of p53(17-29) when it binds to MDM2. The population of cis p53(17-29) was estimated to be 0.8% of the total population in the bound state. The stronger binding of trans p53(17-29) to MDM2 compared to cis may leave a minimal level of p53 available to respond to cellular stress. This study demonstrates that it is feasible to estimate the absolute binding affinity for an intrinsically disordered protein fragment binding to an ordered protein that are in good agreement with experimental results.

Impact of the K24N mutation on the transactivation domain of p53 and its binding to murine double-minute clone 2.

The level of the p53 transcription factor is negatively regulated by the E3 ubiquitin ligase murine double-minute clone 2 (MDM2). The interaction between p53 and MDM2 is essential for the maintenance of genomic integrity for most eukaryotes. Previous structural studies revealed that MDM2 binds to p53 transactivation domain (p53TAD) from residues 17 to 29. The K24N mutation of p53TAD changes a lysine at position 24 to an asparagine. This mutation occurs naturally in the bovine family and is also found in a rare form of human gestational cancer called choriocarcinoma. In this study, we have investigated how the K24N mutation affects the affinity, structure, and dynamics of p53TAD binding to MDM2. Nuclear magnetic resonance studies of p53TAD show that the K24N mutant is more flexible and has less transient helical secondary structure than the wild type. Isothermal titration calorimetry measurements demonstrate that these changes in structure and dynamics do not significantly change the binding affinity for p53TAD-MDM2. Finally, free-energy perturbation and standard molecular dynamic simulations suggest the negligible affinity change is due to a compensating interaction energy between the K24N mutant and the MDM2 when it is bound. Overall, the data suggest that the K24N-MDM2 complex is able to, at least partly, compensate for an increase in the conformational entropy in unbound K24N with an increase in the bound-state electrostatic interaction energy.

Effect of 1-methyladenine on double-helical DNA structures.

Methylation at the N1 site of adenine leads to the formation of cytotoxic 1-methyladenine (m1A). Since the N1 site of adenine is involved in the hydrogen bonding of T.A and A.T Watson-Crick base pairs, it is expected that the pairing interactions will be disrupted upon 1-methylation. In this study, high-resolution NMR investigations were performed to determine the effect of m1A on double-helical DNA structures. Interestingly, instead of disrupting hydrogen bonding, we found that 1-methylation altered the T.A Watson-Crick base pair to T(anti).m1A(syn) Hoogsteen base pair, providing insights into the observed differences in AlkB-repair efficiency between dsDNA and ssDNA.