Search

SUMMARY PARAGRAPHDisease-causing mutations in genes encoding transcription factors (TFs) are a recurrent finding in hematopoietic malignancies and might involve key regulators of lineage adherence and cellular differentiation1–3. Such mutations can affect TF-interactions with their cognate DNA-binding motifs4, 5. Whether and how TF-mutations impact upon the nature of binding to TF composite elements (CE) and influence their interaction with other TFs is unclear. Here, we report a new mechanism of TF alteration in human lymphomas with perturbed B cell identity. It is caused by a recurrent somatic missense mutation c.295T>C (p.Cys99Arg; p.C99R) targeting the center of the DNA-binding domain of Interferon Regulatory Factor 4 (IRF4), a key TF in immune cell-differentiation and -activation6, 7. IRF4-C99R fundamentally alters IRF4 DNA-binding, with loss-of-binding to canonical IRF motifs and neomorphic gain-of-binding to canonical and non-canonical IRF composite elements (CEs). Furthermore, IRF4-C99R thoroughly modifies IRF4 function, by blocking IRF4-dependent plasma cell induction, and up-regulating disease-specific genes in a non-canonical Activator Protein-1 (AP-1)-IRF-CE (AICE)-dependent manner. Our data explain how a single arginine mutation creates a complex switch of TF specificity and gene regulation. These data open the possibility of designing specific inhibitors to block the neomorphic, disease-causing DNA-binding activities of a mutant transcription factor.

Comprehensive genomic profiling using cancer gene panels has been shown to improve treatment options for a variety of cancer types. However, genomic aberrations detected via such gene panels do not necessarily serve as strong predictors of drug sensitivity. In this study, using pharmacogenomics datasets of cell lines, patient-derived xenografts, and ex vivo treated fresh tumor specimens, we demonstrate that utilizing the transcriptome on top of gene panel features substantially improves drug response prediction performance in cancer.

Comprehensive genomic profiling using cancer gene panels has been shown to improve treatment options for a variety of cancer types. However, genomic aberrations detected via such gene panels don’t necessarily serve as strong predictors of drug sensitivity. In this study, using pharmacogenomics datasets of cell lines, patient-derived xenografts, and ex-vivo treated fresh tumor specimens, we demonstrate that utilizing the transcriptome on top of gene panel features substantially improves drug response prediction performance in cancer.

The use of RNA sequencing from wastewater samples is a valuable way for estimating infection dynamics and circulating lineages of SARS-CoV-2. This approach is independent from testing individuals and can therefore become the key tool to monitor this and potentially other viruses. However, it is equally important to develop easily accessible and scalable tools which can highlight critical changes in infection rates and dynamics over time across different locations given sequencing data from wastewater. Here, we provide an analysis of lineage dynamics in Berlin and New York City using wastewater sequencing and present PiGx SARS-CoV-2, a highly reproducible computational analysis pipeline with comprehensive reports. This end-to-end pipeline includes all steps from raw data to shareable reports, additional taxonomic analysis, deconvolution and geospatial time series analyses. Using simulated datasets (in silico generated and spiked-in samples) we could demonstrate the accuracy of our pipeline calculating proportions of Variants of Concern (VOC) from environmental as well as pre-mixed samples (spiked-in). By applying our pipeline on a dataset of wastewater samples from Berlin between February 2021 and January 2022, we could reconstruct the emergence of B.1.1.7(alpha) in February/March 2021 and the replacement dynamics from B.1.617.2 (delta) to BA.1 and BA.2 (omicron) during the winter of 2021/2022. Using data from very-short-reads generated in an industrial scale setting, we could see even higher accuracy in our deconvolution. Lastly, using a targeted sequencing dataset from New York City (receptor-binding-domain (RBD) only), we could reproduce the results recovering the proportions of the so-called cryptic lineages shown in the original study. Overall our study provides an in-depth analysis reconstructing virus lineage dynamics from wastewater. While applying our tool on a wide range of different datasets (from different types of wastewater sample locations and sequenced with different methods), we show that PiGx SARS-CoV-2 can be used to identify new mutations and detect any emerging new lineages in a highly automated and scalable way. Our approach can support efforts to establish continuous monitoring and early-warning projects for detecting SARS-CoV-2 or any other pathogen.

