Epstein-Barr virus transcription factor Zta acts through distal regulatory elements to directly control cellular gene expression

Lytic replication of the human gamma herpes virus Epstein-Barr virus (EBV) is an essential prerequisite for the spread of the virus. Differential regulation of a limited number of cellular genes has been reported in B-cells during the viral lytic replication cycle. We asked whether a viral bZIP transcription factor, Zta (BZLF1, ZEBRA, EB1), drives some of these changes. Using genome-wide chromatin immunoprecipitation coupled to next-generation DNA sequencing (ChIP-seq) we established a map of Zta interactions across the human genome. Using sensitive transcriptome analyses we identified 2263 cellular genes whose expression is significantly changed during the EBV lytic replication cycle. Zta binds 278 of the regulated genes and the distribution of binding sites shows that Zta binds mostly to sites that are distal to transcription start sites. This differs from the prevailing view that Zta activates viral genes by binding exclusively at promoter elements. We show that a synthetic Zta binding element confers Zta regulation at a distance and that distal Zta binding sites from cellular genes can confer Zta-mediated regulation on a heterologous promoter. This leads us to propose that Zta directly reprograms the expression of cellular genes through distal elements

Phylogenetic analysis of metastatic progression in breast cancer using somatic mutations and copy number aberrations.

Several studies using genome-wide molecular techniques have reported various degrees of genetic heterogeneity between primary tumours and their distant metastases. However, it has been difficult to discern patterns of dissemination owing to the limited number of patients and available metastases. Here, we use phylogenetic techniques on data generated using whole-exome sequencing and copy number profiling of primary and multiple-matched metastatic tumours from ten autopsied patients to infer the evolutionary history of breast cancer progression. We observed two modes of disease progression. In some patients, all distant metastases cluster on a branch separate from their primary lesion. Clonal frequency analyses of somatic mutations show that the metastases have a monoclonal origin and descend from a common 'metastatic precursor'. Alternatively, multiple metastatic lesions are seeded from different clones present within the primary tumour. We further show that a metastasis can be horizontally cross-seeded. These findings provide insights into breast cancer dissemination

Regulatory sites for splicing in human basal ganglia are enriched for disease-relevant information

Genome-wide association studies have generated an increasing number of common genetic variants associated with neurological and psychiatric disease risk. An improved understanding of the genetic control of gene expression in human brain is vital considering this is the likely modus operandum for many causal variants. However, human brain sampling complexities limit the explanatory power of brain-related expression quantitative trait loci (eQTL) and allele-specific expression (ASE) signals. We address this, using paired genomic and transcriptomic data from putamen and substantia nigra from 117 human brains, interrogating regulation at different RNA processing stages and uncovering novel transcripts. We identify disease-relevant regulatory loci, find that splicing eQTLs are enriched for regulatory information of neuron-specific genes, that ASEs provide cell-specific regulatory information with evidence for cellular specificity, and that incomplete annotation of the brain transcriptome limits interpretation of risk loci for neuropsychiatric disease. This resource of regulatory data is accessible through our web server, http://braineacv2.inf.um.es/

Regulation of the Na+/K+-ATPase Ena1 Expression by Calcineurin/Crz1 under High pH Stress: A Quantitative Study

