Next Article in Journal
Plinia trunciflora Extract Administration Prevents HI-Induced Oxidative Stress, Inflammatory Response, Behavioral Impairments, and Tissue Damage in Rats
Previous Article in Journal
Clinical Impact of Mediterranean Diet Adherence before and after Bariatric Surgery: A Narrative Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Analysis of Disordered Eating Behavior in the Mexican Population

by José Jaime Martínez-Magaña 1, Sandra Hernandez 2, Ana Rosa Garcia 3, Valeria Cardoso-Barajas 3, Emmanuel Sarmiento 3, Beatriz Camarena 2, Alejandro Caballero 4, Laura Gonzalez 4, Jorge Ameth Villatoro-Velazquez 5, Maria Elena Medina-Mora 5, Marycarmen Bustos-Gamiño 5, Clara Fleiz-Bautista 5, Carlos Alfonso Tovilla-Zarate 6, Isela Esther Juárez-Rojop 7, Humberto Nicolini 1,* and Alma Delia Genis-Mendoza 1,*
1
Laboratorio de Genómica de Enfermedades Psiquiátricas y Neurodegenerativas, Instituto Nacional de Medicina Genómica, Mexico City 14610, Mexico
2
Laboratorio de Farmacogenética, Instituto Nacional de Psiquiatría Ramón de la Fuente Muñiz, Mexico City 14370, Mexico
3
Unidad de Investigación, Hospital Psiquiátrico Infantil Juan N. Navarro, Mexico City 14080, Mexico
4
Unidad de Trastornos Alimenticios, Instituto Nacional de Psiquiatría Ramón de la Fuente Muñiz, Mexico City 14370, Mexico
5
Unidad de Análisis de Datos y Encuestas, Instituto Nacional de Psiquiatría Ramón de la Fuente Muñiz, Mexico City 14370, Mexico
6
División Académica Multidisciplinaria de Comalcalco, Universidad Juárez Autónoma de Tabasco, Comalcalco 86654, Mexico
7
División de Ciencias de la Salud, Universidad Juárez Autónoma de Tabasco, Villahermosa 86100, Mexico
*
Authors to whom correspondence should be addressed.
Submission received: 7 October 2021 / Revised: 4 January 2022 / Accepted: 5 January 2022 / Published: 17 January 2022
(This article belongs to the Section Nutrigenetics and Nutrigenomics)

Abstract

:
Alterations in eating behavior characterized eating disorders (ED). The genetic factors shared between ED diagnoses have been underexplored. The present study performed a genome-wide association study in individuals with disordered eating behaviors in the Mexican population, blood methylation quantitative trait loci (blood-meQTL), summary data-based Mendelian randomization (SMR) analysis, and in silico function prediction by different algorithms. The analysis included a total of 1803 individuals. We performed a genome-wide association study and blood-meQTL analysis by logistic and linear regression. In addition, we analyzed in silico functional variant prediction, phenome-wide, and multi-tissue expression quantitative trait loci. The genome-wide association study identified 44 single-nucleotide polymorphisms (SNP) associated at a nominal value and seven blood-meQTL at a genome-wide threshold. The SNPs show enrichment in genome-wide associations of the metabolic and immunologic domains. In the in silico analysis, the SNP rs10419198 (p-value = 4.85 × 10−5) located on an enhancer mark could change the expression of PRR12 in blood, adipocytes, and brain areas that regulate food intake. Additionally, we found an association of DNA methylation levels of SETBP1 (p-value = 6.76 × 10−4) and SEMG1 (p-value = 5.73 × 10−4) by SMR analysis. The present study supports the previous associations of genetic variation in the metabolic domain with ED.

1. Introduction

Eating disorders (ED) significantly impact individuals’ physical and mental health [1,2,3]. The manifestation of disordered eating behaviors characterizes ED. Disordered eating is the presence of aversive food-related behaviors, including fasting, vomiting, binge eating, and food intake restriction [4,5]. The three primary ED diagnoses are anorexia nervosa, bulimia nervosa, and binge-eating disorder [6]. Individuals diagnosed with ED have different disordered eating patterns, which could change over time [7,8,9]. Even when the diagnosis criteria for ED are well established [6], affected individuals may have a diagnosis crossover [10,11,12]. The diagnosis crossover could result from the close relationship between ED diagnosis criteria and the overlapping of symptoms [8]. Others suggest that this crossover could be an escalation mechanism, mainly in individuals diagnosed with the binge-eating disorder [13].
The etiology of ED is complex, with the interaction between environmental and biological risk factors [14]. Some identified risk factors are genetics, gender, dieting, or weight concerns, occurring early in life or adolescence [15]. In the last decade, different studies have explored the genetic risk factors that underlie many psychiatric disorders, but the analysis of the genetics of ED had been underexplored [16,17]. Furthermore, genome-wide association studies until now had been performed only in individuals diagnosed with AN [18,19,20,21,22]. These studies suggest that AN has a genetic correlation with other psychiatric phenotypes or metabolic traits [18,22]. Nevertheless, the genetic factors shared between ED diagnoses are less explored.
The environment had a role in the etiology of ED. The translation of the environmental risk factors into molecular changes is the focus of epigenetics. DNA methylation is the most studied epigenetic molecular modification, and several studies have reported alterations in DNA methylation in individuals diagnosed with ED [23]. DNA methylation levels depend on the environment and genetic variation, and integrative analysis of genetic variants and DNA methylation levels could identify new biomolecules involved in the etiology of ED. Some authors had developed statistical frameworks to integrate these two sources of genomic information; an example is summary-based Mendelian randomization (SMR) [24]. However, the integration of these two sources of genomic information in individuals diagnosed with ED had not been performed in the Mexican population. The present work aimed to perform a genome-wide association analysis, blood methylation quantitative trait loci (blood-meQTL), summary data-based Mendelian randomization (SMR) analysis, and in silico function prediction in individuals diagnosed with disordered eating behaviors in the Mexican population.

2. Materials and Methods

2.1. Sample Population

