- Research Division, Dasman Diabetes Institute, Dasman, Kuwait
Association studies have implicated single nucleotide polymorphisms (SNPs), particularly rs1421085, from the fat mass and obesity-associated (FTO) gene with body composition phenotypes, obesity, dietary intake, and physical activity in European, East Asian, and African populations. However, the impact of the rs1421085 variant has not been sufficiently tested in ethnic populations (such as Arabs) with high levels of obesity. Further, there is a lack of studies identifying biomarkers that interact with FTO. Therefore, we investigated the association of rs1421085 with obesity and body composition traits and metabolic biomarkers in Arab population. We genotyped rs1421085 SNP in 278 Arab individuals, where multiple biomarkers relating to obesity, inflammation, and other metabolic pathways were quantified. We performed genetic association tests under additive mode of inheritance using linear regression models and found association of rs1421085_C allele with higher levels of body weight, soft lean mass (SLM), and total body water. Examination (using linear regression models under dominant mode of inheritance) of correlation among biomarkers and interaction with genotypes at the variant revealed that measures of these three body composition traits were found mediated by interaction between carrier genotypes (TC+CC) and measures of ghrelin, ApoA1, and ApoB48. Lean body mass (LBM), to which SLM contributes, is an important determinant of physical strength and is a focal point in studies on sarcopenia. Low LBM is known to be associated with higher risk of cardiometabolic disorders. Thus, the finding on the FTO variant as a genetic determinant of SLM via interaction with ghrelin, ApoA1, and ApoB48 is important.
Introduction
The fat mass and obesity-associated (FTO) gene is a well-established obesity-susceptibility locus. Early genome-wide association studies (GWAS) associated FTO genetic variants with obesity traits in Europeans (Dina et al., 2007; Frayling et al., 2007); subsequently, many of these risk variants were replicated in Africans, South and East Asians, and other populations. Furthermore, association of FTO variants with type 2 diabetes mellitus (T2DM), independent of body mass index (BMI), has been demonstrated in East and South Asians (Li et al., 2012) and Scandinavians (Hertel et al., 2011). It is also the case that FTO is associated with body composition traits such as lean body mass, fat free mass, and body fat percentage in people of European ancestry (Zillikens et al., 2017; Hubel et al., 2019).
FTO has been functionally implicated in both adipogenesis and tumorigenesis and on the underlying m6A-dependent demethylase activity which influences several mRNA processing events (Gerken et al., 2007; Bartosovic et al., 2017; Deng et al., 2018). The most notable genetic variants from FTO associated with obesity and body composition traits are a cluster of at least 20 common SNPs from its first intron (Loos and Yeo, 2014). Despite the association of many FTO genetic variants with obesity traits in humans, elucidation of the mechanistic basis for these associations has been elusive. Obesity has reached epidemic proportions, and understanding its genetic basis will help mitigate obesity-associated illness. In this context, there is an urgent need to decipher the causal mechanism for the role of FTO noncoding variants in the aetiology of obesity and its related disorders (Wu and Arora, 2016).
Among the published FTO variants, rs1421085_T > C is of particular interest; previously published global GWAS [as listed in NHGRI-EBI GWAS Catalog (Buniello et al., 2019)] associated the rs1421085 variant with risk for obesity (European) (Wheeler et al., 2013), increased BMI (European) (Felix et al., 2016), risk for T2DM (Japanese, European) (Imamura et al., 2016; Kichaev et al., 2019), increased dietary macronutrient intake (European) (Tanaka et al., 2013), increased body fat percentage (European), decreased hand grip strength (British) (Tikkanen et al., 2018), increased fat-free mass (European) (Hubel et al., 2019), decreased pulse pressure (East Asian) (Takeuchi et al., 2018), and preference for morning sleep (European) (Kichaev et al., 2019). A recent study by Claussnitzer et al. (2015) ascribed a critical role for the rs1421085 variant, in adipocyte-thermogenesis regulation pathway through interaction with ARID5B, IRX3, and IRX5 genes, by utilising an integrative approach combining data on the regulatory role of the variant on its upstream and downstream gene regulators that were identified from genetic and epigenetic profiles of animal models. It therefore indicates that manipulation of this pathway may regulate obesity levels. Similar approaches are needed to delineate the influence of genetic variants on the aetiology of obesity, and this may be achieved through deep phenotyping of biomarkers, known to be involved in related biological processes, and by studying these genetic variants in various cohorts (Price et al., 2015; Szendroedi et al., 2016).
According to the most recent World Health Organization report, worldwide prevalence of obesity has tripled since 1975 (https://www.who.int/news-room/fact-sheets/detail/obesity-and-overweight). The prevalence of obesity is most noticeable in the Gulf region, where nine countries ranked highest in the obesity statistics among adults aged ≥18 years. Despite its high prevalence in the Gulf region, the Arab population is under-represented in global studies of obesity. In this study, we evaluated the influence of the FTO rs1421085 variant on various anthropometric and body composition measures and obesity, inflammatory, and cardiovascular biomarkers in a cohort of Kuwait Arab subjects.
Materials and Methods
Study Participants
The study protocol was reviewed and approved by the Ethical Review Committee of Dasman Diabetes Institute and was conducted according to the guidelines of the Declaration of Helsinki and of the US Federal Policy for the Protection of Human Subjects. Participants were Kuwaitis of Arab ethnicity, which was confirmed through detailed questioning on parental lineage up to three generations. The study population consisted of 278 subjects whose data on age, sex, comorbidities (e.g., diabetes and cardiovascular complications) and baseline characteristics, and vital signs such as height, weight, and waist circumference (WC) were recorded during enrolment. Details on medications taken by the participants for lowering lipid levels and diabetes were collected and used in correction procedures with the association statistics. All participants signed the consent form before participating in the study.
Sample Collection and Processing
Collection of blood samples and measurement of vital signs were performed in accordance with established institutional guidelines. DNA was extracted using the Gentra Puregene® kit (Qiagen, Valencia, CA, USA). Quant-iT™ PicoGreen® dsDNA Assay Kit (Life Technologies, Grand Island, NY, USA) and Epoch Microplate Spectrophotometer (BioTek Instruments) were used to quantify DNA; absorbance values at 260–280 nm were checked for adherence to an optical density range of 1.8–2.1.
Estimation of Plasma Levels of Various Biomarkers
The plasma levels of 46 biomarkers listed in Supplementary Table S1 were determined. Briefly, blood samples were collected in EDTA-treated tubes. Plasma, separated by centrifugation, was aliquoted and stored at −80°C. The levels of apolipoproteins, except for ApoB48, were determined using Luminex assay from Millipore. The levels of inflammatory markers were determined using Luminex immunobead array (R&D Systems). The data were processed using the Bio-Plex Manager Software version 6 (Bio-Rad) and analysed using five-parametric curve fitting. High Sensitivity C-Reactive Protein (HsCRP) level was measured by enzyme-linked immunosorbent assay, as reported earlier (Abu-Farha et al., 2014; Abu-Farha et al., 2016) and adipokines were measured using Luminex assay (R&D Systems). Analytes provided by the manufacturer were used as quality controls to validate specificity and sensitivity.
Bioelectric Impedance Measurements
Body composition of the participants was measured using Bio-impedancemetry Body Composition Analyzer IOI 353 (Jawon Medical Co., South Korea). Total body water (TBW) is the summation of body water that exists inside the cell membrane (intracellular) and water outside the cell membrane (i.e., blood and intestinal fluid, extracellular). Soft lean mass (SLM) is the addition of TBW and proteins in the body and is made up of skeletal and smooth muscle (Bera, 2014). Lean body mass (LBM) is the summation of SLM and minerals. Percentage body fat (PBF) was estimated by subtracting LBM from the total body weight.
Targeted Genotyping
The TaqMan® Genotyping Assay on ABI 7500 Real-Time PCR System from Applied Biosystems (Foster City, CA, USA) was used to perform candidate SNP genotyping. Next, 10 ng of DNA, 5× FIREPol® Master Mix (Solis BioDyne, Estonia), and 1 µl of 20× TaqMan® SNP Genotyping Assay constituted each polymerase chain reaction sample, and thermal cycling conditions were set at 60°C for 1 min and 95°C for 15 min followed by 40 cycles of 95°C for 15 s and 60°C for 1 min. Sanger sequencing, using the BigDye™ Terminator v3.1 Cycle Sequencing on an Applied Biosystems 3730xl DNA Analyzer, was performed for selected cases of homozygotes and heterozygotes to validate genotypes determined by the above techniques.
Power Calculation
We used the Quanto software tool (University of Southern California, Los Angeles, CA, USA) to calculate the power of the study cohort and its ability to delineate quantitative trait variability at a given power (which we set at 80%). We considered the underlying genetic model to be based on additive mode of inheritance; we performed the power calculation considering “gene only” contributions with the underlying genetic model as based on additive mode of inheritance (AA versus AB versus BB) or under “gene-environment (GxE)” contributions (Parnell et al., 2014) with the underlying genetic model as based on dominant mode of inheritance (AA+AB versus BB, where A is the effect allele; see (Zhao et al., 2016) for definitions of genetic models for complex disorders). We set a P-value threshold of ≤ 0.05 for considering a type 1 error as significant. Genetic effect that accounts for at least 0.1%–5% variance in the trait was detected by letting the RG2 (estimate for marginal genetic effect) to assume values in the range of 0.001–0.05 in step of 0.005. We considered the population (mean ± standard deviation) of the quantitative trait and of the environment term (interaction trait) in these calculations.
Quality Assessment of SNP and Trait Measurements and Statistical Association Analysis
We used PLINK (version1.9) to assess quality and statistical association of the rs1421085 SNP. We calculated minor allele frequency (MAF) and Hardy–Weinberg equilibrium for the FTO variant. We considered any quantitative trait values < Q1- 1.5 * IQR or any value > Q3+1.5 * IQR as an outlier and excluded it from further statistical analysis.
A total of 60 traits and biomarkers (Supplementary Table S1) were considered in this study; it is to be noted that there can be dependency among some of these traits. Independent traits were identified from the summary statistics of Pearson correlation analysis among the 60 traits using the matSpD tool (Li and Ji, 2005) (available at https://sites.google.com/site/qutsgel/software/matspd-local-version). The analysis pointed to 23 independent traits.
We performed allele-based statistical association tests for rs1421085 with all 60 traits using linear regression adjusting for regular corrections toward age and sex. We also adjusted for diabetes medication and lipid-lowering medication. Correction for multiple testing was assessed by adjusting the P-value threshold for the 23 independent traits (0.05/effective number of independent traits), which was (0.05/23 = 0.0021). We also examined empirical P-values (Pemp-value), generated using the max(T) permutation procedure available in PLINK, based on 10,000 permutations.
Further, correlations between the rs1421085 variant and levels of biomarkers/traits (such as SLM, TBW and weight), corrected for age/sex, as well as interactions of such correlations with other traits (example: SLM ~ genotypes+Age+sex+ApoA1*ApoB48) were evaluated using multivariable linear regression tests. Summary statistics on the relationships among the traits/biomarkers included percentage variation in response variable (r2), estimates denoting standardised beta-coefficients, standard error, and significance of test (P) for the genotype distributions. P-value threshold was set at 0.05.
Results
Sample Set for Each of the Traits
The quality control procedures for checking outliers in trait measurements and missing information on traits led to varying number of samples for different traits, with the maximum size of 270 for BMI.
Characteristics of the Study Cohort
The summary statistics on the considered 60 traits and biomarkers of the 278 participants are summarised in Supplementary Table S1. The mean age of the cohort was 46.26 ± 12.39 years, with the ratio of males to females as 1:1.2. The cohort comprised mostly class I obese subjects with a mean BMI of 29.93 ± 5.18 kg/m2 and mean WC (99.36 ± 13.36 cm); 48% of the subjects were obese and 41% were diabetic. Mean HbA1c (6.31 ± 1.3%), low-density lipoprotein (3.13 ± 0.97 mmol/L), high-density lipoprotein (HDL) (1.2 ± 0.32 mmol/L), total cholesterol (5.02 ± 1.09 mmol/L), and triglyceride (1.22 ± 0.6 mmol/L) were normal or near optimal. Of the 278 subjects, 104 were taking medication for diabetes and 91 for lowering lipids.
After partitioning based on the genotypes (TT vs TC+CC) at the FTO rs1421085 variant, the cohort exhibited significant differences (P-value ≤ 0.05) in the mean values for weight, BMI, SLM, TBW, leptin, eotaxin, TNFα, MMP-9, MIP1α, IFNγ, and several interleukins (see Supplementary Table S1). The significances of differences were particularly strong for weight (P-value = 0.0047), IL-5 (P-value = 0.00076), IL-6 (P-value = 0.0046), IL-7 (P-value = 0.0043), IFNγ (P-value = 0.0078), and TNFα (P-value = 0.0052). Upon performing the matSpD analysis on these interleukins that showed significant differences, it was found that the effective number of independent variables were 3, suggesting that they were partially correlated.
Association Between FTO rs1421085 Variant and Increased SLM, TBW, and Weight
The summary statistics of associations observed between the FTO variant and the phenotype traits, derived from allele-based association tests based on additive models corrected for age and sex, are listed in Table 1 and Supplementary Table S2. Statistically significant associations (P-value threshold corrected for multiple testing ≤ 0.002) between the variant (with the minor allele C as the effect allele) and three traits (SLM, TBW, and weight) were observed; and values of Pemp-value ≤ 0.05 were observed for these three association signals when the tests were further corrected for medication (Table 1). Though association signals for five other traits were at borderline P-values (C-peptide at P-value = 0.004, IFNγ at P-value = 0.005, eotaxin at P-value = 0.007, TNFα at P-value = 0.007, and IL-5 at P-value = 0.008), Pemp-value did not reach the desired value of ≤ 0.05 even after the tests were corrected for medication (see Supplementary Table S2).
Table 1 Significant results (P-value < 0.002) of association tests for the FTO variant rs1421085 (with C as the effect allele) with phenotype traits, using genetic model based on additive mode of inheritance (TT versus TC versus CC).
Differences in the Mean Values of Weight, SLM, and TBW as per Genotype Distribution at the FTO Variant
The mean values of the three traits (weight, SLM and TBW) associated with the FTO variant showed significant differences (P-value = 0.0047, P-value = 0.019 and P-value = 0.016, respectively) when the cohort was split as per the genotype distributions (TC+CC vs TT) (Figure 1 and Supplementary Table S1). Upon grouping the cohort into non-obese and obese categories, amongs the three traits, only weight in the “obese” category retained statistical significance (P-value = 0.030) between the genotype distributions of (TC+CC) and (TT) (Supplementary Table S3). Thus, the impact due to the effect allele at the FTO variant on SLM, TBW, and weight was likely not due to obese subjects in the cohort. To confirm this, we performed linear regression analysis between each of these three traits and the genotype (TC+CC) with an interaction term involving “obesity status” or “diabetes status” (while correcting for age and sex) and examined the P-values (Supplementary Table S4). The P-values obtained were non-significant, and hence it may be considered that the association signal observed with SLM, TBW, and weight was not due to obesity or diabetes status of the subjects but was entirely due to the genotype of TC+CC, as opposed to TT.
Figure 1 Box-plots displaying data distribution for the six phenotype traits such as (A) SLM, (B) TBW, (C) Weight, (D) Ghrelin, (E) ApoA1, (F) ApoB48 in individuals with genotypes (TC+CC) containing the effect allele or genotypes (TT) homozygous for reference allele at FTO variant rs1421085.
Do Other Traits and Biomarkers Interact With SLM, TBW, and Weight?
Having noted that obesity or diabetes status of the participants does not have an effect on the association signal between SLM, TBW, weight, and the genotypes of TC+CC at the FTO variant, we analyzed the remaining traits and biomarkers of the 60 that we began with using the linear regression model. We found that with TC+CC genotypes (regressed against reference TT genotype), the P-values for interactions of SLM and TBW with ghrelin, ApoA1, and ApoB48 were significant, but for weight, only an interaction with ApoA1 was significant (Table 2). These results illustrate that for TC+CC genotypes at the FTO variant, levels of SLM, TBW, and weight increased and correspondingly levels of ghrelin, ApoA1, and ApoB48 decreased (Figure 2). Though the levels of body weight, total body water, and soft lean mass were seen correlated (at R2 > 0.20) with non-anthropometric traits such as TGL, HDL, FPG, HbA1c, insulin, PBF(%), leptin, ApoC3, glucagon, adiponectin, adipsin, ApoB48, and ghrelin (Supplementary Table S5) in study cohort, the FTO variant with TC+CC genotypes showed interactions only with ghrelin, ApoA1, and ApoB48” in linear regression adjusted for age, sex, and interacting trait.
Table 2 Linear regression model illustrating the link between the FTO variant rs1421085 and interaction between soft lean mass (SLM), total body water (TBW), and weight with ghrelin, ApoA1, and ApoB48.
Figure 2 Significant interactions (negative correlation) observed in soft lean mass (SLM) (A–C), total body water (TBW) (D–F), and weight (G–I) with ghrelin, ApoA1, and ApoB48 with respect to the TC+CC genotype at FTO variant rs1421085.
Results From Power Calculation
Results from Power calculation under both “gene only” hypothesis and “gene-environment hypothesis are presented in Supplementary Table S6. The power calculation indicated that the study cohort had 80% power to detect associations with the FTO variant and could explain 4.1% variance in the traits of SLM, TBW, and weight. The observed effect sizes of 2.249, 1.796, and 3.953 for the association of the FTO variant with SLM, TBW, and weight, respectively (see Table 1) were less than the expected effect sizes of 2.961, 2.302, and 4.746, respectively, from power calculation under “gene-only” hypothesis; in a similar manner, the observed effect sizes were less than expected effect sizes from power calculation under “gene-environment GxE” hypothesis (see Supplementary Table S6).
Discussion
We showed that the FTO variant, rs1421085, was associated with increased body weight, TBW, and SLM in a cohort of Kuwaiti Arab subjects. Further, associations at borderline P-values were observed with lower levels of C-peptide; the inflammatory cytokines, IL-5, TNFα, and IFNγ and the inflammatory chemokine, eotaxin. These observations corroborate the varied associations made by global GWAS between rs1421085 and various physiological processes, including early onset extreme obesity (Wheeler et al., 2013), BMI abnormalities in childhood (Felix et al., 2016), chronotype measurement (morning person) (Jones et al., 2019), percentage body fat (Hubel et al., 2019), lean body mass (fat free mass), type 2 diabetes (Suzuki et al., 2019), energy intake (dietary macronutrient intake) (Tanaka et al., 2013), and hand grip strength (Tikkanen et al., 2018), that were identified in various ethnic groups.
Most notable genetic variants from FTO associated with obesity traits are a cluster of as many as 20 common SNPs from its first intron — particularly notable are rs1421085 (the study variant), rs9939609, rs1121980, rs8050136, rs3751812, rs17817449, and rs8050135. The study variant is in strong linkage disequilibrium with all the other above-listed variants (R2 > 0.90) except for rs8050135 (R2 = 0.0).
In this study, we observed that individuals carrying the rs1421085_C allele exhibited lower levels of inflammatory markers, including IL-4, -5 and -6, TNFα, IFNγ, and MIP1α, among others, while the levels of BMI, TBW, and SLM, as well as leptin and MMP-9, were higher. Associations between FTO SNPs and biomarkers show differing data in literature. For example, Zimmermann et al. (2011) found in young Danish men that the FTO risk allele, rs9939609, which has a strong linkage disequilibrium with rs1421085, was associated with increased CRP level but had no impact on IL-1, IL-6, and TNFα. However, rs9939609 was associated with increased TNFα in Mexican women with gestational diabetes (Saucedo et al., 2017). The impact of rs9939609 on the fasting and postprandial level of IL-6 and TNFα in Brazilian women revealed that, in the fasting stage, levels of IL-6 and TNFα were similar for the A and T risk alleles but in the postprandial period, while TNFα was not altered in individuals with AA genotype compared to TT and TA genotypes, the serum levels of IL-6 and ghrelin were low and those of leptin were high (Magno et al., 2018). Ethnicity is likely to have an effect on the variation in association of FTO with inflammatory markers and investigating these variants may offer a better understanding of their role in the obesity epidemic affecting the world and particularly countries such as Kuwait and others in the Gulf region.
The three body composition traits found associated in this study with the FTO variant are linked to one another: SLM is the summation of TBW and protein; LBM is the summation of SLM and minerals; and weight is the sum total of LBM and body fat mass. The genotypes of CC or TC (as compared to the TT genotype) at the FTO variant upregulate expression of FTO in skeletal muscle (https://gtexportal.org/home/snp/rs1421085), which is a major component of SLM; it is known that an upregulation of FTO can result in increased food intake leading to obesity (Church et al., 2010). Similarly, association of FTO variants with LBM is known (Livshits et al., 2012). TBW has been observed to be higher in obese women than in non-obese women (Sartorio et al., 2005), and very severely obese individuals have high SLM (along with fat-free mass and waist-hip ratio), followed by severe and moderately obese individuals (Jahanlou and Kouzekanani, 2017).
Apart from the observations that exercise and protein intake increase muscle mass, the main determinants of adult muscle mass are yet to be fully characterised. Family and twin studies (Hsu et al., 2005) have suggested that LBM has a high heritability estimate (range, 0.52–0.60). However, to date, only few genes (HSD17B11, VCAN, and ADAMTSL3) have been linked to lean mass (Zillikens et al., 2017) albeit with minor contributions to the trait. The FTO gene has been linked to BMI as it not only affects fat tissue and mass, as established in the past, but also relates to muscle mass as described in this study and by others (Zillikens et al., 2017). The report of Antonio et al. (2018) found that the FTO variant rs1421085 (subject of our study) was significantly associated with PBF, while our study did not find a significant association. This could be attributed to differences in the cohort studied; the study of Antonio et al. compared non-athletes and athletes, who are likely to show more significant differences in PBF (Antonio et al., 2018).
We showed lower levels of ghrelin, ApoA1, and ApoB48 with the effect allele C at the FTO variant in regulating weight, SLM, and TBW. Whereas the association between FTO and ghrelin has been reported (Karra et al., 2013) earlier, for the first time, to the best of our knowledge, we report associations between FTO and ApoA1 and ApoB48.
Ghrelin, termed the “hunger hormone” and leptin, the “satiety hormone”, alter brain responses to food; high plasma levels of ghrelin and low serum levels of leptin have been linked to the regulation of body weight mediated by FTO SNPs (Benedict et al., 2014). ApoA1 is a component of HDL, which is known to transport cholesterol and phospholipids through the bloodstream to the liver. Therefore, our findings were not surprising as FTO may interact with ApoA1 in regulating weight, TBW, and SLM and thus affects obesity and also corroborate studies showing SNPs from ApoA1 associated with lower levels of HDL as obesity risk factors (Rashid and Genest, 2007; Chen et al., 2009). In addition, this is also supported by our previous observations showing interactions between Apo-A1/HDL and a genetic variant in TCN2 regulate WC (Hebbar et al., 2017). Others have demonstrated a positive correlation between ApoB48 and abdominal circumference (Kuo et al., 2019), with obese and hyperlipidemic subjects showing higher levels of ApoB48 (Otokozawa et al., 2009).
SLM contributes to LBM. Low LBM has been shown to be associated with a higher risk of cardiometabolic disorders (Khazem et al., 2018) such as cardiovascular events (Hioki et al., 2015), diabetes (Kalyani et al., 2014), and metabolic syndrome (Kim and Park, 2018). Further, LBM is an important determinant of physical strength, mobility, stamina, and balance (Janssen et al., 2002; Visser et al., 2002; Janssen, 2006) and loss of LBM is associated with sarcopenia (Hansen et al., 1999; Chen et al., 2014; Cruz-Jentoft et al., 2019). It is known that low lean body mass has a strong genetic component, with heritability ranging over 50% (Hsu et al., 2005; Keen-Kim et al., 2006). However, the specific genes underlying the variation in low lean body mass are largely unknown (Urano and Inoue, 2015). Thus, the observation made in our study on FTO variant associated with soft lean mass is novel and important in the field of sarcopenia. Further, physical activity is known to be beneficial in general for the prevention of sarcopenia (Steffl et al., 2017). In this context, it is interesting to note that the studied FTO variant has been implicated in affecting BMI via gene-by-environment interactions such as overall diet and physical activity, alcohol consumption, and deviations from mean sleep duration (Young et al., 2016; Chang et al., 2018). Thus, the finding that the FTO variant rs1421085 is a determinant of soft lean mass (which is primarily made of skeletal and smooth muscle mass) is important.
The MAF of the rs1421085_T > C variant is 0.23 in 1,000 Genomes Project populations (Genomes Project et al., 2015) and varies considerably across the population with the highest value of 0.43 observed in Europeans and the lowest value of 0.06 observed in Africans. Comparatively, allele frequencies of the variant among individuals in Kuwait and Qatar are 0.42 and 0.49, respectively, similar to that seen in Europeans.
Small sample size can be considered as a limitation of the study; however, it is a novel cohort of Arab participants with deep phenotyping. The SNP frequency is very large (40%) and results of power calculation support the notion that the sample size is adequate to support the association of the FTO variant with body weight and lean mass and the interactions with ApoA1, ApoB48, and ghrelin levels.
In conclusion, we report that the FTO variant, rs1421085, is a genetic factor for the body composition traits of TBW, SLM, and body weight in Kuwaiti Arabs. Our findings support an association between the FTO variant and plasma levels of ghrelin, ApoA1, and ApoB48.
Data Availability Statement
The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher. Requests for access to data may be addressed to Dr JA at jehad.abubakr@dasmaninstitute.org.
Ethics Statement
The studies involving human participants were reviewed and approved by Ethical Review Committee of Dasman Diabetes Institute. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
JA, FA-M, and TAT are responsible for the study design. MA-F and AM performed participants recruitment and collection of samples and phenotype. FA and MM performed genotyping experiments. PH carried out data analysis and contributed to interpreting the results. TAT, MA-F, and PH prepared the manuscript. Manuscript was critically reviewed and approved by FA-M. All authors have read and approved the final manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors thank the Tissue Bank Core Facility for processing the samples and extracting the DNA. The Kuwait Foundation for Advancement of the Sciences provided institutional funding that facilitated the study projects (RA 2016-025 and RA 2016-026).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.01411/full#supplementary-material
References
Abu-Farha, M., Behbehani, K., Elkum, N. (2014). Comprehensive analysis of circulating adipokines and hsCRP association with cardiovascular disease risk factors and metabolic syndrome in Arabs. Cardiovasc. Diabetol. 13, 76. doi: 10.1186/1475-2840-13-76
Abu-Farha, M., Abubaker, J., Al-Khairi, I., Cherian, P., Noronha, F., Kavalakatt, S., et al. (2016). Circulating angiopoietin-like protein 8 (betatrophin) association with HsCRP and metabolic syndrome. Cardiovasc. Diabetol. 15, 25. doi: 10.1186/s12933-016-0346-0
Antonio, J., Knafo, S., Kapoor, R., Tartar, J. L. (2018). A fat mass and obesity-associated gene polymorphism influences fat mass in exercise-trained individuals. J. Int. Soc. Sports Nutr. 15 (1), 40. doi: 10.1186/s12970-018-0246-7
Bartosovic, M., Molares, H. C., Gregorova, P., Hrossova, D., Kudla, G., Vanacova, S. (2017). N6-methyladenosine demethylase FTO targets pre-mRNAs and regulates alternative splicing and 3’-end processing. Nucleic Acids Res. 45 (19), 11356–11370. doi: 10.1093/nar/gkx778
Benedict, C., Axelsson, T., Soderberg, S., Larsson, A., Ingelsson, E., Lind, L., et al. (2014). Fat mass and obesity-associated gene (FTO) is linked to higher plasma levels of the hunger hormone ghrelin and lower serum levels of the satiety hormone leptin in older adults. Diabetes 63 (11), 3955–3959. doi: 10.2337/db14-0470
Bera, T. K. (2014). Bioelectrical impedance methods for noninvasive health monitoring: a review. J. Med. Eng. 2014, 381251. doi: 10.1155/2014/381251
Buniello, A., MacArthur, J. A. L., Cerezo, M., Harris, L. W., Hayhurst, J., Malangone, C., et al. (2019). The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 47 (D1), D1005–D1012. doi: 10.1093/nar/gky1120
Chang, J. Y., Park, J. H., Park, S. E., Shon, J., Park, Y. J. (2018). The fat mass- and obesity-associated (FTO) gene to obesity: lessons from mouse models. Obes. (Silver Spring) 26 (11), 1674–1686. doi: 10.1002/oby.22301
Chen, E. S., Mazzotti, D. R., Furuya, T. K., Cendoroglo, M. S., Ramos, L. R., Araujo, L. Q., et al. (2009). Apolipoprotein A1 gene polymorphisms as risk factors for hypertension and obesity. Clin. Exp. Med. 9 (4), 319–325. doi: 10.1007/s10238-009-0051-3
Chen, L. K., Liu, L. K., Woo, J., Assantachai, P., Auyeung, T. W., Bahyah, K. S., et al. (2014). Sarcopenia in Asia: consensus report of the Asian working group for Sarcopenia. J. Am. Med. Dir. Assoc. 15 (2), 95–101. doi: 10.1016/j.jamda.2013.11.025
Church, C., Moir, L., McMurray, F., Girard, C., Banks, G. T., Teboul, L., et al. (2010). Overexpression of FTO leads to increased food intake and results in obesity. Nat. Genet. 42 (12), 1086–1092. doi: 10.1038/ng.713
Claussnitzer, M., Dankel, S. N., Kim, K. H., Quon, G., Meuleman, W., Haugen, C., et al. (2015). FTO obesity variant circuitry and adipocyte browning in humans. N. Engl. J. Med. 373 (10), 895–907. doi: 10.1056/NEJMoa1502214
Cruz-Jentoft, A. J., Bahat, G., Bauer, J., Boirie, Y., Bruyere, O., Cederholm, T., et al. (2019). Sarcopenia: revised European consensus on definition and diagnosis. Age Ageing 48 (1), 16–31. doi: 10.1093/ageing/afy169
Deng, X., Su, R., Stanford, S., Chen, J. (2018). Critical enzymatic functions of FTO in obesity and Cancer. Front. Endocrinol. (Lausanne) 9, 396. doi: 10.3389/fendo.2018.00396
Dina, C., Meyre, D., Gallina, S., Durand, E., Korner, A., Jacobson, P., et al. (2007). Variation in FTO contributes to childhood obesity and severe adult obesity. Nat. Genet. 39 (6), 724–726. doi: 10.1038/ng2048
Felix, J. F., Bradfield, J. P., Monnereau, C., van der Valk, R. J., Stergiakouli, E., Chesi, A., et al. (2016). Genome-wide association analysis identifies three new susceptibility loci for childhood body mass index. Hum. Mol. Genet. 25 (2), 389–403. doi: 10.1093/hmg/ddv472
Frayling, T. M., Timpson, N. J., Weedon, M. N., Zeggini, E., Freathy, R. M., Lindgren, C. M., et al. (2007). A common variant in the FTO gene is associated with body mass index and predisposes to childhood and adult obesity. Science 316 (5826), 889–894. doi: 10.1126/science.1141634
Genomes Project, C., Auton, A., Brooks, L. D., Durbin, R. M., Garrison, E. P., Kang, H. M., et al. (2015). A global reference for human genetic variation. Nature 526 (7571), 68–74. doi: 10.1038/nature15393
Gerken, T., Girard, C. A., Tung, Y. C., Webby, C. J., Saudek, V., Hewitson, K. S., et al. (2007). The obesity-associated FTO gene encodes a 2-oxoglutarate-dependent nucleic acid demethylase. Science 318 (5855), 1469–1472. doi: 10.1126/science.1151710
Hansen, R. D., Raja, C., Aslani, A., Smith, R. C., Allen, B. J. (1999). Determination of skeletal muscle and fat-free mass by nuclear and dual-energy x-ray absorptiometry methods in men and women aged 51-84 y (1-3). Am. J. Clin. Nutr. 70 (2), 228–233. doi: 10.1093/ajcn.70.2.228
Hebbar, P., Alkayal, F., Nizam, R., Melhem, M., Elkum, N., John, S. E., et al. (2017). The TCN2 variant of rs9606756 [Ile23Val] acts as risk loci for obesity-related traits and mediates by interacting with Apo-A1. Obes. (Silver Spring) 25 (6), 1098–1108. doi: 10.1002/oby.21826
Hertel, J. K., Johansson, S., Sonestedt, E., Jonsson, A., Lie, R. T., Platou, C. G., et al. (2011). FTO, type 2 diabetes, and weight gain throughout adult life: a meta-analysis of 41,504 subjects from the Scandinavian HUNT, MDC, and MPP studies. Diabetes 60 (5), 1637–1644. doi: 10.2337/db10-1340
Hioki, H., Miura, T., Motoki, H., Kobayashi, H., Kobayashi, M., Nakajima, H., et al. (2015). Lean body mass index prognostic value for cardiovascular events in patients with coronary artery disease. Heart Asia 7 (2), 12–18. doi: 10.1136/heartasia-2015-010644
Hsu, F. C., Lenchik, L., Nicklas, B. J., Lohman, K., Register, T. C., Mychaleckyj, J., et al. (2005). Heritability of body composition measured by DXA in the diabetes heart study. Obes. Res. 13 (2), 312–319. doi: 10.1038/oby.2005.42
Hubel, C., Gaspar, H. A., Coleman, J. R. I., Finucane, H., Purves, K. L., Hanscombe, K. B., et al. (2019). Genomics of body fat percentage may contribute to sex bias in anorexia nervosa. Am. J. Med. Genet. B. Neuropsychiatr. Genet. 180 (6), 428–438. doi: 10.1002/ajmg.b.32709
Imamura, M., Takahashi, A., Yamauchi, T., Hara, K., Yasuda, K., Grarup, N., et al. (2016). Genome-wide association studies in the Japanese population identify seven novel loci for type 2 diabetes. Nat. Commun. 7, 10531. doi: 10.1038/ncomms10531
Jahanlou, A. S., Kouzekanani, K. (2017). The interaction effect of body mass index and age on fat-free mass, waist-to-hip ratio, and soft lean mass. J. Res. Med. Sci. 22, 86. doi: 10.4103/jrms.JRMS_335_15
Janssen, I., Heymsfield, S. B., Ross, R. (2002). Low relative skeletal muscle mass (sarcopenia) in older persons is associated with functional impairment and physical disability. J. Am. Geriatr. Soc. 50 (5), 889–896. doi: 10.1046/j.1532-5415.2002.50216.x
Janssen, I. (2006). Influence of sarcopenia on the development of physical disability: the cardiovascular health study. J. Am. Geriatr. Soc. 54 (1), 56–62. doi: 10.1111/j.1532-5415.2005.00540.x
Jones, S. E., Lane, J. M., Wood, A. R., van Hees, V. T., Tyrrell, J., Beaumont, R. N., et al. (2019). Genome-wide association analyses of chronotype in 697,828 individuals provides insights into circadian rhythms. Nat. Commun. 10 (1), 343. doi: 10.1038/s41467-018-08259-7
Kalyani, R. R., Tra, Y., Egan, J. M., Ferrucci, L., Brancati, F. (2014). Hyperglycemia is associated with relatively lower lean body mass in older adults. J. Nutr. Health Aging 18 (8), 737–743. doi: 10.1007/s12603-014-0445-0
Karra, E., O’Daly, O. G., Choudhury, A. I., Yousseif, A., Millership, S., Neary, M. T., et al. (2013). A link between FTO, ghrelin, and impaired brain food-cue responsivity. J. Clin. Invest. 123 (8), 3539–3551. doi: 10.1172/JCI44403
Keen-Kim, D., Mathews, C. A., Reus, V. I., Lowe, T. L., Herrera, L. D., Budman, C. L., et al. (2006). Overrepresentation of rare variants in a specific ethnic group may confuse interpretation of association analyses. Hum. Mol. Genet. 15 (22), 3324–3328. doi: 10.1093/hmg/ddl408
Khazem, S., Itani, L., Kreidieh, D., El Masri, D., Tannir, H., Citarella, R., et al. (2018). Reduced lean body mass and cardiometabolic diseases in adult males with overweight and obesity: a pilot study. Int. J. Environ. Res. Public Health 15 (12), 2754. doi: 10.3390/ijerph15122754
Kichaev, G., Bhatia, G., Loh, P. R., Gazal, S., Burch, K., Freund, M. K., et al. (2019). Leveraging polygenic functional enrichment to improve GWAS power. Am. J. Hum. Genet. 104 (1), 65–75. doi: 10.1016/j.ajhg.2018.11.008
Kim, K., Park, S. M. (2018). Association of muscle mass and fat mass with insulin resistance and the prevalence of metabolic syndrome in Korean adults: a cross-sectional study. Sci. Rep. 8 (1), 2703. doi: 10.1038/s41598-018-21168-5
Kuo, I. C., Wu, P. H., Lin, H. Y., Niu, S. W., Huang, J. C., Hung, C. C., et al. (2019). The association of adiponectin with metabolic syndrome and clinical outcome in patients with non-diabetic chronic kidney disease. PloS One 14 (7), e0220158. doi: 10.1371/journal.pone.0220158
Li, J., Ji, L. (2005). Adjusting multiple testing in multilocus analyses using the eigenvalues of a correlation matrix. Hered. (Edinb) 95 (3), 221–227. doi: 10.1038/sj.hdy.6800717
Li, H., Kilpelainen, T. O., Liu, C., Zhu, J., Liu, Y., Hu, C., et al. (2012). Association of genetic variation in FTO with risk of obesity and type 2 diabetes with data from 96,551 East and South Asians. Diabetologia 55 (4), 981–995. doi: 10.1007/s00125-011-2370-7
Livshits, G., Malkin, I., Moayyeri, A., Spector, T. D., Hammond, C. J. (2012). Association of FTO gene variants with body composition in UK twins. Ann. Hum. Genet. 76 (5), 333–341. doi: 10.1111/j.1469-1809.2012.00720.x
Loos, R. J., Yeo, G. S. (2014). The bigger picture of FTO: the first GWAS-identified obesity gene. Nat. Rev. Endocrinol. 10 (1), 51–61. doi: 10.1038/nrendo.2013.227
Magno, F., Guarana, H. C., Fonseca, A. C. P., Cabello, G. M. K., Carneiro, J. R. I., Pedrosa, A. P., et al. (2018). Influence of FTO rs9939609 polymorphism on appetite, ghrelin, leptin, IL6, TNFalpha levels, and food intake of women with morbid obesity. Diabetes Metab. Syndr. Obes. 11, 199–207. doi: 10.2147/DMSO.S154978
Otokozawa, S., Ai, M., Diffenderfer, M. R., Asztalos, B. F., Tanaka, A., Lamon-Fava, S., et al. (2009). Fasting and postprandial apolipoprotein B-48 levels in healthy, obese, and hyperlipidemic subjects. Metabolism 58 (11), 1536–1542. doi: 10.1016/j.metabol.2009.04.040
Parnell, L. D., Blokker, B. A., Dashti, H. S., Nesbeth, P. D., Cooper, B. E., Ma, Y., et al. (2014). CardioGxE, a catalog of gene-environment interactions for cardiometabolic traits. BioData Min. 7, 21. doi: 10.1186/1756-0381-7-21
Price, A. L., Spencer, C. C., Donnelly, P. (2015). Progress and promise in understanding the genetic basis of common diseases. Proc. Biol. Sci. 282 (1821), 20151684. doi: 10.1098/rspb.2015.1684
Rashid, S., Genest, J. (2007). Effect of obesity on high-density lipoprotein metabolism. Obes. (Silver Spring) 15 (12), 2875–2888. doi: 10.1038/oby.2007.342
Sartorio, A., Malavolti, M., Agosti, F., Marinone, P. G., Caiti, O., Battistini, N., et al. (2005). Body water distribution in severe obesity and its assessment from eight-polar bioelectrical impedance analysis. Eur. J. Clin. Nutr. 59 (2), 155–160. doi: 10.1038/sj.ejcn.1602049
Saucedo, R., Valencia, J., Gutierrez, C., Basurto, L., Hernandez, M., Puello, E., et al. (2017). Gene variants in the FTO gene are associated with adiponectin and TNF-alpha levels in gestational diabetes mellitus. Diabetol. Metab. Syndr. 9, 32. doi: 10.1186/s13098-017-0234-0
Steffl, M., Bohannon, R. W., Sontakova, L., Tufano, J. J., Shiells, K., Holmerova, I. (2017). Relationship between sarcopenia and physical activity in older people: a systematic review and meta-analysis. Clin. Interv. Aging 12, 835–845. doi: 10.2147/CIA.S132940
Suzuki, K., Akiyama, M., Ishigaki, K., Kanai, M., Hosoe, J., Shojima, N., et al. (2019). Identification of 28 new susceptibility loci for type 2 diabetes in the Japanese population. Nat. Genet. 51 (3), 379–386. doi: 10.1038/s41588-018-0332-4
Szendroedi, J., Saxena, A., Weber, K. S., Strassburger, K., Herder, C., Burkart, V., et al. (2016). Cohort profile: the German diabetes study (GDS). Cardiovasc. Diabetol. 15, 59. doi: 10.1186/s12933-016-0374-9
Takeuchi, F., Akiyama, M., Matoba, N., Katsuya, T., Nakatochi, M., Tabara, Y., et al. (2018). Interethnic analyses of blood pressure loci in populations of East Asian and European descent. Nat. Commun. 9 (1), 5052. doi: 10.1038/s41467-018-07345-0
Tanaka, T., Ngwa, J. S., van Rooij, F. J., Zillikens, M. C., Wojczynski, M. K., Frazier-Wood, A. C., et al. (2013). Genome-wide meta-analysis of observational studies shows common genetic variants associated with macronutrient intake. Am. J. Clin. Nutr. 97 (6), 1395–1402. doi: 10.3945/ajcn.112.052183
Tikkanen, E., Gustafsson, S., Amar, D., Shcherbina, A., Waggott, D., Ashley, E. A., et al. (2018). Biological insights into muscular strength: genetic findings in the UK biobank. Sci. Rep. 8 (1), 6451. doi: 10.1038/s41598-018-24735-y
Urano, T., Inoue, S. (2015). Recent genetic discoveries in osteoporosis, sarcopenia and obesity. Endocr. J. 62 (6), 475–484. doi: 10.1507/endocrj.EJ15-0154
Visser, M., Kritchevsky, S. B., Goodpaster, B. H., Newman, A. B., Nevitt, M., Stamm, E., et al. (2002). Leg muscle mass and composition in relation to lower extremity performance in men and women aged 70 to 79: the health, aging and body composition study. J. Am. Geriatr. Soc. 50 (5), 897–904. doi: 10.1046/j.1532-5415.2002.50217.x
Wheeler, E., Huang, N., Bochukova, E. G., Keogh, J. M., Lindsay, S., Garg, S., et al. (2013). Genome-wide SNP and CNV analysis identifies common and low-frequency variants associated with severe early-onset obesity. Nat. Genet. 45 (5), 513–517. doi: 10.1038/ng.2607
Wu, C., Arora, P. (2016). Noncoding genome-wide association studies variant for obesity: inroads into mechanism: an overview from the AHA’s council on functional genomics and translational biology. J. Am. Heart Assoc. 5 (7), e003060. doi: 10.1161/JAHA.115.003060
Young, A. I., Wauthier, F., Donnelly, P. (2016). Multiple novel gene-by-environment interactions modify the effect of FTO variants on body mass index. Nat. Commun. 7, 12724. doi: 10.1038/ncomms12724
Zhao, F., Song, M., Wang, Y., Wang, w. (2016). Genetic model. J. Cell Mol. Med. 20 (4), 765. doi: 10.1111/jcmm.12751
Zillikens, M. C., Demissie, S., Hsu, Y. H., Yerges-Armstrong, L. M., Chou, W. C., Stolk, L., et al. (2017). Large meta-analysis of genome-wide association studies identifies five loci for lean body mass. Nat. Commun. 8 (1), 80. doi: 10.1038/s41467-017-00031-7
Keywords: biomarkers, FTO, ghrelin, obesity, soft lean mass, total body water
Citation: Hebbar P, Abu-Farha M, Mohammad A, Alkayal F, Melhem M, Abubaker J, Al-Mulla F and Thanaraj TA (2020) FTO Variant rs1421085 Associates With Increased Body Weight, Soft Lean Mass, and Total Body Water Through Interaction With Ghrelin and Apolipoproteins in Arab Population. Front. Genet. 10:1411. doi: 10.3389/fgene.2019.01411
Received: 14 November 2019; Accepted: 31 December 2019;
Published: 31 January 2020.
Edited by:
Amelia Marti, University of Navarra, SpainReviewed by:
Caren Smith, Jean Mayer USDA Human Nutrition Research Center on Aging at Tufts University, United StatesLarry Parnell, Jean Mayer USDA Human Nutrition Research Center on Aging at Tufts University, United States
Copyright © 2020 Hebbar, Abu-Farha, Mohammad, Alkayal, Melhem, Abubaker, Al-Mulla and Thanaraj. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jehad Abubaker, jehad.abubakr@dasmaninstitute.org; Fahd Al-Mulla, fahd.almulla@dasmaninstitute.org; Thangavel Alphonse Thanaraj, alphonse.thangavel@dasmaninstitute.org
†These authors have contributed equally to this work