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.

TraDIS-Xpress: a high-resolution whole-genome assay identifies novel mechanisms of triclosan action and resistance.

Understanding the genetic basis for a phenotype is a central goal in biological research. Much has been learnt about bacterial genomes by creating large mutant libraries and looking for conditionally important genes. However, current genome-wide methods are largely unable to assay essential genes which are not amenable to disruption. To overcome this limitation, we developed a new version of "TraDIS" (transposon directed insertion-site sequencing) that we term "TraDIS-Xpress" that combines an inducible promoter into the transposon cassette. This allows controlled overexpression and repression of all genes owing to saturation of inserts adjacent to all open reading frames as well as conventional inactivation. We applied TraDIS-Xpress to identify responses to the biocide triclosan across a range of concentrations. Triclosan is endemic in modern life, but there is uncertainty about its mode of action with a concentration-dependent switch from bacteriostatic to bactericidal action unexplained. Our results show a concentration-dependent response to triclosan with different genes important in survival between static and cidal exposures. These genes include those previously reported to have a role in triclosan resistance as well as a new set of genes, including essential genes. Novel genes identified as being sensitive to triclosan exposure include those involved in barrier function, small molecule uptake, and integrity of transcription and translation. We anticipate the approach we show here, by allowing comparisons across multiple experimental conditions of TraDIS data, and including essential genes, will be a starting point for future work examining how different drug conditions impact bacterial survival mechanisms.

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.

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

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.

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.

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.

A genome-wide analysis of Escherichia coli responses to fosfomycin using TraDIS-Xpress reveals novel roles for phosphonate degradation and phosphate transport systems.

BACKGROUND: Fosfomycin is an antibiotic that has seen a revival in use due to its unique mechanism of action and efficacy against isolates resistant to many other antibiotics. In Escherichia coli, fosfomycin often selects for loss-of-function mutations within the genes encoding the sugar importers, GlpT and UhpT. There has, however, not been a genome-wide analysis of the basis for fosfomycin susceptibility reported to date. METHODS: Here we used TraDIS-Xpress, a high-density transposon mutagenesis approach, to assay the role of all genes in E. coli involved in fosfomycin susceptibility. RESULTS: The data confirmed known fosfomycin susceptibility mechanisms and identified new ones. The assay was able to identify domains within proteins of importance and revealed essential genes with roles in fosfomycin susceptibility based on expression changes. Novel mechanisms of fosfomycin susceptibility that were identified included those involved in glucose metabolism and phosphonate catabolism (phnC-M), and the phosphate importer, PstSACB. The impact of these genes on fosfomycin susceptibility was validated by measuring the susceptibility of defined inactivation mutants. CONCLUSIONS: This work reveals a wider set of genes that contribute to fosfomycin susceptibility, including core sugar metabolism genes and two systems involved in phosphate uptake and metabolism previously unrecognized as having a role in fosfomycin susceptibility.

Biocontrol traits of Bacillus licheniformis GL174, a culturable endophyte of Vitis vinifera cv. Glera.

BACKGROUND: Bacillus licheniformis GL174 is a culturable endophytic strain isolated from Vitis vinifera cultivar Glera, the grapevine mainly cultivated for the Prosecco wine production. This strain was previously demonstrated to possess some specific plant growth promoting traits but its endophytic attitude and its role in biocontrol was only partially explored. In this study, the potential biocontrol action of the strain was investigated in vitro and in vivo and, by genome sequence analyses, putative functions involved in biocontrol and plant-bacteria interaction were assessed. RESULTS: Firstly, to confirm the endophytic behavior of the strain, its ability to colonize grapevine tissues was demonstrated and its biocontrol properties were analyzed. Antagonism test results showed that the strain could reduce and inhibit the mycelium growth of diverse plant pathogens in vitro and in vivo. The strain was demonstrated to produce different molecules of the lipopeptide class; moreover, its genome was sequenced, and analysis of the sequences revealed the presence of many protein-coding genes involved in the biocontrol process, such as transporters, plant-cell lytic enzymes, siderophores and other secondary metabolites. CONCLUSIONS: This step-by-step analysis shows that Bacillus licheniformis GL174 may be a good biocontrol agent candidate, and describes some distinguished traits and possible key elements involved in this process. The use of this strain could potentially help grapevine plants to cope with pathogen attacks and reduce the amount of chemicals used in the vineyard.

Comprehensive transcript profiling of two grapevine rootstock genotypes contrasting in drought susceptibility links the phenylpropanoid pathway to enhanced tolerance.

