A novel gene member of the human glycophorin A and B gene family. Molecular cloning and expression

A new gene closely related to the glycophorin A (GPA) and glycophorin B (GPB) genes has been identified in the normal human genome as well as in that of persons with known alterations of GPA and/or GPB expression. This gene, called glycophorin E (GPE), is transcribed into a 0.6-kb message which encodes a 78-amino-acid protein with a putative leader peptide of 19 residues. The first 26 amino acids of the mature protein are identical to those of M-type glycophorin A (GPA), but the C-terminal domain (residues 27-59) differs significantly from those of glycophorins A and B (GPA and GPB). The GPE gene consists of four exons distributed over 30 kb of DNA, and its nucleotide sequence is homologous to those of the GPA and GPB genes in the 5' region, up to exon 3. Because of branch and splice site mutations, the GPE gene contains a large intron sequence partially used as exons in GPA and GPB genes. Compared to its counterpart in the GPB gene, exon 3 of the GPE gene contains several point mutations, an insertion of 24 bp, and a stop codon which shortens the reading frame. Downstream from exon 3, the GPE and the GPB sequences are virtually identical and include the same Alu repeats. Thus, it is likely that the GPE and GPB genes have evolved by a similar mechanism. From the analysis of the GPA, GPB and GPE genes in glycophorin variants [En(a-), S-s-U- and Mk], it is proposed that the three genes are organized in tandem on chromosome 4. Deletion events within this region may remove one or two structural gene(s) and may generate new hybrid structures in which the promoter region of one gene is positioned upstream from the body of another gene of the same family. This model of gene organization provides a basis with which to explain the diversity of the glycophorin gene family.     

Sequence diversification and exon inactivation in the glycophorin A gene family from chimpanzee to human

In humans, the allelic diversity of MNSs glycophorins (GP) occurs mainly through the recombinational modulation of silent exons (pseudoexons) in duplicated genes. To address the origin of such a mechanism, structures of GPA, GPB, and GPE were determined in chimpanzee, the only higher primate known to have achieved a three-gene framework as in humans. Pairwise comparison of the chimpanzee and human genes revealed a high degree of sequence identity and similar exon-intron organization. However, the chimpanzee GPA gene lacks a completely formed M- or N-defining sequence as well as a consensus sequence for the Asn-linked glycosylation. In the case of the GPB gene, exon III is expressed in the chimpanzee but silenced, as a pseudoexon, in the human. Therefore, the protein product in the chimpanzee bears a larger extracellular domain than in the human. For the GPE genes, exon III and exon IV have been inactivated by identical donor splice-site mutations in the two species. Nevertheless, the chimpanzee GPE-like mRNA appeared to be transcribed from a GPB/E composite gene containing no 24-bp insertion sequence in exon V for the transmembrane domain. These results suggest a divergent processing of exonic units from chimpanzee to human in which the inactivation of GPB exon III preserved a limited sequence repertoire for diversification of human glycophorins.Correspondence to: O.O. Blumenfeld     

Molecular analysis of a hybrid gene encoding human glycophorin variant Miltenberger V-like molecule

The genomic structure of a human glycophorin variant, Miltenberger class V-like molecule (MiV*), was examined. Southern blot analysis of total genomic DNA revealed that the 5' half of the MiV* gene derived from glycophorin A (GPA) gene whereas the 3' half derived from glycophorin B (GPB) gene. This structure is reciprocal to another glycophorin variant, Sta, which has a GPB-GPA hybrid structure. The genomic sequences around the crossing-over point were amplified by polymerase chain reaction, and the sequences were determined. Comparison of the nucleotide sequences of the GPA, GPB, and MiV* genes indicates that the crossing-over point is located in the region around the 3' end of intron 3 of the GPA gene. This place is different from the crossing-over point for Sta, which was found to be highly homologous to that for haptoglobin-related genes. However, the nucleotide sequences within the presumptive crossing-over point for the MiV* gene were found to be homologous in a reverse orientation to the crossing-over point proposed for haptoglobin-related genes. These results suggest strongly that homologous recombination through unequal crossing over can be facilitated by specific genomic elements such as those in common for formation of MiV*, Sta, and haptoglobin-related genes. The present study also localized the gene of the third glycophorin, GPE, at chromosome 4, q31.1 band, the same locus as for the GPA and GPB genes. The results indicate that GPE was not involved in generating MiV* or Sta hybrid gene despite the fact that it is localized adjacent to the GPA and GPB genes.     

