The circular RNome of primary breast cancer.

Circular RNAs (circRNAs) are a class of RNAs that is under increasing scrutiny, although their functional roles are debated. We analyzed RNA-seq data of 348 primary breast cancers and developed a method to identify circRNAs that does not rely on unmapped reads or known splice junctions. We identified 95,843 circRNAs, of which 20,441 were found recurrently. Of the circRNAs that match exon boundaries of the same gene, 668 showed a poor or even negative (R < 0.2) correlation with the expression level of the linear gene. In silico analysis showed only a minority (8.5%) of circRNAs could be explained by known splicing events. Both these observations suggest that specific regulatory processes for circRNAs exist. We confirmed the presence of circRNAs of CNOT2, CREBBP, and RERE in an independent pool of primary breast cancers. We identified circRNA profiles associated with subgroups of breast cancers and with biological and clinical features, such as amount of tumor lymphocytic infiltrate and proliferation index. siRNA-mediated knockdown of circCNOT2 was shown to significantly reduce viability of the breast cancer cell lines MCF-7 and BT-474, further underlining the biological relevance of circRNAs. Furthermore, we found that circular, and not linear, CNOT2 levels are predictive for progression-free survival time to aromatase inhibitor (AI) therapy in advanced breast cancer patients, and found that circCNOT2 is detectable in cell-free RNA from plasma. We showed that circRNAs are abundantly present, show characteristics of being specifically regulated, are associated with clinical and biological properties, and thus are relevant in breast cancer.

RUNX1-targeted therapy for AML expressing somatic or germline mutation in RUNX1.

RUNX1 transcription factor regulates normal and malignant hematopoiesis. Somatic or germline mutant RUNX1 (mtRUNX1) is associated with poorer outcome in acute myeloid leukemia (AML). Knockdown or inhibition of RUNX1 induced more apoptosis of AML expressing mtRUNX1 versus wild-type RUNX1 and improved survival of mice engrafted with mtRUNX1-expressing AML. CRISPR/Cas9-mediated editing-out of RUNX1 enhancer (eR1) within its intragenic super-enhancer, or BET protein BRD4 depletion by short hairpin RNA, repressed RUNX1, inhibited cell growth, and induced cell lethality in AML cells expressing mtRUNX1. Moreover, treatment with BET protein inhibitor or degrader (BET-proteolysis targeting chimera) repressed RUNX1 and its targets, inducing apoptosis and improving survival of mice engrafted with AML expressing mtRUNX1. Library of Integrated Network-based Cellular Signatures 1000-connectivity mapping data sets queried with messenger RNA signature of RUNX1 knockdown identified novel expression-mimickers (EMs), which repressed RUNX1 and exerted in vitro and in vivo efficacy against AML cells expressing mtRUNX1. In addition, the EMs cinobufagin, anisomycin, and narciclasine induced more lethality in hematopoietic progenitor cells (HPCs) expressing germline mtRUNX1 from patients with AML compared with HPCs from patients with familial platelet disorder (FPD), or normal untransformed HPCs. These findings highlight novel therapeutic agents for AML expressing somatic or germline mtRUNX1.

Multiregional transcriptomic profiling provides improved prognostic insight in localized non-small cell lung cancer.

