Abstract

Nicotinamide N-methyltransferase (NNMT) catalyzes the methylation of nicotinamide. Our previous works indicate that NNMT is involved in the body mass index and energy metabolism, and recently the association between a SNP (rs694539) of NNMT and a variety of cardiovascular diseases was reported. At present, more than 200 NNMT single nucleotide polymorphisms (SNPs) have been identified in the databases of the human genome projects; however, the association between rs694539 variation and hyperlipidemia has not been reported yet, and whether there are any SNPs in NNMT significantly associated with hyperlipidemia is still unclear. In this paper, we selected 19 SNPs in NNMT as the tagSNPs using Haploview software (Haploview 4.2) first and then performed a case-control study to observe the association between these tagSNPs and hyperlipidemia and finally applied physiological approaches to explore the possible mechanisms through which the NNMT polymorphism induces hyperlipidemia. The results show that a SNP (rs1941404) in NNMT is significantly associated with hyperlipidemia, and the influence of rs1941404 variation on the resting energy expenditure may be the possible mechanism for rs1941404 variation to induce hyperlipidemia.

1. Introduction

Nicotinamide N-methyltransferase (NNMT) catalyzes the methylation of nicotinamide (NAM) using S-adenosylmethionine (SAM) as a methyl donor to generate methylnicotinamide (MNA) [1, 2]. The metabolic process and products of nicotinamide are highly related to various cardiovascular diseases, such as coronary heart disease, stroke, atherosclerosis, and diabetes [3, 4], which are closely related to an energy metabolism imbalance or obesity. Our previous studies showed that the nicotinamide metabolism rate was significantly correlated with the body mass index (BMI) [5] and energy metabolism [6], and these works have recently been supported by multiple reports. Kraus et al. [7] reported that NNMT expression is increased in white adipose tissue (WAT) and liver of obese and diabetic mice, and NNMT knockdown in WAT and liver protects against diet-induced obesity by increasing cellular energy expenditure. Liu et al. [8] found that serum MNA was associated with obesity and diabetes. Hong et al. [9] showed that nicotinamide N-methyltransferase regulated hepatic nutrient metabolism. Additionally, when NNMT catalyzes the methylation of nicotinamide, S-adenosyl homocysteine (SAH) and homocysteine (Hcy) are generated [10]. Hyperhomocysteinemia is one of the independent risk factors for cardiovascular diseases [11]; therefore, the relationship between NNMT and cardiovascular diseases has recently been increasingly reported [1214].

Since the expression of NNMT is directly determined by NNMT, the roles of NNMT in the development of cardiovascular diseases have been reported in recent years. Bubenek et al. [14] reported that the occurrence and development of peripheral arterial occlusive diseases were closely related to NNMT expression and the serum NNMT level, and the NNMT expression level was significantly positively correlated with the low density lipoprotein level and significantly negatively correlated with the high density lipoprotein level. At present, more than 200 NNMT single nucleotide polymorphisms (SNPs) in noncoding region have been identified in the human genome projects, but only one SNP (rs694539) has been reported in the genetic association studies. Souto et al. [11] reported that rs694539 variation is significantly associated with serum Hcy level in a Spanish population. Later, the associations between a variety of cardiovascular diseases and rs694539 variation were reported. van Driel et al. [15] reported that the risk of congenital heart diseases increased by eightfold in rs694539 AG+GG carriers under the conditions of low nicotinamide intake and drug exposure. Giusti et al. [16] reported that rs694539 variation was related to abdominal aortic diseases. de Jonge et al. [17] reported the association of rs694539 variation with lymphoblastic leukemia in children. Sazci et al. [1820] reported the association of rs694539 variant with nonalcoholic steatohepatitis, bipolar disorder, and epilepsy, respectively. However, the association between rs694539 variation and hyperlipidemia has not been reported yet, and whether there are any SNPs in NNMT significantly associated with hyperlipidemia is still unclear.

Generally, hyperlipidemia is thought to be a major factor to induce many cardiovascular diseases. In this paper, we selected 19 SNPs (including rs694539) as the tagSNPs from NNMT DNA sequence in the database of 1000 Genomes Project using Haploview software (Haploview 4.2) first and then performed a case-control study to observe the association between these tagSNPs and hyperlipidemia and finally applied physiological approaches to explore the possible mechanisms through which NNMT polymorphism induces hyperlipidemia.

2. Subjects and Methods

2.1. Subjects

