Evolution of gene families and relationship with organismal evolution: rapid divergence of tissue-specific genes in the early evolution of chordates

To determine a possible relationship between organismal and molecular evolution, the divergence patterns of gene families were examined by taking special notice of functional difference, tissue distribution, and intracellular localization of the members. A phylogenetic analysis of 25 different gene families revealed interesting patterns of divergence of these families: Most gene duplications giving rise to different functions antedate the vertebrates-arthropods separation. On the other hand, in a group of members carrying virtually identical function to one another but differing in tissue distribution (tissue- specific isoform), most gene duplications have occurred independently in each of vertebrates and arthropods after the separation of the two animal groups. In family members encoding molecules localizing in cell compartments (compartmentalized isoforms), the gene duplications antedate the animals-fungi separation. In the cases of the Ca2+ pump and rab subfamilies, the compartmentalized isoforms were shown to have diverged during the early evolution of eukaryotes. A phylogenetic analysis of the tissue-specific isoforms from 26 different subfamilies revealed extensive gene duplications and rapid rates of amino acid substitutions in the early evolution of chordates before the separation of fishes and tetrapods. On the contrary, the genetic variations are relatively low in the later period. This pattern of evolution observed at the molecular level is correlated well with that of tissue evolution based on fossil evidence and morphological data, and thus evolution at the two levels may be related.      

Complex evolutionary history and diverse domain organization of SET proteins suggest divergent regulatory interactions

? Plants and animals possess very different developmental processes, yet share conserved epigenetic regulatory mechanisms, such as histone modifications. One of the most important forms of histone modification is methylation on lysine residues of the tails, carried out by members of the SET protein family, which are widespread in eukaryotes. ? We analyzed molecular evolution by comparative genomics and phylogenetics of the SET genes from plant and animal genomes, grouping SET genes into several subfamilies and uncovering numerous gene duplications, particularly in the Suv, Ash, Trx and E(z) subfamilies. ? Domain organizations differ between different subfamilies and between plant and animal SET proteins in some subfamilies, and support the grouping of SET genes into seven main subfamilies, suggesting that SET proteins have acquired distinctive regulatory interactions during evolution. We detected evidence for independent evolution of domain organization in different lineages, including recruitment of new domains following some duplications. ? More recent duplications in both vertebrates and land plants are probably the result of whole-genome or segmental duplications. The evolution of the SET gene family shows that gene duplications caused by segmental duplications and other mechanisms have probably contributed to the complexity of epigenetic regulation, providing insights into the evolution of the regulation of chromatin structure.     

Protein Tyrosine Phosphatases from Amphioxus, Hagfish, and Ray: Divergence of Tissue-Specific Isoform Genes in the Early Evolution of Vertebrates

Since separation from fungi and plants, multicellular animals evolved a variety of gene families involved in cell-cell communication from a limited number of ancestral precursors by gene duplications in two separate periods of animal evolution. In the very early evolution of animals before the separation of parazoans and eumetazoans, animals underwent extensive gene duplications by which different subtypes (subfamilies) with distinct functions diverged. The multiplicity of members (isoforms) in the same subtype increased by further gene duplications (isoform duplications) in the first half of chordate evolution before the fish-tetrapod split; different isoforms are virtually identical in structure and function but differ in tissue distribution. From cloning and phylogenetic analyses of four subfamilies of the protein tyrosine kinase (PTK) family, we recently showed extensive isoform duplications in a limited period around or just before the cyclostome-gnathostome split. To obtain a reliable estimate for the divergence time of vertebrate isoforms, we have conducted isolation of cDNAs encoding the protein tyrosine phosphatases (PTPs) from Branchiostoma belcheri, an amphioxus, Eptatretus burgeri, a hagfish, and Potamotrygon motoro, a ray. We obtained 33 different cDNAs in total, most of which belong to known PTP subfamilies. The phylogenetic analyses of five subfamilies based on the maximum likelihood method revealed frequent isoform duplications in a period around or just before the gnathostome-cyclostome split. An evolutionary implication was discussed in relation to the Cambrian explosion.     

Six family genes--structure and function as transcription factors and their roles in development

