Update of AMmtDB: a database of multi-aligned Metazoa mitochondrial DNA sequences

The AMmtDB database (http://bighost.area.ba.cnr.it/mitochondriome) has been updated by collecting the multi-aligned sequences of Chordata and Invertebrata mitochondrial genes coding for proteins and tRNAs. Links to the multi-aligned mtDNA intraspecies variants, collected in VarMmtDB at the Mitochondriome web site, have been introduced. The genes coding for proteins are multi-aligned based on the translated sequences and both the nucleotide and amino acid multi-alignments are provided. AMmtDB data selected through SRS can be viewed and managed using GeneDoc or other programs for the management of multi-aligned data depending on the user’s operative system. The multiple alignments have been produced with CLUSTALW and PILEUP programs and then carefully optimized manually.     

Update of AMmtDB: a database of multi-aligned metazoa mitochondrial DNA sequences.

The present paper describes AMmtDB, a database collecting the multi-aligned sequences of vertebrate mitochondrial genes coding for proteins and tRNAs, as well as the multiple alignment of the mammalian mtDNA main regulatory region (D-loop) sequences. The genes coding for proteins are multi-aligned based on the translated sequences and both the nucleotide and amino acid multi-alignments are provided. As far as the genes coding for tRNAs are concerned, the multi-alignments based on the primary and the secondary structures are both provided; for the mammalian D-loop multi-alignments we report the conserved regions of the entire D-loop (CSB1, CSB2, CSB3, the central region, ETAS1 and ETAS2) as defined by Sbisà et al. [ Gene (1997), 205, 125-140). A flatfile format for AMmtDB has been designed allowing its implementation in SRS (http://bio-www.ba.cnr.it:8000/BioWWW/#AMMTDB ). Data selected through SRS can be managed using GeneDoc or other programs for the management of multi-aligned data depending on the user's operative system. The multiple alignments have been produced with CLUSTALV and PILEUP programs and then carefully optimized manually.     

The design of synthetic genes

Computer programs are described that aid in the design of synthetic genes coding for proteins that are targets of a research program in site directed mutagenesis. These programs "reverse-translate" protein sequences into general nucleic acid sequences (those where codons have not yet been selected), map restriction sites into general DNA sequences, identify points in the synthetic gene where unique restriction sites can be introduced, and assist in the design of genes coding for hybrids and evolutionary intermediates between homologous proteins. Application of these programs therefore facilitates the use of modular mutagenesis to create variants of proteins, and the implementation of evolutionary guidance as a strategy for selecting mutants.     

Intrinsic and extrinsic approaches for detecting genes in a bacterial genome.

The unannotated regions of the Escherichia coli genome DNA sequence from the EcoSeq6 database, totaling 1,278 'intergenic' sequences of the combined length of 359,279 basepairs, were analyzed using computer-assisted methods with the aim of identifying putative unknown genes. The proposed strategy for finding new genes includes two key elements: i) prediction of expressed open reading frames (ORFs) using the GeneMark method based on Markov chain models for coding and non-coding regions of Escherichia coli DNA, and ii) search for protein sequence similarities using programs based on the BLAST algorithm and programs for motif identification. A total of 354 putative expressed ORFs were predicted by GeneMark. Using the BLASTX and TBLASTN programs, it was shown that 208 ORFs located in the unannotated regions of the E. coli chromosome are significantly similar to other protein sequences. Identification of 182 ORFs as probable genes was supported by GeneMark and BLAST, comprising 51.4% of the GeneMark 'hits' and 87.5% of the BLAST 'hits'. 73 putative new genes, comprising 20.6% of the GeneMark predictions, belong to ancient conserved protein families that include both eubacterial and eukaryotic members. This value is close to the overall proportion of highly conserved sequences among eubacterial proteins, indicating that the majority of the putative expressed ORFs that are predicted by GeneMark, but have no significant BLAST hits, nevertheless are likely to be real genes. The majority of the putative genes identified by BLAST search have been described since the release of the EcoSeq6 database, but about 70 genes have not been detected so far. Among these new identifications are genes encoding proteins with a variety of predicted functions including dehydrogenases, kinases, several other metabolic enzymes, ATPases, rRNA methyltransferases, membrane proteins, and different types of regulatory proteins.     

Comparative sequence analysis of the complete human sarcomeric myosin heavy chain family: implications for functional diversity.

The conventional myosin motor proteins that drive mammalian skeletal and cardiac muscle contraction include eight sarcomeric myosin heavy chain (MyHC) isoforms. Six skeletal MyHCs are encoded by genes found in tightly linked clusters on human and mouse chromosomes 17 and 11, respectively. The full coding regions of only two out of six mammalian skeletal MyHCs had been sequenced prior to this work. In an effort to assess the extent of sequence diversity within the human MyHC family we present new full-length coding sequences corresponding to four additional human genes: MyHC-IIb, MyHC-extraocular, MyHC-IIa and MyHC-IIx/d. This represents the first opportunity to compare the full coding sequences of all eight sarcomeric MyHC isoforms within a vertebrate organism. Sequence variability has been analyzed in the context of available structure/function data with an emphasis on potential functional diversity within the family. Results indicate that functional diversity among MyHCs is likely to be accomplished by having small pockets of sequence diversity in an otherwise highly conserved molecule.     

