Association of rs2070600 in advanced glycosylation end‐product specific receptor with prognosis of heart failure

Abstract Aims Our objective was to investigate the association of common variants in the coding region of advanced glycosylation end‐product specific receptor (RAGE) and the prognosis of heart failure (HF). Methods and results A total of 3394 HF patients were continuously enrolled from January 2009 to August 2018 with a median follow‐up of 20.4 months. Additionally, 2861 healthy subjects also participated in the study. By sequencing these two groups, we identified a common functional missense variant rs2070600 in the coding region of RAGE, which showed a significant association with the prognosis of HF [hazard ratio = 0.53, 95%, confidence interval (CI) = 0.30–0.94, P = 0.03], but no association with the risk of HF (odds ratio = 0.52, 95%, CI = 0.66–1.04, P = 0.106). A series of functional assays revealed that rs2070600‐A, but not ‐G allele, suppressed the expression of RAGE protein by facilitating the binding of miR‐125a‐3p. Furthermore, the RAGE messenger RNA levels of human peripheral blood lymphocytes were reduced in subjects with the rs2070600‐AA genotype compared with subjects with the rs2070600‐GG or ‐AG genotypes. Additionally, our Western blot results from human heart tissue showed increased RAGE expression in HF samples compared with that in healthy donors. Conclusions Our results demonstrate that the common missense variant rs2070600‐A allele is associated with a reduced risk of cardiovascular death and cardiac transplantation by facilitating the binding of miR‐125a‐3p.


Introduction
Heart failure (HF) is the terminal phase of many cardiovascular diseases, which are often driven by molecular events and genetic causes. Patients with HF have poor prognosis and poor therapeutic options. The optimal development of precision treatment relies on risk stratification and prognostic biomarker identification. Several genome-wide association studies have identified hundreds of loci associated with the risk of cardiomyopathy 1 and coronary artery disease. 2 However, the link between specific genes and HF is largely undescribed.
Advanced glycosylation end product receptor (RAGE, GenBank accession no. NC_000006) localizes on chromosome 6p21 and is a key regulator of endogenous inflammatory effects 3 and reactive oxygen species. 4 The advanced glycosylation end product/RAGE axis is associated with multiple diseases and complications in humans and animal models 5,6 and has served as a key therapeutic target 7 for the synthesis of small-molecule inhibitors. 8 The RAGE gene encodes the advanced glycosylation end product receptor, which belongs to the immunoglobulin superfamily of cell surface receptors because of its special domain 9 that can bind multiple ligands 10 (e.g. advanced glycosylation end products, β-sheets, S100, and HMGB1). This binding allows the receptors to act as nicotinamide adenine dinucleotide phosphate and nuclear factor-κ-gene binding 11 that mediate the pathophysiological process of cardiovascular diseases. [12][13][14][15] Polymorphisms and single-nucleotide polymorphisms (SNPs) of RAGE have been implicated in disease development in previous studies. 16 Notably, a common nonsynonymous SNP, rs2070600, located in a coding region that causes a conversion at position 82 from glycine to serine (G82S), may be linked to rheumatoid arthritis, 17 cancer, 18 lung diseases, 19 diabetic nephropathy, 20 Alzheimer's disease, 21 and coronary artery disease. 22 However, the association of rs2070600 with HF is undefined, 23,24 and further research is needed in particular ethnic groups.
In this report, we hypothesized that the common missense variant of RAGE, rs2070600, is associated with HF outcomes and can be used to predict the prognosis of HF. We evaluated the association between RAGE risk variants and the prognosis of HF. Another case-control study was performed to assess the risk of HF. We then provided functional evidence that the common mutation at the chromosome 6p21.32 risk loci affects RAGE gene expression and investigated the underlying mechanisms at the genetic and cellular levels.

