Genome-wide methylation patterns in papillary thyroid cancer are distinct based on histological subtype and tumor genotype

A diagnostic model based on DNA methylation haplotype block characteristics for identifying papillary thyroid carcinoma from thyroid adenoma

Papillary thyroid carcinoma (PTC) is the most prevalent form of thyroid cancer. Methylation of some genes plays a crucial role in the tendency to malignancy as well as poor prognosis of thyroid cancer, suggesting that methylation features can serve as complementary markers for molecular diagnosis. In this study, we aimed to develop and validate a diagnostic model for PTC based on DNA methylation markers. A total of 142 thyroid nodule tissue samples containing 84 cases of PTC and 58 cases of thyroid adenoma (TA) were collected for reduced representation bisulfite sequencing (RRBS) and subsequent analysis. The diagnostic model was constructed by the logistic regression (LR) method followed by 5-cross validation and based on 94 tissue methylation haplotype block (MHB) markers. The model achieved an area under the receiver operating characteristic curve (AUROC) of 0.974 (95% CI, 0.964-0.981) on 108 training samples and 0.917 (95% CI, 0.864-0.973) on 27 independent testing samples. The diagnostic model scores showed significantly high in males (P = 0.0016), age ≤ 45 years (P = 0.026), high body mass index (BMI) (P = 0.040), lymph node metastasis (P = 0.00052) and larger nodules (P = 0.0017) in the PTC group, and the risk score of this diagnostic model showed significantly high in recurrent PTC group (P = 0.0005). These results suggest that the diagnostic model can be expected to be a powerful tool for PTC diagnosis and there are more potential clinical applications of methylation markers to be excavated.

Fasting regulates mitochondrial function through lncRNA PRKCQ-AS1-mediated IGF2BPs in papillary thyroid carcinoma

Recurring evidence suggests that fasting has extensive antitumor effects in various cancers, including papillary thyroid carcinoma (PTC). However, the underlying mechanism of this relationship with PTC is unknown. In this study, we study the effect of fasting on glycolysis and mitochondrial function in PTC. We find that fasting impairs glycolysis and reduces mitochondrial dysfunction in vitro and in vivo and also fasting in vitro and fasting mimicking diets (FMD) in vivo significantly increase the expression of lncRNA-protein kinase C theta antisense RNA 1 (PRKCQ-AS1), during the inhibition of TPC cell glycolysis and mitochondrial function. Moreover, lncRNA PRKCQ-AS1 was significantly lower in PTC tissues and cells. In addition, PRKCQ-AS1 overexpression increased PTC cell glycolysis and mitochondrial function; PRKCQ-AS1 knockdown has the opposite effect. On further mechanistic analysis, we identified that PRKCQ-AS1 physically interacts with IGF2BPs and enhances protein arginine methyltransferases 7 (PRMT7) mRNA, which is the key player in regulating glycolysis and mitochondrial function in PTC. Hence, PRKCQ-AS1 inhibits tumor growth while regulating glycolysis and mitochondrial functions via IGF2BPs/PRMT7 signaling. These results indicate that lncRNA PRKCQ-AS1 is a key downstream target of fasting and is involved in PTC metabolic reprogramming. Further, the PRKCQ-AS1/IGF2BPs/PRMT7 axis is an ideal therapeutic target for PTC diagnosis and treatment.

Advances and challenges in thyroid cancer: The interplay of genetic modulators, targeted therapies, and AI-driven approaches

Thyroid cancer continues to exhibit a rising incidence globally, predominantly affecting women. Despite stable mortality rates, the unique characteristics of thyroid carcinoma warrant a distinct approach. Differentiated thyroid cancer, comprising most cases, is effectively managed through standard treatments such as thyroidectomy and radioiodine therapy. However, rarer variants, including anaplastic thyroid carcinoma, necessitate specialized interventions, often employing targeted therapies. Although these drugs focus on symptom management, they are not curative. This review delves into the fundamental modulators of thyroid cancers, encompassing genetic, epigenetic, and non-coding RNA factors while exploring their intricate interplay and influence. Epigenetic modifications directly affect the expression of causal genes, while long non-coding RNAs impact the function and expression of micro-RNAs, culminating in tumorigenesis. Additionally, this article provides a concise overview of the advantages and disadvantages associated with pharmacological and non-pharmacological therapeutic interventions in thyroid cancer. Furthermore, with technological advancements, integrating modern software and computing into healthcare and medical practices has become increasingly prevalent. Artificial intelligence and machine learning techniques hold the potential to predict treatment outcomes, analyze data, and develop personalized therapeutic approaches catering to patient specificity. In thyroid cancer, cutting-edge machine learning and deep learning technologies analyze factors such as ultrasonography results for tumor textures and biopsy samples from fine needle aspirations, paving the way for a more accurate and effective therapeutic landscape in the near future.

Classification of Thyroid Tumors Based on DNA Methylation Patterns

Background: Alterations in DNA methylation are stable epigenetic events that can serve as clinical biomarkers. The aim of this study was to analyze methylation patterns among various follicular cell-derived thyroid neoplasms to identify disease subtypes and help understand and classify thyroid tumors. Methods: We employed an unsupervised machine learning method for class discovery to search for distinct methylation patterns among various thyroid neoplasms. Our algorithm was not provided with any clinical or pathological information, relying exclusively on DNA methylation data to classify samples. We analyzed 810 thyroid samples (n = 256 for discovery and n = 554 for validation), including benign and malignant tumors, as well as normal thyroid tissue. Results: Our unsupervised algorithm identified that samples could be classified into three subtypes based solely on their methylation profile. These methylation subtypes were strongly associated with histological diagnosis (p < 0.001) and were therefore named normal-like, follicular-like, and papillary thyroid carcinoma (PTC)-like. Follicular adenomas, follicular carcinomas, oncocytic adenomas, and oncocytic carcinomas clustered together forming the follicular-like methylation subtype. Conversely, classic papillary thyroid carcinomas (cPTC) and tall cell PTC clustered together forming the PTC-like subtype. These methylation subtypes were also strongly associated with genomic drivers: 98.7% BRAFV600E-driven cancers were PTC like, whereas 96.0% RAS-driven cancers had a follicular-like methylation pattern. Interestingly, unlike other diagnoses, follicular variant PTC (FVPTC) samples were split into two methylation clusters (follicular like and PTC like), indicating a heterogeneous group likely to be formed by two distinct diseases. FVPTC samples with a follicular-like methylation pattern were enriched for RAS mutations (36.4% vs. 8.0%; p < 0.001), whereas FVPTC- with PTC-like methylation patterns were enriched for BRAFV600E mutations (52.0% vs. 0%, Fisher exact p = 0.004) and RET fusions (16.0% vs. 0%, Fisher exact p = 0.003). Conclusions: Our data provide novel insights into the epigenetic alterations of thyroid tumors. Since our classification method relies on a fully unsupervised machine learning approach for subtype discovery, our results offer a robust background to support the classification of thyroid neoplasms based on methylation patterns.

Correlation of DNA methylation and lymph node metastasis in papillary thyroid carcinoma

BACKGROUND

Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer with a primarily good prognosis, and its 10-year survival rate is over 90%. However, PTC is prone to early lymph node metastasis.

METHODS

Thyroid cancer tissues from PTC patients with lymphatic metastasis and normal tissues were collected for DNA methylation analysis. Different methylation sites, different methylation regions, gene-enriched pathways, and protein-protein interactions (PPIs) were analyzed.

RESULTS

There were 1004 differentially methylated sites in the PTC group versus the control group; these involved 479 hypermethylated sites in 415 related genes, 525 hypomethylated sites in 482 related genes, 64 differentially methylated regions located in the CpG island region, 34 differentially methylated genes closely related to thyroid cancer, and 17 genes with differentially methylated genes in the DNA promoter region.

CONCLUSION

NDRG4 hypermethylation and FOXO3, ZEB2, and CDK6 hypomethylation were associated with PTC lymph node metastasis.

Association of DNA Promoter Methylation and BRAF Mutation in Thyroid Cancer

