Design, construction, and functional characterization of a tRNA neochromosome in yeast.

Here, we report the design, construction, and characterization of a tRNA neochromosome, a designer chromosome that functions as an additional, de novo counterpart to the native complement of Saccharomyces cerevisiae. Intending to address one of the central design principles of the Sc2.0 project, the ∼190-kb tRNA neochromosome houses all 275 relocated nuclear tRNA genes. To maximize stability, the design incorporates orthogonal genetic elements from non-S. cerevisiae yeast species. Furthermore, the presence of 283 rox recombination sites enables an orthogonal tRNA SCRaMbLE system. Following construction in yeast, we obtained evidence of a potent selective force, manifesting as a spontaneous doubling in cell ploidy. Furthermore, tRNA sequencing, transcriptomics, proteomics, nucleosome mapping, replication profiling, FISH, and Hi-C were undertaken to investigate questions of tRNA neochromosome behavior and function. Its construction demonstrates the remarkable tractability of the yeast model and opens up opportunities to directly test hypotheses surrounding these essential non-coding RNAs.

Dicer promotes transcription termination at sites of replication stress to maintain genome stability.

Nuclear RNAi is an important regulator of transcription and epigenetic modification, but the underlying mechanisms remain elusive. Using a genome-wide approach in the fission yeast S. pombe, we have found that Dcr1, but not other components of the canonical RNAi pathway, promotes the release of Pol II from the 3? end of highly transcribed genes, and, surprisingly, from antisense transcription of rRNA and tRNA genes, which are normally transcribed by Pol I and Pol III. These Dcr1-terminated loci correspond to sites of replication stress and DNA damage, likely resulting from transcription-replication collisions. At the rDNA loci, release of Pol II facilitates DNA replication and prevents homologous recombination, which would otherwise lead to loss of rDNA repeats especially during meiosis. Our results reveal a novel role for Dcr1-mediated transcription termination in genome maintenance and may account for widespread regulation of genome stability by nuclear RNAi in higher eukaryotes.

Retraction: Oncogenic activity of Cdc6 through repression of the INK4/ARF locus.

This corrects the article DOI: 10.1038/nature04585.

Genome-wide search of nucleosome patterns using visual analytics.

MOTIVATION: The Burrows-Wheeler transform (BWT) is widely used for the fast alignment of high-throughput sequence data. This method also has potential applications in other areas of bioinformatics, and it can be specially useful for the fast searching of patterns on coverage data from different sources. RESULTS: We present a nucleosome pattern search method that converts levels of nucleosomal occupancy to a sequence-like format to which BWT searches can be applied. The method is embedded in a nucleosome map browser, 'Nucleosee', an interactive visual tool specifically designed to enhance BWT searches, giving them context and making them suitable for visual discourse analysis of the results. The proposed method is fast, flexible and sufficiently generic for the exploration of data in a broad and interactive way. AVAILABILITY AND IMPLEMENTATION: The proposed algorithm and visual browser are available for testing at http://cpg3.der.usal.es/nucleosee. The source code and installation packages are also available at https://github.com/rodrigoSantamaria/nucleosee. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Nucleosomal signatures impose nucleosome positioning in coding and noncoding sequences in the genome.

In the yeast genome, a large proportion of nucleosomes occupy well-defined and stable positions. While the contribution of chromatin remodelers and DNA binding proteins to maintain this organization is well established, the relevance of the DNA sequence to nucleosome positioning in the genome remains controversial. Through quantitative analysis of nucleosome positioning, we show that sequence changes distort the nucleosomal pattern at the level of individual nucleosomes in three species of Schizosaccharomyces and in Saccharomyces cerevisiae This effect is equally detected in transcribed and nontranscribed regions, suggesting the existence of sequence elements that contribute to positioning. To identify such elements, we incorporated information from nucleosomal signatures into artificial synthetic DNA molecules and found that they generated regular nucleosomal arrays indistinguishable from those of endogenous sequences. Strikingly, this information is species-specific and can be combined with coding information through the use of synonymous codons such that genes from one species can be engineered to adopt the nucleosomal organization of another. These findings open the possibility of designing coding and noncoding DNA molecules capable of directing their own nucleosomal organization.

A conserved role of the RSC chromatin remodeler in the establishment of nucleosome-depleted regions.