The heptarepeats of the C-terminal domain (CTD) of RNA polymerase II (Pol II) are extensively modified throughout the transcription cycle. The CTD coordinates RNA synthesis and processing by recruiting transcription regulation factors as well as RNA capping, splicing and 3’end processing factors. The SPOC domain of PHF3 was recently identified as a new CTD reader domain specifically binding to phosphorylated serine-2 residues in adjacent CTD repeats. Here, we establish the SPOC domains of the human proteins DIDO, SHARP and RBM15 as phosphoserine binding modules that can act as CTD readers but also recognize other phosphorylated binding partners. We report the crystal structure of SHARP (SPEN) SPOC-CTD and identify the molecular determinants for its specific binding to phosphorylated serine-5. PHF3 and DIDO SPOC domains preferentially interact with the Pol II elongation complex, while RBM15 and SHARP SPOC domains engage with the m6A writer and reader proteins. Our findings establish the SPOC domain as a major interface between the transcription machinery and regulators of transcription and co-transcriptional processes.

The use of RNA sequencing from wastewater samples is a valuable way for estimating infection dynamics and circulating lineages of SARS-CoV-2. This approach is independent from testing individuals and can therefore become the key tool to monitor this and potentially other viruses. However, it is equally important to develop easily accessible and scalable tools which can highlight critical changes in infection rates and dynamics over time across different locations given sequencing data from wastewater. Here, we provide an analysis of lineage dynamics in Berlin and New York City using wastewater sequencing and present PiGx SARS-CoV-2, a highly reproducible computational analysis pipeline with comprehensive reports. This end-to-end pipeline includes all steps from raw data to shareable reports, additional taxonomic analysis, deconvolution and geospatial time series analyses. Using simulated datasets (in silico generated and spiked-in samples) we could demonstrate the accuracy of our pipeline calculating proportions of Variants of Concern (VOC) from environmental as well as pre-mixed samples (spiked-in). By applying our pipeline on a dataset of wastewater samples from Berlin between February 2021 and January 2022, we could reconstruct the emergence of B.1.1.7(alpha) in February/March 2021 and the replacement dynamics from B.1.617.2 (delta) to BA.1 and BA.2 (omicron) during the winter of 2021/2022. Using data from very-short-reads generated in an industrial scale setting, we could see even higher accuracy in our deconvolution. Lastly, using a targeted sequencing dataset from New York City (receptor-binding-domain (RBD) only), we could reproduce the results recovering the proportions of the so-called cryptic lineages shown in the original study. Overall our study provides an in-depth analysis reconstructing virus lineage dynamics from wastewater, and that our tool can be used to identify new mutations and to detect any emerging new lineages with different amplification and sequencing methods. Our approach can support efforts to establish continuous monitoring and early-warning projects for detecting SARS-CoV-2 or any other pathogen.

The C-terminal domain (CTD) of the largest subunit of RNA polymerase II (Pol II) is a regulatory hub for transcription and RNA processing. Here, we identify PHD-finger protein 3 (PHF3) as a regulator of transcription and mRNA stability that docks onto Pol II CTD through its SPOC domain. We characterize SPOC as a CTD reader domain that preferentially binds two phosphorylated Serine-2 marks in adjacent CTD repeats. PHF3 drives liquid-liquid phase separation of phosphorylated Pol II, colocalizes with Pol II clusters and tracks with Pol II across the length of genes. PHF3 knock-out or SPOC deletion in human cells results in increased Pol II stalling, reduced elongation rate and an increase in mRNA stability, with marked derepression of neuronal genes. Key neuronal genes are aberrantly expressed in Phf3 knock-out mouse embryonic stem cells, resulting in impaired neuronal differentiation. Our data suggest that PHF3 acts as a prominent effector of neuronal gene regulation by bridging transcription with mRNA decay., Here the authors identify PHF3 SPOC domain as a reader of the phosphorylated RNA polymerase II (Pol II) C-terminal domain. They show that PHF3 clusters with Pol II complexes in cells, drives phase separation of Pol II in vitro, and regulates neuronal gene expression and neuronal differentiation.

