Kenji Satou

Finding microRNA regulatory modules in human genome using rule induction

Background:MicroRNAs (miRNAs) are a class of small non-coding RNA molecules (20–24 nt), which are believed to participate in repression of gene expression. They play important roles in several biological processes (e.g. cell death and cell growth). Both experimental and computational approaches have been used to determine the function of miRNAs in cellular processes. Most efforts have concentrated on identification of miRNAs and their target genes. However, understanding the regulatory mechanism of miRNAs in the gene regulatory network is also essential to the discovery of functions of miRNAs in complex cellular systems. To understand the regulatory mechanism of miRNAs in complex cellular systems, we need to identify the functional modules involved in complex interactions between miRNAs and their target genes.Results:We propose a rule-based learning method to identify groups of miRNAs and target genes that are believed to participate cooperatively in the post-transcriptional gene regulation, so-called miRNA regulatory modules (MRMs). Applying our method to human genes and miRNAs, we found 79 MRMs. The MRMs are produced from multiple information sources, including miRNA-target binding information, gene expression and miRNA expression profiles. Analysis of two first MRMs shows that these MRMs consist of highly-related miRNAs and their target genes with respect to biological processes.Conclusion:The MRMs found by our method have high correlation in expression patterns of miRNAs as well as mRNAs. The mRNAs included in the same module shared similar biological functions, indicating the ability of our method to detect functionality-related genes. Moreover, review of the literature reveals that miRNAs in a module are involved in several types of human cancer.

Inflammation-induced repression of tumor suppressor miR-7 in gastric tumor cells

Inflammation has an important role in cancer development through various mechanisms. It has been shown that dysregulation of microRNAs (miRNAs) that function as oncogenes or tumor suppressors contributes to tumorigenesis. However, the relationship between inflammation and cancer-related miRNA expression in tumorigenesis has not yet been fully understood. Using K19-C2mE and Gan mouse models that develop gastritis and gastritis-associated tumors, respectively, we found that 21 miRNAs were upregulated, and that 29 miRNAs were downregulated in gastric tumors in an inflammation-dependent manner. Among these miRNAs, the expression of miR-7, a possible tumor suppressor, significantly decreased in both gastritis and gastric tumors. Moreover, the expression of miR-7 in human gastric cancer was inversely correlated with the levels of interleukin-1β and tumor necrosis factor-α, suggesting that miR-7 downregulation is related to the severity of inflammatory responses. In the normal mouse stomach, miR-7 expression was at a basal level in undifferentiated gastric epithelial cells, and was induced during differentiation. Moreover, transfection of a miR-7 precursor into gastric cancer cells suppressed cell proliferation and soft agar colony formation. These results suggest that suppression of miR-7 expression is important for maintaining the undifferentiated status of gastric epithelial cells, and thus contributes to gastric tumorigenesis. Although epigenetic changes were not found in the CpG islands around miR-7-1 of gastritis and gastric tumor cells, we found that activated macrophage-derived small molecule(s) (<3 kDa) are responsible for miR-7 repression in gastric cancer cells. Furthermore, the miR-7 expression level significantly decreased in the inflamed gastric mucosa of Helicobacter-infected mice, whereas it increased in the stomach of germfree K19-C2mE and Gan mice wherein inflammatory responses were suppressed. Taken together, these results indicate that downregulation of tumor suppressor miR-7 is a novel mechanism by which the inflammatory response promotes gastric tumorigenesis.

Inferring microRNA and transcription factor regulatory networks in heterogeneous data