In light of ongoing climate changes in wine-growing regions, the selection of drought-tolerant rootstocks is becoming a crucial factor for developing a sustainable viticulture. In this study, M4, a new rootstock genotype that shows tolerance to drought, was compared from a genomic and transcriptomic point of view with the less drought-tolerant genotype 101.14. The root and leaf transcriptome of both 101.14 and the M4 rootstock genotype was analysed, following exposure to progressive drought conditions. Multifactorial analyses indicated that stress treatment represents the main factor driving differential gene expression in roots, whereas in leaves the genotype is the prominent factor. Upon stress, M4 roots and leaves showed a higher induction of resveratrol and flavonoid biosynthetic genes, respectively. The higher expression of VvSTS genes in M4, confirmed by the accumulation of higher levels of resveratrol in M4 roots compared with 101.14, was coupled to an up-regulation of several VvWRKY transcription factors. Interestingly, VvSTS promoter analyses performed on both the resequenced genomes highlighted a significantly higher number of W-BOX elements in the tolerant genotype. It is proposed that the elevated synthesis of resveratrol in M4 roots upon water stress could enhance the plant's ability to cope with the oxidative stress usually associated with water deficit.

A deep survey of alternative splicing in grape reveals changes in the splicing machinery related to tissue, stress condition and genotype.

BACKGROUND: Alternative splicing (AS) significantly enhances transcriptome complexity. It is differentially regulated in a wide variety of cell types and plays a role in several cellular processes. Here we describe a detailed survey of alternative splicing in grape based on 124 SOLiD RNAseq analyses from different tissues, stress conditions and genotypes. RESULTS: We used the RNAseq data to update the existing grape gene prediction with 2,258 new coding genes and 3,336 putative long non-coding RNAs. Several gene structures have been improved and alternative splicing was described for about 30% of the genes. A link between AS and miRNAs was shown in 139 genes where we found that AS affects the miRNA target site. A quantitative analysis of the isoforms indicated that most of the spliced genes have one major isoform and tend to simultaneously co-express a low number of isoforms, typically two, with intron retention being the most frequent alternative splicing event. CONCLUSIONS: As described in Arabidopsis, also grape displays a marked AS tissue-specificity, while stress conditions produce splicing changes to a minor extent. Surprisingly, some distinctive splicing features were also observed between genotypes. This was further supported by the observation that the panel of Serine/Arginine-rich splicing factors show a few, but very marked differences between genotypes. The finding that a part the splicing machinery can change in closely related organisms can lead to some interesting hypotheses for evolutionary adaptation, that could be particularly relevant in the response to sudden and strong selective pressures.

PASS-bis: a bisulfite aligner suitable for whole methylome analysis of Illumina and SOLiD reads.

SUMMARY: The sequencing of bisulfite-treated DNA (Bi-Seq) is becoming a gold standard for methylation studies. The mapping of Bi-Seq reads is complex and requires special alignment algorithms. This problem is particularly relevant for SOLiD color space, where the bisulfite conversion C/T changes two adjacent colors into 16 possible combinations. Here, we present an algorithm that efficiently aligns Bi-Seq reads obtained either from SOLiD or Illumina. An accompanying methylation-caller program creates a genomic view of methylated and unmethylated Cs on both DNA strands. AVAILABILITY AND IMPLEMENTATION: The algorithm has been implemented as an option of the program PASS, freely available at http://pass.cribi.unipd.it. CONTACT: pass@cribi.unipd.it SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Driving through stop signs: predicting stop codon reassignment improves functional annotation of bacteriophages.

The majority of bacteriophage diversity remains uncharacterized, 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 its 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 median) 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 median coding capacity from 66.8% to 90.0% and from 69.0% to 89.8% for UHGV and INPHARED sequences predicted to use translation table 15. 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.

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.

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.

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.

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.

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.

The transcriptional landscape of the deep-sea bacterium Photobacterium profundum in both a toxR mutant and its parental strain.

BACKGROUND: The deep-sea bacterium Photobacterium profundum is an established model for studying high pressure adaptation. In this paper we analyse the parental strain DB110 and the toxR mutant TW30 by massively parallel cDNA sequencing (RNA-seq). ToxR is a transmembrane DNA-binding protein first discovered in Vibrio cholerae, where it regulates a considerable number of genes involved in environmental adaptation and virulence. In P. profundum the abundance and activity of this protein is influenced by hydrostatic pressure and its role is related to the regulation of genes in a pressure-dependent manner. RESULTS: To better characterize the ToxR regulon, we compared the expression profiles of wt and toxR strains in response to pressure changes. Our results revealed a complex expression pattern with a group of 22 genes having expression profiles similar to OmpH that is an outer membrane protein transcribed in response to high hydrostatic pressure. Moreover, RNA-seq allowed a deep characterization of the transcriptional landscape that led to the identification of 460 putative small RNA genes and the detection of 298 protein-coding genes previously unknown. We were also able to perform a genome-wide prediction of operon structure, transcription start and termination sites, revealing an unexpected high number of genes (992) with large 5'-UTRs, long enough to harbour cis-regulatory RNA structures, suggesting a correlation between intergenic region size and UTR length. CONCLUSION: This work led to a better understanding of high-pressure response in P. profundum. Furthermore, the high-resolution RNA-seq analysis revealed several unexpected features about transcriptional landscape and general mechanisms of controlling bacterial gene expression.

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.