Tadashi Imanishi

PCDq: human protein complex database with quality index which summarizes different levels of evidences of protein complexes predicted from H-Invitational protein-protein interactions integrative dataset

BackgroundProteins interact with other proteins or biomolecules in complexes to perform cellular functions. Existing protein-protein interaction (PPI) databases and protein complex databases for human proteins are not organized to provide protein complex information or facilitate the discovery of novel subunits. Data integration of PPIs focused specifically on protein complexes, subunits, and their functions. Predicted candidate complexes or subunits are also important for experimental biologists.DescriptionBased on integrated PPI data and literature, we have developed a human protein complex database with a complex quality index (PCDq), which includes both known and predicted complexes and subunits. We integrated six PPI data (BIND, DIP, MINT, HPRD, IntAct, and GNP_Y2H), and predicted human protein complexes by finding densely connected regions in the PPI networks. They were curated with the literature so that missing proteins were complemented and some complexes were merged, resulting in 1,264 complexes comprising 9,268 proteins with 32,198 PPIs. The evidence level of each subunit was assigned as a categorical variable. This indicated whether it was a known subunit, and a specific function was inferable from sequence or network analysis. To summarize the categories of all the subunits in a complex, we devised a complex quality index (CQI) and assigned it to each complex. We examined the proportion of consistency of Gene Ontology (GO) terms among protein subunits of a complex. Next, we compared the expression profiles of the corresponding genes and found that many proteins in larger complexes tend to be expressed cooperatively at the transcript level. The proportion of duplicated genes in a complex was evaluated. Finally, we identified 78 hypothetical proteins that were annotated as subunits of 82 complexes, which included known complexes. Of these hypothetical proteins, after our prediction had been made, four were reported to be actual subunits of the assigned protein complexes.ConclusionsWe constructed a new protein complex database PCDq including both predicted and curated human protein complexes. CQI is a useful source of experimentally confirmed information about protein complexes and subunits. The predicted protein complexes can provide functional clues about hypothetical proteins. PCDq is freely available at http://h-invitational.jp/hinv/pcdq/.

H-InvDB in 2013: an omics study platform for human functional gene and transcript discovery

