Evaluation of DNA Methylation Profiles of LINE-1, Alu and Ribosomal DNA Repeats in Human Cell Lines Exposed to Radiofrequency Radiation.

A large body of evidence indicates that environmental agents can induce alterations in DNA methylation (DNAm) profiles. Radiofrequency electromagnetic fields (RF-EMFs) are radiations emitted by everyday devices, which have been classified as "possibly carcinogenic"; however, their biological effects are unclear. As aberrant DNAm of genomic repetitive elements (REs) may promote genomic instability, here, we sought to determine whether exposure to RF-EMFs could affect DNAm of different classes of REs, such as long interspersed nuclear elements-1 (LINE-1), Alu short interspersed nuclear elements and ribosomal repeats. To this purpose, we analysed DNAm profiles of cervical cancer and neuroblastoma cell lines (HeLa, BE(2)C and SH-SY5Y) exposed to 900 MHz GSM-modulated RF-EMF through an Illumina-based targeted deep bisulfite sequencing approach. Our findings showed that radiofrequency exposure did not affect the DNAm of Alu elements in any of the cell lines analysed. Conversely, it influenced DNAm of LINE-1 and ribosomal repeats in terms of both average profiles and organisation of methylated and unmethylated CpG sites, in different ways in each of the three cell lines studied.

DLX5/6 GABAergic Expression Affects Social Vocalization: Implications for Human Evolution.

DLX5 and DLX6 are two closely related transcription factors involved in brain development and in GABAergic differentiation. The DLX5/6 locus is regulated by FoxP2, a gene involved in language evolution and has been associated with neurodevelopmental disorders and mental retardation. Targeted inactivation of Dlx5/6 in mouse GABAergic neurons (Dlx5/6VgatCre mice) results in behavioral and metabolic phenotypes notably increasing lifespan by 33%. Here, we show that Dlx5/6VgatCre mice present a hyper-vocalization and hyper-socialization phenotype. While only 7% of control mice emitted more than 700 vocalizations/10 min, 30% and 56% of heterozygous or homozygous Dlx5/6VgatCre mice emitted more than 700 and up to 1,400 calls/10 min with a higher proportion of complex and modulated calls. Hyper-vocalizing animals were more sociable: the time spent in dynamic interactions with an unknown visitor was more than doubled compared to low-vocalizing individuals. The characters affected by Dlx5/6 in the mouse (sociability, vocalization, skull, and brain shape…) overlap those affected in the "domestication syndrome". We therefore explored the possibility that DLX5/6 played a role in human evolution and "self-domestication" comparing DLX5/6 genomic regions from Neanderthal and modern humans. We identified an introgressed Neanderthal haplotype (DLX5/6-N-Haplotype) present in 12.6% of European individuals that covers DLX5/6 coding and regulatory sequences. The DLX5/6-N-Haplotype includes the binding site for GTF2I, a gene associated with Williams-Beuren syndrome, a hyper-sociability and hyper-vocalization neurodevelopmental disorder. The DLX5/6-N-Haplotype is significantly underrepresented in semi-supercentenarians (>105 years of age), a well-established human model of healthy aging and longevity, suggesting their involvement in the coevolution of longevity, sociability, and speech.

A Genome-Wide Association Study of 2304 Extreme Longevity Cases Identifies Novel Longevity Variants.

We performed a genome-wide association study (GWAS) of human extreme longevity (EL), defined as surviving past the 99th survival percentile, by aggregating data from four centenarian studies. The combined data included 2304 EL cases and 5879 controls. The analysis identified a locus in CDKN2B-AS1 (rs6475609, p = 7.13 × 10-8) that almost reached genome-wide significance and four additional loci that were suggestively significant. Among these, a novel rare variant (rs145265196) on chromosome 11 had much higher longevity allele frequencies in cases of Ashkenazi Jewish and Southern Italian ancestry compared to cases of other European ancestries. We also correlated EL-associated SNPs with serum proteins to link our findings to potential biological mechanisms that may be related to EL and are under genetic regulation. The findings from the proteomic analyses suggested that longevity-promoting alleles of significant genetic variants either provided EL cases with more youthful molecular profiles compared to controls or provided some form of protection from other illnesses, such as Alzheimer's disease, and disease progressions.

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.

