Characterization of tumor heterogeneity through segmentation-free representation learning.

The interaction between tumors and their microenvironment is complex and heterogeneous. Recent developments in high-dimensional multiplexed imaging have revealed the spatial organization of tumor tissues at the molecular level. However, the discovery and thorough characterization of the tumor microenvironment (TME) remains challenging due to the scale and complexity of the images. Here, we propose a self-supervised representation learning framework, CANVAS, that enables discovery of novel types of TMEs. CANVAS is a vision transformer that directly takes high-dimensional multiplexed images and is trained using self-supervised masked image modeling. In contrast to traditional spatial analysis approaches which rely on cell segmentations, CANVAS is segmentation-free, utilizes pixel-level information, and retains local morphology and biomarker distribution information. This approach allows the model to distinguish subtle morphological differences, leading to precise separation and characterization of distinct TME signatures. We applied CANVAS to a lung tumor dataset and identified and validated a monocytic signature that is associated with poor prognosis.

MethNet: a robust approach to identify regulatory hubs and their distal targets from cancer data.

Aberrations in the capacity of DNA/chromatin modifiers and transcription factors to bind non-coding regions can lead to changes in gene regulation and impact disease phenotypes. However, identifying distal regulatory elements and connecting them with their target genes remains challenging. Here, we present MethNet, a pipeline that integrates large-scale DNA methylation and gene expression data across multiple cancers, to uncover cis regulatory elements (CREs) in a 1 Mb region around every promoter in the genome. MethNet identifies clusters of highly ranked CREs, referred to as 'hubs', which contribute to the regulation of multiple genes and significantly affect patient survival. Promoter-capture Hi-C confirmed that highly ranked associations involve physical interactions between CREs and their gene targets, and CRISPR interference based single-cell RNA Perturb-seq validated the functional impact of CREs. Thus, MethNet-identified CREs represent a valuable resource for unraveling complex mechanisms underlying gene expression, and for prioritizing the verification of predicted non-coding disease hotspots.

Genome-wide DNA methylation profiles and breast cancer among World Trade Center survivors.

Background: Increased incidence of cancer has been reported among World Trade Center (WTC)-exposed persons. Aberrant DNA methylation is a hallmark of cancer development. To date, only a few small studies have investigated the relationship between WTC exposure and DNA methylation. The main objective of this study was to assess the DNA methylation profiles of WTC-exposed community members who remained cancer free and those who developed breast cancer. Methods: WTC-exposed women were selected from the WTC Environmental Health Center clinic, with peripheral blood collected during routine clinical monitoring visits. The reference group was selected from the NYU Women's Health Study, a prospective cohort study with blood samples collected before 9 November 2001. The Infinium MethylationEPIC array was used for global DNA methylation profiling, with adjustments for cell type composition and other confounders. Annotated probes were used for biological pathway and network analysis. Results: A total of 64 WTC-exposed (32 cancer free and 32 with breast cancer) and 32 WTC-unexposed (16 cancer free and 16 with prediagnostic breast cancer) participants were included. Hypermethylated cytosine-phosphate-guanine probe sites (defined as ß > 0.8) were more common among WTC-exposed versus unexposed participants (14.3% vs. 4.5%, respectively, among the top 5000 cytosine-phosphate-guanine sites). Cancer-related pathways (e.g., human papillomavirus infection, cGMP-PKG) were overrepresented in WTC-exposed groups (breast cancer patients and cancer-free subjects). Compared to the unexposed breast cancer patients, 47 epigenetically dysregulated genes were identified among WTC-exposed breast cancers. These genes formed a network, including Wnt/ß-catenin signaling genes WNT4 and TCF7L2, and dysregulation of these genes contributes to cancer immune evasion. Conclusion: WTC exposure likely impacts DNA methylation and may predispose exposed individuals toward cancer development, possibly through an immune-mediated mechanism.

Inflammation in the tumor-adjacent lung as a predictor of clinical outcome in lung adenocarcinoma.

