Molecular evolution of a developmental pathway: phylogenetic analyses of transforming growth factor-beta family ligands, receptors and Smad signal transducers.

Intercellular signaling by transforming growth factor-beta (TGF-beta) proteins coordinates developmental decisions in many organisms. A receptor complex and Smad signal transducers are required for proper responses to TGF-beta signals. We have taken a phylogenetic approach to understanding the developmental evolutionary history of TGF-beta signaling pathways. We were interested in detecting evolutionary influences among the physically interacting multigene families encoding TGF-beta ligands, receptors, and Smads. Our analyses included new ligands and Smads identified from genomic sequence as well as the newest published family members. From an evolutionary perspective we find that (1) TGF-beta pathways do not predate the divergence of animals, plants, and fungi; (2) ligands of the TGF-beta/activin subfamily likely originated after the divergence of nematodes and arthropods; (3) type I receptors from Caenorhabditis elegans are distinct from other receptors and may reflect an ancestral transitional state between type I and type II receptors; and (4) the Smad family appears to be evolving faster than, and independently of, ligands and receptors. From a developmental perspective we find (1) numerous phylogenetic associations not previously detected in each multigene family; (2) that there are unidentified pathway components that discriminate between type I and type II receptors; (3) that there are more Smads to be discovered in Drosophila and mammals; and (4) that the number of C-terminal serines is the best predictor of a Smad's role in TGF-beta signal transduction. We discuss these findings with respect to the coevolution of physically interacting genes.     

Molecular evolution of the 14-3-3 protein family

Members of the highly conserved and ubiquitous 14-3-3 protein family modulate a wide variety of cellular processes. To determine the evolutionary relationships among specific 14-3-3 proteins in different plant, animal, and fungal species and to initiate a predictive analysis of isoform-specific differences in light of the latest functional and structural studies of 14-3-3, multiple alignments were constructed from forty-six 14-3-3 sequences retrieved from the GenBank and SwissProt databases and a newly identified second 14-3-3 gene fromCaenorhabditis elegans. The alignment revealed five highly conserved sequence blocks. Blocks 2–5 correlate well with the alpha helices 3, 5, 7, and 9 which form the proposed internal binding domain in the three-dimensional structure model of the functioning dimer. Amino acid differences within the functional and structural domains of plant and animal 14-3-3 proteins were identified which may account for functional diversity amongst isoforms. Protein phylogenic trees were constructed using both the maximum parsimony and neighbor joining methods of the PHYLIP(3.5c) package; 14-3-3 proteins fromEntamoeba histolytica, an amitochondrial protozoa, were employed as an outgroup in our analysis. Epsilon isoforms from the animal lineage form a distinct grouping in both trees, which suggests an early divergence from the other animal isoforms. Epsilons were found to be more similar to yeast and plant isoforms than other animal isoforms at numerous amino acid positions, and thus epsilon may have retained functional characteristics of the ancestral protein. The known invertebrate proteins group with the nonepsilon mammalian isoforms. Most of the current 14-3-3 isoform diversity probably arose through independent duplication events after the divergence of the major eukaryotic kingdoms. Divergence of the seven mammalian isoforms beta, zeta, gamma, eta, epsilon, tau, and sigma (stratifin/ HME1) occurred before the divergence of mammalian and perhaps before the divergence of vertebrate species. A possible ancestral 14-3-3 sequence is proposed. Correspondence to: D.C. Shakes     

Rates of nuclear and cytoplasmic mitochondrial DNA sequence divergence in mammals

Differential rates of nucleotide substitution among different gene segments and between distinct evolutionary lineages is well documented among mitochondrial genes and is likely a consequence of locus-specific selective constraints that delimit mutational divergence over evolutionary time. We compared sequence variation of 18 homologous loci (15 coding genes and 3 parts of the control region) among 10 mammalian mitochondrial DNA genomes which allowed us to describe different mitochondrial evolutionary patterns and to produce an estimation of the relative order of gene divergence. The relative rates of divergence of mitochondrial DNA genes in the family Felidae were estimated by comparing their divergence from homologous counterpart genes included in nuclear mitochondrial DNA (Numt, pronounced "new might"), a genomic fossil that represents an ancient transfer of 7.9 kb of mitochondrial DNA to the nuclear genome of an ancestral species of the domestic cat (Felis catus). Phylogenetic analyses of mitochondrial (mtDNA) sequences with multiple outgroup species were conducted to date the ancestral node common to the Numt and the cytoplasmic (Cymt) mtDNA genes and to calibrate the rate of sequence divergence of mitochondrial genes relative to nuclear homologous counterparts. By setting the fastest substitution rate as strictly mutational, an empirical "selective retardation index" is computed to quantify the sum of all constraints, selective and otherwise, that limit sequence divergence of mitochondrial gene sequences over time.      

