Phylogenetic analysis and identification of pseudogenes reveal a progressive loss of zona pellucida genes during evolution of vertebrates

Vertebrate eggs are surrounded by an extracellular matrix with similar functions and conserved individual components: the zona pellucida (ZP) glycoproteins. In mammals, chickens, frogs, and some fish species, we established an updated list of the ZP genes, studied the relationships within the ZP gene family using phylogenetic analysis, and identified ZP pseudogenes. Our study confirmed the classification of ZP genes in six subfamilies: ZPA/ZP2, ZPB/ZP4, ZPC/ZP3, ZP1, ZPAX, and ZPD. The identification of a Zpb pseudogene in the mouse genome, Zp1 pseudogenes in the dog and bovine genomes, and Zpax pseudogenes in the human, chimpanzee, macaque, and bovine genomes showed that the evolution of ZP genes mainly occurs by death of genes. Our study revealed that the extracellular matrix surrounding vertebrate eggs contains three to at least six ZP glycoproteins. Mammals can be classified in three categories. In the mouse, the ZP is composed of three ZP proteins (ZPA/ZP2, ZPC/ZP3, and ZP1). In dog, cattle and, putatively, pig, cat, and rabbit, the zona is composed of three ZP proteins (ZPA/ZP2, ZPB/ZP4, and ZPC/ZP3). In human, chimpanzee, macaque, and rat, the ZP is composed of four ZP proteins (ZPA/ZP2, ZPB/ZP4, ZPC/ZP3, and ZP1). Our review provides new directions to investigate the molecular basis of sperm-egg recognition, a mechanism which is not yet elucidated.     

ZP genes in avian species illustrate the dynamic evolution of the vertebrate egg envelope

The vertebrate egg envelope is composed of a set of related proteins, encoded by the ZP genes. The apparent simplicity of the egg envelope is in contrast to the number of ZP genes identified by conventional cloning and data mining of genome sequences from a number of vertebrates. The vertebrate ZP genes fall into five classes, ZP1, ZP2, ZP3, ZPD and ZPAX. Analysis of chicken genome and EST sequence data has revealed the presence of seven distinct ZP genes, falling into these classes that are expressed in the female reproductive system. Comparison with the repertoire of ZP genes in other vertebrates suggests a major source of diversity in the composition of the egg envelope is a continual process of amplification, diversification and attrition of ZP gene sequences.     

A hatching enzyme substrate in the Xenopus laevis egg envelope is a high molecular weight ZPA homolog

The Xenopus laevis egg envelope is composed of six or more glycoproteins, three of which have been cloned and identified as the mammalian homologs ZPA (ZP2), ZPB (ZP1) and ZPC (ZP3). The remaining glycoproteins are a triplet of high molecular weight components that are selectively hydrolyzed by the hatching enzyme. We have isolated one of these proteins and cloned its cDNA. The mRNA for the protein was found to be expressed only in early stage oocytes, as are other envelope components. From the deduced amino acid sequence, it was indicated to be a secreted glycoprotein with a characteristic ZP domain in the C-terminal half of the molecule. The N-terminal half was unrelated to any known glycoprotein. Comparative sequence analysis of the ZP domain indicated that it was derived from an ancestor of ZPA and ZPB, with the greatest identity to ZPA. This envelope component has been designated ZPAX.     

Identification of the true human orthologue of the mouse Zp1 gene: evidence for greater complexity in the mammalian zona pellucida?

The mammalian zona pellucida is a mixture of glycoproteins, believed to be encoded by three distinct genes, ZP1/ZPB, ZP2/ZPA, and ZP3/ZPC. We have now determined that the true human orthologue of the mouse Zp1 gene is not ZPB, but that there is a distinct human ZP1 gene. Comparison of the human ZP1 and murine Zp1 genes indicates significant conservation of nucleotide and amino acid sequences, of intron-exon size and organisation, and of regulatory sequences. In addition, the mouse and human ZP1 genes are in a region of conserved synteny between human chromosome 11 and mouse chromosome 19.     

Evolution and nomenclature of the zona pellucida gene family

