Identification of differentially expressed genes in papillary thyroid cancers

Construction of novel 7 integrin-related gene signatures in thyroid cancer construction of model based on integrin genes

Advanced and metastatic THCA patients usually have a poor prognosis. Thus, this study aimed to establish a risk model to discriminate the high risk population. The expression and clinical data were obtained from TCGA database. The cluster analysis, lasso, univariate and multivariate cox analyses were used to construct risk model. K-M, ROC and DCA were applied to validate the efficiency and stability of the model. GO, KEGG, and ssGSEA analysis were performed to identify the potential mechanism of signatures. The 7-gene prognosis model was constructed, including FAM27E3, FIGN, GSTM4, BEX5, RBPMS2, PHF13, and DCSTAMP. ROC and DCA results showed our model had a better prognosis prediction performance than other risk models. The high risk score was associated with the poor prognosis of THCA patients with different clinical characteristics. The risk score was closely related to cell cycle. Further, we found that the expressions of signatures were significantly dysregulated in THCA and associated with prognosis. These gene expressions were affected by some clinical characteristics, methylation and CNV. Some signatures played a role in drug sensitivity and pathway activation. We constructed a 7-gene signature model based on the integrin-related genes, which showed a great prognostic value in THCA.

Upregulation of dendrocyte-expressed seven transmembrane protein is associated with unfavorable outcomes in differentiated thyroid cancer

PURPOSE

Dendritic cell infiltrates are increased in thyroid cancer but may have a defective ability to provoke effective immune responses. In this study, we aimed to identify potential thyroid cancer biomarkers linked to dendritic cell development and evaluate their prognostic relevance.

METHODS

Through a bioinformatics search, we identified the dendrocyte-expressed seven transmembrane protein (DCSTAMP) as a prognostic gene involved in dendritic cell differentiation for thyroid cancer. Immunohistochemical analyses of DCSTAMP expression were performed and correlated with clinical outcomes.

RESULTS

DCSTAMP was overexpressed in a variety of types of thyroid cancers, while normal thyroid tissue or benign thyroid lesions exhibited low or undetectable DCSTAMP immunoreactivity. The results of automated quantification were consistent with subjective semiquantitative scoring. Among 144 patients with differentiated thyroid cancer, high DCSTAMP expression was associated with papillary tumor type (p < 0.001), extrathyroidal extension (p = 0.007), lymph node metastasis (p < 0.001), and BRAF V600E mutation (p = 0.029). Patients with tumors showing high DCSTAMP expression had shorter overall (p = 0.027) and recurrence-free (p = 0.042) survival.

CONCLUSION

This study provides the first evidence of DCSTAMP overexpression in thyroid cancer. Apart from the prognostic implications, studies are needed to explore its potential immunomodulatory role in thyroid cancer.

Systems biomarkers for papillary thyroid cancer prognosis and treatment through multi-omics networks

The identification of biomolecules associated with papillary thyroid cancer (PTC) has upmost importance for the elucidation of the disease mechanism and the development of effective diagnostic and treatment strategies. Despite particular findings in this regard, a holistic analysis encompassing molecular data from different biological levels has been lacking. In the present study, a meta-analysis of four transcriptome datasets was performed to identify gene expression signatures in PTC, and reporter molecules were determined by mapping gene expression data onto three major cellular networks, i.e., transcriptional regulatory, protein-protein interaction, and metabolic networks. We identified 282 common genes that were differentially expressed in all PTC datasets. In addition, six proteins (FYN, JUN, LYN, PML, SIN3A, and RARA), two Erb-B2 receptors (ERBB2 and ERBB4), two cyclin-dependent receptors (CDK1 and CDK2), and three histone deacetylase receptors (HDAC1, HDAC2, and HDAC3) came into prominence as proteomic signatures in addition to several metabolites including lactaldehyde and proline at the metabolome level. Significant associations with calcium and MAPK signaling pathways and transcriptional and post-transcriptional activities of 12 TFs and 110 miRNAs were also observed at the regulatory level. Among them, six miRNAs (miR-30b-3p, miR-15b-5p, let-7a-5p, miR-130b-3p, miR-424-5p, and miR-193b-3p) were associated with PTC for the first time in the literature, and the expression levels of miR-30b-3p, miR-15b-5p, and let-7a-5p were found to be predictive of disease prognosis. Drug repositioning and molecular docking simulations revealed that 5 drugs (prochlorperazine, meclizine, rottlerin, cephaeline, and tretinoin) may be useful in the treatment of PTC. Consequently, we report here biomolecule candidates that may be considered as prognostic biomarkers or potential therapeutic targets for further experimental and clinical trials for PTC.

