Chanana, Pritha

Abstract Background: DNA polymerase epsilon (POLE) is encoded by the POLE gene, and POLE-driven tumors are characterized by high mutational rates. POLE-driven tumors are relatively common in endometrial and colorectal cancer, and their presence is increasingly recognized in ovarian cancer (OC) of endometrioid type. POLE-driven cases possess an abundance of TCT>TAT and TCG>TTG somatic mutations characterized by mutational signature 10 from the Catalog of Somatic Mutations in Cancer (COSMIC). By quantifying the contribution of COSMIC mutational signature 10 in RNA sequencing (RNA-seq) we set out to identify POLE-driven tumors in a set of unselected Mayo Clinic OC. Methods: Mutational profiles were calculated using expressed single-nucleotide variants (eSNV) in the Mayo Clinic OC tumors (n=195), The Cancer Genome Atlas (TCGA) OC tumors (n=419), and the Genotype-Tissue Expression (GTEx) normal ovarian tissues (n=84). Non-negative Matrix Factorization (NMF) of the mutational profiles inferred the contribution per sample of four distinct mutational signatures, one of which corresponds to COSMIC mutational signature 10. Results: In the Mayo Clinic OC cohort we identified six tumors with a predicted contribution from COSMIC mutational signature 10 of over five mutations per megabase. These six cases harbored known POLE hotspot mutations (P286R, S297F, V411L, and A456P) and were of endometrioid histotype (P=5e-04). These six tumors were hypermutated with a higher tumor mutation load (mean, 54.02 mutations per megabase) compared to non-POLE endometrioid OC cases (mean, 7.69 mutations per megabase; P=5e-04), and had an early onset (average age of patients at onset, 48.33 years) when compared to non-POLE endometrioid OC cohort (average age at onset, 60.13 years; P=.008). Samples from TCGA and GTEx had a low COSMIC signature 10 contribution (median 0.16 mutations per megabase; maximum 1.78 mutations per megabase) and carried no POLE hotspot mutations.Conclusions: From the largest cohort of RNA-seq from endometrioid OC to date (n=53), we identified six hypermutated samples likely driven by POLE (frequency, 11%). Our result suggests the clinical need to screen for POLE driver mutations in endometrioid OC, which can guide enrollment in immunotherapy clinical trials.

Abstract Introduction: Acquired resistance to platinum remains a major clinical problem in ovarian cancer since systemic therapies have limited efficacy in this setting. There is growing evidence that the tumor microenvironment (TME) contributes to such resistance, as well as the pathogenesis of the disease. The study of the TME is limited by technical challenges of separating the epithelial from the stromal component. To overcome this limitation, we performed RNA sequencing (RNA-Seq) for three pairs of sensitive and induced resistant patient-derived xenograft models since the epithelial transcripts originate from human tissue and the stroma component is murine. Methods: To determine which TME genes might be contributing to chemotherapy resistance, three platinum-sensitive ovarian cancer PDX lines were treated repeatedly with carboplatin/paclitaxel until resistance occurred. Primary platinum-sensitive and induced platinum resistant tumors were characterized histologically by H&E and Masson Trichrome staining. The three pairs of sensitive and resistant PDX models were then subjected to RNA-Seq. Raw sequencing data were processed through Xenome to classify the reads into epithelial carcinoma (human) or TME (mouse) through Mayo Clinic’s internal RNA-Seq pipeline for transcriptome quantification. Differential expression analysis between induced resistant and sensitive models was carried out for both carcinoma and TME components using edgeR paired analysis. Gene set enrichment analysis was used to identify overrepresented biologic functions for differentially expressed genes between the pairs. Top significantly varied, differentially expressed genes between the pairs were validated using species-specific qPCR. Results: Histologic analysis of the matched sensitive/resistant pairs revealed a marked increase of the stromal component in the resistant tumors evident by H&E and Masson trichrome staining. A total number of 95 and 312 genes were differentially expressed in the epithelial and TME component, respectively (FDR<0.05; absolute log2 fold-change >1.5). The stemness marker SOX2 was significantly upregulated in the human component of the resistant tumors as confirmed by qPCR. Genes involved in stroma remodeling were upregulated in the “murine” component of the resistant tumors, and upregulation was confirmed by qPCR (normalized to a murine housekeeper gene). Gene set enrichment analysis showed that the upregulated genes in TME of resistant models were enriched for cytoskeleton related functions, while the downregulated genes were enriched for extracellular matrix-related functions. Conclusion: A bioinformatics strategy able to map human and murine reads separately may identify resistance mechanisms attributed to the two tumor subcompartments and could lead to strategies to potentially target each separately or concurrently. Citation Format: Valentina Zanfagnin, Yuanhang Liu, Pritha Chanana, Xiaonan Pritha, Beverly Long, Joe Pathoulas, Brad Evans, Jaime Davila, Chen Wang, Saravut John Weroha. Transcriptional profiling of tumor stroma using ovarian cancer PDX models with induced platinum-resistance [abstract]. In: Proceedings of the AACR Special Conference on Advances in Ovarian Cancer Research; 2019 Sep 13-16, 2019; Atlanta, GA. Philadelphia (PA): AACR; Clin Cancer Res 2020;26(13_Suppl):Abstract nr B18.

