Single cell transcriptome in aneuploidies reveals mechanisms of gene dosage imbalance.

Aneuploidy is a major source of gene dosage imbalance due to copy number alterations (CNA), and viable human trisomies are model disorders of altered gene expression. We study gene and allele-specific expression (ASE) of 9668 single-cell fibroblasts from trisomy 21 (T21) discordant twins and from mosaic T21, T18, T13 and T8. We examine 928 single cells with deep scRNAseq. Expected and observed overexpression of trisomic genes in trisomic vs. diploid bulk RNAseq is not detectable in trisomic vs. diploid single cells. Instead, for trisomic genes with low-to-average expression, their altered gene dosage is mainly due to the higher fraction of trisomic cells simultaneously expressing these genes, in agreement with a stochastic 2-state burst-like model of transcription. These results, confirmed in a further analysis of 8740 single fibroblasts with shallow scRNAseq, suggest that the specific transcriptional profile of each gene contributes to the phenotypic variability of trisomies. We propose an improved model to understand the effects of CNA and, generally, of gene regulation on gene dosage imbalance.

Bi-allelic Variants in IQSEC1 Cause Intellectual Disability, Developmental Delay, and Short Stature.

We report two consanguineous families with probands that exhibit intellectual disability, developmental delay, short stature, aphasia, and hypotonia in which homozygous non-synonymous variants were identified in IQSEC1 (GenBank: NM_001134382.3). In a Pakistani family, the IQSEC1 segregating variant is c.1028C>T (p.Thr343Met), while in a Saudi Arabian family the variant is c.962G>A (p.Arg321Gln). IQSEC1-3 encode guanine nucleotide exchange factors for the small GTPase ARF6 and their loss affects a variety of actin-dependent cellular processes, including AMPA receptor trafficking at synapses. The ortholog of IQSECs in the fly is schizo and its loss affects growth cone guidance at the midline in the CNS, also an actin-dependent process. Overexpression of the reference IQSEC1 cDNA in wild-type flies is lethal, but overexpression of the two variant IQSEC1 cDNAs did not affect viability. Loss of schizo caused embryonic lethality that could be rescued to 2nd instar larvae by moderate expression of the human reference cDNA. However, the p.Arg321Gln and p.Thr343Met variants failed to rescue embryonic lethality. These data indicate that the variants behave as loss-of-function mutations. We also show that schizo in photoreceptors is required for phototransduction. Finally, mice with a conditional Iqsec1 deletion in cortical neurons exhibited an increased density of dendritic spines with an immature morphology. The phenotypic similarity of the affecteds and the functional experiments in flies and mice indicate that IQSEC1 variants are the cause of a recessive disease with intellectual disability, developmental delay, and short stature, and that axonal guidance and dendritic projection defects as well as dendritic spine dysgenesis may underlie disease pathogenesis.

Bi-allelic Variants in DYNC1I2 Cause Syndromic Microcephaly with Intellectual Disability, Cerebral Malformations, and Dysmorphic Facial Features.

Cargo transport along the cytoplasmic microtubular network is essential for neuronal function, and cytoplasmic dynein-1 is an established molecular motor that is critical for neurogenesis and homeostasis. We performed whole-exome sequencing, homozygosity mapping, and chromosomal microarray studies in five individuals from three independent pedigrees and identified likely-pathogenic variants in DYNC1I2 (Dynein Cytoplasmic 1 Intermediate Chain 2), encoding a component of the cytoplasmic dynein 1 complex. In a consanguineous Pakistani family with three affected individuals presenting with microcephaly, severe intellectual disability, simplification of cerebral gyration, corpus callosum hypoplasia, and dysmorphic facial features, we identified a homozygous splice donor site variant (GenBank: NM_001378.2:c.607+1G>A). We report two additional individuals who have similar neurodevelopmental deficits and craniofacial features and harbor deleterious variants; one individual bears a c.740A>G (p.Tyr247Cys) change in trans with a 374 kb deletion encompassing DYNC1I2, and an unrelated individual harbors the compound-heterozygous variants c.868C>T (p.Gln290∗) and c.740A>G (p.Tyr247Cys). Zebrafish larvae subjected to CRISPR-Cas9 gene disruption or transient suppression of dync1i2a displayed significantly altered craniofacial patterning with concomitant reduction in head size. We monitored cell death and cell cycle progression in dync1i2a zebrafish models and observed significantly increased apoptosis, likely due to prolonged mitosis caused by abnormal spindle morphology, and this finding offers initial insights into the cellular basis of microcephaly. Additionally, complementation studies in zebrafish demonstrate that p.Tyr247Cys attenuates gene function, consistent with protein structural analysis. Our genetic and functional data indicate that DYNC1I2 dysfunction probably causes an autosomal-recessive microcephaly syndrome and highlight further the critical roles of the dynein-1 complex in neurodevelopment.