The role of PI3K signaling pathway and its associated genes in papillary thyroid cancer

BACKGROUND

One of the well-differentiated types of thyroid cancer is papillary thyroid cancer (PTC), often developed by genetic mutations and radiation.

METHODS

In this study, the public NCBI-GEO database was systematically searched. The eligible datasets were the targets for biomarker discovery associated with PI3K signaling pathway.

RESULTS

Only two datasets were suitable and passed the inclusion criteria. The meta-analysis outcomes revealed eleven upregulation and thirteen downregulation genes differentially expressed between PTC and healthy tissues. Moreover, the outcomes for survival and disease-free rates for each gene were illustrated.

CONCLUSIONS

The present research suggests a panel signature of 24 gene biomarkers in diagnosing the PTC.

Cadmium-induced developmental alteration and upregulation of serine-type endopeptidase transcripts in wild freshwater populations of Chironomus plumosus

Cadmium, a toxic heavy metal, is a persistent environmental contaminant with irreversible toxicity to aquatic organisms. Chironomus plumosus, a natural species, is the largest sediment-burrowing aquatic midge in freshwater environments. In this study, we evaluated developmental defects in C. plumosus resulting from Cd exposure. In C. plumosus larvae, Cd exposure induced decreased survival and growth rates, reduction of emergence rate and sex ratio, and delayed emergence, as well as elevating the incidence of split tooth deformities. To identify potential biomarker genes to assess environmental pollutants such as Cd, we identified differentially expressed genes (DEGs) in C. plumosus exposed to various Cd concentrations. Among fourteen characterized DEGs, serine-type endopeptidase (SP) and heat shock protein 70 (HSP70) genes exhibited significant upregulation in C. plumosus larvae after Cd exposure. Therefore, we evaluated SP and HSP70 responses in natural C. plumosus populations collected from three sites of a Korean river and analyzed their correlations with eighteen environmental quality characteristics using principal component analysis. The highest expression of SP and HSP70 transcripts was observed in C. plumosus populations from Yeosu in Korea, which has high concentrations of polluting heavy metals. SP transcript expression was positively correlated with concentrations of Cd, Pb, Al, Fe, NO2, and NO3. These results suggested that environmental pollutants such as Cd can impair proteolytic activity in the digestive system of C. plumosus and may ultimately induce developmental alterations. We therefore suggest SP as a potential biomarker to assess the effects of environmental pollutants in aquatic ecosystems.

Identification of Long Non-Coding RNA Expression Profiles and Co-Expression Genes in Thyroid Carcinoma Based on The Cancer Genome Atlas (TCGA) Database

Background

Thyroid carcinoma is a malignancy with high morbidity and mortality. Genetic alterations play pivot roles in the pathogenesis of thyroid carcinoma, where long noncoding RNA (lncRNA) have been identified to be crucial. This study sought to investigate the biological functions of lncRNA expression profiles in thyroid carcinoma.

Material/Methods

The lncRNAs expression profiles were acquired from The Cancer Genome Atlas (TCGA) database according to 510 thyroid cancer tissues and 58 normal thyroid tissues. By using R package edgeR, differentially expressed RNAs were obtained. Also, an overall survival model was established based on Cox regression and clinical data then testified by Kaplan-Meier plot, receiver operating characteristic (ROC)-curve and C-index analysis. We investigated the co-expressed genes with lncRNAs involved in the prognostic model, as well as Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was conducted R package clusterProfile.

Results

A total of 352 lncRNAs were identified as differentially expressed in thyroid carcinoma, and an overall survival model consisting of 8 signature lncRNAs was proposed (ROC=0.862, C-index=0.893, P<0.05), 3 of which (DOCK9-DT, FAM111A-DT, and LINC01736) represent co-expressed mRNAs. However, as an oncogene, only FAM111A-DT increased the prognostic risk in thyroid carcinoma. Furthermore, we found differential genes LINC01016, LHX1-DT, IGF2-AS, ND MIR1-1HG-AS1, significantly related to lymph node metastasis (P<0.05).

Conclusions

In this study, we clarified the differential lncRNA expression profiles which were related to the tumorigenesis and prognosis in thyroid carcinoma. Our results provide new rationale and understandings to the pathogenesis and regulatory mechanisms of thyroid carcinoma.

