Don M. Matshazi1*
Cecil J. Weale1Rajiv T. Erasmus2Andre P. Kengne3,4
Saarah F. G. Davids1Shanel Raghubeer1Glenda M. Davison1
Tandi E. Matsha1*- 1SAMRC/CPUT/Cardiometabolic Health Research Unit, Department of Biomedical Sciences, Faculty of Health and Wellness Science, Cape Peninsula University of Technology, Cape Town, South Africa
- 2Division of Chemical Pathology, Faculty of Health Sciences, National Health Laboratory Service (NHLS) and Stellenbosch University, Cape Town, South Africa
- 3Non-Communicable Diseases Research Unit, South African Medical Research Council, Cape Town, South Africa
- 4Department of Medicine, University of Cape Town, Cape Town, South Africa
MicroRNAs are non-coding, post-transcriptional regulators of gene expression and their dysregulation has been associated with development of various diseases, including hypertension. Consequently, understanding their role in the pathogenesis and progression of disease is essential. Prior research focusing on microRNAs in disease has provided a basis for understanding disease prognosis and offered possible channels for therapeutic interventions. Herein, we aimed to investigate possible differences in the expression profiles of five microRNAs in the blood of participants grouped on the basis of their hypertension status. This was done to elucidate the possible roles played by these microRNAs in the development of hypertension. Using quantitative reverse transcription polymerase chain reaction, we evaluated the expression levels of miR-126-3p, 30a-5p, 182-5p, 30e-3p, and 1299 in the whole blood of 1456 participants, normotensive (n = 573), screen-detected hypertensive (n = 304) and known hypertensive (n = 579). The expression of miR-126-3p and 182-5p was significantly higher in known hypertensives relative to both screen-detected hypertensives and normotensives, and also in screen-detected hypertensives vs normotensives. A significant association between the expression of miR-126-3p, 182-5p, and 30a-5p and known hypertension was also evident. This study demonstrated dysregulated miR-126-3p, 182-5p, and 30a-5p expression in hypertension, highlighting the possible efficacy of these microRNAs as targets for the diagnosis of hypertension as well as the development of microRNA-based therapies.
Introduction
Hypertension (HPT) is a complex and multifactorial disease responsible for considerable loss of life globally (Campbell et al., 2016a; Mills et al., 2020). It is an important, modifiable risk factor for cardiovascular disease [1, 3], whose prevalence varies globally, but is on an upward trajectory in sub-Saharan Africa (Campbell et al., 2016a; Mills et al., 2016). The prevalence rates of HPT in some sub-Saharan African countries currently rank among the highest on the globe, in stark contrast to the picture from a few decades ago when the region had the lowest blood pressure levels (Danaei et al., 2011). Inroads into understanding the pathophysiology of HPT and advancing treatment options have been made over the years. The pathogenesis has been linked to various biological processes, including endothelial dysfunction, impaired angiogenesis, dysregulation of the renin-angiotensin-aldosterone axis and platelet activation (Taddei et al., 2001; Gkaliagkousi et al., 2010; Touyz et al., 2018). Sizeable financial investments have been made to study the genetic and environmental determinants of HPT and some of these studies have linked the Liddle phenotype (a hereditary disorder characterized by overactivity of the renal tubular epithelial sodium channel as a result of mutations in the SCNN1B and SCNN1G genes), and salt and water retention to the high HPT prevalence levels currently observed in sub-Saharan Africa (Spence and Rayner, 2018). However, in up to 95% of HPT cases, the etiology remains unknown and primary HPT continues to be a leading cause of morbidity and premature mortality globally (Carretero and Oparil, 2000; He and Macgregor, 2007). Described as a “silent killer,” its symptoms manifest in the later stages of the condition when hypertension-mediated target organ damage has possibly taken place (Moore, 2005).
Processes involved in blood pressure homeostasis are tightly regulated by various systems in the body. Amongst others, cellular processes like differentiation, growth and metabolism, are known to be under the control of microRNAs (miRNAs) (Vidigal and Ventura, 2015). These miRNAs are 18-25 nucleotide long, non-coding, post-transcriptional regulators of gene expression and their dysregulation has been linked to the development of cancer, essential HPT, viral disease and endothelial dysfunction (Chen et al., 2008; Li and Kowdley, 2012; Nemecz et al., 2016; Silambarasan et al., 2016). It is plausible that dysregulation of miRNAs may lead to disturbances in the body’s blood pressure regulatory mechanisms and play an important role in the development of HPT. Herein, we investigated circulating levels of miR-1299, miR-30a-5p, miR-30e-3p, miR-126-3p, and miR-182-5p in participants with normal blood pressure, as well as screen-detected and known hypertensives on anti-hypertensive treatment. These miRNAs were chosen as possible targets in our study as next generation sequencing data had previously shown their association with diabetes and HPT (Matsha et al., 2018; Matshazi et al., 2021).
Materials and Methods
Design of the Study and Description of Procedures
This study included males and females from the Vascular and Metabolic Health (VMH) study, a sub-study of the Cape Town Bellville South Study (Matsha et al., 2012). The strategies followed with regards to data collection and the conduction of various procedures in this study were previously described by Matsha et al. (2018). In summary, various anthropometric measurements were taken from each participant and for each variable, the average of three separate readings taken was reported. Body Mass Index was calculated as weight per square meter (kg/m2) where kg is a participant’s weight in kilograms and m2 is the square of their height in meters. World Health Organization guidelines (Chalmers et al., 1999) informed the process of blood pressure measurement. In brief, blood pressure was measured from the right arm of a participant who was in a sitting position and had rested for at least 10-min, using a semi-automatic digital blood pressure monitor (Omron M6 comfort-preformed cuff blood pressure monitor, China). The blood pressure was measured thrice, at three-min intervals, and the lowest systolic blood pressure and corresponding diastolic blood pressure values reported. Participants were then put into their respective blood pressure categories based on; the use of anti-hypertensive medication as known HPT, blood pressure measurement of 140/90 mm Hg and above as screen-detected HPT and normal blood pressure measurement (less than 140/90 mm Hg) as normotensive.
A number of molecular methods were utilized to measure several biochemical parameters in blood samples collected from study participants. These measurements were conducted in an ISO 15189 accredited pathology practice (PathCare Reference Laboratory, Cape Town, South Africa). The molecular methods used included High Performance Liquid Chromatography (BioRad Variant Turbo, BioRad, Hercules, CA, United States), for measurement of glycated hemoglobin, whilst a paramagnetic particle chemiluminescence assay (Beckman DXI, Beckman Coulter, South Africa) was used to measure serum insulin. The Competitive Chemiluminescent (Immulite 2000, Siemens, Munich, Germany) was used to determine levels of cotinine in serum whilst plasma glucose concentration was measured by the enzymatic hexokinase method (Beckman AU, Beckman Coulter, Brea, CA, United States). An enzymatic immunoinhibition—end point assay (Beckman AU, Beckman Coulter, Brea, CA, United States) was used to determine both total cholesterol and high-density lipoprotein cholesterol, whilst triglycerides were measured using a glycerol phosphate oxidase-peroxidase, end point assay (Beckman AU, Beckman Coulter, Brea, CA, United States). The enzymatic selective protection—end point (Beckman AU, Beckman Coulter, Brea, CA, United States) assay was used to determine the levels of low-density lipoprotein cholesterol. In order to determine ultrasensitive C-reactive protein (CRP) levels, a Latex Particle Immunoturbidimetry (Beckman AU, Beckman Coulter, Brea, CA, United States) assay was used. The analysis of miRNA expression through total RNA extraction and quantitative reverse transcription polymerase chain reaction (RT-qPCR) assays was conducted in blood samples that had previously been collected in Tempus Blood RNA tubes (ThermoFisher Scientific, Waltham, MA, United States) and stored at −80°C.
Total RNA Isolation
The MagMax Total RNA isolation kit (ThermoFisher Scientific, Waltham, MA, United States) was used to extract total RNA (including miRNAs) as per the recommendations of the manufacturer. Total RNA extraction was conducted on thawed, 3 mL whole blood samples. As this was a semi-automated procedure, the Kingfisher Flex system (ThermoFisher Scientific, Waltham, MA, United States) conducted the RNA washing and elution steps. Quality checks were conducted on the extracted nucleic acid by determining the RNA concentration and purity on the NanoDrop One spectrophotometer (ThermoFisher Scientific, Waltham, MA, United States). An extracted total RNA sample with a 260/280 value between 1.8 and 2.0, and whose concentration was greater than 20 ng/μl passed the quality check and was therefore used for downstream applications like RT-qPCR.
Quantitative Reverse Transcription PCR
In order for miRNA expression to be determined using RT-qPCR, the extracted RNA had to be converted to cDNA first and this conversion was done using the TaqMan Advanced miRNA cDNA synthesis kit (Applied Biosystems, ThermoFisher Scientific, Waltham, MA, United States) as per the manufacturer’s recommendations. In summary, by sequentially conducting poly-A tailing, adaptor ligation, reverse transcription and miR-Amp steps, we converted 2 μl of total RNA into cDNA, which was the starting material in the succeeding PCR step used to determine miRNA expression levels. This PCR analysis was conducted on the QuantStudio 7 Flex real-time PCR instrument (Life Technologies, Carlsbad, CA, United States) using TaqMan miRNA Assay primers (Applied Biosystems, ThermoFisher Scientific, Waltham, MA, United States). The determination of relative miRNA expression in a sample was done using the 2–ΔCt, whilst fold change differences in miRNA expression between the study groups were computed with the use of the 2–ΔΔCt method. In order to normalize miRNA quantification, miR-16-5p (ThermoFisher Scientific, Waltham, MA, United States) was used as the endogenous control (Livak and Schmittgen, 2001).
Statistical Analysis
The Statistical Product and Service Solutions (SPSS) v.26 software (IBM Corp, United States) was used to conduct data analyses. For normally distributed variables, results are reported as count (and percentages), mean (and standard deviation) whilst for asymmetrically distributed variables, the results are reported as median (25th-75th percentiles). In order to compare the mean and median baseline characteristics across blood pressure groups, the analysis of variance (ANOVA) and Kruskal–Wallis tests were used, respectively. Age, gender and BMI-adjusted Spearman’s partial correlations were used to assess the relationship between miRNAs and other cardiovascular risk profile variables whilst multivariable logistic regression models were used to assess the association of miRNAs with screen-detected and known HPT. A p-value less than 0.05 signified statistically significant findings.
Results
Study Participant Characteristics
A total of 1,456 participants were included in this study. A summary of the participants’ characteristics is shown in Table 1. There were 386 (67.4%) females and 187 (32.6%) males in the normotensive group with an average age of 39.3 ± 13.7 years, 209 (68.8%) females and 95 (31.2%) males in the screen-detected HPT group with an average age of 48.9 ± 13.4 years whilst the known HPT group comprised of 478 (82.6%) females and 101 (17.4%) males whose average age was 58.9 ± 11.0 years. As expected for HPT, there was a significant difference in the average age of participants across the three groups (p < 0.001). Serum gamma-Glutamyltransferase (Gamma GT-S) levels differed significantly between the three groups and the other expected differences (body mass index (BMI), age, systolic blood pressure, waist circumference, high density lipoprotein cholesterol (HDL-c) and low-density lipoprotein (LDL) cholesterol by HPT status in the cardiovascular risk profile were apparent between the participant groups as shown in Table 1.
Table 1. Study participant characteristics based on blood pressure status.
Relative MicroRNA Expression
Normalization of relative miRNA expression was done with reference to the expression of miR-16-5p. Overall, the expression of miR-126-3p and 182-5p differed significantly across all blood pressure groups. Figure 1A shows that in participants with known HPT, the relative expression (2–ΔCt) of miR-126-3p (3.266) was significantly higher than that of the normotensives (1.275) and screen-detected hypertensives (1.573), p < 0.001). For miR-182-5p, the relative expression in known HPT was significantly higher (4.360) compared to both normotensives (1.711) and screen detected hypertensives (2.312), both p < 0.001 as shown in Figure 1B. However, for miR-30a-5p, 1299 and 30e-3p, whilst there were significant differences in their relative expression (2–ΔCt) in known HPT (0.068; 0.011 and 0.014) vs normotensives (0.033; 0.004 and 0.006), respectively, (all p < 0.001), that significance was not seen when screen-detected hypertensives were compared to the normotensives, all p ≥ 0.114 as shown in Figures 2A–C.
Figure 1. Comparison of the relative expression of miR-126-3p and 182-5p in the three blood pressure groups. (A) The expression of miR-126-3p was significantly higher in known hypertensives when compared to screen-detected hypertensives (p < 0.001) and the normotensives (p < 0.001). There was also a significant difference in expression when screen-detected hypertensives were compared to the normotensives (p = 0.001). (B) There was a significantly higher expression of miR-182-5p in known hypertensives when compared to screen-detected hypertensives (p < 0.001) and the normotensives (p < 0.001). When screen-detected hypertensives were compared to the normotensives, that significant difference in expression remained (p = 0.002). The relative miRNA expressions were calculated using the 2– ΔCt method and a one-way ANOVA used to compare differences in miRNA expression between groups. The symbol * denotes a statistically significant p-value.
Figure 2. Comparison of the relative expression of miR-30a-5p, miR-1299 and miR-30e-3p. (A) The expression of the miR-30a-5p was significantly higher in known hypertension compared to screen-detected hypertension (p = 0.001) and the normotensives (p < 0.001). However, there was no significant difference in expression between screen-detected hypertensives and normotensives (p = 0.114). (B) There was a significantly higher expression of miR-1299 in known hypertension when compared to screen-detected hypertension (p = 0.002) and the normotensives (p < 0.001). There was no significant difference in expression when screen-detected hypertensives were compared to the normotensives (p = 0.447). (C) The expression of miR-30e-3p was significantly higher in known hypertension when compared to screen-detected hypertension (p < 0.001) and the normotensives (p < 0.001). There was no significant difference in expression when screen-detected hypertensives were compared to the normotensives (p = 0.407). The relative miRNA expressions were calculated using the 2– ΔCt method and a one-way ANOVA used to compare differences in miRNA expression between groups. The symbol * denotes a statistically significant p-value.
Fold Change Computation
When compared to the normotensives, there was at least a 2.2-fold increase in expression of miR-126-3p, 182-5p, 30a-5p, 1299, and 30e-3p in the known HPT group. MiR-126-3p was the most expressed of the five miRNAs, particularly in known hypertensives versus normotensives (fold change = 2.64) and known hypertensives versus screen-detected hypertensives (fold change = 1.99). However, the fold difference in expression of all five miRNAs in the screen-detected hypertensive participants compared to the normotensives was ≤ 1.65-fold, with the lowest fold difference in expression between these two blood pressure groups observed for miR-30e-3p, whose expression was 1.17-fold higher in screen-detected hypertensives compared to the normotensives.
Correlation of MicroRNA Expression With Anthropometric Measurements
There was a significant positive correlation between the expression of the five miRNAs across the blood pressure groups (r ≥ 0.71, p < 0.001), with the highest correlation in expression seen between miR-126-3p and miR-182-5p (r = 0.983, p < 0.001) as shown in Table 2. Detailed correlations between the expression of each miRNA and biochemical parameters are shown in the Supplementary Tables (see Supplementary Tables 1–5). The expression of all miRNAs correlated negatively with waist circumference, with the highest correlation coefficients seen in the screen-detected HPT group for both miR-126-3p (r = −0.748, p < 0.001) and miR-30e-3p (r = −0.729, p < 0.001). There was a significant, though weak correlation between the expression of both miR-182-5p and miR-30a-5p with systolic blood pressure in the known HPT group. However, there was no significant association between the expression of miR-30e-3p, miR-126-3p and 1299 with systolic blood pressure regardless of blood pressure status. Whilst there was no correlation between the expression of any of the five miRNAs and total cholesterol (TC), there was a significant positive correlation with HDL-c, with the highest correlation coefficients seen in the known HPT group with respect to miR-182-5p (r = 0.629, p = 0.001), miR-30a-5p (r = 0.615, p = 0.002), miR-30e-3p (r = 0.608, p = 0.002) miR-126-3p (r = 0.595, p = 0.003) and finally miR-1299 (r = 0.508, p = 0.013). Gamma GT-S also showed a significantly positive correlation with the expression of all but miR-1299 in the screen-detected HPT group.
Table 2. Partial correlations between miRNA relative expression and anthropometric and biochemical parameters.
Multivariable Regression Analysis
The results of multivariable regression analysis are shown in Table 3. With regards to miR-126-3p, the crude odds (age and gender adjusted only) ratio (OR) was 1.16 (95% confidence interval (CI): 1.05-1.27, ‘p = 0.003) for screen-detected HPT, whilst for known HPT, the OR was 1.58 (95% CI: 1.46-1.71, p < 0.001). For miR-30a-5p, the OR was 1.24 (95% CI: 1.00-1.55, p = 0.053) for screen-detected HPT whilst for known HPT, the OR was 1.63 (95% CI: 1.36-1.95, p < 0.001). The crude odds ratio for miR-182-5p was 1.11 (95% CI: 1.05-1.17, p < 0.001) for screen-detected HPT whilst for known HPT, the OR was 1.24 (95% CI: 1.18-1.30, p < 0.001). The associations between the expression of these miR-126-3p, 182-5p and 30a-5p and HPT (both screen-detected and known) remained significant even when the model was adjusted for BMI, TC, triglycerides (TG) and glycated hemoglobin (HbA1c). When further adjusted for duration of known HPT diagnosis, the associations also remained significant for these three miRNAs. As for miR-30e-3p, the only significant association was with known HPT which had an OR of 1.71 (95% CI: 1.46–2.01, p < 0.001) whilst for screen-detected HPT, the OR was 1.07 (95% CI: 0.91–1.26, p = 0.393). There was no significant association between the expression of miR-1299 and screen-detected or known HPT.
Table 3. Multivariable regression analysis of miRNAs for the presence of screen-detected and known hypertension.
Discussion
This study demonstrated a significantly higher expression of miR-126-3p and miR-182-5p in hypertensives (both screen-detected and known) when compared to the normotensives. However, there was no significant difference in the expression levels of miR-30a-5p, 30e-3p, and 1299 in normotensives, relative to screen-detected hypertensives. Although multivariable logistic regressions showed no association between HPT and the expression of both miR-30e-3p and miR-1299, we observed an association between the expression of miR-126-3p, 182-5p and 30a-5p with screen-detected and known HPT, particularly in the latter. Even after adjustment of the crude model for age, sex, BMI, HbA1c, TG, and TC, the associations remained significant.
Angiogenesis and maintenance of vascular integrity are vital blood pressure regulatory processes in which miR-126 has an essential role (Wang et al., 2008). During hypertensive states, there is loss of endothelial cell function and blood perfusion to the capillaries becomes limited, leading to capillary disappearance. This is described as microvascular rarefaction and is a distinctive characteristic of HPT (Humar et al., 2009). Its effects contribute to HPT-related complications such as organ damage and stroke. It has been reported that miR-126, through the stimulation of proangiogenic activity of vascular endothelial growth factor (VEGF) and fibroblast growth factor (FGF), counteracts microvascular rarefaction by encouraging the formation of blood vessels. This is achieved through the quelling of Spred-1 expression, a known inhibitor of angiogenic signaling (Wang et al., 2008). Our study demonstrated a higher expression of miR-126-3p in hypertensive participants when compared to the normotensives, and this finding has also been reported in another study (Liu et al., 2018). It is plausible that the high expression of miR-126 in hypertensives represents a response to lower the blood pressure through promotion of blood vessel formation, whilst repressing the effects of anti-angiogenic Spred-1. Whilst our study demonstrated a higher expression of miR-126-3p in the whole blood of hypertensives compared to the normotensives, another study demonstrated lower expression of miR-126 in hypertensives when compared to healthy controls, albeit in peripheral blood mononuclear cells (PBMCs) (Kontaraki et al., 2014). There was no significant difference in miR-126 expression between hypertensive and normotensive participants in the study by Chen and co-workers (Chen et al., 2018). These discrepancies could in part be explained by the different sample types and participant recruitment criteria used. It has also been demonstrated that anti-hypertensive medications like nebivolol and atenolol affect miRNA expression (Ye et al., 2013). Kontaraki and colleagues excluded anyone on anti-hypertensive medication from their study, while Chen and colleagues’ hypertensive group included any participants who had been on anti-hypertensive therapy for longer than 3 months (Chen et al., 2018).
The efficacy of miR-182 against glioblastoma multiforme, a therapy-resistant cancer of the brain has been previously reported (Kouri et al., 2015). It has been established as an oncogenic miRNA and its interactions in several types of cancers have been reviewed (Wei et al., 2015). The expression of miR-182-5p was significantly higher in both hypertensive groups when compared to the normotensives in our study. However, there is a paucity of studies indicating a role for miR-182-5p in essential HPT. Nonetheless, in a study evaluating the miRNA expression profiles in the placenta of pregnant participants, the expression of miR-182 was found to be elevated in participants with pre-eclampsia, a mid-term complication of pregnancy characterized by proteinuria and HPT, when compared to normal, pregnant controls (Pineles et al., 2007). Hypertension is a known risk factor for coronary artery disease (CAD) (Weber et al., 2016) and in participants with unprotected left main CAD, the plasma expression of miR-182-5p was significantly higher compared to the non-CAD control group and further analyses indicated its high diagnostic power for uncontrolled left main CAD (Zhu et al., 2019). Whether miR-182-5p plays a role in hemodynamic regulation or indeed the pathogenesis of HPT remains to be elucidated.
In our study, the expression of miR-30a-5p was significantly higher in known hypertensives compared to normotensives. The upregulation of miR-30a in hypertensive participants compared to normal controls was also reported in another study and identified as a possible biomarker target for differentiating white coat HPT from essential HPT and normotensives (Huang et al., 2016). Delta-like ligand 4 (dll4) is mainly expressed in the vascular endothelium and exerts its effects through Notch signaling by attachment to its receptor Notch1. The ligand is a key inhibitor of angiogenesis and the expression of miR-30a works against the expression of dll4 and in doing so, promotes angiogenesis (Jiang et al., 2013). The promotion of angiogenesis through the expression of miR-30a-5p can lead to lowering of the blood pressure and could have been the case in our cohort of known hypertensives.
A limitation of our study was that participants making up the three blood pressure groups were not age and gender-matched. Although these variables were adjusted for in the multivariable regression analysis, their effect may not have been completely eliminated. However, the overall large number of participants involved in the project provided sufficient power to the study as each blood pressure group was adequately represented. Whilst the cross-sectional nature of the study precludes inference about causal relationships between miRNA expression and development of HPT, it provided a basis for setting up longitudinal cohorts in which functional studies can be conducted to further clarify the roles that these non-coding ribonucleic acids (RNAs) play in the pathogenesis of HPT.
In conclusion, we demonstrated for the first time in an African population, the differential expression of miRNAs in the whole blood of participants on the basis of their blood pressure status. These miRNAs could form a panel of biomarker targets for HPT diagnostic purposes. Furthermore, the study validated prior findings on miR-126 and the miR-30 family and highlighted their possible involvement in the pathogenesis of HPT, warranting further investigation into these non-coding RNAs as they could offer potential prognostic and therapeutic avenues for cardiovascular diseases.
Data Availability Statement
Some of the datasets used and/or analyzed during the current study are available from the National Centre for Biotechnology Information (NCBI) Sequence Read Archive on the following link: https://www.ncbi.nlm.nih.gov/sra/PRJNA680302. Other datasets presented are not readily available because of the terms of consent to which participants agreed, but are available from the principal investigator of the main study on reasonable request, which should be directed to TM, bWF0c2hhdEBjcHV0LmFjLnph.
Ethics Statement
This investigation was based on the Cape Town Vascular and Metabolic Health (VMH) study, which was approved by the Research Ethics Committees of the Cape Peninsula University of Technology (CPUT) and Stellenbosch University (respectively, NHREC: REC — 230 408 — 014 and N14/01/003). Ethical approval was also obtained for this cross-sectional sub-study from the CPUT Health and Wellness Sciences Research Ethics Committee (CPUT/HW-REC 2019/H7). The study was conducted as per the provisions of the Declaration of Helsinki. All procedures were explained to the participants in their language of choice. Once the participants fully understood their participation, they signed informed consent forms to allow the collection of blood and anthropometric data.
Author Contributions
TM, RE, and AK: conceptualization and funding acquisition. DM, SR, and CW: methodology. DM and SD: formal analysis. DM and CW: investigation. TM: resources. DM, SR, and SD: data curation. DM: writing—original draft preparation. TM, RE, AK, GD, and SR: writing—review and editing. DM, CW, and SR: validation. DM, SD, and TM: visualization. TM and GD: supervision. TM and SD: project administration. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by the South African Medical Research Council (SAMRC), with funds from National Treasury under its Economic Competitiveness and Support Package (MRC-RFA-UFSP-01-2013/VMH Study) and the South African National Research Foundation (SANRF) (Grant no. 115450).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We would like to thank the Bellville South Ward 009 community for partaking in the study as well as the Bellville South Community Health Forum for supporting the engagement with the Bellville South Community.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.710438/full#supplementary-material
References
Bromfield, S., and Muntner, P. (2013). High Blood Pressure: the Leading Global Burden of Disease Risk Factor and the Need for Worldwide Prevention Programs. Curr. Hypertens. Rep. 15, 134–136. doi: 10.1007/s11906-013-0340-9
Campbell, N. R. C., Bovet, P., Schutte, A. E., Lemogoum, D., and Nkwescheu, A. S. (2016). High Blood Pressure in Sub-Saharan Africa: why Prevention, Detection, and Control are Urgent and Important. J. Clin. Hypertens. 17, 663–667. doi: 10.1111/jch.12599
Carretero, O. A., and Oparil, S. (2000). Essential Hypertension. Part I: definition and Etiology. Circulation 101, 329–335. doi: 10.1161/01.cir.101.3.329
Chalmers, J., MacMahon, S., Mancia, G., Whitworth, J., Beilin, L., Hansson, L., et al. (1999). 1999 World Health Organization-International Society of Hypertension Guidelines for the Management of Hypertension. Clin. Exp. 21, 1009–1060. doi: 10.3109/10641969909061028
Chen, S., Chen, R., Zhang, T., Lin, S., Chen, Z., Zhao, B., et al. (2018). Relationship of Cardiovascular Disease Risk Factors and Noncoding RNAs with Hypertension: a Case-Control Study. BMC Cardiovasc. Disord. 18:58. doi: 10.1186/s12872-018-0795-3
Chen, X., Ba, Y., Ma, L., Cai, X., Yin, Y., Wang, K., et al. (2008). Characterization of microRNAs in Serum: a Novel Class of Biomarkers for Diagnosis of Cancer and Other Diseases. Cell Res. 18, 997–1006.
Danaei, G., Finucane, M. M., Lin, J. K., Singh, G. M., Paciorek, C. J., Cowan, M. J., et al. (2011). National, Regional, and Global Trends in Systolic Blood Pressure Since 1980: systematic Analysis of Health Examination Surveys and Epidemiological Studies with 786 Country-years and 5.4 Million Participants. Lancet 377, 568–577. doi: 10.1016/s0140-6736(10)62036-3
Gkaliagkousi, E., Passacquale, G., Douma, S., Zamboulis, C., and Ferro, A. (2010). Platelet Activation in Essential Hypertension: implications for Antiplatelet Treatment. Am. J. Hypertens. 23, 229–236. doi: 10.1038/ajh.2009.247
He, F. J., and Macgregor, G. A. (2007). Blood Pressure is the Most Important Cause of Death and Disability in the World. Eur. Heart J. Suppl. 9, B23–B28.
Huang, Y. Q., Huang, C., Chen, J. Y., Li, J., and Feng, Y. Q. (2016). The Association of Circulating miR-30a, miR-29 and miR-133 with White-Coat Hypertension. Biomark. Med. 10, 1231–1239. doi: 10.2217/bmm-2016-0215
Humar, R., Zimmerli, L., and Battegay, E. (2009). Angiogenesis and Hypertension: an Update. J. Hum. Hypertens. 23, 773–782. doi: 10.1038/jhh.2009.63
Jiang, Q., Lagos-Quintana, M., Liu, D., Shi, Y., Helker, C., Herzog, W., et al. (2013). MiR-30a Regulates Endothelial Tip Cell Formation and Arteriolar Branching. Hypertension 62, 592–598. doi: 10.1161/hypertensionaha.113.01767
Kontaraki, J. E., Marketou, M. E., Zacharis, E. A., Parthenakis, F. I., and Vardas, P. E. (2014). MicroRNA-9 and microRNA-126 Expression Levels in Patients with Essential Hypertension: potential Markers of Target-Organ Damage. J. Am. Soc. Hypertens. 8, 368–375. doi: 10.1016/j.jash.2014.03.324
Kouri, F. M., Hurley, L. A., Daniel, W. L., Day, E. S., Hua, Y., Hao, L., et al. (2015). MiR-182 Integrates Apoptosis, Growth, and Differentiation Programs in Glioblastoma. Genes Dev. 29, 732–745. doi: 10.1101/gad.257394.114
Li, Y., and Kowdley, K. V. (2012). MicroRNAs in Common Human Diseases. Genomics Proteomics Bioinformatics 10, 246–253. doi: 10.1016/j.gpb.2012.07.005
Liu, J., Liu, J., Shi, L., Zhang, F., Yu, L., Yang, X., et al. (2018). Preliminary Study of microRNA-126 as a Novel Therapeutic Target for Primary Hypertension. Int. J. Mol. Med. 41, 1835–1844.
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2–ΔΔCT Method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Matsha, T. E., Hassan, M. S., Kidd, M., and Erasmus, R. T. (2012). The 30-Year Cardiovascular Risk Profile of South Africans with Diagnosed Diabetes, Undiagnosed Diabetes, Pre-diabetes or Normoglycaemia: the Bellville, South Africa Pilot Study. Cardiovasc. J. Afr. 23, 5–11. doi: 10.5830/cvja-2010-087
Matsha, T. E., Kengne, A. P., Hector, S., Mbu, D. L., Yako, Y. Y., and Erasmus, R. T. (2018). MicroRNA Profiling and Their Pathways in South African Individuals with Prediabetes and Newly Diagnosed Type 2 Diabetes Mellitus. Oncotarget 9, 30485–30498. doi: 10.18632/oncotarget.25271
Matshazi, D. M., Weale, C. J., Erasmus, R. T., Kengne, A. P., Davids, S. F. G., Raghubeer, S., et al. (2021). MicroRNA Profiles in Normotensive and Hypertensive South African Individuals. Front. Cardiovasc. Med. 8:645541. doi: 10.3389/fcvm.2021.645541
Mills, K. T., Bundy, J. D., Kelly, T. N., Reed, J. E., Kearney, P. M., Reynolds, K., et al. (2016). Global Disparities of Hypertension Prevalence and Control: a Systematic Analysis of Population-Based Studies from 90 Countries. Circulation 134, 441–450. doi: 10.1161/circulationaha.115.018912
Mills, K. T., Stefanescu, A., and He, J. (2020). The Global Epidemiology of Hypertension. Nat. Rev. Nephrol. 16, 223–237.
Moore, J. (2005). Hypertension: catching the Silent Killer. Nurse Pract. 30, 16–35. doi: 10.1097/00006205-200510000-00005
Nemecz, M., Alexandru, N., Tanko, G., and Georgescu, A. (2016). Role of MicroRNA in Endothelial Dysfunction and Hypertension. Curr. Hypertens. Rep. 18:87.
Pineles, B. L., Romero, R., Montenegro, D., Tarca, A. L., Han, Y. M., Kim, Y. M., et al. (2007). Distinct Subsets of microRNAs are Expressed Differentially in the Human Placentas of Patients with Preeclampsia. Am. J. Obstet. Gynecol. 196, 261.e1–6. doi: 10.1016/j.ajog.2007.01.008
Silambarasan, M., Tan, J. R., Karolina, D. S., Armugam, A., Kaur, C., and Jeyaseelan, K. (2016). MicroRNAs in Hyperglycemia Induced Endothelial Cell Dysfunction. Int. J. Mol. Sci. 17:518. doi: 10.3390/ijms17040518
Spence, J. D., and Rayner, B. L. (2018). Hypertension in Blacks. Individualized Therapy Based on Renin/Aldosterone Phenotyping. Hypertension 72, 263–269. doi: 10.1161/hypertensionaha.118.11064
Taddei, S., Virdis, A., Ghiadoni, L., Sudano, I., and Salvetti, A. (2001). Endothelial Dysfunction in Hypertension. J. Cardiovasc. Pharmacol. 38, S11–S14.
Touyz, R. M., Alves-Lopes, R., Rios, F. J., Camargo, L. L., Anagnostopoulou, A., Arner, A., et al. (2018). Vascular Smooth Muscle Contraction in Hypertension. Cardiovasc. Res. 114, 529–539.
Vidigal, J. A., and Ventura, A. (2015). The Biological Functions of miRNAs: lessons from in vivo Studies. Trends Cell Biol. 25, 137–147. doi: 10.1016/j.tcb.2014.11.004
Wang, S., Aurora, A. B., Johnson, B. A., Qi, X., McAnally, J., Hill, J. A., et al. (2008). The Endothelial-Specific MicroRNA miR-126 Governs Vascular Integrity and Angiogenesis. Dev. Cell 15, 261–271. doi: 10.1016/j.devcel.2008.07.002
Weber, T., Lang, I., Zweiker, R., Horn, S., Wenzel, R. R., Watschinger, B., et al. (2016). Hypertension and Coronary Artery Disease: epidemiology, Physiology, Effects of Treatment, and Recommendations: a Joint Scientific Statement from the Austrian Society of Cardiology and the Austrian Society of Hypertension. Wien. Klin. Wochenschr. 128, 467–479. doi: 10.1007/s00508-016-0998-5
Wei, Q., Lei, R., and Hu, G. (2015). Roles of miR-182 in Sensory Organ Development and Cancer. Thorac. Cancer 6, 2–9. doi: 10.1111/1759-7714.12164
Ye, H., Ling, S., Castillo, A. C., Thomas, B., Long, B., Qian, J., et al. (2013). Nebivolol Induces Distinct Changes in Profibrosis MicroRNA Expression Compared with Atenolol, in Salt-sensitive Hypertensive Rats. Hypertension 61, 1008–1013. doi: 10.1161/hypertensionaha.111.00892
Keywords: hypertension, microRNA, RT-qPCR, Africa, blood pressure, non-coding RNA
Citation: Matshazi DM, Weale CJ, Erasmus RT, Kengne AP, Davids SFG, Raghubeer S, Davison GM and Matsha TE (2021) Circulating Levels of MicroRNAs Associated With Hypertension: A Cross-Sectional Study in Male and Female South African Participants. Front. Genet. 12:710438. doi: 10.3389/fgene.2021.710438
Received: 16 May 2021; Accepted: 13 August 2021;
Published: 14 September 2021.
Edited by:
Dawn Stephens, Technology Innovation Agency (TIA), South AfricaReviewed by:
Said El Shamieh, Beirut Arab University, LebanonZhipeng Liu, Purdue University, United States
Copyright © 2021 Matshazi, Weale, Erasmus, Kengne, Davids, Raghubeer, Davison and Matsha. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Don M. Matshazi, bWF0c2hhemlkQGdtYWlsLmNvbQ==; Tandi E. Matsha, bWF0c2hhdEBjcHV0LmFjLnph; dGFuZGltYXRzaGFAZ21haWwuY29t