Re-analysis of whole genome sequencing ends a diagnostic odyssey: Case report of an RNU4-2 related neurodevelopmental disorder.

Despite increasing knowledge of disease-causing genes in human genetics, approximately half of the individuals affected by neurodevelopmental disorders remain genetically undiagnosed. Part of this missing heritability might be caused by genetic variants outside of protein-coding genes, which are not routinely diagnostically investigated. A recent preprint identified de novo variants in the non-coding spliceosomal snRNA gene RNU4-2 as a cause of a frequent novel syndromic neurodevelopmental disorder. Here we mined 164 whole genome sequencing (WGS) trios from individuals with neurodevelopmental or multiple congenital anomaly disorders that received diagnostic genomic investigations at our clinic. We identify a recurrent de novo RNU4-2 variant (NR_003137.2(RNU4-2):n.64_65insT) in a 5-year-old girl with severe global developmental delay, hypotonia, microcephaly, and seizures that likely explains her phenotype, given that extensive previous genetic investigations failed to identify an alternative cause. We present detailed phenotyping of the individual obtained during a 5-year follow-up. This includes photographs showing recognizable facial features for this novel disorder, which might allow prioritizing other currently unexplained affected individuals sharing similar facial features for targeted investigations of RNU4-2. This case illustrates the power of re-analysis to solve previously unexplained cases even when a diagnostic genome remains negative.

The role of confined placental mosaicism in fetal growth restriction: A retrospective cohort study.

OBJECTIVE: To evaluate which cytogenetic characteristics of confined placental mosaicism (CPM) detected in the first trimester chorionic villi and/or placentas in terms of chromosome aberration, cell lineage involved and trisomy origin will lead to fetal growth restriction and low birthweight. METHODS: Cohort study using routinely collected perinatal data and cytogenetic data of non-invasive prenatal testing, the first trimester chorionic villi sampling and postnatal placentas. RESULTS: 215 CPM cases were found. Fetal growth restriction (FGR) and low birthweight below the 10th percentile (BW < p10) were seen in 34.0% and 23.1%, respectively. Excluding cases of trisomy 16, 29.1% showed FGR and 17.9% had a BW < p10. The highest rate of FGR and BW < p10 was found in CPM type 3, but differences with type 1 and 2 were not significant. FGR and BW < p10 were significantly more often observed in cases with meiotic trisomies. CONCLUSION: There is an association between CPM and FGR and BW < p10. This association is not restricted to trisomy 16, neither to CPM type 3, nor to CPM involving a meiotic trisomy. Pregnancies with all CPM types and origins should be considered to be at increased risk of FGR and low BW < p10. A close prenatal fetal monitoring is indicated in all cases of CPM.

Early detection of active Human CytomegaloVirus (hCMV) infection in pregnant women using data generated for noninvasive fetal aneuploidy testing.

BACKGROUND: Prenatal hCMV infections can lead to severe embryopathy and neurological sequelae in neonates. Screening during pregnancy is not recommended by global societies, as there is no effective therapy. Recently, several groups showed that maternal-fetal hCMV transmission can be strongly reduced by administering anti-viral agents early in pregnancy. This calls for a screening method to identify at risk pregnancies at an appropriate gestational age, with the possibility for large-scale enrolment. Non-Invasive Prenatal Testing (NIPT) for fetal aneuploidy screening early in pregnancy is already implemented in many countries and performed on a large-scale basis. We investigated the use of whole genome cell-free DNA (cfDNA) sequencing data, generated for the purpose of NIPT, as (pre-)screening tool to identify women with active hCMV-infections, eligible for therapy. METHODS: Coded raw sequencing NIPT data from 204,818 pregnant women from three testing laboratories were analyzed for the presence of hCMV-cfDNA. Samples were stratified by cfDNA-hCMV load. For validation and interpretation, diagnostic hCMV-qPCR and serology testing were performed on a subset of cfDNA-hCMV-positive (n = 112) and -negative (n = 127) samples. FINDINGS: In 1930 samples (0.94%) hCMV fragments were detected. Validation by hCMV-qPCR showed that samples with high cfDNA-hCMV load tested positive and cfDNA-hCMV-negative samples tested negative. In 32/112 cfDNA-hCMV-positive samples (28.6%) the serological profile suggested a recent primary infection: this was more likely in samples with high cfDNA-hCMV load (78.6%) than in samples with low cfDNA-hCMV load (11.0%). In none of the cfDNA-hCMV-negative samples serology was indicative of a recent primary infection. INTERPRETATION: Our study shows that large-scale (pre-)screening for both genetic fetal aberrations and active maternal hCMV infections during pregnancy can be combined in one cfDNA sequencing test, performed on a single blood sample, drawn in the first trimester of pregnancy. FUNDING: This work was partly funded by the Prenatal Screening Foundation Nijmegen, the Netherlands.

