Online ISSN: 1945-4589
Learn about our FREE Post-Publication Promotion Services
Research Paper Volume 15, Issue 17 pp 9193—9216

Identification of a novel competing endogenous RNA network and candidate drugs associated with ferroptosis in aldosterone-producing adenomas

Yu Hanxiao1, *, , Yang Boyun2, *, , Jia Minyue3, , Song Xiaoxiao4, ,

  • 1 Clinical Research Center, The Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China
  • 2 Department of Allergy, The Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China
  • 3 Department of Ultrasound, The Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China
  • 4 Department of Endocrinology and Metabolism, The Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China
* Equal contribution

Received: May 1, 2023       Accepted: August 22, 2023       Published: September 13, 2023      

https://doi.org/10.18632/aging.205028
How to Cite

Copyright: © 2023 Hanxiao et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Aldosterone-producing adenoma (APA), characterized by unilaterally excessive aldosterone production, is a common cause of primary aldosteronism. Ferroptosis, a recently raised iron-dependent mode of programmed cell death, has been involved in the development and therapy of various diseases. This study obtained datasets of the mRNA and lncRNA expression profiles for APA and adjacent adrenal gland (AAG) from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) and lncRNAs (DE lncRNAs) associated with ferroptosis were identified. Enrichment analyses indicated 89 ferroptosis-related DEGs were primarily enriched in ROS related processes and ferroptosis. Two physical cores, and one combined core were identified in the protein–protein interaction (PPI). DEGs and clinical traits were used in conjunction to screen eight hub genes from two hub modules and 89 DEGs. A competitive endogenous RNA (ceRNA) network was constructed via co-express analysis. Thereafter, molecular docking was used to identify potential targets. Two active compounds, QL-X-138 and MK-1775, bound to AURKA and DUOX1, respectively, with the lowest binding energies. Molecular dynamics simulation verified the stability of the two complexes. In summary, our studies identified eight hub genes and a novel ceRNA regulatory network associated with ferroptosis, wherein QL-X-138 and MK-1775 were considered to be potential drugs for treating APA.

Introduction

Primary aldosteronism caused by increased aldosterone levels is well established to be linked to secondary hypertension, with an estimated prevalence of nearly 6% in 1.1 billion hypertensive patients of the worldwide [1, 2], and a prevalence of up to 20% in resistant hypertensive patients [3]. The diagnosis of primary aldosteronism is based on hypertension related to an elevated aldosterone-to-renin ratio and usually hypokalemia [4, 5]. APA as well as idiopathic hyperaldosteronism [4] are the major causes of primary aldosteronism, which is characterized by autonomous production of aldosterone in the adrenal cortex [6]. Although the renin-angiotensin-II-aldosterone system has vigorous feedback regulation, aldosterone production in APA is normally uncontrollable even under low renin conditions [7]. Somatic mutations are known to promote the production of aldosterone, along with the associated genes that mainly encoding ion channels such as KCNJ5 [8], CACNA1D [9, 10], CACNA1H [11], CLCN2 [5, 12] and ATPases (ATP1A1 [9] and ATP2B3 [13]) result in elevating Ca2+ concentration in zona glomerulosa cells [14, 15]. Although most APAs can be correctable for hypertension via adrenalectomy, only 37% of patients achieve complete clinical success [16]. Thus, noninvasive treatment, fewer side effects and improved drug tolerance are necessary to either compete for the mineralocorticoid receptor with aldosterone or suppress the Ca2+ concentration [17]. Novel target drugs can be identified by programmed cell death of tumor cells and regulatory networks combined with the aldosterone production pathway through gene expression profile.

Programmed cell death can be classified into six forms: apoptosis, necroptosis, autophagy, ferroptosis, pyroptosis and necrosis [18]. Since the discovery of the iron-dependent form of cell death termed ferroptosis [19], there has been a rapid increase in research on the same [20]. Ferroptosis is different from other programmed cell deaths in terms of morphology, biochemical features and regulatory genes [19, 21]. It depends on the accumulation of cell ROS and excess lipid peroxidation, which induce oxidative damage in the cell membrane [19]. Excessive iron causes lipid peroxidation through ROS production and the iron-associated Fenton reaction [22]. Moreover, certain proteins that transfer membrane electrons, such as NADPH oxidases (NOX) [19, 23], also have a positive effect on lipid peroxidation via ROS production in tumors [24]. In contrast, growing evidence has approved that some drugs, including ROS-inducing drugs [25, 26], cytokines [27] and ionizing radiation [28], can function as proferroptotic inducers. Conversely, antiferroptotic inhibitors can cause inflammation-related immune suppression [24]. Ferroptosis may be inhibited in malignant tumors [29] and activated in neurodegenerative diseases [30]. Recently Yuhong et al. revealed that the differentially expressed gene BEX1 had a potent effect on protecting adrenocortical cells from ferroptosis in micro-APAs compared to macro-APAs [15]. The sequencing on APA with or without KCNJ5 mutation of spatial metabolomics showed abundant accumulation of proferroptotic metabolites, such as adrenic acid and docosapentaenoic acid [31]. lncRNAs also engage in the ferroptosis process by impacting post-transcriptional alterations and protein stability, as observed in contemporary studies [32]. However, the molecular mechanisms by which ferroptosis-related genes and lncRNAs regulate aldosterone production remain unclear.

Fortunately, with the improvement of sequencing and bioinformatics, large sample sizes data from multiple independent studies can be integrated and explored potential information by researchers. Yang et al. reported the relationship between tumor immune genes and prognosis in ovarian cancer [33] and breast cancer [34, 35] based on the multiple cohorts in the TCGA and GEO databases. In this study, we download mRNA and lncRNA expression profile of APA and AAG tissue with corresponding clinical characteristics from the GEO database. Subsequently, the ferroptosis-related DEGs and DE lncRNAs were screened out, and the enriched pathway was analyzed by Metascape and Gene Set Enrichment Analysis (GSEA). Based on the ferroptosis-related DEGs and DE lncRNAs, a ceRNA network was constructed by predicting target miRNAs. Finally, drugs targeting hub genes related to ferroptosis were anticipated by molecular docking and molecular dynamics simulation.

Results

Identification of DE lncRNAs and ferroptosis-related DEGs

The datasets of GSE64957 and GSE101894 were selected to identify the DEGs and DE lncRNAs in APA tissue compared to that in the AAG tissue, as depicted in the volcano plots (Figure 1A, 1B). In addition, 263 ferroptosis-related genes from FerrDb were used to screen the ferroptosis-related DEGs. Consequently 89 ferroptosis-related DEGs were identified, including 44 upregulated genes and 45 downregulated genes (Figure 1C). The expression of 89 genes in the APA and AAG tissue are presented in the heatmap (Figure 1D).

Screening for the ferroptosis-related DEGs. (A, B) Volcano plots of the DEGs in GSE64957 and GSE101894 between APA and AAG tissue. The red dots represent the upregulated genes, the green dots represent the downregulated genes, and the black dots represent genes with no significant difference. (C) Venn diagram to identify ferroptosis-related DEGs. (D) Expression heatmap of 89 ferroptosis-related DEGs in APA and AAG tissue.

Figure 1. Screening for the ferroptosis-related DEGs. (A, B) Volcano plots of the DEGs in GSE64957 and GSE101894 between APA and AAG tissue. The red dots represent the upregulated genes, the green dots represent the downregulated genes, and the black dots represent genes with no significant difference. (C) Venn diagram to identify ferroptosis-related DEGs. (D) Expression heatmap of 89 ferroptosis-related DEGs in APA and AAG tissue.

Enrichment analysis of 89 ferroptosis-related DEGs

Metascape was used to analyze the functions and pathways of 89 ferroptosis-related DEGs in patients with APA. This tool could identify all statistically enriched terms that were selected in this study, including Gene Ontology biological processes (GO BP), reactome gene sets, KEGG pathways, canonical pathways, and other pathways. Then, the enrichment factors and accumulative hypergeometric P-values were calculated for the filtrate. Based on the Kappa-statistical similarities among the enriched genes, the remaining significant terms were hierarchically clustered into a tree [36]. The top 20 most highly enriched items are presented in Figure 2A. Among these items, the response to oxidative stress, reactive oxygen species metabolic process, superoxide metabolic process and positive cell death regulation of were associated with the production of oxidative stress in the GO BP category. This suggested that these genes enriched in the items played a role in the homeostasis and metabolism of superoxides. Furthermore, KEGG pathway demonstrated that the DEGs were mainly concentrated in the ferroptosis, hypoxia-inducible factor 1 (HIF1) TF (involved in the response to hypoxia), and IL-18 signaling pathways (Figure 2A). Upon selecting the top 20 terms from the clusters of functions and pathways, this tool transformed these terms into a network exhibition (Figure 2B). When the Kappa score was greater than 0.3, the node network could be connected. More specifically, each term was represented by a circle whose size was proportional to the number of genes enriched in that term. In addition, the same function or pathway clusters were indicated by an identical color. These terms were connected to an edge using a similarity score, which determined the thickness of the edge. The GO BP items, response to oxidative stress and reactive oxygen species metabolic processes, had a greater number of genes. Furthermore, they were connected to the ferroptosis pathway and positive regulation of cell death more intensively. To further investigate the relationship between 89 ferroptosis-related DEGs and enriched function as well as pathways, a chord plot was drawn to visualize the DEGs and top 10 enriched terms (Figure 2C). Their responses to oxidative stress, cellular responses to stress and positive regulation of cell death were shown to comprise the major genes. Twenty-three DEGs enriched more than three terms suggesting that they participated in regulating multiple pathways and played a key role in ferroptosis-related functions.

Enriched functions and pathways of ferroptosis-related DEGs. (A) The top 20 most significant enriched terms analyzed from Metascape, ranked according to –log10 (P) value. (B) The interaction network of 20 enriched terms, each identical color circle represents an enriched term. (C) The relationship between ferroptosis-related DEGs and top 10 enriched terms.

Figure 2. Enriched functions and pathways of ferroptosis-related DEGs. (A) The top 20 most significant enriched terms analyzed from Metascape, ranked according to –log10 (P) value. (B) The interaction network of 20 enriched terms, each identical color circle represents an enriched term. (C) The relationship between ferroptosis-related DEGs and top 10 enriched terms.

PPI network analysis, physical and combined cores

