Rutger W. W. Brouwer

Cohesin and CTCF differentially affect chromatin architecture and gene expression in human cells

Significance For the 2m DNA to fit into the tiny cell nucleus, it is wrapped around nucleosomes and folded into loops clustering together in domains. Genome function depends on this 3D-organization, especially on-going dynamic processes like transcription. Techniques studying the network of DNA contacts genome-wide have recently revealed this 3D architecture, but the protein factors behind this are not understood. We study two proteins that are known to help form DNA loops: cohesin and CTCC-binding factor (CTCF). Respective depletion and analysis of DNA contacts genome-wide show that CTCF is required to separate neighboring folding domains and keep cohesin in place, whereas cohesin is important for shaping the domains. Consistently, we observe different changes of gene expression. Recent studies of genome-wide chromatin interactions have revealed that the human genome is partitioned into many self-associating topological domains. The boundary sequences between domains are enriched for binding sites of CTCC-binding factor (CTCF) and the cohesin complex, implicating these two factors in the establishment or maintenance of topological domains. To determine the role of cohesin and CTCF in higher-order chromatin architecture in human cells, we depleted the cohesin complex or CTCF and examined the consequences of loss of these factors on higher-order chromatin organization, as well as the transcriptome. We observed a general loss of local chromatin interactions upon disruption of cohesin, but the topological domains remain intact. However, we found that depletion of CTCF not only reduced intradomain interactions but also increased interdomain interactions. Furthermore, distinct groups of genes become misregulated upon depletion of cohesin and CTCF. Taken together, these observations suggest that CTCF and cohesin contribute differentially to chromatin organization and gene regulation.

Multiplexed chromosome conformation capture sequencing for rapid genome-scale high-resolution detection of long-range chromatin interactions

Chromosome conformation capture (3C) technology is a powerful and increasingly popular tool for analyzing the spatial organization of genomes. Several 3C variants have been developed (e.g., 4C, 5C, ChIA-PET, Hi-C), allowing large-scale mapping of long-range genomic interactions. Here we describe multiplexed 3C sequencing (3C-seq), a 4C variant coupled to next-generation sequencing, allowing genome-scale detection of long-range interactions with candidate regions. Compared with several other available techniques, 3C-seq offers a superior resolution (typically single restriction fragment resolution; approximately 1–8 kb on average) and can be applied in a semi-high-throughput fashion. It allows the assessment of long-range interactions of up to 192 genes or regions of interest in parallel by multiplexing library sequencing. This renders multiplexed 3C-seq an inexpensive, quick (total hands-on time of 2 weeks) and efficient method that is ideal for the in-depth analysis of complex genetic loci. The preparation of multiplexed 3C-seq libraries can be performed by any investigator with basic skills in molecular biology techniques. Data analysis requires basic expertise in bioinformatics and in Linux and Python environments. The protocol describes all materials, critical steps and bioinformatics tools required for successful application of 3C-seq technology.

Clinical significance of immunohistochemistry for detection of BAP1 mutations in uveal melanoma.

Uveal melanoma is a lethal cancer with a strong propensity to metastasize. Limited therapeutic options are available once the disease has disseminated. A strong predictor for metastasis is the loss of chromosome 3. Inactivating mutations in BAP1 encoding the BRCA1-associated protein 1 and located on chromosome 3p21.1, have been described in uveal melanoma and other types of cancer. In this study, we determined the prevalence of somatic BAP1 mutations and examined whether these mutations correlate with the functional expression of BAP1 in uveal melanoma tissue and with other clinical, histopathological and chromosomal parameters. We screened a cohort of 74 uveal melanomas for BAP1 mutations, using different deep sequencing methods. The frequency of BAP1 mutations in our study group was 47%. The expression of BAP1 protein was studied using immunohistochemistry. BAP1 staining was absent in 43% of the cases. BAP1 mutation status was strongly associated with BAP1 protein expression (P<0.001), loss of chromosome 3 (P<0.001), and other aggressive prognostic factors. Patients with a BAP1 mutation and absent BAP1 expression had an almost eightfold higher chance of developing metastases compared with those without these changes (P=0.002). We found a strong correlation between the immunohistochemical and sequencing data and therefore propose that, immunohistochemical screening for BAP1 should become routine in the histopathological work-up of uveal melanoma. Furthermore, our analysis indicates that loss of BAP1 may be particularly involved in the progression of uveal melanoma to an aggressive, metastatic phenotype.