Background Acute inflammation leads to increased myelopoiesis in the bone marrow and recruitment of immune cells to respective sites of injury. Yet, in certain cases of high demand such as acute tissue damage as in acute myocardial infarction (AMI), innate immune cells can be additionally recruited from lymphoid organs such as the spleen. Monocytes are important in both onset and resolution of inflammation comprising a huge repertoire of coding and non-coding RNA molecules regulating a fast-acting immune response. Long non-coding RNAs (lncRNA) emerge as novel regulatory cell-specific molecules, which are able to influence cell physiology in numerous ways. Purpose Delineate monocyte subpopulations from patients with AMI compared to healthy controls based on their non-coding transcriptome, different isoform usage and associated source of origin. Methods and results Next generation bulk RNA sequencing has been performed from human monocytes after pre-sorting for CD14 expression and subsequent separation into three major subpopulations based on their surface markers CD14 and CD16 using a FACS Aria II cell sorter. Non-coding RNA expression shows clear differences between cells from healthy controls versus patients with AMI and among subpopulations. Intensive bioinformatics analysis revealed numerous annotated lncRNA, antisense, pseudogene and circular RNAs displaying significant differences in their expression profiles. Some lncRNA show positive or negative correlations with their respective coding antisense gene or their convergent or divergent sense neighbouring gene, such as linc02207 vs MCTP2, SAP30-AS vs SAP30, BID-AS vs BID, HLX-AS vs HLX or TNFAIP3-AS vs TNFAIP3. We also identified the existence of different isoforms of certain lncRNA and a distinct preference thereof, such as linc02207- 217 isoform being most abundant. Real-time PCR was applied for validation of lncRNA expression profile and their nucleocytoplasmic distribution observing preferential nuclear expression for most lncRNA. In vitro assays for silencing and overexpression approaches as well as treatment with inflammatory stimuli and in silico bioinformatics analysis will help to unravel their functionality. Additionally mapping our bulk RNA data to deconvoluted single cell RNAseq data from publicly available human cell atlas (HCA), is applied to additionally unravel the origin of particular subpopulations within the AMI monocytes, associating their origin to either spleen or bone marrow. Conclusions High read depth bulk RNA sequencing revealed a human monocyte-specific long non-coding RNA transcriptome differentially expressed in three predefined subpopulations and disease stages. Distinct expression of lncRNA and associated coding genes as well as distinct isoform preference implies specific functionalities. Studying these molecules and their interaction partners will unravel functional participants in the development of therapeutic approaches for cardiovascular diseases. Funding Acknowledgement Type of funding source: Public grant(s) – National budget only. Main funding source(s): DZHK, SNF

The coronavirus disease 2019 (COVID-19) pandemic, caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is an ongoing global health threat with more than two million infected people since its emergence in late 2019. Detailed knowledge of the molecular biology of the infection is indispensable for understanding of the viral replication, host responses, and disease progression. We provide gene expression profiles of SARS-CoV and SARS-CoV-2 infections in three human cell lines (H1299, Caco-2 and Calu-3 cells), using bulk and single-cell transcriptomics. Small RNA profiling showed strong expression of the immunity and inflammation-associated microRNA miRNA-155 upon infection with both viruses. SARS-CoV-2 elicited approximately two-fold higher stimulation of the interferon response compared to SARS-CoV in the permissive human epithelial cell line Calu-3, and induction of cytokines such as CXCL10 or IL6. Single cell RNA sequencing data showed that canonical interferon stimulated genes such as IFIT2 or OAS2 were broadly induced, whereas interferon beta (IFNB1) and lambda (IFNL1-4) were expressed only in a subset of infected cells. In addition, temporal resolution of transcriptional responses suggested interferon regulatory factors (IRFs) activities precede that of nuclear factor-κB (NF-κB). Lastly, we identified heat shock protein 90 (HSP90) as a protein relevant for the infection. Inhibition of the HSP90 charperone activity by Tanespimycin/17-N-allylamino-17-demethoxygeldanamycin (17-AAG) resulted in a reduction of viral replication, and of TNF and IL1B mRNA levels. In summary, our study established in vitro cell culture models to study SARS-CoV-2 infection and identified HSP90 protein as potential drug target for therapeutic intervention of SARS-CoV-2 infection.

