Filling annotation gaps in yeast genomes using genome-wide contact maps

Chromosome-scale assembly of the streamlined picoeukaryote Picochlorum sp. SENEW3 genome reveals Rabl-like chromatin structure and potential for C4 photosynthesis

Genome sequencing and assembly of the photosynthetic picoeukaryotic Picochlorum sp. SENEW3 revealed a compact genome with a reduced gene set, few repetitive sequences, and an organized Rabl-like chromatin structure. Hi-C chromosome conformation capture revealed evidence of possible chromosomal translocations, as well as putative centromere locations. Maintenance of a relatively few selenoproteins, as compared to similarly sized marine picoprasinophytes Mamiellales, and broad halotolerance compared to others in Trebouxiophyceae, suggests evolutionary adaptation to variable salinity environments. Such adaptation may have driven size and genome minimization and have been enabled by the retention of a high number of membrane transporters. Identification of required pathway genes for both CAM and C4 photosynthetic carbon fixation, known to exist in the marine mamiellale pico-prasinophytes and seaweed Ulva, but few other chlorophyte species, further highlights the unique adaptations of this robust alga. This high-quality assembly provides a significant advance in the resources available for genomic investigations of this and other photosynthetic picoeukaryotes.

Metagenomes Binning Using Proximity-Ligation Data

Microbial communities are key components of all ecosystems, but characterization of their complete genomic structure remains challenging. Typical analysis tends to elude the complexity of the mixes in terms of species, strains, as well as extrachromosomal DNA molecules. Recently, approaches have been developed that bins DNA contigs into individual genomes and episomes according to their 3D contact frequencies. Those contacts are quantified by chromosome conformation capture experiments (3C, Hi-C), also known as proximity-ligation approaches, applied to metagenomics samples. Here, we present a simple computational pipeline that allows to recover high-quality Metagenomics Assemble Genomes (MAGs) starting from metagenomic 3C or Hi-C datasets and a metagenome assembly.

Chromosome-scale assembly and evolution of the tetraploid Salvia splendens (Lamiaceae) genome

Polyploidization plays a key role in plant evolution, but the forces driving the fate of homoeologs in polyploid genomes, i.e., paralogs resulting from a whole-genome duplication (WGD) event, remain to be elucidated. Here, we present a chromosome-scale genome assembly of tetraploid scarlet sage (Salvia splendens), one of the most diverse ornamental plants. We found evidence for three WGD events following an older WGD event shared by most eudicots (the γ event). A comprehensive, spatiotemporal, genome-wide analysis of homoeologs from the most recent WGD unveiled expression asymmetries, which could be associated with genomic rearrangements, transposable element proximity discrepancies, coding sequence variation, selection pressure, and transcription factor binding site differences. The observed differences between homoeologs may reflect the first step toward sub- and/or neofunctionalization. This assembly provides a powerful tool for understanding WGD and gene and genome evolution and is useful in developing functional genomics and genetic engineering strategies for scarlet sage and other Lamiaceae species.

Computer vision for pattern detection in chromosome contact maps

Chromosomes of all species studied so far display a variety of higher-order organisational features, such as self-interacting domains or loops. These structures, which are often associated to biological functions, form distinct, visible patterns on genome-wide contact maps generated by chromosome conformation capture approaches such as Hi-C. Here we present Chromosight, an algorithm inspired from computer vision that can detect patterns in contact maps. Chromosight has greater sensitivity than existing methods on synthetic simulated data, while being faster and applicable to any type of genomes, including bacteria, viruses, yeasts and mammals. Our method does not require any prior training dataset and works well with default parameters on data generated with various protocols.

Chromatin loops bridging distant loci within chromosomes can be detected by a variety of techniques such as Hi-C. Here the authors present Chromosight, an algorithm applied on mammalian, bacterial, viral and yeast genomes, able to detect various types of pattern in chromosome contact maps, including chromosomal loops.

Implications of the Evolutionary Trajectory of Centromeres in the Fungal Kingdom