The Molecular Evolution of Actin

We have investigated the molecular evolution of plant and nonplant actin genes comparing nucleotide and amino acid sequences of 20 actin genes. Nucleotide changes resulting in amino acid substitutions (replacement substitutions) ranged from 3-7% for all pairwise comparisons of animal actin genes with the following exceptions. Comparisons between higher animal muscle actin gene sequences and comparisons between higher animal cytoplasmic actin gene sequences indicated less than 3% divergence. Comparisons between plant and nonplant actin genes revealed, with two exceptions, 11-15% replacement substitution. In the analysis of plant actins, replacement substitution between soybean actin genes SAc1, SAc3, SAc4 and maize actin gene MAc1 ranged from 8-10%, whereas these members within the soybean actin gene family ranged from 6-9% replacement substitution. The rate of sequence divergence of plant actin sequences appears to be similar to that observed for animal actins. Furthermore, these and other data suggest that the plant actin gene family is ancient and that the families of soybean and maize actin genes have diverged from a single common ancestral plant actin gene that originated long before the divergence of monocots and dicots. The soybean actin multigene family encodes at least three classes of actin. These classes each contain a pair of actin genes that have been designated kappa (SAc1, SAc6), lambda (SAc2, SAc4) and mu (SAc3, SAc7). The three classes of soybean actin are more divergent in nucleotide sequence from one another than higher animal cytoplasmic actin is divergent from muscle actin. The location and distribution of amino acid changes were compared between actin proteins from all sources. A comparison of the hydropathy of all actin sequences, except from Oxytricha, indicated a strong similarity in hydropathic character between all plant and nonplant actins despite the greater number of replacement substitutions in plant actins. These protein sequence comparisons are discussed with respect to the demonstrated and implicated roles of actin in plants and animals, as well as the tissue-specific expression of actin.     

The primary structure of fatty-acid-binding protein from nurse shark liver. Structural and evolutionary relationship to the mammalian fatty-acid-binding protein family.

The primary structure of a fatty-acid-binding protein (FABP) isolated from the liver of the nurse shark (Ginglymostoma cirratum) was determined by high-performance tandem mass spectrometry (employing multichannel array detection) and Edman degradation. Shark liver FABP consists of 132 amino acids with an acetylated N-terminal valine. The chemical molecular mass of the intact protein determined by electrospray ionization mass spectrometry (Mr = 15124 +/- 2.5) was in good agreement with that calculated from the amino acid sequence (Mr = 15121.3). The amino acid sequence of shark liver FABP displays significantly greater similarity to the FABP expressed in mammalian heart, peripheral nerve myelin and adipose tissue (61-53% sequence similarity) than to the FABP expressed in mammalian liver (22% similarity). Phylogenetic trees derived from the comparison of the shark liver FABP amino acid sequence with the members of the mammalian fatty-acid/retinoid-binding protein gene family indicate the initial divergence of an ancestral gene into two major subfamilies: one comprising the genes for mammalian liver FABP and gastrotropin, the other comprising the genes for mammalian cellular retinol-binding proteins I and II, cellular retinoic-acid-binding protein myelin P2 protein, adipocyte FABP, heart FABP and shark liver FABP, the latter having diverged from the ancestral gene that ultimately gave rise to the present day mammalian heart-FABP, adipocyte FABP and myelin P2 protein sequences. The sequence for intestinal FABP from the rat could be assigned to either subfamily, depending on the approach used for phylogenetic tree construction, but clearly diverged at a relatively early evolutionary time point. Indeed, sequences proximately ancestral or closely related to mammalian intestinal FABP, liver FABP, gastrotropin and the retinoid-binding group of proteins appear to have arisen prior to the divergence of shark liver FABP and should therefore also be present in elasmobranchs. The presence in shark liver of an FABP which differs substantially in primary structure from mammalian liver FABP, while being closely related to the FABP expressed in mammalian heart muscle, peripheral nerve myelin and adipocytes, opens a further dimension regarding the question of the existence of structure-dependent and tissue-specific specialization of FABP function in lipid metabolism.     

