Decoding huge phage diversity: a taxonomic classification of Lak megaphages.

High-throughput sequencing for uncultivated viruses has accelerated the understanding of global viral diversity and uncovered viral genomes substantially larger than any that have so far been cultured. Notably, the Lak phages are an enigmatic group of viruses that present some of the largest known phage genomes identified in human and animal microbiomes, and are dissimilar to any cultivated viruses. Despite the wealth of viral diversity that exists within sequencing datasets, uncultivated viruses have rarely been used for taxonomic classification. We investigated the evolutionary relationships of 23 Lak phages and propose a taxonomy for their classification. Predicted protein analysis revealed the Lak phages formed a deeply branching monophyletic clade within the class Caudoviricetes which contained no other phage genomes. One of the interesting features of this clade is that all current members are characterised by an alternative genetic code. We propose the Lak phages belong to a new order, the 'Grandevirales'. Protein and nucleotide-based analyses support the creation of two families, three sub-families, and four genera within the order 'Grandevirales'. We anticipate that the proposed taxonomy of Lak megaphages will simplify the future classification of related viral genomes as they are uncovered. Continued efforts to classify divergent viruses are crucial to aid common analyses of viral genomes and metagenomes.

Nanopore and Illumina sequencing reveal different viral populations from human gut samples.

The advent of viral metagenomics, or viromics, has improved our knowledge and understanding of global viral diversity. High-throughput sequencing technologies enable explorations of the ecological roles, contributions to host metabolism, and the influence of viruses in various environments, including the human intestinal microbiome. However, bacterial metagenomic studies frequently have the advantage. The adoption of advanced technologies like long-read sequencing has the potential to be transformative in refining viromics and metagenomics. Here, we examined the effectiveness of long-read and hybrid sequencing by comparing Illumina short-read and Oxford Nanopore Technology (ONT) long-read sequencing technologies and different assembly strategies on recovering viral genomes from human faecal samples. Our findings showed that if a single sequencing technology is to be chosen for virome analysis, Illumina is preferable due to its superior ability to recover fully resolved viral genomes and minimise erroneous genomes. While ONT assemblies were effective in recovering viral diversity, the challenges related to input requirements and the necessity for amplification made it less ideal as a standalone solution. However, using a combined, hybrid approach enabled a more authentic representation of viral diversity to be obtained within samples.

Investigating the Human Intestinal DNA Virome and Predicting Disease-Associated Virus-Host Interactions in Severe Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS).

Understanding how the human virome, and which of its constituents, contributes to health or disease states is reliant on obtaining comprehensive virome profiles. By combining DNA viromes from isolated virus-like particles (VLPs) and whole metagenomes from the same faecal sample of a small cohort of healthy individuals and patients with severe myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS), we have obtained a more inclusive profile of the human intestinal DNA virome. Key features are the identification of a core virome comprising tailed phages of the class Caudoviricetes, and a greater diversity of DNA viruses including extracellular phages and integrated prophages. Using an in silico approach, we predicted interactions between members of the Anaerotruncus genus and unique viruses present in ME/CFS microbiomes. This study therefore provides a framework and rationale for studies of larger cohorts of patients to further investigate disease-associated interactions between the intestinal virome and the bacteriome.

Investigating Antibody Reactivity to the Intestinal Microbiome in Severe Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS): A Feasibility Study.

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a multisystemic disease of unknown aetiology that is characterised by disabling chronic fatigue and involves both the immune and gastrointestinal (GI) systems. Patients display alterations in GI microbiome with a significant proportion experiencing GI discomfort and pain and elevated blood biomarkers for altered intestinal permeability compared with healthy individuals. To investigate a possible GI origin of ME/CFS we designed a feasibility study to test the hypothesis that ME/CFS pathogenesis is a consequence of increased intestinal permeability that results in microbial translocation and a breakdown in immune tolerance leading to generation of antibodies reactive to indigenous intestinal microbes. Secretory immunoglobulin (Ig) A and serum IgG levels and reactivity to intestinal microbes were assessed in five pairs of severe ME/CFS patients and matched same-household healthy controls. For profiling serum IgG, we developed IgG-Seq which combines flow-cytometry based bacterial cell sorting and metagenomics to detect mucosal IgG reactivity to the microbiome. We uncovered evidence for immune dysfunction in severe ME/CFS patients that was characterised by reduced capacity and reactivity of serum IgG to stool microbes, irrespective of their source. This study provides the rationale for additional studies in larger cohorts of ME/CFS patients to further explore immune-microbiome interactions.

