PRFect: a tool to predict programmed ribosomal frameshifts in prokaryotic and viral genomes

BACKGROUND

One of the stranger phenomena that can occur during gene translation is where, as a ribosome reads along the mRNA, various cellular and molecular properties contribute to stalling the ribosome on a slippery sequence and shifting the ribosome into one of the other two alternate reading frames. The alternate frame has different codons, so different amino acids are added to the peptide chain. More importantly, the original stop codon is no longer in-frame, so the ribosome can bypass the stop codon and continue to translate the codons past it. This produces a longer version of the protein, a fusion of the original in-frame amino acids, followed by all the alternate frame amino acids. There is currently no automated software to predict the occurrence of these programmed ribosomal frameshifts (PRF), and they are currently only identified by manual curation.

RESULTS

Here we present PRFect, an innovative machine-learning method for the detection and prediction of PRFs in coding genes of various types. PRFect combines advanced machine learning techniques with the integration of multiple complex cellular properties, such as secondary structure, codon usage, ribosomal binding site interference, direction, and slippery site motif. Calculating and incorporating these diverse properties posed significant challenges, but through extensive research and development, we have achieved a user-friendly approach. The code for PRFect is freely available, open-source, and can be easily installed via a single command in the terminal. Our comprehensive evaluations on diverse organisms, including bacteria, archaea, and phages, demonstrate PRFect's strong performance, achieving high sensitivity, specificity, and an accuracy exceeding 90%. The code for PRFect is freely available and installs with a single terminal command.

CONCLUSION

PRFect represents a significant advancement in the field of PRF detection and prediction, offering a powerful tool for researchers and scientists to unravel the intricacies of programmed ribosomal frameshifting in coding genes.

Women's perceptions and preferences for cervical cancer screening in light of updated guidelines

BACKGROUND

Updated risk-based guidelines for cervical cancer screening (CCS) allow for individualized screening that minimizes unnecessary testing. However, these guidelines are complex and may not be easily understandable to patients.

PURPOSE

To describe women's perceptions and preferences about CCS in light of recent guideline changes.

METHODOLOGY

This qualitative study conducted in February 2020 used online, text-based focus groups with a sample of US women ( n = 49) ages 27-45. After participants completed a short demographic survey, an experienced moderator used a semistructured guide to solicit participants' perceptions of routine gynecologic care and CCS. We summarized survey data using descriptive statistics. Two authors analyzed transcripts using conventional content analysis and met with other team members to resolve discrepancies and determine final themes.

RESULTS

Most participants were non-Hispanic White (65%), had health insurance (90%), and reported having a routine gynecologic examination in the past year (70%). We identified four common themes: (1) low perceived risk of human papillomavirus (HPV) coupled with low knowledge about HPV as a causative factor for cervical cancer, (2) confusion about-and mistrust of-recent individual risk-based guidelines that determine the indicated timing and type of CCS test (Pap or HPV testing), (3) mixed opinions about performing a self-swab for HPV testing, and (4) conflicting perceptions of trust toward providers and the health care industry.

CONCLUSIONS

Findings highlight women's uncertainty and hesitancy about updated CCS guidelines.

IMPLICATIONS

Provider-patient communication strategies should consider women's gaps in knowledge about HPV, include the rationale for guidelines and types of tests, and build trust between patients and providers.

Host interactions of novel Crassvirales species belonging to multiple families infecting bacterial host, Bacteroides cellulosilyticus WH2

Bacteroides, the prominent bacteria in the human gut, play a crucial role in degrading complex polysaccharides. Their abundance is influenced by phages belonging to the Crassvirales order. Despite identifying over 600 Crassvirales genomes computationally, only few have been successfully isolated. Continued efforts in isolation of more Crassvirales genomes can provide insights into phage-host-evolution and infection mechanisms. We focused on wastewater samples, as potential sources of phages infecting various Bacteroides hosts. Sequencing, assembly, and characterization of isolated phages revealed 14 complete genomes belonging to three novel Crassvirales species infecting Bacteroides cellulosilyticus WH2. These species, Kehishuvirus sp. 'tikkala' strain Bc01, Kolpuevirus sp. 'frurule' strain Bc03, and 'Rudgehvirus jaberico' strain Bc11, spanned two families, and three genera, displaying a broad range of virion productions. Upon testing all successfully cultured Crassvirales species and their respective bacterial hosts, we discovered that they do not exhibit co-evolutionary patterns with their bacterial hosts. Furthermore, we observed variations in gene similarity, with greater shared similarity observed within genera. However, despite belonging to different genera, the three novel species shared a unique structural gene that encodes the tail spike protein. When investigating the relationship between this gene and host interaction, we discovered evidence of purifying selection, indicating its functional importance. Moreover, our analysis demonstrated that this tail spike protein binds to the TonB-dependent receptors present on the bacterial host surface. Combining these observations, our findings provide insights into phage-host interactions and present three Crassvirales species as an ideal system for controlled infectivity experiments on one of the most dominant members of the human enteric virome.