The olfactory receptor gene superfamily: data mining, classification, and nomenclature

The vertebrate olfactory receptor (OR) subgenome harbors the largest known gene family, which has been expanded by the need to provide recognition capacity for millions of potential odorants. We implemented an automated procedure to identify all OR coding regions from published sequences. This led us to the identification of 831 OR coding regions (including pseudogenes) from 24 vertebrate species. The resulting dataset was subjected to neighbor-joining phylogenetic analysis and classified into 32 distinct families, 14 of which include only genes from tetrapodan species (Class II ORs). We also report here the first identification of OR sequences from a marsupial (koala) and a monotreme (platypus). Analysis of these OR sequences suggests that the ancestral mammal had a small OR repertoire, which expanded independently in all three mammalian subclasses. Classification of ``fish-like' (Class I) ORs indicates that some of these ancient ORs were maintained and even expanded in mammals. A nomenclature system for the OR gene superfamily is proposed, based on a divergence evolutionary model. The nomenclature consists of the root symbol `OR', followed by a family numeral, subfamily letter(s), and a numeral representing the individual gene within the subfamily. For example, OR3A1 is an OR gene of family 3, subfamily A, and OR7E12P is an OR pseudogene of family 7, subfamily E. The symbol is to be preceded by a species indicator. We have assigned the proposed nomenclature symbols for all 330 human OR genes in the database. A WWW tool for automated name assignment is provided. Received: / Accepted:     

Evolution of arthropod hemocyanins and insect storage proteins (hexamerins)

Crustacean and cheliceratan hemocyanins (oxygen-transport proteins) and insect hexamerins (storage proteins) are homologous gene products, although the latter do not bind oxygen and do not possess the copper- binding histidines present in the hemocyanins. An alignment of 19 amino acid sequences of hemocyanin subunits and insect hexamerins was made, based on the conservation of elements of secondary structure observed in X-ray structures of two hemocyanin subunits. The alignment was analyzed using parsimony and neighbor-joining methods. Results provide strong indications for grouping together the sequences of the 2 crustacean hemocyanin subunits, the 5 cheliceratan hemocyanin subunits, and the 12 insect hexamerins. Within the insect clade, four methionine- rich proteins, four arylphorins, and two juvenile hormone-suppressible proteins from Lepidoptera, as well as two dipteran proteins, form four separate groups. In the absence of an outgroup sequence, it is not possible to present information about the ancestral state from which these proteins are derived. Although this family of proteins clearly consists of homologous gene products, there remain striking differences in gene organization and site of biosynthesis of the proteins within the cell. Because studies on 18S and 12S rRNA sequences indicate a rather close relationship between insects and crustaceans, we propose that hemocyanin is the ancestral arthropod protein and that insect hexamerins lost their copper-binding capability after divergence of the insects from the crustaceans.      

High-mobility-group family genes from Lampetra japonica reveal their early origin and molecular evolution in the vertebrate lineage