Draft Genome Sequence of a Primate Isolate of Kazachstania pintolopesii.

Kazachstania pintolopesii is an opportunistic mammalian pathobiont from the K. telluris species complex. No draft genomes of this species are currently available. Here, we report the first draft genome sequence of a primate isolate of K. pintolopesii (NCYC 4417).

Role of mucin glycosylation in the gut microbiota-brain axis of core 3 O-glycan deficient mice.

Alterations in intestinal mucin glycosylation have been associated with increased intestinal permeability and sensitivity to inflammation and infection. Here, we used mice lacking core 3-derived O-glycans (C3GnT-/-) to investigate the effect of impaired mucin glycosylation in the gut-brain axis. C3GnT-/- mice showed altered microbial metabolites in the caecum associated with brain function such as dimethylglycine and N-acetyl-L-tyrosine profiles as compared to C3GnT+/+ littermates. In the brain, polysialylated-neural cell adhesion molecule (PSA-NCAM)-positive granule cells showed an aberrant phenotype in the dentate gyrus of C3GnT-/- mice. This was accompanied by a trend towards decreased expression levels of PSA as well as ZO-1 and occludin as compared to C3GnT+/+. Behavioural studies showed a decrease in the recognition memory of C3GnT-/- mice as compared to C3GnT+/+ mice. Combined, these results support the role of mucin O-glycosylation in the gut in potentially influencing brain function which may be facilitated by the passage of microbial metabolites through an impaired gut barrier.

A hybrid and poly-polish workflow for the complete and accurate assembly of phage genomes: a case study of ten przondoviruses.

Bacteriophages (phages) within the genus Przondovirus are T7-like podoviruses belonging to the subfamily Studiervirinae, within the family Autographiviridae, and have a highly conserved genome organisation. The genomes of these phages range from 37 to 42 kb in size, encode 50-60 genes and are characterised by the presence of direct terminal repeats (DTRs) flanking the linear chromosome. These DTRs are often deleted during short-read-only and hybrid assemblies. Moreover, long-read-only assemblies are often littered with sequencing and/or assembly errors and require additional curation. Here, we present the isolation and characterisation of ten novel przondoviruses targeting Klebsiella spp. We describe HYPPA, a HYbrid and Poly-polish Phage Assembly workflow, which utilises long-read assemblies in combination with short-read sequencing to resolve phage DTRs and correcting errors, negating the need for laborious primer walking and Sanger sequencing validation. Our assembly workflow utilised Oxford Nanopore Technologies for long-read sequencing for its accessibility, making it the more relevant long-read sequencing technology at this time, and Illumina DNA Prep for short-read sequencing, representing the most commonly used technologies globally. Our data demonstrate the importance of careful curation of phage assemblies before publication, and prior to using them for comparative genomics.

Draft Genome Sequence of a Preterm Infant-Derived Isolate of Candida parapsilosis.

Candida parapsilosis is a human fungal pathogen of increasing incidence and causes invasive candidiasis, notably in preterm or low-birthweight neonates. Here, we present the genome sequence of C. parapsilosis NCYC 4289, a fecal isolate from a preterm male infant.

Predicting stop codon reassignment improves functional annotation of bacteriophages.

