Liu, Yuanhang

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.

Abstract Objectives Mucinous cystic neoplasm of the liver is characterized by neoplastic mucinous and/or biliary epithelium surrounded by ovarian-type stroma. Immunohistochemical studies have shown that the ovarian-type stroma expresses estrogen receptor, suggesting potential hormonal responsiveness. The molecular biology of mucinous cystic neoplasm of the liver remains poorly studied. Methods Transcriptome sequencing and immunohistochemistry were performed on a series of mucinous cystic neoplasms. Results Mucinous cystic neoplasm of the liver exhibited significantly increased RNA expression of ovarian stromal markers WT1, PR, and ER2 and sex cord stromal markers SF-1, inhibin-α, and calretinin compared with nonneoplastic liver. Immunohistochemistry confirmed the RNA-level data. Evidence for sex hormone biosynthesis was identified by significant overexpression of multiple estrogen biosynthetic enzymes. Expression of 17β-hydroxysteroid dehydrogenase 1 was confirmed immunohistochemically. Pathway analysis also identified significant upregulation of the hedgehog and Wnt pathways and significant downregulation of T-helper 1 and T-helper 2 pathways. Conclusions Mucinous cystic neoplasm of the liver recapitulates ovarian stroma at the morphologic, DNA, RNA, and protein levels. These data support the concept that this tumor likely arises from ectopic primitive gonadal tissue and/or stromal cells with capacity to transdifferentiate to ovarian cortical cells.

Abstract Background: Traditionally, mutational burden and mutational signatures have been assessed by tumor-normal pair DNA sequencing. The requirement of having both normal and tumor samples is not always feasible from a clinical perspective, and led us to investigate the efficacy of using RNA sequencing of only the tumor sample to determine the mutational burden and signatures, and subsequently molecular cause of the cancer. The potential advantages include reducing the cost of testing, and simultaneously providing information on the gene expression profile and gene fusions present in the tumor.Results: In this study, we devised supervised and unsupervised learning methods to determine mutational signatures from tumor RNA-seq data. As applications, we applied the methods to a training set of 587 TCGA uterine cancer RNAseq samples, and examined in an independent testing set of xxx TCGA colorectal cancer RNAseq samples. Both diseases are known associated with microsatellite instable high (MSI-H) and driver defects in DNA polymerase ɛ (POLɛ). From RNAseq called variants, we found majority (>95%) are likely germline variants, leading to C>T enriched germline variants (dbSNP) widely applicable in tumor and normal RNAseq samples. We found significant associations between RNA-derived mutational burdens and MSI/POLɛ status, and insignificant relationship between RNAseq total coverage and derived mutational burdens. Additionally we found that over 80% of variants could be explained by using the COSMIC mutational signature-5, -6 and -10, which are implicated in natural aging, MSI-H, and POLɛ, respectively. For classifying tumor type, within UCEC we achieved a recall of 0.56 and 0.78, and specificity of 0.66 and 0.99 for MSI and POLɛ respectively. By applying learnt RNA signatures from UCEC to COAD, we were able to improve our classification of both MSI and POLɛ. Conclusions: Taken together, our work provides a novel method to detect RNAseq derived mutational signatures with effective procedures to remove likely germline variants. It can leads to accurate classification of underlying driving mechanisms of DNA damage deficiency.

Background: Scattered tubular-like cells (STCs) proliferate and differentiate to support neighboring injured renal tubular cells during recovery from insults. Renal artery stenosis (RAS) induces renal ischemia and hypertension and leads to loss of kidney function, but whether RAS alters renal endogenous repair mechanisms, such as STCs, remains unknown. We hypothesize that RAS in swine modifies the messenger RNA (mRNA) profile of STCs, blunting their in vitro reparative capacity. Methods: CD24+/CD133+ STCs were isolated from pig kidneys after 10-weeks of RAS or sham (n = 3 each) and their gene cargo analyzed using high-throughput mRNAseq. Expression profiles for upregulated and downregulated mRNAs in RAS-STCs were functionally interpreted by gene ontology analysis. STC activation was assessed by counting the total number of STCs in pig kidney sections using flow cytometry, whereas cell proliferation was assessed in vitro. Results: Of all expressed genes, 1430 genes were upregulated and 315 downregulated in RAS- versus Normal-STCs. Expression of selected candidate genes followed the same fold change directions as the mRNAseq findings. Genes upregulated in RAS-STCs were involved in cell adhesion, extracellular matrix remodeling, and kidney development, whereas those downregulated in RAS-STCs are related to cell cycle and cytoskeleton. The percentage of STCs from dissociated kidney cells was higher in RAS versus Normal pigs, but their proliferation rate was blunted. Conclusions: Renal ischemia and hypertension in swine induce changes in the mRNA profile of STCs, associated with increased STC activation and impaired proliferation. These observations suggest that RAS may alter the reparative capacity of STCs.