Biallelic variants in FBXL3 cause intellectual disability, delayed motor development and short stature.

FBXL3 (F-Box and Leucine Rich Repeat Protein 3) encodes a protein that contains an F-box and several tandem leucine-rich repeats (LRR) domains. FBXL3 is part of the SCF (Skp1-Cullin-F box protein) ubiquitin ligase complex that binds and leads to phosphorylation-dependent degradation of the central clock protein cryptochromes (CRY1 and CRY2) by the proteasome and its absence causes circadian phenotypes in mice and behavioral problems. No FBXL3-related phenotypes have been described in humans. By a combination of exome sequencing and homozygosity mapping, we analyzed two consanguineous families with intellectual disability and identified homozygous loss-of-function (LoF) variants in FBXL3. In the first family, from Pakistan, an FBXL3 frameshift variant [NM_012158.2:c.885delT:p.(Leu295Phefs*25)] was the onlysegregating variant in five affected individuals in two family loops (LOD score: 3.12). In the second family, from Lebanon, we identified a nonsense variant [NM_012158.2:c.445C>T:p.(Arg149*)]. In a third patient from Italy, a likely deleterious non-synonymous variant [NM_012158.2:c.1072T>C:p.(Cys358Arg)] was identified in homozygosity. Protein 3D modeling predicted that the Cys358Arg change influences the binding with CRY2 by destabilizing the structure of the FBXL3, suggesting that this variant is also likely to be LoF. The eight affected individuals from the three families presented with a similar phenotype that included intellectual disability, developmental delay, short stature and mild facial dysmorphism, mainly large nose with a bulbous tip. The phenotypic similarity and the segregation analysis suggest that FBXL3 biallelic, LoF variants link this gene with syndromic autosomal recessive developmental delay/intellectual disability.

Extensive cellular heterogeneity of X inactivation revealed by single-cell allele-specific expression in human fibroblasts.

X-chromosome inactivation (XCI) provides a dosage compensation mechanism where, in each female cell, one of the two X chromosomes is randomly silenced. However, some genes on the inactive X chromosome and outside the pseudoautosomal regions escape from XCI and are expressed from both alleles (escapees). We investigated XCI at single-cell resolution combining deep single-cell RNA sequencing with whole-genome sequencing to examine allelic-specific expression in 935 primary fibroblast and 48 lymphoblastoid single cells from five female individuals. In this framework we integrated an original method to identify and exclude doublets of cells. In fibroblast cells, we have identified 55 genes as escapees including five undescribed escapee genes. Moreover, we observed that all genes exhibit a variable propensity to escape XCI in each cell and cell type and that each cell displays a distinct expression profile of the escapee genes. A metric, the Inactivation Score-defined as the mean of the allelic expression profiles of the escapees per cell-enables us to discover a heterogeneous and continuous degree of cellular XCI with extremes represented by "inactive" cells, i.e., cells exclusively expressing the escaping genes from the active X chromosome and "escaping" cells expressing the escapees from both alleles. We found that this effect is associated with cell-cycle phases and, independently, with the XIST expression level, which is higher in the quiescent phase (G0). Single-cell allele-specific expression is a powerful tool to identify novel escapees in different tissues and provide evidence of an unexpected cellular heterogeneity of XCI.

Bi-allelic Loss-of-Function Variants in DNMBP Cause Infantile Cataracts.