BackgroundTranscription factors (TFs) and microRNAs (miRNAs) are primary metazoan gene regulators. Regulatory mechanisms of the two main regulators are of great interest to biologists and may provide insights into the causes of diseases. However, the interplay between miRNAs and TFs in a regulatory network still remains unearthed. Currently, it is very difficult to study the regulatory mechanisms that involve both miRNAs and TFs in a biological lab. Even at data level, a network involving miRNAs, TFs and genes will be too complicated to achieve. Previous research has been mostly directed at inferring either miRNA or TF regulatory networks from data. However, networks involving a single type of regulator may not fully reveal the complex gene regulatory mechanisms, for instance, the way in which a TF indirectly regulates a gene via a miRNA.ResultsWe propose a framework to learn from heterogeneous data the three-component regulatory networks, with the presence of miRNAs, TFs, and mRNAs. This method firstly utilises Bayesian network structure learning to construct a regulatory network from multiple sources of data: gene expression profiles of miRNAs, TFs and mRNAs, target information based on sequence data, and sample categories. Then, in order to produce more meaningful results for further biological experimentation and research, the method searches the learnt network to identify the interplay between miRNAs and TFs and applies a network motif finding algorithm to further infer the network.We apply the proposed framework to the data sets of epithelial-to-mesenchymal transition (EMT). The results elucidate the complex gene regulatory mechanism for EMT which involves both TFs and miRNAs. Several discovered interactions and molecular functions have been confirmed by literature. In addition, many other discovered interactions and bio-markers are of high statistical significance and thus can be good candidates for validation by experiments. Moreover, the results generated by our method are compact, involving a small number of interactions which have been proved highly relevant to EMT.ConclusionsWe have designed a framework to infer gene regulatory networks involving both TFs and miRNAs from multiple sources of data, including gene expression data, target information, and sample categories. Results on the EMT data sets have shown that the proposed approach is able to produce compact and meaningful gene regulatory networks that are highly relevant to the biological conditions of the data sets. This framework has the potential for application to other heterogeneous datasets to reveal the complex gene regulatory relationships.

Phylogenetic reconstruction from non-genomic data

MOTIVATION Recent results related to horizontal gene transfer suggest that phylogenetic reconstruction cannot be determined conclusively from sequence data, resulting in a shift from approaches based on polymorphism information in DNA or protein sequence to studies aimed at understanding the evolution of complete biological processes. The increasing amount of available information on metabolic pathways for several species makes it of greater relevance to understand the similarities and differences among such pathways. These similarities can then be used to infer phylogenetic trees not based exclusively in sequence data, therefore avoiding the previously mentioned problems. RESULTS In this article, we present a method to assess the structural similarity of metabolic pathways for several organisms. Our algorithms work by using one of the three possible enzyme similarity measures (hierarchical, information content, gene ontology), and one of the two clustering methods (neighbor-joining, unweighted pair group method with arithmetic mean), to produce a phylogenetic tree both in Newick and graphic format. The web server implementing our algorithms is optimized to answer queries in linear time. AVAILABILITY The software is available for free public use on a web server, at the address http://www.jaist.ac.jp/~clemente/cgi-bin/phylo.pl. It is available on demand in source code form for research use to educational institutions, non-profit research institutes, government research laboratories and individuals, for non-exclusive use, without the right of the licensee to further redistribute the source code.

Comparative analysis of protein thermostability: Differences in amino acid content and substitution at the surfaces and in the core regions of thermophilic and mesophilic proteins

Abstract In order to investigate the factors responsible for protein thermostability, we performed a comparative analysis. For this study, we prepared a new dataset composed of 47 homologous pairs of thermophilic and mesophilic proteins. It is he largest comparative study dataset ever presented. The frequency and substitution preference of each amino acid type in the dataset were analyzed.Twokinds of residual structural states were considered, i.e. surface (solvent-exposed) and core (buried) regions. On the surface of thermophilic proteins, higher frequencies were observed for Arg, Glu, and Tyr. Analysis of substitution preference also suggests that these often appear by replacement of other amino acid types. The results indicate that Arg, Glu, and Tyr are suitable for location on the surface of thermophilic proteins. On the other hand, at the core of thermophilic proteins, Ala is often appeared. In addition, our t-test analysis provides the first quantitative information about trends in the frequencies and substitution preferences for Cys, Gln, Met, and Ser. The results indicate that Gln and Met on the surface and Cys and Ser in the core are disadvantageous for protein thermostability.

Support vector machines for prediction and analysis of beta and gamma-turns in proteins.

Tight turns have long been recognized as one of the three important features of proteins, together with alpha-helix and beta-sheet. Tight turns play an important role in globular proteins from both the structural and functional points of view. More than 90% tight turns are beta-turns and most of the rest are gamma-turns. Analysis and prediction of beta-turns and gamma-turns is very useful for design of new molecules such as drugs, pesticides, and antigens. In this paper we investigated two aspects of applying support vector machine (SVM), a promising machine learning method for bioinformatics, to prediction and analysis of beta-turns and gamma-turns. First, we developed two SVM-based methods, called BTSVM and GTSVM, which predict beta-turns and gamma-turns in a protein from its sequence. When compared with other methods, BTSVM has a superior performance and GTSVM is competitive. Second, we used SVMs with a linear kernel to estimate the support of amino acids for the formation of beta-turns and gamma-turns depending on their position in a protein. Our analysis results are more comprehensive and easier to use than the previous results in designing turns in proteins.

