Search

Abstract Type II innate lymphoid cells (ILC2s) maintain homeostasis and barrier integrity in mucosal tissues. In both mice and humans, ILC2s poorly reconstitute after allogeneic hematopoietic stem cell transplantation (allo-HSCT). Determining the mechanisms involved in their impaired reconstitution could improve transplant outcomes. By integrating single-cell chromatin and transcriptomic analyses of transplanted ILC2s, we identify a previously unreported population of converted ILC1-like cells in the mouse small intestine post-transplant. Exposure of ILC2s to proinflammatory cytokines resulted in a mixed ILC1-ILC2 phenotype but was able to convert only a small population of ILC2s to ILC1s, which were found post-transplant. Whereas ILC2s protected against acute graft-versus-host disease (aGVHD) mediated mortality, infusion of proinflammatory cytokine-exposed ILC2s accelerated aGvHD. Interestingly, murine ILC2 reconstitution post-HSCT is decreased in the presence of alloreactive T cells. Finally, peripheral blood cells from human patients with aGvHD have an altered ILC2-associated chromatin landscape compared to transplanted controls. These data demonstrate that following transplantation ILC2s convert to a pro-pathogenic population with an ILC1-like chromatin state and provide insights into the contribution of ILC plasticity to the impaired reconstitution of ILC2 cells, which is one of several potential mechanisms for the poor reconstitution of these important cells after allo-HSCT.

Summary: Background: Early-stage triple-negative breast cancer (TNBC) displays clinical and biological diversity. From a biological standpoint, immune infiltration plays a crucial role in TNBC prognosis. Currently, there is a lack of genomic tools aiding in treatment decisions for TNBC. This study aims to assess the effectiveness of a B-cell/immunoglobulin signature (IGG) alone, or in combination with tumor burden, in predicting prognosis and treatment response in patients with TNBC. Methods: Genomic and clinical data were retrieved from 7 cohorts: SCAN-B (N = 874), BrighTNess (n = 482), CALGB-40603 (n = 389), METABRIC (n = 267), TCGA (n = 118), GSE58812 (n = 107), GSE21653 (n = 67). IGG and a risk score integrating IGG with tumor/nodal staging (IGG-Clin) were assessed for event-free survival (EFS) and overall survival (OS) in each cohort. Random effects model was used to derive pooled effect sizes. Association of IGG with pathological complete response (pCR) was assessed in CALGB-40603 and BrighTNess. Immune significance of IGG was estimated through CIBERSORTx and EcoTyper. Findings: IGG was associated with improved EFS (pooled HR = 0.77, [95% CI = 0.70–0.85], I2 = 18%) and OS (pooled HR = 0.79, [0.73–0.85], I2 = 0%) across cohorts, and was predictive of pCR in CALGB-40603 (OR 1.25, [1.10–1.50]) and BrighTNess (OR 1.57 [1.25–1.98]). IGG-Clin was predictive of recurrence (pooled HR = 2.11, [1.75–2.55], I2 = 0%) and death (pooled HR = 1.99, 95% [0.84–4.73], I2 = 79%) across cohorts. IGG was associated with adaptive immune response at CIBERSORTx and EcoTyper analysis. Interpretation: IGG is linked to improved prognosis and pCR in early-stage TNBC. The integration of IGG alongside tumor and nodal staging holds promise as an approach to identify patients benefitting from intensified or de-intensified treatments. Funding: This study received funding from: Associació Beca Marta Santamaria, European Union’s Horizon 2020 research and innovation and Marie Skłodowska–Curie Actions programs, Fundación FERO, Fundación CRIS contra el cáncer, Agència de Gestó d'Ajuts Universitaris i de Recerca, Instituto de Salud Carlos III, Fundación Contigo, Asociación Cáncer de Mama Metastásico IV, Breast Cancer Research Foundation, RESCUER, Fundación científica AECC and FSEOM.

Abstract The identification of prognostic markers in patients receiving neoadjuvant therapy is crucial for treatment optimization in HER2-positive breast cancer, with the immune microenvironment being a key factor. Here, we investigate the complexity of B and T cell receptor (BCR and TCR) repertoires in the context of two phase III trials, NeoALTTO and CALGB 40601, evaluating neoadjuvant paclitaxel with trastuzumab and/or lapatinib in women with HER2-positive breast cancer. BCR features, particularly the number of reads and clones, evenness and Gini index, are heterogeneous according to hormone receptor status and PAM50 subtypes. Moreover, BCR measures describing clonal expansion, namely evenness and Gini index, are independent prognostic factors. We present a model developed in NeoALTTO and validated in CALGB 40601 that can predict event-free survival (EFS) by integrating hormone receptor and clinical nodal status, breast pathological complete response (pCR), stromal tumor-infiltrating lymphocyte levels (%) and BCR repertoire evenness. A prognostic score derived from the model and including those variables, HER2-EveNT, allows the identification of patients with 5-year EFS > 90%, and, in those not achieving pCR, of a subgroup of immune-enriched tumors with an excellent outcome despite residual disease.

Plasma ctDNA is a promising method to determine patient outcome in multiple cancer types. Here, the authors use shallow WGS to create machine learning signatures to identify tumor phenotypes and predict therapy response in patients with metastatic breast cancer.

Obesity promotes triple-negative breast cancer (TNBC), and effective interventions are urgently needed to break the obesity-TNBC link. Epidemiologic studies indicate that bariatric surgery reduces TNBC risk, while evidence is limited or conflicted for weight loss via low-fat diet (LFD) or calorie restriction (CR). Using a murine model of obesity-driven TNBC, we compared the antitumor effects of vertical sleeve gastrectomy (VSG) with LFD, chronic CR, and intermittent CR. Each intervention generated weight and fat loss and suppressed tumor growth relative to obese mice (greatest suppression with CR). VSG and CR regimens exerted both similar and unique effects, as assessed using multiomics approaches, in reversing obesity-associated transcript, epigenetics, secretome, and microbiota changes and restoring antitumor immunity. Thus, in a murine model of TNBC, bariatric surgery and CR each reverse obesity-driven tumor growth via shared and distinct antitumor mechanisms, and CR is superior to VSG in reversing obesity’s procancer effects.