Infantile and childhood-onset cataracts form a heterogeneous group of disorders; among the many genetic causes, numerous pathogenic variants in additional genes associated with autosomal-recessive infantile cataracts remain to be discovered. We identified three consanguineous families affected by bilateral infantile cataracts. Using exome sequencing, we found homozygous loss-of-function variants in DNMBP: nonsense variant c.811C>T (p.Arg271∗) in large family F385 (nine affected individuals; LOD score = 5.18 at θ = 0), frameshift deletion c.2947_2948del (p.Asp983∗) in family F372 (two affected individuals), and frameshift variant c.2852_2855del (p.Thr951Metfs∗41) in family F3 (one affected individual). The phenotypes of all affected individuals include infantile-onset cataracts. RNAi-mediated knockdown of the Drosophila ortholog still life (sif), enriched in lens-secreting cells, affects the development of these cells as well as the localization of E-cadherin, alters the distribution of septate junctions in adjacent cone cells, and leads to a ∼50% reduction in electroretinography amplitudes in young flies. DNMBP regulates the shape of tight junctions, which correspond to the septate junctions in invertebrates, as well as the assembly pattern of E-cadherin in human epithelial cells. E-cadherin has an important role in lens vesicle separation and lens epithelial cell survival in humans. We therefore conclude that DNMBP loss-of-function variants cause infantile-onset cataracts in humans.

Visual impairment and progressive phthisis bulbi caused by recessive pathogenic variant in MARK3.

Developmental eye defects often severely reduce vision. Despite extensive efforts, for a substantial fraction of these cases the molecular causes are unknown. Recessive eye disorders are frequent in consanguineous populations and such large families with multiple affected individuals provide an opportunity to identify recessive causative genes. We studied a Pakistani consanguineous family with three affected individuals with congenital vision loss and progressive eye degeneration. The family was analyzed by exome sequencing of one affected individual and genotyping of all family members. We have identified a non-synonymous homozygous variant (NM_001128918.2: c.1708C > G: p.Arg570Gly) in the MARK3 gene as the likely cause of the phenotype. Given that MARK3 is highly conserved in flies (I: 55%; S: 67%) we knocked down the MARK3 homologue, par-1, in the eye during development. This leads to a significant reduction in eye size, a severe loss of photoreceptors and loss of vision based on electroretinogram (ERG) recordings. Expression of the par-1 p.Arg792Gly mutation (equivalent to the MARK3 variant found in patients) in developing fly eyes also induces loss of eye tissue and reduces the ERG signals. The data in flies and human indicate that the MARK3 variant corresponds to a loss of function. We conclude that the identified mutation in MARK3 establishes a new gene-disease link, since it likely causes structural abnormalities during eye development and visual impairment in humans, and that the function of MARK3/par-1 is evolutionarily conserved in eye development.

The genomic landscape of human cellular circadian variation points to a novel role for the signalosome.

The importance of natural gene expression variation for human behavior is undisputed, but its impact on circadian physiology remains mostly unexplored. Using umbilical cord fibroblasts, we have determined by genome-wide association how common genetic variation impacts upon cellular circadian function. Gene set enrichment points to differences in protein catabolism as one major source of clock variation in humans. The two most significant alleles regulated expression of COPS7B, a subunit of the COP9 signalosome. We further show that the signalosome complex is imported into the nucleus in timed fashion to stabilize the essential circadian protein BMAL1, a novel mechanism to oppose its proteasome-mediated degradation. Thus, circadian clock properties depend in part upon a genetically-encoded competition between stabilizing and destabilizing forces, and genetic alterations in these mechanisms provide one explanation for human chronotype.

Detection of Imprinted Genes by Single-Cell Allele-Specific Gene Expression.

