Improving gene annotation of complete viral genomes

Microviridae bacteriophages influence behavioural hallmarks of food addiction via tryptophan and tyrosine signalling pathways

Food addiction contributes to the obesity pandemic, but the connection between how the gut microbiome is linked to food addiction remains largely unclear. Here we show that Microviridae bacteriophages, particularly Gokushovirus WZ-2015a, are associated with food addiction and obesity across multiple human cohorts. Further analyses reveal that food addiction and Gokushovirus are linked to serotonin and dopamine metabolism. Mice receiving faecal microbiota and viral transplantation from human donors with the highest Gokushovirus load exhibit increased food addiction along with changes in tryptophan, serotonin and dopamine metabolism in different regions of the brain, together with alterations in dopamine receptors. Mechanistically, targeted tryptophan analysis shows lower anthranilic acid (AA) concentrations associated with Gokushovirus. AA supplementation in mice decreases food addiction and alters pathways related to the cycle of neurotransmitter synthesis release. In Drosophila, AA regulates feeding behaviour and addiction-like ethanol preference. In summary, this study proposes that bacteriophages in the gut microbiome contribute to regulating food addiction by modulating tryptophan and tyrosine metabolism.

Tegument proteins of Epstein-Barr virus: Diverse functions, complex networks, and oncogenesis

The tegument is the structure between the envelope and nucleocapsid of herpesvirus particles. Viral (and cellular) proteins accumulate to create the layers of the tegument. Some Epstein-Barr virus (EBV) tegument proteins are conserved widely in Herpesviridae, but others are shared only by members of the gamma-herpesvirus subfamily. As the interface to envelope and nucleocapsid, the tegument functions in virion morphogenesis and budding of the nucleocapsid during progeny production. When a virus particle enters a cell, enzymes such as kinase and deubiquitinase, and transcriptional activators are released from the virion to promote virus infection. Moreover, some EBV tegument proteins are involved in oncogenesis. Here, we summarize the roles of EBV tegument proteins, in comparison to those of other herpesviruses.

Caudovirales bacteriophages are associated with improved executive function and memory in flies, mice, and humans

Growing evidence implicates the gut microbiome in cognition. Viruses, the most abundant life entities on the planet, are a commonly overlooked component of the gut virome, dominated by the Caudovirales and Microviridae bacteriophages. Here, we show in a discovery (n = 114) and a validation cohort (n = 942) that subjects with increased Caudovirales and Siphoviridae levels in the gut microbiome had better performance in executive processes and verbal memory. Conversely, increased Microviridae levels were linked to a greater impairment in executive abilities. Microbiota transplantation from human donors with increased specific Caudovirales (>90% from the Siphoviridae family) levels led to increased scores in the novel object recognition test in mice and up-regulated memory-promoting immediate early genes in the prefrontal cortex. Supplementation of the Drosophila diet with the 936 group of lactococcal Siphoviridae bacteriophages resulted in increased memory scores and upregulation of memory-involved brain genes. Thus, bacteriophages warrant consideration as novel actors in the microbiome-brain axis.

Manual Annotation Studio (MAS): a collaborative platform for manual functional annotation of viral and microbial genomes

Background

Functional genome annotation is the process of labelling functional genomic regions with descriptive information. Manual curation can produce higher quality genome annotations than fully automated methods. Manual annotation efforts are time-consuming and complex; however, software can help reduce these drawbacks.

Results

We created Manual Annotation Studio (MAS) to improve the efficiency of the process of manual functional annotation prokaryotic and viral genomes. MAS allows users to upload unannotated genomes, provides an interface to edit and upload annotations, tracks annotation history and progress, and saves data to a relational database. MAS provides users with pertinent information through a simple point and click interface to execute and visualize results for multiple homology search tools (blastp, rpsblast, and HHsearch) against multiple databases (Swiss-Prot, nr, CDD, PDB, and an internally generated database). MAS was designed to accept connections over the local area network (LAN) of a lab or organization so multiple users can access it simultaneously. MAS can take advantage of high-performance computing (HPC) clusters by interfacing with SGE or SLURM and data can be exported from MAS in a variety of formats (FASTA, GenBank, GFF, and excel).

Conclusions

MAS streamlines and provides structure to manual functional annotation projects. MAS enhances the ability of users to generate, interpret, and compare results from multiple tools. The structure that MAS provides can improve project organization and reduce annotation errors. MAS is ideal for team-based annotation projects because it facilitates collaboration.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08029-8.

A Rare, Virulent Clostridium perfringens Bacteriophage Susfortuna Is the First Isolated Bacteriophage in a New Viral Genus

Background: Clostridium perfringens is a well known swine pathogen. C. perfringens type A is considered the causative agent of enteric diseases in neonatal and weaned piglets. Phage therapy using C. perfringens phages in vivo has previously proved effective. Materials and Methods: Pig fecal samples were used to isolate phages, with Clostridium perfringens type A as host. Complete genome sequencing, comparative genomics, a proteome analysis and electron microscopy were used to characterize the phage. Results: Clostridium phage Susfortuna has a double-stranded DNA genome of 19,046 bp with a G+C% content of 29.2, inverted terminal repeats and 28 predicted coding sequences (CDSs). Putative functions could not be assigned to most of the CDSs (64.3%). Transmission electron microscopy of phage Susfortuna revealed an isometric head and a short protruding tail stub resembling the structure of the Podoviridae family. A proteome analysis of phage Susfortuna identified seven structural proteins, but only one could be assigned with a putative function. Conclusions: Based on the morphology, the inverted terminal repeats and the small genome size, phage Susfortuna belongs to subfamily Picovirinae within the Podoviridae family in the order Caudovirales. Together with C. perfringens bacteriophage CPD7, phage Susfortuna represent a new genus of bacteriophages with very limited DNA sequence similarity to other known C. perfringens phages. Despite the limited DNA sequence similarity, the gene synteny among putative structural genes of phage Susfortuna is conserved among several C. perfringens bacteriophages belonging to the Podoviridae family indicating a common ancestor.

