Interstitial telomere repeats as markers of evolutionary changes in the mammalian karyotype: human chromosome 2

Telomer repeats represented by hexamer (TTAGGG)n at chromosome termini are required for correct function and chromosome stability. At the same time, interstitial telomer sequence (ITS) located far from the chromosome ends are known for several mammalian genomes, including the human genome. It is assumed that these repeats mark the points of fusion or other chromosome reconstructions of ancestors. Exact localization of all interstitial telomer sequences in the genome could greatly improve our understanding of the mechanism of karyotype evolution and species origin. We have developed a software for a search of interstitial telomer sequences in complete sequences of mammalian genomes. We have demonstrated the evolutionary significance of repeats by an example of human chromosome 2. The results and supplementary materials are available at the site of the Institute of Cytology and Genetics: http://www.bionet.nsc.ru/labs/theorylabmain/orlov/telomere/.     

Prediction of eukaryotic mRNA translational properties.

MOTIVATION: It is well known that eukaryotic mRNAs are translated at different levels depending on their sequence characteristics. Evaluation of mRNA translatability is of importance in prediction of the gene expression pattern by computer methods and to improve the recognition of mRNAs within cloned nucleotide sequences. It may also be used in biotechnological experiments to optimize the expression of foreign genes in transgenic organisms. RESULTS: The sets of 5' untranslated region characteristics, significantly different between mRNAs encoding abundant and scarce polypeptides, were determined for mammals, dicot plants and monocot plants, and collected in the LEADER_RNA database. Computer tools for the prediction of mRNA translatability are presented. AVAILABILITY: Programs for mRNA translatability prediction are available at http://wwwmgs.bionet.nsc. ru/programs/acts2/mo_mRNA.htm (for monocots), http://wwwmgs.bionet. nsc.ru/programs/acts2/di_mRNA.htm (for dicots) and http://wwwmgs. bionet.nsc.ru/programs/acts2/ma_mRNA.htm (for mammals). The LEADER_RNA database may be accessed at: http://wwwmgs.bionet.nsc. ru/systems/LeaderRNA/.     

Gene Network and Database for Genes of Wheat’s Resistance to Pathogenic Fungi

