Complete genomic sequence and comparative analysis of the tumorigenic poxvirus Yaba monkey tumor virus

The Yatapoxvirus genus of poxviruses is comprised of Yaba monkey tumor virus (YMTV), Tanapox virus, and Yaba-like disease virus (YLDV), which all have the ability to infect primates, including humans. Unlike other poxviruses, YMTV induces formation of focalized histiocytomas upon infection. To gain a greater understanding of the Yatapoxvirus genus and the unique tumor formation properties of YMTV, we sequenced the 134,721-bp genome of YMTV. The genome of YMTV encodes at least 140 open reading frames, all of which are also found as orthologs in the closely related YLDV. However, 13 open reading frames found in YLDV are completely absent from YMTV. Common to both YLDV and YMTV are the unusually large noncoding regions between many open reading frames. To determine whether any of these noncoding regions might be functionally significant, we carried out a comparative analysis between the putative noncoding regions of YMTV and similar noncoding regions from other poxviruses. This approach identified three new gene poxvirus families, defined as orthologs of YMTV23.5L, YMTV28.5L, and YMTV120.5L, which are highly conserved in virtually all poxvirus species. Furthermore, the comparative analysis also revealed a 40-bp nucleotide sequence at approximately 14,700 bases from the left terminus that was 100% identical in the comparable intergene site within members of the Yatapoxvirus, Suipoxvirus, and Capripoxvirus genera and 95% conserved in the Leporipoxvirus genus. This conserved sequence was shown to function as a poxvirus late promoter element in transfected and infected cells, but other functions, such as an involvement in viral replication or packaging, cannot be excluded. Finally, we summarize the predicted immunomodulatory protein repertoire in the Yatapoxvirus genus as a whole.     

SynBrowse: a synteny browser for comparative sequence analysis

MOTIVATION: The recent efforts of various sequence projects to sequence deeply into various phylogenies provide great resources for comparative sequence analysis. A generic and portable tool is essential for scientists to visualize and analyze sequence comparisons. RESULTS: We have developed SynBrowse, a synteny browser for visualizing and analyzing genome alignments both within and between species. It is intended to help scientists study macrosynteny, microsynteny and homologous genes between sequences. It can also aid with the identification of uncharacterized genes, putative regulatory elements and novel structural features of a species. SynBrowse is a GBrowse (the Generic Genome Browser) family software tool that runs on top of the open source BioPerl modules. It consists of two components: a web-based front end and a set of relational database back ends. Each database stores pre-computed alignments from a focus sequence to reference sequences in addition to the genome annotations of the focus sequence. The user interface lets end users select a key comparative alignment type and search for syntenic blocks between two sequences and zoom in to view the relationships among the corresponding genome annotations in detail. SynBrowse is portable with simple installation, flexible configuration, convenient data input and easy integration with other components of a model organism system. AVAILABILITY: The software is available at http://www.gmod.org CONTACT: vbrendel@iastate.edu     

Acetaminophen: a practical pharmacologic overview.

Acetaminophen is an effective analgesic and antipyretic agent with few adverse effects when used in recommended dosages. The drug is metabolized mainly in the liver, and the several end products have no harmful effects. An intermediate compound in a minor metabolic pathway, however, is toxic; it is normally inactivated by glutathione. In the case of an acetaminophen overdose the hepatic stores of glutathione seem to become depleted, leaving the toxic intermediate free to damage liver tissue. Such damage is unlikely to occur unless the plasma concentration of acetaminophen peaks above 150 micrograms/mL--a level far in excess of the 5 to 20 micrograms/mL achieved with therapeutic doses of the drug. Long-term therapeutic use of acetaminophen does not appear to be associated with liver damage, although some case reports suggest the possibility. Acetaminophen poisoning follows an acute overdose and, if untreated, is manifested clinically by an initial phase of nonspecific signs and symptoms, a latent period in which the liver transaminase levels rise and then, 3 to 5 days after the ingestion, signs of more serious hepatic dysfunction. Most patients do not progress beyond the first or second phase. They and those who survive the third phase recover with no residual injury to the liver. Appropriate antidotal therapy markedly reduces the severity of the initial damage.     

COGNAT: a web server for comparative analysis of genomic neighborhoods

Background

In prokaryotic genomes, functionally coupled genes can be organized in conserved gene clusters enabling their coordinated regulation. Such clusters could contain one or several operons, which are groups of co-transcribed genes. Those genes that evolved from a common ancestral gene by speciation (i.e. orthologs) are expected to have similar genomic neighborhoods in different organisms, whereas those copies of the gene that are responsible for dissimilar functions (i.e. paralogs) could be found in dissimilar genomic contexts. Comparative analysis of genomic neighborhoods facilitates the prediction of co-regulated genes and helps to discern different functions in large protein families.

Aim

We intended, building on the attribution of gene sequences to the clusters of orthologous groups of proteins (COGs), to provide a method for visualization and comparative analysis of genomic neighborhoods of evolutionary related genes, as well as a respective web server.

