Abstract
Tumours with a high mutation burden exhibit considerable neoantigens and tumour-infiltrating lymphocytes. RNA editing by ADAR1 is a source of changes in epitope. However, ADAR1 expression in cancer cells and tumour-infiltrating lymphocyte levels in triple-negative breast cancer have not been well evaluated. We immunohistochemically examined ADAR1 expression in 681 triple-negative breast cancer patients and analysed their clinicopathological characteristics. We also analysed basal-like tumours using The Cancer Genome Atlas data. Among the 681 triple-negative breast cancer patients, 45.8% demonstrated high ADAR1 expression. Tumours with high ADAR1 expression exhibited high tumour-infiltrating lymphocyte levels, considerable CD8 + T lymphocyte infiltration, high histological grade and high expression of interferon-related proteins, including HLA-ABC, MxA and PKR. Among patients with lymph node metastasis, those with high tumour-infiltrating lymphocyte levels and low ADAR1 expression demonstrated the best disease-free survival. The Cancer Genome Atlas data analysis of basal-like tumours revealed significant positive correlation between
Keywords
Introduction
Triple-negative breast cancer (TNBC) accounts for ~15% of all breast cancers and are associated with high histological grade and aggressive clinical behaviour. These cancers are defined by the absence of expression of oestrogen receptor (ER), progesterone receptor (PR) or human epidermal growth factor receptor 2 (HER2). 1 Since TNBC lacks a therapeutic target, patients with this type of breast cancer do not benefit from molecularly targeted treatments such as endocrine therapy or trastuzumab. 2 Recently, a growing body of evidence suggests that these patients are potential candidates for immunotherapy. 3
The use of immunotherapy has markedly altered the field of anti-cancer treatment. 4 Tumours with a high mutation burden are known to have considerable neoantigens and tumour-infiltrating lymphocytes (TILs). 5 Tumour mutation burden is mainly analysed by assessing mutation rates using DNA sequencing data. However, RNA editing is another source of alterations in amino acid sequences and epitopes. The frequency of RNA editing in many types of tumours has been explored using The Cancer Genome Atlas (TCGA) data, and a substantial amount of mutational load due to RNA editing has been detected. 6 The adenosine deaminase acting on RNA (ADAR) family of editases is a group of enzymes responsible for converting an adenosine into an inosine (A-to-I) in structured or double-stranded (ds) RNAs; this is the most prevalent form of RNA editing in mammals. Cellular functions and regulations can be altered by A-to-I editing through changes in messenger RNA (mRNA) splicing, stability, localisation and translation as well as interference with the binding of regulatory RNAs. 7 With the improvement in sequencing and analysis technologies, researchers have studied RNA editing events using RNA sequencing data.6–9
In humans, the ADAR family comprises ADAR1, ADAR2 and ADAR3, which share repeated copies of dsRNA-binding domains and a catalytic adenosine deaminase domain.
10
Among them, ADAR1 is the major contributor to RNA editing in cancer, especially in breast. Anantharaman et al.
11
have demonstrated that breast cancer cells show elevated
Materials and methods
Patient and tissue specimens
We evaluated 681 TNBC patients whose resected tissues were obtained by surgeries conducted between 2004 and 2010 at Asan Medical Center. In this group, 472 of the 681 patients did not present with lymph node metastasis and received four cycles of adjuvant anthracycline and cyclophosphamide (AC: adriamycin, 60 mg/m2 and cyclophosphamide, 600 mg/m2). The remaining 209 patients presented with lymph node metastasis and were treated with four cycles of AC, followed by either four cycles of paclitaxel (175 mg/m2) or four cycles of docetaxel (75 mg/m2). In total, 550 patients (80.8%) received radiotherapy. The median follow-up period was 78.3 months. All patients were chemotherapy-naïve and radiotherapy-naïve preoperatively. Clinical information was obtained from medical records and pathology reports from the surgery. Exemption from informed consent after the de-identification of information was approved by the Institutional Review Board of Asan Medical Center.
Histological evaluation
Haematoxylin and eosin–stained slides were reviewed by two pathologists (H.J.L. and G.G.) and were histopathologically analysed for TILs (defined as the percentage of stroma of invasive carcinoma that is infiltrated by lymphocytes in 10% increments), 14 histological grade, tumour size, pT stage, pN stage and lymphovascular invasion. The histological type was defined based on the 2012 WHO classification criteria, and the histological grade was assessed using the modified Bloom–Richardson classification. 15
Tissue microarray construction and immunohistochemistry
Formalin-fixed, paraffin-embedded tissue samples were evaluated using a tissue-arraying instrument, as previously described. 16 Each sample was arrayed in three 1-mm-diameter cores to minimise tissue loss and mitigate tumour heterogeneity effects. Tissue microarray sections were stained using an automatic immunohistochemical staining device (Benchmark XT; Ventana Medical Systems, Tucson, AZ, USA). Antibodies to target ADAR1 (1:200 dilution; HPA003890; Atlas Antibodies, Stockholm, Sweden), MxA (1:1000; ab95926; Abcam, Cambridge, UK), double-stranded RNA (dsRNA), activated protein kinase (PKR, 1:400; ab32052; Abcam), CD8 (1:200; DAKO, Glostrup, Denmark) and HLA-ABC (1:1600; ab70328, EMR8-5; Abcam) were used. 12 As controls for ADAR1, cell blocks from human breast cancer cell lines were used (Supplementary Figure S1). 13
ADAR1, MxA, PKR and HLA-ABC intensities were evaluated as a four-value intensity score (0, 1, 2 and 3). The percentages of nuclear ADAR expression, membranous and/or cytoplasmic expression of HLA-ABC and cytoplasmic MxA and PKR expression were also evaluated. An ‘immunoreactive score’ was generated as the product of the intensity and percentage of positive cells. Immunoreactive scores for ADAR1, MxA and PKR were dichotomised based on the mean protein expression value. We categorised HLA-ABC expression in tumour cells as one of the three levels, as previously described (strongly positive, expression in >75% of tumour cells; weakly positive, expression in 25%–75% of tumour cells and negative, loss of expression in >75% of tumour cells). 17
The immunostained tissue microarray slides for CD8 were scanned using a digital microscopy scanner (Pannoramic 250 FLASH; 3DHISTECH Ltd., Budapest, Hungary). The numbers of CD8 + cells in tissue microarray cores were counted using the NuclearQuant module of Pannoramic Viewer 1.15.2 (3DHISTECH Ltd.).
TCGA and Cancer Cell Line Encyclopedia data
TCGA data portal was used to download the breast cancer data. 18 In total, 84 cases of basal-like breast cancer were analysed according to the PAM50 predictor. We downloaded level-3 gene expression data derived from an Agilent custom 244K whole-genome microarray and somatic mutation data obtained using exome sequencing from Illumina Genomic Analyzer. Gene expression levels were median-centred by genes.
Gene Set Enrichment Analysis (GSEA) was performed using javaGSEA v2.2.4, with gene sets derived from the Kyoto Encyclopedia of Genes and Genomes pathway database.
19
Basal-like tumours were dichotomised based on the mean
Statistical analysis
All statistical analyses were conducted using SPSS statistical software (version 18; SPSS, Chicago, IL, USA) or R (https://www.R-project.org/). The chi-square test, log-rank test, Spearman’s correlation and forward conditional Cox proportional hazards regression model were used to evaluate the data. All tests were two-sided, and
Results
ADAR1 expression in TNBCs
ADAR1 was expressed in tumour cell nuclei (Figure 1 and Supplementary Figure S1). The breast cell lines with relatively high

Representative figures of ADAR1 expression in breast cancer. (a) Tumour cells show nuclear positivity for ADAR1. (b) Tumour cells are negative for ADRA1 (immunohistochemistry, 400×).
Clinicopathological characteristics of patients with breast cancers according to ADAR1 expression levels.
HLA: human leukocyte antigen; LN: lymph node; TIL: tumour-infiltrating lymphocyte.
Bold values denotes the significance level of 5% (p < 0.05).
Prognostic significance of ADAR1 expression in TNBCs
Because standardised adjuvant chemotherapy differs between patients having tumours with lymph node metastasis and those without lymph node metastasis, we analysed the prognostic significance of ADAR1 expression according to the presence or absence of lymph node metastasis. In the patient group without lymph node metastasis who were treated with adjuvant anthracycline-based chemotherapy, ADAR1 expression exhibited no prognostic significance in disease-free survival (Figure 2(a),

Kaplan–Meier disease-free survival curves according to ADAR1 expression. (a) In the patient group without lymph node metastasis, ADAR1 expression exhibited no prognostic significance in disease-free survival. (b) In the patient group with lymph node metastasis, patients with high ADAR1 expression demonstrated worse disease-free survival, but this was not statistically significant. (c) In the patient group without lymph node metastasis, patients with a high TIL levels tended to demonstrate favourable clinical outcomes, regardless of ADAR1 expression. (d) In the patient group with lymph node metastasis, patients with low TIL levels and high ADAR1 expression demonstrated worst outcomes.
Cox proportional hazard regression model of disease-free survival in the patients with triple-negative breast cancers with lymph node metastasis based on clinicopathological variables.
CI: confidence interval; HR: hazard ratio; LVI: lymphovascular invasion; TIL: tumour-infiltrating lymphocyte.
Bold values denotes the significance level of 5% (
We further analysed the prognostic significance according to the combination of TIL levels and ADAR1 expression. In the patient group without lymph node metastasis, patients with high TIL levels tended to exhibit better clinical outcomes than those with low TIL levels, regardless of ADAR1 expression (Figure 2(c),
Analysis of TCGA data
To explore the relationships among DNA mutation load, RNA editing and TIL levels, we analysed the total number of DNA mutations and

Correlations of
To analyse the biological characteristics of tumours with high
Differentially regulated pathways in basal-like tumours with high
NES: normalised enrichment score; FDR: false discovery rate; RIG-I: retinoic acid-inducible gene I; JAK/STAT: Janus kinase/signal transducer and activator of transcription.
Discussion
The main target of A-to-I editing is
Although ADAR1 expression was positively correlated with TIL levels, patients with high ADAR1 expression demonstrated a worse disease-free survival rate than those with low ADAR1 expression. We have previously shown that high TIL levels and high MxA expression in TNBCs were associated with favourable clinical outcomes.12,22 Conversely, Fumagalli et al.
7
have demonstrated that ADAR silencing in breast cancer cell lines (MDA-MB-231, MCF-7 and BT-474) is associated with a statistically significant decrease in cell proliferation and increase in cell apoptosis, suggesting the role of ADAR in cell proliferation and apoptosis. We performed a GSEA study and demonstrated that tumours with high
A-to-I editing can be an important parameter for cancer treatment strategies because its levels are deregulated during cancer and can promote cancer progression. 21 Analysis using TCGA data demonstrated that RNA editing is abundant in many types of human cancers. Several reports have also revealed that RNA editing rates are higher in tumour samples of breast, thyroid and lung cancer than in normal tissue. However, in chromophobe renal cell carcinoma, renal papillary cell carcinoma and metastatic melanoma, A-to-I editing levels are reduced in tumour tissues.6,10 The downregulation of ADAR1 expression in metastatic melanoma promotes tumour growth and cell proliferation. 10 In oesophageal squamous cell carcinoma, ADAR1 has been identified as an oncogene and negative predictor of overall survival. 10 Therefore, studies exploring the possibility of the modulation of ADAR and downstream molecular activity are warranted to develop novel therapeutic modalities.
The discovery of neoantigens is an important research topic for developing personalised immunotherapies using a vaccine or adoptive antigen-specific T-cell therapy.
23
In general, tumour mutational load is positively correlated with the number of neoantigens, TIL levels and T-cell diversity.24,25 We analysed TCGA data of basal–like breast cancers defined by PAM50 subtype predictors, because basal-like breast cancers are highly concordant with TNBCs defined by immunohistochemistry (concordance rate: ~75%), and basal-like breast cancers that are not triple negative show similar genomic features to those of basal-like TNBCs.
26
However, in the basal-like breast cancers assessed using TCGA data, we could not identify a positive correlation between
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (2015R1C1A1A02036484 and 2016R1D1A1B03935367) and the Korean Health Technology R&D Project, Ministry of Health & Welfare (HI17C0337).
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
