DNA copy number variations characterize benign and malignant thyroid tumors

Coexisting Molecular Alterations Increase the Risk of Malignancy in Thyroid Nodules with Copy Number Alterations

Molecular mutations and alterations play a role in thyroid tumorigenesis. Different alterations are associated with different clinical and pathological characteristics. Copy number alterations (CNAs) are known to be present in some thyroid tumors; however, their idiosyncratic clinicopathological implications are not yet well elucidated. A retrospective chart review was performed to identify patients with CNAs on pre-operative molecular testing results who subsequently underwent surgical treatment between January 2016 and April 2022 at McGill University teaching hospitals. Of the 316 patients with thyroid nodules who opted for molecular testing with ThyroSeqV3 followed by surgery, 67 (21.2%) nodules were positive for CNAs, including 23 Bethesda III, 31 Bethesda IV, 12 Bethesda V and 1 Bethesda VI nodules. On surgical pathology, 29.9% were benign and 70.1% were malignant or non-invasive follicular thyroid neoplasm with papillary-like nuclear features (NIFTP). Among those that were malignant/NIFTP, 17.02% were considered to be aggressive cancers. The presence of other molecular alterations was found to be an independent predictor of malignancy in multivariate analysis (OR = 5.087, 95% C.I. = 1.12-23.04, p = 0.035). No unique factor was correlated with aggressiveness; however, CNA-positive thyroid nodules that were associated with high-risk mutations such as BRAF V600E, TP53, NTRK1/3 fusion, or PTEN mutation with high allele frequency (AF) ended up being aggressive cancers. Most of the CNA-positive thyroid nodules resulted in follicular patterned tumors in 41 (65.2%) cases and oncocytic tumors in 20 (29.9%) cases. This study demonstrates that 70.1% of surgically resected thyroid nodules with CNAs were malignant/NIFTP. Most CNA-positive thyroid nodules were either oncocytic patterned tumors or follicular patterned tumors. Furthermore, CNA-positive thyroid nodules were more likely to be malignant if they were associated with other molecular alterations or mutations.

Copy number variations: A novel molecular marker for papillary thyroid cancer

Background

We aimed to screen tumor-associated functional genes on a large scale through copy number variation (CNV) analysis of papillary thyroid cancer (PTC).

Methods

We analyzed 74 tissue samples from 41 patients with thyroid nodules. The samples were subjected to whole-genome resequencing and then analyzed by the ‘WISECONDOR’ method. Potential chromosome CNV regions were identified between the different sample groups.

Results

Of the 74 samples from 41 patients, 28 were PTC tissue samples, 29 were para-carcinoma tissue samples, 13 were benign tumor tissue samples and 4 were para-benign tumor tissue samples. According to our findings, PTC can be identified by CNVs at the corresponding positions on chromosomes 5, 7, 8, 10, and 17. For carcinoma tissue, the sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), accuracy and area under the curve (AUC) of the test method were 100%, 66.7%, 87.5%, 100.0%, 90.0% and 0.83 (95% confidence interval [CI], 0.67–1.00) and for para-carcinoma tissue, these values were 96.6%, 75.0%, 96.6%, 75.0%, 93.9% and 0.86 (95% CI, 0.60–1.00).

Conclusion

CNV analysis assays involving high-volume sequencing analysis can increase the identification of PTC, potentially avoiding errors caused by position deflection in sampling.

1.

Thyroid nodules can be identified by CNVs at the corresponding positions on chromosomes 5, 7, 8, 10, and 17.

2.

The identification rate of PTC can be greatly increased through high-volume CNV sequencing analysis.

Identification of Three Prognosis-Related Differentially Expressed lncRNAs Driven by Copy Number Variation in Thyroid Cancer

