Molecular characterization of papillary thyroid carcinoma: a potential three-lncRNA prognostic signature

Abstract

Purpose

Papillary thyroid carcinoma (PTC), the most frequent type of malignant thyroid tumor, lacks novel and reliable biomarkers of patients’ prognosis. In the current study, we mined The Cancer Genome Atlas (TCGA) to develop lncRNA signature of PTC.

Patients and methods

The intersection of PTC lncRNAs was obtained from the TCGA database using integrative computational method. By the univariate and multivariate Cox analysis, key lncRNAs were identified to construct the prognostic model. Then, all patients were divided into the high-risk group and low-risk group to perform the Kaplan–Meier (K–M) survival curves and time-dependent receiver operating characteristic (ROC) curve, estimating the prognostic power of the prognostic model. Functional enrichment analysis was also performed. Finally, we verified the results of the TCGA analysis by the Gene Expression Omnibus (GEO) databases and quantitative real-time PCR (qRT-PCR).

Results

After the comprehensive analysis, a three-lncRNA signature (PRSS3P2, KRTAP5-AS1 and PWAR5) was obtained. Interestingly, patients with low-risk scores tended to gain obviously longer survival time, and the area under the time-dependent ROC curve was 0.739. Furthermore, gene ontology (GO) and pathway analysis revealed the tumorigenic and prognostic function of the three lncRNAs. We also found three potential transcription factors to help understand the mechanisms of the PTC-specific lncRNAs. Finally, the GEO databases and qRT-PCR validation were consistent with our TCGA bioinformatics results.

Conclusion

We built a three-lncRNA signature by mining the TCGA database, which could effectively predict the prognosis of PTC.

Keywords:

• The Cancer Genome Atlas
• long noncoding RNAs
• overall survival
• risk score

Introduction

Thyroid cancer (TC) is the most common malignant tumor of the thyroid, accounting for about 3%–4% of newly diagnosed tumors annually.¹ TC includes four pathological types: papillary thyroid carcinoma (PTC), follicular thyroid carcinoma, anaplastic thyroid carcinoma and medullary thyroid carcinoma.² PTC represents about 80% of TC with low malignancy and better prognosis,^3,4 and its incidence is increasing year by year.⁵ PTC develops slowly and can be cured by surgery, thyroid hormone therapy and ¹³¹I isotope therapy.⁶ Most PTC patients have good prognosis with the 10-year survival rate reaching 90%.^7,8 Because of the serious potential side effects of radioactive iodine and patients’ neglect of the PTC prognosis, some PTC patients suffer from recurrence and some PTCs may develop distant metastases with high mortality. Therefore, to improve the long-term survival time, prevent PTC deterioration and reduce the recurrence incidence, the novel and reliable biomarkers are important for proper treatment and the prognosis of PTC patients.

Long noncoding RNAs (lncRNAs) are conservative with more than 200 nucleotides in length, having no significant protein-coding capacity.⁹ Nowadays, lncRNAs have become a hot topic in genetic research and plays an important role in many life activities, such as the dosage compensation effect, epigenetic regulation, cell cycle regulation and cell differentiation regulation. Moreover, increasing evidence indicates that lncRNAs play a potential role as novel biomarkers for prognostic prediction in various cancers.^10,11 It was found that a number of lncRNAs were correlated with the prognosis of PTC. Liao et al¹² reported that BANCR could inhibit tumorigenesis in PTC by performing quantitative real-time PCR (qRT-PCR), cell function and tumor xenograft in nude mice, suggesting that BANCR may serve as a novel prognostic marker. Lan et al¹³ screened the genome-wide PTC to understand the expression profile and potential functions of lncRNAs. Liyanarachchi et al¹⁴ observed the two lncRNAs associated with lymph node metastasis and BRAF mutation in PTC using The Cancer Genome Atlas (TCGA). Another study performed by Zhao et al¹⁵ described a PTC competing endogenous RNA (ceRNA) network and found three survival-related lncRNAs in TCGA. Thus, the TCGA dataset provides a possibility for us to find prognostic biomarkers in PTC.

In the current study, the expression of lncRNAs and related information were obtained from the TCGA database. Then, we built an lncRNA signature for predicting the PTC patients’ survival time to provide a potential biomarker of PTC patients’ prognosis.

Patients and methods

TCGA dataset and sample information