Dissecting the Pre-Columbian Genomic Ancestry of Native Americans along the Andes-Amazonia Divide.

Extensive European and African admixture coupled with loss of Amerindian lineages makes the reconstruction of pre-Columbian history of Native Americans based on present-day genomes extremely challenging. Still open questions remain about the dispersals that occurred throughout the continent after the initial peopling from the Beringia, especially concerning the number and dynamics of diffusions into South America. Indeed, if environmental and historical factors contributed to shape distinct gene pools in the Andes and Amazonia, the origins of this East-West genetic structure and the extension of further interactions between populations residing along this divide are still not well understood. To this end, we generated new high-resolution genome-wide data for 229 individuals representative of one Central and ten South Amerindian ethnic groups from Mexico, Peru, Bolivia, and Argentina. Low levels of European and African admixture in the sampled individuals allowed the application of fine-scale haplotype-based methods and demographic modeling approaches. These analyses revealed highly specific Native American genetic ancestries and great intragroup homogeneity, along with limited traces of gene flow mainly from the Andes into Peruvian Amazonians. Substantial amount of genetic drift differentially experienced by the considered populations underlined distinct patterns of recent inbreeding or prolonged isolation. Overall, our results support the hypothesis that all non-Andean South Americans are compatible with descending from a common lineage, while we found low support for common Mesoamerican ancestors of both Andeans and other South American groups. These findings suggest extensive back-migrations into Central America from non-Andean sources or conceal distinct peopling events into the Southern Continent.

Genetics of Human Longevity Within an Eco-Evolutionary Nature-Nurture Framework.

Human longevity is a complex trait, and to disentangle its basis has a great theoretical and practical consequences for biomedicine. The genetics of human longevity is still poorly understood despite several investigations that used different strategies and protocols. Here, we argue that such rather disappointing harvest is largely because of the extraordinary complexity of the longevity phenotype in humans. The capability to reach the extreme decades of human lifespan seems to be the result of an intriguing mixture of gene-environment interactions. Accordingly, the genetics of human longevity is here described as a highly context-dependent phenomenon, within a new integrated, ecological, and evolutionary perspective, and is presented as a dynamic process, both historically and individually. The available literature has been scrutinized within this perspective, paying particular attention to factors (sex, individual biography, family, population ancestry, social structure, economic status, and education, among others) that have been relatively neglected. The strength and limitations of the most powerful and used tools, such as genome-wide association study and whole-genome sequencing, have been discussed, focusing on prominently emerged genes and regions, such as apolipoprotein E, Forkhead box O3, interleukin 6, insulin-like growth factor-1, chromosome 9p21, 5q33.3, and somatic mutations among others. The major results of this approach suggest that (1) the genetics of longevity is highly population specific; (2) small-effect alleles, pleiotropy, and the complex allele timing likely play a major role; (3) genetic risk factors are age specific and need to be integrated in the light of the geroscience perspective; (4) a close relationship between genetics of longevity and genetics of age-related diseases (especially cardiovascular diseases) do exist. Finally, the urgent need of a global approach to the largely unexplored interactions between the 3 genetics of human body, that is, nuclear, mitochondrial, and microbiomes, is stressed. We surmise that the comprehensive approach here presented will help in increasing the above-mentioned harvest.

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.

A bio-cultural approach to the study of food choice: The contribution of taste genetics, population and culture.

