Next Article in Journal
Relevance of Autophagy in Parenchymal and Non-Parenchymal Liver Cells for Health and Disease
Previous Article in Journal
Zonisamide Administration Improves Fatty Acid β-Oxidation in Parkinson’s Disease
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Role for Circular Non-Coding RNAs in the Pathogenesis of Sporadic Parathyroid Adenomas and the Impact of Gender-Specific Epigenetic Regulation

1
1st Department of Propaedeutic Internal Medicine, National and Kapodistrian University of Athens, School of Medicine, Laikon General Hospital, Athens 11527, Greece
2
Pathology Department, Faculty of Medicine, Aristotle University of Thessaloniki, Thessaloniki, 54124, Greece
3
3rd Department of Surgery, AHEPA Univ. Hospital, Aristotle University of Thessaloniki, Thessaloniki 54124, Greece
4
Laboratory of Clinical and Molecular Endocrinology, 1st Department of Internal Medicine, Aristotle University of Thessaloniki, Thessaloniki 54124, Greece
5
1st Department of Surgery, AHEPA Univ. Hospital, Aristotle University of Thessaloniki, Thessaloniki, 54124, Greece
*
Author to whom correspondence should be addressed.
Cells 2019, 8(1), 15; https://doi.org/10.3390/cells8010015
Submission received: 15 November 2018 / Revised: 22 December 2018 / Accepted: 24 December 2018 / Published: 30 December 2018

Abstract

:
Epigenetic changes, including altered small non-coding RNAs, appear to be implicated in the pathogenesis of sporadic parathyroid adenomas (PAs). In this study, we investigated the circular RNAs (circRNAs) expression profile in sporadic PAs. Sixteen tissue samples of sporadic PAs, and four samples of normal parathyroid tissue (NPT) were investigated. Sample preparation and microarray hybridization were performed based on the Arraystar’s standard protocols, and circRNAs sequences were predicted by bioinformatics tools. We identified 35 circRNAs that were differentially expressed in sporadic PAs compared to NPT; 22 were upregulated, and 13 were downregulated, according to the pre-defined thresholds of fold-change > 2.0 and p < 0.05. In the subgroup analysis of PAs from male patients (n = 7) compared to PAs from female patients (n = 9), we also find a different expression profile. In particular, 19 circRNAs were significantly upregulated, and four circRNAs were significantly downregulated in male patients, compared to female counterparts. We show here for the first time a differential circRNA expression pattern in sporadic PAs compared to NPT, and a different expression profile in PA samples from male compared to female patients, suggesting an epigenetic role in the PA pathogenesis, and also an effect of gender in the epigenetic regulation of PAs.

1. Introduction

Primary hyperparathyroidism (PHP) is considered to be one of the most common endocrinopathies with a prevalence of 0.3% in the general population [1] showing a predilection for postmenopausal women [2]. PHP is characterized by hypercalcemia, hypercalciouria, hypophosphatemia, and high or inappropriately normal levels of serum PTH. The increased secretion of PTH from parathyroid cells in PHP is usually due to a benign monoclonal expansion of parathyroid cells, leading to parathyroid adenoma (PA) or hyperplasia, while malignancy rarely occurs and accounts for less than 1% of the cases [3].
The pathophysiology of parathyroid tumorigenesis involves inactivating germline mutations of suppressive oncogenes in the majority of cases with familial predisposition. In sporadic PAs, several genes, including multiple endocrine neoplasia type 1 (MEN1), cyclin D1/PRAD1, and cyclin-dependent kinase inhibitors (CDKI), such as CDKN1B/p27, have also been implicated in the pathogenesis [4].
Apart from alterations in MEN1, however, genetic abnormalities in other genes appear to occur very rarely, and thus the role of epigenetic changes in parathyroid tumorigenesis is under intensive investigation.
Non-coding RNAs, including long–non-coding RNAs (lncRNAs), short micro-RNAs (miRNAs), and circular RNAs (circRNAs), control various levels of gene expression in physiology and disease, such as chromatin architecture/epigenetic memory, transcription, RNA splicing, editing, translation, and possibly tumorigenesis [5]. CircRNAs are a relatively newly-discovered type of non-coding RNAs that are formed from the covalent linkage of the 3’ and 5’ ends to form a closed loop [6]. As a result of this closed structure, circRNAs have been shown to be highly stable and largely resistant to RNA degradation pathways [7], which suggests that circRNAs may be more useful molecular biomarkers for human diseases than linear non-coding RNAs.
Similar to miRNAs, several circRNAs have shown aberrant expression in human neoplastic diseases. As such, the expression of hsa_circ_0001141 and hsa_circ_001763 has been shown to be significantly silenced in squamous cell carcinoma of the esophagus and colorectal cancer tumors, while the expression of hsa_circ_002059 and hsa_circ_0001649 has been shown to be significantly down-regulated in gastric malignancies, and in hepatocellular carcinoma tumors, respectively [7]. Data are scarce, however, on the expression of circRNAs in parathyroid tissue. To address this issue, we profiled the circRNA tissue expression of sporadic PAs compared to normal parathyroid tissues, in order to gain a better insight on the role of epigenetic changes in the pathogenesis of sporadic PAs.

2. Materials and Methods

2.1. Sample Selection

Parathyroid tissue samples, alongside with clinical and histopathological data, were collected from patients that underwent parathyroidectomy (PTX) from December 2014 to February 2017. Inclusion criteria were: single gland disease, occurring in adults with no family history of primary hyperparathyroidism nor a personal/family history suggestive of multiple endocrine neoplasia or a related syndrome, and no evidence of atypical or malignant features. Cases of parathyroid carcinomas according to Schantz and Castleman’s histologic criteria [8] and parathyroid hyperplasia were excluded from the final analysis. Sixteen patients with PAs that fulfilled the inclusion/exclusion criteria, and normal parathyroid tissue (NPT) samples from four patients, which were obtained at the time of total thyroidectomy for benign multinodular goiter, were finally included in the analysis (Figure 1).
Tissue samples were not suitable for one of the following reasons: no tumor cells present in tumor samples (n = 7). Tissue samples were withdrawn for one of the following reasons: patients with atypical features (n = 3), analysis unsuccessful (n = 4).
We performed an immunohistochemical evaluation of the adenomatous polyposis coli (APC) gene, Ki67, cyclin D1, and parafibromin expression, in order to definitely establish a diagnosis of PA, and to exclude cases of atypical PA and/or parathyroid carcinoma [9,10,11]. The study was approved by the Scientific Committee of the affiliated AHEPA University Hospital.

2.2. Histopathological and Immunohistochemical Analysis

Tissue samples from all patients were collected, fixed in formaldehyde 10% and embedded in paraffin blocks (Formalin Fixed, Paraffin Embeded, FFPE). Sections (3 μm) from all the blocks were taken and stained with Hematoxylin–Eosin (H&E) for morphological evaluation and verification of adenoma or normal PTH, accordingly. In addition, the H&E slides were used for marking the epithelial rich areas for microdissection and tumor microarray (TMA) construction. TMAs constructed after the necessary slides for the below described microdissection were obtained. From each patient’s block, three cores were obtained, each 1 mm in diameter. The cores were taken from areas with high ratio of epithelial to stromal component. From the TMA blocks, one slide was taken and stained with H&E, in order to confirm that there was adequate tissue. Another four sections (3 μm each) from each TMA block was taken and stained for Ki67 (MIB1 clone, DAKO, Glostrup, Denmark), CyclinD1 (EP12 clone, DAKO), APC (C-20 clone, Santa Cruz Biotechnology Inc., Dallax, TX, USA), and parafibromin (2H1 clone, Santa Cruz Biotechnology Inc.). All immunohistochemical stains were performed by using Bond Max and Bond III autostainers (Leica Microsystems, Wezlar, Germany). Parafibromin was evaluated as “positive” for any percentage of nuclear positivity, and “negative” for no nuclear positivity. For APC, cytoplasmic positivity was evaluated as a percentage of positive cells. For CyclinD1 and Ki67, nuclear positivity was evaluated as a percentage of positive cells. Each stain was evaluated in three cores, and the average value was taken into account for statistical analysis.

