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.

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.

Selective adaptation of SARS-CoV-2 Omicron under booster vaccine pressure: a multicentre observational study.

BACKGROUND: High rates of vaccination and natural infection drive immunity and redirect selective viral adaptation. Updated boosters are installed to cope with drifted viruses, yet data on adaptive evolution under increasing immune pressure in a real-world situation are lacking. METHODS: Cross-sectional study to characterise SARS-CoV-2 mutational dynamics and selective adaptation over >1 year in relation to vaccine status, viral phylogenetics, and associated clinical and demographic variables. FINDINGS: The study of >5400 SARS-CoV-2 infections between July 2021 and August 2022 in metropolitan New York portrayed the evolutionary transition from Delta to Omicron BA.1-BA.5 variants. Booster vaccinations were implemented during the Delta wave, yet booster breakthrough infections and SARS-CoV-2 re-infections were almost exclusive to Omicron. In adjusted logistic regression analyses, BA.1, BA.2, and BA.5 had a significant growth advantage over co-occurring lineages in the boosted population, unlike BA.2.12.1 or BA.4. Selection pressure by booster shots translated into diffuse adaptive evolution in Delta spike, contrasting with strong, receptor-binding motif-focused adaptive evolution in BA.2-BA.5 spike (Fisher Exact tests; non-synonymous/synonymous mutation rates per site). Convergent evolution has become common in Omicron, engaging spike positions crucial for immune escape, receptor binding, or cleavage. INTERPRETATION: Booster shots are required to cope with gaps in immunity. Their discriminative immune pressure contributes to their effectiveness but also requires monitoring of selective viral adaptation processes. Omicron BA.2 and BA.5 had a selective advantage under booster vaccination pressure, contributing to the evolution of BA.2 and BA.5 sublineages and recombinant forms that predominate in 2023. FUNDING: The study was supported by NYU institutional funds and partly by the Cancer Center Support Grant P30CA016087 at the Laura and Isaac Perlmutter Cancer Center.

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.

Single-cell RNA sequencing reveals the effects of chemotherapy on human pancreatic adenocarcinoma and its tumor microenvironment.

The tumor microenvironment (TME) in pancreatic ductal adenocarcinoma (PDAC) is a complex ecosystem that drives tumor progression; however, in-depth single cell characterization of the PDAC TME and its role in response to therapy is lacking. Here, we perform single-cell RNA sequencing on freshly collected human PDAC samples either before or after chemotherapy. Overall, we find a heterogeneous mixture of basal and classical cancer cell subtypes, along with distinct cancer-associated fibroblast and macrophage subpopulations. Strikingly, classical and basal-like cancer cells exhibit similar transcriptional responses to chemotherapy and do not demonstrate a shift towards a basal-like transcriptional program among treated samples. We observe decreased ligand-receptor interactions in treated samples, particularly between TIGIT on CD8 + T cells and its receptor on cancer cells, and identify TIGIT as the major inhibitory checkpoint molecule of CD8 + T cells. Our results suggest that chemotherapy profoundly impacts the PDAC TME and may promote resistance to immunotherapy.

Alterations in the cutaneous microbiome of patients with psoriasis and psoriatic arthritis reveal similarities between non-lesional and lesional skin.

OBJECTIVES: To investigate the cutaneous microbiome spanning the entire psoriatic disease spectrum, and to evaluate distinguishing features of psoriasis (PsO) and psoriatic arthritis (PsA). METHODS: Skin swabs were collected from upper and lower extremities of healthy individuals and patients with PsO and PsA. Psoriatic patients contributed both lesional (L) and contralateral non-lesional (NL) samples. Microbiota were analysed using 16S rRNA sequencing. RESULTS: Compared with healthy skin, alpha diversity in psoriatic NL and L skin was significantly reduced (p<0.05) and samples clustered separately in plots of beta diversity (p<0.05). Kocuria and Cutibacterium were enriched in healthy subjects, while Staphylococcus was enriched in psoriatic disease. Microbe-microbe association networks revealed a higher degree of similarity between psoriatic NL and L skin compared with healthy skin despite the absence of clinically evident inflammation. Moreover, the relative abundance of Corynebacterium was higher in NL PsA samples compared with NL PsO samples (p<0.05), potentially serving as a biomarker for disease progression. CONCLUSIONS: These findings show differences in diversity, bacterial composition and microbe-microbe interactions between healthy and psoriatic skin, both L and NL. We further identified bacterial biomarkers that differentiate disease phenotypes, which could potentially aid in predicting the transition from PsO to PsA.