Three subfamilies of genes are acknowledged within the zona pellucida (ZP) gene family. At present, these subfamilies each have two names that are used interchangeably: ZPA or ZP2, ZPB or ZP1, and ZPC or ZP3. The ZPA genes encode the longest protein sequences and the ZPC genes the shortest. Recently, several sequences, which have no clear relationship to the three subfamilies, have been identified. These sequences include two paralogous ZP genes from Xenopus laevis and a single gene from the fish Oryzias latipes. We have conducted extensive phylogenetic analyses of the known ZP genes. As well as establishing the evolutionary relationships among these genes, the analyses make it clear that the dual nomenclature system is no longer feasible, because major paralogous groups are present in the ZPB (ZP1) family of genes of amniotes. We propose a unified system of nomenclature for the ZP gene family that removes the existing ambiguities.     

Cluster of genes encoding the major egg envelope protein of zebrafish

The proteins of fish egg envelopes are encoded by genes that are closely related to the genes for human zona pellucida proteins. A cluster of three genes coding for an egg envelope protein was isolated from the zebrafish, Danio rerio. The three genes, zp2a, zp2b, and zp2c, are located within an 11 kb region and are each comprised of eight exons spanning 1.85 kb. The exon-intron structures of the genes are nearly identical; however, their deduced amino acid sequences diverge at exon 7 (zp2b and zp2c from zp2a) and exon 8 (zp2c from zp2b). Exons 2-7 have a structural organization similar to exons in the carboxy-terminal half of the human zona pellucida ZP1, ZP2, and ZPB genes, suggesting they arose from a common ancestral gene. Sequence comparisons indicate that the deduced zebrafish proteins are most closely related to human ZPB. Zebrafish mRNAs coding for each of the three ZP2 variants have been found either as full-length cDNAs or expressed sequence tags. Distinct from the wf(female) gene of winter flounder which we first reported (Lyons et al., 1993: J Biol Chem 268:21351-21358), expression of the zebrafish zp2 genes was found to be ovary-specific, instead of liver-specific, and the promoter regions of zp2a and zp2b, while different, both contained E-box sequences (CANNTG) that have been demonstrated to be essential for coordination of zona pellucida gene expression in mammalian oocytes. Mixed peptide sequence analysis was used to identify the major polypeptide component of isolated zebrafish egg envelopes as the zp2 gene product.     

Reassessing the molecular biology of sperm-egg recognition with mouse genetics

The zona pellucida is an extracellular coat that surrounds mammalian eggs and early embryos. This insoluble matrix separates germ from somatic cells during folliculogenesis and plays critical roles during fertilization and early development. The mouse and human zona pellucida contain three glycoproteins (ZP1 or ZPB, ZP2, ZP3), the primary structures of which have been deduced by molecular cloning. Targeted mutagenesis of endogenous mouse genes and transgenesis with human homologues provide models to investigate the roles of individual zona components. Collectively, the genetic data indicate that no single mouse zona pellucida protein is obligatory for taxon-specific sperm binding and that two human proteins are not sufficient to support human sperm binding. An observed post-fertilization persistence of mouse sperm binding to "humanized" zona pellucida correlates with uncleaved ZP2. These observations are consistent with a model for sperm binding in which the supramolecular structure of the zona pellucida necessary for sperm binding is modulated by the cleavage status of ZP2.     

The compositional transition between the genomes of cold- and warm-blooded vertebrates: codon frequencies in orthologous genes

The genomes of the ancestors of mammals and birds underwent a compositional change in which the gene-richest regions increased their GC levels. Here we investigated this compositional transition by analyzing the levels of G and C in third codon positions, as well as the codon frequencies of orthologous genes from human, chicken and Xenopus. The results may be summed up as follows: (i) GC-poor genes, that did not undergo the compositional transition, showed only minor differences in orthologous sets from Xenopus, human and chicken; this is remarkable in view of the very many nucleotide substitutions that occurred over the long evolutionary times separating these species; (ii) GC-rich genes, that underwent the compositional transition, showed large differences between Xenopus and warm-blooded vertebrates, but not between chicken and human. In other words, the independent changes that occurred in avian and mammalian genes, on the average, were the same.     

Cloning and sequence comparison of the mouse, human, and chicken engrailed genes reveal potential functional domains and regulatory regions.