Antibody secreting plasma cells (PCs) play an important role in effective humoral immune responses. The low frequency of bone marrow PCs in humans makes it challenging to obtain sufficient numbers of PCs for biologic studies. Previous studies have employed in vitro model systems to generate cells that morphologically, phenotypically, and functionally resemble normal polyclonal PCs. Gene expression profiles of in vitro generated PCs (IVPCs) mirror their normal counterparts, however to date extensive immunoglobulin (Ig) repertoire analysis of IVPCs is lacking. Here, we used a modified 3-step protocol to generate IVPCs and used RNA-seq to explore the transcriptome with emphasis on the Ig repertoire of plasmablasts and PCs. Total B cells were isolated from 3 normal donors and cultured with various cytokines and the B cell activators CpG ODN and CD40L. RNA was obtained from freshly isolated B cells (Day 0; D0) as well as from Day 4 (D4) plasmablasts, and Day 10 (D10) IVPCs. Morphologically, D10 cells exhibited typical PC morphology, including an eccentric nucleus and perinuclear hof. RNA-seq was performed on total RNA from all 3 donors and time points using the Standard TRuSeq v2 library prep and with paired end sequencing on the Illumina HiSeq 4000 platform. Principle component analysis of gene expression data showed that D0, D4 and D10 cells could be clearly segregated across all 3 normal donors. Of importance, transcripts previously described as distinguishing B cells from PCs were found to be differentially expressed including overexpression of CXCR5, CD19, EBF, CD83, PAX5, IRF8 in D0 B cells and overexpression of IRF4, Blimp-1, XBP1, BCMA, SLAMF7, Syndecan-1, CD38 and CD27 in IVPCs, thus validating our in vitro model for generating PCs. Furthermore, expression of cell cycle related transcripts such as CKS1, CDK1, and CCDN2 followed the pattern of low expression in resting B cells, increased expression in plasmablasts, and decreased expression in IVPCs confirming the cells are actively cycling in a manner comparable to cells in vivo. D10 IVPCs also overexpressed transcripts known to be upregulated during the unfolded protein response. As expected from Ig secreting cells, D10 IVPCs had an over-representation of Ig transcripts. At D0, resting B cells had high levels of IgD and IgM heavy chain (HC) transcripts. At D10, IgM transcripts modestly increased with Log2 fold change (FC) = 3 and as expected, IgD levels decreased significantly (Log2 FC = -2.2). IgA and IgG isotype transcripts significantly increased at D10 (Log2 FC > 6.0) with the IgG4 subtype having the greatest Log2 FC at 8.4. Next we focused on the Ig repertoire of D0, D4, and D10 cells. By aligning to known germline Ig sequences in IMGT/V-Quest (www.imgt.org) and then assembling the paired ends of D0, D4 and D10 Ig transcripts, we were able to analyze the Ig repertoire. Since the Ig HC variable (V) region is encoded by V, diversity (D) and joining (J) segments, only fragments that could be confidently determined were considered. All but 3 IGHV transcripts (IGHV3-35, IGHV3-47 and IGHV7-8) and 2 IGHD transcripts (IGHD4-4 and IGHD5-5) were found and all IGHJ segments were represented across the differentiation spectrum. In D0 cells, the number of unique VDJ combinations ranged from 643 to 863 across all 3 normal samples and increased to a range of 2524 to 2867 in D10 IVPCs. When looking at the differential expression of each VDJ combination from D0 to D10, a pairwise t-test for relative frequency showed that there was no significant change greater than 1%, suggesting the repertoire diversity was not skewed, thus proving the conditions for stimulation were not targeting any one starting B cell. Our data also allowed us to track clonal expansions during differentiation as defined by the increasing frequency of sequences with identical nucleotide sequence in the V region and CDR3 (including D and J regions). Hence, a single sequence could be tracked from D0 to D10. Of interest, in a small sampling of the total available sequences, only those B cells with a mutated IGHV region, characteristic of a memory B cell, went on to expand in this system whereas B cells with an unmutated IGHV did not. Our analysis of the Ig repertoire of IVPCs suggests this system provides a functional model to study Ig repertoire along the B cell differentiation process and further delineate the conditions that may result in a clonal expansion, a hallmark of many hematologic malignancies including multiple myeloma. Disclosures No relevant conflicts of interest to declare.

