Abstract:Background: The genes and genetic mechanisms underlying the occurrence and progression of papillary thyroid carcinoma (PTC) are still unknown. This study aimed to find candidate genes related to the pathogenesis and progression of PTC. Methods: RNA sequencing (RNA-seq) data of PTC and normal tissues were downloaded from The Cancer Genome Atlas (TCGA) database with clinical stage data to form a test, validation, and clinical-stage data matrix. We used the test data set to analyze differentially expressed genes (DEGs) and weighted gene co-expression network analysis (WGCNA) to find those gene clusters highly correlated with PTC. We then verified the expression of genes in the interested modules using the validation matrix. The quantitative real-time polymerase chain reaction (qRT-PCR) was used to verify the reliability of the expression of selected genes. Five key genes (GDF15, LCN2, KCNN4, SH3BGRL3, and MMP2) were used to analyze the connection between gene expression and the American Joint Committee on Cancer (AJCC) stage. The upregulated and downregulated DEGs, along with the modules of interest, were subjected to Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment using the Database for Annotation, Visualization, and Integrated Discovery (DAVID). Results: We used WGCNA to find two modules of interest, the yellow module, which was positively associated with PTC, and the blue module, which was negatively correlated with PTC. Four genes (GDF15, LCN2, KCNN4, and SH3BGRL3) from the yellow module were determined to be highly expressed in PTC in the test data matrix and were verified in both the validation data matrix and quantitative real-time PCR, which indicated that these four genes were highly correlated with the occurrence of the PTC. Furthermore, these four genes also had a significantly higher expression in the advanced levels of pathological T, N, and AJCC stage, meaning that they were correlated with the progression of PTC. Genes in the yellow module and upregulated DEGs were significantly enriched in three vital signaling pathways, including focal adhesion, extracellular matrix (ECM)-receptor interaction, and the PI3K-Akt signaling pathway. Conclusions: Four candidate genes (GDF15, LCN2, KCNN4, and SH3BGRL3) may be potential biomarkers for the PTC's pathogenesis and may be useful for predicting the disease stage.

Identification of LINC02454-related key pathways and genes in papillary thyroid cancer by weighted gene coexpression network analysis (WGCNA)

Identification of Specific Long Non-Coding RNA Expression: Profile and Analysis of Association with Clinicopathologic Characteristics and BRAF Mutation in Papillary Thyroid Cancer

Identification of Key Genes and Evaluation of Clinical Outcomes in Lung Squamous Cell Carcinoma Using Integrated Bioinformatics Analysis

Overexpression of Novel Long Intergenic Non-Coding RNA LINC02454 is Associated with a Poor Prognosis in Papillary Thyroid Cancer

Identification of Potential Functional Genes in Papillary Thyroid Cancer by Co-Expression Network Analysis

Identification of candidate genes associated with papillary thyroid carcinoma pathogenesis and progression by weighted gene co-expression network analysis

Identification of Potential LncRNAs As Papillary Thyroid Carcinoma Biomarkers Based on Integrated Bioinformatics Analysis Using TCGA and RNA Sequencing Data

Identification of Key Genes Associated with Papillary Thyroid Microcarcinoma Characteristics by Integrating Transcriptome Sequencing and Weighted Gene Co-Expression Network Analysis.

Identification of novel gene signature predicting lymph node metastasis in papillary thyroid cancer via bioinformatics analysis and in vitro validation

Identification and Validation of Core Genes Involved in the Development of Papillary Thyroid Carcinoma via Bioinformatics Analysis

Integrated Analysis of Long Noncoding RNA Interactions Reveals the Potential Role in Progression of Human Papillary Thyroid Cancer.

Genomic Instability-Related LncRNA Signature Predicts the Prognosis and Highlights LINC01614 is a Tumor Microenvironment-Related Oncogenic Lncrna of Papillary Thyroid Carcinoma.

The Prognostic Value and Potential Immune Mechanisms of Lncrnas Related to Immunogenic Cell Death in Papillary Thyroid Carcinoma

Long Noncoding RNA LINC00284 Facilitates Cell Proliferation in Papillary Thyroid Cancer Via Impairing Mir-3127-5p Targeted E2F7 Suppression.

Integrated bioinformatics analysis and screening of hub genes in papillary thyroid carcinoma

Identification and Validation of Potential Target Genes in Papillary Thyroid Cancer.

Construction and analysis of an aberrant lncRNA-miRNA-mRNA network associated with papillary thyroid cancer

Identification of key immune genes related to lymphatic metastasis of papillary thyroid cancer via bioinformatics analysis and experimental validation

Identification and Analysis of Genes Associated with Papillary Thyroid Carcinoma by Bioinformatics Methods

Long Intergenic Non-Coding RNA 271 is Predictive of a Poorer Prognosis of Papillary Thyroid Cancer

LncRNA ENST00000539653 Acts As an Oncogenic Factor Via MAPK Signalling in Papillary Thyroid Cancer