Compounding Achromobacter Phages for Therapeutic Applications

Achromobacter species colonization of Cystic Fibrosis respiratory airways is an increasing concern. Two adult patients with Cystic Fibrosis colonized by Achromobacter xylosoxidans CF418 or Achromobacter ruhlandii CF116 experienced fatal exacerbations. Achromobacter spp. are naturally resistant to several antibiotics. Therefore, phages could be valuable as therapeutics for the control of Achromobacter. In this study, thirteen lytic phages were isolated and characterized at the morphological and genomic levels for potential future use in phage therapy. They are presented here as the Achromobacter Kumeyaay phage collection. Six distinct Achromobacter phage genome clusters were identified based on a comprehensive phylogenetic analysis of the Kumeyaay collection as well as the publicly available Achromobacter phages. The infectivity of all phages in the Kumeyaay collection was tested in 23 Achromobacter clinical isolates; 78% of these isolates were lysed by at least one phage. A cryptic prophage was induced in Achromobacter xylosoxidans CF418 when infected with some of the lytic phages. This prophage genome was characterized and is presented as Achromobacter phage CF418-P1. Prophage induction during lytic phage preparation for therapy interventions require further exploration. Large-scale production of phages and removal of endotoxins using an octanol-based procedure resulted in a phage concentrate of 1 × 109 plaque-forming units per milliliter with an endotoxin concentration of 65 endotoxin units per milliliter, which is below the Food and Drugs Administration recommended maximum threshold for human administration. This study provides a comprehensive framework for the isolation, bioinformatic characterization, and safe production of phages to kill Achromobacter spp. in order to potentially manage Cystic Fibrosis (CF) pulmonary infections.

Host interactions of novel Crassvirales species belonging to multiple families infecting bacterial host, Bacteroides cellulosilyticus WH2

Bacteroides, the prominent bacteria in the human gut, play a crucial role in degrading complex polysaccharides. Their abundance is influenced by phages belonging to the Crassvirales order. Despite identifying over 600 Crassvirales genomes computationally, only few have been successfully isolated. Continued efforts in isolation of more Crassvirales genomes can provide insights into phage-host-evolution and infection mechanisms. We focused on wastewater samples, as potential sources of phages infecting various Bacteroides hosts. Sequencing, assembly, and characterization of isolated phages revealed 14 complete genomes belonging to three novel Crassvirales species infecting Bacteroides cellulosilyticus WH2. These species, Kehishuvirus sp. 'tikkala' strain Bc01, Kolpuevirus sp. 'frurule' strain Bc03, and 'Rudgehvirus jaberico' strain Bc11, spanned two families, and three genera, displaying a broad range of virion productions. Upon testing all successfully cultured Crassvirales species and their respective bacterial hosts, we discovered that they do not exhibit co-evolutionary patterns with their bacterial hosts. Furthermore, we observed variations in gene similarity, with greater shared similarity observed within genera. However, despite belonging to different genera, the three novel species shared a unique structural gene that encodes the tail spike protein. When investigating the relationship between this gene and host interaction, we discovered evidence of purifying selection, indicating its functional importance. Moreover, our analysis demonstrated that this tail spike protein binds to the TonB-dependent receptors present on the bacterial host surface. Combining these observations, our findings provide insights into phage-host interactions and present three Crassvirales species as an ideal system for controlled infectivity experiments on one of the most dominant members of the human enteric virome.

Impact statement

Bacteriophages play a crucial role in shaping microbial communities within the human gut. Among the most dominant bacteriophages in the human gut microbiome are Crassvirales phages, which infect Bacteroides. Despite being widely distributed, only a few Crassvirales genomes have been isolated, leading to a limited understanding of their biology, ecology, and evolution. This study isolated and characterized three novel Crassvirales genomes belonging to two different families, and three genera, but infecting one bacterial host, Bacteroides cellulosilyticus WH2. Notably, the observation confirmed the phages are not co-evolving with their bacterial hosts, rather have a shared ability to exploit similar features in their bacterial host. Additionally, the identification of a critical viral protein undergoing purifying selection and interacting with the bacterial receptors opens doors to targeted therapies against bacterial infections. Given Bacteroides role in polysaccharide degradation in the human gut, our findings advance our understanding of the phage-host interactions and could have important implications for the development of phage-based therapies. These discoveries may hold implications for improving gut health and metabolism to support overall well-being.