Chromosome segregation during the cell cycle is an evolutionarily conserved, fundamental biological process. Dynamic interaction between spindle microtubules and the kinetochore complex that assembles on centromere DNA is required for faithful chromosome segregation. The first artificial minichromosome was constructed by cloning the centromere DNA of the budding yeast Saccharomyces cerevisiae. Since then, centromeres have been identified in >60 fungal species. The DNA sequence and organization of the sequence elements are highly diverse across these fungal centromeres. In this article, we provide a comprehensive view of the evolution of fungal centromeres. Studies of this process facilitated the identification of factors influencing centromere specification, maintenance, and propagation through many generations. Additionally, we discuss the unique features and plasticity of centromeric chromatin and the involvement of centromeres in karyotype evolution. Finally, we discuss the implications of recurrent loss of RNA interference (RNAi) and/or heterochromatin components on the trajectory of the evolution of fungal centromeres and propose the centromere structure of the last common ancestor of three major fungal phyla-Ascomycota, Basidiomycota, and Mucoromycota.

The Methylotroph Gene Order Browser (MGOB) reveals conserved synteny and ancestral centromere locations in the yeast family Pichiaceae

The yeast family Pichiaceae, also known as the 'methylotrophs clade', is a relatively little studied group of yeasts despite its economic and clinical relevance. To explore the genome evolution and synteny relationships within this family, we developed the Methylotroph Gene Order Browser (MGOB, http://mgob.ucd.ie) similar to our previous gene order browsers for other yeast families. The dataset contains genome sequences from nine Pichiaceae species, including our recent reference sequence of Pichia kudriavzevii. As an example, we demonstrate the conservation of synteny around the MOX1 locus among species both containing and lacking the MOX1 gene for methanol assimilation. We found ancient clusters of genes that are conserved as adjacent between Pichiaceae and Saccharomycetaceae. Surprisingly, we found evidence that the locations of some centromeres have been conserved among Pichiaceae species, and between Pichiaceae and Saccharomycetaceae, even though the centromeres fall into different structural categories-point centromeres, inverted repeats and retrotransposon cluster centromeres.

From dynamic chromatin architecture to DNA damage repair and back

Maintaining the integrity of the genome in the face of DNA damage is crucial to ensure the survival of the cell and normal development. DNA lesions and repair occur in the context of the chromatin fiber, whose 3D organization and movements in the restricted volume of the nucleus are under intense scrutiny. Here, we highlight work from our and other labs that addresses how the dynamic organization of the chromatin fiber affects the repair of damaged DNA and how, conversely, DNA damage and repair affect the structure and dynamics of chromatin in the budding yeast nucleus.

MetaTOR: A Computational Pipeline to Recover High-Quality Metagenomic Bins From Mammalian Gut Proximity-Ligation (meta3C) Libraries

Characterizing the complete genomic structure of complex microbial communities would represent a key step toward the understanding of their diversity, dynamics, and evolution. Current metagenomics approaches aiming at this goal are typically done by analyzing millions of short DNA sequences directly extracted from the environment. New experimental and computational approaches are constantly sought for to improve the analysis and interpretation of such data. We developed MetaTOR, an open-source computational solution that bins DNA contigs into individual genomes according to their 3D contact frequencies. Those contacts are quantified by chromosome conformation capture experiments (3C, Hi-C), also known as proximity-ligation approaches, applied to metagenomics samples (meta3C). MetaTOR was applied on 20 meta3C libraries of mice gut microbiota. We quantified the program ability to recover high-quality metagenome-assembled genomes (MAGs) from metagenomic assemblies generated directly from the meta3C libraries. Whereas nine high-quality MAGs are identified in the 148-Mb assembly generated using a single meta3C library, MetaTOR identifies 82 high-quality MAGs in the 763-Mb assembly generated from the merged 20 meta3C libraries, corresponding to nearly a third of the total assembly. Compared to the hybrid binning softwares MetaBAT or CONCOCT, MetaTOR recovered three times more high-quality MAGs. These results underline the potential of 3C-/Hi-C-based approaches in metagenomic projects.

The Impact of Centromeres on Spatial Genome Architecture

The development of new technologies and experimental techniques is enabling researchers to see what was once unable to be seen. For example, the centromere was first seen as the mediator between spindle fiber and chromosome during mitosis and meiosis. Although this continues to be its most prominent role, we now know that the centromere functions beyond cellular division with important roles in genome organization and chromatin regulation. Here we aim to share the structures and functions of centromeres in various organisms beginning with the diversity of their DNA sequence anatomies. We zoom out to describe their position in the nucleus and ultimately detail the different ways they contribute to genome organization and regulation at the spatial level.

Characterizing meiotic chromosomes' structure and pairing using a designer sequence optimized for Hi-C

In chromosome conformation capture experiments (Hi-C), the accuracy with which contacts are detected varies due to the uneven distribution of restriction sites along genomes. In addition, repeated sequences or homologous regions remain indistinguishable because of the ambiguities they introduce during the alignment of the sequencing reads. We addressed both limitations by designing and engineering 144 kb of a yeast chromosome with regularly spaced restriction sites (Syn-HiC design). In the Syn-HiC region, Hi-C signal-to-noise ratio is enhanced and can be used to measure the shape of an unbiased distribution of contact frequencies, allowing to propose a robust definition of a Hi-C experiment resolution. The redesigned region is also distinguishable from its native homologous counterpart in an otherwise isogenic diploid strain. As a proof of principle, we tracked homologous chromosomes during meiotic prophase in synchronized and pachytene-arrested cells and captured important features of their spatial reorganization, such as chromatin restructuration into arrays of Rec8-delimited loops, centromere declustering, individualization, and pairing. Overall, we illustrate the promises held by redesigning genomic regions to explore complex biological questions.

Rapid Isolation of Centromeres from Scheffersomyces stipitis

Centromeres (CENs) are the chromosomal regions promoting kinetochore formation for faithful chromosome segregation. In yeasts, CENs have been recognized as the essential elements for extra-chromosomal DNA stabilization. However, the epigeneticity of CENs makes their localization on individual chromosomes very challenging, especially in many not well-studied nonconventional yeast species. Previously, we applied a stepwise method to identify a 500-bp CEN5 from Scheffersomyces stipitis chromosome 5 and experimentally confirmed its critical role on improving plasmid stability. Here we report a library-based strategy that integrates in silico GC₃ chromosome scanning and high-throughput functional screening, which enabled the isolation of all eight S. stipitis centromeres with a 16 000-fold reduction in sequence very efficiently. Further identification of a 125-bp CEN core sequence that appears multiple times on each chromosome but all in the unique signature GC₃-valley indicates that CEN location might be accurately discerned by their local GC₃ percentages in a subgroup of yeasts.

3D organization of synthetic and scrambled chromosomes

Although the design of the synthetic yeast genome Sc2.0 is highly conservative with respect to gene content, the deletion of several classes of repeated sequences and the introduction of thousands of designer changes may affect genome organization and potentially alter cellular functions. We report here the Hi-C-determined three-dimensional (3D) conformations of Sc2.0 chromosomes. The absence of repeats leads to a smoother contact pattern and more precisely tractable chromosome conformations, and the large-scale genomic organization is globally unaffected by the presence of synthetic chromosome(s). Two exceptions are synIII, which lacks the silent mating-type cassettes, and synXII, specifically when the ribosomal DNA is moved to another chromosome. We also exploit the contact maps to detect rearrangements induced in SCRaMbLE (synthetic chromosome rearrangement and modification by loxP-mediated evolution) strains.

Proximity ligation scaffolding and comparison of two Trichoderma reesei strains genomes

BACKGROUND

The presence of low complexity and repeated regions in genomes often results in difficulties to assemble sequencing data into full chromosomes. However, the availability of full genome scaffolds is essential to several investigations, regarding for instance the evolution of entire clades, the analysis of chromosome rearrangements, and is pivotal to sexual crossing studies. In non-conventional but industrially relevant model organisms, such as the ascomycete Trichoderma reesei, a complete genome assembly is seldom available.

RESULTS

The chromosome scaffolds of T. reesei QM6a and Rut-C30 strains have been generated using a contact genomic/proximity ligation genomic approach. The original reference assembly, encompassing dozens of scaffolds, was reorganized into two sets of seven chromosomes. Chromosomal contact data also allowed to characterize 10-40 kb, gene-free, AT-rich (76%) regions corresponding to the T. reesei centromeres. Large chromosomal rearrangements (LCR) in Rut-C30 were then characterized, in agreement with former studies, and the position of LCR breakpoints used to assess the likely chromosome structure of other T. reesei strains [QM9414, CBS999.97 (1-1, re), and QM9978]. In agreement with published results, we predict that the numerous chromosome rearrangements found in highly mutated industrial strains may limit the efficiency of sexual reproduction for their improvement.

CONCLUSIONS

The GRAAL program allowed us to generate the karyotype of the Rut-C30 strain, and from there to predict chromosome structure for most T. reesei strains for which sequence is available. This method that exploits proximity ligation sequencing approach is a fast, cheap, and straightforward way to characterize both chromosome structure and centromere sequences and is likely to represent a popular convenient alternative to expensive and work-intensive resequencing projects.

The second decade of 3C technologies: detailed insights into nuclear organization

The relevance of three-dimensional (3D) genome organization for transcriptional regulation and thereby for cellular fate at large is now widely accepted. Our understanding of the fascinating architecture underlying this function is based on microscopy studies as well as the chromosome conformation capture (3C) methods, which entered the stage at the beginning of the millennium. The first decade of 3C methods rendered unprecedented insights into genome topology. Here, we provide an update of developments and discoveries made over the more recent years. As we discuss, established and newly developed experimental and computational methods enabled identification of novel, functionally important chromosome structures. Regulatory and architectural chromatin loops throughout the genome are being cataloged and compared between cell types, revealing tissue invariant and developmentally dynamic loops. Architectural proteins shaping the genome were disclosed, and their mode of action is being uncovered. We explain how more detailed insights into the 3D genome increase our understanding of transcriptional regulation in development and misregulation in disease. Finally, to help researchers in choosing the approach best tailored for their specific research question, we explain the differences and commonalities between the various 3C-derived methods.

Inferring the physical properties of yeast chromatin through Bayesian analysis of whole nucleus simulations

Background

The structure and mechanical properties of chromatin impact DNA functions and nuclear architecture but remain poorly understood. In budding yeast, a simple polymer model with minimal sequence-specific constraints and a small number of structural parameters can explain diverse experimental data on nuclear architecture. However, how assumed chromatin properties affect model predictions was not previously systematically investigated.

Results

We used hundreds of dynamic chromosome simulations and Bayesian inference to determine chromatin properties consistent with an extensive dataset that includes hundreds of measurements from imaging in fixed and live cells and two Hi-C studies. We place new constraints on average chromatin fiber properties, narrowing down the chromatin compaction to ~53–65 bp/nm and persistence length to ~52–85 nm. These constraints argue against a 20–30 nm fiber as the exclusive chromatin structure in the genome. Our best model provides a much better match to experimental measurements of nuclear architecture and also recapitulates chromatin dynamics measured on multiple loci over long timescales.

Conclusion

This work substantially improves our understanding of yeast chromatin mechanics and chromosome architecture and provides a new analytic framework to infer chromosome properties in other organisms.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1199-x) contains supplementary material, which is available to authorized users.

