1. Introduction
Today, metabolic diseases are one of the most important health problems worldwide, with obesity, diabetes, and metabolic syndrome affecting most of the population. Obesity is classified as an abnormal accumulation of fat that causes organic changes. Approximately 10–15% of the world’s population suffers from obesity; in Mexico, 75.2% of the population is obese or overweight, and the state of Sinaloa is among the 10 Mexican states with the highest prevalence of this disease [
1,
2]. Metabolic syndrome is a clinical condition that includes central and abdominal obesity, systemic hypertension, insulin resistance, and atherogenic dyslipidemia, and it is estimated that 34% of people with obesity have metabolic syndrome [
3]. Metabolic diseases are classified as multifactorial diseases in which genetic variants such as single-nucleotide polymorphisms (SNPs) may significantly contribute to the development of these diseases.
The study of the genetics of metabolic diseases began with genome-wide association studies (GWASs), which found a large number of SNPs (more than 120 genes) associated with obesity or metabolic syndrome, mainly in European populations [
4]. For example, the proprotein convertase subtilisin/kexin type 1 (PCSK1) gene is responsible for the activation of several neuropeptides and hormones involved in thermogenesis and feeding. In fact, the PCSK1 SNP is the third most important gene associated with obesity; this gene is one of the first genes recognized as a monogenic cause of obesity [
5,
6,
7]. The zinc finger BTB domain-containing 16 (ZBTB16) gene is involved in the thermogenic process, as well as fatty acid oxidation in brown adipocytes, and genetic variants of ZBTB16 have been associated with changes in anthropometric parameters and lipid levels in the population. They are a risk factor for developing metabolic diseases [
8,
9]. ZPR1 is a zinc finger protein that acts as a transcription factor associated with peroxisome proliferator-activated receptor γ; this protein also interacts with apolipoprotein A5, which helps regulate postprandial triglycerides. SNPs in this gene are associated with elevated triglyceride and glucose levels and altered insulin sensitivity, and they have been implicated in several metabolic disorders (e.g., metabolic syndrome and type 2 diabetes mellitus) [
10,
11,
12]. The peroxisome proliferator-activated receptor γ (
PPARG) gene is a transcription factor that regulates genes involved in metabolic homeostasis and energy expenditure and storage; in addition,
PPARG is the master regulator of adipocyte differentiation, and it modulates metabolism and inflammation in immune cells [
13]. Its polymorphism has been associated with altered energy expenditure, leading to increased body mass index (BMI) and obesity [
14]. The transmembrane protein 18 (TMEM18) gene is a nuclear membrane-localized protein in various brain regions, including the hypothalamus, that controls feeding behavior, and its function is related to the transcriptional repressor that binds to a specific ssDNA sequence. SNPs in this gene are associated with changes in satiety and obesity, but the mechanisms remain unclear [
15,
16].
Studies of allelic variants of PCSK1 rs6235 G allele, TMEM18 rs6548238 C allele, GPX5 rs445870 G allele, ZPR1 rs964184 G allele, ZBTB16 rs7106340 T allele, and PPARG1 rs3856806 T allele are limited or absent in adults from northwestern Mexico, and their association with metabolic diseases is not well known. Furthermore, the possible associations found in this study could help to detect genetic risk factors for the development of metabolic diseases (obesity and/or metabolic syndrome) in our population; moreover, we could lay the groundwork for new specific treatments against metabolic diseases for our population taking into account genetic profiles. Therefore, this work aimed to identify individual associations of six previously described SNPs with anthropometric and biochemical traits and with metabolic diseases such as obesity or metabolic syndrome in a population from northwestern Mexico.
4. Discussion
Obesity and metabolic syndrome are multifactorial diseases in which genetic components may play an important role. In this work, we found that the allelic variants PCSK1 rs6235 G allele, TMEM18 rs6548238 C allele, PPARG1 rs3856806 T allele, and GPX5 rs445870 G allele were associated with anthropometric and biochemical changes in subjects from the state of Sinaloa. In addition, the PCSK1 rs6235 G allele was associated with obesity and metabolic syndrome, and the PPARG1 rs3856806 T allele was associated with obesity.
The allele frequencies of some SNPs In this work were higher compared to other ethnic groups, such as the G allele of
PCSK1 rs6235 (19.99%) and the G allele of
GPX5 rs445870 (33%) [
21,
22], but the frequencies of other SNPs were similar to those of different populations, such as the C allele of
TMEM18 rs6548238 [
23,
24]. The differences in SNP frequencies compared to our study may be because the other studies were conducted in different populations, such as Asians or Europeans.
The presence of some SNPs has been observed to be associated with changes in biochemical or anthropometric traits in some populations, and this study was no exception. The G rs6235 allele of
PCSK1 was associated with increased anthropometric traits, including BMI, weight, waist circumference, and WHR. These results are consistent with other studies; for example, in the Chinese population, an allelic variant of
PCSK1 rs6234 was associated with increased BMI and waist circumference [
25,
26], and, in the young European-American population, these SNPs were associated with increased BMI [
27]. Similar results were found in British and Mexican subjects [
28,
29]. The increase in BMI, weight, waist circumference, and HR associated with the
PCSK1 SNP in this population could be because this SNP is strongly associated with obesity due to alterations in thermogenesis and feeding, among other metabolic changes, increasing the anthropometric characteristics mentioned above [
26].
On the other hand, no previous study has reported an association between the C allele of TMEM18 rs6548238 and increased systolic blood pressure in adults, as observed in this study; this increase in systolic blood pressure may be related to the ability of TMEM18 to cause obesity, but the mechanisms are still unclear.
Regarding biochemical parameters, several studies have shown that
PPARG1 SNPs could increase cholesterol levels in different populations, and our results are consistent with these studies [
30]. This association could be because
PPARG is related to cholesterol efflux in cells, and altering this gene could increase blood cholesterol levels [
31]. Triglyceride level was another trait associated with the C allele of
TMEM18 rs6548238 that was altered in this study, a phenomenon also observed in the Chinese population [
32]; this effect may be related to its ability to induce obesity. As for the G allele of
GPX5 rs445870, there are no previous reports on its association with reduced HDL cholesterol in adults, and its mechanism is unknown.
Regarding metabolic diseases, in this study, we analyzed the individual associations of six SNPs with obesity and metabolic syndrome, among which the G rs6235 allele of
PCSK1 was associated with obesity. This association has also been identified in different populations, such as European subjects. A study conducted on Danish subjects analyzed 3074 obese and 2790 non-obese subjects and found associations between
PCSK1 rs6232 and rs6235 and obesity. In addition, in Swiss subjects, an association between rs6235 and rs6232 and obesity was found in a comparison of 551 obese Swiss class III individuals and 542 randomly selected blood donors. A comparison between 532 non-obese French young adults and 505 obese French children showed an association between
PCSK1 rs6232 (OR: 1.57,
p: 0.009) and rs6235 (OR: 1.50 (1.23,
p: 0.00003) and obesity, results consistent with this study [
33]. Reported associations between
PCSK1 allelic variants (rs6232 and rs6235) in European-American (OR: 1.71,
p: 0.018) and African-American (OR: 1.47,
p: 0.018) individuals and obesity, but not in Hispanic individuals [
27]. In addition, allelic variants of
PCSK1 (rs6232 and rs6235) were also studied in Mexican subjects; 802 non-obese and 404 obese children and 562 non-obese and 614 obese adults from Mexico City were analyzed. The rs6232 SNP was associated with obesity in children (OR: 3.78,
p: 0.0003) and with class III obesity in adults (OR: 2.61,
p: 0.02), consistent with our findings [
21]. Furthermore, associations between
PCSK1 allelic variants and obesity have been reported in Chinese and Taiwanese subjects [
25,
34,
35].
Mutation in the
PCSK1 gene is sufficient to cause obesity, as this gene is part of a family of genes classified as responsible for monogenic obesity; in fact, alterations in this gene are associated with early onset obesity, severe intestinal malabsorption, and hyperphagia [
36]. The strong association between
PCSK1 and obesity is due to its function; mutation in this gene leads to the misprocessing of melanocortin peptides, which alters the role of hormones and neuropeptides important for thermogenesis and feeding and involved in glucose homeostasis, fat oxidation, and energy expenditure [
5,
6,
7,
37].
In addition, the T rs3856806 allele of
PPARG1 was also associated with obesity in adults from northwestern Mexico. Vaisam-Castro G et al. (2021) also found an association between
PPARG rs1801282 and obesity in Brazilians [
38]. In contrast, other studies in Portuguese women [
39] and children [
40] found no association between
PPARG allelic variants and obesity. In a previous study performed by our team on adults from northwestern Mexico, we did not find an individual association with obesity, but the allelic variant analyzed was
PPARG rs1801282 [
17]. However, both
PPARG rs1801282 and rs3856806) are in linkage disequilibrium [
41]; rs1801282 is the most studied, and some work has associated this polymorphism with obesity or obesity-related traits, whereas rs3856806 has been associated with protective activity in dyslipidemia [
38,
42]. Further studies are needed to establish individual associations between
PPARG allelic variants and obesity. The mechanisms of
PPARG1 that were associated with obesity may be due to changes in energy expenditure and energy storage [
42].
On the other hand, the G rs6235 allele of
PCSK1 showed an individual association with this metabolic disorder. The association between the G rs6235 allele of
PCSK1 and metabolic syndrome might be confounded in this population because this allelic variant increased some anthropometric components of metabolic syndrome, such as BMI, weight, waist circumference, and WHR. Interestingly, no previous reports have associated the G rs6235 allele of
PCSK1 with metabolic syndrome, only with some components of metabolic syndrome. For example, in a study analyzing 27,786 European subjects,
PCSK1 rs6235 and rs6232 were associated with increased waist circumference and hip ratio in men and women [
7]. Three
PCSK1 variants (rs10515237, rs6232, rs43636321, and rs3792747) were associated with increased systolic and diastolic blood pressure and risk of hypertension in 7869 European subjects [
43]. In addition,
PCSK1 rs6232 was associated with increased fasting glucose levels [
44,
45]. In addition,
PCSK1 variants may affect HDL cholesterol; one study showed that the livers of
PCSK1 mutant mice had 2.0-fold lower levels of serum apolipoprotein A1, the primary component of HDL, but HDL cholesterol concentration was unaffected [
46]. The effect of
PCSK1 on metabolic syndrome traits could be due to (1) its activity in energy metabolism, as mentioned above, (2) regulation of blood pressure by the renin–angiotensin–aldosterone system (RAAS), as
PCSK1 is involved in the processing of prorenin to renin, and variants in this gene could affect the RAAS [
47], or (3) the involvement of
PCSK1 in glucose and lipid metabolism by modulating insulin and APOA1 production [
48,
49]. On the other hand, the lack of association of the other SNPs (
TMEM18 rs6548238 C allele,
GPX5 rs445870 G allele,
ZPR1 rs964184 G allele, and
ZBTB16 rs7106340 T allele) with metabolic diseases (obesity or metabolic syndrome) in this work could have been influenced by the sample size.
Within SNP research, there are two important approaches: (1) diagnosis, where the more SNPs associated with specific diseases (in this case, obesity and metabolic syndrome) there are in a population, the more genetic profiles can be developed to identify at an early stage subjects with a high probability of developing metabolic diseases; (2) treatment, where having a genetic profile of SNPs associated with a particular disease or trait, we can choose the best therapy to treat the metabolic disease; in this sense, SNP research is in the process of developing personalized therapies according to the genetic profile [
50].
To our knowledge, this is the first report linking the C allele of TMEM18 rs6548238 and the G allele of GPX5 rs445870 with increased systolic blood pressure and decreased HDL cholesterol, respectively, in adults, but the mechanism remains to be elucidated. In addition, this is the first report of an association between the G rs6235 allele of PCSK1 and metabolic syndrome based on international standards.
Limitations of this study included the sample size (with a larger sample, more associations are likely to be found) and the stratification of the population (the same n in each study group of normal weight versus obese and nonmetabolic syndrome versus metabolic syndrome might help to find more associations between SNPs and metabolic syndrome), as the data were not adjusted for admixture (taking into account the geographic origin based on the genetic ancestry of the study subjects might provide different associations). In addition, fewer subjects were analyzed for the G rs1801282 allele of PPARG1 compared with the other SNPs because of technical problems that prevented analysis of the entire sample.