Search

Transcriptome of papillary thyroid cancer (PTC) is well characterized and correlates with some prognostic and genotypic factors, but data addressing the interaction between PTC and tumor microenvironment (TME) are scarce. Therefore, in the present study, we aimed to assess the impact of TME on gene expression profile in PTC. We evaluated the gene expression profile in PTC and normal thyroid cells isolated by laser capture microdissection and in whole tissue slides corresponding to the entire tumor. We included 26 microdissected samples for gene expression analysis (HG-U133 PLUS 2.0, Affymetrix, currently Thermo Fisher Scientific USA): 15 PTC samples, 11 samples of normal thyrocytes, and 30 whole slides (15 PTC and 15 normal thyroid). Transcripts were divided into three groups: differentially expressed both in microdissected and whole slides, transcripts differently expressed in microdissected samples and not changed in whole slides, and transcripts differentially expressed in whole slides and not changed in microdissected samples. Eleven genes were selected for validation in an independent set of samples; among them, four genes differentiated only microdissected PTC and normal cells. Two genes (PTCSC and CTGF) were confirmed. One gene (FOS) was not confirmed by the validation, whereas EGR1 was also significant in whole slide analysis. The other seven genes (TFF3, FN1, MPPED2, MET, KCNJ2, TACSTD2, and GALE) showed differentiated expression in microdissected thyrocytes and in whole tumor slides. Most of identified genes were related to the tumor-microenvironment interaction and confirmed the crosstalk between TME and cancer cells.

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.

Introduction: Somatic RET mutations are detectable in two-thirds of sporadic cases of medullary thyroid cancer (MTC). Recent studies reported a high proportion of RAS somatic mutations in RET negative tumours, which may indicate RAS mutation as a possible alternative genetic event in sporadic MTC tumorigenesis. Thus, the aim of the study was to evaluate the frequency of somatic RAS mutations in sporadic medullary thyroid cancer in the Polish population and to relate the obtained data to the presence of somatic RET mutations. Material and methods: Somatic mutations (RET, RAS genes) were evaluated in 78 snap-frozen MTC samples (57 sporadic and 21 hereditary) by direct sequencing. Next, three randomly selected RET-negative MTC samples were analysed by the next generation sequencing. Results: RAS mutation was detected in 26.5% of 49 sporadic MTC tumours. None of all the analysed samples showed N-RAS mutation. When only RET-negative samples were considered, the prevalence of RAS mutation was 68.7%, compared to 6% observed in RET-positive samples. Most of these mutations were located in H-RAS codon 61 (72%). None of 21 hereditary MTC samples showed any RAS mutations. Conclusions: RAS mutations constitute a frequent molecular event in RET-negative sporadic medullary thyroid carcinoma in Polish patients. However, their role in MTC tumorigenesis remains unclear. (Endokrynol Pol 2015; 66 (2): 121–125)

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.

Epidermal differentiation complex (EDC) comprises a number of genes associated with human skin diseases including psoriasis, atopic dermatitis and hyperkeratosis. These genes have also been linked to numerous cancers, among them skin, gastric, colorectal, lung, ovarian and renal carcinomas. The involvement of EDC components encoding S100 proteins, small proline-rich proteins (SPRRs) and other genes in the tumorigenesis of head and neck squamous cell cancer (HNSCC) has been previously suggested. The aim of the study was to systematically analyze the expression of EDC components on the transcript level in HNSCC. Tissue specimens from 93 patients with HNC of oral cavity and 87 samples from adjacent or distant grossly normal oral mucosawere analyzed. 48 samples (24 tumor and 24 corresponding surrounding tissue) were hybridized to Affymetrix GeneChip Human 1.0 ST Arrays. For validation by quantitative real-time PCR (QPCR) the total RNA from all 180 samples collected in the study was analyzed with Real-Time PCR system and fluorescent amplicon specific-probes. Additional set of samples from 14 patients with laryngeal carcinoma previously obtained by HG-U133 Plus 2.0 microarray was also included in the analyses. The expression of analyzed EDC genes was heterogeneous. Two transcripts (S100A1 and S100A4) were significantly down-regulated in oral cancer when compared to normal mucosa (0.69 and 0.36-fold change, respectively), showing an opposite pattern of expression to the remaining S100 genes. Significant up-regulation in tumors was found for S100A11, S100A7, LCE3D, S100A3 and S100A2 genes. The increased expression of S100A7 was subsequently validated by QPCR, confirming significant differences. The remaining EDC genes, including all encoding SPRR molecules, did not show any differences between oral cancer and normal mucosa. The observed differences were also assessed in the independent set of laryngeal cancer samples, confirming the role of S100A3 and LCE3D transcripts in HNC. In HNC of oral cavity only one family of EDC genes (S100 proteins) showed significant cancer-related differences. A number of other transcripts which showed altered expression in HNC require further validation.