The BRAF V600E mutation and DNA promoter methylation play important roles in the pathogenesis of thyroid cancer (TC). However, the association of these genetic and epigenetic alterations is not clear. In this study, using paired tumor and surrounding normal tissue from the same patients, on a genome-wide scale we tried to identify (a) any association between BRAF mutation and DNA promoter methylation, and (b) if the molecular findings may provide a basis for therapeutic intervention. We included 40 patients with TC (female = 28, male = 12) without distant metastasis. BRAF mutation was present in 18 cases. We identified groups of differentially methylated loci (DML) that are found in (a) both BRAF mutant and wild type, (b) only in BRAF mutant tumors, and (c) only in BRAF wild type. BRAF mutation-specific promoter loci were more frequently hypomethylated, whereas BRAF wild-type-specific loci were more frequently hypermethylated. Common DML were enriched in cancer-related pathways, including the mismatch repair pathway and Wnt-signaling pathway. Wild-type-specific DML were enriched in RAS signaling. Methylation status of checkpoint signaling genes, as well as the T-cell inflamed genes, indicated an opportunity for the potential use of PDL1 inhibitors in BRAF mutant TC. Our study shows an association between BRAF mutation and methylation in TC that may have biological significance.

Genomic and epigenomic profile of thyroid cancer

Thyroid cancer is the most common malignancy of the endocrine system, and its incidence has been steadily increasing. Advances in sequencing have allowed analysis of the entire cancer genome, and has provided new information on the genetic lesions and modifications responsible for the onset, progression, dedifferentiation and metastasis of thyroid carcinomas. Moreover, integrated genomics has advanced our understanding of the development of cancer and its behavior, and has facilitated the identification of new genetic mutations and molecular pathways. The functional analysis of epigenetic modifications, such as DNA methylation, histone acetylation and non-coding RNAs, have contributed to define new regulatory mechanisms that control cell malignancy in thyroid cancer, especially aggressive forms. Here we review the most recent advances in genomics and epigenomics of thyroid cancer, which have resulted in a new classification and interpretation of the initiation and progression of thyroid tumors, providing new tools and opportunities for further investigation and for the clinical development of new treatment strategies.

Induction of BRAF^V600E Mutation in Thyroid Cells Leads to Frequent Hypermethylation

Objectives.

The BRAF^V600E mutation is a major driver mutation in papillary thyroid cancer. The aim of this study was to elucidate the correlation between DNA methylation and gene expression changes induced by the BRAF^V600E mutation in thyroid cells.

Methods.

We used Nthy/BRAF cell lines generated by transfection of Nthy/ori cells with the wild-type BRAF gene (Nthy/WT cells) and the V600E mutant-type BRAF gene (Nthy/V600E cells). We performed gene expression microarray and DNA methylation array analyses for Nthy/WT and Nthy/V600E cells. Two types of array data were integrated to identify inverse correlations between methylation and gene expression. The results were verified in silico using data from The Cancer Genome Atlas (TCGA) and in vivo through pyrosequencing and quantitative real-time polymerase chain reaction (qRT-PCR).

Results.

In the Nthy/V600E cells, 199,821 probes were significantly hypermethylated, and 697 genes showed a “hypermethylation-downregulation” pattern in Nthy/V600E. Tumor suppressor genes and apoptosis-related genes were included. In total, 66,446 probes were significantly hypomethylated, and 227 genes showed a “hypomethylation-upregulation” pattern in Nthy/V600E cells. Protooncogenes and developmental protein-coding genes were included. In the TCGA analysis, 491/697 (70.44%) genes showed a hypermethylation-downregulation pattern, and 153/227 (67.40%) genes showed a hypomethylation-upregulation pattern. Ten selected genes showed a similar methylation-gene expression pattern in pyrosequencing and qRT-PCR.

Conclusion.

Induction of the BRAF^V600E mutation in thyroid cells led to frequent hypermethylation. Anticancer genes, such as those involved in tumor suppression or apoptosis, were downregulated by upstream hypermethylation, whereas carcinogenic genes, such as protooncogenes, were upregulated by hypomethylation. Our results suggest that the BRAF^V600E mutation in thyroid cells modulates DNA methylation and results in cancer-related gene expression.

Identification of 4-methylation driven genes based prognostic signature in thyroid cancer: an integrative analysis based on the methylmix algorithm

Thyroid cancer (TC) is known with a high rate of persistence and recurrence. We aimed to develop a prognostic signature to monitor and assess the survival of TC patients. mRNA expression and methylation data were downloaded from the TCGA database. Then, R package methylmix was applied to construct a mixed model was used to identify methylation-driven genes (MDGs) according to the methylation levels. Furthermore, an MDGs based prognostic signature and predictive nomogram were constructed according to the analysis of univariate and multivariate Cox regression. Totally 62 methylation-driven genes that were mainly enriched in substrate-dependent cell migration, cellular response to mechanical stimulus, et al. were found in TC tissues. aldolase C (AldoC), C14orf62, dishevelled 1 (DVL1), and protein tyrosine phosphatase receptor type C (PTPRC) were identified to be significantly related to patients' survival, and may serve as independent prognostic biomarkers for TC. Additionally, the prognostic methylation signature and a novel prognostic, predictive nomogram was established based on the methylation level of 4 MDGs. In this study, we developed a 4-MDGs based prognostic model, which might be the potential predictors for the survival rate of TC patients, and this findings might provide a novel sight for accurate monitoring and prognosis assessment.

Thyroid Cancer Stem-Like Cells: From Microenvironmental Niches to Therapeutic Strategies

Thyroid cancer (TC) is the most common endocrine malignancy. Recent progress in thyroid cancer biology revealed a certain degree of intratumoral heterogeneity, highlighting the coexistence of cellular subpopulations with distinct proliferative capacities and differentiation abilities. Among those subpopulations, cancer stem-like cells (CSCs) are hypothesized to drive TC heterogeneity, contributing to its metastatic potential and therapy resistance. CSCs principally exist in tumor areas with specific microenvironmental conditions, the so-called stem cell niches. In particular, in thyroid cancer, CSCs' survival is enhanced in the hypoxic niche, the immune niche, and some areas with specific extracellular matrix composition. In this review, we summarize the current knowledge about thyroid CSCs, the tumoral niches that allow their survival, and the implications for TC therapy.

Sex-biased DNA methylation in papillary thyroid cancer

Background: Women have a higher risk of developing papillary thyroid cancer (PTC) than men. DNA methylation is known to differ between men and women. Materials & methods: Using the human methylation 450 BeadChip and RNA-sequence, we profiled the genome-wide DNA methylation patterns of papillary thyroid carcinoma patients and para-carcinoma tissue. Results: We first identified 398 different expression genes (DEGs) between males and females PTC. Then we analyzed the relationship between differentially methylated positions (DMPs) and DEGs at gene promoter regions and identified 39 genes and explored DMP-DEGs' correlation with immune cell infiltration and tumor purity. We also analyzed the relationship between genomic regions and enhancers. Conclusion: Our study identified 39 DMP-DEGs providing some new insights into the mechanisms of methylation-mediated gender differences in PTC.

New Insights into Mechanisms of Endocrine-Disrupting Chemicals in Thyroid Diseases: The Epigenetic Way

In recent years, the presence in the environment of chemical compounds with thyroid-disrupting effects is progressively increased. This phenomenon has risen concern for human health as the preservation of thyroid system homeostasis is essential for fetal development and for maintaining psychological and physiological wellbeing. An increasing number of studies explored the role of different classes of toxicants in the occurrence and severity of thyroid diseases, but large epidemiological studies are limited and only a few animal or in vitro studies have attempted to identify the mechanisms of chemical action. Recently, epigenetic changes such as alteration of methylation status or modification of non-coding RNAs have been suggested as correlated to possible deleterious effects leading to different thyroid disorders in susceptible individuals. This review aims to analyze the epigenetic alterations putatively induced by chemical exposures and involved in the onset of frequent thyroid diseases such as thyroid cancer, autoimmune thyroiditis and disruption of fetal thyroid homeostasis.

DNA methylation methods: Global DNA methylation and methylomic analyses

