Extending the family table: Insights from beyond vertebrates into the regulation of embryonic development by FGFs

Since the discovery of fibroblast growth factors (FGFs) much focus has been placed on elucidating the roles for each vertebrate FGF ligand, receptor, and regulating molecules in the context of vertebrate development, human disorders and cancer. Studies in human, mouse, frog, chick, and zebrafish have made great contributions to our understanding of the role of FGFs in specific processes. However, in recent years, as more genomes are sequenced, information is becoming available from many non‐vertebrate models and a more complete picture of the FGF superfamily as a whole is emerging. In some cases, less redundancy in these FGF signaling systems may allow for more mechanistic insights. Studies in sea anemones have highlighted how ancient FGF signaling is and helped provide insight into the evolution of the FGF gene family. Work in nematodes has shown that different splice forms can be used for functional specificity in invertebrate FGF signaling. Comparing FGFs between urochordates and vertebrates as well as between different insect species reveals important clues into the process of gene loss, duplication and subfunctionalization of FGFs throughout evolution. Finally, comparing all members of the FGF ligand superfamily reveals variability in many properties, which may point to a feature of FGFs as being highly adaptable with regards to protein structure and signaling mechanism. Further studies on FGF signaling outside of vertebrates is likely to continue to complement work in vertebrates by contributing additional insights to the FGF field and providing unexpected information that could be used for medical applications. Birth Defects Research (Part C) 90:214–227, 2010. © 2010 Wiley‐Liss, Inc.     

Of Worms and Men: An Evolutionary Perspective on the Fibroblast Growth Factor (FGF) and FGF Receptor Families

FGFs (fibroblast growth factors) play major roles in a number of developmental processes. Recent studies of several human disorders, and concurrent analysis of gene knock-out and properties of the corresponding recombinant proteins have shown that FGFs and their receptors are prominently involved in the development of the skeletal system in mammals. We have compared the sequences of the nine known mammalian FGFs, FGFs from other vertebrates, and three additional sequences that we extracted from existing databases: two human FGF sequences that we tentatively designated FGF10 and FGF11, and an FGF sequence from C?norhabditis elegans. Similarly, we have compared the sequences of the four FGF receptor paralogs found in chordates with four non-chordate FGF receptors, including one recently identified in C. elegans. The comparison of FGF and FGF receptor sequences in vertebrates and nonvertebrates shows that the FGF and FGF receptor families have evolved through phases of gene duplications, one of which may have coincided with the emergence of vertebrates, in relation with their new system of body scaffold. Received: 6 April 1996 / Accepted: 5 July 1996     

The vestigial olfactory receptor subgenome of odontocete whales: phylogenetic congruence between gene-tree reconciliation and supermatrix methods

The macroevolutionary transition of whales (cetaceans) from a terrestrial quadruped to an obligate aquatic form involved major changes in sensory abilities. Compared to terrestrial mammals, the olfactory system of baleen whales is dramatically reduced, and in toothed whales is completely absent. We sampled the olfactory receptor (OR) subgenomes of eight cetacean species from four families. A multigene tree of 115 newly characterized OR sequences from these eight species and published data for Bos taurus revealed a diverse array of class II OR paralogues in Cetacea. Evolution of the OR gene superfamily in toothed whales (Odontoceti) featured a multitude of independent pseudogenization events, supporting anatomical evidence that odontocetes have lost their olfactory sense. We explored the phylogenetic utility of OR pseudogenes in Cetacea, concentrating on delphinids (oceanic dolphins), the product of a rapid evolutionary radiation that has been difficult to resolve in previous studies of mitochondrial DNA sequences. Phylogenetic analyses of OR pseudogenes using both gene-tree reconciliation and supermatrix methods yielded fully resolved, consistently supported relationships among members of four delphinid subfamilies. Alternative minimizations of gene duplications, gene duplications plus gene losses, deep coalescence events, and nucleotide substitutions plus indels returned highly congruent phylogenetic hypotheses. Novel DNA sequence data for six single-copy nuclear loci and three mitochondrial genes (> 5000 aligned nucleotides) provided an independent test of the OR trees. Nucleotide substitutions and indels in OR pseudogenes showed a very low degree of homoplasy in comparison to mitochondrial DNA and, on average, provided more variation than single-copy nuclear DNA. Our results suggest that phylogenetic analysis of the large OR superfamily will be effective for resolving relationships within Cetacea whether supermatrix or gene-tree reconciliation procedures are used.     

