In silico pattern-based analysis of the human cytomegalovirus genome

More than 200 open reading frames (ORFs) from the human cytomegalovirus genome have been reported as potentially coding for proteins. We have used two pattern-based in silico approaches to analyze this set of putative viral genes. With the help of an objective annotation method that is based on the Bio-Dictionary, a comprehensive collection of amino acid patterns that describes the currently known natural sequence space of proteins, we have reannotated all of the previously reported putative genes of the human cytomegalovirus. Also, with the help of MUSCA, a pattern-based multiple sequence alignment algorithm, we have reexamined the original human cytomegalovirus gene family definitions. Our analysis of the genome shows that many of the coded proteins comprise amino acid combinations that are unique to either the human cytomegalovirus or the larger group of herpesviruses. We have confirmed that a surprisingly large portion of the analyzed ORFs encode membrane proteins, and we have discovered a significant number of previously uncharacterized proteins that are predicted to be G-protein-coupled receptor homologues. The analysis also indicates that many of the encoded proteins undergo posttranslational modifications such as hydroxylation, phosphorylation, and glycosylation. ORFs encoding proteins with similar functional behavior appear in neighboring regions of the human cytomegalovirus genome. All of the results of the present study can be found and interactively explored online (http://cbcsrv.watson.ibm.com/virus/).     

Characterization of Trypanosoma (Duttonella) vivax by isoenzyme analysis

Stocks derived from 10 different primary isolates of T. vivax were subjected to isoenzyme analysis for 34 enzymes by both isoelectric focusing in agarose and electrophoresis in starch gel. Trypanosomes were measured and their morphology examined for comparison with the biochemical data. Thirteen enzymes (14 zymograms) were selected to construct isoenzyme profiles. Nine different zymodemes were identified and only two stocks were identical. Both rodent infectivity and the production of the haemorrhagic syndrome could be correlated with the isoenzyme profiles.     

The fatty acids of Trypanosoma vivax

A study was carried out to determine the lipid composition of the blood-stream form of the African trypanosome. Trypanosoma vivax. Data from thin layer chromatography showed that the major polar lipids were lysophosphatidylcholine, sphingomyelin, phosphatidylcholine, phosphatidylethanolamine, phosphatidylserine and diphosphatidylglycerol. The major neutral lipids were sterol, monoacylglycerol, diacylglycerol, free fatty acid and triacyglycerol. 16:0, 18:0, 18:1 and 18:2 constituted the major fatty acids of both the polar and neutral lipid fractions. The work constituted the first detailed study on the fatty acid composition of this African trypanosome.     

The phosphoglucose isomerases of the bloodstream forms of Trypanosoma brucei and Trypanosoma vivax.

1. The phosphoglucose isomerases (PGI's) of the bloodstream forms of Trypanosoma brucei and T. vivax have been purified some 150-fold, using cellulose ion-exchange chromatography, gel filtration and isoelectric focussing. 2. The two trypanosome enzymes showed many similarities in kinetic properties, but differed from each other somewhat in thermal stability and in isoelectric point. 3. Both trypanosome enzymes differ from PGI's from other sources in having a higher Ki for the competitive inhibitor 6-phosphogluconate.     

Exploring copy number variation in the rabbit (Oryctolagus cuniculus) genome by array comparative genome hybridization

The European rabbit (Oryctolagus cuniculus) is relevant in a large spectrum of fields: it is a livestock, a pet, a biomedical model and a biotechnology tool, a wild resource and a pest. The sequencing of the rabbit genome has opened new perspectives to study this lagomorph at the genome level. We herein investigated for the first time the O. cuniculus genome by array comparative genome hybridization (aCGH) and established a first copy number variation (CNV) genome map in this species comprising 155 copy number variation regions (CNVRs; 95 gains, 59 losses, 1 with both gain and loss) covering ~0.3% of the OryCun2.0 version. About 50% of the 155 CNVRs identified spanned 139 different protein coding genes, 110 genes of which were annotated or partially annotated (including Major Histocompatibility Complex genes) with 277 different gene ontology terms. Many rabbit CNVRs might have a functional relevance that should be further investigated.     

Insights from the GC content analysis of 76genome survey sequences (GSS) from Elaeisoleifera

South American oil-palm (Elaeis oleifera) is not cultivated in tropical countries like Malaysia on large scale due to low yield of palm oil derived from its fruit mesocarp. However, its fruit mesocarp oil contains about 68.6 % oleic acid (C_18:1) which is more than double in comparison to commercially cultivated oilpalm, E. guineensis Jacq Tenera (hybrid of Dura (♀) x Pisifera (♂)). It is also known that E. oleifera is a good source of tocotrienols and carotenoids. Therefore, it is of interest to know the genome sequence of E. oleifera. The objective of this study is to generate genome survey sequences (GSS) to get GC content insight in the E. oleifera genome. The nuclear genomic DNA isolated from young leaf‐tissues was digested with EcoRI and NdeI/DraI restriction enzymes; and three genomic DNA libraries were constructed using Lambda ZAP‐II, pGEM®‐T Easy, and pDONR 222™ as cloning vectors. Generated 76 GSSs were analyzed by using Bioinformatics tools. The analysis result indicates that the adenine, cytosine, guanine and thymine content in generated GSSs are 30%, 20%, 20%, and 30% respectively. In conclusion, based on the precise GC content analysis of the randomly isolated 76 GSSs by using Bioinformatics tools we hypothesize that GC content in E. oleifera genome is 40%. The hypothesized 40% GC content in E. oleifera genome is expected to remain close to the GC content based on the whole genome analysis.^ψThe nucleotide sequence data reported in this paper have been submitted to dbGSS division of the international DNA database (GenBank/DDBJ/EMBL) under accession numbers: {"type":"entrez-nucleotide-range","attrs":{"text":"DX575945- DX575972","start_term":"DX575945","end_term":"DX575972","start_term_id":"94323534","end_term_id":"94323561"}}DX575945- DX575972 and {"type":"entrez-nucleotide-range","attrs":{"text":"EI798032-EI798079","start_term":"EI798032","end_term":"EI798079","start_term_id":"146199019","end_term_id":"146199066"}}EI798032-EI798079.

Abbreviations

gDNA - Nuclear genomic DNA, GSSs - Genome survey sequences K12, SAOP - South American oil‐palm Db1     

Studies on the mechanism of adenosine transport in Trypanosoma vivax

Adenosine uptake in the presence of some metabolic inhibitors and nucleosides has been studied. The uptake of adenosine was inhibited by oubain, phlorizin, iodoacetate and coformycin. Guanosine, on the other hand stimulated adenosine uptake to a considerable extent. Neither thymidine nor inosine caused significant change in adenosine uptake. Results of the time course assay and uptake studies at various concentrations of adenosine suggest that possibly more than one mode of uptake operates in the transport of adenosine in T. Vivax.     

In silico analysis of the restriction fragments length distribution in the human genome.]

The Restriction On Computer (ROC) program (freely available at http://www.mcb.harvard.edu/gilbert/ROC) was developed and used to analyze the restriction fragment length distribution in the human genome. In contrast to other programs searching for restriction sites, ROC simultaneously analyzes several long nucleotide sequences, such as the entire genomes, and in essence simulates electrophoretic analysis of DNA restriction fragments. In addition, this program extracts and analyzes DNA repeats that account for peaks in the restriction fragment length distribution. The ROC analysis data are consistent with the experimental data obtained via in vitro restriction enzyme analysis (taxonomic printing). A difference between the in vitro and in silico results is explained by underrepresentation of tandem DNA repeats in genomic databases. The ROC analysis of individual genome fragments elucidated the nature of several DNA markers, which were earlier revealed by taxonomic printing, and showed that L1 and Alu repeats are nonrandomly distributed in various chromosomes. Another advantage is that the ROC procedure makes it possible to analyze the nonrandom character of a genomic distribution of short DNA sequences. The ROC analysis showed that a low poly(G) frequency is characteristic of the entire human genome, rather than of only coding sequences. The method was proposed for a more complex in silico analysis of the genome. For instance, it is possible to simulate DNA restriction together with blot hybridization and then to analyze the nature of markers revealed.     

Genome-wide comparative in silico analysis of calcium transporters of rice and sorghum

The mechanism of calcium uptake, translocation and accumulation in Poaceae has not yet been fully understood. To address this issue, we conducted genome-wide comparative in silico analysis of the calcium (Ca(2+)) transporter gene family of two crop species, rice and sorghum. Gene annotation, identification of upstream cis-acting elements, phylogenetic tree construction and syntenic mapping of the gene family were performed using several bioinformatics tools. A total of 31 Ca(2+) transporters, distributed on 9 out of 12 chromosomes, were predicted from rice genome, while 28 Ca(2+) transporters predicted from sorghum are distributed on all the chromosomes except chromosome 10 (Chr 10). Interestingly, most of the genes on Chr 1 and Chr 3 show an inverse syntenic relationship between rice and sorghum. Multiple sequence alignment and motif analysis of these transporter proteins revealed high conservation between the two species. Phylogenetic tree could very well identify the subclasses of channels, ATPases and exchangers among the gene family. The in silico cis-regulatory element analysis suggested diverse functions associated with light, stress and hormone responsiveness as well as endosperm- and meristem-specific gene expression. Further experiments are warranted to validate the in silico analysis of the predicted transporter gene family and elucidate the functions of Ca(2+) transporters in various biological processes.     

Dataset of potential targets for Mycobacterium tuberculosis H37Rv through comparative genome analysis

Mycobacterium tuberculosis is the causative agent of the disease, tuberculosis and H37Rv is the most studied clinical strain. We use comparative genome analysis of Mycobacterium tuberculosis H37Rv and human for the identification of potential targets dataset. We used DEG (Database of Essential Genes) to identify essential genes in the H37Rv strain. The analysis shows that 628 of the 3989 genes in Mycobacterium tuberculosis H37Rv were found to be essential of which 324 genes lack similarity to the human genome. Subsequently hypothetical proteins were removed through manual curation. This further resulted in a dataset of 135 proteins with essential function and no homology to human.     

Enzyme-substrate interactions in the purine-specific nucleoside hydrolase from Trypanosoma vivax

Nucleoside hydrolases are key enzymes in the purine salvage pathway of Trypanosomatidae and are considered as targets for drug design. We previously reported the first x-ray structure of an inosine-adenosine-guanosine preferring nucleoside hydrolase (IAG-NH) from Trypanosoma vivax (). Here we report the 2.0-A crystal structure of the slow D10A mutant in complex with the inhibitor 3-deaza-adenosine and the 1.6-A crystal structure of the same enzyme in complex with a genuine substrate inosine. The enzyme-substrate complex shows the substrate bound to the enzyme in a different conformation from 3-deaza-adenosine and provides a snapshot along the reaction coordinate of the enzyme-catalyzed reaction. The chemical groups on the substrate important for binding and catalysis are mapped. The 2'-OH, 3'-OH, and 5'-OH contribute 4.6, 7.5, and 5.4 kcal/mol to k(cat)/K(m), respectively. Specific interactions with the exocyclic groups on the purine ring are not required for catalysis. Site-directed mutagenesis indicates that the purine specificity of the IAG-NHs is imposed by a parallel aromatic stacking interaction involving Trp(83) and Trp(260). The pH profiles of k(cat) and k(cat)/K(m) indicate the existence of one or more proton donors, possibly involved in leaving group activation. However, mutagenesis of the active site residues around the nucleoside base and an alanine scan of a flexible loop near the active site fail to identify this general acid. The parallel aromatic stacking seems to provide the most likely alternative mechanism for leaving group activation.     

Sequence and comparative analysis of the maize NB mitochondrial genome

The NB mitochondrial genome found in most fertile varieties of commercial maize (Zea mays subsp. mays) was sequenced. The 569,630-bp genome maps as a circle containing 58 identified genes encoding 33 known proteins, 3 ribosomal RNAs, and 21 tRNAs that recognize 14 amino acids. Among the 22 group II introns identified, 7 are trans-spliced. There are 121 open reading frames (ORFs) of at least 300 bp, only 3 of which exist in the mitochondrial genome of rice (Oryza sativa). In total, the identified mitochondrial genes, pseudogenes, ORFs, and cis-spliced introns extend over 127,555 bp (22.39%) of the genome. Integrated plastid DNA accounts for an additional 25,281 bp (4.44%) of the mitochondrial DNA, and phylogenetic analyses raise the possibility that copy correction with DNA from the plastid is an ongoing process. Although the genome contains six pairs of large repeats that cover 17.35% of the genome, small repeats (20-500 bp) account for only 5.59%, and transposable element sequences are extremely rare. MultiPip alignments show that maize mitochondrial DNA has little sequence similarity with other plant mitochondrial genomes, including that of rice, outside of the known functional genes. After eliminating genes, introns, ORFs, and plastid-derived DNA, nearly three-fourths of the maize NB mitochondrial genome is still of unknown origin and function.     

A comparative analysis of serpin genes in the silkworm genome

Serine protease inhibitors (serpins) are a superfamily of proteins, most of which control protease-mediated processes by inhibiting their cognate enzymes. Sequencing of the silkworm genome provides an opportunity to investigate serpin structure, function, and evolution at the genome level. There are thirty-four serpin genes in Bombyx mori. Six are highly similar to their Manduca sexta orthologs that regulate innate immunity. Three alternative exons in serpin1 gene and four in serpin28 encode a variable region including the reactive site loop. Splicing of serpin2 pre-mRNA yields variations in serpin2A, 2A′ and 2B. Sequence similarity and intron positions reveal the evolutionary pathway of seven serpin genes in group C. RT-PCR indicates an increase in the mRNA levels of serpin1, 3, 5, 6, 9, 12, 13, 25, 27, 32 and 34 in fat body and hemocytes of larvae injected with bacteria. These results suggest that the silkworm serpins play regulatory roles in defense responses.