[EN] Regulated expression of the Ena1 Na+-ATPase is a crucial event for adaptation to high salt and/or alkaline pH stress in the budding yeast Saccharomyces cerevisiae. ENA1 expression is under the control of diverse signaling pathways, including that mediated by the calcium-regulatable protein phosphatase calcineurin and its downstream transcription factor Crz1. We present here a quantitative study of the expression of Ena1 in response to alkalinization of the environment and we analyze the contribution of Crz1 to this response. Experimental data and mathematical models substantiate the existence of two stress-responsive Crz1-binding sites in the ENA1 promoter and estimate that the contribution of Crz1 to the early response of the ENA1 promoter is about 60%. The models suggest the existence of a second input with similar kinetics, which would be likely mediated by high pH-induced activation of the Snf1 kinase.This work was supported by grants BFU2011-30197-C3-01, BFU2014-54591-C2-1-P and EUI2009-04147 (SysMo2) to JA. (Ministry of Industry and Competitivity, Spain, and Fondo Europeo de Desarrollo Regional [FEDER]). JA is the recipient of an Ajut 2014SGR-4 award (Generalitat de Catalunya). DC was recipient of a predoctoral fellowship from the Spanish Ministry of Education. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Petrezsélyová, S.; López-Malo, M.; Canadell, D.; Roque, A.; Serra-Cardona, A.; Marques Romero, MC.; Vilaprinyó, E.... (2016). Regulation of the Na+/K+-ATPase Ena1 Expression by Calcineurin/Crz1 under High pH Stress: A Quantitative Study. PLoS ONE. 11(6):e0158424-e0158424. https://doi.org/10.1371/journal.pone.0158424Se0158424e015842411

NucTools: analysis of chromatin feature occupancy profiles from high-throughput sequencing data

Background: Biomedical applications of high-throughput sequencing methods generate a vast amount of data in which numerous chromatin features are mapped along the genome. The results are frequently analysed by creating binary data sets that link the presence/absence of a given feature to specific genomic loci. However, the nucleosome occupancy or chromatin accessibility landscape is essentially continuous. It is currently a challenge in the field to cope with continuous distributions of deep sequencing chromatin readouts and to integrate the different types of discrete chromatin features to reveal linkages between them. Results: Here we introduce the NucTools suite of Perl scripts as well as MATLAB- and R-based visualization programs for a nucleosome-centred downstream analysis of deep sequencing data. NucTools accounts for the continuous distribution of nucleosome occupancy. It allows calculations of nucleosome occupancy profiles averaged over several replicates, comparisons of nucleosome occupancy landscapes between different experimental conditions, and the estimation of the changes of integral chromatin properties such as the nucleosome repeat length. Furthermore, NucTools facilitates the annotation of nucleosome occupancy with other chromatin features like binding of transcription factors or architectural proteins, and epigenetic marks like histone modifications or DNA methylation. The applications of NucTools are demonstrated for the comparison of several datasets for nucleosome occupancy in mouse embryonic stem cells (ESCs) and mouse embryonic fibroblasts (MEFs). Conclusions: The typical workflows of data processing and integrative analysis with NucTools reveal information on the interplay of nucleosome positioning with other features such as for example binding of a transcription factor CTCF, regions with stable and unstable nucleosomes, and domains of large organized chromatin K9me2 modifications (LOCKs). As potential limitations and problems we discuss how inter-replicate variability of MNase-seq experiments can be addressed

Expression profiling during arabidopsis/downy mildew interaction reveals a highly-expressed effector that attenuates responses to salicylic acid

Plants have evolved strong innate immunity mechanisms, but successful pathogens evade or suppress plant immunity via effectors delivered into the plant cell. Hyaloperonospora arabidopsidis (Hpa) causes downy mildew on Arabidopsis thaliana, and a genome sequence is available for isolate Emoy2. Here, we exploit the availability of genome sequences for Hpa and Arabidopsis to measure gene-expression changes in both Hpa and Arabidopsis simultaneously during infection. Using a high-throughput cDNA tag sequencing method, we reveal expression patterns of Hpa predicted effectors and Arabidopsis genes in compatible and incompatible interactions, and promoter elements associated with Hpa genes expressed during infection. By resequencing Hpa isolate Waco9, we found it evades Arabidopsis resistance gene RPP1 through deletion of the cognate recognized effector ATR1. Arabidopsis salicylic acid (SA)-responsive genes including PR1 were activated not only at early time points in the incompatible interaction but also at late time points in the compatible interaction. By histochemical analysis, we found that Hpa suppresses SA-inducible PR1 expression, specifically in the haustoriated cells into which host-translocated effectors are delivered, but not in non-haustoriated adjacent cells. Finally, we found a highly-expressed Hpa effector candidate that suppresses responsiveness to SA. As this approach can be easily applied to host-pathogen interactions for which both host and pathogen genome sequences are available, this work opens the door towards transcriptome studies in infection biology that should help unravel pathogen infection strategies and the mechanisms by which host defense responses are overcome

