Genetic expression profile-based screening of genes and pathways associated with papillary thyroid carcinoma

Identification of key biomarkers for thyroid cancer by integrative gene expression profiles

Thyroid cancer is a frequently diagnosed malignancy and the incidence has been increased rapidly in recent years. Despite the favorable prognosis of most thyroid cancer patients, advanced patients with metastasis and recurrence still have poor prognosis. Therefore, the molecular mechanisms of progression and targeted biomarkers were investigated for developing effective targets for treating thyroid cancer. Eight chip datasets from the gene expression omnibus database were selected and the inSilicoDb and inSilicoMerging R/Bioconductor packages were used to integrate and normalize them across platforms. After merging the eight gene expression omnibus datasets, we obtained one dataset that contained the expression profiles of 319 samples (188 tumor samples plus 131 normal thyroid tissue samples). After screening, we identified 594 significantly differentially expressed genes (277 up-regulated genes plus 317 down-regulated genes) between the tumor and normal tissue samples. The differentially expressed genes exhibited enrichment in multiple signaling pathways, such as p53 signaling. By building a protein-protein interaction network and module analysis, we confirmed seven hub genes, and they were all differentially expressed at all the clinical stages of thyroid cancer. A diagnostic seven-gene signature was established using a logistic regression model with the area under the receiver operating characteristic curve (AUC) of 0.967. Seven robust candidate biomarkers predictive of thyroid cancer were identified, and the obtained seven-gene signature may serve as a useful marker for thyroid cancer diagnosis and prognosis.

Identification of key genes of papillary thyroid carcinoma by integrated bioinformatics analysis

BACKGROUND

Papillary thyroid carcinoma (PTC) is one of the fastest-growing malignant tumor types of thyroid cancer. Therefore, identifying the interaction of genes in PTC is crucial for elucidating its pathogenesis and finding more specific molecular biomarkers.

METHODS

Four pairs of PTC tissues and adjacent tissues were sequenced using RNA-Seq, and 3745 differentially expressed genes were screened (P<0.05, |logFC|>1). The enrichment analysis indicated that the vast majority of differentially expressed genes (DEGs) may play a positive role in the development of cancer. Then, the significant modules were analyzed using Cytoscape software in the protein-protein interaction network. Survival analysis, TNM analysis, and immune infiltration analysis of key genes were analyzed. And the expression of ADORA1, APOE, and LPAR5 genes were verified by qPCR in PTC compared with matching adjacent tissues.

RESULTS

Twenty-five genes were identified as hub genes with nodes greater than 10. The expression of 25 genes were verified by the GEPIA database, and the overall survival and disease-free survival analyses were conducted with Kaplan-Meier plotter. We found only three genes were confirmed with our validation and were statistically significant in PTC, namely ADORA1, APOE, and LPAR5. Further analysis found that the mRNA levels and methylation degree of these three genes were significantly correlated with the TNM staging of PTC. And these three genes were related to PTC immune infiltration. Verification of the expression of these three genes by RT-qPCR and Western blot further confirmed the reliability of our results.

CONCLUSION

Our study identified three genes that may play key regulatory roles in the development, metastasis, and immune infiltration of papillary thyroid carcinoma.

Identification of Potential lncRNAs and miRNAs as Diagnostic Biomarkers for Papillary Thyroid Carcinoma Based on Machine Learning

BACKGROUND

Papillary thyroid carcinoma (PTC) accounts for most of the proportion of thyroid cancer (TC). The objective of this study was to identify diagnostic, differentially expressed long noncoding RNAs (lncRNAs) and microRNAs (miRNAs), contributing to understanding the epigenetics mechanism of PTC.

METHODS