Identifying gene modules of thyroid cancer associated with pathological stage by weighted gene co-expression network analysis

Thyroid cancer is the most common type of endocrine tumor. The TNM classification remains a standard for treatment determination and predicting prognosis in thyroid cancer. The genes modules associated with the progression of papillary thyroid carcinoma (PTC) were not clear. We applied a weighted gene co-expression network analysis (WGCNA) and differential expression analysis to systematically identified co-expressed gene modules and hub genes associated with PTC progression based on The Cancer Genome Atlas (TCGA) PTC transcriptome sequencing data. An independent validation cohort, GSE27155, was used to evaluate the preservation of gene modules. We identified two co-expressed genes modules associated with progression of PTC. Enrichment analysis indicated that the two modules were enriched in angiogenesis and extracellular matrix organization. DCN, COL1A1, COL1A2, COL5A2 and COL3A1 were hub genes in the co-expressed network. We systematically identified co-expressed gene modules and hub genes associated with PTC progression for the first time, which provided insights into the mechanisms underlying PTC progression and some potential targets for the treatment of PTC.

Identification of pathways and genes associated with synovitis in osteoarthritis using bioinformatics analyses

Synovitis in osteoarthritis (OA) is a very common condition. However, its underlying mechanism is still not well understood. This study aimed to explore the molecular mechanisms of synovitis in OA. The gene expression profile GSE82107 (downloaded from the Gene Expression Omnibus database) included 10 synovial tissues of the OA patients and 7 synovial tissues of healthy people. Subsequently, differentially expressed gene (DEG) analysis, GO (gene ontology) enrichment analysis, pathway analysis, pathway network analysis, and gene signal network analysis were performed using Gene-Cloud of Biotechnology Information (GCBI). A total of 1,941 DEGs consisting of 1,471 upregulated genes and 470 downregulated genes were determined. Genes such as PSMG3, LRP12 MIA-RAB4B, ETHE1, SFXN1, DAZAP1, RABEP2, and C9orf16 were significantly regulated in synovitis of OA. In particular, the MAPK signalling pathway, apoptosis, and pathways in cancer played the most important roles in the pathway network. The relationships between these pathways were also analysed. Genes such as NRAS, SPHK2, FOS, CXCR4, PLD1, GNAI2, and PLA2G4F were strongly implicated in synovitis of OA. In summary, this study indicated that several molecular mechanisms were implicated in the development and progression of synovitis in OA, thus improving our understanding of OA and offering molecular targets for future therapeutic advances.

Gene expression profile analyze the molecular mechanism of CXCR7 regulating papillary thyroid carcinoma growth and metastasis

Background

To detect genetic expression profile alterations after papillary thyroid carcinoma (PTC) cells transfected with chemokine receptor CXCR7 gene by gene microarray, and gain insights into molecular mechanisms of how CXCR7 regulating PTC growth and metastasis.

Methods

The Human OneArray microarray was used for a complete genome-wide transcript profiling of CXCR7 transfected PTCs (K1-CXCR7 cells), defined as experimental group. Non CXCR7 transfected PTCs (K1 cells) were used as control group. Differential analysis for per gene was performed with a random variance model and t test, p values were adjusted to control the false discovery rate. Gene ontology (GO) on differentially expressed genes to identify the biological processes in modulating the progression of papillary thyroid carcinoma. Pathway analysis was used to evaluate the signaling pathway that differentially expressed genes were involved in. In addition, quantitative real-time polymerase chain reaction (q-PCR) and Western blot were used to verify the top differentially expression genes.

Results

Comparative analysis revealed that the expression level of 1149 genes was changed in response to CXCR7 transfection. After unsupervised hierarchical clustering analysis, 270 differentially expressed genes were filtered, of them 156 genes were up-regulated whereas 114 genes were down-regulated in K1-CXCR7 cells. GO enrichment analysis revealed the differentially expressed genes were mainly involved in biopolymer metabolic process, signal transduction and protein metabolism. Pathway enrichment analysis revealed differentially expressed genes were mainly involved in ECM-receptor interaction, Focal adhesion, MAPK signaling pathway and Cytokine-cytokine receptor interaction pathway. More importantly, the expression level of genes closely associated with tumor growth and metastasis was altered significantly in K1-CXCR7 cells, including up-regulated genes FN1, COL1A1, COL4A1, PDGFRB, LTB, CXCL12, MMP-11, MT1-MMP and down-regulated genes ITGA7, and Notch-1.