DNA methylation provides a pivotal layer of epigenetic regulation in eukaryotes that has significant involvement for numerous biological processes in health and disease. The function of methylation of cytosine bases in DNA was originally proposed as a "silencing" epigenetic marker and focused on promoter regions of genes for decades. Improved technologies and accumulating studies have been extending our understanding of the roles of DNA methylation to various genomic contexts including gene bodies, repeat sequences and transcriptional start sites. The demand for comprehensively describing DNA methylation patterns spawns a diversity of DNA methylation profiling technologies that target its genomic distribution. These approaches have enabled the measurement of cytosine methylation from specific loci at restricted regions to single-base-pair resolution on a genome-scale level. In this review, we discuss the different DNA methylation analysis technologies primarily based on the initial treatments of DNA samples: bisulfite conversion, endonuclease digestion and affinity enrichment, involving methodology evolution, principles, applications, and their relative merits. This review may offer referable information for the selection of various platforms for genome-wide analysis of DNA methylation.

Overexpression of novel long intergenic non‑coding RNA LINC02454 is associated with a poor prognosis in papillary thyroid cancer

It has been revealed from microarray data analysis that long intergenic non-coding RNA 02454 (LINC02454) is highly expressed in papillary thyroid cancer (PTC). The aim of the present study was to explore the potential role of LINC02454 in the tumorigenesis of PTC. The mRNA expression levels of LINC02454 were assessed using data from The Cancer Genome Atlas (TCGA) and the GSE66783 cohort in thyroid cancer, and were validated using reverse transcription-quantitative PCR in 104 patients with PTC recruited in the present study. The association between the LINC02454 mRNA expression levels and the clinicopathological features of the 104 patients with PTC were also analyzed. Functional enrichment analyses were conducted on the differentially expressed genes in the high and low LINC02454 expression groups that were identified from the TCGA cohort. RNA interference, using short interfering (si)RNA against LINC02454, was used to investigate the role of LINC02454 in the biological functions of PTC cells in vitro. The expression level of LINC02454 was significantly increased in PTC tissues (P=0.0011) and was significantly associated with a larger tumor size, T stage, an advanced TNM stage and an increased lymph node metastasis (P<0.05), which was consistent with that in the TCGA and GSE66783 cohort. High expression levels of LINC02454 were observed in patients with PTC that also had BRAF mutations (P<0.001), and were significantly associated with a poorer disease-free survival in the TCGA cohort (P<0.05). Functional enrichment analysis indicated that LINC02454-related genes were significantly enriched in Gene Ontology terms, such as ‘positive regulation of cell proliferation’, ‘positive regulation of cell division’ and ‘cell adhesion’, and the following Kyoto Encyclopedia of Genes and Genomes pathways: ‘Pathways in cancer’ ‘proteoglycans in cancer’ and ‘ECM-receptor interaction’. In vitro, the knockdown of LINC02454 markedly arrested the cells in the G₀/G₁ phase of the cell cycle, and also led to an overall increase in apoptosis, as well as to an unexpected decrease in cell proliferation. LINC02454 may thus potentially function as an oncogene, which inhibits the apoptosis and enhances proliferation of PTC cells. Thus, as suggested by the findings of the present study, LINC02454 may be used as a diagnostic and prognostic biomarker for PTC in the future.

PD-1, CTLA4, PD-L1 and PD-L2 DNA methylation in papillary thyroid carcinoma

Aim: We investigated DNA methylation patterns of immune checkpoint genes PD-1, PD-L1, PD-L2, CTLA4 and an adjacent long noncoding RNA in papillary thyroid carcinoma (PTC). Materials & methods: DNA methylation and mRNA expression were examined in PTCs. DNA methylation was correlated with mRNA expression, BRAF and RAS mutational status, and immune cell infiltration. Results: Inverse correlations between DNA methylation and mRNA expression were observed. Immune checkpoint expression correlated positively, and DNA methylation negatively, with immune cell infiltration. Higher DNA methylation levels accompanied by lower immune checkpoint expression were observed in RAS-mutated tumors. Conclusion: We suggest epigenetic regulation of immune checkpoints in PTC. Methylation was associated with BRAF and RAS mutation status. DNA methylation might be a promising biomarker candidate in the context of immunotherapies in PTC.

Cancer Stem Cells in Thyroid Tumors: From the Origin to Metastasis

Thyroid tumors are extremely heterogeneous varying from almost benign tumors with good prognosis as papillary or follicular tumors, to the undifferentiated ones with severe prognosis. Recently, several models of thyroid carcinogenesis have been described, mostly hypothesizing a major role of the thyroid cancer stem cell (TCSC) population in both cancer initiation and metastasis formation. However, the cellular origin of TCSC is still incompletely understood. Here, we review the principal epigenetic mechanisms relevant to TCSC origin and maintenance in both well-differentiated and anaplastic thyroid tumors. Specifically, we describe the alterations in DNA methylation, histone modifiers, and microRNAs (miRNAs) involved in TCSC survival, focusing on the potential of targeting aberrant epigenetic modifications for developing novel therapeutic approaches. Moreover, we discuss the bidirectional relationship between TCSCs and immune cells. The cells of innate and adaptive response can promote the TCSC-driven tumorigenesis, and conversely, TCSCs may favor the expansion of immune cells with protumorigenic functions. Finally, we evaluate the role of the tumor microenvironment and the complex cross-talk of chemokines, hormones, and cytokines in regulating thyroid tumor initiation, progression, and therapy refractoriness. The re-education of the stromal cells can be an effective strategy to fight thyroid cancer. Dissecting the genetic and epigenetic landscape of TCSCs and their interactions with tumor microenvironment cells is urgently needed to select more appropriate treatment and improve the outcome of patients affected by advanced differentiated and undifferentiated thyroid cancers.

Multi-omics Signatures and Translational Potential to Improve Thyroid Cancer Patient Outcome

Recent advances in high-throughput molecular and multi-omics technologies have improved our understanding of the molecular changes associated with thyroid cancer initiation and progression. The translation into clinical use based on molecular profiling of thyroid tumors has allowed a significant improvement in patient risk stratification and in the identification of targeted therapies, and thereby better personalized disease management and outcome. This review compiles the following: (1) the major molecular alterations of the genome, epigenome, transcriptome, proteome, and metabolome found in all subtypes of thyroid cancer, thus demonstrating the complexity of these tumors and (2) the great translational potential of multi-omics studies to improve patient outcome.

Methylation of tumour suppressor genes associated with thyroid cancer

BACKGROUND

Thyroid carcinoma is the most common endocrine malignancy worldwide. Changes in DNA methylation can cause silencing of normally active genes, especially tumour suppressor genes (TSG) or activation of normally silent genes.

OBJECTIVE

The aim of this study is to evaluate the degree of promoter methylation for a panel of markers for thyroid neoplasms and to establish their relationship with thyroid oncogenesis.

METHODS

To generate a comprehensive DNA methylation signature of TSGs involved in thyroid neoplasia, we use Human TSG EpiTect Methyl II Signature PCR Array-Qiagen for 24 samples (follicular adenomas and papillary thyroid carcinomas) compared with normal thyroid tissue. We extended the evaluation for three TSGs (TP73, WIF1, PDLIM4) using qMS-PCR. Statistical analysis was performed with GraphPad Prism.

RESULTS

We noted four important genes NEUROG1, ESR1, RUNX3, MLH1, which presented methylated promoter in tumour samples compared to normal. We found new characteristic of thyroid tumours: methylation of TP73, WIF1 and PDLIM4 TSGs, which can contribute to thyroid neoplasia. A significant correlation between BRAF V600E mutation and RET/PTC rearrangements with TIMP3 and CDH13, RARB methylation, respectively was observed.

CONCLUSIONS

TSGs promoter hypermethylation is a hallmark of cancer and a test that uses methylation quantification method is suitable for diagnosis and prognosis of thyroid cancer.

Investigating the expression and promoter methylation of RET gene in patients with medullary thyroid cancer with unmutated RET

OBJECTIVE

Thyroid cancer is one of the most common endocrine malignancies. Mutations in the rearranged during transfection (RET) gene, especially in exon 10, 11, and 16, as well as epigenetic modifications, constitute the major underlying molecular events leading to medullary thyroid cancer (MTC). There are few studies on the mutations and epigenetic changes of RET gene in Iranian patients with MTC. In the present study, we aimed to address this question and explore the clinical relevance of such genetic alternations in an Iranian population.

METHODS

Thirty-three patients with confirmed MTC who underwent thyroidectomy surgery in Imam Khomeini Hospital (Tehran, Iran) were enrolled. DNA extracted from cancerous tissues was amplified by polymerase chain reaction (PCR) and then was sequenced for identification of RET mutations. In patients with no identified mutations, the methylation status of RET promoter and its expression were further investigated using methylation-specific PCR and real-time PCR methods, respectively.

