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

Proxi-RIMS-seq2 applied to native microbiomes uncovers hundreds of known and novel m5 C methyltransferase specificities

Methylation patterns in bacteria can be used to study Restriction-Modification (RM) or other defense systems with novel properties. While m4 C and m6 A methylation is well characterized mainly through PacBio sequencing, the landscape of m5 C methylation is under-characterized. To bridge this gap, we performed RIMS-seq2 on microbiomes composed of resolved assemblies of distinct genomes through proximity ligation. This high-throughput approach enables the identification of m5 C methylated motifs and links them to cognate methyltransferases directly on native microbiomes without the need to isolate bacterial strains. Methylation patterns can also be identified on viral DNA and compared to host DNA, strengthening evidence for virus-host interaction. Applied to three different microbiomes, the method unveils over 1900 motifs that were deposited in REBASE. The motifs include a novel 8-base recognition site (CAT m5 CGATG) that was experimentally validated by characterizing its cognate methyltransferase. Our findings suggest that microbiomes harbor arrays of untapped m5 C methyltransferase specificities, providing insights to bacterial biology and biotechnological applications.

Antimicrobial resistance prediction by clinical metagenomics in pediatric severe pneumonia patients

BACKGROUND

Antimicrobial resistance (AMR) is a major threat to children's health, particularly in respiratory infections. Accurate identification of pathogens and AMR is crucial for targeted antibiotic treatment. Metagenomic next-generation sequencing (mNGS) shows promise in directly detecting microorganisms and resistance genes in clinical samples. However, the accuracy of AMR prediction through mNGS testing needs further investigation for practical clinical decision-making.

METHODS

We aimed to evaluate the performance of mNGS in predicting AMR for severe pneumonia in pediatric patients. We conducted a retrospective analysis at a tertiary hospital from May 2022 to May 2023. Simultaneous mNGS and culture were performed on bronchoalveolar lavage fluid samples obtained from pediatric patients with severe pneumonia. By comparing the results of mNGS detection of microorganisms and antibiotic resistance genes with those of culture, sensitivity, specificity, positive predictive value, and negative predictive value were calculated.

RESULTS

mNGS detected bacterial in 71.7% cases (86/120), significantly higher than culture (58/120, 48.3%). Compared to culture, mNGS demonstrated a sensitivity of 96.6% and a specificity of 51.6% in detecting pathogenic microorganisms. Phenotypic susceptibility testing (PST) of 19 antibiotics revealed significant variations in antibiotics resistance rates among different bacteria. Sensitivity prediction of mNGS for carbapenem resistance was higher than penicillins and cephalosporin (67.74% vs. 28.57%, 46.15%), while specificity showed no significant difference (85.71%, 75.00%, 75.00%). mNGS also showed a high sensitivity of 94.74% in predicting carbapenem resistance in Acinetobacter baumannii.

CONCLUSIONS

mNGS exhibits variable predictive performance among different pathogens and antibiotics, indicating its potential as a supplementary tool to conventional PST. However, mNGS currently cannot replace conventional PST.

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.

[Tracking transfers of resistance-carrying bacteria between animals, humans and the environment]

The fight against antibiotic resistance must incorporate the "One Health" concept to be effective. This means having a holistic approach embracing the different ecosystems, human, animal, and environment. Transfers of resistance genes may exist between these three domains and different stresses related to the exposome may influence these transfers. Various targeted or pan-genomic molecular biology techniques can be used to better characterise the dissemination of bacterial clones and to identify exchanges of genes and mobile genetic elements between ecosystems.

Orchestrating chromosome conformation capture analysis with Bioconductor

