Michal Swierniak

Polymorphic mature microRNAs from passenger strand of pre-miR-146a contribute to thyroid cancer

Prior work has shown that heterozygosity G/C of single nucleotide polymorphism (SNP rs2910164) within the precursor of microRNA-146a predisposes to PTC (odds ratio = 1.62, P = 0.000007) although the mechanism was unclear. Here, we show that GC heterozygotes differ from both GG and CC homozygotes by producing 3 mature microRNAs: 1 from the leading strand (miR-146a), and 2 from the passenger strand (miR-146a*G and miR-146a*C), each with its distinct set of target genes. TaqMan analysis of paired tumor/normal samples revealed 1.5- to 2.6-fold overexpression of polymorphic miR-146a* in 7 of 8 tumors compared with the unaffected part of the same gland. The microarray data showed that widely different transcriptomes occurred in the tumors and in unaffected parts of the thyroid from GC and GG patients. The modulated genes are mainly involved in regulation of apoptosis leading to exaggerated DNA-damage response in heterozygotes potentially explaining the predisposition to cancer. We propose that contrary to previously held views transcripts from the passenger strand of miRs can profoundly affect the downstream effects. Heterozygosity for polymorphisms within the premiR sequence can cause epistasis through the production of additional mature miRs. We propose that mature miRs from the passenger strand may regulate many genetic processes.

Next generation sequencing reveals microRNA isoforms in liver cirrhosis and hepatocellular carcinoma

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 functionally 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<0.05, fold change>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 miRNAs, 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.

Epidermal differentiation complex (locus 1q21) gene expression in head and neck cancer and normal mucosa.

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 all180 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.

Molecular differential diagnosis of follicular thyroid carcinoma and adenoma based on gene expression profiling by using formalin-fixed paraffin-embedded tissues

BackgroundDifferential 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.MethodsA 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.ResultsFive 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.ConclusionsThe proposed approach could support histopathological examination: 5-gene classifier may aid in molecular discrimination between FTC and FTA in FFPE material.

Differences in the transcriptome of medullary thyroid cancer regarding the status and type of RET gene mutations

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.

Gene signature of the post-Chernobyl papillary thyroid cancer

PurposeFollowing the nuclear accidents in Chernobyl and later in Fukushima, the nuclear community has been faced with important issues concerning how to search for and diagnose biological consequences of low-dose internal radiation contamination. Although after the Chernobyl accident an increase in childhood papillary thyroid cancer (PTC) was observed, it is still not clear whether the molecular biology of PTCs associated with low-dose radiation exposure differs from that of sporadic PTC.MethodsWe investigated tissue samples from 65 children/young adults with PTC using DNA microarray (Affymetrix, Human Genome U133 2.0 Plus) with the aim of identifying molecular differences between radiation-induced (exposed to Chernobyl radiation, ECR) and sporadic PTC. All participants were resident in the same region so that confounding factors related to genetics or environment were minimized.ResultsThere were small but significant differences in the gene expression profiles between ECR and non-ECR PTC (global test, p < 0.01), with 300 differently expressed probe sets (p < 0.001) corresponding to 239 genes. Multifactorial analysis of variance showed that besides radiation exposure history, the BRAF mutation exhibited independent effects on the PTC expression profile; the histological subset and patient age at diagnosis had negligible effects. Ten genes (PPME1, HDAC11, SOCS7, CIC, THRA, ERBB2, PPP1R9A, HDGF, RAD51AP1, and CDK1) from the 19 investigated with quantitative RT-PCR were confirmed as being associated with radiation exposure in an independent, validation set of samples.ConclusionSignificant, but subtle, differences in gene expression in the post-Chernobyl PTC are associated with previous low-dose radiation exposure.

BRAFV600E-Associated Gene Expression Profile: Early Changes in the Transcriptome, Based on a Transgenic Mouse Model of Papillary Thyroid Carcinoma

Background The molecular mechanisms driving the papillary thyroid carcinoma (PTC) are still poorly understood. The most frequent genetic alteration in PTC is the BRAFV600E mutation–its impact may extend even beyond PTC genomic profile and influence the tumor characteristics and even clinical behavior. Methods In order to identify BRAF-dependent signature of early carcinogenesis in PTC, a transgenic mouse model with BRAFV600E-induced PTC was developed. Mice thyroid samples were used in microarray analysis and the data were referred to a human thyroid dataset. Results Most of BRAF(+) mice developed malignant lesions. Nevertheless, 16% of BRAF(+) mice displayed only benign hyperplastic lesions or apparently asymptomatic thyroids. After comparison of non-malignant BRAF(+) thyroids to BRAF(−) ones, we selected 862 significantly deregulated genes. When the mouse BRAF-dependent signature was transposed to the human HG-U133A microarray, we identified 532 genes, potentially indicating the BRAF signature (representing early changes, not related to developed malignant tumor). Comparing BRAF(+) PTCs to healthy human thyroids, PTCs without BRAF and RET alterations and RET(+), RAS(+) PTCs, 18 of these 532 genes displayed significantly deregulated expression in all subgroups. All 18 genes, among them 7 novel and previously not reported, were validated as BRAFV600E-specific in the dataset of independent PTC samples, made available by The Cancer Genome Atlas Project. Conclusion The study identified 7 BRAF-induced genes that are specific for BRAF V600E-driven PTC and not previously reported as related to BRAF mutation or thyroid carcinoma: MMD, ITPR3, AACS, LAD1, PVRL3, ALDH3B1, and RASA1. The full signature of BRAF-related 532 genes may encompass other BRAF-related important transcripts and require further study.

Somatic mutation profiling of follicular thyroid cancer by next generation sequencing

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.

The prevalence of somatic RAS mutations in medullary thyroid cancer — a Polish population study

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.

Gene Expression (mRNA) Markers for Differentiating between Malignant and Benign Follicular Thyroid Tumours

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.