Purifying selection and birth-and-death evolution in the ubiquitin gene family

Ubiquitin is a highly conserved protein that is encoded by a multigene family. It is generally believed that this gene family is subject to concerted evolution, which homogenizes the member genes of the family. However, protein homogeneity can be attained also by strong purifying selection. We therefore studied the proportion (pS) of synonymous nucleotide differences between members of the ubiquitin gene family from 28 species of fungi, plants, and animals. The results have shown that pS is generally very high and is often close to the saturation level, although the protein sequence is virtually identical for all ubiquitins from fungi, plants, and animals. A small proportion of species showed a low level of pS values, but these values appeared to be caused by recent gene duplication. It was also found that the number of repeat copies of the gene family varies considerably with species, and some species harbor pseudogenes. These observations suggest that the members of this gene family evolve almost independently by silent nucleotide substitution and are subjected to birth-and-death evolution at the DNA level.

The evolution of the vertebrate Dlx gene family.

The vertebrate Dlx gene family consists of homeobox-containing transcription factors distributed in pairs on the same chromosomes as the Hox genes. To investigate the evolutionary history of Dlx genes, we have cloned five new zebrafish family members and have provided additional sequence information for two mouse genes. Phylogenetic analyses of Dlx gene sequences considered in the context of their chromosomal arrangements suggest that an initial tandem duplication produced a linked pair of Dlx genes after the divergence of chordates and arthropods but prior to the divergence of tunicates and vertebrates. This pair of Dlx genes was then duplicated in the chromosomal events that led to the four clusters of Hox genes characteristic of bony fish and tetrapods. It is possible that a pair of Dlx genes linked to the Hoxc cluster has been lost from mammals. We were unable to distinguish between independent duplication and retention of the ancestral state of bony vertebrates to explain the presence of a greater number of Dlx genes in zebrafish than mammals. Determination of the linkage relationship of these additional zebrafish Dlx genes to Hox clusters should help resolve this issue.

Duplication and functional divergence in the chalcone synthase gene family of Asteraceae: evolution with substrate change and catalytic simplification.

Plant-specific polyketide synthase genes constitute a gene superfamily, including universal chalcone synthase [CHS; malonyl-CoA:4-coumaroyl-CoA malonyltransferase (cyclizing) (EC 2.3.1.74)] genes, sporadically distributed stilbene synthase (SS) genes, and atypical, as-yet-uncharacterized CHS-like genes. We have recently isolated from Gerbera hybrida (Asteraceae) an unusual CHS-like gene, GCHS2, which codes for an enzyme with structural and enzymatic properties as well as ontogenetic distribution distinct from both CHS and SS. Here, we show that the GCHS2-like function is encoded in the Gerbera genome by a family of at least three transcriptionally active genes. Conservation within the GCHS2 family was exploited with selective PCR to study the occurrence of GCHS2-like genes in other Asteraceae. Parsimony analysis of the amplified sequences together with CHS-like genes isolated from other taxa of angiosperm subclass Asteridae suggests that GCHS2 has evolved from CHS via a gene duplication event that occurred before the diversification of the Asteraceae. Enzyme activity analysis of proteins produced in vitro indicates that the GCHS2 reaction is a non-SS variant of the CHS reaction, with both different substrate specificity (to benzoyl-CoA) and a truncated catalytic profile. Together with the recent results of Durbin et al. [Durbin, M. L., Learn, G. H., Jr., Huttley, G. A. & Clegg, M. T. (1995) Proc. Natl. Acad. Sci. USA 92, 3338-3342], our study confirms a gene duplication-based model that explains how various related functions have arisen from CHS during plant evolution.

Evolution of the chalcone synthase gene family in the genus Ipomoea.

The evolution of the chalcone synthase [CHS; malonyl-CoA:4-coumaroyl-CoA malonyltransferase (cyclizing), EC 2.3.1.74] multigene family in the genus Ipomoea is explored. Thirteen CHS genes from seven Ipomoea species (family Convolvulaceae) were sequenced--three from genomic clones and the remainder from PCR amplification with primers designed from the 5' flanking region and the end of the 3' coding region of Ipomoea purpurea Roth. Analysis of the data indicates a duplication of CHS that predates the divergence of the Ipomoea species in this study. The Ipomoea CHS genes are among the most rapidly evolving of the CHS genes sequenced to date. The CHS genes in this study are most closely related to the Petunia CHS-B gene, which is also rapidly evolving and highly divergent from the rest of the Petunia CHS sequences.