The study included a total of 1803 individuals of Mexican descent, from three different samples: an adolescent clinical subsample from the Mexican Genomic Database for Cross-Disorder Research (n = 168, MeDaCrosR) [25,26], a clinical sample recruited from the ED Unit of the Instituto Nacional de Psiquiatría Ramón de la Fuente Muñiz (n = 166, INPRFM), and an epidemiological subsample from the Mexican Genomic Database for Addiction Research (n = 1469; MxGDAR) [27,28] (Table 1).
Adolescents from the MeDaCrosR were recruited from the Hospital Psiquiátrico Infantil Juan N Navarro (HPIJNN), from the external consultation area. Eating patterns in MeDaCrosR adolescents were evaluated by children’s specialized psychiatrists through the Eating Attitude Test 26 (EAT-26) [29,30] and the Questionnaire on Eating and Weight Patterns Revised (QWEP-R) in Spanish [31]. The adult sample from the INPRFM was diagnosed according to DSM-IV-TR criteria for ED using the Structured Clinical Interview for Mental Disorders v.2.0. (SCID-II) [32] and evaluated with the Eating Disorder Inventory v.2.0. (EDI-2) for psychological patterns [29,30]. The cases of disordered eating were the groups from MeDaCrosR and INPRFM (n = 333). The population-based controls (n = 1802) were the individuals from the MxGDAR. The Diagnostic Interview for Psychosis and Affective Disorders (DIPAD) screening section was used to evaluate the individuals from MxGDAR [28,33,34,35,36]. The study was performed based in the Helsinki declaration, and was revised by the ethic and investigation committees of the Instituto Nacional de Medicina Genómica (Approval Number CEI/2018/60), Instituto Nacional de Psiquiatría Ramón de la Fuente Muñiz (Approval Number: DGC-279-2008), and the Hospital Psiquiátrico Infantil Juan N Navarro (Approval number II3/01/0913). Every individual fulfilled and signed an informed assent (adolescents) and/or consent (parents/adults).

2.2. Genome-Wide Typing

We extracted DNA from blood and buccal epithelial samples with a modified salting-out method, implemented in the commercial kit Gentra Puregene (Qiagen, Readwood City, CA, USA). According to the manufacturer’s protocol, genotyping was performed in the Instituto Nacional de Medicina Genomica with the commercial Infinium PsychArray Beadchip (Illumina, San Diego, CA, USA). Fluorescent intensities were measured with the iScan (Illumina, San Diego, CA, USA), transformed to genotypes with the GenomeStudio (Illumina, San Diego, CA, USA), and converted to Plink format files. Quality control of genotypes was performed in the Plink software [37,38], based on previously published protocols, briefly, we considered the following criteria: variant calling greater than 95%, a minor allele frequency (MAF) greater than 5%, a Hardy–Weinberg equilibrium Chi-square test p-value greater than 1 × 10−5, and remotion of variants with A/T or G/C alleles (to avoid the flip strand effect). Additionally, we excluded individuals with a genotype call rate of less than 95%. In addition, all individual pairs with an identity-by-state value greater than 1.6 were marked to correct for cryptic relationships, and the individual with the lowest genotype call rate was excluded.

2.3. Statistical Analysis

2.3.1. Population Stratification

We evaluated population stratification using previously reported algorithms and the PC-AiR package [39]. The reference panel was the Human Genome Diversity Project (HGDP) [40]. We included only independent SNPs by filtering with linkage disequilibrium pruning (LD pruning), using the following parameters: a window size of 50 Kb, a step of 2, and a variance inflation factor of 5 implemented in Plink.

2.3.2. Genome-Wide Associations Analysis

We performed genetic associations through multiple logistic regressions, adjusted for age, sex, and ten components of global ancestry as covariables in Plink. We considered a p-value of 5.00 × 10−5 as nominally associated, and a p-value of 5.00 × 10−8 was considered statistically significant on the genome-wide level. After statistical contrasts, we removed all the variants with a MAF lower than 5.0% in cases or controls. The final results included a total of 223,189 SNPs.

2.3.3. Functional Prediction

An in silico functional annotation was also carried out of the associated SNPs using a variant effect predictor (VEP) [41]. A phenome-wide association study was performed in the GWAS atlas of the associated SNPs, and we considered a statistical significance threshold for a p-value < 5 × 10−8 [42]. Pathway analysis was carried out with the online ComPath tool [43]. Finally, we searched the GTEx portal for in silico multi-tissue expression quantitative trait loci analysis (eQTL) results of the associated variants [44,45].

2.3.4. Methylation Quantitative Trait Loci Analysis

We performed blood methylation quantitative trait loci (blood-meQTL) analysis to estimate the possible impact of the associated SNPs in DNA methylation levels. We analyzed a subset of individuals diagnosed with eating disorders of Mexican ascendence from a previously published database [46]. The DNA methylation levels were determined using Illumina MethylationEPIC BeadChips (Illumina, USA). We calculated trans- and cis-meQTL with the -mtscore function of the KING software [47], and a p-value < 5.00 × 10−8 was considered genome-wide statistically significant. We searched the associated SNPs with eQTL and blood-meQTL effect in the regulomeDB for enhancer or promoter marks [48].

2.3.5. Summary Data-Based Mendelian Randomization (SMR)

We applied summary data-based Mendelian randomization analysis (SMR) to identify DNA methylation levels associated with disordered eating by pleiotropy [49]. SMR was developed for quantitative expression loci and uses the principles of Mendelian Randomization (MR) to estimate the pleiotropic association between the level of gene expression and a phenotype. Zhu et al. established that if we denote z as a genetic variant (for example, SNP), x as the expression level of a gene, and y as the trait, then the two-step least-squares estimate of the effect of x on y from an MR analysis is:
bxy = bzy/bzx
where bzy and bzx are the least-squares estimates of y and x on z, respectively, and bxy is interpreted as the effect size of x on y free of confounding from non-genetic factors. After calculating the sampling variance by the Delta method and replacing bzy and bzx by the estimates, we should have:
T SMR = b xy 2 / var ( b xy ) = z zy 2 · z zx 2 / z zy 2 + z zx 2
where zzy and zzx are the z statistics from the GWAS and eQTL study, respectively. The previous method was developed by Zhu et al., and we used the SMR to integrate the GWAS and blood-meQTL summary statistics.

3. Results