The members of the Six gene family were identified as homologues of Drosophila sine oculis which is essential for compound-eye formation. The Six proteins are characterized by the Six domain and the Six-type homeodomain, both of which are essential for specific DNA binding and for cooperative interactions with Eya proteins. Mammals possess six Six genes which can be subdivided into three subclasses, and mutations of Six genes have been identified in human genetic disorders. Characterization of Six genes from various animal phyla revealed the antiquity of this gene family and roles of its members in several different developmental contexts. Some members retain conserved roles as components of the Pax-Six-Eya-Dach regulatory network, which may have been established in the common ancestor of all bilaterians as a toolbox controlling cell proliferation and cell movement during embryogenesis. Gene duplications and cis-regulatory changes may have provided a basis for diverse functions of Six genes in different animal lineages.     

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.     

Concerted gene duplications in the two keratin gene families

Summary Evolutionary trees were derived from the keratin protein sequences using the Phylogeny Analysis Using Parsimony (PAUP) set of programs. Three major unexpected conclusions were derived from the analysis: The smallest keratin protein subunit, K#19 (Moll et al. 1982), is not the most primitive one, but has evolved to fulfill a highly specialized function, presumably to redress the unbalanced synthesis of keratin subunits. Second, the ancestors of keratins expressed in the early embryonic stages, K#8 and K#18, were the first to diverge from the ancestors of all the other keratins. The branches leading to these two keratins are relatively short, indicating a comparatively strong selection against changes in the sequences of these two proteins. Third, the two keratin families show extraodinary parallelism in their patterns of gene duplications. In both families the genes expressed in embryos diverged first, later bursts of gene duplications created the subfamilies expressed in various differentiated cells, and relatively recent gene duplications gave rise to the hair keratin genes and separated the basal cell-specific keratin from those expressed under hyperproliferative conditions. The parallelism of gene duplications in the two keratin gene families implies a mechanism in which duplications in one family influence duplication events in the other family.     

The monosaccharide transporter gene family in Arabidopsis and rice: a history of duplications, adaptive evolution, and functional divergence

Current hypotheses of gene duplicate divergence propose that surviving members of a gene duplicate pair may evolve, under conditions of purifying or nearly neutral selection, in one of two ways: with new function arising in one duplicate while the other retains original function (neofunctionalization [NF]) or partitioning of the original function between the 2 paralogs (subfunctionalization [SF]). More recent studies propose that SF followed by NF (subneofunctionalization [SNF]) explains the divergence of many duplicate genes. In this analysis, we evaluate these hypotheses in the context of the large monosaccharide transporter (MST) gene families in Arabidopsis and rice. MSTs have an ancient origin, predating plants, and have evolved in the seed plant lineage to comprise 7 subfamilies. In Arabidopsis, 53 putative MST genes have been identified, with one subfamily greatly expanded by tandem gene duplications. We searched the rice genome for members of the MST gene family and compared them with the MST gene family in Arabidopsis to determine subfamily expansion patterns and estimate gene duplicate divergence times. We tested hypotheses of gene duplicate divergence in 24 paralog pairs by comparing protein sequence divergence rates, estimating positive selection on codon sites, and analyzing tissue expression patterns. Results reveal the MST gene family to be significantly larger (65) in rice with 2 subfamilies greatly expanded by tandem duplications. Gene duplicate divergence time estimates indicate that early diversification of most subfamilies occurred in the Proterozoic (2500-540 Myr) and that expansion of large subfamilies continued through the Cenozoic (65-0 Myr). Two-thirds of paralog pairs show statistically symmetric rates of sequence evolution, most consistent with the SF model, with half of those showing evidence for positive selection in one or both genes. Among 8 paralog pairs showing asymmetric divergence rates, most consistent with the NF model, nearly half show evidence of positive selection. Positive selection does not appear in any duplicate pairs younger than approximately 34 Myr. Our data suggest that the NF, SF, and SNF models describe different outcomes along a continuum of divergence resulting from initial conditions of relaxed constraint after duplication.     

Voltage-Gated Potassium Channel Genes Are Clustered in Paralogous Regions of the Mouse Genome