Approximately 30% of early-stage lung adenocarcinoma patients present with disease progression after successful surgical resection. Despite efforts of mapping the genetic landscape, there has been limited success in discovering predictive biomarkers of disease outcomes. Here we performed a systematic multi-omic assessment of 143 tumors and matched tumor-adjacent, histologically-normal lung tissue with long-term patient follow-up. Through histologic, mutational, and transcriptomic profiling of tumor and adjacent-normal tissue, we identified an inflammatory gene signature in tumor-adjacent tissue as the strongest clinical predictor of disease progression. Single-cell transcriptomic analysis demonstrated the progression-associated inflammatory signature was expressed in both immune and non-immune cells, and cell type-specific profiling in monocytes further improved outcome predictions. Additional analyses of tumor-adjacent transcriptomic data from The Cancer Genome Atlas validated the association of the inflammatory signature with worse outcomes across cancers. Collectively, our study suggests that molecular profiling of tumor-adjacent tissue can identify patients at high risk for disease progression.

Digital spatial profiling to predict recurrence in grade 3 stage I lung adenocarcinoma.

OBJECTIVE: Early-stage lung adenocarcinoma is treated with local therapy alone, although patients with grade 3 stage I lung adenocarcinoma have a 50% 5-year recurrence rate. Our objective is to determine if analysis of the tumor microenvironment can create a predictive model for recurrence. METHODS: Thirty-four patients with grade 3 stage I lung adenocarcinoma underwent surgical resection. Digital spatial profiling was used to perform genomic (n = 31) and proteomic (n = 34) analyses of pancytokeratin positive and negative tumor cells. K-means clustering was performed on the top 50 differential genes and top 20 differential proteins, with Kaplan-Meier recurrence curves based on patient clustering. External validation of high-expression genes was performed with Kaplan-Meier plotter. RESULTS: There were no significant clinicopathologic differences between patients who did (n = 14) and did not (n = 20) have recurrence. Median time to recurrence was 806 days; median follow-up with no recurrence was 2897 days. K-means clustering of pancytokeratin positive genes resulted in a model with a Kaplan-Meier curve with concordance index of 0.75. K-means clustering for pancytokeratin negative genes was less successful at differentiating recurrence (concordance index 0.6). Genes upregulated or downregulated for recurrence were externally validated using available public databases. Proteomic data did not reach statistical significance but did internally validate the genomic data described. CONCLUSIONS: Genomic difference in lung adenocarcinoma may be able to predict risk of recurrence. After further validation, stratifying patients by this risk may help guide who will benefit from adjuvant therapy.

MethNet: a robust approach to identify regulatory hubs and their distal targets in cancer.

Aberrations in the capacity of DNA/chromatin modifiers and transcription factors to bind non-coding regions can lead to changes in gene regulation and impact disease phenotypes. However, identifying distal regulatory elements and connecting them with their target genes remains challenging. Here, we present MethNet, a pipeline that integrates large-scale DNA methylation and gene expression data across multiple cancers, to uncover novel cis regulatory elements (CREs) in a 1Mb region around every promoter in the genome. MethNet identifies clusters of highly ranked CREs, referred to as 'hubs', which contribute to the regulation of multiple genes and significantly affect patient survival. Promoter-capture Hi-C confirmed that highly ranked associations involve physical interactions between CREs and their gene targets, and CRISPRi based scRNA Perturb-seq validated the functional impact of CREs. Thus, MethNet-identified CREs represent a valuable resource for unraveling complex mechanisms underlying gene expression, and for prioritizing the verification of predicted non-coding disease hotspots.

CCNE1 amplification among metastatic sites in patients with gynecologic high-grade serous carcinoma.