The occupancy of nucleosomes governs access to the eukaryotic genomes and results from a combination of biophysical features and the effect of ATP-dependent remodelling complexes. Most promoter regions show a conserved pattern characterized by a nucleosome-depleted region (NDR) flanked by nucleosomal arrays. The conserved RSC remodeler was reported to be critical to establish NDR in vivo in budding yeast but other evidences suggested that this activity may not be conserved in fission yeast. By reanalysing and expanding previously published data, we propose that NDR formation requires, at least partially, RSC in both yeast species. We also discuss the most prominent biological role of RSC and the possibility that non-essential subunits do not define alternate versions of the complex.

Comparative analysis of methods for genome-wide nucleosome cartography.

Nucleosomes contribute to compacting the genome into the nucleus and regulate the physical access of regulatory proteins to DNA either directly or through the epigenetic modifications of the histone tails. Precise mapping of nucleosome positioning across the genome is, therefore, essential to understanding the genome regulation. In recent years, several experimental protocols have been developed for this purpose that include the enzymatic digestion, chemical cleavage or immunoprecipitation of chromatin followed by next-generation sequencing of the resulting DNA fragments. Here, we compare the performance and resolution of these methods from the initial biochemical steps through the alignment of the millions of short-sequence reads to a reference genome to the final computational analysis to generate genome-wide maps of nucleosome occupancy. Because of the lack of a unified protocol to process data sets obtained through the different approaches, we have developed a new computational tool (NUCwave), which facilitates their analysis, comparison and assessment and will enable researchers to choose the most suitable method for any particular purpose. NUCwave is freely available at http://nucleosome.usal.es/nucwave along with a step-by-step protocol for its use.

RNAi promotes heterochromatic silencing through replication-coupled release of RNA Pol II.

Heterochromatin comprises tightly compacted repetitive regions of eukaryotic chromosomes. The inheritance of heterochromatin through mitosis requires RNA interference (RNAi), which guides histone modification during the DNA replication phase of the cell cycle. Here we show that the alternating arrangement of origins of replication and non-coding RNA in pericentromeric heterochromatin results in competition between transcription and replication in Schizosaccharomyces pombe. Co-transcriptional RNAi releases RNA polymerase II (Pol II), allowing completion of DNA replication by the leading strand DNA polymerase, and associated histone modifying enzymes that spread heterochromatin with the replication fork. In the absence of RNAi, stalled forks are repaired by homologous recombination without histone modification.

Replication dynamics in fission and budding yeasts through DNA polymerase tracking.

The dynamics of eukaryotic DNA polymerases has been difficult to establish because of the difficulty of tracking them along the chromosomes during DNA replication. Recent work has addressed this problem in the yeasts Schizosaccharomyces pombe and Saccharomyces cerevisiae through the engineering of replicative polymerases to render them prone to incorporating ribonucleotides at high rates. Their use as tracers of the passage of each polymerase has provided a picture of unprecedented resolution of the organization of replicons and replication origins in the two yeasts and has uncovered important differences between them. Additional studies have found an overlapping distribution of DNA polymorphisms and the junctions of Okazaki fragments along mononucleosomal DNA. This sequence instability is caused by the premature release of polymerase δ and the retention of non proof-read DNA tracts replicated by polymerase α. The possible implementation of these new experimental approaches in multicellular organisms opens the door to the analysis of replication dynamics under a broad range of genetic backgrounds and physiological or pathological conditions.

Nucleosomal organization of replication origins and meiotic recombination hotspots in fission yeast.

In Schizosaccharomyces pombe, DNA replication origins (ORIs) and meiotic recombination hotspots lack consensus sequences and show a bias towards mapping to large intergenic regions (IGRs). To explore whether this preference depended on underlying chromatin features, we have generated genome-wide nucleosome profiles during mitosis and meiosis. We have found that meiotic double-strand break sites (DSBs) colocalize with nucleosome-depleted regions (NDRs) and that large IGRs include clusters of NDRs that overlap with almost half of all DSBs. By contrast, ORIs do not colocalize with NDRs and they are regulated independently of DSBs. Physical relocation of NDRs at ectopic loci or modification of their genomic distribution during meiosis was paralleled by the generation of new DSB sites. Over 80% of all meiotic DSBs colocalize with NDRs that are also present during mitosis, indicating that the recombination pattern is largely dependent on constitutive properties of the genome and, to a lesser extent, on the transcriptional profile during meiosis. The organization of ORIs and of DSBs regions in S. pombe reveals similarities and differences relative to Saccharomyces cerevisiae.

PTPN13 regulates cellular signalling and ß-catenin function during megakaryocytic differentiation.

