Autotoxin-mediated latecomer killing in yeast communities.

Cellular adaptation to stressful environments such as starvation is essential to the survival of microbial communities, but the uniform response of the cell community may lead to entire cell death or severe damage to their fitness. Here, we demonstrate an elaborate response of the yeast community against glucose depletion, in which the first adapted cells kill the latecomer cells. During glucose depletion, yeast cells release autotoxins, such as leucic acid and L-2keto-3methylvalerate, which can even kill the clonal cells of the ones producing them. Although these autotoxins were likely to induce mass suicide, some cells differentiated to adapt to the autotoxins without genetic changes. If nondifferentiated latecomers tried to invade the habitat, autotoxins damaged or killed the latecomers, but the differentiated cells could selectively survive. Phylogenetically distant fission and budding yeast shared this behavior using the same autotoxins, suggesting that latecomer killing may be the universal system of intercellular communication, which may be relevant to the evolutional transition from unicellular to multicellular organisms.

Chromosome-dependent aneuploid formation in Spo11-less meiosis.

Deficiency in meiotic recombination leads to aberrant chromosome disjunction during meiosis, often resulting in the lethality of gametes or genetic disorders due to aneuploidy formation. Budding yeasts lacking Spo11, which is essential for initiation of meiotic recombination, produce many inviable spores in meiosis, while very rarely all sets of 16 chromosomes are coincidentally assorted into gametes to form viable spores. We induced meiosis in a spo11∆ diploid, in which homolog pairs can be distinguished by single nucleotide polymorphisms and determined whole-genome sequences of their exceptionally viable spores. We detected no homologous recombination in the viable spores of spo11∆ diploid. Point mutations were fewer in spo11∆ than in wild-type. We observed spo11∆ viable spores carrying a complete diploid set of homolog pairs or haploid spores with a complete haploid set of homologs but with aneuploidy in some chromosomes. In the latter, we found the chromosome-dependence in the aneuploid incidence, which was positively and negatively influenced by the chromosome length and the impact of dosage-sensitive genes, respectively. Selection of aneuploidy during meiosis II or mitosis after spore germination was also chromosome dependent. These results suggest a pathway by which specific chromosomes are more prone to cause aneuploidy, as observed in Down syndrome.

Extended TAQing system for large-scale plant genome reorganization.

We previously developed a large-scale genome restructuring technology called the TAQing system. It can induce genomic rearrangements by introducing transient and conditional formation of DNA double-strand breaks (DSBs) via heat activation of a restriction enzyme TaqI, which can cleave DNA at 5'-TCGA-3' sequences in the genome at higher temperatures (37-42°C). Such heat treatment sometimes confers lethal damage in certain plant species and TaqI cannot induce rearrangements in AT-rich regions. To overcome such problems we developed an extended TAQing (Ex-TAQing) system, which enables the use of a wider range of restriction enzymes active at standard plant-growing temperatures. We established the Ex-TAQing system using MseI that can efficiently cleave DNA at room temperature (at temperatures ranging from 22 to 25°C) and the 5'-TTAA-3' sequence which is highly abundant in the Arabidopsis genome. A synthetic intron-spanning MseI gene, which was placed downstream of a heat-shock-inducible promoter, was conditionally expressed upon milder heat treatment (33°C) to generate DSBs in Arabidopsis chromosomes. Genome resequencing revealed various types of genomic rearrangements, including copy number variations, translocation and direct end-joining at MseI cleavage sites. The Ex-TAQing system could induce large-scale rearrangements in diploids more frequently (17.4%, n = 23) than the standard TAQing system. The application of this system to tetraploids generated several strains with chromosomal rearrangements associated with beneficial phenotypes, such as high salinity stress tolerance and hypersensitivity to abscisic acid. We have developed the Ex-TAQing system, allowing more diverse patterns of genomic rearrangements, by employing various types of endonucleases and have opened a way to expand the capacity for artificial genome reorganization.

Conserved HORMA domain-containing protein Hop1 stabilizes interaction between proteins of meiotic DNA break hotspots and chromosome axis.