Up to October 23, 2017, 503 PTC cases and 59 normal thyroid cases were obtained from the TCGA database. Inclusion criteria were as follows: 1) patients who had a pathological diagnosis of PTC and 2) patients whose expression of genes and clinical feature were available. Exclusive criteria were as follows: 1) patients who suffered from other cancers except PTC and 2) patients who did not die of PTC. Overall, this study embraced 485 eligible samples, including 430 PTC tumor cases and 55 normal thyroid cases. On the basis of the seventh American Joint Committee on Cancer (AJCC) TNM staging system, 243, 46, 89 and 52 PTC patients were in the stages I, II, III and IV, respectively. In addition, there were 214 PTC cases with lymphatic metastasis and 216 PTC cases without lymphatic metastasis.

In addition, 32 paired frozen PTC samples were obtained from Nanjing Drum Tower Hospital, including tumor tissues and adjacent normal tissues. These tissues were stored in RNAlater (GenStar BioSolutions, Beijing, People’s Republic of China) at −80°C immediately after exairesis until further use. Pathology reports and quality assessment reports were required to verify the collected samples. Informed consent agreements were signed by all patients in this study. Our study was approved by the ethics committee of Nanjing Drum Tower Hospital.

Differential analysis of expressed lncRNAs

lncRNA expression profiles (level 3) of PTC patients were obtained from the TCGA data portal, and their raw data were processed and normalized by TCGA RNASeqv2 system. Fold change >2 or <0.5, P-value <0.05 and false discovery rate (FDR) <0.05¹⁶ were setup to identify significantly differentially expressed lncRNAs. Then, all samples collected in TCGA database were divided into three groups to select the intersection of lncRNAs, including PTC tumor samples and normal thyroid samples, PTC samples with lymphatic metastasis and without lymphatic metastasis, PTC samples with stage I–II and stage III–IV. Then, the intersection of lncRNAs was obtained for further analysis. Figure S1 shows the flowchart of bioinformatic analysis.

Construction of potential prognostic signatures

To find the prognostic lncRNAs, we performed Cox regression analysis (univariate and multivariate analysis). Then, the risk score model was constructed by corresponding coefficients of the prognostic lncRNAs: risk score = Exp_lncRNA1* β_lncRNA1 + Exp_lncRNA2 * β_lncRNA2 + …+ Exp_lncRNAn * β_lncRNAn, where Exp is the expression level and β is the regression coefficient derived from the multivariate Cox regression model.

The lower the score, the lower the risk of death outcomes. In addition, all samples were divided into the low-risk score group and the high-risk score group according to the median scores. To identify the distinguishing ability of patients’ outcome, the Kaplan–Meier (K–M) survival curves were performed, and the time-dependent receiver operating characteristic (survival ROC)¹⁷ curves were employed to detect the prognostic power of the risk score model. According to all eligible patients’ clinical information, the relationship between the survival time and clinical factor was evaluated using the survival analysis.

The lncRNA signature function

