RNA sequencing identifies multiple fusion transcripts, differentially expressed genes, and reduced expression of immune function genes in BRAF (V600E) mutant vs BRAF wild-type papillary thyroid carcinoma

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

BACKGROUND

Our previous study demonstrated that long intergenic noncoding RNA 02454 (LINC02454) may act as an oncogene to promote the proliferation and inhibit the apoptosis of papillary thyroid cancer (PTC) cells. This study was designed to investigate the mechanisms whereby LINC02454 is related to PTC tumorigenesis.

METHODS

Thyroid cancer RNA sequence data were obtained from The Cancer Genome Atlas (TCGA) database. Weighted gene coexpression network analysis (WGCNA) was applied to identify modules closely associated with PTC. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was used to identify the key pathways, and the maximal clique centrality (MCC) topological method was used to identify the hub genes. The Gene Expression Profiling Interactive Analysis (GEPIA) database was used to compare expression levels of key genes between PTC samples and normal samples and explore the prognostic value of key genes. The key genes were further validated with GEO dataset.

RESULTS

The top 5000 variable genes were investigated, followed by an analysis of 8 modules, and the turquoise module was the most positively correlated with the clinical stage of PTC. KEGG pathway analysis found the top two pathways of the ECM - receptor interaction and MAPK signaling pathway. In addition, five key genes (FN1, LAMB3, ITGA3, SDC4, and IL1RAP) were identified through the MCC algorithm and KEGG analysis. The expression levels of the five key genes were significantly upregulated in thyroid cancer in both TCGA and GEO datasets, and of these five genes, FN1 and ITGA3 were associated with poor disease-free prognosis.

CONCLUSIONS

Our study identified five key genes and two key pathways associated with LINC02454, which might shed light on the underlying mechanism of LINC02454 action in PTC.

WZ-3146 acts as a novel small molecule inhibitor of KIF4A to inhibit glioma progression by inducing apoptosis

BACKGROUND

Glioma is considered the most common primary malignant tumor of the central nervous system. Although traditional treatments have not achieved satisfactory outcomes, recently, targeted therapies for glioma have shown promising efficacy. However, due to the single-target nature of targeted therapy, traditional targeted therapies are ineffective; thus, novel therapeutic targets are urgently needed.

METHODS

The gene expression data for glioma patients were derived from the GEO (GSE4290, GSE50161), TCGA and CGGA databases. Next, the upregulated genes obtained from the above databases were cross-analyzed, finally, 10 overlapping genes (BIRC5, FOXM1, EZH2, CDK1, KIF11, KIF4A, NDC80, PBK, RRM2, and TOP2A) were ultimately screened and only KIF4A expression has the strongest correlation with clinical characteristics in glioma patients. Futher, the TCGA and CGGA database were utilized to explore the correlation of KIF4A expression with glioma prognosis. Then, qRT-PCR and Western blot was used to detect the KIF4A mRNA and protein expression level in glioma cells, respectively. And WZ-3146, the small molecule inhibitor targeting KIF4A, were screened by Cmap analysis. Subsequently, the effect of KIF4A knockdown or WZ-3146 treatment on glioma was measured by the MTT, EdU, Colony formation assay and Transwell assay. Ultimately, GSEA enrichment analysis was performed to find that the apoptotic pathway could be regulated by KIF4A in glioma, in addition, the effect of WZ-3146 on glioma apoptosis was detected by flow cytometry and Western blot.

RESULTS

In the present study, we confirmed that KIF4A is abnormally overexpressed in glioma. In addition, KIF4A overexpression is a key indicator of glioma prognosis; moreover, suppressing KIF4A expression can inhibit glioma progression. We also discovered that WZ-3146, a small molecule inhibitor of KIF4A, can induce apoptosis in glioma cells and exhibit antiglioma effects.

CONCLUSION

In conclusion, these observations demonstrated that targeting KIF4A can inhibit glioma progression. With further research, WZ-3146, a small molecule inhibitor of KIF4A, could be combined with other molecular targeted drugs to cooperatively inhibit glioma progression.

Comparative efficiency of differential diagnostic methods for the identification of BRAF V600E gene mutation in papillary thyroid cancer (Review)

V-Raf murine sarcoma viral oncogene homolog B1 (BRAF) encodes a serine-threonine kinase. The V600E point mutation in the BRAF gene is the most common mutation, predominantly occurring in melanoma, and colorectal, thyroid and non-small cell lung cancer. Particularly in the context of thyroid cancer research, it is routinely employed as a molecular biomarker to assist in diagnosing and predicting the prognosis of papillary thyroid cancer (PTC), and to formulate targeted therapeutic strategies. Currently, several methods are utilized in clinical settings to detect BRAF V600E mutations in patients with PTC. However, the sensitivity and specificity of various detection techniques vary significantly, resulting in diverse detection outcomes. The present review highlights the advantages and disadvantages of the methods currently employed in medical practice, with the aim of guiding clinicians and researchers in selecting the most suitable detection approach for its high sensitivity, reproducibility and potential to develop targeted therapeutic regimens for patients with BRAF gene mutation-associated PTC.

The role of MAPK, notch and Wnt signaling pathways in papillary thyroid cancer: Evidence from a systematic review and meta-analyzing microarray datasets employing bioinformatics knowledge and literature

Papillary thyroid cancer (PTC) is a prevalent kind of thyroid cancer (TC), with the risk of metastasis increasing faster than any other malignancy. So, understanding the role of PTC in pathogenesis requires studying the various gene expressions to find out which particular molecular biomarkers will be helpful. The authors conducted a comprehensive search on the PubMed microarray database and a meta-analysis approach on the remaining ones to determine the differentially expressed genes between PTC and normal tissues, along with the analyses of overall survival (OS) and recurrence-free survival (RFS) rates in patients with PTC. We considered the associated genes with MAPK, Wnt, and Notch signaling pathways. Two GEO datasets have been included in this research, considering inclusion and exclusion criteria. Nineteen genes were found to have higher differences through the meta-analysis procedure. Among them, ten genes were upregulated, and nine genes were downregulated. The expression of 19 genes was examined using the GEPIA2 database, and the Kaplan-Meier plot statistics were used to analyze RFS and the OS rates. We discovered seven significant genes with the validation: PRICKLE1, KIT, RPS6KA5, GADD45B, FGFR2, FGF7, and DTX4. To further explain these findings, it was discovered that the mRNA expression levels of these seven genes and the remaining 12 genes were shown to be substantially linked with the results of the experimental literature investigations on the PTC. Our research found nineteen panels of genes that could be involved in the PTC progression and metastasis and the immune system infiltration of these cancers.

Integration of metabolomics and transcriptomics reveals metformin suppresses thyroid cancer progression via inhibiting glycolysis and restraining DNA replication

The anti-tumoral effects of metformin have been widely studied in several types of cancer, including thyroid cancer; however, the underlying molecular mechanisms remain poorly understood. As an oral hypoglycemic drug, metformin facilitates glucose catabolism and disrupts metabolic homeostasis. Metabolic reprogramming, particularly cellular glucose metabolism, is an important characteristic of malignant tumors. This study aimed to explore the therapeutic effects of metformin in thyroid cancer and the underlying metabolic mechanism. In the present study, it was shown that metformin reduced cell viability, invasion, migration, and EMT, and induced apoptosis and cell cycle G1 phase arrest in thyroid cancer. Transcriptome analysis demonstrated that the differentially expressed genes induced by metformin were involved in several signaling pathways including apoptosis singling pathways, TGF-β signaling, and cell cycle regulation in human thyroid cancer cell lines. In addition, the helicase activity of the CDC45-MCM2-7-GINS complex and DNA replication related genes such as RPA2, RAD51, and PCNA were downregulated in metformin-treated thyroid cancer cells. Moreover, metabolomics analysis showed that metformin-induced significant alterations in metabolic pathways such as glutathione metabolism and polyamine synthesis. Integrative analysis of transcriptomes and metabolomics revealed that metformin suppressed glycolysis by downregulating the key glycolytic enzymes LDHA and PKM2 and upregulating IDH1 expression in thyroid cancer. Furthermore, the anti-tumor role of metformin in thyroid cancer in vivo was shown. Together these results show that metformin plays an anti-tumor role by inhibiting glycolysis and restraining DNA replication in thyroid cancer.

BRAF, TERT and HLA-G Status in the Papillary Thyroid Carcinoma: A Clinicopathological Association Study

