- Open Access
Screening of differentially expressed microRNAs of essential hypertension in Uyghur population
Lipids in Health and Disease volume 18, Article number: 98 (2019)
Essential hypertension can cause many kinds of cardiovascular diseases. The pathogenesis of essential hypertension is very complex, and the mechanism is still unclear. The microRNAs have been identified as novel biomarkers for pre-diagnosis and prognosis of hypertension. However, the kinds of microRNAs that can be used as specific biomarkers for hypertension are unknown.
Methods and results
Plasma samples were isolated from Uyghur subjects with essential hypertension and the healthy individuals. Microarray was used to identify differentially expressed microRNAs. The microarray data were clustered and annotated with online software. The target genes of differentially expressed microRNAs were also analyzed. The microarray results were further verified by quantitative real-time PCR. We identified 257 microRNAs that were differentially expressed between patients with essential hypertension and the healthy individuals. These microRNAs had a total of 6580 target genes. The 47 microRNAs that had target genes, including 24 up-regulated and 23 down-regulated microRNAs, were further screened out to construct a reference set of potential microRNA biomarkers. Most of the 47 microRNAs were located at chromosome 19 (40 microRNAs) and chromosome 1 (45 microRNAs). Their target genes were mainly enriched in metal ion binding, transcription regulation, cell adhesion and junction, indicating that these candidate microRNAs may regulate mineral ion binding and cell communication process of essential hypertension. The quantitative real-time PCR results of miR-198 and miR-1183 (which were the two most significantly up-regulated microRNAs by microarray), and, miR-30e-5p and miR-144-3p (which were the two most significantly down-regulated microRNAs by microarray) were consistent with the microarray results.
A reference set of potential microRNA biomarkers that may be involved in essential hypertension is constructed. Our study may provide experimental evidence for further studying the mechanism of essential hypertension.
Essential hypertension is not only one of the main causes but also one of the clinical manifestations of cardiovascular diseases. It is a highly complex disease that can be induced by both environmental and genetic factors [1, 2]. The mechanism of essential hypertension is still unclear. There are still no precise pre-diagnosis methods or efficient treatment options for essential hypertension [2,3,4]. Thus, it is necessary to identify specific markers for essential hypertension.
MicroRNA is a kind of small non-coding RNA with about 20 bp in length that can regulate gene expression at the post-transcriptional level [4,5,6]. Previous studies demonstrate that microRNAs can be used as diagnostic biomarkers for hypertension and cardiovascular diseases [2, 4]. The abnormal expression of microRNAs in tissues and body fluids has also been reported [2, 4, 6]. It is also found that there were 27 microRNAs differentially expressed in patients with essential hypertension and that the relationship between hmcv-miR-UL122 and hypertension was identified . In addition, the increased expression of miR-505 may play important roles in essential hypertension . Tissue-based studies also illustrate that several microRNAs play remarkable roles in hypertension, such as miR-181a [1, 9], miR-1, miR-21, miR-9 and miR-126 [10, 11].
MicroRNA array and quantitative PCR are two main techniques for detection and examination of microRNAs [2,3,4, 6]. The microRNA profile of hypertension has been investigated. However, different microRNA biomarker sets were obtained due to the diversity of experiment methods and research subjects [1,2,3,4]. More microRNA markers for essential hypertension remain to be identified. In this study, we screened the microRNA profile of essential hypertension in Uyghur population. Our study may provide experimental evidence for investigating the possible role and mechanism of microRNAs in essential hypertension.
Materials and methods
Sample collection and processing
Eight Uyghur subjects aged between 30 and 40 years old were recruited. Four subjects, including 2 males and 2 females, with untreated essential hypertension were included in the hypertension group. Inclusion criteria: 1) blood pressure should meet one of the following criteria: Sitting blood pressure was higher than or equal to 140/90 mill Hg (1 mmHg = 0.133 kPa) measured by Mercury sphygmomanometer; and/or, average daily blood pressure was higher than or equal to 135/85 mill Hg by 24-h ambulatory blood pressure monitor; 2) patients diagnosed with essential hypertension for the first time; 3) patients were untreated. Exclusion criteria: 1) patients with secondary hypertension; 2) patients with cardiac function III, IV, coronary heart disease, severe pulmonary insufficiency, rheumatic diseases, diabetes, endocrine disorders, severe liver or kidney dysfunction, stroke, peripheral arterial disease, or, mental illness; 3) pregnant or lactating women. Four age and gender matched healthy individuals with normal blood pressure (2 males and 2 females) were included in the control group. Healthy subjects did not smoke and had no history of cancer, diabetes, kidney failure, stroke or peripheral arterial disease. For verification of microRNA expression, another 15 patients with essential hypertension (8 males and 7 females) and 15 healthy individuals (4 males and 11 females) were included. The inclusion and exclusion criteria were same as above.
All subjects were fasted before sampling and sampled at the same time in the morning. They were asked to sit still for 10 min before blood collection. Totally, 3 mL of peripheral blood was collected from each subject and the plasma was then isolated.
All subjects or their families signed the informed consent. The study was approved by Ethical committee of Xinjiang Medical University.
RNA extraction and miRNA microarray analysis
Total RNAs were extracted from plasma using mirVana PARIS kit (Ambion-1556) and quantified by NanoDrop ND-2000 (Thermo Scientific). RNA quality was checked by Bioanalyzer 2100 (Agilent, G2939A). MicroRNA expression was detected by Agilent Human miRNA Microarray (Agilent, 8*60 K, ID070156). RNAs were hybridized with the array in hybridization oven (Agilent, G2545A). Then, array signals were scanned by Agilent Scanner (Agilent, G2505C). All processes were performed in accordance with the standard guidelines provided by the manufacturer.
The signals were then analyzed by Feature Extraction software (version10.7.1.1, Agilent Technologies). Raw data was obtained, and normalized with the quartile algorithm using Genespring software (version 12.5; Agilent Technologies) to calculate the microRNA expression. MicroRNAs with fold changes greater than two and a p-value below 0.05 were defined as differentially expressed.
MicroRNAs were clustered using the Hierarchical Clustering method. The chromosome loci of microRNAs were annotated by miRBase (http://www.mirbase.org/) and HGNC (Hugo Gene Nomenclature Committee) (http://www.genenames.org/). Cytogenetic locations of microRNAs and OMIM (Online Mendelian Inheritance in Man) record of essential hypertension were retrieved from OMIM database (http://www.omim.org/).
Target gene prediction, GO annotation and pathway analysis
Target genes of microRNAs were predicted using three databases of Targetscan, (http://www.targetscan.org/), microRNAorg, (http://www.microrna.org/), and pita, (https://genie.weizmann.ac.il/pubs/mir07/), respectively. The resulted genes shared by three databases were selected and then analyzed using Gene Ontology (GO) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway. Furthermore, genes were clustered according to the annotation information by DAVID analysis tool (https://david.ncifcrf.gov/).
Quantitative real-time PCR
Quantitative real-time PCR (qPCR) was conducted to verify microRNA expression. MicroRNAs were extracted using QIAzol Lysis Reagent (QIAGEN, 5346994), miRNeasy Serum/Plasma Spike-In Control (QIAGEN, 219610) and miRNeasy Serum/Plasma Kit (QIAGEN,217184). cDNA was prepared using miScript II RT Kit (QIAGEN, 218161). MicroRNA expression was detected using miScript SYBR Green PCR Kit (QIAGEN, 218073) on ABI7500 (ABI). Four pairs of Primers were used in qPCR assays, including hsa-miR-1183 (CD201–0149, TIANGEN BIOTECH), hsa-miR-198 (CD201–0085, TIANGEN BIOTECH), hsa-miR-30e-5p (CD201–0340, TIANGEN BIOTECH) and hsa-miR-144-3p (CD201–0011, TIANGEN BIOTECH). miRNA39 was used for normalization in the qRT-PCR assays. The relative gene expression was calculated with 2-ΔΔCt method.
Prism statistical software (GraphPad Software Inc) was used for data analysis. Difference in microRNA expression between essential hypertension group and control group was analyzed by student t-test. A value of p ≤ 0.05 was considered to be statistically significant.
The expression profile of microRNAs in the essential hypertension group
Totally, 257 microRNAs that showed differential expression between patients with essential hypertension and the control group were obtained (Table 1). The number of up-regulated microRNAs (161) was higher than that of down-regulated microRNAs (96) (Table 1). Furthermore, fold change distribution of the differentially expressed microRNAs showed distinct pattern. Up-regulated microRNAs changed much more than the down-regulated microRNAs (Fig. 1). The change of hsa-miR-1183 expression (up-regulated, fold change = 44.41) and hsa-miR-30e-5p expression (down-regulated, fold change = 24) were the most significant (Table 2). Cluster analysis classified samples into separate groups based on the microRNA expression in each sample. The microRNA expression of patients in the essential hypertension group was obviously different from that of healthy controls (Fig. 2). Overall, nearly 2/3 microRNAs showed higher level in the essential hypertension samples (Fig. 2).
Target gene identification and filteration of candidate microRNA markers
Target genes of the microRNAs were predicted by microRNAorg, PITA and TargetScan. Target genes that were simutaneously predicted by these three tools were selected (Fig. 3 and Fig. 4). In total, 6580 target genes were obtained (Additional file 1: Figure S1). We further screened out 47 microRNAs that had target genes (Tables 1 and 2). The number of up-regulated and down-regulated microRNAs that had target genes was almost equal (Table 1). Hsa-miR-765, hsa-miR-940, hsa-miR-1183, hsa-miR-557 and hsa-miR-939-5p had over 300 target genes. Hsa-miR-765 had the largest number of target genes (475) (Table 3).
Functional classification and pathway analysis of the target genes
The biological functions of the target genes were analyzed by GO. Most target genes were involved in the regulation of transcription, signal transduction, organismal development and cell adhesion (Additional file 2: Figure S2A). The number of target genes in cytoplasm and nucleus was almost equal (Additional file 2: Figure S2B). Over half of the target genes (53.8%) had protein binding function, and about 40% of the genes may play roles in zinc-ion binding (17.5%) or metal-ion binding (23.9%) (Additional file 2: Figure S2C).
The target genes were classified into different functional groups based on GO-terms. It is well worth noting that most of genes were enriched in ion binding group (ES = 3.46) and transcription regulation group (ES = 3.44). Target genes involved in cell fraction and cell adhesion were enriched with an enrichment score of 5.55 and 4.33, respectively. In addition, pathway analysis was performed to analyze the signal pathways that target genes or microRNAs may participate in. We found that the target genes were mainly involved in cancer pathways, MAPK pathway, regulation of actin cytoskeleton and focal adhesion (Fig. 5).
Cytogenetic location and OMIM annotation of the microRNAs
Cytogenetic location of the 47 microRNA candidates were analyzed. Most of them were located at chromosome 19 (40 microRNAs) and chromosome 1 (45 microRNAs) (Table 3). No microRNA was found on chromosome 12, 15, 18 or Y. Moreover, OMIM records of the patients with essential hypertension were evaluated. Totally, 13 OMIM records about essential hypertension were obtained. The records illustrated the relationship between microRNAs and disease status (Table 3). Hsa-miR-1182, hsa-miR-557 and hsa-miR-144-3p showed a direct correlation with essential hypertension (OMIM ID #145500). Moreover, miR-92b-3p, miR-939-5p, and miR-769-3p were annotated with lipodystrophy, glucocorticoid and apolipoprotein E, which are hypertension-related indicators.
Verification of the potential microRNA markers and microRNA-gene interaction
Four microRNA candidates with significant expression change were verified by qPCR. MiR-198 and miR-1183 showed a higher expression in patients with essential hypertension while miR-30e-5p and miR-144-3p showed lower expression (Fig. 6). Fold change of microRNAs (miR-198, miR-1183, miR-30e-5p and miR-144-3p) measured by qPCR was in accordance with the results of microRNA array (Table 4). In order to explore the relationship between microRNAs and target genes, a microRNA-gene interaction network was constructed (Fig. 7). Target genes of the four microRNAs were selected with a tougher criteria. Target genes were graded using miRBase database firstly and then the genes with scores of more than 80 were selected. Altogether, 33 target genes were obtained. MiR-144-3p shared 21 target genes with miR-30e-5p, while it didn’t share any target gene with miR-198 (Fig. 7). MiR-1183 and miR-198 both regulated ARFGAP2. MiR-1183 and miR-144-3p both targeted MMGT1, ST6GALNAC3 and ZBTB21.
In this study, we identified 257 differentially expressed microRNAs from 4 Uyghur patients (from Xinjiang, China) with essential hypertension and constructed a candidate microRNA biomarker set of 47 microRNAs. To the best of our knowledge, this study is the first report to screen microRNA markers of essential hypertension in Uyghur population of Xinjiang, China. Our results add new information to the potential microRNA markers of essential hypertension.
Among the 47 microRNAs, the number of up-regulated and down-regulated microRNAs was almost equal. However, we found that the expression of up-regulated microRNAs changed more significantly. Compared with the microRNA profiles of hypertension reported by previous studies [2,3,4], many new microRNAs were identified in this study. The difference may be due to the complexity of hypertension or the differences in sampling population or methods [2, 4].
Previous studies have shown that essential hypertension may be caused by genetic factors [1, 2]. Here, we found that the microRNAs were mostly located on chromosome 1 and chromosome 19. Five miRNAs located at different loci of chromosome 1 and seven miRNAs located at chromosome 19, five of which located at the same locus (19q13.42). Further analysis through the OMIM database revealed that miR-1182 was annotated with angiotensinogen, a protein related to essential hypertension [4, 8]. Moreover, miR-557 was related to ATPase and Na+/K+ transport, which play important roles in the blood pressure regulation [12, 13]. miR-769-3p was annotated with apolipoprotein, a molecule that is closely related to cardiac diseases . These results demonstrate that the idenfied microRNAs are related with the metabolic processes of essential hypertension, indicating that these microRNAs may be involved in the development of essential hypertension.
Target genes were clustered into 8 groups by functional annotation and most of them were enriched in metal ion binding and transcription regulation (Table 5). Abnormal Ca2+/Na+ is one of the common phenomena in hypertension and cardiovascular diseases . Annotation analysis both in microRNAs and target genes showed a tendency in ATP and mineral ions metabolism. This suggests that these microRNAs may be used as potential diagnostic markers. Further studies on the role and mechanism of microRNAs in essential hypertension are needed.
In our study, target genes related to cell adhesion and junction were enriched into two clusters. Researchers have found that miR-92a may play a role in mediating cell-to-cell communication [15,16,17]. The miR-1 may affect the genes of gap-junction channels [12, 18]. Our results reveal that microRNAs may participate in regulating cell communications in essential hypertension. Pathway analysis also supported their role in cell to cell communication. Genes in this pathway are involved in focal adhesion and regulation of actin cytoskeleton and adherence junction. However, sources of microRNAs in plasma are various. Apoptotic, necrotic cell death and active secretion all result in detectable microRNAs in plasma .
MiR-198 and miR-1183 were the two most significantly up-regulated microRNAs according to the microarray results, while miR-30e-5p and miR-144-3p were the two most significantly down-regulated microRNAs. QPCR further confirmed the microarray results. However, the fold change of microRNA expression measured by qPCR was not as significant as that by microarray. Li et al. found that miR-1183 was a potential diagnostic marker for rheumatic heart disease . High expression level of miR-198 is reported in acute coronary syndromes . miR-30e-5p is found to suppress the expression of hypertension related gene ADRA2A in human endothelial cells . miR-144-3p regulates cholesterol homeostasis and is identified as a potential therapeutic target of atherosclerosis and acute myocardial infarction . However, how are these four microRNAs involved in essential hypertension are unclear. The significant changes of their expression in this research may pave a way for further study.
Interaction analysis can greatly help us to understand the relationship of microRNAs with target genes. One microRNA may regulate multiple genes, and one gene might be regulated by different microRNAs [3,4,5]. In this study, the target genes of miR-198, miR-1183, miR-30e-5p and miR-144-3p were predicted and a microRNA-gene interaction network was constructed. Totally, 33 target genes were obtained. However, functional analysis and validation of the target genes are not performed and needs further study.
This study has some limitations. First, this study focused on Uyghur subjects from Xinjiang, China. Thus, the results of this study may lack representativeness for the general population. Second, the sample size was relatively small. Third, the number of differentially expressed microRNAs validated by qPCR was small. Fourth, the function of the 33 target genes of miR-1183, miR198, miR144-3p and miR-30e-5p was not further analyzed. Further studies are warranted.
In conclusion, we constructed a reference microRNA pool of essential hypertension and obtained a set of circulating microRNAs, which may be used as potential biomarkers for essential hypertension. These microRNAs showed a remarkable correlation with mineral ion metabolism and cellular communication. Four microRNAs (miR-1183, miR-198, miR-30e-5p and miR-144-3p) that showed most prominent change were further verified. A regulatory network between these microRNAs and their target genes was also constructed, which may provide experimental evidence for further study of the mechanism of essential hypertension.
Hugo Gene Nomenclature Committee
Kyoto Encyclopedia of Genes and Genomes
Online Mendelian Inheritance in Man
Quantitative Real-time PCR
Marques FZ, Campain AE, Tomaszewski M, Zukowska-Szczechowska E, Yang YH, Charchar FJ, et al. Gene expression profiling reveals renin mRNA overexpression in human hypertensive kidneys and a role for microRNAs. Hypertension. 2011;58:1093–8.
Romaine SP, Charchar FJ, Samani NJ, Tomaszewski M. Circulating microRNAs and hypertension--from new insights into blood pressure regulation to biomarkers of cardiovascular risk. Curr Opin Pharmacol. 2016;27:1–7.
Bronze-da-Rocha E. MicroRNAs expression profiles in cardiovascular diseases. Biomed Res Int. 2014;2014:985408.
Marques FZ, Booth SA, Charchar FJ. The emerging role of non-coding RNA in essential hypertension and blood pressure regulation. J Hum Hypertens. 2015;29:459–67.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.
Wahid F, Shehzad A, Khan T, Kim YY. MicroRNAs: synthesis, mechanism, function, and recent clinical trials. Biochim Biophys Acta. 2010;1803:1231–43.
Li S, Zhu J, Zhang W, Chen Y, Zhang K, Popescu LM, et al. Signature microRNA expression profile of essential hypertension and its novel link to human cytomegalovirus infection. Circulation. 2011;124:175–84.
Yang Q, Jia C, Wang P, Xiong M, Cui J, Li L, et al. MicroRNA-505 identified from patients with essential hypertension impairs endothelial cell migration and tube formation. Int J Cardiol. 2014;177:925–34.
Marques FZ, Romaine SP, Denniff M, Eales J, Dormer J, Garrelds IM, et al. Signatures of miR-181a on renal transcriptome and blood pressure. Mol Med. 2015;21(1):739–48.
Kontaraki JE, Marketou ME, Zacharis EA, Parthenakis FI, Vardas PE. MicroRNA-9 and microRNA-126 expression levels in patients with essential hypertension: potential markers of target-organ damage. J Am Soc Hypertens. 2014;8:368–75.
Kontaraki JE, Marketou ME, Zacharis EA, Parthenakis FI, Vardas PE. Differential expression of vascular smooth muscle-modulating microRNAs in human peripheral blood mononuclear cells: novel targets in essential hypertension. J Hum Hypertens. 2014;28:510–6.
Luo X, Yang B, Nattel S. MicroRNAs and atrial fibrillation: mechanisms and translational potential. Nat Rev Cardiol. 2015;12:80–90.
O'Connor DT, Zhu G, Rao F, Taupenot L, Fung MM, Das M, et al. Heritability and genome-wide linkage in US and australian twins identify novel genomic regions controlling chromogranin a: implications for secretion and blood pressure. Circulation. 2008;118:247–57.
Sinha R, Singh R. Role of apolipoprotein AI gene polymorphism (G-75A and C+83T) in essential hypertension in Indian population. Ann Clin Lab Sci. 2014;44:298–303.
Emanueli C, Shearn AI, Angelini GD, Sahoo S. Exosomes and exosomal miRNAs in cardiovascular protection and repair. Vasc Pharmacol. 2015;71:24–30.
Mittelbrunn M, Sanchez-Madrid F. Intercellular communication: diverse structures for exchange of genetic information. Nat Rev Mol Cell Biol. 2012;13:328–35.
Valadi H, Ekstrom K, Bossios A, Sjostrand M, Lee JJ, Lotvall JO. Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat Cell Biol. 2007;9:654–9.
Yang B, Lin H, Xiao J, Lu Y, Luo X, Li B, et al. The muscle-specific microRNA miR-1 regulates cardiac arrhythmogenic potential by targeting GJA1 and KCNJ2. Nat Med. 2007;13:486–91.
van Empel VP, De Windt LJ, da Costa Martins PA. Circulating miRNAs: reflecting or affecting cardiovascular disease? Curr Hypertens Rep. 2012;14:498–509.
Li N, Lian J, Zhao S, Zheng D, Yang X, Huang X, et al. Detection of differentially expressed MicroRNAs in rheumatic heart disease: miR-1183 and miR-1299 as potential diagnostic biomarkers. Biomed Res Int. 2015;2015:524519.
Li M, Zhang J. Circulating MicroRNAs: potential and emerging biomarkers for diagnosis of cardiovascular and cerebrovascular diseases. Biomed Res Int. 2015;2015:730535.
Kriegel AJ, Baker MA, Liu Y, Liu P, Cowley AW Jr, Liang M. Endogenous microRNAs in human microvascular endothelial cells regulate mRNAs encoded by hypertension-related genes. Hypertension. 2015;66:793–9.
Hu YW, Hu YR, Zhao JY, Li SF, Ma X, Wu SG, et al. An agomir of miR-144-3p accelerates plaque formation through impairing reverse cholesterol transport and promoting pro-inflammatory cytokine production. PLoS One. 2014;9:e94997.
This work was supported by the National natural Science Foundation of China [grant number NSFC81660071]. The funding source had no involvement in study design, in the collection, analysis and interpretation of data, in the writing of the report, and in the decision to submit the article for publication.
Availability of data and materials
The datasets during and analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The study was approved by ethical committee of Xinjiang Medical University.
Consent for publication
All patients and their families signed the informed consent.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S2. GO analysis. (A) Target gene coverage of GO biological process term. (B) Target gene coverage of GO cellular component terms. (C) Target gene coverage of GO molecular function terms (TIF 503 kb)