Comparative genomic analysis of 142 bacteriophages infecting Salmonella enterica subsp. enterica

BACKGROUND

Bacteriophages are bacterial parasites and are considered the most abundant and diverse biological entities on the planet. Previously we identified 154 prophages from 151 serovars of Salmonella enterica subsp. enterica. A detailed analysis of Salmonella prophage genomics is required given the influence of phages on their bacterial hosts and should provide a broader understanding of Salmonella biology and virulence and contribute to the practical applications of phages as vectors and antibacterial agents.

RESULTS

Here we provide a comparative analysis of the full genome sequences of 142 prophages of Salmonella enterica subsp. enterica which is the full complement of the prophages that could be retrieved from public databases. We discovered extensive variation in genome sizes (ranging from 6.4 to 358.7 kb) and guanine plus cytosine (GC) content (ranging from 35.5 to 65.4%) and observed a linear correlation between the genome size and the number of open reading frames (ORFs). We used three approaches to compare the phage genomes. The NUCmer/MUMmer genome alignment tool was used to evaluate linkages and correlations based on nucleotide identity between genomes. Multiple sequence alignment was performed to calculate genome average nucleotide identity using the Kalgin program. Finally, genome synteny was explored using dot plot analysis. We found that 90 phage genome sequences grouped into 17 distinct clusters while the remaining 52 genomes showed no close relationships with the other phage genomes and are identified as singletons. We generated genome maps using nucleotide and amino acid sequences which allowed protein-coding genes to be sorted into phamilies (phams) using the Phamerator software. Out of 5796 total assigned phamilies, one phamily was observed to be dominant and was found in 49 prophages, or 34.5% of the 142 phages in our collection. A majority of the phamilies, 4330 out of 5796 (74.7%), occurred in just one prophage underscoring the high degree of diversity among Salmonella bacteriophages.

CONCLUSIONS

Based on nucleotide and amino acid sequences, a high diversity was found among Salmonella bacteriophages which validate the use of prophage sequence analysis as a highly discriminatory subtyping tool for Salmonella. Thorough understanding of the conservation and variation of prophage genomic characteristics will facilitate their rational design and use as tools for bacterial strain construction, vector development and as anti-bacterial agents.

Complete Genome Sequence of Escherichia-Infecting Phage CEC_KAZ_2018, Isolated from Soil

Avian pathogenic Escherichia coli (APEC) bacteria are one of the main problems of the poultry industry. An effective way to combat colibacillosis is to use a phage preparation that lyses the bacteria. Here, we report the isolation of an E. coli-infecting phage, CEC_KAZ_2018, isolated from soil.

Avian pathogenic Escherichia coli (APEC) bacteria are one of the main problems of the poultry industry. An effective way to combat colibacillosis is to use a phage preparation that lyses the bacteria. Here, we report the isolation of an E. coli-infecting phage, CEC_KAZ_2018, isolated from soil.

Analysis of genetic recombination and the pan-genome of a highly recombinogenic bacteriophage species

Bacteriophages are the most prevalent biological entities impacting on the ecosystem and are characterized by their extensive diversity. However, there are two aspects of phages that have remained largely unexplored: genetic flux by recombination between phage populations and characterization of specific phages in terms of the pan-genome. Here, we examined the recombination and pan-genome in Helicobacter pylori prophages at both the genome and gene level. In the genome-level analysis, we applied, for the first time, chromosome painting and fineSTRUCTURE algorithms to a phage species, and showed novel trends in inter-population genetic flux. Notably, hpEastAsia is a phage population that imported a higher proportion of DNA fragments from other phages, whereas the hpSWEurope phages showed weaker signatures of inter-population recombination, suggesting genetic isolation. The gene-level analysis showed that, after parameter tuning of the prokaryote pan-genome analysis program, H. pylori phages have a pan-genome consisting of 75 genes and a soft-core genome of 10 genes, which includes genes involved in the lytic and lysogenic life cycles. Quantitative analysis of recombination events of the soft-core genes showed no substantial variation in the intensity of recombination across the genes, but rather equally frequent recombination among housekeeping genes that were previously reported to be less prone to recombination. The signature of frequent recombination appears to reflect the host–phage evolutionary arms race, either by contributing to escape from bacterial immunity or by protecting the host by producing defective phages.

Comparative Genomic Analysis of 130 Bacteriophages Infecting Bacteria in the Genus Pseudomonas

Bacteria of the genus Pseudomonas are genetically diverse and ubiquitous in the environment. Like other bacteria, those of the genus Pseudomonas are susceptible to bacteriophages which can significantly affect their host in many ways, ranging from cell lysis to major changes in morphology and virulence. Insights into phage genomes, evolution, and functional relationships with their hosts have the potential to contribute to a broader understanding of Pseudomonas biology, and the development of novel phage therapy strategies. Here we provide a broad-based comparative and evolutionary analysis of 130 complete Pseudomonas phage genome sequences available in online databases. We discovered extensive variation in genome size (ranging from 3 to 316 kb), G + C percentage (ranging from 37 to 66%), and overall gene content (ranging from 81–96% of genome space). Based on overall nucleotide similarity and the numbers of shared gene products, 100 out of 130 genome sequences were grouped into 12 different clusters; 30 were characterized as singletons, which do not have close relationships with other phage genomes. For 5/12 clusters, constituent phage members originated from two or more different Pseudomonas host species, suggesting that phage in these clusters can traverse bacterial species boundaries. An analysis of CRISPR spacers in Pseudomonas bacterial genome sequences supported this finding. Substantial diversity was revealed in analyses of phage gene families; out of 4,462 total families, the largest had only 39 members and there were 2,992 families with only one member. An evolutionary analysis of 72 phage gene families, based on patterns of nucleotide diversity at non-synonymous and synonymous sites, revealed strong and consistent signals for purifying selection. Our study revealed highly diverse and dynamic Pseudomonas phage genomes, and evidence for a dominant role of purifying selection in shaping the evolution of genes encoded in them.

The Kaposi's-sarcoma-associated herpesvirus orf35 gene product is required for efficient lytic virus reactivation

Kaposi's sarcoma-associated herpesvirus (KSHV) is implicated in the etiology of several human malignancies. KSHV open reading frame (orf) 35 encodes a conserved gammaherpesvirus protein with an, as yet, unknown function. Employing the bacterial artificial chromosome (BAC) system, we generated a recombinant viral clone that fails to express ORF35 (BAC16-ORF35-stop) but preserves intact adjacent and overlapping reading frames. Using this construct, we studied the role of this previously uncharacterized gene product during lytic reactivation of KSHV. Upon lytic reactivation, the ORF35-stop recombinant virus displayed significantly reduced lytic viral gene expression, viral DNA replication, and progeny virus production as compared to control wild-type virus. Exogenous expression of ORF35-Flag reversed the effects of ORF35 deficiency. These results demonstrate that ORF35 is important for efficient lytic virus reactivation.

Murine gammaherpesvirus 68 ORF35 is required for efficient lytic replication and latency

Murine gammaherpesvirus (MHV) 68, a natural pathogen of field mice, is related to human gammaherpesviruses, Epstein–Barr virus (EBV; human herpesvirus 4) and Kaposi’s sarcoma-associated herpesvirus (KSHV; human herpesvirus 8). The ORF35 of MHV-68 and its homologues of EBV and KSHV are located in the gene cluster composed of ORF34–ORF38 in which each gene overlaps with adjacent genes. Although MHV-68 ORF35 was reported to be an essential gene, its function during infection is presently unknown. In this study, we show, by analysing ORF35-transfected cells, that three serine residues in the C terminus are responsible for the phosphorylation and that the ORF35 protein forms homo-oligomers via a predicted coiled-coil motif. The ORF35 protein expressed by transfection was preferentially located in the cytoplasm of cells uninfected or infected with MHV-68. The recombinant virus lacking ORF35 (35S virus) exhibited genome replication and expression of lytic proteins comparable to those of the WT virus, but reduced levels of virus production, suggesting that the ORF35 protein acts at the virion assembly and/or egress step. Lytic replication in the lung after intranasal infection and the frequency of ex vivo reactivation from latency after intraperitoneal infection were lower in 35S virus-infected mice than in mice infected with the WT or marker-reverted virus. Our results indicate that ORF35 is not essential for MHV-68 lytic replication, but plays an important role in efficient viral replication and reactivation from latency.

Roles of Epstein-Barr virus BGLF3.5 gene and two upstream open reading frames in lytic viral replication in HEK293 cells

The Epstein-Barr virus (EBV) predominantly establishes a latent infection in B lymphocytes, but a small percentage of infected cells switch from the latent state to the lytic cycle, leading to potent viral DNA replication and progeny viruses production. We here focused on a lytic gene BGLF3.5, and first established BGLF3.5 mutants by marker cassette insertion. Unexpectedly, this insertion mutant failed to produce BGLF4 protein and thus progeny production was severely inhibited. Then we carefully made two point mutant viruses (stop codon insertion or frame-shift mutation) and found that BGLF3.5 is not essential for EBV lytic replication processes, such as viral gene expression, DNA replication, or progeny production in the HEK293 cells although its homolog in murine gammaherpesvirus 68 (MHV-68) was reported to be essential. In addition, we examined the roles of two short, upstream open reading frames within the 5'UTR of BGLF3.5 gene in translation of BGLF4.

NCBI viral genomes resource

Recent technological innovations have ignited an explosion in virus genome sequencing that promises to fundamentally alter our understanding of viral biology and profoundly impact public health policy. Yet, any potential benefits from the billowing cloud of next generation sequence data hinge upon well implemented reference resources that facilitate the identification of sequences, aid in the assembly of sequence reads and provide reference annotation sources. The NCBI Viral Genomes Resource is a reference resource designed to bring order to this sequence shockwave and improve usability of viral sequence data. The resource can be accessed at http://www.ncbi.nlm.nih.gov/genome/viruses/ and catalogs all publicly available virus genome sequences and curates reference genome sequences. As the number of genome sequences has grown, so too have the difficulties in annotating and maintaining reference sequences. The rapid expansion of the viral sequence universe has forced a recalibration of the data model to better provide extant sequence representation and enhanced reference sequence products to serve the needs of the various viral communities. This, in turn, has placed increased emphasis on leveraging the knowledge of individual scientific communities to identify important viral sequences and develop well annotated reference virus genome sets.

Analysis of the genetic diversity of ovine herpesvirus 2 in samples from livestock with malignant catarrhal fever

In order to define better virus isolates from animals with malignant catarrhal fever (MCF), segments of three genes of ovine herpesvirus-2 were amplified from diagnostic samples representing MCF cases with a range of clinical presentations in cattle, including head and eye, alimentary and neurological. The variation within each gene segment was estimated by DNA sequencing, which confirmed that the newly-annotated Ov9.5 gene was significantly more polymorphic than either of the other loci tested (segments of ORF50 and ORF75), with alleles that differed at over 60% of nucleotide positions. Despite this, the nine Ov9.5 alleles characterised had identical predicted splicing patterns and could be translated into Ov9.5 polypeptides with at least 49% amino acid identity. This multi-locus approach has potential for use in epidemiological studies and in charactering chains of infection. However there was no association between specific variants of OvHV-2 and the clinical/pathological presentation of MCF in the cattle analysed.

Gene identification in prokaryotic genomes, phages, metagenomes, and EST sequences with GeneMarkS suite

This unit describes how to use several gene-finding programs from the GeneMark line developed for finding protein-coding ORFs in genomic DNA of prokaryotic species, in genomic DNA of eukaryotic species with intronless genes, in genomes of viruses and phages, and in prokaryotic metagenomic sequences, as well as in EST sequences with spliced-out introns. These bioinformatics tools were demonstrated to have state-of-the-art accuracy, and have been frequently used for gene annotation in novel nucleotide sequences. An additional advantage of these sequence-analysis tools is that the problem of algorithm parameterization is solved automatically, with parameters estimated by iterative self-training (unsupervised training).

A novel spliced gene in alcelaphine herpesvirus 1 encodes a glycoprotein which is secreted in vitro

Herpesviruses often contain cryptic, spliced genes that are not obvious from the initial in silico annotation. Alcelaphine herpesvirus 1 (AlHV-1) contains 72 annotated ORFs but there are also a number of gaps between these that may have protein-coding potential. Comparative analysis of coding potential between AlHV-1 and the related ovine herpesvirus 2 (OvHV-2) revealed a putative novel spliced gene that we have termed A9.5. Analysis of cDNA clones from AlHV-1-infected cells revealed three overlapping clones corresponding to A9.5 and the coding sequence was confirmed by reverse transcription PCR of RNA from AlHV-1-infected cattle tissues. The A9.5 gene was predicted to encode a secreted glycoprotein with molecular mass 19 kDa. Empirical analysis showed that a recombinant haemagglutinin-tagged A9.5 fusion protein was secreted from transfected cells and had a molecular mass of 45 kDa, which was reduced to 20 kDa by endoglycosidase F treatment, confirming that A9.5 was a secreted glycoprotein. In situ RNA hybridization showed that A9.5 was expressed in cells associated with malignant catarrhal fever (MCF) lesions in infected cattle. Detailed analysis of the available OvHV-2 sequences revealed an homologous gene (Ov9.5) with conserved splicing signals and predicted amino acid sequence features in both sequenced isolates of this related virus. We have therefore identified a novel spliced gene in two related macaviruses that is expressed in MCF lesions. Future work will determine its importance for the pathogenesis of disease.

Viral infection: an evolving insight into the signal transduction pathways responsible for the innate immune response

The innate immune response is initiated by the interaction of stereotypical pathogen components with genetically conserved receptors for extracytosolic pathogen-associated molecular patterns (PAMPs) or intracytosolic nucleic acids. In multicellular organisms, this interaction typically clusters signal transduction molecules and leads to their activations, thereby initiating signals that activate innate immune effector mechanisms to protect the host. In some cases programmed cell death-a fundamental form of innate immunity-is initiated in response to genotoxic or biochemical stress that is associated with viral infection. In this paper we will summarize innate immune mechanisms that are relevant to viral pathogenesis and outline the continuing evolution of viral mechanisms that suppress the innate immunity in mammalian hosts. These mechanisms of viral innate immune evasion provide significant insight into the pathways of the antiviral innate immune response of many organisms. Examples of relevant mammalian innate immune defenses host defenses include signaling to interferon and cytokine response pathways as well as signaling to the inflammasome. Understanding which viral innate immune evasion mechanisms are linked to pathogenesis may translate into therapies and vaccines that are truly effective in eliminating the morbidity and mortality associated with viral infections in individuals.

Gene identification in prokaryotic genomes, phages, metagenomes, and EST sequences with GeneMarkS suite

This unit describes how to use several gene-finding programs from the GeneMark line developed for finding protein-coding ORFs in genomic DNA of prokaryotic species, in genomic DNA of eukaryotic species with intronless genes, in genomes of viruses and phages, and in prokaryotic metagenomic sequences, as well as in EST sequences with spliced-out introns. These bioinformatics tools were demonstrated to have state-of-the-art accuracy and have been frequently used for gene annotation in novel nucleotide sequences. An additional advantage of these sequence-analysis tools is that the problem of algorithm parameterization is solved automatically, with parameters estimated by iterative self-training (unsupervised training).

Identification, sequencing and molecular analysis of Chp4, a novel chlamydiaphage of Chlamydophila abortus belonging to the family Microviridae

Members of the family Microviridae have been identified in a number of chlamydial species infecting humans (phage CPAR39 in Chlamydophila pneumoniae), other mammals (φCPG1 in Chlamydophila caviae, Chp2 in Chlamydophila abortus and Chp3 in Chlamydophila pecorum) and birds (Chp1 in Chlamydophila psittaci). This study describes the identification and genome sequencing of Chp4, an icosahedral, 4530 bp, ssDNA phage in C. abortus. Chp4 is predicted to contain eight ORFs, six of which could be assigned putative functions based on sequence similarity to characterized bacteriophage. Gene order and content were highly conserved amongst chlamydiaphage, with the highest sequence variability occurring in the IN5 and INS variable regions of the VP1 major coat protein, which has been associated with host cell recognition and binding. Phylogenetic analysis of VP1 indicated that Chp4 is a member of the Chlamydiamicrovirus, and is most closely related to phage φCPG1 and CPAR39.

Novel DNA virus isolated from samples showing endothelial cell necrosis in the Japanese eel, Anguilla japonica

Economic loss due to viral endothelial cell necrosis of eel (VECNE) of Anguilla japonica is a serious problem for the cultured Japanese eel market. However, the viral genome responsible for VECNE is unknown. We recently developed a rapid determination system for viral nucleic acid sequences (RDV) to determine viral genome sequences. In this study, viral DNA fragments were obtained using RDV, and approximately 15-kbp circular full genome sequences were determined using a next-generation sequencing system, overlapping PCR, and Southern blot analysis. One open reading frame (ORF) was homologous to the large T-antigen of polyomavirus; other ORFs have no homology with any nucleic or amino acid sequences of polyomavirus. Therefore, as this DNA virus might comprise a novel virus family, we provisionally named it Japanese eel endothelial cells-infecting virus (JEECV). JEECV was detected in both naturally and experimentally infected eels, suggesting that JEECV potentially causes VECNE.

Towards Viral Genome Annotation Standards, Report from the 2010 NCBI Annotation Workshop