OBJECTIVE: We sought to characterize the variability of CCNE1 amplification among metastatic sites of CCNE1 amplified high grade serous carcinoma (HGSC) cases to investigate the feasibility of targeting this alteration for therapeutic purposes. METHODS: Patients with CCNE1 amplified HGSC who underwent surgical cytoreduction with metastatic sites were identified from institutional molecular profiling reports and a population of HGSC cases screened using digital droplet PCR (ddPCR). Cases with normal CCNE1 copy number were included as controls. Slides from metastatic sites were cut from formalin-fixed paraffin-embedded tissue blocks, dissected for tumor of > 50% purity, and underwent DNA extraction. CCNE1 copy number was determined by ddPCR. Tumor purity was confirmed with mutant TP53 allele fraction from targeted massively parallel sequencing. RESULTS: Four of 15 patients from an institutional database screened by ddPCR were found to have CCNE1 amplification. Three additional patients were identified from a query of institutional commercial clinical reports. Among these 7 CCNE1 amplified cases (2 uterine, 5 ovarian), 5 showed preservation of CCNE1 amplification (copy number > 5) among all metastatic sites. The remaining 2 cases had multiple metastatic sites without preserved CCNE1 amplification. Non-amplified cases had predominantly normal CCNE1 copy number across metastatic sites. CONCLUSIONS: CCNE1 amplification is an early genomic event in HGSC and is preserved in most metastatic sites suggesting a uniform response to pathway targeting therapies.

Dispersal dynamics of SARS-CoV-2 lineages during the first epidemic wave in New York City.

During the first phase of the COVID-19 epidemic, New York City rapidly became the epicenter of the pandemic in the United States. While molecular phylogenetic analyses have previously highlighted multiple introductions and a period of cryptic community transmission within New York City, little is known about the circulation of SARS-CoV-2 within and among its boroughs. We here perform phylogeographic investigations to gain insights into the circulation of viral lineages during the first months of the New York City outbreak. Our analyses describe the dispersal dynamics of viral lineages at the state and city levels, illustrating that peripheral samples likely correspond to distinct dispersal events originating from the main metropolitan city areas. In line with the high prevalence recorded in this area, our results highlight the relatively important role of the borough of Queens as a transmission hub associated with higher local circulation and dispersal of viral lineages toward the surrounding boroughs.

Somatic Focal Copy Number Gains of Noncoding Regions of Receptor Tyrosine Kinase Genes in Treatment-Resistant Epilepsy.

Epilepsy is a heterogenous group of disorders defined by recurrent seizure activity due to abnormal synchronized activity of neurons. A growing number of epilepsy cases are believed to be caused by genetic factors and copy number variants (CNV) contribute to up to 5% of epilepsy cases. However, CNVs in epilepsy are usually large deletions or duplications involving multiple neurodevelopmental genes. In patients who underwent seizure focus resection for treatment-resistant epilepsy, whole genome DNA methylation profiling identified 3 main clusters of which one showed strong association with receptor tyrosine kinase (RTK) genes. We identified focal copy number gains involving epidermal growth factor receptor (EGFR) and PDGFRA loci. The dysplastic neurons of cases with amplifications showed marked overexpression of EGFR and PDGFRA, while glial and endothelial cells were negative. Targeted sequencing of regulatory regions and DNA methylation analysis revealed that only enhancer regions of EGFR and gene promoter of PDGFRA were amplified, while coding regions did not show copy number abnormalities or somatic mutations. Somatic focal copy number gains of noncoding regulatory represent a previously unrecognized genetic driver in epilepsy and a mechanism of abnormal activation of RTK genes. Upregulated RTKs provide a potential avenue for therapy in seizure disorders.

Sequencing identifies multiple early introductions of SARS-CoV-2 to the New York City region.

Effective public response to a pandemic relies upon accurate measurement of the extent and dynamics of an outbreak. Viral genome sequencing has emerged as a powerful approach to link seemingly unrelated cases, and large-scale sequencing surveillance can inform on critical epidemiological parameters. Here, we report the analysis of 864 SARS-CoV-2 sequences from cases in the New York City metropolitan area during the COVID-19 outbreak in spring 2020. The majority of cases had no recent travel history or known exposure, and genetically linked cases were spread throughout the region. Comparison to global viral sequences showed that early transmission was most linked to cases from Europe. Our data are consistent with numerous seeds from multiple sources and a prolonged period of unrecognized community spreading. This work highlights the complementary role of genomic surveillance in addition to traditional epidemiological indicators.

