Next Article in Journal
Rivastigmine–Benzimidazole Hybrids as Promising Multitarget Metal-Modulating Compounds for Potential Treatment of Neurodegenerative Diseases
Next Article in Special Issue
Pathogenic Variants of SLC22A12 (URAT1) and SLC2A9 (GLUT9) in Spanish Patients with Renal Hypouricemia: Founder Effect of SLC2A9 Variant c.374C>T; p.(T125M)
Previous Article in Journal
Calcium Dynamics of the Ventrolateral Preoptic GABAergic Neurons during Spontaneous Sleep-Waking and in Response to Homeostatic Sleep Demands
Previous Article in Special Issue
CLEC16A—An Emerging Master Regulator of Autoimmunity and Neurodegeneration
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 9
10.3390/ijms24098309
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Risk Effects of rs1799945 Polymorphism of the HFE Gene and Intergenic Interactions of GWAS-Significant Loci for Arterial Hypertension in the Caucasian Population of Central Russia

by
Tatiana Ivanova
,
Maria Churnosova
,
Maria Abramova
,
Irina Ponomarenko
,
Evgeny Reshetnikov
,
Inna Aristova
,
Inna Sorokina
and
Mikhail Churnosov
*
Department of Medical Biological Disciplines, Belgorod State National Research University, 308015 Belgorod, Russia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(9), 8309; https://doi.org/10.3390/ijms24098309
Submission received: 11 March 2023 / Revised: 24 April 2023 / Accepted: 25 April 2023 / Published: 5 May 2023
(This article belongs to the Special Issue Genes and Human Diseases)
Browse Figures
Versions Notes

Abstract

:
The aim of this case-control replicative study was to investigate the link between GWAS-impact for arterial hypertension (AH) and/or blood pressure (BP) gene polymorphisms and AH risk in Russian subjects (Caucasian population of Central Russia). AH (n = 939) and control (n = 466) cohorts were examined for ten GWAS AH/BP risk loci. The genotypes/alleles of these SNP and their combinations (SNP–SNP interactions) were tested for their association with the AH development using a logistic regression statistical procedure. The genotype GG of the SNP rs1799945 (C/G) HFE was strongly linked with an increased AH risk (ORrecGG = 2.53; 95%CIrecGG1.03–6.23; ppermGG = 0.045). The seven SNPs such as rs1173771 (G/A) AC026703.1, rs1799945 (C/G) HFE, rs805303 (G/A) BAG6, rs932764 (A/G) PLCE1, rs4387287 (C/A) OBFC1, rs7302981 (G/A) CERS5, rs167479 (T/G) RGL3, out of ten regarded loci, were related with AH within eight SNP–SNP interaction models (<0.001 ≤ pperm-interaction ≤ 0.047). Three polymorphisms such as rs8068318 (T/C) TBX2, rs633185 (C/G) ARHGAP42, and rs2681472 (A/G) ATP2B1 were not linked with AH. The pairwise rs805303 (G/A) BAG6–rs7302981 (G/A) CERS5 combination was a priority in determining the susceptibility to AH (included in six out of eight SNP–SNP interaction models [75%] and described 0.82% AH entropy). AH-associated variants are conjecturally functional for 101 genes involved in processes related to the immune system (major histocompatibility complex protein, processing/presentation of antigens, immune system process regulation, etc.). In conclusion, the rs1799945 polymorphism of the HFE gene and intergenic interactions of BAG6, CERS5, AC026703.1, HFE, PLCE1, OBFC1, RGL3 have been linked with AH risky in the Caucasian population of Central Russia.

1. Introduction

Arterial (essential) hypertension (AH) is known as a disorder with high blood pressure (BP) [1]. AH is one of the most frequent diseases in the world—in 2015, AH was registered in 1.13 billion people [2]. AH appeared in 30–45% adults, with a predominant prevalence (>60%) in individuals aged >60 years [3]. AH is a serious risk factor for the manifestation of coronary heart disease, stroke, chronic kidney disease and dementia [4]. The risk of coronary heart disease and stroke doubles with an increase in systolic BP (SBP) for every 20 mmHg (starting from 115 mmHg) and diastolic BP (DBP) for every 10 mm Hg (starting from 75 mm Hg) [5]. The largest quantity of SBP-related deaths is caused by coronary heart disease, hemorrhagic and ischemic stroke [6].
The role of hereditary determinants in the formation of BP indicators is highly important and is beyond doubt [7,8,9,10,11,12,13]. The materials gained in the twin/family studies reveal a weighty heredity influence on the BP level, which according to different authors, varies on average within 30–55% [12]. In accordance with the catalog of genome-wide studies, GWAS (https://www.ebi.ac.uk/gwas/search?query=hypertension (accessed on 10 September 2022)), there is information on 118 performed GWAS, as a result of which 586 polymorphisms associated with AH were identified (the information is relevant at the end of 2022). If we also take into account the data obtained in full-exome studies (EAWAS), the number of AH-involved polymorphic loci will exceed 1000 [10]. According to Padmanabhan et al. (2021), at the moment, there is publication information of about 1.5 thousand GWAS SNPs linked with various BP phenotypes (systolic/diastolic/mean/pulse BP) [13]. Considering the accessible data on “ordinary” associative studies of AH from the positions of diverse candidate genes [14,15,16,17,18,19,20,21,22,23], the amount of AH-associated polymorphisms can reach several thousand.
With this mega-large-scale amount of genetic information available to research teams dealing with the problem of the hereditary nature of AH, it is very problematic to give an answer to the question, “what are the specific polymorphisms/genes from the currently known more than several thousand such polymorphisms/genes (of which more than 1.5 thousand GWAS/EAWAS-significant), determine the susceptibility to AH in the population of this region (including residents of Central Russia)?”, since not everyone, even GWAS/EAWAS-significant polymorphism, will determine the susceptibility to AH in the study population. In this regard, replication studies are becoming particularly relevant [24,25,26,27,28], aimed at confirming (or, conversely, refuting) the role of GWAS gene polymorphism in the AH formation in residents of a particular territory having their own characteristics of the genetic “constitution”, the action of environmental factors, intergenic and gene–environmental relationships, etc., predetermining and features of the candidate genes involvement in the disease (AH) formation.
This case-control replication study estimated AH risk linked with GWAS impact loci in genes correlated with AH/BP in Caucasian populations of Central Russia.

2. Results

The distribution of genotypes in both AH (p ≥ 0.098) and control (p ≥ 0.251) groups were in compliance with H-We (Table S1). The genotype GG of the SNP rs1799945 (C/G) HFE was strongly linked with an increased AH risk in both Model 1 (ORrecGG = 2.53; 95%CIrecGG1.03–6.23; precGG=0.043; ppermGG = 0.045; power = 89.65%) and Model 2 (ORrecGG = 2.48; 95%CIrecGG1.02–6.07; precGG = 0.045; ppermGG = 0.046; power = 88.37%) (Table 1). The high identical results on the relationship of rs1799945 (C/G) HFE with AH, obtained by us in both Model 1 (covariates list included in the analysis were the following: BMI; TG; TC; HDL-C; LDL-C; blood glucose; smokers) and Model 2 (in addition to the covariates of Model 1, two more factors such as low physical activity and consumption of fatty food with a high fat content were included in the analysis as confounders), may be due to the fact that the AH-impact phenotypic effects of these two additional covariates included in Model 2 have already been taken into account in the AH-significant phenotypic effects of the covariates of Model 1. Like this, BMI and lipid status indicators (TC; TG; LDL-C; HDL-C) (covariates of Model 1) are strongly correlated with the level of fatty food intake (additional covariate of Model 2), and the BMI parameter (covariate of Model 1) will be significantly linked with the level of physical activity (additional covariate of Model 2). These results allow us to use only a list of Model 1 covariates when studying the intergenic interactions that determine susceptibility to AH (the next stage of our study).
Alongside this, seven SNPs, such as rs1173771 (G/A) AC026703.1, rs1799945 (C/G) HFE,rs805303 (G/A) BAG6, rs932764 (A/G) PLCE1, rs4387287 (C/A) OBFC1, rs7302981 (G/A) CERS5, rs167479 (T/G) RGL3, out of ten regarded loci were related with AH within eight SNP–SNP interaction models (two-(pperm-interaction ≤ 0.047), three-(pperm-interaction ≤ 0.006), and four-(pperm-interaction ≤ 0.001) level (SNP) models were considered, and two, three and three models, respectively, were significant) (Table 2). Three polymorphisms such as rs8068318 (T/C) TBX2, rs633185 (C/G) ARHGAP42, and rs2681472 (A/G) ATP2B1 were not linked with AH.
Importantly, firstly, two loci such as rs805303 (G/A) BAG6 and rs7302981 (G/A) CERS5 were engaged in the maximum number (eight (100%) and six (75%), respectively) of AH-involved intergenic interactions models (Table 2). The pairwise SNP–SNP combination between rs805303 (G/A) BAG6 and rs7302981 (G/A) CERS5 has been priority in determining the susceptibility to AH (included in six out eight SNP–SNP interaction models (75%) (Table 2) and described 0.82% AH entropy (Figure 1)) in comparison with both individual loci (only 0.02–0.36% of entropy is determined (Figure 1)) and other loci interactions (defining less than 0.48% of entropy) (Figure 1).
Secondly, four-SNP models such as rs7302981 (G/A) CERS5 × rs805303 (G/A) BAG6 × rs1173771 (G/A) AC026703.1 × rs167479 (T/G) RGL3 have the highest Wald indicators (47.36, pperm-interaction < 0.001), which “show” its leading value in predisposition to AH (Table 2). Thirdly, thirty-six AH-involved genotype combinations were modeled (Table S2), among which the largest effect (the biggest values (+/−) of the regression coefficient beta) has several up-susceptibility and down-susceptibility combinations such as rs7302981-AA CERS5 × rs805303-GG BAG6 × rs1173771-AA AC026703.1 (beta = 2.34 p = 0.034), rs7302981-GA CERS5 × rs805303-GG BAG6×rs1173771-AA AC026703.1 (beta = 2.38 p = 0.022) and rs7302981-GG CERS5 × rs805303-AA BAG6 × rs1173771-GG AC026703.1 × rs167479-GG RGL3 (beta = −2.70 p = 0.021), rs7302981-GG CERS5 × rs1799945-CC HFE × rs805303-GG BAG6 × rs167479-GG RGL3 (beta = −2.15 p = 0.050).

2.1. Functional Annotation of AH-Associated SNPs

2.1.1. Non-Synonymous and Epigenetic-Significant Loci

There have been four missense SNPs presented (rs7302981 (G/A), rs1799945 (C/G), rs167479 (T/G) and rs1046089 (G/A)) among 88 examined polymorphic genetic variants (7 AH-associated loci and 81 linked with them) leading to amino acid (AA) permutations (Cysteine9Arginine, Histidine63Aspartic acid, Proline162Histidine and Arginine1740Histidine respectively) in 4 AH-causal proteins (CERS5, HFE, RGL3 and PRRC2A) (Table S3). These AA substitutions are characterized by SIFT “deleterious” (Proline162Histidine, Arginine1740Histidine), “tolerated” (Cysteine9Arginine, Histidine63Aspartic acid) predictive parameters, “probably damaging” (Histidine63Aspartic acid, Proline162Histidine, Arginine1740Histidine), and “benign” (Cysteine9Arginine) PolyPhen predictive parameters (Table S3).
Among eighty-eight possible AH loci, 11.36% of SNPs (n = 10) were located in conservative nucleotide DNA sequences, 56.81% (n = 50) were placed in introns, 5.68% (n = 5) were localized in five gene exons such as COX14/LASS5, PRRC2A, HFE, RGL (non-synonymous variants), GPD1 (synonymous variant), and 4.55% (n = 4) were disposed in 3′-UTR (COX14/LASS5) and 5′-UTR (OBFC1) gene regions (Table S4). The bioinformatics data depict the presence of analyzed loci in the places of the supposed enhancers (n = 39; 44.32%) and promoters (n = 17; 19.32%), DNA-ase hypersensitive sites (DHsites) (n = 31; 32.33%), the probable areas of binding to transcription factors (TrF) (n = 74; 84.91%) and regulatory proteins (n = 16; 18.18%) (Table S4). In general, 88 possible polymorphic AH loci can potentially participate in the epigenetic regulation of the activity of 22 adjacent genes such as RP11-411N4.1, RGL3, PRRC2A, PLCE1, OBFC1, HFE, NPR3, GPD1, LY6G5B, GPD1, LIMA1, CSNK2B, HIST1H1T, COX14, HIST1H4C, CERS5, HIST1H2BC, AC026703.1, LASS5, C12orf62, BAG6, HIST1H2AC) (Table S4). A weighty number of AH-involved loci exhibit their epigenetic effects in AH target organs, such as the adult/fetal heart, aorta. For instance, in the conjectural promoter/enhancer gene regions in the fetal heart, there are loci rs4387287 (C/A) OBFC1 and rs7302981 (G/A) CERS5 in the right atrium, for left ventricles and aorta, rs4387287 (C/A) OBFC1, and for the right ventricle, rs1173771 (G/A) AC026703.1, rs805303 (G/A) BAG6, rs4387287 (C/A) OBFC1, rs7302981 (G/A) CERS5.
The largest number of strongly coupled loci (n = 34), including polymorphisms with any regulatory effects (n = 32), has been registered for our locus rs7302981 (G/A) CERS5 (Table S4). Two of these loci such as rs78594839 and rs10538142 were located in the DNA regions that interacted with the largest number of TrF (11 each) such as Cart1, Dbx1, Evi-1, Foxp1, HDAC2, HNF1, Hoxa10, Hoxa5, Hoxb13, Sox, Zfp105 and Fox, Foxa, Foxd3, Foxf1, Foxi1, Foxj1, Foxj2, Foxp1, Foxq1, Sox, Zfp105. Three-fourths of all polymorphic loci strongly linked to rs7302981 (G/A) CERS5 (26 of 34, 76.47%) were functionally significant (located in the DNA positions of the presumable promoters/enhancers/DHsites) in the AH target organs (adult/fetal heart, aorta). According to the information on chromatin states (HaploReg data based on 25-state model using 12 imputed marks), several polymorphic loci such as rs836179, rs836180, rs7967705 have the strongly pronounced epigenetic impact in all examined AH target organs (fetal heart, right atrium, left and right ventricles, aorta of adult).