Evaluation of Circulating Tumor DNA as a Liquid Biomarker in Uveal Melanoma.

Purpose: Uveal melanoma (UM) has a high propensity to metastasize. Prognosis is associated with specific driver mutations and copy number variations (CNVs), but limited primary tumor tissue is available for molecular characterization due to eye-sparing irradiation treatment. This study aimed to assess the rise in circulating tumor DNA (ctDNA) levels in UM and evaluate its efficacy for CNV-profiling of patients with UM. Methods: In a pilot study, we assessed ctDNA levels in the blood of patients with UM (n = 18) at various time points, including the time of diagnosis (n = 13), during fractionated stereotactic radiotherapy (fSRT) treatment (n = 6), and upon detection of metastatic disease (n = 13). Shallow whole-genome sequencing (sWGS) combined with in silico size-selection was used to identify prognostically relevant CNVs in patients with UM (n = 26) from peripheral blood retrieved at the time of diagnosis (n = 9), during fSRT (n = 5), during post-treatment follow-up (n = 4), metastasis detection (n = 6), and metastasis follow-up (n = 4). Results: A total of 34 patients had blood analyzed for ctDNA detection (n = 18) and/or CNV analysis (n = 26) at various time points. At the time of diagnosis, 5 of 13 patients (38%) had detectable ctDNA (median = 0 copies/mL). Upon detection of metastatic disease, ctDNA was detected in 10 of 13 patients (77%) and showed increased ctDNA levels (median = 24 copies/mL, P < 0.01). Among the six patients analyzed during fSRT, three (50%) patients had detectable ctDNA at baseline and three of six (50%) patients had undetectable levels of ctDNA. During the fSRT regimen, ctDNA levels remained unchanged (P > 0.05). The ctDNA fractions were undetectable to low in localized disease, and sWGS did not elucidate chromosome 3 status from blood samples. However, in 7 of 10 (70%) patients with metastases, the detection of chromosome 3 loss corresponded to the high metastatic-risk class. Conclusions: The rise in ctDNA levels observed in patients with UM harboring metastases suggests its potential utility for CNV profiling. These findings highlight the potential of using ctDNA for metastasis detection and patient inclusion in therapeutic studies targeting metastatic UM.

Pyruvate metabolism controls chromatin remodeling during CD4+ T cell activation.

Upon antigen-specific T cell receptor (TCR) engagement, human CD4+ T cells proliferate and differentiate, a process associated with rapid transcriptional changes and metabolic reprogramming. Here, we show that the generation of extramitochondrial pyruvate is an important step for acetyl-CoA production and subsequent H3K27ac-mediated remodeling of histone acetylation. Histone modification, transcriptomic, and carbon tracing analyses of pyruvate dehydrogenase (PDH)-deficient T cells show PDH-dependent acetyl-CoA generation as a rate-limiting step during T activation. Furthermore, T cell activation results in the nuclear translocation of PDH and its association with both the p300 acetyltransferase and histone H3K27ac. These data support the tight integration of metabolic and histone-modifying enzymes, allowing metabolic reprogramming to fuel CD4+ T cell activation. Targeting this pathway may provide a therapeutic approach to specifically regulate antigen-driven T cell activation.