The study of food choice, one of the most complex human traits, requires an integrated approach that takes into account environmental, socio-cultural and biological diversity. We recruited 183 volunteers from four geo-linguistic groups and highly diversified in terms of both genetic background and food habits from whom we collected genotypes and phenotypes tightly linked to taste perception. We confirmed previous genetic associations, in particular with stevioside perception, and noted significant differences in food consumption: in particular, broccoli, mustard and beer consumption scores were significantly higher (Adjusted P = 0.02, Adjusted P < 0.0001 and Adjusted P = 0.01, respectively) in North Europeans, when compared to the other groups. Licorice and Parmesan cheese showed lower consumption and liking scores in the Sri Lankan group (Adjusted P = 0.001 and Adjusted P < 0.001, respectively). We also highlighted how rs860170 (TAS2R16) strongly differentiated populations and was associated to salicin bitterness perception. Identifying genetic variants on chemosensory receptors that vary across populations and show associations with taste perception and food habits represents a step towards a better comprehension of this complex trait, aimed at improving the individual health status. This is the first study that concurrently explores the contribution of genetics, population diversity and cultural aspects in taste perception and food consumption.

Assessing the combined effect of extremely low-frequency magnetic field exposure and oxidative stress on LINE-1 promoter methylation in human neural cells.

Extremely low frequency magnetic fields (ELF-MF) have been classified as "possibly carcinogenic", but their genotoxic effects are still unclear. Recent findings indicate that epigenetic mechanisms contribute to the genome dysfunction and it is well known that they are affected by environmental factors. To our knowledge, to date the question of whether exposure to ELF-MF can influence epigenetic modifications has been poorly addressed. In this paper, we investigated whether exposure to ELF-MF alone and in combination with oxidative stress (OS) can affect DNA methylation, which is one of the most often studied epigenetic modification. To this end, we analyzed the DNA methylation levels of the 5'untranslated region (5'UTR) of long interspersed nuclear element-1s (LINE-1 or L1), which are commonly used to evaluate the global genome methylation level. Human neural cells (BE(2)C) were exposed for 24 and 48 h to extremely low frequency pulsed magnetic field (PMF; 50 Hz, 1 mT) in combination with OS. The methylation levels of CpGs located in L1 5'UTR region were measured by MassARRAY EpiTYPER. The results indicate that exposures to the single agents PMF and OS induced weak decreases and increases of DNA methylation levels at different CpGs. However, the combined exposure to PMF and OS lead to significant decrease of DNA methylation levels at different CpG sites. Most of the changes were transient, suggesting that cells can restore homeostatic DNA methylation patterns. The results are discussed and future research directions outlined.

Inferring chronological age from DNA methylation patterns of human teeth.

OBJECTIVE: Current methods to determine chronological age from modern and ancient remains rely on both morphological and molecular approaches. However, low accuracy and the lack of standardized protocols make the development of alternative methods for the estimation of individual's age even more urgent for several research fields, such as biological anthropology, biodemography, forensics, evolutionary genetics, and ancient DNA studies. Therefore, the aim of this study is to identify genomic regions whose DNA methylation level correlates with age in modern teeth. METHODS: We used MALDI-TOF mass spectrometry to analyze DNA methylation levels of specific CpGs located in the ELOVL2, FHL2, and PENK genes. We considered methylation data from cementum, dentin and pulp of 21 modern teeth (from 17 to 77 years old) to construct a mathematical model able to exploit DNA methylation values to predict age of the individuals. RESULTS: The median difference between the real age and that estimated using DNA methylation values is 1.20 years (SD = 1.9) if DNA is recovered from both cementum and pulp of the same modern teeth, 2.25 years (SD = 2.5) if DNA is recovered from dental pulp, 2.45 years (SD = 3.3) if DNA is extracted from cementum and 7.07 years (SD = 7.0) when DNA is recovered from dentin only. DISCUSSION: We propose for the first time the evaluation of DNA methylation at ELOVL2, FHL2, and PENK genes as a powerful tool to predict age in modern teeth for anthropological applications. Future studies are needed to apply this method also to historical and relatively ancient human teeth.

Positive selection of lactase persistence among people of Southern Arabia.

