The EMBL data library

Best practices in bioinformatics training for life scientists

The mountains of data thrusting from the new landscape of modern high-throughput biology are irrevocably changing biomedical research and creating a near-insatiable demand for training in data management and manipulation and data mining and analysis. Among life scientists, from clinicians to environmental researchers, a common theme is the need not just to use, and gain familiarity with, bioinformatics tools and resources but also to understand their underlying fundamental theoretical and practical concepts. Providing bioinformatics training to empower life scientists to handle and analyse their data efficiently, and progress their research, is a challenge across the globe. Delivering good training goes beyond traditional lectures and resource-centric demos, using interactivity, problem-solving exercises and cooperative learning to substantially enhance training quality and learning outcomes. In this context, this article discusses various pragmatic criteria for identifying training needs and learning objectives, for selecting suitable trainees and trainers, for developing and maintaining training skills and evaluating training quality. Adherence to these criteria may help not only to guide course organizers and trainers on the path towards bioinformatics training excellence but, importantly, also to improve the training experience for life scientists.

PSSARD (2.0): a database server for making flexible queries relating amino acid sequences to main-chain secondary structure conformations for proteins of known three-dimensional structure and certain useful applications

We have updated the Protein Sequence-Structure Analysis Relational Database (PSSARD) first published in the Int. J. Biol. Macromol. 36 (2005) 259-262 corresponding to 1573 representative protein chains selected from the Protein Data Bank (PDB). In this, the updated and revised PSSARD (Version 2.0), we have included all proteins in the Protein Data Bank available at the time of developing this database including the NMR PDB entries. The current database corresponds to 22,752 XRAY PDB entries and 3977 NMR PDB entries and is separated accordingly in order to facilitate the appropriate database search. The representative protein chains can also be separately accessed within the current database. We have made a provision to combine more than one field to query the database and the results of any search can be used to carry out further nested searches using a combination of queries. We have provided hyperlinks to the individual PDB entries obtained as the result of any search in PSSARD in order to obtain additional details relevant to the protein structure. Certain applications useful to identify domains and structural motifs are discussed.

PSSARD: protein sequence-structure analysis relational database

We have implemented a relational database comprising a representative dataset of amino acid sequences and their associated secondary structure. The representative amino acid sequences were selected according to the PDB_SELECT program by choosing proteins corresponding to protein crystal structure data deposited in the protein data bank that share less than 25% overall pair-wise sequence identity. The secondary structure was extracted from the protein data bank website. The information content in the database includes the protein description, PDB code, crystal structure resolution, total number of amino acid residues in the protein chain, amino acid sequence, secondary structure conformation and its summary. The database is freely accessible from the website mentioned below and is useful to query on any of the above fields. The database is particularly useful to quickly retrieve amino acid sequences that are compatible to any super-secondary structure conformation from several proteins simultaneously.

A relationship between GC content and coding-sequence length

Since base composition of translational stop codons (TAG, TAA, and TGA) is biased toward a low G+C content, a differential density for these termination signals is expected in random DNA sequences of different base compositions. The expected length of reading frames (DNA segments of sense codons flanked by in-phase stop codons) in random sequences is thus a function of GC content. The analysis of DNA sequences from several genome databases stratified according to GC content reveals that the longest coding sequences-exons in vertebrates and genes in prokaryotes-are GC-rich, while the shortest ones are GC-poor. Exon lengthening in GC-rich vertebrate regions does not result, however, in longer vertebrate proteins, perhaps because of the lower number of exons in the genes located in these regions. The effects on coding-sequence lengths constitute a new evolutionary meaning for compositional variations in DNA GC content.

The structure of adenovirus type 12 DNA integration sites in the hamster cell genome

