Targeted chromatin conformation analysis identifies novel distal neural enhancers of ZEB2 in pluripotent stem cell differentiation

The transcription factor zinc finger E-box binding protein 2 (ZEB2) controls embryonic and adult cell fate decisions and cellular maturation in many stem/progenitor cell types. Defects in these processes in specific cell types underlie several aspects of Mowat–Wilson syndrome (MOWS), which is caused by ZEB2 haplo-insufficiency. Human ZEB2, like mouse Zeb2, is located on chromosome 2 downstream of a ±3.5 Mb-long gene-desert, lacking any protein-coding gene. Using temporal targeted chromatin capture (T2C), we show major chromatin structural changes based on mapping in-cis proximities between the ZEB2 promoter and this gene desert during neural differentiation of human-induced pluripotent stem cells, including at early neuroprogenitor cell (NPC)/rosette state, where ZEB2 mRNA levels increase significantly. Combining T2C with histone-3 acetylation mapping, we identified three novel candidate enhancers about 500 kb upstream of the ZEB2 transcription start site. Functional luciferase-based assays in heterologous cells and NPCs reveal co-operation between these three enhancers. This study is the first to document in-cis Regulatory Elements located in ZEB2’s gene desert. The results further show the usability of T2C for future studies of ZEB2 REs in differentiation and maturation of multiple cell types and the molecular characterization of newly identified MOWS patients that lack mutations in ZEB2 protein-coding exons.

Mild dyserythropoiesis and β-like globin gene expression imbalance due to the loss of histone chaperone ASF1B

The expression of the human β-like globin genes follows a well-orchestrated developmental pattern, undergoing two essential switches, the first one during the first weeks of gestation (ε to γ), and the second one during the perinatal period (γ to β). The γ- to β-globin gene switching mechanism includes suppression of fetal (γ-globin, HbF) and activation of adult (β-globin, HbA) globin gene transcription. In hereditary persistence of fetal hemoglobin (HPFH), the γ-globin suppression mechanism is impaired leaving these individuals with unusual elevated levels of fetal hemoglobin (HbF) in adulthood. Recently, the transcription factors KLF1 and BCL11A have been established as master regulators of the γ- to β-globin switch. Previously, a genomic variant in the KLF1 gene, identified by linkage analysis performed on twenty-seven members of a Maltese family, was found to be associated with HPFH. However, variation in the levels of HbF among family members, and those from other reported families carrying genetic variants in KLF1, suggests additional contributors to globin switching. ASF1B was downregulated in the family members with HPFH. Here, we investigate the role of ASF1B in γ- to β-globin switching and erythropoiesis in vivo. Mouse-human interspecies ASF1B protein identity is 91.6%. By means of knockdown functional assays in human primary erythroid cultures and analysis of the erythroid lineage in Asf1b knockout mice, we provide evidence that ASF1B is a novel contributor to steady-state erythroid differentiation, and while its loss affects the balance of globin expression, it has no major role in hemoglobin switching.

Multifaceted actions of Zeb2 in postnatal neurogenesis from the ventricular-subventricular zone to the olfactory bulb

The transcription factor Zeb2 controls fate specification and subsequent differentiation and maturation of multiple cell types in various embryonic tissues. It binds many protein partners, including activated Smad proteins and the NuRD co-repressor complex. How Zeb2 subdomains support cell differentiation in various contexts has remained elusive. Here, we studied the role of Zeb2 and its domains in neurogenesis and neural differentiation in the young postnatal ventricular-subventricular zone (V-SVZ), in which neural stem cells generate olfactory bulb-destined interneurons. Conditional Zeb2 knockouts and separate acute loss- and gain-of-function approaches indicated that Zeb2 is essential for controlling apoptosis and neuronal differentiation of V-SVZ progenitors before and after birth, and we identified Sox6 as a potential downstream target gene of Zeb2. Zeb2 genetic inactivation impaired the differentiation potential of the V-SVZ niche in a cell-autonomous fashion. We also provide evidence that its normal function in the V-SVZ also involves non-autonomous mechanisms. Additionally, we demonstrate distinct roles for Zeb2 protein-binding domains, suggesting that Zeb2 partners co-determine neuronal output from the mouse V-SVZ in both quantitative and qualitative ways in early postnatal life.

Integrative and perturbation-based analysis of the transcriptional dynamics of TGFβ/BMP system components in transition from embryonic stem cells to neural progenitors

Cooperative actions of extrinsic signals and cell‐intrinsic transcription factors alter gene regulatory networks enabling cells to respond appropriately to environmental cues. Signaling by transforming growth factor type β (TGFβ) family ligands (eg, bone morphogenetic proteins [BMPs] and Activin/Nodal) exerts cell‐type specific and context‐dependent transcriptional changes, thereby steering cellular transitions throughout embryogenesis. Little is known about coordinated regulation and transcriptional interplay of the TGFβ system. To understand intrafamily transcriptional regulation as part of this system's actions during development, we selected 95 of its components and investigated their mRNA‐expression dynamics, gene‐gene interactions, and single‐cell expression heterogeneity in mouse embryonic stem cells transiting to neural progenitors. Interrogation at 24 hour intervals identified four types of temporal gene transcription profiles that capture all stages, that is, pluripotency, epiblast formation, and neural commitment. Then, between each stage we performed esiRNA‐based perturbation of each individual component and documented the effect on steady‐state mRNA levels of the remaining 94 components. This exposed an intricate system of multilevel regulation whereby the majority of gene‐gene interactions display a marked cell‐stage specific behavior. Furthermore, single‐cell RNA‐profiling at individual stages demonstrated the presence of detailed co‐expression modules and subpopulations showing stable co‐expression modules such as that of the core pluripotency genes at all stages. Our combinatorial experimental approach demonstrates how intrinsically complex transcriptional regulation within a given pathway is during cell fate/state transitions.

The dynamic emergence of GATA1 complexes identified in in vitro embryonic stem cell differentiation and in vivo mouse fetal liver