Lung Cancer remains the leading cause of cancer deaths in the USA and worldwide. Non-small cell lung cancer (NSCLC) harbors high transcriptomic intratumor heterogeneity (RNA-ITH) that limits the reproducibility of expression-based prognostic models. In this study, we used multiregional RNA-seq data (880 tumor samples from 350 individuals) from both public (TRACERx) and internal (MDAMPLC) cohorts to investigate the effect of RNA-ITH on prognosis in localized NSCLC at the gene, signature, and tumor microenvironment levels. At the gene level, the maximal expression of hazardous genes (expression negatively associated with survival) but the minimal expression of protective genes (expression positively associated with survival) across different regions within a tumor were more prognostic than the average expression. Following that, we examined whether multiregional expression profiling can improve the performance of prognostic signatures. We investigated 11 gene signatures collected from previous publications and one signature developed in this study. For all of them, the prognostic prediction accuracy can be significantly improved by converting the regional expression of signature genes into sample-specific expression with a simple function-taking the maximal expression of hazardous genes and the minimal expression of protective genes. In the tumor microenvironment, we found a similar rule also seems applicable to immune ITH. We calculated the infiltration levels of major immune cell types in each region of a sample based on expression deconvolution. Prognostic analysis indicated that the region with the lowest infiltration level of protective or highest infiltration level of hazardous immune cells determined the prognosis of NSCLC patients. Our study highlighted the impact of RNA-ITH on the prognostication of NSCLC, which should be taken into consideration to optimize the design and application of expression-based prognostic biomarkers and models. Multiregional assays have the great potential to significantly improve their applications to prognostic stratification.

Intratumoral microbiome of adenoid cystic carcinomas and comparison with other head and neck cancers.

Adenoid cystic carcinoma (ACC) is a rare, usually slow-growing yet aggressive head and neck malignancy. Despite its clinical significance, our understanding of the cellular evolution and microenvironment in ACC remains limited. We investigated the intratumoral microbiomes of 50 ACC tumor tissues and 33 adjacent normal tissues using 16S rRNA gene sequencing. This allowed us to characterize the bacterial communities within the ACC and explore potential associations between the bacterial community structure, patient clinical characteristics, and tumor molecular features obtained through RNA sequencing. The bacterial composition in the ACC was significantly different from that in adjacent normal salivary tissue, and the ACC exhibited diverse levels of species richness. We identified two main microbial subtypes within the ACC: oral-like and gut-like. Oral-like microbiomes, characterized by increased diversity and abundance of Neisseria, Leptotrichia, Actinomyces, Streptococcus, Rothia, and Veillonella (commonly found in healthy oral cavities), were associated with a less aggressive ACC-II molecular subtype and improved patient outcomes. Notably, we identified the same oral genera in oral cancer and head and neck squamous cell carcinomas. In both cancers, they were part of shared oral communities associated with a more diverse microbiome, less aggressive tumor phenotype, and better survival that reveal the genera as potential pancancer biomarkers for favorable microbiomes in ACC and other head and neck cancers. Conversely, gut-like intratumoral microbiomes, which feature low diversity and colonization by gut mucus layer-degrading species, such as Bacteroides, Akkermansia, Blautia, Bifidobacterium, and Enterococcus, were associated with poorer outcomes. Elevated levels of Bacteroides thetaiotaomicron were independently associated with significantly worse survival and positively correlated with tumor cell biosynthesis of glycan-based cell membrane components.

The evolution of lung adenocarcinoma precursors is associated with chromosomal instability and transition from innate to adaptive immune response/evasion.

Studying lung adenocarcinoma (LUAD) early carcinogenesis is challenging, primarily due to the lack of LUAD precursors specimens. We amassed multi-omics data from 213 LUAD and LUAD precursors to identify molecular features underlying LUAD precancer evolution. We observed progressively increasing mutations, chromosomal aberrations, whole genome doubling and genomic instability from precancer to invasive LUAD, indicating aggravating chromosomal instability (CIN). Telomere shortening, a crucial genomic alteration linked to CIN, emerged at precancer stage. Moreover, later-stage lesions demonstrated increasing cancer stemness and decreasing alveolar identity, suggesting epithelial de-differentiation during early LUAD carcinogenesis. The innate immune cells progressively diminished from precancer to invasive LUAD, concomitant with a gradual recruitment of adaptive immune cells (except CD8+ and gamma-delta T cells that decreased in later stages) and upregulation of numerous immune checkpoints, suggesting LUAD precancer evolution is associated with a shift from innate to adaptive immune response and immune evasion mediated by various mechanisms.

The Single-Cell Immunogenomic Landscape of B and Plasma Cells in Early-Stage Lung Adenocarcinoma.