Genomic imprinting results in parental-specific gene expression. Imprinted genes are involved in the etiology of rare syndromes and have been associated with common diseases such as diabetes and cancer. Standard RNA bulk cell sequencing applied to whole-tissue samples has been used to detect imprinted genes in human and mouse models. However, lowly expressed genes cannot be detected by using RNA bulk approaches. Here, we report an original and robust method that combines single-cell RNA-seq and whole-genome sequencing into an optimized statistical framework to analyze genomic imprinting in specific cell types and in different individuals. Using samples from the probands of 2 family trios and 3 unrelated individuals, 1,084 individual primary fibroblasts were RNA sequenced and more than 700,000 informative heterozygous single-nucleotide variations (SNVs) were genotyped. The allele-specific coverage per gene of each SNV in each single cell was used to fit a beta-binomial distribution to model the likelihood of a gene being expressed from one and the same allele. Genes presenting a significant aggregate allelic ratio (between 0.9 and 1) were retained to identify of the allelic parent of origin. Our approach allowed us to validate the imprinting status of all of the known imprinted genes expressed in fibroblasts and the discovery of nine putative imprinted genes, thereby demonstrating the advantages of single-cell over bulk RNA-seq to identify imprinted genes. The proposed single-cell methodology is a powerful tool for establishing a cell type-specific map of genomic imprinting.

HSA21 Single-Minded 2 (Sim2) Binding Sites Co-Localize with Super-Enhancers and Pioneer Transcription Factors in Pluripotent Mouse ES Cells.

The HSA21 encoded Single-minded 2 (SIM2) transcription factor has key neurological functions and is a good candidate to be involved in the cognitive impairment of Down syndrome. We aimed to explore the functional capacity of SIM2 by mapping its DNA binding sites in mouse embryonic stem cells. ChIP-sequencing revealed 1229 high-confidence SIM2-binding sites. Analysis of the SIM2 target genes confirmed the importance of SIM2 in developmental and neuronal processes and indicated that SIM2 may be a master transcription regulator. Indeed, SIM2 DNA binding sites share sequence specificity and overlapping domains of occupancy with master transcription factors such as SOX2, OCT4 (Pou5f1), NANOG or KLF4. The association between SIM2 and these pioneer factors is supported by co-immunoprecipitation of SIM2 with SOX2, OCT4, NANOG or KLF4. Furthermore, the binding of SIM2 marks a particular sub-category of enhancers known as super-enhancers. These regions are characterized by typical DNA modifications and Mediator co-occupancy (MED1 and MED12). Altogether, we provide evidence that SIM2 binds a specific set of enhancer elements thus explaining how SIM2 can regulate its gene network in neuronal features.

Biased allelic expression in human primary fibroblast single cells.

The study of gene expression in mammalian single cells via genomic technologies now provides the possibility to investigate the patterns of allelic gene expression. We used single-cell RNA sequencing to detect the allele-specific mRNA level in 203 single human primary fibroblasts over 133,633 unique heterozygous single-nucleotide variants (hetSNVs). We observed that at the snapshot of analyses, each cell contained mostly transcripts from one allele from the majority of genes; indeed, 76.4% of the hetSNVs displayed stochastic monoallelic expression in single cells. Remarkably, adjacent hetSNVs exhibited a haplotype-consistent allelic ratio; in contrast, distant sites located in two different genes were independent of the haplotype structure. Moreover, the allele-specific expression in single cells correlated with the abundance of the cellular transcript. We observed that genes expressing both alleles in the majority of the single cells at a given time point were rare and enriched with highly expressed genes. The relative abundance of each allele in a cell was controlled by some regulatory mechanisms given that we observed related single-cell allelic profiles according to genes. Overall, these results have direct implications in cellular phenotypic variability.

Tissue-specific effects of genetic and epigenetic variation on gene regulation and splicing.