AMFR dysfunction causes autosomal recessive spastic paraplegia in human that is amenable to statin treatment in a preclinical model.

Hereditary spastic paraplegias (HSP) are rare, inherited neurodegenerative or neurodevelopmental disorders that mainly present with lower limb spasticity and muscle weakness due to motor neuron dysfunction. Whole genome sequencing identified bi-allelic truncating variants in AMFR, encoding a RING-H2 finger E3 ubiquitin ligase anchored at the membrane of the endoplasmic reticulum (ER), in two previously genetically unexplained HSP-affected siblings. Subsequently, international collaboration recognized additional HSP-affected individuals with similar bi-allelic truncating AMFR variants, resulting in a cohort of 20 individuals from 8 unrelated, consanguineous families. Variants segregated with a phenotype of mainly pure but also complex HSP consisting of global developmental delay, mild intellectual disability, motor dysfunction, and progressive spasticity. Patient-derived fibroblasts, neural stem cells (NSCs), and in vivo zebrafish modeling were used to investigate pathomechanisms, including initial preclinical therapy assessment. The absence of AMFR disturbs lipid homeostasis, causing lipid droplet accumulation in NSCs and patient-derived fibroblasts which is rescued upon AMFR re-expression. Electron microscopy indicates ER morphology alterations in the absence of AMFR. Similar findings are seen in amfra-/- zebrafish larvae, in addition to altered touch-evoked escape response and defects in motor neuron branching, phenocopying the HSP observed in patients. Interestingly, administration of FDA-approved statins improves touch-evoked escape response and motor neuron branching defects in amfra-/- zebrafish larvae, suggesting potential therapeutic implications. Our genetic and functional studies identify bi-allelic truncating variants in AMFR as a cause of a novel autosomal recessive HSP by altering lipid metabolism, which may potentially be therapeutically modulated using precision medicine with statins.

Web-accessible application for identifying pathogenic transcripts with RNA-seq: Increased sensitivity in diagnosis of neurodevelopmental disorders.

For neurodevelopmental disorders (NDDs), a molecular diagnosis is key for management, predicting outcome, and counseling. Often, routine DNA-based tests fail to establish a genetic diagnosis in NDDs. Transcriptome analysis (RNA sequencing [RNA-seq]) promises to improve the diagnostic yield but has not been applied to NDDs in routine diagnostics. Here, we explored the diagnostic potential of RNA-seq in 96 individuals including 67 undiagnosed subjects with NDDs. We performed RNA-seq on single individuals' cultured skin fibroblasts, with and without cycloheximide treatment, and used modified OUTRIDER Z scores to detect gene expression outliers and mis-splicing by exonic and intronic outliers. Analysis was performed by a user-friendly web application, and candidate pathogenic transcriptional events were confirmed by secondary assays. We identified intragenic deletions, monoallelic expression, and pseudoexonic insertions but also synonymous and non-synonymous variants with deleterious effects on transcription, increasing the diagnostic yield for NDDs by 13%. We found that cycloheximide treatment and exonic/intronic Z score analysis increased detection and resolution of aberrant splicing. Importantly, in one individual mis-splicing was found in a candidate gene nearly matching the individual's specific phenotype. However, pathogenic splicing occurred in another neuronal-expressed gene and provided a molecular diagnosis, stressing the need to customize RNA-seq. Lastly, our web browser application allowed custom analysis settings that facilitate diagnostic application and ranked pathogenic transcripts as top candidates. Our results demonstrate that RNA-seq is a complementary method in the genomic diagnosis of NDDs and, by providing accessible analysis with improved sensitivity, our transcriptome analysis approach facilitates wider implementation of RNA-seq in routine genome diagnostics.

Epigenetic Analyses of Human Left Atrial Tissue Identifies Gene Networks Underlying Atrial Fibrillation.