Tumor-infiltrating B and plasma cells (TIB) are prevalent in lung adenocarcinoma (LUAD); however, they are poorly characterized. We performed paired single-cell RNA and B-cell receptor (BCR) sequencing of 16 early-stage LUADs and 47 matching multiregion normal tissues. By integrative analysis of ∼50,000 TIBs, we define 12 TIB subsets in the LUAD and adjacent normal ecosystems and demonstrate extensive remodeling of TIBs in LUADs. Memory B cells and plasma cells (PC) were highly enriched in tumor tissues with more differentiated states and increased frequencies of somatic hypermutation. Smokers exhibited markedly elevated PCs and PCs with distinct differentiation trajectories. BCR clonotype diversity increased but clonality decreased in LUADs, smokers, and with increasing pathologic stage. TIBs were mostly localized within CXCL13+ lymphoid aggregates, and immune cell sources of CXCL13 production evolved with LUAD progression and included elevated fractions of CD4 regulatory T cells. This study provides a spatial landscape of TIBs in early-stage LUAD. SIGNIFICANCE: While TIBs are highly enriched in LUADs, they are poorly characterized. This study provides a much-needed understanding of the transcriptional, clonotypic states and phenotypes of TIBs, unraveling their potential roles in the immunopathology of early-stage LUADs and constituting a road map for the development of TIB-targeted immunotherapies for the treatment of this morbid malignancy. This article is highlighted in the In This Issue feature, p. 2483.

Genomic and Single-Cell Landscape Reveals Novel Drivers and Therapeutic Vulnerabilities of Transformed Cutaneous T-cell Lymphoma.

ABSTRACT: Cutaneous T-cell lymphoma (CTCL) is a rare cancer of skin-homing T cells. A subgroup of patients develops large cell transformation with rapid progression to an aggressive lymphoma. Here, we investigated the transformed CTCL (tCTCL) tumor ecosystem using integrative multiomics spanning whole-exome sequencing (WES), single-cell RNA sequencing, and immune profiling in a unique cohort of 56 patients. WES of 70 skin biopsies showed high tumor mutation burden, UV signatures that are prognostic for survival, exome-based driver events, and most recurrently mutated pathways in tCTCL. Single-cell profiling of 16 tCTCL skin biopsies identified a core oncogenic program with metabolic reprogramming toward oxidative phosphorylation (OXPHOS), cellular plasticity, upregulation of MYC and E2F activities, and downregulation of MHC I suggestive of immune escape. Pharmacologic perturbation using OXPHOS and MYC inhibitors demonstrated potent antitumor activities, whereas immune profiling provided in situ evidence of intercellular communications between malignant T cells expressing macrophage migration inhibitory factor and macrophages and B cells expressing CD74. SIGNIFICANCE: Our study contributes a key resource to the community with the largest collection of tCTCL biopsies that are difficult to obtain. The multiomics data herein provide the first comprehensive compendium of genomic alterations in tCTCL and identify potential prognostic signatures and novel therapeutic targets for an incurable T-cell lymphoma. This article is highlighted in the In This Issue feature, p. 1171.

Expansion of Candidate HPV-Specific T Cells in the Tumor Microenvironment during Chemoradiotherapy Is Prognostic in HPV16+ Cancers.

Human papillomavirus (HPV) infection causes 600,000 new cancers worldwide each year. HPV-related cancers express the oncogenic proteins E6 and E7, which could serve as tumor-specific antigens. It is not known whether immunity to E6 and E7 evolves during chemoradiotherapy or affects survival. Using T cells from 2 HPV16+ patients, we conducted functional T-cell assays to identify candidate HPV-specific T cells and common T-cell receptor motifs, which we then analyzed across 86 patients with HPV-related cancers. The HPV-specific clones and E7-related T-cell receptor motifs expanded in the tumor microenvironment over the course of treatment, whereas non-HPV-specific T cells did not. In HPV16+ patients, improved recurrence-free survival was associated with HPV-responsive T-cell expansion during chemoradiotherapy.

