Three amphioxus Wnt genes (AmphiWnt3, AmphiWnt5, and AmphiWnt6) associated with the tail bud: the evolution of somitogenesis in chordates.

The amphioxus tail bud is similar to the amphibian tail bud in having an epithelial organization without a mesenchymal component. We characterize three amphioxus Wnt genes (AmphiWnt3, AmphiWnt5, and AmphiWnt6) and show that their early expression around the blastopore can subsequently be traced into the tail bud; in vertebrate embryos, there is a similar progression of expression domains for Wnt3, Wnt5, and Wnt6 genes from the blastopore lip (or its equivalent) to the tail bud. In amphioxus, AmphiWnt3, AmphiWnt5, and AmphiWnt6 are each expressed in a specific subregion of the tail bud, tentatively suggesting that a combinatorial code of developmental gene expression may help generate specific tissues during posterior elongation and somitogenesis. In spite of similarities within their tail buds, vertebrate and amphioxus embryos differ markedly in the relation between the tail bud and the nascent somites: vertebrates have a relatively extensive zone of unsegmented mesenchyme (i.e., presomitic mesoderm) intervening between the tail bud and the forming somites, whereas the amphioxus tail bud gives rise to new somites directly. It is likely that presomitic mesoderm is a vertebrate innovation made possible by developmental interconversions between epithelium and mesenchyme that first became prominent at the dawn of vertebrate evolution.     

Characterization of amphioxus AmphiWnt8: insights into the evolution of patterning of the embryonic dorsoventral axis

SUMMARY The full-length sequence and developmental expression of an amphioxus Wnt gene ( AmphiWnt8 ) are described. In amphioxus embryos, the expression patterns of AmphiWnt8 suggest patterning roles in the forebrain, in the hindgut, and in the paraxial mesoderm that gives rise to the muscular somites. Phylogenetic analysis indicates that a single Wnt8 subfamily gene in an ancestral chordate duplicated early in vertebrate evolution into a Wnt8 clade and a Wnt8b clade. Coincident with this gene duplication, the functions of the ancestral AmphiWnt8 -like gene appear to have been divided between vertebrate Wnt8b (exclusively neurogenic, especially in the forebrain) and vertebrate Wnt8 (miscellaneous, especially in early somitogenesis). Amphioxus AmphiWnt8 and its vertebrate Wnt8 homologs probably play comparable roles in the early dorsoventral patterning of the embryonic body axis.     

Combining multiple data sets in a likelihood analysis: which models are the best?

Until recently, phylogenetic analyses have been routinely based on homologous sequences of a single gene. Given the vast number of gene sequences now available, phylogenetic studies are now based on the analysis of multiple genes. Thus, it has become necessary to devise statistical methods to combine multiple molecular data sets. Here, we compare several models for combining different genes for the purpose of evaluating the likelihood of tree topologies. Three methods of branch length estimation were studied: assuming all genes have the same branch lengths (concatenate model), assuming that branch lengths are proportional among genes (proportional model), or assuming that each gene has a separate set of branch lengths (separate model). We also compared three models of among-site rate variation: the homogenous model, a model that assumes one gamma parameter for all genes, and a model that assumes one gamma parameter for each gene. On the basis of two nuclear and one mitochondrial amino acid data sets, our results suggest that, depending on the data set chosen, either the separate model or the proportional model represents the most appropriate method for branch length analysis. For all the data sets examined, one gamma parameter for each gene represents the best model for among-site rate variation. Using these models we analyzed alternative mammalian tree topologies, and we describe the effect of the assumed model on the maximum likelihood tree. We show that the choice of the model has an impact on the best phylogeny obtained.     

Isolation of Two Novel WNT Genes, WNT14 and WNT15, One of Which (WNT15) Is Closely Linked to WNT3 on Human Chromosome 17q21