Results

Here we introduce the COmparative Gene Neighborhoods Analysis Tool (COGNAT), a web server for comparative analysis of genomic neighborhoods. The tool is based on the COG database, as well as the Pfam protein families database. As an example, we show the utility of COGNAT in identifying a new type of membrane protein complex that is formed by paralog(s) of one of the membrane subunits of the NADH:quinone oxidoreductase of type 1 (COG1009) and a cytoplasmic protein of unknown function (COG3002).

Reviewers

This article was reviewed by Drs. Igor Zhulin, Uri Gophna and Igor Rogozin.

     

GARSA: genomic analysis resources for sequence annotation

SUMMARY: Growth of genome data and analysis possibilities have brought new levels of difficulty for scientists to understand, integrate and deal with all this ever-increasing information. In this scenario, GARSA has been conceived aiming to facilitate the tasks of integrating, analyzing and presenting genomic information from several bioinformatics tools and genomic databases, in a flexible way. GARSA is a user-friendly web-based system designed to analyze genomic data in the context of a pipeline. EST and GGS data can be analyzed using the system since it accepts (1) chromatograms, (2) download of sequences from GenBank, (3) Fasta files stored locally or (4) a combination of all three. Quality evaluation of chromatograms, vector removing and clusterization are easily performed as part of the pipeline. A number of local and customizable Blast and CDD analyses can be performed as well as Interpro, complemented with phylogeny analyses. GARSA is being used for the analyses of Trypanosoma vivax (GSS and EST), Trypanosoma rangeli (GSS, EST and ORESTES), Bothrops jararaca (EST), Piaractus mesopotamicus (EST) and Lutzomyia longipalpis (EST). AVAILABILITY: The GARSA system is freely available under GPL license (http://www.biowebdb.org/garsa/). For download requests visit http://www.biowebdb.org/garsa/ or contact Dr Alberto Dávila.     

Modeling a minimal ribosome based on comparative sequence analysis

We have determined the three-dimensional organization of ribosomal RNAs and proteins essential for minimal ribosome function. Comparative sequence analysis identifies regions of the ribosome that have been evolutionarily conserved, and the spatial organization of conserved domains is determined by mapping these onto structures of the 30S and 50S subunits determined by X-ray crystallography. Several functional domains of the ribosome are conserved in their three-dimensional organization in the Archaea, Bacteria, Eucaryotic nuclear, mitochondria and chloroplast ribosomes. In contrast, other regions from both subunits have shifted their position in three-dimensional space during evolution, including the L11 binding domain and the alpha-sarcin-ricin loop (SRL). We examined conserved bridge interactions between the two ribosomal subunits, giving an indication of which contacts are more significant. The tRNA contacts that are conserved were also determined, highlighting functional interactions as the tRNA moves through the ribosome during protein synthesis. To augment these studies of a large collection of comparative structural models sampled from all major branches on the phylogenetic tree, Caenorhabditis elegans mitochondrial rRNA is considered individually because it is among the smallest rRNA sequences known. The C.elegans model supports the large collection of comparative structure models while providing insight into the evolution of mitochondrial ribosomes.     

Risk analysis: an overview

Risk analysis increasingly is considered as an integral part of the environmental management decision-making process. Risk, defined as the probability of occurrence of a particular adverse effect on human health or the environment, should not be confounded with hazard, defined as a source of potential injury independent of occurrence. Risk analysis has to be followed by risk management. Some opponents of risk analysis make the reproach that the science used in risk analysis is immature and consequently that the entire process in laden with hidden value judgments. Attempts to overcome these critics are increasingly based on the use of robust biologic data the final considered values system being efficacy-based, efficiency-based or equity-based. Globalization has brought with it new problems, and there is an urgent need to improve risk analysts; to increase its public acceptability and to establish consensus regarding solutions to global environmental problems. In this context biologic-based models and biomarkers hold, the greatest promise for improving risk assessment. These considerations are illustrated by a few examples, also pertaining to low-dose extrapolation and to the problem of thresholds for carcinogenesis. Future directions for development are evoked.     

Canine PKD1 is a single-copy gene: genomic organization and comparative analysis

Most cases of autosomal dominant polycystic kidney disease are caused by mutations in the gene PKD1, encoding polycystin-1. To gain insight into the role of polycystin-1 in tubulogenesis and cystogenesis using the well-characterized canine kidney epithelial cell line MDCK, we have now cloned and characterized the exon/intron structure of the canine gene PKD1. FISH analysis showed that the dog genome lacks the multiple PKD1 homologs present in human. Intron 21 of dog PKD1 lacked the polypyrimidine tract characteristic of the human gene, whereas pyrimidine-rich elements were identified in canine intron 30. Canine polycystin-1 showed a higher degree of homology with the human counterpart and lower homology with mouse and rat. A striking degree of conservation (97% identity) was determined for the leucine-rich repeat domain between dog and human. Also, the homology analysis indicated that 4 of 16 Ig-like repeats (IgIII, IgVII, IgX, and IgXV) are likely to be functionally significant. This is particularly important in light of our recent findings demonstrating that Iglike domains form strong homophilic interactions and can mediate cell-cell adhesion. These data enable detailed analysis of the role of polycystin-1 in cystogenesis and tubulogenesis using the canine MDCK cell line.     

Random sequence analysis of genomic DNA of a hyperthermophile: Aquifex pyrophilus

Aquifex pyrophilus is one of the hyperthermophilic bacteria that can grow at temperatures up to 95°C. To obtain information about its genomic structure, random sequencing was performed on plasmid libraries containing 0.5–2 kb genomic DNA fragments of A. pyrophilus. Comparison of the obtained sequence tags with known proteins revealed that 123 tags showed strong similarity to previously identifed proteins in the PIR or Genebank databases. These included three proteases, two amino acid racemases, and three enzymes utilizing oxygen as substrate. Although the GC ratio of the genome is about 40%, the codon usage of A. pyrophilus showed biased occurrence of G and C at the third position of codons, especially those for amino acids such as asparagine, aspartic acid, cysteine, glutamine, glutamic acid, histidine, lysine, and tyrosine. A higher ratio of positively charged amino acids in A. pyrophilus proteins as compared with proteins from mesophiles suggested that Aquifex proteins might contain increased ion-pair interaction that could help to maintain heat stability. Received: March 1, 1997 / Accepted: May 9, 1997     

CRAC: an integrated approach to the analysis of RNA-seq reads

A large number of RNA-sequencing studies set out to predict mutations, splice junctions or fusion RNAs. We propose a method, CRAC, that integrates genomic locations and local coverage to enable such predictions to be made directly from RNA-seq read analysis. A k-mer profiling approach detects candidate mutations, indels and splice or chimeric junctions in each single read. CRAC increases precision compared with existing tools, reaching 99:5% for splice junctions, without losing sensitivity. Importantly, CRAC predictions improve with read length. In cancer libraries, CRAC recovered 74% of validated fusion RNAs and predicted novel recurrent chimeric junctions. CRAC is available at http://crac.gforge.inria.fr.     

Mapping and sequence analysis of barley hordoindolines

Barley (Hordeum vulgare L.) cultivars vary in traits such as grain hardness and malt quality. However, little is known about the genetic basis of these grain quality traits in barley, while more is known about the basis of grain hardness in wheat (Triticum aestivum L.). Puroindolines are endosperm-specific proteins found in wheat and barley, as well as other members of the Triticeae. In wheat, variation of puroindoline sequence is associated with most of the variability in wheat grain texture. However, no information exists on sequence variation of the barley homologs of puroindolines, the hordoindolines. We have therefore chosen to isolate and characterize the hordoindoline (hin) sequences of eight North American barley cultivars. The barley sequences contain numerous non-conservative amino-acid substitutions relative to their wheat counterparts. However, no significant rearrangements were found in either hinA or hinB of barley. Three hinA and two hinB sequence types were found among the eight barley cultivars examined, indicating substantial allelic variation at this locus. The hinB sequence variability was used to map hinB to the short arm of chromosome 5H in a Steptoe/Morex mapping population, which is coincident with the previously mapped location of hinA and Gsp (grain-softness protein). This chromosomal location also coincides with a small barley malt-extract QTL, suggesting that hordoindoline sequence variation may play a small role in barley grain quality. Efforts to correlate barley seed textural differences and malting quality with hordoindoline sequence type are ongoing. Received: 25 May 2000 / Accepted: 21 September 2000     

Sequence and comparative genomic analysis of actin-related proteins

Actin-related proteins (ARPs) are key players in cytoskeleton activities and nuclear functions. Two complexes, ARP2/3 and ARP1/11, also known as dynactin, are implicated in actin dynamics and in microtubule-based trafficking, respectively. ARP4 to ARP9 are components of many chromatin-modulating complexes. Conventional actins and ARPs codefine a large family of homologous proteins, the actin superfamily, with a tertiary structure known as the actin fold. Because ARPs and actin share high sequence conservation, clear family definition requires distinct features to easily and systematically identify each subfamily. In this study we performed an in depth sequence and comparative genomic analysis of ARP subfamilies. A high-quality multiple alignment of approximately 700 complete protein sequences homologous to actin, including 148 ARP sequences, allowed us to extend the ARP classification to new organisms. Sequence alignments revealed conserved residues, motifs, and inserted sequence signatures to define each ARP subfamily. These discriminative characteristics allowed us to develop ARPAnno (http://bips.u-strasbg.fr/ARPAnno), a new web server dedicated to the annotation of ARP sequences. Analyses of sequence conservation among actins and ARPs highlight part of the actin fold and suggest interactions between ARPs and actin-binding proteins. Finally, analysis of ARP distribution across eukaryotic phyla emphasizes the central importance of nuclear ARPs, particularly the multifunctional ARP4.