OBJECTIVE: Frequency patterns of the lactase persistence (LP)-associated -13,915 G allele and archaeological records pointing to substantial role played by southern regions in the peopling and domestication processes that involved the Arabian Peninsula suggest that Southern Arabia plausibly represented the center of diffusion of such adaptive variant. Nevertheless, a well-defined scenario for evolution of Arabian LP is still to be elucidated and the burgeoning archaeological picture of complex human migrations occurred through the peninsula is not matched by an equivalent high-resolution description of genetic variation underlying this adaptive trait. To fill this gap, we investigated diversity at a wide genomic interval surrounding the LCT gene in different Southern Arabian populations. METHODS: 40 SNPs were genotyped to characterize LCT profiles of 630 Omani and Yemeni individuals to perform population structure, linkage disequilibrium, population differentiation-based and haplotype-based analyses. RESULTS: Typical Arabian LP-related variation was found in Dhofaris and Yemenis, being characterized by private haplotypes carrying the -13,915 G allele, unusual differentiation with respect to northern groups and conserved homozygous haplotype-blocks, suggesting that the adaptive allele was likely introduced in the Arabian gene pool in southern populations and was then subjected to prolonged selective pressure. CONCLUSION: By pointing to Yemen as one of the best candidate centers of diffusion of the Arabian-specific adaptive variant, obtained results indicate the spread of indigenous groups as the main process underlying dispersal of LP along the Arabian Peninsula, supporting a refugia model for Arabian demic movements occurred during the Terminal Pleistocene and Early Holocene.

Inferring the genetic history of lactase persistence along the Italian peninsula from a large genomic interval surrounding the LCT gene.

OBJECTIVE: Although genetic variants related to lactase persistence in European populations were supposed to have firstly undergone positive selection in farmers from the Balkans and Central Europe, demographic and evolutionary dynamics that subsequently shaped the distribution of this adaptive trait across the continent have still to be elucidated. To deepen the knowledge about potential routes of diffusion of lactase persistence to Western Europe we investigated variation at a large genomic region surrounding the LCT gene along the Italian peninsula, a geographical area that played a key role in population movements responsible for Neolithic diffusion across Europe. METHODS: By genotyping 40 highly selected SNPs in more than 400 Italian individuals we described gradients of nucleotide and haplotype variation potentially related to lactase persistence and compared them with those observed in several European and Mediterranean human groups. RESULTS: Multiple migratory events responsible for earlier introduction of the examined alleles in Italy than in Northern European regions could be invoked. Different demic processes occurred along the western and eastern sides of the peninsula were also inferred via linkage disequilibrium and population structure analyses. CONCLUSION: The appreciable genetic continuum observed between people from Northern or Central-Western Italy and Central European populations suggested a local arrival of lactase persistence-related variants mainly via overland routes. On the contrary, diversity of Central-Eastern and Southern Italian groups entailed also gene flow from South-Eastern Mediterranean regions, in accordance to the earlier entrance of the Neolithic in Southern Italy via maritime population movements along the Mediterranean coastlines.

The epigenetic side of human adaptation: hypotheses, evidences and theories.

CONTEXT: Epigenetics represents a still unexplored research field in the understanding of micro- and macro-evolutionary mechanisms, as epigenetic changes create phenotypic diversity within both individuals and populations. OBJECTIVE: The purpose of this review is to dissect the landscape of studies focused on DNA methylation, one of the most described epigenetic mechanisms, emphasizing the aspects that could be relevant in human adaptations. METHODS: Theories and results here considered were collected from the most recent papers published. RESULTS: The matter of DNA methylation inheritance is here described as well as the recent evolutionary theories regarding the role of DNA methylation-and epigenetics in a broader sense-in human evolution. The complex relation between (1) DNA methylation and genetic variability and (2) DNA methylation and the environmental stimuli crucial in shaping genetic and phenotypic variability through the human lineage-such as diet, climate and pathogens exposure-are described. Papers about population epigenetics are also illustrated due to their high relevance in this context. CONCLUSION: Genetic, epigenetic and phenotypic variations of the species, together with cultural ones, are considerably shaped by a vast range of environmental stimuli, thus representing the foundation of all human bio-cultural adaptations.