Thyroid cancer as the malignant tumor with the highest incidence in the endocrine system also shows a fast growth and development. In this work, we developed a new method to identify copy number variation- (CNV-) driven differentially expressed lncRNAs in thyroid cancer for predicting cancer prognosis. The data of RNA sequencing, CNV, methylation, mutation, and clinical details of thyroid cancer were obtained from the Cancer Genome Atlas database (TCGA). Molecular subtypes were clustered by iClusterPlus. Weighted gene co-expression network analysis (WGCNA) was employed to show co-expression modules. DEseq2 was conducted to identify protein coding genes (PCGs) and differentially expressed lncRNAs. CNV was detected using GISTIC 2.0. Three molecular subtypes were identified, and 68 differentially expressed lncRNAs (DElncRNAs) related to cancer were found among different molecular subtypes. CNV of FOXD2-AS1, FAM181A-AS1, and RNF157-AS1 was associated with overall survival and was involved in cancer-related pathways. These three DElncRNAs discovered based on CNV could serve as prognostic biomarkers to predict prognosis for thyroid cancer and new targets to explore molecular drugs.

Copy number variations identified in thyroid FNA specimens are associated with Hürthle cell cytomorphology

BACKGROUND

The fine-needle aspiration (FNA) diagnosis of thyroid Hürthle cell neoplasms (HCNs) remains challenging. This study explored a possible association of copy number variations (CNVs) with Hürthle cell lesions of the thyroid.

METHODS

Thyroid FNA cases that were diagnosed as follicular lesion of undetermined significance (FLUS) or follicular neoplasm (FN)/HCN for which the ThyroSeq version 3 genomic classifier test was performed were retrieved.

RESULTS

A total of 324 thyroid FNA cases (228 FLUS cases, 46 HCN cases, and 50 FN cases) were included in the study. FLUS cases were further classified as Hürthle cell type (follicular lesion of undetermined significance-Hürthle cell type [FLUS-HCT]; 20 cases) or non-Hürthle cell type (follicular lesion of undetermined significance-non-Hürthle cell type [FLUS-NHCT]; 208 cases). HCN and FLUS-HCT cases showed a higher prevalence of CNVs (23 of 66 [35%]) in comparison with those classified as FN or FLUS-NHCT (14 of 258 [5%]; P < .001). A total of 105 patients had histopathologic follow-up. Cases with CNVs were more likely to be neoplastic (18 of 26 [69%]) and associated with Hürthle cell changes (14 of 26 [54%]) in comparison with cases without any molecular alterations (neoplastic, 8 of 24 [33%]; Hürthle cell changes, 2 of 24 [8%]; P < .05). In HCN/FLUS-HCT cases with CNVs (n = 14), Hürthle cell changes (13 of 14 [93%]) and neoplasms (9 of 14 [64%]) were more likely to be seen on surgical follow-up in comparison with the 17 cases without CNVs (Hürthle cell changes, 6 of 17 [35%]; neoplastic, 3 of 17 [18%]; P < .05).

CONCLUSIONS

CNVs identified in thyroid FNA cases are associated with Hürthle cell morphology and are suggestive of a neoplasm with Hürthle cell features in thyroid FNAs classified as FLUS-HCT/HCN. This finding may be helpful in triaging patients who would benefit from surgical management.

Genomic profile of columnar cell variant of papillary thyroid carcinoma

BACKGROUND AND AIMS

Columnar cell variant (CCV) is a rare papillary thyroid carcinoma subtype. The majority of CCV occur in older patients and are large, invasive tumours that pursue an aggressive clinical course. Rare well-circumscribed CCV occur in younger female patients and are comparatively indolent.

METHODS AND RESULTS

