Takehiko Itoh

Efficient de novo assembly of highly heterozygous genomes from whole-genome shotgun short reads

Although many de novo genome assembly projects have recently been conducted using high-throughput sequencers, assembling highly heterozygous diploid genomes is a substantial challenge due to the increased complexity of the de Bruijn graph structure predominantly used. To address the increasing demand for sequencing of nonmodel and/or wild-type samples, in most cases inbred lines or fosmid-based hierarchical sequencing methods are used to overcome such problems. However, these methods are costly and time consuming, forfeiting the advantages of massive parallel sequencing. Here, we describe a novel de novo assembler, Platanus, that can effectively manage high-throughput data from heterozygous samples. Platanus assembles DNA fragments (reads) into contigs by constructing de Bruijn graphs with automatically optimized k-mer sizes followed by the scaffolding of contigs based on paired-end information. The complicated graph structures that result from the heterozygosity are simplified during not only the contig assembly step but also the scaffolding step. We evaluated the assembly results on eukaryotic samples with various levels of heterozygosity. Compared with other assemblers, Platanus yields assembly results that have a larger scaffold NG50 length without any accompanying loss of accuracy in both simulated and real data. In addition, Platanus recorded the largest scaffold NG50 values for two of the three low-heterozygosity species used in the de novo assembly contest, Assemblathon 2. Platanus therefore provides a novel and efficient approach for the assembly of gigabase-sized highly heterozygous genomes and is an attractive alternative to the existing assemblers designed for genomes of lower heterozygosity.

HDAC8 mutations in Cornelia de Lange syndrome affect the cohesin acetylation cycle.

Cornelia de Lange syndrome (CdLS) is a dominantly inherited congenital malformation disorder, caused by mutations in the cohesin-loading protein NIPBL for nearly 60% of individuals with classical CdLS, and by mutations in the core cohesin components SMC1A (∼5%) and SMC3 (<1%) for a smaller fraction of probands. In humans, the multisubunit complex cohesin is made up of SMC1, SMC3, RAD21 and a STAG protein. These form a ring structure that is proposed to encircle sister chromatids to mediate sister chromatid cohesion and also has key roles in gene regulation. SMC3 is acetylated during S-phase to establish cohesiveness of chromatin-loaded cohesin, and in yeast, the class I histone deacetylase Hos1 deacetylates SMC3 during anaphase. Here we identify HDAC8 as the vertebrate SMC3 deacetylase, as well as loss-of-function HDAC8 mutations in six CdLS probands. Loss of HDAC8 activity results in increased SMC3 acetylation and inefficient dissolution of the ‘used’ cohesin complex released from chromatin in both prophase and anaphase. SMC3 with retained acetylation is loaded onto chromatin, and chromatin immunoprecipitation sequencing analysis demonstrates decreased occupancy of cohesin localization sites that results in a consistent pattern of altered transcription seen in CdLS cell lines with either NIPBL or HDAC8 mutations.

Identification of cis-acting sites for condensin loading onto budding yeast chromosomes

Eukaryotic chromosomes reach their stable rod-shaped appearance in mitosis in a reaction dependent on the evolutionarily conserved condensin complex. Little is known about how and where condensin associates with chromosomes. Here, we analyze condensin binding to budding yeast chromosomes using high-resolution oligonucleotide tiling arrays. Condensin-binding sites coincide with those of the loading factor Scc2/4 of the related cohesin complex. The sites map to tRNA and other genes bound by the RNA polymerase III transcription factor TFIIIC, and ribosomal protein and SNR genes. An ectopic B-box element, recognized by TFIIIC, constitutes a minimal condensin-binding site, and TFIIIC and the Scc2/4 complex promote functional condensin association with chromosomes. A similar pattern of condensin binding is conserved along fission yeast chromosomes. This reveals that TFIIIC-binding sites, including tRNA genes, constitute a hitherto unknown chromosomal feature with important implications for chromosome architecture during both interphase and mitosis.

