Specific photoreceptor cell fate pathways are differentially altered in NR2E3-associated diseases

Mutations in NR2E3, a gene encoding an orphan nuclear transcription factor, cause two retinal dystrophies with a distinct phenotype, but the precise role of NR2E3 in rod and cone transcriptional networks remains unclear. To dissect NR2E3 function, we performed scRNA-seq in the retinas of wildtype and two different Nr2e3 mouse models that show phenotypes similar to patients carrying NR2E3 mutations. Our results reveal that rod and cone populations are not homogeneous and can be separated into different sub-classes. We identify a previously unreported cone pathway that generates hybrid cones co-expressing both cone- and rod-related genes. In mutant retinas, this hybrid cone subpopulation is more abundant and includes a subpopulation of rods transitioning towards a cone cell fate. Hybrid photoreceptors with high misexpression of cone- and rod-related genes are prone to regulated necrosis. Overall, our results shed light on the role of NR2E3 in modulating photoreceptor differentiation towards cone and rod fates and explain how different mutations in NR2E3 lead to distinct visual disorders in humans.

Doxorubicin Changes the Spatial Organization of the Genome around Active Promoters

In this study, we delve into the impact of genotoxic anticancer drug treatment on the chromatin structure of human cells, with a particular focus on the effects of doxorubicin. Using Hi-C, ChIP-seq, and RNA-seq, we explore the changes in chromatin architecture brought about by doxorubicin and ICRF193. Our results indicate that physiologically relevant doses of doxorubicin lead to a local reduction in Hi-C interactions in certain genomic regions that contain active promoters, with changes in chromatin architecture occurring independently of Top2 inhibition, cell cycle arrest, and differential gene expression. Inside the regions with decreased interactions, we detected redistribution of RAD21 around the peaks of H3K27 acetylation. Our study also revealed a common structural pattern in the regions with altered architecture, characterized by two large domains separated from each other. Additionally, doxorubicin was found to increase CTCF binding in H3K27 acetylated regions. Furthermore, we discovered that Top2-dependent chemotherapy causes changes in the distance decay of Hi-C contacts, which are driven by direct and indirect inhibitors. Our proposed model suggests that doxorubicin-induced DSBs cause cohesin redistribution, which leads to increased insulation on actively transcribed TAD boundaries. Our findings underscore the significant impact of genotoxic anticancer treatment on the chromatin structure of the human genome.

Early developmental plasticity enables the induction of an intermediate extraembryonic cell state

Two fundamental elements of pre-implantation embryogenesis are cells' intrinsic self-organization program and their developmental plasticity, which allows embryos to compensate for alterations in cell position and number; yet, these elements are still poorly understood. To be able to decipher these features, we established culture conditions that enable the two fates of blastocysts' extraembryonic lineages-the primitive endoderm and the trophectoderm-to coexist. This plasticity emerges following the mechanisms of the first lineage segregation in the mouse embryo, and it manifests as an extended potential for extraembryonic chimerism during the pre-implantation embryogenesis. Moreover, this shared state enables robust assembly into higher-order blastocyst-like structures, thus combining both the cell fate plasticity and self-organization features of the early extraembryonic lineages.

Multiomic atlas with functional stratification and developmental dynamics of zebrafish cis-regulatory elements