Naumovozyma castellii: an alternative model for budding yeast molecular biology

Naumovozyma castellii (Saccharomyces castellii) is a member of the budding yeast family Saccharomycetaceae. It has been extensively used as a model organism for telomere biology research and has gained increasing interest as a budding yeast model for functional analyses owing to its amenability to genetic modifications. Owing to the suitable phylogenetic distance to S. cerevisiae, the whole genome sequence of N. castellii has provided unique data for comparative genomic studies, and it played a key role in the establishment of the timing of the whole genome duplication and the evolutionary events that took place in the subsequent genomic evolution of the Saccharomyces lineage. Here we summarize the historical background of its establishment as a laboratory yeast species, and the development of genetic and molecular tools and strains. We review the research performed on N. castellii, focusing on areas where it has significantly contributed to the discovery of new features of molecular biology and to the advancement of our understanding of molecular evolution. Copyright © 2016 John Wiley & Sons, Ltd.

Neocentromeres Provide Chromosome Segregation Accuracy and Centromere Clustering to Multiple Loci along a Candida albicans Chromosome

Assembly of kinetochore complexes, involving greater than one hundred proteins, is essential for chromosome segregation and genome stability. Neocentromeres, or new centromeres, occur when kinetochores assemble de novo, at DNA loci not previously associated with kinetochore proteins, and they restore chromosome segregation to chromosomes lacking a functional centromere. Neocentromeres have been observed in a number of diseases and may play an evolutionary role in adaptation or speciation. However, the consequences of neocentromere formation on chromosome missegregation rates, gene expression, and three-dimensional (3D) nuclear structure are not well understood. Here, we used Candida albicans, an organism with small, epigenetically-inherited centromeres, as a model system to study the functions of twenty different neocentromere loci along a single chromosome, chromosome 5. Comparison of neocentromere properties relative to native centromere functions revealed that all twenty neocentromeres mediated chromosome segregation, albeit to different degrees. Some neocentromeres also caused reduced levels of transcription from genes found within the neocentromere region. Furthermore, like native centromeres, neocentromeres clustered in 3D with active/functional centromeres, indicating that formation of a new centromere mediates the reorganization of 3D nuclear architecture. This demonstrates that centromere clustering depends on epigenetically defined function and not on the primary DNA sequence, and that neocentromere function is independent of its distance from the native centromere position. Together, the results show that a neocentromere can form at many loci along a chromosome and can support the assembly of a functional kinetochore that exhibits native centromere functions including chromosome segregation accuracy and centromere clustering within the nucleus.