Clinical and genomic signatures of SARS-CoV-2 Delta breakthrough infections in New York.

BACKGROUND: In 2021, Delta became the predominant SARS-CoV-2 variant worldwide. While vaccines have effectively prevented COVID-19 hospitalization and death, vaccine breakthrough infections increasingly occurred. The precise role of clinical and genomic determinants in Delta infections is not known, and whether they contributed to increased rates of breakthrough infections compared to unvaccinated controls. METHODS: We studied SARS-CoV-2 variant distribution, dynamics, and adaptive selection over time in relation to vaccine status, phylogenetic relatedness of viruses, full genome mutation profiles, and associated clinical and demographic parameters. FINDINGS: We show a steep and near-complete replacement of circulating variants with Delta between May and August 2021 in metropolitan New York. We observed an increase of the Delta sublineage AY.25 (14% in vaccinated, 7% in unvaccinated), its spike mutation S112L, and AY.44 (8% in vaccinated, 2% in unvaccinated) with its nsp12 mutation F192V in breakthroughs. Delta infections were associated with younger age and lower hospitalization rates than Alpha. Delta breakthrough infections increased significantly with time since vaccination, and, after adjusting for confounders, they rose at similar rates as in unvaccinated individuals. INTERPRETATION: We observed a modest adaptation of Delta genomes in breakthrough infections in New York, suggesting an improved genomic framework to support Delta's epidemic growth in times of waning vaccine protection despite limited impact on vaccine escape. FUNDING: The study was supported by NYU institutional funds. The NYULH Genome Technology Center is partially supported by the Cancer Center Support Grant P30CA016087 at the Laura and Isaac Perlmutter Cancer Center.

Global DNA Methylation Profiles in Peripheral Blood of WTC-Exposed Community Members with Breast Cancer.

Breast cancer represents the most common cancer diagnosis among World Trade Center (WTC)-exposed community members, residents, and cleanup workers enrolled in the WTC Environmental Health Center (WTC EHC). The primary aims of this study were (1) to compare blood DNA methylation profiles of WTC-exposed community members with breast cancer and WTC-unexposed pre-diagnostic breast cancer blood samples, and (2) to compare the DNA methylation differences among the WTC EHC breast cancer cases and WTC-exposed cancer-free controls. Gene pathway enrichment analyses were further conducted. There were significant differences in DNA methylation between WTC-exposed breast cancer cases and unexposed prediagnostic breast cancer cases. The top differentially methylated genes were Intraflagellar Transport 74 (IFT74), WD repeat-containing protein 90 (WDR90), and Oncomodulin (OCM), which are commonly upregulated in tumors. Probes associated with established tumor suppressor genes (ATM, BRCA1, PALB2, and TP53) were hypermethylated among WTC-exposed breast cancer cases compared to the unexposed group. When comparing WTC EHC breast cancer cases vs. cancer-free controls, there appeared to be global hypomethylation among WTC-exposed breast cancer cases compared to exposed controls. Functional pathway analysis revealed enrichment of several gene pathways in WTC-exposed breast cancer cases including endocytosis, proteoglycans in cancer, regulation of actin cytoskeleton, axon guidance, focal adhesion, calcium signaling, cGMP-PKG signaling, mTOR, Hippo, and oxytocin signaling. The results suggest potential epigenetic links between WTC exposure and breast cancer in local community members enrolled in the WTC EHC program.

