Michał Świerniak

Analysis options for high-throughput sequencing in miRNA expression profiling

BackgroundRecently high-throughput sequencing (HTS) using next generation sequencing techniques became useful in digital gene expression profiling.Our study introduces analysis options for HTS data based on mapping to miRBase or counting and grouping of identical sequence reads. Those approaches allow a hypothesis free detection of miRNA differential expression.MethodsWe compare our results to microarray and qPCR data from one set of RNA samples. We use Illumina platforms for microarray analysis and miRNA sequencing of 20 samples from benign follicular thyroid adenoma and malignant follicular thyroid carcinoma. Furthermore, we use three strategies for HTS data analysis to evaluate miRNA biomarkers for malignant versus benign follicular thyroid tumors.ResultsHigh correlation of qPCR and HTS data was observed for the proposed analysis methods. However, qPCR is limited in the differential detection of miRNA isoforms. Moreover, we illustrate a much broader dynamic range of HTS compared to microarrays for small RNA studies. Finally, our data confirm hsa-miR-197-3p, hsa-miR-221-3p, hsa-miR-222-3p and both hsa-miR-144-3p and hsa-miR-144-5p as potential follicular thyroid cancer biomarkers.ConclusionsCompared to microarrays HTS provides a global profile of miRNA expression with higher specificity and in more detail. Summarizing of HTS reads as isoform groups (analysis pipeline B) or according to functional criteria (seed analysis pipeline C), which better correlates to results of qPCR are promising new options for HTS analysis. Finally, data opens future miRNA research perspectives for HTS and indicates that qPCR might be limited in validating HTS data in detail.

Variants in the ATM-CHEK2-BRCA1 Axis Determine Genetic Predisposition and Clinical Presentation of Papillary Thyroid Carcinoma

The risk of developing papillary thyroid carcinoma (PTC), the most frequent form of thyroid malignancy, is elevated up to 8.6‐fold in first‐degree relatives of PTC patients. The familial risk could be explained by high‐penetrance mutations in yet unidentified genes, or polygenic action of low‐penetrance alleles. Since the DNA‐damaging exposure to ionizing radiation is a known risk factor for thyroid cancer, polymorphisms in DNA repair genes are likely to affect this risk. In a search for low‐penetrance susceptibility alleles we employed Sequenom technology to genotype deleterious polymorphisms in ATM, CHEK2, and BRCA1 in 1,781 PTC patients and 2,081 healthy controls. As a result of the study, we identified CHEK2 rs17879961 (OR = 2.2, P = 2.37e‐10) and BRCA1 rs16941 (odds ratio [OR] = 1.16, P = 0.005) as risk alleles for PTC. The ATM rs1801516 variant modifies the risk associated with the BRCA1 variant by 0.78 (P = 0.02). Both the ATM and BRCA1 variants modify the impact of male gender on clinical variables: T status (P = 0.007), N status (P = 0.05), and stage (P = 0.035). Our findings implicate an important role of variants in the ATM‐ CHEK2‐ BRCA1 axis in modification of the genetic predisposition to PTC and its clinical manifestations.

Association between GWAS-Derived rs966423 Genetic Variant and Overall Mortality in Patients with Differentiated Thyroid Cancer.

Purpose: Five germline genetic variants (rs116909374, rs965513, rs944289, rs966423, and rs2439302) have been associated in genome-wide association studies (GWAS) with increased risk of differentiated thyroid cancer (DTC), but their role in mortality of patients has not been established. Also, no preoperative marker of the clinical outcome of thyroid cancer had yet been identified. The aim of the study was to investigate the relationship between the variants and overall mortality in patients with DTC. Experimental Design: Retrospective study of 1,836 patients (1,643 women, 193 men) with median age at diagnosis of 49 years and overall median follow-up time of 8.7 years after initial treatment at a single comprehensive cancer center between 1990 and 2013. Results: Among 5 variants, rs966423 was associated with increased mortality, which was 6.4% (33 of 518) versus 3.7% (47 of 1,259) in TT carriers versus CC/CT carriers (P = 0.017). The HR of TT versus TC/CC carriers was 1.6 [95% confidence interval (CI), 1.02–2.49; P = 0.038] after adjustment for age at diagnosis and sex. Importantly, the association of rs966423 with mortality remained valid when clinicopathologic risk factors were included in the model (HR, 1.89; 95% CI, 1.14–3.13; P = 0.014). Higher rs966423–associated patient mortality of TT versus CC/CT carriers was also observed in interaction with angioinvasion (adjusted HR, 3.48; 95% CI, 1.67–7.22; P < 0.001), lymph node metastasis (adjusted HR, 3.47; 95% CI, 1.16–10.4; P = 0.018), extrathyroidal invasion (adjusted HR, 2.07; 95% CI, 1.15–3.73; P = 0.013). Conclusions: The presence of the rs966423-TT genotype was associated with a significant increase in overall mortality of patients with DTC. Contrary to BRAF mutation and other somatic changes, the status of germline rs966423 is known before the treatment and might be used in the management of mortality risk by means of modification of therapy. Clin Cancer Res; 22(5); 1111–9. ©2015 AACR.