Genome-wide chromatin conformation capture assays provide formidable insights into the spatial organization of genomes. However, due to the complexity of the data structure, their integration in multi-omics workflows remains challenging. We present data structures, computational methods and visualization tools available in Bioconductor to investigate Hi-C, micro-C and other 3C-related data, in R. An online book ( https://bioconductor.org/books/OHCA/ ) further provides prospective end users with a number of workflows to process, import, analyze and visualize any type of chromosome conformation capture data.

Molecular Dynamics Simulation of a Feather-Boa Model of a Bacterial Chromosome

The chromosome of a bacterium consists of a mega-base pair-long circular DNA, which self-organizes within the micron-sized bacterial cell volume, compacting itself by three orders of magnitude. Unlike eukaryotic chromosomes, it lacks a nuclear membrane and freely floats in the cytosol confined by the cell membrane. It is believed that strong confinement, cross-linking by associated proteins, and molecular crowding all contribute to determine chromosome size and morphology. Modelling the chromosome simply as a circular polymer decorated with closed side loops in a cylindrical confining volume has been shown to already recapture some of the salient properties observed experimentally. Here we describe how a computer simulation can be set up to study structure and dynamics of bacterial chromosomes using this model.

Shedding Light on Bacterial Chromosome Structure: Exploring the Significance of 3C-Based Approaches

The complex architecture of DNA within living organisms is essential for maintaining the genetic information that dictates their functions and characteristics. Among the many complexities of genetic material organization, the folding and arrangement of DNA into chromosomes play a critical role in regulating gene expression, replication, and other essential cellular processes. Bacteria, despite their apparently simple cellular structure, exhibit a remarkable level of chromosomal organization that influences their adaptability and survival in diverse environments. Understanding the three-dimensional arrangement of bacterial chromosomes has long been a challenge due to technical limitations, but the development of Chromosome Conformation Capture (3C) methods revolutionized our ability to explore the hierarchical structure and the dynamics of bacterial genomes. Here, we review the major advances in the field of bacterial chromosome structure using 3C technology over the past decade.

The landscape of the histone-organized chromatin of Bdellovibrionota bacteria

Histone proteins have traditionally been thought to be restricted to eukaryotes and most archaea, with eukaryotic nucleosomal histones deriving from their archaeal ancestors. In contrast, bacteria lack histones as a rule. However, histone proteins have recently been identified in a few bacterial clades, most notably the phylum Bdellovibrionota, and these histones have been proposed to exhibit a range of divergent features compared to histones in archaea and eukaryotes. However, no functional genomic studies of the properties of Bdellovibrionota chromatin have been carried out. In this work, we map the landscape of chromatin accessibility, active transcription and three-dimensional genome organization in a member of Bdellovibrionota (a Bacteriovorax strain). We find that, similar to what is observed in some archaea and in eukaryotes with compact genomes such as yeast, Bacteriovorax chromatin is characterized by preferential accessibility around promoter regions. Similar to eukaryotes, chromatin accessibility in Bacteriovorax positively correlates with gene expression. Mapping active transcription through single-strand DNA (ssDNA) profiling revealed that unlike in yeast, but similar to the state of mammalian and fly promoters, Bacteriovorax promoters exhibit very strong polymerase pausing. Finally, similar to that of other bacteria without histones, the Bacteriovorax genome exists in a three-dimensional (3D) configuration organized by the parABS system along the axis defined by replication origin and termination regions. These results provide a foundation for understanding the chromatin biology of the unique Bdellovibrionota bacteria and the functional diversity in chromatin organization across the tree of life.

Insights in bacterial genome folding

Chromosomes in all domains of life are well-defined structural entities with complex hierarchical organization. The regulation of this hierarchical organization and its functional interplay with gene expression or other chromosome metabolic processes such as repair, replication, or segregation is actively investigated in a variety of species, including prokaryotes. Bacterial chromosomes are typically gene-dense with few non-coding sequences and are organized into the nucleoid, a membrane-less compartment composed of DNA, RNA, and proteins (nucleoid-associated proteins or NAPs). The continuous improvement of imaging and genomic methods has put the organization of these Mb-long molecules at reach, allowing to disambiguate some of their highly dynamic properties and intertwined structural features. Here we review and discuss some of the recent advances in the field of bacterial chromosome organization.

Hybridisation has shaped a recent radiation of grass-feeding aphids

BACKGROUND

Aphids are common crop pests. These insects reproduce by facultative parthenogenesis involving several rounds of clonal reproduction interspersed with an occasional sexual cycle. Furthermore, clonal aphids give birth to live young that are already pregnant. These qualities enable rapid population growth and have facilitated the colonisation of crops globally. In several cases, so-called "super clones" have come to dominate agricultural systems. However, the extent to which the sexual stage of the aphid life cycle has shaped global pest populations has remained unclear, as have the origins of successful lineages. Here, we used chromosome-scale genome assemblies to disentangle the evolution of two global pests of cereals-the English (Sitobion avenae) and Indian (Sitobion miscanthi) grain aphids.

RESULTS

Genome-wide divergence between S. avenae and S. miscanthi is low. Moreover, comparison of haplotype-resolved assemblies revealed that the S. miscanthi isolate used for genome sequencing is likely a hybrid, with one of its diploid genome copies closely related to S. avenae (~ 0.5% divergence) and the other substantially more divergent (> 1%). Population genomics analyses of UK and China grain aphids showed that S. avenae and S. miscanthi are part of a cryptic species complex with many highly differentiated lineages that predate the origins of agriculture. The complex consists of hybrid lineages that display a tangled history of hybridisation and genetic introgression.

CONCLUSIONS

Our analyses reveal that hybridisation has substantially contributed to grain aphid diversity, and hence, to the evolutionary potential of this important pest species. Furthermore, we propose that aphids are particularly well placed to exploit hybridisation events via the rapid propagation of live-born "frozen hybrids" via asexual reproduction, increasing the likelihood of hybrid lineage formation.

Noise reduction strategies in metagenomic chromosome confirmation capture to link antibiotic resistance genes to microbial hosts

The gut microbiota is a reservoir for antimicrobial resistance genes (ARGs). With current sequencing methods, it is difficult to assign ARGs to their microbial hosts, particularly if these ARGs are located on plasmids. Metagenomic chromosome conformation capture approaches (meta3C and Hi-C) have recently been developed to link bacterial genes to phylogenetic markers, thus potentially allowing the assignment of ARGs to their hosts on a microbiome-wide scale. Here, we generated a meta3C dataset of a human stool sample and used previously published meta3C and Hi-C datasets to investigate bacterial hosts of ARGs in the human gut microbiome. Sequence reads mapping to repetitive elements were found to cause problematic noise in, and may importantly skew interpretation of, meta3C and Hi-C data. We provide a strategy to improve the signal-to-noise ratio by discarding reads that map to insertion sequence elements and to the end of contigs. We also show the importance of using spike-in controls to quantify whether the cross-linking step in meta3C and Hi-C protocols has been successful. After filtering to remove artefactual links, 87 ARGs were assigned to their bacterial hosts across all datasets, including 27 ARGs in the meta3C dataset we generated. We show that commensal gut bacteria are an important reservoir for ARGs, with genes coding for aminoglycoside and tetracycline resistance being widespread in anaerobic commensals of the human gut.

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 OMM¹² synthetic bacterial community stably associated within mice gut (gnotobiotic mouse line OMM¹²). 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 OMM¹² 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.

Bacteriophages of the Order Crassvirales: What Do We Currently Know about This Keystone Component of the Human Gut Virome?

The order Crassvirales comprises dsDNA bacteriophages infecting bacteria in the phylum Bacteroidetes that are found in a variety of environments but are especially prevalent in the mammalian gut. This review summarises available information on the genomics, diversity, taxonomy, and ecology of this largely uncultured viral taxon. With experimental data available from a handful of cultured representatives, the review highlights key properties of virion morphology, infection, gene expression and replication processes, and phage-host dynamics.

Longitudinal study of the short- and long-term effects of hospitalisation and oral trimethoprim-sulfadiazine administration on the equine faecal microbiome and resistome

BACKGROUND

Hospitalisation and antimicrobial treatment are common in horses and significantly impact the intestinal microbiota. Antimicrobial treatment might also increase levels of resistant bacteria in faeces, which could spread to other ecological compartments, such as the environment, other animals and humans. In this study, we aimed to characterise the short- and long-term effects of transportation, hospitalisation and trimethoprim-sulfadiazine (TMS) administration on the faecal microbiota and resistome of healthy equids.

METHODS

In a longitudinal experimental study design, in which the ponies served as their own control, faecal samples were collected from six healthy Welsh ponies at the farm (D0-D13-1), immediately following transportation to the hospital (D13-2), during 7 days of hospitalisation without treatment (D14-D21), during 5 days of oral TMS treatment (D22-D26) and after discharge from the hospital up to 6 months later (D27-D211). After DNA extraction, 16S rRNA gene sequencing was performed on all samples. For resistome analysis, shotgun metagenomic sequencing was performed on selected samples.

RESULTS

Hospitalisation without antimicrobial treatment did not significantly affect microbiota composition. Oral TMS treatment reduced alpha-diversity significantly. Kiritimatiellaeota, Fibrobacteres and Verrucomicrobia significantly decreased in relative abundance, whereas Firmicutes increased. The faecal microbiota composition gradually recovered after discontinuation of TMS treatment and discharge from the hospital and, after 2 weeks, was more similar to pre-treatment composition than to composition during TMS treatment. Six months later, however, microbiota composition still differed significantly from that at the start of the study and Spirochaetes and Verrucomicrobia were less abundant. TMS administration led to a significant (up to 32-fold) and rapid increase in the relative abundance of resistance genes sul2, tetQ, ant6-1a, and aph(3")-lb. lnuC significantly decreased directly after treatment. Resistance genes sul2 (15-fold) and tetQ (six-fold) remained significantly increased 6 months later.

CONCLUSIONS

Oral treatment with TMS has a rapid and long-lasting effect on faecal microbiota composition and resistome, making the equine hindgut a reservoir and potential source of resistant bacteria posing a risk to animal and human health through transmission. These findings support the judicious use of antimicrobials to minimise long-term faecal presence, excretion and the spread of antimicrobial resistance in the environment. Video Abstract.

ViralCC retrieves complete viral genomes and virus-host pairs from metagenomic Hi-C data

The introduction of high-throughput chromosome conformation capture (Hi-C) into metagenomics enables reconstructing high-quality metagenome-assembled genomes (MAGs) from microbial communities. Despite recent advances in recovering eukaryotic, bacterial, and archaeal genomes using Hi-C contact maps, few of Hi-C-based methods are designed to retrieve viral genomes. Here we introduce ViralCC, a publicly available tool to recover complete viral genomes and detect virus-host pairs using Hi-C data. Compared to other Hi-C-based methods, ViralCC leverages the virus-host proximity structure as a complementary information source for the Hi-C interactions. Using mock and real metagenomic Hi-C datasets from several different microbial ecosystems, including the human gut, cow fecal, and wastewater, we demonstrate that ViralCC outperforms existing Hi-C-based binning methods as well as state-of-the-art tools specifically dedicated to metagenomic viral binning. ViralCC can also reveal the taxonomic structure of viruses and virus-host pairs in microbial communities. When applied to a real wastewater metagenomic Hi-C dataset, ViralCC constructs a phage-host network, which is further validated using CRISPR spacer analyses. ViralCC is an open-source pipeline available at https://github.com/dyxstat/ViralCC .

Hypothesis-driven probabilistic modelling enables a principled perspective of genomic compartments

The Hi-C method has revolutionized the study of genome organization, yet interpretation of Hi-C interaction frequency maps remains a major challenge. Genomic compartments are a checkered Hi-C interaction pattern suggested to represent the partitioning of the genome into two self-interacting states associated with active and inactive chromatin. Based on a few elementary mechanistic assumptions, we derive a generative probabilistic model of genomic compartments, called deGeco. Testing our model, we find it can explain observed Hi-C interaction maps in a highly robust manner, allowing accurate inference of interaction probability maps from extremely sparse data without any training of parameters. Taking advantage of the interpretability of the model parameters, we then test hypotheses regarding the nature of genomic compartments. We find clear evidence of multiple states, and that these states self-interact with different affinities. We also find that the interaction rules of chromatin states differ considerably within and between chromosomes. Inspecting the molecular underpinnings of a four-state model, we show that a simple classifier can use histone marks to predict the underlying states with 87% accuracy. Finally, we observe instances of mixed-state loci and analyze these loci in single-cell Hi-C maps, finding that mixing of states occurs mainly at the cell level.

Diversity of novel archaeal viruses infecting methanogens discovered through coupling of stable isotope probing and metagenomics

Diversity of viruses infecting non-extremophilic archaea has been grossly understudied. This is particularly the case for viruses infecting methanogenic archaea, key players in the global carbon biogeochemical cycle. Only a dozen of methanogenic archaeal viruses have been isolated so far. In the present study, we implemented an original coupling between stable isotope probing and complementary shotgun metagenomic analyses to identify viruses of methanogens involved in the bioconversion of formate, which was used as the sole carbon source in batch anaerobic digestion microcosms. Under our experimental conditions, the microcosms were dominated by methanogens belonging to the order Methanobacteriales (Methanobacterium and Methanobrevibacter genera). Metagenomic analyses yielded several previously uncharacterized viral genomes, including a complete genome of a head-tailed virus (class Caudoviricetes, proposed family Speroviridae, Methanobacterium host) and several near-complete genomes of spindle-shaped viruses. The two groups of viruses are predicted to infect methanogens of the Methanobacterium and Methanosarcina genera and represent two new virus families. The metagenomics results are in good agreement with the electron microscopy observations, which revealed the dominance of head-tailed virus-like particles and the presence of spindle-shaped particles. The present study significantly expands the knowledge on the viral diversity of viruses of methanogens.

Inferring antibiotic susceptibility from metagenomic data: dream or reality?

BACKGROUND

The diagnosis of bacterial infections continues to rely on culture, a slow process in which antibiotic susceptibility profiles of potential pathogens are made available to clinicians 48h after sampling, at best. Recently, clinical metagenomics (CMg), the metagenomic sequencing of samples with the purpose of identifying microorganisms and determining their susceptibility to antimicrobials, has emerged as a potential diagnostic tool that could prove faster than culture. CMg indeed has the potential to detect antibiotic resistance genes (ARGs) and mutations associated with resistance. Nevertheless, many challenges have yet to be overcome in order to make rapid phenotypic inference of antibiotic susceptibility from metagenomic data a reality.

OBJECTIVES

The objective of this narrative review is to discuss the challenges underlying the phenotypic inference of antibiotic susceptibility from metagenomic data.

SOURCES

We conducted a narrative review using published articles available in the NCBI Pubmed database.

CONTENT

We review the current ARG databases with a specific emphasis on those which now provide associations with phenotypic data. Next, we discuss the bioinformatic tools designed to identify ARGs in metagenomes. We then report on the performance of phenotypic inference from genomic data and the issue predicting the expression of ARGs. Finally, we address the challenge of linking an ARG to this host.

IMPLICATIONS

Significant improvements have recently been made in associating ARG and phenotype, and the inference of susceptibility from genomic data has been demonstrated in pathogenic bacteria such as Staphylococci and Enterobacterales. Resistance involving gene expression is more challenging however, and inferring susceptibility from species such as Pseudomonas aeruginosa remains difficult. Future research directions include the consideration of gene expression via RNA sequencing and machine learning.

HAM-ART: An optimised culture-free Hi-C metagenomics pipeline for tracking antimicrobial resistance genes in complex microbial communities

Shotgun metagenomics is a powerful tool to identify antimicrobial resistance (AMR) genes in microbiomes but has the limitation that extrachromosomal DNA, such as plasmids, cannot be linked with the host bacterial chromosome. Here we present a comprehensive laboratory and bioinformatics pipeline HAM-ART (Hi-C Assisted Metagenomics for Antimicrobial Resistance Tracking) optimised for the generation of metagenome-assembled genomes including both chromosomal and extrachromosomal AMR genes. We demonstrate the performance of the pipeline in a study comparing 100 pig faecal microbiomes from low- and high-antimicrobial use pig farms (organic and conventional farms). We found significant differences in the distribution of AMR genes between low- and high-antimicrobial use farms including a plasmid-borne lincosamide resistance gene exclusive to high-antimicrobial use farms in three species of Lactobacilli. The bioinformatics pipeline code is available at https://github.com/lkalmar/HAM-ART.

HiCBin: binning metagenomic contigs and recovering metagenome-assembled genomes using Hi-C contact maps

Recovering high-quality metagenome-assembled genomes (MAGs) from complex microbial ecosystems remains challenging. Recently, high-throughput chromosome conformation capture (Hi-C) has been applied to simultaneously study multiple genomes in natural microbial communities. We develop HiCBin, a novel open-source pipeline, to resolve high-quality MAGs utilizing Hi-C contact maps. HiCBin employs the HiCzin normalization method and the Leiden clustering algorithm and includes the spurious contact detection into binning pipelines for the first time. HiCBin is validated on one synthetic and two real metagenomic samples and is shown to outperform the existing Hi-C-based binning methods. HiCBin is available at https://github.com/dyxstat/HiCBin.

Diversity and Evolution of Mamiellophyceae: Early-Diverging Phytoplanktonic Green Algae Containing Many Cosmopolitan Species

The genomic revolution has bridged a gap in our knowledge about the diversity, biology and evolution of unicellular photosynthetic eukaryotes, which bear very few discriminating morphological features among species from the same genus. The high-quality genome resources available in the class Mamiellophyceae (Chlorophyta) have been paramount to estimate species diversity and screen available metagenomic data to assess the biogeography and ecological niches of different species on a global scale. Here we review the current knowledge about the diversity, ecology and evolution of the Mamiellophyceae and the large double-stranded DNA prasinoviruses infecting them, brought by the combination of genomic and metagenomic analyses, including 26 metabarcoding environmental studies, as well as the pan-oceanic GOS and the Tara Oceans expeditions.

What makes a megaplasmid?

Naturally occurring plasmids come in different sizes. The smallest are less than a kilobase of DNA, while the largest can be over three orders of magnitude larger. Historically, research has tended to focus on smaller plasmids that are usually easier to isolate, manipulate and sequence, but with improved genome assemblies made possible by long-read sequencing, there is increased appreciation that very large plasmids—known as megaplasmids—are widespread, diverse, complex, and often encode key traits in the biology of their host microorganisms. Why are megaplasmids so big? What other features come with large plasmid size that could affect bacterial ecology and evolution? Are megaplasmids 'just' big plasmids, or do they have distinct characteristics? In this perspective, we reflect on the distribution, diversity, biology, and gene content of megaplasmids, providing an overview to these large, yet often overlooked, mobile genetic elements.

This article is part of the theme issue ‘The secret lives of microbial mobile genetic elements’.

The 3D Organization of Chromatin Colors in Mammalian Nuclei

While many computational methods have been proposed for 3D chromosome reconstruction from chromosomal contact maps, these methods are rarely used for the interpretation of such experimental data, in particular Hi-C data. We posit that this is due to the lack of an easy-to-use implementation of the proposed algorithms, as well as to the important computational cost of most methods. We here give a detailed implementation of the fast ShRec3D algorithm. We provide a tutorial that will enable the reader to reconstruct 3D consensus structures for human chromosomes and to decorate these structures with chromatin epigenetic states. We use this methodology to show that the bivalent chromatin, including Polycomb-rich domains, is spatially segregated and located in between the active and the quiescent chromatin compartments.

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.

Normalization of Chromosome Contact Maps: Matrix Balancing and Visualization

Over the last decade, genomic proximity ligation approaches have reshaped our vision of chromosomes 3D organizations, from bacteria nucleoids to larger eukaryotic genomes. The different protocols (3Cseq, Hi-C, TCC, MicroC [XL], Hi-CO, etc.) rely on common steps (chemical fixation digestion, ligation…) to detect pairs of genomic positions in close proximity. The most common way to represent these data is a matrix, or contact map, which allows visualizing the different chromatin structures (compartments, loops, etc.) that can be associated to other signals such as transcription, protein occupancy, etc. as well as, in some instances, to biological functions.In this chapter we present and discuss the filtering of the events recovered in proximity ligation experiments as well as the application of the balancing normalization procedure on the resulting contact map. We also describe a computational tool for visualizing normalized contact data dubbed Scalogram.The different processes described here are illustrated and supported by the laboratory custom-made scripts pooled into "hicstuff," an open-access python package accessible on github ( https://github.com/koszullab/hicstuff ).

Lysis-Hi-C as a method to study polymicrobial communities and eDNA

Microbes interact in natural communities in a spatially structured manner, particularly in biofilms and polymicrobial infections. While next generation sequencing approaches provide powerful insights into diversity, metabolic capacity, and mutational profiles of these communities, they generally fail to recover in situ spatial proximity between distinct genotypes in the interactome. Hi‐C is a promising method that has assisted in analysing complex microbiomes, by creating chromatin cross‐links in cells, that aid in identifying adjacent DNA, to improve de novo assembly. This study explored a modified Hi‐C approach involving an initial lysis phase prior to DNA cross‐linking, to test whether adjacent cell chromatin can be cross‐linked, anticipating that this could provide a new avenue for study of spatial‐mutational dynamics in structured microbial communities. An artificial polymicrobial mixture of Pseudomonas aeruginosa, Staphylococcus aureus, and Escherichia coli was lysed for 1–18 h, then prepared for Hi‐C. A murine biofilm infection model was treated with sonication, mechanical lysis, or chemical lysis before Hi‐C. Bioinformatic analyses of resulting Hi‐C interspecies chromatin links showed that while microbial species differed from one another, generally lysis significantly increased links between species and increased the distance of Hi‐C links within species, while also increasing novel plasmid‐chromosome links. The success of this modified lysis‐Hi‐C protocol in creating extracellular DNA links is a promising first step toward a new lysis‐Hi‐C based method to recover genotypic microgeography in polymicrobial communities, with potential future applications in diseases with localized resistance, such as cystic fibrosis lung infections and chronic diabetic ulcers.

Prospects for multi-omics in the microbial ecology of water engineering

Advances in high-throughput sequencing technologies and bioinformatics approaches over almost the last three decades have substantially increased our ability to explore microorganisms and their functions - including those that have yet to be cultivated in pure isolation. Genome-resolved metagenomic approaches have enabled linking powerful functional predictions to specific taxonomical groups with increasing fidelity. Additionally, related developments in both whole community gene expression surveys and metabolite profiling have permitted for direct surveys of community-scale functions in specific environmental settings. These advances have allowed for a shift in microbiome science away from descriptive studies and towards mechanistic and predictive frameworks for designing and harnessing microbial communities for desired beneficial outcomes. Water engineers, microbiologists, and microbial ecologists studying activated sludge, anaerobic digestion, and drinking water distribution systems have applied various (meta)omics techniques for connecting microbial community dynamics and physiologies to overall process parameters and system performance. However, the rapid pace at which new omics-based approaches are developed can appear daunting to those looking to apply these state-of-the-art practices for the first time. Here, we review how modern genome-resolved metagenomic approaches have been applied to a variety of water engineering applications from lab-scale bioreactors to full-scale systems. We describe integrated omics analysis across engineered water systems and the foundations for pairing these insights with modeling approaches. Lastly, we summarize emerging omics-based technologies that we believe will be powerful tools for water engineering applications. Overall, we provide a framework for microbial ecologists specializing in water engineering to apply cutting-edge omics approaches to their research questions to achieve novel functional insights. Successful adoption of predictive frameworks in engineered water systems could enable more economically and environmentally sustainable bioprocesses as demand for water and energy resources increases.

Generating High-Resolution Hi-C Contact Maps of Bacteria

During the past decade, Chromosome Conformation Capture (3C/Hi-C)-based methods have been used to probe the 3D structure and organization of bacterial genomes, revealing fundamental aspects of chromosome dynamics. However, the current protocols are expensive, inefficient, and limited in their resolution. Here we present a simple, cost-effective Hi-C approach that is readily applicable to a range of Gram-positive and Gram-negative bacteria.

Genome reconstruction and haplotype phasing using chromosome conformation capture methodologies

Genomic analysis of individuals or organisms is predicated on the availability of high-quality reference and genotype information. With the rapidly dropping costs of high-throughput DNA sequencing, this is becoming readily available for diverse organisms and for increasingly large populations of individuals. Despite these advances, there are still aspects of genome sequencing that remain challenging for existing sequencing methods. This includes the generation of long-range contiguity during genome assembly, identification of structural variants in both germline and somatic tissues, the phasing of haplotypes in diploid organisms and the resolution of genome sequence for organisms derived from complex samples. These types of information are valuable for understanding the role of genome sequence and genetic variation on genome function, and numerous approaches have been developed to address them. Recently, chromosome conformation capture (3C) experiments, such as the Hi-C assay, have emerged as powerful tools to aid in these challenges for genome reconstruction. We will review the current use of Hi-C as a tool for aiding in genome sequencing, addressing the applications, strengths, limitations and potential future directions for the use of 3C data in genome analysis. We argue that unique features of Hi-C experiments make this data type a powerful tool to address challenges in genome sequencing, and that future integration of Hi-C data with alternative sequencing assays will facilitate the continuing revolution in genomic analysis and genome sequencing.

Examining horizontal gene transfer in microbial communities

Bacteria acquire novel DNA through horizontal gene transfer (HGT), a process that enables an organism to rapidly adapt to changing environmental conditions, provides a competitive edge and potentially alters its relationship with its host. Although the HGT process is routinely exploited in laboratories, there is a surprising disconnect between what we know from laboratory experiments and what we know from natural environments, such as the human gut microbiome. Owing to a suite of newly available computational algorithms and experimental approaches, we have a broader understanding of the genes that are being transferred and are starting to understand the ecology of HGT in natural microbial communities. This Review focuses on these new technologies, the questions they can address and their limitations. As these methods are applied more broadly, we are beginning to recognize the full extent of HGT possible within a microbiome and the punctuated dynamics of HGT, specifically in response to external stimuli. Furthermore, we are better characterizing the complex selective pressures on mobile genetic elements and the mechanisms by which they interact with the bacterial host genome.

Interconnecting solvent quality, transcription, and chromosome folding in Escherichia coli

All cells fold their genomes, including bacterial cells, where the chromosome is compacted into a domain-organized meshwork called the nucleoid. How compaction and domain organization arise is not fully understood. Here, we describe a method to estimate the average mesh size of the nucleoid in Escherichia coli. Using nucleoid mesh size and DNA concentration estimates, we find that the cytoplasm behaves as a poor solvent for the chromosome when the cell is considered as a simple semidilute polymer solution. Monte Carlo simulations suggest that a poor solvent leads to chromosome compaction and DNA density heterogeneity (i.e., domain formation) at physiological DNA concentration. Fluorescence microscopy reveals that the heterogeneous DNA density negatively correlates with ribosome density within the nucleoid, consistent with cryoelectron tomography data. Drug experiments, together with past observations, suggest the hypothesis that RNAs contribute to the poor solvent effects, connecting chromosome compaction and domain formation to transcription and intracellular organization.

Generation of gene-level resolution chromosome contact maps in bacteria and archaea

Chromosome conformation capture (Hi-C) has become a routine method for probing the 3D organization of genomes. However, when applied to bacteria and archaea, current protocols are expensive and limited in their resolution. By dissecting the different steps of published eukaryotic and prokaryotic Hi-C protocols, we have developed a cost- and time-effective approach to generate high-resolution (down to 500 bp – 1 kb) contact matrices of both bacteria and archaea genomes.

For complete details on the use and execution of this protocol, please refer to Cockram et al. (2020).

•

Optimized Hi-C protocol for archaeal and bacterial genomes

•

Generation of genome-wide contact maps up to 1 kb resolution

•

Detailed description of steps from cell fixation to sequencing library preparation

•

A cost- and time-effective approach offering gene-level resolution contact maps

Discovering multiple types of DNA methylation from bacteria and microbiome using nanopore sequencing

Bacterial DNA methylation occurs at diverse sequence contexts and plays important functional roles in cellular defense and gene regulation. Existing methods for detecting DNA modification from nanopore sequencing data do not effectively support de novo study of unknown bacterial methylomes. In this work, we observed that a nanopore sequencing signal displays complex heterogeneity across methylation events of the same type. To enable nanopore sequencing for broadly applicable methylation discovery, we generated a training dataset from an assortment of bacterial species and developed a method, named nanodisco ( https://github.com/fanglab/nanodisco ), that couples the identification and fine mapping of the three forms of methylation into a multi-label classification framework. We applied it to individual bacteria and the mouse gut microbiome for reliable methylation discovery. In addition, we demonstrated the use of DNA methylation for binning metagenomic contigs, associating mobile genetic elements with their host genomes and identifying misassembled metagenomic contigs.

Recovery of complete genomes and non-chromosomal replicons from activated sludge enrichment microbial communities with long read metagenome sequencing

New long read sequencing technologies offer huge potential for effective recovery of complete, closed genomes from complex microbial communities. Using long read data (ONT MinION) obtained from an ensemble of activated sludge enrichment bioreactors we recover 22 closed or complete genomes of community members, including several species known to play key functional roles in wastewater bioprocesses, specifically microbes known to exhibit the polyphosphate- and glycogen-accumulating organism phenotypes (namely Candidatus Accumulibacter and Dechloromonas, and Micropruina, Defluviicoccus and Candidatus Contendobacter, respectively), and filamentous bacteria (Thiothrix) associated with the formation and stability of activated sludge flocs. Additionally we demonstrate the recovery of close to 100 circularised plasmids, phages and small microbial genomes from these microbial communities using long read assembled sequence. We describe methods for validating long read assembled genomes using their counterpart short read metagenome-assembled genomes, and assess the influence of different correction procedures on genome quality and predicted gene quality. Our findings establish the feasibility of performing long read metagenome-assembled genome recovery for both chromosomal and non-chromosomal replicons, and demonstrate the value of parallel sampling of moderately complex enrichment communities to obtaining high quality reference genomes of key functional species relevant for wastewater bioprocesses.

Approaches for characterizing and tracking hospital-associated multidrug-resistant bacteria

Hospital-associated infections are a major concern for global public health. Infections with antibiotic-resistant pathogens can cause empiric treatment failure, and for infections with multidrug-resistant bacteria which can overcome antibiotics of "last resort" there exists no alternative treatments. Despite extensive sanitization protocols, the hospital environment is a potent reservoir and vector of antibiotic-resistant organisms. Pathogens can persist on hospital surfaces and plumbing for months to years, acquire new antibiotic resistance genes by horizontal gene transfer, and initiate outbreaks of hospital-associated infections by spreading to patients via healthcare workers and visitors. Advancements in next-generation sequencing of bacterial genomes and metagenomes have expanded our ability to (1) identify species and track distinct strains, (2) comprehensively profile antibiotic resistance genes, and (3) resolve the mobile elements that facilitate intra- and intercellular gene transfer. This information can, in turn, be used to characterize the population dynamics of hospital-associated microbiota, track outbreaks to their environmental reservoirs, and inform future interventions. This review provides a detailed overview of the approaches and bioinformatic tools available to study isolates and metagenomes of hospital-associated bacteria, and their multi-layered networks of transmission.

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.

Essential Amino Acid Enrichment and Positive Selection Highlight Endosymbiont's Role in a Global Virus-Vectoring Pest

Host-associated microbes display remarkable convergence in genome repertoire resulting from selection to supplement missing host functions. Nutritional supplementation has been proposed in the verrucomicrobial endosymbiont Xiphinematobacter sp., which lives within a globally widespread group of plant-parasitic nematodes that vector damaging nepoviruses to plants. Only one genome sequence has been published from this symbiont, leaving unanswered questions about its diversity, host range, role, and selective pressures within its hosts. Because its hosts are exceptionally resistant to culturing, this symbiont is best studied through advanced genomic approaches. To analyze the role of Xiphinematobacter sp. in its host, sequencing was performed on nematode communities, and then genomes were extracted for comparative genomics, gene ontology enrichment tests, polymorphism analysis, de Bruijn-based genome-wide association studies, and tests of pathway- and site-specific selection on genes predicted play a role in the symbiosis. Results showed a closely clustered set of Xiphinematobacter isolates with reduced genomes of ∼917 kbp, for which a new species was proposed. Symbionts shared only 2.3% of genes with outgroup Verrucomicrobia, but comparative analyses showed high conservation of all 10 essential amino acid (EAA) biosynthesis pathways plus several vitamin pathways. These findings were supported by gene ontology enrichment tests and high polymorphisms in these pathways compared with background. Genome-wide association analysis confirmed high between-species fixation of alleles with significant functional enrichment for EAA and thiamine synthesis. Strong positive selection was detected on sites within these pathways, despite several being under increased purifying selection. Together, these results suggest that supplementation of EAAs missing in the host diet may drive this widespread symbiosis.IMPORTANCE Xiphinematobacter spp. are distinctly evolved intracellular symbionts in the phylum Verrucomicrobia, which includes the important human gut-associated microbe Akkermansia muciniphila and many highly abundant free-living soil microbes. Like Akkermansia sp., Xiphinematobacter sp. is obligately associated with the gut of its hosts, which in this case consists of a group of plant-parasitic nematodes that are among the top 10 most destructive species to global agriculture, by vectoring plant viruses. This study examined the hypothesis that the key to this symbiont's stable evolutionary association with its host is through provisioning nutrients that its host cannot make that may be lacking in the nematode's plant phloem diet, such as essential amino acids and several vitamins. The significance of our research is in demonstrating, using population genomics, the signatures of selective pressure on these hypothesized roles to ultimately learn how this independently evolved symbiont functionally mirrors symbionts of phloem-feeding insects.

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.

Euryarchaeal genomes are folded into SMC-dependent loops and domains, but lack transcription-mediated compartmentalization

Hi-C has become a routine method for probing the 3D organization of genomes. However, when applied to prokaryotes and archaea, the current protocols are expensive and limited in their resolution. We develop a cost-effective Hi-C protocol to explore chromosome conformations of these two kingdoms at the gene or operon level. We first validate it on E. coli and V. cholera, generating sub-kilobase-resolution contact maps, and then apply it to the euryarchaeota H. volcanii, Hbt. salinarum, and T. kodakaraensis. With a resolution of up to 1 kb, we explore the diversity of chromosome folding in this phylum. In contrast to crenarchaeota, these euryarchaeota lack (active/inactive) compartment-like structures. Instead, their genomes are composed of self-interacting domains and chromatin loops. In H. volcanii, these structures are regulated by transcription and the archaeal structural maintenance of chromosomes (SMC) protein, further supporting the ubiquitous role of these processes in shaping the higher-order organization of genomes.

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.

A Chromosome-Scale Genome Assembly for the Fusarium oxysporum Strain Fo5176 To Establish a Model Arabidopsis-Fungal Pathosystem

Plant pathogens cause widespread yield losses in agriculture. Understanding the drivers of plant-pathogen interactions requires decoding the molecular dialog leading to either resistance or disease. However, progress in deciphering pathogenicity genes has been severely hampered by suitable model systems and incomplete fungal genome assemblies. Here, we report a significant improvement of the assembly and annotation of the genome of the Fusarium oxysporum (Fo) strain Fo5176. Fo comprises a large number of serious plant pathogens on dozens of plant species with largely unresolved pathogenicity factors. The strain Fo5176 infects Arabidopsis thaliana and, hence, constitutes a highly promising model system. We use high-coverage Pacific Biosciences Sequel long-read and Hi-C sequencing data to assemble the genome into 19 chromosomes and a total genome size of 67.98 Mb. The genome has a N50 of 4 Mb and a 99.1% complete BUSCO score. Phylogenomic analyses based on single-copy orthologs clearly place the Fo5176 strain in the Fo f sp. conglutinans clade as expected. We generated RNAseq data from culture medium and plant infections to train gene predictions and identified ∼18,000 genes including ten effector genes known from other Fo clades. We show that Fo5176 is able to infect cabbage and Brussel sprouts of the Brassica oleracea, expanding the usefulness of the Fo5176 model pathosystem. Finally, we performed large-scale comparative genomics analyses comparing the Fo5176 to 103 additional Fo genomes to define core and accessory genomic regions. In conjunction with the molecular tool sets available for A. thaliana, the Fo5176 genome and annotation provides a crucial step toward the establishment of a highly promising pathosystem.

Recent advances in understanding the ecology of the lung microbiota and deciphering the gut-lung axis

A rapidly expanding new field of lung research has been produced by the emergence of culture-independent next-generation sequencing technologies. While pulmonary microbiome research lags behind the exploration of the microbiome in other organ systems, the field is maturing and has recently produced multiple exciting discoveries. In this mini-review, we will explore recent advances in our understanding of the lung microbiome and the gut-lung axis from an ecological perspective.

Recent advances in understanding the ecology of the lung microbiota and deciphering the gut-lung axis

A rapidly expanding new field of lung research has been produced by the emergence of culture-independent next-generation sequencing technologies. While pulmonary microbiome research lags behind the exploration of the microbiome in other organ systems, the field is maturing and has recently produced multiple exciting discoveries. In this mini-review, we will explore recent advances in our understanding of the lung microbiome and the gut-lung axis from an ecological perspective.

Chromosome-Scale Genome Assembly of Fusarium oxysporum Strain Fo47, a Fungal Endophyte and Biocontrol Agent

Here, we report a chromosome-level genome assembly of Fusarium oxysporum Fo47 (12 pseudomolecules; contig N₅₀: 4.52 Mb), generated using a combination of PacBio long-read, Illumina paired end, and high-throughput chromosome conformation capture sequencing data. Although F. oxysporum causes vascular wilt to over 100 plant species, the strain Fo47 is classified as an endophyte and is widely used as a biocontrol agent for plant disease control. The Fo47 genome carries a single accessory chromosome of 4.23 Mb, compared with the reference genome of F. oxysporum f. sp. lycopersici Fol4287. The high-quality assembly and annotation of the Fo47 genome will be a valuable resource for studying the mechanisms underlying the endophytic interactions between F. oxysporum and plants as well as for deciphering the genome evolution of the F. oxysporum species complex.

Experimental approaches to tracking mobile genetic elements in microbial communities

Horizontal gene transfer is an important mechanism of microbial evolution and is often driven by the movement of mobile genetic elements between cells. Due to the fact that microbes live within communities, various mechanisms of horizontal gene transfer and types of mobile elements can co-occur. However, the ways in which horizontal gene transfer impacts and is impacted by communities containing diverse mobile elements has been challenging to address. Thus, the field would benefit from incorporating community-level information and novel approaches alongside existing methods. Emerging technologies for tracking mobile elements and assigning them to host organisms provide promise for understanding the web of potential DNA transfers in diverse microbial communities more comprehensively. Compared to existing experimental approaches, chromosome conformation capture and methylome analyses have the potential to simultaneously study various types of mobile elements and their associated hosts. We also briefly discuss how fermented food microbiomes, given their experimental tractability and moderate species complexity, make ideal models to which to apply the techniques discussed herein and how they can be used to address outstanding questions in the field of horizontal gene transfer in microbial communities.

Widespread transfer of mobile antibiotic resistance genes within individual gut microbiomes revealed through bacterial Hi-C

The gut microbiome harbors a ‘silent reservoir’ of antibiotic resistance (AR) genes that is thought to contribute to the emergence of multidrug-resistant pathogens through horizontal gene transfer (HGT). To counteract the spread of AR, it is paramount to know which organisms harbor mobile AR genes and which organisms engage in HGT. Despite methods that characterize the overall abundance of AR genes in the gut, technological limitations of short-read sequencing have precluded linking bacterial taxa to specific mobile genetic elements (MGEs) encoding AR genes. Here, we apply Hi-C, a high-throughput, culture-independent method, to surveil the bacterial carriage of MGEs. We compare two healthy individuals with seven neutropenic patients undergoing hematopoietic stem cell transplantation, who receive multiple courses of antibiotics, and are acutely vulnerable to the threat of multidrug-resistant infections. We find distinct networks of HGT across individuals, though AR and mobile genes are associated with more diverse taxa within the neutropenic patients than the healthy subjects. Our data further suggest that HGT occurs frequently over a several-week period in both cohorts. Whereas most efforts to understand the spread of AR genes have focused on pathogenic species, our findings shed light on the role of the human gut microbiome in this process.

Linking antibiotic resistance (AR) in the gut microbiome with their bacterial hosts remains challenging. Here, the authors apply bacterial Hi-C to map mobile genetic elements in metagenomes, and illustrate that genes are present in more diverse taxa in neutropenic patients than healthy subjects.

Exploring neighborhoods in large metagenome assembly graphs using spacegraphcats reveals hidden sequence diversity

Genomes computationally inferred from large metagenomic data sets are often incomplete and may be missing functionally important content and strain variation. We introduce an information retrieval system for large metagenomic data sets that exploits the sparsity of DNA assembly graphs to efficiently extract subgraphs surrounding an inferred genome. We apply this system to recover missing content from genome bins and show that substantial genomic sequence variation is present in a real metagenome. Our software implementation is available at https://github.com/spacegraphcats/spacegraphcats under the 3-Clause BSD License.

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.