A de novo frameshift variant in MED13 gene in a patient with autism spectrum disorder and magnetic resonance imaging abnormalities mimicking tuberous sclerosis.

The mediator complex subunit 13 (MED13) gene is implicated in neurodevelopmental disorders including autism spectrum disorder (ASD), intellectual disability, and speech delay with varying severity and course. Additional, extra central nervous system, features include eye or vision problems, hypotonia, congenital heart abnormalities, and dysmorphisms. We describe a 7-year- and 4-month-old girl evaluated for ASD whose brain magnetic resonance imaging was suggestive of multiple cortical tubers. The exome sequencing (ES - trio analysis) uncovered a unique, de novo, frameshift variant in the MED13 gene (c.4880del, D1627Vfs*17), with a truncating effect on the protein. This case report thus expands the phenotypic spectrum of MED13-related disorders to include brain abnormalities.

A Novel Heterozygous Mutation c.1627G>T (p.Gly543Cys) in the SLC34A1 Gene in a Male Patient with Recurrent Nephrolithiasis and Early Onset Osteopenia: A Case Report.

Serum phosphate concentration is regulated by renal phosphate reabsorption and mediated by sodium-phosphate cotransporters. Germline mutations in genes encoding these cotransporters have been associated with clinical phenotypes, variably characterized by hyperphosphaturia, hypophosphatemia, recurrent kidney stones, skeletal demineralization, and early onset osteoporosis. We reported a 33-year-old male patient presenting a history of recurrent nephrolithiasis and early onset osteopenia in the lumbar spine and femur. He was tested, through next generation sequencing (NGS), by using a customized multigenic panel containing 33 genes, whose mutations are known to be responsible for the development of congenital parathyroid diseases. Two further genes, SLC34A1 and SLC34A3, encoding two sodium-phosphate cotransporters, were additionally tested. A novel germline heterozygous mutation was identified in the SLC34A1 gene, c.1627G>T (p.Gly543Cys), currently not reported in databases of human gene mutations and scientific literature. SLC34A1 germline heterozygous mutations have been associated with the autosomal dominant hypophosphatemic nephrolithiasis/osteoporosis type 1 (NPHLOP1). Consistently, alongside the clinical features of NPHLOP1, our patient experienced recurrent nephrolithiasis and lumbar and femoral osteopenia at a young age. Genetic screening for the p.Gly453Cys variant and the clinical characterization of his first-degree relatives associated the presence of the variant in one younger brother, presenting renal colic and microlithiasis, suggesting p.Gly453Cys is possibly associated with renal altered function in the NPHLOP1 phenotype.

High frequency of heterozygous rare variants of the SLC34A1 and SLC9A3R1 genes in patients with atypical femur fracture.

OBJECTIVE: Atypical femur fractures (AFFs) are rare fragility fractures originating at the lateral cortex of the femur, affecting the subtrochanteric or diaphyseal area of thebone with a transverse morphology. Occurrence of AFF is specifically associated with a small number of rare monogenic congenital metabolic bone disorders, such as hypophosphatasia, and with long-term treatment with antiresorptiondrugs. The exact pathogenesis of these fractures remains poorly understood and, except for cases of diagnosed HPP or other AFF-causing bone diseases, it is not possible to assess which patients are at higher riskof developing AFFs as a consequence of anti-resorption therapy. DESIGN: We genetically screened 25 unrelated patients who had developed at least one AFF. INTERVENTION: Genetic screening was performed through a nextgeneration sequencing analysis with a customized panel containing 76 human genes involved in the regulation of the mineralization processWe genetically screened 25 unrelated patients who had developed at least one AFF. RESULTS: We found a relatively high frequency (32.0%) of heterozygous rare variants inthe SLC34A1 and SLC9A3R1 genes, two genes whose heterozygous inactivating mutations have been respectively associated with autosomal dominant hypophosphatemic nephrolithiasis/osteoporosis types 1 and 2 (NPHLOP1and NPHLOP2). Other heterozygous rare variants were found in the BMPR1B, CYP27B1, FBN1, MEPE, PIGO, and PHOSPHO1 genes, each in a single AFF case (4.0%). CONCLUSIONS AND RELEVANCE: Our findings suggest that rarevariants of SLC34A1 and SLC9A3R1 could represent a possible genetic risk factor for the occurrence of AFFs. On the other hand, AFFs could represent an unsuspected clinical manifestation and/or an anti-resorption therapycorrelatedadverse event in patients with NPHLOP disorders.

