Ubiquitin-specific protease as the underlying gene biomarker for aortic stenosis

Background Aortic stenosis is a common heart valvular disease whose pathological processes include an inflammatory reaction and lipid accumulation. However, its detailed pathogenesis is yet to be completely elucidated. Therefore, it is of great significance to further explore the molecular mechanisms of aortic stenosis. Methods Four datasets were downloaded from the Gene Expression Omnibus (GEO) database. Firstly, the differently expressed genes (DEGs) were screened between control and aortic stenosis samples. Secondly, weighted gene co-expression network analysis (WGCNA) was performed to find the highly relevant gene modules. Enrichment analysis and protein-protein interaction (PPI) networking were also performed, then Cytoscape was used to identify hub genes. Finally, the six participants (3 control participants and 3 patients with aortic stenosis) were recruited at the Tianjin Chest Hospital. In order to verify the expression level of USP14, several molecular experiments were performed, including hematoxylin-eosin (HE) staining, immunohistochemistry, immunofluorescence technology, real time-quantitative polymerase chain reaction (RT-qPCR), and western blotting. Results A total of 9636 DEGs were found between the control and aortic stenosis samples. The DEGs were mainly enriched in the autophagy-animal, cellular lipid catabolic process, apoptosis, and glycoside metabolic process categories. Eleven hub genes were identified via four different algorithms. Following verification of the patient samples, Ubiquitin-specific protease 14 (USP14) was found to be displayed at higher levels in the aortic stenosis samples. Conclusion USP14 might be involved in the occurrence and development of aortic stenosis, so it would be a molecular target for early diagnosis and specific treatment of aortic stenosis. There is a significant association between the high expression of USP14 and aortic stenosis, indicating that this gene may be a genetic risk factor for aortic stenosis.


Background
Aortic stenosis (AS) is a common heart valvular disease [1]. In recent years, the incidence of AS has been on the rise. Apart from its degenerative, congenital (bicuspid valve), and inflammatory causes, AS might also be caused by rheumatic heart disease and other less common diseases, such as Paget's disease, Fabry disease, and lupus erythematosus. The diagnosis of aortic stenosis is based on an echocardiographic assessment [2]. The valve area is 3-4 cm 2 in healthy adults. In general, AS patients do not exhibit clinical symptoms until the valve area is less than 1.0 cm 2 . At that point, the extra work required to open the valve prevents the heart from functioning properly. Currently, age-related degenerative aortic stenosis has become the most common cause of aortic stenosis in adults. The prevalence of aortic stenosis increases with age, averaging 0.2% in people aged 50 to 59 and 9.8% in people aged 80 to 89 [3].
The pathological process of aortic stenosis includes an inflammatory reaction, lipid accumulation, angiotensin converting enzyme activation, and macrophage and T lymphocyte infiltration. The activity of the valve is limited due to the deposition of calcium at the base of the valve, which leads to aortic stenosis. Aortic valve calcification is similar to coronary heart disease and has a high correlation with coronary artery calcification. Hypertension, dyslipidemia, diabetes, and smoking are the risk factors [4,5].
Aortic stenosis is a slowly progressing disease. Patients with aortic stenosis have a long asymptomatic period; those without symptoms need no treatment but should follow-up regularly. However, when symptoms appear or obstruction becomes very severe, rapid diagnosis and treatment is essential [6]. At present, the treatments for aortic stenosis include surgical aortic valve replacement (AVR) and transcatheter aortic valve implantation (TAVI) [7]. In fact, sometimes patients are referred to AVR after the myocardium has started to decompensate, which might cause new tissue damage and lead to serious complications. Thus, early diagnosis and prevention are needed. In addition, TAVI is also associated with paravalvular aortic regurgitation, stroke, vascular complications, and the need for a new pacemaker.
Therefore, it is of great clinical significance to further explore the molecular mechanisms of aortic stenosis, in order to find molecular targets for early diagnosis, prevention, and specific treatment. Microarray technology can use genomic data to explore and mine diseaserelated differentially expressed genes. At present, it is an effective method to explore the pathogenesis of diseases [8]. To meet these goals, this paper uses microarray technology to mine the differentially expressed genes between stenotic valves and normal aortic valves, and to identify the hub genes of aortic stenosis. Results are further verified by molecular experiments.