Next-generation sequencing reveals microRNA markers of adrenocortical tumors malignancy

Background Adrenocortical carcinoma is a rare finding among common adrenocortical tumors, but it is highly aggressive and requires early detection and treatment. Still, the differential diagnosis between benign and malignant lesions is difficult even for experienced pathologists and there is a significant need for novel diagnostic methods. In this study we aimed to reveal a complete set of microRNAs expressed in the adrenal gland and to identify easily detectable, stable and objective biomarkers of adrenocortical malignancy. Methods We employed next-generation sequencing to analyze microRNA profiles in a unique set of 51 samples, assigned to either a learning dataset including 7 adrenocortical carcinomas (ACCs), 8 adrenocortical adenomas (AAs) and 8 control samples (NAs), or a validation dataset including 8 ACCs, 10 AAs and 10 NAs. The results were validated in real-time Q-PCR. Results We detected 411 miRNAs expressed in 1763 length isoforms in the examined samples. Fifteen miRNAs differentiate between malignant (ACC) and non-malignant (AA + NA) tissue in the test set of independent samples. Expression levels of 6 microRNAs, miR-503-5p, miR-483-3p, miR-450a-5p, miR-210, miR-483-5p, miR-421, predict sample status (malignancy/non-malignancy) with at least 95% accuracy in both datasets. The best single-gene malignancy marker, miR-483-3p, has been validated by real-time RT PCR. Conclusions As a result of the study we propose clinically valid and easily detectable biomarkers of adrenocortical malignancy that may significantly facilitate morphological examination. Since microRNAs can be detected in blood, the study brings tools for development of non-invasive diagnostics of adrenocortical carcinomas.

Unsupervised analysis of follicular thyroid tumours transcriptome by oligonucleotide microarray gene expression profiling.

INTRODUCTION Mechanisms driving the invasiveness of follicular thyroid cancer (FTC) are not fully understood. In our study, we undertook an unsupervised analysis of the set of follicular thyroid tumours (adenomas (FTA) and carcinomas) to verify whether the malignant phenotype influences major sources of variability in our dataset. MATERIAL AND METHODS The core set of samples consisted of 52 tumours (27 FTC, 25 FTA). Total RNA was analysed by oligonucleotide microarray (HG-U133 Plus 2.0). Principal Component Analysis (PCA) was applied as a main method of unsupervised analysis. RESULTS An analysis of biological character of genes correlated to the first six PCs was performed. When genes correlated to the first PC were used to cluster FTC and FTA, they appeared in two branches; one, relatively enriched in adenomas, with homogenous expression of subset of genes, and the other containing mainly carcinomas, with down-regulation of these genes and heterogeneous up-regulation in a smaller cluster of transcripts. Genes highly up-regulated in adenomas included some thyroid-specific transcripts. The second cluster of genes, up-regulated in carcinomas, contained mainly immunity-related transcripts. Immune response genes were found in the first, third and sixth principal components, improving the discrimination between carcinomas and adenomas. CONCLUSIONS Our unsupervised analysis indicates that invasiveness of follicular tumours might be considered as the major source of variability in transcriptome analysis. However, the distance between both groups is small and the clusters are overlapping, thus, unsupervised analysis is not sufficient to properly classify them.

Functional analysis of a novel, thyroglobulin-embedded microRNA gene deregulated in papillary thyroid carcinoma.

MicroRNAs, non-coding regulators of gene expression, are known culprits of thyroid cancer. Using next-generation sequencing, we identified a novel microRNA gene, encoded within an important thyroid regulator – thyroglobulin, and analyzed its functionality in the thyroid gland. In vitro and in silico analyses proved that the novel miR-TG is processed from the precursor, and co-expressed with thyroglobulin. Both genes are specific for thyroid tissue and downregulated in papillary thyroid carcinoma by 44% (p = 0.04) and 48% (p = 0.001), respectively. Putative target genes for miR-TG were identified using in silico tools, which pinpointed MAP4K4, an oncogene upregulated in thyroid cancer. Analysis of transcriptome by RNA-seq revealed that overexpression of miR-TG in PTC-derived cell line led to downregulation of several genes, including MAP4K4 (fold change 0,82; p = 0.036). The finding was confirmed by SQ-PCR (fold change 071; p = 0.004). Direct interaction between miR-TG and MAP4K4 was confirmed in the luciferase assay (p = 0.0006). Functional studies showed increase proliferation in K1 cell line transfected with miR-TG. We propose that in normal thyroid miR-TG plays a fine-tuning effect on the maintenance of MAPK pathway, inhibiting the expression of miR’s target MAP4K4. This regulation is disturbed in cancer due to downregulation of the novel, thyroglobulin-embedded microRNA, characterized in this study.

Age at diagnosis and gender modify the risk of 9q22 and 14q13 polymorphisms for papillary thyroid carcinoma

