Ancient divergence of animal protein tyrosine kinase genes demonstrated by a gene family tree including choanoflagellate genes.

Animal-specific gene families involved in cell-cell communication and developmental control comprise many subfamilies with distinct domain structures and functions. They diverged by subfamily-generating duplications and domain shufflings before the parazoan-eumetazoan split. Here, we have cloned 40 PTK cDNAs from choanoflagellates, Monosiga ovata, Stephanoeca diplocostata and Codosiga gracilis, the closest relatives to animals. A phylogeny-based analysis of PTKs revealed that 40 out of 47 subfamilies analyzed have unique domain structures and are possibly generated independently in animal and choanoflagellate lineages by domain shufflings. Seven cytoplasmic subfamilies showed divergence before the animal-choanoflagellate split originated by both duplications and shufflings.

The repertoire of desaturases and elongases reveals fatty acid variations in 56 eukaryotic genomes.

The repertoire of biosynthetic enzymes found in an organism is an important clue for elucidating the chemical structural variations of various compounds. In the case of fatty acids, it is essential to examine key enzymes that are desaturases and elongases, whose combination determine the range of fatty acid structures. We systematically investigated 56 eukaryotic genomes to obtain 275 desaturase and 265 elongase homologs. Phylogenetic and motif analysis indicated that the desaturases consisted of four functionally distinct subfamilies and the elongases consisted of two subfamilies. From the combination of the subfamilies, we then predicted the ability to synthesize six types of fatty acids. Consequently, we found that the ranges of synthesizable fatty acids were often different even between closely related organisms. The reason is that, as well as diverging into subfamilies, the enzymes have functionally diverged within the individual subfamilies. Finally, we discuss how the adaptation to individual environments and the ability to synthesize specific metabolites provides some explanation for the diversity of enzyme functions. This study provides an example of a potent strategy to bridge the gap from genomic knowledge to chemical knowledge.

Multiple receptor-like kinase cDNAs from liverwort Marchantia polymorpha and two charophycean green algae, Closterium ehrenbergii and Nitella axillaris: Extensive gene duplications and gene shufflings in the early evolution of streptophytes.

Plant receptor-like kinases (RLKs) comprise a large family with more than several hundred members in vascular plants. The RLK family is thought to have diverged specifically in the plant kingdom, and no family member has been identified in other lineages except for animals and Plasmodium, both of which have RLK related families of small size. To know the time of divergence of RLK family members by gene duplications and domain shufflings, comprehensive isolations of RLK cDNAs were performed from a nonvascular plant, liverwort Marchantia polymorpha and two charophycean green algae, Closterium ehrenbergii, and Nitella axillaris, thought to be the closest relatives to land plants. We obtained twenty-nine, fourteen, and thirteen RLK related cDNAs from M. polymorpha, C. ehrenbergii, and N. axillaris, respectively. The amino acid sequences of these RLKs were compared with those of vascular plants, and phylogenetic trees were inferred by GAMT, a genetic algorithm-based maximum likelihood (ML) method that outputs multiple trees, together with best one. The inferred ML trees revealed ancient gene duplications generating subfamilies with different domain organizations, which occurred extensively at least before the divergence of vascular and nonvascular plants. Rather it remains possible that the extensive gene duplications occurred during the early evolution of streptophytes. Multicellular-specific somatic embryogenesis receptor kinase (SERK) involved in somatic embryogenesis was found in a unicellular alga C. ehrenbergii, suggesting the evolution of SERK by gene recruitment of a unicellular gene.

Evolutionary history and functional implications of protein domains and their combinations in eukaryotes.