2.3. RNA Extraction from Paraffin-Embedded Tissue

Manual microdissection was performed from each FFPE block to enrich for epithelial tissue before total RNA extraction. In brief, one H&E slide and up to five unstained slides (10 μm thickness) were generated from each FFPE block. The target lesion was marked on the H&E slide for each case, which was then used to guide the removal of non-target tissues (e.g., normal parathyroid tissue) from unstained slides. After an overnight deparaffinization, the tissue area of interest was scraped off the slide into a microcentrifuge tube to facilitate total RNA extraction. RNA extraction was performed by miRNeasy FFPE kit (Qiagen GmbH., Hilden, Germany), according to manufacturer’s instructions. The kit is optimized to isolate RNA molecules that are longer than 18 nucleotides from FFPE tissue samples, and to reverse as much formaldehyde modification as possible without further RNA degradation.
The concentrations of the RNA samples were determined by OD260 by using a NanoDrop ND-1000 instrument. The integrity of RNA was assessed by electrophoresis on a denaturing agarose gel.

2.4. Sample Labeling and Microarray Hybridization

The sample preparation and microarray hybridization were performed based on the Arraystar’s standard protocols (Arraystar Inc., Rockville, MD, USA). All samples were used for microarray analysis.
Briefly, total RNA from each sample was digested with Rnase R (RNase R; Epicenter, Madison, WI, USA) to remove linear RNAs, and to enrich circular RNAs. Then, the enriched circular RNAs were amplified and transcribed into fluorescent complement RNA utilizing a random priming method (Arraystar Super RNA Labeling Kit; Arraystar, Rockville, MD, USA). The fluorescently-labeled cRNAs were quantified with a NanoDrop ND-1000 spectrophotometer (NanoDrop) and hybridized onto the Arraystar Human circRNA Array V2 (8x15K, Arraystar) (Figure S1). After having washed the slides, the arrays were scanned by the Agilent Scanner G2505C. Agilent Feature Extraction software (version 11.0.1.1) was used to analyze acquired array images. CircRNA sequences were predicted by bioinformatic methods [12].

2.5. Quantitative Real Time Reverse Transcription Polymerase Chain Reaction (qRT-PCR)

Following RNA extraction from parathyroid adenomas and normal parathyroid tissue cDNA was synthesized with a reverse transcriptase according to the kit’s instructions (miScript II RT Kit, Qiagen). In brief, we used 1 ug RNA, polyadenylated by poly(A) polymerase and converted it into cDNA by reverse transcriptase with oligo-dT priming using the miScript HiFlex Buffer. The 18S and 28S ribosomal RNAs were used as internal standards. Cycling was performed under standardized conditions with 2× QuantiTect® SYBR Green PCR Master Mix on the QIAGEN Rotor-Gene Q (Corbett Rotor-Gene 6000, Hilden, Germany) real-time PCR cycler. Relative circRNA expression was calculated with the 2−ΔΔCt method.

2.6. Prediction of circRNA–miRNA Interactions

Circular RNAs have been shown to affect miRNA-mediated regulation of gene expression through miRNA sponging [13]. Therefore, the target prediction software, based on Circbase 2.0, TargetScan and Circular RNA Interactome (NIH), was applied to predict putative circRNA/miRNA interactions for the differentially expressed circRNAs in PAs compared to NPT, identified from the microarray hybridization. The Circ-Interactome is a computational tool that enables the prediction and mapping of binding sites for RNA binding proteins (RBP) and miRNAs on reported circRNAs. This tool searches mammalian circRNAs’ names and the genomic position and the best-matching transcripts of the circRNA. Then, it retrieves genomic and mature circRNA sequences and searches RBPs binding to a circRNA, and to sequences upstream/downstream of the circRNA. Finally, it identifies RBPs binding to the circRNA junctions and miRNAs targeting a circRNA [14].

Prediction of circRNA–miRNA–Target Gene Associations

Several online bioinformatics platforms such as TargetScan, CircBase, MirGator, MirBase, CircNet, and NIH Circ interactome were applied in conjunction in order to predict circRNA- miRNA-target gene associations.

2.7. Statistical Analysis

Quantile normalization and subsequent data processing were performed, using the R software limma package. Differentially expressed circRNAs with statistical significance between two groups were identified through Volcano Plot filtering. Differentially expressed circRNAs between two samples were identified through Fold Change filtering. Hierarchical Clustering was performed to show the distinguishable circRNAs expression pattern among samples. When comparing the two groups (PAs versus control), for profile differences, we computed the “fold change’’ between the groups for each circRNA. circRNAs having more than 2-fold changes and p-values < 0.05 were selected as being significantly differentially expressed.

3. Results

Mean age of the included patients was 60 ± 8 years, including nine postmenopausal females and seven males. The control group included samples of NPT from two males and two females mean age 65.7 ± 3.7 years. The demographic and clinical characteristics of the included patients are depicted in Table 1.
As expected, patients with PHP had significantly higher levels of serum calcium and PTH levels compared to controls (sCa: 2.7 ± 0.2 mmol/L vs. 2.28 ± 0.09 mmol/L, p < 0.001, and sPTH: 12.68 ± 3.5 pmol/L vs. 4.42 ± 0.57 pmol/L, p < 0.001, respectively). None of the patients had significant level of hypophosphatemia, although serum phosphate was in the low normal range. After parathyroidectomy, there was a complete biochemical remission of primary hyperparathyroidism in all patients with serum PTH, calcium, and PO4 reaching normal values within one to two months post-surgery.

3.1. Immunohistochemical Findings

Representative images of H&E-stained parathyroid tumor specimens are provided in Figure 2. All tested samples, both PAs and NPT were positive for parafibromin (Figure 2, Table S1). Loss of APC expression was found in one case of PA, while all other cases were positive. Ki67 expression was low in all cases of PAs, reaching up to 1%, and negative in all NPT samples. Expression of cyclinD1 was significantly higher in PAs compared to NPT, as previously described [15,16].
Immunohistochemical findings confirmed the diagnosis of PAs.

3.2. Identification of Differential circRNA Profiles

The Arraystar Human circRNA Microarray was used to sample labeling and microarray hybridization of all the finally selected parathyroid-tissue samples (Figure 3).
We identified 35 circRNAs that were differentially expressed in sporadic PAs compared to normal parathyroid tissue; 22 were upregulated, and 13 were downregulated (Table 2), according to the pre-defined thresholds of fold-change > 2.0 and p < 0.05. Among upregulated circRNAs; hsa_cirRNA_051778, hsa_cirRNA_402533, hsa_cirRNA_406174, and hsa_cirRNA_0008267 showed the highest fold-change (>5-fold), whereas hsa_circRNA_032603 and hsa_circRNA_058097 showed the highest fold-change (>3-fold) among the downregulated circRNAs.

3.3. qRT-PCR Validation of Differential circRNAs

The four significantly upregulated circRNAs and the two significantly downregulated circRNAs with the highest fold change were selected for qRT-PCR validation in the same samples that were used for the microarray analysis. The qRT-PCR results revealed that all six circRNAs but one (namely hsa_circ_RNA_402533) were significantly differentiated between PAs vs. normal parathyroid tissue, with fold-changes mirroring those derived from the microarray data (Figure 4).

3.4. Construction of circRNA-miRNA Network