Data summary

The genomes used in this research are available on Sequence Read Archive (SRA) within the project, PRJNA737576. Bacteroides cellulosilyticus WH2, Kehishuvirus sp. 'tikkala' strain Bc01, Kolpuevirus sp. ' frurule' strain Bc03, and 'Rudgehvirus jaberico' strain Bc11 are all available on GenBank with accessions NZ_CP072251.1 ( B. cellulosilyticus WH2), QQ198717 (Bc01), QQ198718 (Bc03), and QQ198719 (Bc11), and we are working on making the strains available through ATCC. The 3D protein structures for the three Crassvirales genomes are available to download at doi.org/10.25451/flinders.21946034.

PRFect: A tool to predict programmed ribosomal frameshifts in prokaryotic and viral genomes

Background

One of the stranger phenomena that can occur during gene translation is where, as a ribosome reads along the mRNA, various cellular and molecular properties contribute to stalling the ribosome on a slippery sequence, shifting the ribosome into one of the other two alternate reading frames. The alternate frame has different codons, so different amino acids are added to the peptide chain, but more importantly, the original stop codon is no longer in-frame, so the ribosome can bypass the stop codon and continue to translate the codons past it. This produces a longer version of the protein, a fusion of the original in-frame amino acids, followed by all the alternate frame amino acids. There is currently no automated software to predict the occurrence of these programmed ribosomal frameshifts (PRF), and they are currently only identified by manual curation.

Results

Here we present PRFect, an innovative machine-learning method for the detection and prediction of PRFs in coding genes of various types. PRFect combines advanced machine learning techniques with the integration of multiple complex cellular properties, such as secondary structure, codon usage, ribosomal binding site interference, direction, and slippery site motif. Calculating and incorporating these diverse properties posed significant challenges, but through extensive research and development, we have achieved a user-friendly approach. The code for PRFect is freely available, open-source, and can be easily installed via a single command in the terminal. Our comprehensive evaluations on diverse organisms, including bacteria, archaea, and phages, demonstrate PRFect's strong performance, achieving high sensitivity, specificity, and an accuracy exceeding 90%.

Conclusion

PRFect represents a significant advancement in the field of PRF detection and prediction, offering a powerful tool for researchers and scientists to unravel the intricacies of programmed ribosomal frameshifting in coding genes.

The Promise and Pitfalls of Prophages

Phages dominate every ecosystem on the planet. While virulent phages sculpt the microbiome by killing their bacterial hosts, temperate phages provide unique growth advantages to their hosts through lysogenic conversion. Many prophages benefit their host, and prophages are responsible for genotypic and phenotypic differences that separate individual microbial strains. However, the microbes also endure a cost to maintain those phages: additional DNA to replicate and proteins to transcribe and translate. We have never quantified those benefits and costs. Here, we analysed over two and a half million prophages from over half a million bacterial genome assemblies. Analysis of the whole dataset and a representative subset of taxonomically diverse bacterial genomes demonstrated that the normalised prophage density was uniform across all bacterial genomes above 2 Mbp. We identified a constant carrying capacity of phage DNA per bacterial DNA. We estimated that each prophage provides cellular services equivalent to approximately 2.4 % of the cell's energy or 0.9 ATP per bp per hour. We demonstrate analytical, taxonomic, geographic, and temporal disparities in identifying prophages in bacterial genomes that provide novel targets for identifying new phages. We anticipate that the benefits bacteria accrue from the presence of prophages balance the energetics involved in supporting prophages. Furthermore, our data will provide a new framework for identifying phages in environmental datasets, diverse bacterial phyla, and from different locations.

MultiPhATE2: code for functional annotation and comparison of phage genomes

To address a need for improved tools for annotation and comparative genomics of bacteriophage genomes, we developed multiPhATE2. As an extension of multiPhATE, a functional annotation code released previously, multiPhATE2 performs gene finding using multiple algorithms, compares the results of the algorithms, performs functional annotation of coding sequences, and incorporates additional search algorithms and databases to extend the search space of the original code. MultiPhATE2 performs gene matching among sets of closely related bacteriophage genomes, and uses multiprocessing to speed computations. MultiPhATE2 can be re-started at multiple points within the workflow to allow the user to examine intermediate results and adjust the subsequent computations accordingly. In addition, multiPhATE2 accommodates custom gene calls and sequence databases, again adding flexibility. MultiPhATE2 was implemented in Python 3.7 and runs as a command-line code under Linux or MAC operating systems. Full documentation is provided as a README file and a Wiki website.