Zebrafish, a popular organism for studying embryonic development and for modeling human diseases, has so far lacked a systematic functional annotation program akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created a central repository to store and process zebrafish developmental functional genomic data. Our data coordination center ( https://danio-code.zfin.org ) combines a total of 1,802 sets of unpublished and re-analyzed published genomic data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements throughout development, including classes with distinct features dependent on their activity in time and space. We delineated the distinct distance topology and chromatin features between regulatory elements active during zygotic genome activation and those active during organogenesis. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predicted functional relationships between them beyond sequence similarity, thus extending the utility of zebrafish developmental genomics to mammals.

Shear stress switches the association of endothelial enhancers from ETV/ETS to KLF transcription factor binding sites

Endothelial cells (ECs) lining blood vessels are exposed to mechanical forces, such as shear stress. These forces control many aspects of EC biology, including vascular tone, cell migration and proliferation. Despite a good understanding of the genes responding to shear stress, our insight into the transcriptional regulation of these genes is much more limited. Here, we set out to study alterations in the chromatin landscape of human umbilical vein endothelial cells (HUVEC) exposed to laminar shear stress. To do so, we performed ChIP-Seq for H3K27 acetylation, indicative of active enhancer elements and ATAC-Seq to mark regions of open chromatin in addition to RNA-Seq on HUVEC exposed to 6 h of laminar shear stress. Our results show a correlation of gained and lost enhancers with up and downregulated genes, respectively. DNA motif analysis revealed an over-representation of KLF transcription factor (TF) binding sites in gained enhancers, while lost enhancers contained more ETV/ETS motifs. We validated a subset of flow responsive enhancers using luciferase-based reporter constructs and CRISPR-Cas9 mediated genome editing. Lastly, we characterized the shear stress response in ECs of zebrafish embryos using RNA-Seq. Our results lay the groundwork for the exploration of shear stress responsive elements in controlling EC biology.

Lima1 mediates the pluripotency control of membrane dynamics and cellular metabolism

Lima1 is an extensively studied prognostic marker of malignancy and is also considered to be a tumour suppressor, but its role in a developmental context of non-transformed cells is poorly understood. Here, we characterise the expression pattern and examined the function of Lima1 in mouse embryos and pluripotent stem cell lines. We identify that Lima1 expression is controlled by the naïve pluripotency circuit and is required for the suppression of membrane blebbing, as well as for proper mitochondrial energetics in embryonic stem cells. Moreover, forcing Lima1 expression enables primed mouse and human pluripotent stem cells to be incorporated into murine pre-implantation embryos. Thus, Lima1 is a key effector molecule that mediates the pluripotency control of membrane dynamics and cellular metabolism.

Zebrafish transposable elements show extensive diversification in age, genomic distribution, and developmental expression

There is considerable interest in understanding the effect of transposable elements (TEs) on embryonic development. Studies in humans and mice are limited by the difficulty of working with mammalian embryos and by the relative scarcity of active TEs in these organisms. The zebrafish is an outstanding model for the study of vertebrate development, and over half of its genome consists of diverse TEs. However, zebrafish TEs remain poorly characterized. Here we describe the demography and genomic distribution of zebrafish TEs and their expression throughout embryogenesis using bulk and single-cell RNA sequencing data. These results reveal a highly dynamic genomic ecosystem comprising nearly 2000 distinct TE families, which vary in copy number by four orders of magnitude and span a wide range of ages. Longer retroelements tend to be retained in intergenic regions, whereas short interspersed nuclear elements (SINEs) and DNA transposons are more frequently found nearby or within genes. Locus-specific mapping of TE expression reveals extensive TE transcription during development. Although two-thirds of TE transcripts are likely driven by nearby gene promoters, we still observe stage- and tissue-specific expression patterns in self-regulated TEs. Long terminal repeat (LTR) retroelements are most transcriptionally active immediately following zygotic genome activation, whereas DNA transposons are enriched among transcripts expressed in later stages of development. Single-cell analysis reveals several endogenous retroviruses expressed in specific somatic cell lineages. Overall, our study provides a valuable resource for using zebrafish as a model to study the impact of TEs on vertebrate development.

Ronin governs the metabolic capacity of the embryonic lineage for post-implantation development

During implantation, the murine embryo transitions from a “quiet” into an active metabolic/proliferative state, which kick‐starts the growth and morphogenesis of the post‐implantation conceptus. Such transition is also required for embryonic stem cells to be established from mouse blastocysts, but the factors regulating this process are poorly understood. Here, we show that Ronin plays a critical role in the process by enabling active energy production, and the loss of Ronin results in the establishment of a reversible quiescent state in which naïve pluripotency is promoted. In addition, Ronin fine‐tunes the expression of genes that encode ribosomal proteins and is required for proper tissue‐scale organisation of the pluripotent lineage during the transition from blastocyst to egg cylinder stage. Thus, Ronin function is essential for governing the metabolic capacity so that it can support the pluripotent lineage’s high‐energy demands for cell proliferation and morphogenesis.

Chromatin architecture transitions from zebrafish sperm through early embryogenesis

Chromatin architecture mapping in 3D formats has increased our understanding of how regulatory sequences and gene expression are connected and regulated in a genome. The 3D chromatin genome shows extensive remodeling during embryonic development, and although the cleavage-stage embryos of most species lack structure before zygotic genome activation (pre-ZGA), zebrafish has been reported to have structure. Here, we aimed to determine the chromosomal architecture in paternal/sperm zebrafish gamete cells to discern whether it either resembles or informs early pre-ZGA zebrafish embryo chromatin architecture. First, we assessed the higher-order architecture through advanced low-cell in situ Hi-C. The structure of zebrafish sperm, packaged by histones, lacks topological associated domains and instead displays “hinge-like” domains of ∼150 kb that repeat every 1–2 Mbs, suggesting a condensed repeating structure resembling mitotic chromosomes. The pre-ZGA embryos lacked chromosomal structure, in contrast to prior work, and only developed structure post-ZGA. During post-ZGA, we find chromatin architecture beginning to form at small contact domains of a median length of ∼90 kb. These small contact domains are established at enhancers, including super-enhancers, and chemical inhibition of Ep300a (p300) and Crebbpa (CBP) activity, lowering histone H3K27ac, but not transcription inhibition, diminishes these contacts. Together, this study reveals hinge-like domains in histone-packaged zebrafish sperm chromatin and determines that the initial formation of high-order chromatin architecture in zebrafish embryos occurs after ZGA primarily at enhancers bearing high H3K27ac.

Repression of endogenous retroviruses prevents antiviral immune response and is required for mammary gland development

The role of heterochromatin in cell fate specification during development is unclear. We demonstrate that loss of the lysine 9 of histone H3 (H3K9) methyltransferase G9a in the mammary epithelium results in de novo chromatin opening, aberrant formation of the mammary ductal tree, impaired stem cell potential, disrupted intraductal polarity, and loss of tissue function. G9a loss derepresses long terminal repeat (LTR) retroviral sequences (predominantly the ERVK family). Transcriptionally activated endogenous retroviruses generate double-stranded DNA (dsDNA) that triggers an antiviral innate immune response, and knockdown of the cytosolic dsDNA sensor Aim2 in G9a knockout (G9acKO) mammary epithelium rescues mammary ductal invasion. Mammary stem cell transplantation into immunocompromised or G9acKO-conditioned hosts shows partial dependence of the G9acKO mammary morphological defects on the inflammatory milieu of the host mammary fat pad. Thus, altering the chromatin accessibility of retroviral elements disrupts mammary gland development and stem cell activity through both cell-autonomous and non-autonomous mechanisms.

Germ cell differentiation requires Tdrd7-dependent chromatin and transcriptome reprogramming marked by germ plasm relocalization

In many animal models, primordial germ cell (PGC) development depends on maternally deposited germ plasm, which prevents somatic cell fate. Here, we show that PGCs respond to regulatory information from the germ plasm in two distinct phases using two distinct mechanisms in zebrafish. We demonstrate that PGCs commence zygotic genome activation together with the somatic blastocysts with no demonstrable differences in transcriptional and chromatin opening. Unexpectedly, both PGC and somatic blastocysts activate germ-cell-specific genes, which are only stabilized in PGCs by cytoplasmic germ plasm determinants. Disaggregated perinuclear relocalization of germ plasm during PGC migration is regulated by the germ plasm determinant Tdrd7 and is coupled to dramatic divergence between PGC and somatic transcriptomes. This transcriptional divergence relies on PGC-specific cis-regulatory elements characterized by promoter-proximal distribution. We show that Tdrd7-dependent reconfiguration of chromatin accessibility is required for elaboration of PGC fate but not for PGC migration.

•

No evidence for transcriptional activation delay in zebrafish PGCs

•

Germ-plasm-associated post-transcriptional divergence during ZGA

•

Epigenetic reprogramming marks onset of PGC migration

•

Epigenetic reprogramming in PGCs relies on Tdrd7, coupled to germ plasm relocalization

FAN-C: a feature-rich framework for the analysis and visualisation of chromosome conformation capture data

Chromosome conformation capture data, particularly from high-throughput approaches such as Hi-C, are typically very complex to analyse. Existing analysis tools are often single-purpose, or limited in compatibility to a small number of data formats, frequently making Hi-C analyses tedious and time-consuming. Here, we present FAN-C, an easy-to-use command-line tool and powerful Python API with a broad feature set covering matrix generation, analysis, and visualisation for C-like data ( https://github.com/vaquerizaslab/fanc ). Due to its compatibility with the most prevalent Hi-C storage formats, FAN-C can be used in combination with a large number of existing analysis tools, thus greatly simplifying Hi-C matrix analysis.

The Transition from Quiescent to Activated States in Human Hematopoietic Stem Cells Is Governed by Dynamic 3D Genome Reorganization

Lifelong blood production requires long-term hematopoietic stem cells (LT-HSCs), marked by stemness states involving quiescence and self-renewal, to transition into activated short-term HSCs (ST-HSCs) with reduced stemness. As few transcriptional changes underlie this transition, we used single-cell and bulk assay for transposase-accessible chromatin sequencing (ATAC-seq) on human HSCs and hematopoietic stem and progenitor cell (HSPC) subsets to uncover chromatin accessibility signatures, one including LT-HSCs (LT/HSPC signature) and another excluding LT-HSCs (activated HSPC [Act/HSPC] signature). These signatures inversely correlated during early hematopoietic commitment and differentiation. The Act/HSPC signature contains CCCTC-binding factor (CTCF) binding sites mediating 351 chromatin interactions engaged in ST-HSCs, but not LT-HSCs, enclosing multiple stemness pathway genes active in LT-HSCs and repressed in ST-HSCs. CTCF silencing derepressed stemness genes, restraining quiescent LT-HSCs from transitioning to activated ST-HSCs. Hence, 3D chromatin interactions centrally mediated by CTCF endow a gatekeeper function that governs the earliest fate transitions HSCs make by coordinating disparate stemness pathways linked to quiescence and self-renewal.

CHESS enables quantitative comparison of chromatin contact data and automatic feature extraction

Dynamic changes in the three-dimensional (3D) organization of chromatin are associated with central biological processes, such as transcription, replication and development. Therefore, the comprehensive identification and quantification of these changes is fundamental to understanding of evolutionary and regulatory mechanisms. Here, we present Comparison of Hi-C Experiments using Structural Similarity (CHESS), an algorithm for the comparison of chromatin contact maps and automatic differential feature extraction. We demonstrate the robustness of CHESS to experimental variability and showcase its biological applications on (1) interspecies comparisons of syntenic regions in human and mouse models; (2) intraspecies identification of conformational changes in Zelda-depleted Drosophila embryos; (3) patient-specific aberrant chromatin conformation in a diffuse large B-cell lymphoma sample; and (4) the systematic identification of chromatin contact differences in high-resolution Capture-C data. In summary, CHESS is a computationally efficient method for the comparison and classification of changes in chromatin contact data.

Cohesin Disrupts Polycomb-Dependent Chromosome Interactions in Embryonic Stem Cells

How chromosome organization is related to genome function remains poorly understood. Cohesin, loop extrusion, and CCCTC-binding factor (CTCF) have been proposed to create topologically associating domains (TADs) to regulate gene expression. Here, we examine chromosome conformation in embryonic stem cells lacking cohesin and find, as in other cell types, that cohesin is required to create TADs and regulate A/B compartmentalization. However, in the absence of cohesin, we identify a series of long-range chromosomal interactions that persist. These correspond to regions of the genome occupied by the polycomb repressive system and are dependent on PRC1. Importantly, we discover that cohesin counteracts these polycomb-dependent interactions, but not interactions between super-enhancers. This disruptive activity is independent of CTCF and insulation and appears to modulate gene repression by the polycomb system. Therefore, we discover that cohesin disrupts polycomb-dependent chromosome interactions to modulate gene expression in embryonic stem cells.

Endothelial EphB4 maintains vascular integrity and transport function in adult heart

The homeostasis of heart and other organs relies on the appropriate provision of nutrients and functional specialization of the local vasculature. Here, we have used mouse genetics, imaging and cell biology approaches to investigate how homeostasis in the adult heart is controlled by endothelial EphB4 and its ligand ephrin-B2, which are known regulators of vascular morphogenesis and arteriovenous differentiation during development. We show that inducible and endothelial cell-specific inactivation of Ephb4 in adult mice is compatible with survival, but leads to rupturing of cardiac capillaries, cardiomyocyte hypertrophy, and pathological cardiac remodeling. In contrast, EphB4 is not required for integrity and homeostasis of capillaries in skeletal muscle. Our analysis of mutant mice and cultured endothelial cells shows that EphB4 controls the function of caveolae, cell-cell adhesion under mechanical stress and lipid transport. We propose that EphB4 maintains critical functional properties of the adult cardiac vasculature and thereby prevents dilated cardiomyopathy-like defects.

Visualising three-dimensional genome organisation in two dimensions

The three-dimensional organisation of the genome plays a crucial role in developmental gene regulation. In recent years, techniques to investigate this organisation have become more accessible to labs worldwide due to improvements in protocols and decreases in the cost of high-throughput sequencing. However, the resulting datasets are complex and can be challenging to analyse and interpret. Here, we provide a guide to visualisation approaches that can aid the interpretation of such datasets and the communication of biological results.

The Birth of the 3D Genome during Early Embryonic Development

The 3D structure of chromatin in the nucleus is important for the regulation of gene expression and the correct deployment of developmental programs. The differentiation of germ cells and early embryonic development (when the zygotic genome is activated and transcription is taking place for the first time) are accompanied by dramatic changes in gene expression and the epigenetic landscape. Recent studies used Hi-C to investigate the 3D chromatin organization during these developmental transitions, uncovering remarkable remodeling of the 3D genome. Here, we highlight the changes described so far and discuss some of the implications that these findings have for our understanding of the mechanisms and functionality of 3D chromatin architecture.

Generation of Genome-wide Chromatin Conformation Capture Libraries from Tightly Staged Early Drosophila Embryos

Investigating the three-dimensional architecture of chromatin offers invaluable insight into the mechanisms of gene regulation. Here, we describe a protocol for performing the chromatin conformation capture technique in situ Hi-C on staged Drosophila melanogaster embryo populations. The result is a sequencing library that allows the mapping of all chromatin interactions that occur in the nucleus in a single experiment. Embryo sorting is done manually using a fluorescent stereo microscope and a transgenic fly line containing a nuclear marker. Using this technique, embryo populations from each nuclear division cycle, and with defined cell cycle status, can be obtained with very high purity. The protocol may also be adapted to sort older embryos beyond gastrulation. Sorted embryos are used as inputs for in situ Hi-C. All experiments, including sequencing library preparation, can be completed in five days. The protocol has low input requirements and works reliably using 20 blastoderm stage embryos as input material. The end result is a sequencing library for next generation sequencing. After sequencing, the data can be processed into genome-wide chromatin interaction maps that can be analyzed using a wide range of available tools to gain information about topologically associating domain (TAD) structure, chromatin loops, and chromatin compartments during Drosophila development.

Xrp1 genetically interacts with the ALS-associated FUS orthologue caz and mediates its toxicity

Mallik et al. identify Xrp1 as a nuclear chromatin-binding protein involved in gene expression regulation that mediates phenotypes induced by loss of function of cabeza (caz), the Drosophila melanogaster orthologue of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) protein FUS. Knockdown of Xrp1 in motor neurons rescues phenotypes induced by ALS-mutant FUS.

Cabeza (caz) is the single Drosophila melanogaster orthologue of the human FET proteins FUS, TAF15, and EWSR1, which have been implicated in amyotrophic lateral sclerosis (ALS) and frontotemporal dementia. In this study, we identified Xrp1, a nuclear chromatin-binding protein, as a key modifier of caz mutant phenotypes. Xrp1 expression was strongly up-regulated in caz mutants, and Xrp1 heterozygosity rescued their motor defects and life span. Interestingly, selective neuronal Xrp1 knockdown was sufficient to rescue, and neuronal Xrp1 overexpression phenocopied caz mutant phenotypes. The caz/Xrp1 genetic interaction depended on the functionality of the AT-hook DNA-binding domain in Xrp1, and the majority of Xrp1-interacting proteins are involved in gene expression regulation. Consistently, caz mutants displayed gene expression dysregulation, which was mitigated by Xrp1 heterozygosity. Finally, Xrp1 knockdown substantially rescued the motor deficits and life span of flies expressing ALS mutant FUS in motor neurons, implicating gene expression dysregulation in ALS-FUS pathogenesis.

Chromatin Architecture Emerges during Zygotic Genome Activation Independent of Transcription

Chromatin architecture is fundamental in regulating gene expression. To investigate when spatial genome organization is first established during development, we examined chromatin conformation during Drosophila embryogenesis and observed the emergence of chromatin architecture within a tight time window that coincides with the onset of transcription activation in the zygote. Prior to zygotic genome activation, the genome is mostly unstructured. Early expressed genes serve as nucleation sites for topologically associating domain (TAD) boundaries. Activation of gene expression coincides with the establishment of TADs throughout the genome and co-localization of housekeeping gene clusters, which remain stable in subsequent stages of development. However, the appearance of TAD boundaries is independent of transcription and requires the transcription factor Zelda for locus-specific TAD boundary insulation. These results offer insight into when spatial organization of the genome emerges and identify a key factor that helps trigger this architecture.

Dual randomization of oligonucleotides to reduce the bias in ribosome-profiling libraries

Protein translation is at the heart of cellular metabolism and its in-depth characterization is key for many lines of research. Recently, ribosome profiling became the state-of-the-art method to quantitatively characterize translation dynamics at a transcriptome-wide level. However, the strategy of library generation affects its outcomes. Here, we present a modified ribosome-profiling protocol starting from yeast, human cells and vertebrate brain tissue. We use a DNA linker carrying four randomized positions at its 5′ end and a reverse-transcription (RT) primer with three randomized positions to reduce artifacts during library preparation. The use of seven randomized nucleotides allows to efficiently detect library-generation artifacts. We find that the effect of polymerase chain reaction (PCR) artifacts is relatively small for global analyses when sufficient input material is used. However, when input material is limiting, our strategy improves the sensitivity of gene-specific analyses. Furthermore, randomized nucleotides alleviate the skewed frequency of specific sequences at the 3′ end of ribosome-protected fragments (RPFs) likely resulting from ligase specificity. Finally, strategies that rely on dual ligation show a high degree of gene-coverage variation. Taken together, our approach helps to remedy two of the main problems associated with ribosome-profiling data. This will facilitate the analysis of translational dynamics and increase our understanding of the influence of RNA modifications on translation.

TADtool: visual parameter identification for TAD-calling algorithms

Summary: Eukaryotic genomes are hierarchically organized into topologically associating domains (TADs). The computational identification of these domains and their associated properties critically depends on the choice of suitable parameters of TAD-calling algorithms. To reduce the element of trial-and-error in parameter selection, we have developed TADtool: an interactive plot to find robust TAD-calling parameters with immediate visual feedback. TADtool allows the direct export of TADs called with a chosen set of parameters for two of the most common TAD calling algorithms: directionality and insulation index. It can be used as an intuitive, standalone application or as a Python package for maximum flexibility.

Availability and implementation: TADtool is available as a Python package from GitHub (https://github.com/vaquerizaslab/tadtool) or can be installed directly via PyPI, the Python package index (tadtool).

Contact:kai.kruse@mpi-muenster.mpg.de, jmv@mpi-muenster.mpg.de

Supplementary information:Supplementary data are available at Bioinformatics online.

Large-scale nuclear architecture and transcriptional control

Transcriptional regulation is one the most basic mechanisms for controlling gene expression. Over the past few years, much research has been devoted to understanding the interplay between transcription factors, histone modifications and associated enzymes required to achieve this control. However, it is becoming increasingly apparent that the three-dimensional conformation of chromatin in the interphase nucleus also plays a critical role in regulating transcription. Chromatin localisation in the nucleus is highly organised, and early studies described strong interactions between chromatin and sub-nuclear components. Single-gene studies have shed light on how chromosomal architecture affects gene expression. Lately, this has been complemented by whole-genome studies that have determined the global chromatin conformation of living cells in interphase. These studies have greatly expanded our understanding of nuclear architecture and its interplay with different physiological processes. Despite these advances, however, most of the mechanisms used to impose the three-dimensional chromatin structure remain unknown. Here, we summarise the different levels of chromatin organisation in the nucleus and discuss current efforts into characterising the mechanisms that govern it.

Response to Comments on "Drosophila Dosage Compensation Involves Enhanced Pol II Recruitment to Male X-Linked Promoters"

Ferrari et al. and Straub and Becker wrongly claim that an error in the computational analysis calls into question the conclusions of Conrad et al. All the available evidence, including the reanalyzed genomic data, show that the conclusions and the key message of the study remain unchanged: RNA polymerase II recruitment to male X-linked promoters is an important regulatory step during dosage compensation.

DNA-binding specificities of human transcription factors

Although the proteins that read the gene regulatory code, transcription factors (TFs), have been largely identified, it is not well known which sequences TFs can recognize. We have analyzed the sequence-specific binding of human TFs using high-throughput SELEX and ChIP sequencing. A total of 830 binding profiles were obtained, describing 239 distinctly different binding specificities. The models represent the majority of human TFs, approximately doubling the coverage compared to existing systematic studies. Our results reveal additional specificity determinants for a large number of factors for which a partial specificity was known, including a commonly observed A- or T-rich stretch that flanks the core motifs. Global analysis of the data revealed that homodimer orientation and spacing preferences, and base-stacking interactions, have a larger role in TF-DNA binding than previously appreciated. We further describe a binding model incorporating these features that is required to understand binding of TFs to DNA.

Drosophila dosage compensation involves enhanced Pol II recruitment to male X-linked promoters

Through hyperacetylation of histone H4 lysine 16 (H4K16), the male-specific lethal (MSL) complex in Drosophila approximately doubles transcription from the single male X chromosome in order to match X-linked expression in females and expression from diploid autosomes. By obtaining accurate measurements of RNA polymerase II (Pol II) occupancies and short promoter-proximal RNA production, we detected a consistent, genome-scale increase in Pol II activity at the promoters of male X-linked genes. Moreover, we found that enhanced Pol II recruitment to male X-linked promoters is largely dependent on the MSL complex. These observations provide insights into how global modulation of chromatin structure by histone acetylation contributes to the precise control of Pol II function.

m:Explorer: multinomial regression models reveal positive and negative regulators of longevity in yeast quiescence

We developed m:Explorer for identifying process-specific transcription factors (TFs) from multiple genome-wide sources, including transcriptome, DNA-binding and chromatin data. m:Explorer robustly outperforms similar techniques in finding cell cycle TFs in Saccharomyces cerevisiae. We predicted and experimentally tested regulators of quiescence (G0), a model of ageing, over a six-week time-course. We validated nine of top-12 predictions as novel G0 TFs, including Δmga2, Δcst6, Δbas1 with higher viability and G0-essential TFs Tup1, Swi3. Pathway analysis associates longevity to reduced growth, reprogrammed metabolism and cell wall remodeling. m:Explorer (http://biit.cs.ut.ee/mexplorer/) is instrumental in interrogating eukaryotic regulatory systems using heterogeneous data.

The NSL complex regulates housekeeping genes in Drosophila

MOF is the major histone H4 lysine 16-specific (H4K16) acetyltransferase in mammals and Drosophila. In flies, it is involved in the regulation of X-chromosomal and autosomal genes as part of the MSL and the NSL complexes, respectively. While the function of the MSL complex as a dosage compensation regulator is fairly well understood, the role of the NSL complex in gene regulation is still poorly characterized. Here we report a comprehensive ChIP-seq analysis of four NSL complex members (NSL1, NSL3, MBD-R2, and MCRS2) throughout the Drosophila melanogaster genome. Strikingly, the majority (85.5%) of NSL-bound genes are constitutively expressed across different cell types. We find that an increased abundance of the histone modifications H4K16ac, H3K4me2, H3K4me3, and H3K9ac in gene promoter regions is characteristic of NSL-targeted genes. Furthermore, we show that these genes have a well-defined nucleosome free region and broad transcription initiation patterns. Finally, by performing ChIP-seq analyses of RNA polymerase II (Pol II) in NSL1- and NSL3-depleted cells, we demonstrate that both NSL proteins are required for efficient recruitment of Pol II to NSL target gene promoters. The observed Pol II reduction coincides with compromised binding of TBP and TFIIB to target promoters, indicating that the NSL complex is required for optimal recruitment of the pre-initiation complex on target genes. Moreover, genes that undergo the most dramatic loss of Pol II upon NSL knockdowns tend to be enriched in DNA Replication-related Element (DRE). Taken together, our findings show that the MOF-containing NSL complex acts as a major regulator of housekeeping genes in flies by modulating initiation of Pol II transcription.

The MOF chromobarrel domain controls genome-wide H4K16 acetylation and spreading of the MSL complex

The histone H4 lysine 16 (H4K16)-specific acetyltransferase MOF is part of two distinct complexes involved in X chromosome dosage compensation and autosomal transcription regulation. Here we show that the MOF chromobarrel domain is essential for H4K16 acetylation throughout the Drosophila genome and is required for spreading of the male-specific lethal (MSL) complex on the X chromosome. The MOF chromobarrel domain directly interacts with nucleic acids and potentiates MOF's enzymatic activity after chromatin binding, making it a unique example of a chromo-like domain directly controlling acetylation activity in vivo. We also show that the Drosophila-specific N terminus of MOF has evolved to perform sex-specific functions. It modulates nucleosome binding and HAT activity and controls MSL complex assembly, thus regulating MOF function in dosage compensation. We propose that MOF has been especially tailored to achieve tight regulation of its enzymatic activity and enable its dual role on X and autosomes.

How do you find transcription factors? Computational approaches to compile and annotate repertoires of regulators for any genome

Transcription factors (TFs) play an important role in regulating gene expression. The availability of complete genome sequences and associated functional genomic data offer excellent opportunities to understand the transcriptional regulatory system of an entire organism. To do so, however, it is essential to compile a reliable dataset of regulatory components. Here, we review computational methods and publicly accessible resources that help identify TF-coding genes in prokaryotic and eukaryotic genomes. Since the regulatory functions of most TFs remain unknown, we also discuss approaches for combining diverse genomic datasets that will help elucidate their chromosomal organisation, expression, and evolutionary conservation. These analysis methods provide a solid foundation for further investigations of the transcriptional regulatory system.

SpeCond: a method to detect condition-specific gene expression

Transcriptomic studies routinely measure expression levels across numerous conditions. These datasets allow identification of genes that are specifically expressed in a small number of conditions. However, there are currently no statistically robust methods for identifying such genes. Here we present SpeCond, a method to detect condition-specific genes that outperforms alternative approaches. We apply the method to a dataset of 32 human tissues to determine 2,673 specifically expressed genes. An implementation of SpeCond is freely available as a Bioconductor package at http://www.bioconductor.org/packages/release/bioc/html/SpeCond.html.

The nonspecific lethal complex is a transcriptional regulator in Drosophila

Here, we report the biochemical characterization of the nonspecific lethal (NSL) complex (NSL1, NSL2, NSL3, MCRS2, MBD-R2, and WDS) that associates with the histone acetyltransferase MOF in both Drosophila and mammals. Chromatin immunoprecipitation-Seq analysis revealed association of NSL1 and MCRS2 with the promoter regions of more than 4000 target genes, 70% of these being actively transcribed. This binding is functional, as depletion of MCRS2, MBD-R2, and NSL3 severely affects gene expression genome wide. The NSL complex members bind to their target promoters independently of MOF. However, depletion of MCRS2 affects MOF recruitment to promoters. NSL complex stability is interdependent and relies mainly on the presence of NSL1 and MCRS2. Tethering of NSL3 to a heterologous promoter leads to robust transcription activation and is sensitive to the levels of NSL1, MCRS2, and MOF. Taken together, we conclude that the NSL complex acts as a major transcriptional regulator in Drosophila.

Multiplexed massively parallel SELEX for characterization of human transcription factor binding specificities

The genetic code-the binding specificity of all transfer-RNAs--defines how protein primary structure is determined by DNA sequence. DNA also dictates when and where proteins are expressed, and this information is encoded in a pattern of specific sequence motifs that are recognized by transcription factors. However, the DNA-binding specificity is only known for a small fraction of the approximately 1400 human transcription factors (TFs). We describe here a high-throughput method for analyzing transcription factor binding specificity that is based on systematic evolution of ligands by exponential enrichment (SELEX) and massively parallel sequencing. The method is optimized for analysis of large numbers of TFs in parallel through the use of affinity-tagged proteins, barcoded selection oligonucleotides, and multiplexed sequencing. Data are analyzed by a new bioinformatic platform that uses the hundreds of thousands of sequencing reads obtained to control the quality of the experiments and to generate binding motifs for the TFs. The described technology allows higher throughput and identification of much longer binding profiles than current microarray-based methods. In addition, as our method is based on proteins expressed in mammalian cells, it can also be used to characterize DNA-binding preferences of full-length proteins or proteins requiring post-translational modifications. We validate the method by determining binding specificities of 14 different classes of TFs and by confirming the specificities for NFATC1 and RFX3 using ChIP-seq. Our results reveal unexpected dimeric modes of binding for several factors that were thought to preferentially bind DNA as monomers.

Comprehensive reanalysis of transcription factor knockout expression data in Saccharomyces cerevisiae reveals many new targets

Transcription factor (TF) perturbation experiments give valuable insights into gene regulation. Genome-scale evidence from microarray measurements may be used to identify regulatory interactions between TFs and targets. Recently, Hu and colleagues published a comprehensive study covering 269 TF knockout mutants for the yeast Saccharomyces cerevisiae. However, the information that can be extracted from this valuable dataset is limited by the method employed to process the microarray data. Here, we present a reanalysis of the original data using improved statistical techniques freely available from the BioConductor project. We identify over 100,000 differentially expressed genes-nine times the total reported by Hu et al. We validate the biological significance of these genes by assessing their functions, the occurrence of upstream TF-binding sites, and the prevalence of protein-protein interactions. The reanalysed dataset outperforms the original across all measures, indicating that we have uncovered a vastly expanded list of relevant targets. In summary, this work presents a high-quality reanalysis that maximizes the information contained in the Hu et al. compendium. The dataset is available from ArrayExpress (accession: E-MTAB-109) and it will be invaluable to any scientist interested in the yeast transcriptional regulatory system.

Genome-wide analysis reveals MOF as a key regulator of dosage compensation and gene expression in Drosophila

Dosage compensation, mediated by the MSL complex, regulates X-chromosomal gene expression in Drosophila. Here we report that the histone H4 lysine 16 (H4K16) specific histone acetyltransferase MOF displays differential binding behavior depending on whether the target gene is located on the X chromosome versus the autosomes. More specifically, on the male X chromosome, where MSL1 and MSL3 are preferentially associated with the 3' end of dosage compensated genes, MOF displays a bimodal distribution binding to promoters and the 3' ends of genes. In contrast, on MSL1/MSL3 independent X-linked genes and autosomal genes in males and females, MOF binds primarily to promoters. Binding of MOF to autosomes is functional, as H4K16 acetylation and the transcription levels of a number of genes are affected upon MOF depletion. Therefore, MOF is not only involved in the onset of dosage compensation, but also acts as a regulator of gene expression in the Drosophila genome.

BABELOMICS: a systems biology perspective in the functional annotation of genome-scale experiments

We present a new version of Babelomics, a complete suite of web tools for functional analysis of genome-scale experiments, with new and improved tools. New functionally relevant terms have been included such as CisRed motifs or bioentities obtained by text-mining procedures. An improved indexing has considerably speeded up several of the modules. An improved version of the FatiScan method for studying the coordinate behaviour of groups of functionally related genes is presented, along with a similar tool, the Gene Set Enrichment Analysis. Babelomics is now more oriented to test systems biology inspired hypotheses. Babelomics can be found at .

PupaSuite: finding functional single nucleotide polymorphisms for large-scale genotyping purposes

We have developed a web tool, PupaSuite, for the selection of single nucleotide polymorphisms (SNPs) with potential phenotypic effect, specifically oriented to help in the design of large-scale genotyping projects. PupaSuite uses a collection of data on SNPs from heterogeneous sources and a large number of pre-calculated predictions to offer a flexible and intuitive interface for selecting an optimal set of SNPs. It improves the functionality of PupaSNP and PupasView programs and implements new facilities such as the analysis of user's data to derive haplotypes with functional information. A new estimator of putative effect of polymorphisms has been included that uses evolutionary information. Also SNPeffect database predictions have been included. The PupaSuite web interface is accessible through and through .

Next station in microarray data analysis: GEPAS

The Gene Expression Profile Analysis Suite (GEPAS) has been running for more than four years. During this time it has evolved to keep pace with the new interests and trends in the still changing world of microarray data analysis. GEPAS has been designed to provide an intuitive although powerful web-based interface that offers diverse analysis options from the early step of preprocessing (normalization of Affymetrix and two-colour microarray experiments and other preprocessing options), to the final step of the functional annotation of the experiment (using Gene Ontology, pathways, PubMed abstracts etc.), and include different possibilities for clustering, gene selection, class prediction and array-comparative genomic hybridization management. GEPAS is extensively used by researchers of many countries and its records indicate an average usage rate of 400 experiments per day. The web-based pipeline for microarray gene expression data, GEPAS, is available at http://www.gepas.org.

PupasView: a visual tool for selecting suitable SNPs, with putative pathological effect in genes, for genotyping purposes

We have developed a web tool, PupasView, for the selection of single nucleotide polymorphisms (SNPs) with potential phenotypic effect. PupasView constitutes an interactive environment in which functional information and population frequency data can be used as sequential filters over linkage disequilibrium parameters to obtain a final list of SNPs optimal for genotyping purposes. PupasView is the first resource that integrates phenotypic effects caused by SNPs at both the translational and the transcriptional level. PupasView retrieves SNPs that could affect conserved regions that the cellular machinery uses for the correct processing of genes (intron/exon boundaries or exonic splicing enhancers), predicted transcription factor binding sites and changes in amino acids in the proteins for which a putative pathological effect is calculated. The program uses the mapping of SNPs in the genome provided by Ensembl. PupasView will be of much help in studies of multifactorial disorders, where the use of functional SNPs will increase the sensitivity of the identification of the genes responsible for the disease. The PupasView web interface is accessible through http://pupasview.ochoa.fib.es and through http://www.pupasnp.org.

BABELOMICS: a suite of web tools for functional annotation and analysis of groups of genes in high-throughput experiments

We present Babelomics, a complete suite of web tools for the functional analysis of groups of genes in high-throughput experiments, which includes the use of information on Gene Ontology terms, interpro motifs, KEGG pathways, Swiss-Prot keywords, analysis of predicted transcription factor binding sites, chromosomal positions and presence in tissues with determined histological characteristics, through five integrated modules: FatiGO (fast assignment and transference of information), FatiWise, transcription factor association test, GenomeGO and tissues mining tool, respectively. Additionally, another module, FatiScan, provides a new procedure that integrates biological information in combination with experimental results in order to find groups of genes with modest but coordinate significant differential behaviour. FatiScan is highly sensitive and is capable of finding significant asymmetries in the distribution of genes of common function across a list of ordered genes even if these asymmetries were not extreme. The strong multiple-testing nature of the contrasts made by the tools is taken into account. All the tools are integrated in the gene expression analysis package GEPAS. Babelomics is the natural evolution of our tool FatiGO (which analysed almost 22 000 experiments during the last year) to include more sources on information and new modes of using it. Babelomics can be found at .

GEPAS, an experiment-oriented pipeline for the analysis of microarray gene expression data

The Gene Expression Profile Analysis Suite, GEPAS, has been running for more than three years. With >76,000 experiments analysed during the last year and a daily average of almost 300 analyses, GEPAS can be considered a well-established and widely used platform for gene expression microarray data analysis. GEPAS is oriented to the analysis of whole series of experiments. Its design and development have been driven by the demands of the biomedical community, probably the most active collective in the field of microarray users. Although clustering methods have obviously been implemented in GEPAS, our interest has focused more on methods for finding genes differentially expressed among distinct classes of experiments or correlated to diverse clinical outcomes, as well as on building predictors. There is also a great interest in CGH-arrays which fostered the development of the corresponding tool in GEPAS: InSilicoCGH. Much effort has been invested in GEPAS for developing and implementing efficient methods for functional annotation of experiments in the proper statistical framework. Thus, the popular FatiGO has expanded to a suite of programs for functional annotation of experiments, including information on transcription factor binding sites, chromosomal location and tissues. The web-based pipeline for microarray gene expression data, GEPAS, is available at http://www.gepas.org.

Highly specific and accurate selection of siRNAs for high-throughput functional assays

MOTIVATION

Small interfering RNA (siRNA) is widely used in functional genomics to silence genes by decreasing their expression to study the resulting phenotypes. The possibility of performing large-scale functional assays by gene silencing accentuates the necessity of a software capable of the high-throughput design of highly specific siRNA. The main objective sought was the design of a large number of siRNAs with appropriate thermodynamic properties and, especially, high specificity. Since all the available procedures require, to some extent, manual processing of the results to guarantee specific results, specificity constitutes to date, the major obstacle to the complete automation of all the steps necessary for the selection of optimal candidate siRNAs.

RESULT

Here, we present a program that for the first time completely automates the search for siRNAs. In SiDE, the most complete set of rules for the selection of siRNA candidates (including G+C content, nucleotides at determined positions, thermodynamic properties, propensity to form internal hairpins, etc.) is implemented and moreover, specificity is achieved by a conceptually new method. After selecting possible siRNA candidates with the optimal functional properties, putative unspecific matches, which can cause cross-hybridization, are checked in databases containing a unique entry for each gene. These truly non-redundant databases are constructed from the genome annotations (Ensembl). Also intron/exon boundaries, presence of polymorphisms (single nucleotide polymorphisms) specificity for either gene or transcript, and other features can be selected to be considered in the design of siRNAs.

AVAILABILITY

The program is available as a web server at http://side.bioinfo.cnio.es. The program was written under the GPL license.

CONTACT

jdopazo@cnio.es.

PupaSNP Finder: a web tool for finding SNPs with putative effect at transcriptional level

We have developed a web tool, PupaSNP Finder (PupaSNP for short), for high-throughput searching for single nucleotide polymorphisms (SNPs) with potential phenotypic effect. PupaSNP takes as its input lists of genes (or generates them from chromosomal coordinates) and retrieves SNPs that could affect the conserved regions that the cellular machinery uses for the correct processing of genes (intron/exon boundaries or exonic splicing enhancers), predicted transcription factor binding sites (TFBS) and changes in amino acids in the proteins. The program uses the mapping of SNPs in the genome provided by Ensembl. Additionally, user-defined SNPs (not yet mapped in the genome) can be easily provided to the program. Also, additional functional information from Gene Ontology, OMIM and homologies in other model organisms is provided. In contrast to other programs already available, which focus only on SNPs with possible effect in the protein, PupaSNP includes SNPs with possible transcriptional effect. PupaSNP will be of significant help in studies of multifactorial disorders, where the use of functional SNPs will increase the sensitivity of identification of the genes responsible for the disease. The PupaSNP web interface is accessible through http://pupasnp.bioinfo.cnio.es.

DNMAD: web-based diagnosis and normalization for microarray data

SUMMARY

We present a web server for Diagnosis and Normalization of MicroArray Data (DNMAD). DNMAD includes several common data transformations such as spatial and global robust local regression or multiple slide normalization, and allows for detecting several kinds of errors that result from the manipulation and the image analysis of the arrays. This tool offers a user-friendly interface, and is completely integrated within the Gene Expression Pattern Analysis Suite (GEPAS).

AVAILABILITY

The tool is accessible on-line at http://dnmad.bioinfo.cnio.es.

New challenges in gene expression data analysis and the extended GEPAS

Since the first papers published in the late nineties, including, for the first time, a comprehensive analysis of microarray data, the number of questions that have been addressed through this technique have both increased and diversified. Initially, interest focussed on genes coexpressing across sets of experimental conditions, implying, essentially, the use of clustering techniques. Recently, however, interest has focussed more on finding genes differentially expressed among distinct classes of experiments, or correlated to diverse clinical outcomes, as well as in building predictors. In addition to this, the availability of accurate genomic data and the recent implementation of CGH arrays has made mapping expression and genomic data on the chromosomes possible. There is also a clear demand for methods that allow the automatic transfer of biological information to the results of microarray experiments. Different initiatives, such as the Gene Ontology (GO) consortium, pathways databases, protein functional motifs, etc., provide curated annotations for genes. Whereas many resources on the web focus mainly on clustering methods, GEPAS has evolved to cope with the aforementioned new challenges that have recently arisen in the field of microarray data analysis. The web-based pipeline for microarray gene expression data, GEPAS, is available at http://gepas.bioinfo.cnio.es.

GEPAS: A web-based resource for microarray gene expression data analysis

We present a web-based pipeline for microarray gene expression profile analysis, GEPAS, which stands for Gene Expression Profile Analysis Suite (http://gepas.bioinfo.cnio.es). GEPAS is composed of different interconnected modules which include tools for data pre-processing, two-conditions comparison, unsupervised and supervised clustering (which include some of the most popular methods as well as home made algorithms) and several tests for differential gene expression among different classes, continuous variables or survival analysis. A multiple purpose tool for data mining, based on Gene Ontology, is also linked to the tools, which constitutes a very convenient way of analysing clustering results. On-line tutorials are available from our main web server (http://bioinfo.cnio.es).