Improvements in DNA sequencing technologies portend a new era in virology and could possibly lead to a giant leap in our understanding of viral evolution and ecology. Yet, as viral genome sequences begin to fill the world’s biological databases, it is critically important to recognize that the scientific promise of this era is dependent on consistent and comprehensive genome annotation. With this in mind, the NCBI Genome Annotation Workshop recently hosted a study group tasked with developing sequence, function, and metadata annotation standards for viral genomes. This report describes the issues involved in viral genome annotation and reviews policy recommendations presented at the NCBI Annotation Workshop.

Redefining the genetics of murine gammaherpesvirus 68 via transcriptome-based annotation

Viral genetic studies typically focus on large open reading frames (ORFs) identified during genome annotation (ORF-based annotation). Here we describe tools for examining viral gene expression nucleotide by nucleotide across the genome. Using these tools on the 119,450 base pair (bp) genome of murine gammaherpesvirus 68 (gammaHV68) allowed us to establish that gammaHV68 RNA expression was significantly more complex than predicted from ORF-based annotation, including over 73,000 nucleotides of unexpected transcription within 30 expressed genomic regions (EGRs). Approximately 90% of this RNA expression was antisense to genomic regions containing known large ORFs. We verified the existence of previously undefined transcripts in three EGRs and determined which parts of the transcriptome depend on protein or viral DNA synthesis. This study redefines the genetic map of gammaHV68, indicating that herpesviruses contain significantly more genetic complexity than predicted from ORF-based genome annotations, and provides alternative tools and approaches for viral genetic studies.

VIGOR, an annotation program for small viral genomes

Background

The decrease in cost for sequencing and improvement in technologies has made it easier and more common for the re-sequencing of large genomes as well as parallel sequencing of small genomes. It is possible to completely sequence a small genome within days and this increases the number of publicly available genomes. Among the types of genomes being rapidly sequenced are those of microbial and viral genomes responsible for infectious diseases. However, accurate gene prediction is a challenge that persists for decoding a newly sequenced genome. Therefore, accurate and efficient gene prediction programs are highly desired for rapid and cost effective surveillance of RNA viruses through full genome sequencing.

Results

We have developed VIGOR (Viral Genome ORF Reader), a web application tool for gene prediction in influenza virus, rotavirus, rhinovirus and coronavirus subtypes. VIGOR detects protein coding regions based on sequence similarity searches and can accurately detect genome specific features such as frame shifts, overlapping genes, embedded genes, and can predict mature peptides within the context of a single polypeptide open reading frame. Genotyping capability for influenza and rotavirus is built into the program. We compared VIGOR to previously described gene prediction programs, ZCURVE_V, GeneMarkS and FLAN. The specificity and sensitivity of VIGOR are greater than 99% for the RNA viral genomes tested.

Conclusions

VIGOR is a user friendly web-based genome annotation program for five different viral agents, influenza, rotavirus, rhinovirus, coronavirus and SARS coronavirus. This is the first gene prediction program for rotavirus and rhinovirus for public access. VIGOR is able to accurately predict protein coding genes for the above five viral types and has the capability to assign function to the predicted open reading frames and genotype influenza virus. The prediction software was designed for performing high throughput annotation and closure validation in a post-sequencing production pipeline.

Orthopoxvirus genome evolution: the role of gene loss

Poxviruses are highly successful pathogens, known to infect a variety of hosts. The family Poxviridae includes Variola virus, the causative agent of smallpox, which has been eradicated as a public health threat but could potentially reemerge as a bioterrorist threat. The risk scenario includes other animal poxviruses and genetically engineered manipulations of poxviruses. Studies of orthologous gene sets have established the evolutionary relationships of members within the Poxviridae family. It is not clear, however, how variations between family members arose in the past, an important issue in understanding how these viruses may vary and possibly produce future threats. Using a newly developed poxvirus-specific tool, we predicted accurate gene sets for viruses with completely sequenced genomes in the genus Orthopoxvirus. Employing sensitive sequence comparison techniques together with comparison of syntenic gene maps, we established the relationships between all viral gene sets. These techniques allowed us to unambiguously identify the gene loss/gain events that have occurred over the course of orthopoxvirus evolution. It is clear that for all existing Orthopoxvirus species, no individual species has acquired protein-coding genes unique to that species. All existing species contain genes that are all present in members of the species Cowpox virus and that cowpox virus strains contain every gene present in any other orthopoxvirus strain. These results support a theory of reductive evolution in which the reduction in size of the core gene set of a putative ancestral virus played a critical role in speciation and confining any newly emerging virus species to a particular environmental (host or tissue) niche.

Ab initio gene identification in metagenomic sequences

We describe an algorithm for gene identification in DNA sequences derived from shotgun sequencing of microbial communities. Accurate ab initio gene prediction in a short nucleotide sequence of anonymous origin is hampered by uncertainty in model parameters. While several machine learning approaches could be proposed to bypass this difficulty, one effective method is to estimate parameters from dependencies, formed in evolution, between frequencies of oligonucleotides in protein-coding regions and genome nucleotide composition. Original version of the method was proposed in 1999 and has been used since for (i) reconstructing codon frequency vector needed for gene finding in viral genomes and (ii) initializing parameters of self-training gene finding algorithms. With advent of new prokaryotic genomes en masse it became possible to enhance the original approach by using direct polynomial and logistic approximations of oligonucleotide frequencies, as well as by separating models for bacteria and archaea. These advances have increased the accuracy of model reconstruction and, subsequently, gene prediction. We describe the refined method and assess its accuracy on known prokaryotic genomes split into short sequences. Also, we show that as a result of application of the new method, several thousands of new genes could be added to existing annotations of several human and mouse gut metagenomes.

Experimental annotation of channel catfish virus by probabilistic proteogenomic mapping

Experimental identification of expressed proteins by proteomics constitutes the most reliable approach to identify genomic location and structure of protein-coding genes and substantially complements computational genome annotation. Channel catfish herpesvirus (CCV) is a simple comparative model for understanding herpesvirus biology and the evolution of the Herpesviridae. The canonical CCV genome has 76 predicted ORF and only 12 of these have been confirmed experimentally. We describe a modification of a statistical method, which assigns significance measures, q-values, to peptide identifications based on 2-D LC ESI MS/MS, real-decoy database searches and SEQUEST XCorr and DeltaC(n) scores. We used this approach to identify CCV proteins expressed during its replication in cell culture, to determine protein composition of mature virions and, consequently, to refine the canonical CCV genome annotation. To complement trypsin, we used partial proteinase K digestion, which yielded greater proteome coverage. At FDR <5%, for peptide identifications, we identified 25/76 previously predicted ORF using trypsin and 31/76 using proteinase K. Furthermore, we identified 17 novel protein-coding regions (7 potential ATG-initiated ORF). Most of these novel ORF encode small proteins (<100 amino acids). Directed, strand-specific reverse transcription real-time PCR confirmed RNA expression from 6/7 novel ATG-initiated ORF investigated.

Overlapping genes produce proteins with unusual sequence properties and offer insight into de novo protein creation

It is widely assumed that new proteins are created by duplication, fusion, or fission of existing coding sequences. Another mechanism of protein birth is provided by overlapping genes. They are created de novo by mutations within a coding sequence that lead to the expression of a novel protein in another reading frame, a process called "overprinting." To investigate this mechanism, we have analyzed the sequences of the protein products of manually curated overlapping genes from 43 genera of unspliced RNA viruses infecting eukaryotes. Overlapping proteins have a sequence composition globally biased toward disorder-promoting amino acids and are predicted to contain significantly more structural disorder than nonoverlapping proteins. By analyzing the phylogenetic distribution of overlapping proteins, we were able to confirm that 17 of these had been created de novo and to study them individually. Most proteins created de novo are orphans (i.e., restricted to one species or genus). Almost all are accessory proteins that play a role in viral pathogenicity or spread, rather than proteins central to viral replication or structure. Most proteins created de novo are predicted to be fully disordered and have a highly unusual sequence composition. This suggests that some viral overlapping reading frames encoding hypothetical proteins with highly biased composition, often discarded as noncoding, might in fact encode proteins. Some proteins created de novo are predicted to be ordered, however, and whenever a three-dimensional structure of such a protein has been solved, it corresponds to a fold previously unobserved, suggesting that the study of these proteins could enhance our knowledge of protein space.

In silico identification of genes in bacteriophage DNA

One of the most satisfying aspects of a genome sequencing project is the identification of the genes contained within it.These are of two types: those which encode tRNAs and those which produce proteins. After a general introduction on the properties of protein-encoding genes and the utility of the Basic Local Alignment Search Tool (BLASTX) to identify genes through homologs, a variety of tools are discussed by their creators. These include for genome annotation: GeneMark, Artemis, and BASys; and, for genome comparisons: Artemis Comparison Tool (ACT), Mauve, CoreGenes, and GeneOrder.

Malignant catarrhal fever: a review

Malignant catarrhal fever (MCF) is a fatal lymphoproliferative disease of cattle and other ungulates caused by the ruminant gamma-herpesviruses alcelaphine herpesvirus 1 (AlHV-1) and ovine herpesvirus 2 (OvHV-2). These viruses cause inapparent infection in their reservoir hosts (wildebeest for AlHV-1 and sheep for OvHV-2), but fatal lymphoproliferative disease when they infect MCF-susceptible hosts, including cattle, deer, bison, water buffalo and pigs. MCF is an important disease wherever reservoir and MCF-susceptible species mix and currently is a particular problem in Bali cattle in Indonesia, bison in the USA and in pastoralist cattle herds in Eastern and Southern Africa. MCF is characterised by the accumulation of lymphocytes (predominantly CD8(+) T lymphocytes) in a variety of organs, often associated with tissue necrosis. Only a small proportion of these lymphocytes appear to contain virus, although recent results with virus gene-specific probes indicate that more infected cells may be present than previously thought. The tissue damage in MCF is hypothesised to be caused by the indiscriminate activity of MHC-unrestricted cytotoxic T/natural killer cells. The pathogenesis of MCF and the virus life cycle are poorly understood and, currently, there is no effective disease control. Recent sequencing of the OvHV-2 genome and construction of an AlHV-1 bacterial artificial chromosome (BAC) are facilitating studies to understand the pathogenesis of this extraordinary disease. Furthermore, new and improved methods of disease diagnosis have been developed and promising vaccine strategies are being tested. The next few years are likely to be exciting and productive for MCF research.

Transcriptional regulation of the open reading frame 35 encoded by Kaposi's sarcoma-associated herpesvirus

Kaposi's sarcoma-associated herpesvirus (KSHV) is a member of the Gammaherpesvirinae and is causally associated with Kaposi's sarcoma, primary effusion lymphoma, and multicentric Castleman's disease. The KSHV genome encodes over 85 genes; the function of some is entirely unknown. We have characterized the transcriptional regulation of a conserved and uncharacterized Gammaherpesvirinae open reading frame, orf35, which lies in a cluster of several overlapping genes, orf34 to orf38. We identified the transcription start site and analyzed upstream sequences. We found that expression of the KSHV lytic replication and transcription activator (RTA) strongly increased the orf35 promoter activity through a 46-nucleotide region which includes a conserved AP-1 binding site. Electrophoretic mobility shift assay demonstrated direct binding of cJUN and cFOS to the predicted AP-1 binding site. Finally, using a mutated promoter lacking the AP-1 site and dominant-negative cFOS, we established that the RTA-mediated orf35 transactivation is AP-1-dependent.

Annotation of selection strengths in viral genomes

MOTIVATION

Viral genomes tend to code in overlapping reading frames to maximize informational content. This may result in atypical codon bias and particular evolutionary constraints. Due to the fast mutation rate of viruses, there is additional strong evidence for varying selection between intra- and intergenomic regions. The presence of multiple coding regions complicates the concept of K(a)/K(s) ratio, and thus begs for an alternative approach when investigating selection strengths. Building on the paper by McCauley and Hein, we develop a method for annotating a viral genome coding in overlapping reading frames. We introduce an evolutionary model capable of accounting for varying levels of selection along the genome, and incorporate it into our prior single sequence HMM methodology, extending it now to a phylogenetic HMM. Given an alignment of several homologous viruses to a reference sequence, we may thus achieve an annotation both of coding regions as well as selection strengths, allowing us to investigate different selection patterns and hypotheses.

RESULTS

We illustrate our method by applying it to a multiple alignment of four HIV2 sequences, as well as of three Hepatitis B sequences. We obtain an annotation of the coding regions, as well as a posterior probability for each site of the strength of selection acting on it. From this we may deduce the average posterior selection acting on the different genes. Whilst we are encouraged to see in HIV2, that the known to be conserved genes gag and pol are indeed annotated as such, we also discover several sites of less stringent negative selection within the env gene. To the best of our knowledge, we are the first to subsequently provide a full selection annotation of the Hepatitis B genome by explicitly modelling the evolution within overlapping reading frames, and not relying on simple K(a)/K(s) ratios.

Comparative annotation of viral genomes with non-conserved gene structure

MOTIVATION

Detecting genes in viral genomes is a complex task. Due to the biological necessity of them being constrained in length, RNA viruses in particular tend to code in overlapping reading frames. Since one amino acid is encoded by a triplet of nucleic acids, up to three genes may be coded for simultaneously in one direction. Conventional hidden Markov model (HMM)-based gene-finding algorithms may typically find it difficult to identify multiple coding regions, since in general their topologies do not allow for the presence of overlapping or nested genes. Comparative methods have therefore been restricted to likelihood ratio tests on potential regions as to being double or single coding, using the fact that the constrictions forced upon multiple-coding nucleotides will result in atypical sequence evolution. Exploiting these same constraints, we present an HMM based gene-finding program, which allows for coding in unidirectional nested and overlapping reading frames, to annotate two homologous aligned viral genomes. Our method does not insist on conserved gene structure between the two sequences, thus making it applicable for the pairwise comparison of more distantly related sequences.

RESULTS

We apply our method to 15 pairwise alignments of six different HIV2 genomes. Given sufficient evolutionary distance between the two sequences, we achieve sensitivity of approximately 84-89% and specificity of approximately 97-99.9%. We additionally annotate three pairwise alignments of the more distantly related HIV1 and HIV2, as well as of two different hepatitis viruses, attaining results of approximately 87% sensitivity and approximately 98.5% specificity. We subsequently incorporate prior knowledge by 'knowing' the gene structure of one sequence and annotating the other conditional on it. Boosting accuracy close to perfect we demonstrate that conservation of gene structure on top of nucleotide sequence is a valuable source of information, especially in distantly related genomes.

AVAILABILITY

The Java code is available from the authors.

Genetic diversity among five T4-like bacteriophages

BACKGROUND

Bacteriophages are an important repository of genetic diversity. As one of the major constituents of terrestrial biomass, they exert profound effects on the earth's ecology and microbial evolution by mediating horizontal gene transfer between bacteria and controlling their growth. Only limited genomic sequence data are currently available for phages but even this reveals an overwhelming diversity in their gene sequences and genomes. The contribution of the T4-like phages to this overall phage diversity is difficult to assess, since only a few examples of complete genome sequence exist for these phages. Our analysis of five T4-like genomes represents half of the known T4-like genomes in GenBank.

RESULTS

Here, we have examined in detail the genetic diversity of the genomes of five relatives of bacteriophage T4: the Escherichia coli phages RB43, RB49 and RB69, the Aeromonas salmonicida phage 44RR2.8t (or 44RR) and the Aeromonas hydrophila phage Aeh1. Our data define a core set of conserved genes common to these genomes as well as hundreds of additional open reading frames (ORFs) that are nonconserved. Although some of these ORFs resemble known genes from bacterial hosts or other phages, most show no significant similarity to any known sequence in the databases. The five genomes analyzed here all have similarities in gene regulation to T4. Sequence motifs resembling T4 early and late consensus promoters were observed in all five genomes. In contrast, only two of these genomes, RB69 and 44RR, showed similarities to T4 middle-mode promoter sequences and to the T4 motA gene product required for their recognition. In addition, we observed that each phage differed in the number and assortment of putative genes encoding host-like metabolic enzymes, tRNA species, and homing endonucleases.

CONCLUSION

Our observations suggest that evolution of the T4-like phages has drawn on a highly diverged pool of genes in the microbial world. The T4-like phages harbour a wealth of genetic material that has not been identified previously. The mechanisms by which these genes may have arisen may differ from those previously proposed for the evolution of other bacteriophage genomes.

Annotation of unknown yeast ORFs by correlation analysis of microarray data and extensive literature searches

Changes in the expression of genes were used to elucidate the metabolic pathways and regulatory mechanisms that respond to environmental or genetic modifications. Results from previously published chemostat datasets were merged with novel data generated in the present study. ORFs displaying significant changes in expression that correlated with those of other ORFs were analysed using GO mapping tools and supplemented by literature information. The strategy developed was used to propose annotations for ORFs of unknown function. The following ORFs were assigned functions as a result of this study: YMR090w, YGL157w, YGR243w, YLR327c, YER121w, YFR017c, YGR067c, YKL187c, YGR236c (SPG1), YMR107w (SPG4), YMR206w, YER067w, YJL103c, YNL175C (NOP13) YJL200C, YDL070C (FMP16) and YGR173W.

Detecting overlapping coding sequences in virus genomes

Background

Detecting new coding sequences (CDSs) in viral genomes can be difficult for several reasons. The typically compact genomes often contain a number of overlapping coding and non-coding functional elements, which can result in unusual patterns of codon usage; conservation between related sequences can be difficult to interpret – especially within overlapping genes; and viruses often employ non-canonical translational mechanisms – e.g. frameshifting, stop codon read-through, leaky-scanning and internal ribosome entry sites – which can conceal potentially coding open reading frames (ORFs).