Evidence for diversifying selection on erythrocyte-binding antigens of Plasmodium falciparum and P. vivax

Malaria parasite antigens involved in erythrocyte invasion are primary vaccine candidates. The erythrocyte-binding antigen 175K (EBA-175) of Plasmodium falciparum binds to glycophorin A on the human erythrocyte surface via an N-terminal cysteine-rich region (termed region II) and is a target of antibody responses. A survey of polymorphism in a malaria-endemic population shows that nucleotide alleles in eba-175 region II occur at more intermediate frequencies than expected under neutrality, but polymorphisms in the homologous domains of two closely related genes, eba-140 (encoding a second erythrocyte-binding protein) and psieba-165 (a putative pseudogene), show an opposite trend. McDonald-Kreitman tests employing interspecific comparison with the orthologous genes in P. reichenowi (a closely related parasite of chimpanzees) reveal a significant excess of nonsynonymous polymorphism in P. falciparum eba-175 but not in eba-140. An analysis of the Duffy-binding protein gene, encoding a major erythrocyte-binding antigen in the other common human malaria parasite P. vivax, also reveals a significant excess of nonsynonymous polymorphisms when compared with divergence from its ortholog in P. knowlesi (a closely related parasite of macaques). The results suggest that EBA-175 in P. falciparum and DBP in P. vivax are both under diversifying selection from acquired human immune responses.     

Molecular analysis of glycophorin A and B gene structure and expression in homozygous Miltenberger class V (Mi. V) human erythrocytes

In the Miltenberger class V (Mi. V) condition, red cells lack glycophorin A (GPA) and glycophorin B (GPB) but carry instead an unusual glycoprotein thought to be a hybrid molecule produced by the unequal crossing-over between the closely linked genes encoding for GPA and GPB. By Western blot analysis with rabbit anti-GPA antibodies specific for discrete domains of GPA, it was found that the Mi. V glycoprotein (donor F. M.) contains approximately 60 amino acid residues of GPA at its N-terminus. As a preliminary approach to the molecular analysis of this variant the restriction maps of the GPA and GPB genes were established by Southern blot analysis of genomic DNA and from genomic clones isolated from a human leukocyte library constructed in lambda EMBL4. The GPA and GPB genes cover about 30 kb of DNA and are organized into seven exons (A-1-A-7) and five exons (B-1-B-5), respectively. In addition to the normal genes, a third gene (named inv), closely resembling the GPA and GPB genes, was also identified. In the homozygous Mi. V individual the normal GPA and GPB genes were absent, but an unusual form of gene structure was detected by Southern blot analysis. The Mi. V glycoprotein gene was composed of exon B-1 of the GPB gene followed by exons A-2 and A-3 of the GPA gene and the exons B-3, B-4 and B-5 of the GPB gene. Exon B-1 can be distinguished from exon A-1 of GPA since it is located within a different restriction fragment, but both encode the same amino acid sequence (N-terminal region of the signal peptides). Using the polymerase chain reaction, the junction between exon A-3 and exon B-3 was confirmed by amplification of the DNA region where the putative crossing-over has occurred and it was deduced that the Mi. V glycoprotein is a hybrid molecule composed of amino acid residues 1-58 from GPA fused to amino acid residues 27-72 of GPB. In addition, the finding that part of the signal peptide and the 5'-untranslated region are derived from GPB suggests that the genetic background of the Mi. V variant is rather complex and may involve a cascade of recombination or gene conversion events.     

Genetic polymorphism and natural selection in the malaria parasite Plasmodium falciparum.