HORMA domain-containing proteins such as Hop1 play crucial regulatory roles in various chromosomal functions. Here, we investigated roles of the fission yeast Hop1 in the formation of recombination-initiating meiotic DNA double strand breaks (DSBs). Meiotic DSB formation in fission yeast relies on multiple protein-protein interactions such as the one between the chromosome axial protein Rec10 and the DSB-forming complex subunit Rec15. Chromatin immunoprecipitation sequencing demonstrated that Hop1 is colocalized with both Rec10 and Rec15, and we observed physical interactions of Hop1 to Rec15 and Rec10. These results suggest that Hop1 promotes DSB formation by interacting with both axis components and the DSB-forming complex. We also show that Hop1 binding to DSB hotspots requires Rec15 and Rec10, while Hop1 axis binding requires Rec10 only, suggesting that Hop1 is recruited to the axis via Rec10, and to hotspots by hotspot-bound Rec15. Furthermore, we introduced separation-of-function Rec10 mutations, deficient for interaction with either Rec15 or Hop1. These single mutations and hop1Δ conferred only partial defects in meiotic recombination, while the combining the Rec15-binding-deficient rec10 mutation with hop1Δ synergistically reduced meiotic recombination, at least at a model hotspot. Taken together, Hop1 likely functions as a stabilizer for Rec15-Rec10 interaction to promote DSB formation.

A central coupler for recombination initiation linking chromosome architecture to S phase checkpoint.

Higher-order chromosome structure is assumed to control various DNA-templated reactions in eukaryotes. Meiotic chromosomes implement developed structures called "axes" and "loops"; both are suggested to tether each other, activating Spo11 to catalyze meiotic DNA double-strand breaks (DSBs) at recombination hotspots. We found that the Schizosaccharomyces pombe Spo11 homolog Rec12 and its partners form two distinct subcomplexes, DSBC (Rec6-Rec12-Rec14) and SFT (Rec7-Rec15-Rec24). Mde2, whose expression is strictly regulated by the replication checkpoint, interacts with Rec15 to stabilize the SFT subcomplex and further binds Rec14 in DSBC. Rec10 provides a docking platform for SFT binding to axes and can partially interact with DSB sites located in loops depending upon Mde2, which is indicative of the formation of multiprotein-based tethered axis-loop complex. These data lead us to propose a mechanism by which Mde2 functions as a recombination initiation mediator to tether axes and loops, in liaison with the meiotic replication checkpoint.

Replication stress induces accumulation of FANCD2 at central region of large fragile genes.

During mild replication stress provoked by low dose aphidicolin (APH) treatment, the key Fanconi anemia protein FANCD2 accumulates on common fragile sites, observed as sister foci, and protects genome stability. To gain further insights into FANCD2 function and its regulatory mechanisms, we examined the genome-wide chromatin localization of FANCD2 in this setting by ChIP-seq analysis. We found that FANCD2 mostly accumulates in the central regions of a set of large transcribed genes that were extensively overlapped with known CFS. Consistent with previous studies, we found that this FANCD2 retention is R-loop-dependent. However, FANCD2 monoubiquitination and RPA foci formation were still induced in cells depleted of R-loops. Interestingly, we detected increased Proximal Ligation Assay dots between FANCD2 and R-loops following APH treatment, which was suppressed by transcriptional inhibition. Collectively, our data suggested that R-loops are required to retain FANCD2 in chromatin at the middle intronic region of large genes, while the replication stress-induced upstream events leading to the FA pathway activation are not triggered by R-loops.

Local potentiation of stress-responsive genes by upstream noncoding transcription.

It has been postulated that a myriad of long noncoding RNAs (lncRNAs) contribute to gene regulation. In fission yeast, glucose starvation triggers lncRNA transcription across promoter regions of stress-responsive genes including fbp1 (fructose-1,6-bisphosphatase1). At the fbp1 promoter, this transcription promotes chromatin remodeling and fbp1 mRNA expression. Here, we demonstrate that such upstream noncoding transcription facilitates promoter association of the stress-responsive transcriptional activator Atf1 at the sites of transcription, leading to activation of the downstream stress genes. Genome-wide analyses revealed that ∼50 Atf1-binding sites show marked decrease in Atf1 occupancy when cells are treated with a transcription inhibitor. Most of these transcription-enhanced Atf1-binding sites are associated with stress-dependent induction of the adjacent mRNAs or lncRNAs, as observed in fbp1 These Atf1-binding sites exhibit low Atf1 occupancy and high histone density in glucose-rich conditions, and undergo dramatic changes in chromatin status after glucose depletion: enhanced Atf1 binding, histone eviction, and histone H3 acetylation. We also found that upstream transcripts bind to the Groucho-Tup1 type transcriptional corepressors Tup11 and Tup12, and locally antagonize their repressive functions on Atf1 binding. These results reveal a new mechanism in which upstream noncoding transcription locally magnifies the specific activation of stress-inducible genes via counteraction of corepressors.