BACKGROUND: Atrial fibrillation (AF) often arises from structural abnormalities in the left atria (LA). Annotation of the noncoding genome in human LA is limited, as are effects on gene expression and chromatin architecture. Many AF-associated genetic variants reside in noncoding regions; this knowledge gap impairs efforts to understand the molecular mechanisms of AF and cardiac conduction phenotypes. METHODS: We generated a model of the LA noncoding genome by profiling 7 histone post-translational modifications (active: H3K4me3, H3K4me2, H3K4me1, H3K27ac, H3K36me3; repressive: H3K27me3, H3K9me3), CTCF binding, and gene expression in samples from 5 individuals without structural heart disease or AF. We used MACS2 to identify peak regions (P<0.01), applied a Markov model to classify regulatory elements, and annotated this model with matched gene expression data. We intersected chromatin states with expression quantitative trait locus, DNA methylation, and HiC chromatin interaction data from LA and left ventricle. Finally, we integrated genome-wide association data for AF and electrocardiographic traits to link disease-related variants to genes. RESULTS: Our model identified 21 epigenetic states, encompassing regulatory motifs, such as promoters, enhancers, and repressed regions. Genes were regulated by proximal chromatin states; repressive states were associated with a significant reduction in gene expression (P<2×10-16). Chromatin states were differentially methylated, promoters were less methylated than repressed regions (P<2×10-16). We identified over 15 000 LA-specific enhancers, defined by homeobox family motifs, and annotated several cardiovascular disease susceptibility loci. Intersecting AF and PR genome-wide association studies loci with long-range chromatin conformation data identified a gene interaction network dominated by NKX2-5, TBX3, ZFHX3, and SYNPO2L. CONCLUSIONS: Profiling the noncoding genome provides new insights into the gene expression and chromatin regulation in human LA tissue. These findings enabled identification of a gene network underlying AF; our experimental and analytic approach can be extended to identify molecular mechanisms for other cardiac diseases and traits.

Recently Evolved Enhancers Emerge with High Interindividual Variability and Less Frequently Associate with Disease.

Mutations in non-coding regulatory DNA such as enhancers underlie a wide variety of diseases including developmental disorders and cancer. As enhancers rapidly evolve, understanding their function and configuration in non-human disease models can have important clinical applications. Here, we analyze enhancer configurations in tissues isolated from the common marmoset, a widely used primate model for human disease. Integrating these data with human and mouse data, we find that enhancers containing trait-associated variants are preferentially conserved. In contrast, most human-specific enhancers are highly variable between individuals, with a subset failing to contact promoters. These are located further away from genes and more often reside in inactive B-compartments. Our data show that enhancers typically emerge as instable elements with minimal biological impact prior to their integration in a transcriptional program. Furthermore, our data provide insight into which trait variations in enhancers can be faithfully modeled using the common marmoset.

Hominin-specific regulatory elements selectively emerged in oligodendrocytes and are disrupted in autism patients.

Speciation is associated with substantial rewiring of the regulatory circuitry underlying the expression of genes. Determining which changes are relevant and underlie the emergence of the human brain or its unique susceptibility to neural disease has been challenging. Here we annotate changes to gene regulatory elements (GREs) at cell type resolution in the brains of multiple primate species spanning most of primate evolution. We identify a unique set of regulatory elements that emerged in hominins prior to the separation of humans and chimpanzees. We demonstrate that these hominin gains perferentially affect oligodendrocyte function postnatally and are preferentially affected in the brains of autism patients. This preference is also observed for human-specific GREs suggesting this system is under continued selective pressure. Our data provide a roadmap of regulatory rewiring across primate evolution providing insight into the genomic changes that underlie the emergence of the brain and its susceptibility to neural disease.

Multi-contact 4C: long-molecule sequencing of complex proximity ligation products to uncover local cooperative and competitive chromatin topologies.