RESULTS

In MTC patients with no RET mutations, the promoter of the proto-oncogene was hypomethylated. Furthermore, RET gene expression was elevated in patients who revealed no mutations in neither of exon 10, 11, or 16 of the RET gene.

CONCLUSION

Hypomethylation of RET promoter may contribute to MTC pathogenesis. The methylation status of RET promoter could be a new potential prognostic, diagnostic and therapeutic marker in MTC.

Identification of DNA motifs that regulate DNA methylation

DNA methylation is an important epigenetic mark but how its locus-specificity is decided in relation to DNA sequence is not fully understood. Here, we have analyzed 34 diverse whole-genome bisulfite sequencing datasets in human and identified 313 motifs, including 92 and 221 associated with methylation (methylation motifs, MMs) and unmethylation (unmethylation motifs, UMs), respectively. The functionality of these motifs is supported by multiple lines of evidence. First, the methylation levels at the MM and UM motifs are respectively higher and lower than the genomic background. Second, these motifs are enriched at the binding sites of methylation modifying enzymes including DNMT3A and TET1, indicating their possible roles of recruiting these enzymes. Third, these motifs significantly overlap with "somatic QTLs" (quantitative trait loci) of methylation and expression. Fourth, disruption of these motifs by mutation is associated with significantly altered methylation level of the CpGs in the neighbor regions. Furthermore, these motifs together with somatic mutations are predictive of cancer subtypes and patient survival. We revealed some of these motifs were also associated with histone modifications, suggesting a possible interplay between the two types of epigenetic modifications. We also found some motifs form feed forward loops to contribute to DNA methylation dynamics.

New molecular biomarkers in differentiated thyroid carcinoma: Impact of miR-146, miR-221 and miR-222 levels in the evolution of the disease

OBJECTIVE

MircroRNAs (miR) are small, noncoding RNA molecules of 18-25 nucleotides. Their dysregulation has been widely studied in many human tumours including differentiated thyroid cancer (DTC). miRs more frequently associated with these kinds of tumours are miR-146, miR-221 and miR-222. Our objective was to assess the relationship among circulating miR levels and the evolution and outcomes of disease.

DESIGN

We analysed a sample of 60 patients with DTC assigning them to one of three groups according to the dynamic scale of risk (excellent response, incomplete biochemical response and incomplete structural response).

PATIENTS AND MEASUREMENTS

At study inclusion, we determined thyroid-stimulating hormone, thyroxine, thyroglobulin, antithyroglobulin antibodies and plasma levels of miR-146, miR-221 and miR-222.

RESULTS

Male sex and advanced age at diagnosis were associated with the worst disease progression. miR-222 was twofold to threefold higher in tall cell papillary carcinomas (P = 0.038). miR-146 (P = 0.016) and miR-221 (P = 0.050) had a positive correlation with thyroglobulin at the time of sampling. In regression analysis, miR-146 (P = 0.006), miR-221 (P = 0.004) and miR-222 (P = 0.007) predicted more than 70% of the variation in thyroglobulin levels at the time of sampling.

CONCLUSIONS

Elevated miR-222 and miR-146 levels are associated with poorer outcomes of the disease and may have a prognostic value in the management and follow-up of DTC.

DNA methylation in thyroid cancer

In recent years, cancer genomics has provided new insights into genetic alterations and signaling pathways involved in thyroid cancer. However, the picture of the molecular landscape is not yet complete. DNA methylation, the most widely studied epigenetic mechanism, is altered in thyroid cancer. Recent technological advances have allowed the identification of novel differentially methylated regions, methylation signatures and potential biomarkers. However, despite recent progress in cataloging methylation alterations in thyroid cancer, many questions remain unanswered. The aim of this review is to comprehensively examine the current knowledge on DNA methylation in thyroid cancer and discuss its potential clinical applications. After providing a general overview of DNA methylation and its dysregulation in cancer, we carefully describe the aberrant methylation changes in thyroid cancer and relate them to methylation patterns, global hypomethylation and gene-specific alterations. We hope this review helps to accelerate the use of the diagnostic, prognostic and therapeutic potential of DNA methylation for the benefit of thyroid cancer patients.

Association between DNA methylation profile and malignancy in follicular-patterned thyroid neoplasms

Molecular differentiation between benign (follicular thyroid adenoma, FTA) and malignant (follicular thyroid carcinoma, FTC) thyroid neoplasms is challenging. Here, we explored the genome-wide DNA methylation profile of FTA (n.10) and FTC (n.11) compared to normal thyroid (NT) (n.7) tissues. FTC featured 3,564 differentially-methylated CpGs (DMCpG), most (84%) of them hypermethylated, with respect to normal controls. At the principal component analysis (PCA), the methylation profile of FTA occupied an intermediate position between FTC and normal tissue. A large fraction (n. 2,385) of FTC-associated DMCpG were related (intragenic or within 1500 bp from the transcription start site) to annotated genes (n. 1,786). FTC-hypermethylated genes were enriched for targets of the Polycomb transcriptional repressor complex and the specific histone H3 marks (H3K4me2/me3-H3K27me3) found in chromatin domains known as "bivalent". Transcriptome profiling by RNAseq showed that 7.9% of the DMCpGs-associated genes were differentially expressed in FTC compared to NT, suggesting that altered DNA methylation may contribute to their altered expression. Overall, this study suggests that perturbed DNA methylation, in particular hypermethylation, is a component of the molecular mechanisms leading to the formation of FTC and that DNA methylation profiling may help differentiating FTCs from their benign counterpart.

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.

A review on the mechanism of iodide metabolic dysfunction in differentiated thyroid cancer

The incidence of differentiated thyroid cancer (DTC) has been increasing rapidly worldwide, and the risk factors remain unclear. With the growing number of patients with DTC, the related issues have been gradually highlighted. ¹³¹Iodide (¹³¹I) is an important treatment for DTC and has the potential to reduce the risk of recurrence. ¹³¹I is also an effective treatment for distant metastases of thyroid carcinoma. However, iodide metabolism dysfunction in metastatic foci causes patients to lose the opportunity of ¹³¹I treatment. This article reviews the related mechanisms of iodide metabolism dysfunction in DTC cells and summarizes the clinical transformation progression.

Altered Epigenetic Mechanisms in Thyroid Cancer Subtypes

Thyroid carcinoma (TC) is the most frequent malignant neoplasm of the endocrine system. Molecular methods for diagnosis of invasive thyroid disease can be effectively adopted. Epigenetic factors play an important role in the diversity patterns of gene expression and the phenotypic and biological characteristics of TC subtypes. We aimed to review epigenetic changes in the main subtypes of TC, along with a presentation of the methods that have examined these changes, and active clinical trials for the treatment of advanced TCs targeting epigenetic changes. A literature analysis was performed in MEDLINE using PubMed, Elsevier, and Google Scholar for studies published up to 2016, using the keywords: "Epigenetic alterations" OR "Epigenetic changes", "thyroid cancers", "papillary thyroid cancer", "medullary thyroid cancer", "follicular thyroid cancer", and "anaplastic thyroid cancer", which resulted in 310 articles in English. All related abstracts were reviewed and studies were included that were published in English, had available full text, and determined the details of the methods and materials associated with the epigenetic patterns of TC and its subtypes (100 articles). Analysis of epigenetic alterations in TC subtypes helps to identify pathogenesis and can play an important role in the classification and diagnosis of tumors. Epigenetic mechanisms, especially aberrant methylation of DNA and microRNAs (miRs), are likely to play an important role in thyroid tumorigenesis. Further studies are required to elucidate the role of histone modification mechanisms in TC development.

Identification of Tissue-Specific DNA Methylation Signatures for Thyroid Nodule Diagnostics

PURPOSE

Thyroid cancer is frequently difficult to diagnose due to an overlap of cytologic features between malignant and benign nodules. This overlap leads to unnecessary removal of the thyroid in patients without cancer. While providing some improvement over cytopathologic diagnostics, molecular methods frequently fail to provide a correct diagnosis for thyroid nodules. These approaches are based on the difference between cancer and adjacent thyroid tissue and assume that adjacent tissues are the same as benign nodules. However, in contrast to adjacent tissues, benign thyroid nodules can contain genetic alterations that can be found in cancer.Experimental Design: For the development of a new molecular diagnostic test for thyroid cancer, we evaluated DNA methylation in 109 thyroid tissues by using genome-wide single-base resolution DNA methylation analysis. The test was validated in a retrospective cohort containing 65 thyroid nodules.