Sequencing identifies multiple early introductions of SARS-CoV-2 to the New York City Region.

Effective public response to a pandemic relies upon accurate measurement of the extent and dynamics of an outbreak. Viral genome sequencing has emerged as a powerful approach to link seemingly unrelated cases, and large-scale sequencing surveillance can inform on critical epidemiological parameters. Here, we report the analysis of 864 SARS-CoV-2 sequences from cases in the New York City metropolitan area during the COVID-19 outbreak in Spring 2020. The majority of cases had no recent travel history or known exposure, and genetically linked cases were spread throughout the region. Comparison to global viral sequences showed that early transmission was most linked to cases from Europe. Our data are consistent with numerous seeds from multiple sources and a prolonged period of unrecognized community spreading. This work highlights the complementary role of genomic surveillance in addition to traditional epidemiological indicators.

Transposon insertion profiling by sequencing (TIPseq) for mapping LINE-1 insertions in the human genome.

BACKGROUND: Transposable elements make up a significant portion of the human genome. Accurately locating these mobile DNAs is vital to understand their role as a source of structural variation and somatic mutation. To this end, laboratories have developed strategies to selectively amplify or otherwise enrich transposable element insertion sites in genomic DNA. RESULTS: Here we describe a technique, Transposon Insertion Profiling by sequencing (TIPseq), to map Long INterspersed Element 1 (LINE-1, L1) retrotransposon insertions in the human genome. This method uses vectorette PCR to amplify species-specific L1 (L1PA1) insertion sites followed by paired-end Illumina sequencing. In addition to providing a step-by-step molecular biology protocol, we offer users a guide to our pipeline for data analysis, TIPseqHunter. Our recent studies in pancreatic and ovarian cancer demonstrate the ability of TIPseq to identify invariant (fixed), polymorphic (inherited variants), as well as somatically-acquired L1 insertions that distinguish cancer genomes from a patient's constitutional make-up. CONCLUSIONS: TIPseq provides an approach for amplifying evolutionarily young, active transposable element insertion sites from genomic DNA. Our rationale and variations on this protocol may be useful to those mapping L1 and other mobile elements in complex genomes.

Electrophoretic mobility shift assay analysis of NF-κB DNA binding.

The discovery and characterization of the nuclear factor-kappa B (NF-κB) family of transcription factors was predicated on the technical ability to detect protein binding to defined sequences of DNA. Proteins capable of binding to specific sequences of nucleic acid are detected through the use of the electrophoretic mobility shift assay (EMSA), also called a gel shift assay. While newer techniques, including chromatin immunoprecipitation (ChIP), are widely used to assess NF-κB binding to the promoters and enhancers of specific genes, the EMSA remains a powerful experimental tool to quickly test for the presence of NF-κB that is capable of binding DNA. In this way, the EMSA is a useful general readout of the activation state of the NF-κB pathway and an essential tool for the investigation of this important transcription factor family.

Proteasome inhibition by bortezomib increases IL-8 expression in androgen-independent prostate cancer cells: the role of IKKα.

Expression of the proinflammatory and proangiogenic chemokine IL-8, which is regulated at the transcriptional level by NF-κB, is constitutively increased in androgen-independent metastatic prostate cancer and correlates with poor prognosis. Inhibition of NF-κB-dependent transcription was used as an anticancer strategy for the development of the first clinically approved 26S proteasome inhibitor, bortezomib (BZ). Even though BZ has shown remarkable antitumor activity in hematological malignancies, it has been less effective in prostate cancer and other solid tumors; however, the mechanisms have not been fully understood. In this article, we report that proteasome inhibition by BZ unexpectedly increases IL-8 expression in androgen-independent prostate cancer PC3 and DU145 cells, whereas expression of other NF-κB-regulated genes is inhibited or unchanged. The BZ-increased IL-8 expression is associated with increased in vitro p65 NF-κB DNA binding activity and p65 recruitment to the endogenous IL-8 promoter. In addition, proteasome inhibition induces a nuclear accumulation of IκB kinase (IKK)α, and inhibition of IKKα enzymatic activity significantly attenuates the BZ-induced p65 recruitment to IL-8 promoter and IL-8 expression, demonstrating that the induced IL-8 expression is mediated, at least partly, by IKKα. Together, these data provide the first evidence, to our knowledge, for the gene-specific increase of IL-8 expression by proteasome inhibition in prostate cancer cells and suggest that targeting both IKKα and the proteasome may increase BZ effectiveness in treatment of androgen-independent prostate cancer.

