Genetic factors associated with obesity risks in a Kazakhstani population.

OBJECTIVES: There is limited published literature on the genetic risks of chronic inflammatory related disease (eg, obesity and cardiovascular disease) among the Central Asia population. The aim is to determine potential genetic loci as risk factors for obesity for the Kazakhstani population.
SETTING: Kazakhstan. PARTICIPANTS: One hundred and sixty-three Kazakhstani nationals (ethnic groups: both Russians and Kazakhs) were recruited for the cross-sectional study. Linear regression models, adjusted for confounding factors, were used to examine the genetic associations of single nucleotide polymorphisms (SNPs) in 19 genetic loci with obesity (73 obese/overweight individuals and 90 controls).
RESULTS: Overall, logistic regression analyses revealed genotypes C/T in CRP (rs1205), A/C in AGTR1 (rs5186), A/G in CBS (rs234706), G/G in FUT2 (rs602662), A/G in PAI-1 (rs1799889), G/T (rs1801131) and A/G (rs1801133) in MTHFR genes significantly decrease risk of overweight/obesity. After stratification for ethnicity, rs234706 was significantly associated with overweight/obesity in both Russians and Kazakhs, while rs1800871 was significant in Kazakhs only.
CONCLUSIONS: This study revealed that variations in SNPs known to be associated with cardiovascular health can also contribute to the risks of developing obesity in the population of Kazakhstan. © Author(s) (or their employer(s)) 2021. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.

First to study single nucleotide polymorphisms (SNPs) previously associated with overweight/obesity in the Central Asian population.

Different association of SNPs with overweight/obesity in Kazakhs and Russians.

Haplotype analysis of SNPs.

Limited number of participants (n=163).

Participants from the southern and northern regions of Kazakhstan only.

Introduction

Overweight and obesity (OO) are important public health problems worldwide. In 2014, the worldwide proportion of overweight and obese adults (18 years and older) were 39% and 13%, respectively.1 These figures have increased twofold since the 1980s and will continue to grow.1 The prevalence of OO is increasing in the Kazakhstani population as well. Indeed, 55.6% of Kazakhstani population was overweight and 23.7% obese in 2008.2 It has been projected that prevalence of obesity will increase to 55% in 2030.2 Additionally, obesity is also an established risk factor for cardiovascular disease which is the leading cause of deaths in Kazakhstan.3 However, the roles of genetics and lifestyle factors in determining the susceptibility of obesity remains under-reported in the country. The Kazakhstani people have originated from 130 ethnic groups. In 2016, the total population was 17.8 million people, of which 66.5% were Kazakhs and 20.6% were Russians, while the remaining 12.9% included the Uzbeks, Ukrainians, Uighurs, Tatars, Germans, Koreans, Azerbaijanis, Belarusians, Turks, Dungans, Poles and others.4 Many factors have been shown to contribute to OO such as genetics, epigenetics, lack of physical exercise, comorbidities, and psychological and environmental stresses.5 There is growing evidence of the role of genetics on the development of OO. Particularly, numerous studies reported the association of single nucleotide polymorphisms (SNPs) with the risk of getting OO in different ethnic groups.6 For example, the fat mass and obesity-associated gene (FTO) is one of the first genes to be linked with obesity in European, Caucasian American, Hispanic American populations.7–10 The FTO gene polymorphisms were also found to be associated with elevated risk of obesity and obesity-related phenotypes in a North Indian population.11 Kazakhs in Xinjiang, China have the highest prevalence of metabolic syndrome compared with the whole country, namely 26.7% and 16.5%, respectively.12 13 Moreover, Kazakh G/G carriers in rs8057044 of FTO gene increased the odds of having metabolic syndrome compared with the A/A genotype, while mean body mass index (BMI) in the former carriers was significantly higher than in A/A group.14

Meanwhile, both rs9939609 in FTO and rs4994 in adrenoceptor beta 3 (ADRB3) genes were statistically associated with obesity in Russians in Moscow and Sverdlovsk cities, Russia.15 Particularly, obese participants had significantly higher prevalence of A/A or A/T genotypes in rs9939609 of FTO gene than non-obese subjects. Prevalence of rs4994 in the adrenergic gene that has been linked to obesity and metabolic disorders were also statistically different among people with and without obesity.16 Furthermore, this SNP was reported to increase the risk of myocardial infarction, stroke and heart failure in predominantly Caucasian Americans.17 SNPs can exert different effects on people from various ethnic groups. For instance, obesity associated SNPs in NEGR1 and PRKD1 genes elicited heterogeneous effects on obesity risk in people with African and European descent.18 SNP in GBE1 had different effects among Europeans and East-Asians.18 These findings suggest that a huge diversity of gene polymorphisms affect obesity risk differently in different populations. Therefore, in the Kazakhstani population, it would be useful to examine the impact of these SNPs on the risk of obesity and overweight of the different ethnic groups in the country, such as the Kazakhs and the Russians.