As BRAF, TERT, HLA-G, and microRNAs have been individually associated with papillary thyroid carcinoma (PTC), we aimed to evaluate the individual and collaborative role of these markers in PTC in the same patient cohort. HLA-G and BRAF tumor expression was evaluated by immunohistochemistry. Using molecular methods, BRAFV600E and TERT promoter mutations were evaluated in thyroid fine needle aspirates. MicroRNA tumor profiling was investigated using massively parallel sequencing. We observed strong HLA-G (67.96%) while BRAF (62.43%) staining was observed in PTC specimens. BRAF overexpression was associated with poor response to therapy. The BRAFV600E (52.9%) and TERTC228T (13%) mutations were associated with extrathyroidal extension, advanced-age, and advanced-stage cancer. The TERT rs2853669 CC+TC genotypes (38%) were overrepresented in metastatic tumors. Nine modulated microRNAs targeting the BRAF, TERT, and/or HLA-G genes were observed in PTC and involved with cancer-related signaling pathways. The markers were individually associated with PTC features, emphasizing the synergistic effect of BRAFV600E and TERTC228T; however, their collaborative role on PTC outcome was not fully demonstrated. The differentially expressed miRNAs targeting the BRAF and/or HLA-G genes may explain their increased expression in the tumor milieu.

Identifying the Carcinogenic Mechanism of Malignant Struma Ovarii Using Whole-Exome Sequencing and DNA Methylation Analysis

Background: Since malignant struma ovarii is a very rare disease, its carcinogenic mechanism has not been elucidated. Here, we sought to identify the genetic lesions that may have led to the carcinogenesis of a rare case of malignant struma ovarii (follicular carcinoma) with peritoneal dissemination. Methods: DNA was extracted from the paraffin-embedded sections of normal uterine tissues and malignant struma ovarii for genetic analysis. Whole-exome sequencing and DNA methylation analysis were then performed. Results: Germline variants of RECQL4, CNTNAP2, and PRDM2, which are tumor-suppressor genes, were detected by whole-exome sequencing. Somatic uniparental disomy (UPD) was also observed in these three genes. Additionally, the methylation of FRMD6-AS2, SESN3, CYTL1, MIR4429, HIF3A, and ATP1B2, which are associated with tumor growth suppression, was detected by DNA methylation analysis. Conclusions: Somatic UPD and DNA methylation in tumor suppressor genes may be associated with the pathogenesis of malignant struma ovarii. To our knowledge, this is the first report of whole-exome sequencing and DNA methylation analysis in malignant struma ovarii. Genetic and DNA methylation analysis may help elucidate the mechanism of carcinogenesis in rare diseases and guide treatment decisions.

Reveal the potential molecular mechanism of circRNA regulating immune-related mRNA through sponge miRNA in the occurrence and immune regulation of papillary thyroid cancer

BACKGROUND

Papillary thyroid cancer (PTC) is the most common endocrine malignant tumour. The purpose of this study was to explore the potential molecular mechanism of circRNA regulating immune-related mRNA through sponge miRNA in the occurrence and immune regulation of PTC.

METHODS

All data were downloaded from public databases, such as GEO, Immport and TCGA. Differentially expressed (DE) mRNAs (DEmRNAs), DEmiRNAs and DEcircRNAs were identified using metaMA and limma packages. Subsequently, immune-related DEmRNAs were screened, and circRNA-miRNA-mRNA (ceRNA) regulatory network was constructed. In addition, functional annotation, protein-protein interaction (PPI) network construction, immune cell infiltration analysis and Pearson correlation analysis were performed. Finally, qRT-PCR validation and cell experiments were also performed.

RESULTS

In total, 2962 DEmRNAs, 78 DEmiRNAs and 51 DEcircRNAs were obtained. Subsequently, 195 immune-related DEmRNAs were obtained based on Immport database. Cytokine-cytokine receptor interaction was the only signalling pathway obtained in KEGG analysis. Then, 8 hub immune-related DEmRNAs were identified based on PPI network and CytoHubba plug-in. Subsequently, ceRNA sub-network containing hub immune-related DEmRNAs was extracted from ceRNA regulatory network. In ceRNA sub-network, hsa_circ_0082182-hsa-miR-18b-5p-FGF1/PDGFC, hsa_circ_0016404-hsa-miR-1275-FGF1/CTSB/IL13RA1, hsa_circ_0070100-hsa-miR-27a-3p/hsa-miR-27b-3p-TGFBR3, hsa_circ_0060055/hsa_circ_0038718-hsa-miR-150-3p-CXCL14, hsa_circ_0030427/hsa_circ_0002917-hsa-miR-22-3p-BMP7 and hsa_circ_0030427/hsa_circ_0002917-hsa-miR-125a-5p-LIFR axes were identified. Moreover, FGF1, PDGFC, CTSB, IL13RA1, TGFBR3, CXCL14, BMP7, LIFR, hsa-miR-125a-5p, hsa-miR-1275, hsa-miR-150-3p, hsa-miR-18b-5p and hsa-miR-27b-3p were also found to have good diagnostic accuracy and may be potential novel diagnostic markers for PTC. XCell analysis showed that the levels of immune cell infiltration (including Tregs, HSC, DC and Monocytes) were significantly different between the PTC and the control groups. Knockdown of the expression of hsa_circ_0082182 significantly inhibits the activity, proliferation, migration and invasion of TPC-1 cells.

CONCLUSION

Several circRNA-miRNA-mRNA axes identified in this study may be related to the occurrence, progression and survival of PTC. This lays a theoretical foundation for further understanding the molecular mechanism of PTC, and also contributes to clinical management and research.

Berberine inhibits IFN-γ signaling pathway in DSS-induced ulcerative colitis

Aims

The potential signaling pathways and core genes in ulcerative colitis (UC) were investigated in this study. Furthermore, potential mechanisms of BBR in treating UC were also explored.

Methods

Expression profiling by array of UC patients were obtained from Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were determined with the differential analysis. The biological functions of DEGs were analyzed through the Database for Annotation, Visualization and Integrated Discovery (DAVID). The Gene Set Enrichment Analysis (GSEA) was applied to analyze the expression differences between two different phenotype sample sets. Dextran sulfate sodium (DSS) was applied to establish UC model of mice and lipopolysaccharide (LPS) was utilized to induce inflammatory damage of NCM460 cells. Therapeutic effects of berberine (BBR) on disease performance, pathologic changes and serum supernatant indices were analyzed in vivo. To further investigate the potential mechanisms of BBR in treating UC, the expression of genes and proteins in vivo and in vitro were examined by RT-qPCR, immunohistochemical staining and western blotting.

Results

Immune-inflammatory genes were identified and up-regulated significantly in UC patients. In addition, IFN-γ signaling pathway and its core genes were significantly up-regulated in the phenotype of UC. All disease performance and the pathologic changes of UC in mice were evidently ameliorated by BBR treatment. The pro-inflammatory cytokines of serum, including CXCL9, CXCL1, IL-17 and TNF-α, in UC mice were significantly reduced by treatment of BBR. In terms of mechanisms of BBR in treating UC, the pro-inflammatory and immune-related genes, encoding IFN-γ, IRF8, NF-κB and TNF-α decreased significantly in UC mice followed by BBR treatment. Meanwhile, the expression of IFN-γ and its initiated targets, including IRF8, Ifit1, Ifit3, IRF1, were suppressed significantly by BBR treatment in vivo. The blocking of IFN-γ in vitro led to the silence of IFN-γ signaling pathway after exposure to BBR. Furthermore, the blocking of IFN-γ in vitro led to the silence of IFN-γ signaling pathway after exposure to BBR.

Conclusion

BBR holds anti-inflammatory activity and can treat UC effectively. The anti-inflammatory property of BBR is tightly related to the suppression of IFN-γ signaling pathway, which is crucial in immune-inflammatory responses of the colon mucosa.

Identifying key genes of classic papillary thyroid cancer in women aged more than 55 years old using bioinformatics analysis

BACKGROUND

The incidence rate of thyroid carcinoma (THCA) markedly increased in the recent few decades and has been likely over-diagnosed, especially papillary thyroid cancer (PTC) in women. However, the incidence of advanced-stage papillary thyroid cancer is also rising. According to earlier studies, tumors with identical pathology might have different clinical outcomes, which implies some variances in papillary thyroid cancer. Although the mortality of thyroid cancer has remained stable or declined, there is still an important problem in estimating whether it is benign or needs surgery for patients with papillary thyroid cancer.

METHODS

After obtaining data from The Cancer Genome Atlas (TCGA) Project-THCA database by R package TCGA bio links, 18 samples (11 at stage IV as high-risk group and 7 at stage I as low-risk group) were obtained using survival package and edgeR to ensure differential expression; ClusterProfiler package was used to carry on gene set enrichment analysis and searched the possible pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. STRING and Cytoscape were used to construct and modify the protein-protein interaction (PPI) network to get hub genes of differentially expressed genes. Next, the pROC package was used to get the receiver operating characteristic (ROC) curves of hub genes' disease-free survival (DFS). Then, transcription factors (TFs) and miRNAs of key genes were predicted by ENCORI and AnimalTFDB. In the end, TF-target genes-miRNA regulatory network was also constructed by Cytoscape.

RESULTS