Transcriptional Dysregulation in NIPBL and Cohesin Mutant Human Cells

Genome-wide studies using cells from patients with Cornelia de Lange Syndrome reveal a role for cohesin in regulating gene expression in human cells.

Budding Yeast Wpl1(Rad61)-Pds5 Complex Counteracts Sister Chromatid Cohesion-Establishing Reaction

Sister chromatid cohesion, which is mediated by the cohesin complex, is vital for faithful segregation of chromosomes in mitosis and meiosis (reviewed in). Cohesion is established during S phase, and this process requires the function of the acetyltransferase Eco1/Ctf7. The mechanism of the cohesion establishment is, however, still unclear. Here, we describe isolation and identification of genetic suppressors of budding yeast eco1-1 temperature-sensitive mutant. By using a recently described microarray-based method, we successfully mapped 11 intergenic suppressor mutations in two genes, wpl1 (also known as rad61) and pds5. Pds5 is a known accessory factor of cohesin complex, and we show that Wpl1/Rad61 protein forms a complex with Pds5 and colocalizes with cohesin on chromosomes, as its presumed human homolog Wapl. Impaired function of Wpl1-Pds5 complex makes Eco1 dispensable for cell survival. We also provide evidence that Wpl1 is required for efficient association of cohesin with G2 phase chromosomes and that Eco1 promotes dissociation of Wpl1-Pds5 from cohesin via acetylation of Smc3, a cohesin subunit. Taken together, the presented data suggest that Wpl1-Pds5 complex is inhibitory for cohesion establishment and that Eco1 establishes cohesion by hindering the function of Wpl1-Pds5 temporally in S phase.

The defective prophage pool of Escherichia coli O157: prophage-prophage interactions potentiate horizontal transfer of virulence determinants.

Bacteriophages are major genetic factors promoting horizontal gene transfer (HGT) between bacteria. Their roles in dynamic bacterial genome evolution have been increasingly highlighted by the fact that many sequenced bacterial genomes contain multiple prophages carrying a wide range of genes. Enterohemorrhagic Escherichia coli O157 is the most striking case. A sequenced strain (O157 Sakai) possesses 18 prophages (Sp1-Sp18) that encode numerous genes related to O157 virulence, including those for two potent cytotoxins, Shiga toxins (Stx) 1 and 2. However, most of these prophages appeared to contain multiple genetic defects. To understand whether these defective prophages have the potential to act as mobile genetic elements to spread virulence determinants, we looked closely at the Sp1-Sp18 sequences, defined the genetic defects of each Sp, and then systematically analyzed all Sps for their biological activities. We show that many of the defective prophages, including the Stx1 phage, are inducible and released from O157 cells as particulate DNA. In fact, some prophages can even be transferred to other E. coli strains. We also show that new Stx1 phages are generated by recombination between the Stx1 and Stx2 phage genomes. The results indicate that these defective prophages are not simply genetic remnants generated in the course of O157 evolution, but rather genetic elements with a high potential for disseminating virulence-related genes and other genetic traits to other bacteria. We speculate that recombination and various other types of inter-prophage interactions in the O157 prophage pool potentiate such activities. Our data provide new insights into the potential activities of the defective prophages embedded in bacterial genomes and lead to the formulation of a novel concept of inter-prophage interactions in defective prophage communities.

Replisome stability at defective DNA replication forks is independent of S phase checkpoint kinases.

The S phase checkpoint pathway preserves genome stability by protecting defective DNA replication forks, but the underlying mechanisms are still understood poorly. Previous work with budding yeast suggested that the checkpoint kinases Mec1 and Rad53 might prevent collapse of the replisome when nucleotide concentrations are limiting, thereby allowing the subsequent resumption of DNA synthesis. Here we describe a direct analysis of replisome stability in budding yeast cells lacking checkpoint kinases, together with a high-resolution view of replisome progression across the genome. Surprisingly, we find that the replisome is stably associated with DNA replication forks following replication stress in the absence of Mec1 or Rad53. A component of the replicative DNA helicase is phosphorylated within the replisome in a Mec1-dependent manner upon replication stress, and our data indicate that checkpoint kinases control replisome function rather than stability, as part of a multifaceted response that allows cells to survive defects in chromosome replication.