We present the experimental protocol and data analysis toolbox for multi-contact 4C (MC-4C), a new proximity ligation method tailored to study the higher-order chromatin contact patterns of selected genomic sites. Conventional chromatin conformation capture (3C) methods fragment proximity ligation products for efficient analysis of pairwise DNA contacts. By contrast, MC-4C is designed to preserve and collect large concatemers of proximity ligated fragments for long-molecule sequencing on an Oxford Nanopore or Pacific Biosciences platform. Each concatemer of proximity ligation products represents a snapshot topology of a different individual allele, revealing its multi-way chromatin interactions. By inverse PCR with primers specific for a fragment of interest (the viewpoint) and DNA size selection, sequencing is selectively targeted to thousands of different complex interactions containing this viewpoint. A tailored statistical analysis toolbox is able to generate background models and three-way interaction profiles from the same dataset. These profiles can be used to distinguish whether contacts between more than two regulatory sequences are mutually exclusive or, conversely, simultaneously occurring at chromatin hubs. The entire procedure can be completed in 2 w, and requires standard molecular biology and data analysis skills and equipment, plus access to a third-generation sequencing platform.

Bayesian mixture regression analysis for regulation of Pluripotency in ES cells.

BACKGROUND: Observed levels of gene expression strongly depend on both activity of DNA binding transcription factors (TFs) and chromatin state through different histone modifications (HMs). In order to recover the functional relationship between local chromatin state, TF binding and observed levels of gene expression, regression methods have proven to be useful tools. They have been successfully applied to predict mRNA levels from genome-wide experimental data and they provide insight into context-dependent gene regulatory mechanisms. However, heterogeneity arising from gene-set specific regulatory interactions is often overlooked. RESULTS: We show that regression models that predict gene expression by using experimentally derived ChIP-seq profiles of TFs can be significantly improved by mixture modelling. In order to find biologically relevant gene clusters, we employ a Bayesian allocation procedure which allows us to integrate additional biological information such as three-dimensional nuclear organization of chromosomes and gene function. The data integration procedure involves transforming the additional data into gene similarity values. We propose a generic similarity measure that is especially suitable for situations where the additional data are of both continuous and discrete type, and compare its performance with similar measures in the context of mixture modelling. CONCLUSIONS: We applied the proposed method on a data from mouse embryonic stem cells (ESC). We find that including additional data results in mixture components that exhibit biologically meaningful gene clusters, and provides valuable insight into the heterogeneity of the regulatory interactions.

4C-seq from beginning to end: A detailed protocol for sample preparation and data analysis.

Chromosome conformation capture (3C) methods measure DNA contact frequencies based on nuclear proximity ligation, to uncover in vivo genomic folding patterns. 4C-seq is a derivative 3C method, designed to search the genome for sequences contacting a selected genomic site of interest. 4C-seq employs inverse PCR and next generation sequencing to amplify, identify and quantify its proximity ligated DNA fragments. It generates high-resolution contact profiles for selected genomic sites based on limited amounts of sequencing reads. 4C-seq can be used to study multiple aspects of genome organization. It primarily serves to identify specific long-range DNA contacts between individual regulatory DNA modules, forming for example regulatory chromatin loops between enhancers and promoters, or architectural chromatin loops between cohesin- and CTCF- associated domain boundaries. Additionally, 4C-seq contact profiles can reveal the contours of contact domains and can identify the structural domains that co-occupy the same nuclear compartment. Here, we present an improved step-by-step protocol for sample preparation and the generation of 4C-seq sequencing libraries, including an optimized PCR and 4C template purification strategy. In addition, a data processing pipeline is provided which processes multiplexed 4C-seq reads directly from FASTQ files and generates files compatible with standard genome browsers for visualization and further statistical analysis of the data such as peak calling using peakC. The protocols and the pipeline presented should readily allow anyone to generate, visualize and interpret their own high resolution 4C contact datasets.

Chromatin Conformation Links Putative Enhancers in Intracranial Aneurysm-Associated Regions to Potential Candidate Genes.