Electrophoretic mobility shift assay analysis of NFκB transcriptional regulation by nuclear IκBα.

Transcription factor NFκB is a key regulator of genes involved in immune and inflammatory responses, as well as genes regulating cell proliferation and survival. In addition to many inflammatory disorders, NFκB is constitutively activated in a variety of human cancers and leukemia. Thus, inhibition of NFκB DNA binding activity represents an important therapeutic approach for disorders characterized by high levels of constitutive NFκB activity. We have previously shown that NFκB DNA binding activity is suppressed by the nuclear translocation and accumulation of IκBα, which is induced by inhibition of the 26S proteasome. In this chapter, we describe a protocol that uses small inhibitory RNA (si RNA) interference followed by electrophoretic mobility shift assay (EMSA) to analyze the regulation of NFκB DNA binding by nuclear IκBα induced by the proteasome inhibitor MG132. Using this protocol, we show that in human leukemia Hut-78 cells that exhibit high levels of NFκB DNA binding activity, MG132 induces nuclear translocation and accumulation of IκBα, which then specifically inhibits NFκB DNA binding. This protocol uses human leukemia Hut-78 cells; however, it can be easily adapted for other cells exhibiting high levels of constitutive NFκB DNA binding.

Chromatin immunoprecipitation analysis of NFκB transcriptional regulation by nuclear IκBα in human macrophages.

Transcription factor NFκB comprises a family of proteins that serve as crucial regulators of genes involved in host immune and inflammatory responses, cell survival, proliferation, and differentiation. Since transcription of NFκB-dependent genes is increased in numerous inflammatory disorders as well as in many types of cancer and leukemia, inhibition of NFκB-dependent transcription thus represents an important therapeutic target. We have previously shown that in human leukocytes, transcription of NFκB-dependent genes is inhibited by the nuclear translocation and accumulation of IκBα, which can be induced by an inhibitor of CRM1-dependent nuclear export, leptomycin B (LMB). In this chapter, we describe a protocol that uses chromatin immunoprecipitation (ChIP) to analyze the regulation of NFκB recruitment to NFκB-dependent promoters by nuclear IκBα induced by LMB. We show that in lipopolysaccharide (LPS)-stimulated human U-937 macrophages, recruitment of NFκB p65 and p50 proteins to NFκB-dependent promoters of IκBα and cIAP2 genes is suppressed by the LMB-induced nuclear IκBα. Even though in this study we use U-937 macrophages, this protocol should be readily adaptable to analyze the regulation of NFκB recruitment by nuclear IκBα also in other cell types.

Bortezomib induces nuclear translocation of IκBα resulting in gene-specific suppression of NF-κB--dependent transcription and induction of apoptosis in CTCL.