As circRNAs interact with miRNAs via miRNA response elements (MREs), we searched for putative MREs through Arraystar’s circRNA target prediction software, and CircBase, CircNet, and NIH Circ interactome. The predicted miRNAs for the significantly-upregulated and downregulated circRNAs are listed in Table 2, with hsa-miR_1248 being the most frequently targeted miRNA by 5 up-regulated circRNAs (hsa_circRNA_404337, hsa_circRNA_100965, hsa_circRNA_101283, hsa_circ-RNA_007273, hsa_circ-RNA_406174) and one downregulated (hsa_circRNA_021732). Other less frequently targeted miRNAs, include hsa_miR-139-5p, targeted by four upregulated circRNAs (hsa_circRNA_102904, hsa_circRNA_102903, hsa_circRNA_101283, hsa_circRNA_002617), hsa_miR_1264, targeted by three upregulated circular RNAs (hsa -circRNA_404337, hsa_circRNA_100965, hsa_circRNA_002127) and hsa_miR-486-3p targeted by three down-regulated circRNAs (hsa_circRNA_404643, hsa_cirRNA_105038, hsa_cirRNA_032603).
We also searched for the total miRNA binding site (miRBS) that has been reported so far for the three most upregulated and the two most down-regulated circRNAs. We identified 46 miRBS for hsa_circRNA_008627; 41 for hsa_circRNA_058097; 18 for hsa_circRNA_406174; seven miRBS for hsa_circRNA_032603; one for hsa_circRNA_402533.

3.5. circRNA–Genes Interactions

To further explore the underlying molecular mechanism(s) of the differentially expressed circRNAs in sporadic PAs, we used the bioinformatics platforms CircBase, and NIH Circ-Interactome to predict circRNA–genes interactions. The predicted genes associated with the differentially expressed circ-RNAs are depicted in Table 2. Two circRNAs, hsa_circRNA_404337 and hsa_circRNA_051799, interact with genes that have been shown to associate with molecular pathways that are involved in the pathogenesis of sporadic PAs, such as CDKN1B and CDK1, respectively [17,18].
Most of the genes identified are protein-coding, and associated with several cellular processes, such as organization and re-arrangement of the cytoskeleton, cell proliferation, cell metabolism, and transcription (Table S2).

3.6. Parathyroid Adenomas from Male vs. Female Patients

In order to test whether the gender of the patient has an impact in the expression profile of circRNAs in PAs, since the vast majority of PHP due to PAs appear in postmenopausal women, we performed a sub-group analysis of PAs from male patients compared to PAs from female patients. We found a different expression profile between male and female patients. In particular 19 circRNAs were significantly upregulated and four circRNAs were significantly downregulated in male patients compared to female counterparts (Table 3). It is worth noting that the identified differentially expressed circRNAs in the subgroup analysis between male and female patients were different from the ones that we have identified in our initial analysis between PA samples and controls, and there was no overlap in circRNAs expression profiles between the two analyses.
The predicted miRNAs for the significantly upregulated and downregulated circRNAs are also listed in Table 3, with hsa_miR_1184 being the most frequently targeted miRNA by 4 up-regulated circRNAs (hsa_circRNA_008636, hsa_circRNA_000763, hsa_circRNA_092388, hsa_cirRNA_044837) and one down-regulated circRNA (hsa_circRNA_100332). The predicted genes associated with the differentially expressed circRNAs in this analysis are also depicted in Table 3. The majority of these genes are protein coding (Supplemental Table S3) and involves the intracellular ubiquitin proteasome system, nucleic acid and RNA binding, protein heterodimerization and homodimerization activity, vesicle trafficking and alternative splicing of pre-mRNAs. None of the implicated genes, however, is located in the X chromosome (Table 3), and thus we cannot address whether incomplete inactivation of the X chromosome is a possible mechanism, contributing to the altered expression of some circRNAs in males vs. females.
In order to further evaluate our results in the differential expression of sporadic PAs in male vs. female patients, we also conducted a subgroup analysis between samples from males and females in the control group, to test whether the altered expression of circRNAs profile is disease-independent. We did not find any significant difference in this analysis, suggesting that the differences in circRNA expression between male and female patients may be attributed to PA, but the numbers were small.

4. Discussion

In this study, we investigated the circRNA expression profile in sporadic PAs, and we tested the effect of gender.
Thirty-five circRNAs, 22 up-regulated, and 13 down-regulated were identified by microarray analysis to be differentially expressed in sporadic PAs compared to NPT. Applying a set of bioinformatic tools, we constructed a network of circRNA–miRNA and genes interactions, identifying two of the differentially expressed circRNAs, hsa_circRNA_404337 and hsa_circRNA_051799, that interact with genes known to associate with molecular pathways involved in the pathogenesis of sporadic PAs, such as CDKN1B and CDK1 [4,19], respectively. In the sub-analysis within the cohort of sporadic PAs, we also report a differential circRNA expression profile in sporadic PAs from male patients compared to PAs from female counterparts with 19 circRNA being significantly upregulated and four being significantly downregulated. Primary hyperparathyroidism is most frequently present in postmenopausal women, and although the involvement of estrogen has been suggested, and parathyroid tumors are found to express ERβ1 and ERβ2 [20], the exact underlying pathophysiology remains obscure. The differences in circRNA expression profile between male and female patients, but not between male and female controls, point to a potential gender effect on the PA tumorigenesis through epigenetic regulation.
In the past decades an ever-growing list of diverse non-coding RNAs with functional capacity have been identified in eukaryotic cells [21]. CircRNAs are formed by alternative splicing of pre-mRNA, in a process known as ‘backsplicing’, where an upstream splice acceptor is joined to a downstream splice donor [22]. They are expressed in a variety of eukaryotic organisms, demonstrate conservation across mammals, they are abundant, expressed in a regulated manner independent of their cognate linear isoform and in a cell type-specific manner.
Dysregulation of circRNAs profile has been identified in various types of cancers, such as papillary thyroid carcinoma [23], gastric cancer [24], hepatocellular cancer [25,26,27], colorectal cancer [28,29,30], bladder cancer [31], and lung cancer [32], as well as in non-cancer conditions including cardiovascular disease [33], Alzheimer’s disease [34], and diabetes [35,36]. The implication of tissue-specific circRNAs, however, in disease pathogenesis, progression, and therapy resistance remains largely unknown, and is currently under intensive investigation [37].
In cancer, circRNAs appear to have a role in cancer development and progression, by interacting with microRNAs. CircRNAs negatively regulate the inhibitory effects of miRNAs upon their target linear mRNAs by directly binding to miRNAs via MREs, thereby essentially functioning as ‘miRNA sponges’ within the cellular RNA network. Another potential function of circRNAs could be to compete with the splicing of an mRNA [38], while some circRNAs have also been implicated in transcriptional or post-transcriptional gene regulation of their host genes [39] (Figure 5). Compared to mRNAs, circRNAs are more stable, due to their natural resistance to exonucleases, with half-lives greater than 24–48 h, which is significantly longer than the half-life of an average mRNA (8-9 h in human cells) [40,41]; this is why they are considered as being much more promising in understanding the pathophysiology of certain diseases and the development of cancer.
The emerging role of the epigenetic component in tumorigenesis has just now started to be identified in parathyroid neoplasia [6,19].
MiRNA profiles have been investigated in parathyroid tissue before, demonstrating dysregulated miRNA expression profile in parathyroid carcinomas [6,42], and showing a global down-regulation (approximately 80%) of miRNAs compared to normal parathyroid tissue. Other studies have also identified significantly up-regulated miRNAs, such as miR-222 [6] miR-503 [6], and miR-517c [43] in parathyroid carcinomas.
The mature form of miRNA_24-1 was studied in the parathyroid hyperplasia of MEN1. It has been shown that miRNA_24-1 is uniquely expressed in MEN1 parathyroid adenoma tissues that conserve the wild-type allele, while it is downregulated in the MEN1 parathyroid adenoma tissues with the 11q13 LOH, as well as in sporadic parathyroid adenomas [44].
We also search for miRNA response elements (MREs) on the 35 differentially expressed circRNAs, which were then used to identify putative miRNAs that were based on their seed-sequence complementarity [13]. Hsa_miR_1248 was the most frequently encountered miRNA, and although little is known about the potential targets and pathways regulated by this miRNA it has been linked to the pathogenesis of breast [45] and non-small cell lung carcinoma [46], autism spectrum disorders [47], asthma [48], and in the aging process [49]. Genome-wide microarray and in vitro studies [48,49], have shown that miRNA-1248 regulates the expression of genes that play a role in interleukin-related inflammatory pathway, as well as, DNA repair pathways.
None of the miRNAs that were previously reported in parathyroid carcinomas, however, was identified in our analysis through interaction with the reported dysregulated circRNAs.
Apart from the dysregulated expression of non-coding RNAs, the hypermethylation of APC, the gene responsible for familial adenomatous polyposis, occurs frequently in parathyroid adenomas although aberrant APC expression has not been demonstrated [50,51]. In addition, genes that are frequently subject to epigenetic inactivation in human cancers, such as Ras association domain-containing protein 1 (RASSF1A) [52] and hypermethylated in cancer 1 protein (HIC1) [53], are also frequently hypermethylated in parathyroid adenomas, but it is unclear which of these aberrantly methylated genes may be important to the pathogenesis of parathyroid tumors.
Our study has several limitations; (i) the sample size was relatively small, (ii) our analysis was limited to parathyroid tissue, and thus results from other biosamples such as serum or plasma before and after parathyroidectomy that could strengthen our results are missing and (iii) we were not able to perform DNA analysis and investigation of somatic mutations in our samples, mainly due to inadequate material in the total of PA samples used.
Despite these limitations, however, this is the first study investigating the circRNA expression profile in sporadic PAs. We identified 35 differentially expressed circRNAs in PA samples compared to normal parathyroid tissue, suggesting a role of circRNAs in the pathogenesis of these benign tumors. In addition, we reported significant differences in the circRNA expression profile between male and female patients underlying the fact that epigenetics may be involved in the gender specific differences reported in PAs. Furthermore, our results from the bioinformatic analysis reveal a circRNA–miRNA network that could play a key role in the triggering of the tumorigenetic process in the parathyroid cells. As circRNAs regulate the epigenome in multiple ways (Figure 5), future functional studies that can pave the way between ‘circnome’ and PA pathogenesis would significantly increase our knowledge for the molecular pathways in involved.
CircRNA network appears to be cell and tissue specific [37], stable, and with a significant regulatory role in pathogenesis of diseases, providing new strategic approaches to tumorigenesis, and emerge as promising diagnostic and prognostic biomarkers, and as therapeutic targets.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/2073-4409/8/1/15/s1.