2.1.2. Plausible Gene Expression (eQTL) and Splicing (sQTL) Regulatory Potential of AH-Involved SNPs

The four AH-associated loci, such as rs7302981, rs4387287, rs805303, rs1799945, independently and due to strongly linked loci (33 SNPs out of 81 strongly linked loci were cis- and trans-eQTL significant, 40.74%) determine the expression of nine genes in human blood (peripheral) (LY6G5C, AIF1, LIMA1, HCP5, SLK, HSPA1B, LASS5, TRIM38, and ALAS2) (blood eQTL browser materials are presented in Tables S5 and S6).
At that, strongly pronounced gene expression was found with regulatory potential of all 7 AH-correlated loci and 74 coupled SNP (91.36% among all 81 examined LD SNP) in relation to 78 genes in other various organs (ABHD16A, AIF1, APOM, AQP5, ASIC1, ATF1, ATF6B, ATP6V1G2, BAG6, BTN2A3P, C4A, C4B, C6orf48, CCHCR1, CERS5, CLIC1, COX14, CSNK2B, CTC-510F12.3, CYP21A1P, CYP21A2, DDAH2, DIP2B, DXO, GPD1, ANK1, GSTO1, GUSBP2, HCG22, HCP5, HDAC1P1, HFE, HIST1H3E, HLA-B, HLA-DQB1, HLA-DRB5, HLA-DRB6, HLA-S, HSPA1A, LARP4, LIMA1, LY6G5BLY6G5C, LY6G6C, LY6G6D, MICB, LY6G6E, LY6G6F, MIR6891, MPIG6B, NCR3, NEU1, NPR3, POU5F1, PRRC2A, RGL3, RNF5, RP11-457M11.5, RP11-541N10.3, RP3-405J10.3, RP4-605O3.4, SH3PXD2A, SH3PXD2A-AS1, SLC17A1, SLC17A3, SLK, SMARCD1, STK19, STK19B, STN1, TRIM38, U91328.19, UQCRHP1, VWA7, WASF5P, XXbac-BPG248L24.12, ZBTB12, ZNF322) (GTExproject information is showed in Tables S7 and S8). It is important to highlight the serious capability to eQTL regulation of considered polymorphisms, as in the AH target organs such as arterial vessels (VWA7, HFE, ABHD16A, LY6G5B, CYP21A1P, BAG6, ATF1, HLA-DRB5, SMARCD1, LY6G5C, C4A, COX14, RP4-605O3.4, RP4-605O3.4, LY6G5B) and the heart (CSNK2B, C4A, HLA-DRB5, CERS5, CYP21A1P, HIST1H3E, DDAH2, RP4-605O3.4, LY6G5B, STN1), as well as in organs involved in the AH biology such as various parts of the brain (cortex; basal ganglia; hypothalamus; pituitary gland; black substance, etc.) (ABHD16A, CERS5, CYP21A1P, DDAH2, HIST1H3E, LY6G5B, LY6G5C, MPIG6B, RP4-605O3.4, SH3PXD2A), adrenal gland (NPR3, SLC17A1, SLC17A3, GUSBP2, DDAH2, LY6G5B, ABHD16A, RP4-605O3.4, COX14), thyroid (HFE, TRIM38, LY6G5C, LY6G5B, HLA-DRB5, CCHCR1, CYP21A1P, C4A, RNF5, ABHD16A, LY6G6F, RP4-605O3.4, CERS5, LIMA1, COX14, ASIC1), adipose (ABHD16A, ASIC1, BTN2A3P, C4A, C6orf48, CCHCR1, CERS5, COX14, CYP21A1P, CYP21A2, DDAH2, HFE, HLA-B, HLA-DRB5, LY6G5B, LY6G5C, RP11-541N10.3, RP4-605O3.4, STK19B, STN1, U91328.19, VWA7), skeletal muscle (HIST1H3E, LY6G5C, LY6G5B, DDAH2, GPANK1, HLA-DRB5, ATF6B, APOM, C4A, RNF5, CSNK2B, PRRC2A, CYP21A1P, CLIC1, STN1, RP4-605O3.4, COX14, ATF1, DIP2B), blood (HIST1H3E, LY6G5C, LY6G5B, HLA-DRB5, C4B, PRRC2A, CYP21A2, VWA7, C4A, CYP21A1P, AIF1, C6orf48, RP4-605O3.4, ATF1, LIMA1), etc. (Tables S7 and S8).
We observed the connection of the genetic polymorphisms under consideration (7 AH causal loci/65 SNP in LD) with the intron splicing regulation of 32 genes (PRRC2A, AIF1, ATF6B, BAG6, C6orf48, CCHCR1, CYP21A1P, CYP21A2, DDX39B, FLOT1, GPANK1, LSM2, HLA-DQA1, HLA-DRB1, HLA-DRB5, HLA-DRB6, LST1, LY6G5B, LY6G5C, LY6G6C, MICA, STK19, STK19B, VARS, ATF1, CERS5, COX14, FAM186A, RP4-605O3.4, HFE, SMARCD1, SH3PXD2A-AS1) including in disease target organs such as arterial vessels (aorta, coronary artery, etc.) (BAG6, HLA-DRB1, HLA-DRB5, HLA-DRB6LSM2, STK19B, GPANK1, ATF6B), heart (GPANK1, HLA-DRB1, HLA-DRB5, HLA-DRB6, STK19B, LY6G5C) and organs significant for AH pathogenesis: brain (cortex; basal ganglia; pituitary) (LY6G5C, BAG6), adrenal gland (CCHCR1, BAG6, CYP21A2, CYP21A1P), thyroid (GPANK1, HLA-DRB1, STK19, HLA-DRB5, CCHCR1, HLA-DRB6, STK19B, BAG6, FLOT1), adipose (HLA-DRB6, BAG6, HLA-DRB1, AIF1, HLA-DRB5), skeletal muscle (GPANK1, HLA-DRB5, BAG6, HLA-DRB6, HLA-DRB1, HLA-DRB5, CCHCR1), and blood (HLA-DRB5, BAG6, HLA-DRB6, AIF1, GPANK1, HLA-DRB1, LY6G5C) (Tables S9 and S10).

2.1.3. Pathway Analysis of AH-Associated Genes