Syncytial skeletal muscle cells contain hundreds of nuclei in a shared cytoplasm. We investigated nuclear heterogeneity and transcriptional dynamics in the uninjured and regenerating muscle using single-nucleus RNA-sequencing (snRNAseq) of isolated nuclei from muscle fibers. This revealed distinct nuclear subtypes unrelated to fiber type diversity, previously unknown subtypes as well as the expected ones at the neuromuscular and myotendinous junctions. In fibers of the Mdx dystrophy mouse model, distinct subtypes emerged, among them nuclei expressing a repair signature that were also abundant in the muscle of dystrophy patients, and a nuclear population associated with necrotic fibers. Finally, modifications of our approach revealed the compartmentalization in the rare and specialized muscle spindle. Our data identifies nuclear compartments of the myofiber and defines a molecular roadmap for their functional analyses; the data can be freely explored on the MyoExplorer server (https://shiny.mdc-berlin.de/MyoExplorer/)., The transcriptional programs of nuclei in the muscle syncytium were assumed to be homogenous except at the neuromuscular and myotendinous junctions. Here, using single-nucleus transcriptomics, the authors reveal a previously unrecognized diversity and dynamics of myonuclear transcriptional programs.

Neurons and oligodendrocytes are terminally differentiated cells that perform highly specialized functions, which depend on cascades of gene activation and repression to retain homeostatic control over a lifespan. Gene expression is regulated by three-dimensional (3D) genome organisation, from local levels of chromatin compaction to the organisation of topological domains and chromosome compartments. Whereas our understanding of 3D genome architecture has vastly increased in the past decade, it remains difficult to study specialized cells in their native environment without disturbing their activity. To develop the application of Genome Architecture Mapping (GAM) in small numbers of specialized cells in complex tissues, we combined GAM with immunoselection. We applied immunoGAM to map the genome architecture of specific cell populations in the juvenile/adult mouse brain: dopaminergic neurons (DNs) from the midbrain, pyramidal glutamatergic neurons (PGNs) from the hippocampus, and oligodendrocyte lineage cells (OLGs) from the cortex. We integrate 3D genome organisation with single-cell transcriptomics data, and find specific chromatin structures that relate with cell-type specific patterns of gene expression. We discover abundant changes in compartment organisation, especially a strengthening of heterochromatin compartments which establish strong contacts spanning tens of megabases, especially in brain cells. These compartments contain olfactory and taste receptor genes, which are de-repressed in a subpopulation of PGNs with molecular signatures of long-term potentiation (LTP). We also show extensive reorganisation of topological domains where activation of neuronal or oligodendrocyte genes coincides with formation of new TAD borders. Finally, we discover loss of TAD organisation, or ‘TAD melting’, at long (>1Mb) neuronal genes when they are most highly expressed. Our work shows that the 3D organisation of the genome is highly cell-type specific in terminally differentiated cells of the brain, and essential to better understand brain-specific mechanisms of gene regulation.

SUMMARYThe C-terminal domain (CTD) of the largest subunit of RNA polymerase II (Pol II) is a regulatory hub for transcription and RNA processing. Here, we identify PHD-finger protein 3 (PHF3) as a new CTD-binding factor that negatively regulates transcription and mRNA stability. The PHF3 SPOC domain preferentially binds to CTD repeats phosphorylated on Serine-2 and PHF3 tracks with Pol II across the length of genes. PHF3 competes with TFIIS for Pol II binding through its TFIIS-like domain (TLD), thus inhibiting TFIIS-mediated rescue of backtracked Pol II. PHF3 knock-out or PHF3 SPOC deletion in human cells result in gene upregulation and a global increase in mRNA stability, with marked derepression of neuronal genes. Key neuronal genes are aberrantly expressed in Phf3 knock-out mouse embryonic stem cells, resulting in impaired neuronal differentiation. Our data suggest that PHF3 is a prominent effector of neuronal gene regulation at the interface of transcription elongation and mRNA decay.

Cancer is a heterogeneous collection of diseases traditionally classified by the tissue of origin. The diversity of the molecular profiles of cancers has a big impact on the way patients are diagnosed and treated, how they respond to their prescribed treatments, the duration of survival after diagnosis, and factors such as remission, recurrence, or spread (metastasis) of the disease. While such diagnostic and prognostic outcomes are potentially predictable by taking a closer look into the changes of the genome, epigenome, transcriptome, proteome, and various other omics platforms, the contemporary cancer treatments still predominantly don't make the best use of such multi-omics profiling of patient samples. Therefore, multi-omics profiling of cancers holds great potential to define a molecularly coherent subtype definition of cancers in order to achieve the eventual goal of matching the best possible treatment to the subgroup of patients. However, the current subtypes from consortiums such as TCGA have been defined by heterogeneous methods and molecular markers by different teams. A subset of these studies have not attempted to characterize molecular subtypes, but rather taken histopathologically defined subtypes as the gold standard and tried to characterize molecular features of these subtypes. Here we evaluate TCGA cancer subtypes based on the molecular profile coherence score. This novel metric combines survival statistics, pathways information, tumor purity estimates, and mutational signatures. We expect that subtypes that are patient subgroups should display molecular signature homogeneity. We evaluate TCGA subtypes from 21 cancers using these criteria and compare the subtypes with our own definition using multi-omics data in a deep learning framework. We have refined the several subtypes from multiple cancers towards more molecularly coherent patient subgroups. Citation Format: Altuna Akalin, Bora Uyar, Jonathan Ronen, Vedran Franke. Redefining cancer subtypes using multi-omics and deep learning [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 196.

Retrotransposons are “copy-and-paste” insertional mutagens that substantially contribute to mammalian genome content. Retrotransposons often carry long terminal repeats (LTRs) for retrovirus-like reverse transcription and integration into the genome. We report an extraordinary impact of a group of LTRs from the mammalian endogenous retrovirus-related ERVL retrotransposon class on gene expression in the germline and beyond. In mouse, we identified more than 800 LTRs from ORR1, MT, MT2, and MLT families, which resemble mobile gene-remodeling platforms that supply promoters and first exons. The LTR-mediated gene remodeling also extends to hamster, human, and bovine oocytes. The LTRs function in a stage-specific manner during the oocyte-to-embryo transition by activating transcription, altering protein-coding sequences, producing noncoding RNAs, and even supporting evolution of new protein-coding genes. These functions result, for example, in recycling processed pseudogenes into mRNAs or lncRNAs with regulatory roles. The functional potential of the studied LTRs is even higher, because we show that dormant LTR promoter activity can rescue loss of an essential upstream promoter. We also report a novel protein-coding gene evolution—D6Ertd527e—in which an MT LTR provided a promoter and the 5′ exon with a functional start codon while the bulk of the protein-coding sequence evolved through a CAG repeat expansion. Altogether, ERVL LTRs provide molecular mechanisms for stochastically scanning, rewiring, and recycling genetic information on an extraordinary scale. ERVL LTRs thus offer means for a comprehensive survey of the genome's expression potential, tightly intertwining with gene expression and evolution in the germline.

Persistent NF-κB activation is a hallmark of the malignant Hodgkin/Reed-Sternberg (HRS) cells in classical Hodgkin lymphoma (cHL). Genomic lesions, Epstein-Barr virus infection, soluble factors, and tumor–microenvironment interactions contribute to this activation. Here, in an unbiased approach to identify the cHL cell-secreted key factors for NF-κB activation, we have dissected the secretome of cultured cHL cells by chromatography and subsequent mass spectrometry. We identified lymphotoxin-α (LTA) as the causative factor for autocrine and paracrine activation of canonical and noncanonical NF-κB in cHL cell lines. In addition to inducing NF-κB, LTA promotes JAK2/STAT6 signaling. LTA and its receptor TNFRSF14 are transcriptionally activated by noncanonical NF-κB, creating a continuous feedback loop. Furthermore, LTA shapes the expression of cytokines, receptors, immune checkpoint ligands and adhesion molecules, including CSF2, CD40, PD-L1/PD-L2, and VCAM1. Comparison with single-cell gene-activity profiles of human hematopoietic cells showed that LTA induces genes restricted to the lymphoid lineage, as well as those largely restricted to the myeloid lineage. Thus, LTA sustains autocrine NF-κB activation, impacts activation of several signaling pathways, and drives expression of genes essential for microenvironmental interactions and lineage ambiguity. These data provide a robust rationale for targeting LTA as a treatment strategy for cHL patients.

In bioinformatics, as well as other computationally-intensive research fields, there is a need for workflows that can reliably produce consistent output, independent of the software environment or configuration settings of the machine on which they are executed. Indeed, this is essential for controlled comparison between different observations or for the wider dissemination of workflows. Providing this type of reproducibility, however, is often complicated by the need to accommodate the myriad dependencies included in a larger body of software, each of which generally come in various versions. Moreover, in many fields (bioinformatics being a prime example), these versions are subject to continual change due to rapidly evolving technologies, further complicating problems related to reproducibility. Here, we propose a principled approach for building analysis pipelines and managing their dependencies. As a case study to demonstrate the utility of our approach, we present a set of highly reproducible pipelines for the analysis of RNA-seq, ChIP-seq, Bisulfite-seq, and single-cell RNA-seq. All pipelines process raw experimental data, and generate reports containing publication-ready plots and figures, with interactive report elements and standard observables. Users may install these highly reproducible packages and apply them to their own datasets without any special computational expertise beyond the use of the command line. We hope such a toolkit will provide immediate benefit to laboratory workers wishing to process their own data sets or bioinformaticians seeking to automate all, or parts of, their analyses. In the long term, we hope our approach to reproducibility will serve as a blueprint for reproducible workflows in other areas. Our pipelines, along with their corresponding documentation and sample reports, are available at http://bioinformatics.mdc-berlin.de/pigx

In bioinformatics, as well as other computationally intensive research fields, there is a need for workflows that can reliably produce consistent output, from known sources, independent of the software environment or configuration settings of the machine on which they are executed. Indeed, this is essential for controlled comparison between different observations and for the wider dissemination of workflows. However, providing this type of reproducibility and traceability is often complicated by the need to accommodate the myriad dependencies included in a larger body of software, each of which generally comes in various versions. Moreover, in many fields (bioinformatics being a prime example), these versions are subject to continual change due to rapidly evolving technologies, further complicating problems related to reproducibility. Here, we propose a principled approach for building analysis pipelines and managing their dependencies with GNU Guix. As a case study to demonstrate the utility of our approach, we present a set of highly reproducible pipelines called PiGx for the analysis of RNA sequencing, chromatin immunoprecipitation sequencing, bisulfite-treated DNA sequencing, and single-cell resolution RNA sequencing. All pipelines process raw experimental data and generate reports containing publication-ready plots and figures, with interactive report elements and standard observables. Users may install these highly reproducible packages and apply them to their own datasets without any special computational expertise beyond the use of the command line. We hope such a toolkit will provide immediate benefit to laboratory workers wishing to process their own datasets or bioinformaticians seeking to automate all, or parts of, their analyses. In the long term, we hope our approach to reproducibility will serve as a blueprint for reproducible workflows in other areas. Our pipelines, along with their corresponding documentation and sample reports, are available at http://bioinformatics.mdc-berlin.de/pigx

Initiation of zygotic transcription in mammals is poorly understood. In mice, zygotic transcription is first detected shortly after pronucleus formation in 1-cell embryos, but the identity of the transcribed loci and mechanisms regulating their expression are not known. Using total RNA-Seq, we have found that transcription in 1-cell embryos is highly promiscuous, such that intergenic regions are extensively expressed and thousands of genes are transcribed at comparably low levels. Striking is that transcription can occur in the absence of defined core-promoter elements. Furthermore, accumulation of translatable zygotic mRNAs is minimal in 1-cell embryos because of inefficient splicing and 3′ processing of nascent transcripts. These findings provide novel insights into regulation of gene expression in 1-cell mouse embryos that may confer a protective mechanism against precocious gene expression that is the product of a relaxed chromatin structure present in 1-cell embryos. The results also suggest that the first zygotic transcription itself is an active component of chromatin remodeling in 1-cell embryos.

Summary: Biological insights can be obtained through computational integration of genomics data sets consisting of diverse types of information. The integration is often hampered by a large variety of existing file formats, often containing similar information, and the necessity to use complicated tools to achieve the desired results. We have built an R package, genomation, to expedite the extraction of biological information from high throughput data. The package works with a variety of genomic interval file types and enables easy summarization and annotation of high throughput data sets with given genomic annotations. Availability and implementation: The software is currently distributed under MIT artistic license and freely available at http://bioinformatics.mdc-berlin.de/genomation, and through the Bioconductor framework. Contact: dirk.schubeler@fmi.ch, chm2042@med.cornell.edu, altuna.akalin@fmi.ch, or aakalin@gmail.com

High-occupancy target (HOT) regions are the segments of the genome with unusually high number of transcription factor binding sites. These regions are observed in multiple species and thought to have biological importance due to high transcription factor occupancy. Furthermore, they coincide with house-keeping gene promoters and the associated genes are stably expressed across multiple cell types. Despite these features, HOT regions are solemnly defined using ChIP-seq experiments and shown to lack canonical motifs for transcription factors that are thought to be bound there. Although, ChIP-seq experiments are the golden standard for finding genome-wide binding sites of a protein, they are not noise free. Here, we show that HOT regions are likely to be ChIP-seq artifacts and they are similar to previously proposed “hyper-ChIPable” regions. Using ChIP-seq data sets for knocked-out transcription factors, we demonstrate presence of false positive signals on HOT regions. We observe sequence characteristics and genomic features that are discriminatory of HOT regions, such as GC/CpG-rich k-mers and enrichment of RNA-DNA hybrids (R-loops) and DNA tertiary structures (G-quadruplex DNA). The artificial ChIP-seq enrichment on HOT regions could be associated to these discriminatory features. Furthermore, we propose strategies to deal with such artifacts for the future ChIP-seq studies.

The oocyte-to-embryo transition (OET) transforms a differentiated gamete into pluripotent blastomeres. The accompanying maternal-zygotic RNA exchange involves remodeling of the long non-coding RNA (lncRNA) pool. Here, we used next generation sequencing and de novo transcript assembly to define the core population of 1,600 lncRNAs expressed during the OET (lncRNAs). Relative to mRNAs, OET lncRNAs were less expressed and had shorter transcripts, mainly due to fewer exons and shorter 5′ terminal exons. Approximately half of OET lncRNA promoters originated in retrotransposons suggesting their recent emergence. Except for a small group of ubiquitous lncRNAs, maternal and zygotic lncRNAs formed two distinct populations. The bulk of maternal lncRNAs was degraded before the zygotic genome activation. Interestingly, maternal lncRNAs seemed to undergo cytoplasmic polyadenylation observed for dormant mRNAs. We also identified lncRNAs giving rise to trans-acting short interfering RNAs, which represent a novel lncRNA category. Altogether, we defined the core OET lncRNA transcriptome and characterized its remodeling during early development. Our results are consistent with the notion that rapidly evolving lncRNAs constitute signatures of cells-of-origin while a minority plays an active role in control of gene expression across OET. Our data presented here provide an excellent source for further OET lncRNA studies.

In mouse, the oocyte-to-embryo transition entails converting a highly differentiated oocyte to totipotent blastomeres. This transition is driven by degradation of maternal mRNAs, which results in loss of oocyte identity, and reprogramming of gene expression during the course of zygotic gene activation, which occurs primarily during the two-cell stage and confers blastomere totipotency. Full-grown oocytes are transcriptionally quiescent and mRNAs are remarkably stable in oocytes due to the RNA-binding protein MSY2, which stabilizes mRNAs, and low activity of the 5' and 3' RNA degradation machinery. Oocyte maturation initiates a transition from mRNA stability to instability due to phosphorylation of MSY2, which makes mRNAs more susceptible to the RNA degradation machinery, and recruitment of dormant maternal mRNAs that encode for critical components of the 5' and 3' RNA degradation machinery. Small RNAs (miRNA, siRNA, and piRNA) play little, if any, role in mRNA degradation that occurs during maturation. Many mRNAs are totally degraded but a substantial fraction is only partially degraded, their degradation completed by the end of the two-cell stage. Genome activation initiates during the one-cell stage, is promiscuous, low level, and genome wide (and includes both inter- and intragenic regions) and produces transcripts that are inefficiently spliced and polyadenylated. The major wave of genome activation in two-cell embryos involves expression of thousands of new genes. This unique pattern of gene expression is the product of maternal mRNAs recruited during maturation that encode for transcription factors and chromatin remodelers, as well as dramatic changes in chromatin structure due to incorporation of histone variants and modified histones.

The oocyte-to-embryo transition entails genome activation and a dramatic reprogramming of gene expression that is required for continued development. Superimposed on genome activation and reprogramming is development of a transcriptionally repressive state at the level of chromatin structure. Inducing global histone hyperacetylation relieves this repression and histone deacetylases 1 and 2 (HDAC1 and HDAC2) are involved in establishing the repressive state. Because SIN3A is an HDAC1/2-containing complex, we investigated whether it is involved in reprogramming gene expression during the course of genome activation. We find that Sin3a mRNA is recruited during maturation and that inhibiting its recruitment not only inhibits development beyond the 2-cell stage but also compromises the fidelity of reprogramming gene expression. The SIN3A that is synthesized during oocyte maturation reaches a maximum level in the mid-1-cell embryo and is essentially absent by the mid-2-cell stage. Overexpressing SIN3A in 1-cell embryos has no obvious effect on pre- and postimplantation development. These results provide a mechanism by which reprogramming can occur using a maternally inherited transcription machinery, namely, recruitment of mRNAs encoding transcription factors and chromatin remodelers, such as SIN3A.

In mouse, the oocyte-to-embryo transition entails converting a highly differentiated oocyte to totipotent blastomeres. This transition is driven by degradation of maternal mRNAs, which results in loss of oocyte identity, and reprogramming of gene expression during the course of zygotic gene activation, which occurs primarily during the two-cell stage and confers blastomere totipotency. Full-grown oocytes are transcriptionally quiescent and mRNAs are remarkably stable in oocytes due to the RNA-binding protein MSY2, which stabilizes mRNAs, and low activity of the 5' and 3' RNA degradation machinery. Oocyte maturation initiates a transition from mRNA stability to instability due to phosphorylation of MSY2, which makes mRNAs more susceptible to the RNA degradation machinery, and recruitment of dormant maternal mRNAs that encode for critical components of the 5' and 3' RNA degradation machinery. Small RNAs (miRNA, siRNA, and piRNA) play little, if any, role in mRNA degradation that occurs during maturation. Many mRNAs are totally degraded but a substantial fraction is only partially degraded, their degradation completed by the end of the two-cell stage. Genome activation initiates during the one-cell stage, is promiscuous, low level, and genome wide (and includes both inter- and intragenic regions) and produces transcripts that are inefficiently spliced and polyadenylated. The major wave of genome activation in two-cell embryos involves expression of thousands of new genes. This unique pattern of gene expression is the product of maternal mRNAs recruited during maturation that encode for transcription factors and chromatin remodelers, as well as dramatic changes in chromatin structure due to incorporation of histone variants and modified histones.

Aminoacyl-tRNA synthetases (aaRS) are essential enzymes catalyzing the formation of aminoacyl-tRNAs, the immediate precursors for encoded peptides in ribosomal protein synthesis. Previous studies have suggested a link between tRNA aminoacylation and high-molecular-weight cellular complexes such as the cytoskeleton or ribosomes. However, the structural basis of these interactions and potential mechanistic implications are not well understood. To biochemically characterize these interactions we have used a system of two interacting archaeal aaRSs: an atypical methanogenic-type seryl-tRNA synthetase and an archaeal ArgRS. More specifically, we have shown by thermophoresis and surface plasmon resonance that these two aaRSs bind to the large ribosomal subunit with micromolar affinities. We have identified the L7/L12 stalk and the proteins located near the stalk base as the main sites for aaRS binding. Finally, we have performed a bioinformatics analysis of synonymous codons in the Methanothermobacter thermautotrophicus genome that supports a mechanism in which the deacylated tRNAs may be recharged by aaRSs bound to the ribosome and reused at the next occurrence of a codon encoding the same amino acid. These results suggest a mechanism of tRNA recycling in which aaRSs associate with the L7/L12 stalk region to recapture the tRNAs released from the preceding ribosome in polysomes., Nucleic Acids Research, 42 (8), ISSN:1362-4962, ISSN:0301-5610

In mammals, a single Dicer participates in biogenesis of small RNAs in microRNA (miRNA) and RNAi pathways. In mice, endogenous RNAi is highly active in oocytes, but not in somatic cells, which we ascribe here to an oocyte-specific Dicer isoform (Dicer(O)). Dicer(O) lacks the N-terminal DExD helicase domain and has higher cleavage activity than the full-length Dicer in somatic cells (Dicer(S)). Unlike Dicer(S), Dicer(O) efficiently produces small RNAs from long double-stranded (dsRNA) substrates. Expression of the Dicer(O) isoform is driven by an intronic MT-C retrotransposon promoter, deletion of which causes loss of Dicer(O) and female sterility. Oocytes from females lacking the MT-C element show meiotic spindle defects and increased levels of endogenous small interfering RNA (endo-siRNA) targets, phenocopying the maternal Dicer null phenotype. The alternative Dicer isoform, whose phylogenetic origin demonstrates evolutionary plasticity of RNA-silencing pathways, is the main determinant of endogenous RNAi activity in the mouse female germline.

The cohesin complex is crucial for chromosome segregation during mitosis and has recently also been implicated in transcriptional regulation and chromatin architecture. The NIPBL protein is required for the loading of cohesin onto chromatin, but how and where cohesin is loaded in vertebrate cells is unclear. Heterozygous mutations of NIPBL were found in 50% of the cases of Cornelia de Lange Syndrome (CdLS), a human developmental syndrome with a complex phenotype. However, no defects in the mitotic function of cohesin have been observed so far and the links between NIPBL mutations and the observed developmental defects are unclear. We show that NIPBL binds to chromatin in somatic cells with a different timing than cohesin. Further, we observe that high-affinity NIPBL binding sites localize to different regions than cohesin and almost exclusively to the promoters of active genes. NIPBL or cohesin knockdown reduce transcription of these genes differently, suggesting a cohesin-independent role of NIPBL for transcription. Motif analysis and comparison to published data show that NIPBL co-localizes with a specific set of other transcription factors. In cells derived from CdLS patients NIPBL binding levels are reduced and several of the NIPBL-bound genes have previously been observed to be mis-expressed in CdLS. In summary, our observations indicate that NIPBL mutations might cause developmental defects in different ways. First, defects of NIPBL might lead to cohesin-loading defects and thereby alter gene expression and second, NIPBL deficiency might affect genes directly via its role at the respective promoters., Author Summary The cohesin complex is crucial for chromosome segregation during cell divisions but was recently also implicated in transcriptional regulation and chromatin architecture. Cohesin's binding to chromatin depends on NIPBL, a factor that was found to be mutated in 50% of the cases of the human developmental disorder Cornelia de Lange Syndrome (CdLS). To understand the role of NIPBL for cohesin, we need to know when and where the cohesin is loaded onto DNA. Our experiments have identified high-affinity NIPBL binding sites in different cells lines which do not overlap with cohesin-binding, but colocalize with specific transcription factors at active promoters. The activity of the respective genes depends on NIPBL but not cohesin. This is in contrast with other published data showing colocalization of NIPBL and cohesin, and we reveal the existence of different types of NIPBL binding sites that are detected differently by the antibodies used in the different studies. Our observations reveal a dual role for NIPBL in cohesin loading and as potential transcription co-factor, which yields novel insights into how NIPBL defects could cause Cornelia de Lange Syndrome since NIPBL mutations might directly influence developmentally important genes.

With their unique evolutionary position of arguably the simplest extant and earliest branching metazoans, sponges are rapidly gaining focus in genomic research. Only one draft genome sequence of a sponge has been published to date, that of Amphimedon queenslandica. Genomic information from additional taxa is required for comparative genomic approaches within the phylum Porifera to describe their genomic makeup, gene architecture and structure, eventually leading to elucidation of early metazoan gene regulation and development. We sequenced and assembled a preliminary draft of the genome of the freshwater sponge Ephydatia muelleri using next-generation sequencing technology. Paired- end Illumina libraries yielded contigs covering at least 85% of the genome. Addition of SOLiD mate-pair libraries produced scaffolds covering the entire genome. We identified a preliminary gene set through automated annotation and functionally classified them using the STRING Clusters of Orthologous Groups (COG) database. Further sequencing is currently underway to improve coverage and complete the genome.