Cloning of the Drosophila Shaker gene established that a neurological phenotype including locomotor dysfunction can be caused by a mutation in a voltage-gated potassium (K) channel gene. Shaker sequences have been used to isolate a large family of related K channel genes from both flies and mammals. Toward elucidating the evolutionary relationship between loci and the potential causal connection that K channels may have to mammalian genetic disorders, we report here the genetic mapping of 12-16 different murine, voltage-gated K channel genes. We find that multiple genes, in some cases from distantly related K channel subfamilies, occur in clusters in the mouse genome. These mapping results suggest that the K channel gene subfamilies arose through ancient localized gene duplication events, followed by chromosomal duplications and rearrangements as well as further gene duplication. We also note that several neurologic disorders of both mouse and human are associated with the chromosomal regions containing K channel genes.     

Extensive duplications of phototransduction genes in early vertebrate evolution correlate with block (chromosome) duplications

Many gene families in mammals have members that are expressed more or less uniquely in the retina or differentially in specific retinal cell types. We describe here analyses of nine such gene families with regard to phylogenetic relationships and chromosomal location. The families are opsins, G proteins (alpha, beta, and gamma subunits), phosphodiesterases type 6, cyclic nucleotide-gated channels, G-protein-coupled receptor kinases, arrestins, and recoverins. The results suggest that multiple new gene copies arose in all of these families very early in vertebrate evolution during a period with extensive gene duplications. Many of the new genes arose through duplications of large chromosome regions (blocks of genes) or even entire chromosomes, as shown by linkage with other gene families. Some of the phototransduction families belong to the same duplicated regions and were thus duplicated simultaneously. We conclude that gene duplications in early vertebrate evolution probably helped facilitate the specialization of the retina and the subspecialization of different retinal cell types.     

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.     

Evolution of the CYP2ABFGST gene cluster in rat, and a fine-scale comparison among rodent and primate species

The evolution of gene families can be best understood by studying the modern organization and functions of family members, and by comparing parallel families in different species. In this study, the CYP2ABFGST gene cluster has been characterized in rat and compared to the syntenic clusters in mouse and human, providing an interesting example of gene family evolution. In the rat, 18 loci from six subfamilies have been identified by specifically amplifying and sequencing gene fragments from cloned DNA, and have been exactly placed on chromosome 1. The overall organization of the gene cluster in rat is relatively simple, with genes from each subfamily in tandem, and is more similar to the mouse than to the human cluster. We have reconstructed the probable structure of the CYP2ABFGST cluster in the common ancestor of primates and rodents, and inferred a model of the evolution of this gene cluster in the three species. Numerous nontandem and block duplications, inversions, and translocations have occurred entirely inside the cluster, indicating that pairing between duplicate genes is keeping the rearrangements within the cluster region. The initial tandem duplication of a CYP2 gene in an early mammalian ancestor has made this region particularly subject to such localized rearrangements. Even if duplicated genes do not have a large-scale effect on chromosomal rearrangements, on a local level clustered gene families may have contributed significantly to the genomic complexity of modern mammals.     

A comprehensive expression analysis of all members of a gene family encoding cell-wall enzymes allowed us to predict cis-regulatory regions involved in cell-wall construction in specific organs of Arabidopsis.

The Arabidopsis thaliana genome sequencing project has revealed that multigene families, such as those generated by genome duplications, are more abundant among plant genomes than among animal genomes. To gain insight into the evolutionary implications of the multigene families in higher plants, we examined the XTH gene family, a group of genes encoding xyloglucan endotransglucosylase/hydrolase, which are responsible for cell-wall construction in plants. Expression analysis of all members (33 genes) of this family, using quantitative real-time RT-PCR, revealed that most members exhibit distinct expression profiles in terms of tissue specificity and responses to hormonal signals, with some members exhibiting similar expression patterns. By comparing the flanking sequences of individual genes, we identified four sets of large-segment duplications and two sets of solitary gene duplications. In each set of gene duplicates, long nucleotide sequences, ranging from one to two hundred base pairs, are conserved. Furthermore, gene duplicates exhibit similar organ-specific expression profiles. These facts allowed us to predict putative cis-regulatory regions, particularly those responsible for cell-wall construction, and hence for morphogenesis, that are specific for certain organs or tissues in plants.     

Pax6-dependence of Six3, Eya1 and Dach1 expression during lens and nasal placode induction