We retrospectively identified CCV with material available to perform targeted next-generation sequencing and correlated molecular results with clinicopathological features and outcome. Our cohort was comprised of nine CCV. Nearly all were aggressive tumours; however, one was predominantly well-circumscribed and arose in a thyroglossal duct cyst of a 26-year-old woman who had no evidence of disease at last follow-up. Seven (78%) cases demonstrated activating oncogenic driver alterations in BRAF, including BRAF V600E, an activating N486_P490del deletion, and BRAF-AGK fusions. Activating RAS mutations were seen in two (22%) cases. Additionally, three (33%) cases had TERT promoter mutations, four (44%) had loss of the tumour suppressor CDKN2A and one (11%) case had a loss of function TP53 mutation. Most cases (89%) also demonstrated copy number alterations, including recurrent gain of chromosome 1q (five cases) and losses of chromosome 9p (three cases) and 22q (four cases). The one case without secondary pathogenic mutations or copy number alterations was the tumour in the 26-year-old woman.

CONCLUSIONS

We found that CCV is primarily a BRAF-driven tumour, with most also harbouring secondary oncogenic mutations and multiple chromosomal gains and losses. Moreover, our findings suggest that molecular analysis could potentially be used to help risk stratify CCV.

DNA FISH Diagnostic Assay on Cytological Samples of Thyroid Follicular Neoplasms

Simple Summary

Cytopathology cannot distinguish benign from malignant follicular lesions in 20–30% of cases. These indeterminate cases includes the so-called follicular neoplasms (FNs) according to The Bethesda System for Reporting Thyroid Cytopathology. Frozen samples from 66 classic follicular adenomas (cFAs) and carcinomas (cFTCs) studied by array-comparative genomic hybridization identified three specific alterations of cFTCs (losses of 1p36.33-35.1 and 22q13.2-13.31, and gain of whole chromosome X) confirmed by fluorescent in situ hybridization (FISH) in a second independent series of 60 touch preparations from frozen samples of cFAs and cFTCs. In a third independent set of 27 cases of already stained pre-operative fine-needle aspiration cytology samples diagnosed as FNs and histologically verified, FISH analysis using these three markers identified half of cFTCs. Specificity of our assay for identifying cFTCs is higher than 98% which might be comparable with BRAF^600E testing in cases of suspicion of classic papillary thyroid carcinomas.

Abstract

Although fine-needle aspiration cytology (FNAC) is helpful in determining whether thyroid nodules are benign or malignant, this distinction remains a cytological challenge in follicular neoplasms. Identification of genomic alterations in cytological specimens with direct and routine techniques would therefore have great clinical value. A series of 153 cases consisting of 72 and 81 histopathologically confirmed classic follicular adenomas (cFAs) and classic follicular thyroid carcinomas (cFTCs), respectively, was studied by means of different molecular techniques in three different cohorts of patients (pts). In the first cohort (training set) of 66 pts, three specific alterations characterized by array comparative genomic hybridization (aCGH) were exclusively found in half of cFTCs. These structural abnormalities corresponded to losses of 1p36.33-35.1 and 22q13.2-13.31, and gain of whole chromosome X. The second independent cohort (validation set) of 60 pts confirmed these data on touch preparations of frozen follicular neoplasms by triple DNA fluorescent in situ hybridization using selected commercially available probes. The third cohort, consisting of 27 archived cytological samples from an equal number of pts that had been obtained for preoperative FNAC and morphologically classified as and histologically verified to be follicular neoplasms, confirmed our previous findings and showed the feasibility of the DNA FISH (DNA fluorescent in situ hybridization) assay. All together, these data suggest that our triple DNA FISH diagnostic assay may detect 50% of cFTCs with a specificity higher than 98% and be useful as a low-cost adjunct to cytomorphology to help further classify follicular neoplasms on already routinely stained cytological specimens.

Thyroid Carcinomas That Occur in Familial Adenomatous Polyposis Patients Recurrently Harbor Somatic Variants in APC, BRAF, and KTM2D