Data from the GEO database
On October 16, 2019, the key words "(aortic valve) AND calcified" were set to detect the datasets, using a filter of "expression profiling by array." There were five inclusion criteria: a sample number of more than ten per dataset (samples of less than ten were excluded), data from Homo sapiens (data from other species were excluded), a series entry type, expression profiling by array (data using methylation profiling by array were excluded), and a diagnosis of aortic stenosis caused by the aortic valve sclerosis or calcification (data from mitral valve stenosis diagnoses were excluded).
The study obtained the transcriptome expression profiles GSE12644 ( (Table 1).

Screening of DEGs
The four GEO series were merged and normalized by the R packages sva and limma. The DEGs were screened by R package limma; the cut-off criteria were P value < 0.05 and log [Fold Change (FC)] ≥ 100 or ≤ − 100.

WGCNA analysis
WGCNA is an algorithm for mining module information from chip data, which can describe the patterns of genes between microarray samples and find highly relevant gene modules. In this study, WGCNA analysis was conducted by the R package WGCNA. The lowest thresholding power for the scale-free topology fit index was Construction and analysis of the PPI network STRING (http://string-db.org, Version 11.0), an online database, can predict and provide a PPI network after importing the MEdarkgrey model DEGs. For this analysis, all data sources were considered. The minimum required interaction score was medium confidence (0.400). Cytoscape (Version 3.7.2) is an analysis tool used to provide biological network analysis and two-dimensional (2D) visualization for biologists. The algorithms which were used for network analysis included a Betweenness algorithm, Radiality algorithm, EcCentricity algorithm, and Closeness algorithm. In this study, the PPI network and hub genes were constructed and analyzed by the STRING database and Cytoscape.

Identification of hub DEGs associated with cardiovascular diseases
The comparative toxicogenomics database (CTD) (http://ctdbase.org/, released November 6, 2019 [15938]) is a web-based database that can predict the relationship between a gene/protein and disease. Operation instructions are available at http://ctdbase.org/help/. In this study, the relationships between gene products and heart diseases were analyzed by this database.

MiRNA of hub DEGs prediction
TargetScan (www.targetscan.org, released March 2, 2018) is an online database that performs predictive analysis and identifies possible mechanisms for coregulating the expression of hundreds of genes expressed in different cell types. In this study, TargetScan was used to screen for miRNAs that regulate the hub DEGs.
Functional annotation of miRNA DIANA-miRPath v3.0 (http://www.microrna.gr/miR-Pathv3, version 3) is an online analysis tool suite dedicated to conducting functional and pathway enrichment analysis for miRNA. In this study, GO and KEGG pathway enrichment analysis were performed using miRPath (P < 0.05).

Statistical analysis
The analyses in this study were performed using Perl, R (

Participants and ethics
At the Tianjin Chest Hospital, the pathological sections were obtained between 2018 and 2019 from six participants. The inclusion criteria included being between 18 and 80 years old, a diagnosis of aortic stenosis, and a non-surgical history. Patients whose poor cardiac  Independent-samples T test was used to compare continuous data. Pearson's chi-squared test was used to compare the categorical data. The co-morbidities in this research included hypertension, diabetes, and coronary heart disease. "No" of "Co-morbidities" presented that the participants did not suffer from any of the above three co-morbidities. "Yes" of "Co-morbidities" presented that the participant suffered from all above three co-morbidities function, pulmonary function, or liver and kidney function made them unable to tolerate surgery were excluded. The diagnosis of aortic stenosis was based on the preoperative surface echocardiography, intraoperative esophageal echocardiography, and intraoperative findings [6]. The pathological report manifested that a total of 3 participants were diagnosed with aortic stenosis while 3 others were not. The demographic details of the participants are shown in Table 2. The research conformed to the Declaration of Helsinki, and it was authorized by the Human Ethics and Research Ethics Committees of the Tianjin Chest Hospital. Informed consent was provided by all the participants.

Histopathological examination
The aortic valve was dissected and cut into parts. Part of the artery was fixed in a 4% poly-formaldehyde solution and processed for paraffin embedding and sectioning. The aortic valve was cut into 6-μm-thick transverse sections and stained with HE after fixation in a 4% polyformaldehyde solution embedded in paraffin wax. Then, the immunohistochemistry and immunofluorescence assays were performed to detect the expression level of USP14.
The samples of the aortic valve were taken to make paraffin sections. The paraffin sections were deparaffinized to water, sealed with H 2 O 2 , and washed with double distilled water. After antigen retrieval, immunohistochemistry was used to detect the USP14. The specific detection steps were performed according to the instructions of the VECTASTAIN® Elite® ABC Kit (Vector Laboratories Inc., 30 Ingold Road Burlingame, California, USA), as follows.
First, the antigen-fixed paraffin sections were washed with phosphate buffer solution (PBS) 2-3 times (5 min/ time), then blocked with 10% goat serum (TransGen Biotech Co., LTD, Beijing, China) at 37°C for 20 min. Filter paper was used to remove the serum. The mouse anti-human USP14 polyclonal antibody (dilution rate = 1:3000, SAB1406778, Sigma, St. Louis, USA) was added dropwise and incubated overnight at 4°C, washed three times with PBS (5 min/time), and incubated with goat anti-mouse monoclonal antibody (dilution rate = 1:5000, ab205719, Abcam, Cambridge, UK) at 37°C for 1 h. Diaminobenzidine (DAB) was used for color development  (PBS instead of the primary antibody was used as a negative control). Each paraffin section was photographed in 5 fields and counted [9].

RT-qPCR assay
Total RNA was extracted by the RNAiso Plus (Trizol) kit (Thermofisher, Massachusetts, USA), then reverse transcribed to cDNA. RT-qPCR was performed using a Light Cycler® 4800 System with specific primers for USP14. Table 3 presents the primer sequences. The glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as an endogenous control.

Western blotting analysis
Aortic valve tissues were frozen at -80°C. Total protein was isolated in a lysis buffer, resolved by 10% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE), and transferred onto polyvinylidene fluoride (PVDF) membranes by electroblotting. After blocking it with the Odyssey blocking buffer (LI-COR, Lincoln, USA) for 1 h, the membrane was incubated with primary antibodies at 4°C overnight, and then it was incubated with a horseradish peroxidase-conjugated secondary antibody at 37°C for 1 h. Next, the protein bands were scanned and analyzed with the Odyssey scanner [10]. GAPDH was used as an endogenous control. Mouse anti-human GAPDH polyclonal antibody (dilution rate = 1:5000, WH0002597M1, Sigma, St. Louis, USA) was used to detect GAPDH, and the secondary antibody was also goat anti-mouse monoclonal antibody (dilution rate = 1:5000, ab205719, Abcam, Cambridge, UK). USP14 protein was detected using a mouse anti-human USP14 polyclonal antibody (dilution rate = 1:3000, SAB1406778, Sigma, St. Louis, USA). The secondary antibody was goat anti-mouse monoclonal antibody (dilution rate = 1:5000, ab205719,

Identification of DEGs
The merged series contained 5792 up-regulated genes and 3844 down-regulated genes (Fig.1a).

WGCNA analysis and attainment of module DEGs
In this study, WGCNA analysis was conducted using the R package WGCNA. The cluster of patients is shown in Fig. 1b. As shown in Fig. 1c, a power value of 8 was the lowest power for which scale independence was below 0.9, so it was selected to produce a hierarchical clustering tree of the 6301 genes (Fig.  1d). In addition, a dendrogram and heatmap were used to quantify module similarity by correlation (Fig.  1e). Interactions between these modules were then  (Fig. 1f). The associations between clinical traits and the modules were identified based on the correlation between module and clinical traits (Fig.  1g). The heat map of the MEdarkgrey model DEGs' expression is shown in Fig. 1h.

Functional and pathway enrichment analysis of DEGs
GSEA to DEGs elicited an enrichment of the "heart process" and "cardiac muscle contraction" (Fig. 2a-b, Tables 4 and 5). Metscape analysis of DEGs showed an enrichment in the autophagy-animal, cellular lipid catabolic process, apoptosis, and glycoside metabolic process categories (Fig. 2c-e).

Construction and analysis of the protein-protein interaction network
The PPI network of the DEGs was constructed via the STRING online database and analyzed by Cytoscape software (Fig. 3a). Four different algorithms were employed to identify hub genes and 11 common hub genes were obtained (Fig. 3b). A summary of common hub genes is shown in Table 6. The PPI network of common hub genes is shown in Fig. 3c. The heat map of common hub genes is shown in Fig. 3d.

Identification of hub genes
The CTD database showed that common hub genes targeted heart diseases. It also highlighted pathophysiological processes such as fibrosis, inflammation, and necrosis, which are relevant in aortic stenosis pathogenesis (Fig. 4).  Prediction and functional annotation of miRNA associated with hub genes The miRNAs that regulate the hub genes were screened out with TargetScan (Table 7). GO and KEGG analyses of the miRNA were performed by DIANA-miRPath. GO and KEGG enrichment analysis to the screened miRNA demonstrated a preponderance of the cellular lipid metabolic process, apoptotic signaling pathway, nucleic acid binding transcription factor activity, toll-like receptor signaling pathway, Hippo signaling pathway, Arrhythmogenic right ventricular cardiomyopathy, and others (Fig. 5).

Logistic regression analysis and identification of USP14
The ROC curves of common hub genes are shown in Fig. 6. Likewise, the common hub genes and their effect on aortic stenosis, based on univariate logistic proportional regression analysis, are shown in Table 8 and those based on multivariate logistic proportional regression analysis in Table 9. Of the eleven potential biomarkers, USP14 was chosen, partly due to its higher odds radio (OR) calculated after multivariate logistic regression analysis (Table 9).

Histological patterns of the aortic valve
The HE staining demonstrated that the fibrous tissue and muscle cells were arranged regularly in the control aortic valve. However, the fibrous tissue and muscle cells were arranged irregularly in the aortic stenosis tissue (Fig. 7).

The expression level of USP14 in the histology
The immunohistochemistry assay manifested that USP14 was upregulated in the aortic stenosis tissues, compared with the control aortic valve (Fig. 8). Also, an immunofluorescence assay corroborated such finding (Fig. 9).

Verification of the expression of USP14 in the mRNA and protein levels
The RT-qPCR assay showed that the expression of USP14 was higher in the aortic stenosis than the control group (P < 0.05, Fig. 10a). Moreover, in terms of the protein levels, the Western blotting confirmed that USP14

Discussion
In this study, a WGCNA on microarray data from aortic valves predicted that USP14 would be highly expressed in the valves of aortic stenosis patients, which was effectively verified in molecular experiments with samples collected from patients. USP14 may become a target for the diagnosis and treatment of aortic stenosis. Aortic stenosis, a progressive and chronic disease, is one of the most serious cardiac valvular diseases [11]. In   the sclerotic phase, aortic valves exhibit punctate calcification and valve thickening without hemodynamic changes and no obvious clinical symptoms, that is, aortic valve hardening. In the later stage, the valve ring and valve tissue appear to thicken, resulting in aortic valve stenosis and left ventricular outflow tract obstruction, which will eventually lead to serious obstruction of the left ventricular ejection function. The end stage of aortic stenosis is characterized by calcification, sided by the rupture of basement membranes, inflammatory cell infiltration, lipid deposition, endothelial dysfunction, and a series of activation processes caused by the transformation of valvular interstitial cells into osteoblasts [12]. USP14 is the only proteasome-binding deubiquitinating enzyme in the USP family. It can be recruited to 19 s particles to prevent it from being recognized and degraded by proteasome recognition, by cutting the ubiquitin chain on the protein substrate [13]. Accordingly, the inhibition of USP14 could improve the proteasome degradation activity [14]. USP14 is amplified and overexpressed in many types of diseases [15][16][17][18]. Overexpression of USP14 promotes cell proliferation and migration, while down-regulation induces cell apoptosis and inhibits cell proliferation, migration, and invasion. The expression of USP14 plays a pathogenic role in various types of diseases and is negatively correlated with the prognosis of the disease [19].
The coordination of the two reversible processes of ubiquitination and de-ubiquitination can finely regulate the type and degree of ubiquitination of proteins, affect the function and fate of corresponding proteins, and play an important role in many life activities such as cell division, differentiation, and signal transduction [20]. The USP14 protein has a total length of 494 amino acids and is divided into two domains: a ubiquitin-like domain (UBL) and a catalytic domain (CAT) [21]. It is one of the three deubiquitinating enzymes that bind to proteasome and can interact dynamically with proteins through ubiquitin. They regulate the function, metabolism, and degradation of proteins in the body, thus affecting a series of life activities such as the cell cycle, cell differentiation, transcription, DNA repair, immunity, and inflammation. Dishevelled (Dvl) is a key regulator of Wnt signal transduction. USP14 can enhance the signal transduction of Wnt/β-catenin. The inhibition of USP14 can also increase the degree of ubiquitination of Dvl and significantly inhibit the downstream of Wnt signal transduction [18]. Therefore, the consequence of Wnt/β-catenin activation, due to USP14 Dvl de-ubiquitination, is that cell proliferation would be activated greatly, which would cause aortic valve hypertrophy and then lead to aortic stenosis. USP14 promotes the progression of cardiac hypertrophy by increasing GSK-3β phosphorylation Fig. 9 The detection of USP14 content of the aortic valve by immunofluorescence assay [22]. This study found that USP14 is highly expressed in aortic valve tissue, which inhibits the proteasome degradation activity so this research hypothesizes that the inflammatory factors and lipoproteins gathered in the aortic valve cannot be effectively degraded, resulting in aortic stenosis.
Many studies have shown that the development of aortic stenosis is similar to atherosclerosis, which is a lipid inflammatory process [23]. However, the nature of the specific endogenous stimuli triggering aortic stenosis is not yet clear. Excessive activation of macrophages plays a very important role in inflammatory diseases. Under the stimulation of inflammatory factors or bacteria, the activation of macrophages produces more inflammatory cytokines, causing a cascade of inflammation, which eventually becomes uncontrolled inflammation.
Nuclear factor-κB (NF-κB) is an inducible transcriptional regulator of many inflammatory cytokine genes [24]. One study reported that overexpression of USP14 in lung epithelial cells can reduce IκB protein levels and increase cytokine release [25]. In addition, another study showed that USP14 affected NF-κB activation induced by lipopolysaccharide [26]. Based on these results, USP14 can induce the activation of NF-κB and regulate the up-regulation of tumor necrosis factor-α (TNF-α), Fig. 10 Verification of the expression of USP14 in the mRNA and protein level. a Relative expression of USP14 by RT-qPCR analysis. *p < 0.05, compared with control. b, c Western blotting expression of USP14 in the control and aortic stenosis groups interleukin-6 (IL-6), and interleukin-18 (IL-18) expression by activating the NF-κB signal pathway. These inflammatory factors interact to form a positive feedback loop to increase the production of inflammatory factors, aggravating the inflammatory response in the valve, promoting lipid deposition, and ultimately promoting the occurrence of aortic stenosis.

Strength and study limitation
About the strength of the research, bioinformatics and clinical data were combined for analysis in this study. Many algorithms and experimental methods have been used to identify and verify the USP14 as an underlying gene biomarker for aortic stenosis. Although this study has carried out rigorous bioinformatic analysis, there are still some limitations. No animal experiments have been carried out. It is still necessary to further use an animal model to obtain more accurate results, in order to better understand the molecular mechanisms played by USP14 in aortic stenosis.

Conclusions
In summary, bioinformatics technology could be one useful tool to predict the progression of aortic stenosis, by exploring the mechanism of its occurrence and development. DEGs, especially USP14, might be involved in the occurrence and development of aortic stenosis and might eventually become a biomarker and a target to treat this disease. Detection of USP14 in the clinical practice might provide the better evidence to guide the early diagnosis and treatment of aortic stenosis. In the future, diagnostic reagent kit of USP14 and the corresponding drugs should be developed to perform the multicenter randomized controlled clinical trial.