{'fraction_non_alphanumeric': 0.030517455313726148, 'fraction_numerical': 0.014688621348804453, 'mean_word_length': 5.26207475850483, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Gene expression profile in functioning and non-functioning nodules of autonomous multinodular goiter from an area of iodine deficiency: unexpected common characteristics between the two entities

Purpose Toxic multinodular goiter is a heterogeneous disease associated with hyperthyroidism frequently detected in areas with deficient iodine intake, and functioning and non-functioning nodules, characterized by increased proliferation but opposite functional activity, may coexist in the same gland. To understand the distinct molecular pathology of each entity present in the same gland, the gene expression profile was evaluated by using the Affymetrix technology. Methods Total RNA was extracted from nodular and healthy tissues of two patients and double-strand cDNA was synthesized. Biotinylated cRNA was obtained and, after chemical fragmentation, was hybridized on U133A and B arrays. Each array was stained and the acquired images were analyzed to obtain the expression levels of the transcripts. Both functioning and non-functioning nodules were compared versus healthy tissue of the corresponding patient. Results About 16% of genes were modulated in functioning nodules, while in non-functioning nodules only 9% of genes were modulated with respect to the healthy tissue. In functioning nodules of both patients and up-regulation of cyclin D1 and cyclin-dependent kinase inhibitor 1 was observed, suggesting the presence of a possible feedback control of proliferation. Complement components C1s, C7 and C3 were down-regulated in both types of nodules, suggesting a silencing of the innate immune response. Cellular fibronectin precursor was up-regulated in both functioning nodules suggesting a possible increase of endothelial cells. Finally, Frizzled-1 was down-regulated only in functioning nodules, suggesting a role of Wnt signaling pathway in the proliferation and differentiation of these tumors. None of the thyroid-specific gene was deregulated in microarray analysis. Conclusion In conclusion, the main finding from our data is a similar modulation for both kinds of nodules in genes possibly implicated in thyroid growth.


Introduction
Autonomous or toxic multinodular goiter is a heterogeneous disease-producing hyperthyroidism frequently found in iodine-deficient areas [1]. The term "autonomous" multinodular goiter (AMNG) or "toxic" multinodular goiter (TMNG) encompasses a spectrum of different clinical entities, ranging from a single hyperfunctioning nodule within an enlarged thyroid that also contains non-functioning nodules to multiple hyperfunctioning areas scattered throughout the gland, barely distinguishable from non-functioning nodules and extranodular parenchyma [2,3]. Thyroid scintiscan after the administration of tracer doses of radioactive iodine or 99mTc can classify thyroid P. Agretti and G. De Marco contributed equally to this manuscript.

3
nodules as functioning and non-functioning [1][2][3][4]. While functioning nodules are able to trap radioiodine and are also defined "hot" nodules, non-functioning nodules are those that, compared with the normal thyroid tissue, take up little or no radioiodine and, for this peculiarity, assume the typical "cold" appearance at thyroid scintiscan. Thyroid nodules of TMNG or AMNG may be true adenomas, defined by the presence of a well-formed fibrous capsule, or more commonly hyperplastic lesions consisting of aggregates of micro/macro follicles lacking encapsulation and not clearly delimited from the surrounding parenchyma [4]. The development of nodular goiter is a result of long-term exposure of the thyroid gland to proliferative stimuli, such as iodine deficiency, goitrogen substances or congenital errors in thyroid hormone synthesis, resulting in insufficient thyroid hormone production and stimulation of TSH secretion by the pituitary. TSH determines a short-term upregulation of iodine uptake and organification, thyroglobulin synthesis and T3 and T4 secretion, and a long-term proliferation of follicular cells with enlargement of the thyroid gland [3].
Pivotal studies showed that up to 82% of solitary toxic adenomas harbor activating TSH receptor (TSHr) mutations [5][6][7]. Subsequent studies reported that activating TSHr mutations are not only present in most solitary toxic adenomas, but also in hyperfunctioning areas (either adenomas or hyperplastic nodules) within a toxic multinodular goiter [8,9]. In particular, somatic mutations determining the constitutive activation of TSHr are reported in about 60% of functioning nodules, while the residual 40% do not harbor TSHr mutations and the genetic mechanism behind remains poorly understood. On the other hand, while a defective iodine transport and organification are implicated in hypofunctionality of non-functioning nodules, the molecular event accounting for the proliferative advantage in these nodules is poorly understood [10]. No specific genetic mutation has been described so far in nonfunctioning nodules, however, a minority of these nodules harbor gene mutations that are also common to some malignant follicular neoplasms (N-RAS, H-RAS, K-RAS mutations or RET rearrangements) [10].
Thyroid tumors result from changes in gene expression patterns that are important for cellular regulatory processes such as growth, differentiation, DNA duplication, mismatch repair and apoptosis. Classification of human tumors into distinct groups based on their origin and histopathological appearance has been the foundation for diagnosis and treatment. This classification is generally based on cellular architecture, and cell-specific antigens only. In contrast, gene expression assays have the potential to identify thousands of unique characteristics for each tumor type [11]. For this reason, we decided to compare the gene expression profile of healthy tissue and functioning and non-functioning thyroid nodules arising in the same gland of patients with AMNG.

In vitro tests
Serum measurement of thyroid parameters was performed in all patients included in the study. FT4 and FT3 were measured by chemiluminescent immunoassay (Vitros System, Ortho-Clinical Diagnostic, Rochester, NY, USA). Thyrotropin was assessed by a sensitive chemiluminescence assay (Immulite 2000, DPC, Los Angeles, CA, USA).

Patients
Two patients subjected to near-total thyroidectomy for AMNG were included in the study. At diagnosis, the patients showed subclinical hyperthyroidism with normal serum FT4 and FT3 and undetectable TSH. No autoimmunity signs were present in both patients with negative AbTG, AbTPO and TRAB. Physical examination, ultrasound, scintiscan imaging using 131 I and histological examination were used to study thyroid glands.
Surgical tissue specimens were carefully dissected, matching the scintiscan with the whole gland laid on the pathologist tray in its proper anatomic orientation. Functioning and non-functioning nodules identified by scintiscan and healthy normal tissue were isolated and used for histological examination and molecular studies.

Total RNA isolation
Tissue specimens were frozen in liquid nitrogen and stored at -80° C until processed for RNA isolation. Total RNA was isolated from approximately 30-40 mg of frozen tissue using TRIzol reagent (Invitrogen, Carlsbad, CA) according to the manufacturer's instructions. The quality of RNA samples was assessed by electrophoresis through denaturing agarose gels and staining with ethidium bromide to visualize the 18S and 28S RNA bands under UV illumination. The extraction yield was quantified spectrometrically at 260 nm. To produce single-strand cDNA, five micrograms of total  RNA for each sample were reverse transcribed for 1 h at  42 °C in a 20 μl reaction volume using 200 units of Superscript II Rnase H − reverse transcriptase (Invitrogen, Carls

In vitro transcription and microarray analysis
Biotin-labeled cRNA was prepared by in vitro transcription using BioArray High Yield RNA Transcript Labeling Kit (ENZO Life Sciences, Inc., Farmingdale, NY) following the standard Affymetrix protocols. Biotinylated cRNAs were cleanup, quantified and 15 μg were fragmented and added to a hybridization cocktail according to Affymetrix procedures.
After quality determination on Test-3 arrays, the samples were hybridized using the GeneChip Hybridization Oven 640 (Affymetrix, Santa Clara, CA) for 16 h at 45 °C to GeneChip Human Genome U133 set (HG-U133A and HG-U133B). The arrays were washed and stained with streptavidin-phycoerythrin using the Affymetrix antibody amplification protocol for eukaryotic targets (EukGE-WS2 protocol) in the GeneChip Fluidics Station 400 (Affymetrix, Santa Clara, CA). Images were acquired using the GeneArray Scanner (Agilent Technologies, Palo, CA), and row data were collected and analyzed by using the Affymetrix Micro Array Suite (MAS) version 5.0 software.

Gene expression analysis
MAS 5 absolute and comparison expression analysis were performed according to Affymetrix GeneChip Expression Analysis Technical Manual's protocols. An absolute expression analysis determines whether the transcripts represented on the probe array are present, absent or marginal in the sample. The comparison analysis identifies relative changes in the expression level of each transcript represented on the arrays. To compare the signals of several arrays, each experiment has to be scaled to the same target intensity (100) to take into account the inherent differences between the arrays and their hybridization efficiencies. Comparisons were made between the functioning and non-functioning nodules and the healthy non-nodular tissue of each patient on HG-U133A and HG-U133B, using as baseline the values obtained for the healthy non-nodular tissue. The output is a change call of increase, marginal increase, decrease, marginal decrease or no change, a p-value associated with the change call, and the intensity of the difference (signal log ratio).

Linking the expression data to biological pathways
DNA microarray experiments simultaneously measure the expression levels of thousands of genes, generating huge amounts of data. We used GenMAPP (Gene MicroArray Pathway Profiler, downloaded from www. GenMA PP. org) program for viewing and analyzing gene expression data in the context of known biological pathways [12,13].
Excel file of MAS obtained gene expression data, including probe set ID, signal log ratio (SLR), SLR low and high, change call and change p value. Expression data were imported into GenMAPP in a csv format (comma separated values) and converted into a gene expression data set that can be viewed on the different MAPPs with specific color codes. When a MAPP is linked to a gene expression data set, GenMAPP automatically and dynamically color codes the genes on the MAPP based on data and criteria provided by the investigator. In this paper, we defined the following color codes: red corresponding to strong increase (SLR ≥ 1.5 and change = increase), pink corresponding to valid increase (SLR low ≥ 0.1 and p value ≤ 0.04), blue corresponding to strong decrease (SLR ≤ -1.5 and change = decrease), skyblue corresponding to valid decrease (SLR ≤ -0.1 and p value ≥ 0.996).

Reverse transcription
One microgram of total RNA for each sample was reverse transcribed for 1 h at 42 °C in a 20 μl reaction volume using 200 units of Superscript II Rnase H − reverse transcriptase (Invitrogen Life Technologies, Carlsbad, CA, USA) in the presence of 1.5 μM random examers (Pharmacia Biotech, Uppsala, Sweden), 0.01 M DTT and 1 mM dNTP mix.

Determination of mRNA levels by using real-time RT-PCR
Quantitative gene expression studies were performed using TaqMan Gene Expression Assays pre-designed primer and probe sets (Applied Biosystems, Foster City, CA). PCR reaction was carried out in 96-well optical reaction plates using a cDNA equivalent of 0.5 ng total RNA for each sample in a volume of 50 μl using the TaqMan Universal PCR Master Mix (Applied Biosystems, Foster City, CA) according to the manufacturer's instructions. PCR was developed on the ABI Prism 7700 Sequence Detector (Applied Biosystems, Foster City, CA). The thermal cycling conditions comprised an initial denaturation step at 95 °C for 10 min and 40 cycles of two-step PCR, including 15 s of denaturation at 95 °C and 1 min of annealing-elongation at 60 °C, using the standard protocol of the manufacturer. Each sample was assayed in triplicate and the intra-assay coefficient of variation was less than 1%. Experiments were repeated three times. The monitoring of negative controls for each target showed an absence of carryover.
To minimize the errors arising from the variation in the amount of starting RNA among samples, amplification of β-actin mRNA was performed as an internal reference against which other RNA values can be normalized. The primers and the probe for the β-actin RNA were purchased from Applied Biosystems, Foster City, CA and the amplification was started from 0.5 ng total RNA.
Normalized results were expressed as the ratio of the pg RNA of the target gene to the pg RNA of the β-actin gene (mean ± SE of three experiments).

Gene expression analysis
For each human array containing about 22,000 transcripts (probe sets) we deleted all the genes with a call no change, and we selected all the increased genes with a signal log ratio ≥ 1 (2-fold increased with respect to the baseline) and all decreased genes with a signal log ratio ≤ − 1 (2-fold decrease with respect to the baseline). In Table 1 the percentage of changed genes selected for the functioning and non-functioning nodules with respect to the non-nodular healthy tissue in patient #1 and #2 on HG-U133A and B are shown. In both patients, we observed a major number of regulated genes (up-regulated or down-regulated) in the functioning nodules with respect to the non-functioning ones. In particular, the number of increased genes in the functioning nodule of patient #2 was about 5-fold greater with respect to the non-functioning nodule of the same patient ( Table 1).
The software GenMAPP allowed us to introduce selected genes in the context of known biological pathways. In Table 2, 29 metabolic pathways inside of which the software placed a total of 109 selected genes with their name, GeneBank accession number, chromosome localization and signal log ratio relative value are shown.

Linking the expression data to biological pathways
Data analysis prominently revealed a regulation in the expression of genes associated with apoptosis (8 genes), cell cycle (9 genes), classical complement activation (6 genes), electron transport chain (7 genes), inflammatory response (5 genes), GPCRs class A rhodopsin-like (13 genes), G-protein signaling (7 genes), peptide GPCRs (10 genes), TGF β signaling (9 genes) and Wnt signaling (11 genes) in functioning and non-functioning nodules of the two patients. In all the above biological pathways we observed a major number of regulated genes in functioning nodules with respect to the non-functioning ones (Fig. 1). The classical complement activation pathway was the one with the greater number of genes which were simultaneously up-or down-regulated in both functioning and non-functioning nodules ( Fig. 2A, B).

Chromosomal localization of upand down-regulated genes
When analyzing the regulated genes in functioning and non-functioning nodules from patient #1 and #2 for their chromosomal location, it became apparent that multiple chromosomes were more frequently involved. When corrected for chromosome size and each chromosome's individual total predicted gene content, the regulated genes in functioning and non-functioning nodules from the two patients

Up-regulated genes
The cellular fibronectin precursor was the gene with the strongest increase of expression (signal log ratio 3.4) in the functioning nodule of patient #1 (Table 2). Only the cyclin-dependent kinase inhibitor 1 (Cip1/ CDKN1A) gene was up-regulated in both functioning and non-functioning nodules of the two patients, while Cip1, the complement component 1 β chain gene (C1qβ), the tryptophanyl-tRNA synthetase gene (TrpRS), the cellular fibronectin precursor gene and cyclin D1 were up-regulated in functioning nodules of the two patients (Table 2). Only Cip1 gene was up-regulated in non-functioning nodules of both patients. Cip1, TrpRS, uncoupling protein 2 gene (UCP2), putative G-protein coupled receptor gene (RAIG1), cellular fibronectin precursor gene and cyclin D1 gene were up-regulated in functioning and non-functioning nodules of patient #1. Again Cip1 gene was up-regulated in both nodules of patient #2 (Table 2).

Down-regulated genes
The protein kinase C β 2 (PKCβ2) was the gene with the strongest decrease of expression (signal log ratio − 3.5) in non-functioning nodule of patient #2 (Table 2).
C1s and complement protein component C7 were downregulated in both functioning and non-functioning nodules of the two patients, while genes down-regulated in functioning nodules of both patients were C1s, complement protein component C7, long-chain acyl-coenzyme A synthetase (FACL1), G-protein β 5 subunit (Gβ5), cAMP-dependent protein kinase subunit RII β (PRKARIIβ) and Frizzled 1. In Several genes appear more than one time because associated with different biological pathways the non-functioning nodules of the two patients the downregulated genes were C1s, complement protein component C7 and complement component C3 α-and β-subunits (C3). C1s, complement protein component C7 gene and prostacyclin (prostaglandin I2) synthase gene were down-regulated in both nodules from patient #1, while many genes (IAP 1, c-Jun, C1s, complement protein component C7, C3, prostaglandin D2 synthase, protein kinase C β 2, chemokine receptor-4, human Epstein-Barr virus-induced G-protein coupled receptor, Q9BXA0, duffy blood group antigen, lymphocyte-specific protein-tyrosine kinase, interleukin 2 receptor γ chain) were down-regulated in both nodules from patient #2 ( Table 2).

Determination of mRNA levels by using real-time RT-PCR
To validate the microarray expression data with an independent method, we carried out real-time quantitative RT-PCR analysis for a sub-set of 12 genes (Table 3) in functioning, non-functioning and healthy tissues from each patient. The Affymetrix signal intensity data were substantially confirmed by real-time PCR data for the subset of genes ( Table 3). The expression levels of the thyroid-specific genes thyrotropin receptor (TSHr), thyroid peroxidase (TPO), thyroglobulin (Tg) and sodium/iodide symporter (NIS) were not deregulated by the microarray analysis. By real-time  (Table 4), no significative differences in TSHr and Tg signals in functioning and non-functioning nodules with respect to the non-nodular thyroid tissue were shown, and an increase of TPO and NIS signals in functioning nodules were observed (Table 4). Real-time RT-PCR is commonly used to measure gene expression because it is also more sensitive than microarrays in detecting small changes in expression even if it requires more input RNA and is less adaptable to high-throughput studies [14]. This is probably the explanation for inconsistent results with respect to microarray expression data for NIS and TPO genes.

Discussion
Autonomous or toxic nodular goiter is the most common form of hyperthyroidism in iodine deficiency areas where aged patients with long-standing non-toxic goiter experience a progressive increase in size and number of thyroid nodules [1,2,15]. With this process, thyroid function may progress from a fully TSH-dependent condition to autonomy and then to overt thyrotoxicosis [15]. In AMNG or TMNG most of the functioning thyroid nodules coexist with nonfunctioning ones [1,2]. Non-functioning thyroid nodules are inactive and less differentiated with respect to functioning ones, and their molecular etiology is unknown. Moreover, malignant transformation is observed in nearly 5-10% of non-functioning thyroid nodules [16]. While functioning nodules growth is mediated by the activation of cAMP/PKA cascade [17] different metabolic signaling pathways may be implicated in the growth and the loss of the ability to trap iodine in nonfunctioning nodules [10].
Gene expression assays are able to identify thousands of unique characteristics for each tumor type and for this we decided to use this technology to compare the gene expression profile between functioning and non-functioning thyroid nodules with respect to healthy tissue arising in the same glands of two patients with autonomous or toxic nodular goiter coming from an area of iodine deficiency. The number of studied patients is certainly small, but significant for a preliminary and descriptive approach.
We demonstrated that the gene expression profile in both functioning and non-functioning thyroid nodules arising in the same thyroid gland was similar, and most of the modulated genes belonged to the same biological pathways. Nevertheless, in functioning nodules of both patients the total number of genes that changed their expression level (up-regulated and down-regulated genes) was greater with respect to the non-functioning ones.
Using the software GenMAPP we identified a total of 109 selected genes, located in 29 metabolic pathways, that increased or decreased at least 2-fold their expression level in nodular tissues. Most of the transcripts were down-regulated in functioning and non-functioning nodules compared with the surrounding tissue, but we focused our interest particularly on genes that were simultaneously up-regulated or down-regulated in both functioning or in both non-functioning nodules of the two patients.
In functioning nodules of both patients, an up-regulation of the cyclin D1 (BC000076) gene was observed. Cyclin D1 acts by complexing with cyclin-dependent kinase cdk4 or cdk6, promoting phosphorilation and inactivation of the tumor suppressor protein pRb and sequential events including the activation of E2F transcription factor. Overexpression of cyclin D1 contributes to the progression of the cell from G1 to S phase. Interestingly, the cyclin-dependent kinase inhibitor 1 (U09579) gene was up-regulated in the functioning and non-functioning nodules of both patients. We speculated that the simultaneous increase of Cip1 in both entities of the two patients very likely reflects a feedback control mechanism of cell proliferation. The product of Cip1 gene is p21 protein that may be the important intermediate by which p53 mediates its role as an inhibitor of cellular proliferation in response to DNA damage; it may bind to and inhibit cyclin-dependent kinase activity, preventing phosphorylation of critical cyclin-dependent kinase substrates and blocking cell cycle progression. A high expression of this protein may indicate a role of the p21/p53 pathway in the proliferation of thyroid nodules. Different reports described the overexpression of cyclin D1 in thyroid cancer. Wang et al. [18] found cyclin positivity in 63% of the follicular variant of papillary carcinoma and in 60% of follicular adenomas by using immunohistochemistry. An immunohistochemical positivity for p21/Cip1 protein has also been shown to be more frequent in well-differentiated thyroid carcinomas than in follicular adenomas [19].
Complement component C1s (M18767) and complement protein component C7 (J03507) genes were downregulated in both functioning and non-functioning nodules suggesting a possible silencing of the non-specific immune response mediated by the complement activity. C1s combines with C1q to form C1, the first component of the classical pathway of the complement system, while C7 is a constituent of the membrane attack complex. Complement component C3 (K02765) gene was down-regulated in both non-functioning nodules and one functioning nodule of the two patients. This protein contains the C3a anaphylatoxin, a mediator of the local inflammatory process; it induces the contraction of smooth muscle, increases vascular permeability and causes histamine release from mast cells and basophilic leukocytes. Complement system activation has been demonstrated by immunohistochemistry and immunoelectron microscopy in thyroid carcinoma and in thyroid follicular adenomas [20]. It has been shown that complement can be activated by direct C4 binding to the CCP-like module of TPO without any mediation by Ig [21]. The significance of a reduction in complement activity in benign thyroid nodules remains unknown. Recently, an important down regulation of lymphocyte-specific genes in toxic adenomas [22] has been demonstrated using microarray analysis, and the presence of few lymphocytes in these nodules was also described [23]. The decreased complement components expression in functioning and non-functioning nodules likely extends this conclusion to macrophages.
The long-chain acyl-CoA synthetase (L09229) gene was down-regulated only in functioning nodules. This protein determines activation of long-chain fatty acids for both syntheses of cellular lipids and degradation via β-oxidation.
Functioning nodules of AMNG or TMNG are mostly due to chronic activation of the cAMP/PKA pathway [17], so it's not amazing to find modifications in the expression of genes of G-protein signaling in these nodules. G protein β 5 subunit (AF017656) and cAMP-dependent protein kinase subunit RII-β (M31158) genes were down-regulated only in functioning nodules. The first protein is required for the GTPase activity, for replacement of GDP by GTP and for G-protein-effectors interaction, while the second one mediates membrane association by binding to anchoring proteins, including MAP2 kinase. A significant down-regulation of PKCβ2 was observed in both functioning and non-functioning nodules of patient #2. This finding is of particular interest in nodules characterized by increased proliferation and decreased activity. It has been shown that long-term stimulation of the PKC pathway with 12-O-tetradecanoylphobol-13-acetate in ovine, porcine and dog cultured thyroid cells causes a general loss of thyroid-specific functions. In particular, it was demonstrated that PKC inhibited TSHmediated human thyroid cell differentiation [24]. Cellular fibronectin precursor (X02761) gene was up-regulated in functioning nodules of both patients. Fibronectins bind cell-surface proteins such as integrins and various intra and extracellular cell components including collagen, fibrin, heparin, DNA and actin. They are involved in cell adhesion, cell motility, opsonization, wound healing and maintenance of cell shape. The fibronectin gene has been found to be upregulated in papillary thyroid carcinoma compared to normal thyroid [25]. Recently, Prasad et al. [26] demonstrated by immunohistochemistry that fibronectin expression is significantly associated with malignancy and is highly specific for carcinoma compared to adenoma [26]. The role of an increased expression of this gene in a benign pathology as the functioning thyroid nodule remains to be elucidated.
Frizzled 1 (AF072872) gene was down-regulated only in functioning nodules of both patients. Frizzled 7 (AB017365) gene and Frizzled 10 (AB027464) were down-regulated only in non-functioning nodules of patient #1. Recently, it has been demonstrated that elements of Wnt/beta-catenin signaling pathway are expressed in thyroid cells (nodular goiter and normal tissue adjacent to thyroid carcinoma) and are functionally active [27]. The highly conserved Wnt signaling pathway regulates cell proliferation, differentiation and cell fate and might play an important role in proliferation, differentiation and, when dysregulated, in thyroid tumorigenesis.
Summarizing, functioning and non-functioning thyroid nodules are basically tumors with increased proliferation rate and the functional hyperactivity of the functioning ones may be somehow a side characteristic of this specific kind of tumor. Although the activation of TSH receptor downstream effectors and target genes are supposed to be a very specific and prominent phenomenon, they first lead to the tumoral phenotype characterized by increased growth.
In conclusion, to our knowledge this is the first study that compares the gene expression profile in functioning and non-functioning thyroid nodules located within the same gland of patients affected by autonomous thyroid multinodular goiter from an area of iodine deficiency. The most surprising result emerged from this study is represented by the similar modulation in genes implicated in thyroid growth in different entities such as functioning and non-functioning nodules. Although obtained from a small number of patients, our results may represent a hint to understand the molecular pathogenesis of benign thyroid nodules. Further proteomic and metabolomic studies will be necessary to validate discoveries made at the genomic level.
Informed consent Informed consent was obtained from all individual participants included in the study.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.