Liver-specific deletion of de novo DNA methyltransferases protects against glucose intolerance in high-fat diet-fed male mice.

Alterations to gene transcription and DNA methylation are a feature of many liver diseases including fatty liver disease and liver cancer. However, it is unclear whether the DNA methylation changes are a cause or a consequence of the transcriptional changes. It is even possible that the methylation changes are not required for the transcriptional changes. If DNA methylation is just a minor player in, or a consequence of liver transcriptional change, then future studies in this area should focus on other systems such as histone tail modifications. To interrogate the importance of de novo DNA methylation, we generated mice that are homozygous mutants for both Dnmt3a and Dnmt3b in post-natal liver. These mice are viable and fertile with normal sized livers. Males, but not females, showed increased adipose depots, yet paradoxically, improved glucose tolerance on both control diet and high-fat diets (HFD). Comparison of the transcriptome and methylome with RNA sequencing and whole-genome bisulfite sequencing in adult hepatocytes revealed that widespread loss of methylation in CpG-rich regions in the mutant did not induce loss of homeostatic transcriptional regulation. Similarly, extensive transcriptional changes induced by HFD did not require de novo DNA methylation. The improved metabolic phenotype of the Dnmt3a/3b mutant mice may be mediated through the dysregulation of a subset of glucose and fat metabolism genes which increase both glucose uptake and lipid export by the liver. However, further work is needed to confirm this.

A DNA methylation biomarker of alcohol consumption.

The lack of reliable measures of alcohol intake is a major obstacle to the diagnosis and treatment of alcohol-related diseases. Epigenetic modifications such as DNA methylation may provide novel biomarkers of alcohol use. To examine this possibility, we performed an epigenome-wide association study of methylation of cytosine-phosphate-guanine dinucleotide (CpG) sites in relation to alcohol intake in 13 population-based cohorts (ntotal=13 317; 54% women; mean age across cohorts 42-76 years) using whole blood (9643 European and 2423 African ancestries) or monocyte-derived DNA (588 European, 263 African and 400 Hispanic ancestry) samples. We performed meta-analysis and variable selection in whole-blood samples of people of European ancestry (n=6926) and identified 144 CpGs that provided substantial discrimination (area under the curve=0.90-0.99) for current heavy alcohol intake (⩾42 g per day in men and ⩾28 g per day in women) in four replication cohorts. The ancestry-stratified meta-analysis in whole blood identified 328 (9643 European ancestry samples) and 165 (2423 African ancestry samples) alcohol-related CpGs at Bonferroni-adjusted P<1 × 10-7. Analysis of the monocyte-derived DNA (n=1251) identified 62 alcohol-related CpGs at P<1 × 10-7. In whole-blood samples of people of European ancestry, we detected differential methylation in two neurotransmitter receptor genes, the γ-Aminobutyric acid-A receptor delta and γ-aminobutyric acid B receptor subunit 1; their differential methylation was associated with expression levels of a number of genes involved in immune function. In conclusion, we have identified a robust alcohol-related DNA methylation signature and shown the potential utility of DNA methylation as a clinically useful diagnostic test to detect current heavy alcohol consumption.

An epigenome-wide association study meta-analysis of educational attainment.

The epigenome is associated with biological factors, such as disease status, and environmental factors, such as smoking, alcohol consumption and body mass index. Although there is a widespread perception that environmental influences on the epigenome are pervasive and profound, there has been little evidence to date in humans with respect to environmental factors that are biologically distal. Here we provide evidence on the associations between epigenetic modifications-in our case, CpG methylation-and educational attainment (EA), a biologically distal environmental factor that is arguably among the most important life-shaping experiences for individuals. Specifically, we report the results of an epigenome-wide association study meta-analysis of EA based on data from 27 cohort studies with a total of 10 767 individuals. We find nine CpG probes significantly associated with EA. However, robustness analyses show that all nine probes have previously been found to be associated with smoking. Only two associations remain when we perform a sensitivity analysis in the subset of never-smokers, and these two probes are known to be strongly associated with maternal smoking during pregnancy, and thus their association with EA could be due to correlation between EA and maternal smoking. Moreover, the effect sizes of the associations with EA are far smaller than the known associations with the biologically proximal environmental factors alcohol consumption, body mass index, smoking and maternal smoking during pregnancy. Follow-up analyses that combine the effects of many probes also point to small methylation associations with EA that are highly correlated with the combined effects of smoking. If our findings regarding EA can be generalized to other biologically distal environmental factors, then they cast doubt on the hypothesis that such factors have large effects on the epigenome.

708 Common and 2010 rare DISC1 locus variants identified in 1542 subjects: analysis for association with psychiatric disorder and cognitive traits.