The Drosophila eyeless gene plays a central role in fly eye development and controls a subordinate regulatory network consisting of the so, eya and dac genes. All three genes have highly conserved mammalian homologs, suggesting possible conservation of this eye forming regulatory network. sine oculis (so) belongs to the so/Six gene family, and Six3 is prominently expressed in the developing mammalian eye. Eya1 and Dach1 are mammalian homologs of eya and dac, respectively, and although neither Eya1 nor Dach1 knockout mice express prenatal eye defects, possibilities exist for postnatal ocular phenotypes or for functional redundancy between related family members. To examine whether expression relationships analogous to those between ey, so, eya and dac exist in early mammalian oculogenesis, we investigated Pax6, Six3, Eya1 and Dach1 protein expression in murine lens and nasal placode development. Six3 expression in the pre-placode lens ectoderm is initially Pax6-independent, but subsequently both its expression and nuclear localization become Pax6-dependent. Six3, Dach1 and Eya1 nasal expression in pre-placode ectoderm are also initially Pax6-independent, but thereafter become Pax6-dependent. Pax6, Six3, Dach1 and Eya1 are all co-expressed in the developing ciliary marginal zone, a source of retinal stem cells in some vertebrates. An in vitro protein-protein interaction is detected between Six3 and Eya1. Collectively, these findings suggest that the Pax-Eya-Six-Dach network is at best only partly conserved during lens and nasal placode development. However, the findings do not rule out the possibility that such a regulatory network acts at later stages of oculogenesis.     

Expansion of the receptor-like kinase/Pelle gene family and receptor-like proteins in Arabidopsis

Receptor-like kinases (RLKs) are a family of transmembrane proteins with versatile N-terminal extracellular domains and C-terminal intracellular kinases. They control a wide range of physiological responses in plants and belong to one of the largest gene families in the Arabidopsis genome with more than 600 members. Interestingly, this gene family constitutes 60% of all kinases in Arabidopsis and accounts for nearly all transmembrane kinases in Arabidopsis. Analysis of four fungal, six metazoan, and two Plasmodium sp. genomes indicates that the family was represented in all but fungal genomes, indicating an ancient origin for the family with a more recent expansion only in the plant lineages. The RLK/Pelle family can be divided into several subfamilies based on three independent criteria: the phylogeny based on kinase domain sequences, the extracellular domain identities, and intron locations and phases. A large number of receptor-like proteins (RLPs) resembling the extracellular domains of RLKs are also found in the Arabidopsis genome. However, not all RLK subfamilies have corresponding RLPs. Several RLK/Pelle subfamilies have undergone differential expansions. More than 33% of the RLK/Pelle members are found in tandem clusters, substantially higher than the genome average. In addition, 470 of the RLK/Pelle family members are located within the segmentally duplicated regions in the Arabidopsis genome and 268 of them have a close relative in the corresponding regions. Therefore, tandem duplications and segmental/whole-genome duplications represent two of the major mechanisms for the expansion of the RLK/Pelle family in Arabidopsis.     

The Sine oculis/Six class family of homeobox genes in jellyfish with and without eyes: development and eye regeneration

The development of visual organs is regulated in Bilateria by a network of genes where members of the Six and Pax gene families play a central role. To investigate the molecular aspects of eye evolution, we analyzed the structure and expression patterns of cognate members of the Six family genes in jellyfish (Cnidaria, Hydrozoa), representatives of a basal, non-bilaterian phylum where complex lens eyes with spherical lens, an epidermal cornea, and a retina appear for the first time in evolution. In the jellyfish Cladonema radiatum, a species with well-developed lens eyes in the tentacle bulbs, Six1/2-Cr and Six3/6-Cr, are expressed in the eye cup. Six4/5-Cr is mainly expressed in the manubrium, the feeding, and sex organ. All three Six genes are expressed in different subsets of epidermal nerve cells, possibly of the RFamide type which are part of a net connecting the different eyes with each other and the effector organs. Furthermore, expression is found in other tissues, notably in the striated muscle. During eye regeneration, expression of Six1/2-Cr and Six3/6-Cr is upregulated, but not of Six4/5-Cr. In Podocoryne carnea, a jellyfish without eyes, Six1/2-Pc and Six3/6-Pc are also expressed in the tentacle bulbs, Six1/2-Pc additionally in the manubrium and striated muscle, and Six3/6-Pc in the mechanosensory nematocytes of the tentacle. The conserved gene structure and expression patterns of all Cladonema Six genes suggest broad conservation of upstream regulatory mechanisms in eye development.