Ontogeny and Vulnerabilities of Drug-Tolerant Persisters in HER2+ Breast Cancer.

Resistance to targeted therapies is an important clinical problem in HER2-positive (HER2+) breast cancer. "Drug-tolerant persisters" (DTP), a subpopulation of cancer cells that survive via reversible, nongenetic mechanisms, are implicated in resistance to tyrosine kinase inhibitors (TKI) in other malignancies, but DTPs following HER2 TKI exposure have not been well characterized. We found that HER2 TKIs evoke DTPs with a luminal-like or a mesenchymal-like transcriptome. Lentiviral barcoding/single-cell RNA sequencing reveals that HER2+ breast cancer cells cycle stochastically through a "pre-DTP" state, characterized by a G0-like expression signature and enriched for diapause and/or senescence genes. Trajectory analysis/cell sorting shows that pre-DTPs preferentially yield DTPs upon HER2 TKI exposure. Cells with similar transcriptomes are present in HER2+ breast tumors and are associated with poor TKI response. Finally, biochemical experiments indicate that luminal-like DTPs survive via estrogen receptor-dependent induction of SGK3, leading to rewiring of the PI3K/AKT/mTORC1 pathway to enable AKT-independent mTORC1 activation. SIGNIFICANCE: DTPs are implicated in resistance to anticancer therapies, but their ontogeny and vulnerabilities remain unclear. We find that HER2 TKI-DTPs emerge from stochastically arising primed cells ("pre-DTPs") that engage either of two distinct transcriptional programs upon TKI exposure. Our results provide new insights into DTP ontogeny and potential therapeutic vulnerabilities. This article is highlighted in the In This Issue feature, p. 873.

Role of plakophilin-2 expression on exercise-related progression of arrhythmogenic right ventricular cardiomyopathy: a translational study.

AIMS: Exercise increases arrhythmia risk and cardiomyopathy progression in arrhythmogenic right ventricular cardiomyopathy (ARVC) patients, but the mechanisms remain unknown. We investigated transcriptomic changes caused by endurance training in mice deficient in plakophilin-2 (PKP2cKO), a desmosomal protein important for intercalated disc formation, commonly mutated in ARVC and controls. METHODS AND RESULTS: Exercise alone caused transcriptional downregulation of genes coding intercalated disk proteins. The changes converged with those in sedentary and in exercised PKP2cKO mice. PKP2 loss caused cardiac contractile deficit, decreased muscle mass and increased functional/transcriptomic signatures of apoptosis, despite increased fractional shortening and calcium transient amplitude in single myocytes. Exercise accelerated cardiac dysfunction, an effect dampened by pre-training animals prior to PKP2-KO. Consistent with PKP2-dependent muscle mass deficit, cardiac dimensions in human athletes carrying PKP2 mutations were reduced, compared to matched controls. CONCLUSIONS: We speculate that exercise challenges a cardiomyocyte "desmosomal reserve" which, if impaired genetically (e.g., PKP2 loss), accelerates progression of cardiomyopathy.

Antiplatelet Effects of Clopidogrel Vs Aspirin in Virologically Controlled HIV: A Randomized Controlled Trial.

Patients with HIV exhibit platelet activation and increased risk of cardiovascular disease, the prevention of which is not fully known. Fifty-five HIV-positive patients were randomized to clopidogrel, aspirin, or no-treatment for 14 days, and the platelet phenotype and ability to induce endothelial inflammation assessed. Clopidogrel as opposed to aspirin and no-treatment reduced platelet activation (P-selectin and PAC-1 expression). Compared with baseline, platelet-induced proinflammatory transcript expression of cultured endothelial cells were reduced in those assigned to clopidogrel, with no change in the aspirin and no-treatment arms. In HIV, clinical trials of clopidogrel to prevent cardiovascular disease are warranted. (Antiplatelet Therapy in HIV; NCT02559414).