GATA1 is an essential transcriptional regulator of myeloid hematopoietic differentiation towards red blood cells. During erythroid differentiation, GATA1 forms different complexes with other transcription factors such as LDB1, TAL1, E2A and LMO2 ("the LDB1 complex") or with FOG1. The functions of GATA1 complexes have been studied extensively in definitive erythroid differentiation; however, the temporal and spatial formation of these complexes during erythroid development is unknown. We applied proximity ligation assay (PLA) to detect, localize and quantify individual interactions during embryonic stem cell differentiation and in mouse fetal liver (FL) tissue. We show that GATA1/LDB1 interactions appear before the proerythroblast stage and increase in a subset of the CD71⁺/TER119^- cells to activate the terminal erythroid differentiation program in 12.5 day FL. Using Ldb1 and Gata1 knockdown FL cells, we studied the functional contribution of the GATA1/LDB1 complex during differentiation. This shows that the active LDB1 complex appears quite late at the proerythroblast stage of differentiation and confirms the power of PLA in studying the dynamic interaction of proteins in cell differentiation at the single cell level. We provide dynamic insight into the temporal and spatial formation of the GATA1 and LDB1 transcription factor complexes during hematopoietic development and differentiation.

Forces driving the three-dimensional folding of eukaryotic genomes

The last decade has radically renewed our understanding of higher order chromatin folding in the eukaryotic nucleus. As a result, most current models are in support of a mostly hierarchical and relatively stable folding of chromosomes dividing chromosomal territories into A‐ (active) and B‐ (inactive) compartments, which are then further partitioned into topologically associating domains (TADs), each of which is made up from multiple loops stabilized mainly by the CTCF and cohesin chromatin‐binding complexes. Nonetheless, the structure‐to‐function relationship of eukaryotic genomes is still not well understood. Here, we focus on recent work highlighting the biophysical and regulatory forces that contribute to the spatial organization of genomes, and we propose that the various conformations that chromatin assumes are not so much the result of a linear hierarchy, but rather of both converging and conflicting dynamic forces that act on it.

The detailed 3D multi-loop aggregate/rosette chromatin architecture and functional dynamic organization of the human and mouse genomes

BACKGROUND

The dynamic three-dimensional chromatin architecture of genomes and its co-evolutionary connection to its function-the storage, expression, and replication of genetic information-is still one of the central issues in biology. Here, we describe the much debated 3D architecture of the human and mouse genomes from the nucleosomal to the megabase pair level by a novel approach combining selective high-throughput high-resolution chromosomal interaction capture (T2C), polymer simulations, and scaling analysis of the 3D architecture and the DNA sequence.

RESULTS

The genome is compacted into a chromatin quasi-fibre with ~5 ± 1 nucleosomes/11 nm, folded into stable ~30-100 kbp loops forming stable loop aggregates/rosettes connected by similar sized linkers. Minor but significant variations in the architecture are seen between cell types and functional states. The architecture and the DNA sequence show very similar fine-structured multi-scaling behaviour confirming their co-evolution and the above.

CONCLUSIONS

This architecture, its dynamics, and accessibility, balance stability and flexibility ensuring genome integrity and variation enabling gene expression/regulation by self-organization of (in)active units already in proximity. Our results agree with the heuristics of the field and allow "architectural sequencing" at a genome mechanics level to understand the inseparable systems genomic properties.

Exploiting native forces to capture chromosome conformation in mammalian cell nuclei

Mammalian interphase chromosomes fold into a multitude of loops to fit the confines of cell nuclei, and looping is tightly linked to regulated function. Chromosome conformation capture (3C) technology has significantly advanced our understanding of this structure‐to‐function relationship. However, all 3C‐based methods rely on chemical cross‐linking to stabilize spatial interactions. This step remains a “black box” as regards the biases it may introduce, and some discrepancies between microscopy and 3C studies have now been reported. To address these concerns, we developed “i3C”, a novel approach for capturing spatial interactions without a need for cross‐linking. We apply i3C to intact nuclei of living cells and exploit native forces that stabilize chromatin folding. Using different cell types and loci, computational modeling, and a methylation‐based orthogonal validation method, “TALE‐iD”, we show that native interactions resemble cross‐linked ones, but display improved signal‐to‐noise ratios and are more focal on regulatory elements and CTCF sites, while strictly abiding to topologically associating domain restrictions.

Zeb2 Regulates Cell Fate at the Exit from Epiblast State in Mouse Embryonic Stem Cells

In human embryonic stem cells (ESCs) the transcription factor Zeb2 regulates neuroectoderm versus mesendoderm formation, but it is unclear how Zeb2 affects the global transcriptional regulatory network in these cell‐fate decisions. We generated Zeb2 knockout (KO) mouse ESCs, subjected them as embryoid bodies (EBs) to neural and general differentiation and carried out temporal RNA‐sequencing (RNA‐seq) and reduced representation bisulfite sequencing (RRBS) analysis in neural differentiation. This shows that Zeb2 acts preferentially as a transcriptional repressor associated with developmental progression and that Zeb2 KO ESCs can exit from their naïve state. However, most cells in these EBs stall in an early epiblast‐like state and are impaired in both neural and mesendodermal differentiation. Genes involved in pluripotency, epithelial‐to‐mesenchymal transition (EMT), and DNA‐(de)methylation, including Tet1, are deregulated in the absence of Zeb2. The observed elevated Tet1 levels in the mutant cells and the knowledge of previously mapped Tet1‐binding sites correlate with loss‐of‐methylation in neural‐stimulating conditions, however, after the cells initially acquired the correct DNA‐methyl marks. Interestingly, cells from such Zeb2 KO EBs maintain the ability to re‐adapt to 2i + LIF conditions even after prolonged differentiation, while knockdown of Tet1 partially rescues their impaired differentiation. Hence, in addition to its role in EMT, Zeb2 is critical in ESCs for exit from the epiblast state, and links the pluripotency network and DNA‐methylation with irreversible commitment to differentiation. Stem Cells2017;35:611–625

Binding of nuclear factor κB to noncanonical consensus sites reveals its multimodal role during the early inflammatory response