Clonal Hematopoiesis Is Associated with Increased Risk of Severe Neurotoxicity in Axicabtagene Ciloleucel Therapy of Large B-Cell Lymphoma.

To explore the role of clonal hematopoiesis (CH) in chimeric antigen receptor (CAR) T-cell therapy outcomes, we performed targeted deep sequencing on buffy coats collected during the 21 days before lymphodepleting chemotherapy from 114 large B-cell lymphoma patients treated with anti-CD19 CAR T cells. We detected CH in 42 (36.8%) pretreatment samples, most frequently in PPM1D (19/114) and TP53 (13/114) genes. Grade ≥3 immune effector cell-associated neurotoxicity syndrome (ICANS) incidence was higher in CH-positive patients than CH-negative patients (45.2% vs. 25.0%, P = 0.038). Higher toxicities with CH were primarily associated with DNMT3A, TET2, and ASXL1 genes (DTA mutations). Grade ≥3 ICANS (58.9% vs. 25%, P = 0.02) and ≥3 cytokine release syndrome (17.7% vs. 4.2%, P = 0.08) incidences were higher in DTA-positive than in CH-negative patients. The estimated 24-month cumulative incidence of therapy-related myeloid neoplasms after CAR T-cell therapy was higher in CH-positive than CH-negative patients [19% (95% CI, 5.5-38.7) vs. 4.2% (95% CI, 0.3-18.4), P = 0.028]. SIGNIFICANCE: Our study reveals that CH mutations, especially those associated with inflammation (DNMT3A, TET2, and ASXL1), are associated with severe-grade neurotoxicities in lymphoma patients receiving anti-CD19 CAR T-cell therapy. Further studies to investigate the mechanisms and interventions to improve toxicities in the context of CH are warranted. See related content by Uslu and June, p. 382. This article is highlighted in the In This Issue feature, p. 369.

Whole-Genome Sequencing of Common Salivary Gland Carcinomas: Subtype-Restricted and Shared Genetic Alterations.

PURPOSE: Salivary gland carcinomas (SGCs) are pathologically classified into several widely diverse subtypes, of which adenoid cystic carcinoma (ACC), mucoepidermoid carcinoma (MEC), and salivary duct carcinoma (SDC) are the most commonly encountered. A comparative genetic analysis of these subtypes provides detailed information on the genetic alterations that are associated with their tumorigenesis and may lead to the identification of biomarkers to guide tumor-specific clinical trials. EXPERIMENTAL DESIGN: Whole-genome sequencing of 58 common SGCs (20 ACCs, 20 SDCs, and 18 MECs) was performed to catalog structural variations, copy number, rearrangements, and driver mutations. Data were bioinformatically analyzed and correlated with clinicopathologic parameters, and selected targets were validated. RESULTS: Novel and recurrent type-specific and shared genetic alterations were identified within and among 3 subtypes. Mutually exclusive canonical fusion and nonfusion genomic alterations were identified in both ACC and MEC. In ACCs, loss of chromosome 12q was dominant in MYB or MYBL1 fusion-positive tumors and mutations of NOTCH pathway were more common in these fusion negatives. In MECs, CRTC1-MAML2 fusion-positive tumors showed frequent BAP1 mutation, and tumors lacking this fusion were enriched with LRFN1 mutation. SDCs displayed considerable genetic instability, lacked recurrent chromosomal rearrangements, and demonstrated nonoverlapping TP53 mutation and ERBB2 amplification in a subset of tumors. Limited genetic alterations, including focal amplifications of 8q21-q23, were shared by all subtypes and were associated with poor survival. CONCLUSIONS: This study delineates type-specific and shared genetic alterations that are associated with early phenotypic commitment and the biologic progression of common SGCs. These alterations, upon validation, could serve as biomarkers in tumor-specific clinical trials.