Csm3, Tof1, and Mrc1 Form a Heterotrimeric Mediator Complex That Associates with DNA Replication Forks

Mrc1 (mediator of replication checkpoint), Tof1 (topoisomerase I interacting factor), and Csm3 (chromosome segregation in meiosis) are checkpoint-mediator proteins that function during DNA replication and activate the effector kinase Rad53. We reported previously that Mrc1 and Tof1 are constituents of the replication machinery and that both proteins are required for the proper arrest and stabilization of replication forks in the presence of hydroxyurea. In our current study, we show that Csm3 is a component of moving replication forks and that both Tof1 and Csm3 are specifically required for the association of Mrc1 with these structures. In contrast, the deletion of mrc1 did not affect the association of Tof1 and Csm3 with the replication fork complex. In agreement with previous observations in yeast cells, the results of a baculovirus coexpression system showed that these three proteins interact directly with each other to form a mediator complex in the absence of replication forks.

Chromosome length influences replication-induced topological stress

During chromosome duplication the parental DNA molecule becomes overwound, or positively supercoiled, in the region ahead of the advancing replication fork. To allow fork progression, this superhelical tension has to be removed by topoisomerases, which operate by introducing transient DNA breaks. Positive supercoiling can also be diminished if the advancing fork rotates along the DNA helix, but then sister chromatid intertwinings form in its wake. Despite these insights it remains largely unknown how replication-induced superhelical stress is dealt with on linear, eukaryotic chromosomes. Here we show that this stress increases with the length of Saccharomyces cerevisiae chromosomes. This highlights the possibility that superhelical tension is handled on a chromosome scale and not only within topologically closed chromosomal domains as the current view predicts. We found that inhibition of type I topoisomerases leads to a late replication delay of longer, but not shorter, chromosomes. This phenotype is also displayed by cells expressing mutated versions of the cohesin- and condensin-related Smc5/6 complex. The frequency of chromosomal association sites of the Smc5/6 complex increases in response to chromosome lengthening, chromosome circularization, or inactivation of topoisomerase 2, all having the potential to increase the number of sister chromatid intertwinings. Furthermore, non-functional Smc6 reduces the accumulation of intertwined sister plasmids after one round of replication in the absence of topoisomerase 2 function. Our results demonstrate that the length of a chromosome influences the need of superhelical tension release in Saccharomyces cerevisiae, and allow us to propose a model where the Smc5/6 complex facilitates fork rotation by sequestering nascent chromatid intertwinings that form behind the replication machinery.

Rec8 Guides Canonical Spo11 Distribution along Yeast Meiotic Chromosomes

Spo11-mediated DNA double-strand breaks (DSBs) that initiate meiotic recombination are temporally and spatially controlled. The meiotic cohesin Rec8 has been implicated in regulating DSB formation, but little is known about the features of their interplay. To elucidate this point, we investigated the genome-wide localization of Spo11 in budding yeast during early meiosis by chromatin immunoprecipitation using high-density tiling arrays. We found that Spo11 is dynamically localized to meiotic chromosomes. Spo11 initially accumulated around centromeres and thereafter localized to arm regions as premeiotic S phase proceeded. During this stage, a substantial proportion of Spo11 bound to Rec8 binding sites. Eventually, some of Spo11 further bound to both DSB and Rec8 sites. We also showed that such a change in a distribution of Spo11 is affected by hydroxyurea treatment. Interestingly, deletion of REC8 influences the localization of Spo11 to centromeres and in some of the intervals of the chromosomal arms. Thus, we observed a lack of DSB formation in a region-specific manner. These observations suggest that Rec8 would prearrange the distribution of Spo11 along chromosomes and will provide clues to understanding temporal and spatial regulation of DSB formation.