Study design and eligibility
This study was approved by the Ethics Committee of Tongji Hospital and written informed consent was obtained from all participants. The investigation conformed with the principles outlined in the Declaration of Helsinki. Details on sample recruitment, inclusion and exclusion criteria, data collection, and risk factor definition are described in our previous report. 25 The clinical characteristics of the study samples are shown in Table 1. The reference population for the resequencing efforts were 48 healthy, Han Chinese subjects recruited from a group of individuals undergoing routine health examinations at Tongji Hospital in Wuhan, Hubei province. The HF cohort, including 3394 HF patients, was enrolled from hospitalized patients in Tongji Hospital between January 2009 and August 2018. Finally, 3022 patients remained from the 3394 for prognosis analysis, as follow-up and genotype was not completed in 200 and 172 patients, respectively. The primary endpoints were heart transplantation and cardiovascular death. The 2861 healthy, community-based Han Chinese individuals were originally recruited from the general population in 2004, with ages ranging from 23 to 82 years. In addition, 232 patients without coronary artery disease were confirmed with coronary angiography and collected from January to August 2018. Details on inclusion criteria and data collection are provided in the online Supporting Information methods.

Human tissue samples
The heart samples used in this study were obtained from five patients who underwent heart transplantation at Tongji Hospital and three healthy donors that suffered traffic accidents. Samples were frozen in liquid nitrogen and then stored at À80°C until use. The clinical characteristics of the study samples are summarized in Supporting Information, Table S1. This study was approved by the Institutional Ethics Committee of Tongji Hospital and is in accordance with the principles of the Declaration of Helsinki. Written informed consent was obtained from the study participants or their relatives.

Single-nucleotide polymorphism selection and genotyping
DNA of human samples was extracted from peripheral leucocytes according to our previously published method. 26 Human peripheral lymphocytes were isolated using lymphocyte separation medium (Tian Jin Hao Yang Biological Manufacture, China, LTS10770125). The SNP rs2070600 was genotyped in our study using the TaqMan 5′-nuclease assay on the TaqMan 7900HT Sequence Detection System (Applied Biosystems, Foster City, California) under the following conditions: 10 min at 95°C (enzyme activation), 40 cycles at 95°C for 15 s each, and 60°C for 1 min (annealing/extension). The endpoint read was performed for allelic discrimination after amplification. The genotyping procedure was consistent with our previous report. 27 Details regarding primers and probes are provided in Supporting Information, Table S2.

Genetic variation screening
Sequence data were generated by Sanger sequencing of genomic DNA derived from 48 healthy individuals. Polymerase chain reaction (PCR) fragments covering the exons and introns of RAGE (consensus sequence NC_000006.11 GRCh37. p13) were screened using fluorescent dye-terminator cycle, and products were analysed with an Applied Biosystems 3130xl Genetic Analyser (Applied Biosystems, Foster City, CA). The Chromas programme (Technelysium Pty. Ltd., Helensvale, Queensland, Australia) was used to identify putative polymorphisms that were then confirmed by two independent observers. All identified variants were confirmed by repeated sequencing. Details regarding primers are given in Supporting Information , Table S3.

Plasmid construction, cell culture, and transient transfection
Expression plasmids of the predicted nine transcription factors (E2F6, MSC, ZNF263, ID4, TFAP2A, TFAP2C, EGR1, STAT3, and YTHDF1) were constructed with pcDNA3.1 through amplification of human cDNA derived from the above-mentioned human heart tissue. According allele frequency of rs2070600 and ancestral allele in dbSNP, we defined the genotype of rs2070600-GG as wildtype, -AA as mutation. Vectors expressing the RAGE wild-type (rs2070600, -GG allele) or mutant (-AA allele) type were constructed into pcDNA3.1 and flag-pcDNA3.1 using cDNA derived from human peripheral lymphocytes with different rs2070600 genotypes. AC16 and HEK293T cells were used for the following functional assays. Cell culture and transient transfection procedures are described in detail in Supporting Information, Table S4.