Dynamic transition of transcription and chromatin landscape during fission yeast adaptation to glucose starvation.

Shortage of glucose, the primary energy source for all organisms, is one of the most critical stresses influencing cell viability. Glucose starvation promptly induces changes in mRNA and noncoding RNA (ncRNA) transcription. We previously reported that glucose starvation induces long ncRNA (lncRNA) transcription in the 5' segment of a fission yeast gluconeogenesis gene (fbp1+), which leads to stepwise chromatin alteration around the fbp1+ promoter and to subsequent robust gene activation. Here, we analyzed genomewide transcription by strand-specific RNA sequencing, together with chromatin landscape by immunoprecipitation sequencing (ChIP-seq). Clustering analysis showed that distinct mRNAs and ncRNAs are induced at the early, middle and later stages of cellular response to glucose starvation. The starvation-induced transcription depends substantially on the stress-responsive transcription factor Atf1. Using a new computer program that examines dynamic changes in expression patterns, we identified ncRNAs with similar behavior to the fbp1+ lncRNA. We confirmed that there are continuous lncRNAs associated with local reduction of histone density. Overlapping with the regions for transcription of these lncRNAs, antisense RNAs are antagonistically transcribed under glucose-rich conditions. These results suggest that Atf1-dependent integrated networks of mRNA and lncRNA govern drastic changes in cell physiology in response to glucose starvation.

Stress-induced lncRNAs evade nuclear degradation and enter the translational machinery.

Long noncoding RNAs (lncRNAs) play important roles in the regulation of gene expression. In fission yeast, glucose starvation triggers a transcriptional cascade of polyadenylated lncRNAs in the upstream region of the fructose-1,6-bisphosphatase gene (fbp1(+) ), which is correlated with stepwise chromatin remodeling and necessary for the massive induction of fbp1(+) mRNA. Here, we show that these novel metabolic stress-induced lncRNAs (mlonRNAs) are 5'-capped, less stable than fbp1(+) mRNA and sensitive to a certain extent to the nuclear exosome cofactor Rrp6. However, most mlonRNAs seem to escape nuclear degradation and are exported to the cytoplasm, where they localize to polysomes precisely during glucose starvation-induced global translation inhibition. It is likely that ribosomes tend to accumulate in the upstream region of mlonRNAs. Although mlonRNAs contain an unusual amount of upstream AUGs (uAUGs) and small open reading frames (uORFs), they escape Upf1-mediated targeting to the non-sense-mediated decay (NMD) pathway. The deletion of Upf1 had no effect on mlonRNA stability, but considerably destabilized fbp1(+) mRNA, hinting toward a possible novel role of Upf1. Our findings suggest that the stability of mlonRNAs is distinctly regulated from mRNA and previously described noncoding transcripts.

Rapamycin-sensitive mechanisms confine the growth of fission yeast below the temperatures detrimental to cell physiology.

Cells cease to proliferate above their growth-permissible temperatures, a ubiquitous phenomenon generally attributed to heat damage to cellular macromolecules. We here report that, in the presence of rapamycin, a potent inhibitor of Target of Rapamycin Complex 1 (TORC1), the fission yeast Schizosaccharomyces pombe can proliferate at high temperatures that usually arrest its growth. Consistently, mutations to the TORC1 subunit RAPTOR/Mip1 and the TORC1 substrate Sck1 significantly improve cellular heat resistance, suggesting that TORC1 restricts fission yeast growth at high temperatures. Aiming for a more comprehensive understanding of the negative regulation of high-temperature growth, we conducted genome-wide screens, which identified additional factors that suppress cell proliferation at high temperatures. Among them is Mks1, which is phosphorylated in a TORC1-dependent manner, forms a complex with the 14-3-3 protein Rad24, and suppresses the high-temperature growth independently of Sck1. Our study has uncovered unexpected mechanisms of growth restraint even below the temperatures deleterious to cell physiology.