A Targeted Epigenetic Clock for the Prediction of Biological Age.

Epigenetic clocks were initially developed to track chronological age, but accumulating evidence indicates that they can also predict biological age. They are usually based on the analysis of DNA methylation by genome-wide methods, but targeted approaches, based on the assessment of a small number of CpG sites, are advisable in several settings. In this study, we developed a targeted epigenetic clock purposely optimized for the measurement of biological age. The clock includes six genomic regions mapping in ELOVL2, NHLRC1, AIM2, EDARADD, SIRT7 and TFAP2E genes, selected from a re-analysis of existing microarray data, whose DNA methylation is measured by EpiTYPER assay. In healthy subjects (n = 278), epigenetic age calculated using the targeted clock was highly correlated with chronological age (Spearman correlation = 0.89). Most importantly, and in agreement with previous results from genome-wide clocks, epigenetic age was significantly higher and lower than expected in models of increased (persons with Down syndrome, n = 62) and decreased (centenarians, n = 106; centenarians' offspring, n = 143; nutritional intervention in elderly, n = 233) biological age, respectively. These results support the potential of our targeted epigenetic clock as a new marker of biological age and open its evaluation in large cohorts to further promote the assessment of biological age in healthcare practice.

Clinical, Radiometabolic and Immunologic Effects of Olaparib in Locally Advanced Triple Negative Breast Cancer: The OLTRE Window of Opportunity Trial.

INTRODUCTION: Olaparib is effective in metastatic triple negative breast cancer (TNBC) carrying germline mutations in DNA damage repair (DDR) genes BRCA1/2 (gBRCA-mut). The OLTRE window-of-opportunity trial preliminarily investigated potential pathologic, radiometabolic and immune biomarkers of early-response to olaparib in gBRCA-wild-type (wt) TNBC and, as proof-of-concept in gBRCA-mut HER2-negative BC. METHODS: Patients received olaparib for 3 weeks (3w) before standard neoadjuvant chemotherapy and underwent multiple FDG18-PET/CT scan (basal, after olaparib), clinical assessments (basal, every 3w), tumor biopsies and blood samplings (baseline, after olaparib). Clinical and radiometabolic responses were evaluated according to RECIST1.1 and PERCIST criteria. RESULTS: 27 patients with gBRCA-wt TNBC and 8 with gBRCA-mut BC (6 TNBC, 2 HR+/HER2-negative) were enrolled. Three (11.1%) patients showed mutations in non-BRCA1/2 DDR genes and 4 (14.8%) in other genes. 3w olaparib induced 16/35 and 15/27 partial clinical and radiometabolic responses, including in 40.7% and 50.0% gBRCA-wt patients. gBRCA-mut tumors presented numerically higher tumor-infiltrating lymphocytes (TILs) levels and PD-L1 positive tumors. Clinical responders experienced a reduction in T-regs/T-eff ratio (p=0.05), B and NK lymphocytes (p=0.003 both), with an average increase in T-helpers rate (p<0.001) and CD4/CD8 ratio (p=0.02). Ki67% and TILs did not vary significantly (p=0.67 and p=0.77). A numerical increase in PD-L1 positive cases after olaparib was observed, though non-significant (p=0.134). No differences were observed according to gBRCA status and type of response. CONCLUSIONS: Early-stage TNBC might be a target population for olaparib, irrespective of gBRCA mutations. Future trials should combine TILs, PD-L1 and gBRCA status to better identify candidates for escalated/de-escalated treatment strategies including olaparib.

Whole-genome sequencing analysis of semi-supercentenarians.

Extreme longevity is the paradigm of healthy aging as individuals who reached the extreme decades of human life avoided or largely postponed all major age-related diseases. In this study, we sequenced at high coverage (90X) the whole genome of 81 semi-supercentenarians and supercentenarians [105+/110+] (mean age: 106.6 ± 1.6) and of 36 healthy unrelated geographically matched controls (mean age 68.0 ± 5.9) recruited in Italy. The results showed that 105+/110+ are characterized by a peculiar genetic background associated with efficient DNA repair mechanisms, as evidenced by both germline data (common and rare variants) and somatic mutations patterns (lower mutation load if compared to younger healthy controls). Results were replicated in a second independent cohort of 333 Italian centenarians and 358 geographically matched controls. The genetics of 105+/110+ identified DNA repair and clonal haematopoiesis as crucial players for healthy aging and for the protection from cardiovascular events.

Genomic history of the Italian population recapitulates key evolutionary dynamics of both Continental and Southern Europeans.

BACKGROUND: The cline of human genetic diversity observable across Europe is recapitulated at a micro-geographic scale by variation within the Italian population. Besides resulting from extensive gene flow, this might be ascribable also to local adaptations to diverse ecological contexts evolved by people who anciently spread along the Italian Peninsula. Dissecting the evolutionary history of the ancestors of present-day Italians may thus improve the understanding of demographic and biological processes that contributed to shape the gene pool of European populations. However, previous SNP array-based studies failed to investigate the full spectrum of Italian variation, generally neglecting low-frequency genetic variants and examining a limited set of small effect size alleles, which may represent important determinants of population structure and complex adaptive traits. To overcome these issues, we analyzed 38 high-coverage whole-genome sequences representative of population clusters at the opposite ends of the cline of Italian variation, along with a large panel of modern and ancient Euro-Mediterranean genomes. RESULTS: We provided evidence for the early divergence of Italian groups dating back to the Late Glacial and for Neolithic and distinct Bronze Age migrations having further differentiated their gene pools. We inferred adaptive evolution at insulin-related loci in people from Italian regions with a temperate climate, while possible adaptations to pathogens and ultraviolet radiation were observed in Mediterranean Italians. Some of these adaptive events may also have secondarily modulated population disease or longevity predisposition. CONCLUSIONS: We disentangled the contribution of multiple migratory and adaptive events in shaping the heterogeneous Italian genomic background, which exemplify population dynamics and gene-environment interactions that played significant roles also in the formation of the Continental and Southern European genomic landscapes.

Author Correction: Erythropoietin (EPO) haplotype associated with all-cause mortality in a cohort of Italian patients with Type-2 Diabetes.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Erythropoietin (EPO) haplotype associated with all-cause mortality in a cohort of Italian patients with Type-2 Diabetes.

Type-2 Diabetes (T2D), diabetic complications, and their clinical risk factors harbor a substantial genetic component but the genetic factors contributing to overall diabetes mortality remain unknown. Here, we examined the association between genetic variants at 21 T2D-susceptibility loci and all-cause mortality in an elderly cohort of 542 Italian diabetic patients who were followed for an average of 12.08 years. Univariate Cox regression analyses detected age, waist-to-hip ratio (WHR), glycosylated haemoglobin (HbA1c), diabetes duration, retinopathy, nephropathy, chronic kidney disease (CKD), and anaemia as predictors of all-cause mortality. When Cox proportional hazards multivariate models adjusted for these factors were run, three erythropoietin (EPO) genetic variants in linkage disequilibrium (LD) with each other (rs1617640-T/G, rs507392-T/C and rs551238-A/C) were significantly (False Discovery Rate < 0.1) associated with mortality. Haplotype multivariate analysis revealed that patients carrying the G-C-C haplotype have an increased probability of survival, while an opposite effect was observed among subjects carrying the T-T-A haplotype. Our findings provide evidence that the EPO gene is an independent predictor of mortality in patients with T2D. Thus, understanding the mechanisms by which the genetic variability of EPO affects the mortality of T2D patients may provide potential targets for therapeutic interventions to improve the survival of these patients.

Genomic stability, anti-inflammatory phenotype, and up-regulation of the RNAseH2 in cells from centenarians.

Current literature agrees on the notion that efficient DNA repair favors longevity across evolution. The DNA damage response machinery activates inflammation and type I interferon signaling. Both pathways play an acknowledged role in the pathogenesis of a variety of age-related diseases and are expected to be detrimental for human longevity. Here, we report on the anti-inflammatory molecular make-up of centenarian's fibroblasts (low levels of IL-6, type 1 interferon beta, and pro-inflammatory microRNAs), which is coupled with low level of DNA damage (measured by comet assay and histone-2AX activation) and preserved telomere length. In the same cells, high levels of the RNAseH2C enzyme subunit and low amounts of RNAseH2 substrates, i.e. cytoplasmic RNA:DNA hybrids are present. Moreover, RNAseH2C locus is hypo-methylated and RNAseH2C knock-down up-regulates IL-6 and type 1 interferon beta in centenarian's fibroblasts. Interestingly, RNAseH2C locus is hyper-methylated in vitro senescent cells and in tissues from atherosclerotic plaques and breast tumors. Finally, extracellular vesicles from centenarian's cells up-regulate RNAseH2C expression and dampen the pro-inflammatory phenotype of fibroblasts, myeloid, and cancer cells. These data suggest that centenarians are endowed with restrained DNA damage-induced inflammatory response, that may facilitate their escape from the deleterious effects of age-related chronic inflammation.

Molecular Aging of Human Liver: An Epigenetic/Transcriptomic Signature.

The feasibility of liver transplantation from old healthy donors suggests that this organ is able to preserve its functionality during aging. To explore the biological basis of this phenomenon, we characterized the epigenetic profile of liver biopsies collected from 45 healthy liver donors ranging from 13 to 90 years old using the Infinium HumanMethylation450 BeadChip. The analysis indicates that a large remodeling in DNA methylation patterns occurs, with 8,823 age-associated differentially methylated CpG probes. Notably, these age-associated changes tended to level off after the age of 60, as confirmed by Horvath's clock. Using stringent selection criteria, we further identified a DNA methylation signature of aging liver including 75 genomic regions. We demonstrated that this signature is specific for liver compared to other tissues and that it is able to detect biological age-acceleration effects associated with obesity. Finally, we combined DNA methylation measurements with available expression data. Although the intersection between the two omic characterizations was low, both approaches suggested a previously unappreciated role of epithelial-mesenchymal transition and Wnt-signaling pathways in the aging of human liver.

First evidence of association between past environmental exposure to dioxin and DNA methylation of CYP1A1 and IGF2 genes in present day Vietnamese population.

During the Vietnam War, the United States military sprayed over 74 million litres of Agent Orange (AO) to destroy forest cover as a counterinsurgency tactic in Vietnam, Laos and Cambodia. The main ingredient was contaminated by 2,3,7,8-tetrachlorodibenzo-paradioxin (TCDD). DNA methylation (DNAm) differences are potential biomarker of environmental toxicants exposure. The aim of this study was to perform a preliminary investigation of the DNAm levels from peripheral blood of the present-day Vietnamese population, including individuals whose parents, according to historical data, were exposed to AO/TCDD during the war. 94 individuals from heavily sprayed areas (cases) and 94 individuals from non-sprayed areas (controls) were studied, and historical data on alleged exposure of parents collected. 94 cases were analysed considering those whose father/parents participated in the war (N = 29) and considering the place of residence of both parents (64 living in sprayed areas versus 30 in non-contaminated areas). DNAm levels in CYP1A1 and IGF2 genes were measured (MALDI-TOF technology). The analyses showed that: 1) one CpG site in the CYP1A1 and one in the IGF2 gene showed significant differences in DNAm levels between cases and controls; 2) the CYP1A1 region resulted to be hypomethylated (in 9 out of 16 sites/units; p-val<0.01) in 29 individuals whose father/parents participated in the war in the spray zones; 3) we showed that the place of residence of both parents influenced methylation levels of the CYP1A1 and IGF2 genes (p-val<0.05). In conclusion this study indicates that past environmental exposure to dioxin (AO/TCDD) shapes the DNAm profile of CYP1A1 and that the place of living for parents in former spray zones influences DNAm of CYP1A1 and IGF2 genes. These results open the way to new applications of DNAm as potential biomarker(s) of past human exposure to dioxin.