Mammalian cells have developed intricate mechanisms to interpret, integrate, and respond to extracellular stimuli. For example, tumor necrosis factor (TNF) rapidly activates proinflammatory genes, but our understanding of how this occurs against the ongoing transcriptional program of the cell is far from complete. Here, we monitor the early phase of this cascade at high spatiotemporal resolution in TNF-stimulated human endothelial cells. NF-κB, the transcription factor complex driving the response, interferes with the regulatory machinery by binding active enhancers already in interaction with gene promoters. Notably, >50% of these enhancers do not encode canonical NF-κB binding motifs. Using a combination of genomics tools, we find that binding site selection plays a key role in NF-κΒ–mediated transcriptional activation and repression. We demonstrate the latter by describing the synergy between NF-κΒ and the corepressor JDP2. Finally, detailed analysis of a 2.8-Mbp locus using sub-kbp-resolution targeted chromatin conformation capture and genome editing uncovers how NF-κΒ that has just entered the nucleus exploits pre-existing chromatin looping to exert its multimodal role. This work highlights the involvement of topology in cis-regulatory element function during acute transcriptional responses, where primary DNA sequence and its higher-order structure constitute a regulatory context leading to either gene activation or repression.

The core spliceosome as target and effector of non-canonical ATM signalling

In response to DNA damage, tissue homoeostasis is ensured by protein networks promoting DNA repair, cell cycle arrest or apoptosis. DNA damage response signalling pathways coordinate these processes, partly by propagating gene-expression-modulating signals. DNA damage influences not only the abundance of messenger RNAs, but also their coding information through alternative splicing. Here we show that transcription-blocking DNA lesions promote chromatin displacement of late-stage spliceosomes and initiate a positive feedback loop centred on the signalling kinase ATM. We propose that initial spliceosome displacement and subsequent R-loop formation is triggered by pausing of RNA polymerase at DNA lesions. In turn, R-loops activate ATM, which signals to impede spliceosome organization further and augment ultraviolet-irradiation-triggered alternative splicing at the genome-wide level. Our findings define R-loop-dependent ATM activation by transcription-blocking lesions as an important event in the DNA damage response of non-replicating cells, and highlight a key role for spliceosome displacement in this process.

Sox2 regulates the emergence of lung basal cells by directly activating the transcription of Trp63

Lung development is determined by the coordinated expression of several key genes. Previously, we and others have shown the importance of the sex determining region Y-box 2 (Sox2) gene in lung development. Transgenic expression of Sox2 during lung development resulted in cystic airways, and here we show that modulating the timing of ectopic Sox2 expression in the branching regions of the developing lung results in variable cystic lesions resembling the spectrum of the human congenital disorder congenital cystic adenomatoid malformation (CCAM). Sox2 dominantly differentiated naive epithelial cells into the proximal lineage irrespective of the presence of Fgf10. Sox2 directly induced the expression of Trp63, the master switch toward the basal cell lineage and induced the expression of Gata6, a factor involved in the emergence of bronchoalveolar stem cells. We showed that SOX2 and TRP63 are coexpressed in the lungs of human patients with type II CCAM. The combination of premature differentiation toward the proximal cell lineage and the induction of proliferation resulted in the cyst-like structures. Thus, we show that Sox2 is directly responsible for the emergence of two lung progenitor cells: basal cells by regulating the master gene Trp63 and bronchoalveolar stem cells by regulating Gata6.

NLS-tagging: an alternative strategy to tag nuclear proteins

The characterization of transcription factor complexes and their binding sites in the genome by affinity purification has yielded tremendous new insights into how genes are regulated. The affinity purification requires either the use of antibodies raised against the factor of interest itself or by high-affinity binding of a C- or N-terminally added tag sequence to the factor. Unfortunately, fusing extra amino acids to the termini of a factor can interfere with its biological function or the tag may be inaccessible inside the protein. Here, we describe an effective solution to that problem by integrating the ‘tag’ close to the nuclear localization sequence domain of the factor. We demonstrate the effectiveness of this approach with the transcription factors Fli-1 and Irf2bp2, which cannot be tagged at their extremities without loss of function. This resulted in the identification of novel proteins partners and a new hypothesis on the contribution of Fli-1 to hematopoiesis.

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

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.

Identification of multiple subsets of ventral interneurons and differential distribution along the rostrocaudal axis of the developing spinal cord

The spinal cord contains neuronal circuits termed Central Pattern Generators (CPGs) that coordinate rhythmic motor activities. CPG circuits consist of motor neurons and multiple interneuron cell types, many of which are derived from four distinct cardinal classes of ventral interneurons, called V0, V1, V2 and V3. While significant progress has been made on elucidating the molecular and genetic mechanisms that control ventral interneuron differentiation, little is known about their distribution along the antero-posterior axis of the spinal cord and their diversification. Here, we report that V0, V1 and V2 interneurons exhibit distinct organizational patterns at brachial, thoracic and lumbar levels of the developing spinal cord. In addition, we demonstrate that each cardinal class of ventral interneurons can be subdivided into several subsets according to the combinatorial expression of different sets of transcription factors, and that these subsets are differentially distributed along the rostrocaudal axis of the spinal cord. This comprehensive molecular profiling of ventral interneurons provides an important resource for investigating neuronal diversification in the developing spinal cord and for understanding the contribution of specific interneuron subsets on CPG circuits and motor control.

r3Cseq: an R/Bioconductor package for the discovery of long-range genomic interactions from chromosome conformation capture and next-generation sequencing data

The coupling of chromosome conformation capture (3C) with next-generation sequencing technologies enables the high-throughput detection of long-range genomic interactions, via the generation of ligation products between DNA sequences, which are closely juxtaposed in vivo. These interactions involve promoter regions, enhancers and other regulatory and structural elements of chromosomes and can reveal key details of the regulation of gene expression. 3C-seq is a variant of the method for the detection of interactions between one chosen genomic element (viewpoint) and the rest of the genome. We present r3Cseq, an R/Bioconductor package designed to perform 3C-seq data analysis in a number of different experimental designs. The package reads a common aligned read input format, provides data normalization, allows the visualization of candidate interaction regions and detects statistically significant chromatin interactions, thus greatly facilitating hypothesis generation and the interpretation of experimental results. We further demonstrate its use on a series of real-world applications.

In vivo live imaging of RNA polymerase II transcription factories in primary cells