Antiretroviral therapy (ART) is not curative due to the existence of cellular reservoirs of latent HIV-1 that persist during therapy. Current research efforts to cure HIV-1 infection include “shock and kill” strategies to disrupt latency using small molecules or latency-reversing agents (LRAs) to induce expression of HIV-1 enabling cytotoxic immune cells to eliminate infected cells. The modest success of current LRAs urges the field to identify novel drugs with increased clinical efficacy. Aminobisphosphonates (N-BPs) that include pamidronate, zoledronate, or alendronate, are the first-line treatment of bone-related diseases including osteoporosis and bone malignancies. Here, we show the use of N-BPs as a novel class of LRA: we found in ex vivo assays using primary cells from ART-suppressed people living with HIV-1 that N-BPs induce HIV-1 from latency to levels that are comparable to the T cell activator phytohemagglutinin (PHA). RNA sequencing and mechanistic data suggested that reactivation may occur through activation of the activator protein 1 signaling pathway. Stored samples from a prior clinical trial aimed at analyzing the effect of alendronate on bone mineral density, provided further evidence of alendronate-mediated latency reversal and activation of immune effector cells. Decay of the reservoir measured by IPDA was however not detected. Our results demonstrate the novel use of N-BPs to reverse HIV-1 latency while inducing immune effector functions. This preliminary evidence merits further investigation in a controlled clinical setting possibly in combination with therapeutic vaccination.

A joint dimension reduction approach integrates gene expression data from human tumor cohorts, human cancer cell lines, and mouse model tumors to improve prediction of clinical biomarkers from model systems.

Mechanical stretching of skeletal muscle cells induces extensive transcriptional and posttranscriptional changes in genes encoding proteins involved in transcription and muscle cell differentiation. SR proteins also exhibit altered transcription and phosphorylation in response to stretch.

Super-enhancers and their associated transcription factor networks have been shown to influence ovarian cancer biology. Here, based on an integrated set of genomic and epigenomic datasets, the authors identify clinically relevant super-enhancers amplified in ovarian cancer patients and functionally validate their activity.

The clinical approvals of KRAS G12C inhibitors have been a revolutionary advance in precision oncology, but response rates are often modest. To improve patient selection, we developed an integrated model to predict KRAS dependency. By integrating molecular profiles of a large panel of cell lines from the DEMETER2 dataset, we built a binary classifier to predict a tumor's KRAS dependency. Monte Carlo cross validation via ElasticNet within the training set was used to compare model performance and to tune parameters α and λ. The final model was then applied to the validation set. We validated the model with genetic depletion assays and an external dataset of lung cancer cells treated with a G12C inhibitor. We then applied the model to several Cancer Genome Atlas (TCGA) datasets. The final "K20" model contains 20 features, including expression of 19 genes and KRAS mutation status. In the validation cohort, K20 had an AUC of 0.94 and accurately predicted KRAS dependency in both mutant and KRAS wild-type cell lines following genetic depletion. It was also highly predictive across an external dataset of lung cancer lines treated with KRAS G12C inhibition. When applied to TCGA datasets, specific subpopulations such as the invasive subtype in colorectal cancer and copy number high pancreatic adenocarcinoma were predicted to have higher KRAS dependency. The K20 model has simple yet robust predictive capabilities that may provide a useful tool to select patients with KRAS mutant tumors that are most likely to respond to direct KRAS inhibitors.

Therapeutic options for non-small cell lung cancer patients with brain metastases are limited. Here the authors design B7-H3 targeting CAR-T cells engineered to express the chemokine receptor CCR2b, and show improved accumulation in the brain and enhanced anti-tumor activity in preclinical models of lung cancer brain metastases.

Summary: Background: HER2DX is a prognostic and predictive assay in early-stage HER2-positive breast cancer based on clinical features and the expression of 4 gene signatures (immune, proliferation, luminal differentiation and HER2 amplicon), including ERBB2 mRNA levels. Here, we evaluated the ability of HER2DX to predict efficacy of a de-escalated, chemotherapy-free neoadjuvant regimen in HER2-positive/hormone receptor-positive breast cancer. Methods: HER2DX was evaluated on pre-treatment tumour samples from the PerELISA phase II study focused on postmenopausal patients with operable HER2-positive/hormone receptor-positive breast cancer. Patients received 2-weeks of letrozole, and then underwent a re-biopsy for Ki67 evaluation. Patients with endocrine therapy sensitive disease (ESD) (i.e., >20.0% Ki67 relative reduction at week 2) continued letrozole and 5 cycles of trastuzumab and pertuzumab. Primary aim was to test the ability of HER2DX risk-score, HER2DX pCR score and HER2DX ERBB2 mRNA score (as continuous variables and group categories) to predict pathological complete response (pCR) in patients with ESD. Logistic regression and receiver...operator curve (ROC) analysis assessed associations of HER2DX scores with pCR and ESD. Findings: HER2DX was evaluated in 55 patients (86.0%) enrolled in PerELISA and 40 patients (73.0%) had ESD. The pCR rate in patients with ESD was 22.5% (9/40). In this group, HER2DX pCR score and HER2DX ERBB2 mRNA score were significantly associated with pCR (p.ß=.ß0.008 and p.ß=.ß0.003, univariate logistic regression model; area under ROC [AUC].ß=.ß0.803 and 0.896). The pCR rate in low, medium, and high HER2DX pCR score groups was 7.7% (2/26), 46.2% (6/13) and 100.0% (1/1), respectively. The pCR rate in low, medium, and high HER2DX ERBB2 score groups was 0.0% (0/12), 7.7% (1/13) and 53.3% (8/15), respectively. HER2DX pCR score was also significantly associated with Ki-67 response following 2-weeks of letrozole (p.ß=.ß0.002, univariate logistic regression model; AUC.ß=.ß0.775). The rate of ESD in low, medium, and high HER2DX pCR score groups was 89.7% (26/29), 65.0% (13/20) and 16.7% (1/6), respectively. Interpretation: HER2DX predicts response following neoadjuvant letrozole in combination with dual HER2 blockade with trastuzumab and pertuzumab in early-stage HER2-positive/hormone receptor-positive breast cancer. Funding: This study received funding from Reveal Genomics.