A Distinct Contractile Injection System Gene Cluster Found in a Majority of Healthy Adult Human Microbiomes

Many commensal bacteria antagonize each other or their host by producing syringe-like secretion systems called contractile injection systems (CIS). Members of the Bacteroidales family have been shown to produce only one type of CIS-a contact-dependent type 6 secretion system that mediates bacterium-bacterium interactions. Here, we show that a second distinct cluster of genes from Bacteroidales bacteria from the human microbiome may encode yet-uncharacterized injection systems that we term Bacteroidales injection systems (BIS). We found that BIS genes are present in the gut microbiomes of 99% of individuals from the United States and Europe and that BIS genes are more prevalent in the gut microbiomes of healthy individuals than in those individuals suffering from inflammatory bowel disease. Gene clusters similar to that of the BIS mediate interactions between bacteria and diverse eukaryotes, like amoeba, insects, and tubeworms. Our findings highlight the ubiquity of the BIS gene cluster in the human gut and emphasize the relevance of the gut microbiome to the human host. These results warrant investigations into the structure and function of the BIS and how they might mediate interactions between Bacteroidales bacteria and the human host or microbiome.IMPORTANCE To engage with host cells, diverse pathogenic bacteria produce syringe-like structures called contractile injection systems (CIS). CIS are evolutionarily related to the contractile tails of bacteriophages and are specialized to puncture membranes, often delivering effectors to target cells. Although CIS are key for pathogens to cause disease, paradoxically, similar injection systems have been identified within healthy human microbiome bacteria. Here, we show that gene clusters encoding a predicted CIS, which we term Bacteroidales injection systems (BIS), are present in the microbiomes of nearly all adult humans tested from Western countries. BIS genes are enriched within human gut microbiomes and are expressed both in vitro and in vivo Further, a greater abundance of BIS genes is present within healthy gut microbiomes than in those humans with with inflammatory bowel disease (IBD). Our discovery provides a potentially distinct means by which our microbiome interacts with the human host or its microbiome.

Dietary prophage inducers and antimicrobials: toward landscaping the human gut microbiome

The approximately 10¹¹ viruses and microbial cells per gram of fecal matter (dry weight) in the large intestine are important to human health. The responses of three common gut bacteria species, and one opportunistic pathogen, to 117 commonly consumed foods, chemical additives, and plant extracts were tested. Many compounds, including Stevia rebaudiana and bee propolis extracts, exhibited species-specific growth inhibition by prophage induction. Overall, these results show that various foods may change the abundances of gut bacteria by modulating temperate phage and suggests a novel path for landscaping the human gut microbiome.

multiPhATE: bioinformatics pipeline for functional annotation of phage isolates

Summary

To address the need for improved phage annotation tools that scale, we created an automated throughput annotation pipeline: multiple-genome Phage Annotation Toolkit and Evaluator (multiPhATE). multiPhATE is a throughput pipeline driver that invokes an annotation pipeline (PhATE) across a user-specified set of phage genomes. This tool incorporates a de novo phage gene calling algorithm and assigns putative functions to gene calls using protein-, virus- and phage-centric databases. multiPhATE’s modular construction allows the user to implement all or any portion of the analyses by acquiring local instances of the desired databases and specifying the desired analyses in a configuration file. We demonstrate multiPhATE by annotating two newly sequenced Yersinia pestis phage genomes. Within multiPhATE, the PhATE processing pipeline can be readily implemented across multiple processors, making it adaptable for throughput sequencing projects. Software documentation assists the user in configuring the system.

Availability and implementation

multiPhATE was implemented in Python 3.7, and runs as a command-line code under Linux or Unix. multiPhATE is freely available under an open-source BSD3 license from https://github.com/carolzhou/multiPhATE. Instructions for acquiring the databases and third-party codes used by multiPhATE are included in the distribution README file. Users may report bugs by submitting to the github issues page associated with the multiPhATE distribution.

Supplementary information

Supplementary data are available at Bioinformatics online.

PHANOTATE: a novel approach to gene identification in phage genomes

MOTIVATION

Currently there are no tools specifically designed for annotating genes in phages. Several tools are available that have been adapted to run on phage genomes, but due to their underlying design, they are unable to capture the full complexity of phage genomes. Phages have adapted their genomes to be extremely compact, having adjacent genes that overlap and genes completely inside of other longer genes. This non-delineated genome structure makes it difficult for gene prediction using the currently available gene annotators. Here we present PHANOTATE, a novel method for gene calling specifically designed for phage genomes. Although the compact nature of genes in phages is a problem for current gene annotators, we exploit this property by treating a phage genome as a network of paths: where open reading frames are favorable, and overlaps and gaps are less favorable, but still possible. We represent this network of connections as a weighted graph, and use dynamic programing to find the optimal path.

