Mika Hirakawa

KEGG for linking genomes to life and the environment.

KEGG (http://www.genome.jp/kegg/) is a database of biological systems that integrates genomic, chemical and systemic functional information. KEGG provides a reference knowledge base for linking genomes to life through the process of PATHWAY mapping, which is to map, for example, a genomic or transcriptomic content of genes to KEGG reference pathways to infer systemic behaviors of the cell or the organism. In addition, KEGG provides a reference knowledge base for linking genomes to the environment, such as for the analysis of drug-target relationships, through the process of BRITE mapping. KEGG BRITE is an ontology database representing functional hierarchies of various biological objects, including molecules, cells, organisms, diseases and drugs, as well as relationships among them. KEGG PATHWAY is now supplemented with a new global map of metabolic pathways, which is essentially a combined map of about 120 existing pathway maps. In addition, smaller pathway modules are defined and stored in KEGG MODULE that also contains other functional units and complexes. The KEGG resource is being expanded to suit the needs for practical applications. KEGG DRUG contains all approved drugs in the US and Japan, and KEGG DISEASE is a new database linking disease genes, pathways, drugs and diagnostic markers.

From genomics to chemical genomics: new developments in KEGG

The increasing amount of genomic and molecular information is the basis for understanding higher-order biological systems, such as the cell and the organism, and their interactions with the environment, as well as for medical, industrial and other practical applications. The KEGG resource () provides a reference knowledge base for linking genomes to biological systems, categorized as building blocks in the genomic space (KEGG GENES) and the chemical space (KEGG LIGAND), and wiring diagrams of interaction networks and reaction networks (KEGG PATHWAY). A fourth component, KEGG BRITE, has been formally added to the KEGG suite of databases. This reflects our attempt to computerize functional interpretations as part of the pathway reconstruction process based on the hierarchically structured knowledge about the genomic, chemical and network spaces. In accordance with the new chemical genomics initiatives, the scope of KEGG LIGAND has been significantly expanded to cover both endogenous and exogenous molecules. Specifically, RPAIR contains curated chemical structure transformation patterns extracted from known enzymatic reactions, which would enable analysis of genome-environment interactions, such as the prediction of new reactions and new enzyme genes that would degrade new environmental compounds. Additionally, drug information is now stored separately and linked to new KEGG DRUG structure maps.

KEGG for representation and analysis of molecular networks involving diseases and drugs.

Most human diseases are complex multi-factorial diseases resulting from the combination of various genetic and environmental factors. In the KEGG database resource (http://www.genome.jp/kegg/), diseases are viewed as perturbed states of the molecular system, and drugs as perturbants to the molecular system. Disease information is computerized in two forms: pathway maps and gene/molecule lists. The KEGG PATHWAY database contains pathway maps for the molecular systems in both normal and perturbed states. In the KEGG DISEASE database, each disease is represented by a list of known disease genes, any known environmental factors at the molecular level, diagnostic markers and therapeutic drugs, which may reflect the underlying molecular system. The KEGG DRUG database contains chemical structures and/or chemical components of all drugs in Japan, including crude drugs and TCM (Traditional Chinese Medicine) formulas, and drugs in the USA and Europe. This database also captures knowledge about two types of molecular networks: the interaction network with target molecules, metabolizing enzymes, other drugs, etc. and the chemical structure transformation network in the history of drug development. The new disease/drug information resource named KEGG MEDICUS can be used as a reference knowledge base for computational analysis of molecular networks, especially, by integrating large-scale experimental datasets.

JSNP: a database of common gene variations in the Japanese population.

JSNP is a repository of Japanese Single Nucleotide Polymorphism (SNP) data, begun in 2000 and developed through the Prime Ministers Millennium Project. The aim of this undertaking is to identify and collate up to 150 000 SNPs from the Japanese population, located in genes or in adjacent regions that might influence the coding sequence of the genes. The project has been carried out by a collaboration between the Human Genome Center (HGC) in the Institute of Medical Science (IMS) at the University of Tokyo and the Japan Science and Technology Corporation (JST). JSNP serves as both a storage site for the Japanese SNPs obtained from the ongoing project and as a facility for public dissemination to allow researchers access to high quality SNP data. A primary motivation of the project is the construction of a basic data set to identify relationships between polymorphisms and common diseases or the reaction to drugs. As such, emphasis has been placed on the identification of SNPs that lie in candidate regions which may affect phenotype but which would not necessarily directly cause disease. Unrestricted access to JSNP and any associated files is available at http://snp.ims.u-tokyo.ac.jp/.

The KEGG Databases and Tools Facilitating Omics Analysis: Latest Developments Involving Human Diseases and Pharmaceuticals

In this chapter, we demonstrate the usability of the KEGG (Kyoto encyclopedia of genes and genomes) databases and tools, especially focusing on the visualization of the omics data. The desktop application KegArray and many Web-based tools are tightly integrated with the KEGG knowledgebase, which helps visualize and interpret large amount of data derived from high-throughput measurement techniques including microarray, metagenome, and metabolome analyses. Recently developed resources for human disease, drug, and plant research are also mentioned.

A mammalian conserved element derived from SINE displays enhancer properties recapitulating Satb2 expression in early-born callosal projection neurons.

Short interspersed repetitive elements (SINEs) are highly repeated sequences that account for a significant proportion of many eukaryotic genomes and are usually considered “junk DNA”. However, we previously discovered that many AmnSINE1 loci are evolutionarily conserved across mammalian genomes, suggesting that they may have acquired significant functions involved in controlling mammalian-specific traits. Notably, we identified the AS021 SINE locus, located 390 kbp upstream of Satb2. Using transgenic mice, we showed that this SINE displays specific enhancer activity in the developing cerebral cortex. The transcription factor Satb2 is expressed by cortical neurons extending axons through the corpus callosum and is a determinant of callosal versus subcortical projection. Mouse mutants reveal a crucial function for Sabt2 in corpus callosum formation. In this study, we compared the enhancer activity of the AS021 locus with Satb2 expression during telencephalic development in the mouse. First, we showed that the AS021 enhancer is specifically activated in early-born Satb2+ neurons. Second, we demonstrated that the activity of the AS021 enhancer recapitulates the expression of Satb2 at later embryonic and postnatal stages in deep-layer but not superficial-layer neurons, suggesting the possibility that the expression of Satb2 in these two subpopulations of cortical neurons is under genetically distinct transcriptional control. Third, we showed that the AS021 enhancer is activated in neurons projecting through the corpus callosum, as described for Satb2+ neurons. Notably, AS021 drives specific expression in axons crossing through the ventral (TAG1−/NPY+) portion of the corpus callosum, confirming that it is active in a subpopulation of callosal neurons. These data suggest that exaptation of the AS021 SINE locus might be involved in enhancement of Satb2 expression, leading to the establishment of interhemispheric communication via the corpus callosum, a eutherian-specific brain structure.

A SINE-Derived Element Constitutes a Unique Modular Enhancer for Mammalian Diencephalic Fgf8

Transposable elements, including short interspersed repetitive elements (SINEs), comprise nearly half the mammalian genome. Moreover, they are a major source of conserved non-coding elements (CNEs), which play important functional roles in regulating development-related genes, such as enhancing and silencing, serving for the diversification of morphological and physiological features among species. We previously reported a novel SINE family, AmnSINE1, as part of mammalian-specific CNEs. One AmnSINE1 locus, named AS071, showed an enhancer property in the developing mouse diencephalon. Indeed, AS071 appears to recapitulate the expression of diencephalic fibroblast growth factor 8 (Fgf8). Here we established three independent lines of AS071-transgenic mice and performed detailed expression profiling of AS071-enhanced lacZ in comparison with that of Fgf8 across embryonic stages. We demonstrate that AS071 is a distal enhancer that directs Fgf8 expression in the developing diencephalon. Furthermore, enhancer assays with constructs encoding partially deleted AS071 sequence revealed a unique modular organization in which AS071 contains at least three functionally distinct sub-elements that cooperatively direct the enhancer activity in three diencephalic domains, namely the dorsal midline and the lateral wall of the diencephalon, and the ventral midline of the hypothalamus. Interestingly, the AmnSINE1-derived sub-element was found to specify the enhancer activity to the ventral midline of the hypothalamus. To our knowledge, this is the first discovery of an enhancer element that could be separated into respective sub-elements that determine regional specificity and/or the core enhancing activity. These results potentiate our understanding of the evolution of retroposon-derived cis-regulatory elements as well as the basis for future studies of the molecular mechanism underlying the determination of domain-specificity of an enhancer.

Amino Acid Substitution in HCV Core/NS5A Region and Genetic Variation Near IL28B Gene Affect Treatment Efficacy to Interferon plus Ribavirin Combination Therapy

Objective: To evaluate predictive factors of treatment efficacy to interferon (IFN)/ribavirin in patients infected with HCV genotype 1b (HCV-1b). Methods: This study investigated pretreatment predictors, including viral- (aa substitutions in core aa 70/91 and NS5A-ISDR/IRRDR) and host-related factors (genetic variation near IL28B gene), to 48-week IFN/ribavirin in 490 Japanese adults infected with HCV-1b. Results: The proportion of patients who showed end-of-treatment response (ETR), sustained virological response (SVR), and SVR after ETR was 76, 54, and 76%, respectively. There was a significant positive correlation between the number of aa substitutions in ISDR and those in IRRDR. Concerning the substitution of core aa 91, the number of aa substitutions in ISDR/IRRDR of patients with Leu91 was significantly higher than that of patients with Met91. Furthermore, levels of viremia were influenced by aa substitutions in core aa 91 and ISDR/IRRDR. By multivariate analysis, rs8099917 genotype was an important predictor of ETR and SVR. With regard to viral factors, core aa 70/91 was an important predictor of ETR, and SVR after ETR. ISDR was an important predictor of SVR, and SVR after ETR. Conclusion: aa substitution in core/NS5A region and genetic variation near IL28B were important predictors of treatment efficacy to IFN/ribavirin.

Characterization and evolutionary landscape of AmnSINE1 in Amniota genomes

Discovery of a large number of conserved non-coding elements (CNEs) in vertebrate genomes provides a cornerstone to elucidate molecular mechanisms of macroevolution. Extensive comparative genomics has proven that transposons such as short interspersed elements (SINEs) were an important source of CNEs. We recently characterized AmnSINE1, a SINE family in Amniota genomes, some of which are present in CNEs, and demonstrated that two AmnSINE1 loci play an important role in mammalian-specific brain development by functioning as an enhancer (Sasaki et al. Proc. Natl. Acad. Sci. USA 2008). To get more information about AmnSINE1s, we here performed a multi-species search for AmnSINE1, and revealed the distribution and evolutionary history of these SINEs in amniote genomes. The number of AmnSINE1 regions in amniotes ranged from 160 to 1200; the number in the eutherians were under 500 and the largest was that in chicken. Phylogenetic analysis established that each AmnSINE1 locus has evolved uniquely, primarily since the divergence of mammals from reptiles. These results support the notion that AmnSINE1s were amplified as an ancient retroposon in a common ancestor of Amniota and subsequently have survived for 300 Myr because of functions acquired by mutation-coupled exaptation prior mammalian radiation. On the basis of sequence homology and conserved synteny, we detected the orthologs of AmnSINE1 for candidates of further enhancer analysis, which are more conserved than two loci that were shown to have been involved in mammalian brain development. The present work provides a comprehensive data set to test the role of AmnSINE1s, many of which were exapted and contributed to mammalian macroevolution.

Coordinately Co-opted Multiple Transposable Elements Constitute an Enhancer for wnt5a Expression in the Mammalian Secondary Palate

Acquisition of cis-regulatory elements is a major driving force of evolution, and there are several examples of developmental enhancers derived from transposable elements (TEs). However, it remains unclear whether one enhancer element could have been produced via cooperation among multiple, yet distinct, TEs during evolution. Here we show that an evolutionarily conserved genomic region named AS3_9 comprises three TEs (AmnSINE1, X6b_DNA and MER117), inserted side-by-side, and functions as a distal enhancer for wnt5a expression during morphogenesis of the mammalian secondary palate. Functional analysis of each TE revealed step-by-step retroposition/transposition and co-option together with acquisition of a binding site for Msx1 for its full enhancer function during mammalian evolution. The present study provides a new perspective suggesting that a huge variety of TEs, in combination, could have accelerated the diversity of cis-regulatory elements involved in morphological evolution.