The accurate segregation of chromosomes during cell division is essential for maintaining genome integrity. The centromere is the DNA region on each chromosome where assembly of a large protein complex, the kinetochore, is required to maintain proper chromosome segregation. In addition, active centromeres exhibit a specific three-dimensional organization within the nucleus: the centromeres associate with one another in a clustered manner. Neocentromeres, or new centromeres, appear at new places along the chromosome when a native centromere becomes non-functional. We used a yeast model, Candida albicans, and isolated twenty instances in which neocentromeres had formed at different positions. All of these neocentromeres were able to direct chromosome segregation, but some had increased error rates. Like native centromeres, these neocentromeres cluster in the nucleus with the other active centromeres. This implies that formation of a neocentromere leads to reorganization of the three-dimensional structure of the nucleus so that different regions of the chromosome are in closer contact to regions of other chromosomes. Recent work suggests that approximately 3% of cancers may contain chromosomes with neocentromeres. Our observations that many neocentromeres have increased error rates provides insight into genome instability in cancer cells. Changes in chromosome copy number may benefit the cancer cells by increasing numbers of oncogenes and/or drug resistance genes, but may also sensitize the cells to chemotherapy approaches that target chromosome segregation mechanisms.

Evidence for a dual role of actin in regulating chromosome organization and dynamics in yeast