Currently, there is limited information focusing on the genetic determinants in various ethnic groups of Central Asia especially those on overweight/obesity and cardiovascular health. Further research in this area will help to elucidate the role of genes impacting cardiovascular health and their association with obesity. In this study, our aim was to investigate the association of various cardiovascular and endothelial health gene polymorphisms with overweight/obesity in the population living in Kazakhstan.

Materials and methods

Data collection and study population

Between June 2012 and May 2014 with the aim to study environmental, social and genetic factors determining susceptibility to tuberculosis and obesity in Kazakhstan, 1600 participants were recruited into the study.19 20 Participants were recruited from three regions namely Kyzylorda, Kostanay, Almaty oblast and Almaty city with the consideration to include sites in the northern and southern parts of the country, and the most densely populated regions. The informed consent forms were obtained by all participants, who then completed the questionnaire regarding their age, gender, socioeconomic status and habits—later regarded as potential confounders. Exposures were identified as SNPs, outcome—overweight or obesity.

Genetic samples: purification and analysis

Blood samples were collected from each participant using two vacutainer blood collection tubes (BD Vacutainer blood collection tubes) containing the anticoagulant (K2 EDTA) and polymer gel. The tubes were coded with an identifying number. Blood samples were centrifuged to separate plasma and serum, and the blood components were aliquoted into tubes and placed for −800°C storage.

DNA was extracted from the blood samples using the Wizard Genomic DNA Purification Kit (Promega, Madison, Wisconsin, USA) according to the manufacturer’s protocol. The purified DNA was quantitated, and quality ascertained using a NanoDrop spectrophotometer. For SNPs genotyping, an SNP array using quantitative real-time PCR on a Life Technologies Quantstudio 12K Real Time PCR system based on the cycle relative threshold method was employed. The results were entered and analysed with Life Technologies Copy Caller V.2.1 software. The reactions carried out were based on a two-step assay, each with two primers and a Taqman probe, one specific to the target of SNP. Overweight individuals were defined as those with BMI between 25 and 30, while those with BMI higher than 30 were considered as obese according to WHO standards.1

Statistical analysis

The statistical analysis was performed using STATA V.12.0 software (StataCorp). Student’s t-test and χ2 test were applied for bivariate analysis of data. Logistic regression was used to test whether the genotype and the allele in the SNPs were associated with OO. Factors such as age, sex, smoking status and ethnicity were measured for their possible confounding effects on the association between gene polymorphisms and OO in logistic regression. Stratification according to ethnicity (Kazakhs and Russians) was carried out to look at the effect of SNPs in each group. The statistically significant associations were reported according to α=0.05 level. Bonferroni-corrected p values were calculated for pairwise comparisons.

Results

Sociodemographic characteristics

A total of 163 people took part in the study with 73 participants in overweight/obese group and 90 participants in the control group (table 1). Mean BMI among overweight/obese people was equal to 31.81±2.97 kg/m2, while the mean for controls was 20.70±2.40 kg/m2. The mean age of the whole group was 41.47±13.07 years old with 58.02% of women and 41.98% of men. The mean age of the overweight/obese group was 45.39±12.08 years, while it was 38.33±13.05 years in controls. Additionally, proportion of males was significantly smaller (p<0.05) in overweight/obese group (35.29%) than in controls (64.71%). From the ethnicity perspective, there were in total 93 Kazakhs (57.41%) and 69 Russians (42.59%), who participated in the study. Next, smoking habits varied among study groups with higher proportion of smokers being in controls (34.44%) than in overweight/obese participants (19.44%). Finally, SNP frequencies varied significantly among overweight/obese and controls in all gene polymorphisms except those in rs1205 and rs1801133.

Table 1

Sociodemographic characteristics and frequency of studied single nucleotide polymorphisms genotypes

	Total	Overweight/obese	Control	P value