Human USP18 deficiency underlies type 1 interferonopathy leading to severe pseudo-TORCH syndrome

Meuwissen and collaborators define a novel genetic cause of pseudo-TORCH syndrome, which resembles the sequelae of congenital infection and represents a novel type I interferonopathy.

Targeted Chromatin Capture (T2C): a novel high resolution high throughput method to detect genomic interactions and regulatory elements

BackgroundSignificant efforts have recently been put into the investigation of the spatial organization and the chromatin-interaction networks of genomes. Chromosome conformation capture (3C) technology and its derivatives are important tools used in this effort. However, many of these have limitations, such as being limited to one viewpoint, expensive with moderate to low resolution, and/or requiring a large sequencing effort. Techniques like Hi-C provide a genome-wide analysis. However, it requires massive sequencing effort with considerable costs. Here we describe a new technique termed Targeted Chromatin Capture (T2C), to interrogate large selected regions of the genome. T2C provides an unbiased view of the spatial organization of selected loci at superior resolution (single restriction fragment resolution, from 2 to 6 kbp) at much lower costs than Hi-C due to the lower sequencing effort.ResultsWe applied T2C on well-known model regions, the mouse β-globin locus and the human H19/IGF2 locus. In both cases we identified all known chromatin interactions. Furthermore, we compared the human H19/IGF2 locus data obtained from different chromatin conformation capturing methods with T2C data. We observed the same compartmentalization of the locus, but at a much higher resolution (single restriction fragments vs. the common 40 kbp bins) and higher coverage. Moreover, we compared the β-globin locus in two different biological samples (mouse primary erythroid cells and mouse fetal brain), where it is either actively transcribed or not, to identify possible transcriptional dependent interactions. We identified the known interactions in the β-globin locus and the same topological domains in both mouse primary erythroid cells and in mouse fetal brain with the latter having fewer interactions probably due to the inactivity of the locus. Furthermore, we show that interactions due to the important chromatin proteins, Ldb1 and Ctcf, in both tissues can be analyzed easily to reveal their role on transcriptional interactions and genome folding.ConclusionsT2C is an efficient, easy, and affordable with high (restriction fragment) resolution tool to address both genome compartmentalization and chromatin-interaction networks for specific genomic regions at high resolution for both clinical and non-clinical research.

Genome-wide DNA methylation profiling of non-small cell lung carcinomas

BackgroundNon-small cell lung carcinoma (NSCLC) is a complex malignancy that owing to its heterogeneity and poor prognosis poses many challenges to diagnosis, prognosis and patient treatment. DNA methylation is an important mechanism of epigenetic regulation involved in normal development and cancer. It is a very stable and specific modification and therefore in principle a very suitable marker for epigenetic phenotyping of tumors. Here we present a genome-wide DNA methylation analysis of NSCLC samples and paired lung tissues, where we combine MethylCap and next generation sequencing (MethylCap-seq) to provide comprehensive DNA methylation maps of the tumor and paired lung samples. The MethylCap-seq data were validated by bisulfite sequencing and methyl-specific polymerase chain reaction of selected regions.ResultsAnalysis of the MethylCap-seq data revealed a strong positive correlation between replicate experiments and between paired tumor/lung samples. We identified 57 differentially methylated regions (DMRs) present in all NSCLC tumors analyzed by MethylCap-seq. While hypomethylated DMRs did not correlate to any particular functional category of genes, the hypermethylated DMRs were strongly associated with genes encoding transcriptional regulators. Furthermore, subtelomeric regions and satellite repeats were hypomethylated in the NSCLC samples. We also identified DMRs that were specific to two of the major subtypes of NSCLC, adenocarcinomas and squamous cell carcinomas.ConclusionsCollectively, we provide a resource containing genome-wide DNA methylation maps of NSCLC and their paired lung tissues, and comprehensive lists of known and novel DMRs and associated genes in NSCLC.

NARWHAL, a primary analysis pipeline for NGS data