We have isolated and characterized genomic DNA clones for the human and chicken homologues of the mouse En-1 and En-2 genes and determined the genomic structure and predicted protein sequences of both En genes in all three species. Comparison of these vertebrate En sequences with the Xenopus En-2 [Hemmati-Brivanlou et al., 1991) and invertebrate engrailed-like genes showed that the two previously identified highly conserved regions within the En protein ]reviewed in Joyner and Hanks, 1991] can be divided into five distinct subregions, designated EH1 to EH5. Sequences 5' and 3' to the predicted coding regions of the vertebrate En genes were also analyzed in an attempt to identify cis-acting DNA sequences important for the regulation of En gene expression. Considerable sequence similarity was found between the mouse and human homologues both within the putative 5' and 3' untranslated as well as 5' flanking regions. Between the mouse and Xenopus En-2 genes, shorter stretches of sequence similarity were found within the 3' untranslated region. The 5' untranslated regions of the mouse, chicken and Xenopus En-2 genes, however, showed no similarly conserved stretches. In a preliminary analysis of the expression pattern of the human En genes, En-2 protein and RNA were detected in the embryonic and adult cerebellum respectively and not in other tissues tested. These patterns are analogous to those seen in other vertebrates. Taken together these results further strengthen the suggestion that En gene function and regulation has been conserved throughout vertebrate evolution and, along with the five highly conserved regions within the En protein, raise an interesting question about the presence of conserved genetic pathways.     

Regulation of vertebrate eye development by Rx genes

The paired-like homeobox-containing gene Rx has a critical role in the eye development of several vertebrate species including Xenopus, mouse, chicken, medaka, zebrafish and human. Rx is initially expressed in the anterior neural region of developing embryos, and later in the retina and ventral hypothalamus. Abnormal regulation or function of Rx results in severe abnormalities of eye formation. Overexpression of Rx in Xenopus and zebrafish embryos leads to overproliferation of retinal cells. A targeted elimination of Rx in mice results in a lack of eye formation. Mutations in Rx genes are the cause of the mouse mutation eyeless (ey1), the medaka temperature sensitive mutation eyeless (el) and the zebrafish mutation chokh. In humans, mutations in Rx lead to anophthalmia. All of these studies indicate that Rx genes are key factors in vertebrate eye formation. Because these results cannot be easily reconciled with the most popular dogmas of the field, we offer our interpretation of eye development and evolution.     

Evolutionary conservation of expression profiles between human and mouse orthologous genes

Mouse models are often used to study human genes because it is believed that the expression and function are similar for the majority of orthologous genes between the two species. However, recent comparisons of microarray data from thousands of orthologous human and mouse genes suggested rapid evolution of gene expression profiles under minimal or no selective constraint. These findings appear to contradict non-array-based observations from many individual genes and imply the uselessness of mouse models for studying human genes. Because absolute levels of gene expression are not comparable between species when the data are generated by species-specific microarrays, use of relative mRNA abundance among tissues (RA) is preferred to that of absolute expression signals. We thus reanalyze human and mouse genome-wide gene expression data generated by oligonucleotide microarrays. We show that the mean correlation coefficient among expression profiles detected by different probe sets of the same gene is only 0.38 for humans and 0.28 for mice, indicating that current measures of expression divergence are flawed because the large estimation error (discrepancy in expression signal detected by different probe sets of the same gene) is mistakenly included in the between-species divergence. When this error is subtracted, 84% of human-mouse orthologous gene pairs show significantly lower expression divergence than that of random gene pairs. In contrast to a previous finding, but consistent with the common sense, expression profiles of orthologous tissues between species are more similar to each other than to those of nonorthologous tissues. Furthermore, the evolutionary rate of expression divergence and that of coding sequence divergence are found to be weakly, but significantly positively correlated, when RA and the Euclidean distance are used to measure expression-profile divergence. These results highlight the importance of proper consideration of various estimation errors in comparing the microarray data between species.     

Evolution of the vertebrate ABC gene family: analysis of gene birth and death