Pinaceae show elevated rates of gene turnover that are robust to incomplete gene annotation

Gene duplications and gene losses are major determinants of genome evolution and phenotypic diversity. The frequency of gene turnover (gene gains and gene losses combined) is known to vary between organisms. Comparative genomic analyses of gene families can highlight such variation; however, estimates of gene turnover may be biased when using highly fragmented genome assemblies resulting in poor gene annotations. Here, we address potential biases introduced by gene annotation errors in estimates of gene turnover frequencies in a dataset including both well‐annotated angiosperm genomes and the incomplete gene sets of four Pinaceae, including two pine species, Norway spruce and Douglas‐fir. We show that Pinaceae experienced higher gene turnover rates than angiosperm lineages lacking recent whole‐genome duplications. This finding is robust to both known major issues in Pinaceae gene sets: missing gene models and erroneous annotation of pseudogenes. A separate analysis limited to the four Pinaceae gene sets pointed to an accelerated gene turnover rate in pines compared with Norway spruce and Douglas‐fir. Our results indicate that gene turnover significantly contributes to genome variation and possibly to speciation in Pinaceae, particularly in pines. Moreover, these findings indicate that reliable estimates of gene turnover frequencies can be discerned in incomplete and potentially inaccurate gene sets. Because gymnosperms are known to exhibit low overall substitution rates compared with angiosperms, our results suggest that the rate of single‐base pair mutations is uncoupled from the rate of large DNA duplications and deletions associated with gene turnover in Pinaceae.     

Heparin binding preference and structures in the fibroblast growth factor family parallel their evolutionary diversification

The interaction of a large number of extracellular proteins with heparan sulfate (HS) regulates their transport and effector functions, but the degree of molecular specificity underlying protein–polysaccharide binding is still debated. The 15 paracrine fibroblast growth factors (FGFs) are one of the paradigms for this interaction. Here, we measure the binding preferences of six FGFs (FGF3, FGF4, FGF6, FGF10, FGF17, FGF20) for a library of modified heparins, representing structures in HS, and model glycosaminoglycans, using differential scanning fluorimetry. This is complemented by the identification of the lysine residues in the primary and secondary binding sites of the FGFs by a selective labelling approach. Pooling these data with previous sets provides good coverage of the FGF phylogenetic tree, deduced from amino acid sequence alignment. This demonstrates that the selectivity of the FGFs for binding structures in sulfated polysaccharides and the pattern of secondary binding sites on the surface of FGFs follow the phylogenetic relationship of the FGFs, and so are likely to be the result of the natural selection pressures that led to the expansion of the FGF family in the course of the evolution of more complex animal body plans.     

Genome-level and biochemical diversity of the acyl-activating enzyme superfamily in plants

In higher plants, the superfamily of carboxyl-CoA ligases and related proteins, collectively called acyl activating enzymes (AAEs), has evolved to provide enzymes for many pathways of primary and secondary metabolism and for the conjugation of hormones to amino acids. Across the superfamily there is only limited sequence similarity, but a series of highly conserved motifs, including the AMP-binding domain, make it easy to identify members. These conserved motifs are best understood in terms of the unique domain-rotation architecture that allows AAE enzymes to catalyze the two distinct steps of the CoA ligase reaction. Arabidopsis AAE sequences were used to identify the AAE gene families in the sequenced genomes of green algae, mosses, and trees; the size of the respective families increased with increasing degree of organismal cellular complexity, size, and generation time. Large-scale genome duplications and small-scale tandem gene duplications have contributed to AAE gene family complexity to differing extents in each of the multicellular species analyzed. Gene duplication and evolution of novel functions in Arabidopsis appears to have occurred rapidly, because acquisition of new substrate specificity is relatively easy in this class of proteins. Convergent evolution has also occurred between members of distantly related clades. These features of the AAE superfamily make it difficult to use homology searches and other genomics tools to predict enzyme function.     