A single microRNA gene may give rise to several mature products that differ in length, called isomiRs. IsomiRs are known to be tissue specific and functionally relevant. The microRNA sequence heterogeneity of the thyroid gland has yet to be determined.The objective of the study was to provide a comprehensive view of the microRNA transcriptome in normal thyroid and papillary thyroid carcinoma (PTC).We used next-generation deep sequencing to analyze microRNA length heterogeneity and expression profiles of PTC tumors (n = 14), unaffected tissue adjacent to tumors (n = 14), and control, noncancerous thyroid tissue (n = 14). The results were validated with a microarray on an additional set of 9 PTC tumor/normal tissue pairs.Eighty-nine microRNAs were significantly deregulated in PTC compared with normal thyroid tissue (false discovery rate0.05, fold change 0.13-20.7). Top deregulated miRNAs included miR-146b-5p, miR-221-3p, miR-7-3p, miR-551b-3p, miR-486-3p, and miR-144-3p, confirming previous microarray profiling. The expression of miRNAs did not depend on the BRAF mutation status. Interestingly, 85% of the most abundant microRNAs consisted of isoforms that differed from the standard reference sequence deposited in miRBase. Moreover, the reference microRNAs were completely absent in 42.4% and 35.9% of the microRNAs expressed in normal thyroid and PTC tumors, respectively. Numerous isomiRs had altered seed sequences, which led to a different set of target genes. For highly deregulated miR-146b-5p, we detected 6 isoforms (tumor/normal fold change 14.4-28.7, false discovery rate0.002) that varied at their 5' ends with a 1-nt difference that created 2 alternative seeds. The target genes for those 2 seeds overlapped in only 13.1% of genes.Almost all microRNAs exhibit isoforms of variable length and potentially distinct function in thyroid tumorigenesis.

Medullary thyroid cancer (MTC) can be caused by germline mutations of the RET proto-oncogene or occurs as a sporadic form. It is well known that RET mutations affecting the cysteine-rich region of the protein (MEN2A-like mutations) are correlated with different phenotypes than those in the kinase domain (MEN2B-like mutations). Our aim was to analyse the whole-gene expression profile of MTC with regard to the type of RET gene mutation and the cancer genetic background (hereditary vs sporadic). We studied 86 MTC samples. We demonstrated that there were no distinct differences in the gene expression profiles of hereditary and sporadic MTCs. This suggests a homogeneous nature of MTC. We also noticed that the site of the RET gene mutation slightly influenced the gene expression profile of MTC. We found a significant association between the localization of RET mutations and the expression of three genes: NNAT (suggested to be a tumour suppressor gene), CDC14B (involved in cell cycle control) and NTRK3 (tyrosine receptor kinase that undergoes rearrangement in papillary thyroid cancer). This study suggests that these genes are significantly deregulated in tumours with MEN2A-like and MEN2B-like mutations; however, further investigations are necessary to demonstrate any clinical impact of these findings.

