Expression of proteins related to autotaxin–lysophosphatidate signaling in thyroid tumors

Background We aimed to investigate the expression of proteins related with autotaxin (ATX)–lysophosphatidate (LPA) signaling and the clinical implications in primary and metastatic thyroid tumors. Methods We constructed tissue microarrays with 545 primary thyroid tumors [338 papillary thyroid carcinoma (PTC), 111 follicular carcinoma (FC), 69 medullary carcinoma (MC), 23 poorly differentiated carcinoma (PDC), and four anaplastic carcinoma (AC)]. Immunohistochemical stains for proteins related to ATX–LPA signaling (e.g., ATX, LPA1, LPA2, and LPA3) were performed. Results The expression of ATX was highest in MC, while the LPA1 expression was higher in PDC and AC, and the expression of LPA2 and LPA3 was highest in PTC (p < 0.001). Additionally, the expression of ATX, LPA1, and LPA2 was higher in conventional-type PTC than in follicular-variant PTC (p < 0.05). PTC with BRAF V600E mutation showed higher expression of ATX, LPA1, LPA2, and LPA3 than PTC without BRAF V600E mutation (p < 0.001). In univariate analysis, ATX positivity (p = 0.005) and LPA1 positivity (p = 0.014) were correlated with shorter overall survival in PTC. Conclusion Proteins related to the ATX–LPA axis showed different levels of expression in primary thyroid tumors according to subtype.

Thyroid cancer is a relatively common malignancy, affecting about 1% of the population. The most common histologic type of thyroid cancer is papillary thyroid carcinoma (PTC), followed by follicular carcinoma (FC), medullary carcinoma (MC), poorly differentiated carcinoma (PDC), and anaplastic carcinoma (AC). These histologic subtypes of thyroid cancer are known to have different cell origins, clinical manifestations, metastatic patterns, and clinical prognoses. Previously, ATX has been reported to be highly expressed in undifferentiated carcinoma [7], and it is also indicated to be expressed in follicular carcinoma [8]. Furthermore, as ATX expression is higher in thyroid cancer than in benign lesion [9] and is correlated with prognostic factors in thyroid cancer [10], ATX-LPA signaling is likely to have an influence on thyroid cancer biology. However, studies on proteins related to ATX-LPA signaling in recurrent and metastatic thyroid cancer and according to thyroid tumor subtype have Page 2 of 9 Shin and Koo J Transl Med (2019) 17:288 not been reported. Thus, we aimed to assess the expression of proteins related to ATX-LPA signaling in recurrent and metastatic thyroid cancer and according to thyroid tumor subtype.

Patient selection and histologic evaluation
Patients diagnosed with PTC who underwent surgery between January 2012 and December 2013 at Severance Hospital, Sinchon-dong, South Korea, were included in the study. As for the remaining subtypes, patients who underwent surgery between January 2000 and December 2014 at Severance Hospital were included. Those who had received preoperative chemotherapy were excluded. This study was approved by the Institutional Review Board of Yonsei University Severance Hospital. Hematoxylin-and-eosin-stained slides of all cases were retrospectively reviewed by a thyroid pathologist (Koo, JS). Clinicopathologic data were obtained from the patients' medical records and included age at diagnosis, disease recurrence, metastasis, current status, and length of follow-up. Tumor size, location (right or left lobe), extent (confined to the thyroid parenchyme or with extrathyroidal spread), and number of metastatic lymph nodes were also noted from a review of the slides and the surgical pathology reports. Tumor stroma of PTC was classified as follows: desmoplastic type, tumor stroma composed of cellular fibroblast/myofibroblast proliferation; sclerotic type, tumor stroma composed of collagenous component with few cellular components; pauci type, near absence of a stromal reaction; or inflammatory type, tumor stroma composed of inflammatory cells such as lymphocytes.