A balanced t(1;11) translocation that transects the Disrupted in schizophrenia 1 (DISC1) gene shows genome-wide significant linkage for schizophrenia and recurrent major depressive disorder (rMDD) in a single large Scottish family, but genome-wide and exome sequencing-based association studies have not supported a role for DISC1 in psychiatric illness. To explore DISC1 in more detail, we sequenced 528 kb of the DISC1 locus in 653 cases and 889 controls. We report 2718 validated single-nucleotide polymorphisms (SNPs) of which 2010 have a minor allele frequency of <1%. Only 38% of these variants are reported in the 1000 Genomes Project European subset. This suggests that many DISC1 SNPs remain undiscovered and are essentially private. Rare coding variants identified exclusively in patients were found in likely functional protein domains. Significant region-wide association was observed between rs16856199 and rMDD (P=0.026, unadjusted P=6.3 × 10(-5), OR=3.48). This was not replicated in additional recurrent major depression samples (replication P=0.11). Combined analysis of both the original and replication set supported the original association (P=0.0058, OR=1.46). Evidence for segregation of this variant with disease in families was limited to those of rMDD individuals referred from primary care. Burden analysis for coding and non-coding variants gave nominal associations with diagnosis and measures of mood and cognition. Together, these observations are likely to generalise to other candidate genes for major mental illness and may thus provide guidelines for the design of future studies.

Monozygotic twins affected with major depressive disorder have greater variance in methylation than their unaffected co-twin.

Our understanding of major depressive disorder (MDD) has focused on the influence of genetic variation and environmental risk factors. Growing evidence suggests the additional role of epigenetic mechanisms influencing susceptibility for complex traits. DNA sequence within discordant monozygotic twin (MZT) pairs is virtually identical; thus, they represent a powerful design for studying the contribution of epigenetic factors to disease liability. The aim of this study was to investigate whether specific methylation profiles in white blood cells could contribute to the aetiology of MDD. Participants were drawn from the Queensland Twin Registry and comprised 12 MZT pairs discordant for MDD and 12 MZT pairs concordant for no MDD and low neuroticism. Bisulphite treatment and genome-wide interrogation of differentially methylated CpG sites using the Illumina Human Methylation 450 BeadChip were performed in WBC-derived DNA. No overall difference in mean global methylation between cases and their unaffected co-twins was found; however, the differences in females was significant (P=0.005). The difference in variance across all probes between affected and unaffected twins was highly significant (P<2.2 × 10⁻¹6), with 52.4% of probes having higher variance in cases (binomial P-value<2.2 × 10⁻¹6). No significant differences in methylation were observed between discordant MZT pairs and their matched concordant MZT (permutation minimum P=0.11) at any individual probe. Larger samples are likely to be needed to identify true associations between methylation differences at specific CpG sites.

Horn type and horn length genes map to the same chromosomal region in Soay sheep.

The evolution of male weaponry in animals is driven by sexual selection, which is predicted to reduce the genetic variability underlying such traits. Soay sheep have an inherited polymorphism for horn type in both sexes, with males presenting with either large, normal horns or small, deformed horns (scurs). In addition, there is additive genetic variation in horn length among males with normal horns. Given that scurred males cannot win conflicts with normal-horned males, it is unusual that genes conferring scurs should persist in the population. Identifying the genetic basis of these traits should help us in understanding their evolution. We developed microsatellite markers in a targeted region of the Soay sheep genome and refined the location of the Horns locus (Ho) to a approximately 7.4 cM interval on chromosome 10 (LOD=8.78). We then located quantitative trait loci spanning a 34 cM interval with a peak centred close to Ho, which explained the majority of the genetic variation for horn length and base circumference in normal-horned males (LOD=2.51 and LOD=1.04, respectively). Therefore, the genetic variation in both horn type and horn length is attributable to the same chromosomal region. Understanding the maintenance of horn type and length variation will require an investigation of selection on genotypes that (co)determine both traits.

Family-based mitochondrial association study of traits related to type 2 diabetes and the metabolic syndrome in adolescents.

AIMS/HYPOTHESIS: There has been much focus on the potential role of mitochondria in the aetiology of type 2 diabetes and the metabolic syndrome, and many case-control mitochondrial association studies have been undertaken for these conditions. We tested for a potential association between common mitochondrial variants and a number of quantitative traits related to type 2 diabetes in a large sample of >2,000 healthy Australian adolescent twins and their siblings, many of whom were measured on more than one occasion. METHODS: To the best of our knowledge, this is the first mitochondrial association study of quantitative traits undertaken using family data. The maternal inheritance pattern of mitochondria means established association methodologies are unsuitable for analysis of mitochondrial data in families. We present a methodology, implemented in the freely available program Sib-Pair for performing such an analysis. RESULTS: Despite our study having the power to detect variants with modest effects on these phenotypes, only one significant association was found after correction for multiple testing in any of four age groups. This was for mt14365 with triacylglycerol levels (unadjusted p = 0.0006). This association was not replicated in other age groups. CONCLUSIONS/INTERPRETATION: We find little evidence in our sample to suggest that common European mitochondrial variants contribute to variation in quantitative phenotypes related to diabetes. Only one variant showed a significant association in our sample, and this association will need to be replicated in a larger cohort. Such replication studies or future meta-analyses may reveal more subtle effects that could not be detected here because of limitations of sample size.