TheWntgene family consists of at least 15 structurally related genes that encode secreted extracellular signaling factors. Wnt proteins function in a range of critical developmental processes in both vertebrates and invertebrates and are implicated in regulation of cell growth and differentiation in certain adult mammalian tissues, including the mammary gland. We have isolated a number of WNT sequences from human genomic DNA, two of which, designated WNT14 and WNT15, represent novel members of theWntgene family. We also isolated WNT sequences from human mammary cDNA and present evidence that WNT13 is expressed in human breast tissue, in addition to those previously described. WNT14 and WNT15 appear to have originated from an ancestral branch of theWntgene family that also includes theWnt9sequences found in jawless and cartilaginous fishes. AWnt14cDNA was also isolated from chicken and a partialWnt15sequence from mouse. We show that human WNT14 maps to chromosome 1 and that WNT15 maps distal to BRCA1 on chromosome 17q21, where it lies within 125 kb of another WNT family member, WNT3.     

Fluorescent in situ hybridisation to amphioxus chromosomes

We describe an efficient protocol for mapping genes and other DNA sequences to amphioxus chromosomes using fluorescent in situ hybridisation. We apply this method to identify the number and location of ribosomal DNA gene clusters and telomere sequences in metaphase spreads of Branchiostoma floridae. We also describe how the locations of two single copy genes can be mapped relative to each other, and demonstrate this by mapping an amphioxus Pax gene relative to a homologue of the Notch gene. These methods have great potential for performing comparative genomics between amphioxus and vertebrates.     

The amphioxus genome enlightens the evolution of the thyroid hormone signaling pathway

Thyroid hormones (THs) have pleiotropic effects on vertebrate development, with amphibian metamorphosis as the most spectacular example. However, developmental functions of THs in non-vertebrate chordates are largely hypothetical and even TH endogenous production has been poorly investigated. In order to get better insight into the evolution of the thyroid hormone signaling pathway in chordates, we have taken advantage of the recent release of the amphioxus genome. We found amphioxus homologous sequences to most of the genes encoding proteins involved in thyroid hormone signaling in vertebrates, except the fast-evolving thyroglobulin: sodium iodide symporter, thyroid peroxidase, deiodinases, thyroid hormone receptor, TBG, and CTHBP. As only some genes encoding proteins involved in TH synthesis regulation were retrieved (TRH, TSH receptor, and CRH receptor but not their corresponding receptors and ligands), there may be another mode of upstream regulation of TH synthesis in amphioxus. In accord with the notion that two whole genome duplications took place at the base of the vertebrate tree, one amphioxus gene often corresponded to several vertebrate homologs. However, some amphioxus specific duplications occurred, suggesting that several steps of the TH pathway were independently elaborated in the cephalochordate and vertebrate lineages. The present results therefore indicate that amphioxus is capable of producing THs. As several genes of the TH signaling pathway were also found in the sea urchin genome, we propose that the thyroid hormone signaling pathway is of ancestral origin in chordates, if not in deuterostomes, with specific elaborations in each lineage, including amphioxus.     

Molecular evolution of the polyamine oxidase gene family in Metazoa