Conclusions

Gene expression profiling analysis of papillary thyroid carcinoma can further delineate the mechanistic insights on how CXCR7 regulating papillary thyroid carcinoma growth and metastasis. CXCR7 may regulate growth and metastasis of papillary thyroid carcinoma via the activation of PI3K/AKT pathway and its downstream NF-κB signaling, as well as the down-regulation of Notch signaling.

Proliferation rate of stem cells derived from human dental pulp and identification of differentially expressed genes

Stem cells from human exfoliated deciduous teeth (SHED) and dental pulp stem cells (DPSCs) obtained from the dental pulp of human extracted tooth were cultured and characterized to confirm that these were mesenchymal stem cells. The proliferation rate was assessed using AlamarBlue® cell assay. The differentially expressed genes in SHED and DPSCs were identified using the GeneFishing™ technique. The proliferation rate of SHED (P < 0.05) was significantly higher than DPSCs while SHED had a lower multiplication rate and shorter population doubling time (0.01429, 60.57 h) than DPSCs (0.00286, 472.43 h). Two bands were highly expressed in SHED and three bands in DPSCs. Sequencing analysis showed these to be TIMP metallopeptidase inhibitor 1 (TIMP1), and ribosomal protein s8, (RPS8) in SHED and collagen, type I, alpha 1, (COL1A1), follistatin-like 1 (FSTL1), lectin, galactoside-binding, soluble, 1, (LGALS1) in DPSCs. TIMP1 is involved in degradation of the extracellular matrix, cell proliferation and anti-apoptotic function and RPS8 is involved as a rate-limiting factor in translational regulation; COL1A1 is involved in the resistance and elasticity of the tissues; FSTL1 is an autoantigen associated with rheumatoid arthritis; LGALS1 is involved in cell growth, differentiation, adhesion, RNA processing, apoptosis and malignant transformation. This, along with further protein expression analysis, holds promise in tissue engineering and regenerative medicine.

Brahma‑related gene 1-associated expression of 9-27 and IFI-27 is involved in acquired cisplatin resistance of gastric cancer cells

In order to investigate the mechanism of cisplatin resistance, a cisplatin-resistant human gastric cancer cell line was established. Subsequent to the exposure of the YCC-3 gastric cancer cell line to equal concentrations of cis-diammine-dichloroplatinum (II) (cisplatin, CDDP) for 6 months, a cisplatin-resistant cell line was established (YCC-3/R). To determine the molecular mechanism of cisplatin resistance in YCC-3/R cells, differentially expressed genes (DEGs) were investigated between YCC-3 and YCC-3/R by annealing control primer-based reverse transcriptase-polymerase chain reaction (ACP RT-PCR) technology. Eleven DEGs were successfully identified and sequenced. Among them, interferon-induced transmembrane protein 1 (9-27) and interferon α-inducible protein 27 (IFI-27) were markedly increased in YCC-3/R cells. In addition, western blot analysis demonstrated that the Brahma-related gene 1 (BRG1), which was observed to selectively activate 9-27 and IFI-27 genes, was overexpressed in YCC-3/R cells. The results suggested that the BRG1‑associated expression of 9-27 and IFI-27 is involved in cisplatin resistance in gastric cancer cells.

Evaluation of MMP2 and Caspase-3 expression in 107 cases of papillary thyroid carcinoma and its association with prognostic factors

Papillary thyroid carcinoma (PTC), including its variants and widely varying behavior, constitutes about 80% of all thyroid malignancies. Increased knowledge regarding molecular alterations has led to attempts to identify diagnostic or prognostic factors for a reliable preoperative approach to the classification of patients according to risk of recurrence. In this study, 107 cases of PTC with known histological properties, including vascular or capsular invasion, were assessed for expression of MMP2 and Caspase-3 using immunohistochemistry. Considering 10% as a cutoff to discriminate cases with invasive behavior from the non-invasive group, there was no relationship between expression of MMP2 or Caspase-3 in tumor cells and the presence of capsular invasion (p=0.45 and 0.64, respectively), as well as for the expression of Caspase-3 and vascular invasion (p=0.43). In case of MMP2, a borderline correlation was found between the positive reaction of tumor cells with the presence of vascular invasion (p=0.05). So the evaluation of MMP2 in thyroid PTC appears to be of some benefit to the prediction of tumor behavior while Caspase-3 as a marker of prediction seems to be of no use.