Eukaryotic chromosomes undergo movements that are involved in the regulation of functional processes such as DNA repair. To better understand the origin of these movements, we used fluorescence microscopy, image analysis and chromosome conformation capture to quantify the actin contribution to chromosome movements and interactions in budding yeast. We show that both the cytoskeletal and nuclear actin drive local chromosome movements, independently of Csm4, a putative LINC protein. Inhibition of actin polymerization reduces subtelomere dynamics, resulting in more confined territories and enrichment in subtelomeric contacts. Artificial tethering of actin to nuclear pores increased both nuclear pore complex (NPC) and subtelomere motion. Chromosome loci that were positioned away from telomeres exhibited reduced motion in the presence of an actin polymerization inhibitor but were unaffected by the lack of Csm4. We further show that actin was required for locus mobility that was induced by targeting the chromatin-remodeling protein Ino80. Correlated with this, DNA repair by homologous recombination was less efficient. Overall, interphase chromosome dynamics are modulated by the additive effects of cytoskeletal actin through forces mediated by the nuclear envelope and nuclear actin, probably through the function of actin in chromatin-remodeling complexes.

Generation and Analysis of Chromosomal Contact Maps of Yeast Species

Genome-wide derivatives of the chromosome conformation capture (3C) technique are now well-established approaches to study the multiscale average organization of chromosomes from bacteria to mammals. However, the experimental parameters of the protocol have to be optimized for different species, and the downstream experimental products (i.e., pair-end sequences) are influenced by these parameters. Here, we describe a complete pipeline to generate 3C-seq libraries and compute chromosomal contact maps of yeast species.

Spatial reorganization of telomeres in long-lived quiescent cells

Background

The spatiotemporal behavior of chromatin is an important control mechanism of genomic function. Studies in Saccharomyces cerevisiae have broadly contributed to demonstrate the functional importance of nuclear organization. Although in the wild yeast survival depends on their ability to withstand adverse conditions, most of these studies were conducted on cells undergoing exponential growth. In these conditions, as in most eukaryotic cells, silent chromatin that is mainly found at the 32 telomeres accumulates at the nuclear envelope, forming three to five foci.

Results

Here, combining live microscopy, DNA FISH and chromosome conformation capture (HiC) techniques, we report that chromosomes adopt distinct organizations according to the metabolic status of the cell. In particular, following carbon source exhaustion the genome of long-lived quiescent cells undergoes a major spatial re-organization driven by the grouping of telomeres into a unique focus or hypercluster localized in the center of the nucleus. This change in genome conformation is specific to quiescent cells able to sustain long-term viability. We further show that reactive oxygen species produced by mitochondrial activity during respiration commit the cell to form a hypercluster upon starvation. Importantly, deleting the gene encoding telomere associated silencing factor SIR3 abolishes telomere grouping and decreases longevity, a defect that is rescued by expressing a silencing defective SIR3 allele competent for hypercluster formation.

Conclusions

Our data show that mitochondrial activity primes cells to group their telomeres into a hypercluster upon starvation, reshaping the genome architecture into a conformation that may contribute to maintain longevity of quiescent cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0766-2) contains supplementary material, which is available to authorized users.

Analysis methods for studying the 3D architecture of the genome

The rapidly increasing quantity of genome-wide chromosome conformation capture data presents great opportunities and challenges in the computational modeling and interpretation of the three-dimensional genome. In particular, with recent trends towards higher-resolution high-throughput chromosome conformation capture (Hi-C) data, the diversity and complexity of biological hypotheses that can be tested necessitates rigorous computational and statistical methods as well as scalable pipelines to interpret these datasets. Here we review computational tools to interpret Hi-C data, including pipelines for mapping, filtering, and normalization, and methods for confidence estimation, domain calling, visualization, and three-dimensional modeling.

Genome-wide replication landscape of Candida glabrata

Background

The opportunistic pathogen Candida glabrata is a member of the Saccharomycetaceae yeasts. Like its close relative Saccharomyces cerevisiae, it underwent a whole-genome duplication followed by an extensive loss of genes. Its genome contains a large number of very long tandem repeats, called megasatellites. In order to determine the whole replication program of the C. glabrata genome and its general chromosomal organization, we used deep-sequencing and chromosome conformation capture experiments.

Results

We identified 253 replication fork origins, genome wide. Centromeres, HML and HMR loci, and most histone genes are replicated early, whereas natural chromosomal breakpoints are located in late-replicating regions. In addition, 275 autonomously replicating sequences (ARS) were identified during ARS-capture experiments, and their relative fitness was determined during growth competition. Analysis of ARSs allowed us to identify a 17-bp consensus, similar to the S. cerevisiae ARS consensus sequence but slightly more constrained. Megasatellites are not in close proximity to replication origins or termini. Using chromosome conformation capture, we also show that early origins tend to cluster whereas non-subtelomeric megasatellites do not cluster in the yeast nucleus.

Conclusions

Despite a shorter cell cycle, the C. glabrata replication program shares unexpected striking similarities to S. cerevisiae, in spite of their large evolutionary distance and the presence of highly repetitive large tandem repeats in C. glabrata. No correlation could be found between the replication program and megasatellites, suggesting that their formation and propagation might not be directly caused by replication fork initiation or termination.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0177-6) contains supplementary material, which is available to authorized users.

Principles of chromatin organization in yeast: relevance of polymer models to describe nuclear organization and dynamics

Nuclear organization can impact on all aspects of the genome life cycle. This organization is thoroughly investigated by advanced imaging and chromosome conformation capture techniques, providing considerable amount of datasets describing the spatial organization of chromosomes. In this review, we will focus on polymer models to describe chromosome statics and dynamics in the yeast Saccharomyces cerevisiae. We suggest that the equilibrium configuration of a polymer chain tethered at both ends and placed in a confined volume is consistent with the current literature, implying that local chromatin interactions play a secondary role in yeast nuclear organization. Future challenges are to reach an integrated multi-scale description of yeast chromosome organization, which is crucially needed to improve our understanding of the regulation of genomic transaction.

Accurate identification of centromere locations in yeast genomes using Hi-C

Centromeres are essential for proper chromosome segregation. Despite extensive research, centromere locations in yeast genomes remain difficult to infer, and in most species they are still unknown. Recently, the chromatin conformation capture assay, Hi-C, has been re-purposed for diverse applications, including de novo genome assembly, deconvolution of metagenomic samples and inference of centromere locations. We describe a method, Centurion, that jointly infers the locations of all centromeres in a single genome from Hi-C data by exploiting the centromeres’ tendency to cluster in three-dimensional space. We first demonstrate the accuracy of Centurion in identifying known centromere locations from high coverage Hi-C data of budding yeast and a human malaria parasite. We then use Centurion to infer centromere locations in 14 yeast species. Across all microbes that we consider, Centurion predicts 89% of centromeres within 5 kb of their known locations. We also demonstrate the robustness of the approach in datasets with low sequencing depth. Finally, we predict centromere coordinates for six yeast species that currently lack centromere annotations. These results show that Centurion can be used for centromere identification for diverse species of yeast and possibly other microorganisms.

Contact genomics: scaffolding and phasing (meta)genomes using chromosome 3D physical signatures

High-throughput DNA sequencing technologies are fuelling an accelerating trend to assemble de novo or resequence the genomes of numerous species as well as to complete unfinished assemblies. While current DNA sequencing technologies remain limited to reading stretches of a few hundreds or thousands of base pairs, experimental and computational methods are continuously improving with the goal of assembling entire genomes from large numbers of short DNA sequences. However, the algorithms that piece together DNA strands face important limitations due, notably, to the presence of repeated sequences or of multiple haplotypes within one genome, thus leaving many assemblies incomplete. Recently, the realization that the physical contacts experienced by a portion of a DNA molecule could be used as a robust and quantitative assay to determine its genomic position has led to the emerging field of contact genomics, which promises to revolutionize current genome assembly approaches by exploiting the flexible polymer properties of chromosomes. Here we review the current applications of contact genomics to genome scaffolding, haplotyping and metagenomic assembly, then outline the future developments we envision.

Metagenomic chromosome conformation capture (meta3C) unveils the diversity of chromosome organization in microorganisms

Genomic analyses of microbial populations in their natural environment remain limited by the difficulty to assemble full genomes of individual species. Consequently, the chromosome organization of microorganisms has been investigated in a few model species, but the extent to which the features described can be generalized to other taxa remains unknown. Using controlled mixes of bacterial and yeast species, we developed meta3C, a metagenomic chromosome conformation capture approach that allows characterizing individual genomes and their average organization within a mix of organisms. Not only can meta3C be applied to species already sequenced, but a single meta3C library can be used for assembling, scaffolding and characterizing the tridimensional organization of unknown genomes. By applying meta3C to a semi-complex environmental sample, we confirmed its promising potential. Overall, this first meta3C study highlights the remarkable diversity of microorganisms chromosome organization, while providing an elegant and integrated approach to metagenomic analysis.

DOI:http://dx.doi.org/10.7554/eLife.03318.001

Microbial communities play vital roles in the environment and sustain animal and plant life. Marine microbes are part of the ocean's food chain; soil microbes support the turnover of major nutrients and facilitate plant growth; and the microbial communities residing in the human gut support digestion and the immune system, among other roles. These communities are very complex systems, often containing 1000s of different species engaged in co-dependent relationships, and are therefore very difficult to study.

The entire DNA sequence of an organism constitutes its genome, and much of this genetic information is stored in large structures called chromosomes. Examining the genome of a species can provide important clues about its lifestyle and how it evolved. To do this, DNA is extracted from cells and is then usually cut into smaller fragments, amplified, and sequenced. The small stretches of sequence obtained, called reads, are finally assembled, yielding ideally the complete genome of the organism under study.

Metagenomics attempts to interpret the combined genome of all the different species in a microbial community and has been instrumental in deciphering how the different species interact with each other. Metagenomics involves sequencing stretches of the community's DNA and matching these pieces to individual species to ultimately assemble whole genomes. While this may be a relatively straightforward task for communities that contain only a handful of members, the metagenomes derived from complex microbial communities are huge, fragmented, and incomplete. This often makes it very difficult or even nearly impossible to match the inferred DNA stretches to individual species.

A method called chromosome conformation capture (or ‘3C’ for short) can reveal the physical contacts between different regions of a chromosome and between the different chromosomes of a cell. How often each of these chromosomal contacts occurs provides a kind of physical signature to each genome and each individual chromosome within it.

Marbouty et al. took advantage of these interactions to develop a technique that combines metagenomics and chromosome conformation capture—called meta3C—that can analyze the DNA of many different species mixed together. Testing meta3C on artificial mixtures of a few species of yeast or bacteria showed that meta3C can separate the genomes of the different species without any prior knowledge of the composition of the mix. In a single experiment, meta3C can identify individual chromosomes, match each of them to its species of origin, and reveal the three-dimensional structure of each genome in the mix. Further tests showed that meta3C can also interpret more complex communities where the number and types of the species present are not known.

Meta3C holds great promise for understanding how microbial communities work and how the genomes of the species within a community are organized. However, further developments of the technique will be required to investigate communities as diverse as those present in most natural environments.

DOI:http://dx.doi.org/10.7554/eLife.03318.002

High-quality genome (re)assembly using chromosomal contact data

Closing gaps in draft genome assemblies can be costly and time-consuming, and published genomes are therefore often left ‘unfinished.’ Here we show that genome-wide chromosome conformation capture (3C) data can be used to overcome these limitations, and present a computational approach rooted in polymer physics that determines the most likely genome structure using chromosomal contact data. This algorithm—named GRAAL—generates high-quality assemblies of genomes in which repeated and duplicated regions are accurately represented and offers a direct probabilistic interpretation of the computed structures. We first validated GRAAL on the reference genome of Saccharomyces cerevisiae, as well as other yeast isolates, where GRAAL recovered both known and unknown complex chromosomal structural variations. We then applied GRAAL to the finishing of the assembly of Trichoderma reesei and obtained a number of contigs congruent with the know karyotype of this species. Finally, we showed that GRAAL can accurately reconstruct human chromosomes from either fragments generated in silico or contigs obtained from de novo assembly. In all these applications, GRAAL compared favourably to recently published programmes implementing related approaches.

The correct assembly of genomes from sequencing data remains a challenge due to difficulties in correctly assigning the location of repeated DNA elements. Here the authors describe GRAAL, an algorithm that utilizes genome-wide chromosome contact data within a probabilistic framework to produce accurate genome assemblies.

Mating-type switching by chromosomal inversion in methylotrophic yeasts suggests an origin for the three-locus Saccharomyces cerevisiae system

Saccharomyces cerevisiae has a complex system for switching the mating type of haploid cells, requiring the genome to have three mating-type (MAT)-like loci and a mechanism for silencing two of them. How this system originated is unknown, because the three-locus system is present throughout the family Saccharomycetaceae, whereas species in the sister Candida clade have only one locus and do not switch. Here we show that yeasts in a third clade, the methylotrophs, have a simpler two-locus switching system based on reversible inversion of a section of chromosome with MATa genes at one end and MATalpha genes at the other end. In Hansenula polymorpha the 19-kb invertible region lies beside a centromere so that, depending on the orientation, either MATa or MATalpha is silenced by centromeric chromatin. In Pichia pastoris, the orientation of a 138-kb invertible region puts either MATa or MATalpha beside a telomere and represses transcription of MATa2 or MATalpha2. Both species are homothallic, and inversion of their MAT regions can be induced by crossing two strains of the same mating type. The three-locus system of S. cerevisiae, which uses a nonconservative mechanism to replace DNA at MAT, likely evolved from a conservative two-locus system that swapped genes between expression and nonexpression sites by inversion. The increasing complexity of the switching apparatus, with three loci, donor bias, and cell lineage tracking, can be explained by continuous selection to increase sporulation ability in young colonies. Our results provide an evolutionary context for the diversity of switching and silencing mechanisms.