The interaction network of proteins coded by 89 ferroptosis-related DEGs was constructed by Metascape using the MCODE module (Figure 3A). The MCODE module was utilized to screen gene clusters, as per the algorithm described in this method. Three key subnetworks were identified if the network contained between three and five hundred proteins, and one interacted with at least one other member in the list (Figure 3B–3D). The two physical cores contained 11 (Figure 3B) and 9 DEGs (Figure 3C), respectively, suggesting that these DEGs regulated more physical pathways. Furthermore, functions and pathways enriched analyses were independently applied to the three key subnetworks. The DEGs in the first physical core were primarily associated with the cellular response to chemical stress, response to oxidative stress, as well as cellular response to oxidative stress (Figure 3B). DEGs in the second physical core were strongly related to the AGE/RAGE signaling pathway, which played an important role in the production of ROS mediated by NADPH or derived from NOX (Figure 3C) [37, 38]. The combined core, comprising 12 DEGs, had the greatest number of nodes and edges among all the subnetworks (Figure 3D). It was associated with the apoptotic signaling and extrinsic apoptotic signaling pathways, suggesting the involvement of alternative programmed cell death. Moreover, the NOX family (including NOX4, NOX5, DUOX1, and DUOX2 in this combined core) was demonstrated to contribute both to ferroptosis and apoptosis by generating ROS [39]. This combined core might comprise an axis such as the FBW7-NRA41-SCD1 axis [40] that synchronously regulated ferroptosis and apoptosis.

PPI network and three MCODE clusters. (A) A connected interaction network of 89 ferroptosis-related DEGs. (B, C) Two physical cores extracted from the PPI network. (D) The combined core extracted from the PPI network. Each node represents a protein, while each edge represents the relation of protein-protein.

Figure 3. PPI network and three MCODE clusters. (A) A connected interaction network of 89 ferroptosis-related DEGs. (B, C) Two physical cores extracted from the PPI network. (D) The combined core extracted from the PPI network. Each node represents a protein, while each edge represents the relation of protein-protein.

WGCNA and hub genes identification

To screen the hub genes from all genes associated with clinical characteristics, the genes expression profile and clinical traits of patients with APA were combined to construct a gene co-expression network. Weighted gene co-expression network analysis (WGCNA) is a bioinformatic analysis that is used to investigate the relationships between genes and phenotypes efficiencies [41]. The clinical traits, including APA, age, gender, systolic blood pressure (SBP), diastolic blood pressure (DBP), serum K+ (K), plasma aldosterone (ALD), and plasma renin of APA as well as AAG samples were clustered and depicted on a heatmap (Figure 4A). When the soft threshold was set 8, the index of scale-free topological network was 0.9. Hence, the network was closer to the actual network (Figure 4B). The created gene dendrograms and corresponding module colors were drawn in Figure 4C. The genes were divided into 22 modules obtained from the visualized correlations between module eigengenes and clinical characteristics (Figure 4D). APA has a complex pathogenic mechanism of involving serum K+, plasma aldosterone and renin et al. Here, we took multiple clinicopathological factors into consideration and analyzed the correlation between each gene module. The APA, SBP and ALD of clinical characteristics were more positively relevant to the module yellow, and more negatively relevant to the module brown of all modules. Moreover, the most positive and negative correlation coefficient between APA and module yellow as well as brown were 0.85 and -0.78, respectively. This indicated that the genes in module yellow may contribute to the pathogenesis and development of APA, SBP and ALD. In contrast, the genes in module brown may repress the progression of APA, SBP and ALD. The hub genes were identified from the intersection of modules yellow and brown, along with that of ferroptosis-related DEGs (Figure 5A). Remarkably, IL6 and VEGFA were members of physical cores, whereas NOX4, DUOX1, DUOX2, and BNIP3 interacted with other genes in the combined core from PPI network (Figure 3B–3D). This suggested that these genes played an important role in the subnetworks.

WGCNA analysis and identification of hub modules related to the clinical traits. (A) Clustering dendrograms of samples and clinical traits of APA, age, gender, SBP, DBP, K, ALD and renin showed at the bottom. (B) Analysis of scale free network topology for various soft-threshold powers. The image on the left demonstrates the effect of soft-threshold power on the scale free topology model fit index; the right image shows the effect of soft-threshold power on the mean connectivity of the network. (C) Gene clustering dendrograms with dissimilarity based on the topological overlap, assigned with specific module colors. (D) The heatmap of correlation coefficient and p-value between modules and clinical traits. Each row represents a module, and each column corresponds to a clinical trait.

Figure 4. WGCNA analysis and identification of hub modules related to the clinical traits. (A) Clustering dendrograms of samples and clinical traits of APA, age, gender, SBP, DBP, K, ALD and renin showed at the bottom. (B) Analysis of scale free network topology for various soft-threshold powers. The image on the left demonstrates the effect of soft-threshold power on the scale free topology model fit index; the right image shows the effect of soft-threshold power on the mean connectivity of the network. (C) Gene clustering dendrograms with dissimilarity based on the topological overlap, assigned with specific module colors. (D) The heatmap of correlation coefficient and p-value between modules and clinical traits. Each row represents a module, and each column corresponds to a clinical trait.

Identification of eight hub ferroptosis-related DEGs and the expression among mutant KCNJ5 and WT APA samples. (A) Identification of BNIP3, DUOX1, DUOX2, NOX4, HNF4A, AURKA, VEGFA and IL6 from the intersection of hub module yellow and brown, along with 89 ferroptosis-related DEGs. (B) The expression of eight hub genes among the mutant KCNJ5 and WT groups. *P 0.05.

Figure 5. Identification of eight hub ferroptosis-related DEGs and the expression among mutant KCNJ5 and WT APA samples. (A) Identification of BNIP3, DUOX1, DUOX2, NOX4, HNF4A, AURKA, VEGFA and IL6 from the intersection of hub module yellow and brown, along with 89 ferroptosis-related DEGs. (B) The expression of eight hub genes among the mutant KCNJ5 and WT groups. *P < 0.05.

The expression of eight hub genes in mutant KCNJ5 and WT APA samples

Somatic mutations are a common cause of APA. The general prevalence of KCNJ5 mutations was observed to be 43% in patients with APA [42]. The ComBat was used to remove the batch effect to compare the expression of eight hub genes, and the PCA diagram showed the before and after correction (Supplementary Figure 1). The expression of NOX4 and DUOX2, two members of NOX family, was significantly differentiated in mutant KCNJ5 samples compared to that in the WT (Figure 5B). This indicated the impact of the KCNJ5 mutation on the expression of ferroptosis-related genes.

Construction of ceRNA network

It is well known that miRNAs are recognized as targeted regulators by binding mRNAs. The lncRNA-miRNA-mRNA regulatory network across transcriptional talk is formed as a ceRNA [43]. This study identified miRNAs predicted by online tools using gene-miRNA and miRNA-lncRNA target modules. The miRNAs were selected by the targeted intersection of eight hub ferroptosis-related DEGs and correlated DE lncRNAs. The key lncRNA emerged according to the positively correlated to the expression of mRNAs and that located in the cytoplasm. Consequently, based on the interaction of lncRNA-miRNA-mRNA, the ceRNA was constructed with 26 lncRNAs, 132 miRNAs and 4 mRNAs (Figure 6C). Among all the predicted miRNAs, the ones with the highest degree of connection were miR-942, miR-342, miR-877, miR-511, miR-655, miR-4279, miR-3121, miR-339, miR-5571 and miR-1343. This result identified couples of lncRNA-miRNA-mRNA associations, including lnc-RASL11B (ENSG00000248115)-miR-335-AURKA, and lnc-ZBTB18 (ENSG00000226828)-miR-4296/4265/4322-DUOX1 axes, which might be important therapeutic targets for APA.

The analysis of molecular crosstalk. (A) The correlation between eight hub ferroptosis-related DEGs. (B) The correlation analysis between 8 hub ferroptosis-related DEGs and 10 DE lncRNAs. (C) The ceRNA network interaction among lncRNAs, miRNAs, and mRNAs. (D) Two axes ENSG00000248115-hsa-miR-335-AURKA and ENSG00000226828- hsa-miR-4265/hsa-miR-4296/hsa-miR-4322-DUOX1. Green squares represent lncRNAs; blue circles represent miRNAs; red triangles represent mRNAs; the size of the shape represents the degree of connection. *P 0.05; **P 0.01; ***P 0.001.

Figure 6. The analysis of molecular crosstalk. (A) The correlation between eight hub ferroptosis-related DEGs. (B) The correlation analysis between 8 hub ferroptosis-related DEGs and 10 DE lncRNAs. (C) The ceRNA network interaction among lncRNAs, miRNAs, and mRNAs. (D) Two axes ENSG00000248115-hsa-miR-335-AURKA and ENSG00000226828- hsa-miR-4265/hsa-miR-4296/hsa-miR-4322-DUOX1. Green squares represent lncRNAs; blue circles represent miRNAs; red triangles represent mRNAs; the size of the shape represents the degree of connection. *P < 0.05; **P < 0.01; ***P < 0.001.

The impact on co-expression of eight hub ferroptosis-related DEGs and DE lncRNAs

To further elucidate the pathogenesis of APA at the transcriptome level, the co-expression of genes and lncRNAs was analyzed. All eight hub ferroptosis-related DEGs correlated with one another at different levels (P < 0.05) (Figure 6A). Among them, AURKA and BNIP3, AURKA and VEGFA, VEGFA and DUOX2, DUOX1 and NOX4, DUOX1 and HNF4A, DUOX1 and IL6, NOX4 and HNF4A, NOX4 and IL6 were the eight most positively correlated pairs of genes. In contrast, the gene pairs with the most negative correlations were HNF4A and BNIP3, NOX4 and VEGFA, NOX4 and BNIP3, NOX4 and AURKA, DUOX1 and VEGFA. The strong correlation between the eight genes suggested that the genes may impact one another during their expressions in APA. Moreover, to further understand the expression of these eight genes, which are regulated by DE lncRNAs, we selected the top 10 DE lncRNAs with the highest degree of connection using Pearson’s correlation coefficient (Figure 6B). Consequently, a strong positive correlation was observed between TUSC7 and VEGFA. However, the lncRNAs AL109946, AC245060, AC008543, CR786580, MALAT1, AL627309 and genes HNF4A, NEAT1 and NOX4 exhibited a strong negative correlation.

