Locomotion is a defining characteristic of animal life and plays a crucial role in most behaviors. Locomotion involves physical activity, which can have far-reaching effects on physiology and neurobiology, both acutely and chronically. In human populations and in laboratory rodents, higher levels of physical activity are generally associated with positive health outcomes, although excessive exercise can have adverse consequences. Whether and how such relationships occur in wild animals is unknown. Behavioral variation among individuals arises from genetic and environmental factors and their interactions as well as from developmental programming (persistent effects of early-life environment). Although tremendous progress has been made in identifying genetic and environmental influences on individual differences in behavior, early-life effects are not well understood. Early-life effects can in some cases persist across multiple generations following a single exposure and, in principle, may constrain or facilitate the rate of evolution at multiple levels of biological organization. Understanding the mechanisms of such transgenerational effects (e.g., exposure to stress hormones in utero, inherited epigenetic alterations) may prove crucial to explaining unexpected and/or sex-specific responses to selection as well as limits to adaptation. One area receiving increased attention is early-life effects on adult physical activity. Correlational data from epidemiological studies suggest that early-life nutritional stress can (adversely) affect adult human activity levels and associated physiological traits (e.g., body composition, metabolic health). The few existing studies of laboratory rodents demonstrate that both maternal and early-life exercise can affect adult levels of physical activity and related phenotypes. Going forward, rodents offer many opportunities for experimental studies of (multigenerational) early-life effects, including studies that use maternal exposures and cross-fostering designs.

DNA methylation regulates cell type-specific gene expression. Here, in a transgenic mouse model, we show that deletion of the gene encoding DNA methyltransferase Dnmt3a in hypothalamic AgRP neurons causes a sedentary phenotype characterized by reduced voluntary exercise and increased adiposity. Whole-genome bisulfite sequencing (WGBS) and transcriptional profiling in neuronal nuclei from the arcuate nucleus of the hypothalamus (ARH) reveal differentially methylated genomic regions and reduced expression of AgRP neuron-associated genes in knockout mice. We use read-level analysis of WGBS data to infer putative ARH neural cell types affected by the knockout, and to localize promoter hypomethylation and increased expression of the growth factor Bmp7 to AgRP neurons, suggesting a role for aberrant TGF-β signaling in the development of this phenotype. Together, these data demonstrate that DNA methylation in AgRP neurons is required for their normal epigenetic development and neuron-specific gene expression profiles, and regulates voluntary exercise behavior.

Background

Genetic variants can modulate phenotypic outcomes via epigenetic intermediates, for example at methylation quantitative trait loci (mQTL). We present the first large-scale assessment of mQTL at human genomic regions selected for interindividual variation in CpG methylation, which we call correlated regions of systemic interindividual variation (CoRSIVs). These can be assayed in blood DNA and do not reflect interindividual variation in cellular composition.

Results

We use target-capture bisulfite sequencing to assess DNA methylation at 4086 CoRSIVs in multiple tissues from each of 188 donors in the NIH Gene-Tissue Expression (GTEx) program. At CoRSIVs, DNA methylation in peripheral blood correlates with methylation and gene expression in internal organs. We also discover unprecedented mQTL at these regions. Genetic influences on CoRSIV methylation are extremely strong (median R²=0.76), cumulatively comprising over 70-fold more human mQTL than detected in the most powerful previous study. Moreover, mQTL beta coefficients at CoRSIVs are highly skewed (i.e., the major allele predicts higher methylation). Both surprising findings are independently validated in a cohort of 47 non-GTEx individuals. Genomic regions flanking CoRSIVs show long-range enrichments for LINE-1 and LTR transposable elements; the skewed beta coefficients may therefore reflect evolutionary selection of genetic variants that promote their methylation and silencing. Analyses of GWAS summary statistics show that mQTL polymorphisms at CoRSIVs are associated with metabolic and other classes of disease.

Conclusions

A focus on systemic interindividual epigenetic variants, clearly enhanced in mQTL content, should likewise benefit studies attempting to link human epigenetic variation to the risk of disease.

The double burden of under- and overnutrition profoundly affects human health globally. According to the World Health Organization, obesity and diabetes rates have almost doubled worldwide since 1980, and, in 2011, more than 40 million children under 5 years of age were overweight. Ecologic factors, parental genetics and fitness, and the intrauterine environment significantly influence the likelihood of offspring developing the dysmetabolic diathesis of obesity. This report examines the effects of these factors, including preconception, intrauterine and postnatal energy balance affecting programming of transgenerational transmission, and development of chronic diseases later in life-in particular, diabesity and its comorbidities.

Analysis of DNA methylation patterns relies increasingly on sequencing-based profiling methods. The four most frequently used sequencing-based technologies are the bisulfite-based methods MethylC-seq and reduced representation bisulfite sequencing (RRBS), and the enrichment-based techniques methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylated DNA binding domain sequencing (MBD-seq). We applied all four methods to biological replicates of human embryonic stem cells to assess their genome-wide CpG coverage, resolution, cost, concordance and the influence of CpG density and genomic context. The methylation levels assessed by the two bisulfite methods were concordant (their difference did not exceed a given threshold) for 82% for CpGs and 99% of the non-CpG cytosines. Using binary methylation calls, the two enrichment methods were 99% concordant and regions assessed by all four methods were 97% concordant. We combined MeDIP-seq with methylation-sensitive restriction enzyme (MRE-seq) sequencing for comprehensive methylome coverage at lower cost. This, along with RNA-seq and ChIP-seq of the ES cells enabled us to detect regions with allele-specific epigenetic states, identifying most known imprinted regions and new loci with monoallelic epigenetic marks and monoallelic expression.