High-mobility group family (HMG) genes are ubiquitous in vertebrates, including mammals, birds, amphibians and fishes. To elucidate the molecular phylogeny of the HMG genes in the primitive vertebrate, we have cloned three homologues of HMG-box genes, called Lj-HMGB1, Lj-HMGB2 and Lj-HMGBX, from a cDNA library generated from lymphocyte-like cells of the Japanese lamprey (Lampetra japonica), an Agnathan that occupies a critical phylogenetic position between invertebrates and vertebrates. The open reading frames of Lj-HMGB1, Lj-HMGB2 and Lj-HMGBX contained 627 bp, 585 bp and 678 bp, respectively. The analysis of the deduced amino acid sequences indicated that these three putative Lj-HMGB proteins contain four domains: HMG-box A, HMG-box B, an acidic carboxyl-terminal tail and a linker. A phylogenetic analysis revealed that the Lj-HMGB proteins fall outside the vertebrate clade; Lj-HMGBX is descended from a gene ancestral to the mammalian HMGB1/2/3. This discovery implies that there was a gene duplication event in the HMGB1/2/3 gene family that occurred after the divergence of the vertebrates (Cyclostomata) from the Cephalochordata and Urochordata at least 450 million years ago (MYA). The Lj-HMGB1, Lj-HMGB2 and Lj-HMGBX genes were detected in most tissues of the lamprey by RT-PCR. Our findings provide insight into the phylogeny of the HMGB genes in vertebrates.     

Molecular reconstruction and homology modelling of the catalytic domain of the common ancestor of the haemostatic vitamin-K-dependent serine proteinases

The vitamin-K-dependent serine proteinases of coagulation have evolved by a process of gene duplication and divergence, acquiring along the way a considerable degree of functional diversity that has equipped them for their different roles in haemostasis. The cDNA sequences encoding the catalytic domains of the early mammalian ancestors of five vitamin-K-dependent factors (factors VII, IX and X, protein C and prothrombin) were reconstructed by employing cDNA sequence data from a range of extant mammals and by using established phylogenies. The cDNA sequence of the putative common ancestor of these early mammalian proteins was then reconstructed from the five sequences by using a deduced phylogeny that was different in a number of respects from those previously proposed. Factor IX is proposed to have branched off early on, followed by protein C and prothrombin and finally factors VII and X. Significant differences in mutation rates were observed between proteins within a species; factor IX exhibited a lower mutation rate than the other proteins, consistent with its early emergence. Differences in mutation rates were also observed between species for a given protein and these exhibited an inverse correlation with generation time. A biophysically plausible structure for the ancestral vitamin-K-dependent factor protein was constructed by comparative methods. Studies of the functional architecture of this model provide new insights into the evolution of protein-binding specificity in this family of proteins. Received: 3 April 1996 / Revised: 24 April 1996     

Comparison of the beta-like globin gene families of rabbits and humans indicates that the gene cluster 5''-epsilon-gamma-delta-beta-3'' predates the mammalian radiation

The members of the rabbit and human beta-like globin gene families have been compared both by a computer-generated dot matrix graphical analysis of each entire gene and by calculating divergences in the coding regions. The rabbit-human gene pairs beta 4-epsilon, beta 3- gamma, psi beta 2-delta, and beta 1-beta were identified as orthologous on the basis of sequence similarities found in flanking and intervening sequences as well as by quantitative divergence calculations. The orthologous genes are in the same order on the chromosome in each species, which suggests that an ancestral family with the arrangement 5'-epsilon-gamma-delta-beta-3' preceded the mammalian radiation. Descendants of ancestral epsilon have diverged more slowly than other beta-like genes and are expressed only in embryonic life. Descendants of ancestral gamma and beta diverged at a higher rate and are expressed at wider range of developmental times. Descendants of delta have undergone nonreciprocal recombination at a high frequency and are often pseudogenes. Paralogous comparisons among the rabbit beta-like globin genes show that the beta 4-beta 3 and psi beta 2-beta 1 pairs are most similar and that beta 4 and beta 3 are more closely related to beta 1 than to psi beta 2. This fits with a branching pattern where the primordial beta split into ancestral epsilon/gamma and delta/beta genes, which later split into epsilon and gamma or delta and beta, respectively. Rabbit genes beta 4 and beta 1 acquired similar 3' untranslated regions after the epsilon/gamma split but prior to the mammalian radiation, presumably via a gene conversion event. The 5' end of beta 2 apparently converted with beta 1 after the radiation, and afterward it became a pseudogene.      

Evolutionary Analysis of the ErbB Receptor and Ligand Families