Cutaneous T-cell lymphoma (CTCL) is characterized by constitutive activation of nuclear factor κB (NF-κB), which plays a crucial role in the survival of CTCL cells and their resistance to apoptosis. NF-κB activity in CTCL is inhibited by the proteasome inhibitor bortezomib; however, the mechanisms remained unknown. In this study, we investigated mechanisms by which bortezomib suppresses NF-κB activity in CTCL Hut-78 cells. We demonstrate that bortezomib and MG132 suppress NF-κB activity in Hut-78 cells by a novel mechanism that consists of inducing nuclear translocation and accumulation of IκBα (nuclear factor of kappa light polypeptide gene enhancer in B-cells inhibitor, alpha), which then associates with NF-κB p65 and p50 in the nucleus and inhibits NF-κB DNA binding activity. Surprisingly, however, while expression of NF-κB-dependent antiapoptotic genes cIAP1 and cIAP2 is inhibited by bortezomib, expression of Bcl-2 is not suppressed. Chromatin immunoprecipitation indicated that cIAP1 and cIAP2 promoters are occupied by NF-κB p65/50 heterodimers, whereas Bcl-2 promoter is occupied predominantly by p50/50 homodimers. Collectively, our data reveal a novel mechanism of bortezomib function in CTCL and suggest that the inhibition of NF-κB-dependent gene expression by bortezomib is gene specific and depends on the subunit composition of NF-κB dimers recruited to NF-κB-responsive promoters.

Gene-specific repression of proinflammatory cytokines in stimulated human macrophages by nuclear IκBα.

We have previously shown that increased nuclear accumulation of IkappaBalpha inhibits NF-kappaB activity and induces apoptosis in human leukocytes. In this study, we wanted to explore the possibility that the nucleocytoplasmic distribution of IkappaBalpha can be used as a therapeutic target for the regulation of NF-kappaB-dependent cytokine synthesis. Treatment of LPS-stimulated human U937 macrophages with an inhibitor of chromosome region maintenance 1-dependent nuclear export, leptomycin B, resulted in the increased nuclear accumulation of IkappaBalpha and inhibition of NF-kappaB DNA binding activity, caused by the nuclear IkappaBalpha-p65 NF-kappaB interaction. Surprisingly, however, whereas mRNA expression and cellular release of TNF-alpha, the beta form of pro-IL-1 (IL-1beta), and IL-6 were inhibited by the leptomycin B-induced nuclear IkappaBalpha, IL-8 mRNA expression and cellular release were not significantly affected. Analysis of in vivo recruitment of p65 NF-kappaB to NF-kappaB-regulated promoters by chromatin immunoprecipitation in U937 cells and human PBMCs indicated that although the p65 recruitment to TNF-alpha, IL-1beta, and IL-6 promoters was inhibited by the nuclear IkappaBalpha, p65 recruitment to IL-8 promoter was not repressed. Chromatin immunoprecipitation analyses using IkappaBalpha and S536 phosphospecific p65 NF-kappaB Abs demonstrated that although the newly synthesized IkappaBalpha induced by postinduction repression is recruited to TNF-alpha, IL-1beta, and IL-6 promoters but not to the IL-8 promoter, S536-phosphorylated p65 is recruited to IL-8 promoter, but not to TNF-alpha, IL-1beta, or IL-6 promoters. Together, these data indicate that the inhibition of NF-kappaB-dependent transcription by nuclear IkappaBalpha in LPS-stimulated macrophages is gene specific and depends on the S536 phosphorylation status of the recruited p65 NF-kappaB.

Proteasome inhibitors induce apoptosis of prostate cancer cells by inducing nuclear translocation of IkappaBalpha.

Proteasome inhibitors are known to suppress the proteasome-mediated degradation of IkappaBalpha in stimulated cells. This results in the cytoplasmic retention of NFkappaB and its reduced nuclear transcriptional activity. In this study, we show that in the metastatic prostate cancer cells, the proteasome inhibitors exhibit a novel, previously unrecognized effect: they increase the cellular levels of IkappaBalpha, which then translocates to the nucleus, associates with the nuclear p65 NFkappaB, thus inhibiting the constitutive NFkappaB DNA binding activity and inducing apoptosis. The proteasome inhibition-induced nuclear translocation of IkappaBalpha is dependent on de novo protein synthesis, occurs also in other cell types, and does not require IkappaBalpha phosphorylation on Ser-32. Since NFkappaB activity is constitutively increased in many human cancers as well as in inflammatory disorders, the proteasome inhibition-induced nuclear translocation of IkappaBalpha could thus provide a new therapeutic strategy aimed at the specific inhibition of NFkappaB activity by the nuclear IkappaBalpha.