RESULTS

By conducting reduced representation bisulfite sequencing in 109 thyroid specimens, we found significant differences between adjacent tissue, benign nodules, and cancer. These tissue-specific signatures are strongly linked to active enhancers and cancer-associated genes. Based on these signatures, we developed a new epigenetic approach for thyroid diagnostics. According to the validation cohort, our test has an estimated specificity of 97% [95% confidence interval (CI), 81-100], sensitivity of 100% (95% CI, 87-100), positive predictive value of 97% (95% CI, 83-100), and negative predictive value of 100% (95% CI, 86-100).

CONCLUSIONS

These data show that epigenetic testing can provide outstanding diagnostic accuracy for thyroid nodules.See related commentary by Mitmaker et al., p. 457.

Integrated bioinformatics analysis reveals that the expression of cathepsin S is associated with lymph node metastasis and poor prognosis in papillary thyroid cancer

The prognosis of the majority of patients with papillary thyroid cancer (PTC) is excellent, although there are patients who experience disease recurrence and progression. The aim of the present study was to identify potential prognostic risk markers in PTC. Differentially expressed genes (DEGs), identified from four Genome Expression Omnibus cohorts were subjected to functional enrichment analyses with Gene Ontology terms and the Kyoto Encyclopedia of Genes and Genome pathways. Hub genes, filtered from cytoHubba, were validated using the The Cancer Genome Atlas (TCGA) cohort, and their associations with clinicopathological features and prognosis were analyzed. A total of 277 DEGs were identified following data preprocessing. DEGs were primarily enriched in 'small cell lung cancer', 'ECM-receptor interaction', 'pathways in cancer'and 'tyrosine metabolism'. Hub genes [APOE, cathepsin S (CTSS), insulin receptor substrate 1 (IRS1), KIT, LGALS3, RUNX2 and TGFBR1] were extracted from cytoHubba. Their expression in the TCGA cohort was consistent with that in the GEO cohorts. CTSS (P=0.006) and IRS1 (P=0.005) were associated with disease‑free survival, as determined using the Kaplan-Meier analysis. CTSS was an independent risk factor for poor disease‑free survival (HR, 2.649; 95% CI, 1.095-6.409; P=0.031). Patients with high expression of CTSS exhibited different histological types (increased tall-cell subtype and reduced follicular subtype; P<0.001), more frequent lymph node metastasis (P<0.001) and advanced tumor-node-metastasis stages (P=0.049) compared with the low-expression group. High expression of CTSS was independently associated with lymph node metastasis (OR, 2.015; 95% CI, 1.225-3.315; P=0.006). Therefore, CTSS may serve as a predictive risk marker for the progression and prognosis of PTC.

NKX2.5 is expressed in papillary thyroid carcinomas and regulates differentiation in thyroid cells

BACKGROUND

NKX2.5 is a transcription factor transiently expressed during thyroid organogenesis. Recently, several works have pointed out the oncogenic role of NKX2.5 in a variety of tumors. We therefore hypothesized that NKX2.5 could also play a role in thyroid cancer.

METHODS

The validation of NKX2.5 expression was assessed by immunohistochemistry analysis in a Brazilian case series of 10 papillary thyroid carcinoma (PTC) patients. Then, the long-term prognostic value of NKX2.5 and its correlation with clinicopathologic features of 51 PTC patients was evaluated in a cohort with 10-years follow-up (1990-1999). Besides, the effect of NKX2.5 overexpression on thyroid differentiation markers and function was also investigated in a non-tumor thyroid cell line (PCCL3).

RESULTS

NKX2.5 was shown to be expressed in most PTC samples (8/10, case series; 27/51, cohort). Patients who had tumors expressing NKX2.5 showed lower rates of persistence/recurrence (p = 0.013). Overexpression of NKX2.5 in PCCL3 cells led to: 1) downregulation of thyroid differentiation markers (thyrotropin receptor, thyroperoxidase and sodium-iodide symporter); 2) reduced iodide uptake; 3) increased extracellular H₂O₂ generation, dual oxidase 1 mRNA levels and activity of DuOx1 promoter.

CONCLUSIONS

In summary, NKX2.5 is expressed in most PTC samples analyzed and its presence correlates to better prognosis of PTC. In vitro, NKX2.5 overexpression reduces the expression of thyroid differentiation markers and increases ROS production. Thus, our data suggests that NKX2.5 could play a role in thyroid carcinogenesis.

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.

Increased Global DNA Hypomethylation in Distant Metastatic and Dedifferentiated Thyroid Cancer

CONTEXT

Global DNA hypomethylation is a major event for the development and progression of cancer, although the significance in thyroid cancer remains unclear. Therefore, we aimed to investigate its role in thyroid cancer progression and its potential as a prognostic marker.

METHODS

Global hypomethylation of Alu repeats was used as a surrogate marker for DNA global hypomethylation, and was assessed using the Quantification of Unmethylated Alu technique. Mutations in BRAF and RAS were determined by Sanger sequencing.

RESULTS

Ninety primary thyroid tumors were included [28 low-risk differentiated thyroid cancer (DTC), 13 pediatric DTC, 33 distant metastatic DTC, 7 poorly differentiated thyroid cancer (PDTC), and 9 anaplastic thyroid cancer (ATC)], as well as 24 distant metastases and 20 normal thyroid tissues. An increasing hypomethylation was found for distant metastatic DTC [median, 4.0; interquartile range (IQR), 3.1 to 6.2] and PDTC/ATC (median, 9.3; IQR, 7.0 to 12.1) as compared with normal thyroid tissue (median, 2.75; IQR, 2.30 to 3.15), whereas low-risk and pediatric DTC were not affected by hypomethylation. Alu hypomethylation was similar between distant metastases and matched primary tumors. Within distant metastatic DTC, Alu hypomethylation was increased in BRAF vs RAS mutated tumors. Kaplan-Meier and Cox regression analyses showed that thyroid cancer-related and all-cause mortality were associated with tumor hypomethylation, but this association was lost after adjustment for thyroid cancer risk category.

CONCLUSION

Distant metastatic DTC, PDTC, and ATC were increasingly affected by global Alu hypomethylation, suggesting that this epigenetic entity may be involved in thyroid cancer progression and dedifferentiation.

The relationship between RASSF1A promoter methylation and thyroid carcinoma: A meta-analysis of 14 articles and a bioinformatics of 2 databases (PRISMA)

BACKGROUND

DNA promoter methylation can suppresses gene expression and shows an important role in the biological functions of Ras association domain family 1A (RASSF1A). Many studies have performed to elucidate the role of RASSF1A promoter methylation in thyroid carcinoma, while the results were conflicting and heterogeneous. Here, we analyzed the data of databases to determine the relationship between RASSF1A promoter methylation and thyroid carcinoma.

METHODS

We used the data from 14 cancer-normal studies and Gene Expression Omnibus (GEO) database to analyze RASSF1A promoter methylation in thyroid carcinoma susceptibility. The data from the Cancer Genome Atlas project (TCGA) database was used to analyze the relationship between RASSF1A promoter methylation and thyroid carcinoma susceptibility, clinical characteristics, prognosis. Odds ratios were estimated for thyroid carcinoma susceptibility and hazard ratios were estimated for thyroid carcinoma prognosis. The heterogeneity between studies of meta-analysis was explored using H, I values, and meta-regression. We adopted quality criteria to classify the studies of meta-analysis. Subgroup analyses were done for thyroid carcinoma susceptibility according to ethnicity, methods, and primers.

RESULTS

Result of meta-analysis indicated that RASSF1A promoter methylation is associated with higher susceptibility to thyroid carcinoma with small heterogeneity. Similarly, the result from GEO database also showed that a significant association between RASSF1A gene promoter methylation and thyroid carcinoma susceptibility. For the results of TCGA database, we found that RASSF1A promoter methylation is associated with susceptibility and poor disease-free survival (DFS) of thyroid carcinoma. In addition, we also found a close association between RASSF1A promoter methylation and patient tumor stage and age, but not in patients of different genders.

CONCLUSIONS

