Functional linear models for region-based association analysis

Regional association analysis is one of the most powerful tools for gene mapping because instead analysis of individual variants it simultaneously considers all variants in the region. Recent development of the models for regional association analysis involves functional data analysis approach. In the framework of this approach, genotypes of variants within region as well as their effects are described by continuous functions. Such approach allows us to use information about both linkage and linkage disequilibrium and reduce the influence of noise and/or observation errors. Here we define a functional linear mixed model to test association on independent and structured samples. We demonstrate how to test fixed and random effects of a set of genetic variants in the region on quantitative trait. Estimation of statistical properties of new methods shows that type I errors are in accordance with declared values and power is high especially for models with fixed effects of genotypes. We suppose that new functional regression linear models facilitate identification of rare genetic variants controlling complex human and animal traits. New methods are implemented in computer software FREGAT which is available for free download at http://mga.bionet.nsc.ru/soft/FREGAT/.     

An Algorithm of Step-by-step Pedigree Drawing

An algorithm for drawing large, complex pedigrees containing inbred loops and multiple-mate families is presented. The algorithm is based on a step-by-step approach to imaging, when the researcher determines the direction of further extension of the scheme. The algorithm is implemented as the PedigreeQuery software package written in Java. The software has a convenient graphical interface. The software package permits constructing not only whole pedigrees, but also their fragments that are particularly interesting for research. It also allows for adding new information on the phenotypes and genotypes of pedigree members. PedigreeQuery is distributed free of charge; it is available at http://mga.bionet.msc.ru/PedigreeQuery/PedigreeQuery.html and ftp://mga.bionet.nsc.ru/PedigreeQuery/.     

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.     

Breaking loops in large complex pedigrees

For pedigrees with multiple loops, exact likelihoods could not be computed in an acceptable time frame and thus, approximate methods are used. Some of these methods are based on breaking loops and approximations of complex pedigree likelihoods using the exact likelihood of the corresponding zero-loop pedigree. Due to ignoring loops, this method results in a loss of genetic information and a decrease in the power to detect linkage. To minimize this loss, an optimal set of loop breakers has to be selected. In this paper, we present a graph theory based algorithm for automatic selection of an optimal set of loop breakers. We propose using a total relationship between measured pedigree members as a proxy to power. To minimize the loss of genetic information, we suggest selection of such breakers whose duplication in a pedigree would be accompanied by a minimal loss of total relationship between measured pedigree members. We show that our algorithm compares favorably with other existing loop-breaker selection algorithms in terms of conservation of genetic information, statistical power and CPU time of subsequent linkage analysis. We implemented our method in a software package LOOP_EDGE, which is available at http://mga.bionet.nsc.ru/nlru/.     

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/.     

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.     

Interstitial telomeric repeats as markers of evolutionary changes in the mammalian karyotype: Human chromosome 2

The presence of conserved telomeric repeats represented by the hexamer (TTAGGG)_n at the chromosomal termini is necessary for the correct functioning and stability of chromosomes. A number of the genomes of mammals, including human, are known to contain interstitial telomeric sequences located far from the chromosomal termini. It is assumed that these repeats mark the regions of fusions or other rear-rangements of ancestral chromosomes. Exact localization of all interstitial telomeric sequences in the genome could significantly advance the understanding of the mechanisms of karyotype evolution and speciation. In this context, software was developed to search for degenerate interstitial telomeric repeats in complete sequences of mammalian chromosomes. The evolutionary significance of such repeats was demonstrated by the example of human chromosome 2. The results are available at http://www.bionet.nsc.ru/labs/theorylabmain/orlov/telomere/.     

An algorithm of step-by-step pedigree drawing

An algorithm for drawing large, complex pedigrees containing inbred loops and multiple-mate families is presented. The algorithm is based on a step-by-step approach to imaging, when the researcher determines the direction of further extension of the scheme. The algorithm is implemented as the PedigreeQuery software package written in Java. The software has a convenient graphical interface. The software package permits constructing not only whole pedigrees, but also their fragments that are particularly interesting for research. It also allows for adding new information on the phenotypes and genotypes of pedigree members. PedigreeQuery is distributed free of charge; it is available at http://mga.bionet.msc.ru/PedigreeQuery/PedigreeQuery.html and ftp://mga.bionet.msc.ru/PedigreeQuery/.     

Modeling evolution of spatially distributed bacterial communities: a simulation with the haploid evolutionary constructor

Background

Multiscale approaches for integrating submodels of various levels of biological organization into a single model became the major tool of systems biology. In this paper, we have constructed and simulated a set of multiscale models of spatially distributed microbial communities and study an influence of unevenly distributed environmental factors on the genetic diversity and evolution of the community members.

Results

Haploid Evolutionary Constructor software http://evol-constructor.bionet.nsc.ru/ was expanded by adding the tool for the spatial modeling of a microbial community (1D, 2D and 3D versions). A set of the models of spatially distributed communities was built to demonstrate that the spatial distribution of cells affects both intensity of selection and evolution rate.

Conclusion

In spatially heterogeneous communities, the change in the direction of the environmental flow might be reflected in local irregular population dynamics, while the genetic structure of populations (frequencies of the alleles) remains stable. Furthermore, in spatially heterogeneous communities, the chemotaxis might dramatically affect the evolution of community members.

     

Integrated computer system for regulating eukaryotic gene expression

The review describes several modules of the GeneExpress integrated computer system concerning the regulation of gene expression in eukaryotes. Approaches to the presentation of experimental data in databases are considered. The employment of GeneExpress in computer analysis and modeling of the organization and function of genetic systems is illustrated with examples. GeneExpress is available at http://wwwmgs.bionet.nsc.ru/mgs/gnw/.     

Estimation of Allele Frequencies in Ethnically Heterogeneous Populations

A method for reconstructing allele frequencies characteristic of an original ethnically homogeneous population before the start of migration processes is described. Information on both the ethnic group studied and offspring of interethnic marriages is used to estimate the allele frequencies. This makes it possible to increase the informativeness of the sample, which, in the case of ethnic heterogeneity, depends not only on allele frequencies and the total sample size, but also on the ethnic structure of the sample. The problem of estimating allele frequency in an ethnically heterogeneous sample has been solved analytically for diallelic loci. It has been demonstrated that, if offspring of interethnic marriages with the same degree of outbreeding is added to a sample of the ethnic group studied, the sample informativeness does not change. To utilize the information contained in the phenotypes of the offspring of interethnic marriages, representatives of the population from which migration occurs should be included into the sample. The size of the sample ensuring the preassigned accuracy of estimation is minimized at a certain ratio between the numbers of the offspring of interethnic marriages and the “immigrants.” To analyze polyallelic loci, a software package has been developed that allows estimating allele frequencies, determining the errors of these estimates, and planning the sample ensuring the preassigned accuracy of estimation. The package is available free at http://mga.bionet.nsc.ru/PopMixed/PopMixed.html.__________Translated from Genetika, Vol. 41, No. 7, 2005, pp. 990–996.Original Russian Text Copyright © 2005 by Axenovich, Kirichenko.     

An Integrated Computer System for Studying the Regulation of Eukaryotic Gene Expression

The review describes several modules of the GeneExpress integrated computer system concerning the regulation of gene expression in eukaryotes. Approaches to the presentation of experimental data in databases are considered. The employment of GeneExpress in computer analysis and modeling of the organization and function of genetic systems is illustrated with examples. GeneExpress is available at http://wwwmgs.bionet.nsc.ru/mgs/gnw/.     

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.     

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.