We have studied the genetic polymorphism at 10 Plasmodium falciparum loci that are considered potential targets for specific antimalarial vaccines. The polymorphism is unevenly distributed among the loci; loci encoding proteins expressed on the surface of the sporozoite or the merozoite (AMA-1, CSP, LSA-1, MSP-1, MSP-2, and MSP-3) are more polymorphic than those expressed during the sexual stages or inside the parasite (EBA-175, Pfs25, PF48/45, and RAP-1). Comparison of synonymous and nonsynonymous substitutions indicates that natural selection may account for the polymorphism observed at seven of the 10 loci studied. This inference depends on the assumption that synonymous substitutions are neutral, which we test by analyzing codon bias and G+C content in a set of 92 gene loci. We find evidence for an overall trend towards increasing A+T richness, but no evidence for mutation bias. Although the neutrality of synonymous substitutions is not definitely established, this trend towards an A+T rich genome cannot explain the accumulation of substitutions at least in the case of four genes (AMA-1, CSP, LSA-1, and PF48/45) because the Gleft and right arrow C transversions are more frequent than expected. Moreover, the Tajima test manifests positive natural selection for the MSP-1 and, less strongly, MSP-3 polymorphisms; the McDonald-Kreitman test manifests natural selection at LSA-1 and PF48/45. We conclude that there is definite evidence for positive natural selection in the genes encoding AMA-1, CSP, LSA-1, MSP-1, and Pfs48/45. For four other loci, EBA-175, MSP-2, MSP-3, and RAP-1, the evidence is limited. No evidence for natural selection is found for Pfs25.     

Studies on human red-cell membrane glycophorin A and glycophorin B genes in glycophorin-deficient individuals.

1. Genomic DNA derived from individuals who lack glycophorin A (GPA), glycophorin B (GPB) or both of these proteins was subjected to Southern-blot analysis using GPA and GPB cDNA probes. 2. Bands on the Southern blots were assigned to the GPA gene, GPB gene or to a putative pseudogene. 3. Genomic DNA derived from an individual of the Mk phenotype was shown to have deletions in the GPA and GPB genes. The simplest model for the results obtained is that a single deletion spans the GPA and GPB genes in the individual studied.     

一种鉴定MNSs红细胞血型单克隆抗体类型的方法的建立

为鉴定MNSs血型单克隆细胞株6D7C9分泌的抗体类型,通过克隆、亚克隆、细胞转染等分子生物学技术建立了血型糖蛋白GPA、GPB的异源表达系统,并将其作为抗原,通过ELISA、Western 印迹法确定6D7C9分泌的McAb.结果显示,RT－PCR技术成功克隆获得了GPA、GPB血型糖蛋白编码基因,通过分别构建其重组逆转录病毒表达载体pEGZ/GPA及pEGZ/GPB,转染包装细胞293T,再感染L929细胞,经zeocin筛选2周后,RT PCR及流式细胞仪分析证实,L929/GPA和L929/GPB转基因细胞中分别有GPA、GPB目的基因的转录和蛋白表达.用稳定高表达GPA、GPB的转基因细胞通过ELISA和Western 印迹法证实单克隆细胞株6D7C9分泌的是抗GPA/GPB McAb.本研究成功地建立了血型糖蛋白GPA、GPB的异源表达系统,为MNSs血型McAb的检测及GPA、GPB蛋白的功能学研究奠定了基础.     

Pattern of Nucleotide Substitutions in Growth Hormone-Prolactin Gene Family: A Paradigm for Evolution by Gene Duplication

The growth hormone-prolactin gene family in mammals is an interesting example of evolution by gene duplication. Divergence among members of duplicated gene families and among species was examined by using reported gene sequences of growth hormone, prolactin and their receptors. Sequence divergence among species was found to show a general tendency in which a generation-time effect is pronounced for synonymous substitutions but not so for nonsynonymous substitutions. Divergence among duplicated genes is characterized by the relatively high rate of nonsynonymous substitutions, i.e., the rate is close to that of synonymous ones. In view of the stage- and tissue-specific expression of duplicated genes, some of the amino acid substitutions among duplicated genes is likely to be caused by positive Darwinian selection.     

Identification of the crossing-over point of a hybrid gene encoding human glycophorin variant Sta. Similarity to the crossing-over point in haptoglobin-related genes

