Next Article in Journal
Histone Acetylation Changes in Plant Response to Drought Stress
Next Article in Special Issue
Advances and Limitations of Next Generation Sequencing in Animal Diet Analysis
Previous Article in Journal
Common and Unique Genetic Background between Attention-Deficit/Hyperactivity Disorder and Excessive Body Weight
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 12
Issue 9
10.3390/genes12091408
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Years of Schooling Could Reduce Epigenetic Aging: A Study of a Mexican Cohort

by
Juan Carlos Gomez-Verjan
1,*,
Marcelino Esparza-Aguilar
2,
Verónica Martín-Martín
3,
Cecilia Salazar-Perez
2,
Cinthya Cadena-Trejo
1,
Luis Miguel Gutierrez-Robledo
4,
José Jaime Martínez-Magaña
5,
Humberto Nicolini
5 and
Pedro Arroyo
1
1
Dirección de Investigación, Instituto Nacional de Geriatría, Mexico City 10200, Mexico
2
Departamento de Investigación en Epidemiología, Instituto Nacional de Pediatría, Mexico City 04530, Mexico
3
Subdirección de Investigación Médica, Instituto Nacional de Pediatría, Mexico City 04530, Mexico
4
Dirección General, Instituto Nacional de Geriatría, Mexico City 10200, Mexico
5
Laboratorio de Genómica de Enfermedades Psiquiátricas y Neurodegenerativas, Instituto Nacional de Medicina Genómica, Mexico City 04809, Mexico
*
Author to whom correspondence should be addressed.
Genes 2021, 12(9), 1408; https://doi.org/10.3390/genes12091408
Submission received: 25 May 2021 / Revised: 27 August 2021 / Accepted: 3 September 2021 / Published: 13 September 2021
(This article belongs to the Special Issue Nutrigenomics and Cellular Metabolism)
Browse Figures
Review Reports Versions Notes

Abstract

:
Adverse conditions in early life, including environmental, biological and social influences, are risk factors for ill-health during aging and the onset of age-related disorders. In this context, the recent field of social epigenetics offers a valuable method for establishing the relationships among them However, current clinical studies on environmental changes and lifespan disorders are limited. In this sense, the Tlaltizapan (Mexico) cohort, who 52 years ago was exposed to infant malnutrition, low income and poor hygiene conditions, represents a vital source for exploring such factors. Therefore, in the present study, 52 years later, we aimed to explore differences in clinical/biochemical/anthropometric and epigenetic (DNA methylation) variables between individuals from such a cohort, in comparison with an urban-raised sample. Interestingly, only cholesterol levels showed significant differences between the cohorts. On the other hand, individuals from the Tlaltizapan cohort with more years of schooling had a lower epigenetic age in the Horvath (p-value = 0.0225) and PhenoAge (p-value = 0.0353) clocks, compared to those with lower-level schooling. Our analysis indicates 12 differentially methylated sites associated with the PI3-Akt signaling pathway and galactose metabolism in individuals with different durations of schooling. In conclusion, our results suggest that longer durations of schooling could promote DNA methylation changes that may reduce epigenetic age; nevertheless, further studies are needed.

1. Introduction

Biological and social environmental factors influence the development and life expectancy of the average human being [1,2]. The current investigations are centered on the impact of different environmental exposures at different stages of life, including those involved in aging and the manifestation of age-related diseases [3,4,5,6,7]. In this sense, studies have identified a crosslink between early-life adverse conditions and the reduction of life expectancy and “healthspan” [8]. For instance, suffering rheumatic fever and malaria in childhood is associated with a higher prevalence of cardiac disease and increased mortality [9]. Poor socioeconomic status early in life is also predictive of later cognition problems and disability [10]. Another study performed by the Japan Gerontological Evaluation Study Cohort found an association between adverse childhood experiences (parental death, divorce, parental mental illness, family violence, physical abuse, psychological neglect, and psychological abuse) and a greater risk of dementia later in life [11]. In this context, an analysis of the effect of early-life conditions on the disease’s manifestation has two general hypotheses. The antagonistic pleiotropy states that beneficial traits early in life could be detrimental later, due to a decline in the force of natural selection [12]. Furthermore, the Barker hypothesis states that inadequate nutrition in utero may induce an adaptive metabolic program in the fetus that could lead to future disease [13,14]. Some of the most compelling and unfortunate examples of fetal or child detrimental nutrition are historical famines, such as the Chinese Famine or the Hunger Winter [14,15,16,17]. Individuals from those cohorts who were exposed to nutritional deficiency in childhood had a higher risk of age-related illness, such as neurodegenerative diseases, and a meaningful reduction in life expectancy [18,19,20,21].

1.1. Social Epigenetics

As mentioned above, as we grow, we are exposed to different adverse environmental stimuli; however, changes in our lifestyle could modify the risk of disease [22,23,24,25]. In this context, the recent field of social epigenetics seems valuable in creating a network of influences to understand the role of the environment and the further development of diseases [26,27,28,29]. Epigenetics mechanisms are genome dynamics that change with time and environmental conditions [30,31]. The most studied epigenetic mechanism in molecular epidemiology concerns analyzing the genome’s DNA methylation patterns (DNAm). This possesses a strong correlation with different diseases, such as metabolic disorders, neurodegenerative diseases, and aging.
Moreover, different research groups worldwide have developed epigenetic clocks, based on data from DNAm patterns [32,33,34,35,36,37], which have proven to be valuable biomarkers for mortality and to be associated with several age-related diseases, such as cancer, diabetes, cardiovascular disease, and frailty, among others [38,39,40]. Although there is little knowledge available on the effect that changes in lifestyle or social environment could have on epigenetic mechanisms, changes in our environment by lifestyle modification, like the number of years of schooling, physical activity, or diet habits, are linked to changes in aging and age-related disorders [32,35].

1.2. Tlaltizapan Cohort

In Mexico, between 1966 and 1967, Dr. Joaquin Cravioto recruited all the mothers and children born in a rural town in Tlaltizapan, Morelos, Mexico for study. These individuals were exposed to early adverse conditions, such as low birth weight, low income, poor hygiene conditions at birth and in their homes, and malnutrition. The parents were mainly employed in agriculture; therefore, they were usually exposed to organochlorine-organophosphorus pesticides [41,42]. Such individuals were named the Tlaltizapan cohort and represent an essential source for exploring changes in environmental exposures, not only for the volume of information on them but also because they were born under adverse early-life conditions in a rural community in Mexico. Therefore, in the present study, we aimed to explore the differences in clinical/biochemical and anthropometric variables and DNAm variations between individuals from the Tlaltizapan cohort, 52 years later, and compare them with an urban-raised sample.