Efficacy, Safety, and Biomarker Analysis of Combined PD-L1 (Atezolizumab) and VEGF (Bevacizumab) Blockade in Advanced Malignant Peritoneal Mesothelioma.

Malignant peritoneal mesothelioma (MPeM) is a rare but aggressive malignancy with limited treatment options. VEGF inhibition enhances efficacy of immune-checkpoint inhibitors by reworking the immunosuppressive tumor milieu. Efficacy and safety of combined PD-L1 (atezolizumab) and VEGF (bevacizumab) blockade (AtezoBev) was assessed in 20 patients with advanced and unresectable MPeM with progression or intolerance to prior platinum-pemetrexed chemotherapy. The primary endpoint of confirmed objective response rate per RECISTv1.1 by independent radiology review was 40% [8/20; 95% confidence interval (CI), 19.1-64.0] with median response duration of 12.8 months. Six (75%) responses lasted for >10 months. Progression-free and overall survival at one year were 61% (95% CI, 35-80) and 85% (95% CI, 60-95), respectively. Responses occurred notwithstanding low tumor mutation burden and PD-L1 expression status. Baseline epithelial-mesenchymal transition gene expression correlated with therapeutic resistance/response (r = 0.80; P = 0.0010). AtezoBev showed promising and durable efficacy in patients with advanced MPeM with an acceptable safety profile, and these results address a grave unmet need for this orphan disease. SIGNIFICANCE: Efficacy of atezolizumab and bevacizumab vis-à-vis response rates and survival in advanced peritoneal mesothelioma previously treated with chemotherapy surpassed outcomes expected with conventional therapies. Biomarker analyses uncovered epithelial-mesenchymal transition phenotype as an important resistance mechanism and showcase the value and feasibility of performing translationally driven clinical trials in rare tumors.See related commentary by Aldea et al., p. 2674.This article is highlighted in the In This Issue feature, p. 2659.

In vitro differential sensitivity of melanomas to phenothiazines is based on the presence of codon 600 BRAF mutation.

The panel of 60 human cancer cell lines (the NCI-60) assembled by the National Cancer Institute for anticancer drug discovery is a widely used resource. We previously sequenced 24 cancer genes in those cell lines. Eleven of the genes were found to be mutated in three or more of the lines. Using a pharmacogenomic approach, we analyzed the relationship between drug activity and mutations in those 11 genes (APC, RB1, KRAS, NRAS, BRAF, PIK3CA, PTEN, STK11, MADH4, TP53, and CDKN2A). That analysis identified an association between mutation in BRAF and the antiproliferative potential of phenothiazine compounds. Phenothiazines have been used as antipsychotics and as adjunct antiemetics during cancer chemotherapy and more recently have been reported to have anticancer properties. However, to date, the anticancer mechanism of action of phenothiazines has not been elucidated. To follow up on the initial pharmacologic observations in the NCI-60 screen, we did pharmacologic experiments on 11 of the NCI-60 cell lines and, prospectively, on an additional 24 lines. The studies provide evidence that BRAF mutation (codon 600) in melanoma as opposed to RAS mutation is predictive of an increase in sensitivity to phenothiazines as determined by 3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium, inner salt assay (Wilcoxon P = 0.007). That pattern of increased sensitivity to phenothiazines based on the presence of codon 600 BRAF mutation may be unique to melanomas, as we do not observe it in a panel of colorectal cancers. The findings reported here have potential implications for the use of phenothiazines in the treatment of V600E BRAF mutant melanoma.

Linking Associations of Rare Low-Abundance Species to Their Environments by Association Networks.

Studies of microbial communities by targeted sequencing of rRNA genes lead to recovering numerous rare low-abundance taxa with unknown biological roles. We propose to study associations of such rare organisms with their environments by a computational framework based on transformation of the data into qualitative variables. Namely, we analyze the sparse table of putative species or OTUs (operational taxonomic units) and samples generated in such studies, also known as an OTU table, by collecting statistics on co-occurrences of the species and on shared species richness across samples. Based on the statistics we built two association networks, of the rare putative species and of the samples respectively, using a known computational technique, Association networks (Anets) developed for analysis of qualitative data. Clusters of samples and clusters of OTUs are then integrated and combined with metadata of the study to produce a map of associated putative species in their environments. We tested and validated the framework on two types of microbiomes, of human body sites and that of the Populus tree root systems. We show that in both studies the associations of OTUs can separate samples according to environmental or physiological characteristics of the studied systems.