ABSTRACT: BACKGROUND: Polyamine oxidase enzymes catalyze the oxidation of polyamines and acetylpolyamines. Since polyamines are basic regulators of cell growth and proliferation, their homeostasis is crucial for cell life. Members of the polyamine oxidase gene family have been identified in a wide variety of animals, including vertebrates, arthropodes, nematodes, placozoa, as well as in plants and fungi. Polyamine oxidases (PAOs) from yeast can oxidize spermine, N1-acetylspermine, and N1-acetylspermidine, however, in vertebrates two different enzymes, namely spermine oxidase (SMO) and acetylpolyamine oxidase (APAO), specifically catalyze the oxidation of spermine, and N1-acetylspermine/N1-acetylspermidine, respectively. Little is known about the molecular evolutionary history of these enzymes. However, since the yeast PAO is able to catalyze the oxidation of both acetylated and non acetylated polyamines, and in vertebrates these functions are addressed by two specialized polyamine oxidase subfamilies (APAO and SMO), it can be hypothesized an ancestral reference for the former enzyme from which the latter would have been derived. RESULTS: We analysed 36 SMO, 26 APAO and 14 PAO homologue protein sequences from 54 taxa including various vertebrates and invertebrates. The analysis of the full-length sequences and the principal domains of vertebrate and invertebrate PAOs yielded consensus primary protein sequences for vertebrate SMOs and APAOs, and invertebrate PAOs. This analysis, coupled to molecular modeling techniques, also unveiled sequence regions that confer specific structural and functional properties, including substrate specificity, by the different PAO subfamilies. Molecular phylogenetic trees revealed a basal position of all the invertebrates PAO enzymes relative to vertebrate SMOs and APAOs. PAOs from insects constitute a monophyletic clade. Two PAO variants sampled in the amphioxus are basal to the dichotomy between two well supported monophyletic clades including, respectively, all the SMOs and APAOs from vertebrates. The two vertebrate monophyletic clades clustered strictly mirroring the organismal phylogeny of fishes, amphibians, reptiles, birds, and mammals. Evidences from comparative genomic analysis, structural evolution and functional divergence in a phylogenetic framework across Metazoa suggested an evolutionary scenario where the ancestor PAO coding sequence, present in invertebrates as an orthologous gene, has been duplicated in the vertebrate branch to originate the paralogous SMO and APAO genes. A further genome evolution event concerns the SMO gene of placental, but not marsupial and monotremate, mammals which increased its functional variation following an alternative splicing (AS) mechanism. CONCLUSIONS: In this study the explicit integration in a phylogenomic framework of phylogenetic tree construction, structure prediction, and biochemical function data/prediction, allowed inferring the molecular evolutionary history of the PAO gene family and to disambiguate paralogous genes related by duplication event (SMO and APAO) and orthologous genes related by speciation events (PAOs, SMOs/APAOs). Further, while in vertebrates experimental data corroborate SMO and APAO molecular function predictions, in invertebrates the finding of a supported phylogenetic clusters of insect PAOs and the co-occurrence of two PAO variants in the amphioxus urgently claim the need for future structure-function studies.     

青岛文昌鱼过氧化氢酶基因的克隆及生物信息学分析

过氧化氢酶(catalase,CAT)是一类广泛存在于动物、植物和微生物体内的末端氧化酶,是生物体内抗氧化酶体系的重要成员,因此本研究旨在克隆文昌鱼CAT基因并对其进行生物信息学分析。本研究以青岛文昌鱼(Branchiostoma belcheri tsingtauense)为材料,用RACE技术首次克隆了其CAT全长cDNA的基因序列,命名为AmphiCAT(GenBank登录号:KU058636);该序列总长为2640 bp,开放阅读框(open reading frame,ORF)为1533 bp,编码510个氨基酸,含有一个长为19个氨基酸序列的潜在活性位点和一个长为9个氨基酸序列的血红素配体信号,总分子量在线预测约为57.85 kD;经生物软件分析确定该蛋白质无信号肽序列,预测该蛋白质为非分泌性蛋白,属于单功能CAT的clade3分支;该基因的分子进化树表明青岛文昌鱼CAT同软体动物的亲缘关系较近。     

Isolation of pathogen/stress-inducible cDNAs from alfalfa by mRNA differential display

Differential display of mRNA was used to isolate a full-length (SRG1) and a partial (SRG2) alfalfa cDNA induced during infection with the fungal pathogen Colletotrichum trifolii. The deduced amino acid sequences are similar to each other and resemble plant defense-related proteins and tree pollen allergens. SRG1 is a member of a gene family in alfalfa, which may also include the putative defense-related gene PR10. Unlike many defense-related genes described in similar systems, expression of SRG1-like genes does not correlate with resistance to C. trifolii. We speculate SRG1 is induced in response to plant stress.     

Genomic organization and evolution of actin genes in the amphioxus Branchiostoma belcheri and Branchiostoma floridae

