Molecular Biology and Evolution 19:223-229 (2002)
© 2002 Society for Molecular Biology and Evolution
Natural Selection on the Erythrocyte Surface
*Department of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine;
Institute of Biological Anthropology, University of Oxford
 |
Abstract |
|---|
 TOP Abstract IntroductionMaterials and MethodsResultsDiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Surface glycoproteins are principal receptors used by pathogens to invade target cells. It has been suggested that mammalian erythrocyte surface glycoproteins function as decoy receptors attracting pathogens to the anucleated erythrocyte and away from their target tissues. Glycophorin A (GYPA) is solely expressed on the erythrocyte surface where it is the most abundant sialoglycoprotein, although its function is unknown. The pathogen decoy hypothesis may be relevant here, as GYPA has been shown in vitro to bind numerous viruses and bacteria, which do not infect erythrocytes. However, it is also a receptor for erythrocyte invasion by the malarial parasite Plasmodium falciparum. Analyses of gypa sequence variation among six higher primates and within a human population show that there is a large excess of replacement (nonsynonymous) substitutions along each primate lineage (particularly on exons 2–4 encoding the extracellular glycosylated domain of GYPA) and a significant excess of polymorphisms in exon 2 (encoding the terminal portion of the extracellular domain) within humans. These two signatures suggest that there has been exceptionally strong positive selection on this receptor driving GYPA divergence during primate evolution and balancing selection maintaining allelic variation within human populations. The pathogen decoy hypothesis alone is adequate to explain both these signatures of between-species and within-species diversifying selection. This has implications for understanding the functions of erythrocyte surface components and their roles in health and disease.