INTRODUCTION Papillary thyroid cancer (PTC) shows familial occurrence, and some susceptibility single nucleotide polymorphisms (SNPs) have been identified in FOXE1 and near the NKX2-1 locus. The aim of our study was to analyse the association of PTC risk with SNPs in FOXE1 (rs965513, rs1867277, rs1443434) and near the NKX2-1 locus (rs944289) in a Polish population, and, in the second step, the interac-tion between SNPs and patient-related factors (age at diagnosis and gender). MATERIAL AND METHODS A total of 2243 DNA samples from PTC patients and 1160 controls were included in the study. The SNP analysis was performed with the allelic discrimination technique. RESULTS There were significant associations of all SNPs with PTC (rs965513 odds ratio [OR] = 1.72, p = 8 × 10-7; rs1867277 OR = 1.59, p = 1 × 10-6; rs1443434 OR = 1.53, p = 1 × 10-5; rs944289 OR = 1.52, p = 4 × 10-5). Logistic regression analysis revealed an increased PTC risk in the interaction of rs944289 with age at diagnosis (OR = 1.01 per year, p = 6 × 10-4) and a decreased PTC risk in the interaction of male gender with the GGT FOXE1 protective haplotype (OR = 0.69, p = 0.01). CONCLUSIONS the association between PTC and all analysed SNPs was confirmed. It was also shown that patient-related factors modify the predisposition to PTC by increasing the risk for rs944289 per year of age, and by enhancing the protective effect of the FOXE1 GGT haplotype in men.

The Role and Use of MicroRNA in Cancer Diagnosis

MicroRNAs (miRNAs) are small RNA molecules that posttranscriptionally regulate gene expression. The aim of this paper is to present a survey on published application of these molecules in cancer diagnosis including our own results in this area. As an illustrative example of the use of miRNA in cancer diagnosis we present some main points from our study related to thyroid carcinoma.

The rs2910164 Genetic Variant of miR-146a-3p Is Associated with Increased Overall Mortality in Patients with Follicular Variant Papillary Thyroid Carcinoma

Aberrant expression of the sodium-iodide symporter (NIS) and the resistance to post-operative radioactive iodide treatment is a crucial cause of higher mortality of some thyroid cancer patients. In this study, we analyzed the impact of miR-146a on the expression and function of NIS and on the overall survival of thyroid cancer patients. The study included 2441 patients (2163 women; 278 men); including 359 cases with follicular variant of papillary thyroid carcinoma (fvPTC). miR:NIS interactions were analyzed in cell lines using in vivo binding and inhibition assays and radioactive iodine uptake assays. Tumor/blood DNA was used for rs2910164 genotyping. Overall survival was assessed retrospectively. In the results, we showed that miR-146a-3p directly binds to and inhibits NIS. Inhibition of miR-146a-3p restores the expression and function of NIS, increasing radioactive iodine uptake. Rs2910164 functional variant within miR-146a-3p is associated with increased overall mortality among fvPTC female patients. The deaths per 1000 person-years were 29.7 in CC carriers vs. 5.08 in GG/GC-carriers (HR = 6.21, p = 0.006). Higher mortality of CC vs. GG/GC carriers was also observed in patients with lower clinical stage (HR = 22.72, p < 0.001), smaller tumor size (pT1/pT2) (HR = 25.05, p < 0.001), lack of extrathyroidal invasion (HR = 9.03, p = 0.02), lack of nodular invasion (HR = 7.84, p = 0.002), lack of metastases (HR = 6.5, p = 0.005) and older (age at diagnosis >50 years) (HR = 7.8, p = 0.002). MiR-146a-3p underwent somatic mutations in 16.1% of analyzed specimens, mainly towards the deleterious C allele. In this report we propose a novel molecular marker of the clinical outcome of fvPTC patients. Rs2910164 increases the overall mortality with inhibition of NIS and disruption of radioiodine uptake as a possible mechanism.

NGS Reveals Molecular Pathways Affected by Obesity and Weight Loss-Related Changes in miRNA Levels in Adipose Tissue

Both obesity and weight loss may cause molecular changes in adipose tissue. This study aimed to characterize changes in adipose tissue miRNome in order to identify molecular pathways affected by obesity and weight changes. Next generation sequencing (NGS) was applied to identify microRNAs (miRNAs) differentially expressed in 47 samples of visceral (VAT) and subcutaneous (SAT) adipose tissues from normal-weight (N), obese (O) and obese after surgery-induced weight loss (PO) individuals. Subsequently miRNA expression was validated by real-time PCR in 197 adipose tissues and bioinformatics analysis performed to identify molecular pathways affected by obesity-related changes in miRNA expression. NGS identified 344 miRNAs expressed in adipose tissues with ≥5 reads per million. Using >2 and <−2 fold change as cut-offs we showed that the expression of 54 miRNAs differed significantly between VAT-O and SAT-O. Equally, between SAT-O and SAT-N, the expression of 20 miRNAs differed significantly, between SAT-PO and SAT-N the expression of 79 miRNAs differed significantly, and between SAT-PO and SAT-O, the expression of 61 miRNAs differed significantly. Ontological analyses disclosed several molecular pathways regulated by these miRNAs in adipose tissue. NGS-based miRNome analysis characterized changes of the miRNA profile of adipose tissue, which are associated with changes of weight possibly responsible for a differential regulation of molecular pathways in adipose tissue when the individual is obese and after the individual has lost weight.