BACKGROUND: In higher multicellular eukaryotes, complex protein domain combinations contribute to various cellular functions such as regulation of intercellular or intracellular signaling and interactions. To elucidate the characteristics and evolutionary mechanisms that underlie such domain combinations, it is essential to examine the different types of domains and their combinations among different groups of eukaryotes. RESULTS: We observed a large number of group-specific domain combinations in animals, especially in vertebrates. Examples include animal-specific combinations in tyrosine phosphorylation systems and vertebrate-specific combinations in complement and coagulation cascades. These systems apparently underwent extensive evolution in the ancestors of these groups. In extant animals, especially in vertebrates, animal-specific domains have greater connectivity than do other domains on average, and contribute to the varying number of combinations in each animal subgroup. In other groups, the connectivities of older domains were greater on average. To observe the global behavior of domain combinations during evolution, we traced the changes in domain combinations among animals and fungi in a network analysis. Our results indicate that there is a correlation between the differences in domain combinations among different phylogenetic groups and different global behaviors. CONCLUSION: Rapid emergence of animal-specific domains was observed in animals, contributing to specific domain combinations and functional diversification, but no such trends were observed in other clades of eukaryotes. We therefore suggest that the strategy for achieving complex multicellular systems in animals differs from that of other eukaryotes.

Positive selection in the ComC-ComD system of Streptococcal Species.

Competence-stimulating peptide (CSP) and ComD of the streptococcal species are a pheromone and its receptor, respectively, involved in the regulation of competence for natural genetic transformation. We show here that these molecules have undergone positive selection. This study is the first report of positive selection due to competition among bacterial populations.

Identification of rDNA-specific non-LTR retrotransposons in Cnidaria.

Ribosomal RNA genes are abundant repetitive sequences in most eukaryotes. Ribosomal DNA (rDNA) contains many insertions derived from mobile elements including non-long terminal repeat (non-LTR) retrotransposons. R2 is the well-characterized 28S rDNA-specific non-LTR retrotransposon family that is distributed over at least 4 bilaterian phyla. R2 is a large family sharing the same insertion specificity and classified into 4 clades (R2-A, -B, -C, and -D) based on the N-terminal domain structure and the phylogeny. There is no observation of horizontal transfer of R2; therefore, the origin of R2 dates back to before the split between protostomes and deuterostomes. Here, we in silico identified 1 R2 element from the sea anemone Nematostella vectensis and 2 R2-like retrotransposons from the hydrozoan Hydra magnipapillata. R2 from N. vectensis was inserted into the 28S rDNA like other R2, but the R2-like elements from H. magnipapillata were inserted into the specific sequence in the highly conserved region of the 18S rDNA. We designated the Hydra R2-like elements R8. R8 is inserted at 37 bp upstream from R7, another 18S rDNA-specific retrotransposon family. There is no obvious sequence similarity between targets of R2 and R8, probably because they recognize long DNA sequences. Domain structure and phylogeny indicate that R2 from N. vectensis is the member of the R2-D clade, and R8 from H. magnipapillata belongs to the R2-A clade despite its different sequence specificity. These results suggest that R2 had been generated before the split between cnidarians and bilaterians and that R8 is a retrotransposon family that changed its target from the 28S rDNA to the 18S rDNA.

Improvement in the accuracy of multiple sequence alignment program MAFFT.

In 2002, we developed and released a rapid multiple sequence alignment program MAFFT that was designed to handle a huge (up to approximately 5,000 sequences) and long data (approximately 2,000 aa or approximately 5,000 nt) in a reasonable time on a standard desktop PC. As for the accuracy, however, the previous versions (v.4 and lower) of MAFFT were outperformed by ProbCons and TCoffee v.2, both of which were released in 2004, in several benchmark tests. Here we report a recent extension of MAFFT that aims to improve the accuracy with as little cost of calculation time as possible. The extended version of MAFFT (v.5) has new iterative refinement options, G-INS-i and L-INS-i (collectively denoted as [GL]-INS-i in this report). These options use a new objective function combining the weighted sum-of-pairs (WSP) score and a score similar to COFFEE derived from all pairwise alignments. We discuss the improvement in accuracy brought by this extension, mainly using two benchmark tests released very recently, BAliBASE v.3 (for protein alignments) and BRAliBASE (for RNA alignments). According to BAliBASE v.3, the overall average accuracy of L-INS-i was higher than those of other methods successively released in 2004, although the difference among the most accurate methods (ProbCons, TCoffee v.2 and new options of MAFFT) was small. The advantage in accuracy of [GL]-INS-i became greater for the alignments consisting of approximately 50-100 sequences. By utilizing this feature of MAFFT, we also examined another possible approach to improve the accuracy by incorporating homolog information collected from database. The [GL]-INS-i options are applicable to aligning up to approximately 200 sequences, although not applicable to thousands of sequences because of time and space complexities.