The methylation status of RASSF1A promoter is strongly associated with thyroid carcinoma susceptibility and DFS. The RASSF1A promoter methylation test can be applied in the clinical diagnosis of thyroid carcinoma.

Prognostic Classifier Based on Genome-Wide DNA Methylation Profiling in Well-Differentiated Thyroid Tumors

Context

Even though the majority of well-differentiated thyroid carcinoma (WDTC) is indolent, a number of cases display an aggressive behavior. Cumulative evidence suggests that the deregulation of DNA methylation has the potential to point out molecular markers associated with worse prognosis.

Objective

To identify a prognostic epigenetic signature in thyroid cancer.

Design

Genome-wide DNA methylation assays (450k platform, Illumina) were performed in a cohort of 50 nonneoplastic thyroid tissues (NTs), 17 benign thyroid lesions (BTLs), and 74 thyroid carcinomas (60 papillary, 8 follicular, 2 Hürthle cell, 1 poorly differentiated, and 3 anaplastic). A prognostic classifier for WDTC was developed via diagonal linear discriminant analysis. The results were compared with The Cancer Genome Atlas (TCGA) database.

Results

A specific epigenetic profile was detected according to each histological subtype. BTLs and follicular carcinomas showed a greater number of methylated CpG in comparison with NTs, whereas hypomethylation was predominant in papillary and undifferentiated carcinomas. A prognostic classifier based on 21 DNA methylation probes was able to predict poor outcome in patients with WDTC (sensitivity 63%, specificity 92% for internal data; sensitivity 64%, specificity 88% for TCGA data). High-risk score based on the classifier was considered an independent factor of poor outcome (Cox regression, P < 0.001).

Conclusions

The methylation profile of thyroid lesions exhibited a specific signature according to the histological subtype. A meaningful algorithm composed of 21 probes was capable of predicting the recurrence in WDTC.

GABRB2 plays an important role in the lymph node metastasis of papillary thyroid cancer

BACKGROUND

Thyroid cancer is a common malignant tumor of the endocrine system. Its incidence has increased continuously worldwide for the past three decades. With advanced sequencing technology, we discovered that GABRB2 gene is overexpressed in tumor tissues and closely associated with vertebrate nervous systems. However, its role in cancer remains unclear.

METHODS

We conducted a massively parallel whole transcriptome resequencing and a comprehensive analysis of matched papillary thyroid carcinoma (PTC) tumors and normal tissues in 19 patients. Results showed that GABRB2 expression was significantly upregulated in thyroid cancer. Forty-five pairs of tumors and normal tissues were subjected to reverse transcription polymerase chain reaction to validate previous findings. The specific functions of GABRB2 in PTC cell lines (BCPAP, TPC1, and KTC-1) transfected with small interfering RNA were determined through cell colony formation, Cell Counting Kit-8, Transwell migration, Transwell invasion, and apoptosis assays. The effect of DNA demethylation on this gene was also examined.

RESULTS

GABRB2 was remarkably overexpressed in primarily sequenced PTC tumors and validation cohort (T: N = 4.94 ± 3.43:0.83 ± 1.71, P < 0.001), and this observation was consistent with that in the TCGA cohort (T: N = 38.92 ± 35.53:0.30 ± 0.55, P < 0.001). GABRB2 overexpression was correlated with lymph node metastasis in both cohorts (P < 0.01). In vitro experiments revealed that GABRB2 downregulation significantly inhibited the colony formation, migration, and invasion of the three PTC cell lines.

CONCLUSION

GABRB2 plays important tumorigenic functions and acts as a novel oncogene in PTC.

Comparative analysis of MBD-seq and MeDIP-seq and estimation of gene expression changes in a rodent model of schizophrenia

We conducted a comparative study of multiplexed affinity enrichment sequence methodologies (MBD-seq and MeDIP-seq) in a rodent model of schizophrenia, induced by in utero methylazoxymethanol acetate (MAM) exposure. We also examined related gene expression changes using a pooled sample approach. MBD-seq and MeDIP-seq identified 769 and 1771 differentially methylated regions (DMRs) between F2 offspring of MAM-exposed rats and saline control rats, respectively. The assays showed good concordance, with ~56% of MBD-seq-detected DMRs being identified by or proximal to MeDIP-seq DMRs. There was no significant overlap between DMRs and differentially expressed genes, suggesting that DNA methylation regulatory effects may act upon more distal genes, or are too subtle to detect using our approach. Methylation and gene expression gene ontology enrichment analyses identified biological processes important to schizophrenia pathophysiology, including neuron differentiation, prepulse inhibition, amphetamine response, and glutamatergic synaptic transmission regulation, reinforcing the utility of the MAM rodent model for schizophrenia research.

Dual Inhibition of HDAC and Tyrosine Kinase Signaling Pathways with CUDC-907 Inhibits Thyroid Cancer Growth and Metastases

Purpose: There is currently no standard therapy for anaplastic thyroid cancer (ATC) and poorly differentiated thyroid cancer (PDTC), which account for two-thirds of thyroid cancer-related deaths. Driver mutations in the PI3K/AKT and RAF/RAS/MEK/ERK pathways are common in ATC and PDTC. Histone deacetylases (HDAC) regulate cancer initiation and progression. Our aim was to determine the therapeutic efficacy of simultaneously targeting these pathways in thyroid cancer with a single agent and to evaluate biomarkers of treatment response.Experimental Design: CUDC-907 is a first-in-class compound, functioning as a dual inhibitor of HDACs and the PI3K/AKT pathway. We investigated its antiproliferative effect in vitro and in vivoResults: CUDC-907 significantly inhibited cellular proliferation in thyroid cancer cell lines, induced G₂-M arrest with decreased levels of the checkpoint regulators cyclin B1, AURKA, AURKB, PLK1, and increased p21 and p27. Treatment induced apoptosis with increased caspase-3/7 activity and decreased survivin levels and decreased cellular migration and invasion. CUDC-907 treatment caused H3 hyperacetylation and decreased HDAC2 expression. HDAC2 was upregulated in ATC and other thyroid cancer histologic subtypes. CUDC-907 treatment reduced both p-AKT and p-ERK1/2 levels. Finally, CUDC-907 treatment, in a metastatic mouse model of thyroid cancer, showed significant inhibition of growth and metastases, and tumors from treated mice had decreased HDAC2 expression, suggesting that this may be a useful biomarker of response.Conclusions: Dual inhibition of HDAC and the tyrosine kinase signaling pathways with CUDC-907 is a promising treatment strategy for advanced, metastatic thyroid cancer. Clin Cancer Res; 23(17); 5044-54. ©2017 AACR.

Intrinsic LINE-1 Hypomethylation and Decreased Brca1 Expression are Associated with DNA Repair Delay in Irradiated Thyroid Cells

Exposure to ionizing radiation greatly increases the risk of developing papillary thyroid carcinoma (PTC), especially during childhood, mainly due to gradual inactivation of DNA repair genes and DNA damages. Recent molecular characterization of PTC revealed DNA methylation deregulation of several promoters of DNA repair genes. Thus, epigenetic silencing might be a plausible mechanism for the activity loss of tumor suppressor genes in radiation-induced thyroid tumors. Herein, we investigated the impact of ionizing radiation on global methylation and CpG islands within promoter regions of homologous recombination (HR) and non-homologous end joining (NHEJ) genes, as well as its effects on gene expression, using two well-established normal differentiated thyroid cell lines (FRTL5 and PCCL3). Our data reveal that X-ray exposure promoted G₂/M arrest in normal thyroid cell lines. The FRTL5 cells displayed a slower kinetics of double-strand breaks (DSB) repair and a lower long interspersed nuclear element-1 (LINE-1) methylation than the PCCL3 cells. Nevertheless, acute X-ray exposure does not alter the expression of genes involved in HR and NHEJ pathways, apart from the downregulation of Brca1 in thyroid cells. On the other hand, HR and NHEJ gene expressions were upregulated in radiation-induced senescent thyroid cells. Taken together, these data suggest that FRTL5 cells intrinsically have less efficient DNA DSB repair machinery than PCCL3 cells, as well as genomic instability, which could predispose the FRTL5 cells to unrepaired DSB lesions and, therefore, gene mutations.

Integrated data analysis reveals potential drivers and pathways disrupted by DNA methylation in papillary thyroid carcinomas

BACKGROUND