The majority of bacteriophage diversity remains uncharacterised, and new intriguing mechanisms of their biology are being continually described. Members of some phage lineages, such as the Crassvirales, repurpose stop codons to encode an amino acid by using alternate genetic codes. Here, we investigated the prevalence of stop codon reassignment in phage genomes and subsequent impacts on functional annotation. We predicted 76 genomes within INPHARED and 712 vOTUs from the Unified Human Gut Virome catalogue (UHGV) that repurpose a stop codon to encode an amino acid. We re-annotated these sequences with modified versions of Pharokka and Prokka, called Pharokka-gv and Prokka-gv, to automatically predict stop codon reassignment prior to annotation. Both tools significantly improved the quality of annotations, with Pharokka-gv performing best. For sequences predicted to repurpose TAG to glutamine (translation table 15), Pharokka-gv increased the median gene length (median of per genome medians) from 287 to 481 bp for UHGV sequences (67.8% increase) and from 318 to 550 bp for INPHARED sequences (72.9% increase). The re-annotation increased mean coding density from 66.8% to 90.0%, and from 69.0% to 89.8% for UHGV and INPHARED sequences. Furthermore, the proportion of genes that could be assigned functional annotation increased, including an increase in the number of major capsid proteins that could be identified. We propose that automatic prediction of stop codon reassignment before annotation is beneficial to downstream viral genomic and metagenomic analyses.

Genome-Wide Analysis of Innate Susceptibility Mechanisms of Escherichia coli to Colistin.

Colistin is an antibiotic that has seen increasing clinical use for the treatment of human infections caused by Gram-negative pathogens, particularly due to the emergence of multidrug-resistant pathogens. Colistin resistance is also a growing problem and typically results from alterations to lipopolysaccharides mediated by phosphoethanolamine (pETn) transferase enzymes which can be encoded on the chromosome, or plasmids. In this study, we used 'TraDIS-Xpress' (Transposon Directed Insertion site Sequencing with expression), where a high-density transposon mutant library including outward facing promoters in Escherichia coli BW25113 identified genes involved in colistin susceptibility. We examined the genome-wide response of E. coli following exposure to a range of concentrations of colistin. Our TraDIS-Xpress screen confirmed the importance of overexpression of the two-component system basSR (which regulates pETn transferases) but also identified a wider range of genes important for survival in the presence of colistin, including genes encoding membrane associated proteins, DNA repair machinery, various transporters, RNA helicases, general stress response genes, fimbriae and phosphonate metabolism. Validation experiments supported a role in colistin susceptibility for novel candidate genes tested. TraDIS-Xpress is a powerful tool that expands our understanding of the wider landscape of genes involved in response to colistin susceptibility mechanisms.

The Cynomolgus Macaque Intestinal Mycobiome Is Dominated by the Kazachstania Genus and K. pintolopesii Species.

The cynomolgus macaque, Macaca fascicularis, is a non-human primate (NHP) widely used in biomedical research as its genetics, immunology and physiology are similar to those of humans. They may also be a useful model of the intestinal microbiome as their prokaryome resembles that of humans. However, beyond the prokaryome relatively little is known about other constituents of the macaque intestinal microbiome including the mycobiome. Here, we conducted a region-by-region taxonomic survey of the cynomolgus intestinal mycobiota, from duodenum to distal colon, of sixteen captive animals of differing age (from young to old). Using a high-throughput ITS1 amplicon sequencing-based approach, the cynomolgus gut mycobiome was dominated by fungi from the Ascomycota phylum. The budding yeast genus Kazachstania was most abundant, with the thermotolerant species K. pintolopesii highly prevalent, and the predominant species in both the small and large intestines. This is in marked contrast to humans, in which the intestinal mycobiota is characterised by other fungal genera including Candida and Saccharomyces, and Candida albicans. This study provides a comprehensive insight into the fungal communities present within the captive cynomolgus gut, and for the first time identifies K. pintolopesii as a candidate primate gut commensal.

MetaPhage: an Automated Pipeline for Analyzing, Annotating, and Classifying Bacteriophages in Metagenomics Sequencing Data.