A study of archaeal enzymes involved in polar lipid synthesis linking amino acid sequence information, genomic contexts and lipid composition.

Cellular membrane lipids, of which phospholipids are the major constituents, form one of the characteristic features that distinguish Archaea from other organisms. In this study, we focused on the steps in archaeal phospholipid synthetic pathways that generate polar lipids such as archaetidylserine, archaetidylglycerol, and archaetidylinositol. Only archaetidylserine synthase (ASS), from Methanothermobacter thermautotrophicus, has been experimentally identified. Other enzymes have not been fully examined. Through database searching, we detected many archaeal hypothetical proteins that show sequence similarity to members of the CDP alcohol phosphatidyltransferase family, such as phosphatidylserine synthase (PSS), phosphatidylglycerol synthase (PGS) and phosphatidylinositol synthase (PIS) derived from Bacteria and Eukarya. The archaeal hypothetical proteins were classified into two groups, based on the sequence similarity. Members of the first group, including ASS from M. thermautotrophicus, were closely related to PSS. The rough agreement between PSS homologue distribution within Archaea and the experimentally identified distribution of archaetidylserine suggested that the hypothetical proteins are ASSs. We found that an open reading frame (ORF) tends to be adjacent to that of ASS in the genome, and that the order of the two ORFs is conserved. The sequence similarity of phosphatidylserine decarboxylase to the product of the ORF next to the ASS gene, together with the genomic context conservation, suggests that the ORF encodes archaetidylserine decarboxylase, which may transform archaetidylserine to archaetidylethanolamine. The second group of archaeal hypothetical proteins was related to PGS and PIS. The members of this group were subjected to molecular phylogenetic analysis, together with PGSs and PISs and it was found that they formed two distinct clusters in the molecular phylogenetic tree. The distribution of members of each cluster within Archaea roughly corresponded to the experimentally identified distribution of archaetidylglycerol or archaetidylinositol. The molecular phylogenetic tree patterns and the correspondence to the membrane compositions suggest that the two clusters in this group correspond to archaetidylglycerol synthases and archaetidylinositol synthases. No archaeal hypothetical protein with sequence similarity to known phosphatidylcholine synthases was detected in this study.

MAFFT version 5: improvement in accuracy of multiple sequence alignment.

The accuracy of multiple sequence alignment program MAFFT has been improved. The new version (5.3) of MAFFT offers new iterative refinement options, H-INS-i, F-INS-i and G-INS-i, in which pairwise alignment information are incorporated into objective function. These new options of MAFFT showed higher accuracy than currently available methods including TCoffee version 2 and CLUSTAL W in benchmark tests consisting of alignments of >50 sequences. Like the previously available options, the new options of MAFFT can handle hundreds of sequences on a standard desktop computer. We also examined the effect of the number of homologues included in an alignment. For a multiple alignment consisting of approximately 8 sequences with low similarity, the accuracy was improved (2-10 percentage points) when the sequences were aligned together with dozens of their close homologues (E-value < 10(-5)-10(-20)) collected from a database. Such improvement was generally observed for most methods, but remarkably large for the new options of MAFFT proposed here. Thus, we made a Ruby script, mafftE.rb, which aligns the input sequences together with their close homologues collected from SwissProt using NCBI-BLAST.

The evolutionary repertoires of the eukaryotic-type ABC transporters in terms of the phylogeny of ATP-binding domains in eukaryotes and prokaryotes.

