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

Cutting-edge tools for unveiling the dynamics of plasmid-host interactions

The plasmid-mediated transfer of antibiotic resistance genes (ARGs) in complex microbiomes presents a significant global health challenge. This review examines recent technological advancements that have enabled us to move beyond the limitations of culture-dependent detection of conjugation and have enhanced our ability to track and understand the movement of ARGs in real-world scenarios. We critically assess the applications of single-cell sequencing, fluorescence-based techniques and advanced high-throughput chromatin conformation capture (Hi-C) approaches in elucidating plasmid-host interactions at unprecedented resolution. We also evaluate emerging techniques such as CRISPR-based phage engineering and discuss their potential for developing targeted strategies to curb ARG dissemination. Emerging data derived from these technologies have challenged our previous paradigms on plasmid-host compatibility and an awareness of an emerging uncharted realm for ARGs.

Characterization of the diversity, genomic features, host bacteria, and distribution of crAss-like phages in the pig gut microbiome

Phages play an important role in shaping the gut microbiome. CrAss-like phages, which are key members of the gut virome, show high abundance in the human gut and have attracted increasing interest. However, few studies have been found in pigs, and the distribution of crAss-like phages across broader pig populations remains unknown. Here, we obtained 1,251 pig crAss-like phage genomes from 403 metagenomes publicly available and a pig gut virome dataset constructed by ourselves. These crAss-like phage genomes were further clustered into 533 virus operational taxonomic units (vOTUs). Phylogenetic analysis revealed that crAss-like phages in pig guts were distributed across four well-known family-level clusters (Alpha, Beta, Zeta, and Delta) but were absent in the Gamma and Epsilon clusters. Genomic structure analysis identified 149 pig crAss-like phage vOTUs that utilize alternative genetic codes. Gene blocks encoding replication and assembly proteins varied across crAss-like phage clusters. Approximately 64.73% of crAss-like phage genes lacked functional annotations, highlighting a gap in understanding their functional potential. Numerous anti-CRISPR protein genes were identified in crAss-like phage genomes, and CAZymes encoded by these phages were primarily lysozymes. Host prediction indicated that bacterial hosts of pig crAss-like phages primarily belonged to Prevotella, Parabacteroides, and UBA4372. We observed that interactions between crAss-like phages and Prevotella copri might have a possible effect on fat deposition in pigs. Finally, all detected vOTUs exhibited low prevalence across pig populations, suggesting heterogeneity in crAss-like phage compositions. This study provides key resources and novel insights for investigating crAss-like phage-bacteria interactions and benefits research on the effects of crAss-like phages on pig health and production traits.

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

Methylation patterns in bacteria can be used to study restriction-modification or other defense systems with novel properties. While m4C and m6A methylation are well characterized mainly through PacBio sequencing, the landscape of m5C methylation is under-characterized. To bridge this gap, we performed RIMS-seq2 (rapid identification of methyltransferase specificity sequencing) on microbiomes composed of resolved assemblies of distinct genomes through proximity ligation. This high-throughput approach enables the identification of m5C 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 bacteriophage DNA and compared with host DNA, strengthening evidence for phage-host interactions. Applied to three different microbiomes, the method unveiled over 1900 motifs that were deposited in REBASE. The motifs include a novel eight-base recognition site (CATm5CGATG) that was experimentally validated by characterizing its cognate methyltransferase. Our findings suggest that microbiomes harbor arrays of untapped m5C methyltransferase specificities, providing insights into bacterial biology and biotechnological applications.

The chromatin landscape of the histone-possessing Bacteriovorax 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 with 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 (3D) 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 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.

Detection of rare plasmid hosts using a targeted Hi-C approach

Despite the significant role plasmids play in microbial evolution, there is limited knowledge of their ecology, evolution, and transfer in microbial communities. This is partly due to the limitations of current methods in associating a plasmid with its host in microbiomes. To address this knowledge gap, we developed and implemented a novel approach to identify rare plasmid hosts by combining Hi-C, a proximity ligation method, with enrichment for plasmid-specific DNA. We hereafter refer to this approach as Hi-C+. We applied Hi-C and Hi-C+ to soil microbial communities in which we mimicked increasingly rare transfer of an antimicrobial resistance plasmid from a donor to a recipient. This was achieved by inoculating agricultural soil with mixtures of known plasmid-containing and plasmid-free cells at different proportions. We demonstrated that Hi-C can link a plasmid to its host in soil when the relative abundance of that plasmid-host pair is as low as 0.001%. Hi-C+ further improved the detection limit of Hi-C 100-fold and allowed the identification of plasmid hosts at the genus level. As a culture-independent approach, Hi-C+ will significantly improve our understanding of the range and frequency of spread of antibiotic resistance and other genes that are introduced into soil and other microbiomes.

Evolution and function of chromatin domains across the tree of life