2. Materials and Methods

2.1. Sample Population

Cravioto et al. recruited all children born in around one year (February 1966–February 1967) in Tlaltizapan, with a population of 5637 people. A detailed description of the rationale and methodology of the original study was published in a 1969 Monograph [41,43]. In 2018, we undertook a search of all members of the cohort. Out of 336 births, four were stillborn, 26 died in the first year of life, six died before the age of five, and 16 died before reaching adulthood. The Infant Death Rate was 78.3 per thousand live births. Of the remaining 284 cohort members, 40 were residents of the USA. Ten migrated to the other Mexican States, and we were unable to locate 116. It was possible to track down 116 residents of Tlaltizapan; 36 of them refused to participate in the study, but 82 (50 women and 32 men) gave informed consent to be clinically evaluated at the Instituto Nacional de Geriatría (INGER) in Mexico City. A subgroup comprising 32 birth cohort members was selected for the epigenetic study, based on information regarding family characteristics and living conditions at birth. Additionally, seven healthy middle-class older participants who were born and lived in México City were included as urban-raised individuals. All subjects freely agreed to participate and signed written informed consent forms. The protocol was reviewed and approved by the Bioethics and Research Board of the National Institute of Geriatrics (INGER) and Bioethics and Research Board of the National Institute of Pediatrics (INP) under the number INP-INGER 06/2018. All methods were performed based on the Helsinki declaration and the local Ley General de Salud of Mexico. The information was also protected under the Ley Federal de Proteccion de Datos Personales of Mexico.

2.2. Clinical and Anthropometric Measurements

2.2.1. Biochemical Profile