In the genotype-phenotype associations, no genome-wide associations were found; nevertheless, 44 single-nucleotide polymorphisms (SNP) were associated at a nominal level (p-value < 5 × 10−5) (Table 2). The SNPs distributes in 29 cytogenetic bands and 19 coding regions. The protein-coding genes with SNPs associated with disordered eating showed enrichment in the PI3K-Akt signaling pathway (hsa04151, adjusted p-value = 0.0499, FTL3LG and TSC2), vascular muscle contraction (hsa04270, adjusted p-value = 0.0494, KCNMA1 and PRKCE), and cGMP–PKG signaling pathway (hsa04022, adjusted p-value = 0.0494, KCNMA1 and PRKCE). Of the 44 SNPs, 22 were intronic, and 2 were missense variants. The missense variants were the SRMM4 p.Ser243Asn and FLT3LG p.Phe177Leu.
In the PheWAS, we identified the immunological (rs4626924 and rs10419198) and metabolic (rs3205060, rs8041059 and rs10419198) domain with SNPs associated. The SNPs associated with the immunological domains were both intronic, one located in the coding region of LOC107985364 and the other in RCN3. In the metabolic domain, we identified the rs8041059 located in the intron of LIPC and with previous associations with 33 different lipidic traits.
The blood-meQTLs analysis showed seven SNPs associated with DNA methylation levels of five CpG sites (Table 3). The CpG sites map to four genes, including the rs10419198 related to DNA methylation levels in a CpG site annotated to the body of the PRR12. The blood-mQTLs correlated with five CpG sites, annotated to the gene body of SETBP1 (cg12522870) and PRR12 (cg06378142), 200 base pairs upstream of the transcription starting site (TSS200) of the non-coding LOC440704 (cg12412036 and cg09420738) and 1500 upstream of the SEMG1 (cg15921833). In the case of LOC440704, two different CpGs correlated with two SNPs; meanwhile, for the CpGs on SETBP1 and SEMG1, two SNPs connect with the same site.
The multi-tissue eQTL analysis in the GTEx database revealed that the rs10419198 had an eQTL effect on different tissues (frontal cortex, anterior cingulate cortex, cerebellum, putamen, thyroid, visceral adipose tissue, etc.) (Figure 1). In the regulomeDB database, the rs10419198 showed a mark of enhancers in 51 different cell lines.
Of the previous seven SNPs that showed blood-meQTL significant results, we found that two had a pleiotropic effect with DNA methylation levels and disordered eating by SMR analysis (Table 4). The CpG sites were the cg12522870 (SETBP1) and cg15921833 (SEMG1).

4. Discussion

Previous GWAS on ED has shown, mainly in AN, that genetic variants associated with this disorder could significantly impact metabolic pathways [31,32,33,34]. The genetic association found in our study also supports this finding on individuals with disordered eating behavior. In our study, three variants (rs3205060, rs8041059 and rs10419198) associated with disordered eating are previously associated at the genome-wide level with at least one metabolic trait. One associated variant was the rs8041059 which is found in the coding region of the hepatic triglyceride lipase (LIPC) and is associated with differences in high-density lipoproteins (HDL), cholesterol, and triglyceride plasmatic levels [35,36]. LIPC catalyzes the hydrolysis of triglycerides and phospholipids found in circulating lipoproteins [37,38,39]. Some authors hypothesized that this variant regulates triglyceride concentration, affecting the transcription factors’ binding site in the promoter of LIPC. Still, the exact mechanism remains to be explored [40,41,42]. In addition, the association of this SNP could be important for further explorations in the mechanism behind the reported increased HDL levels in individuals diagnosed with AN [43,44].
The other associated variant, rs10419198, is annotated to the intron of RCN3, and could function as an enhancer for PPR12 [36]. Our blood-meQTL and in silico analysis of multi-tissue eQTL showed that this variant modulates the expression of the PRR12. PRR12 is the proline-rich 12 protein, also known as KIAA1205. The exact cellular function of these proteins remains unknown, but it is proposed that it has an essential role in brain development. The role of this gene on the brain is supported by some studies, where individuals with loss-of-function or structural variation on this gene had clinical syndromes characterized by neurodevelopmental and eye abnormalities [45,46,47]. Individuals carriers of high impact variants on PRR12 had neurodevelopmental disorders, for example, autism, intellectual disability, and attention-deficit/hyperactivity disorder (ADHD). ADHD is one of the most frequent comorbidities found in individuals with ED, and the presence of this clinical entity could increase the symptoms of disordered eating [48,49,50,51]. Further analysis of the function of PRR12 could help us understand this high comorbidity.
Environmental factors interact with genetics to the development of ED [48]. A known molecular mark that captures these environmental differences is DNA methylation [50]. DNA methylation is influenced by the environment and genetic, having differences based on the genotype of the individuals [51,52]. Then if the DNA methylation level affected the trait, we could see an effect of differences in the phenotype, known as pleiotropy [49,53]. We performed SMR to capture this in our analysis and found an association of disordered eating with CpG sites on SETBP1 and SEMG1. Mutations on SETBP1 cause the Schinzel Giedion syndrome, individuals affected had a severe intellectual disability. A recent study identified SETBP1 as a key epigenetic regulator of the gene network responsible for visceral organ and brain morphogenesis [54]. On the other hand, the SEMG1 could be a key modulator of energy metabolism on ED. SEMG1 is the semenogelin 1 gene, expressed mainly on the prostate and is suggested to be responsible for asthenozoospermia [55]. In addition, protein pull-down studies suggest that SEMG1 increases energy expenditure by interacting with enzymes responsible for energy regulation, such as lactate dehydrogenase A [56]. Nevertheless, further studies are required to clarify the molecular mechanism of these proteins in ED.
Our study is one of the first to explore genetic associations with ED in the Mexican population at the genome-wide level; nevertheless, we could state some limitations. The main is the small sample size, and also, we do not have a replication cohort. This reduced sample size reduces our statistical power. Another limitation could be the lack of metabolic parameters measured in the sample and that the population-based sample was not screened for disordered eating behavior. On the other hand, integrating different sources of information, the in silico eQTL and blood-meQTL strengthens it.

5. Conclusions