The nucleolar size is associated to the methylation status of ribosomal DNA in breast carcinomas.

BACKGROUND: There is a body of evidence that shows a link between tumorigenesis and ribosome biogenesis. The precursor of mature 18S, 28S and 5.8S ribosomal RNAs is transcribed from the ribosomal DNA gene (rDNA), which exists as 300-400 copies in the human diploid genome. Approximately one half of these copies are epigenetically silenced, but the exact role of epigenetic regulation on ribosome biogenesis is not completely understood. In this study we analyzed the methylation profiles of the rDNA promoter and of the 5' regions of 18S and 28S in breast cancer. METHODS: We analyzed rDNA methylation in 68 breast cancer tissues of which the normal counterpart was partially available (45/68 samples) using the MassARRAY EpiTYPER assay, a sensitive and quantitative method with single base resolution. RESULTS: We found that rDNA locus tended to be hypermethylated in tumor compared to matched normal breast tissues and that the DNA methylation level of several CpG units within the rDNA locus was associated to nuclear grade and to nucleolar size of tumor tissues. In addition we identified a subgroup of samples in which large nucleoli were associated with very limited or absent rDNA hypermethylation in tumor respect to matched normal tissue. CONCLUSIONS: In conclusion, we suggest that rDNA is an important target of epigenetic regulation in breast tumors and that rDNA methylation level is associated to nucleolar size.

Genetic diversity of variants involved in drug response among Tunisian and Italian populations toward personalized medicine.

Adverse drug reactions (ADR) represent a significant contributor to morbidity and mortality, imposing a substantial financial burden. Genetic ancestry plays a crucial role in drug response. The aim of this study is to characterize the genetic variability of selected pharmacogenes involved with ADR in Tunisians and Italians, with a comparative analysis against global populations. A cohort of 135 healthy Tunisians and 737 Italians were genotyped using a SNP array. Variants located in 25 Very Important Pharmacogenes implicated in ADR were extracted from the genotyping data. Distribution analysis of common variants in Tunisian and Italian populations in comparison to 24 publicly available worldwide populations was performed using PLINK and R software. Results from Principle Component and ADMIXTURE analyses showed a high genetic similarity among Mediterranean populations, distinguishing them from Sub-Saharan African and Asian populations. The Fst comparative analysis identified 27 variants exhibiting significant differentiation between the studied populations. Among these variants, four SNPs rs622342, rs3846662, rs7294, rs5215 located in SLC22A1, HMGCR, VKORC1 and KCNJ11 genes respectively, are reported to be associated with ethnic variability in drug responses. In conclusion, correlating the frequencies of genotype risk variants with their associated ADRs would enhance drug outcomes and the implementation of personalized medicine in the studied populations.

Association between Variants of the TRPV1 Gene and Body Composition in Sub-Saharan Africans.