Impact of demography and population dynamics on the genetic architecture of human longevity.

The study of the genetics of longevity has been mainly addressed by GWASs that considered subjects from different populations to reach higher statistical power. The "price to pay" is that population-specific evolutionary histories and trade-offs were neglected in the investigation of gene-environment interactions. We propose a new "diachronic" approach that considers processes occurred at both evolutionary and lifespan timescales. We focused on a well-characterized population in terms of evolutionary history (i.e. Italians) and we generated genome-wide data for 333 centenarians from the peninsula and 773 geographically-matched healthy individuals. Obtained results showed that: (i) centenarian genomes are enriched for an ancestral component likely shaped by pre-Neolithic migrations; (ii) centenarians born in Northern Italy unexpectedly clustered with controls from Central/Southern Italy suggesting that Neolithic and Bronze Age gene flow did not favor longevity in this population; (iii) local past adaptive events in response to pathogens and targeting arachidonic acid metabolism became favorable for longevity; (iv) lifelong changes in the frequency of several alleles revealed pleiotropy and trade-off mechanisms crucial for longevity. Therefore, we propose that demographic history and ancient/recent population dynamics need to be properly considered to identify genes involved in longevity, which can differ in different temporal/spatial settings.

Genes associated with Type 2 Diabetes and vascular complications.

