A Comprehensive Analysis Identified the Key Differentially Expressed Circular Ribonucleic Acids and Methylation-Related Function in Pheochromocytomas and Paragangliomas

We investigated differentially expressed circular RNAs (circRNAs) and their potential functions in pheochromocytomas and paragangliomas (PCC/PGLs). Expression levels of circRNAs in tumor and adjacent normal tissues from seven PCC/PGL patients were analyzed through RNA sequencing. Real-time PCR was conducted to verify the key candidates identified in the sequencing data. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to predict the functions of these circRNAs. A total of 367 circRNAs were found differentially expressed between tumor and normal samples. The top three histone methylation-related circRNAs (hsa_circ_0000567, hsa_circ_0002897, and hsa_circ_0004473) and their target microRNAs (miRNAs) were identified and validated. We then mapped the circRNA-miRNA-messenger RNA (mRNA) coding-noncoding gene co-expression (CNC) networks to show the potential binding relationships between circRNAs and their targets in PCC/PGLs. The top five mRNAs, 88 miRNAs, and 132 circRNAs related to pathogenesis were utilized to map the CNC network, and we observed that the interactions of these candidates with their target miRNAs regulated histone methylation and further mediated PCC/PGL pathogenesis. This study is the first to provide the whole profile of differentially expressed circRNAs in PCC/PGLs. Our data indicate that altered circRNAs may control the pathogenesis of PCC/PGLs by regulating histone methylation processes, highlighting their role as potential biomarkers.

Introduction

Pheochromocytomas and paragangliomas (PCC/PGLs) are highly genetically related, neuroendocrine tumors that have been listed as rare diseases by the World Health Organization (Neumann et al., 2019). Adrenal chromaffin cells secrete catecholamines and give rise to PCC/PGLs (Neumann et al., 2007; Neary et al., 2011; Gill et al., 2011). Approximately 90% of chromaffin cells are present in the adrenal medulla and tumors that occur here are termed PCC (Ferreira et al., 2016). The other 10% of cells are located outside the adrenal medulla, such as in the heart, bladder, and neck, and these tumor cases are called extra-adrenal PCC or PGL (Crona et al., 2017). The major clinical manifestations of PCC/PGLs are unpredictable hypertension, cardiovascular crisis, hyperhidrosis, palpitations, and elevated basal metabolic rate, which are induced by the excessive secretion of catecholamines (Jochmanova and Pacak, 2016; Loper et al., 2016). PCC/PGLs are also involved in some severe syndromes such as von Hippel-Lindau disease (Liu P. et al., 2019), multiple endocrine neoplasia (type I and II), and neurofibromatosis (Rednam et al., 2017). If patients fail to undergo proper treatment in time, their 5-year survival rate drops below 40% due to excessive catecholamine secretion as well as the development of chemotherapy and radiation therapy resistance (Toledo et al., 2016; Taieb and Pacak, 2017; Li et al., 2018). As PCC/PGLs are among the most genetically related human cancers with 60% of patients showing family aggregation (Fishbein et al., 2017; Liu et al., 2018; Job et al., 2019), genetic screening is of great value and significance for the diagnosis and prevention of PCC/PGLs (Calsina et al., 2019). Moreover, over 15 genomic and transcriptomic molecules reported to regulate PCC/PGL development, such as SDHx, VHL, TMEM127, HRAS, FGFR1, ATRX, RET, EPAS1, MAX, EGLN1, and NF1, were found mutated in the germline of PCC/PGLs patients (Gill et al., 2011; Castro-Vega et al., 2015; Rednam et al., 2017; Taieb and Pacak, 2017; Remacha et al., 2017; Pang et al., 2018).

Circular RNAs (circRNAs) are highly evolutionarily conserved, regulatory, noncoding RNAs that play significant roles in various biological and pathological processes, such as gene translation, transcription, and tumor genesis and development; they also act as microRNA (miRNA) sponges. CircRNAs have a stable loop structure that lacks a 5’ cap and 3’ tail and are thus not easily degraded by RNase. Under specific circumstances, some circRNAs can open their loop structure for translation (Memczak et al., 2013; Liu C.X. et al., 2019). Moreover, circRNAs exhibit abnormal expression in numerous cancer types, including liver, stomach, breast, ovarian, and oral cancers; therefore, they have been considered novel biomarkers for tumor screening and are potential therapeutic targets (Hansen et al., 2013; Conn et al., 2015; Meng et al., 2019; Vo et al., 2019).

Considering the high genetic relevance and neuroendocrine characteristics of PCC/PGLs, we speculated that circRNAs may play important roles in regulating the development of PCC/PGLs. Previous studies have identified multiple key mediators for the pathogenesis of PCC/PGLs (Rednam et al., 2017; Remacha et al., 2017; Job et al., 2019), such as genetic alterations of SHDx and VHL, but whether these components are modulated by circRNAs remains unknown. Therefore, we performed RNA sequencing using tumor tissues and para-cancerous normal adrenal medulla tissues from PCC/PGLs patients to examine the aberrant expression of circRNAs and miRNAs. The data were evaluated by comparative analyses of reference sequences and tissue specificity was predicted. circRNAs with significantly altered expression were further screened and their binding sites to miRNAs predicted. In addition, we performed gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses using these differentially expressed circRNAs and mapped their interaction networks.

Materials and Methods

Sample Collection

