Nucleotide substitution at major histocompatibility complex class II loci: evidence for overdominant selection

Molecular evolution of snake toxins: is the functional diversity of snake toxins associated with a mechanism of accelerated evolution?

Recent studies revealed that animal toxins with unrelated biological functions often possess a similar architecture. To tentatively understand the evolutionary mechanisms that may govern this principle of functional prodigality associated with a structural economy, two complementary approaches were considered. One of them consisted of investigating the rates of mutations that occur in cDNAs and/or genes that encode a variety of toxins with the same fold. This approach was largely adopted with phospholipases A2 from Viperidae and to a lesser extent with three-fingered toxins from Elapidae and Hydrophiidae. Another approach consisted of investigating how a given fold can accommodate distinct functional topographies. Thus, a number of topologies by which three-fingered toxins exert distinct functions were investigated either by making chemical modifications and/or mutational analyses or by studying the three-dimensional structure of toxin-target complexes. This review shows that, although the two approaches are different, they commonly indicate that most if not all the surface of a snake toxin fold undergoes natural engineering, which may be associated with an accelerated rate of evolution. The biochemical process by which this phenomenon occurs remains unknown.

Analysis of intron sequences at the class II HLA-DRB1 locus: implications for the age of allelic diversity

Analyses of the coding sequences of HLA class II alleles have revealed high similarity between species, indicating that much of the polymorphism predates the separation of human (Homo) and chimpanzee (Pan), 4-7.4 million years ago. Recent studies of the intron sequences of alleles provide support for a much more recent origin and rapid generation of HLA alleles. At the DRB1 locus, intron analysis indicates that most of the allelic lineages have diverged from each other before the separation of Homo and Pan, consistent with the exon analysis. However, the intron sequences of alleles within lineages are almost identical, indicating that many of the alleles have been generated after the divergence of the Homo and Pan lineages.

Effects of intra-locus recombination of HLA polymorphism

The observed number of alleles (n(a)) at some human major histocompatibility complex (HLA) loci is more than twice as large as the pairwise mean number (KB) of nonsynonymous nucleotide substitutions in the peptide-binding region (PBR), although the KB is about the same as that predicted by a model of balancing selection. In order to assess the joint effect of intra-locus recombination and balancing selection on the HLA polymorphism, large scale computer simulation is conducted. It reveals that (1) both rate and location of break points of recombination are crucial factors to determine n(a), (2) intra-locus recombination tends to decrease KB, and (3) the rate of PBR nonsynonymous substitutions is insensitive to the recombination rate. Although a large number of alleles can be generated by intra-locus recombination, this fact alone is insufficient to account for the HLA polymorphism. It is argued that even patchwork patterns of PBR motifs, presumably owing to intra-locus recombination, must be maintained primarily by balancing selection.

No evidence for major histocompatibility complex-dependent mating patterns in a free-living ruminant population