Multimodal single-cell analysis of cutaneous T-cell lymphoma reveals distinct subclonal tissue-dependent signatures.

Cutaneous T-cell lymphoma (CTCL) is a heterogeneous group of mature T-cell neoplasms characterized by the accumulation of clonal malignant CD4+ T cells in the skin. The most common variant of CTCL, mycosis fungoides (MF ), is confined to the skin in early stages but can be accompanied by extracutaneous dissemination of malignant T cells to the blood and lymph nodes in advanced stages of disease. Sézary syndrome (SS), a leukemic form of disease, is characterized by significant blood involvement. Little is known about the transcriptional and genomic relationship between skin- and blood-residing malignant T cells in CTCL. To identify and interrogate malignant clones in matched skin and blood from patients with leukemic MF and SS, we combine T-cell receptor clonotyping with quantification of gene expression and cell surface markers at the single cell level. Our data reveal clonal evolution at a transcriptional and genetic level within the malignant populations of individual patients. We highlight highly consistent transcriptional signatures delineating skin- and blood-derived malignant T cells. Analysis of these 2 populations suggests that environmental cues, along with genetic aberrations, contribute to transcriptional profiles of malignant T cells. Our findings indicate that the skin microenvironment in CTCL promotes a transcriptional response supporting rapid malignant expansion, as opposed to the quiescent state observed in the blood, potentially influencing efficacy of therapies. These results provide insight into tissue-specific characteristics of cancerous cells and underscore the need to address the patients' individual malignant profiles at the time of therapy to eliminate all subclones.

Molecular analysis of encapsulated papillary carcinoma of the breast with and without invasion.

Encapsulated papillary carcinomas (EPCs) of the breast are a unique variant of papillary carcinoma confined to a cystic space with absent or attenuated myoepithelial cell layer. Although staged as an in situ lesion, it can be associated with invasive ductal carcinoma (IDC). We sought to compare the genomic characteristics of pure EPC and EPC with associated invasive carcinoma (EPCi) at the genomic level. All cases of EPCi harbored recurrent hotspot mutations in PIK3CA. PIK3CA, KMT2A, and CREBBP deleterious somatic events were found across both tumor groups, irrespective of invasion status. At the whole transcriptomic level, EPCi cases displayed remarkably similar mRNA profiles when compared to EPC. When EPCi cases were compared with their corresponding IDC, despite significant overlap, we identified differential gene expression in 39 genes with enrichment of multiple pathways including extracellular matrix regulation, cell adhesion, and collagen fibril organization. Despite morphologic, genotypic, and transcriptomic overlap between pure EPC and EPCi, the latter tumors are likely advanced lesions with PIK3CA activating mutations and enrichment of stromal-related genes implicated in the switch to IDC.

Serial single-cell profiling analysis of metastatic TNBC during Nab-paclitaxel and pembrolizumab treatment.