The present study suggests an association of the rs10419198 enhancer variant of the PRR12 supporting previous reports that have found that disordered eating could imply genetic associations with metabolism.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/nu14020394/s1, summary statistics of the GWAS of disordered eating.

Author Contributions

Conceptualization, J.J.M.-M. and A.D.G.-M.; methodology, A.D.G.-M.; software, J.J.M.-M.; validation, S.H., A.D.G.-M. and H.N.; formal analysis, J.J.M.-M.; investigation, I.E.J.-R., C.A.T.-Z. and A.C.; resources, C.F.-B., M.B.-G., M.E.M.-M., J.A.V.-V., L.G., A.C., B.C. and E.S.; data curation, J.J.M.-M.; writing—original draft preparation, J.J.M.-M. and S.H.; writing—review and editing, all authors; visualization, J.J.M.-M.; supervision, H.N.; project administration, A.D.G.-M.; funding acquisition, H.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Fundación Gonzalo Río Arronte, Grant Number S591, and Instituto Nacional de Medicina Genómica (INMEGEN), Grant Number 06/2018/I.

Institutional Review Board Statement

This study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committees of the Psychiatric Hospital “Juan N. Navarro” (Protocol Code No. II3/01/0913; 11 October 2017), and National Institute of Genomic Medicine (INMEGEN) (Protocol Code No. 06/2018/I; June 2018).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author, omitted due to privacy and ethical issues. The complete summary statistics would be made available as Supplementary Material.

Acknowledgments