We have compared all available deduced protein sequences of the ErbB family of receptors and their ligands. Analysis of the aligned sequences of the receptors indicates that there are some differences in the receptors that are specific to invertebrates. In addition, comparison of the vertebrate ErbB receptors suggest that a gene duplication event generated two ancestral receptors, the ErbB3/ErbB4 precursor and the ErbB1/ErbB2 precursor. Subsequent gene duplications of these precursors generated the four receptors present in mammals. Analysis of the sequences for the known ligands of the ErbB receptors suggests that the vertebrate ligands segregate into the ErbB1 ligands and the ErbB3/ErbB4 ligands, paralleling the evolution of the receptors; however, it is difficult to ascertain any correlation between the invertebrate and the vertebrate ligands. Even though ErbB3 is kinase-impaired, there is significant conservation of the kinase domain within the vertebrate lineage (human, rat, and F. rubripes), suggesting some function for this domain other than kinase activity, such as mediating protein–protein interactions that are involved in receptor dimerization and/or activation of the kinase domain of the heterodimerization partner. To date, no ligand for ErbB2 has been identified, and comparison of the extracellular domains of ErbB2 reveals two regions that are not conserved across the mammalian species. These two regions of divergence align with sequences in ErbB1 that have been shown to be proximal to the amino-terminus and to the carboxyl-terminal region, respectively, of bound EGF. Further, one of these regions contains an insertion, relative to the other members of the mammalian ErbB family, which might affect the ligand binding site and provide a structural basis for this receptor's apparent inability to bind ligand independently. Received: 8 September 1999 / Accepted: 17 January 2000     

Evolutionary implications of the family of 14-3-3 brain protein homologs in Arabidopsis thaliana

The GF14 family of proteins in Arabidopsis thaliana consists of a homologous group of polypeptides ranging in size from 27 kDa to 32 kDa. As a group, GF14 proteins are also homologous to a family of mammalian proteins most commonly referred to as 14-3-3 proteins. Several distinct and different biochemical activities have been historically attributed to the various isoforms of the mammalian 14-3-3 proteins. These data present the possibility that the various activities are performed by functionally distinct lineages of the gene family. Here we present phylogenetic analyses based on the derived amino acid sequences of five GF14 isoforms expressed in Arabidopsis suspension-cultured cells. A high degree of sequence integrity is apparent in the various Arabidopsis isoforms, and the overall structures of the plant forms are quite conserved with regard to the structures of the known mammalian forms. These gene phylogenies indicate no evolutionary conservation of specific isoform lineages within both plants and animals. Rather, the evolutionary history of this protein appears to be characterized by a separate radiation of plant and animal forms from a common ancestral sequence. Even though the plant and animal forms have evolved independently since that ancestral split, large domains are conserved in both major lineages.     

Evolution of mitochondrial uncoupling proteins: novel invertebrate UCP homologues suggest early evolutionary divergence of the UCP family

Current hypothesis about the evolution of uncoupling proteins (UCPs) proposed by suggests that UCP4 is the earliest form of UCP ancestral to all other UCP orthologues. However, this hypothesis is difficult to reconcile with a narrow tissue distribution of UCP4 (which is a brain-specific isoform), suggesting highly specialized rather than anfcestral function for this protein. We searched for UCP2, UCP3, and UCP5 homologues in invertebrate genomes using amplification with degenerate primers designed against UCP2-specific conserved sequences and/or BLASTP search with stringent ad hoc criteria to distinguish between homologues and orthologues of different UCPs. Our study identified invertebrate UCP homologues similar to UCP2 and 3 (which we termed UCP6) and an invertebrate homologue of UCP5. Phylogenetic analysis indicates that there are at least three clades of UCPs in invertebrates, which are closely related to vertebrate UCP1-3, UCP4, and UCP5, respectively, and shows early evolutionary divergence of UCPs, which pre-dates the divergence of protostomes and deuterostomes. It also suggests that the newly identified UCP6 proteins from invertebrates are ancestral to the vertebrate UCP1, UCP2, and UCP3, and that divergence of these three vertebrate orthologues occurred late in evolution of the vertebrates. This study refutes the hypothesis of Hanak and Jezek (2001) that UCP4 is an ancestral form for all UCPs, and shows early evolutionary diversification of this protein family, which corresponds to their proposed functional diversity in regulation of proton leak, antioxidant defense and apoptosis.     

Characterization of multiple members of the HSP70 family in platyfish culture cells: molecular evolution of stress protein HSP70 in vertebrates