Emergence of the ZNF91 Krüppel-associated box-containing zinc finger gene family in the last common ancestor of anthropoidea.

The ZNF91 gene family, a subset of the Krüppel-associated box (KRAB)-containing group of zinc finger genes, comprises more than 40 loci; most reside on human chromosome 19p12-p13.1. We have examined the emergence and evolutionary conservation of the ZNF91 family. ZNF91 family members were detected in all species of great apes, gibbons, Old World monkeys, and New World monkeys examined but were not found in prosimians or rodents. In each species containing the ZNF91 family, the genes were clustered at one major site, on the chromosome(s) syntenic to human chromosome 19. To identify a putative "founder" gene, > 20 murine KRAB-containing zinc finger protein (ZFP) cDNAs were randomly cloned, but none showed sequence similarity to the ZNF91 genes. These observations suggest that the ZNF91 gene cluster is a derived character specific to Anthropoidea, resulting from a duplication and amplification event some 55 million years ago in the common ancestor of simians. Although the ZNF91 gene cluster is present in all simian species, the sequences of the human ZNF91 gene that confer DNA-binding specificity were conserved only in great apes, suggesting that there is not a high selective pressure to maintain the DNA targets of these proteins during evolution.

A third member of the synapsin gene family

Synapsins are a family of neuron-specific synaptic vesicle-associated phosphoproteins that have been implicated in synaptogenesis and in the modulation of neurotransmitter release. In mammals, distinct genes for synapsins I and II have been identified, each of which gives rise to two alternatively spliced isoforms. We have now cloned and characterized a third member of the synapsin gene family, synapsin III, from human DNA. Synapsin III gives rise to at least one protein isoform, designated synapsin IIIa, in several mammalian species. Synapsin IIIa is associated with synaptic vesicles, and its expression appears to be neuron-specific. The primary structure of synapsin IIIa conforms to the domain model previously described for the synapsin family, with domains A, C, and E exhibiting the highest degree of conservation. Synapsin IIIa contains a novel domain, termed domain J, located between domains C and E. The similarities among synapsins I, II, and III in domain organization, neuron-specific expression, and subcellular localization suggest a possible role for synapsin III in the regulation of neurotransmitter release and synaptogenesis. The human synapsin III gene is located on chromosome 22q12–13, which has been identified as a possible schizophrenia susceptibility locus. On the basis of this localization and the well established neurobiological roles of the synapsins, synapsin III represents a candidate gene for schizophrenia.

Amyrel, a paralogous gene of the amylase gene family in Drosophila melanogaster and the Sophophora subgenus

We describe a gene from Drosophila melanogaster related to the alpha-amylase gene Amy. This gene, which exists as a single copy, was named Amyrel. It is strikingly divergent from Amy because the amino acid divergence is 40%. The coding sequence is interrupted by a short intron at position 655, which is unusual in amylase genes. Amyrel has also been cloned in Drosophila ananassae, Drosophila pseudoobscura, and Drosophila subobscura and is likely to be present throughout the Sophophora subgenus, but, to our knowledge, it has not been detected outside. Unexpectedly, there is a strong conservation of 5′ and 3′ flanking regions between Amyrel genes from different species, which is not the case for Amy and which suggests that selection acts on these regions. In contrast to the Amy genes, Amyrel is transcribed in larvae of D. melanogaster but not in adults. However, the protein has not been detected yet. Amyrel evolves about twice as fast as Amy in the several species studied. We suggest that this gene could result from a duplication of Amy followed by accelerated and selected divergence toward a new adaptation.

Three members of a novel small gene-family from Arabidopsis thaliana able to complement functionally an Escherichia coli mutant defective in PAPS reductase activity encode proteins with a thioredoxin-like domain and “APS reductase” activity