Understanding how genetic variation affects distinct cellular phenotypes, such as gene expression levels, alternative splicing and DNA methylation levels, is essential for better understanding of complex diseases and traits. Furthermore, how inter-individual variation of DNA methylation is associated to gene expression is just starting to be studied. In this study, we use the GenCord cohort of 204 newborn Europeans' lymphoblastoid cell lines, T-cells and fibroblasts derived from umbilical cords. The samples were previously genotyped for 2.5 million SNPs, mRNA-sequenced, and assayed for methylation levels in 482,421 CpG sites. We observe that methylation sites associated to expression levels are enriched in enhancers, gene bodies and CpG island shores. We show that while the correlation between DNA methylation and gene expression can be positive or negative, it is very consistent across cell-types. However, this epigenetic association to gene expression appears more tissue-specific than the genetic effects on gene expression or DNA methylation (observed in both sharing estimations based on P-values and effect size correlations between cell-types). This predominance of genetic effects can also be reflected by the observation that allele specific expression differences between individuals dominate over tissue-specific effects. Additionally, we discover genetic effects on alternative splicing and interestingly, a large amount of DNA methylation correlating to alternative splicing, both in a tissue-specific manner. The locations of the SNPs and methylation sites involved in these associations highlight the participation of promoter proximal and distant regulatory regions on alternative splicing. Overall, our results provide high-resolution analyses showing how genome sequence variation has a broad effect on cellular phenotypes across cell-types, whereas epigenetic factors provide a secondary layer of variation that is more tissue-specific. Furthermore, the details of how this tissue-specificity may vary across inter-relations of molecular traits, and where these are occurring, can yield further insights into gene regulation and cellular biology as a whole.

Extrachromosomal driver mutations in glioblastoma and low-grade glioma.

Alteration of the number of copies of double minutes (DMs) with oncogenic EGFR mutations in response to tyrosine kinase inhibitors is a novel adaptive mechanism of glioblastoma. Here we provide evidence that such mutations in DMs, called here amplification-linked extrachromosomal mutations (ALEMs), originate extrachromosomally and could therefore be completely eliminated from the cancer cells. By exome sequencing of seven glioblastoma patients we reveal ALEMs in EGFR, PDGFRA and other genes. These mutations together with DMs are lost by cancer cells in culture. We confirm the extrachromosomal origin of such mutations by showing that wild-type and mutated DMs may coexist in the same tumour. Analysis of 4,198 tumours suggests the presence of ALEMs across different tumour types with the highest prevalence in glioblastomas and low-grade gliomas. The extrachromosomal nature of ALEMs explains the observed drastic changes in the amounts of mutated oncogenes (like EGFR or PDGFRA) in glioblastoma in response to environmental changes.

Frequent cases of RAS-mutated Down syndrome acute lymphoblastic leukaemia lack JAK2 mutations.

Children with Down syndrome (DS) and acute lymphoblastic leukaemia (ALL) have poorer survival and more relapses than non-DS children with ALL, highlighting an urgent need for deeper mechanistic understanding of DS-ALL. Here, using full-exome or cancer genes-targeted sequencing of 42 ALL samples from 39 DS patients, we uncover driver mutations in RAS, (KRAS and NRAS) recurring to a similar extent (15/42) as JAK2 (12/42) mutations or P2RY8-CRLF2 fusions (14/42). RAS mutations are almost completely mutually exclusive with JAK2 mutations (P=0.016), driving a combined total of two-thirds of analysed cases. Clonal architecture analysis reveals that both RAS and JAK2 drove sub-clonal expansions primarily initiated by CRLF2 rearrangements, and/or mutations in chromatin remodellers and lymphocyte differentiation factors. Remarkably, in 2/3 relapsed cases, there is a switch from a primary JAK2- or PTPN11-mutated sub-clone to a RAS-mutated sub-clone in relapse. These results provide important new insights informing the patient stratification strategies for targeted therapeutic approaches for DS-ALL.

Domains of genome-wide gene expression dysregulation in Down's syndrome.

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.

Cardiomyogenesis is controlled by the miR-99a/let-7c cluster and epigenetic modifications.

Understanding the molecular basis of cardiomyocyte development is critical for understanding the pathogenesis of pre- and post-natal cardiac disease. MicroRNAs (miRNAs) are post-transcriptional modulators of gene expression that play an important role in many developmental processes. Here, we show that the miR-99a/let-7c cluster, mapping on human chromosome 21, is involved in the control of cardiomyogenesis by altering epigenetic factors. By perturbing miRNA expression in mouse embryonic stem cells, we find that let-7c promotes cardiomyogenesis by upregulating genes involved in mesoderm specification (T/Bra and Nodal) and cardiac differentiation (Mesp1, Nkx2.5 and Tbx5). The action of let-7c is restricted to the early phase of mesoderm formation at the expense of endoderm and its late activation redirects cells toward other mesodermal derivatives. The Polycomb complex group protein Ezh2 is a direct target of let-7c, which promotes cardiac differentiation by modifying the H3K27me3 marks from the promoters of crucial cardiac transcription factors (Nkx2.5, Mef2c, Tbx5). In contrast, miR-99a represses cardiac differentiation via the nucleosome-remodeling factor Smarca5, attenuating the Nodal/Smad2 signaling. We demonstrated that the identified targets are underexpressed in human Down syndrome fetal heart specimens. By perturbing the expression levels of these miRNAs in embryonic stem cells, we were able to demonstrate that these miRNAs control lineage- and stage-specific transcription factors, working in concert with chromatin modifiers to direct cardiomyogenesis.