H-InvDB (http://www.h-invitational.jp/) is a comprehensive human gene database started in 2004. In the latest version, H-InvDB 8.0, a total of 244 709 human complementary DNA was mapped onto the hg19 reference genome and 43 829 gene loci, including nonprotein-coding ones, were identified. Of these loci, 35 631 were identified as potential protein-coding genes, and 22 898 of these were identical to known genes. In our analysis, 19 309 annotated genes were specific to H-InvDB and not found in RefSeq and Ensembl. In fact, 233 genes of the 19 309 turned out to have protein functions in this version of H-InvDB; they were annotated as unknown protein functions in the previous version. Furthermore, 11 genes were identified as known Mendelian disorder genes. It is advantageous that many biologically functional genes are hidden in the H-InvDB unique genes. As large-scale proteomic projects have been conducted to elucidate the functions of all human proteins, we have enhanced the proteomic information with an advanced protein view and new subdatabase of protein complexes (Protein Complex Database with quality index). We propose that H-InvDB is an important resource for finding novel candidate targets for medical care and drug development.

Association between insulin resistance and plasma amino acid profile in non‐diabetic Japanese subjects

Elevation of the branched‐chain amino acids (BCAAs), valine, leucine and isoleucine; and the aromatic amino acids, tyrosine and phenylalanine, has been observed in obesity‐related insulin resistance. However, there have been few studies on Asians, who are generally less obese and less insulin‐resistant than Caucasian or African‐Americans. In the present study, we investigated the relationship between homeostasis model assessment of insulin resistance (HOMA‐IR) and plasma amino acid concentration in non‐diabetic Japanese participants.

Genetic Relationships among Japanese, Northern Han, Hui, Uygur, Kazakh, Greek, Saudi Arabian, and Italian Populations Based on Allelic Frequencies at Four VNTR (D1S80, D4S43, COL2A1, D17S5) and One STR (ACTBP2) Loci

The genetic polymorphism at four variable number of tandem repeats (D1S80, D4S43, COL2A1, D17S5) and one short tandem repeat (ACTBP2) loci was assessed by polymerase chain reaction analysis of genomic DNA obtained from blood samples of eight human populations (Japanese, Northern Han, Hui, Uygur, Kazakh, Saudi Arabian, Greek, Italian). Allele frequencies at all loci were in the Hardy-Weinberg equilibrium for each population. With the exception of ACTBP2, the allelic distribution patterns for these loci revealed a marked genetic divergence among the eight populations. A dendrogram constructed by the neighbor-joining method based on the allele frequencies of the five loci suggested that the five Asian populations (Japanese, Northern Han, Hui, Uygur, and Kazakh) formed one cluster, whereas the two European populations and one West Asian population (Italian, Greek, and Saudi Arabian) formed another. The genetic relationship among these populations may have been greatly influenced by admixture as a result of the migration of individuals along the Silk Road throughout history.

Reconstructing the Demographic History of the Human Lineage Using Whole-Genome Sequences from Human and Three Great Apes

The demographic history of human would provide helpful information for identifying the evolutionary events that shaped the humanity but remains controversial even in the genomic era. To settle the controversies, we inferred the speciation times (T) and ancestral population sizes (N) in the lineage leading to human and great apes based on whole-genome alignment. A coalescence simulation determined the sizes of alignment blocks and intervals between them required to obtain recombination-free blocks with a high frequency. This simulation revealed that the size of the block strongly affects the parameter inference, indicating that recombination is an important factor for achieving optimum parameter inference. From the whole genome alignments (1.9 giga-bases) of human (H), chimpanzee (C), gorilla (G), and orangutan, 100-bp alignment blocks separated by ≥5-kb intervals were sampled and subjected to estimate τ = μT and θ = 4μgN using the Markov chain Monte Carlo method, where μ is the mutation rate and g is the generation time. Although the estimated τHC differed across chromosomes, τHC and τHCG were strongly correlated across chromosomes, indicating that variation in τ is subject to variation in μ, rather than T, and thus, all chromosomes share a single speciation time. Subsequently, we estimated Ts of the human lineage from chimpanzee, gorilla, and orangutan to be 6.0–7.6, 7.6–9.7, and 15–19 Ma, respectively, assuming variable μ across lineages and chromosomes. These speciation times were consistent with the fossil records. We conclude that the speciation times in our recombination-free analysis would be conclusive and the speciation between human and chimpanzee was a single event.

Computational analysis and functional expression of ancestral copepod luciferase.

We recently reported the cDNA sequences of 11 copepod luciferases from the superfamily Augaptiloidea in the order Calanoida. They were classified into two groups, Metridinidae and Heterorhabdidae/Lucicutiidae families, by phylogenetic analyses. To elucidate the evolutionary processes, we have now further isolated 12 copepod luciferases from Augaptiloidea species (Metridia asymmetrica, Metridia curticauda, Pleuromamma scutullata, Pleuromamma xiphias, Lucicutia ovaliformis and Heterorhabdus tanneri). Codon-based synonymous/nonsynonymous tests of positive selection for 25 identified copepod luciferases suggested that positive Darwinian selection operated in the evolution of Heterorhabdidae luciferases, whereas two types of Metridinidae luciferases had diversified via neutral mechanism. By in silico analysis of the decoded amino acid sequences of 25 copepod luciferases, we inferred two protein sequences as ancestral copepod luciferases. They were expressed in HEK293 cells where they exhibited notable luciferase activity both in intracellular lysates and cultured media, indicating that the luciferase activity was established before evolutionary diversification of these copepod species.

A portable system for rapid bacterial composition analysis using a nanopore-based sequencer and laptop computer

We developed a portable system for 16S rDNA analyses consisting of a nanopore technology-based sequencer, the MinION, and laptop computers, and assessed its potential ability to determine bacterial compositions rapidly. We tested our protocols using a mock bacterial community that contained equimolar 16S rDNA and a pleural effusion from a patient with empyema, for time effectiveness and accuracy. MinION sequencing targeting 16S rDNA detected all 20 of the bacterial species present in the mock bacterial community. Time course analysis indicated that the sequence data obtained during the first 5u2009minutes of sequencing (1,379 bacterial reads) were enough to detect all 20 bacteria in the mock sample and to determine species composition, consistent with results of those obtained from 4u2009hours of sequencing (24,202 reads). Additionally, using a clinical sample extracted from the empyema patient’s pleural effusion, we could identify major bacterial pathogens in that effusion using our rapid sequencing and analysis protocol. All results are comparable to conventional 16S rDNA sequencing results using an IonPGM sequencer. Our results suggest that rapid sequencing and bacterial composition determination are possible within 2u2009hours after obtaining a DNA sample.

Full-length transcriptome-based H-InvDB throws a new light on chromosome-centric proteomics.

H-Invitational Database (H-InvDB; http://hinv.jp/ ) is an integrated database of all human genes and transcripts that started in an international collaborative research project for establishing a functional annotation database of human full-length cDNAs. Because H-InvDB contains an abundance of information for human transcripts, including not only well-characterized protein-coding transcripts but also those without experimental evidence at the protein level, this will be a useful information resource for identifying novel and uncharacterized human proteins (so-called missing proteins). By extending predicted protein data in H-InvDB, we developed the H-Inv Extended Protein Database (H-EPD; http://hinv.jp/hinv/h-epd/ ). From now on, we plan to carry out a database-driven proteome research that makes full use of H-EPD to promote discoveries in the current and future C-HPP. Furthermore, we will push forward with the integration of genome, transcriptome, and proteome databases using a unique tool for connecting distributed databases and would like to develop a knowledge discovery system by incorporating data mining tools.

Selection pressure on human STR loci and its relevance in repeat expansion disease

Short Tandem Repeats (STRs) comprise repeats of one to several base pairs. Because of the high mutability due to strand slippage during DNA synthesis, rapid evolutionary change in the number of repeating units directly shapes the range of repeat-number variation according to selection pressure. However, the remaining questions include: Why are STRs causing repeat expansion diseases maintained in the human population; and why are these limited to neurodegenerative diseases? By evaluating the genome-wide selection pressure on STRs using the database we constructed, we identified two different patterns of relationship in repeat-number polymorphisms between DNA and amino-acid sequences, although both patterns are evolutionary consequences of avoiding the formation of harmful long STRs. First, a mixture of degenerate codons is represented in poly-proline (poly-P) repeats. Second, long poly-glutamine (poly-Q) repeats are favored at the protein level; however, at the DNA level, STRs encoding long poly-Qs are frequently divided by synonymous SNPs. Furthermore, significant enrichments of apoptosis and neurodevelopment were biological processes found specifically in genes encoding poly-Qs with repeat polymorphism. This suggests the existence of a specific molecular function for polymorphic and/or long poly-Q stretches. Given that the poly-Qs causing expansion diseases were longer than other poly-Qs, even in healthy subjects, our results indicate that the evolutionary benefits of long and/or polymorphic poly-Q stretches outweigh the risks of long CAG repeats predisposing to pathological hyper-expansions. Molecular pathways in neurodevelopment requiring long and polymorphic poly-Q stretches may provide a clue to understanding why poly-Q expansion diseases are limited to neurodegenerative diseases.

Evolutionary growth process of highly conserved sequences in vertebrate genomes

Genome sequence comparison between evolutionarily distant species revealed ultraconserved elements (UCEs) among mammals under strong purifying selection. Most of them were also conserved among vertebrates. Because they tend to be located in the flanking regions of developmental genes, they would have fundamental roles in creating vertebrate body plans. However, the evolutionary origin and selection mechanism of these UCEs remain unclear. Here we report that UCEs arose in primitive vertebrates, and gradually grew in vertebrate evolution. We searched for UCEs in two teleost fishes, Tetraodon nigroviridis and Oryzias latipes, and found 554 UCEs with 100% identity over 100 bps. Comparison of teleost and mammalian UCEs revealed 43 pairs of common, jawed-vertebrate UCEs (jUCE) with high sequence identities, ranging from 83.1% to 99.2%. Ten of them retain lower similarities to the Petromyzon marinus genome, and the substitution rates of four non-exonic jUCEs were reduced after the teleost-mammal divergence, suggesting that robust conservation had been acquired in the jawed vertebrate lineage. Our results indicate that prototypical UCEs originated before the divergence of jawed and jawless vertebrates and have been frozen as perfect conserved sequences in the jawed vertebrate lineage. In addition, our comparative sequence analyses of UCEs and neighboring regions resulted in a discovery of lineage-specific conserved sequences. They were added progressively to prototypical UCEs, suggesting step-wise acquisition of novel regulatory roles. Our results indicate that conserved non-coding elements (CNEs) consist of blocks with distinct evolutionary history, each having been frozen since different evolutionary era along the vertebrate lineage.