Analysis of lamprey and hagfish genes reveals a complex history of gene duplications during early vertebrate evolution

It has been proposed that two events of duplication of the entire genome occurred early in vertebrate history (2R hypothesis). Several phylogenetic studies with a few gene families (mostly Hox genes and proteins from the MHC) have tried to confirm these polyploidization events. However, data from a single locus cannot explain the evolutionary history of a complete genome. To study this 2R hypothesis, we have taken advantage of the phylogenetic position of the lamprey to study the history of gene duplications in vertebrates. We selected most gene families that contain several paralogous genes in vertebrates and for which lamprey genes and an out-group are known in databases. In addition, we isolated members of the nuclear receptor superfamily in lamprey. Hagfish genes were also analyzed and found to confirm the lamprey gene analysis. Consistent with the 2R hypothesis, the phylogenetic analysis of 33 selected gene families, dispersed through the whole genome, revealed that one period of gene duplication arose before the lamprey-gnathostome split and this was followed by a second period of gene duplication after the lamprey-gnathostome split. Nevertheless, our analysis suggests that numerous gene losses and other gene-genome duplications occurred during the evolution of the vertebrate genomes. Thus, the complexity of all the paralogy groups present in vertebrates should be explained by the contribution of genome duplications (2R hypothesis), extra gene duplications, and gene losses.     

Functional phylogeny relates LET-756 to fibroblast growth factor 9

Fibroblast growth factors (FGFs) are secreted regulatory proteins involved in various developmental processes. In vertebrates, the FGF superfamily comprises 22 members. In non-vertebrates, six FGF genes have been identified in Ciona intestinalis, three in Drosophila melanogaster, and two (let-756 and egl-17) in Caenorhabditis elegans. The core of LET-756 shares a 30-50% sequence identity with the various members of the superfamily. The relationships between vertebrate and non-vertebrate FGFs are not clear. We made chimeric FGFs by replacing the core region of LET-756 by the cores of various mammalian, fly, and worm FGFs. LET-756 deleted in its core region was no longer able to rescue the lethal phenotype of a let-756 null mutant, and only chimeras containing the cores of FGFs 9, 16, and 20 showed rescue capacity. This core contains an internal motif of six amino acid residues (EFISIA) whose deletion or mutation abolished both the rescue activity and FGF secretion in the supernatant of transfected COS-1 cells. Chimera containing the core of C. intestinalis FGF9/16/20, a potential ortholog of FGF9 lacking the complete EFISIA motif, was not able to rescue the lethal phenotype or be secreted. However, the introduction of the EFISIA motif restored both activities. The data show that the EFISIA motif in the core of LET-756 is essential for its biological activity and that FGFs 9, 16, and 20, which contain that motif, are functionally close to LET-756 and may be evolutionary related. This non-classical mode of secretion using an internal motif is conserved throughout evolution.     

Duplicate retention in signalling proteins and constraints from network dynamics

Duplications are a major driving force behind evolution. Most duplicates are believed to fix through genetic drift, but it is not clear whether this process affects all duplications equally or whether there are certain gene families that are expected to show neutral expansions under certain circumstances. Here, we analyse the neutrality of duplications in different functional classes of signalling proteins based on their effects on response dynamics. We find that duplications involving intermediary proteins in a signalling network are neutral more often than those involving receptors. Although the fraction of neutral duplications in all functional classes increase with decreasing population size and selective pressure on dynamics, this effect is most pronounced for receptors, indicating a possible expansion of receptors in species with small population size. In line with such an expectation, we found a statistically significant increase in the number of receptors as a fraction of genome size in eukaryotes compared with prokaryotes. Although not confirmative, these results indicate that neutral processes can be a significant factor in shaping signalling networks and affect proteins from different functional classes differently.     

Reconstructing ancestral genome content based on symmetrical best alignments and Dollo parsimony