We previously described the cDNA cloning and expression patterns of actin genes from amphioxus Branchiostoma floridae (Kusakabe, R., Kusakabe, T., Satoh, N., Holland, N.D., Holland, L.Z., 1997. Differential gene expression and intracellular mRNA localization of amphioxus actin isoforms throughout development: implications for conserved mechanisms of chordate development. Dev. Genes Evol. 207, 203–215). In the present paper, we report the characterization of cDNA clones for actin genes from a closely related species, Branchiostoma belcheri, and the exon–intron organization of B. floridae actin genes. Each of these two amphioxus species has two types of actin genes, muscle and cytoplasmic. The coding and non-coding regions of each type are well-conserved between the two species. A comparison of nucleotide sequences of muscle actin genes between the two species suggests that a gene conversion may have occurred between two B. floridae muscle actin genes BfMA1 and BfMA2. From the conserved positions of introns between actin genes of amphioxus and those of other deuterostomes, the evolution of deuterostome actin genes can be inferred. Thus, the presence of an intron at codon 328/329 in vertebrate muscle and cytoplasmic actin genes but not in any known actin gene in other deuterostomes suggests that a gene conversion may have occurred between muscle and cytoplasmic actin genes during the early evolution of the vertebrates after separation from other deuterostomes. A Southern blot analysis of genomic DNA revealed that the amphioxus genome contains multiple muscle and cytoplasmic actin genes. Some of these actin genes seem to have arisen from recent duplication and gene conversion. Our findings suggest that the multiple genes encoding muscle and cytoplasmic actin isoforms arose independently in each of the three chordate lineages and that gene duplications and gene conversions established the extant actin multigene family during the evolution of chordates.     

Uninode coding vs gene tree parsimony for phylogenetic reconstruction using duplicate genes

Two different methods of using paralogous genes for phylogenetic inference have been proposed: reconciled trees (or gene tree parsimony) and uninode coding. Gene tree parsimony suffers from 10 serious problems, including differential weighting of nucleotide and gap characters, undersampling which can be misinterpreted as synapomorphy, all of the characters not being allowed to interact, and conflict between gene trees being given equal weight, regardless of branch support. These problems are largely avoided by using uninode coding. The uninode coding method is elaborated to address multiple gene duplications within a single gene tree family and handle problems caused by lack of gene tree resolution. An example of vertebrate phylogeny inferred from nine genes is reanalyzed using uninode coding. We suggest that uninode coding be used instead of gene tree parsimony for phylogenetic inference from paralogous genes.     

Cloning, expression, purification, and characterization of AmphiTip30, a member of short-chain dehydrogenases/reductases family from the amphioxus Branchiostoma belcheri tsingtauense

In this paper we describe the cloning, expression and identification study of the TIP30 gene from amphioxus (Branchiostoma belcheri). The amphioxus TIP30 cDNA is comprised of 1499 bp and is translated in one open-reading frame to give a protein of 237 amino acids, with a predicted 23 amino acids signal peptide, a 147 bp 5'-UTR and a 638 bp 3'-UTR. A multiple alignment of TIP30 from amphioxus with other known TIP30 sequences shows the conservation of most amino acid residues involved in the peculiar structural domains found within TIP30's. Phylogenetic analysis places AmphiTIP30 at the base of the phylogenetic tree, suggesting that AmphiTIP30 is the archetype of the vertebrate TIP30 genes. We express the amphioxus TIP30 gene in Escherichia coli. driven by T7 promoter. The recombinant amphioxus TIP30 protein was purified by HisTrap affinity column. Subsequently, the binding constant and enzyme activity was mensurated. Western blot and immunohistochemistry analysis confirmed that amphioxus has a native molecular mass of approximately 26 kDa, and TIP30 was strongly expressed in ovary. Finally, the initial function of TIP30 is discussed.     

Incorporation of gap characters and lineage-specific regions into phylogenetic analyses of gene families from divergent clades: an example from the kinesin superfamily across eukaryotes