Tissue microarray
A representative area showing tumor and tumor stroma was selected on a hematoxylin-and-eosin-stained slide, and a corresponding spot was marked on the surface of the paraffin block. Using a biopsy needle, the selected area was punched out, and a 3-mm tissue core was transferred to a 6 × 5 recipient block. Two tissue cores from the invasive tumors were extracted to minimize extraction bias. Each tissue core was assigned a unique tissue microarray location number that was linked to a database containing other clinicopathologic data.

Immunohistochemistry
The antibodies used for immunohistochemistry in this study are shown in Table 1. Formalin-fixed, paraffinembedded tissue sections were used for immunohistochemistry. During the experiments, 3-μm-thick tissue sections were deparaffinized and rehydrated in xylene and graded alcohol. We used the Ventana Discovery XT automated stainer (Ventana Medical Systems, Tucson, AZ, USA). Cell Conditioning 1 buffer (citrate buffer, pH: 6.0; Ventana Medical Systems, Tucson, AZ, USA) was used for antigen retrieval. Appropriate positive and negative controls were used for each antibody.

Interpretation of immunohistochemical staining
Immunohistochemical markers were assessed by light microscopy. All stained slides were semi-quantitatively evaluated [11]. Tumor and stromal cell staining was assessed as 0: negative or weak immunostaining in < 1% of the tumor/stroma; 1: focal expression in 1% to 10% of the tumor/stroma; 2: positive in 11% to 50% of the tumor/ stroma; or 3: positive in 51% to 100% of the tumor/ stroma. Entire tumor areas included in the TMA were evaluated in all cases; a grade of 0 was considered negative and grades of 1 to 3 were considered positive.

Statistical analysis
Data were analyzed using Statistical Table S3). Basal characteristics of MC, PDC, and AC are listed in Additional file 1: Table S3.

Expression of proteins related to the ATX-LPA axis in thyroid cancer
The expression of proteins related to the ATX-LPA axis in thyroid cancer showed differences in ATX, LPA1, LPA2, and LPA3 (p < 0.001). ATX was most highly expressed in medullary carcinoma, while LPA1 was predominant in PDC and AC, and LPA2 and LPA3 were predominant in PTC (Table 1 and Fig. 1).
When assessed according to the histologic subtype and BRAF V600E mutation status, ATX (p = 0.002), LPA1 (p = 0.014), and LPA2 (p = 0.001) expression was different according to the histologic subtype, and the expression of ATX, LPA1, LPA2 and LPA3 (p < 0.001) was different according to BRAF V600E mutation status. ATX, LPA1, and LPA2 expression was higher in conventional type than in FVPTC, while PTC with BRAF V600E mutation showed higher expression of ATX, LPA1, LPA2, and LPA3 than did PTC without BRAF V600E mutation ( Table 2 and Fig. 2).

Correlations between clinicopathologic factors and expression of proteins related to the ATX-LPA axis in thyroid cancer
When correlations between the expression of proteins related to the ATX-LPA axis and clinicopathologic factors were analyzed, LPA2 expression was different according to the stromal type (p = 0.009) in PTC, and LPA1 expression was different according to the vascular invasion status in FC (p < 0.001). In PTC, desmoplastic type showed high expression of LPA2. LPA1 negativity was correlated with vascular invasion in FC (Fig. 3).

Impact of the expression of proteins related to the ATX-LPA axis on prognosis of thyroid cancer
The impact of expression of proteins related to the ATX-LPA axis on the prognosis of thyroid cancer was analyzed. As a result, ATX positivity (p = 0.005) and LPA1 positivity (p = 0.014) were two factors that demonstrated correlation with shorter overall survival (Table 5 and Fig. 4) in univariate analysis. Multivariate Cox analysis showed that lymph node metastasis (hazard ratio: 5.978, 95% confidence interval (CI) 1.339-26.69; p = 0.019) was an independent factor associated with shorter DFS, while age (≥ 45 years) (hazard ratio: 16.39, 95% CI 2.129-126.2; p = 0.007) was an independent factor associated with shorter OS (Table 6).