We want to acknowledge the Microarray and Expression Unit of the National Institute of Genomic Medicine for their technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Berkman, N.D.; Lohr, K.N.; Bulik, C.M. Outcomes of Eating Disorders: A Systematic Review of the Literature. Int. J. Eat. Disord. 2007, 40, 293–309. [Google Scholar] [CrossRef]
  2. Smink, F.R.E.; van Hoeken, D.; Hoek, H.W. Epidemiology of Eating Disorders: Incidence, Prevalence and Mortality Rates. Curr. Psychiatry Rep. 2012, 14, 406–414. [Google Scholar] [CrossRef] [Green Version]
  3. Santomauro, D.F.; Melen, S.; Mitchison, D.; Vos, T.; Whiteford, H.; Ferrari, A.J. The Hidden Burden of Eating Disorders: An Extension of Estimates from the Global Burden of Disease Study 2019. Lancet Psychiatry 2021, 8, 320–328. [Google Scholar] [CrossRef]
  4. Jalali-Farahani, S.; Chin, Y.S.; Mohd Nasir, M.T.; Amiri, P. Disordered Eating and Its Association with Overweight and Health-Related Quality of Life Among Adolescents in Selected High Schools of Tehran. Child Psychiatry Hum. Dev. 2015, 46, 485–492. [Google Scholar] [CrossRef]
  5. Zeiler, M.; Waldherr, K.; Philipp, J.; Nitsch, M.; Dür, W.; Karwautz, A.; Wagner, G. Prevalence of Eating Disorder Risk and Associations with Health-Related Quality of Life: Results from a Large School-Based Population Screening: Prevalence of Eating Disorder Risk. Eur. Eat. Disord. Rev. 2016, 24, 9–18. [Google Scholar] [CrossRef]
  6. Vo, M.; Accurso, E.C.; Goldschmidt, A.B.; Le Grange, D. The Impact of DSM-5 on Eating Disorder Diagnoses. Int. J. Eat. Disord. 2017, 50, 578–581. [Google Scholar] [CrossRef] [Green Version]
  7. Pearson, C.M.; Miller, J.; Ackard, D.M.; Loth, K.A.; Wall, M.M.; Haynos, A.F.; Neumark-Sztainer, D. Stability and Change in Patterns of Eating Disorder Symptoms from Adolescence to Young Adulthood. Int. J. Eat. Disord. 2017, 50, 748–757. [Google Scholar] [CrossRef]
  8. Wu, X.Y.; Yin, W.Q.; Sun, H.W.; Yang, S.X.; Li, X.Y.; Liu, H.Q. The Association between Disordered Eating and Health-Related Quality of Life among Children and Adolescents: A Systematic Review of Population-Based Studies. PLoS ONE 2019, 14, e0222777. [Google Scholar] [CrossRef]
  9. Eichen, D.M.; Strong, D.R.; Rhee, K.E.; Rock, C.L.; Crow, S.J.; Epstein, L.H.; Wilfley, D.E.; Boutelle, K.N. Change in Eating Disorder Symptoms Following Pediatric Obesity Treatment. Int. J. Eat. Disord. 2019, 52, 299–303. [Google Scholar] [CrossRef]
  10. Tozzi, F.; Thornton, L.M.; Klump, K.L.; Fichter, M.M.; Halmi, K.A.; Kaplan, A.S.; Strober, M.; Woodside, D.B.; Crow, S.; Mitchell, J.; et al. Symptom Fluctuation in Eating Disorders: Correlates of Diagnostic Crossover. Am. J. Psychiatry 2005, 162, 732–740. [Google Scholar] [CrossRef]
  11. Milos, G.; Spindler, A.; Schnyder, U.; Fairburn, C.G. Instability of Eating Disorder Diagnoses: Prospective Study. Br. J. Psychiatry 2005, 187, 573–578. [Google Scholar] [CrossRef] [Green Version]
  12. Castellini, G.; Lo Sauro, C.; Mannucci, E.; Ravaldi, C.; Rotella, C.M.; Faravelli, C.; Ricca, V. Diagnostic Crossover and Outcome Predictors in Eating Disorders According to DSM-IV and DSM-V Proposed Criteria: A 6-Year Follow-Up Study. Psychosom. Med. 2011, 73, 270–279. [Google Scholar] [CrossRef]
  13. Stice, E.; Marti, C.N.; Rohde, P. Prevalence, Incidence, Impairment, and Course of the Proposed DSM-5 Eating Disorder Diagnoses in an 8-Year Prospective Community Study of Young Women. J. Abnorm. Psychol. 2013, 122, 445–457. [Google Scholar] [CrossRef] [Green Version]
  14. Jacobi, C.; Hayward, C.; de Zwaan, M.; Kraemer, H.C.; Agras, W.S. Coming to Terms with Risk Factors for Eating Disorders: Application of Risk Terminology and Suggestions for a General Taxonomy. Psychol. Bull. 2004, 130, 19–65. [Google Scholar] [CrossRef]
  15. Hilbert, A.; Pike, K.M.; Goldschmidt, A.B.; Wilfley, D.E.; Fairburn, C.G.; Dohm, F.-A.; Walsh, B.T.; Striegel Weissman, R. Risk Factors across the Eating Disorders. Psychiatry Res. 2014, 220, 500–506. [Google Scholar] [CrossRef] [Green Version]
  16. Horwitz, T.; Lam, K.; Chen, Y.; Xia, Y.; Liu, C. A Decade in Psychiatric GWAS Research. Mol. Psychiatry 2019, 24, 378–389. [Google Scholar] [CrossRef]
  17. The Brainstorm Consortium; Anttila, V.; Bulik-Sullivan, B.; Finucane, H.K.; Walters, R.K.; Bras, J.; Duncan, L.; Escott-Price, V.; Falcone, G.J.; Gormley, P.; et al. Analysis of Shared Heritability in Common Disorders of the Brain. Science 2018, 360, eaap8757. [Google Scholar] [CrossRef] [Green Version]
  18. The Price Foundation Collaborative Group; Wang, K.; Zhang, H.; Bloss, C.S.; Duvvuri, V.; Kaye, W.; Schork, N.J.; Berrettini, W.; Hakonarson, H. A Genome-Wide Association Study on Common SNPs and Rare CNVs in Anorexia Nervosa. Mol. Psychiatry 2011, 16, 949–959. [Google Scholar] [CrossRef] [Green Version]
  19. Boraska, V.; Franklin, C.S.; Floyd, J.A.B.; Thornton, L.M.; Huckins, L.M.; Southam, L.; Rayner, N.W.; Tachmazidou, I.; Klump, K.L.; Treasure, J.; et al. A Genome-Wide Association Study of Anorexia Nervosa. Mol. Psychiatry 2014, 19, 1085–1094. [Google Scholar] [CrossRef] [Green Version]
  20. Duncan, L.; Yilmaz, Z.; Walters, R.; Goldstein, J.; Anttila, V.; Bulik-Sullivan, B.; Ripke, S.; Thornton, L.; Hinney, A.; Daly, M.; et al. Genome-Wide Association Study Reveals First Locus for Anorexia Nervosa and Metabolic Correlations. Am. J. Psychiatry 2017, 174, 850–858. [Google Scholar] [CrossRef]
  21. Li, D.; Chang, X.; Connolly, J.J.; Tian, L.; Liu, Y.; Bhoj, E.J.; Robinson, N.; Abrams, D.; Li, Y.R.; Bradfield, J.P.; et al. A Genome-Wide Association Study of Anorexia Nervosa Suggests a Risk Locus Implicated in Dysregulated Leptin Signaling. Sci. Rep. 2017, 7, 3847. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Watson, H.J.; Yilmaz, Z.; Thornton, L.M.; Hübel, C.; Coleman, J.R.; Gaspar, H.A.; Bryois, J.; Hinney, A.; Leppä, V.M.; Mattheisen, M.; et al. Genome-Wide Association Study Identifies Eight Risk Loci and Implicates Metabo-Psychiatric Origins for Anorexia Nervosa. Nat. Genet. 2019, 51, 1207–1214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Hübel, C.; Marzi, S.J.; Breen, G.; Bulik, C.M. Epigenetics in Eating Disorders: A Systematic Review. Mol. Psychiatry 2019, 24, 901–915. [Google Scholar] [CrossRef]
  24. Wu, Y.; Zeng, J.; Zhang, F.; Zhu, Z.; Qi, T.; Zheng, Z.; Lloyd-Jones, L.R.; Marioni, R.E.; Martin, N.G.; Montgomery, G.W.; et al. Integrative Analysis of Omics Summary Data Reveals Putative Mechanisms Underlying Complex Traits. Nat. Commun. 2018, 9, 918. [Google Scholar] [CrossRef]
  25. Martínez-Magaña, J.J.; Genis-Mendoza, A.D.; Villatoro Velázquez, J.A.; Bustos-Gamiño, M.; Juárez-Rojop, I.E.; Tovilla-Zarate, C.A.; Sarmiento, E.; Saucedo, E.; Rodríguez-Mayoral, O.; Fleiz-Bautista, C.; et al. Genome-Wide Association Study of Psychiatric and Substance Use Comorbidity in Mexican Individuals. Sci. Rep. 2021, 11, 6771. [Google Scholar] [CrossRef] [PubMed]
  26. Ruiz-Ramos, D.; Martínez-Magaña, J.J.; García, A.R.; Juarez-Rojop, I.E.; Gonzalez-Castro, T.B.; Tovilla-Zarate, C.A.; Sarmiento, E.; López-Narvaez, M.L.; Nicolini, H.; Genis-Mendoza, A.D. Psychiatric Comorbidity in Mexican Adolescents with a Diagnosis of Eating Disorders Its Relationship with the Body Mass Index. Int. J. Environ. Res. Public Health 2021, 18, 3900. [Google Scholar] [CrossRef]
  27. Martínez-Magaña, J.J.; Genis-Mendoza, A.D.; Villatoro Velázquez, J.A.; Camarena, B.; Martín del Campo Sanchez, R.; Fleiz Bautista, C.; Bustos Gamiño, M.; Reséndiz, E.; Aguilar, A.; Medina-Mora, M.E.; et al. The Identification of Admixture Patterns Could Refine Pharmacogenetic Counseling: Analysis of a Population-Based Sample in Mexico. Front. Pharmacol. 2020, 11, 324. [Google Scholar] [CrossRef] [Green Version]
  28. Nicolini, H.; Martínez-Magaña, J.J.; Genis-Mendoza, A.D.; Villatoro Velázquez, J.A.; Camarena, B.; Fleiz Bautista, C.; Bustos-Gamiño, M.; Aguilar García, A.; Lanzagorta, N.; Medina-Mora, M.E. Cannabis Use in People With Obsessive-Compulsive Symptomatology: Results From a Mexican Epidemiological Sample. Front. Psychiatry 2021, 12, 664228. [Google Scholar] [CrossRef]
  29. Garner, D.M.; Olmsted, M.P.; Bohr, Y.; Garfinkel, P.E. The Eating Attitudes Test: Psychometric Features and Clinical Correlates. Psychol. Med. 1982, 12, 871–878. [Google Scholar] [CrossRef]
  30. García-García, E.; Vázquez-Velázquez, V.; López-Alvarenga, J.C.; Arcila-Martínez, D. [Internal validity and diagnostic utility of the Eating Disorder Inventory in Mexican women]. Salud Publica Mex. 2003, 45, 206–210. [Google Scholar] [CrossRef] [Green Version]
  31. Spitzer, R.L.; Yanovski, S.; Wadden, T.; Wing, R.; Marcus, M.D.; Stunkard, A.; Devlin, M.; Mitchell, J.; Hasin, D.; Horne, R.L. Binge Eating Disorder: Its Further Validation in a Multisite Study. Int. J. Eat. Disord. 1993, 13, 137–153. [Google Scholar] [CrossRef]
  32. Spitzer, R.L.; Williams, J.B.; Gibbon, M.; First, M.B. The Structured Clinical Interview for DSM-III-R (SCID). I: History, Rationale, and Description. Arch. Gen. Psychiatry 1992, 49, 624–629. [Google Scholar] [CrossRef]
  33. Endicott, J.; Spitzer, R.L. A Diagnostic Interview: The Schedule for Affective Disorders and Schizophrenia. Arch. Gen. Psychiatry 1978, 35, 837–844. [Google Scholar] [CrossRef]
  34. Hartz, S.M.; Pato, C.N.; Medeiros, H.; Cavazos-Rehg, P.; Sobell, J.L.; Knowles, J.A.; Bierut, L.J.; Pato, M.T. Comorbidity of Severe Psychotic Disorders With Measures of Substance Use. JAMA Psychiatry 2014, 71, 248. [Google Scholar] [CrossRef] [Green Version]
  35. Docherty, A.R.; Bigdeli, T.B.; Edwards, A.C.; Bacanu, S.; Lee, D.; Neale, M.C.; Wormley, B.K.; Walsh, D.; O’Neill, F.A.; Riley, B.P.; et al. Genome-Wide Gene Pathway Analysis of Psychotic Illness Symptom Dimensions Based on a New Schizophrenia-Specific Model of the OPCRIT. Schizophr. Res. 2015, 164, 181–186. [Google Scholar] [CrossRef] [Green Version]
  36. Pato, M.T.; Sobell, J.L.; Medeiros, H.; Abbott, C.; Sklar, B.M.; Buckley, P.F.; Bromet, E.J.; Escamilla, M.A.; Fanous, A.H.; Lehrer, D.S.; et al. The Genomic Psychiatry Cohort: Partners in Discovery. Am. J. Med. Genet. B Neuropsychiatr. Genet. 2013, 162, 306–312. [Google Scholar] [CrossRef] [Green Version]
  37. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Chang, C.C.; Chow, C.C.; Tellier, L.C.; Vattikuti, S.; Purcell, S.M.; Lee, J.J. Second-Generation PLINK: Rising to the Challenge of Larger and Richer Datasets. GigaScience 2015, 4, s13742-015. [Google Scholar] [CrossRef] [PubMed]
  39. Conomos, M.P.; Miller, M.B.; Thornton, T.A. Robust Inference of Population Structure for Ancestry Prediction and Correction of Stratification in the Presence of Relatedness. Genet. Epidemiol. 2015, 39, 276–293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Harry, D. The Human Genome Diversity Project and Its Implications for Indigenous Peoples. Genewatch Bull. Comm. Responsib. Genet. 1996, 10, 8–9. [Google Scholar]
  41. McLaren, W.; Gil, L.; Hunt, S.E.; Riat, H.S.; Ritchie, G.R.S.; Thormann, A.; Flicek, P.; Cunningham, F. The Ensembl Variant Effect Predictor. Genome Biol. 2016, 17, 122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Tian, D.; Wang, P.; Tang, B.; Teng, X.; Li, C.; Liu, X.; Zou, D.; Song, S.; Zhang, Z. GWAS Atlas: A Curated Resource of Genome-Wide Variant-Trait Associations in Plants and Animals. Nucleic Acids Res. 2020, 48, D927–D932. [Google Scholar] [CrossRef] [Green Version]
  43. Domingo-Fernández, D.; Hoyt, C.T.; Bobis-Álvarez, C.; Marín-Llaó, J.; Hofmann-Apitius, M. ComPath: An Ecosystem for Exploring, Analyzing, and Curating Mappings across Pathway Databases. NPJ SYST. BIOL. APPL. 2018, 4, 43. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Lonsdale, J.; Thomas, J.; Salvatore, M.; Phillips, R.; Lo, E.; Shad, S.; Hasz, R.; Walters, G.; Garcia, F.; Young, N.; et al. The Genotype-Tissue Expression (GTEx) Project. Nat. Genet. 2013, 45, 580–585. [Google Scholar] [CrossRef] [PubMed]
  45. GTEx Consortium Human Genomics. The Genotype-Tissue Expression (GTEx) Pilot Analysis: Multitissue Gene Regulation in Humans. Science 2015, 348, 648–660. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Rodríguez-López, M.L.; Martínez-Magaña, J.J.; Ruiz-Ramos, D.; García, A.R.; Gonzalez, L.; Tovilla-Zarate, C.A.; Sarmiento, E.; Juárez-Rojop, I.E.; Nicolini, H.; Gonzalez-Castro, T.B.; et al. Individuals Diagnosed with Binge-Eating Disorder Have DNA Hypomethylated Sites in Genes of the Metabolic System: A Pilot Study. Nutrients 2021, 13, 1413. [Google Scholar] [CrossRef]
  47. Manichaikul, A.; Mychaleckyj, J.C.; Rich, S.S.; Daly, K.; Sale, M.; Chen, W.-M. Robust Relationship Inference in Genome-Wide Association Studies. Bioinform. Oxf. Engl. 2010, 26, 2867–2873. [Google Scholar] [CrossRef] [Green Version]
  48. Dong, S.; Boyle, A.P. Predicting Functional Variants in Enhancer and Promoter Elements Using RegulomeDB. Hum. Mutat. 2019, 40, 1292–1298. [Google Scholar] [CrossRef] [Green Version]
  49. Zhu, Z.; Zhang, F.; Hu, H.; Bakshi, A.; Robinson, M.R.; Powell, J.E.; Montgomery, G.W.; Goddard, M.E.; Wray, N.R.; Visscher, P.M.; et al. Integration of Summary Data from GWAS and EQTL Studies Predicts Complex Trait Gene Targets. Nat. Genet. 2016, 48, 481–487. [Google Scholar] [CrossRef]
  50. Szyf, M.; Bick, J. DNA Methylation: A Mechanism for Embedding Early Life Experiences in the Genome. Child Dev. 2013, 84, 49–57. [Google Scholar] [CrossRef] [Green Version]
  51. Gibbs, J.R.; van der Brug, M.P.; Hernandez, D.G.; Traynor, B.J.; Nalls, M.A.; Lai, S.-L.; Arepalli, S.; Dillman, A.; Rafferty, I.P.; Troncoso, J.; et al. Abundant Quantitative Trait Loci Exist for DNA Methylation and Gene Expression in Human Brain. PLoS Genet. 2010, 6, e1000952. [Google Scholar] [CrossRef] [Green Version]
  52. Hannon, E.; Spiers, H.; Viana, J.; Pidsley, R.; Burrage, J.; Murphy, T.M.; Troakes, C.; Turecki, G.; O’Donovan, M.C.; Schalkwyk, L.C.; et al. Methylation Quantitative Trait Loci in the Developing Brain and Their Enrichment in Schizophrenia-Associated Genomic Regions. Nat. Neurosci. 2016, 19, 48–54. [Google Scholar] [CrossRef]
  53. Pavlides, J.M.W.; Zhu, Z.; Gratten, J.; McRae, A.F.; Wray, N.R.; Yang, J. Predicting Gene Targets from Integrative Analyses of Summary Data from GWAS and EQTL Studies for 28 Human Complex Traits. Genome Med. 2016, 8, 84. [Google Scholar] [CrossRef] [Green Version]
  54. Piazza, R.; Magistroni, V.; Redaelli, S.; Mauri, M.; Massimino, L.; Sessa, A.; Peronaci, M.; Lalowski, M.; Soliymani, R.; Mezzatesta, C.; et al. SETBP1 Induces Transcription of a Network of Development Genes by Acting as an Epigenetic Hub. Nat. Commun. 2018, 9, 2192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Yu, Q.; Zhou, Q.; Wei, Q.; Li, J.; Feng, C.; Mao, X. SEMG1 May Be the Candidate Gene for Idiopathic Asthenozoospermia. Andrologia 2014, 46, 158–166. [Google Scholar] [CrossRef] [PubMed]
  56. Shuvalov, O.; Kizenko, A.; Petukhov, A.; Fedorova, O.; Daks, A.; Bottrill, A.; Snezhkina, A.V.; Kudryavtseva, A.V.; Barlev, N. SEMG1/2 Augment Energy Metabolism of Tumor Cells. Cell Death Dis. 2020, 11, 1047. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Multi-tissue eQTL (rs10419198) results from the query to GTEx portal. NES is the slope of the linear regression, computed as the effect of the alternative allele relative to the reference allele and m-value is the posterior probability that an effect exists in a tissue.