RESULTS

We compare PHANOTATE to other gene callers by annotating a set of 2133 complete phage genomes from GenBank, using PHANOTATE and the three most popular gene callers. We found that the four programs agree on 82% of the total predicted genes, with PHANOTATE predicting more genes than the other three. We searched for these extra genes in both GenBank's non-redundant protein database and all of the metagenomes in the sequence read archive, and found that they are present at levels that suggest that these are functional protein-coding genes.

AVAILABILITY AND IMPLEMENTATION

https://github.com/deprekate/PHANOTATE.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Global phylogeography and ancient evolution of the widespread human gut virus crAssphage

Microbiomes are vast communities of microorganisms and viruses that populate all natural ecosystems. Viruses have been considered to be the most variable component of microbiomes, as supported by virome surveys and examples of high genomic mosaicism. However, recent evidence suggests that the human gut virome is remarkably stable compared with that of other environments. Here, we investigate the origin, evolution and epidemiology of crAssphage, a widespread human gut virus. Through a global collaboration, we obtained DNA sequences of crAssphage from more than one-third of the world's countries and showed that the phylogeography of crAssphage is locally clustered within countries, cities and individuals. We also found fully colinear crAssphage-like genomes in both Old-World and New-World primates, suggesting that the association of crAssphage with primates may be millions of years old. Finally, by exploiting a large cohort of more than 1,000 individuals, we tested whether crAssphage is associated with bacterial taxonomic groups of the gut microbiome, diverse human health parameters and a wide range of dietary factors. We identified strong correlations with different clades of bacteria that are related to Bacteroidetes and weak associations with several diet categories, but no significant association with health or disease. We conclude that crAssphage is a benign cosmopolitan virus that may have coevolved with the human lineage and is an integral part of the normal human gut virome.

Global phylogeography and ancient evolution of the widespread human gut virus crAssphage

Microbiomes are vast communities of microorganisms and viruses that populate all natural ecosystems. Viruses have been considered to be the most variable component of microbiomes, as supported by virome surveys and examples of high genomic mosaicism. However, recent evidence suggests that the human gut virome is remarkably stable compared with that of other environments. Here, we investigate the origin, evolution and epidemiology of crAssphage, a widespread human gut virus. Through a global collaboration, we obtained DNA sequences of crAssphage from more than one-third of the world's countries and showed that the phylogeography of crAssphage is locally clustered within countries, cities and individuals. We also found fully colinear crAssphage-like genomes in both Old-World and New-World primates, suggesting that the association of crAssphage with primates may be millions of years old. Finally, by exploiting a large cohort of more than 1,000 individuals, we tested whether crAssphage is associated with bacterial taxonomic groups of the gut microbiome, diverse human health parameters and a wide range of dietary factors. We identified strong correlations with different clades of bacteria that are related to Bacteroidetes and weak associations with several diet categories, but no significant association with health or disease. We conclude that crAssphage is a benign cosmopolitan virus that may have coevolved with the human lineage and is an integral part of the normal human gut virome.

A diversity-generating retroelement encoded by a globally ubiquitous Bacteroides phage

BACKGROUND

Diversity-generating retroelements (DGRs) are genetic cassettes that selectively mutate target genes to produce hypervariable proteins. First characterized in Bordetella bacteriophage BPP-1, the DGR creates a hypervariable phage tail fiber that enables host tropism switching. Subsequent surveys for DGRs conclude that the majority identified to date are bacterial or archaeal in origin. This work examines bacteriophage and bacterial genomes for novel phage-encoded DGRs.

RESULTS

This survey discovered 92 DGRs that were only found in phages exhibiting a temperate lifestyle. The majority of phage-encoded DGRs were identified as prophages in bacterial hosts from the phyla Bacteroidetes, Proteobacteria, and Firmicutes. Sequence reads from these previously unidentified prophages were present in viral metagenomes (viromes), indicating these prophages can produce functional viruses. Five phages possessed hypervariable proteins with structural similarity to the tail fiber of BPP-1, whereas the functions of the remaining DGR target proteins were unknown. A novel temperate phage that harbors a DGR cassette targeting a protein of unknown function was induced from Bacteroides dorei. This phage, here named Bacteroides dorei Hankyphage, lysogenizes 13 different Bacteroides species and was present in 34% and 21% of whole-community metagenomes and human-associated viromes, respectively.

CONCLUSIONS