We performed the following standard hematological tests (hemoglobin (g/dL), the percentage of hematocrit, transferring (mg/dL) and the percentage of transferrin saturation, according to Sandler [44], and performed the following standard metabolic tests for glucose (mg/dL) and hemoglobin A1C, according to the World Health Organization (WHO) recommendations [45], and triglycerides (mg/dL), LDL (mg/dL), HDL (mg/dL), cholesterol (mg/dL), creatinine (mg/dL), and uric acid (mg/dL) levels, according to Muñoz et al. [46].

2.2.2. Anthropometric Measurements

We evaluated the following anthropometric measurements: weight (kg), height (cm), knee height (cm), waist circumference (cm), leg circumference (cm), body mass index (BMI—kg/m2 height) and percentage of oxygenation, according to the WHO’s recommendations [47].

2.2.3. Body Composition and Physical Performance

For this parameter, we evaluated the appendicular muscular mass (kg/m2) and the bone mineral density (T-value) using dual X-ray densitometry (DEXA) [48]. Additionally, we evaluate the short physical performance battery (SPPB): tandem movement (s), unipodal movement (s), walking speed (cm/s), rising from a chair (times/30 s), and maximum grip strength (kg) [49].

2.2.4. Health Status

We evaluated the number of comorbidities, and performed a mini-mental evaluation [50], with depression scale, fear of falls scale [51], frailty scale, falls history, and blood pressure (mmHg).

2.3. DNA Methylation

Before clinical and anthropometric evaluation, an EDTA blood sample was collected. According to the manufacturer protocols, DNA was extracted from peripheral blood mononuclear cells using standard methods (Qiagen LLC, Germantown, MD, USA, DNA extraction kit) and was bisulfite-converted using a commercial EZ DNA methylation kit (Zymo Research, Irvine, CA, USA). Then bisulfite converted DNA was processed using the commercial Infinium Methylation EPIC BeadChip (Illumina, San Diego, CA, USA) on the microarray platform located at the Unidad de Microarreglos of the Instituto Nacional de Medicina Genomica (INMEGEN, Mexico City, Mexico), following the manufacturer’s protocol. The EPIC contains 866, 836 CpG sites; this EPIC array has a similar Infinium assay design as the previous 450K array, but it has twice as much coverage of DNA methylation sites (known as CpG sites), which are particularly enriched for promoter and enhancer regulatory regions. Such characteristics provide an increased power and genome coverage to identify novel loci that are relevant for aging and health disparities. The EPIC intensity signals were transformed to idat files by the GenomeStudio software (Version 2.0) (Illumina, CA, USA). Then, the idat files were processed for quality control under the ChAMP package (version 2.22.0). The quality control included the remotion of: (i) 8,250 CpG site with detection p > 0.01, (ii) 2,727 < 3 beads in 5% samples, (iii) 2,969 NoCpGs, (iv) 95,778 SNPs [52], (v) 50 multihit [53,54], and (vi) 16,549 CpGs sites found on X and Y chromosomes. We included 12 technical replicates (two per plate) for quality control. We removed those probes with low intensity, as well as those that could introduce biases in performance.
Additionally, for this analysis, we removed the probes for sexual chromosomes so that all subjects could be more easily compared since sexual chromosome DNA methylation at baseline is inherently more variable in females than in males. After quality control, a total of 740,513 CpG sites remained for the subsequent analysis. Then, a matrix of β-values and M-values was generated.

2.3.1. Epigenetic Clocks

Epigenetic clocks were calculated using the ENmix package (version 1.28.4) after the ChAMP β-mixture quantile normalization (BMIQ) .norm function of the ChAMP package [55], using the β matrix. The Hannum, Horvath, and DNA PhenoAge clock values were calculated. The age-acceleration level was also calculated, regressing the epigenetic age calculated by each clock with the chronological age [37]. Pairwise correlations were also calculated between each clock and the chronological age (for example, chronological age was correlated with PhenoAge). The pairwise correlations were performed in both cohorts (Tlatizapan and urban-raised). Pearson correlations were performed, and a p-value < 0.05 was considered statistically significant.

2.4. Statistical Analysis

2.4.1. Clinical Study, Epigenetic Clocks, and Anthropometric Comparison between the Tlatizapan Cohort and the Urban-Raised Cohort

Continuous variables were represented by means and standard deviations. Meanwhile, categorical variables were represented as both a total and a percentage. The Shapiro–Wilk test evaluated the normal distribution of the variables. All the variables followed the normal distribution. Continuous variables were compared by t-test. T-test with Welch correction, and categorical by chi-squared test. A p-value < 0.05 was considered statistically significant. These analyses were performed using the R-statistical language (version 3.6.3) [56].

2.4.2. Analysis Only on the Tlatizapan Cohort

Epigenetic Clocks

To explore the factors associated with epigenetic aging in the Tlatizapan cohort, we separated the individuals into groups, based on the age acceleration estimated in Section 2.3.1. Individual analyses were performed for each clock (PhenoAge, Horvath, Hannum). Participants were divided based on the mean of their epigenetic age acceleration, i.e., individuals with accelerated aging (those with a higher than zero value for age acceleration) and non-accelerated individuals (those with a value lower than or equal to zero value on age acceleration). Next, we compared all the anthropometric and clinical variables between the groups, as described in Section 2.4.1.

Epigenome-Wide Association Analysis

Epigenome-wide analysis (differentially methylated sites) of DNA methylation was evaluated by linear models implemented in the limma package (version 3.48.2) [57]. In this analysis, we compared only the Tlatizapan cohort between individuals with high (n = 15) and low duration of schooling (n = 17). The cut-off of years of schooling was 9 years (Mexico national mean), with high: ≥ 9 years and low: < 9 years. Single-value deconvolution analysis was performed to explore the confounding variables (all the clinical and anthropometric variables explored in the Tlatizapan cohort) with the champ.SVD function; no variable was found to be correlated. Batch effects, by slide and array position, were removed by the function champ.runCombat. The effect of blood-cell proportion was removed by adjusting the β-values matrix by a reference-based method implemented in the champ.refbase function of the CHAMP package. For statistical contrast, the β-values were transformed to M-values (better adjustment for linear models), and for reporting the differences in methylation levels, the β-values were used (since they were more biologically interpretable) [58]. A p-value < 9 × 10−8 was considered statistically significant because this p-value threshold adequately controls the false positive rate on the EPIC array [59]. After statistical contrast, we annotated the CpG sites to the nearest gene using the IlluminaHumanMethylationEPICanno.ilm10b2.hg19 package (version 0.6.0). Gene enrichment analysis was performed using the online tool Webgestalt, using the KEGG pathway database, and an FDR-corrected p-value < 0.05 was considered statistically significant [60,61].

2.5. Bioethical Considerations

All participating subjects provided written informed consent regarding the study, genetic data, and personal information security. The institutional review boards approved enrollment and consent procedures for this study, given by the Bioethics and Research Board of the National Institute of Geriatrics (INGER) and the Bioethics and Research Board of the National Institute of Pediatrics (INP) under the number INP-INGER 06/2018. All methods employed in the study were performed, following the relevant international guidelines and regulations as well as those from the Ley General de Salud of Mexico. The data are available upon express request addressed to the corresponding author and are currently in safekeeping by INGER, with the approval of the Research and Ethics Committees.

3. Results

3.1. Differences between the Urban-Raised and the Tlatizapan Cohorts

In the analysis of clinical variables between individuals who were urban-raised and individuals in the Tlaltizapan cohort, we found differences only in the serum total cholesterol level (Table 1).

Evaluation of Different Epigenetic Clocks

The lowest correlation with chronological age was found with the Hannum clock. Chronological age was positively correlated with Horvath and negatively correlated with PhenoAge (Figure 1). In the correlation between epigenetic clocks, a higher correlation was between Horvath and Hannum. PhenoAge had a lower correlation with Hannum and Horvath. The epigenetic age of the urban raised was higher than those of the Tlatizapan cohort; nevertheless, this difference could be derived by the higher chronological age of the urban raised individuals.

3.2. Analysis of the Tlatizapan Cohort

3.2.1. Effect of Clinical Variables on Individuals with Accelerated Epigenetic Aging

After the evaluation comparison with the urban-raised individuals, we performed a case-only stratified analysis in the Tlatizapan cohort. The individuals were divided into groups based on the mean epigenetic age acceleration for each clock (i.e., a separated analysis for each epigenetic clock), on individuals with accelerated aging (higher than zero on age acceleration), and non-accelerated aging (lower or equal than zero). In this analysis, we only found differences between the number of years of schooling in the groups constructed with Horvath (t = −2.4070, p-value = 0.0225) and PhenoAge (t = −2.2110, p-value = 0.0353), and no differences were found in Hannum (t = −0.8793, p-value = 0.3863). The individuals with no epigenetic age acceleration had longer schooling than those with epigenetic age acceleration (Figure 2). In addition, we found a correlation between the duration of schooling and age acceleration (r = −0.34, p-value = 0.0452). We also found that individuals of the Tlatizapan cohort with more years of schooling (>9 years) had a reduction in Horvath clock epigenetic aging (mean = −1.88) compared to the urban-raised individuals (mean = 4.77) (t = −2.7684, p-value = 0.0273).

3.2.2. Evaluation of Epigenetic Changes in Individuals with Long and Short Durations of Schooling

Based on the differences found previously, we divided the Tlaltizapan cohort based on the duration of schooling—high (n = 15) and low (n = 17) number of years in education (high: ≥9 years, low: <9 years). In the analysis of differentially methylated sites between individuals, we did not find any statistically significant signal. Nevertheless, we focused on the 12 top CpG sites (p-value < 5 × 10−5) (Table 2; Supplementary Table S1).
The CpG sites that were nominally differentially methylated were found in 6 intergenic regions, five protein-coding regions, and one long non-protein-coding region. Most of the sites (n = 8, 66.67%) were found in open-sea regions. Cg03184819 had a higher difference between groups (Log FC = −0.0966), followed by the CpG site found in the transcription starting site of the lactase gene (Log FC = −0.0892); both sites were hypermethylated in the individuals with a higher number of years of schooling. The genes with differentially methylated sites were enriched in the PI3-Akt signaling pathway (hsa04151, adjusted p-value = 0.0082, MAGI1 and TNC) and galactose metabolism (hsa00052, adjusted p-value = 0.0137, LCT).

4. Discussion

Social epigenetics is a helpful tool for analyzing the effect of exposures to adverse environments during the first years of life. Hence, in the present report, through such a tool, we analyze a subsample that was part of the Tlaltizapan cohort, which was reported as being exposed to early adverse conditions. Interestingly, according to our results, PhenoAge is not positively correlated with chronological age (in comparison to Horvath’s clock); this could be because this biological clock is more associated with pro-inflammatory and DNA damage routes that are perhaps not yet altered in this population, due to the fact that they are still only in their early fifties.
On the other hand, our findings suggest that individuals with age acceleration in the Tlaltizapan cohort had fewer years of schooling. Such results must be interpreted with caution since correlation is not causation. Nevertheless, this is not the first time that education has been correlated with methylation, and several studies also indicated correlations between methylation patterns and education [32,34,35].
Furthermore, McGuiness et al. [62] suggest that in a cohort of 239 men and women from Scotland, each additional year of education was associated with a 2.4% greater global DNA methylation content. Moreover, van Dongen et al. [32] identified 58 CpG sites associated with educational attainment, and that the methylome level of poorly educated people resembles that of smokers. Additionally, Perng et al. [63] found, in a cross-sectional study of 568 children from Bogota, Colombia that changes in global DNA methylation are associated with different concentrations of micronutrients.
The differences in epigenome-wide analysis (influenced by advanced years of schooling) point to the galactose metabolism’s inactivation. D-galactose’s presence has been a clear indicator of aging by promoting mitochondrial dysfunction and cognitive decline; the mechanism of such changes remains unknown, but epigenetics inactivation could be a fundamental point. On the other hand, the PI3K-Akt pathway was enriched among the associated CpG sites, and this pathway is involved in aging and Alzheimer’s disease [64]. Moreover, this pathway targets the modulation of aging, due to its involvement in skin, muscular, and vascular aging [65].
Among our analysis of epigenome-wide differences between individuals with high or low levels of schooling, the genes with nominal associations were LCT, ASB18, MAGI1, ANK2, TNC, LOC100128811. One of the most interesting findings regarding nutritional status is lactase (LCT). A cross-sectional study with 196 malnourished children from Uganda demonstrated that nutritional injury reduces the capacity of the intestinal mucosa to synthesize LCT [66]. Moreover, the genetic expression of lactase was reduced in a study with biopsy specimens from 29 malnourished infants [67].
Although our results could generate novel hypotheses for a possible intervention regarding aging, they must be assessed with caution since there are limitations on our study that must be pointed out. For instance, the sample size reduces the statistical power of our analysis, and the age of the individuals is still too young to indicate age-related diseases (52 years old); additionally, more comparisons with a significantly larger number of urban-raised individuals must be carried out. Nevertheless, as one strength, our study carefully characterizes each variable analyzed. Additionally, ancestry is another interesting point to highlight; the analyzed cohort is Mexican-raised. In this sense, Horvath et al. [68] reported that epigenetic age could be different based on ancestry and that estimations on Latin American and Hispanic populations were needed, so our analysis adds information to this particular field.

5. Conclusions

In the present study, we subsampled a historical cohort who were recruited more than 52 years ago and who were exposed to several adverse conditions in their early development and in infancy, such as low income, low birth weight, malnutrition, and poor hygienic status. Our results indicate that more years of schooling could reduce epigenetic aging and that PI3K-Akt and galactose metabolism are epigenetically altered. It is noted that further studies with larger sample sizes are needed to validate our results.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12091408/s1. Supplementary Table S1 The CpG site association.

Author Contributions

Funding Acquisition: J.C.G.-V., M.E.-A., V.M.-M., L.M.G.-R., P.A. Methodology: J.C.G.-V., M.E.-A., V.M.-M., C.S.-P., C.C.-T., L.M.G.-R., P.A. Data curation: M.E.-A., V.M.-M., J.C.G.-V. Bioinformatic analysis: J.C.G.-V., J.J.M.-M. Writing-original draft: J.C.G.-V., J.J.M.-M., H.N. Writing—review and editing: J.C.G.-V., M.E.-A., L.M.G.-R., P.A. All authors have read and agreed to the published version of the manuscript.

Funding

This project was supported by CONACYT-FOSISS 2017 number 290287 under the name: “Condiciones sociales y de salud al Nacimiento, primera infancia y su relación con las condiciones en la edad adulta y la expresión genómica, como predictores de envejecimiento saludable desde los 50 años de edad”. The publication of this paper was supported by a grant from the Secretaría de Educación, Ciencia, Tecnología e Innovación de la Ciudad de México CM-SECTEI/200/2020 “Red Colaborativa de Investigación Traslacional para el Envejecimiento Saludable de la Ciudad de México” (RECITES).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of the National Institute of Geriatrics (INGER) and Bioethics and Research board of the National Institute of Pediatrics (INP) (protocol code INP-INGER 06/2018, date = 06/2018).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. All subjects freely agree to participate in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to ethical restrictions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gluckman, P.D.; Hanson, M.A.; Spencer, H.G.; Bateson, P. Environmental Influences during Development and Their Later Consequences for Health and Disease: Implications for the Interpretation of Empirical Studies. Proc. R. Soc. B Biol. Sci. 2005, 272, 671–677. [Google Scholar] [CrossRef] [Green Version]
  2. Schell, L.M.; Gallo, M.V.; Ravenscroft, J. Environmental Influences on Human Growth and Development: Historical Review and Case Study of Contemporary Influences. Ann. Hum. Biol. 2009, 36, 459–477. [Google Scholar] [CrossRef]
  3. Bektas, A.; Schurman, S.H.; Sen, R.; Ferrucci, L. Aging, Inflammation and the Environment. Exp. Gerontol. 2018, 105, 10–18. [Google Scholar] [CrossRef]
  4. Christensen, K.; McGue, M. Healthy Ageing, the Genome and the Environment. Nat. Rev. Endocrinol. 2016, 12, 378–380. [Google Scholar] [CrossRef] [Green Version]
  5. Lovely, C.; Rampersad, M.; Fernandes, Y.; Eberhart, J. Gene-Environment Interactions in Development and Disease: Gene-Environment Interactions in Development and Disease. Wiley Interdiscip. Rev. Dev. Biol. 2017, 6, e247. [Google Scholar] [CrossRef]
  6. Vrieze, S.I.; Iacono, W.G.; McGue, M. Confluence of Genes, Environment, Development, and Behavior in a Post Genome-Wide Association Study World. Dev. Psychopathol. 2012, 24, 1195–1214. [Google Scholar] [CrossRef] [Green Version]
  7. Strife, S.; Downey, L. Childhood Development and Access to Nature: A New Direction for Environmental Inequality Research. Organ. Environ. 2009, 22, 99–122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Gomez-Verjan, J.C.; Barrera-Vázquez, O.S.; García-Velázquez, L.; Samper-Ternent, R.; Arroyo, P. Epigenetic Variations Due to Nutritional Status in Early-life and Its Later Impact on Aging and Disease. Clin. Genet. 2020, 98, 313–321. [Google Scholar] [CrossRef] [PubMed]
  9. McEniry, M. Early-Life Conditions and Older Adult Health in Low- and Middle-Income Countries: A Review. J. Dev. Orig. Health Dis. 2013, 4, 10–29. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Britton, A.; Shipley, M.; Singh-Manoux, A.; Marmot, M.G. Successful Aging: The Contribution of Early-Life and Midlife Risk Factors: SUCCESSFUL AGING. J. Am. Geriatr. Soc. 2008, 56, 1098–1105. [Google Scholar] [CrossRef] [Green Version]
  11. Tani, Y.; Fujiwara, T.; Kondo, K. Adverse Childhood Experiences and Dementia: Interactions With Social Capital in the Japan Gerontological Evaluation Study Cohort. Am. J. Prev. Med. 2021, 61, 225–234. [Google Scholar] [CrossRef] [PubMed]
  12. Flatt, T. Survival Costs of Reproduction in Drosophila. Exp. Gerontol. 2011, 46, 369–375. [Google Scholar] [CrossRef] [Green Version]
  13. Hales, C.N.; Barker, D.J.P. Type 2 (Non-Insulin-Dependent) Diabetes Mellitus: The Thrifty Phenotype Hypothesis. Diabetologia 1992, 35, 595–601. [Google Scholar] [CrossRef] [PubMed]
  14. Almond, D.; Currie, J. Killing Me Softly: The Fetal Origins Hypothesis. J. Econ. Perspect. 2011, 25, 153–172. [Google Scholar] [CrossRef] [Green Version]
  15. Tobi, E.W.; Goeman, J.J.; Monajemi, R.; Gu, H.; Putter, H.; Zhang, Y.; Slieker, R.C.; Stok, A.P.; Thijssen, P.E.; Müller, F.; et al. DNA Methylation Signatures Link Prenatal Famine Exposure to Growth and Metabolism. Nat. Commun. 2014, 5, 5592. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Jowett, A.J. The Growth of China’s Population, 1949–1982 (With Special Reference to the Demographic Disaster of 1960–61). Geogr. J. 1984, 150, 155. [Google Scholar] [CrossRef]
  17. Jowett, A.J. The Demographic Responses to Famine: The Case of China 1958? 61. GeoJournal 1991, 23. [Google Scholar] [CrossRef]
  18. Ahmed, F. Epigenetics: Tales of Adversity. Nature 2010, 468, S20. [Google Scholar] [CrossRef] [Green Version]
  19. Schulz, L.C. The Dutch Hunger Winter and the Developmental Origins of Health and Disease. Proc. Natl. Acad. Sci. USA 2010, 107, 16757–16758. [Google Scholar] [CrossRef] [Green Version]
  20. Gensous, N.; Franceschi, C.; Santoro, A.; Milazzo, M.; Garagnani, P.; Bacalini, M.G. The Impact of Caloric Restriction on the Epigenetic Signatures of Aging. Int. J. Mol. Sci. 2019, 20, 2022. [Google Scholar] [CrossRef] [Green Version]
  21. Wang, Z.; Song, J.; Li, Y.; Dong, B.; Zou, Z.; Ma, J. Early-Life Exposure to the Chinese Famine Is Associated with Higher Methylation Level in the INSR Gene in Later Adulthood. Sci. Rep. 2019, 9, 3354. [Google Scholar] [CrossRef]
  22. Cohen, D.A.; Babey, S. Reducing Environmental Risks to Modify Lifestyle Behaviors. Diabetes Spectr. 2015, 28, 2–4. [Google Scholar] [CrossRef] [Green Version]
  23. Anand, P.; Kunnumakara, A.B.; Sundaram, C.; Harikumar, K.B.; Tharakan, S.T.; Lai, O.S.; Sung, B.; Aggarwal, B.B. Cancer Is a Preventable Disease That Requires Major Lifestyle Changes. Pharm. Res. 2008, 25, 2097–2116. [Google Scholar] [CrossRef] [PubMed]
  24. Kolb, H.; Martin, S. Environmental/Lifestyle Factors in the Pathogenesis and Prevention of Type 2 Diabetes. BMC Med. 2017, 15. [Google Scholar] [CrossRef] [PubMed]
  25. Kremers, S.P.J.; Eves, F.F.; Andersen, R.E. Environmental Changes to Promote Physical Activity and Healthy Dietary Behavior. J. Environ. Public Health 2012, 2012, 470858. [Google Scholar] [CrossRef]
  26. Loi, M.; Del Savio, L.; Stupka, E. Social Epigenetics and Equality of Opportunity. Public Health Ethics 2013, 6, 142–153. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Merrill, S.M.; Gladish, N.; Kobor, M.S. Social Environment and Epigenetics. In Behavioral Neurogenomics; Binder, E.B., Klengel, T., Eds.; Springer International Publishing: Cham, Switzerland, 2019; Volume 42, pp. 83–126. ISBN 978-3-030-31264-0. [Google Scholar]
  28. Cunliffe, V.T. The Epigenetic Impacts of Social Stress: How Does Social Adversity Become Biologically Embedded? Epigenomics 2016, 8, 1653–1669. [Google Scholar] [CrossRef] [Green Version]
  29. Curley, J.P.; Jensen, C.L.; Mashoodh, R.; Champagne, F.A. Social Influences on Neurobiology and Behavior: Epigenetic Effects during Development. Psychoneuroendocrinology 2011, 36, 352–371. [Google Scholar] [CrossRef] [Green Version]
  30. Moosavi, A.; Motevalizadeh Ardekani, A. Role of Epigenetics in Biology and Human Diseases. Iran. Biomed. J. 2016, 20, 246–258. [Google Scholar] [CrossRef] [PubMed]
  31. Barouki, R.; Melén, E.; Herceg, Z.; Beckers, J.; Chen, J.; Karagas, M.; Puga, A.; Xia, Y.; Chadwick, L.; Yan, W.; et al. Epigenetics as a Mechanism Linking Developmental Exposures to Long-Term Toxicity. Environ. Int. 2018, 114, 77–86. [Google Scholar] [CrossRef]
  32. BIOS Consortium; van Dongen, J.; Bonder, M.J.; Dekkers, K.F.; Nivard, M.G.; van Iterson, M.; Willemsen, G.; Beekman, M.; van der Spek, A.; van Meurs, J.B.J.; et al. DNA Methylation Signatures of Educational Attainment. NPJ Sci. Learn. 2018, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. BIOS Consortium; Karlsson Linnér, R.; Marioni, R.E.; Rietveld, C.A.; Simpkin, A.J.; Davies, N.M.; Watanabe, K.; Armstrong, N.J.; Auro, K.; Baumbach, C.; et al. An Epigenome-Wide Association Study Meta-Analysis of Educational Attainment. Mol. Psychiatry 2017, 22, 1680–1690. [Google Scholar] [CrossRef] [Green Version]
  34. Pickersgill, M. Epigenetics, Education, and the Plastic Body: Changing Concepts and New Engagements. Res. Educ. 2020, 107, 72–83. [Google Scholar] [CrossRef] [Green Version]
  35. Chen, H.; Hu, H. The Relationship and Mechanism between Education and Functional Health Status Transition among Older Persons in China. BMC Geriatr. 2018, 18, 89. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Hannum, G.; Guinney, J.; Zhao, L.; Zhang, L.; Hughes, G.; Sadda, S.; Klotzle, B.; Bibikova, M.; Fan, J.-B.; Gao, Y.; et al. Genome-Wide Methylation Profiles Reveal Quantitative Views of Human Aging Rates. Mol. Cell 2013, 49, 359–367. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Horvath, S. DNA Methylation Age of Human Tissues and Cell Types. Genome Biol. 2013, 14, R115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Lu, A.T.; Quach, A.; Wilson, J.G.; Reiner, A.P.; Aviv, A.; Raj, K.; Hou, L.; Baccarelli, A.A.; Li, Y.; Stewart, J.D.; et al. DNA Methylation GrimAge Strongly Predicts Lifespan and Healthspan. Aging 2019, 11, 303–327. [Google Scholar] [CrossRef]
  39. Levine, M.E.; Lu, A.T.; Quach, A.; Chen, B.H.; Assimes, T.L.; Bandinelli, S.; Hou, L.; Baccarelli, A.A.; Stewart, J.D.; Li, Y.; et al. An Epigenetic Biomarker of Aging for Lifespan and Healthspan. Aging 2018, 10, 573–591. [Google Scholar] [CrossRef] [Green Version]
  40. Levine, M.E. Assessment of Epigenetic Clocks as Biomarkers of Aging in Basic and Population Research. J. Gerontol. A Biol. Sci. Med. Sci. 2020, 75, 463–465. [Google Scholar] [CrossRef]
  41. Cravioto, J.; Birch, H.; De Licardie, E.; Rosales, L.; Vega, L. The Ecology of Growth and Development in a Mexican Preindustrial Community Report 1: Method and Findings from Birth to One Month of Age. Monogr. Soc. Res. Child Dev. 1969, 34, 1–76. [Google Scholar] [CrossRef]
  42. Galvan-Ramos, R.; Cravioto, J.; Gutierrez, G.; Gomez, F.; Frenk, S. Operation Zacatepec. III. Comparison of an Indirect and a Direct Method in Evaluation of the Nutritional Status of Children of a Rural Community. Bol. Méd. Hosp. Infant. México 1958, 15, 855–873. [Google Scholar]
  43. Vega Franco, L.; Cravioto, J.; Roca de Licardie, E. Ecology of the Size of Newborn Infants in a Preindustrial Community. Bol. Méd. Hosp. Infant. México 1974, 31, 199–206. [Google Scholar]
  44. Sandler, S.G. Management of Patients with Hematologic Malignancies and Aplastic Anemia Who Are Refractory to Platelet Transfusions. Haematol. Budap 1998, 29, 1–11. [Google Scholar]
  45. World Health Organization. Prevention of Diabetes Mellitus. Report of a WHO Study Group; World Health Organization: Geneva, Switzerland, 1994; Volume 844, pp. 1–100. [Google Scholar]
  46. Muñoz, O.M.; García, Á.A.; Fernández-Ávila, D.; Higuera, A.; Ruiz, Á.J.; Aschner, P.; Toro, J.M.; Arteaga, J.M.; Merchán, A.; Sánchez, G.; et al. Guía de práctica clínica para la prevención, detección temprana, diagnóstico, tratamiento y seguimiento de las dislipidemias: Evaluación del riesgo cardiovascular. Rev. Colomb. Cardiol. 2015, 22, 263–269. [Google Scholar] [CrossRef] [Green Version]
  47. World Health Organization. Physical Status: The Use and Interpretation of Anthropometry. Report of a WHO Expert Committee; World Health Organization: Geneva, Switzerland, 1995; Volume 854, pp. 1–452. [Google Scholar]
  48. Shepherd, J.A.; Ng, B.K.; Sommer, M.J.; Heymsfield, S.B. Body Composition by DXA. Bone 2017, 104, 101–105. [Google Scholar] [CrossRef]
  49. Guralnik, J.M.; Simonsick, E.M.; Ferrucci, L.; Glynn, R.J.; Berkman, L.F.; Blazer, D.G.; Scherr, P.A.; Wallace, R.B. A Short Physical Performance Battery Assessing Lower Extremity Function: Association With Self-Reported Disability and Prediction of Mortality and Nursing Home Admission. J. Gerontol. 1994, 49, M85–M94. [Google Scholar] [CrossRef]
  50. Folstein, M.F.; Folstein, S.E.; McHugh, P.R. “Mini-Mental State”. A Practical Method for Grading the Cognitive State of Patients for the Clinician. J. Psychiatr. Res. 1975, 12, 189–198. [Google Scholar] [CrossRef]
  51. Delbaere, K.; Close, J.C.T.; Mikolaizak, A.S.; Sachdev, P.S.; Brodaty, H.; Lord, S.R. The Falls Efficacy Scale International (FES-I). A Comprehensive Longitudinal Validation Study. Age Ageing 2010, 39, 210–216. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Zhou, W.; Laird, P.W.; Shen, H. Comprehensive Characterization, Annotation and Innovative Use of Infinium DNA Methylation BeadChip Probes. Nucleic Acids Res. 2016, 45, e22. [Google Scholar] [CrossRef]
  53. Nordlund, J.; Bäcklin, C.L.; Wahlberg, P.; Busche, S.; Berglund, E.C.; Eloranta, M.-L.; Flaegstad, T.; Forestier, E.; Frost, B.-M.; Harila-Saari, A.; et al. Genome-Wide Signatures of Differential DNA Methylation in Pediatric Acute Lymphoblastic Leukemia. Genome Biol. 2013, 14, r105. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Pidsley, R.; Zotenko, E.; Peters, T.J.; Lawrence, M.G.; Risbridger, G.P.; Molloy, P.; Van Djik, S.; Muhlhausler, B.; Stirzaker, C.; Clark, S.J. Critical Evaluation of the Illumina MethylationEPIC BeadChip Microarray for Whole-Genome DNA Methylation Profiling. Genome Biol. 2016, 17, 208. [Google Scholar] [CrossRef] [Green Version]
  55. Morris, T.J.; Butcher, L.M.; Feber, A.; Teschendorff, A.E.; Chakravarthy, A.R.; Wojdacz, T.K.; Beck, S. ChAMP: 450k Chip Analysis Methylation Pipeline. Bioinforma. Oxf. Engl. 2014, 30, 428–430. [Google Scholar] [CrossRef] [PubMed]
  56. R Core Team. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria, 2020. [Google Scholar]
  57. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. Limma Powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef]
  58. Du, P.; Zhang, X.; Huang, C.-C.; Jafari, N.; Kibbe, W.A.; Hou, L.; Lin, S.M. Comparison of Beta-Value and M-Value Methods for Quantifying Methylation Levels by Microarray Analysis. BMC Bioinform. 2010, 11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Mansell, G.; Gorrie-Stone, T.J.; Bao, Y.; Kumari, M.; Schalkwyk, L.S.; Mill, J.; Hannon, E. Guidance for DNA Methylation Studies: Statistical Insights from the Illumina EPIC Array. BMC Genom. 2019, 20, 366. [Google Scholar] [CrossRef] [Green Version]
  60. Liao, Y.; Wang, J.; Jaehnig, E.J.; Shi, Z.; Zhang, B. WebGestalt 2019: Gene Set Analysis Toolkit with Revamped UIs and APIs. Nucleic Acids Res. 2019, 47, W199–W205. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Kanehisa, M. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  62. McGuinness, D.; McGlynn, L.M.; Johnson, P.C.; MacIntyre, A.; Batty, G.D.; Burns, H.; Cavanagh, J.; Deans, K.A.; Ford, I.; McConnachie, A.; et al. Socio-Economic Status Is Associated with Epigenetic Differences in the PSoBid Cohort. Int. J. Epidemiol. 2012, 41, 151–160. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Perng, W.; Rozek, L.S.; Mora-Plazas, M.; Duchin, O.; Marin, C.; Forero, Y.; Baylin, A.; Villamor, E. Micronutrient Status and Global DNA Methylation in School-Age Children. Epigenetics 2012, 7, 1133–1141. [Google Scholar] [CrossRef]
  64. Chen, Y.; Li, Y.; Hsieh, T.; Wang, C.; Cheng, K.; Wang, L.; Lin, T.; Cheung, C.H.A.; Wu, C.; Chiang, H. Aging-induced Akt Activation Involves in Aging-related Pathologies and Aβ-induced Toxicity. Aging Cell 2019, 18. [Google Scholar] [CrossRef]
  65. Griffin, R.J.; Moloney, A.; Kelliher, M.; Johnston, J.A.; Ravid, R.; Dockery, P.; O’Connor, R.; O’Neill, C. Activation of Akt/PKB, Increased Phosphorylation of Akt Substrates and Loss and Altered Distribution of Akt and PTEN Are Features of Alzheimer’s Disease Pathology. J. Neurochem. 2005, 93, 105–117. [Google Scholar] [CrossRef] [PubMed]
  66. Nyeko, R.; Kalyesubula, I.; Mworozi, E.; Bachou, H. Lactose Intolerance among Severely Malnourished Children with Diarrhoea Admitted to the Nutrition Unit, Mulago Hospital, Uganda. BMC Pediatr. 2010, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Nichols, B.; Dudley, M.; Nichols, V.; Putman, M.; Avery, S.; Fraley, J.; Quaroni, A.; Shiner, M.; Carrazza, F. Effects of Malnutrition on Expression and Activity of Lactase in Children. Gastroenterology 1997, 112, 742–751. [Google Scholar] [CrossRef]
  68. Horvath, S.; Gurven, M.; Levine, M.E.; Trumble, B.C.; Kaplan, H.; Allayee, H.; Ritz, B.R.; Chen, B.; Lu, A.T.; Rickabaugh, T.M.; et al. An Epigenetic Clock Analysis of Race/Ethnicity, Sex, and Coronary Heart Disease. Genome Biol. 2016, 17, 171. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Genes 12 01408 g001 550