Total	163	73	90
Mean age±SD, years	41.47±13.07	45.39±12.08	38.33±13.05	0.005*
Sex, number (%)
Female	94 (58.0%)	48 (51.0%)	46 (49.0%)	0.046*
Male	68 (42.0%)	24 (35.0%)	44 (65.0%)
BMI, kg/m2	31.81±2.97	20.70±2.40	0.000*
Ethnicity, number (%)
Kazakh	93 (57.0%)	33 (35.5%)	60 (64.5%)	0.008*
Russian	69 (43.0%)	39 (56.5%)	30 (43.5%)
Smoking, number (%)
No	117 (72.2%)	58 (49.6%)	59 (50.4%)	0.034*
Yes	45 (27.8%)	14 (31.1%)	31 (68.9%)

Age and BMI are presented with mean±SD. Rest of the variables are given with number and percentage (%). Age and BMI p values were obtained from Student’s t-test with unequal variances. Other variables p values were acquired from the χ2 test.

*Statistically significant p values (<0.05).

BMI, body mass index.

Genes and SNP analysed in the study

In all, 18 genes and 19 SNPs were analysed in this study (table 2). They were categorised into four physiological subgroups relevant to cardiovascular health, namely inflammation, homocysteine and trans-sulfuration, blood pressure regulation and vascular endothelial health. These genes were selected due to previous reports suggesting their role in overweight/obesity or metabolic syndrome. Some of them include CRP (rs1205), MTHFR (rs1801133)—obesity in children or adults, AGTR1 (rs5186), CBS (rs234706), MTHFR (rs1801131, rs1801133)—metabolic syndrome.21–27 SNP in FUT2 gene (rs602662) revealed increased effect on BMI of German cohort of participants.28 The full list of tested genes and SNPs is summarised in table 2.

Table 2

Summary of single nucleotide polymorphisms tested for association with obesity and overweight, their physiological subgroup, genes they belong to and their rs numbers

Genes and physiological subgroup	Gene symbol	Rs number
Inflammation
C reactive protein	CRP	rs1205
Interleukin-10	IL-10	rs1800871
Homocysteine/trans-sulphuration
Methylenetetrahyrofolate reductase	MTHFR	rs1801133
Methylenetetrahyrofolate reductase	MTHFR	rs1801131
Cystathionine-β-synthase	CBS	rs234706
Methionine synthase	MS	rs1805087
Methionine synthase reductase	MSR	rs1801394
Catechol-O-methyltransferase	COMT	rs4680
Betaine homocysteine methyltransferase	BHMT	rs3733890
Fucosyltransferase 2	FUT2	rs602662
Transcobalamine II	TCN2	rs1801198
Blood pressure regulation
Angiotensinogen	AGT	rs699
ACE	ACE	rs4340
Angiotensin II receptor-1	AGTR1	rs5186
Vascular endothelial health
Endothelial nitric oxide synthase	NOS	rs1799983
NADP-oxidase p22 phox	CYBA	rs4673
NADPH-oxidase-encoded by CYBA	CYBA	rs9932581
Plasminogen activator inhibitor-1	PAI-1	rs1799889
Adiponectin	ADIPOQ	rs1501299

Association of SNPs with overweight/obesity in Kazakhstani population

Among the 19 SNPs studied, 8 of them showed statistically significant differences between experimental and control groups. These are presented in table 3A, B. Particularly, genotypes under the category inflammation (rs1205), homocysteine and trans-sulphuration (rs1801131, rs1801133, rs234706, rs602662), blood pressure regulation (rs5186) and vascular endothelial health (rs1799889) SNPs. These SNPs had decreased odds of having overweight/obesity. Most of the significant genotypes were heterozygotes except GG in rs602662. Furthermore, allelic analysis for risk of being overweight/obese revealed a statistical difference only in rs1799889. The G allele decreased odds of overweight/obesity by 0.54 times (95% CI 0.34 to 0.84) compared with A allele carriers. None of the tested potential confounders, age, sex, smoking status and ethnicity, were revealed to have any statistical significance in the logistic regression model.

Association of SNPs with overweight/obesity in Kazakhs and Russians