An Integrated Sequence-Structure Database incorporating matching mRNA sequence, amino acid sequence and protein three-dimensional structure data.

We have constructed a non-homologous database, termed the Integrated Sequence-Structure Database (ISSD) which comprises the coding sequences of genes, amino acid sequences of the corresponding proteins, their secondary structure and straight phi,psi angles assignments, and polypeptide backbone coordinates. Each protein entry in the database holds the alignment of nucleotide sequence, amino acid sequence and the PDB three-dimensional structure data. The nucleotide and amino acid sequences for each entry are selected on the basis of exact matches of the source organism and cell environment. The current version 1.0 of ISSD is available on the WWW at http://www.protein.bio.msu.su/issd/ and includes 107 non-homologous mammalian proteins, of which 80 are human proteins. The database has been used by us for the analysis of synonymous codon usage patterns in mRNA sequences showing their correlation with the three-dimensional structure features in the encoded proteins. Possible ISSD applications include optimisation of protein expression, improvement of the protein structure prediction accuracy, and analysis of evolutionary aspects of the nucleotide sequence-protein structure relationship.     

Update of MmtDB: a Metazoa mitochondrial DNA variants database.

The present paper describes the improvements in MmtDB, a specialised database designed to collect Metazoa mitochondrial DNA variants. Priority in the data collection has been given to Metazoa for which a large amount of variants is available, e.g., for humans. Starting from the sequences available in the Nucleotide Sequence Databases, the redundant sequences have been removed and new sequences from other sources have been added. Value-added information is associated to each variant sequence, e.g., analysed region, experimental method, tissue and cell lines, population data, sex, age, family code and information about the variation events (nucleotide position, involved gene, restriction site gain or loss). Cross-references are introduced to the EMBL Data Library, as well as an internal cross-referencing among MmtDB entries according to tissual, heteroplasmic, familiar and aplotypical correlation. Furthermore MmtDB has a new section, AMmtDB: Aligned Metazoan mitochondrial biosequences. MmtDB can be accessed through the World Wide Web at URL http://WWW.ba.cnr.it/[symbol: see text]areamt08/MmtDBWWW.htm     

The compositional transition of vertebrate genomes: an analysis of the secondary structure of the proteins encoded by human genes

Fluctuations and increments of both C(3) and G(3) levels along the human coding sequences were investigated comparing two sets of Xenopus/human orthologous genes. The first set of genes shows minor differences of the GC(3) levels, the second shows considerable increments of the GC(3) levels in the human genes. In both data sets, the fluctuations of C(3) and G(3) levels along the coding sequences correlated with the secondary structures of the encoded proteins. The human genes that underwent the compositional transition showed a different increment of the C(3) and G(3) levels within and among the structural units of the proteins. The relative synonymous codon usage (RSCU) of several amino acids were also affected during the compositional transition, showing that there exists a correlation between RSCU and protein secondary structures in human genes. The importance of natural selection for the formation of isochore organization of the human genome has been discussed on the basis of these results.     

Analysis of membrane protein genes in a Brazilian isolate of Anaplasma marginale

The sequencing of the complete genome of Anaplasma marginale has enabled the identification of several genes that encode membrane proteins, thereby increasing the chances of identifying candidate immunogens. Little is known regarding the genetic variability of genes that encode membrane proteins in A. marginale isolates. The aim of the present study was to determine the degree of conservation of the predicted amino acid sequences of OMP1, OMP4, OMP5, OMP7, OMP8, OMP10, OMP14, OMP15, SODb, OPAG1, OPAG3, VirB3, VirB9-1, PepA, EF-Tu and AM854 proteins in a Brazilian isolate of A. marginale compared to other isolates. Hence, primers were used to amplify these genes: omp1, omp4, omp5, omp7, omp8, omp10, omp14, omp15, sodb, opag1, opag3, virb3, VirB9-1, pepA, ef-tu and am854. After polimerase chain reaction amplification, the products were cloned and sequenced using the Sanger method and the predicted amino acid sequence were multi-aligned using the CLUSTALW and MEGA 4 programs, comparing the predicted sequences between the Brazilian, Saint Maries, Florida and A. marginale centrale isolates. With the exception of outer membrane protein (OMP) 7, all proteins exhibited 92-100% homology to the other A. marginale isolates. However, only OMP1, OMP5, EF-Tu, VirB3, SODb and VirB9-1 were selected as potential immunogens capable of promoting cross-protection between isolates due to the high degree of homology (over 72%) also found with A. (centrale) marginale.     