Background We previously showed that intracranial aneurysm ( IA )-associated single-nucleotide polymorphisms are enriched in promoters and putative enhancers identified in the human circle of Willis, on which IA s develop, suggesting a role for promoters and enhancers in IAs . We further investigated the role of putative enhancers in the pathogenesis of IA by identifying their potential target genes and validating their regulatory activity. Methods and Results Using our previously published circle of Willis chromatin immunoprecipitation and sequencing data, we selected 34 putative enhancers in IA -associated regions from genome-wide association studies. We then used a chromatin conformation capture technique to prioritize target genes and found that 15 putative enhancers interact with the promoters of 6 target genes: SOX 17, CDKN 2B, MTAP , CNNM 2, RPEL 1, and GATA 6. Subsequently, we assessed the activity of these putative enhancers in vivo in zebrafish embryos and confirmed activity for 8 putative enhancers. Last, we found that all 6 target genes are expressed in the circle of Willis, on the basis of RNA sequencing data and in situ hybridization. Furthermore, in situ hybridization showed that these genes are expressed in multiple cell types in the circle of Willis. Conclusions In 4 of 6 IA -associated genome-wide association study regions, we identified 8 putative enhancers that are active in vivo and interact with 6 nearby genes, suggesting that these genes are regulated by the identified putative enhancers. These genes, SOX 17, CDKN 2B, MTAP , CNNM 2, RPEL 1, and GATA 6, are therefore potential candidate genes involved in IA pathogenesis and should be studied using animal models in the future.

Enhancer hubs and loop collisions identified from single-allele topologies.

Chromatin folding contributes to the regulation of genomic processes such as gene activity. Existing conformation capture methods characterize genome topology through analysis of pairwise chromatin contacts in populations of cells but cannot discern whether individual interactions occur simultaneously or competitively. Here we present multi-contact 4C (MC-4C), which applies Nanopore sequencing to study multi-way DNA conformations of individual alleles. MC-4C distinguishes cooperative from random and competing interactions and identifies previously missed structures in subpopulations of cells. We show that individual elements of the ß-globin superenhancer can aggregate into an enhancer hub that can simultaneously accommodate two genes. Neighboring chromatin domain loops can form rosette-like structures through collision of their CTCF-bound anchors, as seen most prominently in cells lacking the cohesin-unloading factor WAPL. Here, massive collision of CTCF-anchored chromatin loops is believed to reflect 'cohesin traffic jams'. Single-allele topology studies thus help us understand the mechanisms underlying genome folding and functioning.

peakC: a flexible, non-parametric peak calling package for 4C and Capture-C data.

It is becoming increasingly clear that chromosome organization plays an important role in gene regulation. High-resolution methods such as 4C, Capture-C and promoter capture Hi-C (PCHiC) enable the study of chromatin loops such as those formed between promoters and enhancers or CTCF/cohesin binding sites. An important aspect of 4C/Capture-C/PCHiC analyses is the reliable identification of chromatin loops, preferably not based on visual inspection of a DNA contact profile, but on reproducible statistical analysis that robustly scores interaction peaks in the non-uniform contact background. Here, we present peakC, an R package for the analysis of 4C/Capture-C/PCHiC data. We generated 4C data for 13 viewpoints in two tissues in at least triplicate to test our methods. We developed a non-parametric peak caller based on rank-products. Sampling analysis shows that not read depth but template quality is the most important determinant of success in 4C experiments. By performing peak calling on single experiments we show that the peak calling results are similar to the replicate experiments, but that false positive rates are significantly reduced by performing replicates. Our software is user-friendly and enables robust peak calling for one-vs-all chromosome capture experiments. peakC is available at: https://github.com/deWitLab/peakC.

Sensitive Monogenic Noninvasive Prenatal Diagnosis by Targeted Haplotyping.