Abstract Background Preterm birth is the leading cause of neonatal morbidity and mortality, but research efforts in neonatology are complicated due to the unavailability of large volume blood samples. Whole blood assays can be used to overcome this problem by performing both functional and gene expression studies using small amounts of blood. Gene expression studies using RT-qPCR estimate mRNA-levels of target genes normalized to reference genes. The goal of this study was to identify and validate stable reference genes applicable to cord blood samples obtained from developing neonates of different gestational age groups as well as to adult peripheral blood samples. Eight reference gene candidates (ACTB, B2M, GAPDH, GUSB, HPRT, PPIB, RPLP0, RPL13) were analyzed using the three published software algorithms Bestkeeper, GeNorm and NormFinder. Results A normalization factor consisting of ACTB and PPIB allows for comparative expression analyses of neonatal samples from different gestational age groups. Normalization factors consisting of GAPDH and PPIB or ACTB and GAPDH are suitable when samples from preterm and full-term neonates and adults are compared. However, all candidate reference genes except RPLP0 exhibited significant intergroup gene expression variance and a higher gene expression towards an older age which resulted in a small but statistically significant systematic bias. Systematic analysis of RNA-seq data revealed new reference gene candidates with potentially superior stability. Conclusions The current study identified suitable normalization factors and proposed the use of the additional single gene RPLP0 to avoid systematic bias. This combination will enable comparative analyses not only between neonates of different gestational ages, but also between neonates and adults, as it facilitates more detailed investigations of developmental gene expression changes. The use of software algorithms did not prevent unintended systematic bias. This generally highlights the need for careful validation of such results to prevent false interpretation of potential age-dependent changes in gene expression. To identify the most stable reference genes in the future, RNA-seq based global approaches are recommended.

BackgroundTriple negative breast cancer (TNBC) is an aggressive variant of breast cancer that lacks the expression of estrogen and progesterone receptors (ER and PR) and HER2. Nearly 50% of patients with advanced TNBC will develop brain metastases (BrM), commonly with progressive extracranial disease. Immunotherapy has shown promise in the treatment of advanced TNBC; however, the immune contexture of BrM remains largely unknown. We conducted a comprehensive analysis of TNBC BrM and matched primary tumors to characterize the genomic and immune landscape of TNBC BrM to inform the development of immunotherapy strategies in this aggressive disease.MethodsWhole-exome sequencing (WES) and RNA sequencing were conducted on formalin-fixed, paraffin-embedded samples of BrM and primary tumors of patients with clinical TNBC (n = 25, n = 9 matched pairs) from the LCCC1419 biobank at UNC—Chapel Hill. Matched blood was analyzed by DNA sequencing as a comparison for tumor WES for the identification of somatic variants. A comprehensive genomics assessment, including mutational and copy number alteration analyses, neoantigen prediction, and transcriptomic analysis of the tumor immune microenvironment were performed.ResultsPrimary and BrM tissues were confirmed as TNBC (23/25 primaries, 16/17 BrM) by immunohistochemistry and of the basal intrinsic subtype (13/15 primaries and 16/19 BrM) by PAM50. Compared to primary tumors, BrM demonstrated a higher tumor mutational burden. TP53 was the most frequently mutated gene and was altered in 50% of the samples. Neoantigen prediction showed elevated cancer testis antigen- and endogenous retrovirus-derived MHC class I-binding peptides in both primary tumors and BrM and predicted that single-nucleotide variant (SNV)-derived peptides were significantly higher in BrM. BrM demonstrated a reduced immune gene signature expression, although a signature associated with fibroblast-associated wound healing was elevated in BrM. Metrics of T and B cell receptor diversity were also reduced in BrM.ConclusionsBrM harbored higher mutational burden and SNV-derived neoantigen expression along with reduced immune gene signature expression relative to primary TNBC. Immune signatures correlated with improved survival, including T cell signatures. Further research will expand these findings to other breast cancer subtypes in the same biobank. Exploration of immunomodulatory approaches including vaccine applications and immune checkpoint inhibition to enhance anti-tumor immunity in TNBC BrM is warranted.

Abstract Inhibition of the HER2/ERBB2 receptor is a keystone to treating HER2-positive malignancies, particularly breast cancer, but a significant fraction of HER2-positive (HER2+) breast cancers recur or fail to respond. Anti-HER2 monoclonal antibodies, like trastuzumab or pertuzumab, and ATP active site inhibitors like lapatinib, commonly lack durability because of adaptive changes in the tumor leading to resistance. HER2+ cell line responses to inhibition with lapatinib were analyzed by RNAseq and ChIPseq to characterize transcriptional and epigenetic changes. Motif analysis of lapatinib-responsive genomic regions implicated the pioneer transcription factor FOXA1 as a mediator of adaptive responses. Lapatinib in combination with FOXA1 depletion led to dysregulation of enhancers, impaired adaptive upregulation of HER3, and decreased proliferation. HER2-directed therapy using clinically relevant drugs (trastuzumab with or without lapatinib or pertuzumab) in a 7-day clinical trial designed to examine early pharmacodynamic response to antibody-based anti-HER2 therapy showed reduced FOXA1 expression was coincident with decreased HER2 and HER3 levels, decreased proliferation gene signatures, and increased immune gene signatures. This highlights the importance of the immune response to anti-HER2 antibodies and suggests that inhibiting FOXA1-mediated adaptive responses in combination with HER2 targeting is a potential therapeutic strategy.

Abstract Background Contamination of reagents and cross contamination across samples is a long-recognized issue in molecular biology laboratories. While often innocuous, contamination can lead to inaccurate results. Cantalupo et al., for example, found HeLa-derived human papillomavirus 18 (H-HPV18) in several of The Cancer Genome Atlas (TCGA) RNA-sequencing samples. This work motivated us to assess a greater number of samples and determine the origin of possible contaminations using viral sequences. To detect viruses with high specificity, we developed the publicly available workflow, VirDetect, that detects virus and laboratory vector sequences in RNA-seq samples. We applied VirDetect to 9143 RNA-seq samples sequenced at one TCGA sequencing center (28/33 cancer types) over 5 years. Results We confirmed that H-HPV18 was present in many samples and determined that viral transcripts from H-HPV18 significantly co-occurred with those from xenotropic mouse leukemia virus-related virus (XMRV). Using laboratory metadata and viral transcription, we determined that the likely contaminant was a pool of cell lines known as the “common reference”, which was sequenced alongside TCGA RNA-seq samples as a control to monitor quality across technology transitions (i.e. microarray to GAII to HiSeq), and to link RNA-seq to previous generation microarrays that standardly used the “common reference”. One of the cell lines in the pool was a laboratory isolate of MCF-7, which we discovered was infected with XMRV; another constituent of the pool was likely HeLa cells. Conclusions Altogether, this indicates a multi-step contamination process. First, MCF-7 was infected with an XMRV. Second, this infected cell line was added to a pool of cell lines, which contained HeLa. Finally, RNA from this pool of cell lines contaminated several TCGA tumor samples most-likely during library construction. Thus, these human tumors with H-HPV or XMRV reads were likely not infected with H-HPV 18 or XMRV.