The relationship between base composition and codon usage in bacterial genes and its use for the simple and reliable identification of protein-coding sequences

Bacterial genes that code for proteins appear to possess a codon usage characteristic of their overall base composition. This results in different but predictable non-random distributions of nucleotides within codons, permitting the recognition of protein-coding sequences in a wide range of bacterial species. The nature of this distribution depends on the base composition of the coding sequence. The position-specific differences are especially conspicuous in genes of extreme G + C content, allowing the particularly reliable prediction of the reading frame and coding strand of experimentally determined DNA sequences. This fmding has been exploited to identify the coding sequence of the viomycin phosphotransferase (vph) gene of Streptomyces vinaceus. An easily applied computer program (“Frame”) has been written to carry out and display such analyses.     

Structural analysis of Arabidopsis thaliana chromosome 5. X. Sequence features of the regions of 3,076,755 bp covered by sixty P1 and TAC clones.

In our ongoing project to deduce the nucleotide sequence of Arabidopsis thaliana chromosome 5, non-redundant P1 and TAC clones have been sequenced on the basis of the fine physical map, and as of January, 2000, the sequences of 16.6 Mb representing approximately 60% of chromosome 5 have been accumulated and released at our web site. Along with the sequence determination, structural features of the sequenced regions have been analyzed by applying a variety of computer programs, and we already predicted a total of 2697 potential protein coding genes in the 11,166,130 bp regions, which are covered by 159 P1 and TAC clones. In this paper, we describe the structural features of the 3,076,755 bp regions covered by newly analyzed 60 P1 and TAC clones. A total of 715 potential protein coding genes were identified, giving an average density of the genes identified of 1 gene per 4001 bp. Introns were observed in 80% of the genes, and the average number per gene and the average length of the introns were 4.5 and 147 bp, respectively. These sequence features are nearly identical to those in our latest report in which the data were compiled based on a new standard of gene assignment including the computer-predicted hypothetical genes. The regions also contained 12 tRNA genes when searched by similarity to reported tRNA genes and the tRNA scan-SE program. The sequence data and information on the potential genes are available through the World Wide Web database KAOS (Kazusa Arabidopsis data Opening Site) at http://www.kazusa.or.jp/kaos/.     

Mytilus edulis Histone Gene Clusters Containing Only H1 Genes

We isolated five different phage clones containing histone gene clusters with up to five H1 genes per phage clone from a Mytilus edulis genomic library. Among these H1 genes, nine gene types coding for five different H1 proteins have been identified. All H1 histone genes were located on repetitive restriction fragments with only slightly different sizes. The H1 coding regions show highly related sequences, suggesting that the multitude of H1 genes has evolved by gene duplication events. Core histone genes could not be found on these five Mytilus edulis genome fragments. Received: 28 July 1998 / Accepted: 17 May 1999     

Statistical evaluation of the coding capacity of complementary DNA strands

Two independent methods are used to evaluate the protein-coding information content in different classes of DNA sequences. The first method allows to evaluate the statistical relevance of finding unidentified reading frames, longer than 100 codons, on both DNA strands of: a) 117 DNA sequences that code for 142 nuclear proteins; b) 39 stable RNA coding sequences and c) 36 other DNA sequences which include regulatory and as yet unknown function sequences. The finding of 50 reading frames longer than 100 codons (complementary inverted proteins or c.i.p. genes) located on the DNA strand complementary to the protein-coding one is drastically in excess of the number predicted by chance alone. An independent method (testcode) applied to c.i.p. gene sequences, which assigns the probability of coding to a given sequence, predicts that more than 50% of these genes are translated in a functional product. These analyses indicate the existence of a new class of protein-coding genes, located on the DNA sequences complementary to the protein-coding DNA strand.     

Temporal analysis of French Bordetella pertussis isolates by comparative whole-genome hybridization

Bordetella pertussis, a gram-negative beta-proteobacterium, is the agent of whooping cough in humans. Whooping cough remains a public health problem worldwide, despite well-implemented infant/child vaccination programs. It continues to be endemic and is observed cyclically in vaccinated populations. Classical molecular subtyping methods indicate that genome diversity among B. pertussis isolates is limited. Although the whole bacterial genome has been studied by pulsed-field gel electrophoresis, the genes implicated in the diversity have not been identified. We developed a B. pertussis whole-genome DNA microarray representing over 91% of the predicted coding sequences of the sequenced strain Tohama I. Genomic DNA from clinical isolates with various pulsed-field gel electrophoresis profile patterns was competitively hybridized with the DNA microarray and coding sequences were classified as present, absent or duplicated. Our data strongly suggest that the B. pertussis population is dynamic. In France, with highly vaccinated population, the genetic diversity is low and decreasing with time, and clonal expansion correlates with cycles of the disease. This decrease in diversity is essentially due to loss of genes and pseudogenes. The genes deleted are most of the time flanked by insertion sequences.