PTPN13 is a high-molecular weight intracellular phosphatase with several isoforms that exhibits a highly modular structure. Although in recent years different roles have been described for PTPN13, we are still far from understanding its function in cell biology. Here we show that PTPN13 expression is activated during megakaryocytic differentiation at the protein and mRNA level. Our results show that the upregulation of PTPN13 inhibits megakaryocytic differentiation, while PTPN13 silencing triggers differentiation. The ability of PTPN13 to alter megakaryocytic differentiation can be explained by its capacity to regulate ERK and STAT signalling. Interestingly, the silencing of ß-catenin produced the same effect as PTPN13 downregulation. We demonstrate that both proteins coimmunoprecipitate and colocalise. Moreover, we provide evidence showing that PTPN13 can regulate ß-catenin phosphorylation, stability and transcriptional activity. Therefore, the ability of PTPN13 to control megakaryocytic differentiation must be intimately linked to the regulation of ß-catenin function. Moreover, our results show for the first time that PTPN13 is stabilised upon Wnt signalling, which makes PTPN13 an important player in canonical Wnt signalling. Our results show that PTPN13 behaves as an important regulator of megakaryocytic differentiation in cell lines and also in murine haematopoietic progenitors. This importance can be explained by the ability of PTPN13 to regulate cellular signalling, and especially through the regulation of ß-catenin stability and function. Our results hold true for different megakaryocytic cell lines and also for haematopoietic progenitors, suggesting that these two proteins may play a relevant role during in vivo megakaryopoiesis.

Different nucleosomal architectures at early and late replicating origins in Saccharomyces cerevisiae.

BACKGROUND: Eukaryotic genomes are replicated during S phase according to a temporal program. Several determinants control the timing of origin firing, including the chromatin environment and epigenetic modifications. However, how chromatin structure influences the timing of the activation of specific origins is still poorly understood. RESULTS: By performing high-resolution analysis of genome-wide nucleosome positioning we have identified different chromatin architectures at early and late replication origins. These different patterns are already established in G1 and are tightly correlated with the organization of adjacent transcription units. Moreover, specific early and late nucleosomal patterns are fixed robustly, even in rpd3 mutants in which histone acetylation and origin timing have been significantly altered. Nevertheless, higher histone acetylation levels correlate with the local modulation of chromatin structure, leading to increased origin accessibility. In addition, we conducted parallel analyses of replication and nucleosome dynamics that revealed that chromatin structure at origins is modulated during origin activation. CONCLUSIONS: Our results show that early and late replication origins present distinctive nucleosomal configurations, which are preferentially associated to different genomic regions. Our data also reveal that origin structure is dynamic and can be locally modulated by histone deacetylation, as well as by origin activation. These data offer novel insight into the contribution of chromatin structure to origin selection and firing in budding yeast.

Structural diversity and dynamics of genomic replication origins in Schizosaccharomyces pombe.

DNA replication origins (ORI) in Schizosaccharomyces pombe colocalize with adenine and thymine (A+T)-rich regions, and earlier analyses have established a size from 0.5 to over 3 kb for a DNA fragment to drive replication in plasmid assays. We have asked what are the requirements for ORI function in the chromosomal context. By designing artificial ORIs, we have found that A+T-rich fragments as short as 100 bp without homology to S. pombe DNA are able to initiate replication in the genome. On the other hand, functional dissection of endogenous ORIs has revealed that some of them span a few kilobases and include several modules that may be as short as 25-30 contiguous A+Ts capable of initiating replication from ectopic chromosome positions. The search for elements with these characteristics across the genome has uncovered an earlier unnoticed class of low-efficiency ORIs that fire late during S phase. These results indicate that ORI specification and dynamics varies widely in S. pombe, ranging from very short elements to large regions reminiscent of replication initiation zones in mammals.

Asymmetrical nucleosomal DNA signatures regulate transcriptional directionality.

Despite the symmetrical structure of nucleosomes, in vitro studies have shown that transcription proceeds with different efficiency depending on the orientation of the DNA sequence around them. However, it is unclear whether this functional asymmetry is present in vivo and whether it could regulate transcriptional directionality. Here, we report that the proximal and distal halves of nucleosomal DNA contribute differentially to nucleosome stability in the genome. In +1 nucleosomes, this asymmetry facilitates or hinders transcription depending on the orientation of its underlying DNA, and this difference is associated with an asymmetrical interaction between DNA and histones. These properties are encoded in the DNA signature of +1 nucleosomes, since its incorporation in the two orientations into downstream nucleosomes renders them asymmetrically accessible to MNase and inverts the balance between sense and antisense transcription. Altogether, our results show that nucleosomal DNA endows nucleosomes with asymmetrical properties that modulate the directionality of transcription.