The data of lncRNA, miRNA, and mRNA were downloaded from the Cancer Genome Atlas (TCGA) dataset, followed by functional analysis of differentially expressed mRNAs. Optimal diagnostic lncRNA and miRNA biomarkers were identified via random forest. The regulatory network between optimal diagnostic lncRNA and mRNAs and optimal diagnostic miRNA and mRNAs was identified, followed by the construction of ceRNA network of lncRNA-mRNA-miRNA. Expression validation and diagnostic analysis of lncRNAs, miRNAs, and mRNAs were performed. Overexpression of ADD3-AS1 was performed in PTC-UC3 cell lines, and cell proliferation and invasion assay were used for investigating the role of ADD3-AS1 in PTC.

RESULTS

A total of 107 differentially expressed lncRNAs, 81 differentially expressed miRNAs, and 515 differentially expressed mRNAs were identified. 11 lncRNAs and 6 miRNAs were regarded as the optimal diagnostic biomarkers for PTC. The epigenetic modifications via the above diagnostic lncRNAs and miRNAs were identified, including MIR181A2HG-FOXP2-hsa-miR-146b-3p, BLACAT1/ST7-AS1-RPS6KA5-hsa-miR-34a-5p, LBX2-AS1/MIR100HG-CDHR3-hsa-miR-34a-5p, ADD3-AS1-PTPRE-hsa-miR-9-5p, ADD3-AS1-TGFBR1-hsa-miR-214-3p, LINC00506-MMRN1-hsa-miR-4709-3p, and LOC339059-STK32A-hsa-miR-199b-5p. In the functional analysis, MMRN1 and TGFBR1 were involved in cell adhesion and endothelial cell migration, respectively. Overexpression of ADD3-AS1 inhibited cell growth and invasion in PTC cell lines.

CONCLUSION

The identified lncRNAs/miRNAs/mRNA were differentially expressed between normal and cancerous tissues. In addition, identified altered lncRNAs and miRNAs may be potential diagnostic biomarkers for PTC. Additionally, epigenetic modifications via the above lncRNAs and miRNAs may be involved in tumorigenesis of PTC.

Potential five-mRNA signature model for the prediction of prognosis in patients with papillary thyroid carcinoma

Although the mortality rate of papillary thyroid carcinoma (PTC) is relatively low, the recurrence rates of PTC remain high. The high recurrence rates are related to the difficulties in treatment. Gene expression profiles has provided novel insights into potential therapeutic targets and molecular biomarkers of PTC. The aim of the present study was to identify mRNA signatures which may categorize PTCs into high-and low-risk subgroups and aid with the predictions for prognoses. The mRNA expression profiles of PTC and normal thyroid tissue samples were obtained from The Cancer Genome Atlas (TCGA) database. Differentially expressed mRNAs were identified using the ‘EdgeR’ software package. Gene signatures associated with the overall survival of PTC were selected, and enrichment analysis was performed to explore the biological pathways and functions of the prognostic mRNAs using the Database for Visualization, Annotation and Integration Discovery. A signature model was established to investigate a specific and robust risk stratification for PTC. A total of 1,085 differentially expressed mRNAs were identified between the PTC and normal thyroid tissue samples. Among them, 361 mRNAs were associated with overall survival (P<0.05). A 5-mRNA prognostic signature for PTC (ADRA1B, RIPPLY3, PCOLCE, TEKT1 and SALL3) was identified to classify the patients into high-and low-risk subgroups. These prognostic mRNAs were enriched in Gene Ontology terms such as ‘calcium ion binding’, ‘enzyme inhibitor activity’, ‘carbohydrate binding’, ‘transcriptional activator activity’, ‘RNA polymerase II core promoter proximal region sequence-specific binding’ and ‘glutathione transferase activity’, and Kyoto Encyclopedia of Genes and Genomes signaling pathways such as ‘pertussis’, ‘ascorbate and aldarate metabolism’, ‘systemic lupus erythematosus’, ‘drug metabolism-cytochrome P450 and ‘complement and coagulation cascades’. The 5-mRNA signature model may be useful during consultations with patients with PTC to improve the prediction of their prognosis. In addition, the prognostic signature identified in the present study may reveal novel therapeutic targets for patients with PTC.