Tumor and adjacent normal adrenal medulla specimens were collected during surgery from 40 PCC/PGL patients who underwent laparoscopic resection between February 2018 and September 2018 in Xianya Hospital. Detailed patient information is listed in Table 1. Seven PCC/PGL patients were randomly selected for RNA sequencing of their tissues while the other 33 PCC/PGL patients were used to validate the sequencing outcomes. All samples were validated by pathological examination and were frozen and stored immediately in liquid nitrogen until further use. Peripheral blood of 16 PCC/PGL patients and 16 healthy donors was collected and stored immediately at −80°C. All patients were primarily diagnosed with PCC/PGLs via imaging (computerized tomography and color ultrasound) and laboratory examination (urine vanillylmandelic acid), and did not undergo any treatments before surgery. All patients signed informed consent forms before the surgery and study, which was approved by the Ethics Committee of Xiangya Hospital of Central South University.

TABLE 1

Table 1 Characteristics of the seven pheochromocytomas and paragangliomas (PCC/PGL) patients whose samples were used for RNA sequencing.

Total Ribonucleic Acid Extraction

Frozen samples and peripheral blood were taken out of liquid nitrogen storage and immediately mixed with TRIzol reagent (Ambion, Austin, TX). Tissues were homogenized and lysed in TRIzol reagent, after which total RNAs were extracted according to manufacturer’s instructions.

Ribonucleic Acid Quantification and Qualification

Electrophoresis on 2% agarose gels was used to determine whether the total RNAs were degraded or contaminated. A NanoPhotometer^® spectrophotometer (IMPLEN, Westlake Village, CA) and Qubit^® RNA Assay Kit in Qubit^®2.0 Fluorometer (Life Technologies, Carlsbad, CA) were used to detect the purity and concentration of RNA, respectively. The RNA Nano6000 Assay Kit (Agilent Technologies, Santa Clara, CA) was used to assess RNA integrity.

Ribonucleic Acid Library Construction, Clustering, and Sequencing

We used approximately 5 ug of total RNA for each sample and used the Ribo-Zero™ rRNA Removal Kit (Illumina, San Diego, CA) to eliminate ribosomal RNA interference according to manufacturer’s instructions. The residue was purified by two rounds of ethanol precipitation, after which divalent cations were used to splice the remaining RNAs into small fragments under high temperature. The small RNA fragments were reverse-transcribed to complementary DNA (cDNA) and second strand cDNA was synthesized using Escherichia coli DNA polymerase I, deoxyuridine triphosphates (dUTPs), and RNase H. After adding an A-base to the 3ʹ ends of each cDNA fragment and ligating to the NEBNext Adaptor [New England Biolabs (NEB), Ipswich, MA], we purified the library using the AMPureXP system (Beckman Coulter, Brea, CA) and prepared for hybridization. Next, approximately 2 µl USER enzyme (NEB) was added to the cDNA buffer for PCR. TruSeq SR Cluster Kit v3-cBot-HS (Illumia) was used to cluster the samples on the cBot Cluster Generation System. Finally, the Illumina HiSeq 2500/2000 sequencing platform was used to sequence the libraries (150 bp paired-end reads) after cluster generation. Sequencing data were analyzed against the reference genome, after which the circRNAs were filtered out. The software Find_circ and CIRI2 (Zeng et al., 2017) were used to identify circRNAs;

Identification of Differentially Expressed Circular Ribonucleic Acids, Microribonucleic Acids, and Messenger Ribonucleic Acids