CD4 and CD8 T cells are vital components of the immune system. We found that histone deacetylase 3 (HDAC3) is critical for the development of CD4 T cells, as HDAC3-deficient DP thymocytes generate only CD8SP thymocytes in mice. In the absence of HDAC3, MHC Class II-restricted OT-II thymocytes are redirected to the CD8 cytotoxic lineage, which occurs with accelerated kinetics. Analysis of histone acetylation and RNA-seq reveals that HDAC3-deficient DP thymocytes are biased towards the CD8 lineage prior to positive selection. Commitment to the CD4 or CD8 lineage is determined by whether persistent TCR signaling or cytokine signaling predominates, respectively. Despite elevated IL-21R/γc/STAT5 signaling in HDAC3-deficient DP thymocytes, blocking IL-21R does not restore CD4 lineage commitment. Instead, HDAC3 binds directly to CD8-lineage promoting genes. Thus, HDAC3 is required to restrain CD8-lineage genes in DP thymocytes for the generation of CD4 T cells.

Although immune checkpoint inhibitors have resulted in durable clinical benefits in a subset of patients with advanced cancer, some patients who did not respond to initial anti–PD-1 therapy have been found to benefit from the addition of salvage chemotherapy. However, the mechanism responsible for the successful chemoimmunotherapy is not completely understood. Here we show that a subset of circulating CD8 + T cells expressing the chemokine receptor CX3CR1 are able to withstand the toxicity of chemotherapy and are increased in patients with metastatic melanoma who responded to chemoimmunotherapy (paclitaxel and carboplatin plus PD-1 blockade). These CX3CR1 + CD8 + T cells have effector memory phenotypes and the ability to efflux chemotherapy drugs via the ABCB1 transporter. In line with clinical observation, our preclinical models identified an optimal sequencing of chemoimmunotherapy that resulted in an increase of CX3CR1 + CD8 + T cells. Taken together, we found a subset of PD-1 therapy–responsive CD8 + T cells that were capable of withstanding chemotherapy and executing tumor rejection with their unique abilities of drug efflux (ABCB1), cytolytic activity (granzyme B and perforin), and migration to and retention (CX3CR1 and CD11a) at tumor sites. Future strategies to monitor and increase the frequency of CX3CR1 + CD8 + T cells may help to design effective chemoimmunotherapy to overcome cancer resistance to immune checkpoint blockade therapy.