ABC (ATP-binding cassette) transporters play an important role in the communication of various substrates across cell membranes. They are ubiquitous in prokaryotes and eukaryotes, and eukaryotic types (EK-types) are distinguished from prokaryotic types (PK-types) in terms of their genes and domain organizations. The EK-types and PK-types mainly consist of exporters and importers, respectively. Prokaryotes have both the EK-types and the PK-types. The EK-types in prokaryotes are usually called "bacterial multidrug ABC transporters," but they are not well characterized in comparison with the multidrug ABC transporters in eukaryotes. Thus, an exhaustive search of the EK-types among diverse organisms and detailed sequence classification and analysis would elucidate the evolutionary history of EK-types. It would also help shed some light on the fundamental repertoires of the wide variety of substrates through which multidrug ABC transporters in eukaryotes communicate. In this work, we have identified the EK-type ABC transporters in 126 prokaryotes using the profiles of the ATP-binding domain (NBD) of the EK-type ABC transporters from 12 eukaryotes. As a result, 11 clusters were identified from 1,046 EK-types ABC transporters. In particular, two large novel clusters emerged, corresponding to the bacterial multidrug ABC transporters related to the ABCB and ABCC families in eukaryotes, respectively. In the genomic context, most of these genes are located alone or adjacent to genes from the same clusters. Additionally, to detect functional divergences in the NBDs, the Kullback-Leibler divergence was measured among these bacterial multidrug transporters. As a result, several putative functional regions were identified, some corresponding to the predicted secondary structures. We also analyzed a phylogeny of the EK-type ABC transporters in both prokaryotes and eukaryotes, which revealed that the EK-type ABC transporters in prokaryotes have certain repertoires corresponding to the conventional ABC protein groups in eukaryotes. On the basis of these findings, we propose an updated evolutionary hypothesis in which the EK-type ABC transporters in both eukaryotes and prokaryotes consisted of several kinds of ABC transporters in putative ancestor cells before the divergence of eukaryotic and prokaryotic cells.

Identification of cryptochrome DASH from vertebrates.

A new type of cryptochrome, CRY-DASH, has been recently identified. The CRY-DASH proteins constitute the fifth subfamily of the photolyase/cryptochrome family. CRY-DASHs have been identified from Synechocystis sp. PCC 6803, Vibrio cholerae, and Arabidopsis thaliana. The Synechocystis CRY-DASH was the first cryptochrome identified from bacteria, and its biochemical features and tertiary structure have been extensively investigated. To determine how broadly the subfamily is distributed within living organisms, we searched for new CRY-DASH candidates within several databases. We found five sequences as new CRY-DASH candidates, which are derived from four marine bacteria and Neurospora crassa. We also found many CRY-DASH candidates from the EST databases, which included sequences from fish and amphibians. We cloned and sequenced the cDNAs of the zebrafish and Xenopus laevis candidates, based on the EST sequences. The proteins encoded by the two genes were purified and characterized. Both proteins contained folate and flavin cofactors, and have a weak DNA photolyase activity. A phylogenetic analysis revealed that the seven candidates actually belong to the new type of cryptochrome subfamily. This is the first report of the CRY-DASH members from vertebrates and fungi.

Markedly reduced evolutionary rates of transcription factors and cytoplasmic ribosomal RNAs and proteins in higher vertebrates and their evolutionary implications.

To reveal the relationship between organismal evolution and the molecular evolutionary rate, the temporal pattern of evolutionary rates were investigated for various genes during the course of deuterostome evolution. Deuterostome lineage leading to extant mammals was tentatively divided into two periods (the First and the Latter periods) by the time of divergence of bony fishes and mammals. For each of the First and the Latter period, evolutionary rates of 207 gene sets were calculated. In the Latter period, the evolutionary rate was significantly reduced in such informational genes as transcription factors and cytoplasmic ribosomal RNAs and proteins. In contrast, a variety of enzymes and mitochondrial ribosomal proteins evolve at nearly constant rate throughout the First and the Latter periods. The present result suggests that the increase of gene number by extensive gene duplications in the early evolution of vertebrates is responsible for the decrease of evolutionary rate.