Background: Familial adenomatous polyposis (FAP) is a condition typically caused by pathogenic germline mutations in the APC gene. In addition to colon polyps, individuals with FAP have a substantially increased risk of developing papillary thyroid cancer (PTC). Little is known about the events underlying this association, and the prevalence of somatic "second-hit" mutations in APC is controversial. Methods: Whole-genome sequencing was performed on paired thyroid tumor and normal DNA from 12 FAP patients who developed PTC. Somatic mutation profiles were compared with clinical characteristics and previously sequenced sporadic PTC cases. Germline variant profiling was performed to assess the prevalence of variants in genes previously shown to have a role in PTC predisposition. Results: All 12 patients harbored germline mutations in APC, consistent with FAP. Seven patients also had somatic mutations in APC, and seven patients harbored somatic mutations in KMT2D, which encodes a lysine methyl transferase. Mutation of these genes is extremely rare in sporadic PTCs. Notably, only two of the tumors harbored the somatic BRAF p.V600E mutation, which is the most common driver mutation found in sporadic PTCs. Six tumors displayed a cribriform-morular variant of PTC (PTC-CMV) histology, and all six had somatic mutations in APC. Additionally, nine FAP-PTC patients had rare germline variants in genes that were previously associated with thyroid carcinoma. Conclusions: Our data indicate that FAP-associated PTCs typically have distinct mutations compared with sporadic PTCs. Roughly half of the thyroid cancers that arise in FAP patients have somatic "second-hits" in APC, which is associated with PTC-CMV histology. Somatic BRAF p.V600E variants also occur in some FAP patients, a novel finding. We speculate that in carriers of heterozygous pathogenic mutations of tumor suppressor genes such as APC, a cooperating second-hit somatic variant may occur in a different gene such as KTM2D or BRAF, leading to differences in phenotypes. The role of germline variance in genes other than APC (9 of the 12 patients in this series) needs further research.

Coding Molecular Determinants of Thyroid Cancer Development and Progression

Thyroid cancer is the most common endocrine malignancy. Its incidence and mortality rates have increased for patients with advanced-stage papillary thyroid cancer. The characterization of the molecular pathways essential in thyroid cancer initiation and progression has made huge progress, underlining the role of intracellular signaling to promote clonal evolution, dedifferentiation, metastasis, and drug resistance. The discovery of genetic alterations that include mutations (BRAF, hTERT), translocations, deletions (eg, 9p), and copy-number gain (eg, 1q) has provided new biological insights with clinical applications. Understanding how molecular pathways interplay is one of the key strategies to develop new therapeutic treatments and improve prognosis.

INDEX-db: The Indian Exome Reference Database (Phase I)

Deep sequencing-based genetic mapping has greatly enhanced the ability to catalog variants with plausible disease association. Confirming how these identified variants contribute to specific disease conditions, across human populations, poses the next challenge. Differential selection pressure may impact the frequency of genetic variations, and thus detection of association with disease conditions, across populations. To understand genotype to phenotype correlations, it thus becomes important to first understand the spectrum of genetic variation within a population by creating a reference map. In this study, we report the development of phase I of a new database of genetic variations called INDian EXome database (INDEX-db), from the Indian population, with an aim to establish a centralized database of integrated information. This could be useful for researchers involved in studying disease mechanisms at clinical, genetic, and cellular levels.

A panel of four genes accurately differentiates benign from malignant thyroid nodules

Background

Clinicians are confronted with an increasing number of patients with thyroid nodules. Reliable preoperative diagnosis of thyroid nodules remains a challenge because of inconclusive cytological examination of fine-needle aspiration biopsies. Although molecular analysis of thyroid tissue has shown promise as a diagnostic tool in recent years, it has not been successfully applied in routine clinical use, particularly in Chinese patients.

Methods

Whole-transcriptome sequencing of 19 primary papillary thyroid cancer (PTC) samples and matched adjacent normal thyroid tissue (NT) samples were performed. Bioinformatics analysis was carried out to identify candidate diagnostic genes. Then, RT-qPCR was performed to evaluate these candidate genes, and four genes were finally selected. Based on these four genes, diagnostic algorithm was developed (training set: 100 thyroid cancer (TC) and 65 benign thyroid lesions (BTL)) and validated (independent set: 123 TC and 81 BTL) using the support vector machine (SVM) approach.

