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:In this study, 4 pairs of PTC tissues and adjacent tissues were sequenced using RNA-Seq, and 3745 differentially expressed genes (DEGs) were screened. The results of GO and KEGG enrichment analysis indicate that the vast majority of DEGs may play a positive role in the development of cancer. Then, the significant modules were analysed using Cytoscape software in the protein-protein interaction (PPI) network. Survival analysis, TNM analysis, and immune infiltration analysis of key genes are all analyzed. And the expression of ADORA1, APOE and LPAR5 genes was verified by qPCR in papillary thyroid carcinoma compared to their matching adjacent tissues.Results: A total of 25 genes were identified as hub genes with nodes greater than 10. The expression of 25 key genes in PTC were verified by the GEPIA database, and the overall survival and disease free survival analyses of these key genes were conducted with Kaplan–Meier plots. We found that 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 are significantly correlated with the TNM staging of PTC, and these three genes are related to PTC immune infiltration. Verification of the expression of these three genes by RT-qPCR 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.Key words :RNA-Seq, papillary thyroid carcinoma, key gene, bioinformatics