Effect of BRAFmutational status on expression profiles in conventional papillary thyroid carcinomas
BMC Genomics volume 16, Article number: S6 (2015)
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.
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.
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.
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.
Over the last decades the incidence rate of thyroid cancer increases worldwide . In Saudi Arabia, thyroid carcinoma (TC) is considered the second most common cancer in young women . About 80% of all TCs are PTCs. The majority of PTCs are histologically classified as conventional PTCs. The follicular variant of PTC (FVPTC) represents the largest subtype and accounts for about 30% of all PTCs . Minor and rare subtypes include Hurthle cell variant PTC and insular PTC which bears an aggressive clinical behavior . Conventional PTCs are characterized on the molecular level by a moderate to high frequency of BRAF mutations (40 % - 70 %) that distinguishes them from FVPTCs (10 % - 20 %) .
BRAF is a cytoplasmic receptor serine/threonine kinase and a key molecule in the mitogen activated protein kinase (MAPK) pathway. BRAF is mutated in diverse human malignancies although frequency and clinical presentation varies considerably between different types of cancers . Over 90 % of all BRAF mutations are a valine by glycine substitution at codon 600 (V600E) in exon 15. Other BRAF mutations affect commonly codons adjunct to codon 600. Although the impact of BRAF mutations in PTC is controversially discussed, many studies found an association of BRAFmut PTCs with unfavorable clinical features including larger tumor size, advanced tumor stage, vessel invasion, capsular invasion, tumor extension, higher risk for lymph node (LN) involvement, distant metastasis, and poor prognosis [7–10]. Consistent with this, patients with a BRAFmut PTC are considerably older than those with a BRAFwt PTC [5, 10]. Whereas a BRAF mutation represents are valuable target for molecular therapy in advanced solid tumors including PTCs, molecular profiles of BRAFwt PTCs are less known and genetic screening for valuable target genes is primarily limited to research studies . The major deregulated key genes in the BRAFwt group are RET and RAS. RAS consists of the highly related genes for HRAS, KRAS, and NRAS. Within the MAPK pathway the RAS molecules transmit signals to the downstream target BRAF. The most common RET/PTC fusions are paracentric fusions with the gene entitled coiled-coil domain containing 6 (CCDC6) contributing to ~80 % and with the nuclear receptor coactivator 4 (NCOA4) contributing to ~10 % of all known RET/PTC rearrangements. The frequency of RAS mutations and RET/PTC rearrangements differs between the populations studied and depends in part on the inclusion/exclusion criteria for the different histological PTC subtypes . RAS mutations and RET rearrangements are unlikely to act as molecular drivers for onset of malignancy as they are already present in benign thyroid neoplasms [13, 14]. This distinguishes them from BRAF mutations which are virtually absent in precursor lesions of PTCs.
Until now only a few studies compared expression profiles between BRAFwt and BRAFmut PTCs using array technologies [15–18]. The relevance of expressional screening in PTC according to their BRAF mutational status is in part related to the different clinical behavior of both PTC groups with the necessity to identify appropriate biomarkers and in part related to the different tumor biology of both groups which is not thoroughly understood . In our study we took advantage of current state-of-the-art technology to screen and analyze a case series of BRAFwt and BRAFmut PTCs for detecting new biomarkers which could become useful to distinguish both groups on the molecular level. We did not include FVPTCs and other smaller histological subtypes of PTC in our screening to minimize expressional bias which might be related to a different histology.
We examined 25 specimens from PTCs and seven normal thyroid samples (TN) from patients who were treated surgically in the period between November 2008 and February 2013 at the King Abdulaziz University Hospital, Jeddah, and the King Faisal Specialist Hospital & Research Center (KFSH&RC), Jeddah. In two BRAFmut PTCs, specimens were derived from a recurrence or a local metastasis and in one BRAFwt PTC from an LN metastasis. Normal thyroid specimens were derived from histopathologically unaffected normal thyroid tissue in the course of lobo- or thyroidetomies of thyroid lesions (4 goiters, 1 hyperthyroiditis, 1 PTC, and 1 FVPTC). Histopathological diagnosis and staging of thyroid lesions was performed by an experienced oncologic pathologist (JM) according to established criteria [19, 20]. Five BRAFwt PTCs were classified as micro PTCs (≤ 1 cm). This study was approved by the Research Ethics Committee of the King Abdulaziz University, Faculty of Medicine, #358-10, and the Institutional Review Board of the KFSH&RC, #IRB2010-07, and included written informed consent provided by the participants.
DNA extraction and BRAFmutational screening
Genomic DNA from fresh-frozen samples was extracted using the QIAmp DNA tissue kit (Qiagen, Hilden, Germany). DNA concentration was measured with the Nanodrop ND-1000 spectrophotometer (Thermo Scientific, Wilmington, DE). Screening of the BRAF mutational hotspot region of exon 15 was performed as described earlier involving direct sequencing of PCR products spanning the region .
RNA sample and array processing
Total RNA was extracted from freshly preserved thyroid tissue specimens using the Qiagen RNeasy Mini Kit (Qiagen, Hilden, Germany) including an on-column DNAse treatment according to manufacturer’s recommendations. Quality of the purified RNA was verified on an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA). RNA integrity number for all evaluated samples was at least 5.0. RNA concentrations were determined using a NanoDrop ND-1000 spectrophotometer. Samples containing each 250 ng of RNA were processed using the Ambion WT Expression Kit (Life Technologies, Austin, TX) and the GeneChip WT Terminal Labeling and Controls Kit (Affymetrix, Santa Clara, CA) according to the manufacturers` recommendations. Fragmentation and endlabeling of samples were monitored by electrophoresis on 3 % agarose gels. Affymetrix GeneChip hybridization, wash and stain kits were utilized in subsequent processing steps. Hybridization mixtures containing each 5500 ng of cDNA were hybridized at 45°C for 17 hrs and 60 rpm to Affymetrix Human Gene 1.0 ST GeneChip arrays. Subsequent to wash and staining at the GeneChip Fluidics Station 450, the arrays were scanned with the GeneChip Scanner 3000 7G. Probe cell intensity data (CEL files) were generated using the GeneChip Command Console Software (AGCC). Human Gene 1.0 ST GeneChip arrays interrogate in total with a set of 764,885 probes 36.079 reference sequences (NCBI build 36).
Gene Expression Analysis
CEL files were imported to Partek Genomics Suite version 6.6 (Partek Inc., MO) and a log-transformed data set of robust multi-array averaged (RMA), background-adjusted, and normalized values was generated. Non-annotated genes and multiple transcripts generated from the same gene were excluded from further analysis. Principal component analysis (PCA) was performed to assess quality as well as overall variance in gene expression between groups of samples. Analysis of Variance (ANOVA) was applied to generate a list of differentially expressed genes using a p-value with a false discovery rate (FDR) (Step up method) < 0.05 and a fold change > 2.0. Two dimensional average linkage hierarchical clustering was performed using Spearman’s correlation as a similarity matrix. The generated array data set complies with MIAME  and was submitted to NCBI’s Gene Expression Omnibus (GEO), accession number GSE54958.
Functional network and pathway analysis
To define molecular networks and canonical pathways in differentially regulated gene sets, pathway analyses were performed by using Ingenuity Pathways Analysis (IPA) software (Ingenuity Systems, Redwood City, CA). For this purpose, statistically differentially expressed genes and their corresponding probe set ID, gene symbol as clone identifier, p-value and fold change were imported into IPA. The program identifies with its functional algorithms those biological functions, interacting drugs and/or diseases that are most significantly related to a differentially expressed gene set. The canonical pathway analysis identifies pathways that are most significantly related to the data set.
Demographic data of patients and histopathological criteria of BRAFwt and BRAFmut PTC are listed in Table 1. A gender shift towards females was observed in the BRAFwt group and mean age was considerable lower in the BRAFwt than in the BRAFmut group (30.9 years vs. 45.9 years). Histopathological criteria including tumor size, LN involvement, tumor focality and tumor stage were comparably more unfavorable in BRAFmut PTCs (Table 1).
Expression BRAFwtvs. BRAFmutPTCs
Employing whole-transcript microarrays (HuGene 1.0 ST) we compared expression profiles of 11 BRAFwt PTCs with 14 BRAFmut PTCs. Seven TN specimens were used as a reference for normal thyroid tissue expression which allowed us to identify differentially expressed genes between both PTC groups and TN samples. Three-D presentation of the PCA displays clustering of BRAFwt PTCs, BRAFmut PTCs and TN samples (Figure 1). We identified, after excluding non-annotated genes and multiple transcriptional isoforms, 237 annotated genes that were significantly differentially expressed (p-value with FDR <0.05 and a fold change >2) between BRAFwt and BRAFmut PTCs (Additional file 1). Of these genes, 127 were up-, and 110 were downregulated in BRAFwt compared to BRAFmut PTCs. The most significantly upregulated genes in BRAFwt were inositol 1,4,5-triphosphate receptor, type 1 (ITPR1), hepatic leukemia factor (HLF), potassium voltage-gated channel, shaker-related subfamily, beta member 1 (KCNAB1), engulfment and cell motility 1 (ELMO1), Rho GTPase activating protein 24 (ARHGAP24), thyroid peroxidase (TPO) and solute carrier family 4, sodium bicarbonate cotransporter, member 4 (SLC4A4). The most significantly downregulated genes in BRAFwt PTCs were dendrocyte expressed seven transmembrane protein (DCSTAMP), ladinin 1 (LAD1), keratin 19 (KRT19), chromosome 19 open reading frame 33 (C19orf33), poliovirus receptor-related 4 (PVRL4), EPH receptor A10 (EPHA10), and TBC1 domain family, member 2 (TBC1D2). The list of 237 genes contained one microRNA (MIR492) that was downregulated and one microRNA (MIR32) that was upregulated in BRAFwtvs. BRAFmut PTCs. Additional file 1 contains also the cellular/extracellular location, function of gene products as well as a selection of drugs known to interact with a gene product. Hierarchical cluster analysis of the 237 genes separates both PTC groups according to their differential gene expression (Figure 2). TN samples were included in this data set to provide normal expression values.
Expression PTCs vs. TN samples
Comparison of genes differentially expressed between BRAFwt and TN samples revealed a set of 8249 genes and between BRAFmut and TN samples a set of 8836 genes. To identify genes of interests, e.g. conceivable immunohistochemistry markers, we selected those genes which are differentially expressed between all three groups (Table 2). In sum, 32 genes met the statistical criteria and the vast majority revealed the comparatively highest expression in BRAFmut PTCs.
Gene networks and canonical pathways
Most significant network functions identified by IPA algorithms and associated with the set of 237 differentially expressed genes were involved in cell signaling, cancer, and cellular development (Figure 3). This comprehensive network includes 12 molecules which are overexpressed and 16 molecules which are underexpressed in BRAFwt compared to BRAFmut PTCs. The most significantly associated canonical pathway related to the differentially expressed set of 237 genes is the proteinkinase A (PKA) signaling cascade (Additional file 2). The PKA pathway is involved in second messenger signaling and it is stimulated by upstream cascades including the G protein-coupled receptor pathway (Figure 4). The most significant pathways with thyroid specific functions were the thyroid hormone metabolic and thyronamine/iodothyronamine metabolic pathways.
We performed one of the first studies using whole transcript microarrays to compare expression profiles solely in conventional BRAFwt and BRAFmut PTCs. Comparison to studies including different histological subtypes of PTCs may result in detecting a lower number of common genes. An expression array study including histological variants of PTC identified, on basis of enhanced stringent threshold criteria, over 80 up- and downregulated genes in the BRAFmut group in comparison to PTCs with either a RAS mutation or a RET/PTC rearrangement  and the 40 most up- and downregulated genes have an overlap of ~40 % to our list of 237 differentially expressed genes.
Selected differentially expressed molecules
DCSTAMP, also known as TM7SF4, has been previously identified as one of the most overexpressed genes in BRAFmut PTC compared to BRAFwt PTC as well in PTC with undetermined mutational status compared to normal thyroid tissue [15, 22, 23]. It has been supposed that DCSTAMP expression is an immune response related to BRAFmut tumors . DCSTAMP contains signature motifs owned by the family of transmembrane serine proteases and it exhibits degradation activity against extracellular matrix proteins. One of the functions of the hepatic leukemia factor (HLF) as a transcription factor is to mediate thyroid hormone activation from the thyroid hormone receptor/retinoid X receptor heterodimer to the hypoxia-inducible factor (HIF-1alpha) . The function of C19orf33, also known as H2RSP (Hepatocyte growth factor activator inhibitor type 2-related small protein) is virtually unknown. An enhanced expression with higher levels in LN positive tumors has been observed for C19orf33 in colorectal adenocarcinoma cells at the invasive front . The voltage gated channel molecule KCNAB1 exhibits diverse functions including electrolyte transport, and insulin secretion. Downregulation of KCNAB1 expression has been identified in follicular thyroid carcinomas compared to benign follicular adenomas . ELMO1 functionally interacts with dedicator of cytokinesis 1 (DOCK1) that promotes Rac guanine exchange factor (GEF) activity for Rac proteins of the Rho GTPases family. GTP-loaded Rac proteins initiate downstream pathways that promote cell elongation, migration, and cytoskeleton remodeling [27, 28]. The active ELMO1/DOCK1 complex is anchored via phosphoinositides to the membrane. In our study, the DOCK1 related DOCK5 gene was identified as a 2-fold upregulated molecule in BRAFwt compared to BRAFmut PTCs. The leucine-rich repeat and immunoglobulin-like domain-containing nogo receptor-interacting protein 2 (LINGO2) encodes a single-pass type I membrane protein that is primarily expressed during development in cells adjunct to the epithelial lining of the olfactory pit and in adult brain . In our study, LINGO2 was downregulated in BRAFwt PTCs and even more in BRAFmut PTCs compared to TN samples which may imply a tumor suppressor function for this molecule. Higher expression levels of the ERK1/2-specific cytoplasmic dual specificity phosphatase 6 (DUSP6) in comparison to benign and normal thyroid cells has been previously associated with PTC, especially with advanced thyroid carcinomas [30, 31]. In our study, DUSP6 was 2.1 times higher expressed in BRAFmut than in BRAFwt PTCs assuming that the ERK1/2 related pathway is frequently more utilized in the BRAFmut group. MET overexpression in thyroid cancer has been identified in a number of studies and this molecule was 2.3-fold higher expressed in our study in BRAFmut compared to BRAFwt PTCs which is in accordance with another survey  (Table 2).
Second messenger molecules
ITPR1 is an intracellular receptor for inositol 1,4,5-trisphosphate (IP3) and implicated in the thyroid hormone synthesis pathway. The receptor mediates calcium release from the endoplasmic reticulum upon stimulation by IP3. Downregulation of ITPR1 has been demonstrated in thyroid cancer in comparison to non-malignant thyroid tissue in a number of studies [32–35]. Among the members of the phospholipase C (PLC) family, PLCH1 was upregulated in BRAFwt PTCs whereas PLCD3 was downregulated. PLC molecules hydrolyze phosphatidylinositol 4,5-bisphosphate to generate the second messengers diacylglycerol and IP3.
Thyroid hormone pathway molecules
LAD1 has been identified as an overexpressed gene in BRAFmut thyroid carcinomas compared to those with a RET/PTC rearrangement . Lad1 expression is regulated by the glucocorticoid receptor (GR) and requires for induction the GR coactivators thyroid hormone receptor associated protein 220 (MED1) and the thyroid hormone receptor associated protein 170 (MED14) . Downregulation of TPO has been considered in a number of studies comparing thyroid carcinomas with benign thyroid tumors or normal thyroid samples [32, 33, 37, 38]. However, downregulation of TPO and the sodium iodine symporter genes as been previously associated with BRAFmut PTCs in comparison to PTCs with a RET/PTC rearrangement . This is in line with observations that BRAFmut tumors are refractory for radioactive iodine ablation due to downregulation of thyroid hormone biosynthesis pathways . TPO was in our study 13-fold downregulated in BRAFmut compared to BRAFwt PTCs. In addition, other genes involved in thyroid hormone biosynthesis pathway including the solute carrier family-5 member-8 (SLC5A8), solute carrier family 26 (anion exchanger), member 4 (SLC26A4), deiodinase, iodothyronine, type I (DIO1) and deiodinase, iodothyronine, type II (DIO2) were downregulated in our set of BRAFmut PTCs. Other members of the SLC family that were downregulated in BRAFmut PTCs were SLC1A3, SLC4A4, SLC16A2, and SLC26A7 whereas SLC22A31, SLC30A2, and SLC34A2 were downregulated in BRAFwt PTCs.
ERBB3 and ERRB4genes
Two of the four structurally related receptor tyrosine kinases ErbB, namely ERBB3 and ERBB4, were differentially expressed in the two PTC groups, i.e. ERBB3 was downregulated and ERBB4 upregulated in BRAFwt PTCs (Additional file 1). Of notice, ERBB4 was also lower expressed in BRAFmut PTCs compared to TN samples (p = 0.0035) (data not shown). Oncogenic functions of ERBB3 and ERBB4, which can form heterodimers and signal through the PI3K/AKT signaling pathway have not been elucidated yet in detail in relation to the BRAF mutational status in thyroid cancer [41–43]; however, decreased expression of ERBB4 in PTC vs. normal thyroid tissue has been demonstrated in a RT-PCR study whereas ERBB2 and ERBB3 expression was shown to be increased . On the protein level, a tissue microarray study in proliferative thyroid lesions found a correlation of ERBB3 expression with LN metastasis whereas ERBB4 expression correlated with lower tumor stage . A possible link of ERBB3 to apoptosis can be deduced from a functional in vivo study wherein deletion of ERBB3 in mouse intestinal epithelium induced tumor-specific cell death .
The most significantly downregulated microRNA in BRAFwt compared to BRAFmut PTCs was mir492 and the most upregulated one was mir32. Overexpression of mir492 has been previously linked to progressive hepatoblastoma and tumorigenesis of retinoblastoma [47, 48]. Of notice, mir492 is processed from the KRT19 gene  and both are higher expressed. A possible target for miR492 is the 3`UTR of KRT19  which is 4.4-fold higher expressed in BRAFmut PTCs compared to BRAFwt PTCs suggesting an accumulation of the KRT19 mRNA. Another possible target of miR492 is the thyroid hormone receptor-associated protein 3 (THRAP3) mRNA that harbors two predicted target sites for miR492  and that is significantly downregulated in our case series by approximately 3-fold in both, BRAFwt and BRAFmut PTCs compared to TN samples. Upregulation of mir32 in thyroid cancer vs. benign thyroid tumors has been detected in a microarray study; however functional implication of this microRNA in PTC is not known yet .
In summary, our microarray expression study provides a detailed overview of differentially expressed genes, networks, and pathways between BRAFwt and BRAFmut PTCs that gain interest for basic molecular genetics and translational studies in PTCs.
- BRAF wt :
BRAF wild type
- BRAF mut BRAF :
follicular variant of PTC
papillary thyroid carcinoma
Pellegriti G, Frasca F, Regalbuto C, Squatrito S, Vigneri R: Worldwide increasing incidence of thyroid cancer: update on epidemiology and risk factors. Journal of cancer epidemiology. 2013, 2013 (965212):
Hussain F, Iqbal S, Mehmood A, Bazarbashi S, Elhassan T, Chaudhri N: Incidence of thyroid cancer in the Kingdom of Saudi Arabia, 2000-2010. Hematol Oncol Stem Cell Ther. 2013, 6 (2): 58-64. 10.1016/j.hemonc.2013.05.004.
Yu XM, Schneider DF, Leverson G, Chen H, Sippel RS: Follicular Variant of Papillary Thyroid Carcinoma is a Unique Clinical Entity: A Population-Based Study of 10,740 Cases. Thyroid : official journal of the American Thyroid Association. 2013, 23 (10): 1263-1268. 10.1089/thy.2012.0453.
Roman S, Sosa JA: Aggressive variants of papillary thyroid cancer. Current opinion in oncology. 2013, 25: 33-38. 10.1097/CCO.0b013e32835b7c6b.
Schulten HJ, Salama S, Al-Mansouri Z, Alotibi R, Al-Ghamdi K, Al-Hamour OA, Sayadi H, Al-Aradati H, Al-Johari A, Huwait E, et al: BRAF mutations in thyroid tumors from an ethnically diverse group. Hereditary cancer in clinical practice. 2012, 10: 10-10.1186/1897-4287-10-10.
El-Osta H, Falchook G, Tsimberidou A, Hong D, Naing A, Kim K, Wen S, Janku F, Kurzrock R: BRAF mutations in advanced cancers: clinical characteristics and outcomes. PloS one. 2011, 6: e25806-10.1371/journal.pone.0025806.
Lee JH, Lee ES, Kim YS: Clinicopathologic significance of BRAF V600E mutation in papillary carcinomas of the thyroid: a meta-analysis. Cancer. 2007, 110: 38-46. 10.1002/cncr.22754.
Tufano RP, Teixeira GV, Bishop J, Carson KA, Xing M: BRAF mutation in papillary thyroid cancer and its value in tailoring initial treatment: a systematic review and meta-analysis. Medicine. 2012, 91: 274-286. 10.1097/MD.0b013e31826a9c71.
Xing M, Alzahrani AS, Carson KA, Viola D, Elisei R, Bendlova B, Yip L, Mian C, Vianello F, Tuttle RM, et al: Association between BRAF V600E mutation and mortality in patients with papillary thyroid cancer. JAMA : the journal of the American Medical Association. 2013, 309: 1493-1501. 10.1001/jama.2013.3190.
Kim TH, Park YJ, Lim JA, Ahn HY, Lee EK, Lee YJ, Kim KW, Hahn SK, Youn YK, Kim KH, et al: The association of the BRAF(V600E) mutation with prognostic factors and poor clinical outcome in papillary thyroid cancer: a meta-analysis. Cancer. 2012, 118: 1764-1773. 10.1002/cncr.26500.
Morris V, Kopetz S: BRAF inhibitors in clinical oncology. F1000prime reports. 2013, 5: 11-
Patel KN, Singh B: Genetic considerations in thyroid cancer. Cancer control: journal of the Moffitt Cancer Center. 2006, 13 (2): 111-118.
Schulten HJ, Salama S, Al-Ahmadi A, Al-Mansouri Z, Mirza Z, Al-Ghamdi K, Al-Hamour OA, Huwait E, Gari M, Al-Qahtani MH, Al-Maghrabi J: Comprehensive Survey of HRAS, KRAS, and NRAS Mutations in Proliferative Thyroid Lesions from An Ethnically Diverse Population. Anticancer Res. 2013, 33: 4779-4784.
Ishizaka Y, Kobayashi S, Ushijima T, Hirohashi S, Sugimura T, Nagao M: Detection of retTPC/PTC transcripts in thyroid adenomas and adenomatous goiter by an RT-PCR method. Oncogene. 1991, 6: 1667-1672.
Giordano TJ, Kuick R, Thomas DG, Misek DE, Vinco M, Sanders D, Zhu Z, Ciampi R, Roh M, Shedden K, et al: Molecular classification of papillary thyroid carcinoma: distinct BRAF, RAS, and RET/PTC mutation-specific gene expression profiles discovered by DNA microarray analysis. Oncogene. 2005, 24: 6646-6656. 10.1038/sj.onc.1208822.
Durand S, Ferraro-Peyret C, Joufre M, Chave A, Borson-Chazot F, Selmi-Ruby S, Rousset B: Molecular characteristics of papillary thyroid carcinomas without BRAF mutation or RET/PTC rearrangement: relationship with clinico-pathological features. Endocr Relat Cancer. 2009, 16: 467-481. 10.1677/ERC-08-0081.
Frattini M, Ferrario C, Bressan P, Balestra D, De Cecco L, Mondellini P, Bongarzone I, Collini P, Gariboldi M, Pilotti S, et al: Alternative mutations of BRAF, RET and NTRK1 are associated with similar but distinct gene expression patterns in papillary thyroid cancer. Oncogene. 2004, 23: 7436-7440. 10.1038/sj.onc.1207980.
Rusinek D, Szpak-Ulczok S, Jarzab B: Gene expression profile of human thyroid cancer in relation to its mutational status. Journal of molecular endocrinology. 2011, 47: R91-103. 10.1530/JME-11-0023.
Thyroid. AJCC Cancer Staging Manual. Edited by: Edge SB BD, Compton CC, Fritz AG, Greene FL, Trotti A. 2010, New York: Springer, 87-96.
DeLellis RA, Williams ED: Thyroid and parathyroid tumors. Pathology and genetics of tumours of endocrine organs (IARC WHO classification of tumours). Edited by: DeLellis RA, Lloyd RV, Heitz PU, Eng C. 2004, Lyons: IARC Press, 51-56.
Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, et al: Minimum information about a microarray experiment (MIAME)-toward standards for microarray data. Nature genetics. 2001, 29: 365-371. 10.1038/ng1201-365.
Kim HS, Kim do H, Kim JY, Jeoung NH, Lee IK, Bong JG, Jung ED: Microarray analysis of papillary thyroid cancers in Korean. The Korean journal of internal medicine. 2010, 25: 399-407. 10.3904/kjim.2010.25.4.399.
Siraj AK, Bavi P, Abubaker J, Jehan Z, Sultana M, Al-Dayel F, Al-Nuaim A, Alzahrani A, Ahmed M, Al-Sanea O, et al: Genome-wide expression analysis of Middle Eastern papillary thyroid cancer reveals c-MET as a novel target for cancer therapy. J Pathol. 2007, 213: 190-199. 10.1002/path.2215.
Otto T, Fandrey J: Thyroid hormone induces hypoxia-inducible factor 1alpha gene expression through thyroid hormone receptor beta/retinoid x receptor alpha-dependent activation of hepatic leukemia factor. Endocrinology. 2008, 149: 2241-2250. 10.1210/en.2007-1238.
Uchiyama S, Itoh H, Naganuma S, Nagaike K, Fukushima T, Tanaka H, Hamasuna R, Chijiiwa K, Kataoka H: Enhanced expression of hepatocyte growth factor activator inhibitor type 2-related small peptide at the invasive front of colon cancers. Gut. 2007, 56: 215-226. 10.1136/gut.2005.084079.
Pfeifer A, Wojtas B, Oczko-Wojciechowska M, Kukulska A, Czarniecka A, Eszlinger M, Musholt T, Stokowy T, Swierniak M, Stobiecka E, et al: Molecular differential diagnosis of follicular thyroid carcinoma and adenoma based on gene expression profiling by using formalin-fixed paraffin-embedded tissues. BMC medical genomics. 2013, 6: 38-10.1186/1755-8794-6-38.
Hanawa-Suetsugu K, Kukimoto-Niino M, Mishima-Tsumagari C, Akasaka R, Ohsawa N, Sekine S, Ito T, Tochio N, Koshiba S, Kigawa T, et al: Structural basis for mutual relief of the Rac guanine nucleotide exchange factor DOCK2 and its partner ELMO1 from their autoinhibited forms. Proceedings of the National Academy of Sciences of the United States of America. 2012, 109: 3305-3310. 10.1073/pnas.1113512109.
Patel M, Pelletier A, Cote JF: Opening up on ELMO regulation: New insights into the control of Rac signaling by the DOCK180/ELMO complex. Small GTPases. 2011, 2: 268-275. 10.4161/sgtp.2.5.17716.
Haines BP, Rigby PW: Expression of the Lingo/LERN gene family during mouse embryogenesis. Gene Expr Patterns. 2008, 8: 79-86. 10.1016/j.modgep.2007.10.003.
Lee JU, Huang S, Lee MH, Lee SE, Ryu MJ, Kim SJ, Kim YK, Kim SY, Joung KH, Kim JM, et al: Dual specificity phosphatase 6 as a predictor of invasiveness in papillary thyroid cancer. European journal of endocrinology / European Federation of Endocrine Societies. 2012, 167: 93-101. 10.1530/EJE-12-0010.
Degl'Innocenti D, Romeo P, Tarantino E, Sensi M, Cassinelli G, Catalano V, Lanzi C, Perrone F, Pilotti S, Seregni E, et al: DUSP6/MKP3 is overexpressed in papillary and poorly differentiated thyroid carcinoma and contributes to neoplastic properties of thyroid cancer cells. Endocr Relat Cancer. 2013, 20: 23-37. 10.1530/ERC-12-0078.
Griffith OL, Melck A, Jones SJ, Wiseman SM: Meta-analysis and meta-review of thyroid cancer gene expression profiling studies identifies important diagnostic biomarkers. J Clin Oncol. 2006, 24: 5043-5051. 10.1200/JCO.2006.06.7330.
Huang Y, Prasad M, Lemon WJ, Hampel H, Wright FA, Kornacker K, LiVolsi V, Frankel W, Kloos RT, Eng C, et al: Gene expression in papillary thyroid carcinoma reveals highly consistent profiles. Proceedings of the National Academy of Sciences of the United States of America. 2001, 98: 15044-15049. 10.1073/pnas.251547398.
Finley DJ, Zhu B, Barden CB, Fahey TJ: Discrimination of benign and malignant thyroid nodules by molecular profiling. Ann Surg. 2004, 240: 425-436. 10.1097/01.sla.0000137128.64978.bc. discussion 436-427
Barden CB, Shister KW, Zhu B, Guiter G, Greenblatt DY, Zeiger MA, Fahey TJ: Classification of follicular thyroid tumors by molecular signature: results of gene profiling. Clinical cancer research : an official journal of the American Association for Cancer Research. 2003, 9: 1792-1800.
Wang Z, Chen W, Kono E, Dang T, Garabedian MJ: Modulation of glucocorticoid receptor phosphorylation and transcriptional activity by a C-terminal-associated protein phosphatase. Molecular endocrinology. 2007, 21: 625-634. 10.1210/me.2005-0338.
Prasad NB, Somervell H, Tufano RP, Dackiw AP, Marohn MR, Califano JA, Wang Y, Westra WH, Clark DP, Umbricht CB, et al: Identification of genes differentially expressed in benign versus malignant thyroid tumors. Clinical cancer research : an official journal of the American Association for Cancer Research. 2008, 14: 3327-3337. 10.1158/1078-0432.CCR-07-4495.
Tanaka T, Umeki K, Yamamoto I, Sugiyama S, Noguchi S, Ohtaki S: Immunohistochemical loss of thyroid peroxidase in papillary thyroid carcinoma: strong suppression of peroxidase gene expression. J Pathol. 1996, 179 (1): 89-94. 10.1002/(SICI)1096-9896(199605)179:1<89::AID-PATH546>3.0.CO;2-R.
Romei C, Ciampi R, Faviana P, Agate L, Molinaro E, Bottici V, Basolo F, Miccoli P, Pacini F, Pinchera A, Elisei R: BRAFV600E mutation, but not RET/PTC rearrangements, is correlated with a lower expression of both thyroperoxidase and sodium iodide symporter genes in papillary thyroid cancer. Endocr Relat Cancer. 2008, 15: 511-520. 10.1677/ERC-07-0130.
Xing M: Molecular pathogenesis and mechanisms of thyroid cancer. Nature reviews Cancer. 2013, 13: 184-199. 10.1038/nrc3431.
Kondo T, Ezzat S, Asa SL: Pathogenetic mechanisms in thyroid follicular-cell neoplasia. Nature reviews Cancer. 2006, 6: 292-306. 10.1038/nrc1836.
Mincione G, Di Marcantonio MC, Tarantelli C, D'Inzeo S, Nicolussi A, Nardi F, Donini CF, Coppa A: EGF and TGF-beta1 Effects on Thyroid Function. Journal of thyroid research. 2011, 2013 (431718):
Montero-Conde C, Ruiz-Llorente S, Dominguez JM, Knauf JA, Viale A, Sherman EJ, Ryder M, Ghossein RA, Rosen N, Fagin JA: Relief of feedback inhibition of HER3 transcription by RAF and MEK inhibitors attenuates their antitumor effects in BRAF-mutant thyroid carcinomas. Cancer discovery. 2013, 3: 520-533. 10.1158/2159-8290.CD-12-0531.
Kato S, Kobayashi T, Yamada K, Nishii K, Sawada H, Ishiguro H, Itoh M, Funahashi H, Nagasaka A: Expression of erbB receptors mRNA in thyroid tissues. Biochim Biophys Acta. 2004, 1673: 194-200. 10.1016/j.bbagen.2004.04.016.
Wiseman SM, Griffith OL, Melck A, Masoudi H, Gown A, Nabi IR, Jones SJ: Evaluation of type 1 growth factor receptor family expression in benign and malignant thyroid lesions. American journal of surgery. 2008, 195: 667-673. 10.1016/j.amjsurg.2008.01.008. discussion 673
Lee D, Yu M, Lee E, Kim H, Yang Y, Kim K, Pannicia C, Kurie JM, Threadgill DW: Tumor-specific apoptosis caused by deletion of the ERBB3 pseudo-kinase in mouse intestinal epithelium. The Journal of clinical investigation. 2009, 119: 2702-2713. 10.1172/JCI36435.
von Frowein J, Pagel P, Kappler R, von Schweinitz D, Roscher A, Schmid I: MicroRNA-492 is processed from the keratin 19 gene and up-regulated in metastatic hepatoblastoma. Hepatology. 2011, 53: 833-842. 10.1002/hep.24125.
Zhao JJ, Yang J, Lin J, Yao N, Zhu Y, Zheng J, Xu J, Cheng JQ, Lin JY, Ma X: Identification of miRNAs associated with tumorigenesis of retinoblastoma by miRNA microarray analysis. Child's nervous system : ChNS : official journal of the International Society for Pediatric Neurosurgery. 2009, 25: 13-20. 10.1007/s00381-008-0701-x.
Hsu SD, Chu CH, Tsou AP, Chen SJ, Chen HC, Hsu PW, Wong YH, Chen YH, Chen GH, Huang HD: miRNAMap 2.0: genomic maps of microRNAs in metazoan genomes. Nucleic acids research. 2008, 36: D165-169.
Kitano M, Rahbari R, Patterson EE, Xiong Y, Prasad NB, Wang Y, Zeiger MA, Kebebew E: Expression profiling of difficult-to-diagnose thyroid histologic subtypes shows distinct expression profiles and identify candidate diagnostic microRNAs. Annals of surgical oncology. 2011, 18: 3443-3452. 10.1245/s10434-011-1766-4.
We thank Nadia Bagtian, Ohoud Subhi, and Shireen Hussain for excellent technical assistance. This study was supported by King Abdulaziz City for Science and Technology (KACST) grants 13-BIO2289-03 and 09-BIO820-03.
Publication charges for this article have been funded by the Center of Excellence in Genomic Medicine Research (CEGMR), King Abdulaziz University, Jeddah, Kingdom of Saudi Arabia.
This article has been published as part of BMC Genomics Volume 16 Supplement 1, 2015: Selected articles from the 2nd International Genomic Medical Conference (IGMC 2013): Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/16/S1
The authors declare that they have no competing interests.
JM, EH, MG, and MHQ made substantial contributions to conception and design of the study. RA, AA, and MA processed expression arrays, performed BRAF mutational analysis, and were involved in data interpretation. KG, FM, and OAH were responsible for surgeries, oversight of clinical databases and contributed to the conception and design of the study. JM performed histological examinations. SK and HJS performed data analysis. HJS had general oversight of the study. HJS, JM, and MHQ interpreted data and drafted the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: The 237 most differentially expressed genes in BRAFwt vs. BRAFmut papillary thyroid carcinomas.(DOCX 47 KB)
Additional file 2: The canonical PKA pathway is a second messenger cascade and involved in diverse functions as growth, development, metabolism, DNA replication/recombination, DNA repair and cellular organization. A number of molecules including members of PKA, ryanodine receptors (RYR), inositol trisphosphate receptors (IP3R), and lymphoid enhancing factors/T-cell factors (TCF/LEF) are upregulated (red) and number of molecules including members of phospholipases C (PLC), 14-3-3 proteins, and protein tyrosine phosphatases (PTP) are downregulated (green) in BRAFwt compared to BRAFmut PTCs.(JPG 6 MB)
About this article
Cite this article
Schulten, HJ., Alotibi, R., Al-Ahmadi, A. et al. Effect of BRAFmutational status on expression profiles in conventional papillary thyroid carcinomas. BMC Genomics 16 (Suppl 1), S6 (2015). https://doi.org/10.1186/1471-2164-16-S1-S6