Papillary thyroid carcinoma (PTC) is a common endocrine neoplasm with a recent increase in incidence in many countries. Although PTC has been explored by gene expression and DNA methylation studies, the regulatory mechanisms of the methylation on the gene expression was poorly clarified. In this study, DNA methylation profile (Illumina HumanMethylation 450K) of 41 PTC paired with non-neoplastic adjacent tissues (NT) was carried out to identify and contribute to the elucidation of the role of novel genic and intergenic regions beyond those described in the promoter and CpG islands (CGI). An integrative and cross-validation analysis were performed aiming to identify molecular drivers and pathways that are PTC-related.

RESULTS

The comparisons between PTC and NT revealed 4995 methylated probes (88% hypomethylated in PTC) and 1446 differentially expressed transcripts cross-validated by the The Cancer Genome Atlas data. The majority of these probes was found in non-promoters regions, distant from CGI and enriched by enhancers. The integrative analysis between gene expression and DNA methylation revealed 185 and 38 genes (mainly in the promoter and body regions, respectively) with negative and positive correlation, respectively. Genes showing negative correlation underlined FGF and retinoic acid signaling as critical canonical pathways disrupted by DNA methylation in PTC. BRAF mutation was detected in 68% (28 of 41) of the tumors, which presented a higher level of demethylation (95% hypomethylated probes) compared with BRAF wild-type tumors. A similar integrative analysis uncovered 40 of 254 differentially expressed genes, which are potentially regulated by DNA methylation in BRAFV600E-positive tumors. The methylation and expression pattern of six selected genes (ERBB3, FGF1, FGFR2, GABRB2, HMGA2, and RDH5) were confirmed as altered by pyrosequencing and RT-qPCR.

CONCLUSIONS

DNA methylation loss in non-promoter, poor CGI and enhancer-enriched regions was a significant event in PTC, especially in tumors harboring BRAFV600E. In addition to the promoter region, gene body and 3'UTR methylation have also the potential to influence the gene expression levels (both, repressing and inducing). The integrative analysis revealed genes potentially regulated by DNA methylation pointing out potential drivers and biomarkers related to PTC development.

CXCL12 methylation-mediated epigenetic regulation of gene expression in papillary thyroid carcinoma

Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer, and its incidence rate is rapidly growing. It is necessary to understand the pathogenesis of PTC to develop effective diagnosis methods. Promoter methylation has been recognized to contribute to the alterations in gene expression observed in tumorigenesis. Our RNA-seq data identified 1191 differentially expressed mRNAs and 147 differentially expressed lncRNAs in PTC. Next, promoter methylation of these genes was detected by reduced representation bisulfite sequencing (RRBS) technology and comprehensively analyzed to identify differential methylation. In total, 14 genes (13 mRNAs and 1 lncRNA), in which methylation was intimately involved in regulating gene expression, were proposed as novel diagnostic biomarkers. To gain insights into the relationships among these 14 genes, a core co-function network was constructed based on co-expression, co-function and co-methylation data. Notably, CXCL12 was identified as an essential gene in the network that was closely connected with the other genes. These data suggested that CXCL12 down-regulation in PTC may be caused by promoter hypermethylation. Our study was the first to perform an RRBS analysis for PTC and suggested that CXCL12 may contribute to PTC development by methylation-mediated epigenetic regulation of gene expression.

Nodular thyroid disease in the elderly: novel molecular approaches for the diagnosis of malignancy

Epithelial thyroid cancers (TC) comprise two differentiated histotypes (DTC), the papillary (PTC) and the follicular (FTC) thyroid carcinomas which, following dedifferentiation, are assumed to give rise to the poorly differentiated thyroid carcinomas and the rare, but highly aggressive and invariably fatal, anaplastic thyroid carcinomas. Although thyroid cancer mortality has not been changed, its annual incidence has increased over the last two decades, mainly because of the improved ability to diagnose malignant transformation in small non-palpable thyroid nodules. Despite DTC patients have a favorable prognosis, aggressive disease is more frequently observed in the elderly showing a higher disease-specific mortality. Of relevance is the high prevalence of nodular thyroid disease in aged patients being higher than 90%, in women older than 60 year, and 60% in men older than 80 year. This implies a careful evaluation of thyroid nodules in this group of patients in order to exclude malignancy. In fact, despite the tremendous progress in the comprehension of the underlying molecular mechanisms deregulated in DTC progression, several aspects of their clinical management remain to be solved and novel diagnostic strategies are sorely needed. Here, we will attempt to review new molecular approaches, which are currently being exploited in order to ameliorate the diagnosis of thyroid nodules.

Epigenetic modulators of thyroid cancer

There are some well known factors involved in the etiology of thyroid cancer, including iodine deficiency, radiation exposure at early ages, or some genetic changes. However, epigenetic modulators that may contribute to development of these tumors and be helpful to for both their diagnosis and treatment have recently been discovered. The currently known changes in DNA methylation, histone modifications, and non-coding RNAs in each type of thyroid carcinoma are reviewed here.

RASSF1A promoter methylation is associated with increased risk of thyroid cancer: a meta-analysis

Objective

Previous studies have reported that Ras-associated domain family 1A (RASSF1A), the most commonly silenced tumor suppressor via promoter methylation, played vital roles in the development of carcinogenesis. The purpose of this meta-analysis was to determine whether RASSF1A promoter methylation increased the risk of thyroid cancer.

Methods

PubMed, Embase, ISI Web of Knowledge, and Chinese National Knowledge Infrastructure databases were searched to obtain eligible studies. The pooled odds ratios (ORs) and 95% confidence intervals (CIs) were calculated to estimate the strength of the associations, using Stata 12.0 software. The methodological quality of included studies was evaluated using Newcastle–Ottawa scale table. Egger’s test and Begg’s test were applied to detect publication biases. TSA 0.9 software was used to calculate the required information size and whether the result was conclusive.

Results

A total of 10 articles with 12 studies that included 422 thyroid cancer patients, identifying the association of RASSF1A promoter methylation with thyroid cancer risk, were collected in this meta-analysis. Overall, RASSF1A promoter methylation significantly increased the risk of thyroid cancer (total, OR=8.27, CI=4.38–15.62, P<0.05; Caucasian, OR=9.25, CI=3.97–21.56, P<0.05; Asian, OR=7.01, CI=2.68–18.38, P<0.05). In the subgroup analysis based on sample type, a significant association between thyroid cancer group and control group was found (normal tissue, OR=9.55, CI=4.21–21.67, P<0.05; adjacent tissue, OR=6.80, CI=2.49–18.56, P<0.05). The frequency of RASSF1A promoter methylation in follicular thyroid carcinoma was higher than in control group (OR=11.88, CI=5.80–24.32, P<0.05). In addition, the results indicated that the RASSF1A promoter methylation was correlated with papillary thyroid carcinoma in Caucasians and Asians (total, OR=8.07, CI=3.54–18.41, P<0.05; Caucasian, OR=11.35, CI=2.39–53.98, P<0.05; Asian, OR=6.67, CI=2.53–17.64, P<0.05). On the basis of the trial sequential analysis, the significant association of RASSF1A promoter methylation with thyroid cancer risk was found, and there was no need to perform further studies.

Conclusion

This meta-analysis confirms that RASSF1A promoter methylation is a risk factor for thyroid tumor.

Epigenetic Alterations and Canonical Pathway Disruption in Papillary Thyroid Cancer: A Genome-wide Methylation Analysis

BACKGROUND

Alterations in DNA methylation have been demonstrated in a variety of malignancies, including papillary thyroid cancer (PTC). The full extent of dysregulation in PTC and the downstream affected pathways remains unclear. Here we report a genome-wide analysis of PTC methylation, the dysregulation of various canonical pathways, and assess its potential as a diagnostic test.

METHODS

A discovery set utilized 49 PTCs and matched normal controls from The Cancer Genome Atlas. Another set of 16 PTCs and 13 normal controls were used as a replication set. Genome-wide methylation analysis was done using Illumina 450 K methylation chips. Differentially methylated loci (DML) were identified by comparing PTC and matched normal tissues. DML were defined as false-discovery rate p < 0.05 and absolute Δβ ≥ 0.2. DML were then analyzed for pathway and disease commonalities using Qiagen Ingenuity Pathway Analysis.

RESULTS