Analysis of genotype and allelic effect on the risk of obesity/overweight in Kazakhs and Russians separately revealed the following results (table 4). After stratification, rs1205, rs1801131 and rs1801133 lost their statistical significance in Kazakhs and Russians. Meanwhile, two SNPs, rs234706 and rs602662, were still statistically significant in both ethnic groups. AG genotype in rs234706 had a similar effect on overweight/obesity risk among Russians and Kazakhs, namely ORs were equal to 0.128 (95% CI 0.02 to 0.78) and 0.197 (95% CI 0.05 to 0.78), respectively. Furthermore, G allele in rs602662 decreased the probability of being overweight/obese by 58% in Kazakhs (95% CI 0.22 to 0.80) and by 50% in Russians (95% CI 0.05 to 0.78). Effect on either one of two ethnicities was found in rs5186 (Kazakhs), rs1799889 (Russians) and rs1800871 (Kazakhs). Genotype AC in rs5186 decreased the risk of obesity and overweight by 0.26 times (95% CI 0.09 to 0.78) when compared with AA genotype carriers, while AG genotype in rs1800871 lowered it by 70% (95% CI 0.09 to 0.93) compared with the reference (AA) group. Meanwhile, analysis of rs1799889 revealed statistical significance of two genotypes in Russians, AG and GG. Individuals with the former genotype had 0.27-fold (95% CI 0.07 to 0.97) decreased risk of being obese and overweight than AA carriers, while OR value in those with the latter genotype was smaller, 0.21 (95% CI 0.05 to 0.91). Analysis of frequency distribution of genotypes alleles in two ethnicities revealed statistical difference in rs1205 (both), rs5186 (allele), rs234706 (both), rs602662 (allele) and rs1800871 (both).

Table 4

The genotype and allele associations with overweight/obesity in Kazakhs and Russians and the frequencies of the genetic variation

SNP (gene)	Genotype/allele	Kazakhs		Russians		Kazakhs	Russians	P value
		OR (95% CI)	P value	OR (95% CI)	P value	n (%)	n (%)
rs1205 (CRP)	Genotype
	CC	1		1		25 (27.1)	35 (50.7)	0.020*
	CT	0.46 (0.17 to 1.26)	0.132	0.67 (0.25 to 1.78)	0.420	49 (53.3)	30 (43.5)
	TT	1.02 (0.30 to 3.45)	0.977	1.99 (0.19 to 21.23)	0.565	18 (19.6)	4 (5.8)
	Allele
	C	1		1		99 (53.8)	100 (72.5)	0.001*
	T	0.95 (0.52 to 1.74)	0.861	0.93 (0.44 to 1.98)	0.854	85 (46.2)	38 (27.5)
rs5186 (AGTR1)	Genotype
	AA	1		1		60 (65.2)	34 (49.3)	0.096
	AC	0.26 (0.09 to 0.78)	0.016*	0.51 (0.19 to 1.37)	0.182	30 (32.6)	31 (44.9)
	CC	1.31 (0.08 to 21.9)	0.852	Absent		2 (2.2)	4 (5.8)
	Allele
	A	1		1		150 (81.5)	99 (71.7)	0.038*
	C	0.42 (0.17 to 1.03)	0.059	0.99 (0.47 to 2.10)	0.987	34 (18.5)	39 (28.3)
rs234706 (CBS)	Genotype
	AA	1		1		6 (6.6)	16 (23.2)	0.000*
	AG	0.128 (0.02 to 0.78)	0.026*	0.197 (0.05 to 0.78)	0.020*	59 (64.8)	50 (72.5)
	GG	0.80 (0.12 to 5.20)	0.815	Absent		26 (28.6)	3 (4.3)
	Allele
	A	1		1		71 (39.1)	82 (59.4)	0.000*
	G	1.67 (0.88 to 3.18)	0.115	0.72 (0.36 to 1.44)	0.354	111 (60.9)	56 (40.6)
rs602662 (FUT2)	Genotype
	AA	1		1		12 (13.1)	17 (24.6)	0.062
	AG	0.51 (0.14 to 1.92)	0.320	1.09 (0.31 to 3.81)	0.892	36 (39.1)	30 (43.5)
	GG	0.21 (0.05 to 0.81)	0.023*	0.31 (0.08 to 1.17)	0.084	44 (47.8)	22 (31.9)
	Allele
	A	1		1		60 (32.6)	64 (46.4)	0.012*
	G	0.42 (0.22 to 0.80)	0.008*	0.50 (0.25 to 0.99)	0.046*	124 (67.4)	74 (53.6)
rs1799889 (PAI-1)	Genotype
	AA	1		1		21 (22.8)	19 (27.5)	0.714
	AG	0.6 (0.21 to 1.68)	0.332	0.27 (0.07 to 0.97)	0.045*	51 (55.4)	34 (49.3)
	GG	0.28 (0.07 to 1.10)	0.069	0.21 (0.05 to 0.91)	0.037*	20 (21.8)	16 (23.2)
	Allele
	A	1		1		93 (50.5)	72 (52.2)	0.772
	G	0.58 (0.31 to 1.07)	0.081	0.47 (0.24 to 0.93)	0.031*	91 (49.5)	66 (47.8)