Fission Yeast TORC1 Promotes Cell Proliferation through Sfp1, a Transcription Factor Involved in Ribosome Biogenesis.

Target of rapamycin complex 1 (TORC1) is activated in response to nutrient availability and growth factors, promoting cellular anabolism and proliferation. To explore the mechanism of TORC1-mediated proliferation control, we performed a genetic screen in fission yeast and identified Sfp1, a zinc-finger transcription factor, as a multicopy suppressor of temperature-sensitive TORC1 mutants. Our observations suggest that TORC1 phosphorylates Sfp1 and protects Sfp1 from proteasomal degradation. Transcription analysis revealed that Sfp1 positively regulates genes involved in ribosome production together with two additional transcription factors, Ifh1/Crf1 and Fhl1. Ifh1 physically interacts with Fhl1, and the nuclear localization of Ifh1 is regulated in response to nutrient levels in a manner dependent on TORC1 and Sfp1. Taken together, our data suggest that the transcriptional regulation of the genes involved in ribosome biosynthesis by Sfp1, Ifh1, and Fhl1 is one of the key pathways through which nutrient-activated TORC1 promotes cell proliferation.

Imputation-free reconstructions of three-dimensional chromosome architectures in human diploid single-cells using allele-specified contacts.

Single-cell Hi-C analysis of diploid human cells is difficult because of the lack of dense chromosome contact information and the presence of homologous chromosomes with very similar nucleotide sequences. Thus here, we propose a new algorithm to reconstruct the three-dimensional (3D) chromosomal architectures from the Hi-C dataset of single diploid human cells using allele-specific single-nucleotide variations (SNVs). We modified our recurrence plot-based algorithm, which is suitable for the estimation of the 3D chromosome structure from sparse Hi-C datasets, by newly incorporating a function of discriminating SNVs specific to each homologous chromosome. Here, we eventually regard a contact map as a recurrence plot. Importantly, the proposed method does not require any imputation for ambiguous segment information, but could efficiently reconstruct 3D chromosomal structures in single human diploid cells at a 1-Mb resolution. Datasets of segments without allele-specific SNVs, which were considered to be of little value, can also be used to validate the estimated chromosome structure. Introducing an additional mathematical measure called a refinement further improved the resolution to 40-kb or 100-kb. The reconstruction data supported the notion that human chromosomes form chromosomal territories and take fractal structures where the dimension for the underlying chromosome structure is a non-integer value.

TAQing2.0 for genome reorganization of asexual industrial yeasts by direct protein transfection.

Genomic rearrangements often generate phenotypic diversification. We previously reported the TAQing system where genomic rearrangements are induced via conditional activation of a restriction endonuclease in yeast and plant cells to produce mutants with marked phenotypic changes. Here we developed the TAQing2.0 system based on the direct delivery of endonucleases into the cell nucleus by cell-penetrating peptides. Using the optimized procedure, we introduce a heat-reactivatable endonuclease TaqI into an asexual industrial yeast (torula yeast), followed by a transient heat activation of TaqI. TAQing2.0 leads to generation of mutants with altered flocculation and morphological phenotypes, which exhibit changes in chromosomal size. Genome resequencing suggested that torula yeast is triploid with six chromosomes and the mutants have multiple rearrangements including translocations having the TaqI recognition sequence at the break points. Thus, TAQing2.0 is expected as a useful method to obtain various mutants with altered phenotypes without introducing foreign DNA into asexual industrial microorganisms.

lncRNA transcription induces meiotic recombination through chromatin remodelling in fission yeast.

Noncoding RNAs (ncRNAs) are involved in various biological processes, including gene expression, development, and disease. Here, we identify a novel consensus sequence of a cis-element involved in long ncRNA (lncRNA) transcription and demonstrate that lncRNA transcription from this cis-element activates meiotic recombination via chromatin remodeling. In the fission yeast fbp1 gene, glucose starvation induces a series of promoter-associated lncRNAs, referred to as metabolic-stress-induced lncRNAs (mlonRNAs), which contribute to chromatin remodeling and fbp1 activation. Translocation of the cis-element required for mlonRNA into a well-characterized meiotic recombination hotspot, ade6-M26, further stimulates transcription and meiotic recombination via local chromatin remodeling. The consensus sequence of this cis-element (mlon-box) overlaps with meiotic recombination sites in the fission yeast genome. At one such site, the SPBC24C6.09c upstream region, meiotic double-strand break (DSB) formation is induced in an mlon-box-dependent manner. Therefore, mlonRNA transcription plays a universal role in chromatin remodeling and the regulation of transcription and recombination.