Endometrial stromal tumors include translocation-associated low- and high-grade endometrial stromal sarcomas (ESS) and highly malignant undifferentiated uterine sarcomas (UUS). UUS is considered a poorly defined group of aggressive tumors and is often seen as a diagnosis of exclusion after ESS and leiomyosarcoma (LMS) have been ruled out. We performed a comprehensive analysis of gene expression, copy number variation, point mutations, and immune cell infiltrates in the largest series to date of all major types of uterine sarcomas to shed light on the biology of UUS and to identify potential novel therapeutic targets. We show that UUS tumors have a distinct molecular profile from LMS and ESS. Gene expression and immunohistochemical analyses revealed the presence of high numbers of tumor-associated macrophages (TAMs) in UUS, which makes UUS patients suitable candidates for therapies targeting TAMs. Our results show a high genomic instability of UUS and downregulation of several TP53-mediated tumor suppressor genes, such as NDN, CDH11, and NDRG4. Moreover, we demonstrate that UUS carry somatic mutations in several oncogenes and tumor suppressor genes implicated in RAS/PI3K/AKT/mTOR, ERBB3, and Hedgehog signaling.

Medullary thyroid cancer (MTC) can be caused by germline mutations of the RET proto-oncogene or occurs as a sporadic form. It is well known that RET mutations affecting the cysteine-rich region of the protein (MEN2A-like mutations) are correlated with different phenotypes than those in the kinase domain (MEN2B-like mutations). Our aim was to analyse the whole-gene expression profile of MTC with regard to the type of RET gene mutation and the cancer genetic background (hereditary vs sporadic). We studied 86 MTC samples. We demonstrated that there were no distinct differences in the gene expression profiles of hereditary and sporadic MTCs. This suggests a homogeneous nature of MTC. We also noticed that the site of the RET gene mutation slightly influenced the gene expression profile of MTC. We found a significant association between the localization of RET mutations and the expression of three genes: NNAT (suggested to be a tumour suppressor gene), CDC14B (involved in cell cycle control) and NTRK3 (tyrosine receptor kinase that undergoes rearrangement in papillary thyroid cancer). This study suggests that these genes are significantly deregulated in tumours with MEN2A-like and MEN2B-like mutations; however, further investigations are necessary to demonstrate any clinical impact of these findings.

a b s t r a c t Hepatocellular carcinoma (HCC) represents the major histological subtype of liver cancer. Tumorigenic changes in hepatic cells potentially result from aberrant expression of microRNAs (miRNAs). Individual microRNA gene may give rise to miRNAs of different length, named isomiRNAs that proved to be func- tionally relevant. Since microRNA length heterogeneity in hepatic tissue has not been described before, we employed next-generation sequencing to comprehensively analyze microRNA transcriptome in HCC tumors (n = 24) and unaffected tissue adjacent to tumors (n = 24), including samples with (n = 15) and without cirrhosis (n = 9). We detected 374 microRNAs expressed in liver, including miR-122-5p that constituted over 39% of the hepatic miRnome. Among the liver expressed miRs, the levels of 64 significantly differed between tumor and control samples (FDR 2). Top deregulated miRNAs included miR-1269a (T/N = 22.95), miR-3144-3p (T/N = 5.24), miR-183-5p (T/N = 4.63), miR-10b-5p (T/N = 3.87), miR-490-3p (T/N = 0.13), miR-199a-5p (T/N = 0.17), miR-199a-3p/miR-199b-3p (T/N = 0.19), miR-214-5p (T/N = 0.20) and miR-214-3p (T/N = 0.21). Almost all miRNA genes produced several mature molecules differing in length (isomiRNAs). The reference sequence was not the most prevalent in 38.6% and completely absent in 10.5% of isomiRNAs. Over 26.1% of miRNAs produced isoforms carrying ≥ 2 alternative seed regions, of which 35.5% constituted novel, previously unknown seeds. This fact sheds new light on the percentage of the human genome regulated by microRNAs and their variants. Among the most deregulated miR- NAs, miR-199a-3p/miR-199b-3p (T/N fold change = 0.18, FDR = 0.005) was expressed in 9 isoforms with 3 different seeds, concertedly leading to upregulation of TGF-beta signaling pathway (OR = 1.99; p = 0.004). In conclusion, the study reveals the comprehensive miRNome of hepatic tissue and provides new tools for investigation of microRNA-dependent pathways in cirrhotic liver and hepatocellular carcinoma. This article is part of a Directed Issue entitled: Rare Cancers.