The kinesin superfamily across eukaryotes was used to examine how incorporation of gap characters scored from conserved regions shared by all members of a gene family and incorporation of amino acid and gap characters scored from lineage‐specific regions affect gene‐tree inference of the gene family as a whole. We addressed these two questions in the context of two different densities of sequence sampling, four alignment programs, and two methods of tree construction. Taken together, our findings suggest the following. First, gap characters should be incorporated into gene‐tree inference, even for divergent sequences. Second, gene regions that are not conserved among all or most sequences sampled should not be automatically discarded without evaluation of potential phylogenetic signal that may be contained in gap and/or sequence characters. Third, among the four alignment programs evaluated using their default alignment parameters, Clustal may be expected to output alignments that result in the greatest gene‐tree resolution and support. Yet, this high resolution and support should be regarded as optimistic, rather than conservative, estimates. Fourth, this same conclusion regarding resolution and support holds for Bayesian gene‐tree analyses relative to parsimony‐jackknife gene‐tree analyses. We suggest that a more conservative approach, such as aligning the sequences using DIALIGN‐T or MAFFT, analyzing the appropriate characters using parsimony, and assessing branch support using the jackknife, is more appropriate for inferring gene trees of divergent gene families. © The Willi Hennig Society 2007.     

Elucidation of subfamily segregation and intramolecular coevolution of the olfactomedin-like proteins by comprehensive phylogenetic analysis and gene expression pattern assessment

The categorization of genes by structural distinctions relevant to biological characteristics is very important for understanding of gene functions and predicting functional implications of uncharacterized genes. It was absolutely necessary to deploy an effective and efficient strategy to deal with the complexity of the large olfactomedin-like (OLF) gene family sharing sequence similarity but playing diversified roles in many important biological processes, as the simple highest-hit homology analysis gave incomprehensive results and led to inappropriate annotation for some uncharacterized OLF members. In light of evolutionary information that may facilitate the classification of the OLF family and proper association of novel OLF genes with characterized homologs, we performed phylogenetic analysis on all 116 OLF proteins currently available, including two novel members cloned by our group. The OLF family segregated into seven subfamilies and members with similar domain compositions or functional properties all fell into relevant subfamilies. Furthermore, our Northern blot analysis and previous studies revealed that the typical human OLF members in each subfamily exhibited tissue-specific expression patterns, which in turn supported the segregation of the OLF subfamilies with functional divergence. Interestingly, the phylogenetic tree topology for the OLF domains alone was almost identical with that of the full-length tree representing the unique phylogenetic feature of full-length OLF proteins and their particular domain compositions. Moreover, each of the major functional domains of OLF proteins kept the same phylogenetic feature in defining similar topology of the tree. It indicates that the OLF domain and the various domains in flanking non-OLF regions have coevolved and are likely to be functionally interdependent. Expanded by a plausible gene duplication and domain couplings scenario, the OLF family comprises seven evolutionarily and functionally distinct subfamilies, in which each member shares similar structural and functional characteristics including the composition of coevolved and interdependent domains. The phylogenetically classified and preliminarily assessed subfamily framework may greatly facilitate the studying on the OLF proteins. Furthermore, it also demonstrated a feasible and reliable strategy to categorize novel genes and predict the functional implications of uncharacterized proteins based on the comprehensive phylogenetic classification of the subfamilies and their relevance to preliminary functional characteristics.     

The pregnancy-specific glycoprotein family of the immunoglobulin superfamily: Identification of new members and estimation of family size

The members of the carcinoembryonic antigen (CEA)/pregnancy-specific glycoprotein (PSG) gene family have a characteristic N-terminal domain that is homologous to the immunoglobulin variable region. We have estimated the size of the PSG subfamily by identification of N-domain exons from isolated genomic clones and from total genomic DNA through PCR amplification and DNA sequence determination. The PSG subfamily contains at least 11 different genes. For 7 of these, two DNA sequences differing from each other in 1 to 4 nucleotides were detected. Most likely, they represent different alleles. They are PSG1, PSG2, PSG3, PSG4, PSG5, PSG6, PSG7, PSG8, PSG11, PSG12, and PSG13. Six of the N-domain sequences described here are new. All of the PSGs except PSG1, PSG4, and PSG8 contained the arginine-glycine-aspartic acid sequence at position 93-95 corresponding to the complementarity determining region 3 of immunoglobulin. Parsimony analysis of 24 CEA and PSG sequences using 12 members of the immunoglobulin gene superfamily as outgroups to root the family tree shows that the N-domain of the CEA group genes evolved in one major branch and the PSG group genes in the other.