Effect of BRAF mutational status on expression profiles in conventional papillary thyroid carcinomas.

BACKGROUND: Whereas 40 % to 70 % of papillary thyroid carcinomas (PTCs) are characterized by a BRAF mutation (BRAFmut), unified biomarkers for the genetically heterogeneous group of BRAF wild type (BRAFwt) PTCs are not established yet. Using state-of-the-art technology we compared RNA expression profiles between conventional BRAFwt and BRAFmut PTCs. METHODS: Microarrays covering 36,079 reference sequences were used to generate whole transcript expression profiles in 11 BRAFwt PTCs including five micro PTCs, 14 BRAFmut PTCs, and 7 normal thyroid specimens. A p-value with a false discovery rate (FDR) < 0.05 and a fold change > 2 were used as a threshold of significance for differential expression. Network and pathway utilities were employed to interpret significance of expression data. BRAF mutational status was established by direct sequencing the hotspot region of exon 15. RESULTS: We identified 237 annotated genes that were significantly differentially expressed between BRAFwt and BRAFmut PTCs. Of these, 110 genes were down- and 127 were upregulated in BRAFwt compared to BRAFmut PTCs. A number of molecules involved in thyroid hormone metabolism including thyroid peroxidase (TPO) were differentially expressed between both groups. Among cancer-associated molecules were ERBB3 that was downregulated and ERBB4 that was upregulated in BRAFwt PTCs. Two microRNAs were significantly differentially expressed of which miR492 bears predicted functions relevant to thyroid-specific molecules. The protein kinase A (PKA) and the G protein-coupled receptor pathways were identified as significantly related signaling cascades to the gene set of 237 genes. Furthermore, a network of interacting molecules was predicted on basis of the differentially expressed gene set. CONCLUSIONS: The expression study focusing on affected genes that are differentially expressed between BRAFwt and BRAFmut conventional PTCs identified a number of molecules which are connected in a network and affect important canonical pathways. The identified gene set adds to our understanding of the tumor biology of BRAFwt and BRAFmut PTCs and contains genes/biomarkers of interest.

Comparison of microarray expression profiles between follicular variant of papillary thyroid carcinomas and follicular adenomas of the thyroid.

BACKGROUND: Follicular variant of papillary thyroid carcinoma (FVPTC) and follicular adenoma (FA) are histologically closely related tumors and differential diagnosis remains challenging. RNA expression profiling is an established method to unravel molecular mechanisms underlying the histopathology of diseases. METHODS: BRAF mutational status was established by direct sequencing the hotspot region of exon 15 in six FVPTCs and seven FAs. Whole-transcript arrays were employed to generate expression profiles in six FVPTCs, seven FAs and seven normal thyroid tissue samples. The threshold of significance for differential expression on the gene and exon level was a p-value with a false discovery rate (FDR) < 0.05 and a fold change cutoff > 2. Two dimensional average linkage hierarchical clustering was generated using differentially expressed genes. Network, pathway, and alternative splicing utilities were employed to interpret significance of expression data on the gene and exon level. RESULTS: Expression profiling in FVPTCs and FAs, all of which were negative for a BRAF mutation, revealed 55 transcripts that were significantly differentially expressed, 40 of which were upregulated and 15 downregulated in FVPTCs vs. FAs. Amongst the most significantly upregulated genes in FVPTCs were GABA B receptor, 2 (GABBR2), neuronal cell adhesion molecule (NRCAM), extracellular matrix protein 1 (ECM1), heparan sulfate 6-O-sulfotransferase 2 (HS6ST2), and retinoid X receptor, gamma (RXRG). The most significantly downregulated genes in FVPTCs included interaction protein for cytohesin exchange factors 1 (IPCEF1), G protein-coupled receptor 155 (GPR155), Purkinje cell protein 4 (PCP4), chondroitin sulfate N-acetylgalactosaminyltransferase 1 (CSGALNACT1), and glutamate receptor interacting protein 1 (GRIP1). Alternative splicing analysis detected 87 genes, 52 of which were also included in the list of 55 differentially expressed genes. Network analysis demonstrated multiple interactions for a number of differentially expressed molecules including vitamin D (1,25- dihydroxyvitamin D3) receptor (VDR), SMAD family member 9 (SMAD9), v-kit Hardy-Zuckerman 4 feline sarcoma viral oncogene homolog (KIT), and RXRG. CONCLUSIONS: This is one of the first studies using whole-transcript expression arrays to compare expression profiles between FVPTCs and FAs. A set of differentially expressed genes has been identified that contains valuable candidate genes to differentiate both histopathologically related tumor types on the molecular level.

Comprehensive survey of HRAS, KRAS, and NRAS mutations in proliferative thyroid lesions from an ethnically diverse population.