Transcription steps are marked by different modifications of the C-terminal domain of RNA polymerase II (RNAPII). Phosphorylation of Ser5 and Ser7 by cyclin-dependent kinase 7 (CDK7) as part of TFIIH marks initiation, whereas phosphorylation of Ser2 by CDK9 marks elongation. These processes are thought to take place in localized transcription foci in the nucleus, known as "transcription factories," but it has been argued that the observed clusters/foci are mere fixation or labeling artifacts. We show that transcription factories exist in living cells as distinct foci by live-imaging fluorescently labeled CDK9, a kinase known to associate with active RNAPII. These foci were observed in different cell types derived from CDK9-mCherry knock-in mice. We show that these foci are very stable while highly dynamic in exchanging CDK9. Chromatin immunoprecipitation (ChIP) coupled with deep sequencing (ChIP-seq) data show that the genome-wide binding sites of CDK9 and initiating RNAPII overlap on transcribed genes. Immunostaining shows that CDK9-mCherry foci colocalize with RNAPII-Ser5P, much less with RNAPII-Ser2P, and not with CDK12 (a kinase reported to be involved in the Ser2 phosphorylation) or with splicing factor SC35. In conclusion, transcription factories exist in living cells, and initiation and elongation of transcripts takes place in different nuclear compartments.

KLF10 gene expression is associated with high fetal hemoglobin levels and with response to hydroxyurea treatment in β-hemoglobinopathy patients

AIM

In humans, fetal hemoglobin (HbF) production is controlled by many intricate mechanisms that, to date, remain only partly understood.

PATIENTS & METHODS

Pharmacogenomic analysis of the effects of hydroxyurea (HU) on HbF production was undertaken in a collection of Hellenic β-thalassemia and sickle cell disease (SCD) compound heterozygotes and a collection of healthy and KLF1-haploinsufficient Maltese adults, to identify genomic signatures that follow high HbF patterns.

RESULTS

KLF10 emerged as a top candidate. Moreover, genotype analysis of β-thalassemia major and intermedia patients and an independent cohort of β-thalassemia/SCD compound heterozygous patients that do or do not respond to HU treatment showed that the homozygous mutant state of a tagSNP in the KLF10 3'UTR is not present in β-thalassemia intermedia patients and is underrepresented in β-thalassemia/SCD compound heterozygous patients that respond well to HU treatment.

CONCLUSION

These data suggest that KLF10 may constitute a pharmacogenomic marker to discriminate between response and nonresponse to HU treatment.

Directed migration of cortical interneurons depends on the cell-autonomous action of Sip1

GABAergic interneurons mainly originate in the medial ganglionic eminence (MGE) of the embryonic ventral telencephalon (VT) and migrate tangentially to the cortex, guided by membrane-bound and secreted factors. We found that Sip1 (Zfhx1b, Zeb2), a transcription factor enriched in migrating cortical interneurons, is required for their proper differentiation and correct guidance. The majority of Sip1 knockout interneurons fail to migrate to the neocortex and stall in the VT. RNA sequencing reveals that Sip1 knockout interneurons do not acquire a fully mature cortical interneuron identity and contain increased levels of the repulsive receptor Unc5b. Focal electroporation of Unc5b-encoding vectors in the MGE of wild-type brain slices disturbs migration to the neocortex, whereas reducing Unc5b levels in Sip1 knockout slices and brains rescues the migration defect. Our results reveal that Sip1, through tuning of Unc5b levels, is essential for cortical interneuron guidance.

The mammalian gene function resource: the International Knockout Mouse Consortium

In 2007, the International Knockout Mouse Consortium (IKMC) made the ambitious promise to generate mutations in virtually every protein-coding gene of the mouse genome in a concerted worldwide action. Now, 5 years later, the IKMC members have developed high-throughput gene trapping and, in particular, gene-targeting pipelines and generated more than 17,400 mutant murine embryonic stem (ES) cell clones and more than 1,700 mutant mouse strains, most of them conditional. A common IKMC web portal (www.knockoutmouse.org) has been established, allowing easy access to this unparalleled biological resource. The IKMC materials considerably enhance functional gene annotation of the mammalian genome and will have a major impact on future biomedical research.

Super-resolution imaging reveals three-dimensional folding dynamics of the β-globin locus upon gene activation

The chromatin architecture is constantly changing because of cellular processes such as proliferation, differentiation and changes in the expression profile during gene activation or silencing. Unravelling the changes that occur in the chromatin structure during these processes has been a topic of interest for many years. It is known that gene activation of large gene loci is thought to occur by means of an active looping mechanism. It was also shown for the β-globin locus that the gene promoter interacts with an active chromatin hub by means of an active looping mechanism. This means that the locus changes in three-dimensional (3D) nuclear volume and chromatin shape. As a means of visualizing and measuring these dynamic changes in chromatin structure of the β-globin locus, we used a 3D DNA-FISH method in combination with 3D image acquisition to volume render fluorescent signals into 3D objects. These 3D chromatin structures were geometrically analysed, and results prior to and after gene activation were quantitatively compared. Confocal and super-resolution imaging revealed that the inactive locus occurs in several different conformations. These conformations change in shape and surface structure upon cell differentiation into a more folded and rounded structure that has a substantially smaller size and volume. These physical measurements represent the first non-biochemical evidence that, upon gene activation, an actively transcribing chromatin hub is formed by means of additional chromatin looping.

A new function of ROD1 in nonsense-mediated mRNA decay

RNA-binding proteins play a crucial role in the post-transcriptional regulation of gene expression. Polypyrimidine tract binding protein (PTB in humans) has been extensively characterized as an important splicing factor, and has additional functions in 3' end processing and translation. ROD1 is a PTB paralog containing four RRM (RNA recognition motif) domains. Here, we discover a function of ROD1 in nonsense-mediated mRNA decay (NMD). We show that ROD1 and the core NMD factor UPF1 interact and co-regulate an extensive number of target genes. Using a reporter system, we demonstrate that ROD1, similarly to UPF1 and UPF2, is required for the destabilization of a known NMD substrate. Finally, we show through RIP-seq that ROD1 and UPF1 associate with a significant number of common transcripts.

Enhancers and silencers: an integrated and simple model for their function

Regulatory DNA elements such as enhancers, silencers and insulators are embedded in metazoan genomes, and they control gene expression during development. Although they fulfil different roles, they share specific properties. Herein we discuss some examples and a parsimonious model for their function is proposed. All are transcription units that tether their target promoters close to, or distant from, transcriptional hot spots (or 'factories').