The genome of all organisms is spatially organized to function efficiently. The advent of genome-wide chromatin conformation capture (Hi-C) methods has revolutionized our ability to probe the three-dimensional (3D) organization of genomes across diverse species. In this Review, we compare 3D chromatin folding from bacteria and archaea to that in mammals and plants, focusing on topology at the level of gene regulatory domains. In doing so, we consider systematic similarities and differences that hint at the origin and evolution of spatial chromatin folding and its relation to gene activity. We discuss the universality of spatial chromatin domains in all kingdoms, each encompassing one to several genes. We also highlight differences between organisms and suggest that similar features in Hi-C matrices do not necessarily reflect the same biological process or function. Furthermore, we discuss the evolution of domain boundaries and boundary-forming proteins, which indicates that structural maintenance of chromosome (SMC) proteins and the transcription machinery are the ancestral sculptors of the genome. Architectural proteins such as CTCF serve as clade-specific determinants of genome organization. Finally, studies in many non-model organisms show that, despite the ancient origin of 3D chromatin folding and its intricate link to gene activity, evolution tolerates substantial changes in genome organization.

MatP local enrichment delays segregation independently of tetramer formation and septal anchoring in Vibrio cholerae

Vibrio cholerae harbours a primary chromosome derived from the monochromosomal ancestor of the Vibrionales (ChrI) and a secondary chromosome derived from a megaplasmid (ChrII). The coordinated segregation of the replication terminus of both chromosomes (TerI and TerII) determines when and where cell division occurs. ChrI encodes a homologue of Escherichia coli MatP, a protein that binds to a DNA motif (matS) that is overrepresented in replication termini. Here, we use a combination of deep sequencing and fluorescence microscopy techniques to show that V. cholerae MatP structures TerI and TerII into macrodomains, targets them to mid-cell during replication, and delays their segregation, thus supporting that ChrII behaves as a bona fide chromosome. We further show that the extent of the segregation delay mediated by MatP depends on the number and local density of matS sites, and is independent of its assembly into tetramers and any interaction with the divisome, in contrast to what has been previously observed in E. coli.

Imputing Metagenomic Hi-C Contacts Facilitates the Integrative Contig Binning Through Constrained Random Walk with Restart

Metagenomic Hi-C (metaHi-C) has shown remarkable potential for retrieving high-quality metagenome-assembled genomes from complex microbial communities. Nevertheless, existing metaHi-C-based contig binning methods solely rely on Hi-C interactions between contigs, disregarding crucial biological information such as the presence of single-copy marker genes. To overcome this limitation, we introduce ImputeCC, an integrative contig binning tool optimized for metaHi-C datasets. ImputeCC integrates both Hi-C interactions and the discriminative power of single-copy marker genes to group marker-gene-containing contigs into preliminary bins. It also introduces a novel constrained random walk with restart algorithm to enhance Hi-C connectivity among contigs. Comprehensive assessments using both mock and real metaHi-C datasets from diverse environments demonstrate that ImputeCC consistently outperforms other Hi-C-based contig binning tools. A genus-level analysis of the sheep gut microbiota reconstructed by ImputeCC underlines its capability to recover key species from dominant genera and identify previously unknown genera.

Efficient and easy-to-use capturing three-dimensional metagenome interactions with GutHi-C

Hi-C can obtain three-dimensional chromatin structure information and is widely used for genome assembly. We constructed the GutHi-C technology. As shown in the graphical abstract, it is a highly efficient and quick-to-operate method and can be widely used for human, livestock, and poultry gut microorganisms. It provides a reference for the Hi-C methodology of the microbial metagenome. DPBS, Dulbecco's phosphate-buffered saline; Hi-C, high-through chromatin conformation capture; LB, Luria-Bertani; NGS, next-generation sequencing; PCR, polymerase chain reaction; QC, quality control.

A gap-free genome assembly of Fusarium oxysporum f. sp. conglutinans, a vascular wilt pathogen

Fusarium oxysporum is an asexual filamentous fungus that causes vascular wilt in hundreds of crop plants and poses a threat to public health through Fusariosis. F. oxysporum f. sp. conglutinans strain Fo5176, originally isolated from Brassica oleracea, is pathogenic to Arabidopsis, making it a model pathosystem for dissecting the molecular mechanisms underlying host-pathogen interactions. Assembling the F. oxysporum genome is notoriously challenging due to the presence of repeat-rich accessory chromosomes. Here, we report a gap-free genome assembly of Fo5176 using PacBio HiFi and Hi-C data. The 69.56 Mb assembly contained 18 complete chromosomes, including all centromeres and most telomeres (20/36), representing the first gap-free genome sequence of a pathogenic F. oxysporum strain. In total, 21,460 protein-coding genes were annotated, a 26.3% increase compared to the most recent assembly. This high-quality reference genome for F. oxysporum f. sp. conglutinans Fo5176 provides a valuable resource for further research into fungal pathobiology and evolution.

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

Methylation patterns in bacteria can be used to study Restriction-Modification (RM) or other defense systems with novel properties. While m4C and m6A methylation is well characterized mainly through PacBio sequencing, the landscape of m5C 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 m5C 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 (CATm5CGATG) that was experimentally validated by characterizing its cognate methyltransferase. Our findings suggest that microbiomes harbor arrays of untapped m5C 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 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.

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.