Figure 1. Correlation matrix of the epigenetic clocks and the chronological age. Note: The correlation was performed in both cohorts, Tlaltizapan and urban-raised, and collapsed. Blue colors represent a positive correlation, and red colors represent a negative correlation. The square size signifies the p-value of the correlation (the significance is directly proportional to the square size; all the correlations were statistically significant (p-value < 0.05).
Figure 1. Correlation matrix of the epigenetic clocks and the chronological age. Note: The correlation was performed in both cohorts, Tlaltizapan and urban-raised, and collapsed. Blue colors represent a positive correlation, and red colors represent a negative correlation. The square size signifies the p-value of the correlation (the significance is directly proportional to the square size; all the correlations were statistically significant (p-value < 0.05).
Genes 12 01408 g001
Genes 12 01408 g002 550
Figure 2. Differences between individuals with epigenetic age acceleration and years of schooling. Note: Blue colors represent individuals with no acceleration on epigenetic clocks, and red represents individuals with accelerated epigenetic aging.
Figure 2. Differences between individuals with epigenetic age acceleration and years of schooling. Note: Blue colors represent individuals with no acceleration on epigenetic clocks, and red represents individuals with accelerated epigenetic aging.
Genes 12 01408 g002
Table
Table 1. Clinical, biochemical, and epigenetic variables.
Table 1. Clinical, biochemical, and epigenetic variables.
Tlaltizapan Cohort (n = 32)Urban Raised (n = 7)Stat (p-Value)
Age52.26 (0.42)62.41 (0.35)3.98 (0.0105)
Gender
Female19 (59.37)5 (71.43)0.03 (0.8690)
Male13 (40.63)2 (28.57)0.03 (0.8690)
Years of schooling10.31 (3.55)13.86 (8.99)1.03 (0.3421)
BMI29.35 (4.61)29.12 (4.61)−0.12 (0.9063)
Visceral fat3.09 (1.01)3.17 (0.86)0.21 (0.8379)
Biochemical variables
Glucose108.34 (49.87)125.43 (69.01)0.62 (0.5535)
Triglycerides212.59 (210.05)230.71 (190.92)0.22 (0.8281)
Total Cholesterol196.75 (29.94)236.29 (33.14)2.91 (0.0189)
HDL44.07 (12.84)45.81 (12.09)0.34 (0.7400)
LDL111.98 (38.60)144.30 (44.98)1.76 (0.1155)
Creatinine0.64 (0.18)0.71 (0.14)1.15 (0.2763)
Uric Acid5.62 (0.65)5.21 (1.42)−1.17 (0.2535)
Prealbumine24.92 (9.53)28.82 (5.61)1.37 (0.1970)
Reactive Protein C0.40 (0.17)0.38 (0.08)−0.55 (0.5935)
Transferrine247.76 (87.39)276.40 (42.01)1.24 (0.2341)
Glycosylated Haemoglobin6.45 (1.77)5.73 (1.22)−1.22 (0.2525)
Fibrinogen330.75 (50.13)353.83 (24.79)1.71 (0.1078)
Iron Fixing357.06 (60.39)388.67 (78.15)0.94 (0.3828)
Mental Health
Mini-mental25.56 (2.65)26.42 (1.81)1.04 (0.3167)
Depression6 (18.75)2 (28.57)0.00 (0.9472)
Drugs9 (28.13)5 (71.43)5.18 (0.1591)
Tobacco Index1.50 (5.45)0.86 (1.18)−0.61 (0.5477)
Alcohol28 (87.50)4 (57.14)1.83 (0.1763)
Epigenetic Clocks
Horvath58.37 (3.61)63.84 (4.56)2.96 (0.0138)
Horvath acceleration−0.14 (3.56)0.71 (4.24)0.47 (0.6579)
Hannum56.77 (4.02)60.16 (5.40)1.75 (0.1096)
Hannum acceleration0.05 (3.60)−0.25 (4.45)−0.15 (0.8814)
DNAm PhenoAge49.95 (3.62)54.74 (4.87)2.06 (0.0591)
PhenoAge acceleration0.11 (6.64)−0.58 (4.02)−0.34 (0.7383)
Note. Bold represents statistically significant differences (not adjusted for the multiple tests).
Table
Table 2. Epigenome-wide nominal associations between high and low duration of schooling.
Table 2. Epigenome-wide nominal associations between high and low duration of schooling.
1 PositionCpG Site2 LogFCp-ValueHigh AvgLow AvgGene3 Gene LocCGI 4
2:38496264cg19269093−0.04532.1563 × 10−5 0.87470.8294 IGROpenSea
2:136595281cg04750100−0.08922.5258 × 10−50.41710.3279LCTTSS1500OpenSea
2:237163447cg25305153−0.07793.3160 × 10−50.81260.7347ASB18BodyOpenSea
3:65561644cg05244979−0.05103.4089 × 10−50.68690.6359MAGI1BodyOpenSea
4:113970506cg028151710.11654.5334 × 10−50.42180.5383ANK2TSS1500OpenSea
7:156716133cg03184819−0.09074.9140 × 10−50.25290.1622 IGROpenSea
9:4435234cg08538646−0.07174.2326 × 10−50.85850.7868 IGROpenSea
9:117818174cg07712264−0.05554.0462 × 10−50.70590.6504TNCBodyOpenSea
10:25460855cg15018193−0.04391.1161 × 10−50.75830.7144LOC100128811BodySelf
11:102124935cg06458665−0.03972.2954 × 10−50.86350.8237 IGROpenSea
12:26451968cg13537590−0.05001.0817 × 10−50.71140.6614 IGROpenSea
22:50710746cg22416596−0.04044.9783 × 10−50.68840.6481 IGRShore
Notes. 1 Human genome position of CpG site (GRCh37/hg19). 2 LogFC = logarithm of fold change. 3 Gene Loc = location of CpG site relative to coding gene. 4 CGI = CpG island. 5 IGR = Intergenic region.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gomez-Verjan, J.C.; Esparza-Aguilar, M.; Martín-Martín, V.; Salazar-Perez, C.; Cadena-Trejo, C.; Gutierrez-Robledo, L.M.; Martínez-Magaña, J.J.; Nicolini, H.; Arroyo, P. Years of Schooling Could Reduce Epigenetic Aging: A Study of a Mexican Cohort. Genes 2021, 12, 1408. https://doi.org/10.3390/genes12091408

AMA Style

Gomez-Verjan JC, Esparza-Aguilar M, Martín-Martín V, Salazar-Perez C, Cadena-Trejo C, Gutierrez-Robledo LM, Martínez-Magaña JJ, Nicolini H, Arroyo P. Years of Schooling Could Reduce Epigenetic Aging: A Study of a Mexican Cohort. Genes. 2021; 12(9):1408. https://doi.org/10.3390/genes12091408

Chicago/Turabian Style

Gomez-Verjan, Juan Carlos, Marcelino Esparza-Aguilar, Verónica Martín-Martín, Cecilia Salazar-Perez, Cinthya Cadena-Trejo, Luis Miguel Gutierrez-Robledo, José Jaime Martínez-Magaña, Humberto Nicolini, and Pedro Arroyo. 2021. "Years of Schooling Could Reduce Epigenetic Aging: A Study of a Mexican Cohort" Genes 12, no. 9: 1408. https://doi.org/10.3390/genes12091408

APA Style

Gomez-Verjan, J. C., Esparza-Aguilar, M., Martín-Martín, V., Salazar-Perez, C., Cadena-Trejo, C., Gutierrez-Robledo, L. M., Martínez-Magaña, J. J., Nicolini, H., & Arroyo, P. (2021). Years of Schooling Could Reduce Epigenetic Aging: A Study of a Mexican Cohort. Genes, 12(9), 1408. https://doi.org/10.3390/genes12091408

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