Here, the number of known DGR-containing phages is increased from four to 92. All of these phages exhibit a temperate lifestyle, including a cosmopolitan human-associated phage. Targeted hypervariation by temperate phages may be a ubiquitous mechanism underlying phage-bacteria interaction in the human microbiome.

Phage Genome Annotation Using the RAST Pipeline

Phages are complex biomolecular machineries that have to survive in a bacterial world. Phage genomes show many adaptations to their lifestyle such as shorter genes, reduced capacity for redundant DNA sequences, and the inclusion of tRNAs in their genomes. In addition, phages are not free-living, they require a host for replication and survival. These unique adaptations provide challenges for the bioinformatics analysis of phage genomes. In particular, ORF calling, genome annotation, noncoding RNA (ncRNA) identification, and the identification of transposons and insertions are all complicated in phage genome analysis. We provide a road map through the phage genome annotation pipeline, and discuss the challenges and solutions for phage genome annotation as we have implemented in the rapid annotation using subsystems (RAST) pipeline.

Computational approaches to predict bacteriophage-host relationships

Metagenomics has changed the face of virus discovery by enabling the accurate identification of viral genome sequences without requiring isolation of the viruses. As a result, metagenomic virus discovery leaves the first and most fundamental question about any novel virus unanswered: What host does the virus infect? The diversity of the global virosphere and the volumes of data obtained in metagenomic sequencing projects demand computational tools for virus–host prediction. We focus on bacteriophages (phages, viruses that infect bacteria), the most abundant and diverse group of viruses found in environmental metagenomes. By analyzing 820 phages with annotated hosts, we review and assess the predictive power of in silico phage–host signals. Sequence homology approaches are the most effective at identifying known phage–host pairs. Compositional and abundance-based methods contain significant signal for phage–host classification, providing opportunities for analyzing the unknowns in viral metagenomes. Together, these computational approaches further our knowledge of the interactions between phages and their hosts. Importantly, we find that all reviewed signals significantly link phages to their hosts, illustrating how current knowledge and insights about the interaction mechanisms and ecology of coevolving phages and bacteria can be exploited to predict phage–host relationships, with potential relevance for medical and industrial applications.

New viruses infecting bacteria are increasingly being discovered in many environments through sequence-based explorations. To understand their role in microbial ecosystems, computational tools are indispensable to prioritize and guide experimental efforts. This review assesses and discusses a range of bioinformatic approaches to predict bacteriophage–host relationships when all that is known is their genome sequence.

Altered hippocampal plasticity by prenatal kynurenine administration, kynurenine-3-monoxygenase (KMO) deletion or galantamine

Glutamate receptors sensitive to N-methyl-D-aspartate (NMDA) are involved in embryonic brain development but their activity may be modulated by the kynurenine pathway of tryptophan metabolism which includes an agonist (quinolinic acid) and an antagonist (kynurenic acid) at these receptors. Our previous work has shown that prenatal inhibition of the pathway produces abnormalities of brain development. In the present study kynurenine and probenecid (both 100mg/kg, doses known to increase kynurenic acid levels in the brain) were administered to female Wistar rats on embryonic days E14, E16 and E18 of gestation and the litter was allowed to develop to post-natal day P60. Western blotting revealed no changes in hippocampal expression of several proteins previously found to be altered by inhibition of the kynurenine pathway including the NMDA receptor subunits GluN1, GluN2A and GluN2B, as well as doublecortin, Proliferating Cell Nuclear Antigen (PCNA), sonic hedgehog and unco-ordinated (unc)-5H1 and 5H3. Mice lacking the enzyme kynurenine-3-monoxygenase (KMO) also showed no changes in hippocampal expression of several of these proteins or the 70-kDa and 100-kDa variants of Disrupted in Schizophrenia-1 (DISC1). Electrical excitability of pyramidal neurons in the CA1 region of hippocampal slices was unchanged, as was paired-pulse facilitation and inhibition. Long-term potentiation was decreased in the kynurenine-treated rats and in the KMO(-/-) mice, but galantamine reversed this effect in the presence of nicotinic receptor antagonists, consistent with evidence that it can potentiate glutamate at NMDA receptors. It is concluded that interference with the kynurenine pathway in utero can have lasting effects on brain function of the offspring, implying that the kynurenine pathway is involved in the regulation of early brain development.

GenomePeek-an online tool for prokaryotic genome and metagenome analysis