Author Contributions

Conceptualization, M.P.Y. and J.G.Y.; methodology, M.P.Y., C.P., N.M., A.G., S.C., S.M., T.P., E.D.; validation, A.G., E.S. and E.D..; resources, K.K., P.Z.; writing—original draft preparation, M.P.Y. and C.P.; writing—review and editing, P.Z., I.K., J.G.Y. supervision, J.G.Y.

Funding

This research received no external funding.

Acknowledgments

We would like to acknowledge Olga Tsave, for her valuable assistance at all steps of the laboratory work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Adami, S.; Marcocci, C.; Gatti, D. Epidemiology of primary hyperparathyroidism in Europe. J. Bone Miner. Res. 2002, 17 (Suppl. S2), N18–N23. [Google Scholar]
  2. Lundgren, E.; Hagstrom, E.G.; Lundin, J.; Winnerback, K.; Roos, J.; Ljunghall, S.; Rastad, J. Primary hyperparathyroidism revisited in menopausal women with serum calcium in the upper normal range at population-based screening 8 years ago. World J. Surg. 2002, 26, 931–936. [Google Scholar] [CrossRef] [PubMed]
  3. Fernandez-Ranvier, G.G.; Khanafshar, E.; Jensen, K.; Zarnegar, R.; Lee, J.; Kebebew, E.; Duh, Q.Y.; Clark, O.H. Parathyroid carcinoma, atypical parathyroid adenoma, or parathyromatosis? Cancer 2007, 110, 255–264. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Marini, F.; Cianferotti, L.; Giusti, F.; Brandi, M.L. Molecular genetics in primary hyperparathyroidism: the role of genetic tests in differential diagnosis, disease prevention strategy, and therapeutic planning. A 2017 update. Clin. Cases Miner. Bone Metab. 2017, 14, 60–70. [Google Scholar] [CrossRef] [PubMed]
  5. Ling, H.; Fabbri, M.; Calin, G.A. MicroRNAs and other non-coding RNAs as targets for anticancer drug development. Nat. Rev. Drug Discov. 2013, 12, 847–865. [Google Scholar] [CrossRef] [Green Version]
  6. Verdelli, C.; Forno, I.; Vaira, V.; Corbetta, S. Epigenetic alterations in human parathyroid tumors. Endocrine 2015, 49, 324–332. [Google Scholar] [CrossRef] [PubMed]
  7. Chen, Y.; Li, C.; Tan, C.; Liu, X. Circular RNAs: A new frontier in the study of human diseases. J. Med. Genet. 2016, 53, 359–365. [Google Scholar] [CrossRef]
  8. Schantz, A.; Castleman, B. Parathyroid carcinoma. A study of 70 cases. Cancer 1973, 31, 600–605. [Google Scholar] [CrossRef]
  9. Juhlin, C.C.; Nilsson, I.L.; Johansson, K.; Haglund, F.; Villablanca, A.; Hoog, A.; Larsson, C. Parafibromin and APC as screening markers for malignant potential in atypical parathyroid adenomas. Endocr. Pathol. 2010, 21, 166–177. [Google Scholar] [CrossRef]
  10. Kruijff, S.; Sidhu, S.B.; Sywak, M.S.; Gill, A.J.; Delbridge, L.W. Negative parafibromin staining predicts malignant behavior in atypical parathyroid adenomas. Ann. Surg. Oncol. 2014, 21, 426–433. [Google Scholar] [CrossRef]
  11. Hosny Mohammed, K.; Siddiqui, M.T.; Willis, B.C.; Zaharieva Tsvetkova, D.; Mohamed, A.; Patel, S.; Sharma, J.; Weber, C.; Cohen, C. Parafibromin, APC, and MIB-1 Are Useful Markers for Distinguishing Parathyroid Carcinomas from Adenomas. Appl. Immunohistochem. Mol. Morphol. 2017, 25, 731–735. [Google Scholar] [CrossRef] [PubMed]
  12. Salzman, J.; Gawad, C.; Wang, P.L.; Lacayo, N.; Brown, P.O. Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types. PLoS ONE 2012, 7, e30733. [Google Scholar] [CrossRef] [PubMed]
  13. Sand, M.; Bechara, F.G.; Gambichler, T.; Sand, D.; Bromba, M.; Hahn, S.A.; Stockfleth, E.; Hessam, S. Circular RNA expression in cutaneous squamous cell carcinoma. J. Dermatol. Sci. 2016, 83, 210–218. [Google Scholar] [CrossRef] [PubMed]
  14. Dudekula, D.B.; Panda, A.C.; Grammatikakis, I.; De, S.; Abdelmohsen, K.; Gorospe, M. CircInteractome: A web tool for exploring circular RNAs and their interacting proteins and microRNAs. RNA Biol. 2016, 13, 34–42. [Google Scholar] [CrossRef] [PubMed]
  15. Vasef, M.A.; Brynes, R.K.; Sturm, M.; Bromley, C.; Robinson, R.A. Expression of cyclin D1 in parathyroid carcinomas, adenomas, and hyperplasias: A paraffin immunohistochemical study. Mod. Pathol. 1999, 12, 412–416. [Google Scholar] [PubMed]
  16. Hsi, E.D.; Zukerberg, L.R.; Yang, W.I.; Arnold, A. Cyclin D1/PRAD1 expression in parathyroid adenomas: An immunohistochemical study. J. Clin. Endocrinol. Metab. 1996, 81, 1736–1739. [Google Scholar] [CrossRef] [PubMed]
  17. Costa-Guda, J.; Soong, C.P.; Parekh, V.I.; Agarwal, S.K.; Arnold, A. Germline and somatic mutations in cyclin-dependent kinase inhibitor genes CDKN1A, CDKN2B, and CDKN2C in sporadic parathyroid adenomas. Horm. Cancer 2013, 4, 301–307. [Google Scholar] [CrossRef] [PubMed]
  18. Arya, A.K.; Bhadada, S.K.; Singh, P.; Sachdeva, N.; Saikia, U.N.; Dahiya, D.; Behera, A.; Bhansali, A.; Rao, S.D. Promoter hypermethylation inactivates CDKN2A, CDKN2B and RASSF1A genes in sporadic parathyroid adenomas. Sci. Rep. 2017, 7, 3123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Costa-Guda, J.; Arnold, A. Genetic and epigenetic changes in sporadic endocrine tumors: Parathyroid tumors. Mol. Cell. Endocrinol. 2014, 386, 46–54. [Google Scholar] [CrossRef] [Green Version]
  20. Haglund, F.; Ma, R.; Huss, M.; Sulaiman, L.; Lu, M.; Nilsson, I.L.; Hoog, A.; Juhlin, C.C.; Hartman, J.; Larsson, C. Evidence of a functional estrogen receptor in parathyroid adenomas. J. Clin. Endocrinol. Metab. 2012, 97, 4631–4639. [Google Scholar] [CrossRef]
  21. Morris, K.V.; Mattick, J.S. The rise of regulatory RNA. Nat. Rev. Genet. 2014, 15, 423–437. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Barrett, S.P.; Salzman, J. Circular RNAs: Analysis, expression and potential functions. Development 2016, 143, 1838–1847. [Google Scholar] [CrossRef] [PubMed]
  23. Peng, N.; Shi, L.; Zhang, Q.; Hu, Y.; Wang, N.; Ye, H. Microarray profiling of circular RNAs in human papillary thyroid carcinoma. PLoS ONE 2017, 12, e0170287. [Google Scholar] [CrossRef] [PubMed]
  24. Li, P.; Chen, S.; Chen, H.; Mo, X.; Li, T.; Shao, Y.; Xiao, B.; Guo, J. Using circular RNA as a novel type of biomarker in the screening of gastric cancer. Clin. Chim. Acta 2015, 444, 132–136. [Google Scholar] [CrossRef] [PubMed]
  25. Qin, M.; Liu, G.; Huo, X.; Tao, X.; Sun, X.; Ge, Z.; Yang, J.; Fan, J.; Liu, L.; Qin, W. Hsa_circ_0001649: A circular RNA and potential novel biomarker for hepatocellular carcinoma. Cancer Biomark. 2016, 16, 161–169. [Google Scholar] [CrossRef] [Green Version]
  26. Shang, X.; Li, G.; Liu, H.; Li, T.; Liu, J.; Zhao, Q.; Wang, C. Comprehensive circular RNA profiling reveals that hsa_circ_0005075, a new circular RNA biomarker, is involved in hepatocellular crcinoma development. Medicine 2016, 95, e3811. [Google Scholar] [CrossRef] [PubMed]
  27. Yu, L.; Gong, X.; Sun, L.; Zhou, Q.; Lu, B.; Zhu, L. The Circular RNA Cdr1as act as an oncogene in hepatocellular carcinoma through targeting miR-7 Expression. PLoS ONE 2016, 11, e0158347. [Google Scholar] [CrossRef] [PubMed]
  28. Dou, Y.; Cha, D.J.; Franklin, J.L.; Higginbotham, J.N.; Jeppesen, D.K.; Weaver, A.M.; Prasad, N.; Levy, S.; Coffey, R.J.; Patton, J.G.; et al. Circular RNAs are down-regulated in KRAS mutant colon cancer cells and can be transferred to exosomes. Sci. Rep. 2016, 6, 37982. [Google Scholar] [CrossRef] [Green Version]
  29. Huang, G.; Zhu, H.; Shi, Y.; Wu, W.; Cai, H.; Chen, X. cir-ITCH plays an inhibitory role in colorectal cancer by regulating the Wnt/beta-catenin pathway. PLoS ONE 2015, 10, e0131225. [Google Scholar] [CrossRef]
  30. Tang, Q.; Zou, Z.; Zou, C.; Zhang, Q.; Huang, R.; Guan, X.; Li, Q.; Han, Z.; Wang, D.; Wei, H.; et al. MicroRNA-93 suppress colorectal cancer development via Wnt/beta-catenin pathway downregulating. Tumour. Biol. 2015, 36, 1701–1710. [Google Scholar] [CrossRef]
  31. Zhong, Z.; Lv, M.; Chen, J. Screening differential circular RNA expression profiles reveals the regulatory role of circTCF25-miR-103a-3p/miR-107-CDK6 pathway in bladder carcinoma. Sci. Rep. 2016, 6, 30919. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Wan, L.; Zhang, L.; Fan, K.; Cheng, Z.X.; Sun, Q.C.; Wang, J.J. Circular RNA-ITCH suppresses lung cancer proliferation via inhibiting the Wnt/beta-catenin pathway. Biomed. Res. Int. 2016, 2016, 1579490. [Google Scholar] [CrossRef] [PubMed]
  33. Zhao, Z.; Li, X.; Gao, C.; Jian, D.; Hao, P.; Rao, L.; Li, M. Peripheral blood circular RNA hsa_circ_0124644 can be used as a diagnostic biomarker of coronary artery disease. Sci. Rep. 2017, 7, 39918. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Lukiw, W.J. Circular RNA (circRNA) in Alzheimer’s disease (AD). Front. Genet. 2013, 4, 307. [Google Scholar] [CrossRef] [PubMed]
  35. Zhao, Z.; Li, X.; Jian, D.; Hao, P.; Rao, L.; Li, M. Hsa_circ_0054633 in peripheral blood can be used as a diagnostic biomarker of pre-diabetes and type 2 diabetes mellitus. Acta Diabetol. 2017, 54, 237–245. [Google Scholar] [CrossRef] [PubMed]
  36. Xu, H.; Guo, S.; Li, W.; Yu, P. The circular RNA Cdr1as, via miR-7 and its targets, regulates insulin transcription and secretion in islet cells. Sci. Rep. 2015, 5, 12453. [Google Scholar] [CrossRef] [Green Version]
  37. Greene, J.; Baird, A.M.; Brady, L.; Lim, M.; Gray, S.G.; McDermott, R.; Finn, S.P. Circular RNAs: Biogenesis, Function and Role in Human Diseases. Front. Mol. Biosci. 2017, 4, 38. [Google Scholar] [CrossRef] [PubMed]
  38. Ashwal-Fluss, R.; Meyer, M.; Pamudurti, N.R.; Ivanov, A.; Bartok, O.; Hanan, M.; Evantal, N.; Memczak, S.; Rajewsky, N.; Kadener, S. circRNA biogenesis competes with pre-mRNA splicing. Mol. Cell. 2014, 56, 55–66. [Google Scholar] [CrossRef]
  39. Hansen, T.B.; Wiklund, E.D.; Bramsen, J.B.; Villadsen, S.B.; Statham, A.L.; Clark, S.J.; Kjems, J. miRNA-dependent gene silencing involving Ago2-mediated cleavage of a circular antisense RNA. EMBO J. 2011, 30, 4414–4422. [Google Scholar] [CrossRef] [Green Version]
  40. Garneau, N.L.; Wilusz, J.; Wilusz, C.J. The highways and byways of mRNA decay. Nat. Rev. Mol. Cell Biol. 2007, 8, 113–126. [Google Scholar] [CrossRef]
  41. Neff, A.T.; Lee, J.Y.; Wilusz, J.; Tian, B.; Wilusz, C.J. Global analysis reveals multiple pathways for unique regulation of mRNA decay in induced pluripotent stem cells. Genome Res. 2012, 22, 1457–1467. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Rahbari, R.; Holloway, A.K.; He, M.; Khanafshar, E.; Clark, O.H.; Kebebew, E. Identification of differentially expressed microRNA in parathyroid tumors. Ann. Surg. Oncol. 2011, 18, 1158–1165. [Google Scholar] [CrossRef] [PubMed]
  43. Vaira, V.; Elli, F.; Forno, I.; Guarnieri, V.; Verdelli, C.; Ferrero, S.; Scillitani, A.; Vicentini, L.; Cetani, F.; Mantovani, G.; et al. The microRNA cluster C19MC is deregulated in parathyroid tumours. J. Mol. Endocrinol. 2012, 49, 115–124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Luzi, E.; Marini, F.; Giusti, F.; Galli, G.; Cavalli, L.; Brandi, M.L. The negative feedback-loop between the oncomir Mir-24-1 and menin modulates the Men1 tumorigenesis by mimicking the “Knudson’s second hit”. PLoS ONE 2012, 7, e39767. [Google Scholar] [CrossRef] [PubMed]
  45. Tanic, M.; Yanowski, K.; Gomez-Lopez, G.; Rodriguez-Pinilla, M.S.; Marquez-Rodas, I.; Osorio, A.; Pisano, D.G.; Martinez-Delgado, B.; Benitez, J. MicroRNA expression signatures for the prediction of BRCA1/2 mutation-associated hereditary breast cancer in paraffin-embedded formalin-fixed breast tumors. Int. J. Cancer 2015, 136, 593–602. [Google Scholar] [CrossRef] [PubMed]
  46. Xu, J.; Tian, S.; Yin, Z.; Wu, S.; Liu, L.; Qian, Y.; Pei, D.; Gao, W.; Xu, J.; Yin, Y.; et al. MicroRNA-binding site SNPs in deregulated genes are associated with clinical outcome of non-small cell lung cancer. Lung Cancer 2014, 85, 442–448. [Google Scholar] [CrossRef] [PubMed]
  47. Jyonouchi, H.; Geng, L.; Streck, D.L.; Dermody, J.J.; Toruner, G.A. MicroRNA expression changes in association with changes in interleukin-1ss/interleukin10 ratios produced by monocytes in autism spectrum disorders: Their association with neuropsychiatric symptoms and comorbid conditions (observational study). J. Neuroinflammation 2017, 14, 229. [Google Scholar] [CrossRef] [PubMed]
  48. Panganiban, R.P.; Pinkerton, M.H.; Maru, S.Y.; Jefferson, S.J.; Roff, A.N.; Ishmael, F.T. Differential microRNA epression in asthma and the role of miR-1248 in regulation of IL-5. Am. J. Clin. Exp. Immunol. 2012, 1, 154–165. [Google Scholar] [PubMed]
  49. Noren Hooten, N.; Fitzpatrick, M.; Wood, W.H., 3rd; De, S.; Ejiogu, N.; Zhang, Y.; Mattison, J.A.; Becker, K.G.; Zonderman, A.B.; Evans, M.K. Age-related changes in microRNA levels in serum. Aging (Albany NY) 2013, 5, 725–740. [Google Scholar] [CrossRef] [Green Version]
  50. Juhlin, C.C.; Haglund, F.; Villablanca, A.; Forsberg, L.; Sandelin, K.; Branstrom, R.; Larsson, C.; Hoog, A. Loss of expression for the Wnt pathway components adenomatous polyposis coli and glycogen synthase kinase 3-beta in parathyroid carcinomas. Int. J. Oncol. 2009, 34, 481–492. [Google Scholar]
  51. Svedlund, J.; Auren, M.; Sundstrom, M.; Dralle, H.; Akerstrom, G.; Bjorklund, P.; Westin, G. Aberrant WNT/beta-catenin signaling in parathyroid carcinoma. Mol. Cancer 2010, 9, 294. [Google Scholar] [CrossRef] [PubMed]
  52. Juhlin, C.C.; Kiss, N.B.; Villablanca, A.; Haglund, F.; Nordenstrom, J.; Hoog, A.; Larsson, C. Frequent promoter hypermethylation of the APC and RASSF1A tumour suppressors in parathyroid tumours. PLoS ONE 2010, 5, e9472. [Google Scholar] [CrossRef] [PubMed]
  53. Svedlund, J.; Koskinen Edblom, S.; Marquez, V.E.; Akerstrom, G.; Bjorklund, P.; Westin, G. Hypermethylated in cancer 1 (HIC1), a tumor suppressor gene epigenetically deregulated in hyperparathyroid tumors by histone H3 lysine modification. J. Clin. Endocrinol. Metab. 2012, 97, E1307–E1315. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Flow chart detailing the selection of parathyroid adenomas (PA) tissue samples.