A shift from 28 to 37 degrees C in the incubation temperature of a culture of the platyfish fibroblast cell line, EHS cells (platyfish fibroblast cell line), induced a set of stress proteins. A two-dimensional polyacrylamide gel electrophoresis (2D-PAGE) analysis showed that the cells expressed three genetically distinct forms of heat-shock protein 70 (HSP70) family proteins: heat-inducible forms of HSP70, the constitutively expressed heat-shock cognate protein 70 (HSC70) and its phosphorylated isoform, and the glucose-regulated protein 78 (GRP78). Three different clones encoding two major isoforms of heat-inducible HSP70, platyfish HSP70-1 and HSP70-2, and of the HSC70 were isolated from a platyfish cDNA library. We compared the deduced amino acid sequences of the platyfish HSP70 and HSC70 proteins with those of other vertebrates. Phylogenetic analysis showed that vertebrate HSP70 could be classified into four cluster groups: (a) fish HSP70, with two isoforms of heat-inducible HSP70 in fish, fish HSP70-1 and HSP70-2; (b) the mammalian testis-specific HSP70-related protein HST70; (c) the mammalian heat-inducible HSP70B'; and (d) the mammalian major histocompatibility complex (MHC)-linked HSP70, including the MHC-linked heat-inducible HSP70 and the testis-specific HSP70-related protein. These findings suggest that vertebrate HSP70 was derived from a single ancestral HSP70 gene during vertebrate evolution and that multiple copies of heat-inducible HSP70 were probably evolved during genetic divergence in fish and higher vertebrates.     

Differential evolution of voltage-gated sodium channels in tetrapods and teleost fishes

The voltage-gated sodium channel (SCN) alpha subunits are large proteins with central roles in the generation of action potentials. They consist of approximately 2,000 amino acids encoded by 24-27 exons. Previous evolutionary studies have been unable to reconcile the proposed gene duplication schemes with the species distribution and molecular phylogeny of the genes. We have carefully annotated the complete SCN gene sequences, correcting numerous database errors, for a broad range of vertebrate species and analyzed their phylogenetic relationships. We have also compared the chromosomal positions of the SCN genes relative to adjacent gene families. Our studies show that the ancestor of the vertebrates probably had a single sodium channel gene with two characteristic AT-AC introns, the second of which is unique to vertebrate SCN genes. This ancestral gene, located close to a HOX gene cluster, was quadrupled along with HOX in the two rounds of basal vertebrate tetraploidizations to generate the ancestors of the four channels SCN1A, SCN4A, SCN5A, and SCN8A. The third tetraploidization in the teleost fish ancestor doubled this set of genes and all eight are still present in at least three of four investigated teleost fish genomes. In tetrapods, the gene family expanded by local duplications before the radiation of amniotes, generating the cluster SCN5A, SCN10A, and SCN11A on one chromosome and the cluster SCN1A, SCN2A, SCN3A, and SCN9A on a different chromosome. In eutherian mammals, a tenth gene, SCN7A, arose in a local duplication in the SCN1A gene cluster. The SCN7A gene has undergone rapid evolution and has lost the ability to cause action potentials-instead, it functions as a sodium sensor. The three genes in the SCN5A cluster were translocated from the HOX-bearing chromosome in a mammalian ancestor along with several adjacent genes. This evolutionary scenario is supported by the adjacent TGF-β receptor superfamily (comprised of five distinct families) and the cysteine-serine-rich nuclear protein gene family as well as the HOX clusters. The independent expansions of the SCN repertoires in tetrapods and teleosts suggest that the functional diversification may differ between the two lineages.     

Evolution and diversity of axon guidance Robo receptor family genes