UNLABELLED The NARWHAL software pipeline has been developed to automate the primary analysis of Illumina sequencing data. This pipeline combines a new and flexible de-multiplexing tool with open-source aligners and automated quality assessment. The entire pipeline can be run using only one simple sample-sheet for diverse sequencing applications. NARWHAL creates a sample-oriented data structure and outperforms existing tools in speed. AVAILABILITY https://trac.nbic.nl/narwhal/.

Targeted Next Generation Sequencing reveals previously unidentified TSC1 and TSC2 mutations.

BackgroundTuberous sclerosis complex (TSC) is an autosomal dominant disorder caused by mutations in TSC1 and TSC2. Conventional DNA diagnostic screens identify a TSC1 or TSC2 mutation in 75 - 90% of individuals categorised with definite TSC. The remaining individuals either have a mutation that is undetectable using conventional methods, or possibly a mutation in another as yet unidentified gene.MethodsHere we apply a targeted Next Generation Sequencing (NGS) approach to screen the complete TSC1 and TSC2 genomic loci in 7 individuals fulfilling the clinical diagnostic criteria for definite TSC in whom no TSC1 or TSC2 mutations were identified using conventional screening methods.ResultsWe identified and confirmed pathogenic mutations in 3 individuals. In the remaining individuals we identified variants of uncertain clinical significance. The identified variants included mosaic changes, changes located deep in intronic sequences and changes affecting promoter regions that would not have been identified using exon-only based analyses.ConclusionsTargeted NGS of the TSC1 and TSC2 loci is a suitable method to increase the yield of mutations identified in the TSC patient population.

A novel complex, RUNX1-MYEF2, represses hematopoietic genes in erythroid cells

ABSTRACT RUNX1 is known to be an essential transcription factor for generating hematopoietic stem cells (HSC), but much less is known about its role in the downstream process of hematopoietic differentiation. RUNX1 has been shown to be part of a large transcription factor complex, together with LDB1, GATA1, TAL1, and ETO2 (N. Meier et al., Development 133:4913–4923, 2006) in erythroid cells. We used a tagging strategy to show that RUNX1 interacts with two novel protein partners, LSD1 and MYEF2, in erythroid cells. MYEF2 is bound in undifferentiated cells and is lost upon differentiation, whereas LSD1 is bound in differentiated cells. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) and microarray expression analysis were used to show that RUNX1 binds approximately 9,000 target sites in erythroid cells and is primarily active in the undifferentiated state. Functional analysis shows that a subset of the target genes is suppressed by RUNX1 via the newly identified partner MYEF2. Knockdown of Myef2 expression in developing zebrafish results in a reduced number of HSC.

Sp1/Sp3 transcription factors regulate hallmarks of megakaryocyte maturation and platelet formation and function.

Sp1 and Sp3 belong to the specificity proteins (Sp)/Krüppel-like transcription factor family. They are closely related, ubiquitously expressed, and recognize G-rich DNA motifs. They are thought to regulate generic processes such as cell-cycle and growth control, metabolic pathways, and apoptosis. Ablation of Sp1 or Sp3 in mice is lethal, and combined haploinsufficiency results in hematopoietic defects during the fetal stages. Here, we show that in adult mice, conditional pan-hematopoietic (Mx1-Cre) ablation of either Sp1 or Sp3 has minimal impact on hematopoiesis, whereas the simultaneous loss of Sp1 and Sp3 results in severe macrothrombocytopenia. This occurs in a cell-autonomous manner as shown by megakaryocyte-specific (Pf4-Cre) double-knockout mice. We employed flow cytometry, cell culture, and electron microscopy and show that although megakaryocyte numbers are normal in bone marrow and spleen, they display a less compact demarcation membrane system and a striking inability to form proplatelets. Through megakaryocyte transcriptomics and platelet proteomics, we identified several cytoskeleton-related proteins and downstream effector kinases, including Mylk, that were downregulated upon Sp1/Sp3 depletion, providing an explanation for the observed defects in megakaryopoiesis. Supporting this notion, selective Mylk inhibition by ML7 affected proplatelet formation and stabilization and resulted in defective ITAM receptor-mediated platelet aggregation.