Phages are the most abundant biological entities on the planet, and they play an important role in controlling density, diversity, and network interactions among bacterial communities through predation and gene transfer. To date, a variety of bacteriophage identification tools have been developed that differ in the phage mining strategies used, input files requested, and results produced. However, new users attempting bacteriophage analysis can struggle to select the best methods and interpret the variety of results produced. Here, we present MetaPhage, a comprehensive reads-to-report pipeline that streamlines the use of multiple phage miners and generates an exhaustive report. The report both summarizes and visualizes the key findings and enables further exploration of key results via interactive filterable tables. The pipeline is implemented in Nextflow, a widely adopted workflow manager that enables an optimized parallelization of tasks in different locations, from local server to the cloud; this ensures reproducible results from containerized packages. MetaPhage is designed to enable scalability and reproducibility; also, it can be easily expanded to include new miners and methods as they are developed in this continuously growing field. MetaPhage is freely available under a GPL-3.0 license at https://github.com/MattiaPandolfoVR/MetaPhage. IMPORTANCE Bacteriophages (viruses that infect bacteria) are the most abundant biological entities on earth and are increasingly studied as members of the resident microbiota community in many environments, from oceans to soils and the human gut. Their identification is of great importance to better understand complex bacterial dynamics and microbial ecosystem function. A variety of metagenome bacteriophage identification tools have been developed that differ in the phage mining strategies used, input files requested, and results produced. To facilitate the management and the execution of such a complex workflow, we developed MetaPhage (MP), a comprehensive reads-to-report pipeline that streamlines the use of multiple phage miners and generates an exhaustive report. The pipeline is implemented in Nextflow, a widely adopted workflow manager that enables an optimized parallelization of tasks. MetaPhage is designed to enable scalability and reproducibility and offers an installation-free, dependency-free, and conflict-free workflow execution.

Absence of Bacteria Permits Fungal Gut-To-Brain Translocation and Invasion in Germfree Mice but Ageing Alone Does Not Drive Pathobiont Expansion in Conventionally Raised Mice.

Age-associated changes in the structure of the intestinal microbiome and in its interaction with the brain via the gut-brain axis are increasingly being implicated in neurological and neurodegenerative diseases. Intestinal microbial dysbiosis and translocation of microbes and microbial products including fungal species into the brain have been implicated in the development of dementias such as Alzheimer's disease. Using germ-free mice, we investigated if the fungal gut commensal, Candida albicans, an opportunistic pathogen in humans, can traverse the gastrointestinal barrier and disseminate to brain tissue and whether ageing impacts on the gut mycobiome as a pre-disposing factor in fungal brain infection. C. albicans was detected in different regions of the brain of colonised germ-free mice in both yeast and hyphal cell forms, often in close association with activated (Iba-1+) microglial cells. Using high-throughput ITS1 amplicon sequencing to characterise the faecal gut fungal composition of aged and young SPF mice, we identified several putative gut commensal fungal species with pathobiont potential although their abundance was not significantly different between young and aged mice. Collectively, these results suggest that although some fungal species can travel from the gut to brain where they can induce an inflammatory response, ageing alone is not correlated with significant changes in gut mycobiota composition which could predispose to these events. These results are consistent with a scenario in which significant disruptions to the gut microbiota or intestinal barrier, beyond those which occur with natural ageing, are required to allow fungal escape and brain infection.

Fecal microbiota transfer between young and aged mice reverses hallmarks of the aging gut, eye, and brain.

BACKGROUND: Altered intestinal microbiota composition in later life is associated with inflammaging, declining tissue function, and increased susceptibility to age-associated chronic diseases, including neurodegenerative dementias. Here, we tested the hypothesis that manipulating the intestinal microbiota influences the development of major comorbidities associated with aging and, in particular, inflammation affecting the brain and retina. METHODS: Using fecal microbiota transplantation, we exchanged the intestinal microbiota of young (3 months), old (18 months), and aged (24 months) mice. Whole metagenomic shotgun sequencing and metabolomics were used to develop a custom analysis workflow, to analyze the changes in gut microbiota composition and metabolic potential. Effects of age and microbiota transfer on the gut barrier, retina, and brain were assessed using protein assays, immunohistology, and behavioral testing. RESULTS: We show that microbiota composition profiles and key species enriched in young or aged mice are successfully transferred by FMT between young and aged mice and that FMT modulates resulting metabolic pathway profiles. The transfer of aged donor microbiota into young mice accelerates age-associated central nervous system (CNS) inflammation, retinal inflammation, and cytokine signaling and promotes loss of key functional protein in the eye, effects which are coincident with increased intestinal barrier permeability. Conversely, these detrimental effects can be reversed by the transfer of young donor microbiota. CONCLUSIONS: These findings demonstrate that the aging gut microbiota drives detrimental changes in the gut-brain and gut-retina axes suggesting that microbial modulation may be of therapeutic benefit in preventing inflammation-related tissue decline in later life. Video abstract.