rs1800871 (IL-10)	Genotype
	AA	1		1		20 (22.2)	7 (10.4)	0.018*
	AG	0.30 (0.09 to 0.93)	0.038*	0.17 (0.02 to 1.58)	0.119	43 (47.8)	26 (38.8)
	GG	0.69 (0.21 to 2.20)	0.529	0.21 (0.02 to 1.95)	0.170	27 (30)	34 (50.8)
	Allele
	A	1		1		83 (46.1)	40 (29.9)	0.004*
	G	0.87 (0.47 to 1.61)	0.657	0.71 (0.33 to 1.52)	0.379	97 (53.9)	94 (70.1)
rs1801131 (MTHFR)	Genotype
	GG	1		1		9 (9.9)	13 (19.1)	0.083
	GT	0.33 (0.08 to 1.41)	0.132	0.29 (0.07 to 1.21)	0.089	45 (49.5)	23 (33.8)
	TT	0.43 (0.09 to 1.89)	0.267	0.74 (0.19 to 2.94)	0.670	37 (40.6)	32 (47.1)
	Allele
	G	1		1		63 (34.6)	49 (36.0)	0.794
	T	0.85 (0.44 to 1.61)	0.613	1.05 (0.52 to 2.13)	0.891	119 (65.4)	87 (64.0)
rs1801133 (MTHFR)	Genotype
	AA	1		Absent		9 (9.7)	0	0.122
	AG	0.33 (0.07 to 1.60)	0.168	0.64 (0.24 to 1.72)	0.375	29 (31.2)	35 (51.5)
	GG	0.89 (0.22 to 3.71)	0.882	1		55 (59.1)	33 (48.5)
	Allele
	A	Reference		Reference		47 (25.3)	42 (30.9)	0.266
	G	1.41 (0.69 to 2.88)	0.346	0.81 (0.39 to 1.69)	0.568	139 (74.7)	94 (69.1)

Genotypes GG (rs1801133), GG (rs234706) and CC (rs5186) were absent in Russian group. GG genotype was chosen as the reference group for rs1801133 in Russians.

*Statistically significant p values (<0.05).

SNP, single nucleotide polymorphism.

Haplotype analysis of rs1801131 and rs1801133 found on MTHFR gene and overweight/obesity

Table 5 depicts haplotype analysis of rs1801131 and rs1801133 located on the MTHFR gene. Combination of T on rs1801131 and G rs1801133 had a 0.27-times (95% CI 0.11–0.70, p value, 0.007) protective effect against obesity and overweight than the reference G and A alleles on respective SNPs.

Table 5

Haplotype analysis of rs1801131 and rs1801133 found on MTHFR gene with association with overweight/obesity

	Rs1801131	Rs1801133	Overweight/obese, n (%)	Control, n (%)	OR (95% CI)	P value
Haplotype 1	G	A	8 (26.7)	22 (73.3)	Reference
Haplotype 2	T	G	46 (38.7)	73 (61.3)	0.27 (0.11 to 0.70)	0.007*

*Statistically significant p values (<0.05).

Discussion

Findings from this study suggested that SNPs which are known for their role in cardiovascular health are also associated with risk of overweight/obesity as observed in a Kazakhstani population. Particularly, those SNPs belonging to genes in the homocysteine/trans-sulphuration pathway. All of the SNPs frequencies were statistically different when compared between controls and overweight/obese people except rs1205 of the CRP gene and rs1801133 of the MTHFR gene. Further logistic regression analysis of association of genetic polymorphisms with overweight/obesity showed that heterozygotes had decreased risk for OO including rs1205 in CRP, rs5186 in AGTR1, rs234706 in CBS, rs1800871 in interleukin 10 (IL-10), and both rs1801131 and rs1801133 in MTHFR. Meanwhile, homozygote GG in rs602662 in FUT2 and both homozygote GG and heterozygote in rs1799889 of the PAI-1 gene, were shown to be protective against OO. Moreover, after adjusting for age, sex and ethnicity, all SNPs retained their statistical significance except for rs1205, rs1800871 and rs1801133. In contrast, adjustment resulted in a new statistical significance for the G allele in SNP rs602662 of the FUT2 gene. Furthermore, when stratification analysis for ethnicity was carried out it showed that certain SNPs can have statistical significance on the risk of OO in Kazakhs only (rs5186, rs1800871), in Russians only (rs1799889) and in both ethnic groups (rs234706 and rs602662). Finally, haplotype analysis of MTHFR gene polymorphisms revealed a novel finding on having T allele in rs1801131 and G in rs1801133 can protect their carriers from OO.