In humans, the transient receptor potential vanilloid 1 (TRPV1) gene is activated by exogenous (e.g., high temperatures, irritating compounds such as capsaicin) and endogenous (e.g., endocannabinoids, inflammatory factors, fatty acid metabolites, low pH) stimuli. It has been shown to be involved in several processes including nociception, thermosensation, and energy homeostasis. In this study, we investigated the association between TRPV1 gene variants, sensory perception (to capsaicin and PROP), and body composition (BMI and bioimpedance variables) in human populations. By comparing sequences deposited in worldwide databases, we identified two haplotype blocks (herein referred to as H1 and H2) that show strong stabilizing selection signals (MAF approaching 0.50, Tajima's D > +4.5) only in individuals with sub-Saharan African ancestry. We therefore studied the genetic variants of these two regions in 46 volunteers of sub-Saharan descent and 45 Italian volunteers (both sexes). Linear regression analyses showed significant associations between TRPV1 diplotypes and body composition, but not with capsaicin perception. Specifically, in African women carrying the H1-b and H2-b haplotypes, a higher percentage of fat mass and lower extracellular fluid retention was observed, whereas no significant association was found in men. Our results suggest the possible action of sex-driven balancing selection at the non-coding sequences of the TRPV1 gene, with adaptive effects on water balance and lipid deposition.

Epigenetics of the far northern Yakutian population.

BACKGROUND: Yakuts are one of the indigenous populations of the subarctic and arctic territories of Siberia characterized by a continental subarctic climate with severe winters, with the regular January average temperature in the regional capital city of Yakutsk dipping below - 40 °C. The epigenetic mechanisms of adaptation to such ecologies and environments and, in particular, epigenetic age acceleration in the local population have not been studied before. RESULTS: This work reports the first epigenetic study of the Yakutian population using whole-blood DNA methylation data, supplemented with the comparison to the residents of Central Russia. Gene set enrichment analysis revealed, among others, geographic region-specific differentially methylated regions associated with adaptation to climatic conditions (water consumption, digestive system regulation), aging processes (actin filament activity, cell fate), and both of them (channel activity, regulation of steroid and corticosteroid hormone secretion). Further, it is demonstrated that the epigenetic age acceleration of the Yakutian representatives is significantly higher than that of Central Russia counterparts. For both geographic regions, we showed that epigenetically males age faster than females, whereas no significant sex differences were found between the regions. CONCLUSIONS: We performed the first study of the epigenetic data of the Yakutia cohort, paying special attention to region-specific features, aging processes, age acceleration, and sex specificity.

Epigenetic aging differences between Wichí and Criollos from Argentina: Insights from genomic history and ecology.

Background and objectives: Epigenetic estimators based on DNA methylation levels have emerged as promising biomarkers of human aging. These estimators exhibit natural variations across human groups, but data about indigenous populations remain underrepresented in research. This study aims to investigate differences in epigenetic estimators between two distinct human populations, both residing in the Gran Chaco region of Argentina, the Native-American Wichí, and admixed Criollos who are descendants of intermarriages between Native Americans and the first European colonizers, using a population genetic approach. Methodology: We analyzed 24 Wichí (mean age: 39.2 ± 12.9 yo) and 24 Criollos (mean age: 41.1 ± 14.0 yo) for DNA methylation levels using the Infinium MethylationEPIC (Illumina) to calculate 16 epigenetic estimators. Additionally, we examined genome-wide genetic variation using the HumanOmniExpress BeadChip (Illumina) to gain insights into the genetic history of these populations. Results: Our results indicate that Native-American Wichí are epigenetically older compared to Criollos according to five epigenetic estimators. Analyses within the Criollos population reveal that global ancestry does not influence the differences observed, while local (chromosomal) ancestry shows positive associations between specific SNPs located in genomic regions over-represented by Native-American ancestry and measures of epigenetic age acceleration (AgeAccelHannum). Furthermore, we demonstrate that differences in population ecologies also contribute to observed epigenetic differences. Conclusions and implications: Overall, our study suggests that while the genomic history may partially account for the observed epigenetic differences, non-genetic factors, such as lifestyle and ecological factors, play a substantial role in the variability of epigenetic estimators, thereby contributing to variations in human epigenetic aging.

Evolutionary Implications of Environmental Toxicant Exposure.

Homo sapiens have been exposed to various toxins and harmful compounds that change according to various phases of human evolution. Population genetics studies showed that such exposures lead to adaptive genetic changes; while observing present exposures to different toxicants, the first molecular mechanism that confers plasticity is epigenetic remodeling and, in particular, DNA methylation variation, a molecular mechanism proposed for medium-term adaptation. A large amount of scientific literature from clinical and medical studies revealed the high impact of such exposure on human biology; thus, in this review, we examine and infer the impact that different environmental toxicants may have in shaping human evolution. We first describe how environmental toxicants shape natural human variation in terms of genetic and epigenetic diversity, and then we describe how DNA methylation may influence mutation rate and, thus, genetic variability. We describe the impact of these substances on biological fitness in terms of reproduction and survival, and in conclusion, we focus on their effect on brain evolution and physiology.