BamToCov: an efficient toolkit for sequence coverage calculations.

MOTIVATION: Many genomics applications require the computation of nucleotide coverage of a reference genome or the ability to determine how many reads map to a reference region. RESULTS: BamToCov is a toolkit for rapid and flexible coverage computation that relies on the most memory efficient algorithm and is designed for integration in pipelines, given its ability to read alignment files from streams. The tools in the suite can process sorted BAM or CRAM files, allowing the user to extract coverage information via different filtering approaches and to save the output in different formats (BED, Wig or counts). The BamToCov algorithm can also handle strand-specific and/or physical coverage analyses. AVAILABILITY AND IMPLEMENTATION: This program, accessory utilities and their documentation are freely available at https://github.com/telatin/BamToCov. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Chemical biology-whole genome engineering datasets predict new antibacterial combinations.

Trimethoprim and sulfamethoxazole are used commonly together as cotrimoxazole for the treatment of urinary tract and other infections. The evolution of resistance to these and other antibacterials threatens therapeutic options for clinicians. We generated and analysed a chemical-biology-whole-genome data set to predict new targets for antibacterial combinations with trimethoprim and sulfamethoxazole. For this we used a large transposon mutant library in Escherichia coli BW25113 where an outward-transcribing inducible promoter was engineered into one end of the transposon. This approach allows regulated expression of adjacent genes in addition to gene inactivation at transposon insertion sites, a methodology that has been called TraDIS-Xpress. These chemical genomic data sets identified mechanisms for both reduced and increased susceptibility to trimethoprim and sulfamethoxazole. The data identified that over-expression of FolA reduced trimethoprim susceptibility, a known mechanism for reduced susceptibility. In addition, transposon insertions into the genes tdk, deoR, ybbC, hha, ldcA, wbbK and waaS increased susceptibility to trimethoprim and likewise for rsmH, fadR, ddlB, nlpI and prc with sulfamethoxazole, while insertions in ispD, uspC, minC, minD, yebK, truD and umpG increased susceptibility to both these antibiotics. Two of these genes' products, Tdk and IspD, are inhibited by AZT and fosmidomycin respectively, antibiotics that are known to synergise with trimethoprim. Thus, the data identified two known targets and several new target candidates for the development of co-drugs that synergise with trimethoprim, sulfamethoxazole or cotrimoxazole. We demonstrate that the TraDIS-Xpress technology can be used to generate information-rich chemical-genomic data sets that can be used for antibacterial development.

Evaluating the potential for respiratory metagenomics to improve treatment of secondary infection and detection of nosocomial transmission on expanded COVID-19 intensive care units.

BACKGROUND: Clinical metagenomics (CMg) has the potential to be translated from a research tool into routine service to improve antimicrobial treatment and infection control decisions. The SARS-CoV-2 pandemic provides added impetus to realise these benefits, given the increased risk of secondary infection and nosocomial transmission of multi-drug-resistant (MDR) pathogens linked with the expansion of critical care capacity. METHODS: CMg using nanopore sequencing was evaluated in a proof-of-concept study on 43 respiratory samples from 34 intubated patients across seven intensive care units (ICUs) over a 9-week period during the first COVID-19 pandemic wave. RESULTS: An 8-h CMg workflow was 92% sensitive (95% CI, 75-99%) and 82% specific (95% CI, 57-96%) for bacterial identification based on culture-positive and culture-negative samples, respectively. CMg sequencing reported the presence or absence of ß-lactam-resistant genes carried by Enterobacterales that would modify the initial guideline-recommended antibiotics in every case. CMg was also 100% concordant with quantitative PCR for detecting Aspergillus fumigatus from 4 positive and 39 negative samples. Molecular typing using 24-h sequencing data identified an MDR-K. pneumoniae ST307 outbreak involving 4 patients and an MDR-C. striatum outbreak involving 14 patients across three ICUs. CONCLUSION: CMg testing provides accurate pathogen detection and antibiotic resistance prediction in a same-day laboratory workflow, with assembled genomes available the next day for genomic surveillance. The provision of this technology in a service setting could fundamentally change the multi-disciplinary team approach to managing ICU infections. The potential to improve the initial targeted treatment and rapidly detect unsuspected outbreaks of MDR-pathogens justifies further expedited clinical assessment of CMg.