Type 2 Diabetes (T2D) is a chronic disease associated with a number of micro- and macrovascular complications that increase the morbidity and mortality of patients. The risk of diabetic complications has a strong genetic component. To this end, we sought to evaluate the association of 40 single nucleotide polymorphisms (SNPs) in 21 candidate genes with T2D and its vascular complications in 503 T2D patients and 580 healthy controls. The genes were chosen because previously reported to be associated with T2D complications and/or with the aging process. We replicated the association of T2D risk with IGF2BP rs4402960 and detected novel associations with TERT rs2735940 and rs2736098. The addition of these SNPs to a model including traditional risk factors slightly improved risk prediction. After stratification of patients according to the presence/absence of vascular complications, we found significant associations of variants in the CAT, FTO, and UCP1 genes with diabetic retinopathy and nephropathy. Additionally, a variant in the ADIPOQ gene was found associated with macrovascular complications. Notably, these genes are involved in some way in mitochondrial biology and reactive oxygen species regulation. Hence, our findings strongly suggest a potential link between mitochondrial oxidative homeostasis and individual predisposition to diabetic vascular complications.

Epigenetically induced ectopic expression of UNCX impairs the proliferation and differentiation of myeloid cells.

We here describe a leukemogenic role of the homeobox gene UNCX, activated by epigenetic modifications in acute myeloid leukemia (AML). We found the ectopic activation of UNCX in a leukemia patient harboring a t(7;10)(p22;p14) translocation, in 22 of 61 of additional cases [a total of 23 positive patients out of 62 (37.1%)], and in 6 of 75 (8%) of AML cell lines. UNCX is embedded within a low-methylation region (canyon) and encodes for a transcription factor involved in somitogenesis and neurogenesis, with specific expression in the eye, brain, and kidney. UNCX expression turned out to be associated, and significantly correlated, with DNA methylation increase at its canyon borders based on data in our patients and in archived data of patients from The Cancer Genome Atlas. UNCX-positive and -negative patients displayed significant differences in their gene expression profiles. An enrichment of genes involved in cell proliferation and differentiation, such as MAP2K1 and CCNA1, was revealed. Similar results were obtained in UNCX-transduced CD34+ cells, associated with low proliferation and differentiation arrest. Accordingly, we showed that UNCX expression characterizes leukemia cells at their early stage of differentiation, mainly M2 and M3 subtypes carrying wild-type NPM1 We also observed that UNCX expression significantly associates with an increased frequency of acute promyelocytic leukemia with PML-RARA and AML with t(8;21)(q22;q22.1); RUNX1-RUNX1T1 classes, according to the World Health Organization disease classification. In summary, our findings suggest a novel leukemogenic role of UNCX, associated with epigenetic modifications and with impaired cell proliferation and differentiation in AML.