Of 485,577 CpG sites analyzed, 1226 DML were identified in our discovery and replication sets, and 1061 (86.5 %) DML showed hypomethylation when comparing tumor with normal tissue. Support vector machine classification was able to differentiate benign from malignant tissue in 107 (94.7 %) of 113 tested samples, including 15 (83.3 %) of 18 samples lacking a clearly deleterious mutation. Statistically significant associations with multiple canonical pathways, diseases, and biofunctions were observed including PI3K, PTEN, wnt/β-catenin, and p53.

CONCLUSIONS

Epigenetic dysregulation of multiple canonical pathways are associated with the development of PTC. This methylation signature shows promise as a future adjunctive screening test for thyroid nodules.

Distinct methylation profiles characterize fusion-positive and fusion-negative rhabdomyosarcoma

Rhabdomyosarcoma comprises two major subtypes, fusion positive (PAX3-FOXO1 or PAX7-FOXO1) and fusion negative. To investigate the significance of DNA methylation in these subtypes, we analyzed methylation profiles of 37 rhabdomyosarcoma tumors and 10 rhabdomyosarcoma cell lines, as well as 8 normal tissues. Unsupervised clustering of DNA methylation clearly distinguished the fusion-positive and fusion-negative subsets. The fusion-positive tumors showed substantially lower overall levels of methylation compared with fusion-negative tumors. Comparison with the methylation pattern of normal skeletal muscle and bone marrow indicates that fusion-negative rhabdomyosarcoma is more similar to these normal tissues compared with fusion-positive rhabdomyosarcoma, and suggests that many of the methylation differences between these subtypes arise from 'aberrant' hyper- and hypomethylation events in fusion-positive rhabdomyosarcoma. Integrative methylation and gene expression analysis revealed that methylation differences between fusion-positive and fusion-negative tumors could either be positively or negatively associated with mRNA expression. There was no significant difference in the distribution of PAX3-FOXO1-binding sites between genes with and without differential methylation. However, the finding that PAX3-FOXO1-binding sites were enriched among genes that were both differentially methylated and differentially expressed suggests that the fusion protein interacts with DNA methylation to regulate target gene expression. An 11-gene DNA methylation signature, classifying the rhabdomyosarcoma tumors into fusion-positive and fusion-negative subsets, was established and validated by pyrosequencing assays. Notably, EMILIN1 (part of the 11-gene signature) showed higher methylation and lower mRNA expression in fusion-positive compared with fusion-negative tumors, and demonstrated demethylation and re-expression in multiple fusion-positive cell lines after treatment with 5-aza-2'-deoxycytidine. In conclusion, our study demonstrates that fusion-positive and fusion-negative rhabdomyosarcoma tumors possess characteristic methylation profiles that contribute to the expression differences between these fusion subtypes. These findings indicate an important relationship between fusion status and epigenetic changes in rhabdomyosarcoma, present a novel approach for ascertaining fusion status, and may identify new therapeutic targets in rhabdomyosarcoma.

Detection of global hypermethylation in well-differentiated thyroid neoplasms by immunohistochemical (5-methylcytidine) analysis

PURPOSE

While global hypomethylation of DNA has been found in several malignancies, studies on thyroid tumours have shown controversial results using different techniques. To help resolve this issue, we assessed methylation status using two different techniques in papillary thyroid carcinomas (PTC) and follicular adenomas (FA) and carcinomas (FTC), comparing adjacent non-neoplastic thyroid tissue.

METHODS

A series of 15 FA, 18 FTC and 17 PTC were assessed by: (1) measurement of methylation levels of long interspersed nuclear elements (LINE-1) using a combined bisulfite restriction analysis polymerase chain reaction protocol and (2) immunostaining with an anti-5-methylcytidine antibody that detects methylated DNA regardless of the DNA sequence. Immunostaining was scored by image analysis.

RESULTS

Methylation levels of LINE-1 in FA, FTC and PTC were not significantly different from adjacent normal tissue. There was no significant difference in methylation levels of LINE-1 between FA, FTC and PTC (p = 0.44). By immunohistochemical staining for methylation, the 5-methylcytidine score was significantly higher in tumours than in normal tissue counterparts, for FA (p < 0.001), FTC (p = 0.04) and PTC (p = 0.02). PTC showed the highest 5-methylcytidine expression amongst all tumours which was significantly different from FTC (p = 0.015), but not FA (p = 0.09). There was no correlation in methylation level between LINE-1 and 5-methylcytidine scores for each group and overall.

CONCLUSIONS

Well-differentiated thyroid neoplasms (FA, FTC and PTC) were not found by two independent methods to undergo global hypomethylation as part of an oncogenic sequence from normal tissue to carcinoma. Instead, hypermethylation was detected in all types of tumours, implying that this epigenetic event may contribute to oncogenic development of thyroid neoplasms (both benign and malignant).

A novel 96well-formatted micro-gap plate enabling drug response profiling on primary tumour samples

Drug-based treatments are the most widely used interventions for cancer management. Personalized drug response profiling remains inherently challenging with low cell count harvested from tumour sample. We present a 96well-formatted microfluidic plate with built-in micro-gap that preserves up to 99.2% of cells during multiple assay/wash operation and only 9,000 cells needed for a single reagent test (i.e. 1,000 cells per test spot x 3 selected concentration x triplication), enabling drug screening and compatibility with conventional automated workstations. Results with MCF7 and MDA-MB-231 cell lines showed that no statistical significance was found in dose-response between the device and conventional 96-well plate control. Primary tumour samples from breast cancer patients tested in the device also showed good IC₅₀ prediction. With drug screening of primary cancer cells must consider a wide range of scenarios, e.g. suspended/attached cell types and rare/abundant cell availability, the device enables high throughput screening even for suspended cells with low cell count since the signature microfluidic cell-trapping feature ensures cell preservation in a multiple solution exchange protocol.

Methylation of promoter and expression silencing of GPX3 gene in hepatocellular carcinoma tissue

BACKGROUND/AIMS

Accumulating evidence has identified transcriptional silencing by aberrant methylation of CpG islands as a potential mechanism for the inactivation of tumor suppressor genes. The role of aberrant methylation of the GPX3 promoter in hepatocellular carcinoma (HCC) is not yet clear. We investigated the association of the status of GPX3 promoter methylation and GPX3 protein expression with the clinicopathological progression of HCC.

METHODOLOGY

Sixty HCC tumor and matched non-cancerous tissues were included in this study, and methylation was examined using MSP. GPX3 mRNA and protein levels were evaluated using RT-PCR and western blot analysis, respectively. Clinicopathological data were compiled for correlation analyses.

RESULTS

Among the 60 HCC cases, 76.7% (46/60) showed at elevated DNA methylation and displayed significantly lower levels of GPX3 mRNA and protein expression. Low or undetectable GPX3 protein expression was observed in 10 of 60 tumors. GPX3 promoter methylation was detected in 46 of 60 (76.7%) tumors, while no GPX3 gene promoter methylation was observed in the matched non-cancerous specimens. There was a negative correlation between promoter methylation and GPX3 mRNA levels (P<0.05). Analysis of clinicopathological data revealed that both mRNA and protein were significantly associated with portal tumor thrombosis, metastasis and differentiation. In additional, GPX3 methylation showed a relationship with portal tumor thrombosis, metastasis and differentiation.

CONCLUSIONS

Our results suggest that promoter methylation may be a mechanism for inactivation of GPX3, possibly leading to subsequent carcinogenesis and progression of HCC.

Integrated genomic characterization of papillary thyroid carcinoma

Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer. Here, we describe the genomic landscape of 496 PTCs. We observed a low frequency of somatic alterations (relative to other carcinomas) and extended the set of known PTC driver alterations to include EIF1AX, PPM1D, and CHEK2 and diverse gene fusions. These discoveries reduced the fraction of PTC cases with unknown oncogenic driver from 25% to 3.5%. Combined analyses of genomic variants, gene expression, and methylation demonstrated that different driver groups lead to different pathologies with distinct signaling and differentiation characteristics. Similarly, we identified distinct molecular subgroups of BRAF-mutant tumors, and multidimensional analyses highlighted a potential involvement of oncomiRs in less-differentiated subgroups. Our results propose a reclassification of thyroid cancers into molecular subtypes that better reflect their underlying signaling and differentiation properties, which has the potential to improve their pathological classification and better inform the management of the disease.