FUT2 protein is crucial for synthesis of H antigen that is needed for the A and B antigens production pathway. FUT2 has been found to have a role in the composition of gut microbiota, which, in turn, can affect the risk of obesity due to inadequate interaction between microbes and human organism.29–31 This supports our finding that there is a protective association of homozygotes (GG), which is also known as 5G/5G genotype, against OO providing sufficient representation of H antigens in gut epithelium. The recent study in German and Danish populations also revealed that G allele had a positive association with BMI, though our study’s analysis did not show any statistical significance for the G allele.28

Though the exact role of PAI-1 in obesity is still unclear, numerous studies show that adipose tissue produces higher amounts of PAI-1 in obese people than in controls, while knockdown of PAI-1 in mice resulted in fat loss; suggesting a role of PAI-1 in fat mass control.32–36 Research in Scandinavian countries found a statistically significant difference in the distribution of heterozygote (AG) among obese and control patients, which supports the results of this study.37 Moreover, AG and AA genotypes are associated with increased obesity in Turkish population, which shares a common ancestry with Kazakhs.38 Association from a different perspective was shown in Mexican study, which revealed increased risk for OO for patients with haplogenotype of several SNPs, namely rs1799889 (G)+rs2227631 (A)+rs757716 (C).39 However, Argentinean researchers did not find a statistically significant association of this SNP with obesity.40 Additionally, study in Mexico did not find any significant results in terms of PAI-1 SNP rs1799889 and the risk of being obese.41

IL-10 is an important anti-inflammatory cytokine, and low levels of IL-10 are associated with obesity, and thus leading to impaired immune system and other consequences including cardiovascular diseases.42 43 In our study AA genotype in rs1800871 of IL-10 had the highest risk of getting obese in Kazakhs only. Meanwhile, research of Indian patients with metabolic syndrome and high BMI, showed increased risk of obesity for the A allele carriers, though our study did not find statistical difference for this allele.44 We hypothesise that patients with the AA genotype could have the lowest levels of IL-10, and hence the increased risk of being OO.

MTHFR participates in the methionine to cysteine conversion process, and it is believed that rs1801131 and rs1801133 polymorphisms can cause overexpression of leptin gene, and thus leading to growth of adipose tissue.45 In this study, the GT genotype in SNP rs1801133 of MTHFR showed protective features against overweight/obesity. Meanwhile, similar research using logistic regression revealed positive results only in study of Chinese patients with type 2 diabetes and AA genotype.46 In contrast, studies in Lithuania, China, Italy and Ecuador, which also used logistic regression, did not find any statistically significant link between this SNP and overweight/obesity.45 47–49 In the studies among Chinese children, A allele was found to have association with obesity, though in our study it did not reach significance possibly due to smaller sample size.22 27

SNP rs1801131 on the same gene, MTHFR, has been studied less extensively, positive association of this SNP from our research was not supported by studies carried out on Chinese, Ecuadorian and European populations.48–50 However, results from an Italian research, that studied both rs1801131 and rs1801133, found increased risk for carriers with the AC genotype in rs1801131 for having an overweight/obesity trait.45 As discussed earlier, GT genotype was also implied to play a role in metabolic syndrome (hypertension) among Chinese elderly.26 However, this finding was not observed in the longitudinal study and in our study.

In this study CBS SNP rs234706’s is found to have a significant association with obesity, CBS like MTHFR, also participates in the pathway of conversion of methionine to cysteine, was not confirmed to be significantly associated with obesity by neither the large European study in nine countries (Germany, Spain, Italy, Hungary, Austria, Sweden, France, Greece and the UK) nor the study carried out in India.50 51 In the GWAS study A allele of rs234706 was found to have significant dysregulation of CBS, suggesting its role in impaired lipid metabolism.25 Though in our study A allele did not reach statistical significance, only AG genotype. Furthermore, there is insufficient number of studies on humans to elucidate the role of CBS in obesity risk. However, the study in rodents suggested a role for CBS in cardiac lipotoxicity, in which deficiency of at least one CBS gene (CBS +/−) was associated with increased fatty acid uptake, heart weight and fasting insulin levels compared with the wild type mice (CBS +/+).52