Acceleration of leukocytes' epigenetic age as an early tumor and sex-specific marker of breast and colorectal cancer.

Changes in blood epigenetic age have been associated with several pathological conditions and have recently been described to anticipate cancer development. In this work, we analyze a publicly available leukocytes methylation dataset to evaluate the relation between DNA methylation age and the prospective development of specific types of cancer. We calculated DNA methylation age acceleration using five state-of-the-art estimators (three multi-site: Horvath, Hannum, Weidner; and two CpG specific: ELOV2 and FHL2) in a cohort including 845 subjects from the EPIC-Italy project and we compared 424 samples that remained cancer-free over the approximately ten years of follow-up with 235 and 166 subjects who developed breast and colorectal cancer, respectively. We show that the epigenetic age estimated from blood DNA methylation data is statistically significantly associated to future breast and male colorectal cancer development. These results are corroborated by survival analysis that shows significant association between age acceleration and cancer incidence suggesting that the chance of developing age-related diseases may be predicted by circulating epigenetic markers, with a dependence upon tumor type, sex and age estimator. These are encouraging results towards the non-invasive and perspective usage of epigenetic biomarkers.

The epigenetic landscape of age-related diseases: the geroscience perspective.

In this review, we summarize current knowledge regarding the epigenetics of age-related diseases, focusing on those studies that have described DNA methylation landscape in cardio-vascular diseases, musculoskeletal function and frailty. We stress the importance of adopting the conceptual framework of "geroscience", which starts from the observation that advanced age is the major risk factor for several of these pathologies and aims at identifying the mechanistic links between aging and age-related diseases. DNA methylation undergoes a profound remodeling during aging, which includes global hypomethylation of the genome, hypermethylation at specific loci and an increase in inter-individual variation and in stochastic changes of DNA methylation values. These epigenetic modifications can be an important contributor to the development of age-related diseases, but our understanding on the complex relationship between the epigenetic signatures of aging and age-related disease is still poor. The most relevant results in this field come from the use of the so called "epigenetics clocks" in cohorts of subjects affected by age-related diseases. We report these studies in final section of this review.