Computational discovery of miR-TF regulatory modules in human genome.

MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression at the post-transcriptional level. They play an important role in several biological processes such as cell development and differentiation. Similar to transcription factors (TFs), miRNAs regulate gene expression in a combinatorial fashion, i.e., an individual miRNA can regulate multiple genes, and an individual gene can be regulated by multiple miRNAs. The functions of TFs in biological regulatory networks have been well explored. And, recently, a few studies have explored miRNA functions in the context of gene regulation networks. However, how TFs and miRNAs function together in the gene regulatory network has not yet been examined. In this paper, we propose a new computational method to discover the gene regulatory modules that consist of miRNAs, TFs, and genes regulated by them. We analyzed the regulatory associations among the sets of predicted miRNAs and sets of TFs on the sets of genes regulated by them in the human genome. We found 182 gene regulatory modules of combinatorial regulation by miRNAs and TFs (miR-TF modules). By validating these modules with the Gene Ontology (GO) and the literature, it was found that our method allows us to detect functionally-correlated gene regulatory modules involved in specific biological processes. Moreover, our miR-TF modules provide a global view of coordinated regulation of target genes by miRNAs and TFs.

Using inductive logic programming for predicting protein-protein interactions from multiple genomic data

Protein-protein interactions play an important role in many fundamental biological processes. Computational approaches for predicting protein-protein interactions are essential to infer the functions of unknown proteins, and to validate the results obtained of experimental methods on protein-protein interactions. We have developed an approach using Inductive Logic Programming (ILP) for protein-protein interaction prediction by exploiting multiple genomic data including protein-protein interaction data, SWISS-PROT database, cell cycle expression data, Gene Ontology, and InterPro database. The proposed approach demonstrates a promising result in terms of obtaining high sensitivity/specificity and comprehensible rules that are useful for predicting novel protein-protein interactions. We have also applied our method to a number of protein-protein interaction data, demonstrating an improvement on the expression profile reliability (EPR) index.

Computational discovery of transcriptional regulatory rules

MOTIVATION Even in a simple organism like yeast Saccharomyces cerevisiae, transcription is an extremely complex process. The expression of sets of genes can be turned on or off by the binding of specific transcription factors to the promoter regions of genes. Experimental and computational approaches have been proposed to establish mappings of DNA-binding locations of transcription factors. However, although location data obtained from experimental methods are noisy owing to imperfections in the measuring methods, computational approaches suffer from over-prediction problems owing to the short length of the sequence motifs bound by the transcription factors. Also, these interactions are usually environment-dependent: many regulators only bind to the promoter region of genes under specific environmental conditions. Even more, the presence of regulators at a promoter region indicates binding but not necessarily function: the regulator may act positively, negatively or not act at all. Therefore, identifying true and functional interactions between transcription factors and genes in specific environment conditions and describing the relationship between them are still open problems. RESULTS We developed a method that combines expression data with genomic location information to discover (1) relevant transcription factors from the set of potential transcription factors of a target gene; and (2) the relationship between the expression behavior of a target gene and that of its relevant transcription factors. Our method is based on rule induction, a machine learning technique that can efficiently deal with noisy domains. When applied to genomic location data with a confidence criterion relaxed to P-value = 0.005, and three different expression datasets of yeast S.cerevisiae, we obtained a set of regulatory rules describing the relationship between the expression behavior of a specific target gene and that of its relevant transcription factors. The resulting rules provide strong evidence of true positive gene-regulator interactions, as well as of protein-protein interactions that could serve to identify transcription complexes. AVAILABILITY Supplementary files are available from http://www.jaist.ac.jp/~h-pham/regulatory-rules

The superstructure toward open bioinformatics grid

The grid design strongly depends on not only a network infrastructure but also a superstructure, that is, a social structure of virtual organizations where people trust each other, share resources and work together. Open Bioinformatics Grid (OBIGrid) is a grid aimed at building a cooperative bioinformatics environment for computer sicentists and biologists. In October 2003, OBIGrid consisted of 293 nodes with 492 CPUs provided by 27 sites at universities, laboratories and other enterprises, connected by a virtual private network over the Internet. So many organizations have participated because OBIGrid has been conscious of constructing a superstructure on a grid as well as a grid infrastructure. For the benefit of OBIGrid participants, we have developed a series of life science application services: an open bioinformatics environment (OBIEnv), a scalable genome database (OBISgd), a genome annotation system (OBITco), a biochemical network simulator (OBIYagns), and to name a few.