MOTIVATION: Gene duplications and losses (GDLs) are important events in genome evolution. They result in expansion or contraction of gene families, with a likely role in phenotypic evolution. As more genomes become available and their annotations are improved, software programs capable of rapidly and accurately identifying the content of ancestral genomes and the timings of GDLs become necessary to understand the unique evolution of each lineage. RESULTS: We report EvolMAP, a new algorithm and software that utilizes a species tree-based gene clustering method to join all-to-all symmetrical similarity comparisons of multiple gene sets in order to infer the gene composition of multiple ancestral genomes. The algorithm further uses Dollo parsimony-based comparison of the inferred ancestral genes to pinpoint the timings of GDLs onto evolutionary intervals marked by speciation events. Using EvolMAP, first we analyzed the expansion of four families of G-protein coupled receptors (GPCRs) within animal lineages. Additional to demonstrating the unique expansion tree for each family, results also show that the ancestral eumetazoan genome contained many fewer GPCRs than modern animals, and these families expanded through concurrent lineage-specific duplications. Second, we analyzed the history of GDLs in mammalian genomes by comparing seven proteomes. In agreement with previous studies, we report that the mammalian gene family sizes have changed drastically through their evolution. Interestingly, although we identified a potential source of duplication for 75% of the gained genes, remaining 25% did not have clear-cut sources, revealing thousands of genes that have likely gained their distinct sequence identities within the descent of mammals. AVAILABILITY: Query server, source code and executable are available at http://kosik-web.mcdb.ucsb.edu/evolmap/index.htm .     

Prevalence of intron gain over intron loss in the evolution of paralogous gene families

The mechanisms and evolutionary dynamics of intron insertion and loss in eukaryotic genes remain poorly understood. Reconstruction of parsimonious scenarios of gene structure evolution in paralogous gene families in animals and plants revealed numerous gains and losses of introns. In all analyzed lineages, the number of acquired new introns was substantially greater than the number of lost ancestral introns. This trend held even for lineages in which vertical evolution of genes involved more intron losses than gains, suggesting that gene duplication boosts intron insertion. However, dating gene duplications and the associated intron gains and losses based on the molecular clock assumption showed that very few, if any, introns were gained during the last ~100 million years of animal and plant evolution, in agreement with previous conclusions reached through analysis of orthologous gene sets. These results are generally compatible with the emerging notion of intensive insertion and loss of introns during transitional epochs in contrast to the relative quiet of the intervening evolutionary spans.     

Diversity, classification and function of the plant protein kinase superfamily

Eukaryotic protein kinases belong to a large superfamily with hundreds to thousands of copies and are components of essentially all cellular functions. The goals of this study are to classify protein kinases from 25 plant species and to assess their evolutionary history in conjunction with consideration of their molecular functions. The protein kinase superfamily has expanded in the flowering plant lineage, in part through recent duplications. As a result, the flowering plant protein kinase repertoire, or kinome, is in general significantly larger than other eukaryotes, ranging in size from 600 to 2500 members. This large variation in kinome size is mainly due to the expansion and contraction of a few families, particularly the receptor-like kinase/Pelle family. A number of protein kinases reside in highly conserved, low copy number families and often play broadly conserved regulatory roles in metabolism and cell division, although functions of plant homologues have often diverged from their metazoan counterparts. Members of expanded plant kinase families often have roles in plant-specific processes and some may have contributed to adaptive evolution. Nonetheless, non-adaptive explanations, such as kinase duplicate subfunctionalization and insufficient time for pseudogenization, may also contribute to the large number of seemingly functional protein kinases in plants.     

Eleven ancestral gene families lost in mammals and vertebrates while otherwise universally conserved in animals

Background  

Gene losses played a role which may have been as important as gene and genome duplications and rearrangements, in modelling today species' genomes from a common ancestral set of genes. The set and diversity of protein-coding genes in a species has direct output at the functional level. While gene losses have been reported in all the major lineages of the metazoan tree of life, none have proposed a focus on specific losses in the vertebrates and mammals lineages. In contrast, genes lost in protostomes (i.e. arthropods and nematodes) but still present in vertebrates have been reported and extensively detailed. This probable over-anthropocentric way of comparing genomes does not consider as an important phenomena, gene losses in species that are usually described as "higher". However reporting universally conserved genes throughout evolution that have recently been lost in vertebrates and mammals could reveal interesting features about the evolution of our genome, particularly if these losses can be related to losses of capability.     

Two patterns of genome organization in mammals: the chromosomal distribution of duplicate genes in human and mouse