Figure 1. Flow chart detailing the selection of parathyroid adenomas (PA) tissue samples.
Cells 08 00015 g001
Figure 2. Immunohistochemical staining of sporadic PA and control samples. Representative images of H&E-stained 40× (A) PA samples, (B) control samples, parafibromin-stained 400×, (C) PA samples, (D) control samples, APC-stained 400× (E) PA samples, (F) control samples, Cyclin D1-stained 400×, (G) PA samples, (H) control samples, ki67-stained 400× (I) PA samples, (J) control samples. PA, parathyroid adenomas; H&E, Hematoxylin-Eosin; APC, adenomatous polyposis coli.
Figure 2. Immunohistochemical staining of sporadic PA and control samples. Representative images of H&E-stained 40× (A) PA samples, (B) control samples, parafibromin-stained 400×, (C) PA samples, (D) control samples, APC-stained 400× (E) PA samples, (F) control samples, Cyclin D1-stained 400×, (G) PA samples, (H) control samples, ki67-stained 400× (I) PA samples, (J) control samples. PA, parathyroid adenomas; H&E, Hematoxylin-Eosin; APC, adenomatous polyposis coli.
Cells 08 00015 g002
Figure 3. Scatter plots and volcano plots to identify differentially-expressed circular RNAs (circRNAs). Scatter plots used to identify differentially-expressed circRNAs in (A) PAs vs. NPT and (B) PAs from F patients vs. PAs from M patients. The y-axis represents the mean normalized circRNA signal values for each comparator group (log2-scaled). The green fold-change lines represent 2.0× fold-changes, so the circRNAs lying above and below these green lines displayed greater than a 2.0-fold upregulation or downregulation. Volcano plots used to identify differentially-expressed circRNAs in (C) PAs vs. NPT and (D) PAs from F patients vs. PAs from M patients. The x-axis represents fold-change values (log2-scaled), while the y-axis represents p-values (−log10-scaled). The green vertical lines correspond to 2.0× upregulation and downregulation, respectively, while the green horizontal line corresponds to a p-value of 0.05. On this basis, the red rectangles represent the differentially-expressed circRNAs of statistical significance. NSV, normalized signal value; PAs, parathyroid adenomas; NPT, normal parathyroid tissue; F, female; M, male.
Figure 3. Scatter plots and volcano plots to identify differentially-expressed circular RNAs (circRNAs). Scatter plots used to identify differentially-expressed circRNAs in (A) PAs vs. NPT and (B) PAs from F patients vs. PAs from M patients. The y-axis represents the mean normalized circRNA signal values for each comparator group (log2-scaled). The green fold-change lines represent 2.0× fold-changes, so the circRNAs lying above and below these green lines displayed greater than a 2.0-fold upregulation or downregulation. Volcano plots used to identify differentially-expressed circRNAs in (C) PAs vs. NPT and (D) PAs from F patients vs. PAs from M patients. The x-axis represents fold-change values (log2-scaled), while the y-axis represents p-values (−log10-scaled). The green vertical lines correspond to 2.0× upregulation and downregulation, respectively, while the green horizontal line corresponds to a p-value of 0.05. On this basis, the red rectangles represent the differentially-expressed circRNAs of statistical significance. NSV, normalized signal value; PAs, parathyroid adenomas; NPT, normal parathyroid tissue; F, female; M, male.
Cells 08 00015 g003
Figure 4. qRT-PCR validation of differentially expressed circRNAs in PAs. qRT-PCR was used to validate the dysregulated expression of the four significantly upregulated and the two significantly downregulated circRNAs that showed the highest fold-expression in PAs (n = 16) vs. normal parathyroid tissue (n = 4) in the microarray analysis. The same samples used for microarray analysis were also used for qPCR validation. The y-axis represents the relative fold change compared to controls. Data are displayed as means ± SEMs. p values for Microarray and qRT-PCR analysis are indicated in the figure.
Figure 4. qRT-PCR validation of differentially expressed circRNAs in PAs. qRT-PCR was used to validate the dysregulated expression of the four significantly upregulated and the two significantly downregulated circRNAs that showed the highest fold-expression in PAs (n = 16) vs. normal parathyroid tissue (n = 4) in the microarray analysis. The same samples used for microarray analysis were also used for qPCR validation. The y-axis represents the relative fold change compared to controls. Data are displayed as means ± SEMs. p values for Microarray and qRT-PCR analysis are indicated in the figure.
Cells 08 00015 g004
Figure 5. Stochastic approach of the potential interactions between genes and non-coding RNAs, in the pathogenesis of parathyroid adenomas. Some circRNAs play critical roles in fine-tuning the level of miRNA-mediated regulation of gene expression by sequestering the miRNAs. miR: microRNAs; circ: circular RNAs; EMT: epithelial–mesenchymal transition; RBPs: RNA-binding proteins, QKI: Quaking homolog, KH domain RNA binding (RNA binding protein).
Figure 5. Stochastic approach of the potential interactions between genes and non-coding RNAs, in the pathogenesis of parathyroid adenomas. Some circRNAs play critical roles in fine-tuning the level of miRNA-mediated regulation of gene expression by sequestering the miRNAs. miR: microRNAs; circ: circular RNAs; EMT: epithelial–mesenchymal transition; RBPs: RNA-binding proteins, QKI: Quaking homolog, KH domain RNA binding (RNA binding protein).
Cells 08 00015 g005
Table 1. Demographic and biochemical profile of included patients and healthy controls.
Table 1. Demographic and biochemical profile of included patients and healthy controls.
Patient CodeAgeBaseline Values *Post-PTX Values #
s CasPO4iPTHs CasPO4iPTH
NR: 2.10–2.50 mmol/LNR: 1.0–1.5 mmol/LNR: 1.4–6.8 pmol/LNR: 2.1–2.5 mmol/LNR: 1.0–1.5 mmol/LNR: 1.4–6.8 pmol/L
Di_M_27543.111.218.42.631.90.6
Di_M_28642.721.47.62.6521.2
Di_M_32702.710.918.72.91.51.5
Di_M_35482.851.2311.62.31.20.6
Di_M_36483.081.0912.92.311.431.7
Di_M_37782.631.118.51.881.671
Di_M_38562.681.411.52.231.82.2
Di_F_9653.271.512.72.631.390.6
Di_F_10652.751.412.12.311.440.6
Di_F_11582.740.916.22.451.234.2
Di_F_12582.691.210.32.21.33.5
Di_F_14672.681.129.51.931.30.6
Di_F_15672.481.1613.41.91.63.4
Di_F_20582.641.2210.82.341.482.8
Di_F_21572.751.18.62.141.651.4
Di_F_34522.631.0910.22.091.34.3
Con_1M652.231.44.3///
Con_2M712.411.453.8///
Con_3F652.191.324.4///
Con_4F622.291.255.2///
* Baseline values: were measured within two months prior to parathyroidectomy. # PostPTX values: were measured 1–2 months after parathyroidectomy. M, male; F, female; sCa, serum calcium; sPO4, serum phosphate; iPTH, intact parathyroid hormone; NR, normal range; PTX, parathyroidectomy.
Table 2. Circular RNAs with significantly altered expression in parathyroid adenomas compared to normal parathyroid.
Table 2. Circular RNAs with significantly altered expression in parathyroid adenomas compared to normal parathyroid.
CircRNACircRNA ID ‡Fold Changep ValueFDRChr.Gene/Symbolhsa-miRNA targets ‡‡
Up-regulated CircRNAs
hsa-CircRNA_102904hsa_Circ_00084594.60.0010.126KANSL1L/C2orf67(Chr.2)139-5p, 520f, 144, 1206, 146b-3p
hsa-CircRNA_404337hsa_Circ_01400942.140.0010.148RPS6KA3 (Chr. X)1248, 1261,1264, 331-3p, 623, 644
hsa-CircRNA_043614hsa_Circ_00436142.180.0460.482KRT14 (chr.17)1184, 1205, 338-3p, 361-3p, 1224-3p, 766
hsa-CircRNA_102903hsa_Circ_00026172.260.0060.23KANSL1L/C2orf67 (Chr.2)127-5p, 1290, 139-5p, 149, 409-3p
hsa-CircRNA_100965hsa_Circ_00245012.210.0020.153DDX6 (Chr.11)1248, 1251, 1264, 577, 607, 660
hsa-CircRNA_001379hsa_Circ_00005162.110.0070.234RPPH1 (Chr.14)1296
hsa-CircRNA_402533 **-5.580.0250.387OVOL2 (Chr.20)500b-3p
hsa-CircRNA_003949hsa_Circ_00039492.830.0060.23PTN (Chr.7)1256, 145, 326, 330-5p, 620
hsa-CircRNA_102116hsa_Circ_00032583.790.0260.392ZNF652 (Chr.17)502-5p, 203, 1278, 653, 766
hsa-CircRNA_404504hsa_Circ_01415862.430.0330.424STIL (Chr.1)-
hsa-CircRNA_051778hsa_Circ_00517785.340.0480.487BCAT2 (Chr.19)-
hsa-CircRNA_400537hsa_Circ_00934502.080.0050.223MASTL (Chr.10)1236, 1279, 1825,223, 346, 520
hsa-CircRNA_101283hsa_Circ_00305862.510.0430.465ABCC4 (Chr.13)1238, 1248, 217, 651, 767-5p, 139-5p
hsa-CircRNA_007273hsa_Circ_00072732.190.0010.146RAB11FIP5 (chr.2)1243, 1248, 1288, 361-3p, 876-3p
hsa-CircRNA_401782hsa_Circ_00431384.710.0460.479TAF15 (Chr.17)1257, 556-5p, 758, 1282, 767-3p
hsa-CircRNA_002122hsa_Circ_00021224.070.0250.391CTCF (Chr.16)1264, 136, 561, 644, 874
hsa-CircRNA_002617hsa_Circ_00026172.130.020.36KANSL1L/C2orf67 (Chr.2)127-5p, 1290, 139-5p, 149, 409,
hsa-CircRNA_406174 ***hsa_Circ_00056225.840.0360.436PITPNB (Chr.22)604, 589, 1248, 593, 1248
hsa-CircRNA_104269hsa_Circ_00787682.160.0060.234FAM120B (Chr.6)1292, 661,1182, 1272, 1322
hsa-CircRNA_004286hsa_Circ_00042862.530.0020.162SUSD1 (Chr.9)1299, 182,31, 330-3p, 384
hsa-CircRNA_406165-3.620.0450.477COL6A2 (Chr.21)-
hsa-CircRNA_008267hsa_Circ_00082676.360.0250.387LINC00969/SDHAP2(Chr.3)145, 1304, 198,498, 581
Down-regulated CircRNAs
hsa-CircRNA_404643hsa_Circ_01117462.310.010.272PIK3C2B (Chr.1)1224-3p, 1272, 1280, 486-3p, 586
hsa-CircRNA_105038hsa_Circ_01050382.840.0050.229FLNA (Chr. X)1225-3p, 1247, 149, 486-3p, 661
hsa-CircRNA_101474hsa_Circ_01014742.030.0020.153THBS1 (Chr.15)1322, 507, 557, 578, 615-5p, 663b
hsa-CircRNA_405040hsa_Circ_00995042.40.0010.153PLXNC1 Chr.12570, 1229, 513a-5p, 873, 1322
hsa-CircRNA_101051hsa_Circ_00261432.270.0410.458TROAP Chr.12637, 887, 1286, 661, 1247
hsa-CircRNA_032603hsa_Circ_00326033.010.0060.234LTBP2 Chr.14338-3p, 766, 940, 657, 486-3p
hsa-CircRNA_058097hsa_Circ_00580973.610.0030.198FN1 Chr.2876-3p, 1231, 1238, 515-5p, 558
hsa-CircRNA_074530hsa_Circ_00745302.080.0040.2CD74 Chr.5136, 145, 519e, 515-3p,767-3p
hsa-CircRNA_092437hsa_Circ_00007412.40.0440.472POLR2A Chr.17574-5p, 637, 604, 7, 490-5p
hsa-CircRNA_021732hsa_Circ_01117462.820.0020.153CD44 Chr.111248, 1236, 145, 370, 636
hsa-CircRNA_101017hsa_circ_00255062.010.0020.162GPRC5A Chr.121252, 296-5p,583, 873, 940
hsa-CircRNA_001766hsa_circ_00017662.090.0080.247PDIA4 Chr.7182, 1203, 136, 182, 766
hsa-CircRNA_102368hsa_circ_00476632.040.0430.467RPL17 Chr.181204, 1208, 1228, 586
‡: The experimentally assayed circRNAs were cross-checked for the genomic position and the best transcript in circBase (http://circbase.org). ‡‡: The circRNA -corresponding hsa-miRNA targets with the best context+score percentile (NIH Circular RNA interactome). All identified circularRNAs are exonic arising from the exons of the linear transcript except for: **: Antisense, where their gene locus overlaps with the linear RNA, but is transcribed from the opposite strand. ***: Sense overlapping, transcribed from same gene locus as the linear transcript, but not classified into "exonic" and "intronic". has: Homo sapiens; Circ: circular; chr: chromosome.
Table 3. Circular RNA expression profile in parathyroid adenomas from male vs. female patients.
Table 3. Circular RNA expression profile in parathyroid adenomas from male vs. female patients.
CircRNACircRNA ID‡Fold Changep ValueFDRGene Symbol/hsa-miRNA targets‡‡
Chr.
Up-regulated CircRNAs
hsa-CircRNA_401977hsa_circ_01087032.670.0010.09NEDD4L1231,1298, 149, 330-3p, 421, 561
Chr.18
hsa-CircRNA_103224hsa_circ_00633292.520.0010.07DDX17191, 571, 605, 622, 640
Chr.22
hsa-CircRNA_102850hsa_circ_00010812.150.0030.12DCAF17558,624,653,1200,1224-3p
Ch.2
hsa-CircRNA_008636 ***hsa_circ_00086362.230.0170.21AAGAB1256, 330-3p, 526b, 1182,1184
Chr.15
hsa-CircRNA_000763 **hsa_circ_00076302.640.0250.24MBNL11184, 1263, 1305, 192, 215
Chr.3
hsa-CircRNA_026646hsa_circ_00266462.280.0010.08PCBP21225-3p, 184, 1225-5p, 1233, 1245
Chr.12
hsa-CircRNA_406841 **-2.060.0050.13MIR5695-
Chr.19
hsa-CircRNA_092388 ***hsa_circ_00003153.260.0160.2INCENP767-3p, 942, 1184, 578, 622
Chr.11
hsa-CircRNA_089762hsa_circ_00897622.89<0.0010.04JA760602384, 554, 1249, 1283, 21
hsa-CircRNA_044837hsa_circ_00448372.970.0240.23YPEL21184, 1205, 1257, 1264, 1273
Chr.17
hsa-CircRNA_061170hsa_circ_00611702.170.0370.27RTEL11289, 338-3p, 370, 377, 409-3p
Chr.20
hsa-CircRNA_000711 *hsa_circ_00014432.140.0010.08MAML3873
Chr.4
hsa-CircRNA_001678 ***hsa_circ_00005172.390.0090.16RPPH11225-3p, 1233, 1258, 1292, 1296
Chr.14
hsa-CircRNA_051799hsa_circ_00517993.310.0010.07BAX1827, 198, 383, 486-3p, 1180
Chr.19
hsa-CircRNA_000365hsa_circ_00003652.070.0050.13TBCEL1233, 1304, 339-3p, 421, 510
Chr.11
hsa-CircRNA_404879 ***hsa_circ_00957842.220.0260.24PRR5L1179, 1197, 1205, 1208,1229
Chr.11
hsa-CircRNA_000167***hsa_circ_00005182.380.0110.17RPPH11204, 1225-3p,1233, 1258, 146b-3p
Chr.14
hsa-CircRNA_089763hsa_circ_00897632.740.0010.07JA7606001179, 1238, 1248, 1206, 1228
hsa-CircRNA_002082***hsa_circ_00020822.450.0450.29MALAT11252, 217, 1208, 1224-3p, 1238
Chr.11
Down-regulated CircRNAs
hsa-CircRNA_100332hsa_circ_0014130/0.0260.24PIP5K1A1205, 1246, 1182, 1183, 1184
Chr.1
hsa-CircRNA_004286hsa_circ_0042860/<0.0010.05CRLF31827, 192, 215, 346, 516b
Chr.17
hsa-CircRNA_059571hsa_circ_0059571/0.0090.17RALGAPA21250, 1263, 136, 1183, 1200
Chr.20
hsa-CircRNA_005310hsa_circ_0005310/0.0410.28SBNO11246, 1251, 1264, 1290, 375
Chr.12
All identified circularRNAs are exonic, arising from the exons of the linear transcript, except for *: intronic, arising from an intron of the linear transcript. **: Antisense, where their gene locus overlap with the linear RNA, but transcribed from the opposite strand. ***: Sense overlapping, transcribed from same gene locus as the linear transcript, but not classified into "exonic" and "intronic". hsa: Homo sapiens, Circ: circular, chr: chromosome.

Share and Cite

MDPI and ACS Style

Yavropoulou, M.P.; Poulios, C.; Michalopoulos, N.; Gatzou, A.; Chrisafi, S.; Mantalovas, S.; Papavramidis, T.; Daskalaki, E.; Sofou, E.; Kotsa, K.; et al. A Role for Circular Non-Coding RNAs in the Pathogenesis of Sporadic Parathyroid Adenomas and the Impact of Gender-Specific Epigenetic Regulation. Cells 2019, 8, 15. https://doi.org/10.3390/cells8010015

AMA Style

Yavropoulou MP, Poulios C, Michalopoulos N, Gatzou A, Chrisafi S, Mantalovas S, Papavramidis T, Daskalaki E, Sofou E, Kotsa K, et al. A Role for Circular Non-Coding RNAs in the Pathogenesis of Sporadic Parathyroid Adenomas and the Impact of Gender-Specific Epigenetic Regulation. Cells. 2019; 8(1):15. https://doi.org/10.3390/cells8010015

Chicago/Turabian Style

Yavropoulou, Maria P., Christos Poulios, Nickos Michalopoulos, Ariadni Gatzou, Sofia Chrisafi, Stylianos Mantalovas, Theodosis Papavramidis, Emily Daskalaki, Electra Sofou, Kalliopi Kotsa, and et al. 2019. "A Role for Circular Non-Coding RNAs in the Pathogenesis of Sporadic Parathyroid Adenomas and the Impact of Gender-Specific Epigenetic Regulation" Cells 8, no. 1: 15. https://doi.org/10.3390/cells8010015

APA Style

Yavropoulou, M. P., Poulios, C., Michalopoulos, N., Gatzou, A., Chrisafi, S., Mantalovas, S., Papavramidis, T., Daskalaki, E., Sofou, E., Kotsa, K., Kesisoglou, I., Zebekakis, P., & Yovos, J. G. (2019). A Role for Circular Non-Coding RNAs in the Pathogenesis of Sporadic Parathyroid Adenomas and the Impact of Gender-Specific Epigenetic Regulation. Cells, 8(1), 15. https://doi.org/10.3390/cells8010015

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop