RNA-seq analysis identifies transcriptomic profiles associated with anal cancer recurrence among people living with HIV

Abstract Background Chemoradiation therapy (CRT) is the standard of care for squamous cell carcinoma of the anus (SCCA), the most common type of anal cancer. However, approximately one fourth of patients still relapse after CRT. Methods We used RNA-sequencing technology to characterize coding and non-coding transcripts in tumor tissues from CRT-treated SCCA patients and compare them between 9 non-recurrent and 3 recurrent cases. RNA was extracted from FFPE tissues. Library preparations for RNA-sequencing were created using SMARTer Stranded Total RNA-Seq Kit. All libraries were pooled and sequenced on a NovaSeq 6000. Function and pathway enrichment analysis was performed with Metascape and enrichment of gene ontology (GO) was performed with Gene Set Enrichment Analysis (GSEA). Results There were 449 differentially expressed genes (DEGs) observed (390 mRNA, 12 miRNA, 17 lincRNA and 18 snRNA) between the two groups. We identified a core of upregulated genes (IL4, CD40LG, ICAM2, HLA-I (HLA-A, HLA-C) and HLA-II (HLA-DQA1, HLA-DRB5) in the non-recurrent SCCA tissue enriching to the gene ontology term ‘allograft rejection’, which suggests a CD4+ T cell driven immune response. Conversely, in the recurrent tissues, keratin (KRT1, 10, 12, 20) and hedgehog signaling pathway (PTCH2) genes involved in ‘Epidermis Development,’, were significantly upregulated. We identified miR-4316, that inhibit tumor proliferation and migration by repressing vascular endothelial growth factors, as being upregulated in non-recurrent SCCA. On the contrary, lncRNA-SOX21-AS1, implicated in the progression of many other cancers, was also found to be more common in our recurrent compared to non-recurrent SCCA. Conclusions Our study identified key host factors which may drive the recurrence of SCCA and warrants further studies to understand the mechanism and evaluate their potential use in personalized treatment. Key Message Our study used RNA sequencing (RNA-seq) to identify pivotal factors in coding and non-coding transcripts which differentiate between patients at risk for recurrent anal cancer after treatment. There were 449 differentially expressed genes (390 mRNA, 12 miRNA, 17 lincRNA and 18 snRNA) between 9 non-recurrent and 3 recurrent squamous cell carcinoma of anus (SCCA) tissues. The enrichment of genes related to allograft rejection was observed in the non-recurrent SCCA tissues, while the enrichment of genes related to epidermis development was positively linked with recurrent SCCA tissues.

[1]  R. Pfeiffer,et al.  Impact of HIV on Anal Squamous Cell Carcinoma Rates in the United States, 2001-2015. , 2022, Journal of the National Cancer Institute.

[2]  A. Marzi,et al.  Transcriptional Analysis of Infection With Early or Late Isolates From the 2013–2016 West Africa Ebola Virus Epidemic Does Not Suggest Attenuated Pathogenicity as a Result of Genetic Variation , 2021, Frontiers in Microbiology.

[3]  M. Fernández-Zapico,et al.  Hedgehog/GLI Signaling Pathway: Transduction, Regulation, and Implications for Disease , 2021, Cancers.

[4]  R. Lin,et al.  Identification of ferroptosis genes in immune infiltration and prognosis in thyroid papillary carcinoma using network analysis , 2021, BMC Genomics.

[5]  Neil A. Miller,et al.  Factors affecting migration to GRCh38 in laboratories performing clinical next-generation sequencing. , 2021, The Journal of molecular diagnostics : JMD.

[6]  Fangkun Liu,et al.  CD96 Correlates With Immune Infiltration and Impacts Patient Prognosis: A Pan-Cancer Analysis , 2021, Frontiers in Oncology.

[7]  Fengfeng Wang,et al.  Long non-coding RNA SOX21-AS1 modulates lung cancer progress upon microRNA miR-24-3p/PIM2 axis , 2021, Bioengineered.

[8]  Shu-juan Xie,et al.  RNA sequencing: new technologies and applications in cancer research , 2020, Journal of Hematology & Oncology.

[9]  Lanzhen Li,et al.  Long non‐coding RNA SOX21‐AS1 enhances the stemness of breast cancer cells via the Hippo pathway , 2020, FEBS open bio.

[10]  S. H. Park,et al.  The Role of Autophagy in Skin Fibroblasts, Keratinocytes, Melanocytes, and Epidermal Stem Cells. , 2020, The Journal of investigative dermatology.

[11]  Dan-dan Zhao,et al.  Identification of differentially expressed metastatic genes and their signatures to predict the overall survival of uveal melanoma patients by bioinformatics analysis. , 2020, International journal of ophthalmology.

[12]  D. Tibboel,et al.  SOX21 modulates SOX2-initiated differentiation of epithelial cells in the extrapulmonary airways , 2020, bioRxiv.

[13]  Chunlei Yuan,et al.  SNORA71B promotes breast cancer cells across blood–brain barrier by inducing epithelial-mesenchymal transition , 2020, Breast Cancer.

[14]  Qiu-ming He,et al.  LncRNA SNHG17 aggravated prostate cancer progression through regulating its homolog SNORA71B via a positive feedback loop , 2020, Cell Death & Disease.

[15]  Yue Xu,et al.  MicroRNA-4316 inhibits gastric cancer proliferation and migration via directly targeting VEGF-A , 2020, Cancer Cell International.

[16]  Linjia Su,et al.  Arf6-driven endocytic recycling of CD147 determines HCC malignant phenotypes , 2019, Journal of Experimental & Clinical Cancer Research.

[17]  J. Palefsky,et al.  Recent trends in squamous cell carcinoma of the anus incidence and mortality in the United States, 2001-2015. , 2019, Journal of the National Cancer Institute.

[18]  Yan Zhou,et al.  Long noncoding RNA SOX21‐AS1 promotes cervical cancer progression by competitively sponging miR‐7/VDAC1 , 2019, Journal of cellular physiology.

[19]  Steven L Salzberg,et al.  Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype , 2019, Nature Biotechnology.

[20]  Alireza Hadj Khodabakhshi,et al.  Metascape provides a biologist-oriented resource for the analysis of systems-level datasets , 2019, Nature Communications.

[21]  Yu Yang,et al.  LncRNA CALML3-AS1 promotes tumorigenesis of bladder cancer via regulating ZBTB2 by suppression of microRNA-4316. , 2018, Biochemical and biophysical research communications.

[22]  F. Wang,et al.  Circular RNA circMYO9B facilitates breast cancer cell proliferation and invasiveness via upregulating FOXP4 expression by sponging miR-4316. , 2018, Archives of biochemistry and biophysics.

[23]  Benno Schwikowski,et al.  The Cytoscape Automation app article collection , 2018, F1000Research.

[24]  Y. Chung,et al.  Whole-exome sequencing identified mutational profiles of squamous cell carcinomas of anus. , 2018, Human pathology.

[25]  Harrison Pielke-Lombardo,et al.  GSEA-InContext: identifying novel and common patterns in expression experiments , 2018, bioRxiv.

[26]  R. Pfeiffer,et al.  Anal Cancer Risk Among People With HIV Infection in the United States. , 2018, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[27]  D. Hussein,et al.  Comprehensive molecular biomarker identification in breast cancer brain metastases , 2017, Journal of Translational Medicine.

[28]  Jing Wang,et al.  Comprehensive Genomic Profiling of Metastatic Squamous Cell Carcinoma of the Anal Canal , 2017, Molecular Cancer Research.

[29]  Toby C. Cornish,et al.  Association of HIV Status With Local Immune Response to Anal Squamous Cell Carcinoma: Implications for Immunotherapy , 2017, JAMA oncology.

[30]  E. Bahassi,et al.  Comparative molecular profiling of HPV‐induced squamous cell carcinomas , 2017, Cancer medicine.

[31]  P. Krawczyk,et al.  Evaluation of miR-506 and miR-4316 expression in early and non-invasive diagnosis of colorectal cancer , 2017, International Journal of Colorectal Disease.

[32]  W. Dougall,et al.  TIGIT and CD96: new checkpoint receptor targets for cancer immunotherapy , 2017, Immunological reviews.

[33]  Z. Lu,et al.  RNA Biomarkers: Frontier of Precision Medicine for Cancer , 2017, Non-coding RNA.

[34]  Brian C. Jackson,et al.  Regulation of Head and Neck Squamous Cancer Stem Cells by PI3K and SOX2 , 2017, Journal of the National Cancer Institute.

[35]  P. Agostinis,et al.  DAMP—Induced Allograft and Tumor Rejection: The Circle Is Closing , 2016, American journal of transplantation : official journal of the American Society of Transplantation and the American Society of Transplant Surgeons.

[36]  Y. Shiue,et al.  Aberrant DNA hypermethylation-silenced SOX21-AS1 gene expression and its clinical importance in oral cancer , 2016, Clinical Epigenetics.

[37]  G. Getz,et al.  Genomic Evolution after Chemoradiotherapy in Anal Squamous Cell Carcinoma , 2016, Clinical Cancer Research.

[38]  H. Katki,et al.  Cancer risk among the HIV-infected elderly in the United States , 2016, AIDS.

[39]  James X. Sun,et al.  Comprehensive genomic profiling of anal squamous cell carcinoma reveals distinct genomically defined classes. , 2016, Annals of oncology : official journal of the European Society for Medical Oncology.

[40]  Z. Gatalica,et al.  Comprehensive multiplatform biomarker analysis of 199 anal squamous cell carcinomas , 2015, Oncotarget.

[41]  R. Pfeiffer,et al.  Excess cancers among HIV-infected people in the United States. , 2015, Journal of the National Cancer Institute.

[42]  H. Kourie,et al.  Anal cancer treatment: current status and future perspectives. , 2015, World journal of gastroenterology.

[43]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[44]  T. Fujiwara,et al.  Peptide intra-tumor injection for cancer immunotherapy , 2013, Human vaccines & immunotherapeutics.

[45]  G. Srikrishna S100A8 and S100A9: New Insights into Their Roles in Malignancy , 2011, Journal of Innate Immunity.

[46]  J. Ledermann,et al.  Chemoradiation for the treatment of epidermoid anal cancer: 13-year follow-up of the first randomised UKCCCR Anal Cancer Trial (ACT I) , 2010, British Journal of Cancer.

[47]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[48]  W. Chazin,et al.  S100A8/A9: a Janus-faced molecule in cancer therapy and tumorgenesis. , 2009, European journal of pharmacology.

[49]  S. Cole,et al.  A Meta-Analysis of the Incidence of Non-AIDS Cancers in HIV-Infected Individuals , 2009, Journal of acquired immune deficiency syndromes.

[50]  D. Foell,et al.  The endogenous Toll–like receptor 4 agonist S100A8/S100A9 (calprotectin) as innate amplifier of infection, autoimmunity, and cancer , 2009, Journal of leukocyte biology.

[51]  R. Biggar,et al.  Risk of human papillomavirus-associated cancers among persons with AIDS. , 2009, Journal of the National Cancer Institute.

[52]  J. Bendell,et al.  Anal cancer: an overview. , 2007, The oncologist.

[53]  P. Angel,et al.  S100A8 and S100A9 in inflammation and cancer. , 2006, Biochemical pharmacology.

[54]  M. Katoh,et al.  Hedgehog signaling pathway and gastric cancer , 2005, Cancer biology & therapy.

[55]  M. Stamos,et al.  Patterns of recurrence in anal canal carcinoma. , 2000, Archives of surgery.

[56]  R. Hollingsworth,et al.  PTCH2, a novel human patched gene, undergoing alternative splicing and up-regulated in basal cell carcinomas. , 1999, Cancer research.

[57]  H. Bartelink,et al.  Concomitant radiotherapy and chemotherapy is superior to radiotherapy alone in the treatment of locally advanced anal cancer: results of a phase III randomized trial of the European Organization for Research and Treatment of Cancer Radiotherapy and Gastrointestinal Cooperative Groups. , 1997, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[58]  R. James Epidermoid anal cancer: results from the UKCCCR randomised trial of radiotherapy alone versus radiotherapy, 5-fluorouracil, and mitomycin , 1996, The Lancet.

[59]  G. Tanum [Anal cancer. A review]. , 1996, Tidsskrift for den Norske laegeforening : tidsskrift for praktisk medicin, ny raekke.