The EMBL data library

Back to the Colorectal Cancer Consensus Molecular Subtype Future

PURPOSE OF REVIEW

This review seeks to provide an informed prospective on the advances in molecular profiling and analysis of colorectal cancer (CRC). The goal is to provide a historical context and current summary on how advances in gene and protein sequencing technology along with computer capabilities led to our current bioinformatic advances in the field.

RECENT FINDINGS

An explosion of knowledge has occurred regarding genetic, epigenetic, and biochemical alterations associated with the evolution of colorectal cancer. This has led to the realization that CRC is a heterogeneous disease with molecular alterations often dictating natural history, response to treatment, and outcome. The consensus molecular subtypes (CMS) classification classifies CRC into four molecular subtypes with distinct biological characteristics, which may form the basis for clinical stratification and subtype-based targeted intervention. This review summarizes new developments of a field moving "Back to the Future." CRC molecular subtyping will better identify key subtype specific therapeutic targets and responses to therapy.

A brief history of bioinformatics

It is easy for today’s students and researchers to believe that modern bioinformatics emerged recently to assist next-generation sequencing data analysis. However, the very beginnings of bioinformatics occurred more than 50 years ago, when desktop computers were still a hypothesis and DNA could not yet be sequenced. The foundations of bioinformatics were laid in the early 1960s with the application of computational methods to protein sequence analysis (notably, de novo sequence assembly, biological sequence databases and substitution models). Later on, DNA analysis also emerged due to parallel advances in (i) molecular biology methods, which allowed easier manipulation of DNA, as well as its sequencing, and (ii) computer science, which saw the rise of increasingly miniaturized and more powerful computers, as well as novel software better suited to handle bioinformatics tasks. In the 1990s through the 2000s, major improvements in sequencing technology, along with reduced costs, gave rise to an exponential increase of data. The arrival of ‘Big Data’ has laid out new challenges in terms of data mining and management, calling for more expertise from computer science into the field. Coupled with an ever-increasing amount of bioinformatics tools, biological Big Data had (and continues to have) profound implications on the predictive power and reproducibility of bioinformatics results. To overcome this issue, universities are now fully integrating this discipline into the curriculum of biology students. Recent subdisciplines such as synthetic biology, systems biology and whole-cell modeling have emerged from the ever-increasing complementarity between computer science and biology.

SINA: accurate high-throughput multiple sequence alignment of ribosomal RNA genes

Motivation: In the analysis of homologous sequences, computation of multiple sequence alignments (MSAs) has become a bottleneck. This is especially troublesome for marker genes like the ribosomal RNA (rRNA) where already millions of sequences are publicly available and individual studies can easily produce hundreds of thousands of new sequences. Methods have been developed to cope with such numbers, but further improvements are needed to meet accuracy requirements.

Results: In this study, we present the SILVA Incremental Aligner (SINA) used to align the rRNA gene databases provided by the SILVA ribosomal RNA project. SINA uses a combination of k-mer searching and partial order alignment (POA) to maintain very high alignment accuracy while satisfying high throughput performance demands.

SINA was evaluated in comparison with the commonly used high throughput MSA programs PyNAST and mothur. The three BRAliBase III benchmark MSAs could be reproduced with 99.3, 97.6 and 96.1 accuracy. A larger benchmark MSA comprising 38 772 sequences could be reproduced with 98.9 and 99.3% accuracy using reference MSAs comprising 1000 and 5000 sequences. SINA was able to achieve higher accuracy than PyNAST and mothur in all performed benchmarks.

Availability: Alignment of up to 500 sequences using the latest SILVA SSU/LSU Ref datasets as reference MSA is offered at http://www.arb-silva.de/aligner. This page also links to Linux binaries, user manual and tutorial. SINA is made available under a personal use license.

Contact:epruesse@mpi-bremen.de

Supplementary information:Supplementary data are available at Bioinformatics online.

Aetiological diagnosis of brain abscesses and spinal infections: application of broad range bacterial polymerase chain reaction analysis

OBJECTIVE

To evaluate the usefulness of the broad range bacterial rDNA polymerase chain reaction (PCR) method combined with DNA sequencing in the aetiological diagnosis of intracranial or spinal infections in neurosurgical patients.

METHODS

In addition to conventional methods, the broad range bacterial PCR approach was applied to examine pus or tissue specimens from cerebral or spinal lesions in patients treated in a neurosurgical unit for a clinical or neuroradiological suspicion of bacterial brain abscess or spondylitis.

RESULTS

Among the 44 patients with intracranial or spinal lesions, the final diagnosis suggested bacterial disease in 25 patients, among whom the aetiological agent was identified in 17. A causative bacterial species was identified only by the rDNA PCR method in six cases, by both the PCR methodology and bacterial culture in six cases, and by bacterial culture alone in five. All samples in which a bacterial aetiology was identified only by the PCR approach were taken during antimicrobial treatment, and in three patients the method yielded the diagnosis even after >/= 12 days of parenteral treatment. One case also identified by the PCR approach alone involved a brain abscess caused by Mycoplasma hominis, which is not readily cultured by routine methods.

CONCLUSIONS

In patients with brain abscesses and spinal infections, the broad range bacterial rDNA PCR approach may be the only method to provide an aetiological diagnosis when the patient is receiving antimicrobial treatment, or when the causative agent is fastidious.

Database on the structure of large subunit ribosomal RNA