Sarcomas of the uterus are derived from uterine smooth muscle (leiomyosarcomas) or from endometrial stromal cells. The latter group includes three subtypes: the low-grade (LG) and high-grade (HG) endometrial stromal sarcomas (ESS), each defined by characteristic chromosomal translocations, and undifferentiated uterine sarcoma (UUS), the most aggressive form. The majority of UUS patients present with high-stage disease and even patients with stage I tumors often die within 2 years from diagnosis. Adjuvant radio- and chemotherapy do not improve the clinical outcome. No prior comprehensive molecular analysis of all three types of endometrial stromal tumors has been reported. We performed genomic analysis (whole exome sequencing and aCGH), gene expression profiling (RNA-seq and microarrays), immunohistochemistry (IHC) and FISH on 31 endometrial stromal tumors, including 17 UUS and 14 LG and HG ESS. All ESS cases carried characteristic JAZF1, PHF1, MBTD1 or YWHAE rearrangements. Selected findings were validated by IHC, qRT-PCR and Sanger sequencing. Genomic analysis revealed that UUS are characterized by complex chromosomal aberrations. The most frequent somatic mutations (in TP53, KRAS, FBXW7, PIK3CA and ERBB3) and copy number changes (e.g. gains of 19q, 8q11.1-q24.3 and 3q26.2-q29) in UUS resemble those seen in uterine carcinosarcomas and carcinomas. UUS over-expressed numerous genes encoding M2 macrophage-specific markers, including CD163, CCL18, mannose receptor, stabilin-1, and macrophage galactose-type C-type lectin 2. Immunohistochemistry for CD163 and CD68 confirmed high tumor-associated macrophages (TAM) counts in UUS tissues compared to the ESS. In UUS, we also observed significantly higher mRNA expression for genes implicated in angiogenesis (CD34, MMP9), immunosuppression and tumor invasion (e.g. CCR2, IL10RA, IRAK3, CCL13, CXCL9, CXCL10). TAMs are associated with progression, resistance to cytotoxic therapy and metastatic spread in most human tumor types. In animal models, TAMs were shown to induce increased angiogenesis. CSF1 is a major regulator of macrophage function and is implicated in these tumor-promoting effects. Clinical trials have shown an effect of inhibitors of the CSF1 pathway in tenosynovial giant cell tumors and our data indicate that these inhibitors may be considered for the treatment of UUS. Our study also showed significantly elevated expression in UUS of two additional therapeutic targets involving macrophage pathways: CCR2 (C-C chemokine receptor type 2) that can be targeted by monoclonal antibodies, and BTK (Bruton's tyrosine kinase) that can be targeted by ibrutinib. In conclusion, our findings demonstrate abundant presence of TAMs and over-expression of TAM-associated markers in UUS compared to the other tumors derived from endometrial stroma. Moreover, our results indicate novel potential therapeutic targets for UUS patients management. Citation Format: Joanna Przybyl, Magdalena Kowalewska, Anna Quattrone, Barbara Dewaele, Vanessa Vanspauwen, Sushama Varma, Robert T. Sweeney, Michal Swierniak, Elwira Bakula-Zalewska, Janusz A. Siedlecki, Mariusz Bidzinski, Jan Cools, Matt van de Rijn, Maria Debiec-Rychter. Tumor associated macrophages in undifferentiated uterine sarcoma: association with angiogenesis and therapeutic implications. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 3182.

Introduction Endometrial stromal tumors (EST) are the second most prevalent uterine sarcomas. EST are divided into three histosubtypes: benign endometrial stromal nodules (ESN), low-grade endometrial stromal sarcomas (ESS) and the most malignant undifferentiated stromal sarcomas (UES). Moreover, due to clinicopathological features intermediate between classical low-grade ESS and UES, some tumors have been described in the literature as “high-grade ESS”. The histological distinction between these three subtypes has important prognostic implications. The molecular biology of these tumors remains poorly understood because of their low incidence. Patients and Methods In the present study, we performed gene expression profiling of 22 EST using microarray, RNA-Seq and qRT-PCR techniques. Our cohort included 10 classical low-grade ESS (8 with ESS-associated gene fusions), 8 UES and 4 tumors which were classified as high-grade ESS based on the presence of YWHAE-FAM22A/B fusion and histological features. Results Using gene expression microrarrays, we have identified 81 entities differentially expressed in analyzed samples (fold change > 2, p < 0.005). Hierarchical clustering analysis revealed that high-grade ESS formed a separate sub-cluster with a distinct molecular profile as compared with the low-grade ESS and UES cases. Low-grade ESS cases formed 2 distinct sub-clusters, each containing 5 cases, one of which presented gene expression profile partially corresponding to high-grade ESS cases. GO analysis of genes strongly up-regulated in UES specimens showed over-representation of immune response (GO:0006955) and immune system process (GO:0002376) pathways (MARCO, MYO1F, FCN1, GPR183, HLA-DPB1, FCGR2B, NCF1, CCL18, CD28, PTAFR, CCL13, STK4, CD14, NCF4 and CTLA4 genes) (p < 0.05). UES cases were also characterized by significant down-regulation of WT1 tumor suppressor gene as compared to the low- and high-grade ESS. Moreover, both high-grade ESS and UES cases were characterized by significantly lower expression level of PEG3 gene, which has been reported as a tumor suppressor gene. All low-grade ESS cases presented significant over-expression of estrogen receptor ESR1 gene. Microarray results were validated by RNA-Seq and qRT-PCR experiments. Conclusions This study provides gene expression profiles that distinguish between EST subtypes and sheds light on the biology of these tumors. Our results point to immune system involvement in the pathogenesis of UES and support the notion that the current EST classification should discriminate high-grade ESS as a distinct subtype. Citation Format: Joanna Przybyl, Magdalena Kowalewska, Anna Quattrone, Barbara Dewaele, Vanessa Vanspauwen, Julio Finalet-Ferreiro, Michal Swierniak, Elwira Bakula-Zalewska, Janusz A. Siedlecki, Mariusz Bidzinski, Jan Cools, Maria Debiec-Rychter. Gene expression profiling distinguishes between endometrial stromal tumors subtypes. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 4698. doi:10.1158/1538-7445.AM2014-4698

Background Differential diagnosis between malignant follicular thyroid cancer (FTC) and benign follicular thyroid adenoma (FTA) is a great challenge for even an experienced pathologist and requires special effort. Molecular markers may potentially support a differential diagnosis between FTC and FTA in postoperative specimens. The purpose of this study was to derive molecular support for differential post-operative diagnosis, in the form of a simple multigene mRNA-based classifier that would differentiate between FTC and FTA tissue samples. Methods A molecular classifier was created based on a combined analysis of two microarray datasets (using 66 thyroid samples). The performance of the classifier was assessed using an independent dataset comprising 71 formalin-fixed paraffin-embedded (FFPE) samples (31 FTC and 40 FTA), which were analysed by quantitative real-time PCR (qPCR). In addition, three other microarray datasets (62 samples) were used to confirm the utility of the classifier. Results Five of 8 genes selected from training datasets (ELMO1, EMCN, ITIH5, KCNAB1, SLCO2A1) were amplified by qPCR in FFPE material from an independent sample set. Three other genes did not amplify in FFPE material, probably due to low abundance. All 5 analysed genes were downregulated in FTC compared to FTA. The sensitivity and specificity of the 5-gene classifier tested on the FFPE dataset were 71% and 72%, respectively. Conclusions The proposed approach could support histopathological examination: 5-gene classifier may aid in molecular discrimination between FTC and FTA in FFPE material.