NARWHAL, a primary analysis pipeline for NGS data

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/.

Parallel high-performance grid computing: capabilities and opportunities of a novel demanding service and business class allowing highest resource efficiency

Especially in the life-science and the health-care sectors the huge IT requirements are imminent due to the large and complex systems to be analysed and simulated. Grid infrastructures play here a rapidly increasing role for research, diagnostics, and treatment, since they provide the necessary large-scale resources efficiently. Whereas grids were first used for huge number crunching of trivially parallelizable problems, increasingly parallel high-performance computing is required. Here, we show for the prime example of molecular dynamic simulations how the presence of large grid clusters including very fast network interconnects within grid infrastructures allows now parallel high-performance grid computing efficiently and thus combines the benefits of dedicated super-computing centres and grid infrastructures. The demands for this service class are the highest since the user group has very heterogeneous requirements: i) two to many thousands of CPUs, ii) different memory architectures, iii) huge storage capabilities, and iv) fast communication via network interconnects, are all needed in different combinations and must be considered in a highly dedicated manner to reach highest performance efficiency. Beyond, advanced and dedicated i) interaction with users, ii) the management of jobs, iii) accounting, and iv) billing, not only combines classic with parallel high-performance grid usage, but more importantly is also able to increase the efficiency of IT resource providers. Consequently, the mere "yes-we-can" becomes a huge opportunity like e.g. the life-science and health-care sectors as well as grid infrastructures by reaching higher level of resource efficiency.

Visualization, analysis, and design of COMBO-FISH probes in the grid-based GLOBE 3D genome platform

The genome architecture in cell nuclei plays an important role in modern microscopy for the monitoring of medical diagnosis and therapy since changes of function and dynamics of genes are interlinked with changing geometrical parameters. The planning of corresponding diagnostic experiments and their imaging is a complex and often interactive IT intensive challenge and thus makes high-performance grids a necessity. To detect genetic changes we recently developed a new form of fluorescence in situ hybridization (FISH) - COMBinatorial Oligonucleotide FISH (COMBO-FISH) - which labels small nucleotide sequences clustering at a desired genomic location. To achieve a unique hybridization spot other side clusters have to be excluded. Therefore, we have designed an interactive pipeline using the grid-based GLOBE 3D Genome Viewer and Platform to design and display different labelling variants of candidate probe sets. Thus, we have created a grid-based virtual "paper" tool for easy interactive calculation, analysis, management, and representation for COMBO-FISH probe design with many an advantage: Since all the calculations and analysis run in a grid, one can instantly and with great visual ease locate duplications of gene subsequences to guide the elimination of side clustering sequences during the probe design process, as well as get at least an impression of the 3D architectural embedding of the respective chromosome region, which is of major importance to estimate the hybridization probe dynamics. Beyond, even several people at different locations could work on the same process in a team wise manner. Consequently, we present how a complex interactive process can profit from grid infrastructure technology using our unique GLOBE 3D Genome Platform gateway towards a real interactive curative diagnosis planning and therapy monitoring.

GRIMP: a web- and grid-based tool for high-speed analysis of large-scale genome-wide association using imputed data

The current fast growth of genome-wide association studies (GWAS) combined with now common computationally expensive imputation requires the online access of large user groups to high-performance computing resources capable of analyzing rapidly and efficiently millions of genetic markers for ten thousands of individuals. Here, we present a web-based interface--called GRIMP--to run publicly available genetic software for extremely large GWAS on scalable super-computing grid infrastructures. This is of major importance for the enlargement of GWAS with the availability of whole-genome sequence data from the 1000 Genomes Project and for future whole-population efforts.

Fine-structured multi-scaling long-range correlations in completely sequenced genomes--features, origin, and classification

The sequential organization of genomes, i.e. the relations between distant base pairs and regions within sequences, and its connection to the three-dimensional organization of genomes is still a largely unresolved problem. Long-range power-law correlations were found using correlation analysis on almost the entire observable scale of 132 completely sequenced chromosomes of 0.5 × 10⁶ to 3.0 × 10⁷ bp from Archaea, Bacteria, Arabidopsis thaliana, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Drosophila melanogaster, and Homo sapiens. The local correlation coefficients show a species-specific multi-scaling behaviour: close to random correlations on the scale of a few base pairs, a first maximum from 40 to 3,400 bp (for Arabidopsis thaliana and Drosophila melanogaster divided in two submaxima), and often a region of one or more second maxima from 10⁵ to 3 × 10⁵ bp. Within this multi-scaling behaviour, an additional fine-structure is present and attributable to codon usage in all except the human sequences, where it is related to nucleosomal binding. Computer-generated random sequences assuming a block organization of genomes, the codon usage, and nucleosomal binding explain these results. Mutation by sequence reshuffling destroyed all correlations. Thus, the stability of correlations seems to be evolutionarily tightly controlled and connected to the spatial genome organization, especially on large scales. In summary, genomes show a complex sequential organization related closely to their three-dimensional organization.

Exportin 4 mediates a novel nuclear import pathway for Sox family transcription factors

SRY and other Sox-type transcription factors are important developmental regulators with various implications in human disease. In this study, we identified Exp4 (exportin 4) as an interaction partner of Sox2 in mouse embryonic stem cells and neural progenitors. We show that, besides its established function in nuclear export, Exp4 acts as a bona fide nuclear import receptor for Sox2 and SRY. Thus, Exp4 is an example of a nuclear transport receptor carrying distinct cargoes into different directions. In contrast to a published study, we observed that the import activity of Imp-alpha (importin-a) isoforms toward Sox2 is negligible. Instead, we found that Imp9 and the Imp-beta/7 heterodimer mediate nuclear import of Sox2 in parallel to Exp4. Import signals for the three pathways overlap and include conserved residues in the Sox2 high-mobility group (HMG) box domain that are also critical for DNA binding. This suggests that nuclear import of Sox proteins is facilitated by several parallel import pathways.

Exportin 4 mediates a novel nuclear import pathway for Sox family transcription factors