Cytological results of ultrasound-guided fine-needle aspiration cytology for thyroid nodules: emphasis on correlation with sonographic findings

PURPOSE

To compare the cytological results of ultrasound-guided fine-needle aspiration (US-FNA) cytology of thyroid nodules to sonographic findings and determine whether US findings are helpful in the interpretation of cytological results.

MATERIALS AND METHODS

Among the thyroid nodules that underwent US-FNA cytology, we included the 819 nodules which had a conclusive diagnosis. Final diagnosis was based on pathology from surgery, repeated FNA cytology or follow-up of more than one year. Cytological results were divided into five groups: benign, indeterminate (follicular or Hurthle cell neoplasm), suspicious for malignancy, malignant, and inadequate. US findings were categorized as benign or suspicious. Cytological results and US categories were analyzed.

RESULTS

Final diagnosis was concluded upon in 819 nodules based on pathology (n=311), repeated FNA cytology (n=204) and follow-up (n=304), of which 634 were benign and 185 were malignant. There were 560 benign nodules, 141 malignant nodules, 49 nodules with inadequate results, 21 with indeterminate results, and 48 that were suspicious for malignancy. The positive and negative predictive values of the US categories were 59.1% and 97.0%, and those of the cytological results were 93.7% and 98.9%. The US categories were significantly correlated with final diagnosis in the benign (p=0.014) and suspicious for malignancy (p<0.001) cytological result groups, but not in the inadequate and indeterminate cytological results groups. The false positive and negative rates of cytological results were 1.9% and 3.2%.

CONCLUSION

Sonographic findings can be useful when used alongside cytological results, especially in nodules with cytological results that are benign or suspicious for malignancy.

Differential expression of a set of genes in follicular and classic variants of papillary thyroid carcinoma

Fine-needle aspiration biopsy (FNA) is currently the best initial diagnostic test for evaluation of a thyroid nodule. FNA cytology cannot discriminate between benign and malignant thyroid nodules in up to 30% of thyroid nodules. Therefore, an adjunct to FNA is needed to clarify these lesions as benign or malignant. Using differential display-polymerase chain reaction method, the gene expression differences between follicular and classic variants of papillary thyroid carcinoma (PTC) and benign thyroid nodules were evaluated in a group of 42 patients. Computational gene function analyses via Cytoscape, FuncBASE, and GeneMANIA led us to a functional network of 17 genes in which a core sub-network of five genes coexists. Although the exact mechanisms underlying in thyroid cancer biogenesis are not currently known, our data suggest that the pattern of transformation from healthy cells to cancer cells of PTC is different in follicular variant than in classic variant.

Microarray analysis of papillary thyroid cancers in Korean

BACKGROUND/AIMS

Papillary thyroid cancer (PTC) is the most common malignancy of the thyroid gland. It involves several molecular mechanisms. The BRAF V600E mutation has been identified as the most common genetic abnormality in PTC. Moreover, it is known to be more prevalent in Korean PTC patients than in patients from other countries. We investigated distinct genetic profiles in Korean PTC through cDNA microarray analysis.

METHODS

Transcriptional profiles of five PTC samples and five paired normal thyroid tissue samples were generated using cDNA microarrays. The tumors were genotyped for BRAF mutations. The results of the cDNA microarray gene expression analysis were confirmed by real-time PCR and immunohistochemistry analysis of 35 PTC patients.

RESULTS

Four of the five patients whose PTC tissues were subjected to microarray analysis were found to carry the BRAF V600E mutation. Microarrays analysis of the five PTC tissue samples showed the expression of 96 genes to be increased and that of 16 genes decreased. Real-time reverse transcription-polymerase chain reaction (RT-PCR) confirmed increased expression of SLC34A2, TM7SF4, COMP, KLK7, and KCNJ2 and decreased expression of FOXA2, SLC4A4, LYVE-1, and TFCP2L1 in PTC compared with normal tissue. Of these genes, TFCP2L1, LYVE-1, and KLK7 were previously unidentified in PTC microarray analysis. Notably, Foxa2 activity in PTC was reduced, as shown by its cytoplasmic localization, in immunohistochemical analyses.

CONCLUSIONS

These findings demonstrate both similarities and differences between our results and previous reports. In Korean cases of PTC, Foxa2 activity was reduced with its cytoplasmic accumulation. Further studies are needed to confirm the relationship between FOXA2 and BRAF mutations in Korean cases of PTC.