BACKGROUND The KRAS proto-oncogene is among the most frequently mutated genes in cancer, yet for 40 years it remained an elusive therapeutic target. Recently, allosteric inhibitors that covalently bind to KRAS G12C mutations have been approved for use in lung adenocarcinomas. Although responses are observed, they are often short-lived, thus making in-depth characterization of the mechanisms of resistance of paramount importance.METHODS Here, we present a rapid-autopsy case of a patient who had a KRASG12C-mutant lung adenocarcinoma who initially responded to a KRAS G12C inhibitor but then rapidly developed resistance. Using deep-RNA and whole-exome sequencing comparing pretreatment, posttreatment, and matched normal tissues, we uncover numerous mechanisms of resistance to direct KRAS inhibition.RESULTS In addition to decreased KRAS G12C–mutant allele frequency in refractory tumors, we also found reactivation of the MAPK pathway despite no new mutations in KRAS or its downstream mediators. Tumor cell–intrinsic and non–cell autonomous mechanisms included increased complement activation, coagulation, and tumor angiogenesis, and several lines of evidence of immunologic evasion.CONCLUSION Together, our findings reveal numerous mechanisms of resistance to current KRAS G12C inhibitors through enrichment of clonal populations, KRAS-independent downstream signaling, and diverse remodeling of the tumor microenvironment.FUNDING Richard and Fran Duley, Jimmy and Kay Mann, the NIH, and the North Carolina Biotechnology Center.