An overview describing a gene network that controls the formation of plant responses to diseases caused by pathogenic fungi (http://wwwmgs.bionet.nsc.ru/mgs/gnw/genenet//viewer/Plant%20fungus%20pathogen.html) is presented. The gene network represents the coordinated interactions of genes, proteins, and regulatory molecules, including integrated defense mechanisms that prevent the development of infection, localize the lesion, and minimize damage. The gene network was reconstructed on the basis of literature data, and the elements of the gene network were associated with the records of the PGR database (Pathogenesis-Related Genes, http://srs6.bionet.nsc.ru/srs6bin/cgi-bin/wgetz?-page+top+-newId), where information on plant genes resistant to pathogenic fungi is accumulated. Reconstruction of the gene network allows us to formalize, visualize, and systematize possible mechanisms for the response of plant cells to fungal infection, which may be useful for the planning of experiments and interpretation of experimental data in this field of science.     

Correlation of codon biases and potential secondary structures with mRNA translation efficiency in unicellular organisms

Gene expression is known to correlate with the degree of codon bias in many unicellular organisms. However, such a correlation is not observed in some organisms. It was demonstrated that inverted complementary repeats within coding DNA sequences (ORFs) should be considered for proper estimation of the translation efficiency because they can form secondary structures that obstruct ribosome movement. A program was developed for estimating the potential expression of ORFs in unicellular organisms on the basis of their genome sequences. The program computes the elongation efficiency index (EEI) and takes into account three key factors: codon bias, the average number of inverted complementary repeats, and the free energies of potential stem-loop structures formed by these repeats. The influence of these factors on translation was numerically estimated. Their optimal ratio was computed for each organism. EEIs of 384 unicellular organisms (351 bacteria, 28 archaea, and 5 eukaryotes) were computed using the annotated genomes available from GenBank. Five potential evolutionary strategies of translational optimization were determined in the organisms studied. A considerable difference in preferential translational strategies was observed between bacteria and archaea. Significant correlations between EEIs and gene expression levels were shown for two species (Saccharomyces cerevisiae and Helicobacter pylori), using the available microarray data. The method allows the numerical estimation of the translation efficiency of an ORF and optimization of the nucleotide composition of heterologous genes in specified unicellular organisms. The program is available at http://wwwmgs.bionet.nsc.ru/mgs/programs/eei-calculator.     

Prediction of contact numbers of amino acid residues using a neural network regression algorithm

The profile of contact numbers of amino acid residues in proteins contains important information about the protein structure and is connected with the accessibility of residues to solvent. Here we propose a method for predicting the profile of contact numbers of residues in protein from its amino acid sequence. The method is based on regression using a neural network algorithm. The algorithm predicts two types of profiles, namely, the total number of contacts and the number of close contacts with the neighbors in the chain. The Pearson coefficient of correlation between the actual and predicted values of total contact numbers amounted to 0.526–0.703. As for the number of close contacts, this coefficient was higher (0.662–0.743) for all the considered threshold contact distances (6, 8, 10, and 12 Å). The program for prediction of contact numbers CONNP is available at http://wwwmgs2.bionet.nsc.ru/reloaded.     

Prediction of contact numbers of amino acid residues using a neural network regression algorithm

The profile of contact numbers of amino acid residues in proteins contains important information about the protein structure and is connected with the accessibility of residues to solvent. Here we propose a method for predicting the profile of contact numbers of residues in protein from its amino acid sequence. The method is based on regression using a neural network algorithm. The algorithm predicts two types of profiles, namely, the total number of contacts and the number of close contacts with the neighbors in the chain. The Pearson coefficient of correlation between the actual and predicted values of total contact numbers amounted to 0.526–0.703. As for the number of close contacts, this coefficient was higher (0.662–0.743) for all the considered threshold contact distances (6, 8, 10, and 12 Å). The program for prediction of contact numbers CONNP is available at http://wwwmgs2.bionet.nsc.ru/reloaded.

     

Telomere sequence localization and karyotype evolution in higher plants

Data for chromosomal localization of theArabidopsis-type of telomeric sequence repeats (TTTAGGG)_n are compiled for 44 species belonging to 14 families of angiosperms, gymnosperms and bryophytes. For 23 species and seven families this is the first report. Species of all families, except theAlliaceae, revealed these sequences at their chromosome termini. This indicates thatArabidopsis-type telomeric repeats are highly conserved. It is inferred that they represent the basic telomere sequence of higher plant phyla. In theAlliaceae, a deviating sequence (and mechanism?) for the stabilization of chromosome termini has possibly evolved secondarily. Nine species revealed interstitial telomeric sequences in addition to the terminal ones, in three species (Vicia faba, Pinus elliottii, P. sylvestris) also at centromeric positions. Interstitial telomeric sequences may indicate karyotype reconstructions, in particular alterations of chromosome numbers by chromosome fusion — or inversions with one breakpoint within the terminal array of repeats. They may contribute to stabilization of chromosome breaks, especially centric fissions, and increase the frequency of meiotic and illegitimate recombination.     

Intrachromosomal location of the telomeric repeat (TTAGGG)n

Eukaryotic telomeres are specialized DNA-protein structures that are thought to ensure chromosomal stability and complete replication of the chromosome ends. All telomeres which have been studied consist of a tandem array of G-rich repeats which seem to be sufficient for telomere function. Originally, the human telomeric repeat (TTAGGG)_n was assumed to be exclusively located at the very end of all human chromosomes. More recent evidence, however, suggests an extension into proterminal regions. Very little is known about the interstitial distribution of telomeric repeats. Here we present evidence for the presence of (TTAGGG)_n repeats in internal loci on the long and short arms of different human chromosomes. In addition, we studied the genomic organization of these repeats in more detail and discuss possible functions of interstitial telomeric repeats in the human genome.     

A systems approach to morphogenesis in Arabidopsis thaliana: I. AGNS database

In systems biology, study of a complex and multicomponent system, such as morphogenesis, comprises accumulation of data on morphogenetic processes in databases, classification and logical analysis of this information, and computer simulation of the processes in question using the data accumulated and the results of their analysis. This paper describes realization of the first steps in a systems study of morphogenesis (annotating research papers, compiling information in a database, data systematization, and their logical analysis) by the example of Arabidopsis thaliana, a model object in plant molecular biology. The database AGNS (Arabidopsis GeneNet Supplementary; http://wwwmgs.bionet.nsc.ru/agns) contains the experimentally confirmed information from published papers on specific features of gene expression and phenotypes of wild-type, mutant, and transgenic A. thaliana plants. AGNS queries and logical data analysis with the aid of specially developed software makes it possible to model various morphogenetic processes from gene expression to functioning of gene networks and their contribution to the development of certain traits.

     

SELEX_DB: an activated database on selected randomized DNA/RNA sequences addressed to genomic sequence annotation

SELEX_DB is a novel curated database on selected randomized DNA/RNA sequences designed for accumulation of experimental data on functional site sequences obtained by using SELEX and SELEX-like technologies from the pools of random sequences. This database also contains the programs for DNA/RNA functional site recognition within arbitrary nucleotide sequences. The first release of SELEX_DB has been installed under SRS and is available through the WWW at http://wwwmgs.bionet.nsc.ru/mgs/systems/selex/     

Nucleosome formation potential of exons, introns, and Alu repeats.

A program for constructing nucleosome formation potential profile was applied for investigation of exons, introns, and repetitive sequences. The program is available at http://wwwmgs.bionet.nsc.ru/mgs/programs/recon/. We have demonstrated that introns and repetitive sequences exhibit higher nucleosome formation potentials than exons. This fact may be explained by functional saturation of exons with genetic code, hindering the localization of efficient nucleosome positioning sites.     

A systems approach to morphogenesis in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>: I. AGNS database

In systems biology, study of a complex and multicomponent system, such as morphogenesis, comprises accumulation of data on morphogenetic processes in databases, classification and logical analysis of this information, and computer simulation of the processes in question using the data accumulated and the results of their analysis. This paper describes realization of the first steps in a systems study of morphogenesis (annotating research papers, compiling information in a database, data systematization, and their logical analysis) by the example of Arabidopsis thaliana, a model object in plant molecular biology. The database AGNS (Arabidopsis GeneNet Supplementary; http://wwwmgs.bionet.nsc.ru/agns) contains the experimentally confirmed information from published papers on specific features of gene expression and phenotypes of wild-type, mutant, and transgenic A. thaliana plants. AGNS queries and logical data analysis with the aid of specially developed software makes it possible to model various morphogenetic processes from gene expression to functioning of gene networks and their contribution to the development of certain traits.     

Investigating extended regulatory regions of genomic DNA sequences.

MOTIVATION: Despite the growing volume of data on primary nucleotide sequences, the regulatory regions remain a major puzzle with regard to their function. Numerous recognising programs considering a diversity of properties of regulatory regions have been developed. The system proposed here allows the specific contextual, conformational and physico-chemical properties to be revealed based on analysis of extended DNA regions. RESULTS: The Internet-accessible computer system RegScan, designed to analyse the extended regulatory regions of eukaryotic genes, has been developed. The computer system comprises the following software: (i) programs for classification dividing a set of promoters into TATA-containing and TATA-less promoters and promoters with and without CpG islands; (ii) programs for constructing (a) nucleotide frequency profiles, (b) sequence complexity profiles and (c) profiles of conformational and physico-chemical properties; (iii) the program for constructing the sets of degenerate oligonucleotide motifs of a specified length; and (iv) the program searching for and visualising repeats in nucleotide sequences. The system has allowed us to demonstrate the following characteristic patterns of vertebrate promoter regions: the TATA box region is flanked by regions with an increased G+C content and increased bending stiffness, the TATA box content is asymmetric and promoter regions are saturated with both direct and inverted repeats. AVAILABILITY: The computer system RegScan is available via the Internet at http://www.mgs.bionet.nsc. ru/Systems/RegScan, http://www.cbil.upenn.edu/mgs/systems/r egscan/.