The Antwerp database on large subunit ribosomal RNA now contains 607 complete or nearly complete aligned sequences. The alignment incorporates secondary structure information for each sequence. Other information about the sequences, such as literature references, accession numbers and taxonomic information is also available. Information from the database can be downloaded or searched on the rRNA WWW Server at URL http://rrna.uia.ac.be/

The hrp23 protein in the balbiani ring pre-mRNP particles is released just before or at the binding of the particles to the nuclear pore complex

Balbiani ring (BR) pre-mRNP particles reside in the nuclei of salivary glands of the dipteran Chironomus tentans and carry the message for giant-sized salivary proteins. In the present study, we identify and characterize a new protein component in the BR ribonucleoprotein (RNP) particles, designated hrp23. The protein with a molecular mass of 20 kD has a single RNA-binding domain and a glycine-arginine-serine-rich auxiliary domain. As shown by immunoelectron microscopy, the hrp23 protein is added to the BR transcript concomitant with transcription, is still present in the BR particles in the nucleoplasm, but is absent from the BR particles that are bound to the nuclear pore complex or are translocating through the central channel of the complex. Thus, hrp23 is released just before or at the binding of the particles to the nuclear pore complex. It is noted that hrp23 behaves differently from two other BR RNP proteins earlier studied: hrp36 and hrp45. These proteins both reach the nuclear pore complex, and hrp36 even accompanies the RNA into the cytoplasm. It is concluded that each BR RNA-binding protein seems to have a specific flow pattern, probably related to the particular role of the protein in gene expression.

The relationship between synonymous codon usage and protein structure

The hypothesis that synonymous codon usage is related to protein three-dimensional structure is examined by investigating the correlation between synonymous codon usage and protein secondary structure. All except two codons in E. coli show the same secondary structural preference for alpha-helix, beta-strand or coil as that of amino acids to be encoded by the respective codons, while 17 codons show secondary structural bias in mammalian proteins. The results indicate that there is no significant correlation between synonymous codon usage and protein secondary structure in E. coli, but there is a correlation in mammals. It could be deduced that synonymous codons carry much less structural information in prokaryotes than in eukaryotes due to their divergent evolutionary mechanism.

Diagnosis of meningococcal meningitis by broad-range bacterial PCR with cerebrospinal fluid

We used broad-range bacterial PCR combined with DNA sequencing to examine prospectively cerebrospinal fluid (CSF) samples from patients with suspected meningitis. Fifty-six CSF samples from 46 patients were studied during the year 1995. Genes coding for bacterial 16S and/or 23S rRNA genes could be amplified from the CSF samples from five patients with a clinical picture consistent with acute bacterial meningitis. For these patients, the sequenced PCR product shared 98.3 to 100% homology with the Neisseria meningitidis sequence. For one patient, the diagnosis was initially made by PCR alone. Of the remaining 51 CSF samples, for 50 (98.0%) samples the negative PCR findings were in accordance with the negative findings by bacterial culture and Gram staining, as well as with the eventual clinical diagnosis for the patient. However, the PCR test failed to detect the bacterial rRNA gene in one CSF sample, the culture of which yielded Listeria monocytogenes. These results invite new research efforts to be focused on the application of PCR with broad-range bacterial primers to improve the etiologic diagnosis of bacterial meningitis. In a clinical setting, Gram staining and bacterial culture still remain the cornerstones of diagnosis.

Sequence and gene content in 35 kb genomic clone mapping in the human Xq27.1 region

This paper presents detailed analysis of the entire sequence of a cosmid clone, 26H7, containing 35 kb of human DNA. This cosmid resides on the q27.1 region of the human X chromosome between, DXS1232 and DXS119 loci. Novel potential small exons were detected for which conventional gene identification strategies (Northern blot analysis and extensive cDNA library screening) proved to be inefficient. Of the standard repetitive elements we found: 8 Alu's making up 6.2% of the sequence; 10 MIR segments (4.1%); 5 LINE1 elements (4.8%), 3 MIR2 (1.0%); 2 MLT (2.9%), and 1 MSTA (0.7%) representing about 20% of the total sequence. The overall GC content was rather low, only 42% and no CpG island was detected using rare restriction enzymes. However, a CpG-rich region was identified. Computer aided analysis of the sequence inferred the presence of three possible genes: one of them was found to be homologous to the U7 RNA family elements; a second is reported in this paper, however at the moment no significant homology has been found in the data bank. The third predicted gene has not as yet been found to be detectable by RT-PCR. We also report in this paper the identification of X-chromosome specific repeated sequences.

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 EMBL nucleotide sequence database

The EMBL Nucleotide Sequence Database (http://www.ebi.ac.uk/embl. html ) constitutes Europe's primary nucleotide sequence resource. DNA and RNA sequences are directly submitted from researchers and genome sequencing groups and collected from the scientific literature and patent applications (Fig. 1). In collaboration with DDBJ and GenBank the database is produced, maintained and distributed at the European Bioinformatics Institute. Database releases are produced quarterly and are distributed on CD-ROM. EBI's network services allow access to the most up-to-date data collection via Internet and World Wide Web interface, providing database searching and sequence similarity facilities plus access to a large number of additional databases.

Analysis of donor splice sites in different eukaryotic organisms

We present here a new algorithm for functional site analysis. It is based on four main assumptions: each variation of nucleotide composition makes a different contribution to the overall binding free energy of interaction between a functional site and another molecule; nonfunctioning site-like regions (pseudosites) are absent or rare in genomes; there may be errors in the sample of sites; and nucleotides of different site positions are considered to be mutually dependent. In this algorithm, the site set is divided into subsets, each described by a certain consensus. Donor splice sites of the human protein-coding genes were analyzed. Comparing the results with other methods of donor splice site prediction has demonstrated a more accurate prediction of consensus sequences AG/GU(A,G), G/GUnAG, /GU(A,G)AG, /GU(A,G)nGU, and G/GUA than is achieved by weight matrix and consensus (A,C)AG/GU(A,G)AGU with mismatches. The probability of the first type error, E1, for the obtained consensus set was about 0.05, and the probability of the second type error, E2, was 0.15. The analysis demonstrated that accuracy of the functional site prediction could be improved if one takes into account correlations between the site positions. The accuracy of prediction by using human consensus sequences was tested on sequences from different organisms. Some differences in consensus sequences for the plant Arabidopsis sp., the invertebrate Caenorhabditis sp., and the fungus Aspergillus sp. were revealed. For the yeast Saccharomyces sp. only one conservative consensus, /GUA(U,A,C)G(U,A,C), was revealed (E1 = 0.03, E2 = 0.03). Yeast is a very interesting model to use for analysis of molecular mechanisms of splicing.

Identification and quantification of Bifidobacterium species isolated from food with genus-specific 16S rRNA-targeted probes by colony hybridization and PCR

A Bifidobacterium genus-specific target sequence in the V9 variable region of the 16S rRNA has been elaborated and was used to develop a hybridization probe. The specificity of this probe, named lm3 (5'-CGGGTGCTI*CCCACTTTCATG-3'), was used to identify all known type strains and distinguish them from other bacteria. All of the 30 type strains of Bifidobacterium which are available at the German culture collection Deutsche Sammlung von Mikroorganismen und Zellkulturen, 6 commercially available production strains, and 34 closely related relevant strains (as negative controls) were tested. All tested bifidobacteria showed distinct positive signals by colony hybridization, whereas all negative controls showed no distinct dots except Gardnerella vaginalis DSM4944 and Propionibacterium freudenreichii subsp. shermanii DSM4902, which gave slight signals. Furthermore, we established a method for isolation and identification of bifidobacteria from food by using a PCR assay without prior isolation of DNA but breaking the cells with proteinase K. By this method, all Bifidobacterium strains lead to a DNA product of the expected size. We also established a quick assay to quantitatively measure Bifidobacterium counts in food and feces by dilution plating and colony hybridization. We were able to demonstrate that 2.1 x 10(6) to 2.3 x 10(7) colonies/g of sour milk containing bifidobacteria hybridized with the specific nucleotide probe. With these two methods, genus-specific colony hybridization and genus-specific PCR, it is now possible to readily and accurately detect any bifidobacteria in food and fecal samples and to discriminate between them and members of other genera.

Seventy-five percent accuracy in protein secondary structure prediction

In this study we present an accurate secondary structure prediction procedure by using an query and related sequences. The most novel aspect of our approach is its reliance on local pairwise alignment of the sequence to be predicted with each related sequence rather than utilization of a multiple alignment. The residue-by-residue accuracy of the method is 75% in three structural states after jack-knife tests. The gain in prediction accuracy compared with the existing techniques, which are at best 72%, is achieved by secondary structure propensities based on both local and long-range effects, utilization of similar sequence information in the form of carefully selected pairwise alignment fragments, and reliance on a large collection of known protein primary structures. The method is especially appropriate for large-scale sequence analysis of efforts such as genome characterization, where precise and significant multiple sequence alignments are not available or achievable.

Analysis of the distribution of binding sites for a tissue-specific transcription factor in the vertebrate genome

Hepatocyte nuclear factor 1 (HNF1) is a dimeric homeoprotein expressed in hepatocytes and in a few other epithelial cells where it helps regulate the expression of a specific subset of genes. In an attempt to identify novel target genes for HNF1 and to assess the distribution of its target sites within the vertebrate genome, we performed a computer-assisted search within the available databases using a weighted matrix. Several hundred potential target sequences were identified within the GenBank and EMBL data banks. DNA binding assays demonstrated that more than 95%, of the new sites tested (52 sites among 54) bound HNF1. Surprisingly many HNF1 target sites were found in genes that are transcribed in cell types that do not contain the protein. On the other hand these sites are 2.5 to five times more frequent in hepatic genes than expected. It seems that the presence of HNF1 sites in liver-specific genes was favoured, but that no counter-selection occurred within the rest of the genome. HNF1 binding sites in liver genes are more often associated in clusters with sites for other transcription factors and the enrichment is more pronounced in promoter regions. We identified more than 100 liver specific genes that are potentially regulated by HNF1.

The translational signal database, TransTerm: more organisms, complete genomes

TransTerm is a database of initiation and termination sequence contexts from more than 250 organisms listed in GenBank, including the four complete genomes:Haemophilus influenzae, Methanococcus jannaschii, Mycoplasma genitalium,and Saccharomyces cerevisiae. For the current release, more than 60 000 coding sequences were analysed. The tabulated data include initiation and termination contexts organised by species along with quantitative parameters about individual coding sequences (length, %GC, GC3, Nc and CAI). There are also tables of initiation- and termination-region nucleotide-frequencies, codon usage tables and summaries of stop signal usage. TransTerm is available on the World Wide Web at: http://biochem.otago.ac.nz:800/Transterm/homepage.h tml

The future of protein secondary structure prediction accuracy

BACKGROUND

The accuracy of secondary structure prediction for a protein from knowledge of its sequence has been significantly improved by about 7% to the 70-75% range by inclusion of information residing in sequences similar to the query sequence. The scientific literature has been inconsistent, if not negative, regarding chances for further improvement from the vast knowledge to be provided by genome sequencing efforts.

RESULTS

By applying a prediction technique that is particularly sensitive to added sequence information to a standard set of query sequences with related primary structures taken from chronologically successive releases of the SWISS-PROT database, it is shown that prediction accuracy can be expected to reach 80-85% with a large 10-fold increase in present sequence knowledge.

CONCLUSIONS

Even with present prediction approaches, improvement in prediction accuracy can still be expected, albeit limited to no more than 10%.

GIF-DB, a WWW database on gene interactions involved in Drosophila melanogaster development

GIF-DB (Gene Interactions in the Fly Database) is a new WWW database (http://www-biol.univ-mrs.fr/ approximately lgpd/GIFTS_home_page. html ) describing gene molecular interactions involved in the process of embryonic pattern formation in the flyDrosophila melanogaster. The detailed information is distributed in specific lines arranged into an EMBL- (or SWISS-PROT-) like format. GIF-DB achieves a high level of integration with other databases such as FlyBase, EMBL and SWISS-PROT through numerous hyperlinks. The original concept of interaction databases examplified by GIF-DB could be extended to other biological subjects and organisms so as to study gene regulatory networks in an evolutionary perspective.

The EMBL Nucleotide Sequence Database

The EMBL Nucleotide Sequence Database is a comprehensive database of DNA and RNA sequences directly submitted from researchers and genome sequencing groups and collected from the scientific literature and patent applications. In collaboration with DDBJ and GenBank the database is produced, maintained and distributed at the European Bioinformatics Institute (EBI) and constitutes Europe's primary nucleotide sequence resource. Database releases are produced quarterly and are distributed on CD-ROM. EBI's network services allow access to the most up-to-date data collection via Internet and World Wide Web interface, providing database searching and sequence similarity facilities plus access to a large number of additional databases.

Database on the structure of large ribosomal subunit RNA

The latest release of the large ribosomal subunit RNA database contains 429 sequences. All these sequences are aligned, and incorporate secondary structure information. The rRNA WWW Server at URL http://rrna.uia.ac.be/ provides researchers with an easily accessible resource to obtain the data in this database in a number of computer-readable formats. A new query interface has been added to the server. If necessary, the data can also be obtained by anonymous ftp from the same site.

MmtDB: a Metazoa mitochondrial DNA variants database

The present paper describes the structure of MmtDB-a specialized database designed to collect Metazoa mitochondrial DNA variants. Priority in the data collection is given to the Metazoa species for which a large amount of variants is available, as it is the case for human variants. Starting from the sequences available in the Nucleotide Sequence Databases, the redundant sequences are removed and new sequences from other sources are added. Value-added information are 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's gain or loss). Cross-references are introduced to the EMBL Data Library, as well as an internal cross-referencing among MmtDB entries according to their tissual, heteroplasmic, familiar and aplotypical correlation. MmtDB can be accessed through the World Wide Web at URL [see text].

A novel algorithm for the search of 5S rRNA genes in DNA databases: comparison with other methods and identification of new potential 5S rRNA genes

We report here a new algorithm for the identification of 5S rRNA genes in DNA databases. Based on an improved version of the general weight matrix method, this search procedure relies on the recognition of three informative regions within 5S rRNA genes, and on the weighted evaluation of the distance between them. As an additional step, the algorithm extends the weight matrix analysis to the full-length 5S rRNA sequence. This combined strategy, which includes a fast, but poorly selective, preliminary search procedure and an auxiliary step, that is slow but highly selective, strongly reduces the number of false positive instances, yielding a total false positive rate of 0.00076%. On the other hand, 97.5% of the 1045 known 5S rRNA genes were correctly recognized by this algorithm, and 29 previously unidentified potential 5S rRNA sequences were uncovered. A detailed analysis of these candidate sequences, including prediction of 5S rRNA secondary structure and checking for the presence of transcriptional termination signals, showed that eight of them correspond to authentic 5S rRNA genes. The performance of this specialized algorithm for the detection of 5S rRNA genes was compared with that of the general hidden Markov model search procedure. Due to their utilization of different filtering rules, the two approaches proved to be highly complementary. Their combined use will thus provide a very effective tool for the detection of dispersed 5S rRNA genes, either active or inactive, in the vertebrate genome.

A protein of the SR family of splicing factors binds extensively to exonic Balbiani ring pre-mRNA and accompanies the RNA from the gene to the nuclear pore

We report on the molecular cloning and intracellular localization of a heterogeneous nuclear ribonucleoprotein (hnRNP), Ct-hrp45, one of the major components of pre-mRNP particles in Chironomus tentans. It is shown that hrp45 belongs to the SR family of splicing factors and exhibits high sequence similarity to Drosophila SRp55/B52 and human SF2/ASF. The distribution of hrp45 within the C. tentans salivary gland cells is studied by immunocytology. The hrp45 protein is found to be abundant in the nucleus, whereas it is undetectable in the cytoplasm. The fate of hrp45 in specific pre-mRNP particles, the Balbiani ring (BR) granules, is revealed by immunoelectron microscopy. It is observed that hrp45 is associated with the growing BR pre-mRNP particles and is being added continuously concomitant with the growth of the transcript, indicating that hrp45 is bound extensively to exon 4, which comprises 80-90% of the primary transcript. Furthermore, hrp45 remains bound to the BR RNP particles in the nucleoplasm and is not released until the particles translocate through the nuclear pore. Thus, hrp45 behaves as an hnRNP protein linked to exon RNA (and perhaps also to the introns) rather than as a spliceosome component connected to the assembly and disassembly of spliceosomes. It seems that hrp45, and possibly also other SR family proteins, is playing an important role in the structural organization of pre-mRNP particles and is perhaps participating not only in splicing but also in other intranuclear events.

Friedreich's ataxia protein: phylogenetic evidence for mitochondrial dysfunction

Friedreich's ataxia is the most common inherited spinocerebellar ataxia. A decade of linkage and physical mapping studies have culminated in the identification of the Friedreich's ataxia gene. The presence of homologues in purple bacterial genomes, but not in other bacteria, allows us to infer a mitochondrial location for frataxin (Friedreich's ataxia protein) on the basis of bacterial phylogeny. Frataxin possesses a non-globular N-terminus domain providing a candidate mitochondrial targeting peptide. Clues to the function of frataxin are provided by the mitochondrial location, a clinically similar ataxia with vitamin E deficiency, and certain neuropathies with mitochondrial DNA instability caused by mutations in nuclear genes.

Protein secondary structural types are differentially coded on messenger RNA

Tricodon regions on messenger RNAs corresponding to a set of proteins from Escherichia coli were scrutinized for their translation speed. The fractional frequency values of the individual codons as they occur in mRNAs of highly expressed genes from Escherichia coli were taken as an indicative measure of the translation speed. The tricodons were classified by the sum of the frequency values of the constituent codons. Examination of the conformation of the encoded amino acid residues in the corresponding protein tertiary structures revealed a correlation between codon usage in mRNA and topological features of the encoded proteins. Alpha helices on proteins tend to be preferentially coded by translationally fast mRNA regions while the slow segments often code for beta strands and coil regions. Fast regions correspondingly avoid coding for beta strands and coil regions while the slow regions similarly move away from encoding alpha helices. Structural and mechanistic aspects of the ribosome peptide channel support the relevance of sequence fragment translation and subsequent conformation. A discussion is presented relating the observation to the reported kinetic data on the formation and stabilization of protein secondary structural types during protein folding. The observed absence of such strong positive selection for codons in non-highly expressed genes is compatible with existing theories that mutation pressure may well dominate codon selection in non-highly expressed genes.

Ribosome-mediated translational pause and protein domain organization

Because regions on the messenger ribonucleic acid differ in the rate at which they are translated by the ribosome and because proteins can fold cotranslationally on the ribosome, a question arises as to whether the kinetics of translation influence the folding events in the growing nascent polypeptide chain. Translationally slow regions were identified on mRNAs for a set of 37 multidomain proteins from Escherichia coli with known three-dimensional structures. The frequencies of individual codons in mRNAs of highly expressed genes from E. coli were taken as a measure of codon translation speed. Analysis of codon usage in slow regions showed a consistency with the experimentally determined translation rates of codons; abundant codons that are translated with faster speeds compared with their synonymous codons were found to be avoided; rare codons that are translated at an unexpectedly higher rate were also found to be avoided in slow regions. The statistical significance of the occurrence of such slow regions on mRNA spans corresponding to the oligopeptide domain termini and linking regions on the encoded proteins was assessed. The amino acid type and the solvent accessibility of the residues coded by such slow regions were also examined. The results indicated that protein domain boundaries that mark higher-order structural organization are largely coded by translationally slow regions on the RNA and are composed of such amino acids that are stickier to the ribosome channel through which the synthesized polypeptide chain emerges into the cytoplasm. The translationally slow nucleotide regions on mRNA possess the potential to form hairpin secondary structures and such structures could further slow the movement of ribosome. The results point to an intriguing correlation between protein synthesis machinery and in vivo protein folding. Examination of available mutagenic data indicated that the effects of some of the reported mutations were consistent with our hypothesis.

Submission of nucleotide sequence data to EMBL/GenBank/DDBJ

This review outlines the various methods available for submitting sequence data to the EMBL Nucleotide Sequence Database. Depending on the type of sequence data and the facilities available to the submitter, one method may be more suitable than another. Recent developments have been the World Wide Web submission tool, procedures for bulk submissions and genome projects.

Bacterial 16S rDNA polymerase chain reaction in the detection of intra-amniotic infection

OBJECTIVE

Bacterial polymerase chain reaction (PCR) was used to detect early subclinical intraamhiotic infection. We used universal primers which amplify a DNA fragment of 16S ribosomal DNA (rDNA) from all known bacteria and sequenced the positive samples to identify the bacterial species.

DESIGN

Transabdominally obtained amniotic fluid samples from 20 pregnant women with prelabour rupture of the fetal membranes (PROM), showing no signs of clinical infection, and 16 control samples were analysed with universal bacterial PCR. In addition, routine bacterial culture and amniotic fluid glucose were studied.

RESULTS

Out of 20 PROM patients, five were positive in the PCR. PCR detected Ureaplasma urealyticum in two cases, Haemophilus influenzae in one case, Streptococcus oralis in one case and Fusobacterium sp. in one case. Only two of these were positive in a routine bacterial culture. Both were multibacterial infections, which caused discrepancies between the PCR and culture results. Two patients developed infectious complications: both were identified with the PCR assay. Amniotic fluid glucose was lower in PCR positive patients compared with PCR negative patients.

CONCLUSION

Bacterial 16S rDNA PCR, in properly controlled conditions, promises to be a fast and reliable test for early intra-amniotic infection especially concerning Ureaplasma urealyticum.

Conservation of aconitase residues revealed by multiple sequence analysis. Implications for structure/function relationships

Aconitases have recently regained much attention, because one member of this family, iron regulatory protein-1 (IRP-1), has been found to play a dual role as a cytoplasmic aconitase and a regulatory RNA-binding protein. This finding has highlighted a novel role for Fe-S clusters as post-translational regulatory switches. We have aligned 28 members of the Fe-S isomerase family, identified highly conserved amino acid residues, and integrated this information with data on the crystallographic structure of mammalian mitochondrial aconitase. We propose structural and/or functional roles for the previously unrecognized conserved residues. Our findings illustrate the value of detailed protein sequence analysis when high-resolution crystallographic data are already available.

Lambda clone B22 contains a 7676 bp genomic fragment of Saccharomyces cerevisiae chromosome VII spanning the VAM7-SPM2 intergenic region and containing three novel transcribed open reading frames

A genomic clone of 7676 bp designated B22 from Saccharomyces cerevisiae has been sequenced. The 5' end matches the previously described gene, VAM7, and the 3' end matches the previously described gene, SPM2, both of which have been assigned to the left arm of chromosome VII. The intergenic region contains three transcribed open reading frames (ORFs). The first is related to an uncharacterized ORF of Bacillus subtilis and more weakly to MesJ in Escherichia coli; this is found as a single transcript of 1.1 kb by Northern blotting. The second ORF encodes a small ras-like GTPase of 222 residues with strong homology to yeast Ypt8p and to mammalian Rab11; this is found as a single transcript of 1.1 kb by Northern blotting. The third ORF generates a transcript of 1.6 kb and encodes a protein of 382 residues including a perfect match to the consensus sequence of a C2H2 zinc finger domain; it shares a strong homology with yeast Mig1p and Cre-A from Aspergillus, Emericella and E. coli. This ORF also has a striking similarity to a putative 43 kDa zinc finger protein encoded by an ORF (YEL8) immediately downstream of YPT8, raising the possibility that a region between VAM7 and SPM2 on chromosome VII arose as a duplication of the YPT8-YEL8 region of chromosome V, followed by a translocation.

The sequence of a 21.3 kb DNA fragment from the left arm of yeast chromosome XIV reveals LEU4, MET4, POL1, RAS2, and six new open reading frames

The nucleotide sequence of a fragment of 21 308 bp from the left arm of Saccharomyces cerevisiae chromosome XIV has been determined. Analysis of the sequence revealed 13 open reading frames (ORFs) longer than 300 bp, four of which correspond to the previously identified genes LEU4, MET4, POL1 and RAS2. One putative protein, N2160, shares considerable homology (32% identity) with a hypothetical protein encoded by a gene located on chromosome XV as well as with human OCRL protein (36% identity), involved in Lowe's syndrome. N2185 contains ten predicted transmembrane segments and is similar to another putative protein (YKL146) from yeast.

A pre-mRNA-binding protein accompanies the RNA from the gene through the nuclear pores and into polysomes

In the larval salivary glands of C. tentans, it is possible to visualize by electron microscopy how Balbiani ring (BR) pre-mRNA associates with proteins to form pre-mRNP particles, how these particles move to and through the nuclear pore, and how the BR RNA is engaged in the formation of giant polysomes in the cytoplasm. Here, we study C. tentans hrp36, an abundant protein in the BR particles, and establish that it is similar to the mammalian hnRNP A1. By immuno-electron microscopy it is demonstrated that hrp36 is added to BR RNA concomitant with transcription, remains in nucleoplasmic BR particles, and is translocated through the nuclear pore still associated with BR RNA. It appears in the giant BR RNA-containing polysomes, where it remains as an abundant protein in spite of ongoing translation.

Information services of the European Bioinformatics Institute

The scope of the EBI is focused on providing better services to the scientific community. Technological advancements in the hardware area provide EBI with means of producing data much faster than before, and with greater accuracy since there is now a better technical ability to produce more exhaustive searches through larger indices. Hand in hand with the technological developments, research and development work is continuing on better indexing systems and more efficient ways of establishing and maintaining the future databases. The existing links of communication between EBI and the user community are exploited to study the needs of the scientific community, to provide better services, and to enhance the quality of databases by interpreting user feedback and updates. A very important goal is to enhance the awareness of the scientific (and, maybe even more, the nonscientific) public of the importance of the modern field of bioinformatics and to introduce special meetings and courses, in which more specific subjects will be studied in depth. Another aspect of this goal is to help in constructing special bioinformatics programs in university faculties. In such programs, in contrast to the existing layout, students will pursue studies in a combined environment that provides basic training in biology and in computation. Currently, one of the main problems in the field is that scientists are either biologists, who are self-educated in the field of computers and programming, or computer scientists without sufficient knowledge of biology. It is hoped that a combined program will provide a high level of education in both fields of interest at the appropriate ratios. Building an efficient and friendly interface between the EBI and the user community is the basis for any future development. This aim is achieved by using the most modern server systems while continuously researching newer and better systems and interfaces. This task can never be complete without involvement of the user community by providing feedback to any of EBI's services. A better bioinformatics community is a necessity for any future development of the biological research aiming at a better society.

TransTerm: a database of translational signals

The TransTerm database of sequence contexts of stop and start codons has been expanded to include approximately 50% more species than last year's release. It now contains 148 organisms and >39 500 coding sequences; it is now available on the World Wide Web. The database includes: (i) initiation and termination sequence contexts organized by species; (ii) summary parameters about the individual sequences (sequence length, GC%, GC3, Nc, CAI) in addition to tables of base frequencies for each species' stop and start codon sequence context; (iii) species codon usage tables; and (iv) summary tables of stop signal frequency.

The European Bioinformatics Institute (EBI) databases

The European Bioinformatics Institute (EBI) maintains and distributes the EMBL Nucleotide Sequence database, Europe's primary nucleotide sequence data resource. The EBI also maintains and distributes the SWISS-PROT Protein Sequence database, in collaboration with Amos Bairoch of the University of Geneva. Over fifty additional specialist molecular biology databases, as well as software and documentation of interest to molecular biologists are available. The EBI network services include database searching and sequence similarity searching facilities.

Database on the structure of large ribosomal subunit RNA

Our database on large ribosomal subunit RNA contained 334 sequences in July, 1995. All sequences in the database are aligned, taking into account secondary structure. The aligned sequences are provided, together with incorporated secondary structure information, in several computer-readable formats. These data can easily be obtained through the World Wide Web. The files in the database are also available via anonymous ftp.

Compilation of DNA sequences of Escherichia coli K12 (ECD and ECDC; update 1995)

We have compiled the DNA sequence data for Escherichia coli available from the GenBank and EMBL data libraries and independently from the literature. Unlike the previous updates of our E.coli databases, we provide the most recent version preferentially via the World Wide Web System (use URL: http://susi.bio.unigiessen.de/usr/local/www++ +/html/ecdc.html). Our database includes an assembled set of contiguous sequences. Each of these contigs compiles all available sequence information, including those derived from a variety of elder sequences. The organization of the database allows one to find the exact physical location of each individual gene or regulatory region, even regarding discrepancies in nomenclature. The WWW program allows access into the original EMBL and SWISSPROT datafiles. A FASTA and BLAST search may be performed online. Besides the WWW format a flat file version may be obtained via ftp. The complete compilation, including a full set of genetic map data and the E.coli protein index, can be obtained in machine readable form from the EMBL data library as a part of the CD-ROM issue of the EMBL sequence database, released and updated every three months. After deletion of all detected overlaps a total of 3 333 878 individual bp was determined by the end of September 1995. This corresponds to a total of 71.71% of the entire E.coli chromosome consisting of about 4720 kbp. About 94 kbp (2%) are available additionally, but have not yet been definitely mapped.

DNA and ascomycete systematics

Present knowledge of different types of RNAs as phylogenetic markers among the fungi is discussed, and examples of phylogenetically informative 18S rRNA signature sequences are given. Such signatures give phylogenetic information that is not provided by parsimony or distance analyses of longer gene sequences. A single signature cannot be used as a decisive criterion for defining taxa, but signature sequences give invaluable hints on phylogenetic relationships and can be included in data matrices as morphological criteria when using parsimony analysis. Key words: DNA, ascomycete phylogeny, signature sequences.

The chromo superfamily: new members, duplication of the chromo domain and possible role in delivering transcription regulators to chromatin

Using computer methods for detecting conserved amino acid sequence motifs, we show that the chromatin organization modifier (chromo) domain that has been previously identified in several proteins involved in transcription down-regulation is present in a much larger group of (putative) chromatin-binding proteins, some of which are positive rather than negative regulators of transcription. The most interesting new members of the chromo superfamily are Drosophila male-specific lethal (MSL-3) protein involved in the X chromosome gene dosage compensation in the males and human retinoblastoma-binding protein RBP-1. We show that the chromo domain is duplicated in several chromatin-binding proteins and use this observation to interpret recent results on chromatin binding obtained with chimeric chromo domain-containing proteins. We hypothesize that the chromo domain may be a vehicle that delivers both positive and negative transcription regulators to the sites of their action on chromatin.

Evolution of energetic metabolism: the respiration-early hypothesis

The main energy-transducing metabolic systems originated and diversified very early in the evolution of life. This makes it difficult to unravel the precise steps in the evolution of the proteins involved in these processes. Recent molecular data suggest that homologous proteins of aerobic respiratory chains can be found in Bacteria and Archaea, which points to a common ancestor that possessed these proteins. Other molecular data predict that this ancestor was unlikely to perform oxygenic photosynthesis. This evidence, that aerobic respiration has a single origin and may have evolved before oxygen was released to the atmosphere by photosynthetic organisms, is contrary to the textbook viewpoint.

Quantitative fluorescence in situ hybridization of Bifidobacterium spp. with genus-specific 16S rRNA-targeted probes and its application in fecal samples

Three 16S rRNA hybridization probes were developed and tested for genus-specific detection of Bifidobacterium species in the human fecal flora. Variable regions V2, V4, and V8 of the 16S rRNA contained sequences unique to this genus and proved applicable as target sites for oligodeoxynucleotide probes. Determination of the genus specificity of the oligonucleotides was performed by whole-cell hybridization with fluorescein isothiocyanate-labelled probes. To this end, cells were fixed on glass slides, hybridized with the probes, and monitored by videomicroscopy. In combination with image analysis, this allowed quantification of the fluorescence per cell and objective evaluation of hybridization experiments. One of the probes developed was used to determine the population of Bifidobacterium spp. in human fecal samples. A comparison was made with results obtained by cultural methods for enumeration. Since both methods gave similar population estimates, it was concluded that all bifidobacteria in feces were culturable. However, since the total culturable counts were only a fraction of the total microscopic counts, the contribution of bifidobacteria to the total intestinal microflora was overestimated by almost 10-fold when cultural methods were used as the sole method for enumeration.

A new obligately chemolithoautotrophic, nitrite-oxidizing bacterium, Nitrospira moscoviensis sp. nov. and its phylogenetic relationship

A gram-negative, non-motile, non-marine, nitrite-oxidizing bacterium was isolated from an enrichment culture initiated with a sample from a partially corroded area of an iron pipe of a heating system in Moscow, Russia. The cells were 0.9-2.2 microns x 0.2-0.4 microns in size. They were helical- to vibroid-shaped and often formed spirals with up to three turns 0.8-1.0 micron in width. The organism possessed an enlarged periplasmic space and lacked intracytoplasmic membranes and carboxysomes. The cells tended to excrete extracellular polymers, forming aggregates. The bacterium grew optimally at 39 degrees C and pH 7.6-8.0 in a mineral medium with nitrite as sole energy source and carbon dioxide as sole carbon source. The optimal nitrite concentration was 0.35 mM. Nitrite was oxidized to nitrate stoichiometrically. The doubling time was 12 h in a mineral medium with 7.5 mM nitrite. The cell yield was low; only 0.9 mg protein/l was formed during oxidation of 7.5 mM nitrite. Under anoxic conditions, hydrogen was used as electron donor with nitrate as electron acceptor. Organic matter (yeast extract, meat extract, peptone) supported neither mixotrophic nor heterotrophic growth. At concentrations as low as 0.75 g organic matter/l or higher, growth of nitrite-oxidizing cells was inhibited. The cells contained cytochromes of the b- and c-type. The G+C content of DNA was 56.9 +/- 0.4 mol%. The chemolithoautotrophic nitrite-oxidizer differed from the terrestrial members of the genus Nitrobacter with regard to morphology and substrate range and equaled Nitrospira marina in both characteristics. The isolated bacterium is designated as a new species of the genus Nitrospira.(ABSTRACT TRUNCATED AT 250 WORDS)

Culture and phylogenetic characterization of Trichomitus trypanoides Duboscq & Grassè 1924, n. comb.: a trichomonad flagellate isolated from the hindgut of the termite Reticulitermes santonensis Feytaud

A trichomonad flagellate strain R1 was isolated from the hindgut contents of the termite Reticulitermes santonensis Feytaud. The flagellate was cultivated at 28 degrees C in anaerobic medium containing yeast extract, minerals and vitamins. The isolate fed on living bacteria. It showed the typical morphological and ultrastructural features of the trichomonads, closely resembling Trichomitus trypanoides. In order to determine its phylogenetic position the small subunit ribosomal DNA (SSU rDNA) of the flagellate was amplified in vitro using the polymerase chain reaction (PCR), cloned in a plasmid vector and sequenced. Comparison of the obtained sequence with so far available SSU rRNA/rDNA sequences showed strongest similarity (89%) to the sequence of Tritrichomonas foetus. The phylogenetic analysis with parsimony and distance matrix methods placed Trichomitus trypanoides strain R1 near by the root of the phylogenetically so far analyzed eukaryotic organisms. This confirms that termites harbour hindgut symbionts, which originate from very early evolved eukaryotes.

Alternative mRNA splicing and differential promoter utilization determine tissue-specific expression of the apolipoprotein B mRNA-editing protein (Apobec1) gene in mice. Structure and evolution of Apobec1 and related nucleoside/nucleotide deaminases

Apolipoprotein (apo) B mRNA editing consists of a C-->U conversion involving the first base of the codon CAA, encoding Gln 2153, to UAA, a stop codon. Editing occurs in the intestine only in most mammals, and in both the liver and intestine in a few mammalian species including mouse. We have cloned the cDNA for the mouse apoB mRNA editing protein, apobec1. Expression of mouse apobec1 cDNA in HepG2 cells results in the editing of the intracellular apoB mRNA. The cDNA predicts a 229-amino acid protein showing 92, 66, and 70% identity to the rat, rabbit, and human proteins, respectively. Based on the estimated values of divergence of apobec1 sequences in terms of the numbers of synonymous and non-synonymous suhstitutions per site, we found that apobec1 is a fairly rapidly evolving protein. Sequence comparison among mammalian apobec1 sequences has permitted the identification of seven conserved regions that may be functionally important for editing activity. We present a phylogenetic tree relating apobec1 sequences to double-stranded RNA adenosine deaminase and other nucleotide/nucleoside deaminases. Northern blot analysis indicates that apobec1 mRNA exists in two different sizes, a approximately 2.2-kilobase (kb) form in small intestine and a approximately 2.4-kb form in liver, spleen, kidney, lung, muscle, and heart. To study the molecular basis for the different sized apobec1 mRNAs, we cloned the apobec1 gene and characterized its exon-intron organization together with the sequences expressed in the hepatic and intestinal mRNA. The mouse apobec1 gene contains 8 exons and spans approximately 25 kb, and is located in chromosome 6. The major hepatic mRNA contains all 8 exons, whereas the major small intestinal mRNA misses the first 3 exons and its transcription is initiated in exon 4. The intestinal mRNA also contains at its 5' end a unique 102-nucleotide piece that is absent in the liver mRNA. We also identified two alternatively spliced hepatic apobec1 mRNAs with different acceptor sites in exon 4. Transient expression studies using promoter-reporter gene constructs in HeLa, Hepa, and Caco-2 cells indicate that the 5'-flanking sequences of the liver mRNA (i.e. upstream of exon 1) have predominantly hepatic promoter activity and the 5'-flanking sequences of the major small intestine mRNA (i.e. upstream of exon 4) have preferential intestinal promoter activity.(ABSTRACT TRUNCATED AT 400 WORDS)