Infinium HumanMethylation 450 BeadChip Arrays by Illumina (Illumina HM450K) are among the most popular CpG microarray platforms widely used in biological and medical research. Several recent studies highlighted the potentially confounding impact of the genomic variation on the results of methylation studies performed using Illumina’s Infinium methylation probes. However, the complexity of SNPs impact on the methylation level measurements (β values) has not been comprehensively described. In our comparative study of European and Asian populations performed using Illumina HM450K, we found that the majority of Infinium probes, which differentiated two examined groups, had SNPs in their target sequence. Characteristic tri-modal or bi-modal patterns of β values distribution among individual samples were observed for CpGs with SNPs in the first and second position, respectively. To better understand how SNPs affect methylation readouts, we investigated their impact in the context of SNP position and type, and of the Illumina probe type (Infinium I or II). Our study clearly demonstrates that SNP variation existing in the genome, if not accounted for, may lead to false interpretation of the methylation signal differences suggested by some of the Illumina Infinium probes. In addition, it provides important practical clues for discriminating between differences due to the methylation status and to the genomic polymorphisms, based on the inspection of methylation readouts in individual samples. This approach is of special importance when Illumina Infinium assay is used for any comparative population studies, whether related to cancer, disease, ethnicity where SNP frequencies differentiate the studied groups.

Distinguishing between follicular thyroid cancer (FTC) and follicular thyroid adenoma (FTA) constitutes a long-standing diagnostic problem resulting in equivocal histopathological diagnoses. There is therefore a need for additional molecular markers. To identify molecular differences between FTC and FTA, we analyzed the gene expression microarray data of 52 follicular neoplasms. We also performed a meta-analysis involving 14 studies employing high throughput methods (365 follicular neoplasms analyzed). Based on these two analyses, we selected 18 genes differentially expressed between FTA and FTC. We validated them by quantitative real-time polymerase chain reaction (qRT-PCR) in an independent set of 71 follicular neoplasms from formaldehyde-fixed paraffin embedded (FFPE) tissue material. We confirmed differential expression for 7 genes (CPQ, PLVAP, TFF3, ACVRL1, ZFYVE21, FAM189A2, and CLEC3B). Finally, we created a classifier that distinguished between FTC and FTA with an accuracy of 78%, sensitivity of 76%, and specificity of 80%, based on the expression of 4 genes (CPQ, PLVAP, TFF3, ACVRL1). In our study, we have demonstrated that meta-analysis is a valuable method for selecting possible molecular markers. Based on our results, we conclude that there might exist a plausible limit of gene classifier accuracy of approximately 80%, when follicular tumors are discriminated based on formalin-fixed postoperative material.

The molecular etiology of follicular thyroid tumors is largely unknown, rendering the diagnostics of these tumors challenging. The somatic alterations present in these tumors apart from RAS gene mutations and PAX8/PPARG translocations are not well described. To evaluate the profile of somatic alteration in follicular thyroid tumors, a total of 82 thyroid tissue samples derived from 48 patients were subjected to targeted Illumina HiSeq next generation sequencing of 372 cancer-related genes. New somatic alterations were identified in oncogenes (MDM2, FLI1), transcription factors and repressors (MITF, FLI1, ZNF331), epigenetic enzymes (KMT2A, NSD1, NCOA1, NCOA2), and protein kinases (JAK3, CHEK2, ALK). Single nucleotide and large structural variants were most and least frequently identified, respectively. A novel translocation in DERL/COX6C was detected. Many somatic alterations in non-coding gene regions with high penetrance were observed. Thus, follicular thyroid tumor somatic alterations exhibit complex patterns. Most tumors contained distinct somatic alterations, suggesting previously unreported heterogeneity.

TERT promoter (TERTp) mutations are important factors in papillary thyroid carcinomas (PTCs). They are associated with tumor aggressiveness, recurrence, and disease-specific mortality and their use in risk stratification of PTC patients has been proposed. In this study we investigated the prevalence of TERTp mutations in a cohort of Polish patients with PTCs and the association of these mutations with histopathological factors, particularly in coexistence with the BRAF V600E mutation. A total of 189 consecutive PTC specimens with known BRAF mutational status were evaluated. TERTp mutations were detected in 8.5% of cases (16/189) with the C228T mutation being the most frequent. In six of the PTC specimens (3.2%), four additional TERTp alterations were found, which included one known polymorphism (rs2735943) and three previously unreported alterations. The association analysis revealed that the TERTp hotspot mutations were highly correlated with the presence of the BRAF V600E mutation and their coexistence was significantly associated with gender, advanced patient age, advanced disease stage, presence of lymph node metastases, larger tumor size, and tumor-capsule infiltration. While correlations were identified, the possibility of TERTp mutations being key molecular modulators responsible for PTC aggressiveness requires further studies.

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.