Our research obtained the top 9 candidate genes from the whole network (IFNA1, MRC1, LGALS3, LOX, POSTN, TIMP1, CD276, SDC4, and TLR2). According to the ROC results, TIMP1, LOX, CD276, IFNA1, TLR2, and POSTN were considered to play a more critical role in malignant papillary thyroid cancer or immature cancer of papillary thyroid cancer. Our analysis concludes that TIMP1, LOX, CD276, IFNA1, TLR2, and POSTN are identified as thyroid cancer biomarkers, which lead to the different clinical courses of a woman older than 55 years old with papillary thyroid cancer. Especially CD276, POSTN, and IFNA1 may be considered as new biomarkers associated with the prognosis of thyroid cancer.

CONCLUSIONS

TIMP1, LOX, CD276, IFNA1, TLR2, and POSTN have different expressions in PTCs, which lead to the various clinical courses of a woman older than 55 years old with papillary thyroid cancer. Especially CD276, POSTN, and IFNA1 may be considered as new potential biomarkers associated with the prognosis of thyroid cancer. In addition, TF-miRNA-target gene regulatory network may help further reach for PTC.

Papillary Thyroid Carcinoma Landscape and Its Immunological Link With Hashimoto Thyroiditis at Single-Cell Resolution

The tumor microenvironment heterogeneity of papillary thyroid cancer (PTC) is poorly characterized. The relationship between PTC and Hashimoto thyroiditis (HT) is also in doubt. Here, we used single-cell RNA sequencing to map the transcriptome landscape of PTC from eight PTC patients, of which three were concurrent with HT. Predicted copy number variation in epithelial cells and mesenchymal cells revealed the distinct molecular signatures of carcinoma cells. Carcinoma cells demonstrated intertumoral heterogeneity based on BRAF V600E mutation or lymph node metastasis, and some altered genes were identified to be correlated with disease-free survival in The Cancer Genome Atlas datasets. In addition, transcription factor regulons of follicular epithelial cells unveil the different transcription activation state in PTC patients with or without concurrent HT. The immune cells in tumors exhibited distinct transcriptional states, and the presence of tumor-infiltrating B lymphocytes was predominantly linked to concurrent HT origin. Trajectory analysis of B cells and plasma cells suggested their migration potential from HT adjacent tissues to tumor tissues. Furthermore, we revealed diverse ligand–receptor pairs between non-immune cells, infiltrating myeloid cells, and lymphocytes. Our results provided a single-cell landscape of human PTC. These data would deepen the understanding of PTC, as well as the immunological link between PTC and HT.

In-depth proteomic profiling captures subtype-specific features of craniopharyngiomas

Craniopharyngiomas are rare epithelial tumors derived from pituitary gland embryonic tissue. This epithelial tumor can be categorized as an adamantinomatous craniopharyngioma (ACP) or papillary craniopharyngioma (PCP) subtype with histopathological and genetic differences. Genomic and transcriptomic profiles of craniopharyngiomas have been investigated; however, the proteomic profile has yet to be elucidated and added to these profiles. Recent improvements in high-throughput quantitative proteomic approaches have introduced new opportunities for a better understanding of these diseases and the efficient discovery of biomarkers. We aimed to confirm subtype-associated proteomic changes between ACP and PCP specimens. We performed a system-level proteomic study using an integrated approach that combines mass spectrometry-based quantitative proteomic, statistical, and bioinformatics analyses. The bioinformatics analysis showed that differentially expressed proteins between ACP and PCP were significantly involved in mitochondrial organization, fatty acid metabolic processes, exocytosis, the inflammatory response, the cell cycle, RNA splicing, cell migration, and neuron development. Furthermore, using network analysis, we identified hub proteins that were positively correlated with ACP and PCP phenotypes. Our findings improve our understanding of the pathogenesis of craniopharyngiomas and provide novel insights that may ultimately translate to the development of craniopharyngioma subtype-specific therapeutics.

Deregulation of HLA-I in cancer and its central importance for immunotherapy

It is now well accepted that many tumors undergo a process of clonal selection which means that tumor antigens arising at various stages of tumor progression are likely to be represented in just a subset of tumor cells. This process is thought to be driven by constant immunosurveillance which applies selective pressure by eliminating tumor cells expressing antigens that are recognized by T cells. It is becoming increasingly clear that the same selective pressure may also select for tumor cells that evade immune detection by acquiring deficiencies in their human leucocyte antigen (HLA) presentation pathways, allowing important tumor antigens to persist within cells undetected by the immune system. Deficiencies in antigen presentation pathway can arise by a variety of mechanisms, including genetic and epigenetic changes, and functional antigen presentation is a hard phenomenon to assess using our standard analytical techniques. Nevertheless, it is likely to have profound clinical significance and could well define whether an individual patient will respond to a particular type of therapy or not. In this review we consider the mechanisms by which HLA function may be lost in clinical disease, we assess the implications for current immunotherapy approaches using checkpoint inhibitors and examine the prognostic impact of HLA loss demonstrated in clinical trials so far. Finally, we propose strategies that might be explored for possible patient stratification.

The Potential Role of IL1RAP on Tumor Microenvironment-Related Inflammatory Factors in Stomach Adenocarcinoma

This study was performed to investigate the role of interleukin-1 receptor accessory protein (IL1RAP) in stomach carcinoma in vitro and in vivo, determine whether IL1RAP knockdown could regulate the development of stomach carcinoma, and elucidate the relationship between IL1RAP knockdown and inflammation by tumor microenvironment-related inflammatory factors in stomach carcinoma. We first used TCGA and GEPIA systems to predict the potential function of IL1RAP. Second, western blot and RT-PCR were used to analyze the expression, or mRNA level, of IL1RAP at different tissue or cell lines. Third, the occurrence and development of stomach carcinoma in vitro and in vivo were observed by using IL1RAP knockdown lentivirus. Finally, the inflammation of stomach carcinoma in vitro and in vivo was observed. Results show that in GEPIA and TCGA systems, IL1RAP expression in STAD tumor tissue was higher than normal, and high expression of IL1RAP in STAD patients had a worse prognostic outcome. Besides, GSEA shown IL1RAP was negative correlation of apopopsis, TLR4 and NF-κB signaling pathway. We also predicted that IL1RAP may related to IL-1 s, IL-33, and IL-36 s in STAD. The IL1RAP expression and mRNA level in tumor, or MGC803, cells were increased. Furthermore, IL1RAP knockdown by lentivirus could inhibit stomach carcinoma development in vitro and in vivo through weakening tumor cell proliferation, migration, invasion, therefore reducing tumor volume, weight, and biomarker levels, and increasing apoptotic level. Finally, we found IL1RAP knockdown could increase inflammation of tumor microenvironment-related inflammatory factors of stomach carcinoma, in vitro and in vivo. Our study demonstrates that IL1RAP is possibly able to regulate inflammation and apoptosis in stomach carcinoma. Furthermore, TLR4, NF-κB, IL-1 s, IL-33, and IL-36 s maybe the downstream target factor of IL1RAP in inflammation. These results may provide a new strategy for stomach carcinoma development by regulating inflammation.

Somatic Mutation Profiling of Papillary Thyroid Carcinomas by Whole-exome Sequencing and Its Relationship with Clinical Characteristics

The incidence of papillary thyroid carcinomas (PTCs) has increased rapidly during the past several decades. Until now, the mechanisms underlying the tumorigenesis of PTCs have remained largely unknown. Next-generation-sequencing (NGS) provides new ways to investigate the molecular pathogenesis of PTCs. To characterize the somatic alterations associated with PTCs, we performed whole-exome sequencing (WES) of PTCs from 23 Chinese patients. This study revealed somatic mutations in genes with relevant functions for tumorigenesis, such as BRAF, BCR, CREB3L2, DNMT1, IRS2, MSH6, and TP53. We also identified novel somatic gene alterations which may be potentially involved in PTC progression. Gene set enrichment analysis revealed that the cellular response to hormone stimulus, epigenetic modifications, such as protein/histone methylation and protein alkylation, as well as MAPK, PI3K-AKT, and FoxO/mTOR signaling pathways, were significantly altered in the PTCs studied here. Moreover, Protein-Protein Interaction (PPI) network analysis of our mutated gene selection highlighted EP300, KRAS, PTEN, and TP53 as major core genes. The correlation between gene mutations and clinicopathologic features of the PTCs defined by conventional ultrasonography (US) and contrast-enhanced ultrasonography (CEUS) were assessed. These analyses established significant associations between subgroups of mutations and respectively taller-than-wide, calcified, and peak time iso- or hypo-enhanced and metastatic PTCs. In conclusion, our study supplements the genomic landscape of PTCs and identifies new actionable target candidates and clinicopathology-associated mutations. Extension of this study to larger cohorts will help define comprehensive genomic aberrations in PTCs and validate target candidates. These new targets may open methods of individualized treatments adapted to the clinicopathologic specifics of the patients.

Malignant struma ovarii: next-generation sequencing of six cases revealed Nras, Braf, and Jak3 mutations

PURPOSE