Sequencing analysis of BRAF mutations in human cancers.

Cancers arise because of the accumulation of mutations in critical genes that alter normal programs of cell proliferation, differentiation, and death. The RAS-RAF-MEK-ERK-MAP kinase pathway mediates cellular responses to growth signals. RAS is mutated to an oncogenic form in approximately 15% of human cancer. The three RAF genes code for cytoplasmic serine/threonine kinases that are regulated by binding RAS. ARAF and c-RAF are infrequently mutated in human cancer. However, BRAF is mutated in a wide range of human cancers. Most mutations are within the kinase domain, with a single amino acid substitution (V600E) accounting for most mutations.

CASK mutations are frequent in males and cause X-linked nystagmus and variable XLMR phenotypes.

Mutations of the calcium/calmodulin-dependent serine protein kinase (CASK) gene have recently been associated with X-linked mental retardation (XLMR) with microcephaly, optic atrophy and brainstem and cerebellar hypoplasia, as well as with an X-linked syndrome having some FG-like features. Our group has recently identified four male probands from 358 probable XLMR families with missense mutations (p.Y268H, p.P396S, p.D710G and p.W919R) in the CASK gene. Congenital nystagmus, a rare and striking feature, was present in two of these families. We screened a further 45 probands with either nystagmus or microcephaly and mental retardation (MR), and identified two further mutations, a missense mutation (p.Y728C) and a splice mutation (c.2521-2A>T) in two small families with nystagmus and MR. Detailed clinical examinations of all six families, including an ophthalmological review in four families, were undertaken to further characterise the phenotype. We report on the clinical features of 24 individuals, mostly male, from six families with CASK mutations. The phenotype was variable, ranging from non-syndromic mild MR to severe MR associated with microcephaly and dysmorphic facial features. Carrier females were variably affected. Congenital nystagmus was found in members of four of the families. Our findings reinforce the CASK gene as a relatively frequent cause of XLMR in females and males. We further define the phenotypic spectrum and demonstrate that affected males with missense mutations or in-frame deletions in CASK are frequently associated with congenital nystagmus and XLMR, a striking feature not previously reported.

Evaluation of linkage of breast cancer to the putative BRCA3 locus on chromosome 13q21 in 128 multiple case families from the Breast Cancer Linkage Consortium.

The known susceptibility genes for breast cancer, including BRCA1 and BRCA2, only account for a minority of the familial aggregation of the disease. A recent study of 77 multiple case breast cancer families from Scandinavia found evidence of linkage between the disease and polymorphic markers on chromosome 13q21. We have evaluated the contribution of this candidate "BRCA3" locus to breast cancer susceptibility in 128 high-risk breast cancer families of Western European ancestry with no identified BRCA1 or BRCA2 mutations. No evidence of linkage was found. The estimated proportion (alpha) of families linked to a susceptibility locus at D13S1308, the location estimated by Kainu et al. [(2000) Proc. Natl. Acad. Sci. USA 97, 9603-9608], was 0 (upper 95% confidence limit 0.13). Adjustment for possible bias due to selection of families on the basis of linkage evidence at BRCA2 did not materially alter this result (alpha = 0, upper 95% confidence limit 0.18). The proportion of linked families reported by Kainu et al. (0.65) is excluded with a high degree of confidence in our dataset [heterogeneity logarithm of odds (HLOD) at alpha = 0.65 was -11.0]. We conclude that, if a susceptibility gene does exist at this locus, it can only account for a small proportion of non-BRCA1/2 families with multiple cases of early-onset breast cancer.