Three different cDNAs, Prh-19, Prh-26, and Prh-43 [3′-phosphoadenosine-5′-phosphosulfate (PAPS) reductase homolog], have been isolated by complementation of an Escherichia coli cysH mutant, defective in PAPS reductase activity, to prototrophy with an Arabidopsis thaliana cDNA library in the expression vector λYES. Sequence analysis of the cDNAs revealed continuous open reading frames encoding polypeptides of 465, 458, and 453 amino acids, with calculated molecular masses of 51.3, 50.5, and 50.4 kDa, respectively, that have strong homology with fungal, yeast, and bacterial PAPS reductases. However, unlike microbial PAPS reductases, each PRH protein has an N-terminal extension, characteristic of a plastid transit peptide, and a C-terminal extension that has amino acid and deduced three-dimensional homology to thioredoxin proteins. Adenosine 5′-phosphosulfate (APS) was shown to be a much more efficient substrate than PAPS when the activity of the PRH proteins was tested by their ability to convert 35S-labeled substrate to acid-volatile 35S-sulfite. We speculate that the thioredoxin-like domain is involved in catalytic function, and that the PRH proteins may function as novel “APS reductase” enzymes. Southern hybridization analysis showed the presence of a small multigene family in the Arabidopsis genome. RNA blot hybridization with gene-specific probes revealed for each gene the presence of a transcript of ≈1.85 kb in leaves, stems, and roots that increased on sulfate starvation. To our knowledge, this is the first report of the cloning and characterization of plant genes that encode proteins with APS reductase activity and supports the suggestion that APS can be utilized directly, without activation to PAPS, as an intermediary substrate in reductive sulfate assimilation.

Mendel-GFDb and Mendel-ESTS: databases of plant gene families and ESTs annotated with gene family numbers and gene family names