Gene duplication occurs repeatedly in the evolution of genomes, and the rearrangement of genomic segments has also occurred repeatedly over the evolution of eukaryotes. We studied the interaction of these two factors in mammalian evolution by comparing the chromosomal distribution of multigene families in human and mouse. In both species, gene families tended to be confined to a single chromosome to a greater extent than expected by chance. The average number of families shared between chromosomes was nearly 60% higher in mouse than in human, and human chromosomes rarely shared large numbers of gene families with more than one or two other chromosomes, whereas mouse chromosomes frequently did so. A higher proportion of duplicate gene pairs on the same chromosome originated from recent duplications in human than in mouse, whereas a higher proportion of duplicate gene pairs on separate chromosomes arose from ancient duplications in human than in mouse. These observations are most easily explained by the hypotheses that (1) most gene duplications arise in tandem and are subsequently separated by segmental rearrangement events, and (2) that the process of segmental rearrangement has occurred at a higher rate in the lineage of mouse than in that of human.     

Choosing among alternative trees of multigene families

Estimation of gene trees is the first step in testing alternative hypotheses about the evolution of multigene families. The standard practice for inferring gene family history is to construct trees that meet some objective criteria based on the fit of the character state changes (nucleotide or amino acid changes) to the gene tree. Unfortunately, analysis of character state data can be misleading. In addition, this approach ignores information about the relationships of the species from which the genes have been sampled. In this paper I explore using statistics of fit between the character data and gene trees and the reconciliation of the gene and species trees for choosing among alternative evolutionary hypotheses of gene families. In particular, I advocate a two-pronged strategy for choosing among alternative gene trees. First, the character data are used to define a set of acceptable gene trees (i.e., trees that are not significantly different from the minimum length tree). Next, the set of acceptable gene trees is reconciled with a known species tree, and the gene tree requiring the fewest number of gene duplications and losses is adopted as the best estimate of evolutionary history. The approach is illustrated using three gene families: BMP, EGR, and LDH.     

Lineage-specific expansions provide genomic complexity among sea urchin GTPases

In every organism, GTP-binding proteins control many aspects of cell signaling. Here, we examine in silico several GTPase families from the Strongylocentrotus purpuratus genome: the monomeric Ras superfamily, the heterotrimeric G proteins, the dynamin superfamily, the SRP/SR family, and the "protein biosynthesis" translational GTPases. Identified were 174 GTPases, of which over 90% are expressed in the embryo as shown by tiling array and expressed sequence tag data. Phylogenomic comparisons restricted to Drosophila, Ciona, and humans (protostomes, urochordates, and vertebrates, respectively) revealed both common and unique elements in the expected composition of these families. Galpha and dynamin families contain vertebrate expansions, consistent with whole genome duplications, whereas SRP/SR and translational GTPases are highly conserved. Unexpectedly, Ras superfamily analyses revealed several large (5+) lineage-specific expansions in the sea urchin. For Rho, Rab, Arf, and Ras subfamilies, comparing total human gene numbers to the number of sea urchin genes with vertebrate orthologs suggests reduced genomic complexity in the sea urchin. However, gene duplications in the sea urchin increase overall numbers such that total sea urchin gene numbers approximate vertebrate gene numbers for each monomeric GTPase family. These findings suggest that lineage-specific expansions may be an important component of genomic evolution in signal transduction.     

Orthologs, paralogs and genome comparisons

During the past decade, ancient gene duplications were recognized as one of the main forces in the generation of diverse gene families and the creation of new functional capabilities. New tools developed to search data banks for homologous sequences, and an increased availability of reliable three-dimensional structural information led to the recognition that proteins with diverse functions can belong to the same superfamily. Analyses of the evolution of these superfamilies promises to provide insights into early evolution but are complicated by several important evolutionary processes. Horizontal transfer of genes can lead to a vertical spread of innovations among organisms, therefore finding a certain property in some descendants of an ancestor does not guarantee that it was present in that ancestor. Complete or partial gene conversion between duplicated genes can yield phylogenetic trees with several, apparently independent gene duplications, suggesting an often surprising parallelism in the evolution of independent lineages. Additionally, the breakup of domains within a protein and the fusion of domains into multifunctional proteins makes the delineation of superfamilies a task that remains difficult to automate.     