The decreased risk for OO that we obtained for CRP SNP rs1205, however, is not supported by studies in Brazilians and Mexican-Americans.53 54 CRP is an important inflammatory marker, often induced by adipose tissue via proinflammatory cytokines (Tumour Necrosis Factor -α, IL-6).55 56 It has been reported to be associated with various cardiovascular diseases such as myocardial infarction, stroke, peripheral vascular disease.57 58 Also, elevated CRP levels have been linked with increased risk of obesity.55–59 In the Brazilian study of children and adolescents T allele had increased odds for hypercholesterolaemia, although in our study it did not get statistical significance.21 Furthermore, it was found that CC genotype in rs1205 of CRP results in high levels of the corresponding protein in Taiwanese, and our study provides possible additional evidence for this association.60

AGTR1 protein takes part in blood pressure regulation, namely the renin angiotensin system. Though there is limited information on the exact pathway in which it participate in promoting obesity, there are several speculations regarding its role, including the expression of AGTR1 in adipose tissue, its participation in the production of various proteins and molecules by adipocytes (prostacyclin, norepinephrine, nitric oxide).61–63 As a result, mature adipocytes release cytokines, which lead to inflammation, vessel damage, hypertension.64 These studies provide support and explanation for the association of rs1205 in AGTR1 and OO found in our study. However, while our results show that carriers with AC in AGTR1 SNP rs5186 had protective effect against OO, in Romanian population carriers with either AC or CC genotype had the opposite effect, particularly it was associated with an increased risk of obesity.65 However, Egyptian and Tunisian studies revealed no statistically significant results.66 67

Though statistically significant differences in both allele and genotype distribution among Kazakh and Russians were found in CRP, CBS and IL-10 gene polymorphisms, stratification analysis revealed ethnic-specific associations regarding their risk for OO. In particular, rs234706 (CBS) should be considered a potential risk factor for OO in both ethnicities, while rs1800871 (IL-10) is statistically significant for Kazakhs but not for Russians.

Our analysis indicated that carriers of T in rs1801131 and G in rs1801133 in MTHFR had decreased OR for being OO. And this is an important novel finding as no research reports studying the influence of rs1801131 and rs1801133 haplotype on obesity risk was found up to February 2019. However, there are a number of studies showing the association of just MTHFR rs1801133, C677T with obesity. It will be interesting and important for further studies with bigger sample size to confirm the association of the MTHFR haplotype to OO in future.

This study has several major limitations. First, the total number of participants studied is small (n=163), and they were predominantly from one region of Kazakhstan, Kostanay city. Therefore, other regions remain under-represented. Second, there is a lack of measured clinical parameters such as blood pressure, waist and hip circumference, and other comorbidities including hypertension, metabolic syndrome. Third, the lack of similar studies from the Central Asian regions do not allow proper comparison of these new findings. Therefore, further research in these populations would be useful for elucidating the role of these SNPs and the risk of overweight/obesity. Finally, multiple SNPs interaction, which might also influence risk of OO, was not analysed in this study and this will be a useful step in our future research.

Conclusions

This study found associations of SNPs involved in cardiovascular health with overweight/obesity trait in Kazakhstani population. Genetic polymorphisms from the two major ethnicities, Kazakhs and Russians, were analysed, and associations in one or both were found. These findings will be helpful in predicting genetic risk for overweight/obesity and provide further understanding and guide treatment strategies for this condition in the country.

Table 3A

List of SNPs, genotypes and frequency

SNP (gene)	Total (%)	Overweight/obese (%)	Control (%)	P value
rs1205 (CRP)
CC	60 (37.2)	32 (53.3)	28 (46.7)	0.091
CT	79 (49.1)	28 (35.4)	51 (64.6)
TT	22 (13.7)	11 (50.0)	11 (50.0)
rs5186 (AGTR1)
AA	94 (58.4)	47 (50.0)	47 (50.0)	0.010*
AC	61 (37.9)	19 (31.1)	42 (68.9)
CC	6 (3.7)	5 (83.3)	1 (16.7)
rs234706 (CBS)
AA	22 (13.8)	17 (77.3)	5 (22.7)	0.000*
AG	109 (68.1)	35 (32.1)	74 (67.9)
GG	29 (18.1)	19 (65.5)	10 (34.5)
rs602662 (FUT2)
AA	29 (18.0)	18 (62.1)	11 (37.9)	0.001*
AG	66 (41.0)	35 (53.0)	31 (47.0)
GG	66 (41.0)	18 (27.3)	48 (72.7)
rs1799889 (PAI-1)
AA	40 (24.8)	25 (62.5)	15 (37.5)	0.015*
AG	85 (52.8)	35 (41.2)	50 (58.8)
GG	36 (22.4)	11 (30.6)	25 (69.4)
rs1800871 (IL-10)
AA	27 (17.2)	16 (59.3)	11 (40.7)	0.041*
AG	69 (43.9)	23 (33.3)	46 (66.7)
GG	61 (38.9)	30 (49.2)	31 (50.8)
rs1801131 (MTHFR)
GG	22 (13.9)	14 (63.6)	8 (36.4)	0.022*
GT	68 (42.7)	22 (32.4)	46 (67.6)
TT	69 (43.4)	33 (47.8)	36 (52.2)
rs1801133 (MTHFR)
AA	14 (8.7)	9 (64.3)	5 (35.7)	0.079
AG	61 (37.9)	21 (34.4)	40 (65.6)
GG	86 (53.4)	41 (47.7)	45 (52.3)