Struma ovarii (SO) is a highly specialized ovarian teratoma, consisting of thyroid tissue. Rarely, carcinomas histologically identical to their thyroid counterparts may occur, and are comprehensively defined as malignant struma ovarii (MSO). Their optimal management is controversial, and the molecular profile of the malignant counterpart in the ovary is incompletely known. In this study, the clinicopathological and molecular features of six MSO from different Italian Institutions were analysed, to explore genetic profiles of potential therapeutic interest.

METHODS

The histopathological features and immunoprofile (according to the known markers Galectin-3, HBME1, cytokeratin 19 and CD56) were reviewed. In addition, all cases underwent genetic analysis with a next-generation sequencing (NGS) hot spot cancer panel detecting mutations in 50 genes involved in cancerogenesis. RET/PTC rearrangements and TERT promoter alterations were also evaluated.

RESULTS

Papillary carcinoma in all similar to its thyroid counterpart was found in five of six cases, including classical (two tumors) and follicular variant (three tumors) types. The last case was a poorly differentiated carcinoma. An activating gene mutation, was detected in five of six cases, including two NRAS, two BRAF, and one JAK3 oncogene mutations. No alterations were found in the other panel genes, nor in TERT promoter, or in RET chromosomal regions.

CONCLUSIONS

MSO is a rare condition. Papillary carcinoma is the predominant malignant type, sharing both histomorphological and molecular features of its thyroid counterpart. Interestingly, the single case of poorly differentiated carcinoma displayed a JAK3 mutation. The presence of such driving mutation could be of potential interest in guiding postoperative treatment.

Bioinformatics analysis identified shared differentially expressed genes as potential biomarkers for Hashimoto's thyroiditis-related papillary thyroid cancer

Background: Although the etiology of Hashimoto's thyroiditis (HT), a common autoimmune endocrine disease, is unknown, studies suggest a potential association with genetic factors and environmental conditions inducing excessive iodine intake. Additionally, HT patients have a high risk of papillary thyroid cancer (PTC), which is probably related to the chronic inflammation and autoimmune pathologic process occurring in HT, as it is thought to be associated with neoplastic transformation. Methods: Bioinformatics approaches can identify differentially expressed genes (DEGs) and analyze DEG functions in diseases. R software was used in this study to identify DEGs in HT and PTC using data in Gene Expression Omnibus (GEO). The online tools DAVID, Reactome, and AmiGO were employed for annotation, visualization, and integration of DEGs related to HT and PTC, and the STRING database and Cytoscape software were applied to predict and visualize protein-protein networks (PPIs) for DEG-encoded proteins. Coexpressed DEGs in HT and PTC were validated by reverse transcription PCR (RT-PCR). Results: In total, 326, 231, and 210 DEGs in HT specimens and samples of central PTC and PTC invasive areas, respectively, were detected. According to the PPI network, PTPN6, HLA-A, C3AR1, LCK and ITGB2 are hub genes among HT-DEGs, whereas FN1, CDH2, SERPINA1, and CYR61 are PTC-DEG hub genes. The shared DEGs LTF and CCL21 were validated by RT-PCR. Both bioinformatics and RT-PCR analyses showed LTF and CCL21 to be upregulated in HT tissues and downregulated in PTC tissues. Conclusions: We identified that expression of LTF and CCL21 are significantly different in HT and PTC, suggesting an underlying association between HT and PTC.

Emerging BRAF Mutations in Cancer Progression and Their Possible Effects on Transcriptional Networks

Gene mutations can induce cellular alteration and malignant transformation. Development of many types of cancer is associated with mutations in the B-raf proto-oncogene (BRAF) gene. The encoded protein is a component of the mitogen-activated protein kinases/extracellular signal-regulated kinases (MAPK/ERK) signaling pathway, transmitting information from the outside to the cell nucleus. The main function of the MAPK/ERK pathway is to regulate cell growth, migration, and proliferation. The most common mutations in the BRAF gene encode the V600E mutant (class I), which causes continuous activation and signal transduction, regardless of external stimulus. Consequently, cell proliferation and invasion are enhanced in cancer patients with such mutations. The V600E mutation has been linked to melanoma, colorectal cancer, multiple myeloma, and other types of cancers. Importantly, emerging evidence has recently indicated that new types of mutations (classes II and III) also play a paramount role in the development of cancer. In this minireview, we discuss the influence of various BRAF mutations in cancer, including aberrant transcriptional gene regulation in the affected tissues.

Immune Cell Confrontation in the Papillary Thyroid Carcinoma Microenvironment

Background

Papillary thyroid cancer has been associated with chronic inflammation. A systematic understanding of immune cell infiltration in PTC is essential for subsequent immune research and new diagnostic and therapeutic strategies.

Methods

Three different algorithms, single-sample gene set enrichment analysis (ssGSEA), immune cell marker and CIBERSORT, were used to evaluate immune cell infiltration levels (abundance and proportion) in 10 data sets (The Cancer Genome Atlas [TCGA], GSE3467, GSE3678, GSE5364, GSE27155, GSE33630, GSE50901, GSE53157, GSE58545, and GSE60542; a total of 799 PTC and 194 normal thyroid samples). Consensus unsupervised clustering divided PTC patients into low-immunity and high-immunity groups. Weighted gene coexpression network analysis (WGCNA) and gene set enrichment analysis (GSEA) were used to analyze the potential mechanisms causing differences in the immune response.

Results

Compared with normal tissues, PTC tissues had a higher overall immune level and higher abundance levels and proportions of M2 macrophages, Tregs, monocytes, neutrophils, dendritic cells (DCs), mast cells (MCs), and M0 macrophages. Compared with early PTC, advanced PTC showed higher immune infiltration and higher abundance levels and proportions of M2 macrophages, Tregs, monocytes, neutrophils, DCs, MCs, and M0 macrophages. Compared to the low-immunity group, the high-immunity group exhibited more advanced stages, larger tumor sizes, greater lymph node metastases, higher tall-cell PTCs, lower follicular PTC proportions, more BRAF mutations, and fewer RAS mutations. Epstein-Barr virus (EBV) infection was the most significantly enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway for key module genes.

Conclusions

In human PTC, M2 macrophages, Tregs, monocytes, neutrophils, DCs, MCs, and M0 macrophages appear to play a tumor-promoting role, while M1 macrophages, CD8+ T cells, B cells, NK cells, and T follicular helper (TFH) cells (including eosinophils, γδ T cells, and Th17 cells with weak supporting evidence) appear to play an antitumor role. During the occurrence and development of PTC, the overall immune level was increased, and the abundance and proportion of tumor-promoting immune cells were significantly increased, indicating that immune escape had been aggravated. Finally, we speculate that EBV may play an important role in changing the immune microenvironment of PTC tumors.

NanoString in the screening of genetic abnormalities associated with thyroid cancer

In the setting of cancer pathology, molecular characterization of tumors providing diagnostic and predictive information is acquiring more and more relevance. Moreover, the advent of innovative technologies continuously improves the knowledge of the molecular landscape of tumors and strengthens the links between clinics, tumor pathology and molecular features. In the clinical management of patients with thyroid nodules and thyroid tumors, the aid of molecular testing is encouraged but still not strongly recommended by current guidelines. Also for this reason this field of study is attracting much interest. The nCounter system is a relatively new technology based on a direct hybridization of fluorescent probes to specific nucleic acid targets, followed by digital measurement of signals; the reaction is highly multiplexable and results are robust and reproducible. This review reports and discusses the available data related to the application of this specific technique to thyroid nodules and thyroid tumors samples. The available data indicate that nCounter system represents a solid approach for the research of relevant diagnostic and prognostic biomarkers in thyroid pathology.

Differences in Gene Expression Profile of Primary Tumors in Metastatic and Non-Metastatic Papillary Thyroid Carcinoma-Do They Exist?

Molecular mechanisms of distant metastases (M1) in papillary thyroid cancer (PTC) are poorly understood. We attempted to analyze the gene expression profile in PTC primary tumors to seek the genes associated with M1 status and characterize their molecular function. One hundred and twenty-three patients, including 36 M1 cases, were subjected to transcriptome oligonucleotide microarray analyses: (set A-U133, set B-HG 1.0 ST) at transcript and gene group level (limma, gene set enrichment analysis (GSEA)). An additional independent set of 63 PTCs, including 9 M1 cases, was used to validate results by qPCR. The analysis on dataset A detected eleven transcripts showing significant differences in expression between metastatic and non-metastatic PTC. These genes were validated on microarray dataset B. The differential expression was positively confirmed for only two genes: IGFBP3, (most significant) and ECM1. However, when analyzed on an independent dataset by qPCR, the IGFBP3 gene showed no differences in expression. Gene group analysis showed differences mainly among immune-related transcripts, indicating the potential influence of tumor immune infiltration or signal within the primary tumor. The differences in gene expression profile between metastatic and non-metastatic PTC, if they exist, are subtle and potentially detectable only in large datasets.

Expression based biomarkers and models to classify early and late-stage samples of Papillary Thyroid Carcinoma