Identification and expression profile analysis of the protein kinase gene superfamily in maize development

Eukaryotic protein kinases (ePKs) evolved as a family of highly dynamic molecular switches that serve to orchestrate the activity of almost all cellular processes. Some of the functionally characterized ePKs from plants have been found to be components of signaling networks, such as those for the perception of biotic agents, light quality and quantity, plant hormones, and various adverse environmental conditions. To date, only a tiny fraction of plant ePKs have been functionally identified, and even fewer have been identified in maize [Zea mays (Zm)]. In this study, we have identified 1,241 PK-encoding genes in the maize genome. Phylogenetic analyses identified eight gene groups with considerable conservation among groups, and each group could be further divided into multiple families and/or subfamilies. Similar intron/exon structural patterns were observed in the same families/subfamilies, strongly supporting their close evolutionary relationship. Chromosome distribution and genetic analysis revealed that tandem duplications and segmental/whole-genome duplications might represent two of the major mechanisms contributing to the expansion of the PK superfamily in maize. The dynamic expression patterns of ZmPK genes across the 60 different developmental stages of 11 organs showed that some members of this superfamily exhibit tissue-specific expression, whereas others are more ubiquitously expressed, indicative of their important roles in performing diverse developmental and physiological functions during the maize life cycle. Furthermore, RNA-sequence-based gene expression profiling of PKs along a leaf developmental gradient and in mature bundle sheath and mesophyll cells indicated that ZmPK genes are involved in various physiological processes, such as cell-fate decisions, photosynthetic differentiation, and regulation of stomatal development. Our results provide new insights into the function and evolution of maize PKs and will be useful in studies aimed at revealing the global regulatory network of maize development, thereby contributing to the maize molecular breeding with enhanced quality traits.     

啮齿目核糖酸酶5基因(RNase5)的分子进化

RNase5是RNASE A基因超家族中的一个重要成员,是分子进化研究的理想模型之一。基于基因组水平,我们对啮齿目的3个进化枝10科17个物种开展RNase5的分子进化研究。利用TBlastN及BlastN方法鉴定每个基因组的RNase5基因,发现该基因在啮齿目的Ctenohystrica所有物种发生丢失,时间是在Ctenohystrica形成之后;邻接法和最大似然法构建的系统发育树均支持RNase5在“与小家鼠相关的进化枝”的小家鼠、褐家鼠和拉布拉多白足鼠发生三次独立基因复制事件;利用PAML软件的枝模型、位点模型及枝-位点模型计算选择压力,均检测到RNase5基因受到强烈的正选择作用。总之,我们的研究深入系统开展了RNase5在啮齿目中的分子进化,增加了该基因研究的多样性,为进一步系统认识该基因在动物的适应性进化遗传机制奠定了基础。     

Evolution of cytochrome P-450 proteins [published erratum appears in Mol Biol Evol 1988 Mar;5(2):199]

Thirty-four cytochrome P-450 sequences from one bacterial and six vertebrate species have been aligned with the aid of a computer alignment algorithm. Phylogenetic trees were constructed using the unweighted-pair-group and neighbor-joining methods. The two trees differed at only a single branch point near the base of the tree. The cytochrome P-450 superfamily of proteins clustered into eight families and contained 16 gene-duplication events. The first gene duplication occurred approximately 1,360 Myr before the present (Mybp) and gave rise to cytochrome P-450s found in two different cellular organelles, the mitochondria and the endoplasmic reticulum. Both groups utilize cholesterol or its metabolites as substrates, implying that cholesterol existed greater than 1,360 Mybp. The fourth gene duplication (approximately 900 Mybp) gave rise to the drug-metabolizing P-450s. These proteins aid in the detoxification of foreign chemicals, as opposed to the metabolism of endogenous compounds. The importance of the capacity to metabolize drugs is reflected in 11 further gene duplications occurring in this lineage. The first occurred approximately 800 Mybp and gave rise to the two major P-450 families, the phenobarbital and 3-methylcholanthrene families. An apparent increase in the rate of cytochrome P-450 evolution is noted between the bird-mammal divergence (300 Mybp) and the mammalian radiation (75 Mybp).