Cases and controls were all recruited from unrelated Chinese Han ethnicity volunteers. The investigation was approved by the local ethics committee at Jiangxi Normal University, and all participants gave the written informed consent. This study conforms to the latest revision of the Declaration of Helsinki.

The demographic and clinical characteristics of the 395 hyperlipidemic patients and the 316 controls are reported in Table 1. Usually, levels of plasma lipids increase with age, but after the age of 60 years, the trend is reversed. Levels of plasma lipids decrease with the increasing of age in the people aged ≥ 60 years. If a person aged ≥ 60 years has never had hyperlipidemia, then he/she has little chance to be of hyperlipidemia in the future, which means those people are the real nonsusceptible population to hyperlipidemia. Therefore, we recruited the controls from the volunteers aged ≥ 60 years. Meantime, simple hyperlipidemia usually occurs before the age of 60 years, and most hyperlipidemic patients aged ≥ 60 years have one or more related diseases, such as hypertension and diabetes. To enroll simple hyperlipidemic patients, we recruited the cases from the hyperlipidemic patients aged < 60 years. Finally, 395 patients aged < 60 years with simple hyperlipidemia (physician diagnosis, a recorded serum total cholesterol level > 220 mg/dL, and/or serum triglyceride level > 150 mg/dL) and 316 controls who were aged ≥ 60 years and had never been diagnosed as hyperlipidemia were enrolled. Exclusion criteria for the cases were poorly controlled diabetes mellitus (blood glucose > 6.67 mmol/L and/or glycosylated hemoglobin > 6.0%) and a possible secondary dyslipidemia, including thyroid and liver disease, renal failure, and proteinuria [21]; exclusion criteria for the controls were diabetes mellitus (hypoglycemic treatment and/or fasting blood glucose > 7 mmol/L), hypertension (hypertension treatment and/or blood pressure > 140/60 mmHg), or current treatment with lipid-affecting drugs [21].

To explore the roles of NNMT SNPs in the occurrence of hyperlipidemia, after the case-control studies, in which the SNPs significantly associated with hyperlipidemia and the susceptible genotypes were observed, we recruited an additional 150 healthy Chinese Han ethnicity male college students (aged 17–23 years and without any diagnosed diseases) to genotype the same SNPs. Then, these 150 subjects were divided into the susceptible group and the nonsusceptible group according to their genotypes and carried out the physiological studies.

2.2. Database Query and Determination of Investigation Sites