Clustered regulatory elements at nucleosome-depleted regions punctuate a constant nucleosomal landscape in Schizosaccharomyces pombe.

BACKGROUND: Nucleosomes facilitate the packaging of the eukaryotic genome and modulate the access of regulators to DNA. A detailed description of the nucleosomal organization under different transcriptional programmes is essential to understand their contribution to genomic regulation. RESULTS: To visualize the dynamics of individual nucleosomes under different transcriptional programmes we have generated high-resolution nucleosomal maps in Schizosaccharomyces pombe. We show that 98.5% of the genome remains almost invariable during mitosis and meiosis while remodelling is limited to approximately 1100 nucleosomes in the promoters of a subset of meiotic genes. These inducible nucleosome-depleted regions (NDR) and also those constitutively present in the genome overlap precisely with clusters of binding sites for transcription factors (TF) specific for meiosis and for different functional classes of genes, respectively. Deletion of two TFs affects only a small fraction of all the NDRs to which they bind in vivo, indicating that TFs collectively contribute to NDR maintenance. CONCLUSIONS: Our results show that the nucleosomal profile in S. pombe is largely maintained under different physiological conditions and patterns of gene expression. This relatively constant landscape favours the concentration of regulators in constitutive and inducible NDRs. The combinatorial analysis of binding motifs in this discrete fraction of the genome will facilitate the definition of the transcriptional regulatory networks.

Gcn5 facilitates Pol II progression, rather than recruitment to nucleosome-depleted stress promoters, in Schizosaccharomyces pombe.

In the fission yeast, the MAP kinase Sty1 and the transcription factor Atf1 regulate up to 400 genes in response to environmental signals, and both proteins have been shown to bind to their promoters in a stress-dependent manner. In a genetic search, we have isolated the histone H3 acetyltransferase Gcn5, a component of the SAGA complex, as being essential for oxidative stress survival and activation of those genes. Upon stress, Gcn5 is recruited to promoters and coding sequences of stress genes in a Sty1- and Atf1-dependent manner, causing both an enhanced acetylation of histone H3 and nucleosome eviction. Unexpectedly, recruitment of RNA polymerase II (Pol II) is not impaired in Δgcn5 cells. We show here that stress genes display a 400-bp long nucleosome depleted region upstream of the transcription start site even prior to activation. Stress treatment does not alter promoter nucleosome architecture, but induces eviction of the downstream nucleosomes at stress genes, which is not observed in Δgcn5 cells. We conclude that, while Pol II is recruited to nucleosome-free stress promoters in a transcription factor dependent manner, Gcn5 mediates eviction of nucleosomes positioned downstream of promoters, allowing efficient Pol II progression along the genes.

Oncogenic activity of Cdc6 through repression of the INK4/ARF locus.

The INK4/ARF locus encodes three tumour suppressors (p15(INK4b), ARF and p16(INK4a)) and is among the most frequently inactivated loci in human cancer. However, little is known about the mechanisms that govern the expression of this locus. Here we have identified a putative DNA replication origin at the INK4/ARF locus that assembles a multiprotein complex containing Cdc6, Orc2 and MCMs, and that coincides with a conserved noncoding DNA element (regulatory domain RD(INK4/ARF)). Targeted and localized RNA-interference-induced heterochromatinization of RD(INK4/ARF) results in transcriptional repression of the locus, revealing that RD(INK4/ARF) is a relevant transcriptional regulatory element. Cdc6 is overexpressed in human cancers, where it might have roles in addition to DNA replication. We have found that high levels of Cdc6 result in RD(INK4/ARF)-dependent transcriptional repression, recruitment of histone deacetylases and heterochromatinization of the INK4/ARF locus, and a concomitant decrease in the expression of the three tumour suppressors encoded by this locus. This mechanism is reminiscent of the silencing of the mating-type HM loci in yeast by replication factors. Consistent with its ability to repress the INK4/ARF locus, Cdc6 has cellular immortalization activity and neoplastic transformation capacity in cooperation with oncogenic Ras. Furthermore, human lung carcinomas with high levels of Cdc6 are associated with low levels of p16(INK4a). We conclude that aberrant expression of Cdc6 is oncogenic by directly repressing the INK4/ARF locus through the RD(INK4/ARF) element.