Stressful situation if CENP-A not front and CENter

The exclusive localization of the histone H3 variant CENP-A to centromeres is essential for accurate chromosome segregation. Ubiquitin-mediated proteolysis helps to ensure that CENP-A does not mislocalize to euchromatin, which can lead to genomic instability. Consistent with this, overexpression of the budding yeast CENP-A(Cse4) is lethal in cells lacking Psh1, the E3 ubiquitin ligase that targets CENP-A(Cse4) for degradation. To identify additional mechanisms that prevent CENP-A(Cse4) misincorporation and lethality, we analyzed the genome-wide mislocalization pattern of overexpressed CENP-A(Cse4) in the presence and absence of Psh1 by chromatin immunoprecipitation followed by high throughput sequencing. We found that ectopic CENP-A(Cse4) is enriched at promoters that contain histone H2A.Z(Htz1) nucleosomes, but that H2A.Z(Htz1) is not required for CENP-A(Cse4) mislocalization. Instead, the INO80 complex, which removes H2A.Z(Htz1) from nucleosomes, promotes the ectopic deposition of CENP-A(Cse4). Transcriptional profiling revealed gene expression changes in the psh1Δ cells overexpressing CENP-A(Cse4). The down-regulated genes are enriched for CENP-A(Cse4) mislocalization to promoters, while the up-regulated genes correlate with those that are also transcriptionally up-regulated in an htz1Δ strain. Together, these data show that regulating centromeric nucleosome localization is not only critical for maintaining centromere function, but also for ensuring accurate promoter function and transcriptional regulation

Advancing brain barriers RNA sequencing: guidelines from experimental design to publication

Background: RNA sequencing (RNA-Seq) in its varied forms has become an indispensable tool for analyzing differential gene expression and thus characterization of specific tissues. Aiming to understand the brain barriers genetic signature, RNA seq has also been introduced in brain barriers research. This has led to availability of both, bulk and single-cell RNA-Seq datasets over the last few years. If appropriately performed, the RNA-Seq studies provide powerful datasets that allow for significant deepening of knowledge on the molecular mechanisms that establish the brain barriers. However, RNA-Seq studies comprise complex workflows that require to consider many options and variables before, during and after the proper sequencing process.Main body: In the current manuscript, we build on the interdisciplinary experience of the European PhD Training Network BtRAIN (https://www.btrain-2020.eu/) where bioinformaticians and brain barriers researchers collaborated to analyze and establish RNA-Seq datasets on vertebrate brain barriers. The obstacles BtRAIN has identified in this process have been integrated into the present manuscript. It provides guidelines along the entire workflow of brain barriers RNA-Seq studies starting from the overall experimental design to interpretation of results. Focusing on the vertebrate endothelial blood–brain barrier (BBB) and epithelial blood-cerebrospinal-fluid barrier (BCSFB) of the choroid plexus, we provide a step-by-step description of the workflow, highlighting the decisions to be made at each step of the workflow and explaining the strengths and weaknesses of individual choices made. Finally, we propose recommendations for accurate data interpretation and on the information to be included into a publication to ensure appropriate accessibility of the data and reproducibility of the observations by the scientific community.Conclusion: Next generation transcriptomic profiling of the brain barriers provides a novel resource for understanding the development, function and pathology of these barrier cells, which is essential for understanding CNS homeostasis and disease. Continuous advancement and sophistication of RNA-Seq will require interdisciplinary approaches between brain barrier researchers and bioinformaticians as successfully performed in BtRAIN. The present guidelines are built on the BtRAIN interdisciplinary experience and aim to facilitate collaboration of brain barriers researchers with bioinformaticians to advance RNA-Seq study design in the brain barriers community