SRY and other Sox-type transcription factors are important developmental regulators with various implications in human disease. In this study, we identified Exp4 (exportin 4) as an interaction partner of Sox2 in mouse embryonic stem cells and neural progenitors. We show that, besides its established function in nuclear export, Exp4 acts as a bona fide nuclear import receptor for Sox2 and SRY. Thus, Exp4 is an example of a nuclear transport receptor carrying distinct cargoes into different directions. In contrast to a published study, we observed that the import activity of Imp-α (importin-a) isoforms toward Sox2 is negligible. Instead, we found that Imp9 and the Imp-β/7 heterodimer mediate nuclear import of Sox2 in parallel to Exp4. Import signals for the three pathways overlap and include conserved residues in the Sox2 high-mobility group (HMG) box domain that are also critical for DNA binding. This suggests that nuclear import of Sox proteins is facilitated by several parallel import pathways.

e-Human Grid Ecology - understanding and approaching the inverse tragedy of the commons in the e-Grid society

With ever-new technologies emerging also the amount of information to be stored and processed is growing exponentially and is believed to be always at the limit. In contrast, however, huge resources are available in the IT sector alike e.g. the renewable energy sector, which are often even not at all used. This under-usage bares any rational especially in the IT sector where e.g. virtualisation and grid approaches could be fast implemented due to the great technical and fast turnover opportunities. Here, we describe this obvious paradox for the first time as the Inverse Tragedy of the Commons, in contrast to the Classical Tragedy of the Commons where resources are overexploited. From this perspective the grid IT sector attempting to share resources for better efficiency, reveals two challenges leading to the heart of the paradox: i) From a macro perspective all grid infrastructures involve not only mere technical solutions but also dominantly all of the autopoietic social sub-systems ranging from religion to policy. ii) On the micro level the individual players and their psychology and risk behaviour are of major importance for acting within the macro autopoietic framework. Thus, the challenges of grid implementation are similar to those of e.g. climate protection. This is well described by the classic Human Ecology triangle and our extension to a rectangle: invironment-individual-society-environment. Extension of this classical interdisciplinary field of basic and applied research to an e-Human Grid Ecology rational, allows the Inverse Tragedy of the Commons of the grid sector to be understood and approached better and implies obvious guidelines in the day-to-day management for grid and other (networked) resources, which is of importance for many fields with similar paradoxes as in (e-)society.

The GLOBE 3D Genome Platform - towards a novel system-biological paper tool to integrate the huge complexity of genome organization and function

Genomes are tremendous co-evolutionary holistic systems for molecular storage, processing and fabrication of information. Their system-biological complexity remains, however, still largely mysterious, despite immense sequencing achievements and huge advances in the understanding of the general sequential, three-dimensional and regulatory organization. Here, we present the GLOBE 3D Genome Platform a completely novel grid based virtual "paper" tool and in fact the first system-biological genome browser integrating the holistic complexity of genomes in a single easy comprehensible platform: Based on a detailed study of biophysical and IT requirements, every architectural level from sequence to morphology of one or several genomes can be approached in a real and in a symbolic representation simultaneously and navigated by continuous scale-free zooming within a unique three-dimensional OpenGL and grid driven environment. In principle an unlimited number of multi-dimensional data sets can be visualized, customized in terms of arrangement, shape, colour, and texture etc. as well as accessed and annotated individually or in groups using internal or external data bases/facilities. Any information can be searched and correlated by importing or calculating simple relations in real-time using grid resources. A general correlation and application platform for more complex correlative analysis and a front-end for system-biological simulations both using again the huge capabilities of grid infrastructures is currently under development. Hence, the GLOBE 3D Genome Platform is an example of a grid based approach towards a virtual desktop for genomic work combining the three fundamental distributed resources: i) visual data representation, ii) data access and management, and iii) data analysis and creation. Thus, the GLOBE 3D Genome Platform is the novel system-biology oriented information system urgently needed to access, present, annotate, and to simulate the holistic genome complexity in a unique gateway towards a real understanding, educative presentation and curative manipulation planning of this tremendous evolutionary information grail - genomes.

The 3D structure of the immunoglobulin heavy-chain locus: implications for long-range genomic interactions

The immunoglobulin heavy-chain (Igh) locus is organized into distinct regions that contain multiple variable (V(H)), diversity (D(H)), joining (J(H)) and constant (C(H)) coding elements. How the Igh locus is structured in 3D space is unknown. To probe the topography of the Igh locus, spatial distance distributions were determined between 12 genomic markers that span the entire Igh locus. Comparison of the distance distributions to computer simulations of alternative chromatin arrangements predicted that the Igh locus is organized into compartments containing clusters of loops separated by linkers. Trilateration and triple-point angle measurements indicated the mean relative 3D positions of the V(H), D(H), J(H), and C(H) elements, showed compartmentalization and striking conformational changes involving V(H) and D(H)-J(H) elements during early B cell development. In pro-B cells, the entire repertoire of V(H) regions (2 Mbp) appeared to have merged and juxtaposed to the D(H) elements, mechanistically permitting long-range genomic interactions to occur with relatively high frequency.

Increased gamma-globin gene expression in beta-thalassemia intermedia patients correlates with a mutation in 3'HS1

We report a novel set of genetic markers in the DNaseI hypersensitive sites comprising the human beta-globin locus chromatin hub (CH), namely HS-111 and 3'HS1. The HS-111 (-21 G>A) and 3'HS1 (+179 C>T) transitions form CH haplotypes, which occur at different frequencies in beta-thalassemia intermedia and major patients and normal (nonthalassemic) individuals. We also show that the 3'HS1 (+179 C>T) variation results in a GATA-1 binding site and correlates with increased fetal hemoglobin production in beta-thalassemia intermedia patients. In contrast, the HS-111 (+126 G>A) transition, found in three normal chromosomes, is simply a rare polymorphism. We conclude that the CH haplotypes are useful genetic determinants for beta-thalassemia major and intermedia patients, while the 3'HS1 (+179 C>T) mutation may have functional consequences in gamma-globin genes expression.

Gata-3 is an essential regulator of mammary-gland morphogenesis and luminal-cell differentiation