Genome-wide sequencing analysis of Sgs1, Exo1, Rad51, and Srs2 in DNA repair by homologous recombination.

Homologous recombination is essential to maintain genome stability in response to DNA damage. Here, we have used genome-wide sequencing to quantitatively analyze at nucleotide resolution the dynamics of DNA end resection, re-synthesis, and gene conversion at a double-strand break. Resection initiates asymmetrically in an MRX-independent manner before proceeding steadily in both directions. Sgs1, Exo1, Rad51, and Srs2 differently regulate the rate and symmetry of early and late resection. Exo1 also ensures the coexistence of resection and re-synthesis, while Srs2 guarantees a constant and symmetrical DNA re-polymerization. Gene conversion is MMR independent, spans only a minor fraction of the resected region, and its unidirectionality depends on Srs2. Finally, these repair factors prevent the development of alterations remote from the DNA lesion, such as subtelomeric instability, duplication of genomic regions, and over-replication of Ty elements. Altogether, this approach allows a quantitative analysis and a direct genome-wide visualization of DNA repair by homologous recombination.

A role for the Saccharomyces cerevisiae Rtt109 histone acetyltransferase in R-loop homeostasis and associated genome instability.

The stability of the genome is occasionally challenged by the formation of DNA-RNA hybrids and R-loops, which can be influenced by the chromatin context. This is mainly due to the fact that DNA-RNA hybrids hamper the progression of replication forks, leading to fork stalling and, ultimately, DNA breaks. Through a specific screening of chromatin modifiers performed in the yeast Saccharomyces cerevisiae, we have found that the Rtt109 histone acetyltransferase is involved in several steps of R-loop-metabolism and their associated genetic instability. On the one hand, Rtt109 prevents DNA-RNA hybridization by the acetylation of histone H3 lysines 14 and 23 and, on the other hand, it is involved in the repair of replication-born DNA breaks, such as those that can be caused by R-loops, by acetylating lysines 14 and 56. In addition, Rtt109 loss renders cells highly sensitive to replication stress in combination with R-loop-accumulating THO-complex mutants. Our data evidence that the chromatin context simultaneously influences the occurrence of DNA-RNA hybrid-associated DNA damage and its repair, adding complexity to the source of R-loop-associated genetic instability.

Analysis of DNA strand-specific differential expression with high density tiling microarrays.

BACKGROUND: DNA microarray technology allows the analysis of genome structure and dynamics at genome-wide scale. Expression microarrays (EMA) contain probes for annotated open reading frames (ORF) and are widely used for the analysis of differential gene expression. By contrast, tiling microarrays (TMA) have a much higher probe density and provide unbiased genome-wide coverage. The purpose of this study was to develop a protocol to exploit the high resolution of TMAs for quantitative measurement of DNA strand-specific differential expression of annotated and non-annotated transcripts. RESULTS: We extensively filtered probes present in Affymetrix Genechip Yeast Genome 2.0 expression and GeneChip S. pombe 1.0FR tiling microarrays to generate custom Chip Description Files (CDF) in order to compare their efficiency. We experimentally tested the potential of our approach by measuring the differential expression of 4904 genes in the yeast Schizosaccharomyces pombe growing under conditions of oxidative stress. The results showed a Pearson correlation coefficient of 0.943 between both platforms, indicating that TMAs are as reliable as EMAs for quantitative expression analysis. A significant advantage of TMAs over EMAs is the possibility of detecting non-annotated transcripts generated only under specific physiological conditions. To take full advantage of this property, we have used a target-labelling protocol that preserves the original polarity of the transcripts and, therefore, allows the strand-specific differential expression of non-annotated transcripts to be determined. By using a segmentation algorithm prior to generating the corresponding custom CDFs, we identified and quantitatively measured the expression of 510 transcripts longer than 180 nucleotides and not overlapping previously annotated ORFs that were differentially expressed at least 2-fold under oxidative stress. CONCLUSIONS: We show that the information derived from TMA hybridization can be processed simultaneously for high-resolution qualitative and quantitative analysis of the differential expression of well-characterized genes and of previously non-annotated and antisense transcripts. The consistency of the performance of TMA, their genome-wide coverage and adaptability to updated genome annotations, and the possibility of measuring strand-specific differential expression makes them a tool of choice for the analysis of gene expression in any organism for which TMA platforms are available.