Summary: Background: Both clinical and genomic data independently predict survival and treatment response in early-stage HER2-positive breast cancer. Here we present the development and validation of a new HER2DX risk score, and a new HER2DX pathological complete response (pCR) score, both based on a 27-gene expression plus clinical feature-based classifier. Methods: HER2DX is a supervised learning algorithm incorporating tumour size, nodal staging, and 4 gene expression signatures tracking immune infiltration, tumour cell proliferation, luminal differentiation, and the expression of the HER2 amplicon, into a single score. 434 HER2-positive tumours from the Short-HER trial were used to train a prognostic risk model; 268 cases from an independent cohort were used to verify the accuracy of the HER2DX risk score. In addition, 116 cases treated with neoadjuvant anti-HER2-based chemotherapy were used to train a predictive model of pathological complete response (pCR); two independent cohorts of 91 and 67 cases were used to verify the accuracy of the HER2DX pCR likelihood score. Five publicly available independent datasets with >1,000 patients with early-stage HER2-positive disease were also analysed. Findings: In Short-HER, HER2DX variables were associated with good risk outcomes (i.e., immune, and luminal) and poor risk outcomes (i.e., proliferation, and tumour and nodal staging). In an independent cohort, continuous HER2DX risk score was significantly associated with disease-free survival (DFS) (p=0·002); the 5-year DFS in the low-risk group was 97·4% (94·4-100·0%). For the neoadjuvant pCR predictor training cohort, HER2DX variables were associated with pCR (i.e., immune, proliferation and HER2 amplicon) and non-pCR (i.e., luminal, and tumour and nodal staging). In both independent test set cohorts, continuous HER2DX pCR likelihood score was significantly associated with pCR (p

Effective precision medicine strategies rely on the ability to predict tumour behaviour based on molecular characteristics. Here, the authors build models to predict multiple distinct gene expression patterns using DNA copy number alterations

Abstract Background In the USA, the breast cancer mortality rate is 41% higher for African-American women than non-Hispanic White women. While numerous gene expression studies have classified biological features that vary by race and may contribute to poorer outcomes, few studies have experimentally tested these associations. CRYβB2 gene expression has drawn particular interest because of its association with overall survival and African-American ethnicity in multiple cancers. Several reports indicate that overexpression of the CRYβB2 pseudogene, CRYβB2P1, and not CRYβB2 is linked with race and poor outcome. It remains unclear whether either or both genes are linked to breast cancer outcomes. This study investigates CRYβB2 and CRYβB2P1 expression in human breast cancers and breast cancer cell line models, with the goal of elucidating the mechanistic contribution of CRYβB2 and CRYβB2P1 to racial disparities. Methods Custom scripts for CRYβB2 or CRYβB2P1 were generated and used to identify reads that uniquely aligned to either gene. Gene expression according to race and tumor subtype were assessed using all available TCGA breast cancer RNA sequencing alignment samples (n = 1221). In addition, triple-negative breast cancer models engineered to have each gene overexpressed or knocked out were developed and evaluated by in vitro, biochemical, and in vivo assays to identify biological functions. Results We provide evidence that CRYβB2P1 is expressed at higher levels in breast tumors compared to CRYβB2, but only CRYβB2P1 is significantly increased in African-American tumors relative to White American tumors. We show that independent of CRYβB2, CRYβB2P1 enhances tumorigenesis in vivo via promoting cell proliferation. Our data also reveal that CRYβB2P1 may function as a non-coding RNA to regulate CRYβB2 expression. A key observation is that the combined overexpression of both genes was found to suppress cell growth. CRYβB2 overexpression in triple-negative breast cancers increases invasive cellular behaviors, tumor growth, IL6 production, immune cell chemoattraction, and the expression of metastasis-associated genes. These data underscore that both CRYβB2 and CRYβB2P1 promote tumor growth, but their mechanisms for tumor promotion are likely distinct. Conclusions Our findings provide novel data emphasizing the need to distinguish and study the biological effects of both CRYβB2 and CRYβB2P1 as both genes independently promote tumor progression. Our data demonstrate novel molecular mechanisms of two understudied, disparity-linked molecules.

The prognosis of cancer patients with lymph node (LN) metastasis is worse than those without. Here, the authors report that while histone deacetylase 11 (HDAC11) inhibition suppresses tumor growth within the LN, it also promotes cancer cell migration out of the LN to form distant metastasis, and therefore suggest caution with HDAC inhibitors.

Polymerase theta is a widely conserved DNA polymerase that mediates Theta Mediated End Joining. Here authors present a synthetic lethal CRISPR screen to identify DDR gene mutations that induce cellular addiction to Pol theta.

Abstract Background Measures of the adaptive immune response have prognostic and predictive associations in melanoma and other cancer types. Specifically, intratumoral T cell density and function have considerable prognostic and predictive value in skin cutaneous melanoma (SKCM). Less is known about the significance of tumor-infiltrating B cells in SKCM. Our goal was to understand the prognostic and predictive value of B cell phenotypic subsets in SKCM using RNA sequencing. Methods We used our previously published algorithm, V’DJer, to assemble B cell receptor (BCR) repertoires and estimate diversity from short-read RNA sequencing (RNA-seq). We applied machine learning-based cellular phenotype classifiers to measure relative similarity of bulk tumor sample gene expression profiles and different B cell phenotypes. We assessed these aspects of B cell biology in 473 SKCM from the Cancer Genome Atlas Project (TCGA) as well as in RNA-seq data corresponding to tumor samples procured from patients who received CTLA-4 and PD-1 inhibitors for metastatic SKCM. Results We found that the BCR repertoire was associated with different clinical factors, such as tumor tissue site and sex. However, increased clonality of the BCR repertoire was favorably prognostic in SKCM and was prognostic even after first conditioning on various clinical factors. Mutation burden was not correlated with any BCR measurement, and no specific mutation had an altered BCR repertoire. Lack of an assembled BCR in pre-treatment tumor tissues was associated with a lack of anti-tumor response to a CTLA-4 inhibitor in metastatic SKCM. Conclusions These findings suggest an important prognostic and predictive role for B cell characteristics in SKCM. This has implications for melanoma immunobiology and potential development of immunogenomics features to predict survival and response to immunotherapy.

Abstract Background The seed and soil hypothesis was proposed over a century ago to describe why cancer cells (seeds) grow in certain organs (soil). Since then, the genetic properties that define the cancer cells have been heavily investigated; however, genomic mediators within the organ microenvironment that mediate successful metastatic growth are less understood. These studies sought to identify cancer- and organ-specific genomic programs that mediate metastasis. Methods In these studies, a set of 14 human breast cancer patient-derived xenograft (PDX) metastasis models was developed and then tested for metastatic tropism with two approaches: spontaneous metastases from mammary tumors and intravenous injection of PDX cells. The transcriptomes of the cancer cells when growing as tumors or metastases were separated from the transcriptomes of the microenvironment via species-specific separation of the genomes. Drug treatment of PDX spheroids was performed to determine if genes activated in metastases may identify targetable mediators of viability. Results The experimental approaches that generated metastases in PDX models were identified. RNA sequencing of 134 tumors, metastases, and normal non-metastatic organs identified cancer- and organ-specific genomic properties that mediated metastasis. A common genomic response of the liver microenvironment was found to occur in reaction to the invading PDX cells. Genes within the cancer cells were found to be either transiently regulated by the microenvironment or permanently altered due to clonal selection of metastatic sublines. Gene Set Enrichment Analyses identified more than 400 gene signatures that were commonly activated in metastases across basal-like PDXs. A Src signaling signature was found to be extensively upregulated in metastases, and Src inhibitors were found to be cytotoxic to PDX spheroids. Conclusions These studies identified that during the growth of breast cancer metastases, there were genomic changes that occurred within both the cancer cells and the organ microenvironment. We hypothesize that pathways upregulated in metastases are mediators of viability and that simultaneously targeting changes within different cancer cell pathways and/or different tissue compartments may be needed for inhibition of disease progression.

Abstract We propose a statistical boosting method, termed I-Boost, to integrate multiple types of high-dimensional genomics data with clinical data for predicting survival time. I-Boost provides substantially higher prediction accuracy than existing methods. By applying I-Boost to The Cancer Genome Atlas, we show that the integration of multiple genomics platforms with clinical variables improves the prediction of survival time over the use of clinical variables alone; gene expression values are typically more prognostic of survival time than other genomics data types; and gene modules/signatures are at least as prognostic as the collection of individual gene expression data.

Summary: Mucoepidermoid carcinoma (MEC) is a life-threatening salivary gland cancer that is driven primarily by a transcriptional coactivator fusion composed of cyclic AMP-regulated transcriptional coactivator 1 (CRTC1) and mastermind-like 2 (MAML2). The mechanisms by which the chimeric CRTC1/MAML2 (C1/M2) oncoprotein rewires gene expression programs that promote tumorigenesis remain poorly understood. Here, we show that C1/M2 induces transcriptional activation of the non-canonical peroxisome proliferator-activated receptor gamma coactivator-1 alpha (PGC-1α) splice variant PGC-1α4, which regulates peroxisome proliferator-activated receptor gamma (PPARγ)-mediated insulin-like growth factor 1 (IGF-1) expression. This mitogenic transcriptional circuitry is consistent across cell lines and primary tumors. C1/M2-positive tumors exhibit IGF-1 pathway activation, and small-molecule drug screens reveal that tumor cells harboring the fusion gene are selectively sensitive to IGF-1 receptor (IGF-1R) inhibition. Furthermore, this dependence on autocrine regulation of IGF-1 transcription renders MEC cells susceptible to PPARγ inhibition with inverse agonists. These results yield insights into the aberrant coregulatory functions of C1/M2 and identify a specific vulnerability that can be exploited for precision therapy.

Small cell carcinoma of the ovary, hypercalcemic type (SCCOHT) is a rare and aggressive form of ovarian cancer. SCCOHT tumors have inactivating mutations in SMARCA4 (BRG1), one of the two mutually exclusive ATPases of the SWI/SNF chromatin remodeling complex. To address the role that BRG1 loss plays in SCCOHT tumorigenesis, we performed integrative multi-omic analyses in SCCOHT cell lines +/- BRG1 reexpression. BRG1 reexpression induced a gene and protein signature similar to an epithelial cell and gained chromatin accessibility sites correlated with other epithelial originating TCGA tumors. Gained chromatin accessibility and BRG1 recruited sites were strongly enriched for transcription-factor-binding motifs of AP-1 family members. Furthermore, AP-1 motifs were enriched at the promoters of highly upregulated epithelial genes. Using a dominant-negative AP-1 cell line, we found that both AP-1 DNA-binding activity and BRG1 reexpression are necessary for the gene and protein expression of epithelial genes. Our study demonstrates that BRG1 reexpression drives an epithelial-like gene and protein signature in SCCOHT cells that depends upon by AP-1 activity.

BACKGROUND Surgery remains the frontline therapy for patients with localized clear cell renal cell carcinoma (ccRCC); however, 20%–40% recur. Angiogenesis inhibitors have improved survival in metastatic patients and may result in responses in the neoadjuvant setting. The impact of these agents on the tumor genetic heterogeneity or the immune milieu is largely unknown. This phase II study was designed to evaluate safety, response, and effect on tumor tissue of neoadjuvant pazopanib.METHODS ccRCC patients with localized disease received pazopanib (800 mg daily; median 8 weeks), followed by nephrectomy. Five tumors were examined for mutations by whole exome sequencing from samples collected before therapy and at nephrectomy. These samples underwent RNA sequencing; 17 samples were available for posttreatment assessment.RESULTS Twenty-one patients were enrolled. The overall response rate was 8 of 21 (38%). No patients with progressive disease. At 1-year, response-free survival and overall survival was 83% and 89%, respectively. The most frequent grade 3 toxicity was hypertension (33%, 7 of 21). Sequencing revealed strong concordance between pre- and posttreatment samples within individual tumors, suggesting tumors harbor stable core profiles. However, a reduction in private mutations followed treatment, suggesting a selective process favoring enrichment of driver mutations.CONCLUSION Neoadjuvant pazopanib is safe and active in ccRCC. Future genomic analyses may enable the segregation of driver and passenger mutations. Furthermore, tumor infiltrating immune cells persist during therapy, suggesting that pazopanib can be combined with immune checkpoint inhibitors without dampening the immune response.FUNDING Support was provided by Novartis and GlaxoSmithKline as part of an investigator-initiated study.

ABSTRACT Following inhalation, Yersinia pestis rapidly colonizes the lung to establish infection during primary pneumonic plague. Although several adhesins have been identified in Yersinia spp., the factors mediating early Y. pestis adherence in the lung remain unknown. To identify genes important for Y. pestis adherence during primary pneumonic plague, we used transposon insertion sequencing (Tn-seq). Wild-type and capsule mutant (Δcaf1) Y. pestis transposon mutant libraries were serially passaged in vivo to enrich for nonadherent mutants in the lung using a mouse model of primary pneumonic plague. Sequencing of the passaged libraries revealed six mutants that were significantly enriched in both the wild-type and Δcaf1 Y. pestis backgrounds. The enriched mutants had insertions in genes that encode transcriptional regulators, chaperones, an endoribonuclease, and YPO3903, a hypothetical protein. Using single-strain infections and a transcriptional analysis, we identified a significant role for YPO3903 in Y. pestis adherence in the lung and showed that YPO3903 regulated transcript levels of psaA, which encodes a fimbria previously implicated in Y. pestis adherence in vitro. Deletion of psaA had a minor effect on Y. pestis adherence in the lung, suggesting that YPO3903 regulates other adhesins in addition to psaA. By enriching for mutations in genes that regulate the expression or assembly of multiple genes or proteins, we obtained screen results indicating that there may be not just one dominant adhesin but rather several factors that contribute to early Y. pestis adherence during primary pneumonic plague. IMPORTANCE Colonization of the lung by Yersinia pestis is a critical first step in establishing infection during primary pneumonic plague, a disease characterized by high lethality. However, the mechanisms by which Y. pestis adheres in the lung after inhalation remain elusive. Here, we used Tn-seq to identify Y. pestis genes important for adherence early during primary pneumonic plague. Our mutant enrichment strategy resulted in the identification of genes important for regulation and assembly of genes and proteins rather than adhesin genes themselves. These results reveal that there may be multiple Y. pestis adhesins or redundancy among adhesins. Identifying the adhesins regulated by the genes identified in our enrichment screen may reveal novel therapeutic targets for preventing Y. pestis adherence and the subsequent development of pneumonic plague.

Summary: The Mre11-Rad50-Nbs1 complex is a DNA double-strand break sensor that mediates a tumor-suppressive DNA damage response (DDR) in cells undergoing oncogenic stress, yet the mechanisms underlying this effect are poorly understood. Using a genetically inducible primary mammary epithelial cell model, we demonstrate that Mre11 suppresses proliferation and DNA damage induced by diverse oncogenic drivers through a p53-independent mechanism. Breast tumorigenesis models engineered to express a hypomorphic Mre11 allele exhibit increased levels of oncogene-induced DNA damage, R-loop accumulation, and chromosomal instability with a characteristic copy number loss phenotype. Mre11 complex dysfunction is identified in a subset of human triple-negative breast cancers and is associated with increased sensitivity to DNA-damaging therapy and inhibitors of ataxia telangiectasia and Rad3 related (ATR) and poly (ADP-ribose) polymerase (PARP). Thus, deficiencies in the Mre11-dependent DDR drive proliferation and genome instability patterns in p53-deficient breast cancers and represent an opportunity for therapeutic exploitation. : The origins of genome instability in cancer remain poorly understood. Fagan-Solis et al. reveal a p53-independent genome integrity checkpoint pathway mediated by Mre11 that protects against genome instability in breast cancer. Mre11 dysfunction in breast cancer models induces a genomic loss signature and vulnerability to PARP and ATR inhibitors. Keywords: breast cancer, genome instability, chromosomal instability, DNA damage response, oncogenic stress, Mre11, R loops, genomic scar, replication stress

BackgroundRecent molecular characterization of urothelial cancer (UC) has suggested potential pathways in which to direct treatment, leading to a host of targeted therapies in development for UC. In parallel, gene expression profiling has demonstrated that high-grade UC is a heterogeneous disease. Prognostic basal-like and luminal-like subtypes have been identified and an accurate transcriptome BASE47 classifier has been developed. However, these phenotypes cannot be broadly investigated due to the lack of a clinically viable diagnostic assay. We sought to develop and evaluate a diagnostic classifier of UC subtype with the goal of accurate classification from clinically available specimens.MethodsTumor samples from 52 patients with high-grade UC were profiled for BASE47 genes concurrently by RNAseq as well as NanoString. After design and technical validation of a BASE47 NanoString probeset, results from the RNAseq and NanoString were used to translate diagnostic criteria to the Nanostring platform. Evaluation of repeatability and accuracy was performed to derive a final Nanostring based classifier. Diagnostic classification resulting from the NanoString BASE47 classifier was validated on an independent dataset (n = 30). The training and validation datasets accurately classified 87% and 93% of samples, respectively.ResultsHere we have derived a NanoString-platform BASE47 classifier that accurately predicts basal-like and luminal-like subtypes in high grade urothelial cancer. We have further validated our new NanoString BASE47 classifier on an independent dataset and confirmed high accuracy when compared with our original Transcriptome BASE47 classifier.ConclusionsThe NanoString BASE47 classifier provides a faster turnaround time, a lower cost per sample to process, and maintains the accuracy of the original subtype classifier for better clinical implementation.

Summary: We studied 137 primary testicular germ cell tumors (TGCTs) using high-dimensional assays of genomic, epigenomic, transcriptomic, and proteomic features. These tumors exhibited high aneuploidy and a paucity of somatic mutations. Somatic mutation of only three genes achieved significance—KIT, KRAS, and NRAS—exclusively in samples with seminoma components. Integrated analyses identified distinct molecular patterns that characterized the major recognized histologic subtypes of TGCT: seminoma, embryonal carcinoma, yolk sac tumor, and teratoma. Striking differences in global DNA methylation and microRNA expression between histology subtypes highlight a likely role of epigenomic processes in determining histologic fates in TGCTs. We also identified a subset of pure seminomas defined by KIT mutations, increased immune infiltration, globally demethylated DNA, and decreased KRAS copy number. We report potential biomarkers for risk stratification, such as miRNA specifically expressed in teratoma, and others with molecular diagnostic potential, such as CpH (CpA/CpC/CpT) methylation identifying embryonal carcinomas. : Shen et al. identify molecular characteristics that classify testicular germ cell tumor types, including a separate subset of seminomas defined by KIT mutations. This provides a set of candidate biomarkers for risk stratification and potential therapeutic targeting. Keywords: The Cancer Genome Atlas, testicular germ cell tumors, seminoma, nonseminoma, DNA methylation, exome sequencing, KIT, copy number, miR-375

Summary: Self-renewal and differentiation of adult stem cells are tightly regulated partly through configuration of chromatin structure by chromatin remodelers. Using knockout mice, we here demonstrate that bromodomain PHD finger transcription factor (BPTF), a component of the nucleosome remodeling factor (NURF) chromatin-remodeling complex, is essential for maintaining the population size of hematopoietic stem/progenitor cells (HSPCs), including long-term hematopoietic stem cells (HSCs). Bptf-deficient HSCs are defective in reconstituted hematopoiesis, and hematopoietic-specific knockout of Bptf caused profound defects including bone marrow failure and anemia. Genome-wide transcriptome profiling revealed that BPTF loss caused downregulation of HSC-specific gene-expression programs, which contain several master transcription factors (Meis1, Pbx1, Mn1, and Lmo2) required for HSC maintenance and self-renewal. Furthermore, we show that BPTF potentiates the chromatin accessibility of key HSC “stemness” genes. These results demonstrate an essential requirement of the chromatin remodeler BPTF and NURF for activation of “stemness” gene-expression programs and proper function of adult HSCs. : Wang and colleagues show that a chromatin remodeler, BPTF, sustains appropriate functions of hematopoietic stem/progenitor cells (HSPCs). BPTF loss causes bone marrow failure and anemia. The authors further define a BPTF-dependent gene-expression program in HSPCs, which contains key HSC stemness factors. These results demonstrate an essential requirement of the BPTF-associated chromatin remodelers for HSC functionality and adult hematopoiesis. Keywords: Bptf, hematopoietic stem cells, chromatin remodeler, Meis1, Pbx1, Mn1, DNA accessibility, NURF, AP1 complex

In metastatic breast cancer, subclonal evolution can drive drug resistance. Here, the authors genetically and transcriptionally follow the evolution of four breast cancers over time and treatment, and suggest a phenotype-targeted treatment strategy to adapt to cancer as it evolves.

Triple negative breast cancer: Obesity and metabolism fuel disease spread Metabolic changes contribute to the metastatic potential of triple negative breast cancer (TNBC), a mouse study shows. Stephen Hursting and colleagues from the University of North Carolina at Chapel Hill, USA, established metastatic mouse TNBC cells driven by Wnt-1, a signaling protein that’s highly active in this aggressive subtype of breast cancer. In a lab dish, these cells showed signs of increased invasiveness; and when transplanted into mice, the cells readily formed tumors that metastasized to the lungs. Obese mice experienced more aggressive tumor growth and spread than normal-weight animals. Gene expression analyses revealed that TNBC cells with metastatic potential have an energetic leg-up over their non-metastatic counterparts in the face of obesity-induced metabolic changes, suggesting that targeting metabolic perturbations could help reduce the burden of metastatic TNBC, particularly for obese women.

Background & Aims: Celiac disease (CeD) provides an opportunity to study autoimmunity and the transition in immune cells as dietary gluten induces small intestinal lesions. Methods: Seventy-three celiac disease patients on a long-term, gluten-free diet ingested a known amount of gluten daily for 6 weeks. A peripheral blood sample and intestinal biopsy specimens were taken before and 6 weeks after initiating the gluten challenge. Biopsy results were reported on a continuous numeric scale that measured the villus-heightâtoâcrypt-depth ratio to quantify gluten-induced intestinal injury. Pooled B and T cells were isolated from whole blood, and RNA was analyzed by DNA microarray looking for changes in peripheral B- and T-cell gene expression that correlated with changes in villus height to crypt depth, as patients maintained a relatively healthy intestinal mucosa or deteriorated in the face of a gluten challenge. Results: Gluten-dependent intestinal damage from baseline to 6 weeks varied widely across all patients, ranging from no change to extensive damage. Genes differentially expressed in B cells correlated strongly with the extent of intestinal damage. AÂ relative increase in B-cell gene expression correlated with a lack of sensitivity to gluten whereas their relative decrease correlated with gluten-induced mucosal injury. A core B-cell gene module, representing a subset of B-cell genes analyzed, accounted for the correlation with intestinal injury. Conclusions: Genes comprising the core B-cell module showed a net increase in expression from baseline to 6 weeks in patients with little to no intestinal damage, suggesting that these individuals may have mounted a B-cell immune response to maintain mucosal homeostasis and circumvent inflammation. DNA microarray data were deposited at the GEO repository (accession number: GSE87629; available: https://www.ncbi.nlm.nih.gov/geo/). Keywords: Oral Tolerance, Mucosal Immunity, Autoimmunity, Regulatory B Cell

Renal cell carcinoma (RCC) is a common and aggressive malignancy. Here, the authors generate two mouse models of the most common RCC subtypes: the human papillary RCC throughMYC activation and clear cell RCC through MYC activation combined with Vhl and Cdkn2adeletion.

Oesophageal squamous-cell carcinoma (ESCC) is a leading cause of cancer death, and half of ESCC cases occur in China. Here, the authors provide an in depth genomic landscape for this disease and identify specific mutation signatures—one of which is linked to alcohol intake.

Abstract Changes in the quantity of genetic material, known as somatic copy number alterations (CNAs), can drive tumorigenesis. Many methods exist for assessing CNAs using microarrays, but considerable technical issues limit current CNA calling based upon DNA sequencing. We present SynthEx, a novel tool for detecting CNAs from whole exome and genome sequencing. SynthEx utilizes a “synthetic-normal” strategy to overcome technical and financial issues. In terms of accuracy and precision, SynthEx is highly comparable to array-based methods and outperforms sequencing-based CNA detection tools. SynthEx robustly identifies CNAs using sequencing data without the additional costs associated with matched normal specimens.

Summary: Alpha-1 adrenergic receptors (α1-ARs) play adaptive and protective roles in the heart. Dabuzalgron is an oral selective α1A-AR agonist that was well tolerated in multiple clinical trials of treatment for urinary incontinence, but has never been used to treat heart disease in humans or animal models. In this study, the authors administered dabuzalgron to mice treated with doxorubicin (DOX), a widely used chemotherapeutic agent with dose-limiting cardiotoxicity that can lead to heart failure (HF). Dabuzalgron protected against DOX-induced cardiotoxicity, likely by preserving mitochondrial function. These results suggest that activating cardiac α1A-ARs with dabuzalgron, a well-tolerated oral agent, might represent a novel approach to treating HF. Key Words: alpha adrenergic receptors, anthracyclines, cardioprotection, catecholamines, heart failure

Summary: The incidence of human papilloma virus (HPV)-associated head and neck squamous cell carcinoma (HNSCC) is increasing and implicated in more than 60% of all oropharyngeal carcinomas (OPSCCs). Although whole-genome, transcriptome, and proteome analyses have identified altered signaling pathways in HPV-induced HNSCCs, additional tools are needed to investigate the unique pathobiology of OPSCC. Herein, bioinformatics analyses of human HPV(+) HNSCCs revealed that all tumors express full-length E6 and identified molecular subtypes based on relative E6 and E7 expression levels. To recapitulate the levels, stoichiometric ratios, and anatomic location of E6/E7 expression, we generated a genetically engineered mouse model whereby balanced expression of E6/E7 is directed to the oropharyngeal epithelium. The addition of a mutant PIK3CAE545K allele leads to the rapid development of pre-malignant lesions marked by immune cell accumulation, and a subset of these lesions progress to OPSCC. This mouse provides a faithful immunocompetent model for testing treatments and investigating mechanisms of immunosuppression. : Carper et al. present the “iKHP” mouse, in which HPV16 oncogenes are inducibly activated in vivo in a tissue-specific and temporal manner. Oropharyngeal-specific expression of E6/E7 with PIK3CAE545K in these mice promotes the development of premalignant lesions marked by immune cell infiltration, but only a subset spontaneously convert to OPSCC. Keywords: genetically engineered mouse model, GEMM, human papilloma virus, HPV, HPV16, head and neck squamous cell carcinoma, HNSCC, oropharyngeal squamous cell carcinoma, OPSCC, immunocompetent, tumorigenesis, oncogene, tumor suppressor, LumiFluor

Cigarette smoke is well recognized to cause injury to the airways and the alveolar walls over time. This injury usually requires many years of exposure, suggesting that the lungs may rapidly develop responses that initially protect it from this repetitive injury. Our studies tested the hypotheses that smoke induces an inflammatory response and changes in mRNA profiles that are dependent on sex and the health status of the lung, and that the response of the lungs to smoke differs after 1 day compared to 5 days of exposure. Male and female wildtype (WT) and Scnn1b-transgenic (βENaC) mice, which have chronic bronchitis and emphysematous changes due to dehydrated mucus, were exposed to cigarette smoke or sham air conditions for 1 or 5 days. The inflammatory response and gene expression profiles were analyzed in lung tissue. Overall, the inflammatory response to cigarette smoke was mild, and changes in mediators were more numerous after 1 than 5 days. βENaC mice had more airspace leukocytes than WT mice, and smoke exposure resulted in additional significant alterations. Many genes and gene sets responded similarly at 1 and 5 days: genes involved in oxidative stress responses were upregulated while immune response genes were downregulated. However, certain genes and biological processes were regulated differently after 1 compared to 5 days. Extracellular matrix biology genes and gene sets were upregulated after 1 day but downregulated by 5 days of smoke compared to sham exposure. There was no difference in the transcriptional response to smoke between WT and βENaC mice or between male and female mice at either 1 or 5 days. Taken together, these studies suggest that the lungs rapidly alter gene expression after only one exposure to cigarette smoke, with few additional changes after four additional days of repeated exposure. These changes may contribute to preventing lung damage.