Abstract Background: Archived formalin fixed paraffin embedded (FFPE) samples are valuable clinical resources to examine clinically relevant morphology features and also to study genetic changes. However, DNA quality and quantity of FFPE samples are sub-optimal, and resulting NGS-based genetics variant detections are prone to false positives. Evaluations of wet-lab and bioinformatics approaches are needed to optimize variant detection from FFPE samples. Results: As a pilot study, we designed within-subject triplicate samples of DNA derived from paired FFPE and fresh frozen breast tissues to highlight FFPE-specific artifacts. For FFPE samples, we tested two FFPE DNA extraction methods to determine impact of wet-lab procedures on variant calling: QIAGEN QIAamp DNA Mini Kit ("QA"), and QIAGEN GeneRead DNA FFPE Kit ("QGR"). We also used negative-control (NA12891) and positive control samples (Horizon Discovery Reference Standard FFPE). All DNA sample libraries were prepared for NGS according to the QIAseq Human Breast Cancer Targeted DNA Panel protocol and sequenced on the HiSeq 4000. Variant calling and filtering were performed using QIAGEN Gene Globe Data Portal. Detailed variant concordance comparisons and mutational signature analysis were performed to investigate effects of FFPE samples compared to paired fresh frozen samples, along with different DNA extraction methods. In this study, we found that five times or more variants were called with FFPE samples, compared to their paired fresh-frozen tissue samples even after applying molecular barcoding error-correction and default bioinformatics filtering recommended by the vendor. We also found that QGR as an optimized FFPE-DNA extraction approach leads to much fewer discordant variants between paired fresh frozen and FFPE samples. Approximately 92 % of the uniquely called FFPE variants were of low allelic frequency range (<5%), and collectively shared a “C>T|G>A” mutational signature known to be representative of FFPE artifacts resulting from cytosine deamination. Based on control samples and FFPE-frozen replicates, we derived an effective filtering strategy with associated empirical false-discovery estimates. Conclusions: Through this study, we demonstrated feasibility of calling and filtering genetic variants from FFPE tissue samples using a combined strategy with molecular barcodes, optimized DNA extraction, and bioinformatics methods incorporating genomics context such as mutational signature and variant allelic frequency.

Abstract Background Archived formalin fixed paraffin embedded (FFPE) samples are valuable clinical resources to examine clinically relevant morphology features and also to study genetic changes. However, DNA quality and quantity of FFPE samples are often sub-optimal, and resulting NGS-based genetics variant detections are prone to false positives. Evaluations of wet-lab and bioinformatics approaches are needed to optimize variant detection from FFPE samples. Results As a pilot study, we designed within-subject triplicate samples of DNA derived from paired FFPE and fresh frozen breast tissues to highlight FFPE-specific artifacts. For FFPE samples, we tested two FFPE DNA extraction methods to determine impact of wet-lab procedures on variant calling: QIAGEN QIAamp DNA Mini Kit ("QA"), and QIAGEN GeneRead DNA FFPE Kit ("QGR"). We also used negative-control (NA12891) and positive control samples (Horizon Discovery Reference Standard FFPE). All DNA sample libraries were prepared for NGS according to the QIAseq Human Breast Cancer Targeted DNA Panel protocol and sequenced on the HiSeq 4000. Variant calling and filtering were performed using QIAGEN Gene Globe Data Portal. Detailed variant concordance comparisons and mutational signature analysis were performed to investigate effects of FFPE samples compared to paired fresh frozen samples, along with different library preparations. In this study, we found that five times or more variants were called with FFPE samples, compared to their paired fresh-frozen tissue samples even after applying molecular barcoding error-correction and default bioinformatics filtering recommended by the vendor. We also found that QGR as an optimized FFPE-DNA extraction approach leads to much fewer discordant variants between paired fresh frozen and FFPE samples. Approximately 92% of the uniquely called FFPE variants were of low allelic frequency range (<5%), and collectively shared a “C>T|G>A” mutational signature known to be representative of FFPE artifacts resulting fromcytosine deamination. Based on control samples and FFPE-frozen replicates, we derived an effective filtering strategy with associated empirical false-discovery estimates. Conclusions Through this study, we demonstrated feasibility of calling and filtering genetic variants from FFPE tissue samples using a combined strategy with molecular barcodes, optimized DNA extraction, and bioinformatics methods incorporating genomics context such as mutational signature and variant allelic frequency.

Abstract Background Archived formalin fixed paraffin embedded (FFPE) samples are valuable clinical resources to examine clinically relevant morphology features and also to study genetic changes. However, DNA quality and quantity of FFPE samples are often sub-optimal, and resulting NGS-based genetics variant detections are prone to false positives. Evaluations of wet-lab and bioinformatics approaches are needed to optimize variant detection from FFPE samples. Results As a pilot study, we designed within-subject triplicate samples of DNA derived from paired FFPE and fresh frozen breast tissues to highlight FFPE-specific artifacts. For FFPE samples, we tested two FFPE DNA extraction methods to determine impact of wet-lab procedures on variant calling: QIAGEN QIAamp DNA Mini Kit ("QA"), and QIAGEN GeneRead DNA FFPE Kit ("QGR"). We also used negative-control (NA12891) and positive control samples (Horizon Discovery Reference Standard FFPE). All DNA sample libraries were prepared for NGS according to the QIAseq Human Breast Cancer Targeted DNA Panel protocol and sequenced on the HiSeq 4000. Variant calling and filtering were performed using QIAGEN Gene Globe Data Portal. Detailed variant concordance comparisons and mutational signature analysis were performed to investigate effects of FFPE samples compared to paired fresh frozen samples, along with different library preparations. In this study, we found that five times or more variants were called with FFPE samples, compared to their paired fresh-frozen tissue samples even after applying molecular barcoding error-correction and default bioinformatics filtering recommended by the vendor. We also found that QGR as an optimized FFPE-DNA extraction approach leads to much fewer discordant variants between paired fresh frozen and FFPE samples. Approximately 92% of the uniquely called FFPE variants were of low allelic frequency range (<5%), and collectively shared a “C>T|G>A” mutational signature known to be representative of FFPE artifacts resulting fromcytosine deamination. Based on control samples and FFPE-frozen replicates, we derived an effective filtering strategy with associated empirical false-discovery estimates. Conclusions Through this study, we demonstrated feasibility of calling and filtering genetic variants from FFPE tissue samples using a combined strategy with molecular barcodes, optimized DNA extraction, and bioinformatics methods incorporating genomics context such as mutational signature and variant allelic frequency.

Ex vivo CRISPR/Cas9‐mediated gene editing in hepatocytes using homology‐directed repair (HDR) is a potential alternative curative therapy to organ transplantation for metabolic liver disease. However, a major limitation of this approach in quiescent adult primary hepatocytes is that nonhomologous end‐joining is the predominant DNA repair pathway for double‐strand breaks (DSBs). This study explored the hypothesis that ex vivo hepatocyte culture could reprogram hepatocytes to favor HDR after CRISPR/Cas9‐mediated DNA DSBs. Quantitative PCR (qPCR), RNA sequencing, and flow cytometry demonstrated that within 24 hours, primary mouse hepatocytes in ex vivo monolayer culture decreased metabolic functions and increased expression of genes related to mitosis progression and HDR. Despite the down‐regulation of hepatocyte function genes, hepatocytes cultured for up to 72 hours could robustly engraft in vivo. To assess functionality long‐term, primary hepatocytes from a mouse model of hereditary tyrosinemia type 1 bearing a single‐point mutation were transduced ex vivo with two adeno‐associated viral vectors to deliver the Cas9 nuclease, target guide RNAs, and a 1.2‐kb homology template. Adeno‐associated viral Cas9 induced robust cutting at the target locus, and, after delivery of the repair template, precise correction of the point mutation occurred by HDR. Edited hepatocytes were transplanted into recipient fumarylacetoacetate hydrolase knockout mice, resulting in engraftment, robust proliferation, and prevention of liver failure. Weight gain and biochemical assessment revealed normalization of metabolic function. Conclusion: The results of this study demonstrate the potential therapeutic effect of ex vivo hepatocyte‐directed gene editing after reprogramming to cure metabolic disease in a preclinical model of hereditary tyrosinemia type 1.