PURPOSE: Immunotherapy has recently been shown to improve outcomes for advanced PD-L1-positive triple-negative breast cancer (TNBC) in the Impassion130 trial, leading to FDA approval of the first immune checkpoint inhibitor in combination with taxane chemotherapy. To further develop predictive biomarkers and improve therapeutic efficacy of the combination, interrogation of the tumor immune microenvironment before therapy as well as during each component of treatment is crucial. Here we use single-cell RNA sequencing (scRNA-seq) on tumor biopsies to assess immune cell changes from two patients with advanced TNBC treated in a prospective trial at predefined serial time points, before treatment, on taxane chemotherapy and on chemo-immunotherapy. METHODS: Both patients (one responder and one progressor) received the trial therapy, in cycle 1 nab-paclitaxel given as single agent, in cycle 2 nab-paclitaxel in combination with pembrolizumab. Tumor core biopsies were obtained at baseline, 3 weeks (after cycle 1, chemotherapy alone) and 6 weeks (after cycle 2, chemo-immunotherapy). Single-cell RNA sequencing (scRNA-seq) of both cancer cells and infiltrating immune cells isolated were performed from fresh tumor core biopsy specimens by 10 × chromium sequencing. RESULTS: ScRNA-seq analysis showed significant baseline heterogeneity of tumor-infiltrating immune cell populations between the two patients as well as modulation of the tumor microenvironment by chemotherapy and immunotherapy. In the responding patient there was a population of PD-1high-expressing T cells which significantly decreased after nab-paclitaxel plus pembrolizumab treatment as well as a presence of tissue-resident memory T cells (TRM). In contrast, tumors from the patient with rapid disease progression showed a prevalent and persistent myeloid compartment. CONCLUSIONS: Our study provides a deep cellular analysis of on-treatment changes during chemo-immunotherapy for advanced TNBC, demonstrating not only feasibility of single-cell analyses on serial tumor biopsies but also the heterogeneity of TNBC and differences in on-treatment changes in responder versus progressor.

Lower Airway Dysbiosis Affects Lung Cancer Progression.

In lung cancer, enrichment of the lower airway microbiota with oral commensals commonly occurs, and ex vivo models support that some of these bacteria can trigger host transcriptomic signatures associated with carcinogenesis. Here, we show that this lower airway dysbiotic signature was more prevalent in the stage IIIB-IV tumor-node-metastasis lung cancer group and is associated with poor prognosis, as shown by decreased survival among subjects with early-stage disease (I-IIIA) and worse tumor progression as measured by RECIST scores among subjects with stage IIIB-IV disease. In addition, this lower airway microbiota signature was associated with upregulation of the IL17, PI3K, MAPK, and ERK pathways in airway transcriptome, and we identified Veillonella parvula as the most abundant taxon driving this association. In a KP lung cancer model, lower airway dysbiosis with V. parvula led to decreased survival, increased tumor burden, IL17 inflammatory phenotype, and activation of checkpoint inhibitor markers. SIGNIFICANCE: Multiple lines of investigation have shown that the gut microbiota affects host immune response to immunotherapy in cancer. Here, we support that the local airway microbiota modulates the host immune tone in lung cancer, affecting tumor progression and prognosis.See related commentary by Zitvogel and Kroemer, p. 224.This article is highlighted in the In This Issue feature, p. 211.

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.

Genome-Wide DNA Methylation Profiles in Community Members Exposed to the World Trade Center Disaster.

The primary goal of this pilot study was to assess feasibility of studies among local community members to address the hypothesis that complex exposures to the World Trade Center (WTC) dust and fumes resulted in long-term epigenetic changes. We enrolled 18 WTC-exposed cancer-free women from the WTC Environmental Health Center (WTC EHC) who agreed to donate blood samples during their standard clinical visits. As a reference WTC unexposed group, we randomly selected 24 age-matched cancer-free women from an existing prospective cohort who donated blood samples before 11 September 2001. The global DNA methylation analyses were performed using Illumina Infinium MethylationEpic arrays. Statistical analyses were performed using R Bioconductor package. Functional genomic analyses were done by mapping the top 5000 differentially expressed CpG sites to the Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway database. Among cancer-free subjects, we observed substantial methylation differences between WTC-exposed and unexposed women. The top 15 differentially methylated gene probes included BCAS2, OSGIN1, BMI1, EEF1A2, SPTBN5, CHD8, CDCA7L, AIDA, DDN, SNORD45C, ZFAND6, ARHGEF7, UBXN8, USF1, and USP12. Several cancer-related pathways were enriched in the WTC-exposed subjects, including endocytosis, mitogen-activated protein kinase (MAPK), viral carcinogenesis, as well as Ras-associated protein-1 (Rap1) and mammalian target of rapamycin (mTOR) signaling. The study provides preliminary data on substantial differences in DNA methylation between WTC-exposed and unexposed populations that require validation in further studies.