Transplantation of autologous mesenchymal stem cells (MSCs) may be a viable option for treatment of several diseases. MSCs efficacy depends on adequate function of their mitochondria, which might be impaired in a noxious milieu. We hypothesized that obesity compromises MSCs mitochondrial structure and function, possibly via micro‐RNA (miRNA)‐based mechanisms. MSCs were collected from swine abdominal adipose tissue after 16 weeks of Lean or Obese diet (n = 7 each). Mitochondrial structure was assessed by electron microscopy and function by membrane potential and cytochrome‐c oxidase (COX)‐IV activity. Oxidative stress was assessed by Mito‐SOX and dihydroethidium staining. Next‐generation sequencing (RNA‐seq) was performed to identify miRNAs expression in MSCs, and predicted mitochondrial target genes were then identified (MitoCarta). Compared to Lean‐MSCs, mitochondria from Obese‐MSCs were smaller and showed cristae remodeling and loss. Mitochondrial membrane potential and COX‐IV activity decreased in Obese‐MSCs, associated with increased mitochondrial oxidative stress. RNA‐seq generated reads for 413 miRNAs, of which 5 miRNAs were upregulated in Obese‐MSCs (fold change >2, p < 0.05) and found to target 43 specific mitochondrial genes. Obesity impairs MSC mitochondrial structure and function, possibly mediated partly through miRNA‐induced mitochondrial gene regulation, leading to increased oxidative stress. Importantly, these alterations may limit the therapeutic use of autologous MSCs in subjects with obesity.

AbstractAutologous transplantation of mesenchymal stem cells (MSCs) is a viable option for the treatment of several diseases. Evidence indicates that MSCs release extracellular vesicles (EVs) and that EVs shuttle miRNAs to damaged parenchymal cells to activate an endogenous repair program. We hypothesize that comorbidities may interfere with the packaging of cargo in MSC‐derived EVs. Therefore, we examined whether metabolic syndrome (MetS) modulates the miRNA content packed within MSC‐derived EVs. MSCs were collected from swine abdominal adipose tissue after 16 weeks of lean or obese diet (n = 7 each). Next‐generation RNA sequencing of miRNAs (miRNA‐seq) was performed to identify miRNAs enriched in MSC‐derived EVs and their predicted target genes. Functional pathway analysis of the top 50 target genes of the top 4 miRNAs enriched in each group was performed using gene ontology analysis. Lean‐ and MetS‐EVs were enriched in, respectively, 14 and 8 distinct miRNAs. Target genes of miRNAs enriched in MetS‐EVs were implicated in the development of MetS and its complications, including diabetes‐related pathways, validated transcriptional targets of AP1 family members Fra1 and Fra2, Class A/1 (Rhodopsin‐like receptors), and Peptide ligand‐binding receptors. In contrast, miRNAs enriched in Lean EVs target primarily EphrinA‐EPHA and the Rho family of GTPases. MetS alters the miRNA content of EVs derived from porcine adipose tissue MSCs. These alterations could impair the efficacy and limit the therapeutic use of autologous MSCs in subjects with MetS. Our findings may assist in developing adequate regenerative strategies to preserve the reparative potency of MSCs in individuals with MetS. © 2017 International Society for Advancement of Cytometry

AbstractIntroductionWe hypothesized that common Alzheimer's disease (AD)‐associated variants within the triggering receptor expressed on myeloid (TREM) gene cluster influence disease through gene expression.MethodsExpression microarrays on temporal cortex and cerebellum from ∼400 neuropathologically diagnosed subjects and two independent RNAseq replication cohorts were used for expression quantitative trait locus analysis.ResultsA variant within a DNase hypersensitive site 5′ of TREM2, rs9357347‐C, associates with reduced AD risk and increased TREML1 and TREM2 levels (uncorrected P = 6.3 × 10−3 and 4.6 × 10−2, respectively). Meta‐analysis on expression quantitative trait locus results from three independent data sets (n = 1006) confirmed these associations (uncorrected P = 3.4 × 10−2 and 3.5 × 10−3, Bonferroni‐corrected P = 6.7 × 10−2 and 7.1 × 10−3, respectively).DiscussionOur findings point to rs9357347 as a functional regulatory variant that contributes to a protective effect observed at the TREM locus in the International Genomics of Alzheimer's Project genome‐wide association study meta‐analysis and suggest concomitant increase in TREML1 and TREM2 brain levels as a potential mechanism for protection from AD.