Vertebrate evolution has been largely driven by the duplication of genes that allow for the acquisition of new functions. The ATP-binding cassette (ABC) proteins constitute a large and functionally diverse family of membrane transporters. The members of this multigene family are found in all cellular organisms, most often engaged in the translocation of a wide variety of substrates across lipid membranes. Because of the diverse function of these genes, their large size, and the large number of orthologs, ABC genes represent an excellent tool to study gene family evolution. We have identified ABC proteins from the sea squirt (Ciona intestinalis), zebrafish (Danio rerio), and chicken (Gallus gallus) and, using phylogenetic analysis, identified those genes with a one-to-one orthologous relationship to human ABC proteins. All ABC protein subfamilies found in Ciona and zebrafish correspond to the human subfamilies, with the exception of a single ABCH subfamily gene found only in zebrafish. Multiple gene duplication and deletion events were identified in different lineages, indicating an ongoing process of gene evolution. As many ABC genes are involved in human genetic diseases, and important drug transport phenotypes, the understanding of ABC gene evolution is important to the development of animal models and functional studies.     

Identification of the carboxyl termini of porcine zona pellucida glycoproteins ZPB and ZPC

The extracellular matrix surrounding mammalian oocytes plays important roles in fertilization and is known as the zona pellucida (ZP). The ZP consists of three glycoproteins, ZPA, ZPB, and ZPC, which contain homologous regions known as ZP domains. The ZP domain is also found in many other secretory glycoproteins. Putative transmembrane domains present at the C-termini of ZP glycoprotein precursors are removed as the proteins proceed through the secretory pathway. However, the details of this processing have been unclear. In particular, the precise locations of the C-termini of mammalian zona proteins have not yet been determined. In this study, the C-terminal residues of porcine ZPB and ZPC were identified as Ala-462 and Ser-332, respectively, by mass spectrometry of C-terminal polypeptide fragments of these proteins. These results suggest that ZPB is processed at its furin consensus site, whereas ZPC is processed N-terminal to the furin consensus site. In addition, the analyses of porcine ZPB and ZPC fragments revealed that disulfide bonds within the ZP domains are divided into two groups, suggesting that the ZP domain consists of two subdomains.     

Low rates of expression profile divergence in highly expressed genes and tissue-specific genes during mammalian evolution

Evolutionary rates provide important information about the pattern and mechanism of evolution. Although the rate of gene sequence evolution has been well studied, the rate of gene expression evolution is poorly understood. In particular, it is unclear whether the gene expression level and tissue specificity influence the divergence of expression profiles between orthologous genes. Here we address this question using a microarray data set comprising the expression signals of 10,607 pairs of orthologous human and mouse genes from over 60 tissues per species. We show that the level of gene expression and the degree of tissue specificity are generally conserved between the human and mouse orthologs. The rate of gene expression profile change during evolution is negatively correlated with the level of gene expression, measured by either the average or the highest level among all tissues examined. This is analogous to the observation that the rate of gene (or protein) sequence evolution is negatively correlated with the gene expression level. The impacts of the degree of tissue specificity on the evolutionary rate of gene sequence and that of expression profile, however, are opposite. Highly tissue-specific genes tend to evolve rapidly at the gene sequence level but slowly at the expression profile level. Thus, different forces and selective constraints must underlie the evolution of gene sequence and that of gene expression.     

A comparative map of macrochromosomes between chicken and Japanese quail based on orthologous genes

In order to develop a comparative map between chicken and quail, we identified orthologous gene markers based on chicken genomic sequences and localized them on the Japanese quail Kobe-NIBS linkage map, which had previously been constructed with amplified fragment length polymorphisms. After sequencing the intronic regions of 168 genes located on chicken chromosomes 1-8, polymorphisms among Kobe-NIBS quail family parents were detected in 51 genes. These orthologous markers were mapped on eight Japanese quail linkage groups (JQG), and they allowed the comparison of JQG to chicken macrochromosomes. The locations of the genes and their orders were quite similar between the two species except within a previously reported inversion on quail chromosome 2. Therefore, we propose that the respective quail linkage groups are macrochromosomes and designated as quail chromosomes CJA 1-8.     

Structural conservation of mouse and rat zona pellucida glycoproteins. Probing the native rat zona pellucida proteome by mass spectrometry