Statistical analysis
Data were analysed using SPSS version 24.0 (SPSS, Inc., Chicago, Illinois) for Windows (Microsoft Corp., Redmond, WA). Linkage disequilibrium was calculated using Haploview version 4.1. The polymorphisms were tested for Hardy-Weinberg equilibrium using the χ 2 test. We performed multivariate logistic regression analyses based on different genetic models with adjustment for traditional risk factors to test the association between SNPs and HF risk. The genetic models utilized were as follows: (i) additive model (AA vs. AG vs. GG); (ii) recessive model, AA vs. (AG + GG); and (iii) dominant model (AA + AG) vs. GG. Analysis of clinical prognosis of HF was performed using the Cox proportional hazards model. Significant differences were assessed by either one-way ANOVA followed by Bonferroni's post hoc test or unpaired or paired two-tailed Student's t-test, where appropriate. All biostatistics calculations were performed using Prism (GraphPad). Data are expressed as mean ± SEM of n experiments. All probability values were two-sided, and P < 0.05 was considered significant.

DNA resequencing results
Through RAGE gene resequencing performed in 48 unrelated age-matched and sex-matched, healthy controls from the Han Chinese population, we identified 17 polymorphisms in RAGE.
Of these, three were located in the exon region, and 14 were located in the intron region ( Table 2). All 17 polymorphisms were in Hardy-Weinberg equilibrium in our resequencing population. The results were confirmed using 1000 Genomes Browser (https://www.internationalgenome.org/1000-genomes-browsers). Because the SNPs located in the coding region were more likely to affect the function of proteins by changing amino acids, we chose to evaluate the effect of common missense mutations in RAGE on the prognosis of HF. Rs2070600 in the RAGE coding region leads to the conversion of glycine to serine (G82S), as reported in the Human Gene Mutation Database (http://www.hgmd.cf.ac.uk/ac/index. php).
Region of advanced glycosylation end-product specific receptor SNP rs2070600 G > A had a minor allele frequency (MAF) of 0.236 in the resequencing population (Table 1), 0.225 in HF patients, 0.221 in controls, and 0.233 in non-coronary artery disease patients ( Table 4). According to the dbSNP and 1000 Genomes Browser, we identified guanine 'G' coding G82 to be the major allele (wildtype) and adenine 'A' as the minor allele (mutation) throughout the report.

In silico predictions of functional single-nucleotide polymorphisms
Among the 17 variants identified, 14 were located in the intronic region while only one common missense variant, rs2070600 (MAF > 0.05), was located in an exonic region. This SNP was in exon 3 and showed strong linkage disequilibrium with rs9391855, which is located in intron 8 (D' = 1, r 2 = 1) (Supporting Information, Figure S1).
Region of advanced glycosylation end-product specific receptor rs2070600 variant affects the prognosis of heart failure Among the 3022 HF patients, cardiovascular death or cardiac transplantation occurred in 295 patients (16.7%) in the GG genotype group (n = 1775), 169 patients (15.3%, n = 1091) in the AG genotype group, and 12 patients (7.7%) in the AA genotype group (n = 156) and beta blocker users in 1389 objectors (46%). Cox proportional hazards analysis showed that the rs2070600-A allele was significantly associated with a decreased risk of cardiovascular death and cardiac transplantation in crude [hazard ratio (HR) = 0.49, 95% confidence interval (CI) = 0.28-0.87, P = 0.015] and multivariable-adjusted (HR = 0.53, 95% CI = 0.30-0.94; P = 0.030) recessive genetic models ( Table 3) ( Figure 1A). However, under dominant models, the rs2070600 genotype did not significantly influence the prognosis of HF patients in the crude model (HR = 0.84, 95% CI = 0.70-1.02, P = 0.074) or multivariable model (HR = 0.86, 95% CI = 0.71-1.04, P = 0.113) after adjusting for sex, age, hypertension, diabetes, hyperlipidaemia, smoking, and beta-blocker use. In the additive models, the rs2070600-AA homozygotes showed a higher event-free survival rate compared with survival rate of rs2070600-GG homozygotes in crude (AAvsGG; HR = 0.53, 95% CI = 0.30-0.92; P = 0.024), but not multivariable (HR = 0.59, 95% CI = 0.34-1.03; P = 0.064) models ( Figure 1). After grouped by beta-blocker use, the association between rs2070600 and the prognosis of HF showed no statistical difference in group without beta-blocker use but a trend in group with beta-blocker use in recessive model (Supporting Information, Table S8, Figure S11). Considering the fact that only 46% of the HF patients used beta-blockers, the number of patients in group with beta-blocker use may be too small to detect a statistically significant difference, which needed further investigation. In addition, we surveyed the association between the SNP and HF risk. We carried out a case-control study consisting of 3022 HF patients and 2861 control subjects. Our results showed that the rs2070600 did not significantly influence the susceptibility to HF in recessive model (AA vs. AG + GG, OR = 0.52, 95% CI = 0.66-1.04, P = 0.106) and any of other models ( Table 4).

Increased expression of RAGE in the failing heart
To explore the expression of RAGE in failing heart, we reanalysed the RNA sequencing data from GEO database (GSE57338). The expression of RAGE is higher in the HF group   Figure S9). Furthermore, we detected RAGE expression in heart tissue samples from HF patients (five cases) and healthy donors (three controls) by western blot analysis. The results showed that the protein levels of RAGE were significantly higher in HF patients compared to levels in control subjects. To further confirm the association between the rs2070600 genotype and RAGE protein expression, the DNA was extracted and sequenced. As shown in Figure 2A, there were four major allele homozygotes (rs2070600-GG) and one heterozygote in HF patient group compared with three major allele homozygotes in the healthy donor group. Western blot analysis revealed that RAGE expression increased by 48% in rs2070600 GG allele homozygote HF tissues (P = 0.041).

In vivo effects of rs2070600 variant on messenger RNA levels
To further investigate the hypothesis that rs2070600 could influence gene expression, genotype-dependent RAGE expression was assessed in RNA extracts of peripheral blood lymphocytes from 232 normal individuals. The rs2070600-AA genotype showed significantly lower RAGE mRNA levels compared with the rs2070600-GG genotype ( Figure 2B).
To confirm the influence of rs2070600 on protein expression, plasmids harbouring G or A alleles of rs2070600 were transiently transfected into AC16 cells and HEK293T cells. Reduced RAGE expression was detected in AC16 cells transfected with rs2070600-A allele compared with expression after transfection with rs2070600-G allele ( Figure 2C). This result was reiterated in 293 T cell lines ( Figure 2D) and HUVEC cell lines (Supporting Information, Figure S3). To eliminate the possibility that rs2070600 could affect the binding affinity of the antibody (ABclonal, China, A13264), which contains a sequence corresponding to amino acids 41-340 of human AGER (NP_001127.1) (https://abclonal. com.cn/catalog/A13264), a pcDNA3.1 vector containing Flag-pcDNA3.1-RAGE-rs2070600-G or -A was transfected into AC16 and 293 T, and an anti-FLAG antibody was used to detect the level of RAGE protein. The results showed that RAGE expression of the rs2070600-A allele was lower than that of the -G allele, which was consistent with the above results (Supporting Information, Figure S6).

Micro RNAs modulate RAGE expression
To further investigate the mechanism by which rs2070600 could influence the expression of RAGE, bioinformatics analysis and cell experiments were performed. Bioinformatics Figure 1 The association between genotypes of rs2070600 and the outcome of heart failure in three genetic models. (A-C) Survival analysis of different genotypes on the prognosis of heart failure using Cox proportional hazards analysis after adjusted for traditional risk factors (sex, age, hypertension, diabetes, hyperlipidemia, smoking status, and beta-blocker treatment) in recessive (A), dominant (B), and additive model (C), respectively. Cox proportional hazards analysis showed the association of genotypes of rs2070600-A allele with cardiovascular deaths or cardiac transplantation in recessive model.   Mir-125a-3p directly targets the rs2070600-A allele and mediates allelic expression. (A) Schema graph to illustrate the loci of rs2070600 G/A variant occurs at the miR-125a-3p binding site. Mir-125a-3p co-action rs2070600-A allele which in the coding region sequence of the gene to regulate the expression of RAGE. T base-paired with U in the Watson-Crick mode (solid line), whereas allele G did not (shown without line). (B) miR-125a-3p negatively regulates the protein level of RAGE in AC16 analysed by Western blotting. (C-D) PcDNA3.1-RAGE-rs2070600-G (C) or -rs2070600-A (D) were co-transfected with miR-125a-3p-mimics, negative control miRNA (miR NC), miR-125a-3p-inhibitor, or inhibitor negative control (Inhibitor NC) into AC16. Overexpression of MiR-125a-3p significantly reduced the level of RAGE in pcDNA3.1-RAGE-rs2070600-A transfected cell, while the expression of RAGE increased after inhibition of MiR-125a-3p. No obvious differences were observed in AC16 transfected with pcDNA3.1-RAGE-rs2070600-G.

3568
S. Li et al.
analysis revealed that rs2070600 was located in the potential binding site of several transcription factors including E2F6,  MSC, ZNF263, ID4, TFAP2A, TFAP2C, EGR1, STAT3, and  YTHDF1. However, no significant changes in RAGE expression were observed in AC16 cells transfected with the above nine transcription factors (Supporting Information, Figure S4).
Considering that a common variant in the coding region could affect the binding affinity of microRNA, 28 we then focused on searching for the potential microRNA that could provide the interpretation. Seven candidate miRNAs were predicted to be influenced by rs2070600 using bioinformatics analysis (Supporting Information, Table S5).
To identify the functional miRNAs, the pcDNA3.1-RAGE-rs2070600-A construct was co-transfected with the predicted miRNAs into AC16 cells. As shown in Figure 3A, B, miR125a-3p displayed efficient suppression of RAGE expression, while the other six miRNAs did not significantly decrease the expression of RAGE (Supporting Information, Figure S5). We further explored the effects of the association between miRNA-125a-3p overexpression and inhibition of RAGE expression. The PcDNA3.1-RAGE-rs2070600-G o r pcDNA3.1-RAGE-rs2070600-A constructs were co-transfected with either 100 nmol of has-miR-125a-3p, has-miR-125a-3p inhibitor(s), or miR-negative control (RIBOBIO Co., Ltd, Guangzhou, China). Overexpression of miRNA-125a-3p downregulated RAGE expression level, while inhibition of miR-125a-3p significantly increased the expression of RAGE in AC16 cells transfected with pcDNA3.1-RAGE-rs2070600-A. However, no effects of miRNA-125a-3p or inhibitor on RAGE expression were observed in AC16 cells transfected with pcDNA3.1-RAGE-rs2070600-G ( Figure 3B-C). Additionally, the relative level of miR-125a-3p was detected in peripheral blood lymphocytes, and we found no difference among participants carrying the rs2070600-GG, AG, or AA genotypes (Supporting Information, Figure S6). Furthermore, compared with health control, an increased tendency expression of mir-125a-3p in failing heart was found, but it was no statistical difference (Supporting Information, Figure S7A). To investigate how inflammatory stimuli influence the miR-125a-3p expressing in different cells types, we used 100 ng/L lipopolysaccharide to treat different cells (AC16, 293 T, HUVEC). We had a consistent result that lipopolysaccharide can induce mir-125a-3p high expression in intervention groups compared to controls. The details were supplied in Figure S7B-D. Additionally, the DNA sequences proximal to rs2070600 in AC16(A) and 293 T(B) cell lines were offered in Figure S8.

Discussion
We identified a common SNP (rs2070600) in the coding region sequence (CDS) of RAGE significantly associated with the prognosis of HF in a prospective, observational, single-centre study consisting of 3394 HF patients and 2861 control participants. The SNP affects the expression of RAGE through interacting with miR-125a-3p.
Chronic inflammation contributes to the progression of HF. 29 RAGE is a key mediator of inflammatory response. Under normal conditions, it plays a small role in normal physiological development and in healthy homeostasis. After acute infection or massive injury, RAGE is active and releases proinflammatory factors, 30 directly recruits leucocytes, and detects pathogen DNA. 31 RAGE aggravated the inflammatory response and tissue damage in HF following induction of amassing ligands. 32 Persistent activation of RAGE/nuclear factor-κ-gene binding modulates pro-inflammatory cytokine expression, such as interleukin (IL)-1, IL-6, and tumor necrosis factor-α, which leads to the activation of cardiomyocyte death pathways. 33 Plasma soluble RAGE levels were obviously higher in the rs2070600-GG genotype of Korean subjects compared with the AA and AG genotypes. 34 As a protein biomarker, soluble RAGE was associated with all-cause mortality in the Framingham Heart Study (n = 3523). 35 From the above studies, we could infer that the rs2070600 haplotype may be linked to the outcomes of cardiovascular diseases. Our study further uncovered the direct relationship between rs2070600 genotype and prognosis of HF. Based on our case-control study, rs2070600 was not associated with the risk of HF, which was in line with the Framingham offspring study 36 where no association was observed between rs2070600 and the incidence of cardiovascular disease in diabetic or non-diabetic subjects (n = 1632). In brief, although the common missense variant (rs2070600) of RAGE does not cause HF, it may play an important role in the progression of HF.
We identified an exclusive missense mutation, rs2070600, in the RAGE CDS region that has a MAF of 0.232. Further bioinformatics analysis strongly suggested that rs2070600 is a functional variant. By querying the GTEX database (https:// gtexportal.org), it was confirmed that rs2070600 is the expression quantitative trait loci. In a study exploring the association between rs2070600 and lung function, Miller et al. 19 showed that the major allele of rs2070600 is linked to higher FEV1 and FEV1/FVC in 1024 U.K. participants. Rs2070600 (G82S) is a non-synonymous polymorphism that changes the acid sequence of the encoded protein from glycine to serine. A series of assays were carried out to investigate the underlying mechanisms. In vivo gene expression analysis of human lymphocytes revealed that the rs2070600-AA genotype decreased RAGE expression compared with RAGE expression resultant from the rs2070600-GG or rs2070600-AG genotypes, which was consistent with the results from cell experiments that RAGE expression was lower in cells transfected with pcDNA3.1-rs2070600-A compared with the expression in cells transfected with pcDNA3.1-rs2070600-T allele. Taken together, we confirmed the predicted result of variant rs2070600 that the minor allele could decrease RAGE protein levels.
Multiple studies have shown that mRNAs can bind to CDSs and effectively inhibit translation. 37,38 The current study revealed a possible underlying mechanism by which miR-125a-3p directly targets the rs2070600-A allele and decreases the expression of RAGE. Our result was consistent with a recent report by Fu et al. 39 that the characteristic of miRNAs that directly targeted the CDS were prone to extensive base-pairing at the 3′ seed. However, the specific role of miR-125a-3p in the pathogenesis of HF should be investigated in future studies. To our knowledge, we have discovered the RAGE gene variant locus of a disease-associated miR-SNP interaction found within the CDS of RAGE. Furthermore, assessing whether rs2070600-AA genotype risk is mingled by environmental factors or other common RAGE SNPs is necessary to explore the regulatory controls that impact RAGE expression.
Several limitations of our study should be mentioned. The first limitation of our study is that this is an observational study and that the positive association between RAGE variants and HF needs to be confirmed in prospective cohort studies. Second, the HF patients in our study were recruited mainly from Wuhan Tongji Hospital, and the results need to be validated in further multicentre studies including larger HF populations. Third, we cannot rule out the possibility that other variants in linkage disequilibrium with rs2070600 may account for the association with the prognosis of HF. Finally, other regulatory factors may be involved in the regulation of RAGE gene expression.

Conclusions
The genetic and functional evidence indicated that rs2070600 in the coding region of RAGE is associated with the prognosis of HF in the Han Chinese population. The rs2070600-AA genotype shows a favourable prognosis for HF compared with the rs2070600-G allele carriers, which is dependent on beta-blocker use. A mechanistic study demonstrated that miR-125a-3p could directly target the rs2070600-A allele to decrease RAGE expression. Targeting RAGE is an attractive way to decrease HF-associated mortality and improve the prognosis of HF in the future.

Supporting information
Additional supporting information may be found online in the Supporting Information section at the end of the article.  Table S8. Stratified analysis of beta blocker user for analysis the association rs2070600 with heart failure outcomes Figure S1. Linkage disequilibrium (LD) structure and haplotype blocks of the RAGE gene. LD (D') for identified polymorphisms in RAGE, as generated by Haploview version 4.1 from the genotype data of 48 random controls. Haplotype blocks derived from these genotypes using the solid spine LD setting are outlined in black. Figure S2. The effect of rs2070600 genotype on local mRNA structure. RNA structure of~250-nts of RAGE mRNA flanking rs2070600 (red open-circle) was predicted using UNAfold. Results revealed a genotype-specific structure of RAGE mRNA. The table summarizes unfold microRNA targeting prediction results which indicate that rs2070600 A allele is energetically more favorable. Figure S3. Rs2070600 genotype influences expression of RAGE in HUVEC cells. The rs2070600 mutant A allele displayed significantly reduced protein expression of RAGE in HUVEC cells. The data are presented as mean ± SD from three independent experiments. Figure S4. Interaction between transcription factors and RAGE expression in AC16 cells. The predicted transcription factors didn't significantly influence the expression of RAGE. AC16 cells were transfected with pcDNA3.1 empty vector or pcDNA3.1-TFAP2C construct, pcDNA3.1-ZNF263 construct, pcDNA3.1-YTFDH1 construct (A), or pcDNA3.1-ID4 construct, pcDNA3.1-E2F6 construct (B), or pcDNA3.1-MSC construct (C), or pcDNA3.1-TFAP2A construct (D), or pcDNA3.1-STAT3 construct (E), or pcDNA3.1-EGR1 construct (F). The data are presented as mean±SD from three independent experiments. N.S., no significant. Figure S5. Interaction between miRNA and RAGE expression in AC16 cells. Six of the predicted miRNAs could not significantly suppress the expression of RAGE. AC16 cells were transfected pcDNA3.1-RAGE-rs2070600-A 300ng with 100 nmol miR-Negative Control (RIBOBIO Co., Ltd, Guangzhou, China) and other 8 miRNAs. The data are presented as mean ± SD from three independent experiments. NS, no significant. Figure S6 Flag label plasmids of pcDNA-rs2070600-G (RAGE-WT) and -rs2070600-A (RAGE-MU) transfected into AC16 and 293T cells, and the relative level of miR-125a-3p in peripheral blood lymphocytes. Figure S7 Relative miR-125a-3p expressed of vivo and vitro (A) Relative miR-125a-3p expression in failing heart and health controls measured by real-time PCR. Data are expressed as mean ± SEM. LPS 100ng/L induced miR-125a-3p expressed in vitro. Relative miR-125a-3p expression in AC16 (B), 293T(C) and HUVES(D). Figure S8. DNA sequences proximal to rs2070600 in AC16(A) and 293T(B) cell lines. Figure S9 Differential of RAGE/Endogenous RAGE ligands gene expression analysis in GSE57338 Figure S10 Relative RAGE expression in vitro normalize by GAPDH, RAGE mRNA expression levels was detected by qPCR, and GAPDH was used as endogenous control to mRNA. (A) AC16,(B)293T. Figure S11 Effects of rs2070600 on the prognosis of HF patients grouped by β-blocker use under different models. By Cox proportional hazards analysis after adjusted for traditional risk factors (sex, age, hypertension, diabetes, hyperlipidemia, smoking status), the link of rs2070600 with heart failure outcomes showed no statistical significance in group with and without β-blocker use in recessive(A-B), dominant (C-D) and additive models(E-F).