During pregnancy, cell-free DNA (cfDNA) in maternal blood encompasses a small percentage of cell-free fetal DNA (cffDNA), an easily accessible source for determination of fetal disease status in risk families through non-invasive procedures. In case of monogenic heritable disease, background maternal cfDNA prohibits direct observation of the maternally inherited allele. Non-invasive prenatal diagnostics (NIPD) of monogenic diseases therefore relies on parental haplotyping and statistical assessment of inherited alleles from cffDNA, techniques currently unavailable for routine clinical practice. Here, we present monogenic NIPD (MG-NIPD), which requires a blood sample from both parents, for targeted locus amplification (TLA)-based phasing of heterozygous variants selectively at a gene of interest. Capture probes-based targeted sequencing of cfDNA from the pregnant mother and a tailored statistical analysis enables predicting fetal gene inheritance. MG-NIPD was validated for 18 pregnancies, focusing on CFTR, CYP21A2, and HBB. In all cases we could predict the inherited alleles with >98% confidence, even at relatively early stages (8 weeks) of pregnancy. This prediction and the accuracy of parental haplotyping was confirmed by sequencing of fetal material obtained by parallel invasive procedures. MG-NIPD is a robust method that requires standard instrumentation and can be implemented in any clinic to provide families carrying a severe monogenic disease with a prenatal diagnostic test based on a simple blood draw.

Enhancers reside in a unique epigenetic environment during early zebrafish development.

BACKGROUND: Enhancers, not promoters, are the most dynamic in their DNA methylation status throughout development and differentiation. Generally speaking, enhancers that are primed to or actually drive gene expression are characterized by relatively low levels of DNA methylation (hypo-methylation), while inactive enhancers display hyper-methylation of the underlying DNA. The direct functional significance of the DNA methylation state of enhancers is, however, unclear for most loci. RESULTS: In contrast to conventional epigenetic interactions at enhancers, we find that DNA methylation status and enhancer activity during early zebrafish development display very unusual correlation characteristics: hypo-methylation is a unique feature of primed enhancers whereas active enhancers are generally hyper-methylated. The hypo-methylated enhancers that we identify (hypo-enhancers) are enriched close to important transcription factors that act later in development. Interestingly, hypo-enhancers are de-methylated shortly before the midblastula transition and reside in a unique epigenetic environment. Finally, we demonstrate that hypo-enhancers do become active at later developmental stages and that they are physically associated with the transcriptional start site of target genes, irrespective of target gene activity. CONCLUSIONS: We demonstrate that early development in zebrafish embodies a time window characterized by non-canonical DNA methylation-enhancer relationships, including global DNA hypo-methylation of inactive enhancers and DNA hyper-methylation of active enhancers.

Cause and Consequence of Tethering a SubTAD to Different Nuclear Compartments.

Detailed genomic contact maps have revealed that chromosomes are structurally organized in megabase-sized topologically associated domains (TADs) that encompass smaller subTADs. These domains segregate in the nuclear space to form active and inactive nuclear compartments, but cause and consequence of compartmentalization are largely unknown. Here, we combined lacO/lacR binding platforms with allele-specific 4C technologies to track their precise position in the three-dimensional genome upon recruitment of NANOG, SUV39H1, or EZH2. We observed locked genomic loci resistant to spatial repositioning and unlocked loci that could be repositioned to different nuclear subcompartments with distinct chromatin signatures. Focal protein recruitment caused the entire subTAD, but not surrounding regions, to engage in new genomic contacts. Compartment switching was found uncoupled from transcription changes, and the enzymatic modification of histones per se was insufficient for repositioning. Collectively, this suggests that trans-associated factors influence three-dimensional compartmentalization independent of their cis effect on local chromatin composition and activity.

Epigenomic annotation of gene regulatory alterations during evolution of the primate brain.

Although genome sequencing has identified numerous noncoding alterations between primate species, which of those are regulatory and potentially relevant to the evolution of the human brain is unclear. Here we annotated cis-regulatory elements (CREs) in the human, rhesus macaque and chimpanzee genomes using chromatin immunoprecipitation followed by sequencing (ChIP-seq) in different anatomical regions of the adult brain. We found high similarity in the genomic positioning of rhesus macaque and human CREs, suggesting that the majority of these elements were already present in a common ancestor 25 million years ago. Most of the observed regulatory changes between humans and rhesus macaques occurred before the ancestral separation of humans and chimpanzees, leaving a modest set of regulatory elements with predicted human specificity. Our data refine previous predictions and hypotheses on the consequences of genomic changes between primate species and allow the identification of regulatory alterations relevant to the evolution of the brain.