The mammalian zona pellucida is an egg extracellular matrix to which sperm bind. Mouse zonae are composed of three glycoproteins (ZP1, ZP2, and ZP3), while rat zonae contain four (ZP1, ZP2, ZP3, and ZP4/ZPB). Mouse sperm bind to zonae comprised solely of mouse ZP2 and ZP3. In this report, we show that rat sperm also bind to these zonae, indicating that ZP2 and ZP3 contain a "minimum structure(s)" to which rodent sperm can bind, and ZP1 and ZP4/ZPB are dispensable in these two rodents. These data are consistent with our mass spectrometric analysis of the native rat zona pellucida proteome (defined as the fraction of the total rat proteome to which the zonae glycoproteins contribute) demonstrating that the rat zonae glycoproteins share a high degree of conservation of structural features with respect to their mouse counterparts. The primary sequences of the rat zonae proteins have been deduced from cDNA. Each zona protein undergoes extensive co- and post-translational modification prior to its secretion and incorporation into an extracellular zona matrix. Each has a predicted N-terminal signal peptide that is cleaved off once protein translation begins and an anchoring C-terminal transmembrane domain from which the mature protein is released. Mass spectrometric analysis with a limited amount of native material allowed determination of the mature N-termini of rat ZP1 and ZP3, both of which are characterized by cyclization of glutamine to pyroglutamate; the N-terminus of ZP2 was identified by Edman degradation. The mature C-termini of ZP1 and ZP3 end two amino acids upstream of a conserved dibasic residue that is part of, but distinct from, the consensus furin cleavage sequence, while the C-terminus of ZP2 was not determined. Each zona protein contains a "zona domain" with eight conserved cysteine residues that is thought to play a role in the polymerization of the zona proteins into matrix filaments. Partial disulfide bond assignment indicates that the intramolecular disulfide patterns in rat ZP1, ZP2, and ZP3 are identical to those of their corresponding mouse counterparts. Last, nearly all potential N-glycosylation sites are occupied in the rat zonae glycoproteins (three of three for ZP1, six or seven of seven for ZP2, and four or five of six for ZP3). In comparison, potential O-glycosylation sites are numerous (59-83 Ser/Thr residues), but only two regions were observed to carry O-glycans in rat ZP3.     

Isolation of novel GRO genes and a phylogenetic analysis of the CXC chemokine subfamily in mammals

Approximately 15 different alpha, or CXC, chemokines have thus far been isolated from 11 species of mammals. Among the best studied chemokines are the 12 human proteins that are encoded by 11 paralogous genes. In order to better understand the evolution and function of this group of genes, we isolated and characterized six novel GRO and GRO-related cDNA sequences from the cow (Bos taurus), the sheep (Ovis aries), the rabbit (Oryctolagus cuniculus), and the guinea pig (Cavia porcellus). The amino acid sequence of the diverged guinea pig GRO or KC gene is only 50%-60% similar to presumed orthologs from other species, while the sheep and cow GRO proteins are 90%-99% similar to each other. The presence of multiple GRO genes in the cow, the rabbit, and the sheep is consistent with what has been observed for humans. Phylogenetic analyses of amino acid sequences from 44 proteins indicate that genes orthologous to many of the 11 known from humans exist in other species. One such gene, interleukin 8, or IL8, has been isolated from nine species, including the rodent guinea pig; however, this gene is absent in the rat and the mouse, indicating a unique gene loss event in the rat/mouse (muroid rodent) lineage. The KC (or MIP2) gene of rodents appears to be orthologous to the GRO gene found in other taxonomic orders. Combined evidence from different sources suggests that IP10 and MIG share sister taxon relationships on the evolutionary tree, while the remaining paralogous genes represent independent lineages, with limited evidence for kinship between them. This observation indicates that these genes originated nearly contemporaneously via a series of gene duplication events. Relative-rate tests for synonymous and nonsynonymous nucleotide substitutions in the KC and IL8 genes did not detect rate heterogeneity; however, there are several notable features regarding the IL8 genes. For example, the IL8 proteins from two Old World monkeys are as similar to one another as they are to the IL8 protein from humans, and all observed nucleotide differences between the IL8 genes of the two monkeys cause amino acid changes; in other words, there are no synonymous differences between them.