GSEA analysis of the top 10 lncRNAs

In contrast to the gene enrichment of functions and pathways, the top 10 DE lncRNAs enriched pathways that showed the highest degree of connection with GSEA were also analyzed. The results revealed the majority of 10 lncRNAs regulated immune and aldosterone-biosynthesis related pathway, such as the renin-angiotensin system, steroid biosynthesis, calcium signaling pathway, primary immunodeficiency, chemokine signaling pathway, and natural killer cell mediated cytotoxicity (Figure 7).

Six most significant enriched KEGG pathway of 10 lncRNAs. (A) Renin-angiotensin system (NES=-1.913, P 0.05). (B) Steroid biosynthesis (NES=-1.999, P 0.05). (C) Calcium signaling pathway (NES=-2.162, P 0.05). (D) Primary immunodeficiency (NES=-1.884, P 0.05). (E) Chemokine signaling pathway (NES=-2.475, P 0.05). (F) Natural killer cell mediated cytotoxicity (NES=-2.118, P 0.05). NES: normalized enrichment score.

Figure 7. Six most significant enriched KEGG pathway of 10 lncRNAs. (A) Renin-angiotensin system (NES=-1.913, P < 0.05). (B) Steroid biosynthesis (NES=-1.999, P < 0.05). (C) Calcium signaling pathway (NES=-2.162, P < 0.05). (D) Primary immunodeficiency (NES=-1.884, P < 0.05). (E) Chemokine signaling pathway (NES=-2.475, P < 0.05). (F) Natural killer cell mediated cytotoxicity (NES=-2.118, P < 0.05). NES: normalized enrichment score.

Prediction of candidate drugs through CMAP

To further enhance the precision of drugs screening, 44 upregulated and 45 downregulated ferroptosis-related genes were queried on the CMAP. The detailed results showed that the top ranked drugs, including fludrocortisone and fluocinolone, were positively associated with APA suggesting that they promoted APA progression or aldosterone production (Table 1). Conversely, the bottom ranked drugs that were negatively associated with APA included MK-1775, GW-843682X, QL-X-138 and OM-137, which were the candidate drugs. Drugs with potent adverse reactions or potential toxicities were excluded from this study.

Table 1. The positive and negative small molecular compounds for APA obtained from the connectivity map (CMAP) database.