To investigate the function of the risk score model, we first identified the genes that were highly related to the lncRNA signature (Pearson |R|>0.6) in the TCGA. Then, the gene ontology (GO) biological processes and Kyoto Encyclopedia of Genes and Genomes pathways were analyzed using Database for Annotation, Visualization and Integrated Discovery, and we were only interested in the significant level (P-value <0.05 and FDR <0.05). Moreover, the protein–protein interaction (PPI) network was performed to analyze the inter-function among co-related genes using GENets (https://string-db.org/).¹⁸

Screening of potential transcription factors

In the current study, to identify the transcription factors of the lncRNA-related genes, FunRich, an analysis tool of transcription factors, was used (http://www.funrich.org).¹⁹

The lncRNA signature verification

To confirm the reliability and validity of the risk score model, we first used other cohorts from Gene Expression Omnibus (GEO) databases. The lncRNA data of PTC (GSE35570, GSE60542 and GSE33630) was obtained from the GEO database, which was based on the Platform GPL570.

Second, we used qRT-PCR to examine their actual expression in 32 pairs of PTC samples. Their total RNA was extracted using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s protocol. After the purity was detected, reverse transcription reactions and RT-PCR were performed as per our previous studies.^20,21 Reverse transcription reactions using StarScript II First-strand cDNA Synthesis Mix (GenStar BioSolutions) were conducted in two steps according to its specification. qRT-PCR was carried out to detect the expression levels of candidate lncRNAs with the Step One Plus™ PCR System (Applied Biosystems) using RealStar Green Fast Mixture (GenStar BioSolutions) according to the manufacturer’s protocol. The primer sequences were as follows: GAPDH: forward 5-GGGAGCCAAAAGGGTCATCA-3, reverse 5-TGATG-GCATGGACTGTGGTC-3, product size 203 bp; human PRSS3P2: forward 5-TGGGGTGGCGGATGATCTTGG-3, reverse 5-ACAGTGGGTGGTGTCAGCAGG-3, product size 128 bp; human KRTAP5-AS1: forward 5-GCCCT-GCCTTCAGCCTCCTCA-3, reverse 5-GAAACGTCGGCCTGGCCTTGT-3, product size 186 bp; human PWAR5: forward 5-GAGACTCAAAGGCAAGAACTA-3, reverse 5-TGATGTGGGTGTTGATACTGT-3, product size 462 bp. Results were normalized to the expression of GAPDH. The results of qRT-PCR were calculated by the 2^−∆∆Ct method: ∆Ct = Ct_lncRNAs – Ct_GAPDH and ∆∆Ct = ∆Ct_{tumor tissues} – ∆Ct_{adjacent non-tumor tissues}.

Statistical analysis

All data are expressed as mean ± SD. The threshold of FDR was set as 0.05 to reduce the false positives in multiple tests in bioinformatics analysis. Moreover, the threshold of P-value was set as 0.05 to evaluate the null hypothesis. The key lncRNAs were selected by univariate/multivariate Cox regression model to build the risk score model, and the predicted power was identified by K–M survival curves and survival ROC. The Student’s t-test and the paired t-test were performed to calculate the results of GEO datasets and qRT-PCR. The abovementioned analyses were conducted using SPSS 21.0 (IBM Corporation, Armonk, NY, USA).

Results

Patient characteristics

According to the inclusion and exclusion criteria, 430 patients were included in the current study. These PTC patients were 47.682±15.574 years old. Among them, there were 111 male (25.8%) and 319 (74.2%) female patients. The overall survival (OS) time was 1,195.479±945.363 days, and eleven of 430 (2.6%) PTC patients had died.

Identification of key survival-related lncRNAs

As shown in the flowchart of Figure S1, 83 key lncRNAs were selected through the comparison between the four groups (tissues in stages I–II with non-lymphatic metastases, stages III–IV with non-lymphatic metastases, stages I–II with lymph node metastasis and stages III–IV with lymph node metastasis vs adjacent non-tumorous tissues; Figure 1 and Table S1).

Figure 1 Cluster analysis of consistent differential lncRNA expression.

Note: Red indicates that the lncRNA has higher expression level; green indicates that the lncRNA has lower expression.

According to the correlation between 83 key lncRNAs and clinical features of 430 PTC patients, nine survival-related lncRNAs were obtained with univariate Cox regression model (Table 1). To improve the robustness of the candidate lncRNAs, three lncRNAs (PRSS3P2, KRTAP5-AS1 and PWAR5) were identified for further analysis after performing multivariate Cox regression analysis, revealing that the three lncRNAs were independent prognostic signatures (Table 2).

Table 1 Prognostic value of the differentially expressed lncRNAs by univariate Cox regression analysis

lncRNA	Estimate	Chi-squared	P-value	HR (95% CI)
PRSS3P2	–1.662	5.641	0.018	0.190 (0.041–0.879)
CYP1B1-AS1	–1.485	4.308	0.038	0.226 (0.049–1.050)
LINC00602	1.614	5.254	0.022	5.022 (1.084–23.273)
SNORD116-4	1.672	5.724	0.017	5.324 (1.148–24.678)
RPSAP52	–1.596	5.123	0.024	0.203 (0.044–0.939)
DYP19L2P4	1.448	5.317	0.021	4.225 (1.118–16.185)
GGT3P	–1.511	4.476	0.034	0.221 (0.048–1.024)
KRTAP5-AS1	–1.490	4.351	0.037	0.225 (0.049–1.044)
PWAR5	1.628	5.373	0.020	5.096 (1.100–23.611)

Table 2 Prognostic value of the differentially expressed lncRNAs by multivariate Cox regression analysis

lncRNA	β	HR	HR 95% CI	P-value
PRSS3P2	–1.620	0.198	0.042–0.924	0.039a
KRTAP5-AS1	–1.742	0.175	0.038–0.814	0.026a
PWAR5	1.707	5.514	1.184–25.681	0.030a

Note:

a P<0.05.

Weighted by their corresponding coefficients, a three-lncRNA prognostic signature was constructed: risk score = Exp_PRSS3P2 * (−1.620) + Exp_KRTAP5-AS1 * (−1.742) + Exp_PWAR5* 1.707.

As shown in the risk score model, the coefficient of PRSS3P2 and KRTAP5-AS1 was negative, indicating that they were positively related to the OS of PTC patients while the expression of PWAR5 was a negative factor.

To evaluate the predicted power of the risk model, the patients were divided into two groups (high- and low-risk score groups) using the median risk score value as the cutoff (Figure 2). The K–M analysis confirmed that the OS time of the low-risk score patients was 1,216.495±934.724 days and the high-risk score patients was 1,174.657±959.669 days, indicating that the risk score was significantly associated with OS (P=0.022; Figure 3A). Furthermore, the three-lncRNA signature exhibited the strong predicted power of 5-year survival of the PTC patients (the area under the ROC curve was 0.739; Figure 3B), suggesting that the risk score model performed good sensitivity and specificity in predicting the PTC patients’ survival.

Figure 2 Risk score analysis of the differentially expressed lncRNA signature of PTC.

Note: Survival status and duration of cases and risk score of lncRNA signature.

Figure 3 The three differentially expressed lncRNA signatures of PTC for the outcome.

Notes: (A) The risk score is shown by the time-dependent ROC curve for predicting 5-year survival. (B) The K–M test of the risk score for the OS.

Abbreviations: AUC, area under the curve; K–M, Kaplan–Meier; OS, overall survival; PTC, papillary thyroid carcinoma; ROC, receiver operating characteristic.

Correlation between clinical characteristics and survival time

To examine the association between the three-lncRNA signature with clinical features and survival time in PTC patients, survival analysis was performed. The results indicated that tumor stage, residual tumor, age and neoplasm cancer status could predict the poorer survival of PTC patients (Figure S2). In addition, they had prognostic values, including tumor stage (area under the curve [AUC] =0.766, P=0.003), residual tumor (AUC =0.703, P=0.029), age (AUC =0.836, P=0.000) and neoplasm cancer status (AUC =0.807, P=0.010; Figure S2).

Functional assessment of the three-lncRNA signature

A total of 1,098 genes identified in TCGA database co-expressed with lncRNAs of the risk score model (PRSS3P2, KRTAP5-AS1 and PWAR5; |R|>0.6), including 411 genes with PRSS3P2, 90 genes with KRTAP5-AS1 and 597 genes with PWAR5, respectively. A total of 367 GO terms and 73 pathways (P-value <0.05 and an enrichment score of >1.5) were identi fied. The results showed that the top GO biological process was regulation of transcription, DNA dependent (GO:0006351) and gene expression (GO:0010467; Figure 4A), and the results suggested that these genes were mainly enriched in ribosome and metabolic pathways (Figure 4B). The relationship and function of these genes were revealed in the PPI network (Figure 5).

Figure 4 Enrichment of KEGG pathways and GO terms for co-expressed mRNAs, and the map represents the PPI network of co-expressed genes.

Note: (A and B) KEGG and GO analysis of the related genes.

Abbreviations: GO, gene ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; PPI, protein–protein interaction; NAFLD, nonalcoholic fatty liver disease.

Figure 5 PPI network of mRNAs.

Abbreviation: PPI, protein–protein interaction.

Screening of potential transcription factors

Based on the data from FunRich, the transcription factors for lncRNA co-expressed genes were ELK1, ETV7 and GABPA, which may be helpful to understand the mechanisms of the co-expressed genes of specific lncRNAs (Figure 6).

Figure 6 Enriched transcription factors by lncRNA co-expressed genes.

GEO dataset and qRT-PCR verification

The expression level of these three lncRNAs is summarized in Table 3. The expression level of PRASS3P2 and PWAR5 was evidently different between the PTC tissues and the normal tissues as shown in GSE35570, GSE60542 and GSE33630. The KRTAP5-AS1 expression in PTC was significantly higher than normal tissues as shown in GSE35570 and GSE33630. The results of GEO dataset were consistent with the abovementioned TCGA results. Unfortunately, no survival information of PTC lncRNA was available in the GEO datasets.

Table 3 Expression data of three lncRNAs in GEO datasets

GEO datasets	Year	Country	Platform	Sample	n	PRSS3P2	P-value	KRTAP5-AS1	P-value	PWAR5	P-value
GSE35570	2015	Poland	GPL570	PTC	65	4.10±1.36	<0.01a	2.45±0.24	0.022a	5.59±1.26	<0.01a
				Normal	51	2.82±0.25		2.36±0.08		7.03±1.19
GSE60542	2015	Belgium	GPL570	PTC	33	6.94±0.86	<0.01a	4.43±0.16	0.556	4.94±0.80	0.002a
				Normal	30	5.96±0.17		4.40±0.14		5.96±1.05
GSE33630	2012	Belgium	GPL570	PTC	47	6.52±1.08	<0.01a	5.49±0.21	<0.01a	4.62±0.74	<0.01a
				Normal	47	5.50±0.35		5.33±0.17		5.49±0.85

Note:

a lncRNAs have significantly different expression (P<0.05).

Abbreviations: GEO, Gene Expression Omnibus; PTC, papillary thyroid carcinoma.

Finally, we performed the qRT-PCR to measure the actual expression of PRSS3P2, KRTAP5-AS1 and PWAR5. The results showed that PWAR5 was downregulated while PRSS3P2 and KRTAP5-AS1 were upregulated in the PTC tissues, which were consistent with the TCGA analysis (Figure 7).

Figure 7 qRT-PCR validation of the three lncRNAs.

Notes: Comparison of fold change (2−^∆∆Ct) of lncRNAs between TCGA and qRT-PCR results. Fold change (2−^∆∆Ct) of lncRNAs of qRT-PCR results: ∆Ct = Ct_lncRNAs – Ct_GAPDH and ∆∆Ct = ∆Ct_{tumor tissues} – ∆C_{tadjacent non-tumor tissues}.

Abbreviations: qRT-PCR, quantitative real time PCR; TCGA, The Cancer Genome Atlas.

Discussion

TC is a common malignant endocrine disorder worldwide with increasing incidence in the past few decades. The “cancer statistics 201,821”²² presented the estimated number of new cases of TC expected in America in 2018, pointing out that about 53,990 cases of TC would be expected to be diagnosed and about 2,060 patients would die of TC. PTC is the most common subtype, constituting 80%–85% of all TC cases.²³ In general, PTC has a relatively favorable prognosis with a high 5-year survival rate.²⁴ However, there are still about 5%–10% of patients suffering from the PTC recurrence and facing the aggressive and deathful outcomes.²⁵ Therefore, scholars are committed to find novel survival signatures of PTC. Liu et al²⁶ found that Yes-activated Protein-1 could play as the biomarker and high-risk indicator in PTC. Yi et al²⁷ analyzed the PTC data in TCGA database, suggesting that the ESR1 and ESR ratio (ESR1/ESR) could predict the female patients’ survival and has the potential to be a therapeutic target. A prospective study was performed to show that miR-146b-5p and miR-222 were significantly associated with central lymph node metastases in PTC, indicating that they were the independent predictors of PTC prognosis.²⁸ In the present study, we mined the TCGA data to obtain PTC-specific lncRNAs for predicting the patients’ survival.

To detect the PTC-specific lncRNAs, all cases were divided into four groups (stage I–II with lymphatic metastasis, stage I–II without lymphatic metastasis, stage III–IV with lymphatic metastasis and stage IV without lymphatic metastasis) of patient tumor tissues, and then they were compared with adjacent non-tumor tissues. Then, 83 key lncRNAs were selected for further studies. Among them, several studies that are closely associated with cancers have been reported. Jones et al²⁹ studied in a large sample to find that the telomerase RNA component variation promoted the increased risk of the colorectal cancer. lncRNA FER1L4 was downregulated in tumor tissues and could suppress tumor proliferation in endometrial carcinoma and gastric cancer.^30,31 Su et al³² performed the feedback loop study to find that the lncRNA FOXD2-AS1 was correlated with progression and recurrence in bladder cancer. Drawing inspiration from these cancer-specific signatures, we were interested in whether the PTC-specific lncRNAs could predict the patients’ prognosis.

After the Cox regression analysis, three prognosis-related lncRNAs (PRSS3P2, KRTAP5-AS1 and PWAR5) were obtained. It was reported that PRSS3P2, trypsinogen 6, was closely related to chronic pancreatitis.³³ Song et al³⁴ constructed ceRNA network of gastric cancer and claimed that lncRNA KRTAP5-AS1 could play the role of ceRNAs to regulate Claudin-4, a biomarker associated with patients’ prognosis. PWAR5, also known as PAR5, could regulate the proliferation and progression of glioma by binding to EZH2.³⁵ In addition, to validate the analyzed results from TCGA data, the expression of the abovementioned three key lncRNAs was measured using qRT-PCR. The results indicated that TCGA analysis and qRT-PCR results from 32 PTC patients were in 100% agreement. Combined with the literature and experimental results, these three PTC-specific lncRNAs had the potential to be prognostic indicators.

Then, we used these three survival-related lncRNAs to construct a three-lncRNA signature, which proved to be a prognostic biomarker of PTC. Based on the risk score model, the high-risk patients had poorer prognosis and the three-lncRNA signature could largely predict the 5-year survival time (AUC =0.739). Similarly, several studies also used the TCGA dataset to build the risk score model to explore the potential signature of PTC prognosis and diagnosis. Han et al³⁶ mined the data in the TCGA database, and then 19 significant gene pairs with high diagnostic ability were detected. Luo et al³⁷ identified a three-lncRNA signature (AC079630.2, CRNDE and CTD-2171N6.1) to predict the PTC patients’ survival. Li et al³⁸ also provided a four-lncRNA signature (RP11-536N17.1, RP11-508M8.1, AC026150.8 and CTD-2139B15.2) to improve the prognosis prediction of PTC patients. However, compared with the previous studies (Luo et al’s and Li et al’s studies), we had some advantages. First, we have considered the tumor stage and lymphatic metastasis which were closely associated with PTC prognosis to obtain the key lncRNAs for further study. Second, if the expression of genes was “0” whose patients were more than 10% of all data, they would be eliminated, which could avoid bias. Third, the power of 5-year survival prediction was detected by the survival ROC, and the risk score model had high prognostic power (AUC =0.739). Fourth, we used “FDR <0.05 and P<0.05” as the inclusion criteria.³⁹ Finally, we used GEO dataset and qRT-PCR to verify the results of TCGA analysis. We also wanted to validate the prognostic value of the three-lncRNA signature using the GEO datasets. Unfortunately, no survival information of PTC lncRNA was available in the datasets.³⁷ Therefore, the candidate lncRNAs in the risk sore model were different in our studies, and the three-lncRNA prognostic signature was credible to be a potential biomarker.

In addition, we investigated the function of the three-lncRNA signature. As shown in Figure 4, the co-expressed genes were mainly associated with the DNA-dependent, gene expression, ribosome and metabolic pathways. Abnormal regulation of these functions is involved in various cancers, such as lung cancer, breast cancer, bladder cancer and colorectal cancer.^40–43 Thus, the abnormal regulation of gene function may play a crucial role in the genesis and progress of PTC. Moreover, we identified three important transcription factors (ELK1, ETV7 and GABPA) that may significantly influence the PTC. Bullock et al⁴⁴ claimed that FOXE1 and ELK1 established a new regulatory pathway to regulate the TCs. Jacques et al⁴⁵ found that GABPA was closely associated with human thyroid tumors. However, there were no reports about the relation between ETV7 and TCs. In the future study, we can study on these three transcription factors by further research.

Finally, we analyzed the correlation between clinical characteristics and survival time. The results indicated that tumor stage, residual tumor, age and neoplasm cancer status were the risk factors of poor survival. A plenty of evidence supported our results that tumor stage, residual tumor and neoplasm cancer status were tightly bound to the poor prognosis.^16,46,47 According to AJCC, 45 years old is a common age-cutoff in PTC. (https://emedicine.medscape.com/article/2006643-overview).^48,49 In the current study, the PTC patients aged <45 years have longer survival time, indicating that the younger the PTC patients are, the better the prognosis they had.

Conclusion

By mining the TCGA data, we built a three-lncRNA signature, which could effectively predict the prognosis of PTC. Nevertheless, the present study has some limitations. First, there were only eleven dead PTC cases in the TCGA dataset, which might cause bias to obtain the correct results. Second, the samples of TCGA mostly were white people, and the results should be validated in a proof test of Chinese cohort in the future. Ultimately, we hope that the three-lncRNA signature can help to predict PTC prognosis in the future.

Acknowledgments

This study was supported by the Technology Development Project of Nanjing (YKK10080) and the National Natural Science Foundation of China (81673132 and 81472939). The authors appreciate Mr Cheng Donglin for his technical assistance.

Supplementary materials

Figure S1 Flowchart of bioinformatic analysis.

Figure S2 The K–M curves and predictive value of the risk score for clinical features.

Notes: The prognostic value of different clinical features for the OS of PTC patients. ROC curve is used to predict different clinical features.

Abbreviations: AUC, area under the curve; K–M, Kaplan–Meier; OS, overall survival; PTC, papillary thyroid carcinoma; ROC, receiver operating characteristic.

Table S1 Abnormal expression of the intersection of lncRNAs in PTC

lncRNA	Style	Fold change*
AADACP1	Down	0.092
LINC01139	Down	0.12
TPTE2P1	Down	0.15
LOC100130238	Down	0.15
LINC00473	Down	0.17
SLC26A4-AS1	Down	0.23
HAND2-AS1	Down	0.24
FAM167A-AS1	Down	0.24
LINC00092	Down	0.24
B3GALT5-AS1	Down	0.26
LOC143666	Down	0.28
MIR9-3HG	Down	0.28
LINC01257	Down	0.32
PWAR5	Down	0.33
TDH	Down	0.33
LINC01140	Down	0.33
SNORD116-4	Down	0.33
MIR4697HG	Down	0.35
TERC	Down	0.36
ATP6V0E2-AS1	Down	0.37
DPY19L2P4	Down	0.38
ST7-AS1	Down	0.39
LINC01550	Down	0.4
LINC00982	Down	0.4
LINC00602	Down	0.4
ZNF826P	Down	0.4
ANKRD20A8P	Down	0.4
SNORD116-20	Down	0.4
GOLGA8IP	Down	0.42
VLDLR-AS1	Down	0.42
PGM5-AS1	Down	0.42
FAM95B1	Down	0.42
LINC00652	Down	0.43
LINC00261	Down	0.43
LINC01126	Down	0.43
FAM181A-AS1	Down	0.44
LINC00936	Down	0.44
AGPAT4-IT1	Down	0.46
MIR22HG	Down	0.46
GVINP1	Down	0.47
FAR2P1	Down	0.48
LRRC37A6P	Down	0.49
KRTAP5-AS1	Up	2.01
GGT3P	Up	2.07
EGOT	Up	2.07
MBL1P	Up	2.11
PLEKHA8P1	Up	2.14
LINC00152	Up	2.19
CYP2B7P	Up	2.23
ESPNP	Up	2.23
DLEU2	Up	2.3
LINC01366	Up	2.5
FCGR1CP	Up	2.5
LBX2-AS1	Up	2.55
LOC285629	Up	2.6
GGT8P	Up	2.62
ABHD11-AS1	Up	2.76
SMIM10L2A	Up	2.76
FOXD2-AS1	Up	2.79
LOC100126784	Up	2.82
KRTAP5-AS1	Up	2.85
LOC93429	Up	2.95
FLJ23867	Up	2.95
ABCC6P2	Up	3.04
MIR31HG	Up	3.05
PP14571	Up	3.05
FER1L4	Up	3.09
HCG22	Up	3.17
LOC441204	Up	3.34
NR2F1-AS1	Up	3.52
CYP1B1-AS1	Up	3.53
LINC00887	Up	3.58
PRSS3P2	Up	3.7
MIR924HG	Up	3.83
SFTA1P	Up	3.84
FLJ16779	Up	3.88
YWHAEP7	Up	3.94
FBLL1	Up	4.36
EGFEM1P	Up	4.54
LINC00284	Up	4.65
TINCR	Up	4.79
ABCC6P1	Up	5.32
RPSAP52	Up	9.8
LOC400794	Up	20.13

Note:

* These lncRNAs have significantly different expression (fold change ≥2 or ≤0.5, and P-value <0.01).

Abbreviations: PTC, papillary thyroid carcinoma.

Disclosure

The authors report no conflicts of interest in this work.