Foreign DNA can integrate into the genomes of mammalian cells, and this process plays major roles in viral oncogenesis and in the generation of transgenic organisms and will be important in evolving regimens for human somatic gene therapy. In the present study, the insertion sites of adenovirus type 12 (Ad12) DNA genomes have been analyzed in detail in the Ad12-transformed hamster cell line T637, its revertants, which have lost most of the >20 Ad12 genome equivalents integrated chromosomally in cell line T637, and in the Ad12-induced tumor T191. Some of these junction sites have been molecularly cloned, and the nucleotide sequences at the sites of transition between viral and cellular DNAs have been determined. The sites of linkage between the hamster cellular and the foreign (viral) DNA are characterized by the frequent occurrence of patch homologies between the recombination partners. The cellular junction sites investigated here are not transcriptionally active. One of the cellular DNA sequences abutting the right Ad12 DNA terminus in cell line T637 (os2) is represented only once in the hamster genome and has a strikingly low abundance of 5'-CG-3' dinucleotide sequences. One 5'-GCGC-3' sequence close to the Ad12 DNA integration site is heavily methylated in normal cells, Ad12-transformed cells, and Ad12-induced tumor cells. The second such sequence is more remote from the junction site, is partly methylated in BHK21 hamster cells, and shows differences in methylation in different Ad12-transformed cell lines. This site is unmethylated in liver DNA. The cellular DNA sequence at the site of Ad12 linkage in the tumor T191 exhibits homologies to highly repetitive sequences of the Alu family and to an origin of hamster DNA replication containing an Alu element. A number of junction sites between Ad12 DNA and hamster or mouse DNA in Ad12-transformed cell lines or Ad12-induced tumor cell lines, investigated here and previously, are characterized by stem-loop structures encompassing the junction sites.

Human T-cell receptor variable gene segment families

Multiple DNA and protein sequence alignments have been constructed for the human T-cell receptor alpha/delta, beta, and gamma (TCRA/D, B, and G) variable (V) gene segments. The traditional classification into subfamilies was confirmed using a much larger pool of sequences. For each sequence, a name was derived which complies with the standard nomenclature. The traditional numbering of V gene segments in the order of their discovery was continued and changed when in conflict with names of other segments. By discriminating between alleles at the same locus versus genes from different loci, we were able to reduce the number of more than 150 different TCRBV sequences in the database to a repertoire of only 47 functional TCRBV gene segments. An extension of this analysis to the over 100 TCRAV sequences results in a predicted repertoire of 42 functional TCRAV gene segments. Our alignment revealed two residues that distinguish between the highly homologous V delta and V alpha, one at a site that in VH contacts the constant region, the other at the interface between immunoglobulin VH and VL. This site may be responsible for restricted pairing between certain V delta and V gamma chains. On the other hand, V beta and V gamma appear to be related by the fact that their CDR2 length is increased by four residues as compared with that of V alpha/delta peptides.

Mouse T-cell receptor variable gene segment families

All mouse T-cell receptor alpha/delta, beta, and gamma variable (Tcra/d-, b-, and g-V) gene segments were aligned to compare the sequences with one another, to group them into subfamilies, and to derive a name which complies with the standard nomenclature. It was necessary to change the names of some V gene segments because they conflicted with those of other segments. The traditional classification into subfamilies was re-evaluated using a much larger pool of sequences. In the mouse, most V gene segments can be grouped into subfamilies of closely related genes with significantly less similarity between different subfamilies.

Tandemly repeated pentanucleotides in DNA sequences of eucaryotes

Genetic sequence data banks were scanned in order to retrieve tandemly repeated pentanucleotides (pnts). It was found that among 102 (=(1024-4)/2/5) possible distinct pnts roughly each fourth is involved in tandem repeats. It is shown that tandemly repeated pnts are composed of frequently occurring di- and trinucleotides and that those pnts which occur frequently in the form of mono- or di-pnts form also tandem repeats either in the form of satellites or in the form of shorter tandem repeats. Human satellite III is taken as a specific example. It is shown that the first guanine within GG-AAT pnt exhibits the highest mutability. Sequential distribution of base changes gives evidence that the mutations do not occur at random positions but in a correlated fashion so that long stretches of original pnts remain intact. It is found that pnts related to the satellite III are present in introns and flanking regions of some structural genes, but are not preserved between orthologous genes of related species. The results corroborate the most plausible mechanism of their evolution--rapid amplification followed by successive divergence of repeat units by various mutational processes.

A deductive database system for analyzing human nucleotide sequence data