Figure 1. Multi-tissue eQTL (rs10419198) results from the query to GTEx portal. NES is the slope of the linear regression, computed as the effect of the alternative allele relative to the reference allele and m-value is the posterior probability that an effect exists in a tissue.
Nutrients 14 00394 g001
Table 1. Overview of the characteristics of the samples.
Table 1. Overview of the characteristics of the samples.
CharacteristicMeDaCrosR (n = 168)INPRFM (n = 166)MxGDAR (n = 1469)
Age, mean (s.d)13.96 (1.94)19.32 (4.82)35.86 (15.77)
Gender
 Male, n(%)42 (0.25)16 (9.64)388 (26.41)
 Female, n(%)126 (0.75)150 (90.36)1081 (73.59)
Table 2. Associated SNPs to disordered eating in Mexican population.
Table 2. Associated SNPs to disordered eating in Mexican population.
SNPBandPositionA1/A2MAF CasesMAF ControlsORL95U95p-ValueGeneEffect
rs170301291p36.311:7059150A/G0.37870.31451.6851.3252.1412.03 × 10−5CAMTA1Intron
rs11120813 1:7062993A/G0.44140.37201.6461.3062.0752.45 × 10−5
rs6690584 1:7078434G/T0.44010.36451.7181.3592.1705.87 × 10−5
rs75212041p36.131:19138295T/C0.53290.41931.6721.3302.1001.03 × 10−5Intergenic-
rs120247381q31.11:190694813A/G0.52850.42381.5771.2671.9644.65 × 10−5LINC01720Intron
rs46269241q42.31:234909298C/T0.22600.28620.58610.45360.75734.38 × 10−5LOC107985364
rs8672862p212:45982030A/G0.45050.38211.6551.3212.0741.18 × 10−5PRKCEIntron
rs116771962p122:75830221A/G0.29490.37540.59470.46880.75461.83 × 10−5Intergenic-
rs32050602q31.12:175425346G/A0.42490.34101.6571.3182.0841.57 × 10−5WIPF13′-UTR
rs75694392q352:220590633C/T0.30910.37120.570.44720.72665.67 × 10−6Intergenic-
rs35542515 4:161798045A/C0.27330.20631.931.4722.5291.91 × 10−6
rs27489916p12.26:52596516C/T0.42340.30991.6621.3032.1184.16 × 10−5
rs38012207p14.17:42247876G/A0.54940.45061.7291.3792.1672.11 × 10−6GLI3Intron
rs3801232 7:42253313T/C0.52840.42821.7781.4122.2389.70 × 10−7
rs4724100 7:42264679C/T0.52540.43161.7261.3712.1743.50 × 10−6
rs45077688q13.38:70642018A/G0.12720.17030.50270.36350.69523.23 × 10−5SLCO5A1Intron
rs101148819q21.139:76676071T/C0.52540.42981.6281.2932.0493.34 × 10−5Intergenic-
rs1224151410p12.3110:21602923A/G0.12570.20880.45250.32930.62191.02 × 10−6
rs186502010q22.310:78688976C/T0.45660.37641.6341.3012.0522.45 × 10−5KCNMA1Intron
rs791807410q26.310:134277154A/G0.23800.15471.9221.4482.5516.20 × 10−6LOC105378569
rs10870311 10:134290526A/C0.32280.22791.7511.3472.2752.77 × 10−5Intergenic-
rs1077247112p13.212:11600364A/G0.38020.27541.661.3012.1174.58 × 10−5LOC440084Intron
rs729760612q24.312:119568596A/G0.18860.16051.9181.4152.5992.66 × 10−5SRRM4Missense (p.Ser243Asn)
rs4075945 12:119569784T/C0.18860.16091.9151.4132.5992.78 × 10−5 Intron
rs12809631 12:131045190A/C0.14670.19540.53410.39920.71462.41 × 10−5RIMBP2
rs214406714q32.3114:101952406T/C0.21560.23300.55470.41980.73303.42 × 10−5Intergenic-
rs1007904 14:101955905A/G0.23800.25890.57200.43700.74884.80 × 10−5
rs716346815q1215:26587077T/C0.20360.12431.9151.3992.6214.96 × 10−5
rs3922665 15:26590830G/A0.24250.15521.8851.4132.5141.60 × 10−5
rs804105915q21.315:58743709T/C0.27100.19991.7321.3292.2564.72 × 10−5LIPCIntron
rs1107366515q25.315:87295120G/A0.40270.32811.6261.2952.0412.78 × 10−5AGBL1
rs1713576416p13.316:2111779T/C0.24400.31440.54550.42490.70031.99 × 10−6TSC2Intron
rs1186272916p13.1216:14146098G/A0.23950.17891.8091.3642.40003.87 × 10−5Intergenic-
rs1245476318q12.318:42434615A/G0.41020.33411.6731.3282.1081.26 × 10−5SETBP1Intron
rs991014 18:42439886A/G0.40960.33491.7051.3502.1547.49 × 10−6
rs104212219q13.319:49989424C/T0.27690.35670.56770.44470.72465.46 × 10−6FLT3LGMissense (p.Phe177Leu)
rs10419198 19:50038017T/C0.30840.38330.60540.47980.76384.85 × 10−5RCN3Intron
rs607417020p12.220:10671078A/G0.41620.35011.59401.27301.99704.85 × 10−5Intergenic-
rs4813048 20:11169603T/C0.25750.18211.84401.38702.45202.55 × 10−5
rs604368420p12.120:16023836A/C0.40960.30661.63301.29802.05602.88 × 10−5MACROD2Intron
rs610408220q13.1220:43897362C/T0.34230.28271.67501.30602.14704.70 × 10−5LOC105372630
rs282400621q21.121:18099779C/T0.43860.35841.59301.27601.99004.00 × 10−5Intergenic-
rs2824065 21:18187408C/T0.39970.29681.69101.33802.13701.09 × 105
rs71330155 21:22059184A/C0.15870.22340.52760.39220.70982.38 × 10−5
Note: SNP = single-nucleotide polymorphism, Band = cytogeneitc band, Position = genomic coordinates, A1/A2 = minor allele/major allele frequency, MAF = minor allele frequency, OR = odds ratio, L95 = low 95% confidence interval, U95 = upper 95% confidence interval, Effect = in-silico variant effect prediction.
Table 3. Blood methylation quantitative loci (blood-meQTL).
Table 3. Blood methylation quantitative loci (blood-meQTL).
SNPCpGGeneLocationBetaSEp-Value
rs12024738cg12412036LOC440704TSS200−0.08690.01142.4216 × 10−14
rs12245880cg09420738 0.10410.01551.7351 × 10−11
rs12454763cg12522870SETBP1Body−0.07750.01091.6124 × 10−12
rs991014 −0.07750.01091.6124 × 10−12
rs10419198cg06378142PRR12Body−0.37680.05026.2212 × 10−14
rs2233903cg15921833SEMG1TSS15000.21240.01535.1650 × 10−44
rs6104082 0.18950.01952.3470 × 10−22
Note. SNP = single-nucleotide polymorphism, CpG = cytosine to guanine nucleotides, SE = standard error.
Table 4. Results of the summary data-based Mendelian randomization.
Table 4. Results of the summary data-based Mendelian randomization.
SNPCpGBeta SMRSE SMRp-Value SMR
rs991014cg12522870−1.98890.58526.7589 × 10−4
rs6104082cg159218330.71800.20855.7263 × 10−4
Note. SNP = single-nucleotide polymorphism, CpG = cytosine to guanine nucleotides SE = standard error.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Martínez-Magaña, J.J.; Hernandez, S.; Garcia, A.R.; Cardoso-Barajas, V.; Sarmiento, E.; Camarena, B.; Caballero, A.; Gonzalez, L.; Villatoro-Velazquez, J.A.; Medina-Mora, M.E.; et al. Genome-Wide Analysis of Disordered Eating Behavior in the Mexican Population. Nutrients 2022, 14, 394. https://0-doi-org.brum.beds.ac.uk/10.3390/nu14020394

AMA Style

Martínez-Magaña JJ, Hernandez S, Garcia AR, Cardoso-Barajas V, Sarmiento E, Camarena B, Caballero A, Gonzalez L, Villatoro-Velazquez JA, Medina-Mora ME, et al. Genome-Wide Analysis of Disordered Eating Behavior in the Mexican Population. Nutrients. 2022; 14(2):394. https://0-doi-org.brum.beds.ac.uk/10.3390/nu14020394

Chicago/Turabian Style

Martínez-Magaña, José Jaime, Sandra Hernandez, Ana Rosa Garcia, Valeria Cardoso-Barajas, Emmanuel Sarmiento, Beatriz Camarena, Alejandro Caballero, Laura Gonzalez, Jorge Ameth Villatoro-Velazquez, Maria Elena Medina-Mora, and et al. 2022. "Genome-Wide Analysis of Disordered Eating Behavior in the Mexican Population" Nutrients 14, no. 2: 394. https://0-doi-org.brum.beds.ac.uk/10.3390/nu14020394

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop