DNA methylation is an essential epigenetic mark whose role in gene regulation and its dependency on genomic sequence and environment are not fully understood. In this study we provide novel insights into the mechanistic relationships between genetic variation, DNA methylation and transcriptome sequencing data in three different cell-types of the GenCord human population cohort. We find that the association between DNA methylation and gene expression variation among individuals are likely due to different mechanisms from those establishing methylation-expression patterns during differentiation. Furthermore, cell-type differential DNA methylation may delineate a platform in which local inter-individual changes may respond to or act in gene regulation. We show that unlike genetic regulatory variation, DNA methylation alone does not significantly drive allele specific expression. Finally, inferred mechanistic relationships using genetic variation as well as correlations with TF abundance reveal both a passive and active role of DNA methylation to regulatory interactions influencing gene expression.DOI: http://dx.doi.org/10.7554/eLife.00523.001

Schizophrenia is a common disorder characterized by psychotic symptoms; diagnostic criteria have been established. Family, twin and adoption studies suggest that both genetic and environmental factors influence susceptibility (heritability is approximately 71%; ref. 2), however, little is known about the aetiology of schizophrenia. Clinical and family studies suggest aetiological heterogeneity. Previously, we reported that regions on chromosomes 22, 3 and 8 may be associated with susceptibility to schizophrenia, and collaborations provided some support for regions on chromosomes 8 and 22 (refs 9-13). We present here a genome-wide scan for schizophrenia susceptibility loci (SSL) using 452 microsatellite markers on 54 multiplex pedigrees. Non-parametric linkage (NPL) analysis provided significant evidence for an SSL on chromosome 13q32 (NPL score=4.18; P=0.00002), and suggestive evidence for another SSL on chromosome 8p21-22 (NPL=3.64; P=0.0001). Parametric linkage analysis provided additional support for these SSL. Linkage evidence at chromosome 8 is weaker than that at chromosome 13, so it is more probable that chromosome 8 may be a false positive linkage. Additional putative SSL were noted on chromosomes 14q13 (NPL=2.57; P=0.005), 7q11 (NPL=2.50, P=0.007) and 22q11 (NPL=2.42, P=0.009). Verification of suggestive SSL on chromosomes 13q and 8p was attempted in a follow-up sample of 51 multiplex pedigrees. This analysis confirmed the SSL in 13q14-q33 (NPL=2.36, P=0.007) and supported the SSL in 8p22-p21 (NPL=1.95, P=0.023).

Trisomy 21 is the most frequent genetic cause of cognitive impairment. To assess the perturbations of gene expression in trisomy 21, and to eliminate the noise of genomic variability, we studied the transcriptome of fetal fibroblasts from a pair of monozygotic twins discordant for trisomy 21. Here we show that the differential expression between the twins is organized in domains along all chromosomes that are either upregulated or downregulated. These gene expression dysregulation domains (GEDDs) can be defined by the expression level of their gene content, and are well conserved in induced pluripotent stem cells derived from the twins' fibroblasts. Comparison of the transcriptome of the Ts65Dn mouse model of Down's syndrome and normal littermate mouse fibroblasts also showed GEDDs along the mouse chromosomes that were syntenic in human. The GEDDs correlate with the lamina-associated (LADs) and replication domains of mammalian cells. The overall position of LADs was not altered in trisomic cells; however, the H3K4me3 profile of the trisomic fibroblasts was modified and accurately followed the GEDD pattern. These results indicate that the nuclear compartments of trisomic cells undergo modifications of the chromatin environment influencing the overall transcriptome, and that GEDDs may therefore contribute to some trisomy 21 phenotypes.

We performed exome sequencing to detect somatic mutations in protein-coding regions in seven melanoma cell lines and donor-matched germline cells. All melanoma samples had high numbers of somatic mutations, which showed the hallmark of UV-induced DNA repair. Such a hallmark was absent in tumor sample-specific mutations in two metastases derived from the same individual. Two melanomas with non-canonical BRAF mutations harbored gain-of-function MAP2K1 and MAP2K2 (MEK1 and MEK2, respectively) mutations, resulting in constitutive ERK phosphorylation and higher resistance to MEK inhibitors. Screening a larger cohort of individuals with melanoma revealed the presence of recurring somatic MAP2K1 and MAP2K2 mutations, which occurred at an overall frequency of 8%. Furthermore, missense and nonsense somatic mutations were frequently found in three candidate melanoma genes, FAT4, LRP1B and DSC1.

Down syndrome (DS) is one of the most frequent congenital birth defects, and the most common genetic cause of mental retardation. In most cases, DS results from the presence of an extra copy of chromosome 21. DS has a complex phenotype, and a major goal of DS research is to identify genotype -phenotype correlations. Cases of partial trisomy 21 and other HSA21 rearrangements associated with DS features could identify genomic regions associated with specific phenotypes. We have developed a BAC array spanning HSA21q and used array comparative genome hybridization (aCGH) to enable high-resolution mapping of pathogenic partial aneuploidies and unbalanced translocations involving HSA21. We report the identification and mapping of 30 pathogenic chromosomal aberrations of HSA21 consisting of 19 partial trisomies and 11 partial monosomies for different segments of HSA21. The breakpoints have been mapped to within B85 kb. The majority of the breakpoints (26 of 30) for the partial aneuploidies map within a 10-Mb region. Our data argue against a single DS critical region. We identify susceptibility regions for 25 phenotypes for DS and 27 regions for monosomy 21. However, most of these regions are still broad, and more cases are needed to narrow down the phenotypic maps to a reasonable number of candidate genomic elements per phenotype.