INTRODUCTION

Recently, the rise in the incidences of thyroid cancer worldwide renders it to be the sixth most common cancer among women. Commonly, Fine Needle Aspiration biopsy predominantly facilitates the diagnosis of the nature of thyroid nodules. However, it is inconsiderable in determining the tumor's state, i.e., benign or malignant. This study aims to identify the key RNA transcripts that can segregate the early and late-stage samples of Thyroid Carcinoma (THCA) using RNA expression profiles.

MATERIALS AND METHODS

In this study, we used the THCA RNA-Seq dataset of The Cancer Genome Atlas, consisting of 500 cancer and 58 normal (adjacent non-tumorous) samples obtained from the Genomics Data Commons (GDC) data portal. This dataset was dissected to identify key RNA expression features using various feature selection techniques. Subsequently, samples were classified based on selected features employing different machine learning algorithms.

RESULTS

Single gene ranking based on the Area Under the Receiver Operating Characteristics (AUROC) curve identified the DCN transcript that can classify the early-stage samples from late-stage samples with 0.66 AUROC. To further improve the performance, we identified a panel of 36 RNA transcripts that achieved F1 score of 0.75 with 0.73 AUROC (95% CI: 0.62-0.84) on the validation dataset. Moreover, prediction models based on 18-features from this panel correctly predicted 75% of the samples of the external validation dataset. In addition, the multiclass model classified normal, early, and late-stage samples with AUROC of 0.95 (95% CI: 0.84-1), 0.76 (95% CI: 0.66-0.85) and 0.72 (95% CI: 0.61-0.83) on the validation dataset. Besides, a five protein-coding transcripts panel was also recognized, which segregated cancer and normal samples in the validation dataset with F1 score of 0.97 and 0.99 AUROC (95% CI: 0.91-1).

CONCLUSION