There is no control over the information provided with sequences when they are deposited in the sequence databases. Consequently mistakes can seed the incorrect annotation of other sequences. Grouping genes into families and applying controlled annotation overcomes the problems of incorrect annotation associated with individual sequences. Two databases (http://www.mendel.ac.uk) were created to apply controlled annotation to plant genes and plant ESTs: Mendel-GFDb is a database of plant protein (gene) families based on gapped-BLAST analysis of all sequences in the SWISS-PROT family of databases. Sequences are aligned (ClustalW) and identical and similar residues shaded. The families are visually curated to ensure that one or more criteria, for example overall relatedness and/or domain similarity relate all sequences within a family. Sequence families are assigned a ‘Gene Family Number’ and a unified description is developed which best describes the family and its members. If authority exists the gene family is assigned a ‘Gene Family Name’. This information is placed in Mendel-GFDb. Mendel-ESTS is primarily a database of plant ESTs, which have been compared to Mendel-GFDb, completely sequenced genomes and domain databases. This approach associated ESTs with individual sequences and the controlled annotation of gene families and protein domains; the information being placed in Mendel-ESTS. The controlled annotation applied to genes and ESTs provides a basis from which a plant transcription database can be developed.

The Arabidopsis CBF Gene Family Is Composed of Three Genes Encoding AP2 Domain-Containing Proteins Whose Expression Is Regulated by Low Temperature but Not by Abscisic Acid or Dehydration1

We have identified two genes from Arabidopsis that show high similarity with CBF1, a gene encoding an AP2 domain-containing transcriptional activator that binds to the low-temperature-responsive element CCGAC and induces the expression of some cold-regulated genes, increasing plant freezing tolerance. These two genes, which we have named CBF2 and CBF3, also encode proteins containing AP2 DNA-binding motifs. Furthermore, like CBF1, CBF2 and CBF3 proteins also include putative nuclear-localization signals and potential acidic activation domains. The CBF2 and CBF3 genes are linked to CBF1, constituting a cluster on the bottom arm of chromosome IV. The high level of similarity among the three CBF genes, their tandem organization, and the fact that they have the same transcriptional orientation all suggest a common origin. CBF1, CBF2, and CBF3 show identical expression patterns, being induced very rapidly by low-temperature treatment. However, in contrast to most of the cold-induced plant genes characterized, they are not responsive to abscisic acid or dehydration. Taken together, all of these data suggest that CBF2 and CBF3 may function as transcriptional activators, controlling the level of low-temperature gene expression and promoting freezing tolerance through an abscisic acid-independent pathway.

Differential Expression of Three Members of the 1-Aminocyclopropane-1-Carboxylate Synthase Gene Family in Carnation1

We investigated the expression patterns of three 1-aminocyclopropane-1-carboxylate (ACC) synthase genes in carnation (Dianthus caryophyllus cv White Sim) under conditions previously shown to induce ethylene biosynthesis. These included treatment of flowers with 2,4-dichlorophenoxyacetic acid, ethylene, LiCl, cycloheximide, and natural and pollination-induced flower senescence. Accumulation of ACC synthase transcripts in leaves following mechanical wounding and treatment with 2,4-dichlorophenoxyacetic acid or LiCl was also determined by RNA gel-blot analysis. As in other species, the carnation ACC synthase genes were found to be differentially regulated in a tissue-specific manner. DCACS2 and DCACS3 were preferentially expressed in styles, whereas DCACS1 mRNA was most abundant in petals. Cycloheximide did not induce increased accumulation of ACC synthase transcripts in carnation flowers, whereas the expression of ACC synthase was up-regulated by auxin, ethylene, LiCl, pollination, and senescence in a floral-organ-specific manner. Expression of the three ACC synthases identified in carnation did not correspond to elevated ethylene biosynthesis from wounded or auxin-treated leaves, and there are likely additional members of the carnation ACC synthase gene family responsible for ACC synthase expression in vegetative tissues.

Human tumor suppressor EXT gene family members EXTL1 and EXTL3 encode α1,4- N-acetylglucosaminyltransferases that likely are involved in heparan sulfate/ heparin biosynthesis

The tumor suppressors EXT1 and EXT2 are associated with hereditary multiple exostoses and encode bifunctional glycosyltransferases essential for chain polymerization of heparan sulfate (HS) and its analog, heparin (Hep). Three highly homologous EXT-like genes, EXTL1–EXTL3, have been cloned, and EXTL2 is an α1,4-GlcNAc transferase I, the key enzyme that initiates the HS/Hep synthesis. In the present study, truncated forms of EXTL1 and EXTL3, lacking the putative NH2-terminal transmembrane and cytoplasmic domains, were transiently expressed in COS-1 cells and found to harbor α-GlcNAc transferase activity. EXTL3 used not only N-acetylheparosan oligosaccharides that represent growing HS chains but also GlcAβ1–3Galβ1-O-C2H4NH-benzyloxycarbonyl (Cbz), a synthetic substrate for α-GlcNAc transferase I that determines and initiates HS/Hep synthesis. In contrast, EXTL1 used only the former acceptor. Neither EXTL1 nor EXTL3 showed any glucuronyltransferase activity as examined with N-acetylheparosan oligosaccharides. Heparitinase I digestion of each transferase-reaction product showed that GlcNAc had been transferred exclusively through an α1,4-configuration. Hence, EXTL3 most likely is involved in both chain initiation and elongation, whereas EXTL1 possibly is involved only in the chain elongation of HS and, maybe, Hep as well. Thus, their acceptor specificities of the five family members are overlapping but distinct from each other, except for EXT1 and EXT2 with the same specificity. It now has been clarified that all of the five cloned human EXT gene family proteins harbor glycosyltransferase activities, which probably contribute to the synthesis of HS and Hep.

A gene family required for human germ cell development evolved from an ancient meiotic gene conserved in metazoans

The Deleted in AZoospermia (DAZ) genes encode potential RNA-binding proteins that are expressed exclusively in prenatal and postnatal germ cells and are strong candidates for human fertility factors. Here we report the identification of an additional member of the DAZ gene family, which we have called BOULE. With the identification of this gene, it is clear that the human DAZ gene family contains at least three members: DAZ, a Y-chromosome gene cluster that arose 30–40 million years ago and whose deletion is linked to infertility in men; DAZL, the “father” of DAZ, a gene that maps to human chromosome 3 and has homologs required for both female and male germ cell development in other organisms; and BOULE, a gene that we propose is the “grandfather” of DAZ and maps to human chromosome 2. Human and mouse BOULE resemble the invertebrate meiotic regulator Boule, the proposed ortholog of DAZ, in sequence and expression pattern and hence likely perform a similar meiotic function. In contrast, the previously identified human DAZ and DAZL are expressed much earlier than BOULE in prenatal germ stem cells and spermatogonia; DAZL also is expressed in female germ cells. These data suggest that homologs of the DAZ gene family can be grouped into two subfamilies (BOULE and DAZL) and that members of the DAZ family evolved from an ancestral meiotic regulator, Boule, to assume distinct, yet overlapping, functions in germ cell development.

Three Drought-Responsive Members of the Nonspecific Lipid-Transfer Protein Gene Family in Lycopersicon pennellii Show Different Developmental Patterns of Expression1

Genomic clones of two nonspecific lipid-transfer protein genes from a drought-tolerant wild species of tomato (Lycopersicon pennellii Corr.) were isolated using as a probe a drought- and abscisic acid (ABA)-induced cDNA clone (pLE16) from cultivated tomato (Lycopersicon esculentum Mill.). Both genes (LpLtp1 and LpLtp2) were sequenced and their corresponding mRNAs were characterized; they are both interrupted by a single intron at identical positions and predict basic proteins of 114 amino acid residues. Genomic Southern data indicated that these genes are members of a small gene family in Lycopersicon spp. The 3′-untranslated regions from LpLtp1 and LpLtp2, as well as a polymerase chain reaction-amplified 3′-untranslated region from pLE16 (cross-hybridizing to a third gene in L. pennellii, namely LpLtp3), were used as gene-specific probes to describe expression in L. pennellii through northern-blot analyses. All LpLtp genes were exclusively expressed in the aerial tissues of the plant and all were drought and ABA inducible. Each gene had a different pattern of expression in fruit, and LpLtp1 and LpLtp2, unlike LpLtp3, were both primarily developmentally regulated in leaf tissue. Putative ABA-responsive elements were found in the proximal promoter regions of LpLtp1 and LpLtp2.

Conserved promoter elements in the CYP6B gene family suggest common ancestry for cytochrome P450 monooxygenases mediating furanocoumarin detoxification.

Despite the fact that Papilio glaucus and Papilio polyxenes share no single hostplant species, both species feed to varying extents on hostplants that contain furanocoumarins. P. glaucus contains two nearly identical genes, CYP6B4v2 and CYP6B5v1, and P. polyxenes contains two related genes, CYP6B1v3 and CYP6B3v2. Except for CYP6B3v2, the substrate specificity of which has not yet been defined, each of the encoded cytochrome P450 monooxygenases (P450s) metabolizes an array of linear furanocoumarins. All four genes are transcriptionally induced in larvae by exposure to the furanocoumarin xanthotoxin; several are also induced by other furanocoumarins. Comparisons of the organizational structures of these genes indicate that all have the same intron/exon arrangement. Sequences in the promoter regions of the P. glaucus CYP6B4v2/CYP6B5v1 genes and the P. polyxenes CYP6B3v2 gene are similar but not identical to the -146 to -97 region of CYP6B1v3 gene, which contains a xanthotoxin-responsive element (XRE-xan) important for basal and xanthotoxin-inducible transcription of CYP6B1v3. Complements of the xenobiotic-responsive element (XRE-AhR) in the dioxin-inducible human and rat CYP1A1 genes also exist in all four promoters, suggesting that these genes may be regulated by dioxin. Antioxidant-responsive elements (AREs) in mouse and rat glutathione S-transferase genes and the Barbie box element (Bar) in the bacterial CYP102 gene exist in the CYP6B1v3, CYP6B4v2, and CYP6B5v1 promoters. Similarities in the protein sequences, intron positions, and xanthotoxin- and xenobiotic-responsive promoter elements indicate that these insect CYP6B genes are derived from a common ancestral gene. Evolutionary comparisons between these P450 genes are the first available for a group of insect genes transcriptionally regulated by hostplant allelochemicals and provide insights into the process by which insects evolve specialized feeding habits.