Linear Regression Links Transcriptomic Data and Cellular Raman Spectra.

Raman microscopy is an imaging technique that has been applied to assess molecular compositions of living cells to characterize cell types and states. However, owing to the diverse molecular species in cells and challenges of assigning peaks to specific molecules, it has not been clear how to interpret cellular Raman spectra. Here, we provide firm evidence that cellular Raman spectra and transcriptomic profiles of Schizosaccharomyces pombe and Escherichia coli can be computationally connected and thus interpreted. We find that the dimensions of high-dimensional Raman spectra and transcriptomes measured by RNA sequencing can be reduced and connected linearly through a shared low-dimensional subspace. Accordingly, we were able to predict global gene expression profiles by applying the calculated transformation matrix to Raman spectra, and vice versa. Highly expressed non-coding RNAs contributed to the Raman-transcriptome linear correspondence more significantly than mRNAs in S. pombe. This demonstration of correspondence between cellular Raman spectra and transcriptomes is a promising step toward establishing spectroscopic live-cell omics studies.

Phenotypic diversification by enhanced genome restructuring after induction of multiple DNA double-strand breaks.

DNA double-strand break (DSB)-mediated genome rearrangements are assumed to provide diverse raw genetic materials enabling accelerated adaptive evolution; however, it remains unclear about the consequences of massive simultaneous DSB formation in cells and their resulting phenotypic impact. Here, we establish an artificial genome-restructuring technology by conditionally introducing multiple genomic DSBs in vivo using a temperature-dependent endonuclease TaqI. Application in yeast and Arabidopsis thaliana generates strains with phenotypes, including improved ethanol production from xylose at higher temperature and increased plant biomass, that are stably inherited to offspring after multiple passages. High-throughput genome resequencing revealed that these strains harbor diverse rearrangements, including copy number variations, translocations in retrotransposons, and direct end-joinings at TaqI-cleavage sites. Furthermore, large-scale rearrangements occur frequently in diploid yeasts (28.1%) and tetraploid plants (46.3%), whereas haploid yeasts and diploid plants undergo minimal rearrangement. This genome-restructuring system (TAQing system) will enable rapid genome breeding and aid genome-evolution studies.

Correlation of Meiotic DSB Formation and Transcription Initiation Around Fission Yeast Recombination Hotspots.

Meiotic homologous recombination, a critical event for ensuring faithful chromosome segregation and creating genetic diversity, is initiated by programmed DNA double-strand breaks (DSBs) formed at recombination hotspots. Meiotic DSB formation is likely to be influenced by other DNA-templated processes including transcription, but how DSB formation and transcription interact with each other has not been understood well. In this study, we used fission yeast to investigate a possible interplay of these two events. A group of hotspots in fission yeast are associated with sequences similar to the cyclic AMP response element and activated by the ATF/CREB family transcription factor dimer Atf1-Pcr1. We first focused on one of those hotspots, ade6-3049, and Atf1. Our results showed that multiple transcripts, shorter than the ade6 full-length messenger RNA, emanate from a region surrounding the ade6-3049 hotspot. Interestingly, we found that the previously known recombination-activation region of Atf1 is also a transactivation domain, whose deletion affected DSB formation and short transcript production at ade6-3049 These results point to a possibility that the two events may be related to each other at ade6-3049 In fact, comparison of published maps of meiotic transcripts and hotspots suggested that hotspots are very often located close to meiotically transcribed regions. These observations therefore propose that meiotic DSB formation in fission yeast may be connected to transcription of surrounding regions.

Three-dimensional reconstruction of single-cell chromosome structure using recurrence plots.

Single-cell analysis of the three-dimensional (3D) chromosome structure can reveal cell-to-cell variability in genome activities. Here, we propose to apply recurrence plots, a mathematical method of nonlinear time series analysis, to reconstruct the 3D chromosome structure of a single cell based on information of chromosomal contacts from genome-wide chromosome conformation capture (Hi-C) data. This recurrence plot-based reconstruction (RPR) method enables rapid reconstruction of a unique structure in single cells, even from incomplete Hi-C information.