Exome sequencing identifies putative drivers of progression of transient myeloproliferative disorder to AMKL in infants with Down syndrome.

Some neonates with Down syndrome (DS) are diagnosed with self-regressing transient myeloproliferative disorder (TMD), and 20% to 30% of those progress to acute megakaryoblastic leukemia (AMKL). We performed exome sequencing in 7 TMD/AMKL cases and copy-number analysis in these and 10 additional cases. All TMD/AMKL samples contained GATA1 mutations. No exome-sequenced TMD/AMKL sample had other recurrently mutated genes. However, 2 of 5 TMD cases, and all AMKL cases, showed mutations/deletions other than GATA1, in genes proven as transformation drivers in non-DS leukemia (EZH2, APC, FLT3, JAK1, PARK2-PACRG, EXT1, DLEC1, and SMC3). One patient at the TMD stage revealed 2 clonal expansions with different GATA1 mutations, of which 1 clone had an additional driver mutation. Interestingly, it was the other clone that gave rise to AMKL after accumulating mutations in 7 other genes. Data suggest that GATA1 mutations alone are sufficient for clonal expansions, and additional driver mutations at the TMD stage do not necessarily predict AMKL progression. Later in infancy, leukemic progression requires "third-hit driver" mutations/somatic copy-number alterations found in non-DS leukemias. Putative driver mutations affecting WNT (wingless-related integration site), JAK-STAT (Janus kinase/signal transducer and activator of transcription), or MAPK/PI3K (mitogen-activated kinase/phosphatidylinositol-3 kinase) pathways were found in all cases, aberrant activation of which converges on overexpression of MYC.

The complex SNP and CNV genetic architecture of the increased risk of congenital heart defects in Down syndrome.

Congenital heart defect (CHD) occurs in 40% of Down syndrome (DS) cases. While carrying three copies of chromosome 21 increases the risk for CHD, trisomy 21 itself is not sufficient to cause CHD. Thus, additional genetic variation and/or environmental factors could contribute to the CHD risk. Here we report genomic variations that in concert with trisomy 21, determine the risk for CHD in DS. This case-control GWAS includes 187 DS with CHD (AVSD = 69, ASD = 53, VSD = 65) as cases, and 151 DS without CHD as controls. Chromosome 21-specific association studies revealed rs2832616 and rs1943950 as CHD risk alleles (adjusted genotypic P-values <0.05). These signals were confirmed in a replication cohort of 92 DS-CHD cases and 80 DS-without CHD (nominal P-value 0.0022). Furthermore, CNV analyses using a customized chromosome 21 aCGH of 135K probes in 55 DS-AVSD and 53 DS-without CHD revealed three CNV regions associated with AVSD risk (FDR ≤ 0.05). Two of these regions that are located within the previously identified CHD region on chromosome 21 were further confirmed in a replication study of 49 DS-AVSD and 45 DS- without CHD (FDR ≤ 0.05). One of these CNVs maps near the RIPK4 gene, and the second includes the ZBTB21 (previously ZNF295) gene, highlighting the potential role of these genes in the pathogenesis of CHD in DS. We propose that the genetic architecture of the CHD risk of DS is complex and includes trisomy 21, and SNP and CNV variations in chromosome 21. In addition, a yet-unidentified genetic variation in the rest of the genome may contribute to this complex genetic architecture.

Passive and active DNA methylation and the interplay with genetic variation in gene regulation.

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.