One of the human glycophorin variants, Stones (Sta), has been shown to be the product of a hybrid gene of which the 5'-half derived from the glycophorin B (GPB) gene whereas the 3'-half derived from the glycophorin A (GPA) gene. The present study reveals the crossing-over point of this hybrid gene from the analysis of polymerase chain reaction products. The genomic sequences encompassing the region corresponding to exon 3 to exon 4 of GPA were amplified by polymerase chain reaction with oligonucleotide primers synthesized according to GPA and GPB genomic sequences (Kudo, S., and Fukuda, M. (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 4619-4623). After subcloning the products, the nucleotide sequences derived from GPA, GPB, and putative Sta genes were determined. Comparison of the nucleotide sequences of GPA, GPB, and Sta genes indicate that the crossing-over took place 200 base pairs upstream from the first nucleotide of exon 4. Intriguingly, the nucleotide sequence surrounding the putative crossing-over point is homologous to the crossing-over point proposed for haptoglobin genes (Maeda, N., McEvoy, S.M., Harris, H.F., Huisman, T.H.J., and Smithies, O. (1986) Proc. Natl. Acad. Sci. U.S.A. 83, 7395-7399). These results suggest strongly that homologous recombination through unequal crossing-over can be facilitated by specific genomic elements, such as those in common in these two crossing-over events. The present study also revealed that this Sta individual has a variant GPA gene; substitution of adenine for guanine at the nucleotide for codon 39 results in substitution of lysine for arginine at amino acid 39, and loss of an SstI restriction site.     

Molecular evolution of the major capsid protein VP1 of enterovirus 70.

Nucleotide sequences of the genome RNA encoding capsid protein VP1 (918 nucleotides) of 18 enterovirus 70 (EV70) isolates collected from various parts of the world in 1971 to 1981 were determined, and nucleotide substitutions among them were studied. The genetic distances between isolates were calculated by the pairwise comparison of nucleotide difference. Regression analysis of the genetic distances against time of isolation of the strains showed that the synonymous substitution rate was very high at 21.53 x 10(-3) substitution per nucleotide per year, while the nonsynonymous rate was extremely low at 0.32 x 10(-3) substitution per nucleotide per year. The rate estimated by the average value of synonymous and nonsynonymous substitutions (W.-H. Li, C.-C. Wu, and C.-C. Luo, Mol. Biol. Evol. 2:150-174, 1985) was 5.00 x 10(-3) substitution per nucleotide per year. Taking the average value of synonymous and nonsynonymous substitutions as genetic distances between isolates, the phylogenetic tree was inferred by the unweighted pairwise grouping method of arithmetic average and by the neighbor-joining method. The tree indicated that the virus had evolved from one focal place, and the time of emergence was estimated to be August 1967 +/- 15 months, 2 years before first recognition of the pandemic of acute hemorrhagic conjunctivitis. By superimposing every nucleotide substitution on the branches of the phylogenetic tree, we analyzed nucleotide substitution patterns of EV70 genome RNA. In synonymous substitutions, the proportion of transitions, i.e., C<==>U and G<==>A, was found to be extremely frequent in comparison with that reported on other viruses or pseudogenes. In addition, parallel substitutions (independent substitutions at the same nucleotide position on different branches, i.e., different isolates, of the tree) were frequently found in both synonymous and nonsynonymous substitutions. These frequent parallel substitutions and the low nonsynonymous substitution rate despite the very high synonymous substitution rate described above imply a strong restriction on nonsynonymous substitution sites of VP1, probably due to the requirement for maintaining the rigid icosahedral conformation of the virus.     

New procedures for glycophorin A purification with high yield and high purity

Glycophorin A (GPA) is the major glycoprotein of the human erythrocyte membrane. It is known to form, in SDS gels as well as in a membrane environment, homodimers, and also heterodimers with the homologous molecule Glycophorin B (GPB). It is shown in this report that the propensity of GPA to dimerize with GPB precludes satisfactory preparation with high yield of pure GPA using classical techniques including SEC and RPLC. It was demonstrated using multiple angle light scattering that GPA is eluted from RPLC columns as dimers. A convenient procedure was devised which allowed us to get pure GPA with high yield. This procedure consists of selectively blocking GPA–GPB heterodimer formation by selective modification of Cysteine 50 of GPB before RPLC.     

A new method for estimating synonymous and nonsynonymous rates of nucleotide substitution considering the relative likelihood of nucleotide and codon changes

A new method is proposed for estimating the number of synonymous and nonsynonymous nucleotide substitutions between homologous genes. In this method, a nucleotide site is classified as nondegenerate, twofold degenerate, or fourfold degenerate, depending on how often nucleotide substitutions will result in amino acid replacement; nucleotide changes are classified as either transitional or transversional, and changes between codons are assumed to occur with different probabilities, which are determined by their relative frequencies among more than 3,000 changes in mammalian genes. The method is applied to a large number of mammalian genes. The rate of nonsynonymous substitution is extremely variable among genes; it ranges from 0.004 X 10(-9) (histone H4) to 2.80 X 10(-9) (interferon gamma), with a mean of 0.88 X 10(-9) substitutions per nonsynonymous site per year. The rate of synonymous substitution is also variable among genes; the highest rate is three to four times higher than the lowest one, with a mean of 4.7 X 10(-9) substitutions per synonymous site per year. The rate of nucleotide substitution is lowest at nondegenerate sites (the average being 0.94 X 10(-9), intermediate at twofold degenerate sites (2.26 X 10(-9)). and highest at fourfold degenerate sites (4.2 X 10(-9)). The implication of our results for the mechanisms of DNA evolution and that of the relative likelihood of codon interchanges in parsimonious phylogenetic reconstruction are discussed.     

Mitochondrial DNA polymorphisms as risk factors for Parkinson’s disease and Parkinson’s disease dementia

The activity of complex I of the mitochondrial respiratory chain has been found to be decreased in patients with Parkinsons disease (PD), but no mutations have been identified in genes encoding complex I subunits. Recent studies have suggested that polymorphisms in mitochondrial DNA (mtDNA)-encoded complex I genes (MTND) modify susceptibility to PD. We hypothesize that the risk of PD is conveyed by the total number of nonsynonymous substitutions in the MTND genes in various mtDNA lineages rather than by single mutations. To test this possibility, we determined the number of nonsynonymous substitutions of the seven MTND genes from 183 Finns. The differences in the total number of nonsynonymous substitutions and the nonsynonymous to synonymous substitution rate ratio (K_a/K_s) of MTND genes between the European mtDNA haplogroup clusters (HV, JT, KU, IWX) were analysed by using a statistical approach. Patients with PD (n=238) underwent clinical examination together with mtDNA haplogroup analysis and the clinical features between patient groups defined by the number of nonsynonymous substitutions were compared. Our analysis revealed that the haplogroup clusters HV and KU had a lower average number of amino acid replacements and a lower K_a/K_s ratio in the MTND genes than clusters JT and IWX. Supercluster JTIWX with the highest number of amino acid replacements was more frequent among PD patients and even more frequent among patients with PD who developed dementia. Our results suggest that a relative excess of nonsynonymous mutations in MTND genes in supercluster JTWIX is associated with an increased risk of PD and the disease progression to dementia.     

Distinct regions of human glycophorin A enhance human red cell anion exchanger (band 3; AE1) transport function and surface trafficking

Human red cell glycophorin A (GPA) enhances the expression of band 3 anion transport activity at the cell surface of Xenopus oocytes. This effect of GPA could occur in two ways, enhancement of band 3 anion transport function or enhancement of band 3 trafficking to the cell surface. We have examined the GPA effect using GPA mutants. We compared the sequences of GPA and its homolog glycophorin B (GPB; which does not facilitate band 3 cell-surface activity or trafficking) to identify candidate regions of GPA for study. We constructed several GPA or GPB mutants, including naturally occurring GPA/GPB hybrid molecules and insertion, deletion, and substitution mutants. We analyzed the effects of the mutant proteins on band 3-specific chloride transport and surface presentation using co-expression in Xenopus oocytes. We find that the C-terminal cytoplasmic tail of GPA enhances trafficking of band 3 to the cell surface, whereas the extracellular residues 68-70 increase the specific anion transport activity of band 3. In addition, examination of the oligomerization of GPA mutants showed that single amino acid substitutions N-terminal to the transmembrane domain greatly reduce SDS-stable GPA dimer formation, implying that regions outside the transmembrane domain of GPA are important for GPA dimer formation.     

Evolution of the Sry genes

Existing DNA sequence data on the Sry gene, the mammalian sex- determining locus in the Y chromosome, were analyzed for primates, rodents, and bovids. In all three taxonomic groups, the terminal sequences evolved faster than the HMG (high mobility group) boxes, and this applies both to synonymous (Ks) and nonsynonymous (Ka) nucleotide substitutions. Similar intragenic correlation between synonymous and nonsynonymous substitution rates was not found either in other mammalian genes that contain a conservative box (Sox, Msx) or in the MADS-box genes of plants. The rate of nonsynonymous substitutions exceeds significantly that of synonymous substitutions in the terminal Sry sequences of apes. We did not find good support for the hypothesis that the high evolutionary rate of Sry would be associated with a promiscuous mating system.