The analysis of the human genome is one of the most significant topics in both biology and medical science. There is a growing need for a well-designed database system for searching and analyzing the human genome data. We developed a deductive database system to search and analyze nucleotide sequence data derived from the GenBank primates data. A deductive database system is a next generation one and it contains an inference mechanism that can handle problems beyond the capabilities of classical database systems. Database queries are described in logical rules. These rules are simple even for molecular biologists who are not experts in computer programs because they are declarative and do not require the procedural commands that are usually used in computer programs. Furthermore, queries based on logical rules are powerful enough to express complicated biological problems. Particularly, recursive rules are suitable for examining secondary structures of nucleotide sequences. In our analysis of TfR's IRE, we noted five stem-and-loop structures.

Barley beta-glucosidase: expression during seed germination and maturation and partial amino acid sequences

Unlike most of the hydrolytic enzymes that participate in endosperm mobilization, beta-glucosidase of barley (Hordeum vulgare) seeds does not increase during germination, even in the presence of exogenously added gibberellic acid. However, the germination process affects the physical properties of beta-glucosidase in terms of charge and apparent molecular weight. Analysis of developing barley grains shows that the enzyme is synthesized two weeks before maturation and is stored in the endosperm of the dry dormant seed. Partial amino acid sequencing of the purified beta-glucosidase demonstrates significant similarity between the barley enzyme and beta-glycosidases that belong to family 1 of glycosyl hydrolases.

Nucleotide sequence of human adenovirus type 12 DNA: comparative functional analysis

A fresh inoculum of human adenovirus type 12 (Ad12) was obtained from the American Type Culture Collection and passaged once on human embryonic kidney cells, and Ad12 DNA was prepared from the first-passage yield to avoid higher passages which might have generated host-virus DNA recombinants. The 18 PstI fragments of Ad12 DNA were cloned into the pBluescript KS vector, and the entire nucleotide sequence of both strands from all 18 fragments was determined by using successive oligodeoxyribonucleotide primers. Ad12 DNA extends over 34,125 nucleotide pairs, and its molecular weight is calculated to be about 22 x 10(6). The nucleotide sequence of Ad12 DNA was subjected to computer analyses that determined possible open reading of frames on the two strands, the leader sequences, the position of the virus-associated RNA coding region, possible TATA, and polyadenylation signals. The distribution of the Ad12 open reading frames was similar to that in the previously sequenced Ad2 DNA, but there were also distinct differences. Ad12 DNA has an inverted terminal redundancy of 161 nucleotides, compared with 102 nucleotides in Ad2 DNA. There were stretches of sequence identity between Ad2 and Ad12 DNAs at both termini; the overall sequence similarity between the two viral genomes ranged between 59% (polypeptide IX) and 77% (in the E2 region), with high homology also in the sequences for the adenovirus DNA polymerase.

Dinucleotides and G+C content in human genes: opposite behavior of GpG, GpC, and TpC at II-III codon positions and in introns

We have studied the behavior of the dinucleotide preferences under G+C content variation in human genes. The doublet preferences for each dinucleotide were compared between two functionally distinct zones in genes, the II-III codon positions, and the introns. The 16 dinucleotides have been tentatively classified in three groups: AA, AC, CC, CT, and GA, doublets showing no difference between introns and II-III codon positions in the full range of G+C variation TG and TA, which differ in the full range of G+C variation AT, AG, GT, TC, TT, GG, GC, CG, and CA, which show differences in regions over 50% G+C A remarkable pattern observed concerns the behavior of GG, GC, and TC, which showed opposite trends in II-III codon positions and in introns. If codon positions and introns are under the same structural requirements and the same mutational bias, our results indicate that the differences observed could be related to post-transcriptional constraints acting on mRNA.

A pathogen-induced gene of barley encodes a HSP90 homologue showing striking similarity to vertebrate forms resident in the endoplasmic reticulum

The full-length nucleotide sequence of a barley (Hordeum vulgare L.) leaf mRNA, found to increase rapidly in amount during infection attempts by the powdery mildew fungus (Erysiphe graminis DC. ex Mérat), is reported. The mRNA encodes a polypeptide of 809 amino acid residues which, by sequence comparison, was identified as a member of the 90 kDa heat shock protein (HSP90) family. The encoded protein most resembles the endoplasmic reticulum (ER) resident HSP90 protein, the 94 kDa glucose-regulated protein (GRP94) of vertebrates, as it possesses both the characteristic N-terminal domain including a signal peptide sequence and the C-terminal ER retention signal (Lys-Asp-Glu-Leu). A transcript cross-hybridizing at high stringency accumulated rapidly in leaves upon heat shock treatment. Genomic DNA blot analysis indicated the presence of a family of related genes in the barley genome.