We identified 36 important RNA transcripts whose expression segregated early and late-stage samples with reasonable accuracy. The models and dataset used in this study are available from the webserver CancerTSP (http://webs.iiitd.edu.in/raghava/cancertsp/).

The multifaceted anti-cancer effects of BRAF-inhibitors

The BRAF gene is commonly involved in normal processes of cell growth and differentiation. The BRAF (V600E) mutation is found in several human cancer, causing an increase of cell proliferation due to a modification of the ERK/MAPK-signal cascade. In particular, BRAFV600E mutation is found in those melanoma or thyroid cancer refractory to the common therapy and with a more aggressive phenotype. BRAF V600E was found to influence the composition of the so-called tumour microenvironment modulating both solid (immune-cell infiltration) and soluble (chemokines) mediators, which balance characterize the ultimate behaviour of the tumour, making it more or less aggressive. In particular, the presence of BRAFV600E mutation would be associated with a change of this balance to a more aggressive phenotype of the tumour and a worse prognosis. The investigation of the possible modulation of those components of tumour microenvironment is nowadays object of several studies as a new potential target therapy in those more complicated cases. At present several clinical trials both in melanoma and thyroid cancer are using BRAF-inhibitors with encouraging results, which are derived also from numerous in vitro pre-clinical studies aimed at evaluate the possible modulation of immune-cell density and of specific pro-tumorigenic chemokine secretion (CXCL8 and CCL2) by several BRAF-inhibitors in the context of melanoma and thyroid cancer. This review will encompass in vitro and in vivo studies which investigated the modulation of the tumour microenvironment by BRAF-inhibitors, highlighting also the most recent clinical trials with a specific focus on melanoma and thyroid cancer.

Coding Molecular Determinants of Thyroid Cancer Development and Progression

Thyroid cancer is the most common endocrine malignancy. Its incidence and mortality rates have increased for patients with advanced-stage papillary thyroid cancer. The characterization of the molecular pathways essential in thyroid cancer initiation and progression has made huge progress, underlining the role of intracellular signaling to promote clonal evolution, dedifferentiation, metastasis, and drug resistance. The discovery of genetic alterations that include mutations (BRAF, hTERT), translocations, deletions (eg, 9p), and copy-number gain (eg, 1q) has provided new biological insights with clinical applications. Understanding how molecular pathways interplay is one of the key strategies to develop new therapeutic treatments and improve prognosis.

Novel TG-FGFR1 and TRIM33-NTRK1 transcript fusions in papillary thyroid carcinoma

Papillary thyroid carcinoma (PTC) is most common among all thyroid cancers. Multiple genomic alterations occur in PTC, and gene rearrangements are one of them. Here we screened 14 tumors for novel fusion transcripts by RNA‐Seq. Two samples harboring RET/PTC1 and RET/PTC3 rearrangements were positive controls whereas the remaining ones were negative regarding the common PTC alterations. We used Sanger sequencing to validate potential fusions. We detected 2 novel potentially oncogenic transcript fusions: TG‐FGFR1 and TRIM33‐NTRK1. We detected 4 novel fusion transcripts of unknown significance accompanying the TRIM33‐NTRK1 fusion: ZSWIM5‐TP53BP2, TAF4B‐WDR1, ABI2‐MTA3, and ARID1B‐PSMA1. Apart from confirming the presence of RET/PTC1 and RET/PTC3 in positive control samples, we also detected known oncogenic fusion transcripts in remaining samples: TFG‐NTRK1, ETV6‐NTRK3, MKRN1‐BRAF, EML4‐ALK, and novel isoform of CCDC6‐RET.

Comprehensive Genetic Characterization of Human Thyroid Cancer Cell Lines: A Validated Panel for Preclinical Studies

PURPOSE

Thyroid cancer cell lines are valuable models but have been neglected in pancancer genomic studies. Moreover, their misidentification has been a significant problem. We aim to provide a validated dataset for thyroid cancer researchers.

EXPERIMENTAL DESIGN

We performed next-generation sequencing (NGS) and analyzed the transcriptome of 60 authenticated thyroid cell lines and compared our findings with the known genomic defects in human thyroid cancers.

RESULTS

Unsupervised transcriptomic analysis showed that 94% of thyroid cell lines clustered distinctly from other lineages. Thyroid cancer cell line mutations recapitulate those found in primary tumors (e.g., BRAF, RAS, or gene fusions). Mutations in the TERT promoter (83%) and TP53 (71%) were highly prevalent. There were frequent alterations in PTEN, PIK3CA, and of members of the SWI/SNF chromatin remodeling complex, mismatch repair, cell-cycle checkpoint, and histone methyl- and acetyltransferase functional groups. Copy number alterations (CNA) were more prevalent in cell lines derived from advanced versus differentiated cancers, as reported in primary tumors, although the precise CNAs were only partially recapitulated. Transcriptomic analysis showed that all cell lines were profoundly dedifferentiated, regardless of their derivation, making them good models for advanced disease. However, they maintained the BRAF^V600E versus RAS-dependent consequences on MAPK transcriptional output, which correlated with differential sensitivity to MEK inhibitors. Paired primary tumor-cell line samples showed high concordance of mutations. Complete loss of p53 function in TP53 heterozygous tumors was the most prominent event selected during in vitro immortalization.

CONCLUSIONS

This cell line resource will help inform future preclinical studies exploring tumor-specific dependencies.

Detection of BRAF V600E mutation in fine-needle aspiration fluid of papillary thyroid carcinoma by droplet digital PCR

OBJECTIVES

Papillary thyroid carcinoma (PTC) accounts for 85% of thyroid carcinoma, which is the most common endocrine tumor. For the diagnosis of PTC, ultrasound-guided fine needle aspiration (FNA) with pathological evaluation is the standard test and BRAF V600E mutation is the most common molecular marker associated with the occurrence, progression and poor clinicopathological characteristics of PTC. However, because of the small amount of the tumor cells obtained by FNA for pathological evaluation or BRAF V600E mutation detection, more sensitive and accurate methods are required. Our study aimed to investigate the performance of droplet digital PCR (ddPCR) in detecting BRAF V600E mutation in FNA samples from PTC patients.

METHODS

One hundred and sixty suspected thyroid cancer patients were enrolled, including 146 PTC patients, 2 follicular thyroid carcinoma (FTC) and 12 benign patients, identified by FNA biopsy according to the NCCN clinical practice guidelines of Thyroid Carcinoma. ddPCR and amplification-refractory mutation system (ARMS, AmoyDx) were used to detect BRAFV600E mutation and the results were compared.

RESULTS

ddPCR had high reproducibility (CV0.1% = 22.82% and CV10% = 4.85%) and the detection sensitivity can reach 1–2 copies/μl (0.01%). Among the 160 patients, 128 BRAF V600E mutations were detected, including 4 ARMS negative patients and 3 benign cases [corrected].

CONCLUSIONS

Our results demonstrated that ddPCR could be used in detecting BRAF V600E mutation from FNA fluid samples with higher sensitivity and accuracy than ARMS.

Network Analyses of Integrated Differentially Expressed Genes in Papillary Thyroid Carcinoma to Identify Characteristic Genes

Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer. Identifying characteristic genes of PTC are of great importance to reveal its potential genetic mechanisms. In this paper, we proposed a framework, as well as a measure named Normalized Centrality Measure (NCM), to identify characteristic genes of PTC. The framework consisted of four steps. First, both up-regulated genes and down-regulated genes, collectively called differentially expressed genes (DEGs), were screened and integrated together from four datasets, that is, GSE3467, GSE3678, GSE33630, and GSE58545; second, an interaction network of DEGs was constructed, where each node represented a gene and each edge represented an interaction between linking nodes; third, both traditional measures and the NCM measure were used to analyze the topological properties of each node in the network. Compared with traditional measures, more genes related to PTC were identified by the NCM measure; fourth, by mining the high-density subgraphs of this network and performing Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, several meaningful results were captured, most of which were demonstrated to be associated with PTC. The experimental results proved that this network framework and the NCM measure are useful for identifying more characteristic genes of PTC.

Immune Gene Signature Delineates a Subclass of Papillary Thyroid Cancer with Unfavorable Clinical Outcomes

Papillary thyroid carcinoma (PTC) represents a heterogeneous disease with diverse clinical outcomes highlighting a need to identify robust biomarkers with clinical relevance. We applied non-negative matrix factorization-based deconvolution to publicly available gene expression profiles of thyroid cancers in the Cancer Genome Atlas (TCGA) consortium. Among three metagene signatures identified, two signatures were enriched in canonical BRAF-like and RAS-like thyroid cancers with up-regulation of genes involved in oxidative phosphorylation and cell adhesions, respectively. The third metagene signature representing up-regulation of immune-related genes further segregated BRAF-like and RAS-like PTCs into their respective subgroups of immunoreactive (IR) and immunodeficient (ID), respectively. BRAF-IR PTCs showed enrichment of tumor infiltrating immune cells, tall cell variant PTC, and shorter recurrence-free survival compared to BRAF-ID PTCs. RAS-IR and RAS-ID PTC subtypes included majority of normal thyroid tissues and follicular variant PTC, respectively. Immunopathological features of PTC subtypes such as immune cell fraction, repertoire of T cell receptors, cytolytic activity, and expression level of immune checkpoints such as and PD-L1 and CTLA-4 were consistently observed in two different cohorts. Taken together, an immune-related metagene signature can classify PTCs into four molecular subtypes, featuring the distinct histologic type, genetic and transcriptional alterations, and potential clinical significance.

The Transcription Factor ETV5 Mediates BRAFV600E-Induced Proliferation and TWIST1 Expression in Papillary Thyroid Cancer Cells

The ETS family of transcription factors is involved in several normal remodeling events and pathological processes including tumor progression. ETS transcription factors are divided into subfamilies based on the sequence and location of the ETS domain. ETV5 (Ets variant gene 5; also known as ERM) is a member of the PEA3 subfamily. Our meta-analysis of normal, benign, and malignant thyroid samples demonstrated that ETV5 expression is upregulated in papillary thyroid cancer and was predominantly associated with BRAF V600E or RAS mutations. However, the precise role of ETV5 in these lesions is unknown. In this study, we used the KTC1 cell line as a model for human advanced papillary thyroid cancer (PTC) because the cells harbor the heterozygous BRAF (V600E) mutation together with the C250T TERT promoter mutation. The role of ETV5 in PTC proliferation was tested using RNAi followed by high-throughput screening. Signaling pathways driving ETV5 expression were identified using specific pharmacological inhibitors. To determine if ETV5 influences the expression of epithelial-to-mesenchymal (EMT) markers in these cells, an EMT PCR array was used, and data were confirmed by qPCR and ChIP-qPCR. We found that ETV5 is critical for PTC cell growth, is expressed downstream of the MAPK pathway, and directly upregulates the transcription factor TWIST1, a known marker of intravasation and metastasis. Increased ETV5 expression could therefore be considered as a marker for advanced PTCs and a possible future therapeutic target.

Improving the value of public RNA-seq expression data by phenotype prediction

Publicly available genomic data are a valuable resource for studying normal human variation and disease, but these data are often not well labeled or annotated. The lack of phenotype information for public genomic data severely limits their utility for addressing targeted biological questions. We develop an in silico phenotyping approach for predicting critical missing annotation directly from genomic measurements using well-annotated genomic and phenotypic data produced by consortia like TCGA and GTEx as training data. We apply in silico phenotyping to a set of 70 000 RNA-seq samples we recently processed on a common pipeline as part of the recount2 project. We use gene expression data to build and evaluate predictors for both biological phenotypes (sex, tissue, sample source) and experimental conditions (sequencing strategy). We demonstrate how these predictions can be used to study cross-sample properties of public genomic data, select genomic projects with specific characteristics, and perform downstream analyses using predicted phenotypes. The methods to perform phenotype prediction are available in the phenopredict R package and the predictions for recount2 are available from the recount R package. With data and phenotype information available for 70,000 human samples, expression data is available for use on a scale that was not previously feasible.

Kallikreins Stepwise Scoring Reveals Three Subtypes of Papillary Thyroid Cancer with Prognostic Implications

BACKGROUND

Papillary thyroid cancer (PTC) is the most common type of thyroid cancer. Unlike most cancers, its incidence has dramatically increased in the last decades mainly due to increased diagnosis of indolent PTCs. Adequate risk stratification is crucial to avoid the over-treatment of low-risk patients, as well as the under-treatment of high-risk patients, but the currently available markers are still insufficient. Kallikreins (KLKs) are emergent biomarkers in cancer, but their involvement in PTC is unknown.

METHODS

This study analyzed DNA methylation (HumanMethylation arrays) and gene expression (RNA-Seq) of KLKs, BRAF and RAS mutations, and clinical data from four published thyroid cancer data sets including normal and tumor tissues (n = 73, n = 475, n = 20, and n = 82) as discovery, training, and validation series. The C4.5 classification algorithm was used to generate a decision tree. Disease-free survival was estimated using Kaplan-Meier and Cox approaches. Specific analyses were performed using real-time polymerase chain reaction and immunohistochemistry.

RESULTS

The entire KLK family was deregulated in PTC, displaying a specific epigenetic and transcriptional profile strongly associated with BRAF^V600E or RAS mutations. Thus, a decision-tree algorithm was developed based on three KLKs with >80% sensitivity and >95% specificity, identifying BRAF- and RAS-mutated tumors. Notably, tumors lacking these mutations were classified as BRAF- or RAS-like. Most importantly, the KLK algorithm uncovered a novel PTC subtype showing favorable prognostic features.

CONCLUSIONS

The KLK algorithm could lead to a new clinically applicable strategy with important implications for the risk stratification of PTC and the management of patients.

BRAFwild papillary thyroid carcinoma has two distinct mRNA expression patterns with different clinical behaviors

BACKGROUND

Using a large set of genomic data from The Cancer Genome Atlas (TCGA), we classified BRAF^wild papillary thyroid carcinomas (PTCs) into 2 subtypes with distinct molecular patterns and different clinical behaviors. We also suggested gene signatures (RAS-score) to predict molecular subtypes and clinical behaviors of BRAF^wild PTC.

METHOD

Integrated genomic analysis was done using all genomic data of PTC in TCGA data portal (https://tcga-data.nci.nih.gov) and cancer browser (https://genome-cancer.ucsc.edu). Using Gene Ontology and a logistic regression test, we selected gene signatures (RAS-score) and applied this prediction model to the validation cohort (GSE60542).

RESULT

When we performed multiplatform genomic analysis, BRAF^wild PTCs were divided into 2 molecular subtypes. Each subtype showed distinct molecular patterns and clinical behaviors. Gene signatures successfully predicted molecular subtype in another validation cohort.

CONCLUSION

We found that BRAF^wild PTCs were divided into 2 molecular subtypes and each subtype showed distinct molecular patterns, different activated pathways, and different clinical behaviors.

CRABP1, C1QL1 and LCN2 are biomarkers of differentiated thyroid carcinoma, and predict extrathyroidal extension

Background

The prognostic variability of thyroid carcinomas has led to the search for accurate biomarkers at the molecular level. Follicular thyroid carcinoma (FTC) is a typical example of differentiated thyroid carcinomas (DTC) in which challenges are faced in the differential diagnosis.

Methods

We used high-throughput paired-end RNA sequencing technology to study four cases of FTC with different degree of capsular invasion: two minimally invasive (mFTC) and two widely invasive FTC (wFTC). We searched by genes differentially expressed between mFTC and wFTC, in an attempt to find biomarkers of thyroid cancer diagnosis and/or progression. Selected biomarkers were validated by real-time quantitative PCR in 137 frozen thyroid samples and in an independent dataset (TCGA), evaluating the diagnostic and the prognostic performance of the candidate biomarkers.

Results

We identified 17 genes significantly differentially expressed between mFTC and wFTC. C1QL1, LCN2, CRABP1 and CILP were differentially expressed in DTC in comparison with normal thyroid tissues. LCN2 and CRABP1 were also differentially expressed in DTC when compared with follicular thyroid adenoma. Additionally, overexpression of LCN2 and C1QL1 were found to be independent predictors of extrathyroidal extension in DTC.

Conclusions

We conclude that the underexpression of CRABP1 and the overexpression of LCN2 may be useful diagnostic biomarkers in thyroid tumours with questionable malignity, and the overexpression of LCN2 and C1QL1 may be useful for prognostic purposes.

Electronic supplementary material

The online version of this article (10.1186/s12885-017-3948-3) contains supplementary material, which is available to authorized users.

ETV6-NTRK3 and STRN-ALK kinase fusions are recurrent events in papillary thyroid cancer of adult population

OBJECTIVE

PTC-specific analysis identified novel fusions involving RET, BRAF, NTRK1, NTRK3, AGK and ALK genes in adults and pediatric PTCs. Although many novel fusions are PTC-specific events and, therefore, are ideal for diagnosis purposes, validation across additional and larger patient cohorts is essential for introducing these potential diagnostic or prognostic biomarkers into the clinical practice. As most of the BRAF, NTRK3 and ALK fusions were initially found in pediatric PTC or in more aggressive thyroid carcinomas, and there is a great disparity across population, in this study, we screened a large set of adult-sporadic PTC cases for the most prevalent kinase fusion lately described in the TCGA.

DESIGN AND METHODS

The prevalence of the fusions was determined by RT-PCR in 71 classical PTC, 45 follicular variants of PTC (FVPTC), 19 follicular thyroid adenomas (FTAs) and 22 follicular thyroid carcinomas (FTCs).

RESULTS

ETV6-NTRK3 was exclusively found in FVPTC, in both encapsulated and infiltrative variants, but was not found in FTAs and FTCs. STRN-ALK was found in both classical PTC and FVPTC. No AGK-BRAF fusion was identified in this series, endorsing that AGK-BRAF is a genetic event mainly associated with pediatric PTCs.

CONCLUSIONS

The identification of kinase fusions in thyroid carcinomas helps to expand our knowledge about the landscape of oncogenic alterations in PTC. As ETV6-NTRK3 and STRN-ALK are recurrent and not identified in benign lesions, they can certainly help with diagnosis of thyroid nodules. Further analysis is needed to define if they can also be useful for prognosis and guiding therapy.

Investigating the mechanisms of papillary thyroid carcinoma using transcriptome analysis

As the predominant thyroid cancer, papillary thyroid cancer (PTC) accounts for 75–85% of thyroid cancer cases. This research aimed to investigate transcriptomic changes and key genes in PTC. Using RNA-sequencing technology, the transcriptional profiles of 5 thyroid tumor tissues and 5 adjacent normal tissues were obtained. The single nucleotide polymorphisms (SNPs) were identified by SAMtools software and then annotated by ANNOVAR software. After differentially expressed genes (DEGs) were selected by edgR software, they were further investigated by enrichment analysis, protein domain analysis, and protein-protein interaction (PPI) network analysis. Additionally, the potential gene fusion events were predicted using FusionMap software. A total of 70,172 SNPs and 2,686 DEGs in the tumor tissues, as well as 83,869 SNPs in the normal tissues were identified. In the PPI network, fibronectin 1 (FN1; degree=31) and transforming growth factor β receptor 1 (TGFβR1; degree=22) had higher degrees. A total of 7 PPI pairs containing the non-synonymous risk SNP loci in the interaction domains were identified. Particularly, the interaction domains involved in the interactions of FN1 and 5 other proteins (such as FN1-tenascin C, TNC) had non-synonymous risk SNP loci. Furthermore, 11 and 4 gene fusion events were identified in all of the tumor tissues and normal tissues, respectively. Additionally, the NK2 homeobox 1-surfactant associated 3 (NKX2-1-SFTA3) gene fusion was identified in both tumor and normal tissues. These results indicated that TGFβR1 and the NKX2-1-SFTA3 gene fusion may be involved in PTC. Furthermore, FN1 and TNC containing the non-synonymous risk SNP loci might serve a role in PTC by interacting with each other.

Overexpression of teneurin transmembrane protein 1 is a potential marker of disease progression in papillary thyroid carcinoma

Although papillary thyroid cancer is a relatively indolent malignancy, its progression may be associated with dedifferentiation and resistance to radioactive iodine treatment. In this study, patterns of differentially expressed genes in association with disease progression were systemically evaluated. We firstly performed transcriptome analyses for four matched cancerous and noncancerous tissue pairs of the classical subtype of papillary thyroid cancer. Among the upregulated and downregulated genes, the expression of 164 and 183 genes increased and decreased, respectively, from stage I to stage IV. Functional enrichment and pathway analysis showed that angiogenesis pathway was upregulated, whereas oxidation-reduction and metabolism of reactive oxygen species were downregulated. Teneurin transmembrane protein 1 (TENM1) expression was highly upregulated in cancerous tissues and negative in benign thyroid tissues. By immunohistochemistry, TENM1 expression in papillary thyroid cancer was associated with the classical subtype (p = 0.018), extrathyroidal invasion (p = 0.001), BRAF V600E mutation (p < 0.001), and an advanced stage (p = 0.019). Taken together, our results indicate that distinct pathways are involved in papillary thyroid cancer progression, and TENM1 is a potential marker of cancer progression.

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

BACKGROUND

In recent years, long non-coding RNAs (lncRNAs) have been shown to play a critical regulatory role in cancer biology. However, the contribution of lncRNAs to papillary thyroid cancer (PTC) remains largely unknown.

METHODS

RNA sequencing and quantitative reverse transcription polymerase chain reaction were used to detect and verify changes to the transcriptome profile in 12 PTC tissues compared to paired normal adjacent tissues. The statistical correlation between differentially expressed lncRNAs and clinicopathologic characteristics was analyzed, and potential lncRNA functions were predicted by examining annotations for the co-expressed mRNAs. Furthermore, the specific subgroup patterns of the PTC transcriptome remodeled by BRAF mutations were also analyzed.

RESULTS

A total of 188 lncRNAs and 505 mRNAs were differentially expressed in 50% or more of the PTC tissues (fold change >2; p < 0.05) as assessed by RNA-sequencing compared with paired normal adjacent tissues. Forty-seven lncRNAs and 39 mRNAs were verified in 31 pairs of PTC specimens using quantitative reverse transcription polymerase chain reaction, and the results were consistent with the RNA sequencing data. The lncRNAs NONHSAT076747 and NONHSAT122730 were associated with lymph node metastasis, and NONHSAG051968 expression was negatively correlated with tumor size. A co-expression network between differentially expressed lncRNAs and protein-coding RNAs was constructed and analyzed, and functional analysis suggested that the differentially expressed genes mainly participate in ECM-receptor interactions and the focal adhesion pathway. Furthermore, a specific PTC transcriptome subtype pattern stratified by BRAF mutation was also uncovered. The p53 signaling pathway was the most highly enriched pathway among the BRAF mutation-related genes.

CONCLUSIONS

This study reveals specific changes to the lncRNA profile associated with PTC, and provides new insight into its pathogenesis. The PTC-associated lncRNAs NONHSAG051968, NONHSAT076747, and NONHSAT122730 might be potential diagnostic and therapeutic targets for PTC patients, and lncRNAs with subtype-specific expression stratified by BRAF mutation might be significant in individual molecular subtypes.

Next-generation sequencing in thyroid cancer

Next-generation sequencing (NGS) in thyroid cancer allows for simultaneous high-throughput sequencing analysis of variable genetic alterations and provides a comprehensive understanding of tumor biology. In thyroid cancer, NGS offers diagnostic improvements for fine needle aspiration (FNA) cytology of thyroid with indeterminate features. It also contributes to patient management, providing risk stratification of patients based on the risk of malignancy. Furthermore, NGS has been adopted in cancer research. It is used in molecular tumor classification, and molecular prediction of recurrence and metastasis in papillary thyroid carcinoma. This review covers previous NGS analyses in variable types of thyroid cancer, where samples including FNA cytology, fresh frozen tissue, and formalin-fixed, paraffin-embedded tissues were used. This review also focuses on the clinical and research implications of using NGS to study and treat thyroid cancer.

Fusion proteins in head and neck neoplasms: Clinical implications, genetics, and future directions for targeting

Fusion proteins resulting from chromosomal rearrangements are known to drive the pathogenesis of a variety of hematological and solid neoplasms such as chronic myeloid leukemia and non-small-cell lung cancer. Efforts to elucidate the role they play in these malignancies have led to important diagnostic and therapeutic triumphs, including the famous development of the tyrosine kinase inhibitor dasatinib targeting the BCR-ABL fusion. Until recently, there has been a paucity of research investigating fusion proteins harbored by head and neck neoplasms. The discovery and characterization of novel fusion proteins in neoplasms originating from the thyroid, nasopharynx, salivary glands, and midline head and neck structures offer substantial contributions to our understanding of the pathogenesis and biological behavior of these neoplasms, while raising new therapeutic and diagnostic opportunities. Further characterization of these fusion proteins promises to facilitate advances on par with those already achieved with regard to hematologic malignancies in the precise, molecularly guided diagnosis and treatment of head and neck neoplasms. The following is a subsite specific review of the clinical implications of fusion proteins in head and neck neoplasms and the future potential for diagnostic targeting.

Recurrent hormone-binding domain truncated ESR1 amplifications in primary endometrial cancers suggest their implication in hormone independent growth

The estrogen receptor alpha (ERα) is highly expressed in both endometrial and breast cancers, and represents the most prevalent therapeutic target in breast cancer. However, anti-estrogen therapy has not been shown to be effective in endometrial cancer. Recently it has been shown that hormone-binding domain alterations of ERα in breast cancer contribute to acquired resistance to anti-estrogen therapy. In analyses of genomic data from The Cancer Genome Atlas (TCGA), we observe that endometrial carcinomas manifest recurrent ESR1 gene amplifications that truncate the hormone-binding domain encoding region of ESR1 and are associated with reduced mRNA expression of exons encoding the hormone-binding domain. These findings support a role for hormone-binding alterations of ERα in primary endometrial cancer, with potentially important therapeutic implications.

The "next-generation" knowledge of papillary thyroid carcinoma

The application of Next-Generation Sequencing for studying the genetics of papillary thyroid carcinomas (PTC) has recently revealed new somatic mutations and gene fusions as potential new tumor-initiating events in patients without any known driver lesion. Gene and miRNA expression analyses defined clinically relevant subclasses correlated to tumor progression. In addition, it has been shown that tumor driver mutations in BRAF, and RET rearrangements - altogether termed "BRAF-like" carcinomas - have a very similar expression pattern and constitute a distinct category. Conversely, "RAS-like" carcinomas have a different genomic, epigenomic, and proteomic profile. These findings justify the need to reconsider PTC classification schemes.

Novel fusion transcripts in bladder cancer identified by RNA-seq

Urothelial carcinoma (UC) is the most common type of bladder cancer and is the second most frequently diagnosed genitourinary tumor. The identification of fusion genes in bladder cancer might provide new perspectives for its classification and significance. In this study, we present a thorough search on three UC samples for novel fusion transcripts in bladder cancer using high-throughput RNA sequencing. We used stringent requirements for 819 fusion candidates and nominated 10 candidate fusion transcripts. Among them four novel fusion genes SEPT9/CYHR, IGF1R/TTC23, SYT8/TNNI2 and CASZ1/DFFA were validated and characterized in 48 formalin-fixed paraffin-embedded (FFPE) specimens of bladder cancer. Chromosomal rearrangements of regions 17q25, 15q26.3 and 1p36.22 resulting in the fusion transcripts SEPT9/CYHR, IGF1R/TTC23 and CASZ1/DFFA, appeared to be rare or unique events because they were not detected in the 48 UC samples. In contrast, the SYT8/TNNI2 fusion transcript resulting from transcription-induced chimerism by read-through mechanisms was a rather common and tumor-specific event occurring in 37.5% (18/48) of the UC specimens. Further investigation of functional and clinical relevance of novel fusion genes remains to be elucidated to reveal their role in bladder carcinogenesis.

RNA sequencing identifies crucial genes in papillary thyroid carcinoma (PTC) progression

PURPOSE

The study aims to uncover molecular mechanisms of PTC (papillary thyroid carcinoma) progression and provide therapeutic biomarkers.

METHODS

The paired tumor and control tissues were obtained from 5 PTC patients. RNA was extracted and cDNA libraries were constructed. RNA-sequencing (RNA-seq) was performed on the Illumina HiSeq2000 platform using paired-end method. After preprocessing of the RNA-seq data, gene expression value was calculated by RPKM. Then the differentially expressed genes (DEGs) were identified with edgeR. Functional enrichment and protein-protein interaction (PPI) network analyses were conducted for the DEGs. Module analysis of the PPI network was also performed. Transcription factors (TFs) of DEGs were predicted.

RESULTS

A cohort of 496 up-regulated DEGs mainly correlating with the ECM degradation pathways, and 440 down-regulated DEGs predominantly enriching in transmembrane transport process were identified. Hub nodes in the PPI network were RRM2 and a set of collagens (COL1A1, COL3A1 and COL5A1), which were also remarkable in module 3 and module 5, respectively. Genes in module 3 were associated with cell cycle pathways, while in module 5 were related to ECM degradation pathways. PLAU, PSG1 and EGR2 were the crucial TFs with higher transcriptional activity in PTC than in control.

CONCLUSION

Several genes including COL1A1, COL3A1, RRM2, PLAU, and EGR2 might be used as biomarkers of PTC therapy. Among them, COL1A1 and COL3A1 might exert their functions via involving in ECM degradation pathway, while RRM2 through cell cycle pathway. PLAU might be an active TF, whereas EGR2 might be a tumor suppressor.

Deconstructing innate immune signaling in myelodysplastic syndromes

Overexpression of immune-related genes is widely reported in myelodysplastic syndromes (MDSs), and chronic immune stimulation increases the risk for developing MDS. Aberrant innate immune activation, such as that caused by increased toll-like receptor (TLR) signaling, in MDS can contribute to systemic effects on hematopoiesis, in addition to cell-intrinsic defects on hematopoietic stem/progenitor cell (HSPC) function. This review will deconstruct aberrant function of TLR signaling mediators within MDS HSPCs that may contribute to cell-intrinsic consequences on hematopoiesis and disease pathogenesis. We will discuss the contribution of chronic TLR signaling to the pathogenesis of MDS based on evidence from patients and mouse genetic models.

Intratumor heterogeneity and clonal evolution in an aggressive papillary thyroid cancer and matched metastases

The contribution of intratumor heterogeneity to thyroid metastatic cancers is still unknown. The clonal relationships between the primary thyroid tumors and lymph nodes (LN) or distant metastases are also poorly understood. The objective of this study was to determine the phylogenetic relationships between matched primary thyroid tumors and metastases. We searched for non-synonymous single-nucleotide variants (nsSNVs), gene fusions, alternative transcripts, and loss of heterozygosity (LOH) by paired-end massively parallel sequencing of cDNA (RNA-Seq) in a patient diagnosed with an aggressive papillary thyroid cancer (PTC). Seven tumor samples from a stage IVc PTC patient were analyzed by RNA-Seq: two areas from the primary tumor, four areas from two LN metastases, and one area from a pleural metastasis (PLM). A large panel of other thyroid tumors was used for Sanger sequencing screening. We identified seven new nsSNVs. Some of these were early events clonally present in both the primary PTC and the three matched metastases. Other nsSNVs were private to the primary tumor, the LN metastases and/or the PLM. Three new gene fusions were identified. A novel cancer-specific KAZN alternative transcript was detected in this aggressive PTC and in dozens of additional thyroid tumors. The PLM harbored an exclusive whole-chromosome 19 LOH. We have presented the first, to our knowledge, deep sequencing study comparing the mutational spectra in a PTC and both LN and distant metastases. This study has yielded novel findings concerning intra-tumor heterogeneity, clonal evolution and metastases dissemination in thyroid cancer.

MHC class I loss is a frequent mechanism of immune escape in papillary thyroid cancer that is reversed by interferon and selumetinib treatment in vitro

PURPOSE

To evaluate MHC class I expression on papillary thyroid cancer (PTC) and analyze changes in MHC expression and associated immune activation with current and experimental treatments for thyroid cancer using in vitro PTC cell lines.

EXPERIMENTAL DESIGN

MHC class I expression and assessment of tumor-infiltrating leukocyte populations were evaluated by immunohistochemistry. PTC cell lines were analyzed for HLA-ABC expression by flow cytometry following tyrosine kinase inhibitor, IFNα or IFNγ, or radiation treatment. Functional changes in antigenicity were assessed by coculture of allogeneic donor peripheral blood leukocytes (PBL) with pretreated or untreated PTC cell lines and measurement of T-cell activation and cytokine production.

RESULTS

Both MHC class I and β2-microglobulin expression was reduced or absent in 76% of PTC specimens and was associated with reduced tumor-infiltrating immune cells, including effector (CD3(+), CD8(+), CD16(+)) and suppressor (FoxP3(+)) populations. Treatment of PTC cell lines with the MEK1/2 inhibitor selumetinib or IFN increased HLA-ABC expression. This phenotypic change was associated with increased T-cell activation (%CD25(+) of CD3(+)) and IL2 production by PBL cocultured with treated PTC cell lines. Additive effects were seen with combination selumetinib and IFN treatment.

CONCLUSIONS

MHC class I expression loss is frequent in human PTC specimens and represents a significant mechanism of immune escape. Increased antigenicity following selumetinib and IFN treatment warrants further study for immunotherapy of progressive PTC.