Results

We discovered four genes, namely fibronectin 1 (FN1), gamma-aminobutyric acid type A receptor beta 2 subunit (GABRB2), neuronal guanine nucleotide exchange factor (NGEF) and high-mobility group AT-hook 2 (HMGA2). A SVM model with these four genes performed with 97.0 % sensitivity, 93.8 % specificity, 96.0 % positive predictive value (PPV), and 95.3 % negative predictive value (NPV) in training set. For additional independent validation, it also showed good performance (92.7 % sensitivity, 90.1 % specificity, 93.4 % PPV, and 89.0 % NPV).

Conclusions

Our diagnostic panel can accurately distinguish benign from malignant thyroid nodules using a simple and affordable method, which may have daily clinical application in the near future.

Electronic supplementary material

The online version of this article (doi:10.1186/s13046-016-0447-3) contains supplementary material, which is available to authorized users.

Genome-wide microarray expression and genomic alterations by array-CGH analysis in neuroblastoma stem-like cells

Neuroblastoma has a very diverse clinical behaviour: from spontaneous regression to a very aggressive malignant progression and resistance to chemotherapy. This heterogeneous clinical behaviour might be due to the existence of Cancer Stem Cells (CSC), a subpopulation within the tumor with stem-like cell properties: a significant proliferation capacity, a unique self-renewal capacity, and therefore, a higher ability to form new tumors. We enriched the CSC-like cell population content of two commercial neuroblastoma cell lines by the use of conditioned cell culture media for neurospheres, and compared genomic gains and losses and genome expression by array-CGH and microarray analysis, respectively (in CSC-like versus standard tumor cells culture). Despite the array-CGH did not show significant differences between standard and CSC-like in both analyzed cell lines, the microarray expression analysis highlighted some of the most relevant biological processes and molecular functions that might be responsible for the CSC-like phenotype. Some signalling pathways detected seem to be involved in self-renewal of normal tissues (Wnt, Notch, Hh and TGF-β) and contribute to CSC phenotype. We focused on the aberrant activation of TGF-β and Hh signalling pathways, confirming the inhibition of repressors of TGF-β pathway, as SMAD6 and SMAD7 by RT-qPCR. The analysis of the Sonic Hedgehog pathway showed overexpression of PTCH1, GLI1 and SMO. We found overexpression of CD133 and CD15 in SIMA neurospheres, confirming that this cell line was particularly enriched in stem-like cells. This work shows a cross-talk among different pathways in neuroblastoma and its importance in CSC-like cells.

Interleukin (IL)-21 promoter polymorphism increases the risk of thyroid cancer in Chinese population

Polymorphisms in Interleukin (IL)-21 have been researched in several cancers, but the association between IL-21 polymorphisms and thyroid cancer remains unclarified. This case-control study explored the role of five tagSNPs (rs12508721C>T, rs907715G>A, rs13143866G>A, rs2221903A>G and rs4833837A>G) in IL-21 gene in thyroid cancer development. IL-21 genotypes were examined in 615 thyroid cancer patients and 600 controls in Chinese population, and the associations with the risk of thyroid cancer were estimated by logistic regression. Moreover, the potential role of rs12508721C>T in thyroid cancer was further explored by biochemical assays. Compared with the rs12508721CC genotype, CT genotype presented a significantly decreased risk of thyroid cancer (adjusted odds ratios [OR]=0.72; 95%CI=0.57-0.94), the TT carriers had a further decreased risk of thyroid cancer (OR=0.56; 95%CI=0.41-0.87). Furthermore, our quantitative real-time PCR and Enzyme-linked immunosorbent assay (ELISA) results demonstrated that the presence of rs12508721T allele led to more IL-21 expression. However, no significant difference was found in genotype frequencies for other four sites between cases and controls. These findings suggested that rs12508721 polymorphism in IL-21 might be a genetic modifier for the development of thyroid cancer.