DESeq R package(1.8.3) (https://www.bioconductor.org) was used to perform read-count analysis of differentially expressed RNAs [circRNAs, miRNAs, and messenger RNAs (mRNAs)]. An adjusted p value of < 0.05 was considered statistically significant and related RNAs were recognized as differentially expressed. miRanda software (https://omictools.com › miranda-tool) was used to predict the target miRNA binding sites of circRNAs.

Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Pathway Enrichment Analysis

Gene Ontology (GO) enrichment analysis (http://www.webgestalt.org), an online software, was used to predict the molecular function, biological process, cellular components, and location of the target miRNAs of differentially expressed circRNAs. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis (http://kobas.cbi.pku.edu.cn/anno_iden.php) was used to predict the possible metabolic pathways of some differentially expressed circRNAs Enrichment with p < 0.05 was regarded as statistically significant.

Real Time Polymerase Chain Reaction Analysis

Reverse transcription of total RNA was performed using the Superscript IV Reverse Transcriptase (Thermo Fisher Scientific, Waltham, MA) at 50°C for 15 min to obtain cDNA. Real-time PCR was then performed using SYBR Green Reagent on an ABI PRISM 7500 Fast Real-Time PCR System (Thermo Fisher Scientific, Waltham, MA). Primers used for qPCR were synthesized at Sango Biotech (Shanghai, China) and are listed in Supplementary Table 6. Relative expression of the different genes was normalized to β-actin expression, and arbitrary units were used to display normalized gene expression. The data were analyzed using the 2^−ΔΔCt method.

Circular Ribonucleic Acid-Microribonucleic Acid-Messenger Ribonucleic Acid Coding-Noncoding Gene Co-Expression Network Analysis

Cytoscape software (Shannon et al., 2003) was used to construct the circRNA-miRNA-mRNA networks.

Statistical Analysis

All data were analyzed using SPSS version 23 (IBM, Armonk, NY) and GraphPad Prism 5 (GraphPad Software Inc., La Jolla, CA). A paired Student’s t-test was used to test significance of the results and p values < 0.05 were considered statistically significant.

Results

Circular Ribonucleic Acid Expression Profiles of Human Pheochromocytomas and Paragangliomas Tumor Tissues

The density distribution of circRNAs on the chromosome is shown in Figure 1A. (Chromosome 1–9 and X were shown in the panel, details of circRNAs on all the chromosome can be seen in Supplementary Table 1 the density distribution of circRNAs on all the chromosome). Different sources (exons, introns, and intergenic regions) of circRNAs were counted (Figure 1B) because circRNAs can be acquired from the splicing of exons and introns. We thus speculated that most circRNAs in PCC/PGLs are obtained from the splicing of exons in both tumor and adjacent normal tissues. Next, we analyzed the length distribution of circRNAs in different patient samples (Figure 1C) and found that the majority of circRNAs were between 0 and 2,000 nt. Additionally, we mapped the transcription per million (TPM) distribution graph for an overall examination of the gene expression pattern in different samples (Figure 1D).

FIGURE 1

Figure 1 Profiles of circular RNAs (circRNAs) identified in pheochromocytomas and paragangliomas (PCC/PGLs). (A) Density distribution of circRNAs on the chromosome. (B) The different sources of circRNAs (exons, introns, and intergenic regions). (C) Length distribution graph of circRNAs in different patient samples. (D) Transcription per million (TPM) distribution graph. (RNASEQC: tumor samples of PCC/PGLs patients for RNA sequencing, RNASEQN: normal adrenal medulla tissues of PCC/PGLs patients for RNA sequencing).

Identification and Validation of Differentially Expressed Circular Ribonucleic Acid in Pheochromocytomas and Paragangliomas

A total of 367 differentially expressed circRNAs were identified after RNA sequencing and qualification, among which 255 were upregulated and 112 were downregulated. The top 20 differentially expressed circRNAs are listed in Table 2. We next performed hierarchical cluster analysis of differentially expressed circRNAs and acquired the volcano map (log2FoldChange > 1; adjusted p < 0.05) and heat map for these samples (Figures 2A, B). We selected the top 11 differentially expressed circRNAs and conducted real-time qPCR using the other 33 pairs of PCC/PGLs samples to test if these circRNAs show the same relative expression as in the RNA sequencing data; information and genotype of the 33 patients is available in Supplementary Table 2 information and genotype of the 33 patients and the primer sequences for each circRNA is listed in Supplementary Table 6. As shown in Figure 2C–E, expression levels of hsa_circ_0000567, hsa_circ_0002897, hsa_circ_0004473, hsa_circ_0000972, hsa_circ_0000825, hsa_circ_0003265, and hsa_circ_0007279 were significantly higher in tumor tissues than in para-cancerous tissues, whereas the expression of hsa_circ_0056892, hsa_circ_0019773, and hsa_circ_0007444 was markedly downregulated in tumors tissues compared with adjacent tissues. However, there was no difference in hsa_circ_0001573 expression levels via real-time PCR. Overall, these data indicate that the sequencing outcomes were reliable and that these altered circRNAs may play important roles in tumor genesis or pathological processes and can serve as potential biomarkers or therapeutic targets.

TABLE 2

Table 2 The top 20 differentially expressed circular ribonucleic acids (circRNAs) in pheochromocytomas and paragangliomas.

FIGURE 2

Figure 2 Identification and validation of differentially expressed circular RNAs (circRNAs) in pheochromocytomas and paragangliomas (PCC/PGLs). (A) Heatmap of differentially expressed circRNAs. (RNASEQC, tumor samples of PCC/PGLs patients for RNA sequencing; RNASEQN, normal adrenal medulla tissues of PCC/PGLs patients for RNA sequencing). (B) Volcano map of differentially expressed circRNAs; 255 were upregulated and 122 were downregulated. (PHEOC, tumor tissues of PCC/PGLs patients; PHEON, normal adrenal medulla tissues of PCC/PGLs patients.) (C–E) Real time PCR validation of the top 11 differentially expressed circRNAs in 33 pairs of PCC/PGL tumor tissues and adjacent tissues.

Gene Ontology Enrichment and Analysis of Genes Associated With Differentially Expressed Circular Ribonucleic Acids

After obtaining the differentially expressed circRNAs, we performed GO analysis on the source genes in each group based on the correspondence between circRNAs and their binding miRNAs and source genes (Figure 3A). We found that the main biological process terms were related to cell metabolism and transport, such as regulation of cholesterol metabolic process (GO:0090181), positive regulation of ion transport (GO:0043270), regulation of transport (GO:0051049), and regulation of hydrolase activity (GO:0051336). The highest enriched cellular component terms were associated with microtubule (GO:0005874) and microtubule cytoskeleton (GO:0015630). The most enriched molecular function terms were related to histone methyltransferase activity (H3-K36 specific; GO:0046975), transfer RNA (tRNA) (cytosine-5-)-methyltransferase activity (GO:0016428), and tRNA (cytosine) methyltransferase activity (GO:0016427).

GO enrichment results showed that the degree of enrichment of histone methylation-related circRNAs occupied the top positions of molecular function, suggesting that methylation regulation of target genes by circRNAs may play important roles in the pathogenesis of PCC/PGLs. We further performed real-time PCR to analyze the top three histone methylation-related circRNAs in the peripheral blood of 16 PCC/PGL patients (see Supplementary Table 3 information and genotype of the 16 patients) and 16 healthy donors, and found that circRNAs has_circ_0000567, hsa_circ_0002897, ahashsa_circ_0004473 also exhibited relatively high expression levels in PCC/PGL patients compared with healthy donors (Figure 3B).

FIGURE 3

Figure 3 Gene ontology prediction and validation of histone methylation-related circular RNAs (circRNAs) in peripheral blood. (A) Gene ontology analysis of the top enriched molecular function of the differently expressed circRNAs. (B) Validation of the top three differentially expressed circRNAs by real time PCR using the peripheral blood of PCC/PGL patients and healthy donors.

Prediction and Analysis of the Histone Methylation-Related Circular Ribonucleic Acid-Microribonucleic Acid-Messenger Ribonucleic Acid Coding-Noncoding Gene Co-Expression Network

Next, we selected the three histone methylation-related circRNAs (hsa_circ_0000567, hsa_circ_0002897, and hsa_circ_0004473) with the most significant differences to predict and analyze their target miRNAs and mRNAs. The top five target miRNAs and related mRNAs were used to map the coding-noncoding gene co-expression (CNC) network (Figure 4A). Notably, the predicted binding miRNAs of the three upregulated circRNAs in PCC/PGL tumor tissues were all downregulated according to the miRNA sequencing data (Figure 4B and Supplementary Table 4 the circRNAs-miRNAs binding relationship). Our findings suggest that these histone methylation-related circRNAs act as RNA sponges and regulate their target miRNAs, highlighting their potential as novel therapeutic targets for PCC/PGLs.

FIGURE 4

Figure 4 Prediction of histone methylation-related messenger RNAs (miRNAs). (A) The coding-noncoding gene co-expression (CNC) prediction network of the top three methylation-related circRNAs (hsa_circ_0000567, hsa_circ_0002897, and hsa_circ_0004473) and their target miRNAs and mRNAs (top five miRNAs and mRNAs are shown on the map). (B) Heatmap of differentially expressed miRNAs by RNA sequencing. (RNASEQC, tumor samples of PCC/PGLs patients for RNA sequencing; RNASEQN, normal adrenal medulla tissues of PCC/PGLs patients for RNA sequencing).

Kyoto Encyclopedia of Genes and Genomes Pathway Enrichment and Analysis of Differentially Expressed Circular Ribonucleic Acids

KEGG pathway analysis revealed that the most enriched pathways were related to insulin secretion, axon guidance, endocytosis, bile secretion, glutamatergic synapse, serotonergic synapse, lysine degradation, and steroid biosynthesis. The top 20 enriched pathways were selected and displayed in a KEGG enrichment scatter plot diagram (Figure 5).

FIGURE 5

Figure 5 Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the differently expressed circular RNAs (circRNAs) and their target gene-related enriched pathways.

Prediction and Analysis of the Relationship Between Known Pheochromocytoma and Paraganglioma Susceptibility Genes and the Circular Ribonucleic Acid-Microribonucleic Acid-Messenger Ribonucleic Acid Network

The currently known PCC/PGL-related susceptibility genes include SDHx, VHL, TMEM127, HRAS, FGFR1, ATRX, RET, EPAS1, MAX, EGLN1, and NF1 (Castro-Vega et al., 2015). We thus examined whether the circRNAs exhibit regulatory relationships with the mRNA transcripts of these susceptibility genes. After analyzing and screening the mRNA sequencing data (Supplementary Table 5 mRNA sequencing data of the 7 patients), we selected five differently expressed mRNAs, 88 miRNAs, and 132 circRNAs to map the CNC network. These findings support the role of circRNAs as RNA sponges and in regulating the expression levels of PCC/PGL susceptibility genes (Figure 6).

FIGURE 6

Figure 6 The competing endogenous RNA network of the known pheochromocytoma and paraganglioma (PCC/PGL) susceptibility genes associated with the circular ribonucleic acid-microribonucleic acid-messenger ribonucleic acid (circRNA-miRNA-mRNA) relationship. Five mRNAs (red), 88 miRNAs (yellow), and 132 circRNAs (green) were included in the analysis.

Discussion

CircRNA is a special type of highly evolutionarily conserved, regulatory, noncoding RNA, which was originally considered noise or a redundancy of the transcription process (Memczak et al., 2013). However, an increasing number of studies demonstrated that circRNAs play important roles in various pathological or physiological processes (Zeng et al., 2017; Meng et al., 2019; Kolling et al., 2019; Kristensen et al., 2019; Li et al., 2019; Liang et al., 2019). Studies showed that circRNAs are mainly produced by RNA cleavage, are widely present throughout eukaryotic cells, and have high stability, species conservation, and tissue specificity. CircRNAs are also involved in intracellular RNA regulation networks and are closely associated with the development of many diseases, including systemic sclerosis and several cancer types, such as liver, lung, stomach, oral, nasopharyngeal, ovarian, prostate, uterine, kidney, and bladder cancers (Liu H. et al., 2019; Wang et al., 2019; Wu et al., 2019; Yao et al., 2019; Yu et al., 2019; Zhang et al., 2019; Zhi et al., 2019). Moreover, it has been reported that circRNAs act as RNA sponges, binding to their corresponding miRNAs to regulate the expression of target genes (Memczak et al., 2013). Nevertheless, the aberrant expression of circRNAs on a transcriptional level as well as their regulatory functions have not been explored in PCC/PGL patients. Therefore, in this study, we conducted a comprehensive analysis of differentially expressed circRNAs and their potential epigenetic role in PCC/PGLs, subsequently constructing a circRNA-miRNA CNC network that shows the systematic regulatory function of circRNAs in the pathogenesis of PCC/PGLs.

PCC/PGLs are highly genetically related tumors, and thus research on the regulatory function of its transcription processes may further our understanding of PCC/PGL development and better contribute to the diagnosis and treatment of the disease (Neary et al., 2011; Liu et al., 2015; Zhu et al., 2019; Liu Y. et al., 2019). In the last decades, numerous studies have screened the genetic alterations of PCC/PGLs in an attempt to elucidate the major mutant genes and the disease subtypes caused by those genes (Neumann et al., 2007; Neary et al., 2011; Gill et al., 2011; Taieb and Pacak, 2017; Yi et al., 2018; Abrard et al., 2019). As PCC/PGLs are closely related to epigenetic regulation, comprehensive regulatory network studies are more conducive for explaining their pathogenesis than are studies conducted on a few molecules. Thus, we performed transcriptome analysis of seven PCC/PGL patients via RNA sequencing and subsequently identified 3927 mRNAs, 283 miRNAs, and 367 circRNAs that were differentially expressed. We then selected the top 11 differently expressed circRNAs for real-time PCR validation in 33 pairs of tumor tissues and adjacent normal tissues of PCC/PGL patients, and found that the cross-regulated target genes were consistent with the sequencing results.

Some studies have shown that circRNAs can regulate DNA methylation in certain diseases. Wang et al. (Morley-Smith and Lyon, 2016) reported that CircIBTK can inhibit DNA demethylation and inactivate the AKT pathway in systemic lupus erythematosus patients. Our GO and KEGG pathway analysis of the sequencing results also showed that the major differentially expressed circRNAs were enriched in methylation-related molecular functions such as histone methyltransferase activity and tRNA (cytosine-5-)-methyltransferase activity, which suggests that these circRNAs may affect the development and progression of PCC/PGLs by participating in the modification of target gene methylation. Further analysis of our miRNA sequencing data showed that the histone methylation-related miRNAs predicted by the circRNA-miRNA CNC network were indeed downregulated in tumor tissues, supporting the notion that these circRNAs and miRNAs interact with each other and regulate histone methylation. Moreover, in the KEGG pathway-enriched data, genes with higher enrichment were mainly related to insulin secretion and axon guidance, which coincides with the pathological mechanism of PCC/PGLs patients, wherein neuroendocrine tumors secrete a large amount of catecholamines to disrupt the endocrine and metabolism balance, leading to a series of clinical symptoms such as hypertension and obesity (Neary et al., 2011). Therefore, these differentially expressed circRNAs may affect the secretion and metabolism processes of PCC/PGL patients by modulating their target genes. Our findings thus provide new strategies and targets for the diagnosis and treatment of such patients.

It is well known that circRNAs are stable, specific, and abundant (Wang et al., 2018; Yu et al., 2019; Zhi et al., 2019), making them suitable for detection as circulating biomarkers in liquid biopsies. Indeed, we identified some differentially expressed circRNAs in PCC/PGL patients, such as hsa_circ_0000567, hsa_circ_0002897, and hsa_circ_0004473, that bind to miRNA and affect histone methylation; these three circRNAs were also found differentially expressed between the peripheral blood of PCC/PGL patients and healthy donors. Previous studies have shown that the hsa-miR-183/182/96 cluster, hsa-miR-21-3p, hsa-miR-551b-3p, and hsa-miR-202-5p are associated with metastatic risk and progression of PCC/PGLs (Castro-Vega et al., 2015; Calsina et al., 2019). Notably, our study showed that hsa_circ_0000567 may bind to hsa-miR-96-3p and regulate histone methylation, highlighting hsa_circ_0000567 as a diagnostic and prognostic biomarker for PCC/PGL patients. Meanwhile, SDHx, MEN1, RET, and HF susceptibility genes were identified as source genes according to our mRNA sequencing data, after which the circRNA-miRNA-mRNA CNC network of PCC/PGL patients was mapped, providing a basis as well as novel ideas for future PCC/PGL research and diagnostic biomarker discovery.

In conclusion, the present study demonstrated that differentially expressed circRNAs such as hsa_circ_0000567, hsa_circ_0002897, and hsa_circ_0004473 can interact with their target miRNAs and regulate histone methylation during the pathogenesis of PCC/PGLs. KEGG pathway analysis revealed that the differently expressed circRNAs were involved in endocrine-related pathways, consistent with the clinical manifestations of PCC/PGLs. In addition, we identified the circRNAs associated with PCC/PGL susceptibility genes and mapped their regulatory network, which may provide potential prognostic biomarkers and therapeutic targets for PCC/PGL patients.

Data Availability Statement

The accession number for the data deposited in the GEO database is GSE139518.

Ethics Statement

All patients signed the informed consent forms before the surgery and study, which was approved by the Ethics Committee of Xiangya Hospital of Central South University.

Author Contributions

AY and ML designed the study and conducted the experiment and wrote down the manuscript, DC, CW, and QX collected the samples, CX help analyze the result and gave writing suggestions. LZ, YW, XZ, and YP help review the manuscript. LL guided the study. All the authors approved the final version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (grant number 81400773), Innovation-Driven Project of Central South University (No.020CX046), and Funds for the Shenghua Yuying Talents Program of Central South University.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors would like to thank Novogene company (Beijing, China) for the RNA sequencing.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00015/full#supplementary-material

Supplementary Table 1 | CircRNAs on all the chromosome.

Supplementary Table 2 | The information of the 33 patients.

Supplementary Table 3 | The information of the 16 PCC peripheral blood.

Supplementary Table 4 | PHEOCvsPHEON.differentially expressed miRNA.

Supplementary Table 5 | PHEOCvsPHEON.differentially expressed Mrna.

Supplementary Table 6 | Primer sequences for identification of circRNA by qRT-PCR.

References

Abrard, S., Parfait, E., Fouquet, O., Rineau, E. (2019). Intra-aortic counterpulsation therapy allowing the diagnosis of a pheochromocytoma. Am. J. Respir. Crit. Care Med. 200, e89–e90. doi: 10.1164/rccm.201903-0617IM

PubMed Abstract | CrossRef Full Text | Google Scholar

Calsina, B., Castro-Vega, L. J., Torres-Perez, R., Inglada-Perez, L., Curras-Freixes, M., Roldan-Romero, J. M., et al. (2019). Integrative multi-omics analysis identifies a prognostic miRNA signature and a targetable miR-21-3p/TSC2/mTOR axis in metastatic pheochromocytoma/paraganglioma. Theranostics. 9, 4946–4958. doi: 10.7150/thno.35458

PubMed Abstract | CrossRef Full Text | Google Scholar

Castro-Vega, L. J., Letouze, E., Burnichon, N., Buffet, A., Disderot, P. H., Khalifa, E., et al. (2015). Multi-omics analysis defines core genomic alterations in pheochromocytomas and paragangliomas. Nat. Commun. 6, 6044. doi: 10.1038/ncomms7044

PubMed Abstract | CrossRef Full Text | Google Scholar

Conn, S. J., Pillman, K. A., Toubia, J., Conn, V. M., Salmanidis, M., Phillips, C. A., et al. (2015). The RNA binding protein quaking regulates formation of circRNAs. Cell 160, 1125–1134. doi: 10.1016/j.cell.2015.02.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Crona, J., Taieb, D., Pacak, K. (2017). New perspectives on pheochromocytoma and paraganglioma: toward a molecular classification. Endocr. Rev. 38, 489–515. doi: 10.1210/er.2017-00062

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferreira, V. M., Marcelino, M., Piechnik, S. K., Marini, C., Karamitsos, T. D., Ntusi, N. A. B., et al. (2016). Pheochromocytoma is characterized by catecholamine-mediated myocarditis, focal and diffuse myocardial fibrosis, and myocardial dysfunction. J. Am. Coll. Cardiol. 67, 2364–2374. doi: 10.1016/j.jacc.2016.03.543

PubMed Abstract | CrossRef Full Text | Google Scholar

Fishbein, L., Leshchiner, I., Walter, V., Danilova, L., Robertson, A. G., Johnson, A. R., et al. (2017). Comprehensive molecular characterization of pheochromocytoma and paraganglioma. Cancer Cell 31, 181–193. doi: 10.1016/j.ccell.2017.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Gill, A. J., Pachter, N. S., Clarkson, A., Tucker, K. M., Winship, I. M., Benn, D. E., et al. (2011). Renal tumors and hereditary pheochromocytoma-paraganglioma syndrome type 4. N. Engl. J. Med. 364, 885–886. doi: 10.1056/NEJMc1012357

PubMed Abstract | CrossRef Full Text | Google Scholar

Hansen, T. B., Jensen, T. I., Clausen, B. H., Bramsen, J. B., Finsen, B., Damgaard, C. K., et al. (2013). Natural RNA circles function as efficient microRNA sponges. Nature. 495, 384–388. doi: 10.1038/nature11993

PubMed Abstract | CrossRef Full Text | Google Scholar

Job, S., Draskovic, I., Burnichon, N., Buffet, A., Cros, J., Lepine, C., et al. (2019). Telomerase activation and ATRX mutations are independent risk factors for metastatic pheochromocytoma and paraganglioma. Clin. Cancer Res. 25, 760–770. doi: 10.1158/1078-0432.CCR-18-0139

PubMed Abstract | CrossRef Full Text | Google Scholar

Jochmanova, I., Pacak, K. (2016). Pheochromocytoma: the first metabolic endocrine cancer. Clin. Cancer Res. 22, 5001–5011. doi: 10.1158/1078-0432.CCR-16-0606

PubMed Abstract | CrossRef Full Text | Google Scholar

Kolling, M., Haddad, G., Wegmann, U., Kistler, A., Bosakova, A., Seeger, H., et al. (2019). Circular RNAs in urine of kidney transplant patients with acute T cell-mediated allograft rejection. Clin. Chem. 65, 1287–1294. doi: 10.1373/clinchem.2019.305854

PubMed Abstract | CrossRef Full Text | Google Scholar

Kristensen, L. S., Andersen, M. S., Stagsted, L. V. W., Ebbesen, K. K., Hansen, T. B., Kjems, J. (2019). The biogenesis, biology and characterization of circular RNAs. Nat. Rev. Genet. 20, 675–691. doi: 10.1038/s41576-019-0158-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M., Liu, Q., Liu, P., Yi, X., Guan, X., Yu, A., et al. (2018). Long-term follow-up and novel splice donor mutation in MEN1 in a Chinese family. Oncotarget. 9, 1577–1586. doi: 10.18632/oncotarget.23100

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Xu, J. D., Fang, X. H., Zhu, J. N., Yang, J., Pan, R., et al. (2019). Circular RNA circRNA_000203 aggravates cardiac hypertrophy via suppressing miR26b-5p and miR-140-3p binding to Gata4. Cardiovasc. Res. doi: 10.1093/cvr/cvz215

CrossRef Full Text | Google Scholar

Liang, Y., Song, X., Li, Y., Ma, T., Su, P., Guo, R., et al. (2019). Targeting the circBMPR2/miR-553/USP4 axis as a potent therapeutic approach for breast cancer. Mol. Ther. Nucleic Acids 17, 347–361. doi: 10.1016/j.omtn.2019.05.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H., Li, W. Z., Wang, X. Y., Pei, Y. G., Long, X. Y., Chen, C. Y., et al. (2015). A rare case of extra-adrenal pheochromocytoma localized to the ovary and detected via abdominal computed tomography angiography. Oncol. Lett. 9, 774–776. doi: 10.3892/ol.2014.2718

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, P., Li, M., Guan, X., Yu, A., Xiao, Q., Wang, C., et al. (2018). Clinical syndromes and genetic screening strategies of pheochromocytoma and paraganglioma. J. Kidney Cancer VHL. 5, 14–22. doi: 10.15586/jkcvhl.2018.113

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, P., Zhu, F., Li, M., Dube, D. A., Liu, Q., Wang, C., et al. (2019). Von Hippel-Lindau “Black Forest” mutation inherited in a large Chinese family. Gland Surg. 8, 343–353. doi: 10.21037/gs.2019.08.03

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C. X., Li, X., Nan, F., Jiang, S., Gao, X., Guo, S. K., et al. (2019). Structure and degradation of circular RNAs regulate PKR activation in innate immunity. Cell. 177, 865–880.e821. doi: 10.1016/j.cell.2019.03.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H., Liu, Y., Bian, Z., Zhang, J., Zhang, R., Chen, X., et al. (2019). Correction to: circular RNA YAP1 inhibits the proliferation and invasion of gastric cancer cells by regulating the miR-367-5p/p27 Kip1 axis. Mol. Cancer 18, 117. doi: 10.1186/s12943-019-1045-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Liu, L., Zhu, F. (2019). Therapies targeting the signal pathways of pheochromocytoma and paraganglioma. Onco. Targets Ther. 12, 7227–7241. doi: 10.2147/OTT.S219056

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, P., Li, X., Guo, X., Chen, J., Li, C., Chen, M., et al. (2019). Circular RNA DOCK1 promotes bladder carcinoma progression via modulating circDOCK1/hsa-miR-132-3p/Sox5 signalling pathway. Cell Prolif. 52, e12614. doi: 10.1111/cpr.12614

PubMed Abstract | CrossRef Full Text | Google Scholar

Loper, H., Kant, R., Olson, J., Streeten, E., Munir, K. (2016). Invasive phaeochromocytoma presenting as a right atrial mass. Lancet Diabetes Endocrinol. 4, 286. doi: 10.1016/S2213-8587(15)00466-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Memczak, S., Jens, M., Elefsinioti, A., Torti, F., Krueger, J., Rybak, A., et al. (2013). Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 495, 333–338. doi: 10.1038/nature11928

PubMed Abstract | CrossRef Full Text | Google Scholar

Meng, S., Zhou, H., Feng, Z., Xu, Z., Tang, Y., Wu, M. (2019). Epigenetics in neurodevelopment: emerging role of circular RNA. Front. Cell Neurosci. 13, 327. doi: 10.3389/fncel.2019.00327

PubMed Abstract | CrossRef Full Text | Google Scholar

Morley-Smith, A. C., Lyon, A. R. (2016). Stressing the importance of cardiac assessment in pheochromocytoma. J. Am. Coll. Cardiol. 67, 2375–2377. doi: 10.1016/j.jacc.2016.03.510

PubMed Abstract | CrossRef Full Text | Google Scholar

Neary, N. M., King, K. S., Pacak, K. (2011). Drugs and pheochromocytoma–don’t be fooled by every elevated metanephrine. N. Engl. J. Med. 364, 2268–2270. doi: 10.1056/NEJMc1101502

PubMed Abstract | CrossRef Full Text | Google Scholar

Neumann, H. P., Vortmeyer, A., Schmidt, D., Werner, M., Erlic, Z., Cascon, A., et al. (2007). Evidence of MEN-2 in the original description of classic pheochromocytoma. N. Engl. J. Med. 357, 1311–1315. doi: 10.1056/NEJMoa071407

PubMed Abstract | CrossRef Full Text | Google Scholar

Neumann, H. P. H., Young, W. F., Jr., Eng, C. (2019). Pheochromocytoma and Paraganglioma. N. Engl. J. Med. 381, 552–565. doi: 10.1056/NEJMra1806651

PubMed Abstract | CrossRef Full Text | Google Scholar

Pang, Y., Lu, Y., Caisova, V., Liu, Y., Bullova, P., Huynh, T. T., et al. (2018). Targeting NAD(+)/PARP DNA repair pathway as a novel therapeutic approach to SDHB-mutated cluster I pheochromocytoma and paraganglioma. Clin. Cancer Res. 24, 3423–3432. doi: 10.1158/1078-0432.CCR-17-3406

PubMed Abstract | CrossRef Full Text | Google Scholar

Rednam, S. P., Erez, A., Druker, H., Janeway, K. A., Kamihara, J., Kohlmann, W. K., et al. (2017). Von Hippel-Lindau and hereditary pheochromocytoma/paraganglioma syndromes: clinical features, genetics, and surveillance recommendations in childhood. Clin. Cancer Res. 23, e68–e75. doi: 10.1158/1078-0432.CCR-17-0547

PubMed Abstract | CrossRef Full Text | Google Scholar

Remacha, L., Comino-Mendez, I., Richter, S., Contreras, L., Curras-Freixes, M., Pita, G., et al. (2017). Targeted exome sequencing of Krebs cycle genes reveals candidate cancer-predisposing mutations in pheochromocytomas and paragangliomas. Clin. Cancer Res. 23, 6315–6324. doi: 10.1158/1078-0432.CCR-16-2250

PubMed Abstract | CrossRef Full Text | Google Scholar

Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

Taieb, D., Pacak, K. (2017). New insights into the nuclear imaging phenotypes of cluster 1 pheochromocytoma and paraganglioma. Trends Endocrinol. Metab. 28, 807–817. doi: 10.1016/j.tem.2017.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Toledo, R. A., Qin, Y., Cheng, Z. M., Gao, Q., Iwata, S., Silva, G. M., et al. (2016). Recurrent mutations of chromatin-remodeling genes and kinase receptors in pheochromocytomas and paragangliomas. Clin. Cancer Res. 22, 2301–2310. doi: 10.1158/1078-0432.CCR-15-1841

PubMed Abstract | CrossRef Full Text | Google Scholar

Vo, J. N., Cieslik, M., Zhang, Y., Shukla, S., Xiao, L., Zhang, Y., et al. (2019). The landscape of circular RNA in cancer. Cell. 176, 869–881.e813. doi: 10.1016/j.cell.2018.12.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Zhang, C., Wu, Z., Chen, Y., Shi, W. (2018). CircIBTK inhibits DNA demethylation and activation of AKT signaling pathway via miR-29b in peripheral blood mononuclear cells in systemic lupus erythematosus. Arthritis Res. Ther. 20, 118. doi: 10.1186/s13075-018-1618-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Liu, J., Ma, J., Sun, T., Zhou, Q., Wang, W., et al. (2019). Exosomal circRNAs: biogenesis, effect and application in human diseases. Mol. Cancer 18, 116. doi: 10.1186/s12943-019-1041-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, F., Han, B., Wu, S., Yang, L., Leng, S., Li, M., et al. (2019). Circular RNA TLK1 aggravates neuronal injury and neurological deficits after ischemic stroke via miR-335-3p/TIPARP. J. Neurosci. 39, 7369–7393. doi: 10.1523/JNEUROSCI.0299-19.2019

PubMed Abstract | CrossRef Full Text | Google Scholar

Yao, Y., Chen, X., Yang, H., Chen, W., Qian, Y., Yan, Z., et al. (2019). Hsa_circ_0058124 promotes papillary thyroid cancer tumorigenesis and invasiveness through the NOTCH3/GATAD2A axis. J. Exp. Clin. Cancer Res. 38, 318. doi: 10.1186/s13046-019-1321-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Yi, X., Guan, X., Chen, C., Zhang, Y., Zhang, Z., Li, M., et al. (2018). Adrenal incidentaloma: machine learning-based quantitative texture analysis of unenhanced CT can effectively differentiate sPHEO from lipid-poor adrenal adenoma. J. Cancer 9, 3577–3582. doi: 10.7150/jca.26356

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, J., Yang, M., Zhou, B., Luo, J., Zhang, Z., Zhang, W., et al. (2019). CircRNA-104718 acts as competing endogenous RNA and promotes hepatocellular carcinoma progression through microRNA-218-5p/TXNDC5 signaling pathway. Clin. Sci. 33, 1487–1503. doi: 10.1042/CS20190394

CrossRef Full Text | Google Scholar

Zeng, X., Lin, W., Guo, M., Zou, Q. (2017). A comprehensive overview and evaluation of circular RNA detection tools. PloS Comput. Biol. 13, e1005420. doi: 10.1371/journal.pcbi.1005420

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Wang, Q., Wang, Q., Liu, A., Qin, F., Sun, Q., et al. (2019). Circular RNA expression profiling in the nucleus accumbens: Effects of electroacupuncture treatment on morphine-induced conditioned place preference. Addict. Biol., e12794. doi: 10.1111/adb.12794

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhi, X., Zhang, J., Cheng, Z., Bian, L., Qin, J. (2019). circLgr4 drives colorectal tumorigenesis and invasion through Lgr4-targeting peptide. Int. J. Cancer. doi: 10.1002/ijc.32549

CrossRef Full Text | Google Scholar

Zhu, Y. J., Zheng, B., Luo, G. J., Ma, X. K., Lu, X. Y., Lin, X. M., et al. (2019). Circular RNAs negatively regulate cancer stem cells by physically binding FMRP against CCAR1 complex in hepatocellular carcinoma. Theranostics. 9, 3526–3540. doi: 10.7150/thno.32796

PubMed Abstract | CrossRef Full Text | Google Scholar