RankCMAP nameBelongs toCanonical SMILES
1ApicidinHDAC1 inhibitorCCC(C)C1C(=O)N2CCCCC2C(=O)NC(C(=O)NC(C(=O)N1)CC3=CN(C4=CC=CC=C43)OC)CCCCCC(=O)CC
2EntinostatHDAC inhibitorC1=CC=C(C(=C1)N)NC(=O)C2=CC=C(C=C2)CNC(=O)OCC3=CN=CC=C3
3FludrocortisoneGlucocorticoid receptor agonistCC12CCC(=O)C=C1CCC3C2(C(CC4(C3CCC4(C(=O)CO)O)C)O)F
4FluocinoloneNACC12CC(C3(C(C1CC(C2(C(=O)CO)O)O)CC(C4=CC(=O)C=CC43C)F)F)O
5ClobetasolNACC1CC2C3CCC4=CC(=O)C=CC4(C3(C(CC2(C1(C(=O)CCl)O)C)O)F)C
2831OM-137Aurora kinase inhibitorCC1=C(SC(=N1)N)C(=O)NN=CC2=CC(=C(C=C2)O)OC
2832QL-X-138mTOR inhibitorCC1=C(C=C(C=C1)N2C(=O)C=CC3=CN=C4C=CC(=CC4=C32)C5=CNN=C5)NC(=O)C=C
2833GW-843682XNACOC1=C(C=C2C(=C1)N=CN2C3=CC(=C(S3)C(=O)N)OCC4=CC=CC=C4C(F)(F)F)OC
2834CycloheximideNACC1CC(C(=O)C(C1)C(CC2CC(=O)NC(=O)C2)O)C
2835EmetineProtein synthesis inhibitorCCC1CN2CCC3=CC(=C(C=C3C2CC1CC4C5=CC(=C(C=C5CCN4)OC)OC)OC)OC
2836GestrinoneAndrogen receptor modulatorCCC12C=CC3=C4CCC(=O)C=C4CCC3C1CCC2(C#C)O
2837MK-1775NACC(C)(C1=NC(=CC=C1)N2C3=NC(=NC=C3C(=O)N2CC=C)NC4=CC=C(C=C4)N5CCN(CC5)C)O
1HDAC, histone deacetylase.

Identification of candidate drugs by molecular docking and molecular dynamics simulation

A preliminarily investigation was performed to verify whether there was direct targeting of the eight hub ferroptosis-related proteins of candidate drugs. The molecular docking results evaluated with binding energy have been presented in Table 2 (excluding the NOX4 and DUOX2 as unspecified protein structures). The small compound, QL-X-138, had the highest affinity but the lowest binding energy when targeted towards AURKA (-9.8 kcal/mol). Similarly, MK-1775 showed the highest affinity but the lowest binding energy when targeted towards DUOX1 (-9.3 kcal/mol). This analysis revealed the formation of two stable conformations with the lowest binding energies, among the 24 docking results. Next, the interactions between small compounds and amino acid residues were performed by Pymol and Ligplot. The ligand QL-X-138 formed a hydrogen bond with residue Lys162 and a pi-cation bond with Trp277 of protein AURKA (Figure 8A, 8C). While binding to DUOX1, MK-1775 formed a hydrogen bond with residue His1319 (Figure 8B, 8D). Two axes related to the two proteins in the ceRNA network were obtained (Figure 6D). The miRNAs and lncRNAs selected depend on the basis of their degree of connection might regulate AURKA and DUOX1 through ENSG00000248115-hsa-miR-335-AURKA and ENSG00000226828-hsa-miR-4265/hsa-miR-4296/hsa-miR-4322-DUOX1 axes at the transcriptome level, respectively.

Table 2. The predicted binding energy between candidate drugs and hub genes.

GenePDB identifierTest compoundBinding energy (kcal/mol)
AURKA1MUOQL-X-138-9.8
DUOX17D3FMK-1775-9.3
AURKA1MUOMK-1775-9.1
AURKA1MUOGW-843682X-9.1
DUOX17D3FGW-843682X-8.9
DUOX17D3FQL-X-138-8.6
VEGFA1BJ1MK-1775-8.4
BNIP32KA2QL-X-138-8.4
BNIP32KA2MK-1775-8.0
HNF4A1PZLGW-843682X-8.0
VEGFA1BJ1QL-X-138-7.8
VEGFA1BJ1GW-843682X-7.5
IL65FUCGW-843682X-7.3
BNIP32KA2GW-843682X-7.3
IL65FUCMK-1775-7.2
HNF4A1PZLOM-137-7.0
IL65FUCQL-X-138-6.8
AURKA1MUOOM-137-6.8
VEGFA1BJ1OM-137-6.4
DUOX17D3FOM-137-6.3
IL65FUCOM-137-5.9
BNIP32KA2OM-137-5.9
HNF4A1PZLMK-1775-5.4
HNF4A1PZLQL-X-138-3.0
The molecular docking mode of small compounds and proteins. (A) The three-dimensional structure of interaction between QL-X-138 and AURKA. (B) MK-1775 and DUOX1. (C) The two-dimensional structure of interaction between QL-X-138 and AURKA. (D) MK-1775 and DUOX1. The yellow dotted line represents the hydrogen bond, while the red dotted line represents the pi-cation bond.

Figure 8. The molecular docking mode of small compounds and proteins. (A) The three-dimensional structure of interaction between QL-X-138 and AURKA. (B) MK-1775 and DUOX1. (C) The two-dimensional structure of interaction between QL-X-138 and AURKA. (D) MK-1775 and DUOX1. The yellow dotted line represents the hydrogen bond, while the red dotted line represents the pi-cation bond.

Moreover, the conformational stability of these two complexes was also verified by molecular dynamics simulation. The protein Calpha root-mean-square deviation (RMSD) data of QL-X-138 and AURKA, as well as MK-1775 and DUOX1, were obtained in terms of protein mainchains and protein-ligands. In comparison with the first frame in Figure 9A, 9D, the RMSD of the protein backbone of QL-X-138 and AURKA system became stable after 10 ns, with a deviation of 0.19. While the MK-1775 and DUOX1 system tended to become stable after 40 ns, with a deviation of 0.57. As for the protein-ligand RMSD, the system of QL-X-138 and AURKA, MK-1775 and DUOX1 both achieved equilibrium at 10 ns and 40 ns, with a deviation of 0.25 and 0.61, respectively (Figure 9B, 9E). Meanwhile, the RMSD of ligand QL-X-138 and MK-1775 was also stable at 5 ns and 36 ns as shown in Figure 9C, 9F, with a deviation of 0.09 and 0.14, respectively. Additionally, after reaching stability, the RMSD of the protein backbone, protein-ligand and ligand of the QL-X-138 and AURKA system maintained up to 45 ns. For all systems, the RMSD became stable after 40 ns. These data provided computational insight into the contribution of molecular dynamics simulation to the stable conformation of the systems of QL-X-138 and AURKA, MK-1775 and DUOX1.

The RMSD plot in the process of molecular dynamics simulation. (A) The RMSD of protein AURKA. (B) The RMSD of protein-ligand QL-X-138 and AURKA system. (C) The RMSD of ligand QL-X-138. (D) The RMSD of protein DUOX1. (E) The RMSD of protein-ligand MK-1775 and DUOX1 system. (F) The RMSD of ligand MK-1775.

Figure 9. The RMSD plot in the process of molecular dynamics simulation. (A) The RMSD of protein AURKA. (B) The RMSD of protein-ligand QL-X-138 and AURKA system. (C) The RMSD of ligand QL-X-138. (D) The RMSD of protein DUOX1. (E) The RMSD of protein-ligand MK-1775 and DUOX1 system. (F) The RMSD of ligand MK-1775.

Discussion

Since its discovery, ferroptosis has been associated with various human diseases, including cancer [24], neurological diseases [44], diabetes mellitus [45], and retinal degeneration [46] et al. However, to our knowledge, little attention has been paid to the involvement of ferroptosis in the pathogenesis of APA. In our study, we first identified a novel ceRNA regulatory network based on ferroptosis-related DEGs and DE lncRNAs. Furthermore, we explored candidate drugs that could target the eight hub genes for the APA treatment with molecular docking and molecular dynamics simulation. These findings revealed the pivotal genes involved in the regulation and therapy through ferroptosis signaling pathways. Although ferroptosis and immune infiltration have been implicated in pathogenesis of various tumors with abundant studies, little is known about the latent effect on APA [20, 33, 34]. Additionally, using transcriptomic and microarray sequencing analyses to identify the new biomarkers can improve APA diagnosis and treatment as they have been used in cancer [35].

We identified 89 ferroptosis-related DEGs associated with APA. The enrichment analysis of GO BP and KEGG pathways indicated that these genes mainly regulated cellular functions including responses to oxidative stress, reactive oxygen species metabolic process, superoxide metabolic process and HIF1 TFpathway. Morphological analysis revealed that ferroptosis occurs with elevated bilayer membrane density and decreased mitochondrial volume, resulting in cell membrane intact. Biochemical analyses have shown that glutathione depletion and lipid peroxide accumulation with reduced activity of GPX4 occurred in the process of ferroptosis [47]. Multiple signal molecules and metabolic pathways are involved in the occurrence and development of ferroptosis, such as HIF-1α, NADPH oxidase, GPX4, iron and lipid metabolism. Recently, a study in kidney injury demonstrated that ferroptosis induced by folic acid was decreased with FG-4592 pretreatment, an inhibitor of prolyl hydroxylase of HIF, through the Akt/GSK-3β/Nrf2 pathway [48]. Ferroptosis, enhanced by excess ROS production, injured renal tubules, and promoted diabetic nephropathy via HIF-1α/HO-1 pathway in db/db mice [49]. Among the eight hub genes, NOX4, a major ROS producer, catalyzes the conversion of oxygen to superoxide using NADPH [50]. Min et al. observed that NOX4 increased the astrocytic ferroptosis level and ferroptosis-dependent cytotoxicity through lipid peroxidation and impaired mitochondrial metabolism in Alzheimer's diseases [51].

The accumulation of ROS, a major inducer of ferroptosis, triggers autophagy when ferroptosis is induced by erastin [52]. Autophagy and ferroptosis are two types of programmed cell death that remain controversial regarding relationship, and ROS may become a ‘bridge’ connecting them. Damage to the cell membrane by ROS may be initiated from lipid peroxidation chain reaction which leads to the ferroptosis, autophagy, or apoptosis [39], indicating that ROS played a critical role in concatenating ferroptosis and other programmed cell deaths. AURKA is a cell cycle-regulating kinase that regulates the formation of mitotic spindle during chromosome segregation [53]. However, AURKA has been reported to be upregulated in several tumors. A recent study revealed that the suppression of AURKA caused cell death via GPX4, a well-known ferroptosis inhibitor [54]. The inhibition of GPX4 resulted in the accumulation of lipid-based ROS, thereby inducing ferroptosis [55], suggesting a strong link between AURKA and ferroptosis. In this study, we identified the lnc-RASL11B (ENSG00000248115)-miR-335-AURKA regulated network.

The proinflammatory cytokine IL6 participates in the IL6/Jak2/Stat3 pathway to regulate iron and maintain iron homeostasis [56]. However, high expression of IL6 disturbed iron homeostasis and subsequently induced ferroptosis by producing lipid peroxidation [57]. One transcription factor, HNF4A, identified as a ferroptosis suppressor mediates one category of genes called ferroptosis downregulated factors (FDF) by affecting the biosynthesis of glutathione in liver cancer. Additionally, disintegration of KAT2B preferentially prevented HNF4A from binding to FDF promoter and facilitated HIC1 binding to ferroptosis upregulated factors [58]. Besides ROS-related genes, pleiotropic cytokine, serine/threonine kinase and even gene regulated FDF are involved in ferroptosis.

In general, 26 DE lncRNAs constituting the ceRNA network were associated with ferroptosis. The top 10 lncRNAs regulate KEGG pathways, including renin-angiotensin system, steroid biosynthesis, calcium signaling pathway, primary immunodeficiency, chemokine signaling pathway, and natural killer cell mediated cytotoxicity. One study revealed that the overexpression of lncRNA KCNQ1OT1 mediated the upregulation of L-type Ca2+ channels, CACNA1C, by binding miR-384 in angiotensin-II induced atrial fibrillation [59]. In cultured HeLa cells, under low oxygen conditions, lncRNA MALAT1 was upregulated by the CaMKK/AMPK/HIF-1α axis [60]. Increasing evidence has demonstrated the effects of mutant genes that encoding calcium ion channels and pumps on the production of aldosterone. Elevated intracellular calcium concentration in zona glomerulosa cells of the adrenal cortex induces CYP11B2 overexpression and autonomous production of aldosterone [61]. Additionally, the overexpressed lncRNA MEG3 interacted with miR-125a-5p and suppressed its expression, whereas the downregulation of MEG3 increased the expression of FOXP3 via miR-125a-5p in the CD4+ T cells in immune thrombocytopenic purpura [62]. Silencing of lncRNA NEAT1 by siRNA repressed lipid uptake and the inflammatory response by targeting miR-342-3p in THP-1 macrophage, providing a strategy to treat atherosclerosis [63].

KCNJ5, which encodes the potassium channel GIRK4, is the most frequent mutation in APA. The prevalence of APA induced by KCNJ5 mutation was 34% in Europe [64], while 73% in Japan [65] and 70.7 % in China [66]. In patients with APA, KCNJ5 achieved a better prognosis than those with WT based on the PASO [16, 66]. In our results, NOX4 and DUOX2 showed significantly differentially expressed in KCNJ5 mutant samples compared to WT samples that the two ferroptosis-related genes may be involved in pathogenesis and prognosis. The inhibition of PIEZO1, a calcium channel, could decrease ferroptosis-related process such as ROS and lipid peroxidation in pulmonary endothelial cells, which is accompanied by calcium efflux [67], suggesting that calcium ion channel may impact the ferroptosis via calcium influx/efflux. ROS produced by NOX4 decreased the expression of Kv1.5 (a potassium channel) by activating P-Smad2/3 and P-ERK ½ [68], indicating that NOX4 may influence potassium channel by increasing ROS levels. However, there is little evidence regarding the effect of KCNJ5 mutation on the ferroptosis in APA, and NOX4 and DUOX2 may be novel targets for APA investigation and therapy.

Growing evidence demonstrated that various complicated diseases can scarcely be explained by one or a few genomic variations. Importantly, constructing a ceRNA network could help us understand the entire protein-coding genes dimension through the cross talk of non-coding RNAs. The interaction between non-coding RNAs based on the ceRNA can regulate the post-transcription levels of genes. According to the ceRNA hypothesis, miRNA response elements act as letters of “RNA language” to communicate with mRNAs and lncRNAs to regulate their expression. RNAs crosstalk is efficient when they share numerous miRNA response elements [43]. In this study, we identified the key DEGs and DE lncRNAs associated with ferroptosis in APA. The miRWalk 2.0 database predicted 132 miRNAs that targeted four hub genes and regulated 26 DE lncRNAs to construct the ceRNA network. In addition, well-defined lncRNAs can be functionalized by ceRNA involved in diverse mechanisms to balance the relationship between mRNAs and shared miRNAs. KCNQ1OT1 functioned as a ceRNA and stimulated the expression of HDAC3 by competitively binding to miR-452-3p, thereby repressing ABCA1 expression and lipid accumulation in macrophage [69]. The lncRNA LINC00336 acted as a sponge to inhibit miR-6852 and regulated ferroptosis by cystathionine-β-synthase. LINC00336 inhibited ferroptosis by binding to the RNA-binding protein ELAVL1 [70], implying that lncRNAs may play multiple roles in crosstalk in the ceRNA network.

In considering new insights regarding APA treatment, we used CMAP to screen for possible candidate drugs; it was developed by CMAP team of the Broad community, and the data are publicly available [71]. In the CMAP database, drug signatures fit over one million gene expression profiles obtained from approximately 50 thousand unique perturbagens and are identified by computing algorithms, thereby making the discovery of new drugs or identification of novel uses for existing drugs more efficient. Similar to the expression change between parbendazole and osteogenic hMSCs, Andrea et al. identified a novel bone anabolic compound that caused osteogenic differentiation by altering the activity between BMP-2 and cytoskeleton [72]. CMAP analysis revealed different enrichment for CX3CR+ and CX3CR- gene sets of single lung interstitial macrophage subpopulations [73]. We uploaded 44 upregulated genes and 45 downregulated genes in APA to query CMAP, and the differential expression profile was fit to the drugs identified using the Kolmogorov-Smirnov statistical test. Fludrocortisone, a synthetic mineralocorticoid, is involved in the top ranked drugs that leads to modest hyperaldosteronism [74]. The utilization of this category of drugs may cause clinical characteristics similar to those of APA. Therefore, the bottom ranked drugs were selected to determine the targeted treatment of APA through ferroptosis, in silico.

To avoid repetitive experiments and increase the research efficiency, bioinformatics technology is a highly efficient approach for predicting novel therapeutic drugs for diseases in a low-cost manner. In our study, QL-X-138 and MK-1775 bound to the AURKA and DUOX1, respectively, with the lowest binding energies in molecular docking. As MK-1775 is known as a WEE1 inhibitor, we analyzed its binding energy as a positive control to measure the affinity of docking results. The binding energies were -9.8 kcal/mol for the QL-X-138 and AURKA complex, -9.3 kcal/mol for the MK-1775 and DUOX1 complex and close to -10.3 kcal/mol for the MK-1775 and WEE1 complex. This indicated that QL-X-138 and MK-1775 had good binding activities with their targeted proteins and could be candidate drugs for further experimental verification. In addition, according to the molecular dynamics simulation, the RMSD stability is indispensable to deduce a good binding affinity [75]. The complexes of QL-X-138 and AURKA, MK-1775 and DUOX1 revealed stable RMSD in the simulated results. Molecular dynamics simulation is frequently used to predict how biomolecules, especially proteins, will respond at the atomic level when perturbed [76]. Along with its more popular application in recent years, molecular dynamics simulation has helped to design new drugs such as opioid analgesics with reduced side effects [77] and the development of existing drugs as selective agonists compared to the traditional opioid fentanyl [78]. Drug resistance is a barrier to chemotherapy in tumor related therapies. The ceRNA network of drug-resistant constructed in non-small cell lung cancer offered a novel approach to solving this problem [79]. Given that the current drugs have severe side effect in treatment of APA, ceRNA construction could help us to find novel candidate drugs with bioinformatic analyses that have been used in cancers. In this study, AURKA and DUOX1 acted as drug targets and the regulatory axis of ceRNA could be a promising therapeutic strategy for further studies.

In the present study, we provided insights into the pathogenesis of APA at the transcriptome level, explored the relationship between ferroptosis and APA, and established a methodological reference for further confirmation of experiments. Nevertheless, the sample set size for analysis was relatively small, and a large set of APA samples would be necessary to conduct in vivo and in vitro experiments to verify our findings.

Materials and Methods

Data collection

We obtained microarray data associated with APA and AAG tissue from GEO database. The datasets for mRNA expression profiling were screened as follows: (1) Homo sapiens datasets; (2) datasets contain more than 10 APA samples; (3) number of gene was more than 10000; (4) datasets extra contain lncRNA or KCNJ5 mutation samples expression profiles. At last, datasets GSE64957 and GSE101894 including mRNA and lncRNA expression profiling were collected to analysis DEGs and DE lncRNAs. The dataset GSE64957 contained transcriptome profile of zona glomerulosa (ZG), zona fasciculata (ZF) and APA, and ZG was chosen for comparison as the AAG, since the ZG was known as the zone of aldosterone synthesis and secretion in adrenal gland [17]. Datasets GSE156931 and GSE64957 including KCNJ5 mutation and wild type (WT) samples were collected to analysis hub genes expression between KCNJ5 mutation and WT.

We obtained 263 ferroptosis-related genes from FerrDb [80] which was an integrated and up-to-date database involving ferroptosis-related genes, molecules and diseases.

Data normalization and identification of ferroptosis-related DEGs

Differential expression analysis of GSE64957 was performed on the GEO2R online tool; the original files of GSE101894 and GSE156931 were downloaded from GEO and normalized by Sangerbox tools (http://www.sangerbox.com/tool) as described before [81]. The expression profiles of APA and AAG samples were compared to identified DEGs and DE lncRNAs. The screening criteria were | log2 (fold-change) | > 0.5 and adjusted P < 0.05. We selected the most significant ones when the DEGs or DE lncRNAs were duplicated. A venn diagram of DEGs and ferroptosis-related genes was drawn by Venny 2.1 (http://bioinfogp.cnb.csic.es/tools/venny/index.html).

Volcano plot and heatmap plot analysis

To better visualize the ferroptosis-related DEGs, the volcano plot and heatmap were created by Excel and Helm software, respectively.

Enrichment analysis

Metascape (http://metascape.org) is a reliable platform for gene list enrichment analysis [36]. Based on the function, we used the Custom Expression module to verify the enriched and related neighbor pathway of ferroptosis-related DEGs. The chord plot was drawn by Sangerbox tools.

GSEA software was used to assess the correlation of the lncRNAs in a predefined set of RNA table to other genes in regulating Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Statistical significance was set at P < 0.01 and false discovery rate q < 0.25.

PPI network construction

The ferroptosis-related DEGs in APA were uploaded to Metascape tool with Custom Expression module to obtain the PPI network. MCODE algorithm was then applied to this network to identify neighborhoods where proteins were densely connected. GO enrichment analysis was applied to each MCODE network to assign “meanings” to the network component. We can identify the physical and combined cores with this algorithm.

Screening for hub ferroptosis-related DEGs

The normalized data was analyzed by WGCNA package in R-Studio 3.6.1 software [41]. The application of WGCNA could group genes into modules based on their co-expression patterns. Thus, it can establish a connection between the characteristics of the samples (such as disease status or treatment) and changed expression levels of genes. This can help researchers identify potential disease-related genes or pathways. Also, it offered a systemic insight into the regulatory network and signaling pathway which may be related to the clinical traits. So clinical traits-related modules were constructed to identify the hub genes with WGCNA [82]. Topological overlap matrix (TOM) was created by transforming adjacency matrix, so that constructed network accorded with the real biological network state. Based on the TOM dissimilarity, genes were divided into several modules. Then the scale- free network was constructed to identify the gene co-expression modules with similar expression patterns. And the modules were defined by cutting the clustering tree into branches with a dynamic tree algorithm. Different colors were given for better visualization [83]. All modules were assembled by module eigengenes, and the principal component of the modules that were figured as a synthetic gene stand for the gene expression profiles [83]. Here, we set the soft-thresholding power as 8, and minimal module size as 50. Correlation between clinical traits and module eigengenes were calculated [82]. The hub ferroptosis-related DEGs were verified from the most correlative modules eigengenes and ferroptosis-related DEGs.

Expression of hub ferroptosis-related DEGs in mutant samples

To compare the expression level of hub ferroptosis-related DEGs in KCNJ5 mutation and WT APA samples, box-plots were drawn with datasets of GSE156931 and GSE64957 by Sangerbox tools. The datasets GSE156931 and GSE64957 were integrated with R package Siliconmerging and corrected of batch effect with ComBat in R package SVA to compare the expression of eight hub genes [81]. The principal component analysis (PCA) was used with R package Stats.

Correlation analysis between hub ferroptosis-related DEGs and DE lncRNAs

The relationship between hub ferroptosis-related DEGs and DE lncRNAs was assessed using Pearson correlation. The association was considered significant if the correlation coefficient | R2 | > 0.5 at P < 0.001. The correlative plot was drawn with the top 10 most correlative and highest connecting degree DE lncRNAs and hub ferroptosis-related DEGs.

Prediction of target miRNAs

miRDB, PITA, RNA22, miRWalk, miRanda, RNAhybrid, and Targetscan databases in miRWalk 2.0 online tools [84] were used to predict the target miRNA of hub ferroptosis-related DEGs. We selected miRNAs that were found in at least five databases. LncBase online tool [85] was used to predict miRNAs of DE lncRNAs to identify the lncRNAs with the highest degree of connection.

Construction of ceRNA network

The intersections of predicted miRNAs by hub ferroptosis-related DEGs and DE lncRNAs were selected to construct lncRNA-miRNA-mRNA ceRNA network. Based on the ceRNA network theory [43], the lncRNA positively correlated to the mRNA would be selected. The lncLocator was used as a subcellular localization predictor to further screen lncRNA [86]. Cytoscape software was used to depict network plot and found the mRNAs and lncRNAs with the highest degree of connection.

CMAP analysis for candidate drugs

CMAP (https://clue.io/) is a database that explores the relationships between cell physiology, diseases and therapies, and it contains over one million transcriptional expression profiles perturbed by small molecular compounds [71]. Based on this database, we can find the positive and negative drugs associated with diseases which contains a similar expression profile. We uploaded the ferroptosis-related DEGs consisted with 44 upregulated genes and 45 downregulated genes. The bottom ranked drugs resulted in opposite expression profiles compared to those by diseases, and might have a therapeutic effect.

Molecular docking

The protein 3D structure of eight hub ferroptosis-related DEGs were downloaded from protein data bank (https://www.rcsb.org/). And the structures of drugs targeted to these proteins were downloaded from ZINC (http://zinc.docking.org/). We first removed ligand and water macromolecule, and added hydrogen atoms with Pymol 2.5. The proteins were set to rigid and determined as root. Second, we set twisted key to flexible saved as pdbqt file format by AutoDock Tools 1.5.6. At last, Autodock Vina 1.1.2 was used to computing the molecular docking. The conformations with lowest binding energy were selected for further analysis of their interactions and binding patterns through Pymol 2.5 and Ligplot 2.2.

Molecular dynamics simulation

Following the molecular docking, we tested the stability of the drugs binding to the proteins with lowest binding energy. The force field for ligand and the partial charges were obtained by AMBER 18 with the antechamber module. The standard AMBER force field (ff99SB) [87] for protein and the general AMBER force field (GAFF) [88] for ligand were used in the molecular dynamics simulation by GROMACS 5.1.2 software. Each complex system was immersed in a rectangular box of TIP3P with water molecules, which expanded 12 Å away from any solute atoms. Next, we added proper numbers of Na+ or Cl- to neutralize the complex systems. After minimizing energy, each complex system was heated in from 0 to 300K within 300 ps. Under stationary temperature of 300K for complex, a 50ns molecular dynamics simulation with a step of 2 fs was performed. In the process of sample, the coordinates were saved every 2 ps for subsequently free energy calculation and free energy decomposition analysis.

Statistical analysis

Student's t-test was used to compare the gene expression of KCNJ5 mutation and WT APA samples. Pearson correlation was used to analysis the relationship between mRNAs and lncRNAs. These were performed on the Sangerbox tools.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation. Some codes are available on GitHub (https://github.com/aldosteronelab/WGCNA).

Supplementary Materials

Supplementary Figure 1

Abbreviations

APA: aldosterone-producing adenoma; AAG: adjacent adrenal gland; GEO: Gene Expression Omnibus; DEGs: differentially expressed genes; PPI: protein–protein interaction; ceRNA: competitive endogenous RNA; NOX: NADPH oxidases; GSEA: Gene Set Enrichment Analysis; GO BP: Gene Ontology biological processes; HIF1: hypoxia-inducible factor 1; SBP: systolic blood pressure; DBP: diastolic blood pressure; ALD: plasma aldosterone; RMSD: root-mean-square deviation; FDF: ferroptosis downregulated factors; ZG: zona glomerulosa; ZF: zona fasciculata; WT: wild type; KEGG: Kyoto Encyclopedia of Genes and Genomes; WGCNA: weighted gene co-expression network analysis; TOM: topological overlap matrix; PCA: principal component analysis.

Author Contributions

HY conceptualised the study, analyzed and interpreted all data, and wrote the manuscript. BY collected and assembled data. MJ analyzed the data. XS conceptualised the study, interpreted study data and revised the manuscript. All authors contributed to the article and approved the submitted version.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Ethical Statement

As GEO is a public database, the present study exempts the approval of the local ethics committee and follows GEO access policies and publication guidelines.

Funding

This study was supported by Science and technology of Zhejiang Province Project (LGF21H070003), Zhejiang Provincial Medical and Health Technology Project (2020380946, 2022502078). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

  • 1. Monticone S, Burrello J, Tizzani D, Bertello C, Viola A, Buffolo F, Gabetti L, Mengozzi G, Williams TA, Rabbia F, Veglio F, Mulatero P. Prevalence and Clinical Manifestations of Primary Aldosteronism Encountered in Primary Care Practice. J Am Coll Cardiol. 2017; 69:1811–20. https://doi.org/10.1016/j.jacc.2017.01.052 [PubMed]
  • 2. NCD Risk Factor Collaboration (NCD-RisC). Worldwide trends in blood pressure from 1975 to 2015: a pooled analysis of 1479 population-based measurement studies with 19·1 million participants. Lancet. 2017; 389:37–55. https://doi.org/10.1016/S0140-6736(16)31919-5 [PubMed]
  • 3. Calhoun DA. Hyperaldosteronism as a common cause of resistant hypertension. Annu Rev Med. 2013; 64:233–47. https://doi.org/10.1146/annurev-med-042711-135929 [PubMed]
  • 4. Funder JW, Carey RM, Mantero F, Murad MH, Reincke M, Shibata H, Stowasser M, Young WF Jr. The Management of Primary Aldosteronism: Case Detection, Diagnosis, and Treatment: An Endocrine Society Clinical Practice Guideline. J Clin Endocrinol Metab. 2016; 101:1889–916. https://doi.org/10.1210/jc.2015-4061 [PubMed]
  • 5. Fernandes-Rosa FL, Daniil G, Orozco IJ, Göppner C, El Zein R, Jain V, Boulkroun S, Jeunemaitre X, Amar L, Lefebvre H, Schwarzmayr T, Strom TM, Jentsch TJ, Zennaro MC. A gain-of-function mutation in the CLCN2 chloride channel gene causes primary aldosteronism. Nat Genet. 2018; 50:355–61. https://doi.org/10.1038/s41588-018-0053-8 [PubMed]
  • 6. Zennaro MC, Boulkroun S, Fernandes-Rosa F. Genetic Causes of Functional Adrenocortical Adenomas. Endocr Rev. 2017; 38:516–37. https://doi.org/10.1210/er.2017-00189 [PubMed]
  • 7. Ye P, Mariniello B, Mantero F, Shibata H, Rainey WE. G-protein-coupled receptors in aldosterone-producing adenomas: a potential cause of hyperaldosteronism. J Endocrinol. 2007; 195:39–48. https://doi.org/10.1677/JOE-07-0037 [PubMed]
  • 8. Choi M, Scholl UI, Yue P, Björklund P, Zhao B, Nelson-Williams C, Ji W, Cho Y, Patel A, Men CJ, Lolis E, Wisgerhof MV, Geller DS, et al. K+ channel mutations in adrenal aldosterone-producing adenomas and hereditary hypertension. Science. 2011; 331:768–72. https://doi.org/10.1126/science.1198785 [PubMed]
  • 9. Azizan EA, Poulsen H, Tuluc P, Zhou J, Clausen MV, Lieb A, Maniero C, Garg S, Bochukova EG, Zhao W, Shaikh LH, Brighton CA, Teo AE, et al. Somatic mutations in ATP1A1 and CACNA1D underlie a common subtype of adrenal hypertension. Nat Genet. 2013; 45:1055–60. https://doi.org/10.1038/ng.2716 [PubMed]
  • 10. Scholl UI, Goh G, Stölting G, de Oliveira RC, Choi M, Overton JD, Fonseca AL, Korah R, Starker LF, Kunstman JW, Prasad ML, Hartung EA, Mauras N, et al. Somatic and germline CACNA1D calcium channel mutations in aldosterone-producing adenomas and primary aldosteronism. Nat Genet. 2013; 45:1050–4. https://doi.org/10.1038/ng.2695 [PubMed]
  • 11. Nanba K, Blinder AR, Rege J, Hattangady NG, Else T, Liu CJ, Tomlins SA, Vats P, Kumar-Sinha C, Giordano TJ, Rainey WE. Somatic CACNA1H Mutation As a Cause of Aldosterone-Producing Adenoma. Hypertension. 2020; 75:645–9. https://doi.org/10.1161/HYPERTENSIONAHA.119.14349 [PubMed]
  • 12. Scholl UI, Stölting G, Schewe J, Thiel A, Tan H, Nelson-Williams C, Vichot AA, Jin SC, Loring E, Untiet V, Yoo T, Choi J, Xu S, et al. CLCN2 chloride channel mutations in familial hyperaldosteronism type II. Nat Genet. 2018; 50:349–54. https://doi.org/10.1038/s41588-018-0048-5 [PubMed]
  • 13. Beuschlein F, Boulkroun S, Osswald A, Wieland T, Nielsen HN, Lichtenauer UD, Penton D, Schack VR, Amar L, Fischer E, Walther A, Tauber P, Schwarzmayr T, et al. Somatic mutations in ATP1A1 and ATP2B3 lead to aldosterone-producing adenomas and secondary hypertension. Nat Genet. 2013; 45:440–4. https://doi.org/10.1038/ng.2550 [PubMed]
  • 14. Oki K, Plonczynski MW, Luis Lam M, Gomez-Sanchez EP, Gomez-Sanchez CE. Potassium channel mutant KCNJ5 T158A expression in HAC-15 cells increases aldosterone synthesis. Endocrinology. 2012; 153:1774–82. https://doi.org/10.1210/en.2011-1733 [PubMed]
  • 15. Yang Y, Tetti M, Vohra T, Adolf C, Seissler J, Hristov M, Belavgeni A, Bidlingmaier M, Linkermann A, Mulatero P, Beuschlein F, Reincke M, Williams TA. BEX1 Is Differentially Expressed in Aldosterone-Producing Adenomas and Protects Human Adrenocortical Cells From Ferroptosis. Hypertension. 2021; 77:1647–58. https://doi.org/10.1161/HYPERTENSIONAHA.120.16774 [PubMed]
  • 16. Williams TA, Lenders JWM, Mulatero P, Burrello J, Rottenkolber M, Adolf C, Satoh F, Amar L, Quinkler M, Deinum J, Beuschlein F, Kitamoto KK, Pham U, et al, and Primary Aldosteronism Surgery Outcome (PASO) investigators. Outcomes after adrenalectomy for unilateral primary aldosteronism: an international consensus on outcome measures and analysis of remission rates in an international cohort. Lancet Diabetes Endocrinol. 2017; 5:689–99. https://doi.org/10.1016/S2213-8587(17)30135-3 [PubMed]
  • 17. Shaikh LH, Zhou J, Teo AE, Garg S, Neogi SG, Figg N, Yeo GS, Yu H, Maguire JJ, Zhao W, Bennett MR, Azizan EA, Davenport AP, et al. LGR5 Activates Noncanonical Wnt Signaling and Inhibits Aldosterone Production in the Human Adrenal. J Clin Endocrinol Metab. 2015; 100:E836–44. https://doi.org/10.1210/jc.2015-1734 [PubMed]
  • 18. Moujalled D, Strasser A, Liddell JR. Molecular mechanisms of cell death in neurological diseases. Cell Death Differ. 2021; 28:2029–44. https://doi.org/10.1038/s41418-021-00814-y [PubMed]
  • 19. Dixon SJ, Lemberg KM, Lamprecht MR, Skouta R, Zaitsev EM, Gleason CE, Patel DN, Bauer AJ, Cantley AM, Yang WS, Morrison B 3rd, Stockwell BR. Ferroptosis: an iron-dependent form of nonapoptotic cell death. Cell. 2012; 149:1060–72. https://doi.org/10.1016/j.cell.2012.03.042 [PubMed]
  • 20. Tang Y, Li C, Zhang YJ, Wu ZH. Ferroptosis-Related Long Non-Coding RNA signature predicts the prognosis of Head and neck squamous cell carcinoma. Int J Biol Sci. 2021; 17:702–11. https://doi.org/10.7150/ijbs.55552 [PubMed]
  • 21. Xie Y, Hou W, Song X, Yu Y, Huang J, Sun X, Kang R, Tang D. Ferroptosis: process and function. Cell Death Differ. 2016; 23:369–79. https://doi.org/10.1038/cdd.2015.158 [PubMed]
  • 22. Stockwell BR, Friedmann Angeli JP, Bayir H, Bush AI, Conrad M, Dixon SJ, Fulda S, Gascón S, Hatzios SK, Kagan VE, Noel K, Jiang X, Linkermann A, et al. Ferroptosis: A Regulated Cell Death Nexus Linking Metabolism, Redox Biology, and Disease. Cell. 2017; 171:273–85. https://doi.org/10.1016/j.cell.2017.09.021 [PubMed]
  • 23. Yang WH, Ding CC, Sun T, Rupprecht G, Lin CC, Hsu D, Chi JT. The Hippo Pathway Effector TAZ Regulates Ferroptosis in Renal Cell Carcinoma. Cell Rep. 2019; 28:2501–8.e4. https://doi.org/10.1016/j.celrep.2019.07.107 [PubMed]
  • 24. Chen X, Kang R, Kroemer G, Tang D. Broadening horizons: the role of ferroptosis in cancer. Nat Rev Clin Oncol. 2021; 18:280–96. https://doi.org/10.1038/s41571-020-00462-0 [PubMed]
  • 25. Gorrini C, Harris IS, Mak TW. Modulation of oxidative stress as an anticancer strategy. Nat Rev Drug Discov. 2013; 12:931–47. https://doi.org/10.1038/nrd4002 [PubMed]
  • 26. Louandre C, Ezzoukhry Z, Godin C, Barbare JC, Mazière JC, Chauffert B, Galmiche A. Iron-dependent cell death of hepatocellular carcinoma cells exposed to sorafenib. Int J Cancer. 2013; 133:1732–42. https://doi.org/10.1002/ijc.28159 [PubMed]
  • 27. Wang W, Green M, Choi JE, Gijón M, Kennedy PD, Johnson JK, Liao P, Lang X, Kryczek I, Sell A, Xia H, Zhou J, Li G, et al. CD8+ T cells regulate tumour ferroptosis during cancer immunotherapy. Nature. 2019; 569:270–4. https://doi.org/10.1038/s41586-019-1170-y [PubMed]
  • 28. Lei G, Zhang Y, Koppula P, Liu X, Zhang J, Lin SH, Ajani JA, Xiao Q, Liao Z, Wang H, Gan B. The role of ferroptosis in ionizing radiation-induced cell death and tumor suppression. Cell Res. 2020; 30:146–62. https://doi.org/10.1038/s41422-019-0263-3 [PubMed]
  • 29. Liang C, Zhang X, Yang M, Dong X. Recent Progress in Ferroptosis Inducers for Cancer Therapy. Adv Mater. 2019; 31:e1904197. https://doi.org/10.1002/adma.201904197 [PubMed]
  • 30. Tang M, Chen Z, Wu D, Chen L. Ferritinophagy/ferroptosis: Iron-related newcomers in human diseases. J Cell Physiol. 2018; 233:9179–90. https://doi.org/10.1002/jcp.26954 [PubMed]
  • 31. Gong S, Sun N, Meyer LS, Tetti M, Koupourtidou C, Krebs S, Masserdotti G, Blum H, Rainey WE, Reincke M, Walch A, Williams TA. Primary Aldosteronism: Spatial Multiomics Mapping of Genotype-Dependent Heterogeneity and Tumor Expansion of Aldosterone-Producing Adenomas. Hypertension. 2023; 80:1555–67. https://doi.org/10.1161/HYPERTENSIONAHA.123.20921 [PubMed]
  • 32. Wang Z, Chen X, Liu N, Shi Y, Liu Y, Ouyang L, Tam S, Xiao D, Liu S, Wen F, Tao Y. A Nuclear Long Non-Coding RNA LINC00618 Accelerates Ferroptosis in a Manner Dependent upon Apoptosis. Mol Ther. 2021; 29:263–74. https://doi.org/10.1016/j.ymthe.2020.09.024 [PubMed]
  • 33. Yang L, Wang S, Zhang Q, Pan Y, Lv Y, Chen X, Zuo Y, Hao D. Clinical significance of the immune microenvironment in ovarian cancer patients. Mol Omics. 2018; 14:341–51. https://doi.org/10.1039/c8mo00128f [PubMed]
  • 34. Wang S, Zhang Q, Yu C, Cao Y, Zuo Y, Yang L. Immune cell infiltration-based signature for prognosis and immunogenomic analysis in breast cancer. Brief Bioinform. 2021; 22:2020–31. https://doi.org/10.1093/bib/bbaa026 [PubMed]
  • 35. Wang S, Xiong Y, Zhang Q, Su D, Yu C, Cao Y, Pan Y, Lu Q, Zuo Y, Yang L. Clinical significance and immunogenomic landscape analyses of the immune cell signature based prognostic model for patients with breast cancer. Brief Bioinform. 2021; 22:bbaa311. https://doi.org/10.1093/bib/bbaa311 [PubMed]
  • 36. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, Benner C, Chanda SK. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019; 10:1523. https://doi.org/10.1038/s41467-019-09234-6 [PubMed]
  • 37. Chen YH, Chen ZW, Li HM, Yan XF, Feng B. AGE/RAGE-Induced EMP Release via the NOX-Derived ROS Pathway. J Diabetes Res. 2018; 2018:6823058. https://doi.org/10.1155/2018/6823058 [PubMed]
  • 38. Pathomthongtaweechai N, Chutipongtanate S. AGE/RAGE signaling-mediated endoplasmic reticulum stress and future prospects in non-coding RNA therapeutics for diabetic nephropathy. Biomed Pharmacother. 2020; 131:110655. https://doi.org/10.1016/j.biopha.2020.110655 [PubMed]
  • 39. Su LJ, Zhang JH, Gomez H, Murugan R, Hong X, Xu D, Jiang F, Peng ZY. Reactive Oxygen Species-Induced Lipid Peroxidation in Apoptosis, Autophagy, and Ferroptosis. Oxid Med Cell Longev. 2019; 2019:5080843. https://doi.org/10.1155/2019/5080843 [PubMed]
  • 40. Ye Z, Zhuo Q, Hu Q, Xu X, Mengqi Li, Zhang Z, Xu W, Liu W, Fan G, Qin Y, Yu X, Ji S. FBW7-NRA41-SCD1 axis synchronously regulates apoptosis and ferroptosis in pancreatic cancer cells. Redox Biol. 2021; 38:101807. https://doi.org/10.1016/j.redox.2020.101807 [PubMed]
  • 41. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008; 9:559. https://doi.org/10.1186/1471-2105-9-559 [PubMed]
  • 42. Lenzini L, Rossitto G, Maiolino G, Letizia C, Funder JW, Rossi GP. A Meta-Analysis of Somatic KCNJ5 K(+) Channel Mutations In 1636 Patients With an Aldosterone-Producing Adenoma. J Clin Endocrinol Metab. 2015; 100:E1089–95. https://doi.org/10.1210/jc.2015-2149 [PubMed]
  • 43. Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell. 2011; 146:353–8. https://doi.org/10.1016/j.cell.2011.07.014 [PubMed]
  • 44. Alim I, Caulfield JT, Chen Y, Swarup V, Geschwind DH, Ivanova E, Seravalli J, Ai Y, Sansing LH, Ste Marie EJ, Hondal RJ, Mukherjee S, Cave JW, et al. Selenium Drives a Transcriptional Adaptive Program to Block Ferroptosis and Treat Stroke. Cell. 2019; 177:1262–79.e25. https://doi.org/10.1016/j.cell.2019.03.032 [PubMed]
  • 45. Li D, Jiang C, Mei G, Zhao Y, Chen L, Liu J, Tang Y, Gao C, Yao P. Quercetin Alleviates Ferroptosis of Pancreatic β Cells in Type 2 Diabetes. Nutrients. 2020; 12:2954. https://doi.org/10.3390/nu12102954 [PubMed]
  • 46. Tang Z, Ju Y, Dai X, Ni N, Liu Y, Zhang D, Gao H, Sun H, Zhang J, Gu P. HO-1-mediated ferroptosis as a target for protection against retinal pigment epithelium degeneration. Redox Biol. 2021; 43:101971. https://doi.org/10.1016/j.redox.2021.101971 [PubMed]
  • 47. Li J, Cao F, Yin HL, Huang ZJ, Lin ZT, Mao N, Sun B, Wang G. Ferroptosis: past, present and future. Cell Death Dis. 2020; 11:88. https://doi.org/10.1038/s41419-020-2298-2 [PubMed]
  • 48. Li X, Zou Y, Xing J, Fu YY, Wang KY, Wan PZ, Zhai XY. Pretreatment with Roxadustat (FG-4592) Attenuates Folic Acid-Induced Kidney Injury through Antiferroptosis via Akt/GSK-3β/Nrf2 Pathway. Oxid Med Cell Longev. 2020; 2020:6286984. https://doi.org/10.1155/2020/6286984 [PubMed]
  • 49. Feng X, Wang S, Sun Z, Dong H, Yu H, Huang M, Gao X. Ferroptosis Enhanced Diabetic Renal Tubular Injury via HIF-1α/HO-1 Pathway in db/db Mice. Front Endocrinol (Lausanne). 2021; 12:626390. https://doi.org/10.3389/fendo.2021.626390 [PubMed]
  • 50. Bedard K, Krause KH. The NOX family of ROS-generating NADPH oxidases: physiology and pathophysiology. Physiol Rev. 2007; 87:245–313. https://doi.org/10.1152/physrev.00044.2005 [PubMed]
  • 51. Park MW, Cha HW, Kim J, Kim JH, Yang H, Yoon S, Boonpraman N, Yi SS, Yoo ID, Moon JS. NOX4 promotes ferroptosis of astrocytes by oxidative stress-induced lipid peroxidation via the impairment of mitochondrial metabolism in Alzheimer’s diseases. Redox Biol. 2021; 41:101947. https://doi.org/10.1016/j.redox.2021.101947 [PubMed]
  • 52. Park E, Chung SW. ROS-mediated autophagy increases intracellular iron levels and ferroptosis by ferritin and transferrin receptor regulation. Cell Death Dis. 2019; 10:822. https://doi.org/10.1038/s41419-019-2064-5 [PubMed]
  • 53. Otto T, Sicinski P. Cell cycle proteins as promising targets in cancer therapy. Nat Rev Cancer. 2017; 17:93–115. https://doi.org/10.1038/nrc.2016.138 [PubMed]
  • 54. Gomaa A, Peng D, Chen Z, Soutto M, Abouelezz K, Corvalan A, El-Rifai W. Epigenetic regulation of AURKA by miR-4715-3p in upper gastrointestinal cancers. Sci Rep. 2019; 9:16970. https://doi.org/10.1038/s41598-019-53174-6 [PubMed]
  • 55. Yang WS, Stockwell BR. Ferroptosis: Death by Lipid Peroxidation. Trends Cell Biol. 2016; 26:165–76. https://doi.org/10.1016/j.tcb.2015.10.014 [PubMed]
  • 56. Yang L, Wang H, Yang X, Wu Q, An P, Jin X, Liu W, Huang X, Li Y, Yan S, Shen S, Liang T, Min J, Wang F. Auranofin mitigates systemic iron overload and induces ferroptosis via distinct mechanisms. Signal Transduct Target Ther. 2020; 5:138. https://doi.org/10.1038/s41392-020-00253-0 [PubMed]
  • 57. Bin S, Xin L, Lin Z, Jinhua Z, Rui G, Xiang Z. Targeting miR-10a-5p/IL-6R axis for reducing IL-6-induced cartilage cell ferroptosis. Exp Mol Pathol. 2021; 118:104570. https://doi.org/10.1016/j.yexmp.2020.104570 [PubMed]
  • 58. Zhang X, Du L, Qiao Y, Zhang X, Zheng W, Wu Q, Chen Y, Zhu G, Liu Y, Bian Z, Guo S, Yang Y, Ma L, et al. Ferroptosis is governed by differential regulation of transcription in liver cancer. Redox Biol. 2019; 24:101211. https://doi.org/10.1016/j.redox.2019.101211 [PubMed]
  • 59. Shen C, Kong B, Liu Y, Xiong L, Shuai W, Wang G, Quan D, Huang H. YY1-induced upregulation of lncRNA KCNQ1OT1 regulates angiotensin II-induced atrial fibrillation by modulating miR-384b/CACNA1C axis. Biochem Biophys Res Commun. 2018; 505:134–40. https://doi.org/10.1016/j.bbrc.2018.09.064 [PubMed]
  • 60. Sallé-Lefort S, Miard S, Nolin MA, Boivin L, Paré MÈ, Debigaré R, Picard F. Hypoxia upregulates Malat1 expression through a CaMKK/AMPK/HIF-1α axis. Int J Oncol. 2016; 49:1731–6. https://doi.org/10.3892/ijo.2016.3630 [PubMed]
  • 61. Zennaro MC, Boulkroun S, Fernandes-Rosa FL. Pathogenesis and treatment of primary aldosteronism. Nat Rev Endocrinol. 2020; 16:578–89. https://doi.org/10.1038/s41574-020-0382-4 [PubMed]
  • 62. Li JQ, Hu SY, Wang ZY, Lin J, Jian S, Dong YC, Wu XF, Da, Cao LJ. Long non-coding RNA MEG3 inhibits microRNA-125a-5p expression and induces immune imbalance of Treg/Th17 in immune thrombocytopenic purpura. Biomed Pharmacother. 2016; 83:905–11. https://doi.org/10.1016/j.biopha.2016.07.057 [PubMed]
  • 63. Wang L, Xia JW, Ke ZP, Zhang BH. Blockade of NEAT1 represses inflammation response and lipid uptake via modulating miR-342-3p in human macrophages THP-1 cells. J Cell Physiol. 2019; 234:5319–26. https://doi.org/10.1002/jcp.27340 [PubMed]
  • 64. Boulkroun S, Beuschlein F, Rossi GP, Golib-Dzib JF, Fischer E, Amar L, Mulatero P, Samson-Couterie B, Hahner S, Quinkler M, Fallo F, Letizia C, Allolio B, et al. Prevalence, clinical, and molecular correlates of KCNJ5 mutations in primary aldosteronism. Hypertension. 2012; 59:592–8. https://doi.org/10.1161/HYPERTENSIONAHA.111.186478 [PubMed]
  • 65. Nanba K, Yamazaki Y, Bick N, Onodera K, Tezuka Y, Omata K, Ono Y, Blinder AR, Tomlins SA, Rainey WE, Satoh F, Sasano H. Prevalence of Somatic Mutations in Aldosterone-Producing Adenomas in Japanese Patients. J Clin Endocrinol Metab. 2020; 105:e4066–73. https://doi.org/10.1210/clinem/dgaa595 [PubMed]
  • 66. Zhang C, Wu L, Jiang L, Su T, Zhou W, Zhong X, Xie J, Sun F, Zhu Y, Jiang Y, Wang W. KCNJ5 Mutation Contributes to Complete Clinical Success in Aldosterone-Producing Adenoma: A Study From a Single Center. Endocr Pract. 2021; 27:736–42. https://doi.org/10.1016/j.eprac.2021.01.007 [PubMed]
  • 67. Guo XW, Zhang H, Huang JQ, Wang SN, Lu Y, Cheng B, Dong SH, Wang YY, Li FS, Li YW. PIEZO1 Ion Channel Mediates Ionizing Radiation-Induced Pulmonary Endothelial Cell Ferroptosis via Ca2+/Calpain/VE-Cadherin Signaling. Front Mol Biosci. 2021; 8:725274. https://doi.org/10.3389/fmolb.2021.725274 [PubMed]
  • 68. Lu G, Xu C, Tang K, Zhang J, Li Q, Peng L, Wang Y, Huang Z, Gao X. H2S inhibits angiotensin II-induced atrial Kv1.5 upregulation by attenuating Nox4-mediated ROS generation during atrial fibrillation. Biochem Biophys Res Commun. 2017; 483:534–40. https://doi.org/10.1016/j.bbrc.2016.12.110 [PubMed]
  • 69. Yu XH, Deng WY, Chen JJ, Xu XD, Liu XX, Chen L, Shi MW, Liu QX, Tao M, Ren K. LncRNA kcnq1ot1 promotes lipid accumulation and accelerates atherosclerosis via functioning as a ceRNA through the miR-452-3p/HDAC3/ABCA1 axis. Cell Death Dis. 2020; 11:1043. https://doi.org/10.1038/s41419-020-03263-6 [PubMed]
  • 70. Wang M, Mao C, Ouyang L, Liu Y, Lai W, Liu N, Shi Y, Chen L, Xiao D, Yu F, Wang X, Zhou H, Cao Y, et al. Long noncoding RNA LINC00336 inhibits ferroptosis in lung cancer by functioning as a competing endogenous RNA. Cell Death Differ. 2019; 26:2329–43. https://doi.org/10.1038/s41418-019-0304-y [PubMed]
  • 71. Subramanian A, Narayan R, Corsello SM, Peck DD, Natoli TE, Lu X, Gould J, Davis JF, Tubelli AA, Asiedu JK, Lahr DL, Hirschman JE, Liu Z, et al. A Next Generation Connectivity Map: L1000 Platform and the First 1,000,000 Profiles. Cell. 2017; 171:1437–52.e17. https://doi.org/10.1016/j.cell.2017.10.049 [PubMed]
  • 72. Brum AM, van de Peppel J, van der Leije CS, Schreuders-Koedam M, Eijken M, van der Eerden BC, van Leeuwen JP. Connectivity Map-based discovery of parbendazole reveals targetable human osteogenic pathway. Proc Natl Acad Sci USA. 2015; 112:12711–6. https://doi.org/10.1073/pnas.1501597112 [PubMed]
  • 73. Chakarov S, Lim HY, Tan L, Lim SY, See P, Lum J, Zhang XM, Foo S, Nakamizo S, Duan K, Kong WT, Gentek R, Balachander A, et al. Two distinct interstitial macrophage populations coexist across tissues in specific subtissular niches. Science. 2019; 363:eaau0964. https://doi.org/10.1126/science.aau0964 [PubMed]
  • 74. Tu W, Eckert GJ, Hannon TS, Liu H, Pratt LM, Wagner MA, Dimeglio LA, Jung J, Pratt JH. Racial differences in sensitivity of blood pressure to aldosterone. Hypertension. 2014; 63:1212–8. https://doi.org/10.1161/HYPERTENSIONAHA.113.02989 [PubMed]
  • 75. Doniach S, Eastman P. Protein dynamics simulations from nanoseconds to microseconds. Curr Opin Struct Biol. 1999; 9:157–63. https://doi.org/10.1016/S0959-440X(99)80022-0 [PubMed]
  • 76. Hollingsworth SA, Dror RO. Molecular Dynamics Simulation for All. Neuron. 2018; 99:1129–43. https://doi.org/10.1016/j.neuron.2018.08.011 [PubMed]
  • 77. Manglik A, Lin H, Aryal DK, McCorvy JD, Dengler D, Corder G, Levit A, Kling RC, Bernat V, Hübner H, Huang XP, Sassano MF, Giguère PM, et al. Structure-based discovery of opioid analgesics with reduced side effects. Nature. 2016; 537:185–90. https://doi.org/10.1038/nature19112 [PubMed]
  • 78. Spahn V, Del Vecchio G, Labuz D, Rodriguez-Gaztelumendi A, Massaly N, Temp J, Durmaz V, Sabri P, Reidelbach M, Machelska H, Weber M, Stein C. A nontoxic pain killer designed by modeling of pathological receptor conformations. Science. 2017; 355:966–9. https://doi.org/10.1126/science.aai8636 [PubMed]
  • 79. Kong X, Hu S, Yuan Y, Du Y, Zhu Z, Song Z, Lu S, Zhao C, Yan D. Analysis of lncRNA, miRNA and mRNA-associated ceRNA networks and identification of potential drug targets for drug-resistant non-small cell lung cancer. J Cancer. 2020; 11:3357–68. https://doi.org/10.7150/jca.40729 [PubMed]
  • 80. Zhou N, Bao J. FerrDb: a manually curated resource for regulators and markers of ferroptosis and ferroptosis-disease associations. Database (Oxford). 2020; 2020:baaa021. https://doi.org/10.1093/database/baaa021 [PubMed]
  • 81. Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007; 8:118–27. https://doi.org/10.1093/biostatistics/kxj037 [PubMed]
  • 82. Hu X, Bao M, Huang J, Zhou L, Zheng S. Identification and Validation of Novel Biomarkers for Diagnosis and Prognosis of Hepatocellular Carcinoma. Front Oncol. 2020; 10:541479. https://doi.org/10.3389/fonc.2020.541479 [PubMed] Erratum in: Front Oncol. 2020; 10:617539. https://doi.org/10.3389/fonc.2020.617539 [PubMed]
  • 83. Wang M, Wang L, Pu L, Li K, Feng T, Zheng P, Li S, Sun M, Yao Y, Jin L. LncRNAs related key pathways and genes in ischemic stroke by weighted gene co-expression network analysis (WGCNA). Genomics. 2020; 112:2302–8. https://doi.org/10.1016/j.ygeno.2020.01.001 [PubMed]
  • 84. Dweep H, Gretz N. miRWalk2.0: a comprehensive atlas of microRNA-target interactions. Nat Methods. 2015; 12:697. https://doi.org/10.1038/nmeth.3485 [PubMed]
  • 85. Paraskevopoulou MD, Vlachos IS, Karagkouni D, Georgakilas G, Kanellos I, Vergoulis T, Zagganas K, Tsanakas P, Floros E, Dalamagas T, Hatzigeorgiou AG. DIANA-LncBase v2: indexing microRNA targets on non-coding transcripts. Nucleic Acids Res. 2016; 44:D231–8. https://doi.org/10.1093/nar/gkv1270 [PubMed]
  • 86. Cao Z, Pan X, Yang Y, Huang Y, Shen HB. The lncLocator: a subcellular localization predictor for long non-coding RNAs based on a stacked ensemble classifier. Bioinformatics. 2018; 34:2185–94. https://doi.org/10.1093/bioinformatics/bty085 [PubMed]
  • 87. Duan Y, Wu C, Chowdhury S, Lee MC, Xiong G, Zhang W, Yang R, Cieplak P, Luo R, Lee T, Caldwell J, Wang J, Kollman P. A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations. J Comput Chem. 2003; 24:1999–2012. https://doi.org/10.1002/jcc.10349 [PubMed]
  • 88. Wang J, Wolf RM, Caldwell JW, Kollman PA, Case DA. Development and testing of a general amber force field. J Comput Chem. 2004; 25:1157–74. https://doi.org/10.1002/jcc.20035 [PubMed]
Download PDF
Corresponding Author

Song Xiaoxiao
xsong103@zju.edu.cn

Keywords
aldosterone-producing adenomas ferroptosis competing endogenous RNA candidate drugs differentially expressed genes
Paper Sections