The methylation of nuclear and mitochondrial DNA in ageing phenotypes and longevity.

An increasing body of data is progressively indicating that the comprehension of the epigenetic landscape, actively integrated with the genetic elements, is crucial to delineate the molecular basis of the inter-individual complexity of ageing process. Indeed, it has emerged that DNA methylation changes occur during ageing, consisting mainly in a progressive process of genome demethylation, in a hypermethylation of gene-specific CpG dinucleotides, as well as in an inter-individual divergence of the epigenome due to stochastic events and environmental exposures throughout life, namely as epigenetic drift. Additionally, it has also come to light an implication of the mitochondrial genome in the regulation of the intracellular epigenetic landscape, as demonstrated by the being itself object of epigenetic modifications. An overview of DNA methylation changes occurring during ageing process at both nuclear and mitochondrial level will be described in this review, also taking into account the recent and promising data available on the 5-hydroxymethylcytosine.

Systemic Age-Associated DNA Hypermethylation of ELOVL2 Gene: In Vivo and In Vitro Evidences of a Cell Replication Process.

Epigenetic remodeling is one of the major features of the aging process. We recently demonstrated that DNA methylation of ELOVL2 and FHL2 CpG islands is highly correlated with age in whole blood. Here we investigated several aspects of age-associated hypermethylation of ELOVL2 and FHL2. We showed that ELOVL2 methylation is significantly different in primary dermal fibroblast cultures from donors of different ages. Using epigenomic data from public resources, we demonstrated that most of the tissues show ELOVL2 and FHL2 hypermethylation with age. Interestingly, ELOVL2 hypermethylation was not found in tissues with very low replication rate. We demonstrated that ELOVL2 hypermethylation is associated with in vitro cell replication rather than with senescence. We confirmed intra-individual hypermethylation of ELOVL2 and FHL2 in longitudinally assessed participants from the Doetinchem Cohort Study. Finally we showed that, although the methylation of the two loci is not associated with longevity/mortality in the Leiden Longevity Study, ELOVL2 methylation is associated with cytomegalovirus status in nonagenarians, which could be informative of a higher number of replication events in a fraction of whole-blood cells. Collectively, these results indicate that ELOVL2 methylation is a marker of cell divisions occurring during human aging.

Epigenetic Variability across Human Populations: A Focus on DNA Methylation Profiles of the KRTCAP3, MAD1L1 and BRSK2 Genes.

Natural epigenetic diversity has been suggested as a key mechanism in microevolutionary processes due to its capability to create phenotypic variability within individuals and populations. It constitutes an important reservoir of variation potentially useful for rapid adaptation in response to environmental stimuli. The analysis of population epigenetic structure represents a possible tool to study human adaptation and to identify external factors that are able to naturally shape human DNA methylation variability. The aim of this study is to investigate the dynamics that create epigenetic diversity between and within different human groups. To this end, we first used publicly available epigenome-wide data to explore population-specific DNA methylation changes that occur at macro-geographic scales. Results from this analysis suggest that nutrients, UVA exposure and pathogens load might represent the main environmental factors able to shape DNA methylation profiles. Then, we evaluated DNA methylation of candidate genes (KRTCAP3, MAD1L1, and BRSK2), emerged from the previous analysis, in individuals belonging to different populations from Morocco, Nigeria, Philippines, China, and Italy, but living in the same Italian city. DNA methylation of the BRSK2 gene is significantly different between Moroccans and Nigerians (pairwise t-test: CpG 6 P-value = 5.2*10 (-) (3); CpG 9 P-value = 2.6*10 (-) (3); CpG 10 P-value = 3.1*10 (-) (3); CpG 11 P-value = 2.8*10 (-) (3)). Comprehensively, these results suggest that DNA methylation diversity is a source of variability in human groups at macro and microgeographical scales and that population demographic and adaptive histories, as well as the individual ancestry, actually influence DNA methylation profiles.