DDBJ with new system and face.

DDBJ (http://www.ddbj.nig.ac.jp) collected and released 1 880 115 entries or 1 134 086 245 bases in the period from July 2006 to June 2007. The released data contains the high-throughput cDNAs of cricket and high-quality draft genome of medaka among others. Our computer system has been upgraded since March 2007. Another new aspect is an efficient data retrieval tool that has recently been equipped and served at DDBJ. It is called All-round Retrieval for Sequence and Annotation, which enables the user to search for keywords also in the Feature/Qualifier of the International Nucleotide Sequence Database Collaboration (http://www.insdc.org/). We will also replace our home page with a more efficient one by the end of 2007.

Comparative genomics study reveals Red Sea Bacillus with characteristics associated with potential microbial cell factories (MCFs).

Recent advancements in the use of microbial cells for scalable production of industrial enzymes encourage exploring new environments for efficient microbial cell factories (MCFs). Here, through a comparison study, ten newly sequenced Bacillus species, isolated from the Rabigh Harbor Lagoon on the Red Sea shoreline, were evaluated for their potential use as MCFs. Phylogenetic analysis of 40 representative genomes with phylogenetic relevance, including the ten Red Sea species, showed that the Red Sea species come from several colonization events and are not the result of a single colonization followed by speciation. Moreover, clustering reactions in reconstruct metabolic networks of these Bacillus species revealed that three metabolic clades do not fit the phylogenetic tree, a sign of convergent evolution of the metabolism of these species in response to special environmental adaptation. We further showed Red Sea strains Bacillus paralicheniformis (Bac48) and B. halosaccharovorans (Bac94) had twice as much secreted proteins than the model strain B. subtilis 168. Also, Bac94 was enriched with genes associated with the Tat and Sec protein secretion system and Bac48 has a hybrid PKS/NRPS cluster that is part of a horizontally transferred genomic region. These properties collectively hint towards the potential use of Red Sea Bacillus as efficient protein secreting microbial hosts, and that this characteristic of these strains may be a consequence of the unique ecological features of the isolation environment.

Molecular cloning and characterization of a novel glycoprotein, gp34, that is specifically induced by the human T-cell leukemia virus type I transactivator p40tax.

We have cloned and sequenced a cDNA encoding gp34, a novel glycoprotein expressed in cells bearing human T-cell leukemia virus type I (HTLV-I). HTLV-I has a trans-acting transcriptional activator, p40tax, that is thought to be implicated in leukemogenesis through the activation of cellular enhancers. With a subline (JPX-9) of the human T-cell line Jurkat, in which p40tax is inducible, gp34 was shown to be of cellular origin and to be transcriptionally activated by p40tax. It was also demonstrated that two species of mRNA are generated from one copy of the gp34 gene and that these mRNAs encode the identical gp34 product and differ in the 3' untranslated region. Analysis of the deduced amino acid sequence of gp34 showed that it lacks typical signal peptides; however, it has a hydrophobic stretch for membrane anchoring and four possible N-linked glycosylation sites at the carboxy-terminal portion, indicating that it belongs to the family of membrane proteins whose carboxy-terminal portion protrudes out of the cell. The gp34 gene displayed relatively delayed induction compared with other genes activated by p40tax. Taken together with the observation of the dependence of gp34 expression on HTLV-I p40tax, unlike other p40tax-dependent genes such as those for the interleukin-2 receptor alpha chain and c-fos, which are expressed or induced under physiological conditions, we predict that the mechanism involved in the induction of gp34 expression by p40tax is distinct from and more intricate than those for the previously characterized genes.

Precise switching of DNA replication timing in the GC content transition area in the human major histocompatibility complex.

The human genome is composed of long-range G+C% (GC%) mosaic structures thought to be related to chromosome bands. We previously reported a boundary of megabase-sized GC% mosaic domains at the junction area between major histocompatibility complex (MHC) classes II and III, proposing it as a possible chromosome band boundary. DNA replication timing during the S phase is known to be correlated cytogenetically with chromosome band zones, and thus the band boundaries have been predicted to contain a switch point for DNA replication timing. In this study, to identify to the nucleotide sequence level the replication switch point during the S phase, we determined the precise DNA replication timing for MHC classes II and III, focusing on the junction area. To do this, we used PCR-based quantitation of nascent DNA obtained from synchronized human myeloid leukemia HL60 cells. The replication timing changed precisely in the boundary region with a 2-h difference between the two sides, supporting the prediction that this region may be a chromosome band boundary. We supposed that replication fork movement terminates (pauses) or significantly slows in the switch region, which contains dense Alu clusters; polypurine/polypyrimidine tracts; di-, tri-, or tetranucleotide repeats; and medium-reiteration-frequency sequences. Because the nascent DNA in the switch region was recovered at low efficiency, we investigated whether this region is associated with the nuclear scaffold and found three scaffold-associated regions in and around the switch region.

DDBJ in collaboration with mass-sequencing teams on annotation.

In the past year, we at DDBJ (DNA Data Bank of Japan; http://www.ddbj.nig.ac.jp) collected and released 1,066,084 entries or 718,072,425 bases including the whole chromosome 22 of chimpanzee, the whole-genome shotgun sequences of silkworm and various others. On the other hand, we hosted workshops for human full-length cDNA annotation and participated in jamborees of mouse full-length cDNA annotation. The annotated data are made public at DDBJ. We are also in collaboration with a RIKEN team to accept and release the CAGE (Cap Analysis Gene Expression) data under a new category, MGA (Mass Sequences for Genome Annotation). The data will be useful for studying gene expression control in many aspects.

Genomic annotation of 15,809 ESTs identified from pooled early gestation human eyes.

To complement cDNA libraries from the human eye at early gestation and to discover candidate genes associated with early ocular development, we used freshly dissected human eyeballs from week 9-14 of gestation to construct the early human fetal eye cDNA library. A total of 15,809 clones were isolated and sequenced from the unamplified and unnormalized library. We screened 11,246 good-quality ESTs, leading to the identification of 5,534 nonredundant clusters. Among them, 4,010 (72%) genes matched in the human protein database (Ensembl). The remaining 28% (1,524) corresponded to potentially novel or previously unidentified ESTs. We used BLASTX to compare our EST data with eight organisms and found common expression of a high portion of genes: Caenorhabditis briggsae (26%), Caenorhabditis elegans (27%), Anopheles gambiae (37%), Drosophila melanogaster (32%), Danio rerio (42%), Fugu rubripes (49%), Rattus norvegicusvalitus (52%), and Mus musculus (59%). Nevertheless, 48% (2,680 of 5,534) of the genes expressed in the early developing eye were not shared with current NEIBank human eye cDNA data. In addition, eight known retinal disease genes existed in our ESTs. Among them, six (COL11A1, BBS5, PDE6B, OAT, VMD2, and PGK1) were conserved among the genomes of other organisms, indicating that our annotated EST set provides not only a valuable resource for gene discovery and functional genomic analysis but also for phylogenetic analysis. Our foremost early gestation human eye cDNA library could provide detailed comparisons across species to identify physiological functions of genes and to elucidate evolutionary mechanisms.

Molecular characterization of the developmental gene in eyes: through data-mining on integrated transcriptome databases.

OBJECTIVES: Our aim was to utilize publicly available and proprietary sources to discover candidate genes important for ocular development. DESIGN AND METHODS: The collated information on our 5092 non-redundant clusters was grouped and functional annotation was conducted using gene ontology (FatiGO) for categorizing them with respect to molecular function. The web-based viewer technological platform (H-InvDB) was employed for transcription analyses of in-house high quality fetal eye Expressed Sequence Tags (ESTs). Eye-specific ESTs were also analyzed across species by using EMBEST. RESULTS: According to adult eye cDNA libraries, nucleic acid binding and cell structure/cytoskeletal protein genes were the most abundant among the ESTs of fetal eyes. Using cDNA assembly in H-InvDB, 20 (80%) of the 25 most commonly expressed genes in the human eye are also expressed in extraocular tissues. The crystalline gamma S gene is highly expressed in the eye, but not in other tissues. We used EMBEST to compare human fetal eye and octopus eye ESTs and the expression similarity was low (1.6%). This indicated that our fetal eye library contains genes necessary for the developmental process and biological function of the eye, which may not be expressed in the fully developed octopus eyes. The human fetal eye cDNA library also contained highly abundant eye tissue genes, including alphaA-crystallin, eukaryotic translation elongation factor 1 alpha 1 (EEF1A1), bestrophin (VMD2), cystatin C, and transforming growth factor, beta-induced (BIGH3). CONCLUSIONS: Our annotated EST set provides a valuable resource for gene discovery and functional genomic analysis. This display will help to appreciate the strengths and weaknesses of the different technological platforms, so that in future studies the maximum amount of beneficial information can be derived from the appropriate use of each method.

DNA Data Bank of Japan (DDBJ) in XML.

The DNA Data Bank of Japan (DDBJ, http://www.ddbj.nig.ac.jp) has collected and released more entries and bases than last year. This is mainly due to large-scale submissions from Japanese sequencing teams on mouse, rice, chimpanzee, nematoda and other organisms. The contributions of DDBJ over the past year are 17.3% (entries) and 10.3% (bases) of the combined outputs of the International Nucleotide Sequence Databases (INSD). Our complete genome sequence database, Genome Information Broker (GIB), has been improved by incorporating XML. It is now possible to perform a more sophisticated database search against the new GIB than the ordinary BLAST or FASTA search.

DDBJ in the stream of various biological data.

In the past year we at DDBJ (http://www.ddbj.nig. ac.jp) have made a steady increase in the number of data submissions with a 50.6% increment in the number of bases or 46.5% increment in the number of entries. Among them the genome data of man, ascidian and rice hold the top three. Our activity has extended to providing a tool that enables sequence retrieval using regular expressions, and to launching our SOAP server and web services to facilitate the acquisition of proper data and tools from a huge number of biological data resources on websites worldwide. We have also opened our public gene expression database, CIBEX.

DNA Data Bank of Japan (DDBJ) for genome scale research in life science.

The DNA Data Bank of Japan (DDBJ, http://www.ddbj.nig.ac.jp) has made an effort to collect as much data as possible mainly from Japanese researchers. The increase rates of the data we collected, annotated and released to the public in the past year are 43% for the number of entries and 52% for the number of bases. The increase rates are accelerated even after the human genome was sequenced, because sequencing technology has been remarkably advanced and simplified, and research in life science has been shifted from the gene scale to the genome scale. In addition, we have developed the Genome Information Broker (GIB, http://gib.genes.nig.ac.jp) that now includes more than 50 complete microbial genome and Arabidopsis genome data. We have also developed a database of the human genome, the Human Genomics Studio (HGS, http://studio.nig.ac.jp). HGS provides one with a set of sequences being as continuous as possible in any one of the 24 chromosomes. Both GIB and HGS have been updated incorporating newly available data and retrieval tools.

Codon substitution in evolution and the "saturation" of synonymous changes.

A mathematical model for codon substitution is presented, taking into account unequal mutation rates among different nucleotides and purifying selection. This model is constructed by using a 61 X 61 transition probability matrix for the 61 nonterminating codons. Under this model, a computer simulation is conducted to study the numbers of silent (synonymous) and amino acid-altering (nonsynonymous) nucleotide substitutions when the underlying mutation rates among the four kinds of nucleotides are not equal. It is assumed that the substitution rates are constant over evolutionary time, the codon frequencies being in equilibrium, and, thus, the numbers of synonymous and nonsynonymous substitutions both increase linearly with evolutionary time. It is shown that, when the mutation rates are not equal, the estimate of synonymous substitutions obtained by F. Perler, A. Efstratiadis, P. Lomedico, W. Gilbert, R. Kolodner and J. Dodgson's "Percent Corrected Divergence" method increases nonlinearly, although the true number of synonymous substitutions increases linearly. It is, therefore, possible that the "saturation" of synonymous substitutions observed by Perler et al. is due to the inefficiency of their method to detect all synonymous substitutions.

Rates of synonymous substitution and base composition of nuclear genes in Drosophila.

We compared the rates of synonymous (silent) substitution among various genes in a number of species of Drosophila. First, we found that even for a particular gene, the rate of synonymous substitution varied considerably with Drosophila lineages. Second, we showed a large variation in synonymous substitution rates among nuclear genes in Drosophila. These rates of synonymous substitution were correlated negatively with C content and positively with A content at the third codon positions. Nucleotide sequences were also compared between pseudogenes and their functional homologs. The C content of the pseudogenes was lower than that of the functional genes and the A content of the former was higher than that of the latter. Because the synonymous substitution for functional genes and the nucleotide substitution for pseudogenes are exempted from any selective constraint at the protein level, these observations could be explained by a biased pattern of mutation in the Drosophila nuclear genome. Such a bias in the mutation pattern may affect the molecular clock (local clock) of each nuclear gene of each species. Finally, we obtained the average rates of synonymous substitution for three gene groups in Drosophila; 11.0 x 10(-9), 17.5 x 10(-9) and 27.1 x 10(-9)/site/year.

Establishment of a phylogenetic survey system for AIDS-related lentiviruses and demonstration of a new HIV-2 subgroup.

We designed a universal primer (UNIPOL) for DNA amplification of AIDS-related viruses. The phylogenetic tree constructed from the presumed sequences amplified with UNIPOL was representative of the tree calculated from whole pol gene sequences so far reported. UNIPOL was able to amplify the sequences of all four major groups of primate lentiviruses and also that of a distinct virus from a Ghanaian patient with an AIDS-related complex, designated GH-2. This strain scarcely hybridizes with known HIV/simian immunodeficiency virus (SIV) DNA probes. Sequence analysis of the only amplified fragment revealed rapidly that GH-2 was quite similar to the recently reported HIV-2ALT(D205) and that these two viruses form a new subgroup distint from known HIV-2 and SIVmac/SIVsm in the large HIV-2 group. This system will be useful for further phylogenetic study of various primate lentiviruses.

How can human and simian immunodeficiency viruses utilize chemokine receptors as their coreceptors?

Several chemokine receptors (CKRs) act as coreceptors of human immunodeficiency virus type 1 (HIV-1), type 2 (HIV-2) and simian immunodeficiency virus (SIV). These CKRs interact with the V3 domain of the envelope (env) protein of HIV/SIV. In this study, we found that the amino acid sequences of two chemokines (SDF-1beta and RANTES), whose receptors (CXCR4 and CCR5) act as major coreceptors for HIV-1, HIV-2 or SIV, showed statistically significant similarity to those of the region containing the third variable (V3) and the third conserved (C3) domains (the V3--C3 domain) of the env protein of HIV-1 and HIV-2. We made a multiple alignment of amino acid sequences for 24 chemokines and the region encompassing the second conserved (C2), V3 and C3 domains (the C2--V3--C3 region) of 10 strains of HIV/SIV. Surprisingly, the hydropathic profile and several important amino acids for protein conformation, such as cysteine and tryptophan, are remarkably conserved between chemokines and the V3--C3 region of HIV/SIV. Moreover, hydrophobic amino acids, such as leucine, isoleucine and valine, are found to be clustered both in the amino-terminal region of chemokines and the C2 domain of HIV/SIV. Thus, chemokines have significantly similar profiles of amino acid properties to those of the C2--V3--C3 region of the env protein of HIV/SIV. These findings raise a hypothesis that chemokines and the C2--V3--C3 region have a common origin. Namely, the HIV/SIV ancestor incorporated a chemokine gene into its env gene. The captured chemokine gene has rapidly diverged by frequent mutations specific to the retroviral genome, and thereby obtained the ability to interact with various CKRs in a short period of time. This paper proposes that the capture of a ligand gene of the host cells into the viral genome may be one of the important mechanisms of viral evolution to expand its host range and generate new viral species.

Positively selected amino acid sites in the entire coding region of hepatitis C virus subtype 1b.

To predict the amino acid sites important for the clearance of hepatitis C virus (HCV) subtype 1b in vivo, positively selected amino acid sites were detected by analyzing the sequence data collected from the international DNA databank. The rate of nonsynonymous substitutions per nonsynonymous site was compared with that of synonymous substitutions per synonymous site for each codon site in the entire coding region. As a result, 13 out of 3010 amino acid sites were found to be positively selected. Among the 13 positively selected amino acid sites, eight were located in the structural proteins and five were in the nonstructural proteins. Moreover, eight were located in B-cell epitopes and two were in T-cell epitopes. These observations suggest that both the antibody and the cytotoxic T lymphocyte are involved in the clearance of HCV subtype 1b in vivo. These positively selected amino acid sites represent candidate vaccination targets for HCV subtype 1b.

Significant differences between the G+C content of synonymous codons in orthologous genes and the genomic G+C content.

The relationship between the overall G+C content of the genome (GC) and the GC content at the third codon positions (GC3) of genes, which we refer to as a GC3-plot, was examined using 15 currently available complete genome sequences. A remarkably linear relationship was found between these two quantities, confirming previous observations of a strong positive correlation in the GC3-plot. In order to conduct a more detailed analysis of the GC3-plot, we examined the GC3 content by separating orthologous codons into three categories: synonymously different codons (namely identical amino acids, IA), different amino acids (DA), and identical codons (IC), for a pairwise comparison of two closely related species. When we took pairwise species comparisons between Mycoplasma genitalium (Mg) and Mycoplasma pneumoniae (Mp) and between Mycobacterium tuberculosis (Mt) and Mycobacterium leprae (Ml) as examples, we found that for Mp and Ml, the GC3 for IA deviated the most from the linear expectation in the GC3-plot, whereas for Mg and Mt the deviation was minimal. These findings suggest that the major changes of GC content took place in Mp and Ml, but not in Mg and Mt. This analysis also enables us to predict the future direction of the evolutionary changes of the genomic GC content.

Evolution of heteromorphic sex chromosomes in the order Aulopiformes.

The fish order Aulopiformes contains both synchronously hermaphroditic and gonochoristic species. From the cytogenetic viewpoint, few reports show that gonochoristic Aulopiformes have heteromorphic sex chromosomes. Because fish in this order give us a unique opportunity to elucidate the evolution of sex chromosomes, it is important to examine a phylogenetic relationship in Aulopiformes by both molecular evolutionary and cytogenetic methods. Thus, we conducted molecular phylogenetic and cytogenetic studies of six Aulopiform species. Our results suggested that hermaphroditic species were evolutionarily derived from gonochoristic species. It follows that the hermaphroditic species might have lost the heteromorphic sex chromosomes during evolution. Here, we suggest a possibility that heteromorphic sex chromosomes can disappear from the genome, even if they have appeared once in evolution. Taking into account Ohno's hypothesis that heteromorphic sex chromosomes might have emerged from autosomes, we propose the hypothesis that heteromorphic sex chromosomes may have undergone repeated events of appearance and disappearance during the course of fish evolution.

Bacterial features in the genome of Methanococcus jannaschii in terms of gene composition and biased base composition in ORFs and their surrounding regions.

As a result of genome projects, the complete nucleotide sequence of the entire genome of an archaeon, Methanococcus jannaschii, was recently determined as well as other complete sequences of bacterial and eucaryal genomes. When all the 1680 predicted protein-coding [corrected] genes of M. jannaschii were classified on the basis of sequence similarity, it was found that this archaeon had a chimeric set of 1016 bacterial-type, 471 eucaryal-type and 193 species- or archaebacteria-specific genes. However, most of the genes predicted to be involved in translation and transcription pathways including RNA genes were of the eucaryal-type with only a few exceptions such as 16S ribosomal RNA and some translation factor-like genes. This appeared curious since previous studies indicated that methanogens have bacterial features in gene organization and expression. To understand the apparent inconsistency between physiological observations and the result of the classification of genes for transcription and translation, we examined the structural relatedness of the genome of M. jannaschii to those of other species. In practice, we compared base compositional patterns in ORFs and their surrounding regions. This made it possible to reveal the relationships among the translation- and transcription-related structures in genomes. In this study, we conducted a statistical test called 'G-test' to evaluate the base biases around the boundaries of ORFs. We then found that M. jannaschii possesses more bacterial features in base biases than eucaryal ones, e.g. strong G biases at the positions corresponding to the Shine-Dalgarno site. This indicates that the few exceptional bacterial genes for translation, such as 16S ribosomal RNA and translation factor-like genes, play crucial roles in the translation pathway in M. jannaschii. The possibility that the genome structure in the last common ancestor of all present species was bacterial is discussed.

Evolutionary significance of intra-genome duplications on human chromosomes.

Phylogenetic analyses indicated that a series of paralogous gene pairs, found in two extensive regions on human chromosomal bands 6p21.3 and 9q33-34, were created by at least two independent duplications. The duplicated genes on chromosomal band 6p21.3 include the genes for type 11 collagen alpha2 subunit (COL11A2), NOTCH4 (mouse int-3 homologue), 70 kDa heat shock protein (HSPA1A, HSPA1B, and HSPA1L), valyl-tRNA synthetase 2 (VARS2), complement components (C2 and C4), pre-B cell leukemia transcription factor 2 (PBX2), retinoid X receptor beta (RXRB), NAT/RING3, and four other proteins. Their paralogous genes on chromosomal band 9q33-34 are genes for type 5 collagen alpha1 subunit (COL5A1), NOTCH1, 78 kDa glucose-regulated protein (HSPA5), valyl-tRNA synthetase 1 (VARS1), complement component V (C5), PBX3, retinoid X receptor alpha (RXRA), ORFX/RING3L, and others. Among these, the genes for collagen, complement components, NAT/RING3, PBX, and RXR appear to have been duplicated around the time of vertebrate emergence, supporting the idea that they were duplicated simultaneously at that time. Another group of genes that includes NOTCH and HSP appear to have diverged long before that time. A comparison of the physical maps of these two regions revealed that the genes which duplicated in the same period were arranged in almost the same order in the two regions, with the assumption of a few chromosomal rearrangements. We propose a possible model for the evolution of these regions, taking into account the molecular mechanisms of regional duplication, gene duplication, translocation, and inversion. We also propose that a comparative mapping of paralogous genes within the human genome would be useful for identifying new genes.

Nucleotide sequence and molecular evolution of mouse retrovirus-like IAP elements.

We determined the nucleotide (nt) sequences of cDNA and genomic clones for murine intracisternal type A particle (IAP) elements, which are retrovirus-like repetitive sequences in rodent genomes. The nucleotide sequence of the cDNA resembled that of retrovirus RNA genomes in its lack of the U5 sequence within the 3' long terminal repeat. By sequence comparison of our clones with reported rodent IAP elements, we located the probable gag, pol and env gene regions. The sequences for the pol, env and the 3' two-thirds of the gag region were conserved among the IAP elements. In the regions, synonymous substitutions occurred more frequently than non-synonymous ones, which suggested that the regions in question were functionally constrained until fairly recently. The rate of nucleotide substitutions in the regions was estimated to be 6-10 X 10(-9) nt per site per year, and significantly higher than that of the cellular genes. These rates may exemplify a characteristic of the nucleotide substitutions for an endogenous retrovirus. The sequence homology between the IAP element and IgE-binding factor gene is discussed.