The hypervariable DXS255 locus contains a LINE-1 repetitive element with a CpG island that is extensively methylated only on the active X chromosome

The DXS255 locus at Xp11.22 is highly polymorphic due to a 26-bp variable number of tandem repeats (VNTR) motif. In previous studies, one of the MspI sites flanking the VNTR manifested a correlation between methylation and X chromosome inactivation. Here we show, by DNA sequence analysis, that this MspI site is located within the CpG island at the 5' end of a LINE-1 element, which is 2.5 kb from the VNTR. The methylation status of the CpG island was assessed in Southern blotting experiments using the methylation-sensitive enzymes HpaII, HhaI, and BssHII. All these sites were completely methylated on active X chromosomes, consistent with previously reported findings of full methylation of LINE-1 elements throughout the genome. However, on inactive X chromosomes these sites were predominantly unmethylated, although patterns were found to be heterogeneous. The results suggest that LINE-1 elements on the inactive X chromosome are not suppressed by full methylation of their CpG islands. The differential methylation of the DXS255 CpG island provides the basis for a highly informative X inactivation analysis system.

A standardized format for handling data on plasmids, viruses and transposons: The PVT database format

The PVT format described here has been designed to store and retrieve genetic data on plasmids, viruses or transposons with special focus on their applications. Both naturally-occurring and engineered elements can be included in it. A variety of data can be accommodated in fields that are grouped in blocks: name and type of element, database administration, element administration, history, propagation, selection and host, biological properties, cloned insert and applications. The number of fields, now 157, can be expanded as required. Most properties can be described in simple logical fields. The format is organized to permit rapid searches and to facilitate communication between database and user; connection with culture and/or DNA collections is also envisaged and adequate fields for these tasks have been provided. The format allows cross-reference with that originated by the Microbial Information Network Europe (MINE) for computer storage and handling of bacterial or fungal strain data.

Characterization of a highly polymorphic region near the first exon of the human MAOA gene containing a GT dinucleotide and a novel VNTR motif

The genes encoding the A and B forms of the human monoamine oxidase enzymes (MAOA and MAOB) are localized at Xp11.23-Xp11.4. We report the characterization of a highly informative polymorphic region within a 2.9-kb cloned fragment containing the first exon of the MAOA gene. The polymorphic region consists of a GT microsatellite directly adjacent to an imperfectly duplicated novel 23-bp VNTR motif. DNA sequencing within and flanking the repeated segment allowed the design of specific amplification primers. In 56 unrelated females, 15 different alleles were identified with sizes ranging from 285 to 388 bp. The alleles differed in both the number of dinucleotide and the number of VNTR repeats, yielding a highly informative polymorphic marker locus with a calculated heterozygosity value of 75%.

The Ribosomal Database Project

The Ribosomal Database Project (RDP) complies ribosomal sequences and related data, and redistributes them in aligned and phylogenetically ordered form to its user community. It also offers various software packages for handling, analyzing and displaying sequences. In addition, the RDP offers (or will offer) certain analytic services. At present the project is in an intermediate stage of development.

Molecular biological databases — present and future

The importance of databases as a research tool in molecular biology is growing steadily, and a wide range of databases relevant to genome research is currently available. However, the design of current databases is inadequate for accurate representation and analysis of the results of large-scale genome mapping and sequencing projects. A new generation of databases is required to master the challenges of the future.

Molecular approaches for analyzing differential gene expression: differential cDNA library construction and screening

Complementary deoxyribonucleic acid (cDNA) libraries can be used as a means to isolate and identify cell-specific messenger ribonucleic acid (mRNA) sequences. The basic elements of cDNA library construction and screening are reviewed in the context of analyzing differentially expressed mRNAs. A brief overview of the recombinant DNA systems applied to cDNA library construction and the principles of screening cDNA libraries by plaque hybridization are provided. Methods for comparing mRNA populations by differential screening and by competition hybridization are discussed, and methods for constructing subtracted cDNA libraries, enriched in differentially expressed sequences, are presented. Also reviewed are the analysis of differentially expressed cDNAs by Southern and Northern hybridization, RNase protection, polymerase chain reaction, and sequencing.