Based on the 101-gene (previously, these genes functionally related to 88 disorder-involved polymorphic genetic variants (7 AH-associated/81 linked with them) were identified) enrichment analysis results, an extremely large number (about 150) of biological pathways was discovered (Table S11). Among the pathways, processes related to the immune system prevailed; the greatest statistical significance was shown by such pathways as major histocompatibility complex (MHC) protein (PC00149, pfdr = 4.18 × 10−11), antigen processing and presentation (GO:0019882, pfdr = 4.74 × 10−10), positive regulation of immune system process (GO:0002684, pfdr = 3.66 × 10−7), adaptive immune response based on somatic recombination of immune receptors built from immunoglobulin superfamily domains (GO:0002460, pfdr = 2.46 × 10−6), lymphocyte mediated immunity (GO:0002449, pfdr = 2.73 × 10−6), etc.
The simulated network of intergenic interactions at predisposition to AH (Figure 2) is based on the following hypothetical “mechanisms”: co-expression (48.87%), physical interactions (30.06%), common protein domains (7.78%), joint localization (7.56%), and forecast interactions (5.72%). Among the “main” interacting genes “additionally” included in the genetic network (besides the 101 AH-recognized genes), the first rank positions are occupied by the two genes LSM3 and CSNK2A1. The most pronounced interactions (the physical interactions weight index is 1) were discovered for such gene pairs as HLA-DQB1/HLA-DQA1, LSM3/LSM2, CSNK2A1/CSNK2B (Table S12).
When investigating AH-involved protein–protein interaction (PPI) networks and functional enrichment analysis performed using the STRING online resource, the following data were obtained. Figure 3 shows that the AH-involved proteins significantly interact with each other (PPI enrichment p value < 1.0 × 10−16) and such PPI as Ly-6 antigen/uPA receptor-like/acetylcholine receptor regulator activity (pFDR = 3.16 × 10−16), Tenascin-X/Proline-rich protein 3 (pFDR = 6.23 × 10−8) and proline-rich protein 3/TRIM10/RING-HC finger (pFDR - = 6.23 × 10−8) are of paramount importance in these processes. These PPI are carried out on the basis of protein domains such asimmunoglobulin C1-set domain (pFDR = 9.59 × 10−8), Class I histocompatibility antigen, domains alpha 1 and 2 (pFDR = 0.003) and Class II histocompatibility antigen, beta domain (pFDR - = 0.013).
Among the PPI, three clusters of functionally related proteins can be distinguished (Figure 4D, clusters are indicated by different colors, red, green, and blue). The first PPI cluster (Figure 4A (indicated in red) includes 30 proteins such as HLA-B, HLA-DQB1, HLA-DRB5, HLA-DRB6, HLA-S, HSPA1A, LY6G5BLY6G5C, LY6G6C, LY6G6D, MICB, ZBTB12, AIF1, LST1, MICA, etc.) was represented mainly by interactions of proteins associated with the innate/adaptive immune responses: Interferon-gamma-mediated signaling pathway (GO:0060333, pFDR = 1.97 × 10−5), Antigen processing and presentation of peptide antigen (GO:0048002, pFDR = 5.19 × 10−5), Regulation of immune system process (GO:0002682, pFDR = 0.002), Immune response-activating cell surface receptor signaling pathway (GO:0002429, pFDR = 0.009), Cytokine-mediated signaling pathway (GO:0019221, pFDR = 0.02), etc. The second PPI cluster (Figure 4B (indicated in green) also includes 30 proteins such asHSPA1B, HSPA1A, C4A, C4B, CLIC1, BAG6, FLOT1, etc.) was characterized by interactions of heat shock proteins, complement systems proteins, etc., in which such pathways as Blood microparticle (GO:0072562, pFDR = 0.002), Extracellular exosome (GO:0070062, pFDR = 0.004), Misfolded protein binding (GO:0051787, pFDR = 0.04), etc., are involved. The third PPI cluster (Figure 4C (indicated in blue) includes 15 proteins such as HIST1H2BC, HIST1H2AC, HIST1H4C, etc.) is distinguished by the interaction of various histone proteins based on such Reactome pathways as recognition and association of DNA glycosylase with sites containing an affected pyrimidine (HSA-110328, pFDR = 0.005) and purine (HSA-110330, pFDR = 0.005), Cleavage of the damaged pyrimidine (HSA-110329, pFDR = 0.005) and purine (HSA-110331, pFDR = 0.005), Meiotic synapsis (HSA-1221632, pFDR = 0.005), Packaging Of Telomere Ends (HSA-171306, pFDR = 0.005), Pre-NOTCH Transcription and Translation (HSA-1912408, pFDR = 0.005), Formation of the beta-catenin:TCF transactivating complex (HSA-201722, pFDR = 0.005), PRC2 methylates histones and DNA (HSA-212300, pFDR = 0.005), Condensation of Prophase Chromosomes (HSA-2299718, pFDR = 0.005), etc.

3. Discussion