Results

In a previous paper we introduced a new statistic – MLOGD (Maximum Likelihood Overlapping Gene Detector) – for detecting and analysing overlapping CDSs. Here we present (a) an improved MLOGD statistic, (b) a greatly extended suite of software using MLOGD, (c) a database of results for 640 virus sequence alignments, and (d) a web-interface to the software and database. Tests show that, from an alignment with just 20 mutations, MLOGD can discriminate non-overlapping CDSs from non-coding ORFs with a typical accuracy of up to 98%, and can detect CDSs overlapping known CDSs with a typical accuracy of 90%. In addition, the software produces a variety of statistics and graphics, useful for analysing an input multiple sequence alignment.

Conclusion

MLOGD is an easy-to-use tool for virus genome annotation, detecting new CDSs – in particular overlapping or short CDSs – and for analysing overlapping CDSs following frameshift sites. The software, web-server, database and supplementary material are available at .

ZCURVE_V: a new self-training system for recognizing protein-coding genes in viral and phage genomes

Background

It necessary to use highly accurate and statistics-based systems for viral and phage genome annotations. The GeneMark systems for gene-finding in virus and phage genomes suffer from some basic drawbacks. This paper puts forward an alternative approach for viral and phage gene-finding to improve the quality of annotations, particularly for newly sequenced genomes.

Results

The new system ZCURVE_V has been run for 979 viral and 212 phage genomes, respectively, and satisfactory results are obtained. To have a fair comparison with the currently available software of similar function, GeneMark, a total of 30 viral genomes that have not been annotated by GeneMark are selected to be tested. Consequently, the average specificity of both systems is well matched, however the average sensitivity of ZCURVE_V for smaller viral genomes (< 100 kb), which constitute the main parts of viral genomes sequenced so far, is higher than that of GeneMark. Additionally, for the genome of Amsacta moorei entomopoxvirus, probably with the lowest genomic GC content among the sequenced organisms, the accuracy of ZCURVE_V is much better than that of GeneMark, because the later predicts hundreds of false-positive genes. ZCURVE_V is also used to analyze well-studied genomes, such as HIV-1, HBV and SARS-CoV. Accordingly, the performance of ZCURVE_V is generally better than that of GeneMark. Finally, ZCURVE_V may be downloaded and run locally, particularly facilitating its utilization, whereas GeneMark is not downloadable. Based on the above comparison, it is suggested that ZCURVE_V may serve as a preferred gene-finding tool for viral and phage genomes newly sequenced. However, it is also shown that the joint application of both systems, ZCURVE_V and GeneMark, leads to better gene-finding results. The system ZCURVE_V is freely available at: .

Conclusion

ZCURVE_V may serve as a preferred gene-finding tool used for viral and phage genomes, especially for anonymous viral and phage genomes newly sequenced.

GeneMark: web software for gene finding in prokaryotes, eukaryotes and viruses

The task of gene identification frequently confronting researchers working with both novel and well studied genomes can be conveniently and reliably solved with the help of the GeneMark web software (http://opal.biology.gatech.edu/GeneMark/). The website provides interfaces to the GeneMark family of programs designed and tuned for gene prediction in prokaryotic, eukaryotic and viral genomic sequences. Currently, the server allows the analysis of nearly 200 prokaryotic and >10 eukaryotic genomes using species-specific versions of the software and pre-computed gene models. In addition, genes in prokaryotic sequences from novel genomes can be identified using models derived on the spot upon sequence submission, either by a relatively simple heuristic approach or by the full-fledged self-training program GeneMarkS. A database of reannotations of >1000 viral genomes by the GeneMarkS program is also available from the web site. The GeneMark website is frequently updated to provide the latest versions of the software and gene models.

Identification of proteins associated with murine cytomegalovirus virions

Proteins associated with the murine cytomegalovirus (MCMV) viral particle were identified by a combined approach of proteomic and genomic methods. Purified MCMV virions were dissociated by complete denaturation and subjected to either separation by sodium dodecyl sulfate-polyacrylamide gel electrophoresis and in-gel digestion or treated directly by in-solution tryptic digestion. Peptides were separated by nanoflow liquid chromatography and analyzed by tandem mass spectrometry (LC-MS/MS). The MS/MS spectra obtained were searched against a database of MCMV open reading frames (ORFs) predicted to be protein coding by an MCMV-specific version of the gene prediction algorithm GeneMarkS. We identified 38 proteins from the capsid, tegument, glycoprotein, replication, and immunomodulatory protein families, as well as 20 genes of unknown function. Observed irregularities in coding potential suggested possible sequence errors in the 3'-proximal ends of m20 and M31. These errors were experimentally confirmed by sequencing analysis. The MS data further indicated the presence of peptides derived from the unannotated ORFs ORF(c225441-226898) (m166.5) and ORF(105932-106072). Immunoblot experiments confirmed expression of m166.5 during viral infection.

Complete genomic nucleotide sequence of the temperate bacteriophage Aa Phi 23 of Actinobacillus actinomycetemcomitans

The entire double-stranded DNA genome of the Actinobacillus actinomycetemcomitans bacteriophage Aa Phi 23 was sequenced. Linear DNA contained in the phage particles is circularly permuted and terminally redundant. Therefore, the physical map of the phage genome is circular. Its size is 43,033 bp with an overall molar G+C content of 42.5 mol%. Sixty-six potential open reading frames (ORFs) were identified, including an ORF resulting from a translational frameshift. A putative function could be assigned to 23 of them. Twenty-three other ORFs share homologies only with hypothetical proteins present in several bacteria or bacteriophages, and 20 ORFs seem to be specific for phage Aa Phi 23. The organization of the phage genome and several genetic functions share extensive similarities to that of the lambdoid phages. However, Aa Phi 23 encodes a DNA adenine methylase, and the DNA packaging strategy is more closely related to the P22 system. The attachment sites of Aa Phi 23 (attP) and several A. actinomycetemcomitans hosts (attB) are 49 bp long.