The diversity of axon guidance (AG) receptors reflects gains in complexity of the animal nervous system during evolution. Members of the Roundabout (Robo) family of receptors interact with Slit proteins and play important roles in many developmental processes, including AG and neural crest cell migration. There are four members of the Robo gene family. However, the evolutionary history of Robo family genes remain obscure. We analyzed the distribution of Robo family members in metazoan species ranging in complexity from hydras to humans. We undertook a phylogenetic analysis in metazoans, synteny analysis, and ancestral chromosome mapping in vertebrates, and detected selection pressure and functional divergence among four mammalian Robo paralogs. Based on our analysis, we proposed that the ancestral Robo gene could have undergone a tandem duplication in the vertebrate ancestor; then one round of whole genome duplication events occurred before the divergence of ancestral lamprey and gnathostome, generating four paralogs in early vertebrates. Robo4 paralog underwent segmental loss in the following evolutionary process. Our results showed that Robo3 paralog is under more powerful purifying selection pressure compared with other three paralogs, which could correlate with its unique expression pattern and function. Furthermore, we found four sites under positive selection pressure on the Ig1‐2 domains of Robo4 that might interfere with its binding to Slits ligand. Diverge analysis at the amino acid level showed that Robo4 paralog have relatively greater functional diversifications than other Robo paralogs. This coincides with the fact that Robo4 predominantly functions in vascular endothelial cells but not the nervous system.     

Three novel mammalian toll-like receptors: gene structure, expression, and evolution

We describe three novel genes, encoding members of the Toll-like receptor (Tlr) family (TLR7, TLR8, and TLR9). These Tlr family members, unlike others reported to date, were identified within a genomic database. TLR7 and TLR8 each have three exons, two of which have coding function, and lie in close proximity to one another at Xp22, alongside a pseudogene. The remaining gene (TLR9) resides at 3p21.3 (in linkage with the MyD88 gene), and is expressed in at least two splice forms, one of which is monoexonic and one of which is biexonic, the latter encoding a protein with 57 additional amino acids at the N-terminus. The novel Tlrs comprise a cluster as nearest phylogenetic neighbors. Combining all sequence data related to Toll-like receptors, we have drawn several inferences concerning the phylogeny of vertebrate and invertebrate Tlrs. According to our best estimates, mammalian TLRs 1 and 6 diverged from a common mammalian ancestral gene 95 million years ago. TLR4, which encodes the endotoxin sensor in present-day mammals, emerged as a distinct entity 180 million years ago. TLRs 3 and 5 diverged from a common ancestral gene approximately 150 million years ago, as did Tlr7 and Tlr8. Very likely, fewer Tlrs existed during early vertebrate evolution: at most three or four were transmitted with the primordial vertebrate line. Phylogenetic data that we have adduced in the course of this work also suggest the existence of a Drosophila equivalent of MyD88, and indicate that the plasma membrane protein SIGIRR is close functional relative of MyD88 in mammals. Finally, a single present-day representative of the Toll-like proteins in Drosophila has striking cytoplasmic domain homology to mammalian Tlrs within the cluster that embraces TLRs 1, 2, 4, and 6. This would suggest that an ancestral (pre-vertebrate) Tlr may have adopted a pro-inflammatory function 500 million years ago.     

Comparative analysis of methylthioalkylmalate synthase (MAM) gene family and flanking DNA sequences in Brassica oleracea and Arabidopsis thaliana

Gene BoGSL-PRO is associated with presence of 3-carbon side-chain glucosinolates (GSL). This gene is a member of the methylthioalkylmalate synthase (MAM) gene family. A BAC clone of Brassica oleracea, B21F5, containing this gene, was sequenced, annotated and compared to its corresponding region in Arabidopsis thaliana. Twelve protein-coding genes and 10 transposable elements were found in this clone. The corresponding region in A. thaliana chromosome I has 14 genes and no transposable elements. Analysis of MAM gene family in both species, which also include genes controlling 4-carbon side-chain GSL, separated the genes in two groups based on exon numbers and function. Phylogenetic analysis of the amino acid sequences encoded by these genes suggest that these two groups were produced by a duplication that must have occurred before the divergence of the Rosid and Asterid lineages of angiosperms. Comparison with putative orthologs from several prokaryotes further suggest that the members of the gene family with 10 exons, which encode proteins involved in 4-carbon side-chain GSL biosynthesis, were derived via truncation of the 3′ end from ancestral genes more similar in length to those with 12 exons, which encode proteins involved in 3-carbon side-chain GSL biosynthesis. Lower gene density in B. oleracea compared to A. thaliana is due in part to presence of transposable elements (TE) mostly in inter-genic regions.