Insights into early animal evolution form the genome of the xenacoelomorph worm Xenoturbella bocki.

The evolutionary origins of Bilateria remain enigmatic. One of the more enduring proposals highlights similarities between a cnidarian-like planula larva and simple acoel-like flatworms. This idea is based in part on the view of the Xenacoelomorpha as an outgroup to all other bilaterians which are themselves designated the Nephrozoa (protostomes and deuterostomes). Genome data can provide important comparative data and help to understand the evolution and biology of enigmatic species better. Here we assemble and analyse the genome of the simple, marine xenacoelomorph Xenoturbella bocki, a key species for our understanding of early bilaterian evolution. Our highly contiguous genome assembly of X. bocki has a size of ~111 Mbp in 18 chromosome like scaffolds, with repeat content and intron, exon and intergenic space comparable to other bilaterian invertebrates. We find X. bocki to have a similar number of genes to other bilaterians and to have retained ancestral metazoan synteny. Key bilaterian signalling pathways are also largely complete and most bilaterian miRNAs are present. Overall, we conclude that X. bocki has a complex genome typical of bilaterians, which does not reflect the apparent simplicity of its body plan that has been so important to proposals that the Xenacoelomorpha are the simple sister group of the rest of the Bilateria.

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.

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.

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.

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.

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.

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-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.

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 .

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.

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.

instaGRAAL: chromosome-level quality scaffolding of genomes using a proximity ligation-based scaffolder.

Hi-C exploits contact frequencies between pairs of loci to bridge and order contigs during genome assembly, resulting in chromosome-level assemblies. Because few robust programs are available for this type of data, we developed instaGRAAL, a complete overhaul of the GRAAL program, which has adapted the latter to allow efficient assembly of large genomes. instaGRAAL features a number of improvements over GRAAL, including a modular correction approach that optionally integrates independent data. We validate the program using data for two brown algae, and human, to generate near-complete assemblies with minimal human intervention.

Chromosome organization by a conserved condensin-ParB system in the actinobacterium Corynebacterium glutamicum.

Higher-order chromosome folding and segregation are tightly regulated in all domains of life. In bacteria, details on nucleoid organization regulatory mechanisms and function remain poorly characterized, especially in non-model species. Here, we investigate the role of DNA-partitioning protein ParB and SMC condensin complexes in the actinobacterium Corynebacterium glutamicum. Chromosome conformation capture reveals SMC-mediated long-range interactions around ten centromere-like parS sites clustered at the replication origin (oriC). At least one oriC-proximal parS site is necessary for reliable chromosome segregation. We use chromatin immunoprecipitation and photoactivated single-molecule localization microscopy to show the formation of distinct, parS-dependent ParB-nucleoprotein subclusters. We further show that SMC/ScpAB complexes, loaded via ParB at parS sites, mediate chromosomal inter-arm contacts (as previously shown in Bacillus subtilis). However, the MukBEF-like SMC complex MksBEFG does not contribute to chromosomal DNA-folding; instead, this complex is involved in plasmid maintenance and interacts with the polar oriC-tethering factor DivIVA. Our results complement current models of ParB-SMC/ScpAB crosstalk and show that some condensin complexes evolved functions that are apparently uncoupled from chromosome folding.

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.

Generation of a Metagenomics Proximity Ligation 3C Library of a Mammalian Gut Microbiota.

Microbial species thrive in very diverse environments and play fundamental roles in their equilibrium and dynamics. Metagenomics consists in extracting, sequencing, and studying the DNA present in ecosystems to better understand their regulation. Ideally, the maximal amount of information would be gathered from the full sequences of the genomes, episomes, and phages present in the microbial communities. Current high-throughput DNA sequencing produces reads ranging in size from a few dozen base pairs for the most commonly used technologies to several kb for emerging single-molecule real-time sequencing techniques. Although valuable information can be extracted from processing these DNA sequences into contigs, reconstructing full genomes remains a difficult task. Clustering contigs according to their similarities or read coverage covariations gives some insights on these genomes, but remains limited as viral sequences, or recent horizontal gene transfers, often differ from their host genomes. We recently developed meta3C, a proximity ligation approach that bins contigs in a sequence-independent way by quantifying and exploiting their tridimensional collisions frequencies in vivo. This technique has demonstrated a great potential to reconstruct genomes as well as to assign plasmids and phages to their hosts. It nevertheless requires a specific processing of the microbial samples before sequencing, which has to be carefully planned.

Tridimensional infiltration of DNA viruses into the host genome shows preferential contact with active chromatin.

Whether non-integrated viral DNAs distribute randomly or target specific positions within the higher-order architecture of mammalian genomes remains largely unknown. Here we use Hi-C and viral DNA capture (CHi-C) in primary human hepatocytes infected by either hepatitis B virus (HBV) or adenovirus type 5 (Ad5) virus to show that they adopt different strategies in their respective positioning at active chromatin. HBV contacts preferentially CpG islands (CGIs) enriched in Cfp1 a factor required for its transcription. These CGIs are often associated with highly expressed genes (HEG) and genes deregulated during infection. Ad5 DNA interacts preferentially with transcription start sites (TSSs) and enhancers of HEG, as well as genes upregulated during infection. These results show that DNA viruses use different strategies to infiltrate genomic 3D networks and target specific regions. This targeting may facilitate the recruitment of transcription factors necessary for their own replication and contribute to the deregulation of cellular gene expression.