|
Introduction |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Host-pathogen interactions are major determinants in evolution, with cell surface receptors likely to be under constant selection to protect tissues from infection. Erythrocytes express an unexplained density and diversity of glycoprotein structures on their surface (Mourant, Kopec, and Domaniewska-Sobczak 1978
, Pp. 1–22) which cannot be explained by their primary function of transporting oxygen. Glycophorin A (GYPA) is the most abundant erythrocyte surface sialoglycoprotein (Cartron and London 1992
; Huang and Blumenfeld 1995
) with 0.5–1.0 million copies on each cell. The gypa gene is located on chromosome 4q28–31, in tandem with two paralogous genes (gypb and gype) which arose by duplication during primate evolution (Rearden et al. 1993
). It is expressed solely on cells of the erythroid lineage (Huang and Blumenfeld 1995
), and variation in its extracellular portion in humans determines the two major alleles of the MN blood group system, with further genetic changes accounting for rare blood group variants (Huang and Blumenfeld 1995
). However, its function is still unknown. The dense negative charge that the sugar groups give the erythrocyte surface may play a role in preventing adhesion between cells and vascular surfaces (Huang and Blumenfeld 1995
) having implications for blood viscosity; alternatively there are suggestions that it may be involved in cell ageing (Aminoff 1988
). The extremely rare individuals who are GYPA-negative show that the null phenotype is viable (Cartron and London 1992
), although its rarity would seem to indicate some cost associated with an absence of the receptor.
Many pathogens use sugar groups of surface glycoproteins, in particular sialic acid residues, as receptors for invasion (Karlsson 1995
). It has therefore been suggested that erythrocyte glycoproteins may function as decoy receptors (Gagneux and Varki 1999
) attracting pathogens to the erythrocyte, which as a mature cell lacks a nucleus and DNA replication, and away from target tissues (Gagneux and Varki 1999
). This represents a plausible function for GYPA which carries many complex sugar groups (Tomita and Marchesi 1975
; Pisano et al. 1993
) and is supported by studies in which GYPA has been shown to bind numerous viruses (Paul and Lee 1987
; Nishimura et al. 1988
; Tavakkol and Burness 1990
; Wybenga et al. 1996
) and bacteria (Baseman, Banai, and Kahane 1982
; Brooks et al. 1989
; Saada et al. 1991
). GYPA is, however, also known to be a principal receptor for the malarial parasite Plasmodium falciparum, which infects erythrocytes (Pasvol, Wainscoat, and Weatherall 1982
; Sim et al. 1994
). If GYPA generally functions as a decoy receptor, then it is likely to be under selection for targeting the diversity of potential pathogens. It is expected that this would cause diversifying selection on the gypa gene sequence.
To investigate the selective forces affecting the evolution of the GYPA receptor, a survey of sequence variation in the gypa gene was carried out. Tests for departures from neutrality were performed on pairwise comparisons of synonymous and nonsynonymous differences in the region encoding the mature GYPA product among six primate species and on genetic diversity in the region encoding the extracellular portion of GYPA in a human population.
|
Materials and Methods |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Sequences of gypa and Other Nuclear Genes from Different Primates
The available gypa coding sequences (exons 2–6) of six primates were studied: human, chimpanzee, gorilla, orangutan, gibbon, and macaque (GenBank accession numbers AJ309706, AJ309708, AF015174, AF023467, AF023468, and AF023469). Exon sequences for nine other protein-coding loci in the same six higher primates were obtained from the following GenBank accession numbers—CC chemokine receptor-5 (ccr5): (HSU83326, AF177894, AF177901, AF177900, AF177884, AF177890).
(1,2) fucosyltransferase (H) (fut1): (M35531, AF080603, AF080605, AF111935, AF045545, AF080607). Lysozyme (lyz) (J03801, U76912–5, X60236). A
globin and G globin (A, G): (U01317, M92294–6, J05174, X53419). Protamine 2 (prm2): (NM_002762, X72968, X71336–9). Rhesus blood group associated glycoprotein (rhag): (NM_000324, AF177621–5). Tyrosinase (tyr): (M27160, AF183588–92, AF183599–603, AF183604–8, AF183629–33, AF183644–48). Zinc finger protein, Y linked (zfy): (NM_003411, AB041908, AB041910, AB041912, AB041914, AB041918).
Allele Sequence of gypa Exons 2–4 Sampled from a Human Population
Purified genomic DNA from 33 unrelated individuals (Yoruba ethnicity, southwestern Nigeria) and a single chimpanzee (sample kindly provided by N. Mundy) were used as templates for amplification of a 3.6-kb region of gypa, from exon 2 to the boundary of exon 5 (see fig. 1
). gypa-specific amplification primers GYPAIfwd5'-GCTTAGCTCAGGGACTGGAGG-3' and GYPAIIIrev5'-CACCTTGCCTTTTAATAGAAAGC-3' were designed to ensure no amplification of a paralogous sequence in gypb or gype, which share strong homology with gypa. PCR products were ligated into pGem®-T Easy Vector (Promega) and transformed by heat shock into JM 109 E. Coli High Efficiency Competent Cells. Isolated plasmids were prepared using a Plasmid Miniprep column (QIAGEN). A single clone representing one allele was isolated and sequenced from each individual to avoid problems associated with determination of haplotypes in heterozygotes. Each template was sequenced using internal sequencing primers (available on request) with 3' BIG DYETM dye terminator cycle sequencing premix kit (Applied Biosystems). Sequencing was carried out on a Perkin-Elmer ABI PrismTM 377 DNA Sequencer (Applied Biosystems) and sequences were checked and assembled using Sequence Navigator Version 1.0.1 (Applied Biosystems). Polymorphic sites were identified, and all singletons in the data set were confirmed by reamplification from genomic DNA and sequencing to exclude any errors introduced by PCR or cloning (deposited human sequences: GenBank accession numbers: AJ309828–45, AJ311318–32; chimpanzee AJ309708).
|
Statistical Analyses of Between-species and Within-species Diversity
The ratio of nonsynonymous substitions per nonsynonymous site to synonymous substitutions per synonymous site (dN/dS) was calculated for all pairs of the six primate species (pairwise analysis) using the method of Nei and Gojobori (Nei and Gojobori 1986
) by the program MEGA 2.1 (http://www.megasoftware.net/) (Kumar, Tamura, and Nei 1994
), with standard errors for dN and dS calculated using a bootstrap method in the MEGA program. Maximum likelihood estimates of dN/dS (the 
parameter) for each lineage in the phylogeny (see fig. 2
) were derived by a method employing different evolutionary models, using the CODEML program of PAML Version 3.0b (http://abacus.gene.ucl.ac.uk/software/paml.html) (Nielsen and Yang 1998
; Yang 1998
; Yang et al. 2000
). Investigation of evolutionary rates between lineages was carried out using alternative likelihood models, one with a single dN/dS ratio (M0, one-ratio) (Goldman and Yang 1994
) estimated for all branches and another that allowed independent dN/dS ratios for each branch (FR, free-ratio) (Yang 1998
). Statistical testing of differences between the different models (Nielsen and Yang 1998
; Yang 1998
; Yang et al. 2000
) was carried out using twice the log likelihood difference (2
l) which conforms to a chi-square distribution (the Likelihood Ratio test, LRT), with the degrees of freedom based on the difference between the number of parameters estimated from the models.
|
Three tests of neutrality (Hudson, Kreitman, and Aguade 1987
; Tajima 1989
; Fu and Li 1993
) based on analyses of patterns of intraspecific and interspecific nucleotide differences, using the human population data set and a chimpanzee sequence out-group were performed using the program DNAsp Version 3.5 (http://www.bio.ub.es/
julio/DnaSP.html) (Rozas and Rozas 1999
).
Prediction of O-linked Glycosylation Sites in Primate GYPA Sequences
Published primate gypa sequences (as above) were translated into amino acid sequences and entered into NetOGlyc Version 2.0 (Hansen et al. 1998
) (www.cbs.dtu.dk/services/NetOGlyc). This predicts O-glycosylation sites on serine and threonine residues in peptides based on the amino acid sequence context.
|
Results |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Analysis of Sequence Differences in the gypa Gene Between Primate Species
The sequences encoding the mature GYPA product (378 bp, exons 2–6; fig. 1 ) in humans and five other species (chimpanzee, gorilla, orangutan, gibbon, and macaque) were studied. All the 15 pairwise ratios of dN/dS (number of nonsynonymous differences per nonsynonymous site vs. the number of synonymous differences per synonymous site) (Nei and Gojobori 1986
), among the six species were greater than 1.0 (mean 2.1, minimum 1.4, and maximum 3.3) (table 1
). This indicates an excess of amino acid replacement changes compared with neutral expectations.
|
To further examine this excess of nonsynonymous change across primate lineages, a maximum likelihood approach, which takes into account transition or transversion rate bias or nonuniform codon usage (Nielsen and Yang 1998
; Yang 1998
; Yang et al. 2000
), was used to estimate the dN/dS ratio for each separate lineage in the phylogeny. Different likelihood models were tested, including one with a fixed dN/dS ratio (one-ratio) across all branches (Goldman and Yang 1994
) and a second with independent ratios (free-ratio) for each branch (Yang 1998
) (fig. 2
). Both models gave dN/dS ratios well in excess of 1.0 (one-ratio = 2.6, see fig. 2
for the free-ratios which are >1.0 for individual branches). The likelihood of these two models did not differ significantly (log likelihood value for one-ratio model = -1064.0, for free-ratio model = -1060.2; 
2df = 8 = 8.88, P = 0.47), indicating that positive selection has operated broadly to cause divergence of the GYPA protein in different primate lineages. Furthermore, the dN/dS ratio for the part of the sequence encoding the extracellular region is particularly high (exons 2–4: maximum likelihood dN/dS = 4.7, average pairwise dN/dS = 2.4) compared with the transmembrane and intracellular region (exon 5–6: maximum likelihood dN/dS value = 1.4, average pairwise dN/dS = 1.6).
Analysis of Nine Other Nuclear Genes in Primates
A comparable analysis of the nine other nuclear genes, for which coding sequences in the same six primates were available (CCR5, fut1, A, G, lyz, prm2, rhag, tyr, and zfy), indicated the exceptional nature of GYPA diversity in primates. None of these nine genes have an estimated excess of nonsynonymous changes in all lineages. All have average pairwise dN/dS ratios of less than 1.0 (except for rhag for which dN/dS = 1.00) (fig. 3a
). As indicated in figure 3b
the average dN pairwise difference for gypa (0.113) is exceptionally high (mean for the other genes = 0.018). Hence, the high dN/dS ratio is caused by a high rate of nonsynonymous substituion and not by a low rate of synonymous substitution (the average dS pairwise difference for gypa (0.055) is typical; mean for the other genes = 0.059). Together these data highlight the distinctive signature of diversifying selection on gypa across the evolution of higher primates. Comparable sequences are not yet available for other mammalian orders, except for the mouse homologue that is so divergent from primate sequences that a proper alignment for analysis is not possible (Matsui, Natori, and Obinata 1989
).
|
Analysis of Sequence Variation in the gypa Gene Within a Human Population
To investigate whether intraspecific diversifying selection also operates within humans, the nucleotide sequence of a 3.6-kb region spanning the extracellular encoding part of the gypa gene (exons 2–4 and intron regions; fig. 1 ) was determined for one allele of each of the 33 unrelated Nigerian individuals. Among the 33 chromosomes, there were 25 different alleles of the complete sequence and three different alleles when considering only the exon coding sequence (fig. 4 ). There are five polymorphic sites in the coding sequence (all in exon 2), two of which are the nonsynonymous substitutions that determine the MN polymorphism (Huang and Blumenfeld 1995
).
|
is a polymorphic site within humans, 
is a fixed difference between species. Glycosylation sites (Pisano et al. 1993
for O-linked oligosaccharides and for N-linked oligosaccharidesAnalysis of variation among the 33 sampled chromosomes (table 2 ) showed a marked excess of nucleotide diversity in exon 2 (
= 0.024) compared with other exons ( = 0) and introns ( = 0.003). Statistical departures from neutral expectations were first tested using Tajima's (Tajima 1989
) and Fu and Li's (Fu and Li 1993
) methods. Tajima's D test (Tajima 1989
) tests for a departure from neutrality as measured by the difference between (observed average pairwise diversity) and 
(expected nucleotide diversity under neutrality derived from the number of segregating sites, S). Under balancing selection, rare alleles are selected for and maintained at intermediate frequencies, thus elevating and, therefore, making the value of the test statistic (D) more positive. Fu and Li's test (Fu and Li 1993
) reveals an excess or lack of singletons by comparing estimates of based on the number of mutations in external branches of a phylogeny versus an estimate from S (the D index) or (the F index). Absence of singletons under balancing selection will result in a positive value for each statistic. The value of Tajima's D for exon 2 is significantly positive (D = 2.54; table 2 ), indicating that alleles tend toward intermediate frequencies in the population. Fu and Li's F value for exon 2 is also significantly positive (F = 1.86; table 2
), reflecting the absence of singletons among the nucleotides in this exon. The positive values of both these statistics support a hypothesis of selection, maintaining diversity within humans in exon 2 and contrast with the negative values seen elsewhere in the gene (table 2
) and at other human genetic loci (Kaessmann et al. 2001)
.
|
Comparison of Human gypa Diversity with Divergence from Chimpanzee gypa
Alignment of the human alleles with the chimpanzee sequence for exons 2–4 (fig. 4 ) reveals that the human nucleotide polymorphisms are limited to the first half of exon 2 (up to codon 12). By marked contrast, the fixed differences between the two species all occur from codon 24 onwards (fig. 4 ). Heterogeneity between exon 2 and exons 3–4 was tested by measuring interspecific divergence (human vs. chimp) and intraspecific polymorphism (in humans) (table 3 ) using the HKA test (Hudson, Kreitman, and Aguade 1987
). The HKA test (Hudson, Kreitman, and Aguade 1987
) tests for heterogeneity between two genetic loci by measuring their divergence between two sister species and comparing this with polymorphism within the same loci in one of the species. Under neutrality, the degree of polymorphism and divergence for the two loci is expected to be correlated. Here, the hypothesis that exon 2 is evolving at the same rate as exons 3–4 is rejected (P = 0.02; 2 test). This indication of non-neutrality in exon 2, in itself and in combination with Tajima's and Fu and Li's tests above, indicates that there is selective maintenance of diversity in exon 2 encoding the outermost external portion of GYPA in humans.
|
|
Discussion |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Molecular diversity in cell surface receptors may be adaptive, if it provides a means to prevent the binding of specific ligands to the cell (escape variation) or a means to actively bind specific ligands (target-seeking). Given that the human malarial parasite P. falciparum uses GYPA as a receptor for erythrocyte invasion, it could be an agent of selection, but it is not likely to be the cause of the non-neutral pattern shown here, for two reasons. Firstly, the nonhuman primate malarial parasites do not use GYPA as a receptor (where studied primate malarial parasites use the Duffy antigen as a principal receptor for invasion [Miller et al. 1977
]), and therefore, could not apply selection broadly in primates. Secondly, if P. falciparum did discriminate allelic structures of GYPA in humans (evidence suggests it does not [Pasvol, Wainscoat, and Weatherall 1982
; Binks et al. 2001
]), it is more likely that this would give an advantage of one allele over another rather than the selective maintenance of multiple alleles as seen here. An alternative hypothesis is that GYPA may be exploited by pathogens as a Trojan horse to carry them on the erythrocyte in the circulation to target cells. This would be analogous to HIV binding to dendritic cells, which transmit the virus directly to the susceptible cells in the lymph nodes, the major site of HIV replication (Pohlmann et al. 2001a, 2001b
). However, such a hypothesis is again unable to explain the selective maintenance of alleles in humans.
The joint signatures of between-species and within-species diversifying selection can, however, be explained by the ability of GYPA to bind nonerythrocytic pathogens to the erythrocyte surface (Baseman, Banai, and Kahane 1982
; Paul and Lee 1987
; Nishimura et al. 1988
; Brooks et al. 1989
; Tavakkol and Burness 1990
; Saada et al. 1991
; Wybenga et al. 1996
) by acting as a decoy receptor. Surface glycoproteins containing sialic acids, are principal receptors used by pathogens to invade their target tissues (Karlsson 1995
). Given its abundance on the erythrocyte surface and its heavy glycosylation (including sialic acid groups), GYPA could function as a decoy receptor attracting pathogens that bind cell surface sugar groups to the anucleated erythrocyte and away from more vital target tissues. Such a flypaper strategy would explain why the GYPA sequences have diverged so greatly among species, with GYPA in each species adapting specifically to its own pathogens.
This hypothesis also explains the maintenance of GYPA variation within humans by balancing selection, because heterozygotes that express two structural forms of GYPA would be potentially able to target more pathogens. In addition to the primate sequence data given here, there is also general evidence of variation in the glycosylation found on the erythrocyte surface both within humans (Gardner et al. 1989
) and among different mammals (Gagneux and Varki 1999
). This is supported by comparing the predicted O-linked glycosylation sites of primate GYPA sequences, which demonstrate how changes in amino acid sequence elsewhere in the protein can influence the predicted glycosylation of conserved serine or threonine residues through indirect effects such as protein folding (fig. 5
). These differences in amino acid sequence that affect the folding and tertiary structure of surface proteins may also affect pathogen binding. Consistent with this, it has been observed that human erythrocytes with the GYPA allelic blood group M antigen are preferentially bound by certain strains of Escherichia coli, compared with cells with only N (Vaisanen et al. 1982
; Brooks et al. 1989
), where the difference between the two blood groups is defined by amino acid differences and not by glycosylation differences. Functional studies on the interaction between GYPA and nonerythrocytic pathogens and the effect that sequence variation has on glycosylation and pathogen binding will be important for testing the pathogen decoy hypothesis further. The strong non-neutral patterns in the gypa gene identified in this study highlight the erythrocyte surface as an important area for evolutionary studies.
|
|
Supplementary Material |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
New sequences described in this study: Human GYPA ex2–6 (AJ309706) Chimp GYPA ex2–6 (AJ309708). Human GYPA ex2–4 (AJ309828–45, AJ311318–32).
|
Acknowledgements |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
This work was supported by the Wellcome Trust (Prize Studentship for J.B.). We are grateful to Spencer Polley for advice on statistical analyses and Robin Weiss for helpful comments and suggestions.
|
Footnotes |
|---|
Howard Ochman, Reviewing Editor
Abbreviations: GYPA, Glycophorin A; dN, number of nonsynonymous substitutions per nonsynonymous site; dS, number of synonymous substitutions per synonymous site. 
Keywords: natural selection
erythrocyte
glycophorin A
infectious disease
Address for correspondence and reprints: Jake Baum, Department of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, Keppel Street, London WC1E 7HT, United Kingdom. jacob.baum{at}lshtm.ac.uk
|
References |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Aminoff D., 1988 The role of sialoglycoconjugates in the aging and sequestration of red cells from circulation Blood Cells 14:229-257[Web of Science][Medline]
Baseman J. B., M. Banai, I. Kahane, 1982 Sialic acid residues mediate Mycoplasma pneumoniae attachment to human and sheep erythrocytes Infect. Immun 38:389-391
Binks R. H., J. Baum, A. M. J. Oduola, D. E. Arnot, H. A. Babiker, G. Kremsner, C. Roper, B. M. Greenwood, D. J. Conway, 2001 Population genetic analysis of the Plasmodium falciparum erythrocyte binding antigen-175 (eba-175) gene Mol. Biochem. Parasitol 114:63-70[Web of Science][Medline]
Brooks D. E., J. Cavanagh, D. Jayroe, J. Janzen, R. Snoek, T. J. Trust, 1989 Involvement of the MN blood group antigen in shear-enhanced hemagglutination induced by the Escherichia coli F41 adhesin Infect. Immun 57:377-383
Cartron J., J. London, 1992 The protein and gene structure of red cell glycophorins Pp. 101–151 in P. Agre and J. Cartron, eds. Protein blood group antigens of the human red cell: structure, function, and clinical significance. Johns Hopkins University Press, Baltimore
Fu Y. X., W. H. Li, 1993 Statistical tests of neutrality of mutations Genetics 133:693-709[Abstract]
Gagneux P., A. Varki, 1999 Evolutionary considerations in relating oligosaccharide diversity to biological function Glycobiology 9:747-755
Gardner B., S. F. Parsons, A. H. Merry, D. J. Anstee, 1989 Epitopes on sialoglycoprotein alpha: evidence for heterogeneity in the molecule Immunology 68:283-289[Web of Science][Medline]
Goldman N., Z. Yang, 1994 A codon-based model of nucleotide substitution for protein-coding DNA sequences Mol. Biol. Evol 11:725-736[Abstract]
Hansen J. E., O. Lund, N. Tolstrup, A. A. Gooley, K. L. Williams, S. Brunak, 1998 NetOGlyc: prediction of mucin type O-glycosylation sites based on sequence context and surface accessibility Glycoconj. J 15:115-130[Web of Science][Medline]
Huang C. H., O. O. Blumenfeld, 1995 MNSs blood groups and major glycophorins: molecular basis for allelic variation Pp. 153–188 in J. P. Cartron and P. Rouger, eds. Blood cell biochemistry: molecular basis of major human blood group antigens. Plenum Press, New York
Hudson R. R., M. Kreitman, M. Aguade, 1987 A test of neutral molecular evolution based on nucleotide data Genetics 116:153-159
Kaessmann H., V. Wiebe, G. Weiss, S. Paabo, 2001 Great ape DNA sequences reveal a reduced diversity and an expansion in humans Nat. Genet 27:155-156[Web of Science][Medline]
Karlsson K. A., 1995 Microbial recognition of target-cell glycoconjugates Curr. Opin. Struct. Biol 5:622-635[Web of Science][Medline]
Kudo S., M. Fukuda, 1989 Structural organization of glycophorin A and B genes: glycophorin B gene evolved by homologous recombination at Alu repeat sequences Proc. Natl. Acad. Sci. USA 86:4619-4623
Kumar S., K. Tamura, M. Nei, 1994 MEGA: Molecular Evolutionary Genetics Analysis software for microcomputers Comput. Appl. Biosci 10:189-91
Matsui Y., S. Natori, M. Obinata, 1989 Isolation of the cDNA clone for mouse glycophorin, erythroid-specific membrane protein Gene 77:325-332[Web of Science][Medline]
Miller L. H., J. D. Haynes, F. M. McAuliffe, T. Shiroishi, J. R. Durocher, M. H. McGinniss, 1977 Evidence for differences in erythrocyte surface receptors for the malarial parasites, Plasmodium falciparum and Plasmodium knowlesi J. Exp. Med 146:277-281
Mourant A. E., A. C. Kopec, K. Domaniewska-Sobczak, 1978 Blood groups and diseases Oxford University Press, Oxford
Nei M., T. Gojobori, 1986 Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions Mol. Biol. Evol 3:418-426[Abstract]
Nielsen R., Z. Yang, 1998 Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene Genetics 148:929-936
Nishimura H., K. Sugawara, F. Kitame, K. Nakamura, 1988 Attachment of influenza C virus to human erythrocytes J. Gen. Virol 69:2545-2553
Pasvol G., J. S. Wainscoat, D. J. Weatherall, 1982 Erythrocytes deficiency in glycophorin resist invasion by the malarial parasite Plasmodium falciparum Nature 297:64-66[Medline]
Paul R. W., P. W. Lee, 1987 Glycophorin is the reovirus receptor on human erythrocytes Virology 159:94-101[Web of Science][Medline]
Pisano A., J. W. Redmond, K. L. Williams, A. A. Gooley, 1993 Glycosylation sites identified by solid-phase Edman degradation: O-linked glycosylation motifs on human glycophorin A Glycobiology 3:429-435
Pohlmann S., F. Baribaud, B. Lee, G. J. Leslie, M. D. Sanchez, K. Hiebenthal-Millow, J. Munch, F. Kirchhoff, R. W. Doms, 2001a. DC-SIGN interactions with human immunodeficiency virus type 1 and 2 and simian immunodeficiency virus J. Virol 75:4664-4672
Pohlmann S., E. J. Soilleux, F. Baribaud, G. J. Leslie, L. S. Morris, J. Trowsdale, B. Lee, N. Coleman, R. W. Doms, 2001b. DC-SIGNR, a DC-SIGN homologue expressed in endothelial cells, binds to human and simian immunodeficiency viruses and activates infection in trans Proc. Natl. Acad. Sci. USA 98:2670-2675
Rearden A., A. Magnet, S. Kudo, M. Fukuda, 1993 Glycophorin B and glycophorin E genes arose from the glycophorin A ancestral gene via two duplications during primate evolution J. Biol. Chem 268:2260-2267
Rozas J., R. Rozas, 1999 DnaSP version 3: an integrated program for molecular population genetics and molecular evolution analysis Bioinformatics 15:174–175.
Saada A. B., Y. Terespolski, A. Adoni, I. Kahane, 1991 Adherence of Ureaplasma urealyticum to human erythrocytes Infect. Immun 59:467-469
Sim B. K. L., C. E. Chitnis, K. Wasniowska, T. J. Hadley, L. H. Miller, 1994 Receptor and ligand domains for invasion of erythrocytes by Plasmodium falciparum Science 264:1941-1944
Tajima F., 1989 Statistical method for testing the neutral mutation hypothesis by DNA polymorphism Genetics 123:585-595
Tavakkol A., A. T. Burness, 1990 Evidence for a direct role for sialic acid in the attachment of encephalomyocarditis virus to human erythrocytes Biochemistry 29:10684-10690[Medline]
Tomita M., V. T. Marchesi, 1975 Amino-acid sequence and oligosaccharide attachment sites of human erythrocyte glycophorin Proc. Natl. Acad. Sci. USA 72:2964-2968
Vaisanen V., T. K. Korhonen, M. Jokinen, C. G. Gahmberg, C. Ehnholm, 1982 Blood group M specific haemagglutinin in pyelonephritogenic Escherichia coli Lancet 1:1192.
Wybenga L. E., R. F. Epand, S. Nir, J. W. Chu, F. J. Sharom, T. D. Flanagan, R. M. Epand, 1996 Glycophorin as a receptor for Sendai virus Biochemistry 35:9513-9518[Medline]
Yang Z., 1998 Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution Mol. Biol. Evol 15:568-573[Abstract]
Yang Z., R. Nielsen, N. Goldman, A. M. Pedersen, 2000 Codon-substitution models for heterogeneous selection pressure at amino acid sites Genetics 155:431-449
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  M. Fumagalli, R. Cagliani, U. Pozzoli, S. Riva, G. P. Comi, G. Menozzi, N. Bresolin, and M. Sironi Widespread balancing selection and pathogen-driven selection at blood group antigen genes Genome Res., February 1, 2009; 19(2): 199 - 212. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 B. Nystedt, A. C. Frank, M. Thollesson, and S. G. E. Andersson Diversifying Selection and Concerted Evolution of a Type IV Secretion System in Bartonella Mol. Biol. Evol., February 1, 2008; 25(2): 287 - 300. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
A. Varki and T. K. Altheide Comparing the human and chimpanzee genomes: Searching for needles in a haystack Genome Res., December 1, 2005; 15(12): 1746 - 1758. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 C. Layez, P. Nogueira, V. Combes, F. T. M. Costa, I. Juhan-Vague, L. H. P. da Silva, and J. Gysin Plasmodium falciparum rhoptry protein RSP2 triggers destruction of the erythroid lineage Blood, November 15, 2005; 106(10): 3632 - 3638. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 P. K. Shin, P. Pawar, K. Konstantopoulos, and J. M. Ross Characteristics of new Staphylococcus aureus-RBC adhesion mechanism independent of fibrinogen and IgG under hydrodynamic shear conditions Am J Physiol Cell Physiol, September 1, 2005; 289(3): C727 - C734. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 E. J. Vallender and B. T. Lahn Positive selection on the human genome Hum. Mol. Genet., October 1, 2004; 13(suppl_2): R245 - R254. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 T. Angata, E. H. Margulies, E. D. Green, and A. Varki Large-scale sequencing of the CD33-related Siglec gene cluster in five mammalian species reveals rapid evolution by multiple mechanisms PNAS, September 7, 2004; 101(36): 13251 - 13256. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
L. C. Filip and N. I. Mundy Rapid Evolution by Positive Darwinian Selection in the Extracellular Domain of the Abundant Lymphocyte Protein CD45 in Primates Mol. Biol. Evol., August 1, 2004; 21(8): 1504 - 1511. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
G. H. Pogson and K. A. Mesa Positive Darwinian Selection at the Pantophysin (Pan I) Locus in Marine Gadid Fishes Mol. Biol. Evol., January 1, 2004; 21(1): 65 - 75. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 P. Gagneux, M. Cheriyan, N. Hurtado-Ziola, E. C. M. B. van der Linden, D. Anderson, H. McClure, A. Varki, and N. M. Varki Human-specific Regulation of {alpha}2-6-linked Sialic Acids J. Biol. Chem., November 28, 2003; 278(48): 48245 - 48250. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 J. Lu, W.-H. Li, and C.-I Wu Comment on "Chromosomal Speciation and Molecular Divergence-Accelerated Evolution in Rearranged Chromosomes" Science, November 7, 2003; 302(5647): 988b - 988b. [Full Text] [PDF]
|
|
|
|
|
|
N. I. Mundy and S. Cook Positive Selection During the Diversification of Class I Vomeronasal Receptor-like (V1RL) Genes, Putative Pheromone Receptor Genes, in Human and Primate Evolution Mol. Biol. Evol., November 1, 2003; 20(11): 1805 - 1810. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
H.-Y. Wang, H. Tang, C.-K. J. Shen, and C.-I Wu Rapidly Evolving Genes in Human. I. The Glycophorins and Their Possible Role in Evading Malaria Parasites Mol. Biol. Evol., November 1, 2003; 20(11): 1795 - 1804. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J. M. Akey, G. Zhang, K. Zhang, L. Jin, and M. D. Shriver Interrogating a High-Density SNP Map for Signatures of Natural Selection Genome Res., December 1, 2002; 12(12): 1805 - 1814. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||