0.000* indicates p value less than 0.0005.

*Statistically significant p values (<0.05).

SNP, single nucleotide polymorphism.

Table 3B

The genotype and allele associations with overweight/obesity

SNP (gene)	Genotype/allele	ORunadjusted (95% CI)	P value	ORadjusted (95% CI)	P value
rs1205 (CRP)	Genotype
	CC	1		1
	CT	0.48 (0.24 to 0.95)	0.036*	0.51 (0.24 to 1.07)	0.156
	TT	0.88 (0.33 to 0.44)	0.789	0.96 (0.32 to 2.85)	1.000
	Allele
	C	1		1
	T	0.80 (0.51 to 1.26)	0.327	0.83 (0.51 to 1.37)	0.488
rs5186 (AGTR1)	Genotype
	AA	1		1
	AC	0.45 (0.23 to 0.89)	0.021*	0.34 (0.15 to 0.72)	0.010*
	CC	4.99 (0.56 to 44.44)	0.149	4.59 (0.41 to 51.35)	0.430
	Allele
	A	1		1
	C	0.79 (0.47 to 1.35)	0.393	0.64 (0.36 to 1.14)	0.134
rs234706 (CBS)	Genotype
	AA	1		1
	AG	0.14 (0.05 to 0.41)	0.000*	0.21 (0.06 to 0.66)	0.016*
	GG	0.56 (0.16 to 1.96)	0.364	1.31 (0.32 to 5.39)	1.000
	Allele
	A	1		1
	G	0.95 (0.61 to 1.47)	0.803	1.20 (0.74 to 1.94)	0.454
rs602662 (FUT2)	Genotype
	AA	1		1
	AG	0.69 (0.28 to 1.68)	0.415	0.73 (0.28 to 1.89)	1.000
	GG	0.23 (0.09 to 0.58)	0.002*	0.26 (0.09 to 0.71)	0.018*
	Allele
	A	1		1
	G	0.93 (0.53 to 1.63)	0.796	0.45 (0.27 to 0.74)	0.002*
rs1799889 (PAI-1)	Genotype
	AA	1		1
	AG	0.42 (0.19 to 0.91)	0.028*	0.44 (0.19 to 1.01)	0.112
	GG	0.26 (0.10 to 0.69)	0.006*	0.23 (0.08 to 0.66)	0.012*
	Allele
	A	1		1
	G	0.54 (0.34 to 0.84)	0.006*	0.51 (0.31 to 0.83)	0.007*
rs1800871 (IL-10)	Genotype
	AA	1		1
	AG	0.34 (0.14 to 0.86)	0.022*	0.34 (0.12 to 0.94)	0.076
	GG	0.67 (0.27 to 1.66)	0.384	0.66 (0.24 to 1.85)	0.876
	Allele
	A	1		1
	G	0.95 (0.60 to 1.50)	0.826	0.95 (0.58 to 1.57)	0.860
rs1801131 (MTHFR)	Genotype
	GG	1		1
	GT	0.27 (0.09 to 0.75)	0.012*	0.26 (0.08 to 0.79)	0.034*
	TT	0.52 (0.19 to 1.41)	0.200	0.50 (0.17 to 1.45)	0.408
	Allele
	G	1		1
	T	0.93 (0.58 to 1.47)	0.741	0.91 (0.55 to 1.49)	0.715
rs1801133 (MTHFR)	Genotype
	AA	1		1
	AG	0.29 (0.09 to 0.98)	0.047*	0.29 (0.08 to 1.10)	0.070
	GG	0.51 (0.16 to 1.63)	0.255	0.66 (0.18 to 2.36)	0.533
	Allele
	A	1		1
	G	1.02 (0.62 to 1.66)	0.060	1.21 (0.71 to 2.06)	0.469

*Statistically significant p values (<0.05).

SNP, single nucleotide polymorphism.