Discussion
We investigated the expression of proteins related to the ATX-LPA axis in thyroid cancer. First, ATX was more highly expressed in medullary carcinoma than in other subtypes. Previous studies have shown that ATXrelated protein expression was higher in thyroid cancer or metastatic thyroid cancer than in normal thyroid tissue or benign thyroid neoplasm [9,10], but there has not been any report published to date on ATX-related protein expression in thyroid medullary carcinoma.  The possible mechanism for the high expression of ATX in thyroid medullary carcinoma is the JAK/ STAT3 pathway. In previous studies on breast cancer, activation of the JAK/STAT3 pathway was reported to increase the expression of ATX [12], which is suggested to represent a putative STAT3 target gene. Meanwhile, it has been reported that the JAK/STAT3 pathway is also activated in thyroid medullary carcinoma [13,14]; consequently, it is only plausible that ATX expression is somehow associated with the JAK/STAT3 pathway. Our study showed that LPA1 was highly expressed in PDC and AC, and that LPA2 and LPA3 were similarly highly expressed in PTC. A previous study on thyroid cancer reported that LPA1 messenger RNA expression showed no difference among normal thyroid tissue, benign thyroid nodule, and thyroid cancer [15]. However, the study had included PTC and FC for thyroid cancer; therefore, its direct comparison with our study, which includes PDC and AC in addition to PTC and FC for thyroid cancer, is not possible. Separately, another previous study suggested CD97 as a dedifferentiation marker in thyroid cancer [16], the expression of which is reported to be correlated with LPA receptor in thyroid cancer [17]. Such a result is related with the high expression of LPA1 in PDC and AC in our study. LPA2 Fig. 2 Expression of proteins related to ATX-LPA axis in TPC. ATX, LPA1, LPA2, and LPA3 expression is higher in conventional-type PTC than in FVPTC, and PTC with BRAF V600E mutation shows higher expression of those than does PTC without BRAF V600E mutation   17:288 has been previously reported to be highly expressed in PTC, which is in concordance with our results [15]. On the other hand, our study showed low expression of proteins related to the ATX-LPA axis in PTC and FVPTC without BRAF V600E mutation, a result that disagreed with those of the previous study that reported that ATX protein expression is not correlated with BRAF mutation status [9]. Usually, most FVPTC cases do not harbor BRAF V600E mutation; thus, it is necessary to assess the correlation between BRAF V600E mutation and the ATX-LPA axis. The most plausible mechanism is that BRAF V600E mutation activates cytokine release, which in turn increases ATX protein expression. Tumor cells showing BRAF V600E mutation tend to have increased cytokine secretion such as that of interleukin (IL)-1β, IL-6, and IL-8 [18], because cytokines such as tumor necrosis factor alpha or IL-1β increase ATX and LPA levels. In our study, ATX and LPA1 expression was a poor prognosis factor in PTC, which is in concordance with previous study results that had identified ATX as a poor prognosis factor in breast cancer [19] and in prostate cancer [20] and, separately, LPA2 as a poor prognosis factor in breast cancer [21].  The clinical implication of our study is that the ATX-LPA signaling axis may be a possible therapeutic target in thyroid cancer. Ki16425, a nonlipid competitive inhibitor of LPA1 and LPA3, actually suppressed bone metastasis in a mouse model [22]. Debio0719, an R-stereoisomer of Ki16425, was found to inhibit distant metastasis [23,24]. BrP-LPA, which is a dual ATX and pan-LPAR inhibitor, inhibited cell migration and invasion [25]. In addition, an ATX inhibitor, ONO-843050, decreased tumor volume to 50% to 60% in a mouse model of PTC [9]. Therefore, the ATX-LPA axis may be an effective treatment target for thyroid cancer, and further clinical trials are necessary.