Nineteen NNMT tagSNPs were selected from the 208 known genotyped SNPs occurring within the Chinese Han population (CHB+CHS), contained within the database of the 1000 Genomes Project (http://browser.1000genomes.org/), with Haploview software (Haploview 4.2) after configuring the criteria (MAF > 0.10 and ) in this study. The positions of these 19 tagSNPs in relation to NNMT are shown in Figure 1.

2.3. Gene Polymorphism Detection
2.3.1. Main Instruments and Reagents

The main instruments included a PCR machine (Norwalk, CT. 06859 USA), an electrophoresis apparatus (Beijing Junyi Electrophoresis Co., Ltd.), an automatic UV-Visible spectrometer, a biological electrophoresis image analysis system (Shanghai Furi Science and Technology Co., Ltd.), and a sequencer (ABI). The main reagents included PCR primers, dNTPs (Shanghai Hanyu Biological Engineering Co., Ltd.), a Taq enzyme system, ddH2O (Shanghai Biowing Applied Biotechnology Co., Ltd.), 1.5 mL Eppendorf tubes, pipette tips of all sizes, and 96-well PCR plates (Haimen Yonghui Experimental Equipment Co., Ltd.).

2.3.2. DNA Extraction and Primer and Probe Design

Venous blood was collected, and genomic DNA was extracted using a DNA extraction kit (Promega, USA) and then stored at −80°C for future use. The NNMT gene sequence was downloaded from the NCBT database. Primers and probes were designed using Primer 3 online (Version 0.4.0, http://bioinfo.ut.ee/primer3-0.4.0/).

2.3.3. SNP Genotyping

Polymerase chain reaction-ligase detection reaction (PCR-LDR) was used to detect the genotypes of each SNP. The principle of PCR-LDR is to use a high-temperature ligase to recognize gene polymorphism loci. PCR-LDR first utilizes a multiplex PCR to obtain gene fragments that contain the mutation sites to be measured, after which a multiplex LDR is performed, and finally, the fluorescent products of LDR were differentiated by ABI sequencer (PRISM 3730). To check the reliability of the genotyping results, 10% of PCR-LDR reactions were carried out in duplicate, and more than 99.5% had matching results. Additionally, 50 samples were randomly selected to genotype the most significantly associated SNP (rs1941404) identified in the association analysis using Sanger sequencing method, and 100% matched up perfectly with the results obtained with PCR-LDR method.

2.4. Physiological Tests
2.4.1. Body Fat Percentage (BF%)

BF% was measured with bioimpedance measurement using an X-SCAN PLUS body composition analyzer (X-SCAN PLUS II, Jawon Medical Co., Ltd., Republic of Korea). Tests were carried out in the morning on an empty stomach.

2.4.2. Resting Energy Expenditure (REE), Resting Energy Expenditure per Unit Body Surface Area (REEU), and Respiratory Quotient (RQ)

150 subjects were randomly selected from Chinese Han male college students to carry out a respiratory gas exchange analysis using a computerized indirect calorimetry with a ventilated-hood system (Metalyzer 3B, CORTEX Biophysik GmbH, Leipzig, Germany), as described previously [22], but one of them quitted for personal reason, and only 149 of them completed the test. The measurements were performed in a quiet room at 24°C, between 8:00 and 9:00 a.m., after an overnight fast of ≥12 h. The subjects sat in chairs quietly for 30 min before measurements were started and kept their position until measurements were completed. REE was calculated from measured oxygen consumption and carbon dioxide production according to the Weir equation: REE = 3.9 × O2 used + 1.1 × CO2 produced [23]. REEU was calculated as REE divided by body surface area (BSA), which was calculated using the Stevenson equation [24]: BSA (m2) = 0.0061 × height (cm) + 0.0128 × weight (kg) − 0.1529. RQ was calculated as VCO2/VO2.

2.5. Statistics

The distributions of allele and genotype frequencies, Hardy-Weinberg equilibrium (HWE), were conducted online using http://analysis.bio-x.cn/ [25]. Two classification logistic regressions were performed in the analyses of genotype effects and genetic models. The values of continuous variables were expressed as the mean value ± standard deviation (SD) and compared with independent -test using IBM SPSS Statistics 20.0 (SPSS Inc., Chicago, IL, USA). All values were two-sided, and values < 0.05 were considered to be statistically significant. Bonferroni correction was performed for multiple comparisons, and the corrected significance level was .

3. Results and Discussion

3.1. NNMT Genotyping and Allele Frequency Analysis

The distributions of the genotype and allele frequencies of the 19 SNPs between the case and control groups are shown in Table 2. There are 5 SNPs (rs694539, rs12285641, rs11214926, rs2244175, and rs4646335) showing in the HWE tests on the control group, which means those 5 SNPs are not in HWE and should not be considered in genetic association studies. Among the other 14 SNPs, which are in HWE (), 3 of them (rs10891644, rs4646335, and rs1941404) show in the comparisons of genotype frequencies. However, after Bonferroni correction, only rs1941404 shows the corrected significant differences between the groups (), which indicates that rs1941404 is the most significantly associated tagSNP in NNMT with hyperlipidemia. At this locus, the case group exhibits a higher frequency of allele C and a higher frequency of genotype CC than those of the control group, respectively.

Usually, when a locus is not in HWE, then this suggests deviations from random mating, population substructure, migration, natural selection, small population sizes, or genotyping errors. Considering that the genotyping results of this study have been verified with 10% samples tested in duplicate and validated using Sanger sequencing method, genotyping errors can be excluded, but what exactly makes the departure from HWE at those 5 SNPs is still unclear. Therefore, those 5 SNPs were not considered for further analysis. As mentioned previously, rs694539 is the most studied locus in NNMT and has been reported to be associated with many cardiovascular diseases. However, in this study, the distributions of genotypes and allele frequencies of rs694539 do not show significant differences () between the case and control groups (Table 2). The distribution of SNP allele frequencies has obvious population differences. Existing reports on the association between rs694539 variation and the diseases all employed a Spanish, Japanese, or Israeli population as the research subjects, whereas this study recruited a Chinese Han population. The differences in the study subjects might be the main reason for the different results obtained in the existing reports and this study. Moreover, although many diseases have been reported to be associated with rs694539, the association between rs694539 and hyperlipidemia has not been reported yet, suggesting that possibly the rs694539 is not the SNP which has the most significant regulatory effects on hyperlipidemia. Therefore, we do not analyze rs694539 further, and the subsequent analyses are primarily focused on rs1941404.

3.2. Genotype Effects of the rs1941404 Variation on Hyperlipidemia

Genotype effects usually are analyzed using the genetic models, recessive, dominant, additive, and codominant. However, these genetic models are not independent; sometimes several genetic models show the significance () at the same time. To avoid the hash of possible genetic models and precisely determine the inheritance mode, Zintzaras and Santos [26] provided a whole solution and introduced the degree of dominance (). As shown in Table 3, the chances of the subjects with the genotypes CC, CT, and TT to have hyperlipidemia are 70%, 48%, and 50%, which demonstrates that the homozygote CC has the greatest risk for hyperlipidemia and the risks of the genotypes TT and CT are almost equal, thus indicating that the genotype effects of the rs1941404 variation are a recessive inheritance mode. The results of two classification logistic regression analyses (Table 3) show that the recessive model and additive model are significant (), which further confirms that the genotype effects of the rs1941404 variation are a recessive inheritance mode, because the differences between the homozygote CC and TT carriers are significant () and the OR values show that the risk to be hyperlipidemia for the CC is 2.432 and 2.316 times higher than those for the TT+CT and the TT, respectively. Table 3 also shows that the codominant model is significant () which suggests that the heterozygote CT may have a significantly different risk of being hyperlipidemia from the homozygote CC+TT, and the values of (−1 < < 0) show that the risk of disease for CT is in the middle of the two homozygotes and close to TT. Therefore, genotype effects of rs1941404 variation basically belong to a recessive inheritance mode, and the CC carriers are the hyperlipidemia susceptible population and the CT or TT carriers are the nonsusceptible population, and therefore in the following exploration of the mechanisms, the differences between the CC and the CT+TT will be focused on.

3.3. Physiological Exploration of the Possible Mechanisms

Medically, hyperlipidemia is the one symptom of metabolic syndrome, which is thought to be caused by an underlying disorder of energy utilization and storage [27], so that BF% and the characters of the resting energy expenditure (REE) may be able to offer some clues for the exploration of the possible mechanisms, through which rs1941404 variation induces the hyperlipidemia. As shown in Figure 2, BF% of the CC carriers is substantially equivalent to that of the (CT+TT) carriers (). However, the REE and REEU of the CC genotype are significantly lower than those of the (CT+TT) genotype () and the RQ is significantly higher than that of the (CT+TT) genotypes ().

Theoretically, it is not surprising that NNMT variation has a great impact on an organism’s energy metabolism. The physiological function of NNMT is to transfer a methyl group from SAM and catalyze the methylation of NAM. Firstly, NAM is the precursor of NAD+. NAM methylation directly affects the level of NAD+, which is the key coenzyme for energy metabolism and is directly involved in the entire process of aerobic oxidation. Therefore, changes in NAD+ level will directly affect the energy metabolism process. Secondly, the changes of NNMT level may indirectly affect the expression of various related proteins. For example, Zhang et al. analyzed the SW480 cell lines before and after NNMT transfection [28] and found that more than 30 genes demonstrated greater than twofold differences in gene expression before and after transfection, of which most genes were related to glucose and lipid metabolism and the oxidative respiratory chain.

Epidemiological studies have shown that the hyperlipidemia has a significant association with obesity or energy metabolism imbalance. The results of the present study suggest that although the BF% of the CC genotype subjects was not higher than that of the CT+TT genotype subjects (Figure 2(a)), the probability of the CC genotype subjects to be fat in the future is higher than that of the CT+TT genotype subjects, because the CC genotype subjects have the lower REE and REEU values (Figures 2(b) and 2(c)), which are the most important factors leading to obesity. The RQ value reflects the proportion of the energy supply from different energy materials. The RQ value of the CC genotype subjects is higher and very close to 1, which indicates that before the appearance of the hyperlipidemia these subjects utilize a relatively high proportion of glucose to generate energy other than lipids. Therefore, the effects of rs1941404 variations on REE may be a primary physiological mechanism leading to the development of the hyperlipidemia.

4. Conclusion

The rs1941404 variation in NNMT is significantly associated with hyperlipidemia and the CC is the risk genotype (recessive inherence mode). The influences of rs1941404 variation on the resting energy expenditure may be the possible physiological mechanisms for rs1941404 variation to induce hyperlipidemia.

Competing Interests

The authors declare that there is no conflict of interests.

Acknowledgments

This work was supported by the National Science Foundation of China (NSFC 21365013), by the Science Project from the Jiangxi Provincial Education Department (GJJ 13240), and by the Open Project Program of Key Lab of Training, Monitoring and Intervention of Aquatic Sports of General Administration of Sport of China, Jiangxi Normal University (no. 201604).