In the present study associations have been replicated of the rs1799945 (C/G) HFE and intergenic interactions of seven GWAS-significant loci for AH/BP (rs1173771 (G/A) AC026703.1, rs1799945 (C/G) HFE, rs805303 (G/A) BAG6, rs932764 (A/G) PLCE1, rs4387287 (C/A) OBFC1, rs7302981 (G/A) CERS5, rs167479 (T/G) RGL3) with AH in the Caucasian population of Central Russia. AH risk SNPs (with more 80 proxy variants) are supposedly functionally efficient with respect to the 101 genes implicated in various immune system pathways. Three studied SNPs (rs8068318 (T/C) TBX2, rs633185 (C/G) ARHGAP42, rs2681472 (A/G) ATP2B1) did not confirm the association with AH.
Among the 10 GWAS AH/BP polymorphisms studied, we confirmed independent associations with AH for only one locus—rs1799945 (C/G) HFE gene (OR = 2.53 for genotype GG). There are quite numerous literature data (based on the results of GWAS) convincingly indicating the involvement of the SNP rs1799945 (C/G) HFE in the BP levels formation [29,30,31,32,33,34,35] and susceptibility to AH development [29,31].
In three GWAS, a high BP level and an increased risk of AH have been marked by the genetic variant G rs1799945 [29,31,34], and in three other GWAS, lower BP values (systolic/diastolic/mean/pulse blood pressure) were marked by a reference allele C of this polymorphism [30,31,32,33]. In our work, the risk influence on the AH development is exerted by the genotype GG rs1799945 (C/G) HFE (OR = 2.53), which fully corresponds to the above-mentioned literature GWAS data on the role of allelic variants of this locus (G-risky vs. C-protective) in the formation of BP and AH in various populations of the world.
Rather interesting results about the biological significance of the rs1799945 (C/G) HFE in humans of European ancestry, which can to some extent give a biomedical explanation of the relationship of this polymorphism with susceptibility to AH, were obtained in the study Gill et al. [36].
The authors, on the one hand, on the basis of GWAS data, showed a strong link between the locus rs1799945 (C/G) HFE and the human organism iron status (these were analyzed such serum parameters as levels of the iron, transferrin, ferritin, and saturation of transferrin), on the other hand, using a phenome-wide association study with mendelian randomization (MR-PheWAS analyses), established a substantial influence (causal effect) of the status of iron in organisms at risk of anemia and hypercholesterolemia development [36]. A direct link between hypercholesterolemia and a high AH risk (and in general, a high risk of morbidity/mortality associated with cardiovascular diseases) has been known for a long time and is currently beyond doubt (including in the countries of Eastern Europe) [2,37]. High levels of total cholesterol (TC), triglycerides (TG), low-density lipoprotein cholesterol (LDL-C), and low levels of high-density lipoprotein cholesterol (HDL-C) in AH patients in comparison with the control were also found in the studied sample (Table 3), which, firstly, confirms the significant role of hypercholesterolemia in AH development (risk factor); secondly, it gives reason to assume that one of the significant mechanisms underlying the involvement of rs1799945 (C/G) HFE in AH formation in the population studied by us may be iron status and hypercholesterolemia associated with it.
The data derived by us in silico also demonstrate a serious functional potential of the rs1799945 (C/G) HFE in the organism (it and 7 loci strongly linked to it have been functionally important for 15 genes, including replacement of amino acid Histidine by Aspartic acid in the 63 position of the HFE protein, epigenetic changes in 5 genes (HIST1H2BC; HIST1H2AC; HIST1H1T; HIST1H4C; HFE), eQTL/sQTL influences on 11/1 genes (ALAS2; BTN2A3P; GUSBP2; HFE; HIST1H3E; RP11-457M11.5; SLC17A1; SLC17A3; TRIM38; U91328.19; ZNF322/HFE)), which can also justify its involvement in AH susceptibility. For example, histone genes, such as H2A, H3E, H1T, H2B, H4C, etc., are functionally associated with the rs1799945 locus and are of paramount importance in the regulation of chromatin structure resulting in the modification of DNA “activity” (suppress/activate gene transcription): H2A variants hold the positions of entry/exit along the nucleosomal DNA wrap and thus control the “availability” of DNA; various fractions of H3 histone proteins and their modifications are known “markers” of functionally active DNA regions (promoters/enhancers); the H1T variant is a linker histone and “coordinates” chromatin packaging [38]. In another example, using the cDNA library of the heart of an early human embryo, a new zinc finger gene called ZNF322 was isolated, which, as the authors showed, is expressed both at various embryonic stages (from 80 days to 24 weeks; thus, it is in numerous tissues of an adult [39]. The authors have shown that ZNF322 is a “transcription activator” (via MAPK signaling pathways) of reporter genes such as SRE, AP-1, which have been important for AH pathobiology [39].
Using the in silico approach widely utilized in modern genetic research, we established the multifarious pleiotropic functional actions (chromatin changes; eQTL/sQTL) of 7 AH-involved loci (and 81 linked SNPs) in relation to 101 genes. Furthermore, for a considerable number of genes (more 20 genes), these influences were registered by us in AH target organs such as arterial vessels (aorta, coronary artery, etc.) (OBFC1, HLA-DRB5, HLA-DRB1, HLA-DRB6, LY6G5B, VWA7, C4A, HFE, ABHD16A, CYP21A1P, SMARCD1, LY6G5C, COX14, GPANK1, RP4-605O3.4, RP4-605O3.4, BAG6, ATF1, LSM2, STK19B, ATF6B) and adult/fetal heart (OBFC1, HLA-DRB1, HLA-DRB6, HLA-DRB5, CSNK2B, C4A, CERS5, CYP21A1P, HIST1H3E, DDAH2, RP4-605O3.4, STN1, LY6G5B, STK19B, LY6G5C, GPANK1). The phenotypic functionality of the aforementioned genes may be responsible for the pathophysiology of AH. Thus, for instance, the OBFC1 gene (also called STN1 subunit of CST complex) encodes one of the subunits of alpha accessory factor that stimulates the enzyme activity and thereby initiates DNA replication; as well, this protein is important for telomere-linked complexes and telomere regulation mechanisms (it provides telomere length homeostasis by inhibiting telomerase, recruits and activates the corresponding DNA polymerase, facilitates repair and initiates DNA replication, etc.) [40,41]. Said et al. detected the strongest association of genetically determined telomere length (including OBFC1 gene) with cardiovascular disease and hypertension (based on Mendelian randomization data of UK Biobank 134,773 individuals) [42]. The strong effect of expression-involved OBFC1 polymorphic variants (including eQTL in adventitial tissue of aorta) on cardiovascular disease risk has been shown in previous publication data [43]. Considering our in silico data, several genes of the body immune system (HLA-DRB1, HLA-DRB6, HLA-DRB5) and the extensive pathways related with them associated with the immune responses (MHC protein-involved reactions, immune system process regulation, antigen processing and presentation, etc.) are among the causal factors of hypertension. Pronounced changes in the immune system of AH patients (increased proinflammatory interleukins and other cytokines (tumor necrosis factors, chemokines, etc.) plasma levels, heat shock proteins overexpression, etc.), a substantial contribution of various immune-dependent processes to the development/aggravation of the disease (innate/adaptive immunity activation, changes of the immune reactivity, B- and T- cells pro- and anti- inflammatory responses, manifestations of autoantigenic reactions (due to heat shock proteins, isoketal-modified proteins, etc.), cascade of cytokine responses, etc.), and the need to correct the immune imbalance during hypertension therapy have been noted in numerous previous studies [44,45,46,47]. Rodriguez-Iturbe et al. considers the formation and exacerbation of hypertension as a step-by-step process of involving a variety of mechanisms of the immune system (episodic formation of danger-linked molecular substations and Toll-like receptor expression, activation of the innate immune system and the appearance of inflammatory reactions (infiltration, etc.) in target organs (kidney, vessels, etc.), involvement adaptive immunity reactions in the renal and vascular inflammation, the appearance of an imbalance in inflammatory/anti-inflammatory responses, etc.) and indicate the need to identify genetic factors/traits of hypertension related to the immune response [44]. Our work confirms the paramount importance of genetic factors (including certain AH-impact GWAS polymorphisms) associated with the multiple immune system pathways in AH formation in the Caucasian population of Central Russia.
It is extremely interesting that in a previously conducted genetic research aimed at finding associations of GWAS AH candidate genes with pre-eclampsia in a sample that included 452 pregnant women with pre-eclampsia and 498 pregnant women without pre-eclampsia (the same ten GWAS AH/BP risk loci panel was studied in the same Caucasian population of Central Russia), the significant role of rs1799945 (C/G) HFE (OR = 2.24) in the development of this pregnancy complication was detected [48]. Herewith, it is important to note the complete coincidence of the data obtained for AH (this study) and pre-eclampsia (obtained by us earlier [48]) in orientation (the risk value of the polymorphic variant G) and the associated strength of the polymorphism rs1799945 (C/G) HFE (OR indicators were approximately the same in value and amounted to OR = 2.53 for AH and OR = 2.24 for pre-eclampsia). These facts may indicate, on the one hand, the “reliability” of the identified associations obtained from completely different samples in the same ethno-territorial group of the Russian population and, on the other hand, the proven risk role of GWAS-impact rs1799945 (C/G) HFE in the development of hypertensive conditions (AH as an independent disease, and arterial hypertension is the main symptom of pre-eclampsia) in the Caucasian population of Central Russia, which opens up good prospects for further use of this polymorphism with a prognostic purpose in the practical medicine in this territory of Russia. Meanwhile, there is also an obvious need for additional evidence of the effectiveness of the rs1799945 (C/G) HFE polymorphism in other human diseases correlated with hypertensive manifestations (for example, kidney disease, metabolic syndrome, etc.) in the Caucasian population of Central Russia, which may be the subject of further research.

4. Materials and Methods

4.1. Study Subjects

When planning this work, the number of samples (patient/control) necessary for the study were determined using the Genetic Association Study (GAS) Power Calculator software (online source: http://csg.sph.umich.edu/abecasis/gas_power_calculator/ (accessed on 18 November 2022)) (a multiplicative model of the disease was considered). Taking into account the prevalence of AH among adults, according to the literature [1], on average, about 30–45% with the required research power of 80% and 5% error of the 1st kind (α = 0.05) to identify differences in the frequencies of polymorphic variants between patients and control at the level of OR = 1.23–1.30, the total sample number (patients and control) should be at least 1100 subjects with the prevalence of polymorphic genetic variants among the population ≥10%.
The sample of the present case-control study consisted of 1405 unrelated Caucasian participants (Russian origin; born in Central Russia [49,50]) including 939 AH patients and 466 AH free individuals. The studied subjects were recruited during the 2013–2016 period at the Cardiology Department of St. Joasaph Belgorod Regional Clinical Hospital. All procedures in this study were performed following the tenets of the Declaration of Helsinki. The study protocol was approved by the of the Human Investigation Committee (Ethics Committee) of Belgorod State University. All the participants were fully informed of the purpose and procedures, and written consent was obtained from each participant. The AH diagnosis was verified by certified cardiologists according to the WHO/ESC/ESH recommendation [1] (this was described in detail by us earlier [21]). AH was defined as SBP ≥ 140 mmHg and/or DBP ≥ 90 mmHg (office parameters) [1]. All AH patients had ≥1 year AH clinical anamnesis, and 81.79% received antiAH drugs. Control subjects were AH free (SBP < 140mmHg and DBP < 90mmHg), did not receive antiAH drugs and did not have pronounced metabolic (type 2 diabetes mellitus) and cardiovascular (coronary artery disease) disorders. Individuals with severe autoimmune/allergic/oncological/hematological disorder were excluded [51]. Data characterizing diet and lifestyle were obtained for all subjects. In accordance with WHO/FAO Expert recommendations, “low fruit/vegetable consumption” was considered as daily consumption of less than 400 g of fruits and vegetables, “low physical activity” was evaluated as average weekly physical activity of medium intensity (total at home and at work) less than 2.5 h, “high fatty foods consumption” was estimated to be the share of the average daily consumption of fatty foods from the total amount of food consumed ≥10%, and “high sodium consumption” was considered to be the daily consumption of salt in the volume of a teaspoon or more (≥5 g) [52].
Table 3 presents phenotypic characteristics of the AH and AH-free participants previously received by us [21]: AH patients had high indicators of BMI, blood glucose, TC, TG, LDL-C, smokers and low parameters of HDL-C vs. AH-free individuals (p < 0.001). These data gave us a reason to use the above indicators as confounders in genetic analysis (Model 1). Besides this, the percentage of individuals with low physical activity and preferring fatty foods was higher among AH patients (p < 0.001) (Table 3). These two AH-significant risk factors were additionally included as covariates in Model 2 (together with all covariates of Model 1).

4.2. Experimental Genetic Analysis (DNA Isolation; SNPs Selection; SNPs Genotyping)

Five milliliters of the venous blood was drawn from the ulnar vein and collected into tubes containing 100 μL of 10% EDTA [53]. High molecular weight genomic DNA was extracted from peripheral blood leukocytes, using the standard (phenol/chloroform) protocols [54].
For this study, we selected 10 loci taking into account special criteria such as (1) previously GWAS-linked with BP/AH in Caucasians; (2) possessed impact functionality (evaluated in silico by HaploReg programme [55]), and (3) had significant polymorphism (the frequency of the minor allele was at least 10% among Europeans (HaploReg data [55] and the data of previously conducted studies [48,56,57] were taken into account)). The ten common SNPs (rs167479 of RGL3, rs8068318 of TBX2, rs2681472 of ATP2B1, rs7302981 of CERS5, rs633185 of ARHGAP42, rs4387287 of OBFC1, rs932764 of PLCE1, rs805303 of BAG6, rs1799945 of HFE, rs1173771 of AC026703.1) were chosen for this study based on the above special criteria. All ten SNPs were BP-associated in Europeans, and all 10 SNPs were AH-linked: eight loci were AH-correlated in Europeans and two SNPs (rs2681472 ATP2B1; rs4387287 OBFC1) were AH-significant in the cohort subjects with European prevalence (more 85% out all participations) (Table S13). Nine polymorphic variants among studied 10 SNPs where AH/BP was associated in at least two GWAS (only the locus rs4387287 OBFC1 was AH/BP-associated in one GWAS) (Table S13). All 10 loci selected for the study were found with a frequency of 14% or more among Europeans (Table S14) and were functional (Table S14).
Well known and widely used in genetic research, the TaqMan probe method of Polymerase Chain Reaction was used for genotyping DNA on the examined SNPs. The allele-detection process was performed on a CFX96 Real-Time System (Applied Biosystems, Foster City, CA, USA) to determine the allelic discrimination [58]. The genotyping procedure was carried out at the department of Medical Biological Disciplines (Laboratory of “Human Molecular Genetics”) of Belgorod State National Research University. We paid special attention to the quality control of genotyping: duplicate samples of DNA of participants (≈4–6% out all samples) (independent internal positive control) and samples containing the reaction mixture but not the DNA content (independent internal negative control) were interspersed throughout the plates used for genotype analysis for quality control purposes [59,60]. Importantly, the status of “AH/AH free” of the subjects and the positions of positive/negative independent internal controls on the plates were “hidden” for laboratory personnel during the entire procedure of experimental genetic research. The concordance rates for quality control samples (positive/negative internal controls) were 100% for all assays. An additional measure of quality control of genotyping was the assessment of the compliance of the observed distribution of genotypes with the expected parameters when fulfilling the Hardy–Weinberg law (we performed this procedure at the next stage of our study—statistical analysis of genetic data). The implementation of the Hardy–Weinberg law for the loci under consideration was an additional argument indicating the sufficient quality of the genotyping performed.

4.3. Statistical Analysis of Genetic Data

The samples included in the study (939 AH and 466 AH free) can identify differences in the frequencies of minor alleles of studied SNPs (with an estimated prevalence of the SNPs minor allele among Europeans of 14–46% (Table S14) (HaploReg data [55]) between the AH and AH free groups at OR = 1.25–1.37 (additive model), OR = 1.42–1.44 (dominant model), OR = 1.47–3.10 (recessive model) (calculations were performed in the Quanto program (http://hydra.usc.edu/gxe, 2006 (accessed on 18 November 2022)), at power = 80%; α = 0.05 for 2-sided test).
For each SNP, the Hardy–Weinberg equilibrium was assessed by Pearson’s goodness-of-fit χ2 statistical test [61,62]. The degrees of significance of differences in polymorphic variants (genotype/allele and genotype combinations frequencies) of studied genes (four inheritance models were considered for individual SNPs (allelic; recessive; additive; dominant) [63]; for the general SNP list, intergenic interactions were analyzed [64]) between AH and AH-free cohorts were calculated using gPLINK [65], MB-MDR [66,67], MDR [68,69] packages. The odds ratios (ORgenetic model) and their 95% confidence intervals (CIgenetic model) were obtained by logistic regression [70,71] while adjusting for multiple comparisons (applied permutation testing [72]) and confounding factors (listed above for Model 1 and Model 2 in the “Study subjects” section according to the data in Table 3). Statistical significance was set at 5%, or pperm (pperm-interaction) ≤ 0.05. For AH-associated SNP rs1799945 (C/G) HFE was calculated for statistical power with Quanto tool (http://hydra.usc.edu/gxe, 2006 (accessed on 18 November 2022)).

4.4. SNPs/Gene Predict Functionality/Functions

The well-known bioinformatics information (PolyPhen [73]; SIFT [74], HaploReg [55]; Blood eQTL browser [75]; GTExConsortium [76]; Gene Ontology [77]; GeneMANIA [78]; STRING [79]) was applied to examine in silico an association at the AH-involved loci/genes and high-linked SNPs (used parameter r2 equal or more 0.80 [80,81]) (according to the HaploReg database [55]) with functional prediction effects [82,83].

5. Conclusions

The GWAS-impact for AH/BP polymorphic locus rs1799945 of the HFE gene and intergenic interactions of BAG6, CERS5, AC026703.1, HFE, PLCE1, OBFC1, RGL3 are associated with the risk of developing AH in the Caucasian population of Central Russia. Alongside this, three studied SNPs such as rs8068318 (T/C) TBX2, rs633185 (C/G) ARHGAP42, rs2681472 (A/G) ATP2B1, did not confirm the association with AH.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24098309/s1. References [84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100] are cited in Supplementary Materials.

Author Contributions

Conceptualization, T.I. and I.P.; Data curation, M.A. and I.S.; Formal analysis, E.R., I.A. and I.P.; Project administration, M.C. (Mikhail Churnosov); Writing—original draft, T.I., M.A. and M.C. (Maria Churnosova); Writing—review and editing, M.C. (Mikhail Churnosov), E.R. and I.A. All authors provided final approval of the version to be published. All authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This study was conducted according to the guidelines of the Declaration of Helsinki and was approved by the Local Ethical Committee of the Belgorod State University (21 January 2013, No. 1).

Informed Consent Statement

Informed consent was obtained from all subjects involved in this study.

Data Availability Statement

The data generated in the present study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Williams, B.; Mancia, G.; Spiering, W.; Agabiti Rosei, E.; Azizi, M.; Burnier, M.; Clement, D.L.; Coca, A.; de Simone, G.; Dominiczak, A.; et al. 2018 ESC/ESH Guidelines for the management of arterial hypertension. Eur. Heart J. 2018, 39, 3021–3104. [Google Scholar] [CrossRef] [PubMed]
  2. NCD Risk Factor Collaboration. 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. [Google Scholar] [CrossRef] [PubMed]
  3. Chow, C.K.; Teo, K.K.; Rangarajan, S.; Islam, S.; Gupta, R.; Avezum, A.; Bahonar, A.; Chifamba, J.; Dagenais, G.; Diaz, R.; et al. Prevalence, awareness, treatment, and control of hypertension in rural and urban communities in high-, middle-, and low-income countries. JAMA 2013, 310, 959–968. [Google Scholar] [CrossRef]
  4. Zhou, B.; Perel, P.; Mensah, G.A.; Ezzati, M. Global epidemiology, health burden and effective interventions for elevated blood pressure and hypertension. Nat. Rev. Cardiol. 2021, 18, 785–802. [Google Scholar] [CrossRef] [PubMed]
  5. Lewington, S.; Clarke, R.; Qizilbash, N.; Peto, R.; Collins, R.; Prospective Studies Collaboration. Age-specific relevance of usual blood pressure to vascular mortality: A meta-analysis of individual data for one million adults in 61 prospective studies. Lancet 2002, 360, 1903–1913. [Google Scholar] [CrossRef] [PubMed]
  6. Forouzanfar, M.H.; Liu, P.; Roth, G.A.; Ng, M.; Biryukov, S.; Marczak, L.; Alexander, L.; Estep, K.; Hassen Abate, K.; Akinyemiju, T.F.; et al. Global Burden of Hypertension and Systolic Blood Pressure of at Least 110 to 115 mm Hg, 1990–2015. JAMA 2017, 317, 165–182. [Google Scholar] [CrossRef]
  7. Ehret, G.B. Genome-wide association studies: Contribution of genomics to understanding blood pressure and essential hypertension. Curr. Hypertens. Rep. 2010, 12, 17–25. [Google Scholar] [CrossRef] [PubMed]
  8. Ehret, G.B.; Caulfield, M.J. Genes for blood pressure: An opportunity to understand hypertension. Eur. Heart J. 2013, 34, 951–961. [Google Scholar] [CrossRef]
  9. Evangelou, E.; Warren, H.R.; Mosen-Ansorena, D.; Mifsud, B.; Pazoki, R.; Gao, H.; Ntritsos, G.; Dimou, N.; Cabrera, C.P.; Karaman, I.; et al. Genetic analysis of over 1 million people identifies 535 new loci associated with blood pressure traits. Nat. Genet. 2018, 50, 1412–1425. [Google Scholar] [CrossRef]
  10. Surendran, P.; Feofanova, E.V.; Lahrouchi, N.; Ntalla, I.; Karthikeyan, S.; Cook, J.; Chen, L.; Mifsud, B.; Yao, C.; Kraja, A.T.; et al. Discovery of rare variants associated with blood pressure regulation through meta-analysis of 1.3 million individuals. Nat. Genet. 2020, 52, 1314–1332. [Google Scholar] [CrossRef]
  11. Lip, S.; Padmanabhan, S. Genomics of Blood Pressure and Hypertension: Extending the Mosaic Theory Toward Stratification. Can. J. Cardiol. 2020, 36, 694–705. [Google Scholar] [CrossRef] [PubMed]
  12. Chen, J.; Wang, W.; Li, Z.; Xu, C.; Tian, X.; Zhang, D. Heritability and genome-wide association study of blood pressure in Chinese adult twins. Mol. Genet. Genom. Med. 2021, 9, e1828. [Google Scholar] [CrossRef] [PubMed]
  13. Padmanabhan, S.; Dominiczak, A.F. Genomics of hypertension: The road to precision medicine. Nat. Rev. Cardiol. 2021, 18, 235–250. [Google Scholar] [CrossRef] [PubMed]
  14. Bushueva, O.; Solodilova, M.; Churnosov, M.; Ivanov, V.; Polonikov, A. The Flavin-containing monooxygenase 3 gene and essential hypertension: The joint effect of polymorphism E158K and cigarette smoking on disease susceptibility. Int. J. Hypertens. 2014, 2014, 712169. [Google Scholar] [CrossRef]
  15. Polonikov, A.V.; Ushachev, D.V.; Ivanov, V.P.; Churnosov, M.I.; Freidin, M.B.; Ataman, A.V.; Harbuzova, V.Y.; Bykanova, M.A.; Bushueva, O.Y.; Solodilova, M.A. Altered erythrocyte membrane protein composition mirrors pleiotropic effects of hypertension susceptibility genes and disease pathogenesis. J. Hypertens. 2015, 33, 2265–2277. [Google Scholar] [CrossRef]
  16. Polonikov, A.V.; Bushueva, O.Y.; Bulgakova, I.V.; Freidin, M.B.; Churnosov, M.I.; Solodilova, M.A.; Shvetsov, Y.D.; Ivanov, V.P. A comprehensive contribution of genes for aryl hydrocarbon receptor signaling pathway to hypertension susceptibility. Pharmacogenet. Genom. 2017, 27, 57–69. [Google Scholar] [CrossRef]
  17. Fan, W.; Qu, X.; Li, J.; Wang, X.; Bai, Y.; Cao, Q.; Ma, L.; Zhou, X.; Zhu, W.; Liu, W.; et al. Associations between polymorphisms of the ADIPOQ gene and hypertension risk: A systematic and meta-analysis. Sci. Rep. 2017, 9, 41683. [Google Scholar] [CrossRef]
  18. Polonikov, A.; Bykanova, M.; Ponomarenko, I.; Sirotina, S.; Bocharova, A.; Vagaytseva, K.; Stepanov, V.; Churnosov, M.; Bushueva, O.; Solodilova, M.; et al. The contribution of CYP2C gene subfamily involved in epoxygenase pathway of arachidonic acids metabolism to hypertension susceptibility in Russian population. Clin. Exp. Hypertens. 2017, 39, 306–311. [Google Scholar] [CrossRef]
  19. Moskalenko, M.I.; Milanova, S.N.; Ponomarenko, I.V.; Polonikov, A.V.; Churnosov, M.I. Study of associations of polymorphism of matrix metalloproteinases genes with the development of arterial hypertension in men. Kardiologiia 2019, 59, 31–39. [Google Scholar] [CrossRef]
  20. Zhang, H.; Pushkarev, B.; Zhou, J.; Mu, Y.; Bolshakova, O.; Shrestha, S.; Wang, N.; Jian, B.; Jin, M.; Zhang, K.; et al. CACNA1C rs1006737 SNP increases the risk of essential hypertension in both Chinese Han and ethnic Russian people of Northeast Asia. Medicine 2021, 100, e24825. [Google Scholar] [CrossRef]
  21. Moskalenko, M.; Ponomarenko, I.; Reshetnikov, E.; Dvornyk, V.; Churnosov, M. Polymorphisms of the matrix metalloproteinase genes are associated with essential hypertension in a Caucasian population of Central Russia. Sci. Rep. 2021, 11, 5224. [Google Scholar] [CrossRef] [PubMed]
  22. Takahashi, Y.; Yamazaki, K.; Kamatani, Y.; Kubo, M.; Matsuda, K.; Asai, S. A genome-wide association study identifies a novel candidate locus at the DLGAP1 gene with susceptibility to resistant hypertension in the Japanese population. Sci. Rep. 2021, 11, 19497. [Google Scholar] [CrossRef] [PubMed]
  23. Alsamman, A.M.; Almabrazi, H.; Zayed, H. Whole-Genome Sequencing of 100 Genomes Identifies a Distinctive Genetic Susceptibility Profile of Qatari Patients with Hypertension. J. Pers. Med. 2022, 12, 722. [Google Scholar] [CrossRef] [PubMed]
  24. Jamshidi, J.; Asnaashari, A.; Alipoor, R.; Mohammadi, S.; Roostaei, S.; Samadian, M.M.; Honarmand Aliabadi, S.; Bahramali, E. ATP2B1 rs2681472 and STK39 rs35929607 polymorphisms and risk of Hypertension in Iranian Population. Med. J. Islam. Repub. Iran 2018, 32, 14. [Google Scholar] [CrossRef] [PubMed]
  25. Lule, S.A.; Mentzer, A.J.; Namara, B.; Muwenzi, A.G.; Nassanga, B.; Kizito, D.; Akurut, H.; Lubyayi, L.; Tumusiime, J.; Zziwa, C.; et al. A genome-wide association and replication study of blood pressure in Ugandan early adolescents. Mol. Genet. Genom. Med. 2019, 7, e00950. [Google Scholar] [CrossRef]
  26. Soltész, B.; Pikó, P.; Sándor, J.; Kósa, Z.; Ádány, R.; Fiatal, S. The genetic risk for hypertension is lower among the Hungarian Roma population compared to the general population. PLoS ONE 2020, 15, e0234547. [Google Scholar] [CrossRef]
  27. Kaur, H.; Crawford, D.C.; Liang, J.; Benchek, P.; Zhu, X.; Kallianpur, A.R.; Bush, W.S.; COGENT BP Consortium. Replication of European hypertension associations in a case-control study of 9534 African Americans. PLoS ONE 2021, 16, e0259962. [Google Scholar] [CrossRef]
  28. Althwab, S.A.; Ahmed, A.A.; Rasheed, Z.; Alkhowailed, M.; Hershan, A.; Alsagaby, S.; Alblihed, M.A.; Alaqeel, A.; Alrehaili, J.; Alhumaydhi, F.A.; et al. ATP2B1 genotypes rs2070759 and rs2681472 polymorphisms and risk of hypertension in Saudi population. Nucleosides Nucleotides Nucleic Acids 2021, 40, 1075–1089. [Google Scholar] [CrossRef]
  29. International Consortium for Blood Pressure Genome-Wide Association Studies; Ehret, G.B.; Munroe, P.B.; Rice, K.M.; Bochud, M.; Johnson, A.D.; Chasman, D.I.; Smith, A.V.; Tobin, M.D.; Verwoert, G.C.; et al. Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk. Nature 2011, 478, 103–109. [Google Scholar] [CrossRef]
  30. Surendran, P.; Drenos, F.; Young, R.; Warren, H.; Cook, J.P.; Manning, A.K.; Grarup, N.; Sim, X.; Barnes, D.R.; Witkowska, K.; et al. Trans-ancestry meta-analyses identify rare and common variants associated with blood pressure and hypertension. Nat. Genet. 2016, 48, 1151–1161. [Google Scholar] [CrossRef]
  31. Liu, C.; Kraja, A.T.; Smith, J.A.; Brody, J.A.; Franceschini, N.; Bis, J.C.; Rice, K.; Morrison, A.C.; Lu, Y.; Weiss, S.; et al. Meta-analysis identifies common and rare variants influencing blood pressure and overlapping with metabolic trait loci. Nat. Genet. 2016, 48, 1162–1170. [Google Scholar] [CrossRef] [PubMed]
  32. Ehret, G.B.; Ferreira, T.; Chasman, D.I.; Jackson, A.U.; Schmidt, E.M.; Johnson, T.; Thorleifsson, G.; Luan, J.; Donnelly, L.A.; Kanoni, S.; et al. The genetics of blood pressure regulation and its target organs from association studies in 342,415 individuals. Nat. Genet. 2016, 48, 1171–1184. [Google Scholar] [CrossRef] [PubMed]
  33. Hoffmann, T.J.; Ehret, G.B.; Nandakumar, P.; Ranatunga, D.; Schaefer, C.; Kwok, P.Y.; Iribarren, C.; Chakravarti, A.; Risch, N. Genome-wide association analyses using electronic health records identify new loci influencing blood pressure variation. Nat. Genet. 2017, 49, 54–64. [Google Scholar] [CrossRef] [PubMed]
  34. Wain, L.V.; Vaez, A.; Jansen, R.; Joehanes, R.; van der Most, P.J.; Erzurumluoglu, A.M.; O’Reilly, P.F.; Cabrera, C.P.; Warren, H.R.; Rose, L.M.; et al. Novel Blood Pressure Locus and Gene Discovery Using Genome-Wide Association Study and Expression Data Sets from Blood and the Kidney. Hypertension 2017, 70, e4–e19. [Google Scholar] [CrossRef] [PubMed]
  35. Sung, Y.J.; Winkler, T.W.; de Las Fuentes, L.; Bentley, A.R.; Brown, M.R.; Kraja, A.T.; Schwander, K.; Ntalla, I.; Guo, X.; Franceschini, N.; et al. A Large-Scale Multi-ancestry Genome-wide Study Accounting for Smoking Behavior Identifies Multiple Significant Loci for Blood Pressure. Am. J. Hum. Genet. 2018, 102, 375–400. [Google Scholar] [CrossRef]
  36. Gill, D.; Benyamin, B.; Moore, L.S.P.; Monori, G.; Zhou, A.; Koskeridis, F.; Evangelou, E.; Laffan, M.; Walker, A.P.; Tsilidis, K.K.; et al. Associations of genetically determined iron status across the phenome: A mendelian randomization study. PLoS Med. 2019, 16, e1002833. [Google Scholar] [CrossRef]
  37. Lu, W.; Pikhart, H.; Tamosiunas, A.; Kubinova, R.; Capkova, N.; Malyutina, S.; Pająk, A.; Bobak, M. Prevalence, awareness, treatment and control of hypertension, diabetes and hypercholesterolemia, and associated risk factors in the Czech Republic, Russia, Poland and Lithuania: A cross-sectional study. BMC Public Health 2022, 22, 883. [Google Scholar] [CrossRef]
  38. Talbert, P.B.; Henikoff, S. Histone variants at a glance. J. Cell Sci. 2021, 134, jcs244749. [Google Scholar] [CrossRef]
  39. Li, Y.; Wang, Y.; Zhang, C.; Yuan, W.; Wang, J.; Zhu, C.; Chen, L.; Huang, W.; Zeng, W.; Wu, X.; et al. ZNF322, a novel human C2H2 Kruppel-like zinc-finger protein, regulates transcriptional activation in MAPK signaling pathways. Biochem. Biophys. Res. Commun. 2004, 325, 1383–1392. [Google Scholar] [CrossRef]
  40. GeneCards: The Human Gene Database. Available online: https://www.genecards.org/cgi-bin/carddisp.pl?gene=STN1 (accessed on 16 October 2022).
  41. Wan, M.; Qin, J.; Songyang, Z.; Liu, D. OB fold-containing protein 1 (OBFC1), a human homolog of yeast Stn1, associates with TPP1 and is implicated in telomere length regulation. J. Biol. Chem. 2009, 284, 26725–26731. [Google Scholar] [CrossRef]
  42. Said, M.A.; Eppinga, R.N.; Hagemeijer, Y.; Verweij, N.; van der Harst, P. Telomere Length and Risk of Cardiovascular Disease and Cancer. J. Am. Coll. Cardiol. 2017, 70, 506–507. [Google Scholar] [CrossRef] [PubMed]
  43. Maubaret, C.G.; Salpea, K.D.; Romanoski, C.E.; Folkersen, L.; Cooper, J.A.; Stephanou, C.; Li, K.W.; Palmen, J.; Hamsten, A.; Neil, A.; et al. Association of TERC and OBFC1 haplotypes with mean leukocyte telomere length and risk for coronary heart disease. PLoS ONE 2013, 8, e83122. [Google Scholar] [CrossRef] [PubMed]
  44. Rodriguez-Iturbe, B.; Pons, H.; Johnson, R.J. Role of the Immune System in Hypertension. Physiol. Rev. 2017, 97, 1127–1164. [Google Scholar] [CrossRef] [PubMed]
  45. Drummond, G.R.; Vinh, A.; Guzik, T.J.; Sobey, C.G. Immune mechanisms of hypertension. Nat. Rev. Immunol. 2019, 19, 517–532. [Google Scholar] [CrossRef]
  46. Carnagarin, R.; Matthews, V.; Zaldivia, M.T.K.; Peter, K.; Schlaich, M.P. The bidirectional interaction between the sympathetic nervous system and immune mechanisms in the pathogenesis of hypertension. Br. J. Pharmacol. 2019, 176, 1839–1852. [Google Scholar] [CrossRef]
  47. Mikolajczyk, T.P.; Guzik, T.J. Adaptive Immunity in Hypertension. Curr. Hypertens. Rep. 2019, 21, 68. [Google Scholar] [CrossRef]
  48. Churnosov, M.; Abramova, M.; Reshetnikov, E.; Lyashenko, I.V.; Efremova, O.; Churnosova, M.; Ponomarenko, I. Polymorphisms of hypertension susceptibility genes as a risk factors of preeclampsia in the Caucasian population of central Russia. Placenta 2022, 129, 51–61. [Google Scholar] [CrossRef]
  49. Litovkina, O.; Nekipelova, E.; Dvornyk, V.; Polonikov, A.; Efremova, O.; Zhernakova, N.; Reshetnikov, E.; Churnosov, M. Genes involved in the regulation of vascular homeostasis determine renal survival rate in patients with chronic glomerulo-nephritis. Gene 2014, 546, 112–116. [Google Scholar] [CrossRef]
  50. Reshetnikov, E.A.; Akulova, L.Y.; Dobrodomova, I.S.; Dvornyk, V.Y.; Polonikov, A.V.; Churnosov, M.I. The insertion-deletion polymorphism of the ACE gene is associated with increased blood pressure in women at the end of pregnancy. J. Renin Angiotensin Aldosterone Syst. 2015, 16, 623–632. [Google Scholar] [CrossRef]
  51. Ivanova, T.A. Polymorphic loci of AC026703.1 and HFE genes are associated with severe hypertension. Res. Results Biomed. 2023, 9, 22–38. [Google Scholar] [CrossRef]
  52. World Health Organization. Diet, Nutrition, and the Prevention of Chronic Diseases: Report of a Joint WHO/FAO Expert 802 Consultation; World Health Organization: Geneva, Switzerland, 2003; 148p. [Google Scholar]
  53. Ponomarenko, I.; Reshetnikov, E.; Altuchova, O.; Polonikov, A.; Sorokina, I.; Yermachenko, A.; Dvornyk, V.; Golovchenko, O.; Churnosov, M. Association of genetic polymorphisms with age at menarche in Russian women. Gene 2019, 686, 228–236. [Google Scholar] [CrossRef] [PubMed]
  54. Eliseeva, N.; Ponomarenko, I.; Reshetnikov, E.; Dvornyk, V.; Churnosov, M. LOXL1 gene polymorphism candidates for exfoliation glaucoma are also associated with a risk for primary open-angle glaucoma in a Caucasian population from central Russia. Mol. Vis. 2021, 27, 262–269. [Google Scholar] [PubMed]
  55. Ward, L.D.; Kellis, M. HaploReg v4: Systematic mining of putative causal variants, cell types, regulators and target genes for human complex traits and disease. Nucleic Acids Res. 2016, 44, D877–D881. [Google Scholar] [CrossRef]
  56. Abramova, M.; Churnosova, M.; Efremova, O.; Aristova, I.; Reshetnikov, E.; Polonikov, A.; Churnosov, M.; Ponomarenko, I. Effects of pre-pregnancy over-weight/obesity on the pattern of association of hypertension susceptibility genes with preeclampsia. Life 2022, 12, 2018. [Google Scholar] [CrossRef] [PubMed]
  57. Abramova, M.Y. Genetic markers of severe preeclampsia. Res. Results Biomed. 2022, 8, 305–316. [Google Scholar] [CrossRef]
  58. Reshetnikov, E.; Zarudskaya, O.; Polonikov, A.; Bushueva, O.; Orlova, V.; Krikun, E.; Dvornyk, V.; Churnosov, M. Genetic markers for inherited thrombophilia are associated with fetal growth retardation in the population of Central Russia. J. Obstet. Gynaecol. Res. 2017, 43, 1139–1144. [Google Scholar] [CrossRef]
  59. Golovchenko, O.; Abramova, M.; Ponomarenko, I.; Reshetnikov, E.; Aristova, I.; Polonikov, A.; Dvornyk, V.; Churnosov, M. Functionally significant polymorphisms of ESR1and PGR and risk of intrauterine growth restriction in population of Central Russia. Eur. J. Obstet. Gynecol. Reprod. Biol. 2020, 253, 52–57. [Google Scholar] [CrossRef]
  60. Starikova, D.; Ponomarenko, I.; Reshetnikov, E.; Dvornyk, V.; Churnosov, M. Novel Data about Association of the Functionally Significant Polymorphisms of the MMP9 Gene with Exfoliation Glaucoma in the Caucasian Population of Central Russia. Ophthalmic Res. 2021, 64, 458–464. [Google Scholar] [CrossRef]
  61. Sirotina, S.; Ponomarenko, I.; Kharchenko, A.; Bykanova, M.; Bocharova, A.; Vagaytseva, K.; Stepanov, V.; Churnosov, M.; Solodilova, M.; Polonikov, A. A Novel Polymorphism in the Promoter of the CYP4A11 Gene Is Associated with Susceptibility to Coronary Artery Disease. Dis. Markers 2018, 2018, 5812802. [Google Scholar] [CrossRef]
  62. Tikunova, E.; Ovtcharova, V.; Reshetnikov, E.; Dvornyk, V.; Polonikov, A.; Bushueva, O.; Churnosov, M. Genes of tumor necrosis factors and their receptors and the primary open angle glaucoma in the population of Central Russia. Int. J. Ophthalmol. 2017, 10, 1490–1494. [Google Scholar] [CrossRef]
  63. Pavlova, N.; Demin, S.; Churnosov, M.; Reshetnikov, E.; Aristova, I.; Churnosova, M.; Ponomarenko, I. The Modifying Effect of Obesity on the Association of Matrix Metalloproteinase Gene Polymorphisms with Breast Cancer Risk. Biomedicines 2022, 10, 2617. [Google Scholar] [CrossRef] [PubMed]
  64. Ponomarenko, I.; Reshetnikov, E.; Polonikov, A.; Verzilina, I.; Sorokina, I.; Elgaeva, E.E.; Tsepilov, Y.A.; Yermachenko, A.; Dvornyk, V.; Churnosov, M. Candidate genes for age at menarche are associated with endometriosis. Reprod. Biomed. Online 2020, 41, 943–956. [Google Scholar] [CrossRef] [PubMed]
  65. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef]
  66. Calle, M.L.; Urrea, V.; Malats, N.; Van Steen, K. Mbmdr: An R package for exploring gene-gene interactions associated with binary or quantitative traits. Bioinformatics 2010, 26, 2198–2199. [Google Scholar] [CrossRef] [PubMed]
  67. Ponomarenko, I.V. Using the method of Multifactor Dimensionality Reduction (MDR) and its modifications for analysis of gene-gene and gene-environment interactions in genetic-epidemiological studies (review). Res. Results Biomed. 2019, 5, 4–21. [Google Scholar] [CrossRef]
  68. Sourceforge. Available online: http://sourceforge.net/projects/mdr (accessed on 22 November 2022).
  69. Pavlova, N.; Demin, S.; Churnosov, M.; Reshetnikov, E.; Aristova, I.; Churnosova, M.; Ponomarenko, I. Matrix Metalloproteinase Gene Polymorphisms Are Associated with Breast Cancer in the Caucasian Women of Russia. Int. J. Mol. Sci. 2022, 23, 12638. [Google Scholar] [CrossRef]
  70. Reshetnikov, E.; Ponomarenko, I.; Golovchenko, O.; Sorokina, I.; Batlutskaya, I.; Yakunchenko, T.; Dvornyk, V.; Polonikov, A.; Churnosov, M. The VNTR polymorphism of the endothelial nitric oxide synthase gene and blood pressure in women at the end of pregnancy. Taiwan J. Obstet. Gynecol. 2019, 58, 390–395. [Google Scholar] [CrossRef]
  71. Minyaylo, O.; Ponomarenko, I.; Reshetnikov, E.; Dvornyk, V.; Churnosov, M. Functionally significant polymorphisms of the MMP-9 gene are associated with peptic ulcer disease in the Caucasian population of Central Russia. Sci. Rep. 2021, 11, 13515. [Google Scholar] [CrossRef]
  72. Che, R.; Jack, J.R.; Motsinger-Reif, A.A.; Brown, C.C. An adaptive permutation approach for genome-wide association study: Evaluation and recommendations for use. BioData Min. 2014, 7, 9. [Google Scholar] [CrossRef]
  73. Adzhubei, I.; Jordan, D.M.; Sunyaev, S.R. Predicting functional effect of human missense mutations using PolyPhen-2. Curr. Protoc. Hum. Genet. 2013, 76, 7–20. [Google Scholar] [CrossRef]
  74. Kumar, P.; Henikoff, S.; Ng, P.C. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat. Protoc. 2009, 7, 1073–1081. [Google Scholar] [CrossRef] [PubMed]
  75. Westra, H.J.; Peters, M.J.; Esko, T.; Yaghootkar, H.; Schurmann, C.; Kettunen, J.; Christiansen, M.W.; Fairfax, B.P.; Schramm, K.; Powell, J.E.; et al. Systematic identification of trans eQTLs as putative drivers of known disease associations. Nat. Genet. 2013, 45, 1238–1243. [Google Scholar] [CrossRef] [PubMed]
  76. GTEx Consortium. The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science 2020, 369, 1318–1330. [Google Scholar] [CrossRef] [PubMed]
  77. Gene Ontology Consortium. The Gene Ontology resource: Enriching a GOld mine. Nucleic Acids Res. 2021, 49, D325–D334. [Google Scholar] [CrossRef] [PubMed]
  78. Franz, M.; Rodriguez, H.; Lopes, C.; Zuberi, K.; Montojo, J.; Bader, G.D.; Morris, Q. GeneMANIA update 2018. Nucleic Acids Res. 2018, 46, W60–W64. [Google Scholar] [CrossRef] [PubMed]
  79. Szklarczyk, D.; Kirsch, R.; Koutrouli, M.; Nastou, K.; Mehryary, F.; Hachilif, R.; Gable, A.L.; Fang, T.; Doncheva, N.T.; Pyysalo, S.; et al. The STRING database in 2023: Protein-protein association networks and functional enrichment analyses for any sequenced genome of interest. Nucleic Acids Res. 2023, 51, D638–D646. [Google Scholar] [CrossRef]
  80. Ponomarenko, I.; Reshetnikov, E.; Polonikov, A.; Verzilina, I.; Sorokina, I.; Yermachenko, A.; Dvornyk, V.; Churnosov, M. Candidate genes for age at menarche are associated with uterine leiomyoma. Front. Genet. 2021, 11, 512940. [Google Scholar] [CrossRef]
  81. Golovchenko, I.; Aizikovich, B.; Golovchenko, O.; Reshetnikov, E.; Churnosova, M.; Aristova, I.; Ponomarenko, I.; Churnosov, M. Sex Hormone Candidate Gene Polymorphisms Are Associated with Endometriosis. Int. J. Mol. Sci. 2022, 23, 13691. [Google Scholar] [CrossRef]
  82. Ponomarenko, I.; Reshetnikov, E.; Polonikov, A.; Sorokina, I.; Yermachenko, A.; Dvornyk, V.; Churnosov, M. Candidate genes for age at menarche are associated with endometrial hyperplasia. Gene 2020, 757, 4933. [Google Scholar] [CrossRef]
  83. Novakov, V.; Novakova, O.; Churnosova, M.; Sorokina, I.; Aristova, I.; Polonikov, A.; Reshetnikov, E.; Churnosov, M. Intergenic Interactions of SBNO1, NFAT5 and GLT8D1 Determine the Susceptibility to Knee Osteoarthritis among Europeans of Russia. Life 2023, 13, 405. [Google Scholar] [CrossRef]
  84. Feitosa, M.F.; Kraja, A.T.; Chasman, D.I.; Sung, Y.J.; Winkler, T.W.; Ntalla, I.; Guo, X.; Franceschini, N.; Heng, C.K.; Sim, X.; et al. Novel genetic associations for blood pressure identified via gene-alcohol interaction in up to 570K individuals across multipleancestries. PLoS ONE 2018, 13, e0198166. [Google Scholar] [CrossRef]
  85. German, C.A.; Sinsheimer, J.S.; Klimentidis, Y.C.; Zhou, H.; Zhou, J.J. Ordered multinomial regression for genetic association analysis of ordinal phenotypes at Biobank scale. Genet. Epidemiol. 2019, 44, 248–260. [Google Scholar] [CrossRef]
  86. Giri, A.; Hellwege, J.N.; Keaton, J.M.; Park, J.; Qiu, C.X.; Warren, H.R.; Torstenson, E.S.; Kovesdy, C.P.; Sun, Y.V.; Wilson, O.D.; et al. Trans-ethnic association study of blood pressure determinants in over 750,000 individuals. Nat. Genet. 2019, 51, 51–62. [Google Scholar] [CrossRef]
  87. Hartiala, J.A.; Han, Y.; Jia, Q.; Hilser, J.R.; Huang, P.; Gukasyan, J.; Schwartzman, W.S.; Cai, Z.; Biswas, S.; Trégouët, D.A.; et al. Genome-wide analysis identifies novel susceptibility loci for myocardial infarction. Eur. Heart J. 2021, 42, 919–933. [Google Scholar] [CrossRef] [PubMed]
  88. Jeong, H.; Jin, H.-S.; Kim, S.-S.; Shin, D. Identifying Interactions between Dietary Sodium, Potassium, Sodium–Potassium Ratios, and FGF5 rs16998073 Variants and Their Associated Risk for Hypertension in Korean Adults. Nutrients 2020, 12, 2121. [Google Scholar] [CrossRef] [PubMed]
  89. Kato, N.; Loh, M.; Takeuchi, F.; Verweij, N.; Wang, X.; Zhang, W.; Kelly, T.N.; Saleheen, D.; Lehne, B.; Leach, I.M.; et al. Trans-ancestry genome-wide association study identifies 12 genetic loci influencing blood pressure and implicates a role for DNA methylation. Nat. Genet. 2015, 47, 1282–1293. [Google Scholar] [CrossRef] [PubMed]
  90. Koyama, S.; Ito, K.; Terao, C.; Akiyama, M.; Horikoshi, M.; Momozawa, Y.; Matsunaga, H.; Ieki, H.; Ozaki, K.; Onouchi, Y.; et al. Population-specific and trans-ancestry genome-wide analyses identify distinct and shared genetic risk loci for coronary artery disease. Nat. Genet. 2020, 52, 1169–1177. [Google Scholar] [CrossRef]
  91. Levy, D.; Ehret, G.B.; Rice, K.; Verwoert, G.C.; Launer, L.J.; Dehghan, A.; Glazer, N.L.; Morrison, A.C.; Johnson, A.D.; Aspelund, T.; et al. Genome-wide association study of blood pressure and hypertension. Nat. Genet. 2009, 41, 677–687. [Google Scholar] [CrossRef]
  92. Nelson, C.P.; Goel, A.; Butterworth, A.S.; Kanoni, S.; Webb, T.R.; Marouli, E.; Zeng, L.; Ntalla, I.; Lai, F.Y.; Hopewell, J.C.; et al. Association analyses based on false discovery rate implicate new loci for coronary artery disease. Nat. Genet. 2017, 49, 1385–1391. [Google Scholar] [CrossRef]
  93. Nikpay, M.; Goel, A.; Won, H.H.; Hall, L.M.; Willenborg, C.; Kanoni, S.; Saleheen, D.; Kyriakou, T.; Nelson, C.P.; Hopewell, J.C.; et al. A comprehensive 1000 Genomes-based genome-wide association meta-analysis of coronary artery disease. Nat. Genet. 2015, 47, 1121–1130. [Google Scholar] [CrossRef]
  94. Sakaue, S.; Kanai, M.; Tanigawa, Y.; Karjalainen, J.; Kurki, M.; Koshiba, S.; Narita, A.; Konuma, T.; Yamamoto, K.; Akiyama, M.; et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nat. Genet. 2021, 53, 1415–1424. [Google Scholar] [CrossRef] [PubMed]
  95. Shungin, D.; Winkler, T.W.; Croteau-Chonka, D.C.; Ferreira, T.; Locke, A.E.; Mägi, R.; Strawbridge, R.J.; Pers, T.H.; Fischer, K.; Justice, A.E.; et al. New genetic loci link adipose and insulin biology to body fat distribution. Nature 2015, 518, 187–196. [Google Scholar] [CrossRef] [PubMed]
  96. Tachmazidou, I.; Süveges, D.; Min, J.L.; Ritchie, G.R.S.; Steinberg, J.; Walter, K.; Iotchkova, V.; Schwartzentruber, J.; Huang, J.; Memari, Y.; et al. Whole-Genome Sequencing Coupled to Imputation Discovers Genetic Signals for Anthropometric Traits. Am. J. Hum. Genet. 2017, 100, 865–884. [Google Scholar] [CrossRef] [PubMed]
  97. Takeuchi, F.; Akiyama, M.; Matoba, N.; Katsuya, T.; Nakatochi, M.; Tabara, Y.; Narita, A.; Saw, W.Y.; Moon, S.; Spracklen, C.N.; et al. Interethnic analyses of blood pressure loci in populations of East Asian and European descent. Nat. Commun. 2018, 9, 5052. [Google Scholar] [CrossRef]
  98. Wain, L.V.; Verwoert, G.C.; O’Reilly, P.F.; Shi, G.; Johnson, T.; Johnson, A.D.; Bochud, M.; Rice, K.M.; Henneman, P.; Smith, A.V.; et al. Genome-wide association study identifies six new loci influencing pulse pressure and mean arterial pressure. Nat. Genet. 2011, 43, 1005–1011. [Google Scholar] [CrossRef]
  99. Zhou, W.; Nielsen, J.B.; Fritsche, L.G.; Dey, R.; Gabrielsen, M.E.; Wolford, B.N.; LeFaive, J.; VandeHaar, P.; Gagliano, S.A.; Gifford, A.; et al. Efficiently controlling for case-control imbalance and sample relatedness in large-scale genetic association studies. Nat. Genet. 2018, 50, 1335–1341. [Google Scholar] [CrossRef]
  100. Zhu, Z.; Wang, X.; Li, X.; Lin, Y.; Shen, S.; Liu, C.L.; Hobbs, B.D.; Hasegawa, K.; Liang, L.; International COPD Genetics Consortium; et al. Genetic overlap of chronic obstructive pulmonary disease and cardiovascular disease-related traits: A large-scale genome-wide cross-trait analysis. Respir. Res. 2019, 20, 64. [Google Scholar] [CrossRef]
Ijms 24 08309 g001 550
Figure 1. The entropy graph of the SNP × SNP interactions with AH based on the MDR analysis. Positive values of entropy indicate synergistic interactions, while the negative values indicate redundancy. The red and orange colors denote strong and moderate synergism, respectively, brown color denotes the independent effect, and green and blue colors denote moderate and strong antagonism.
Figure 1. The entropy graph of the SNP × SNP interactions with AH based on the MDR analysis. Positive values of entropy indicate synergistic interactions, while the negative values indicate redundancy. The red and orange colors denote strong and moderate synergism, respectively, brown color denotes the independent effect, and green and blue colors denote moderate and strong antagonism.
Ijms 24 08309 g001
Ijms 24 08309 g002 550
Figure 2. The interaction networks of the candidate genes for the AH in various tissues/organs inferred using GeneMANIA (http://genemania.org (accessed on 16 October 2022)). The candidate genes are cross-shaded.
Figure 2. The interaction networks of the candidate genes for the AH in various tissues/organs inferred using GeneMANIA (http://genemania.org (accessed on 16 October 2022)). The candidate genes are cross-shaded.
Ijms 24 08309 g002
Ijms 24 08309 g003 550
Figure 3. AH-involved protein–protein interaction networks inferred using STRING resource.
Figure 3. AH-involved protein–protein interaction networks inferred using STRING resource.
Ijms 24 08309 g003
Ijms 24 08309 g004 550
Figure 4. AH-involved protein–protein interaction clusters inferred using STRING resource (three groups of PPI clusters are highlighted in color: cluster 1, red (A); cluster 2, green (B); cluster 3, blue (C); summary of three clusters (D)).
Figure 4. AH-involved protein–protein interaction clusters inferred using STRING resource (three groups of PPI clusters are highlighted in color: cluster 1, red (A); cluster 2, green (B); cluster 3, blue (C); summary of three clusters (D)).
Ijms 24 08309 g004
Table
Table 1. Associations of the studied gene polymorphisms with AH.
Table 1. Associations of the studied gene polymorphisms with AH.
Gene (SNP, Major/Minor Alleles)nAllelic Model Additive Model Dominant ModelRecessive Model
OR95% CIpOR95% CIpOR95% CIpOR95% CIp
L95U95L95U95L95U95L95U95
Model 1
AC026703.1 (rs1173771,G/A)13170.900.771.060.2160.890.691.140.3540.750.511.100.1401.030.651.620.905
HFE (rs1799945,C/G)13730.940.771.150.5501.040.781.400.7810.900.631.280.5592.531.036.230.043
BAG6 (rs805303,G/A)13490.970.821.140.6830.870.681.110.2640.950.671.340.7500.650.401.040.075
PLCE1 (rs932764,A/G)13190.870.741.020.0940.810.631.040.0960.650.431.020.0560.880.581.320.525
OBFC1 (rs4387287,C/A)12600.900.671.220.5500.900.651.250.5420.900.611.330.6050.790.311.990.611
ARHGAP42 (rs633185,C/G)13771.020.851.220.8131.090.841.420.5251.080.761.520.6731.250.682.320.472
CERS5 (rs7302981,G/A)13021.030.871.220.7110.940.731.210.6151.030.711.480.8820.760.471.220.249
ATP2B1 (rs2681472,A/G)13291.040.821.310.7621.170.821.670.3841.170.791.740.4371.510.415.500.532
TBX2 (rs8068318,T/C)12921.100.921.330.2971.140.861.520.3561.170.821.660.3981.250.612.550.541
RGL3 (rs167479,T/G)13330.930.791.090.3670.820.641.050.1100.860.571.290.4600.690.471.020.061
Model 2
AC026703.1 (rs1173771,G/A) 0.890.691.150.3740.760.511.110.1531.030.651.640.893
HFE (rs1799945,C/G) 1.010.751.360.9460.870.611.250.4652.481.026.070.045
BAG6 (rs805303,G/A) 0.890.691.140.3430.980.691.390.9060.650.401.060.085
PLCE1 (rs932764,A/G) 0.810.631.040.0930.640.421.020.0540.880.581.330.535
OBFC1 (rs4387287,C/A) 0.880.631.220.4480.870.591.290.4900.780.302.010.608
ARHGAP42 (rs633185,C/G) 1.110.851.440.4551.010.781.560.5921.280.692.360.441
CERS5 (rs7302981,G/A) 0.930.721.210.6001.010.701.460.9480.760.471.240.275
ATP2B1 (rs2681472,A/G) 1.150.801.660.4511.170.781.760.4481.200.324.520.792
TBX2 (rs8068318,T/C) 1.150.861.540.3481.170.841.720.3241.140.552.360.716
RGL3 (rs167479,T/G) 0.800.631.030.0780.810.541.220.3200.680.461.010.057
Note: For Model 2, calculations of the allelic model were not performed because their results are identical to those of Model 1 (covariates are not used when calculating the allelic model). All results were obtained after adjustment for covariates. List covariates for Model 1: BMI, TC, TG, LDL-C, HDL-C, blood glucose, smokers. List covariates for Model 2: BMI, TC, TG, LDL-C, HDL-C, blood glucose, smokers, low physical activity, high fatty foods consumption. OR—odds ratio; 95% CI—95% confidence interval; pperm values ≤ 0.05 are shown in bold.
Table
Table 2. SNP × SNP interactions significantly associated with AH.
Table 2. SNP × SNP interactions significantly associated with AH.
NSNP × SNP Interaction ModelsNHbetaHWHNLbetaLWLpperm
Two-order interaction models
1rs7302981 CERS5 × rs805303 BAG6 20.3728.023−0.4339.260.037
2rs805303 BAG6 × rs1173771 AC026703.111.3187.992−0.4989.440.047
Three-order interaction models
1rs932764 PLCE1 × rs7302981 CERS5 × rs805303 BAG6 10.7663.742−0.79922.250.001
2rs932764 PLCE1 × rs805303 BAG6 × rs4387287 OBFC120.86114.335−0.66820.700.006
3rs7302981 CERS5 × rs805303 BAG6 × rs1173771 AC026703.1 30.57511.103−1.08520.180.006
Four-order interaction models
1rs7302981 CERS5 × rs805303 BAG6 × rs1173771 AC026703.1 × rs167479 RGL310.6596.6410−1.08947.36<0.001
2rs7302981 CERS5 × rs1799945 HFE × rs805303 BAG6 × rs167479 RGL320.68210.039−1.08543.66<0.001
3rs932764 PLCE1 × rs7302981 CERS5 × rs805303 BAG6 × rs167479 RGL330.89511.666−1.08138.100.001
Note: The results were obtained using the MB-MDR method with adjustment for covariates (Model 1); NH—number of significant high risk genotypes in the interaction; betaH—regression coefficient for high risk exposition in the step2 analysis; WH—Wald statistic for high risk category; NL—number of significant low- risk genotypes in the interaction; betaL—regression coefficient for low risk exposition in the step2 analysis; WL—Wald statistic for low risk category; pperm—permutation p value for the interaction model (1.000 permutations).
Table
Table 3. Phenotypic characteristics of the study participants.
Table 3. Phenotypic characteristics of the study participants.
ParametersAH,
Mean ± SD, % (n)		Controls,
Mean ± SD, % (n)		p
N939466
Gender (Male/
Female)		60.06/39.94
(564/375)		55.15/44.85
(257/209)		0.09
Age (years)58.08 ± 8.9157.82 ± 9.520.77
BMI (kg/m2)30.78 ± 5.0824.94 ± 3.14<0.001
SBP (mmHg)182.48 ± 28.26122.58 ± 11.49<0.001
DBP (mmHg)105.84 ± 13.4777.65 ± 6.93<0.001
TC (mM)5.71 ± 1.295.26 ± 1.04<0.001
HDL-C (mM)1.34 ± 0.421.52 ± 0.42<0.001
LDL-C (mM)3.78 ± 1.113.22 ± 0.74<0.001
TG (mM)1.92 ± 1.031.22 ± 0.71<0.001
BG (mM)5.92 ± 1.684.88 ± 0.95<0.001
Smoking 38.33 (353)19.76 (84)<0.001
Alcohol abuse 5.79 (53)3.12 (13)0.051
Low physical activity58.68 (551)27.47 (128)<0.001
Low fruit/vegetable consumption11.39 (107)8.15 (38)0.074
High fatty foods consumption24.71 (232)10.30 (48)<0.001
High sodium consumption16.72 (157)13.30 (62)0.113
Note: Clinical characteristics of age, BMI, SBP, DBP, HDL-C, LDL-C, TG and TC are given as means ± SD and other values as number of individuals; BMI—body mass index; BG—blood glucose; SBP—systolic blood pressure; DBP—diastolic blood pressure; TC—total cholesterol; HDL-C—high-density lipoprotein cholesterol; TG—triglycerides; LDL-C—low-density lipoprotein cholesterol.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ivanova, T.; Churnosova, M.; Abramova, M.; Ponomarenko, I.; Reshetnikov, E.; Aristova, I.; Sorokina, I.; Churnosov, M. Risk Effects of rs1799945 Polymorphism of the HFE Gene and Intergenic Interactions of GWAS-Significant Loci for Arterial Hypertension in the Caucasian Population of Central Russia. Int. J. Mol. Sci. 2023, 24, 8309. https://doi.org/10.3390/ijms24098309

AMA Style

Ivanova T, Churnosova M, Abramova M, Ponomarenko I, Reshetnikov E, Aristova I, Sorokina I, Churnosov M. Risk Effects of rs1799945 Polymorphism of the HFE Gene and Intergenic Interactions of GWAS-Significant Loci for Arterial Hypertension in the Caucasian Population of Central Russia. International Journal of Molecular Sciences. 2023; 24(9):8309. https://doi.org/10.3390/ijms24098309

Chicago/Turabian Style

Ivanova, Tatiana, Maria Churnosova, Maria Abramova, Irina Ponomarenko, Evgeny Reshetnikov, Inna Aristova, Inna Sorokina, and Mikhail Churnosov. 2023. "Risk Effects of rs1799945 Polymorphism of the HFE Gene and Intergenic Interactions of GWAS-Significant Loci for Arterial Hypertension in the Caucasian Population of Central Russia" International Journal of Molecular Sciences 24, no. 9: 8309. https://doi.org/10.3390/ijms24098309

APA Style

Ivanova, T., Churnosova, M., Abramova, M., Ponomarenko, I., Reshetnikov, E., Aristova, I., Sorokina, I., & Churnosov, M. (2023). Risk Effects of rs1799945 Polymorphism of the HFE Gene and Intergenic Interactions of GWAS-Significant Loci for Arterial Hypertension in the Caucasian Population of Central Russia. International Journal of Molecular Sciences, 24(9), 8309. https://doi.org/10.3390/ijms24098309

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop