Multiscale Structuring of the E. coli Chromosome by Nucleoid-Associated and Condensin Proteins.

As in eukaryotes, bacterial genomes are not randomly folded. Bacterial genetic information is generally carried on a circular chromosome with a single origin of replication from which two replication forks proceed bidirectionally toward the opposite terminus region. Here, we investigate the higher-order architecture of the Escherichia coli genome, showing its partition into two structurally distinct entities by a complex and intertwined network of contacts: the replication terminus (ter) region and the rest of the chromosome. Outside of ter, the condensin MukBEF and the ubiquitous nucleoid-associated protein (NAP) HU promote DNA contacts in the megabase range. Within ter, the MatP protein prevents MukBEF activity, and contacts are restricted to ∼280 kb, creating a domain with distinct structural properties. We also show how other NAPs contribute to nucleoid organization, such as H-NS, which restricts short-range interactions. Combined, these results reveal the contributions of major evolutionarily conserved proteins in a bacterial chromosome organization.

Condensin- and Replication-Mediated Bacterial Chromosome Folding and Origin Condensation Revealed by Hi-C and Super-resolution Imaging.

Chromosomes of a broad range of species, from bacteria to mammals, are structured by large topological domains whose precise functional roles and regulatory mechanisms remain elusive. Here, we combine super-resolution microscopies and chromosome-capture technologies to unravel the higher-order organization of the Bacillus subtilis chromosome and its dynamic rearrangements during the cell cycle. We decipher the fine 3D architecture of the origin domain, revealing folding motifs regulated by condensin-like complexes. This organization, along with global folding throughout the genome, is present before replication, disrupted by active DNA replication, and re-established thereafter. Single-cell analysis revealed a strict correspondence between sub-cellular localization of origin domains and their condensation state. Our results suggest that the precise 3D folding pattern of the origin domain plays a role in the regulation of replication initiation, chromosome organization, and DNA segregation.

Metagenome Analysis Exploiting High-Throughput Chromosome Conformation Capture (3C) Data.

Microbial communities are complex and constitute important parts of our environment. Genomic analysis of these populations is a dynamic research area but remains limited by the difficulty in assembling full genomes of individual species. Recently, a new method for metagenome assembly/analysis based on chromosome conformation capture has emerged (meta3C). This approach quantifies the collisions experienced by DNA molecules to identify those sharing the same cellular compartments, allowing the characterization of genomes present within complex mixes of species. The exploitation of these chromosome 3D signatures holds promising perspectives for genome sequencing of discrete species in complex populations. It also has the potential to assign correctly extra-chromosomal elements, such as plasmids, mobile elements and phages, to their host cells.

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.

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.

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

MOTIVATIONS: De novo sequencing of genomes is followed by annotation analyses aiming at identifying functional genomic features such as genes, non-coding RNAs or regulatory sequences, taking advantage of diverse datasets. These steps sometimes fail at detecting non-coding functional sequences: for example, origins of replication, centromeres and rDNA positions have proven difficult to annotate with high confidence. Here, we demonstrate an unconventional application of Chromosome Conformation Capture (3C) technique, which typically aims at deciphering the average 3D organization of genomes, by showing how functional information about the sequence can be extracted solely from the chromosome contact map. RESULTS: Specifically, we describe a combined experimental and bioinformatic procedure that determines the genomic positions of centromeres and ribosomal DNA clusters in yeasts, including species where classical computational approaches fail. For instance, we determined the centromere positions in Naumovozyma castellii, where these coordinates could not be obtained previously. Although computed centromere positions were characterized by conserved synteny with neighboring species, no consensus sequences could be found, suggesting that centromeric binding proteins or mechanisms have significantly diverged. We also used our approach to refine centromere positions in Kuraishia capsulata and to identify rDNA positions in Debaryomyces hansenii. Our study demonstrates how 3C data can be used to complete the functional annotation of eukaryotic genomes. AVAILABILITY AND IMPLEMENTATION: The source code is provided in the Supplementary Material. This includes a zipped file with the Python code and a contact matrix of Saccharomyces cerevisiae. CONTACT: romain.koszul@pasteur.fr SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Purification of G1 daughter cells from different Saccharomycetes species through an optimized centrifugal elutriation procedure.

Centrifugal elutriation discriminates cells according to their sedimentation coefficients, generating homogeneous samples well suited for genomic comparative approaches. It can, for instance, isolate G1 daughter cells from a Saccharomyces cerevisiae unsynchronized population, alleviating ageing and cell-cycle biases when conducting genome-wide/single-cell studies. The present report describes a straightforward and robust procedure to determine whether a cell population of virtually any yeast species can be efficiently elutriated, while offering solutions to optimize success. This approach was used to characterize elutriation parameters and S-phase progression of four yeast species (S. cerevisiae, Candida glabrata, Lachancea kluyveri and Pichia sorbitophila) and could theoretically be applied to any culture of single, individual cells.

Phage-inducible chromosomal minimalist islands (PICMIs), a novel family of small marine satellites of virulent phages.

Phage satellites are bacterial genetic elements that co-opt phage machinery for their own dissemination. Here we identify a family of satellites, named Phage-Inducible Chromosomal Minimalist Islands (PICMIs), that are broadly distributed in marine bacteria of the family Vibrionaceae. A typical PICMI is characterized by reduced gene content, does not encode genes for capsid remodelling, and packages its DNA as a concatemer. PICMIs integrate in the bacterial host genome next to the fis regulator, and encode three core proteins necessary for excision and replication. PICMIs are dependent on virulent phage particles to spread to other bacteria, and protect their hosts from other competitive phages without interfering with their helper phage. Thus, our work broadens our understanding of phage satellites and narrows down the minimal number of functions necessary to hijack a tailed phage.

Chromosome-level genome assembly of the European green woodpecker Picus viridis.

The European green woodpecker, Picus viridis, is a widely distributed species found in the Western Palearctic region. Here, we assembled a highly contiguous genome assembly for this species using a combination of short- and long-read sequencing and scaffolded with chromatin conformation capture (Hi-C). The final genome assembly was 1.28 Gb and features a scaffold N50 of 37 Mb and a scaffold L50 of 39.165 Mb. The assembly incorporates 89.4% of the genes identified in birds in OrthoDB. Gene and repetitive content annotation on the assembly detected 15,805 genes and a ∼30.1% occurrence of repetitive elements, respectively. Analysis of synteny demonstrates the fragmented nature of the P. viridis genome when compared to the chicken (Gallus gallus). The assembly and annotations produced in this study will certainly help for further research into the genomics of P. viridis and the comparative evolution of woodpeckers. Five historical and seven contemporary samples have been resequenced and may give insights on the population history of this species.

Transcription-induced domains form the elementary constraining building blocks of bacterial chromosomes.

Transcription generates local topological and mechanical constraints on the DNA fiber, leading to the generation of supercoiled chromosome domains in bacteria. However, the global impact of transcription on chromosome organization remains elusive, as the scale of genes and operons in bacteria remains well below the resolution of chromosomal contact maps generated using Hi-C (~5-10 kb). Here we combined sub-kb Hi-C contact maps and chromosome engineering to visualize individual transcriptional units. We show that transcriptional units form discrete three-dimensional transcription-induced domains that impose mechanical and topological constraints on their neighboring sequences at larger scales, modifying their localization and dynamics. These results show that transcriptional domains constitute primary building blocks of bacterial chromosome folding and locally impose structural and dynamic constraints.

Chromosome folding and prophage activation reveal specific genomic architecture for intestinal bacteria.

BACKGROUND: Bacteria and their viruses, bacteriophages, are the most abundant entities of the gut microbiota, a complex community of microorganisms associated with human health and disease. In this ecosystem, the interactions between these two key components are still largely unknown. In particular, the impact of the gut environment on bacteria and their associated prophages is yet to be deciphered. RESULTS: To gain insight into the activity of lysogenic bacteriophages within the context of their host genomes, we performed proximity ligation-based sequencing (Hi-C) in both in vitro and in vivo conditions on the 12 bacterial strains of the OMM12 synthetic bacterial community stably associated within mice gut (gnotobiotic mouse line OMM12). High-resolution contact maps of the chromosome 3D organization of the bacterial genomes revealed a wide diversity of architectures, differences between environments, and an overall stability over time in the gut of mice. The DNA contacts pointed at 3D signatures of prophages leading to 16 of them being predicted as functional. We also identified circularization signals and observed different 3D patterns between in vitro and in vivo conditions. Concurrent virome analysis showed that 11 of these prophages produced viral particles and that OMM12 mice do not carry other intestinal viruses. CONCLUSIONS: The precise identification by Hi-C of functional and active prophages within bacterial communities will unlock the study of interactions between bacteriophages and bacteria across conditions (healthy vs disease). Video Abstract.

Normalization of a chromosomal contact map.

BACKGROUND: Chromatin organization has been increasingly studied in relation with its important influence on DNA-related metabolic processes such as replication or regulation of gene expression. Since its original design ten years ago, capture of chromosome conformation (3C) has become an essential tool to investigate the overall conformation of chromosomes. It relies on the capture of long-range trans and cis interactions of chromosomal segments whose relative proportions in the final bank reflect their frequencies of interactions, hence their spatial proximity in a population of cells. The recent coupling of 3C with deep sequencing approaches now allows the generation of high resolution genome-wide chromosomal contact maps. Different protocols have been used to generate such maps in various organisms. This includes mammals, drosophila and yeast. The massive amount of raw data generated by the genomic 3C has to be carefully processed to alleviate the various biases and byproducts generated by the experiments. Our study aims at proposing a simple normalization procedure to minimize the influence of these unwanted but inevitable events on the final results. RESULTS: Careful analysis of the raw data generated previously for budding yeast S. cerevisiae led to the identification of three main biases affecting the final datasets, including a previously unknown bias resulting from the circularization of DNA molecules. We then developed a simple normalization procedure to process the data and allow the generation of a normalized, highly contrasted, chromosomal contact map for S. cerevisiae. The same method was then extended to the first human genome contact map. Using the normalized data, we revisited the preferential interactions originally described between subsets of discrete chromosomal features. Notably, the detection of preferential interactions between tRNA in yeast and CTCF, PolII binding sites in human can vary with the normalization procedure used. CONCLUSIONS: We quantitatively reanalyzed the genomic 3C data obtained for S. cerevisiae, identified some of the biases inherent to the technique and proposed a simple normalization procedure to analyse them. Such an approach can be easily generalized for genomic 3C experiments in other organisms. More experiments and analysis will be necessary to reach optimal resolution and accuracies of the maps generated through these approaches. Working with cell population presenting highest levels of homogeneity will prove useful in this regards.

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.

Bacterial Membrane Vesicles as a Novel Strategy for Extrusion of Antimicrobial Bismuth Drug in Helicobacter pylori.

Bacterial antibiotic resistance is a major threat to human health. A combination of antibiotics with metals is among the proposed alternative treatments. Only one such combination is successfully used in clinics; it associates antibiotics with the metal bismuth to treat infections by Helicobacter pylori. This bacterial pathogen colonizes the human stomach and is associated with gastric cancer, killing 800,000 individuals yearly. The effect of bismuth in H. pylori treatment is not well understood in particular for sublethal doses such as those measured in the plasma of treated patients. We addressed this question and observed that bismuth induces the formation of homogeneously sized membrane vesicles (MVs) with unique protein cargo content enriched in bismuth-binding proteins, as shown by quantitative proteomics. Purified MVs of bismuth-exposed bacteria were strongly enriched in bismuth as measured by inductively coupled plasma optical emission spectrometry (ICP-OES), unlike bacterial cells from which they originate. Thus, our results revealed a novel function of MVs in bismuth detoxification, where secreted MVs act as tool to discard bismuth from the bacteria. Bismuth also induces the formation of intracellular polyphosphate granules that are associated with changes in nucleoid structure. Nucleoid compaction in response to bismuth was established by immunogold electron microscopy and refined by the first chromosome conformation capture (Hi-C) analysis of H. pylori. Our results reveal that even low doses of bismuth induce profound changes in H. pylori physiology and highlight a novel defense mechanism that involves MV-mediated bismuth extrusion from the bacteria and a probable local DNA protective response where polyphosphate granules are associated with nucleoid compaction. IMPORTANCE Bacterial resistance to antibiotics is a major threat to human health. Treatments combining antibiotics with metals were proposed to circumvent this hurdle. Only one such combination is successfully used in clinics associating antibiotics with the metal bismuth to treat infections by the human pathogen Helicobacter pylori. H. pylori causes 800,000 deaths by gastric cancer yearly. How bismuth impacts H. pylori and its response to this toxic metal were ill defined. We discovered that upon bismuth exposure, H. pylori secretes membrane vesicles that are enriched in bismuth. Bismuth also induces the formation of intracellular polyphosphate granules associated with compaction of the chromosome. Upon bismuth exposure, H. pylori displays both defense and protection mechanisms, with bismuth extrusion by vesicles and shielding of the chromosome.

Chromosomal assembly of the flat oyster (Ostrea edulis L.) genome as a new genetic resource for aquaculture.

The European flat oyster (Ostrea edulis L.) is a native bivalve of the European coasts. Harvest of this species has declined during the last decades because of the appearance of two parasites that have led to the collapse of the stocks and the loss of the natural oyster beds. O. edulis has been the subject of numerous studies in population genetics and on the detection of the parasites Bonamia ostreae and Marteilia refringens. These studies investigated immune responses to these parasites at the molecular and cellular levels. Several genetic improvement programs have been initiated especially for parasite resistance. Within the framework of a European project (PERLE 2) that aims to produce genetic lines of O. edulis with hardiness traits (growth, survival, resistance) for the purpose of repopulating natural oyster beds in Brittany and reviving the culture of this species in the foreshore, obtaining a reference genome becomes essential as done recently in many bivalve species of aquaculture interest. Here, we present a chromosome-level genome assembly and annotation for the European flat oyster, generated by combining PacBio, Illumina, 10X linked, and Hi-C sequencing. The finished assembly is 887.2 Mb with a scaffold-N50 of 97.1 Mb scaffolded on the expected 10 pseudochromosomes. Annotation of the genome revealed the presence of 35,962 protein-coding genes. We analyzed in detail the transposable element (TE) diversity in the flat oyster genome, highlighted some specificities in tRNA and miRNA composition, and provided the first insight into the molecular response of O. edulis to M. refringens. This genome provides a reference for genomic studies on O. edulis to better understand its basic physiology and as a useful resource for genetic breeding in support of aquaculture and natural reef restoration.

Closed and High-Quality Bacterial Genome Sequences of the Oligo-Mouse-Microbiota Community.

The Oligo-Mouse-Microbiota (OMM12) gnotobiotic murine model is an increasingly popular model in microbiota studies. However, following Illumina and PacBio sequencing, the genomes of the 12 strains could not be closed. Here, we used genomic chromosome conformation capture (Hi-C) data to reorganize, close, and improve the quality of these 12 genomes.

MetaHiC phage-bacteria infection network reveals active cycling phages of the healthy human gut.

Bacteriophages play important roles in regulating the intestinal human microbiota composition, dynamics, and homeostasis, and characterizing their bacterial hosts is needed to understand their impact. We applied a metagenomic Hi-C approach on 10 healthy human gut samples to unveil a large infection network encompassing more than 6000 interactions bridging a metagenomic assembled genomes (MAGs) and a phage sequence, allowing to study in situ phage-host ratio. Whereas three-quarters of these sequences likely correspond to dormant prophages, 5% exhibit a much higher coverage than their associated MAG, representing potentially actively replicating phages. We detected 17 sequences of members of the crAss-like phage family, whose hosts diversity remained until recently relatively elusive. For each of them, a unique bacterial host was identified, all belonging to different genus of Bacteroidetes. Therefore, metaHiC deciphers infection network of microbial population with a high specificity paving the way to dynamic analysis of mobile genetic elements in complex ecosystems.

High-throughput identification of viral termini and packaging mechanisms in virome datasets using PhageTermVirome.

Viruses that infect bacteria (phages) are increasingly recognized for their importance in diverse ecosystems but identifying and annotating them in large-scale sequence datasets is still challenging. Although efficient scalable virus identification tools are emerging, defining the exact ends (termini) of phage genomes is still particularly difficult. The proper identification of termini is crucial, as it helps in characterizing the packaging mechanism of bacteriophages and provides information on various aspects of phage biology. Here, we introduce PhageTermVirome (PTV) as a tool for the easy and rapid high-throughput determination of phage termini and packaging mechanisms using modern large-scale metagenomics datasets. We successfully tested the PTV algorithm on a mock virome dataset and then used it on two real virome datasets to achieve the rapid identification of more than 100 phage termini and packaging mechanisms, with just a few hours of computing time. Because PTV allows the identification of free fully formed viral particles (by recognition of termini present only in encapsidated DNA), it can also complement other virus identification softwares to predict the true viral origin of contigs in viral metagenomics datasets. PTV is a novel and unique tool for high-throughput characterization of phage genomes, including phage termini identification and characterization of genome packaging mechanisms. This software should help researchers better visualize, map and study the virosphere. PTV is freely available for downloading and installation at https://gitlab.pasteur.fr/vlegrand/ptv .

Chromosome-level genome assembly reveals homologous chromosomes and recombination in asexual rotifer Adineta vaga.

Bdelloid rotifers are notorious as a speciose ancient clade comprising only asexual lineages. Thanks to their ability to repair highly fragmented DNA, most bdelloid species also withstand complete desiccation and ionizing radiation. Producing a well-assembled reference genome is a critical step to developing an understanding of the effects of long-term asexuality and DNA breakage on genome evolution. To this end, we present the first high-quality chromosome-level genome assemblies for the bdelloid Adineta vaga, composed of six pairs of homologous (diploid) chromosomes with a footprint of paleotetraploidy. The observed large-scale losses of heterozygosity are signatures of recombination between homologous chromosomes, either during mitotic DNA double-strand break repair or when resolving programmed DNA breaks during a modified meiosis. Dynamic subtelomeric regions harbor more structural diversity (e.g., chromosome rearrangements, transposable elements, and haplotypic divergence). Our results trigger the reappraisal of potential meiotic processes in bdelloid rotifers and help unravel the factors underlying their long-term asexual evolutionary success.

Chromosomal scale assembly of parasitic wasp genome reveals symbiotic virus colonization.

Endogenous viruses form an important proportion of eukaryote genomes and a source of novel functions. How large DNA viruses integrated into a genome evolve when they confer a benefit to their host, however, remains unknown. Bracoviruses are essential for the parasitism success of parasitoid wasps, into whose genomes they integrated ~103 million years ago. Here we show, from the assembly of a parasitoid wasp genome at a chromosomal scale, that bracovirus genes colonized all ten chromosomes of Cotesia congregata. Most form clusters of genes involved in particle production or parasitism success. Genomic comparison with another wasp, Microplitis demolitor, revealed that these clusters were already established ~53 mya and thus belong to remarkably stable genomic structures, the architectures of which are evolutionary constrained. Transcriptomic analyses highlight temporal synchronization of viral gene expression without resulting in immune gene induction, suggesting that no conflicts remain between ancient symbiotic partners when benefits to them converge.