The transcription factor Gata-3 is a defining marker of the 'luminal' subtypes of breast cancer. To gain insight into the role of Gata-3 in breast epithelial development and oncogenesis, we have explored its normal function within the mammary gland by conditionally deleting Gata-3 at different stages of development. We report that Gata-3 has essential roles in the morphogenesis of the mammary gland in both the embryo and adult. Through the discovery of a novel marker (beta3-integrin) of luminal progenitor cells and their purification, we demonstrate that Gata-3 deficiency leads to an expansion of luminal progenitors and a concomitant block in differentiation. Remarkably, introduction of Gata-3 into a stem cell-enriched population induced maturation along the alveolar luminal lineage. These studies provide evidence for the existence of an epithelial hierarchy within the mammary gland and establish Gata-3 as a critical regulator of luminal differentiation.

Transcriptome and phenotypic analysis reveals Gata3-dependent signalling pathways in murine hair follicles

The transcription factor Gata3 is crucially involved in epidermis and hair follicle differentiation. Yet, little is known about how Gata3 co-ordinates stem cell lineage determination in skin, what pathways are involved and how Gata3 differentially regulates distinct cell populations within the hair follicle. Here, we describe a conditional Gata3-/- mouse (K14-Gata3-/-) in which Gata3 is specifically deleted in epidermis and hair follicles. K14-Gata3-/- mice show aberrant postnatal growth and development, delayed hair growth and maintenance, abnormal hair follicle organization and irregular pigmentation. After the first hair cycle, the germinative layer surrounding the dermal papilla was not restored; instead, proliferation was pronounced in basal epidermal cells. Transcriptome analysis of laser-dissected K14-Gata3-/- hair follicles revealed mitosis, epithelial differentiation and the Notch, Wnt and BMP signaling pathways to be significantly overrepresented. Elucidation of these pathways at the RNA and protein levels and physiologic endpoints suggests that Gata3 integrates diverse signaling networks to regulate the balance between hair follicle and epidermal cell fates.

Multiple interactions between regulatory regions are required to stabilize an active chromatin hub

The human beta-globin locus control region (LCR) is required for the maintenance of an open chromatin configuration of the locus. It interacts with the genes and the hypersensitive regions flanking the locus to form an active chromatin hub (ACH) transcribing the genes. Proper developmental control of globin genes is largely determined by gene proximal regulatory sequences. Here, we provide the first functional evidence of the role of the most active sites of the LCR and the promoter of the beta-globin gene in the maintenance of the ACH. When the human beta-globin gene promoter is deleted in the context of a full LCR, the ACH is maintained with the beta-globin gene remaining in proximity. Additional deletion of hypersensitive site HS3 or HS2 of the LCR shows that HS3, but not HS2, in combination with the beta-globin promoter is crucial for the maintenance of the ACH at the definitive stage. We conclude that multiple interactions between the LCR and the beta-globin gene are required to maintain the appropriate spatial configuration in vivo.

Hearing loss following Gata3 haploinsufficiency is caused by cochlear disorder

Patients with HDR syndrome suffer from hypoparathyroidism, deafness, and renal dysplasia due to a heterozygous deletion of the transcription factor GATA3. Since GATA3 is prominently expressed in both the inner ear and different parts of the auditory nervous system, it is not clear whether the deafness in HDR patients is caused by peripheral and/or central deficits. Therefore, we have created and examined heterozygous Gata3 knockout mice. Auditory brainstem response (ABR) thresholds of alert heterozygous Gata3 mice, analyzed from 1 to 19 months of age, showed a hearing loss of 30 dB compared to wild-type littermates. Neither physiological nor morphological abnormalities were found in the brainstem, cerebral cortex, the outer or the middle ear. In contrast, cochleae of heterozygous Gata3 mice showed significant progressive morphological degeneration starting with the outer hair cells (OHCs) at the apex and ultimately affecting all hair cells and supporting cells in the entire cochlea. Together, these findings indicate that hearing loss following Gata3 haploinsufficiency is peripheral in origin and that this defect is detectable from early postnatal development and maintains through adulthood.

Fluorescence in situ hybridization analysis of transcript dynamics in cells

Since the first description of in situ hybridization in 1969 the technique has advanced to allow sensitive detection of DNA and mRNA molecules at the cellular and subcellular levels. In particular fluorescence in situ hybridization (FISH) has become a frequently used tool in basic and applied biomedical research since detection is sensitive and allows discrimination of multiple targets in the same sample. By using RNA-FISH we have been able to detect primary transcripts of the human embryonic, fetal, and adult globins in erythroid cells to study the competitive transcription mechanism or variegated expression patterns of the human beta-globin locus. We have correlated such expression patterns with other parameters such as cell type, cell cycle, replication, and stage of differentiation by simultaneous detection of, e.g., incorporated BrdUTPs, proteins (e.g., cyclins A and E, PCNA, histones), and globin (primary) transcripts and/or locus integration sites. Thus a combination of FISH and immunofluorescence methods allow the visualization of different processes taking place in the nucleus relative to each other in terms of three-dimensional space and structure and time (development, cell cycle).

Modification of human beta-globin locus PAC clones by homologous recombination in Escherichia coli

We report here modifications of human beta-globin PAC clones by homologous recombination in Escherichia coli DH10B, utilising a plasmid temperature sensitive for replication, the recA gene and a wild-type copy of the rpsL gene which allows for an efficient selection for plasmid loss in this host. High frequencies of recombination are observed even with very small lengths of homology and the method has general utility for introducing insertions, deletions and point mutations. No rearrangements were detected with the exception of one highly repetitive genomic sequence when either the E.COLI: RecA- or the lambdoid phage encoded RecT and RecE-dependent recombination systems were used.

Cloning and sequencing of the CRABP-I locus from chicken and pufferfish: analysis of the promoter regions in transgenic mice