Nitrogen-assimilating enzymes in land plants and algae: phylogenic and physiological perspectives.

An important biochemical feature of autotrophs, land plants and algae, is their incorporation of inorganic nitrogen, nitrate and ammonium, into the carbon skeleton. Nitrate and ammonium are converted into glutamine and glutamate to produce organic nitrogen compounds, for example proteins and nucleic acids. Ammonium is not only a preferred nitrogen source but also a key metabolite, situated at the junction between carbon metabolism and nitrogen assimilation, because nitrogen compounds can choose an alternative pathway according to the stages of their growth and environmental conditions. The enzymes involved in the reactions are nitrate reductase (EC 1.6.6.1-2), nitrite reductase (EC 1.7.7.1), glutamine synthetase (EC 6.3.1.2), glutamate synthase (EC 1.4.1.13-14, 1.4.7.1), glutamate dehydrogenase (EC 1.4.1.2-4), aspartate aminotransferase (EC 2.6.1.1), asparagine synthase (EC 6.3.5.4), and phosphoenolpyruvate carboxylase (EC 4.1.1.31). Many of these enzymes exist in multiple forms in different subcellular compartments within different organs and tissues, and play sometimes overlapping and sometimes distinctive roles. Here, we summarize the biochemical characteristics and the physiological roles of these enzymes. We also analyse the molecular evolution of glutamine synthetase, glutamate synthase and glutamate dehydrogenase, and discuss the evolutionary relationships of these three enzymes.

MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform.

A multiple sequence alignment program, MAFFT, has been developed. The CPU time is drastically reduced as compared with existing methods. MAFFT includes two novel techniques. (i) Homo logous regions are rapidly identified by the fast Fourier transform (FFT), in which an amino acid sequence is converted to a sequence composed of volume and polarity values of each amino acid residue. (ii) We propose a simplified scoring system that performs well for reducing CPU time and increasing the accuracy of alignments even for sequences having large insertions or extensions as well as distantly related sequences of similar length. Two different heuristics, the progressive method (FFT-NS-2) and the iterative refinement method (FFT-NS-i), are implemented in MAFFT. The performances of FFT-NS-2 and FFT-NS-i were compared with other methods by computer simulations and benchmark tests; the CPU time of FFT-NS-2 is drastically reduced as compared with CLUSTALW with comparable accuracy. FFT-NS-i is over 100 times faster than T-COFFEE, when the number of input sequences exceeds 60, without sacrificing the accuracy.

Effect of environmental antigens on the Ig diversification and the selection of productive V-J joints in the bursa.

In chickens, a single set of unique functional segments of both Ig H and L chain genes is rearranged during early embryogenesis to generate a pool of B cell progenitors that will be diversified in the bursa by gene conversion, forming the preimmune repertoire. After hatching, bursal cells are exposed to environmental Ags in the bursal lumen. We prepared B cells from each single bursal follicle and used PCR-directed Ig L chain gene analysis to study the differentiation of B cells and the effect of antigenic stimulation from the bursal lumen on the neonatal chicken B cell repertoire formation. Selective amplification of B cell clones with a productive V-J joint was observed during the late embryonic stage, possibly by the interaction with ligands expressed on the bursal stroma and further accelerated in the neonatal chicken. Administration of the artificial Ags into the bursal lumen before the isolation of bursa by bursal duct ligation in the embryo caused a significant increase in lymphocytes with a productive V-J joint in the neonatal chicken bursa compared with the isolated bursa. Intra- and interclonal diversity of a complementarity-determining region measured by an evolutionary distance increased during bursal development. Clonal diversification did not require stimulation by artificial Ags from the bursal lumen. Thus, the preimmune repertoire in the bursa is generated by gene conversion during Ag-independent B cell proliferation, and antigenic stimulation from the bursal epithelium to bursal B cells plays roles in the selection of clones with a productive V-J joint.