As more and more prokaryotic sequencing takes place, a method to quickly and accurately analyze this data is needed. Previous tools are mainly designed for metagenomic analysis and have limitations; such as long runtimes and significant false positive error rates. The online tool GenomePeek (edwards.sdsu.edu/GenomePeek) was developed to analyze both single genome and metagenome sequencing files, quickly and with low error rates. GenomePeek uses a sequence assembly approach where reads to a set of conserved genes are extracted, assembled and then aligned against the highly specific reference database. GenomePeek was found to be faster than traditional approaches while still keeping error rates low, as well as offering unique data visualization options.

A highly abundant bacteriophage discovered in the unknown sequences of human faecal metagenomes

Metagenomics, or sequencing of the genetic material from a complete microbial community, is a promising tool to discover novel microbes and viruses. Viral metagenomes typically contain many unknown sequences. Here we describe the discovery of a previously unidentified bacteriophage present in the majority of published human faecal metagenomes, which we refer to as crAssphage. Its ~97 kbp genome is six times more abundant in publicly available metagenomes than all other known phages together; it comprises up to 90% and 22% of all reads in virus-like particle (VLP)-derived metagenomes and total community metagenomes, respectively; and it totals 1.68% of all human faecal metagenomic sequencing reads in the public databases. The majority of crAssphage-encoded proteins match no known sequences in the database, which is why it was not detected before. Using a new co-occurrence profiling approach, we predict a Bacteroides host for this phage, consistent with Bacteroides-related protein homologues and a unique carbohydrate-binding domain encoded in the phage genome.

Metagenomic studies of microbial communities often report DNA sequences from unidentified viruses. Here, Dutilh et al. analyse metagenomic data to reveal the complete genome of an abundant, ubiquitous virus from human faeces, and predict that the virus infects bacteria of the Bacteroides group.

Multivariate analysis of functional metagenomes

Metagenomics is a primary tool for the description of microbial and viral communities. The sheer magnitude of the data generated in each metagenome makes identifying key differences in the function and taxonomy between communities difficult to elucidate. Here we discuss the application of seven different data mining and statistical analyses by comparing and contrasting the metabolic functions of 212 microbial metagenomes within and between 10 environments. Not all approaches are appropriate for all questions, and researchers should decide which approach addresses their questions. This work demonstrated the use of each approach: for example, random forests provided a robust and enlightening description of both the clustering of metagenomes and the metabolic processes that were important in separating microbial communities from different environments. All analyses identified that the presence of phage genes within the microbial community was a predictor of whether the microbial community was host-associated or free-living. Several analyses identified the subtle differences that occur with environments, such as those seen in different regions of the marine environment.

PHACTS, a computational approach to classifying the lifestyle of phages

Motivation: Bacteriophages have two distinct lifestyles: virulent and temperate. The virulent lifestyle has many implications for phage therapy, genomics and microbiology. Determining which lifestyle a newly sequenced phage falls into is currently determined using standard culturing techniques. Such laboratory work is not only costly and time consuming, but also cannot be used on phage genomes constructed from environmental sequencing. Therefore, a computational method that utilizes the sequence data of phage genomes is needed.

Results: Phage Classification Tool Set (PHACTS) utilizes a novel similarity algorithm and a supervised Random Forest classifier to make a prediction whether the lifestyle of a phage, described by its proteome, is virulent or temperate. The similarity algorithm creates a training set from phages with known lifestyles and along with the lifestyle annotation, trains a Random Forest to classify the lifestyle of a phage. PHACTS predictions are shown to have a 99% precision rate.

Availability and implementation: PHACTS was implemented in the PERL programming language and utilizes the FASTA program (Pearson and Lipman, 1988) and the R programming language library ‘Random Forest’ (Liaw and Weiner, 2010). The PHACTS software is open source and is available as downloadable stand-alone version or can be accessed online as a user-friendly web interface. The source code, help files and online version are available at http://www.phantome.org/PHACTS/.

Contact:katelyn@rohan.sdsu.edu; redwards@sciences.sdsu.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

Global changes in the hippocampal proteome following exposure to an enriched environment