Conventionally, the extraordinary diversity of the vertebrate major histocompatibility complex (MHC is thought to have evolved in response to parasites and pathogens affecting fitness. More recently, reproductive mechanisms such as disassortative mating have been suggested as alternative mechanisms maintaining MHC diversity. A large unmanaged population of Soay sheep (Ovis aries L.) was used to investigate reproductive mechanisms in the maintenance of MHC diversity. Animals were sampled as new-born lambs and between 887 and 1209 individuals were typed at each of five microsatellite markers located either within or flanking the ovine MHC. All loci were in Hardy-Weinberg proportions. A novel likelihood-based approach was developed to analyse mating patterns using paternity data. No evidence for non-random mating with respect to MHC markers was found using this technique. We conclude that MHC diversity in the St Kildan Soay sheep population is unlikely to be maintained by mating preferences and that, in contrast with evidence from experimental mice populations, MHC variation plays no role in the mating structure of this population.

Nucleotide sequence evolution at the kappa-casein locus: evidence for positive selection within the family Bovidae

Kappa-casein is a mammalian milk protein involved in a number of important physiological processes. In the gut, the ingested protein is split into an insoluble peptide (para kappa-casein) and a soluble hydrophilic glycopeptide (caseinomacropeptide). Caseinomacropeptide is responsible for increased efficiency of digestion, prevention of neonate hypersensitivity to ingested proteins, and inhibition of gastric pathogens. Variation within this peptide has significant effects associated with important traits such as milk production. The nucleotide sequences for regions of kappa-casein exon and intron four were determined for representatives of the artiodactyl family Bovidae. The pattern of nucleotide substitution in kappa-casein sequences for distantly related bovid taxa demonstrates that positive selection has accelerated their divergence at the amino acid sequence level. This selection has differentially influenced the molecular evolution of the two kappa-casein split peptides and is focused within a 34-codon region of caseinomacropeptide.

Strong balancing selection at HLA loci: evidence from segregation in South Amerindian families

The genotypic proportions for major histocompatibility complex loci, HLA-A and HLA-B, of progeny in families in 23 South Amerindian tribes in which segregation for homozygotes and heterozygotes could occur are examined. Overall, there is a large deficiency of homozygotes compared with Mendelian expectations (for HLA-A, 114 observed and 155.50 expected and for HLA-B 110 observed and 144.75 expected), consistent with strong balancing selection favoring heterozygotes. There is no evidence that these deficiencies were associated with particular alleles or with the age of the individuals sampled. When these families were divided into four mating types, there was strong selection against homozygotes, averaging 0.462 for three of the mating types over the two loci. For the other mating type in which the female parent is homozygous and shares one allele with the heterozygous male parent, there was no evidence of selection against homozygotes. A theoretical model incorporating these findings surprisingly does not result in a stable polymorphism for two alleles but does result in an excess of heterozygotes and a minimum fitness at intermediate allele frequencies. However, for more than two alleles, balancing selection does occur and the model approaches the qualities of the symmetrical heterozygote advantage model as the number of alleles increases.

The atypical codon usage of the plant psbA gene may be the remnant of an ancestral bias

The psbA gene of the chloroplast genome has a codon usage that is unusual for plant chloroplast genes. In the present study the evolutionary status of this codon usage is tested by reconstructing putative ancestral psbA sequences to determine the pattern of change in codon bias during angiosperm divergence. It is shown that the codon biases of the ancestral genes are much stronger than all extant flowering plant psbA genes. This is related to previous work that demonstrated a significant increase in synonymous substitution in psbA relative to other chloroplast genes. It is suggested, based on the two lines of evidence, that the codon bias of this gene currently is not being maintained by selection. Rather, the atypical codon bias simply may be a remnant of an ancestral codon bias that now is being degraded by the mutation bias of the chloroplast genome, in other words, that the psbA gene is not at equilibrium. A model for the evolution of selective pressure on the codon usage of plant chloroplast genes is discussed.

Molecular evolution of the vertebrate immune system

Adaptive immunity is unique to the vertebrates, and the molecules involved (including immunoglobulins, T cell receptors and the major histocompatibility complex molecules) seem to have diversified very rapidly early in vertebrate history. Reconstruction of gene phylogenies has yielded insights into the evolutionary origin of a number of molecular systems, including the complement system and the major histocompatibility complex (MHC). These analyses have indicated that the C5 component of complement arose by gene duplication prior to the divergence of C3 and C4, which suggests that the alternative complement pathway was the first to evolve. In the case of the MHC, phylogenetic analysis supports the hypothesis that MHC class II molecules evolved before class I molecules. The fact that the MHC-linked proteasome components that specifically produce peptides for presentation by class I MHC appear to have originated before the separation of jawed and jawless vertebrates suggests that the MHC itself may have been present at this time. Immune system gene families have evolved by gene duplication, interlocus recombination and (in some cases) positive Darwinian selection favoring diversity at the amino acid level.

Dissimilar evolution of B-locus versus A-locus and class II loci of the HLA region in South American Indian tribes

Native American populations have a limited HLA polymorphism compared with other ethnic groups. In spite of this, many novel HLA-B locus alleles, not observed in other populations, have been identified in South American tribes, and rapid evolution of this locus has been suggested. We have studied unrelated subjects of the Toba (TOB n = 116), Wichi (WIC n = 46) and Pilaga (PIL n = 14) tribes from northeastern Argentina to investigate the extent of the HLA polymorphism and obtain clues of selective forces that may have acted in these populations. In these tribes the number of HLA alleles is small at all loci except HLA-B, which presents 22 alleles. Seven novel alleles were characterized including 5 of HLA-B (B*35092, B*3518, B*3519, B*4009, B*4803) 1 at HLA-A (A*0219) and 1 at DRB1 (DRB1*0417). All these variants may have arisen by gene conversion events. Some of the novel variants represent the most frequent alleles of these populations (B*4803 in TOB and PIL; B*3519 in WIC) or are the most frequent subtypes in their lineages. HLA-A, B, DRB1,DQA1 and DQB1, but not DPB1, display relatively similar gene frequencies. This results in high heterozygosity in all the tribes for all the loci studied except HLA-DPB1. The larger polymorphism and the generation and maintenance of novel alleles at the HLA-B locus suggests a more specialized response of this locus to evolutionary forces. These effects may be related to the nature of the polymorphism, to the number of founder alleles and to the functional characteristics of the individual alleles.

Evolution and selection of primate T cell antigen receptor BV8 gene subfamily

The set of potential T cell receptor specificities is highly diverse. The relative contributions of T cell receptor (TCR) V beta gene segment polymorphisms, duplications, deletions, and gene conversions to this final T cell receptor protein diversity are unknown. To study these mechanisms, we sequenced and compared closely related primate TCR gene segments from BV8S1, S2, and S5. Interspecies comparisons show that these gene segments have sustained multiple duplication, gene conversion, and deletion events during the last 35 million years of anthropoid primate evolution. BV8 coding sequences are generally conserved with respect to their flanking noncoding sequences, but we find no evidence for positive or negative selection in sequences coding for the first two putative complementarity-determining (ligand-binding) regions. Sequences of TCRBV8 gene segments from unrelated humans demonstrate no nonsynonymous substitutions in nonleader regions of either the BV8S1 or S2 gene segments. We conclude that gene duplication, deletion, and conversion mechanism contribute in a substantial way to the overall diversity of the TCRBV8 gene segment repertoire in primate evolution and that germline substitutions and consequent polymorphisms in CDRs 1 and 2 of these gene segments probably do not play an active role in generating TCR beta chain protein variation.

High-resolution HLA-DRB typing using denaturing gradient gel electrophoresis and direct sequencing

High-resolution HLA-DRB typing is required for bone marrow transplantation between unrelated donors and recipients and also for identification of novel HLA-DRB alleles. Here we describe a method for the unambiguous identification of HLA-DRB alleles using the polymerase chain reaction (PCR), denaturing gradient gel electrophoresis (DGGE) and direct sequencing. The highly variable second exon of all HLA-DRB1, -DRB3, -DRB4, -DRB5, -DRB6 and -DRB7 alleles was amplified using a single pair of generic DRB-specific primers and alleles were separated by DGGE. DNA was then reamplified from plugs removed from the gel and the sequences of these alleles were determined using fluorescent-based sequencing and allele-assignment software. The validity of this typing procedure was confirmed by identification of HLA-DRB alleles for 17 individuals previously characterized by PCR-SSP and/or cloning and sequencing techniques. We identified 34 different HLA-DRB alleles in these 17 unrelated individuals. Importantly, our analysis revealed HLA-DRB1 alleles which had not been identified using the PCR-SSP typing technique. Additionally, alleles from the HLA-DRB3, -DRB4 and -DRB5 loci were identified. Whereas traditional HLA-DRB typing methods provide limited information or require the use of multiple oligonucleotide primers or probes, our technique provides a reliable, specific and relatively rapid way of identifying all HLA-DRB alleles for high-resolution tissue typing.

Evolution of the mitochondrial DNA control region in the mbuna (Cichlidae) species flock of Lake Malawi, East Africa

Considerable controversy has surrounded the application of mitochondrial DNA data to reconstruction of evolutionary relationships among the endemic cichlids of Lake Malawi. Central to this debate has been the issue of whether lineage sorting is complete, and thus whether these data actually reflect species phylogeny, or simply gene genealogy. Review of all mtDNA control region sequences available for members of one monophyletic subset of this species flock, the Malawi rockfishes, or mbuna, strongly indicates that lineage sorting is incomplete: Character-based analyses of these sequences reconstruct gene, not species, interrelationships. Analysis of the pattern of nucleotide substitutions differentiating these mtDNA alleles suggests that pyrimidine residues undergo transition substitutions more often than do purines. Estimation of the magnitude of derived sequence differentiation in light of the reconstructed gene genealogy suggests that the mbuna may be of considerably more recent vintage than previous molecular characterizations have indicated.

Genetic drift can dominate short-term human immunodeficiency virus type 1 nef quasispecies evolution in vivo

The evolution of human immunodeficiency virus (HIV) type 1 nef quasispecies in a patient clonally infected with a contaminated batch of blood clotting factor IX was monitored. nef sequences were derived at 11, 25, and 41 months postinfection from infected peripheral blood mononuclear cells after molecular cloning of PCR-amplified proviral DNA. The phylogenetic relationships among a total of 41 informative sequences were established by split decomposition analysis and used as a basis to establish a substitution matrix and to score synonymous (s) and nonsynonymous (ns) substitutions. The number of observed in-phase stop codons within the nef sequences was comparable to that expected on a random basis. Similarly, the numbers of observed s and ns substitutions did not differ significantly from expected values. No codon position was preferentially mutated. The maximum sequence divergence increased in a linear manner, with approximately 4.4 nucleotide and approximately 3.2 amino acid changes per year. It appears that stochastic processes strongly influence short-term HIV nef quasispecies evolution in vivo.

Ancestral polymorphism of Mhc class II genes in mice: implications for balancing selection and the mammalian molecular clock

To investigate the evolutionary dynamics at Mhc class II DR genes of mice (genus Mus), we sequenced the peptide binding regions (PBRs) of 41 DRB (= E beta) genes and eight DRA (= E alpha) genes from 15 strains representing eight species. As expected trees of these PBR sequences imply extensive maintenance of ancestral DRB alleles across species. We use a coalescent simulation model to show that the number of interspecific coalescent events (c) observed on these trees was higher than the number expected for neutral genealogies and similar sample sizes and is more consistent with balancing selection that with neutrality. Patterns of ancestral polymorphism in mouse DRB alleles were also used to examine the tempo of synonymous substitution in the PBR of mouse class II genes. Both absolute and relative rate tests on DRA and DRB genes imply increased substitution rates at two- and fourfold degenerate sites of mice and rats relative to primates, and decreased rates for the DRB genes of primates relative to ungulate and carnivore relatives. Thus rates of synonymous substitution at Mhc DR genes in mammals appear to be subject to generation time effects in ways similar to those found at other mammalian genes.

Structures of genes encoding phospholipase A2 inhibitors from the serum of Trimeresurus flavoviridis snake

Inhibitors (PLIs) against snake venom gland phospholipases A2 (PLA2s) have been found in their sera. A cDNA encoding a PLI from Trimeresurus flavoviridis (Tf, habu snake, Crotalinae) serum, cPLI-A, was isolated from the Tf liver cDNA library and sequenced. Northern blot analysis with cPLI-A showed that PLIs are expressed only in liver. Genes for PLIs, gPLI-A and gPLI-B, were isolated from the Tf genomic DNA library and their nucleotide (nt) sequences were determined. The genes consisted of four exons and three introns, and exon 4 encoded the carbohydrate recognition domain (CRD)-like motif. Comparison of the nt sequences between gPLI-A and gPLI-B showed that these genes are highly homologous, including introns, except that exon 3 is rich in nonsynonymous nt substitutions which are almost four times as frequent as synonymous nt substitutions. This evolutionary feature of PLI genes is different from that of venom gland PLA2 isozyme genes in which nonsynonymous nt substitutions are spread over the entire mature protein-coding region.

MHC-DRB allelic sequences incorporate distinct intragenic trans-specific segments

The second exon of primate MHC-DRB genes encodes discrete areas of allelic hypervariability (HVR), which are used as the basis for lineage assignments to determine genetic and evolutionary relationships. Comparisons of these regions have led to the "trans-species hypothesis", which proposes that certain MHC alleles from one species are more closely related to those from other species than they are to each other; i.e., that allelic lineages are ancestral in origin. We evaluated this paradigm in an analysis of macaque and baboon MHC-DRB genes using oligotyping and sequencing of 87 new nonhuman primate DRB alleles. A remarkable conservation of sequence motifs in the HVRIII region (codon 60-79) was observed, detected both by hybridization and by sequencing; some of these motifs were found in species such as prosimians that have diverged from the human lineage 50 MYA. However, these fixed HVRIII motif sequences nevertheless occur on a background of diverse lineages suggesting that it is the segmental motif, rather than the allele per se which is trans-specific in origin. Sequences within the first hypervariable region (codons 7-14) identified lineage assignments to several DRB loci (DRB1, DRB3, DRB4, DRB5, DRB6 and DRB7), although a large number of DRB nucleotide sequences did not correspond to a defined allelic motif, suggesting that many of the nonhuman sequences lack human HVRI homologs and have accumulated additional intraspecies variation subsequent to speciation. While there are certain allelic lineages in HVRI that show trans-species conservation, other sequence motifs seem purely species-specific. These differences suggest that HVRI and HVRIII regions have distinct mechanisms for maintenance of trans-specific sequence elements, with different evolutionary histories for segmental nucleotide conservation.

How large was the founding population of Darwin's finches?

A key assumption of many allopatric speciation models is that evolution in peripheral or isolated populations is facilitated by drastic reductions in population size. Population bottlenecks are believed to lead to rapid changes in gene frequencies through genetic drift, to facilitate rapid emergence of novel phenotypes, and to enhance reproductive isolation via genetic revolutions. For such effects to occur, founding populations must be very small, and remain small for some time after founding. This assumption has, however, rarely been tested in nature. One approach is to exploit the polymorphism of the major histocompatibility complex (Mhc ) to obtain information about the founding population. Here, we use the Mhc polymorphism to estimate the size of the founding population of Darwin's finches in the Galápagos Archipelago. The results indicate that the population could not have been smaller than 30 individuals.

Mutation rate varies among alleles at a microsatellite locus: phylogenetic evidence

The understanding of the mutational mechanism that generates high levels of variation at microsatellite loci lags far behind the application of these genetic markers. A phylogenetic approach was developed to study the pattern and rate of mutations at a dinucleotide microsatellite locus tightly linked to HLA-DQB1 (DQCAR). A random Japanese population (n = 129) and a collection of multiethnic samples (n = 941) were typed at the DQB1 and DQCAR loci. The phylogeny of DQB1 alleles was then reconstructed and DQCAR alleles were superimposed onto the phylogeny. This approach allowed us to group DQCAR alleles that share a common ancestor. The results indicated that the DQCAR mutation rate varies drastically among alleles within this single microsatellite locus. Some DQCAR alleles never mutated during a long period of evolutionary time. Sequencing of representative DQCAR alleles showed that these alleles lost their ability to mutate because of nucleotide substitutions that shorten the length of uninterrupted CA repeat arrays; in contrast, all mutating alleles had relatively longer perfect CA repeat sequences.

Human immunodeficiency virus type 1 molecular evolution and the measure of selection

Human immunodeficiency virus (HIV) envelope genes are highly variable between and often within individuals. Part of this variability is thought to be the result of immune-mediated positive selection for sequence diversity. To measure positive selection it has become customary in HIV research to calculate the ratio of the proportions of synonymous (ds) and nonsynonymous (dn) substitutions per potential synonymous or nonsynonymous site, respectively. However, another measure that can be used is the difference between ds and dn, delta d. We show, by example, that using the ratio, ds/dn, or the difference, delta d, may lead us to different conclusions regarding the existence of positive selection pressure. We conclude by noting that until we understand the processes that mediate nucleotide variation in a host selective environment, inferences based on summary statistics characterizing types of nucleotide substitutions should be made with caution.

Functional significance of polymorphism among MHC class II gene promoters

The functional significance of polymorphism among MHC class II promoters in man and mouse is here reviewed, mainly in terms of the hypothesis of differential expression. The hypothesis proposes that differences between antigen-presenting cells in MHC class II expression exert a co-dominant effect on the Th1-Th2 cytokine balance, such that class II molecules of one type come to control to a greater extent the production of one group of cytokines, and those of another type the production of the alternative group. The survey deals with the influence of signal strength and antigen-presenting cell type on T-cell subset differentiation; functional differences between MHC class II molecules not obviously related to determinant selection; disease protection mediated by HLA alleles; mechanisms possibly responsible for allotypic and isotypic bias; overdominance (heterozygous advantage) in selection for expression of class II alleles; MHC class II promoter structure and function; inter-locus and inter-allele variability within human MHC class II gene upstream regulatory regions; a comparison of these polymorphisms in mouse and man; read-out of class II promoter function; and a comparison with expression of MHC class I. We conclude that the evidence that this variation is functionally active (i.e. controls expression) is increasing, but is not yet compelling. The crucial test still to come, we suggest, is whether or not the biological effects attributable to this polymorphism will line up with molecular studies on expression.

Phylogenetic analysis in molecular evolutionary genetics

Recent developments of statistical methods in molecular phylogenetics are reviewed. It is shown that the mathematical foundations of these methods are not well established, but computer simulations and empirical data indicate that currently used methods such as neighbor joining, minimum evolution, likelihood, and parsimony methods produce reasonably good phylogenetic trees when a sufficiently large number of nucleotides or amino acids are used. However, when the rate of evolution varies extensively from branch to branch, many methods may fail to recover the true topology. Solid statistical tests for examining the accuracy of trees obtained by neighbor joining, minimum evolution, and least-squares method are available, but the methods for likelihood and parsimony trees are yet to be refined. Parsimony, likelihood, and distance methods can all be used for inferring amino acid sequences of the proteins of ancestral organisms that have become extinct.

Accelerated evolution of crotalinae snake venom gland serine proteases

Eight cDNAs encoding serine proteases isolated from Trimeresurus flavoviridis (habu snake) and T. gramineus (green habu snake) venom gland cDNA libraries showed that nonsynonymous nucleotide substitutions have accumulated in the mature protein-coding regions to cause amino acid changes. Southern blot analysis of T. flavoviridis genomic DNAs using two proper probes indicated that venom gland serine protease genes form a multigene family in the genome. These observations suggest that venom gland serine proteases have diversified their amino acid sequences in an accelerating manner. Since a similar feature has been previously discovered in crotalinae snake venom gland phospholipase A2 (PLA2) isozyme genes, accelerated evolution appears to be universal in plural isozyme families of crotalinae snake venom gland.

Compensatory aspects of allele diversity at immunoglobulin loci: gene correlations in rabbit populations devoid of light chain diversity (Oryctolagus cuniculus L.; Kerguelen Islands)

Is there a selective advantage of increased diversity at one immunoglobulin locus when diversity at another locus is low? A previous paper demonstrated excess heterozygosity at the rabbit light chain b locus when heterozygosity was low at the heavy chain constant region e locus. Here we consider the reverse situation by analyzing allele distributions at heavy chain loci in populations fixed for the light chain b locus. We analyzed the a locus that encodes the predominantly expressed heavy chain variable region, and the d and e loci that control different parts of the Ig gamma class constant region. While there was excess heterozygosity, genetic differentiation between localities was extensive and was most pronounced for females. This was in marked contrast with observations in areas where b-locus diversity was important and confirms a negative correlation between e- and b-locus heterozygosity. Trigenic disequilibria corresponded to a significant negative correlation between e- and a-locus heterozygosity due mainly to strong variation among localities within the context of pronounced (digenic) linkage disequilibria. Although substantial, the average increase in a/e-locus single heterozygosity implemented by higher order disequilibria within localities was not significant.

Molecular analysis of DLA-DRBB1 polymorphism

The polymorphism of the canine major histocompatibility complex class II DRB gene, DRBB1, was analyzed. Exon 2 that encodes the polymorphic beta 1 domain was amplified by the polymerase chain reaction (PCR). The PCR product from 250 dogs was cloned and sequenced. Eighteen alleles were identified. Most of the variation in amino acid composition occurred at positions in the peptide binding site. Inheritance of these sequences showed Mendelian segregation with one or two alleles per dog. Cluster analysis of the nucleotide and predicted amino acid sequences subdivided the canine DRBB1 alleles into three major allelic groups. The number of nonsynonymous changes was higher than the number of synonymous changes in the putative antigen recognition sites indicative of positive selection. The data generated can serve as a basis for developing a typing assay for the canine DRBB1 gene.

Peptide binding function and the paradox of HLA disease associations

Many studies have found associations between HLA loci and susceptibility or resistance to both infectious and autoimmune diseases, but often subsequent studies fail to find the same association. Such inconsistencies are not surprising if we consider what is known of the population biology and evolution of HLA genes, especially the consequences of natural selection favouring heterozygosity in the peptide binding regions (PBR) of HLA molecules. Because of past recombination event and/or convergent evolution, alleles that are not closely related in overall sequence may come to resemble each other in the PBR. Because peptide binding is likely to be important for the role of HLA molecules in both infectious and autoimmune disease, a strategy of searching for HLA disease associations that groups alleles in functional categories based on PBR motifs may prove more successful than conventional strategies. Likewise, in developing approaches for molecular typing, it may be advisable to develop methods that group alleles on the basis of shared PBR motifs rather than simply on the basis of shared primer sites in less functionally important regions.

Allelic diversity and gene genealogy at the self-incompatibility locus in the Solanaceae

The self-incompatibility (S) locus of flowering plants offers an example of extreme polymorphism maintained by balancing selection. Estimates of recent and long-term effective population size (Ne) were determined for two solanaceous species by examination of S-allele diversity. Estimates of recent Ne in two solanaceous species differed by an order of magnitude, consistent with differences in the species' ecology. In one species, the evidence was consistent with historical population restriction despite a large recent Ne. In the other, no severe bottleneck was indicated over millions of years. Bottlenecks are integral to founder-event speciation, and loci that are subject to balancing selection can be used to evaluate the frequency of this mode of speciation.

Variation and evolution of class I Mhc in sexual and parthenogenetic geckos

We present the first Mhc class I sequences in geckos. We compared Mhc variation in gekkonid species that reproduce sexually (Hemidactylus frenatus, Lepidodactylus aureolineatus, L. moestus, L. sp. Arno, L. sp. Takapoto) to others reproducing parthenogenetically (H. garnotii, L. lugubris). These comparisons include the known maternal (L. moestus) and paternal (L. sp. Arno) ancestors of the asexual L. lugubris. Sequences similar to other vertebrate species were obtained from both nuclear and cDNA templates indicating that these sequences are derived from expressed class I Mhc loci. Southern blot analysis using gecko class I probes, revealed that parthenogenetic clonal lineages of independent evolutionary origin have no within-clone band variation at class I loci and that no detectable recombination between restriction sites had taken place. Variability in the sexual species was similar to mammalian taxa, i.e. class I genes are highly variable in outbreeding sexual populations. Sequence analysis of the alpha-2 domain of class I genes identified point mutations in a clonal lineage of L. lugubris which led to amino acid substitutions. Potential transspecific allelic lineages were also observed. The persistence of asexual lineages with little or no class I diversification over thousands of generations seems to argue against strong selection for Mhc multi-allelism caused by pathogen-Mhc allele specificity. On the other hand, the high level of heterozygosity in the parthenogenetic species (a consequence of their hybrid origin) may provide clonal lineages with adequate antigen presenting diversity to survive and compete with sexual relatives.

Stem-loop potential in MHC genes: a new way of evaluating positive Darwinian selection?

The domains of polymorphic major histocompatibility complex (MHC) proteins which interact with peptides and T-cell receptors are considered to have been under positive evolutionary selection pressure. Evidence for this is a high ratio of non-synonymous to synonymous mutations in the corresponding genomic domains. By this criterion snake venom phospholipase A2 genes have also been under positive selection pressure. Recent studies of the latter genes indicate that positive selection has overridden an evolutionary pressure on base order which normally promotes the potential to extrude single-strand stem-loops from supercoiled duplex DNA ( fold pressure ). This has resulted in base order-dependent stem-loop potential being shifted to introns, which are highly conserved between species. It is now shown that, like snake venom phospholipase A2 genes, the domains of polymorphic MHC genes which appear to have responded to positive selection pressure have decreased base order-dependent stem-loop potential. The evolutionary pressure to generate stem-loop potential (believed to be important for recombination) has been overridden less in exons under negative purifying selection than in exons under positive Darwinian selection. Thus, base order-dependent stem-loop potential shows promise as an independent indicator of positive selection.

Ancestral major histocompatibility complex DRB genes beget conserved patterns of localized polymorphisms

Genes within the major histocompatibility complex (MHC) are characterized by extensive polymorphism within species and also by a remarkable conservation of contemporary human allelic sequences in evolutionarily distant primates. Mechanisms proposed to account for strict nucleotide conservation in the context of highly variable genes include the suggestion that intergenic exchange generates repeated sets of MHC DRB polymorphisms [Gyllensten, U. B., Sundvall, M. & Erlich, H. A. (1991) Proc. Natl. Acad. Sci. USA 88, 3686-3690; Lundberg, A. S. & McDevitt, H. 0. (1992) Proc. Natl. Acad. Sci. USA 89, 6545-6549]. We analyzed over 50 primate MHC DRB sequences, and identified nucleotide elements within macaque and baboon DRB6-like sequences with deletions corresponding to specific exon 2 hypervariable regions, which encode a discrete alpha helical segment of the MHC antigen combining site. This precisely localized deletion provides direct evidence implicating segmental exchange of MHC-encoded DRB gene fragments as one of the evolutionary mechanisms both generating and maintaining MHC diversity. Intergenic exchange at this site may be fundamental to the diversification of immune protection in populations by permitting alteration in the specificity of the MHC that determines the repertoire of antigens bound.

The genetic structure of a province as revealed by surnames and HLA genes: potential utility in transplantation policy

The combined use of surnames, which simulate selectively neutral genetic markers, and HLA genes may constitute a useful tool for the genetic survey of a small area. We found a coincidence between HLA genetic structure and surname "selectively neutral' pattern, in an Italian province, although HLA genes indicate a more pronounced genetic isolation for one particular subregion of the province. The patchiness of HLA allele distribution that is sometimes created by drift and/or selection can be used as the basis both of anthropological and of epidemiological studies. The creation of genetic maps of relatively small areas may also constitute a tool for the selection of bone marrow donors.

Positive selection and rates of evolution in immunodeficiency viruses from humans and chimpanzees

Evolutionary theory predicts the recent spread of primate immunodeficiency viruses (PIVs) to new human populations to be accompanied by positive selection in response to new host environments and/or by random genetic drift. I assess evidence for positive selection in human and chimpanzee PIVs type I (PIV1s), using ratios of synonymous to nonsynonymous nucleotide change based on branch lengths and outgroup rooting. Ratios are smaller for PIV1s from humans than for PIV1 from a chimpanzee for the pol, gag, and env glycoprotein 120 (gp120) regions, indicating greater effects of positive selection in PIV1s from humans. Parsimony-based relative rate tests for amino acid changes showed significant differences between PIV1s from humans and chimpanzees in 18 of 48 pairwise comparisons, with all 18 showing faster rates of change in PIV1s from humans. This study indicates that in some instances, the recent evolution of human PIV1s follows a speciational pattern, in which increased diversification of taxa is correlated with greater amounts of character change appearing and being maintained through time. This extends the generality of the speciational pattern to a group of organisms (viruses) having the fastest known rates of anagenetic change for nucleotide characters and indicates that comprehensive understanding of PIV1 evolution requires consideration of both anagenetic change within viral lineages and the relative historical success of different viral clades. Phylogenetic analyses show that neither PIV1s infecting humans nor those infecting chimpanzees represent monophyletic groups and suggest multiple host-species shifts for PIV1s.

Population biology of antigen presentation by MHC class I molecules

In principle, the function of major histocompatibility complex (MHC) molecules is simple: to bind a peptide and engage a T cell. In practice, placing this function within the context of the immune response begs questions of population biology; How does the immune response emerge from the interactions among populations of peptides, T cells and MHC molecules? Within a population of vertebrates, how does MHC polymorphism stamp individuality on the response? Does polymorphism confer differential advantages in responding to parasites? How are the pressures on the MHC reflected in turnover of alleles? The role of mutation, recombination, selection, and drift in the generation and maintenance of MHC class 1 polymorphism are considered.

Evolution in a chronic RNA virus infection: selection on HTLV-I tax protein differs between healthy carriers and patients with tropical spastic paraparesis

HTLV-I causes T-cell leukemia and tropical spastic paraparesis (TSP) in a minority of infected people, whereas the majority remain healthy. The virus differs little in sequence between isolates but has been shown to have a quasispecies structure. Using the Nei and Gojobori algorithm, we have shown that the proportion of nonsynonymous to synonymous changes in HTLV-I proviral tax gene sequences from healthy seropositive subjects (Dn/Ds = 0.9 to 1.3) is significantly higher than those from TSP patients (Dn/Ds = 0.3 to 0.6). Here we show that the distinction between healthy seropositives and TSP patients can only be seen with proviral tax sequences, but not with cDNA, the amino-terminal or carboxy-terminal half of tax, or the rex gene. The Dn/Ds ratio of proviral tax sequences was used to analyze two TSP patients with atypical features and to investigate the influence of cytotoxic T cells (CTL) on the viral quasispecies.

Rapid allelic diversification and intensified selection at antigen recognition sites of the Mhc class II DPB1 locus during hominoid evolution

The evolution of polymorphism at the Mhc class II DPB1 locus was studied by comparison of chimpanzee (Pan troglodytes), pygmy chimpanzee (Pan paniscus), gorilla (Gorilla gorilla) and human DPB1 alleles. Extensive polymorphism was found in all hominoids. The clustering of sequences in the phylogenetic tree is consistent with rapid generation of the DPB1 polymorphism. Analysis of the substitution pattern for human alleles shows an excess of non-synonymous changes to synonymous changes at antigen recognition sites, indicating that the amino acid polymorphism at these sites is being maintained by selection. By contrast, no excess of nonsynonymous changes was found at the antigen recognition sites of nonhuman hominoid species. Thus, it appears that diversifying selection on the DPB1 polymorphism has intensified in the lineage leading to humans. No evidence was found for the existence of ancient allelic lineages predating the divergence of the hominoid species. The number of synonymous differences among DPB1 alleles is lower than among DQB1 and DRB1 alleles, indicative of a more recent origin for the DPB1 polymorphism and consistent with the more rapid evolution suggested by the phylogenetic tree.

Trans-species polymorphism of class II Mhc loci in danio fishes

A characteristic feature of the major histocompatibility complex (Mhc) polymorphism in mammals is the existence of allelic lineages shared by related species. This trans-species polymorphism has thus far been documented only in primates, rodents, and artiodactyls. In this communication we provide evidence that it also exists in cyprinid (bony) fishes at the class II A and B loci coding for the alpha and beta polypeptide chains of the class II alpha:beta heterodimers. The study has focused on three species of the family Cyprinidae, subfamily Rasborinae: the zebrafish (Danio rerio), the giant danio (D. malabaricus), and the pearl danio (D. albolineatus). The polymerase chain reaction was used to amplify and then sequence intron 1 and exon 2 of the class II B loci and exon 2 of the class II A loci in these species. Phylogenetic analysis of the sequences revealed the existence of allelic lineages whose divergence predates the divergence of the three species at both the A and B loci. The lineages at the B locus in particular are separated by large genetic distances. The polymorphism is concentrated in the peptide-binding region sites and is apparently maintained by balancing selection. Sharing of this unique Mhc feature by both bony fishes and mammals suggests that the main function of the Mhc (presentation of peptides to T lymphocytes) has not changed during the last 400 million years of its evolution.

Characterization of Mhc genes in a multigenerational family of ring-necked pheasants

Little is known about the major histocompatibility (Mhc) genes of birds in different taxonomic groups or about how Mhc genes may be organized in avian species divergent by evolution or habitat. Yet it seems likely that much might be learned from birds about the evolution, organization, and function of this intricate complex of polymorphic genes. In this study a close relative of the chicken, the ring-necked pheasant (Phasianus colchicus), was examined for the presence and organization of Mhc B-G genes. The patterns of restriction fragments revealed by chicken B-G probes in Southern hybridizations and the patterns of pheasant erythrocyte polypeptides revealed in immunoblots by antisera raised against chicken B-G polypeptides provide genetic, molecular, and biochemical data confirming earlier serological evidence for the presence of B-G genes in the pheasant, and hence, the presence of a family of B-G genes in at least a second species of birds. The high polymorphism exhibited by the pheasant B-G gene family allowed genetic differences among individuals within the small experimental population in this study to be detected easily by restriction fragment patterns. Further evidence was found for the organization of the pheasant Mhc class I and class II genes into genetically independent clusters. Whether these gene clusters are fully comparable to the B and Rfp-Y systems in the chicken or whether yet another organization of Mhc genes has been encountered in the pheasant remains to be determined.

Dynamics of Mhc evolution in birds and crocodilians: amplification of class II genes with degenerate primers

Genes of the major histocompatibility complex (Mhc) are the most polymorphic functional loci in mammalian populations, but little is known of Mhc variability in natural populations of nonmammalian vertebrates. To help extend such studies to birds and relatives, we present a pair of degenerate primers that amplify polymorphic segments of one chain (the beta chain) of the class II genes from the major histocompatibility complex (Mhc) of archosaurs (birds+crocodilians). The primers target two conserved regions lying within portions of the antigen-binding site (ABS) encoded by the second exon and amplify multiple genes from both genomic DNA and cDNA. The pattern of nucleotide substitution in ABS codons of 51 sequences amplified and cloned from five species of passerine birds and an alligator (Alligator mississippiensis) indicates that archosaurian class II beta genes are subject to selective forces similar to those operating in mammalian populations. Hybridization of a genomic clone generated by the primers revealed highly polymorphic bands in a sample of Florida scrub jays (Aphelocoma coerulescens coerulescens). Because the primers amplify only part of the ABS from multiple class II genes, they will be useful primarily for generating species specific clones, thereby providing a critical inroad to more detailed structural and evolutionary studies.

The evolution of the type I interferon gene family in mammals

A phylogenetic analysis of mammalian type I interferon (IFN) genes showed: (1) that the three main subfamilies of these genes in mammals (IFN-beta, IFN-alpha, and IFN-omega) diverged after the divergence of birds and mammals but before radiation of the eutherian orders and (2) that IFN-beta diverged first. Although apparent cases of interlocus recombination among mouse IFN-alpha genes were identified, the hypothesis that coding regions of IFN-alpha genes have been homogenized within species by interlocus recombination was not supported. Flanking regions as well as coding regions of IFN-alpha were more similar within human and mouse than between these species; and reconstruction of the pattern of nucleotide substitution in IFN-alpha coding regions of four mammalian species by the maximum parsimony method suggested that parallel substitutions have occurred far more frequently between species than within species. Therefore, it seems likely that IFN-alpha genes have duplicated independently within different eutherian orders. In general, type I IFN genes are subject to purifying selection, which in the case of IFN-alpha and IFN-beta is strongest in the putative receptor-binding domains. However, analysis of the pattern of nucleotide substitution among IFN-omega genes suggested that positive Darwinian selection may have acted in some cases to diversify members of this subfamily at the amino acid level.

An evolutionary perspective on signaling in behavior and immunology

Signaling and communication are important at different levels of biological organization. Signals exchanged between cells of the immune system initiate and coordinate the immune response; signals exchanged between individuals often coordinate social behavior. Behavioral ecologists interested in the evolution and functional design of signals exchanged between individuals have produced a theory of signaling and communication that stresses the importance of cooperation and conflict; if a conflict of interest between signaler and receiver is great enough, signals evolve towards greater reliability. We suggest that the application of signaling system theory to signals exchanged between cells within an individual will allow for a better understanding of immunity and intra-individual communication in general, with potential for novel approaches to the treatment of disease.