A localized negative genetic correlation constrains microevolution of coat color in wild sheep.

The evolutionary changes that occur over a small number of generations in natural populations often run counter to what is expected on the basis of the heritability of traits and the selective forces acting upon them. In Soay sheep, dark coat color is associated with large size, which is heritable and positively correlated with fitness, yet the frequency of dark sheep has decreased. This unexpected microevolutionary trend is explained by genetic linkage between the causal mutation underlying the color polymorphism and quantitative trait loci with antagonistic effects on size and fitness. As a consequence, homozygous dark sheep are large, but have reduced fitness relative to phenotypically indistinguishable dark heterozygotes and light sheep. This result demonstrates the importance of understanding the genetic basis of fitness variation when making predictions about the microevolutionary consequences of selection.

Compelling evidence that a single nucleotide substitution in TYRP1 is responsible for coat-colour polymorphism in a free-living population of Soay sheep.

Identifying the genes that underlie phenotypic variation in natural populations is a central objective of evolutionary genetics. Here, we report the identification of the gene and causal mutation underlying coat colour variation in a free-living population of Soay sheep (Ovis aries). We targeted tyrosinase-related protein 1 (TYRP1), a positional candidate gene based on previous work that mapped the Coat colour locus to an approximately 15cM window on chromosome 2. We identified a non-synonymous substitution in exon IV that was perfectly associated with coat colour. This polymorphism is predicted to cause the loss of a cysteine residue that is highly evolutionarily conserved and likely to be of functional significance. We eliminated the possibility that this association is due to the presence of strong linkage disequilibrium with an unknown regulatory mutation by demonstrating that there is no difference in relative TYRP1 expression between colour morphs. Analysis of this putative causal mutation in a complex pedigree of more than 500 sheep revealed almost perfect co-segregation with coat colour (chi2-test, p<0.0001, LOD=110.20), and very tight linkage between Coat colour and TYRP1 (LOD=29.50).

Bone density in sheep: genetic variation and quantitative trait loci localisation.

Bone density (BD) is an important factor in osteoporotic fracture risk in humans. However, BD is a complex trait confounded by environmental influences and polygenic inheritance. Sheep provide a potentially useful model for studying differences in BD, as they provide a means of circumventing complex environmental factors and are a similar weight to humans. The aims of this study were to establish whether there is genetic variation in BD in sheep and then to localise quantitative trait loci (QTLs) associated with this variation. We also aimed to evaluate the relationship between fat and muscle body components and BD in sheep. Results showed that there was significant (P < 0.01) genetic variation among Coopworth sheep sires for BD. This genetic difference was correlated (P < 0.01) with body weight and muscle mass. A number of QTLs exceeding the suggestive threshold were identified (nine in total). Of these, two (chromosomes 1, P < 0.05; chromosome 24, P < 0.01) were significant using genome-wide permutation significance thresholds (2000 iterations). The position of the QTL on chromosome 24 coincided with a number of other body composition QTLs, indicating possible pleiotropic effects or the presence of multiple genes affecting body composition at that site. This study shows that sheep are potentially a useful model for studying the genetics of BD.

Linkage disequilibrium in domestic sheep.

The last decade has seen a dramatic increase in the number of livestock QTL mapping studies. The next challenge awaiting livestock geneticists is to determine the actual genes responsible for variation of economically important traits. With the advent of high density single nucleotide polymorphism (SNP) maps, it may be possible to fine map genes by exploiting linkage disequilibrium between genes of interest and adjacent markers. However, the extent of linkage disequilibrium (LD) is generally unknown for livestock populations. In this article microsatellite genotype data are used to assess the extent of LD in two populations of domestic sheep. High levels of LD were found to extend for tens of centimorgans and declined as a function of marker distance. However, LD was also frequently observed between unlinked markers. The prospects for LD mapping in livestock appear encouraging provided that type I error can be minimized. Properties of the multiallelic LD coefficient D' were also explored. D' was found to be significantly related to marker heterozygosity, although the relationship did not appear to unduly influence the overall conclusions. Of potentially greater concern was the observation that D' may be skewed when rare alleles are present. It is recommended that the statistical significance of LD is used in conjunction with coefficients such as D' to determine the true extent of LD.