Exposure to an enriched environment promotes neurochemical, structural and neurophysiological changes in the brain and is associated with enhanced synaptic plasticity and improved hippocampal-dependent learning. Using a global proteomics-based approach we have now been able to reveal the altered expression of a diverse range of hippocampal proteins following exposure to an enriched environment. Male Hooded Lister rats (8 weeks) were subjected to a 6-week regimen in which they were housed in either non-enriched (open field) or enriched conditions (toys, wheels etc.). Whole protein extracts from stratum pyramidale and stratum radiatum of area CA1 were then isolated and subjected to differential gel electrophoresis [McNair K, Davies CH, Cobb SR (2006) Plasticity-related regulation of the hippocampal proteome. Eur J Neurosci 23(2):575-580]. Of the 2469 resolvable protein spots detected in this study, 42 spots (1.7% of the detectable proteome) derived from predominantly somatic fractions and 32 proteins spots from dendritic fractions (1.3% of detectable proteome) were significantly altered in abundance following exposure to an enriched environment (somatic: 14 increased/28 decreased abundance, range -1.5 to +1.4-fold change; dendritic: 16 increased, 16 decreased abundance, range -1.6 to +3.0-fold change). Following in-gel tryptic digestion and Maldi-Tof/Q-star mass spectrometry, database searching revealed the identity of 50 protein spots displaying environmental enrichment-related modulation of expression. Identified proteins belonged to a variety of functional classes with gene ontology analysis revealing the majority (>70%) of regulated proteins to be part of the energy metabolism, cytoplasmic organization/biogenesis and signal transduction processes.

Regulation and function of spinal and peripheral neuronal B1 bradykinin receptors in inflammatory mechanical hyperalgesia

Activation of either B1 or B2 bradykinin receptors by kinins released from damaged tissues contributes to the development and maintenance of inflammatory hyperalgesia. Whereas B2 agonists activate sensory neurones directly, B1 agonists were thought only to have indirect actions on sensory neurones. The recent discovery of constitutive B1 receptor expression in the rat nervous system lead us to re-investigate the role of neuronal B1 receptors in inflammatory hyperalgesia. Therefore we have examined B1 bradykinin receptor regulation in rat dorsal root ganglia in a model of inflammatory hyperalgesia, and correlated it with hyperalgesic behaviour. Twenty-four hours after injection of Freund's complete adjuvant into one hindpaw, there was a significant increase in B1 protein expression (measured by immunohistochemistry) in both ipsilateral and contralateral dorsal root ganglion neurones, whereas axotomy resulted in reduction of B1 protein in ipsilateral dorsal root ganglia. In behavioural experiments, the B1 antagonist desArg10HOE140, administered by either intrathecal or systemic routes, attenuated Freund's complete adjuvant-induced mechanical hyperalgesia in the inflamed paw, but did not affect mechanical allodynia. The B1 agonist, desArg9BK, did not affect paw withdrawal thresholds in nai;ve rats following intraplantar administration into the paw, whilst intrathecal administration elicited mechanical hyperalgesia. However, after Freund's complete adjuvant-induced inflammation, desArg9BK caused a marked mechanical hyperalgesia, by either route, of the contralateral, uninflamed hindpaw, correlating with the observed contralateral and ipsilateral increases in receptor levels. Our results suggest a functional role for B1 receptors expressed both in the periphery and in the spinal cord, in mechanical hyperalgesia during inflammation.

The role of central and peripheral Cannabinoid1 receptors in the antihyperalgesic activity of cannabinoids in a model of neuropathic pain

We have examined the effects of cannabinoid agonists on hyperalgesia in a model of neuropathic pain in the rat and investigated the possible sites of action. The antihyperalgesic activity of the cannabinoids was compared with their ability to elicit behavioural effects characteristic of central cannabinoid activity. WIN55,212-2 (0.3-10 mg kg(-1)), CP-55,940 (0.03-1 mg kg(-1)) and HU-210 (0.001-0.03 mg kg(-1)) were all active in a 'tetrad' of tests consisting of tail-flick, catalepsy, rotarod and hypothermia following subcutaneous administration, with a rank order of potency in each of HU-210 > CP-55,940 > WIN55,212-2. The effects of WIN55,212-2 in each assay were blocked by the Cannabinoid1 (CB1) antagonist SR141716A. In the partial sciatic ligation model of neuropathic pain WIN55,212-2, CP-55,940 and HU-210 produced complete reversal of mechanical hyperalgesia within 3 h of subcutaneous administration with D50 values of 0.52, 0.08 and 0.005 mg kg(-1), respectively. In this model WIN55,212-2 was also effective against thermal hyperalgesia and mechanical allodynia. WIN55,212-2 produced pronounced reversal of mechanical hyperalgesia following intrathecal administration that was blocked by the CB1 antagonist SR141716A. Following intraplantar administration into the ipsilateral hindpaw, WIN55,212-2 produced up to 70% reversal of mechanical hyperalgesia, although activity was also observed at high doses following injection into the contralateral paw. The antihyperalgesic effect of WIN55,212-2 injected into the ipsilateral paw was blocked by subcutaneously administered SR141716A, but was not affected by intrathecally administered SR141716A. These data show that cannabinoids are highly potent and efficacious antihyperalgesic agents in a model of neuropathic pain. This activity is likely to be mediated via an action in both the CNS and in the periphery.