BACKGROUND: The distribution and kind of rat sarcoma viral oncogenes homolog (RAS) mutations, as well as their clinical impact on different types of thyroid lesions, vary widely among the different populations studied. We performed a comprehensive mutational survey in the highly related RAS genes HRAS, KRAS, and NRAS in a case series of proliferative thyroid lesions with known BRAF mutational status, originating from an ethnically diverse group. MATERIALS AND METHODS: Mutational hotspot regions encompassing codons 12, 13, and 61 of the RAS genes were directly sequenced in 381 cases of thyroid lesions. In addition, the putative NRAS hotspot region encompassing codon 97 was sequenced in 36 thyroid lesions. The case series included lesions of Hashimoto's thyroiditis (HT), nodular goiters, hyperplastic nodules, follicular adenomas (FAs), Hurthle cell variants of FA, papillary thyroid carcinomas (PTCs), follicular variants of PTC (FVPTCs), microcarcinomas of PTC (micro PTCs; tumor size ≤1 cm), follicular TCs (FTCs), Hurthle cell variants of FTC, and non-well-differentiated TCs (NWDTCs). RESULTS: We identified RAS mutations in 16 out of 57 (28.1%) FAs, 2 out of 8 (25%) NWDTCs, 8 out of 42 (19.0%) FVPTCs, 2 out of 10 (20.0%) FTCs, 1 out of 12 (8.3%) Hurthle cell variants of FA, 3 out of 46 (6.5%) goiters, 1 out of 18 (5.6%) hyperplastic nodules, 3 out of 56 (5.4%) micro PTCs, 2 out of 115 (1.7%) PTCs, 0 out of 7 (0%) Hurthle cell variants of FTC, and 0 out of 10 (0%) HT lesions. NRAS codon 61 mutation was the predominant form, followed by HRAS codon 61 mutation. Only three mutations affected RAS codons 12 and 13, two of which were identified in goiters. No codon 97 mutation was detected in the examined FVPTCs. An as yet undescribed deletion of KRAS codon 59 was identified in one FA. DISCUSSION: RAS mutations in our case series were commonly associated with follicular-patterned thyroid lesions. Our data suggest that FAs with a RAS mutation may constitute precursor lesions for TC with follicular histology. The newly-discovered KRAS codon 59 deletion is one of the first reported codon deletions in a RAS hotspot region.

BRAF mutations in thyroid tumors from an ethnically diverse group.

BACKGROUND: The molecular etiology of thyroid carcinoma (TC) and other thyroid diseases which may present malignant precursor lesions is not fully explored yet. The purpose of this study was to estimate frequency, type and clinicopathological value of BRAF exon 15 mutations in different types of cancerous and non-cancerous thyroid lesions originating in an ethnically diverse population. METHODS: BRAF exon 15 was sequenced in 381 cases of thyroid lesions including Hashimoto´s thyroiditis, nodular goiters, hyperplastic nodules, follicular adenomas (FA), papillary TC (PTC), follicular variant PTC (FVPTC), microcarcinomas of PTC (micro PTC; tumor size ≤ 1 cm), follicular TC (FTC), and non-well differentiated TC (non-WDTC). RESULTS: We identified BRAF mutations in one of 69 FA, 72 of 115 (63%) PTC, seven of 42 (17%) FVPTC, 10 of 56 (18%) micro PTC, one of 17 (6%) FTC, and one of eight (13%) non-WDTC. Most of the cases showed the common V600E mutation. One case each of PTC, FVPTC, and FTC harbored a K601E mutation. A novel BRAF mutation was identified in a FA leading to deletion of threonine at codon 599 (p.T599del). A rare 3-base pair insertion was detected in a stage III PTC resulting in duplication of threonine at codon 599 (p.T599dup). Patients with PTC harboring no BRAF mutation (BRAFwt) were on average younger than those with a BRAF mutation (BRAFmut) in the PTC (36.6 years vs. 43.8 years). Older age (≥ 45 years) in patients with PTC was significantly associated with tumor size ≥ 4 cm (P = 0.018), vessel invasion (P = 0.004), and distant metastasis (P = 0.001). Lymph node (LN) involvement in PTC significantly correlated with tumor size (P = 0.044), and vessel invasion (P = 0.013). Of notice, taken the whole TC group, family history of thyroid disease positively correlated with capsular invasion (P = 0.025). CONCLUSIONS: Older age is manifold associated with unfavorable tumor markers in our series. The K601E identified in a PTC, FVPTC, and FTC seems to be more distributed among different histological types of TC than previously thought. The T599del is a yet undescribed mutation and the rare T599dup has not been reported as a mutation in PTC so far.