The Origin of Plasma-Derived Bacterial Extracellular Vesicles in Healthy Individuals and Patients with Inflammatory Bowel Disease: A Pilot Study.

The gastrointestinal tract harbors the gut microbiota, structural alterations of which (dysbiosis) are linked with an increase in gut permeability ("leaky gut"), enabling luminal antigens and bacterial products such as nanosized bacterial extracellular vesicles (BEVs) to access the circulatory system. Blood-derived BEVs contain various cargoes and may be useful biomarkers for diagnosis and monitoring of disease status and relapse in conditions such as inflammatory bowel disease (IBD). To progress this concept, we developed a rapid, cost-effective protocol to isolate BEV-associated DNA and used 16S rRNA gene sequencing to identify bacterial origins of the blood microbiome of healthy individuals and patients with Crohn's disease and ulcerative colitis. The 16S rRNA gene sequencing successfully identified the origin of plasma-derived BEV DNA. The analysis showed that the blood microbiota richness, diversity, or composition in IBD, healthy control, and protocol control groups were not significantly distinct, highlighting the issue of 'kit-ome' contamination in low-biomass studies. Our pilot study provides the basis for undertaking larger studies to determine the potential use of blood microbiota profiling as a diagnostic aid in IBD.

Comparison of PCR versus PCR-Free DNA Library Preparation for Characterising the Human Faecal Virome.

The human intestinal microbiota is abundant in viruses, comprising mainly bacteriophages, occasionally outnumbering bacteria 10:1 and is termed the virome. Due to their high genetic diversity and the lack of suitable tools and reference databases, the virome remains poorly characterised and is often referred to as "viral dark matter". However, the choice of sequencing platforms, read lengths and library preparation make study design challenging with respect to the virome. Here we have compared the use of PCR and PCR-free methods for sequence-library construction on the Illumina sequencing platform for characterising the human faecal virome. Viral DNA was extracted from faecal samples of three healthy donors and sequenced. Our analysis shows that most variation was reflecting the individually specific faecal virome. However, we observed differences between PCR and PCR-free library preparation that affected the recovery of low-abundance viral genomes. Using three faecal samples in this study, the PCR library preparation samples led to a loss of lower-abundance vOTUs evident in their PCR-free pairs (vOTUs 128, 6202 and 8364) and decreased the alpha-diversity indices (Chao1 p-value = 0.045 and Simpson p-value = 0.044). Thus, differences between PCR and PCR-free methods are important to consider when investigating "rare" members of the gut virome, with these biases likely negligible when investigating moderately and highly abundant viruses.

Dadaist2: A Toolkit to Automate and Simplify Statistical Analysis and Plotting of Metabarcoding Experiments.

The taxonomic composition of microbial communities can be assessed using universal marker amplicon sequencing. The most common taxonomic markers are the 16S rDNA for bacterial communities and the internal transcribed spacer (ITS) region for fungal communities, but various other markers are used for barcoding eukaryotes. A crucial step in the bioinformatic analysis of amplicon sequences is the identification of representative sequences. This can be achieved using a clustering approach or by denoising raw sequencing reads. DADA2 is a widely adopted algorithm, released as an R library, that denoises marker-specific amplicons from next-generation sequencing and produces a set of representative sequences referred to as 'Amplicon Sequence Variants' (ASV). Here, we present Dadaist2, a modular pipeline, providing a complete suite for the analysis that ranges from raw sequencing reads to the statistics of numerical ecology. Dadaist2 implements a new approach that is specifically optimised for amplicons with variable lengths, such as the fungal ITS. The pipeline focuses on streamlining the data flow from the command line to R, with multiple options for statistical analysis and plotting, both interactive and automatic.