Retinoic acid (RA), a derivative of vitamin A, is an important molecule for development and homeostasis of vertebrate organisms. The intracellular retinoic acid binding protein CRABP-I has a high affinity for RA, and is thought to be involved in the mechanism of RA signalling. CRABP-I is well conserved in evolution and shows a specific expression pattern during development, but mice made deficient for the protein by gene targeting appear normal. However, the high degree of homology with CRABP-I from other species indicates that the protein has been subject to strong selective conservation, indicative of an important biological function. In this paper we have compared the conservation in the expression pattern of the mouse, chicken and pufferfish CRABP-I genes to substantiate this argument further. First we cloned and sequenced genes and promoter regions of the CRABP-I genes from chicken and the Japanese pufferfish, Fugu rubripes. Sequence comparison with the mouse gene did not show any large blocks of homology in the promoter regions. Nevertheless, the promoter of the chicken gene directed expression to a subset of the tissues that show expression with the promoter from the mouse gene. The pattern observed with the pufferfish promoter is even more restricted, essentially to rhombomere 4 only, indicating that this region may be functionally the most important for CRABP-I expression in the developing embryo.

Frameshift mutants of beta amyloid precursor protein and ubiquitin-B in Alzheimer's and Down patients

The cerebral cortex of Alzheimer's and Down syndrome patients is characterized by the presence of protein deposits in neurofibrillary tangles, neuritic plaques, and neuropil threads. These structures were shown to contain forms of beta amyloid precursor protein and ubiquitin-B that are aberrant (+1 proteins) in the carboxyl terminus. The +1 proteins were not found in young control patients, whereas the presence of ubiquitin-B+1 in elderly control patients may indicate early stages of neurodegeneration. The two species of +1 proteins displayed cellular colocalization, suggesting a common origin, operating at the transcriptional level or by posttranscriptional editing of RNA. This type of transcript mutation is likely an important factor in the widely occurring nonfamilial early- and late-onset forms of Alzheimer's disease.

The chromosome make-up of mouse embryonic stem cells is predictive of somatic and germ cell chimaerism

Mouse pluripotent embryonic stem (ES) cells, once reintroduced into a mouse blastocyst, can contribute to the formation of all tissues, including the germline, of an organism referred to as a chimaeric. However, the reasons why this contribution often appears erratic are poorly understood. We have tested the notion that the chromosome make-up may be important in contributing both to somatic cell chimaerism and to germ line transmission. We found that the percentage of chimaerism of ES cell-embryo chimaeras, the absolute number of chimaeras and the ratio of chimaeras to total pups born all correlate closely with the percentage of euploid metaphases in the ES cell clones injected into the murine blastocyst. The majority of the ES cell clones that we tested, which were obtained from different gene targeting knockout experiments and harboured 50 to 100% euploid metaphases, did transmit to the germline; in contrast, none of the ES cell clones with more than 50% of chromosomally abnormal metaphases transmitted to the germline. Euploid ES cell clones cultured in vitro for more than 20 passages rapidly became severely aneuploid, and again this correlated closely with the percentage of chimaerism and with the number of ES cell-embryo chimaeras obtained per number of blastocysts injected. At the same time, the ability of these clones to contribute to the germline was lost when the proportion of euploid cells dropped below 50%. This study suggests that aneuploidy, rather than 'loss of totipotency', in ES cells, is the major cause of failure in obtaining contributions to all tissues of the adult chimaera, including the germline. Because euploidy is predictive of germline transmission, karyotype analysis is crucial and time/cost saving in any gene-targeting experiment.

GATA transcription factors associate with a novel class of nuclear bodies in erythroblasts and megakaryocytes

The nuclear distribution of GATA transcription factors in murine haemopoietic cells was examined by indirect immunofluorescence. Specific bright foci of GATA-1 fluorescence were observed in erythroleukaemia cells and primary murine erythroblasts and megakaryocytes, in addition to diffuse nucleoplasmic localization. These foci, which were preferentially found adjacent to nucleoli or at the nuclear periphery, did not represent sites of active transcription or binding of GATA-1 to consensus sites in the beta-globin loci. Immunoelectron microscopy demonstrated the presence of intensely labelled structures likely to represent the GATA-1 foci seen by immunofluorescence. The GATA-1 nuclear bodies differed from previously described nuclear structures and there was no co-localization with nuclear antigens involved in RNA processing or other ubiquitous (Spl, c-Jun and TBP) or haemopoietic (NF-E2) transcription factors. Interestingly, GATA-2 and GATA-3 proteins also localized to the same nuclear bodies in cell lines co-expressing GATA-1 and -2 or GATA-1 and -3 gene products. This pattern of distribution is, thus far, unique to the GATA transcription factors and suggests a protein-protein interaction with other components of the nuclear bodies via the GATA zinc finger domain.

Mouse model for the lysosomal disorder galactosialidosis and correction of the phenotype with overexpressing erythroid precursor cells

The lysosomal storage disorder galactosialidosis results from a primary deficiency of the protective protein/cathepsin A (PPCA), which in turn affects the activities of beta-galactosidase and neuraminidase. Mice homozygous for a null mutation at the PPCA locus present with signs of the disease shortly after birth and develop a phenotype closely resembling human patients with galactosialidosis. Most of their tissues show characteristic vacuolation of specific cells, attributable to lysosomal storage. Excessive excretion of sialyloligosaccharides in urine is diagnostic of the disease. Affected mice progressively deteriorate as a consequence of severe organ dysfunction, especially of the kidney. The deficient phenotype can be corrected by transplanting null mutants with bone marrow from a transgenic line overexpressing human PPCA in erythroid precursor cells. The transgenic bone marrow gives a more efficient and complete correction of the visceral organs than normal bone marrow. Our data demonstrate the usefulness of this animal model, very similar to the human disease, for experimenting therapeutic strategies aimed to deliver the functional protein or gene to affected organs. Furthermore, they suggest the feasibility of gene therapy for galactosialidosis and other disorders, using bone marrow cells engineered to overexpress and secrete the correcting lysosomal protein.

Targeted disruption of the GATA3 gene causes severe abnormalities in the nervous system and in fetal liver haematopoiesis

GATA-3 is one member of a growing family of related transcription factors which share a strongly conserved expression pattern in all vertebrate organisms. In order to elucidate GATA-3 function using a direct genetic approach, we have disrupted the murine gene by homologous recombination in embryonic stem cells. Mice heterozygous for the GATA3 mutation are fertile and appear in all respects to be normal, whereas homozygous mutant embryos die between days 11 and 12 postcoitum (p.c.) and display massive internal bleeding, marked growth retardation, severe deformities of the brain and spinal cord, and gross aberrations in fetal liver haematopoiesis.