Functional implications of Major Histocompatibility (MH) variation using estuarine fish populations

Recently, there has been a dramatic expansion of studies ofmajor histocompatibility complex (MHC) variation aimed at discoveringfunctional differences in immunity across wild populations ofdiverse vertebrate species. Some species with relatively lowgenetic diversity or under strong directional selection by pathogenshave revealed fascinating cases of MHC allelic disease linkage.More generally in genetically diverse species, however, theselinkages may be hard to find. In this paper, we review approachesfor assessing functional variation in MHC and discuss theirpotential use for discovering smaller-scale intraspecific spatialand temporal patterns of MHC variation. Then, we describe andillustrate an approach using the structural model to producea population composite of variation in antigen-binding regionsby mapping population-specific substitutions onto functionalregions of the molecule. We are producing models of variationin major histocompatibility (MH) loci for populations of non-migratoryfish (killifish, Fundulus heteroclitus) resident at sites thatvary dramatically in environmental quality. We discuss the goalof relating MH population variation to functional differencesin disease susceptibility such as those inferred by observationsof parasitic infection and direct measurement of bacterial challengesin the laboratory. Our study has focused on relatively well-studiedkillifish populations, including those resident in a highlydisturbed, chemically contaminated estuary and nearby less contaminatedsites. Population-specific genetic changes at MHC antigen-bindingloci are described, and evidence relevant to functional implicationsof these changes is reviewed. Population-specific patterns ofvariation in antigen-binding regions in combination with a rangeof assessments of immune function will provide a powerful newapproach to reveal functional changes in MHC.     

Overdispersed molecular clock at the major histocompatibility complex loci.

The extent of amino acid differences of major histocompatibility complex molecules within species is unusually high, consistent with the finding that some pairs of alleles have persisted for more than ten million years and the view that the polymorphism has been maintained by natural selection. The disparity between synonymous and non-synonymous substitutions in the antigen recognition site, however, suggests that some non-synonymous sites have undergone a number of substitutions whereas others have little or none. To describe statistically such an overdispersed underlying process, commonly used Poisson processes are inadequate. An alternative process leads to the surprising conclusion that each non-synonymous site has accumulated as many as 2.6 substitutions, on the average, in the two lineages leading to humans and mice. The standard deviation is also very large (6.6) and the dispersion index (the ratio of the variance to the mean) is at least 17. The substitution process thus inferred qualitatively agrees with the disposition (a boomerang pattern) of substitutions between HLA-A2 and Aw68 alleles, and quantitatively agrees well with that expected where the evolution of major histocompatibility complex molecules has long been driven mostly by balancing selection.     

A method for estimating the intensity of overdominant selection from the distribution of allele frequencies

A method is proposed for estimating the intensity of overdominant selection scaled by the effective population size, S = 2Ns, from allele frequencies. The method is based on the assumption that, with strong overdominant selection, allele frequencies are nearly at their deterministic equilibrium values and that, to a first approximation, deviations depend only on S. Simulations verify that reasonably accurate estimates of S can be obtained for realistic sample sizes. The method is applied to data from several loci in the major histocompatibility complex (Mhc) in numerous human populations. For alleles distinguished by both serological typing and the sequence of the peptide-binding region, our estimates of S are comparable to those obtained by analysis of DNA sequences in showing that selection is strongest on HLA-B and weaker on HLA-A, HLA-DRB1, and HLA-DQA1. The intensity of selection on HLA-B varied considerably among populations. Two populations, Native American and Inuit, showed an excess rather than a deficiency in homozygosity. Comparable estimates of S were obtained for alleles at Mhc class II loci distinguished by serological reactions (serotyping) and by differences in the amino acid sequences of the peptide-binding region (molecular typing). A comparison of two types of data for DQA1 and DRB1 showed that serotyping led to generally lower estimates of S.     

Polymorphism and Balancing Selection at Major Histocompatibility Complex Loci

Amino acid replacements in the peptide-binding region (PBR) of the functional major histocompatibility complex (Mhc) genes appear to be driven by balancing selection. Of the various types of balancing selection, we have examined a model equivalent to overdominance that confers heterozygote advantage. As discussed by A. Robertson, overdominance selection tends to maintain alleles that have more or less the same degree of heterozygote advantage. Because of this symmetry, the model makes various testable predictions about the genealogical relationships among different alleles and provides ways of analyzing DNA sequences of Mhc alleles. In this paper, we analyze DNA sequences of 85 alleles at the HLA-A, -B, -C, -DRB1 and -DQB1 loci with respect to the number of alleles and extent of nucleotide differences at the PBR, as well as at the synonymous (presumably neutral) sites. Theory suggests that the number of alleles that differ at the sites targeted by selection (presumably the nonsynonymous sites in the PBR) should be equal to the mean number of nucleotide substitutions among pairs of alleles. We also demonstrate that the nucleotide substitution rate at the targeted sites relative to that of neutral sites may be much larger than 1. The predictions of the presented model are in surprisingly good agreement with the actual data and thus provide means for inferring certain population parameters. For overdominance selection in a finite population at equilibrium, the product of selection intensity (s) against homozygotes and the effective population size (N) is estimated to be 350-3000, being largest at the B locus and smallest at the C locus. We argue that N is of the order of 10(5) and s is several percent at most, if the mutation rate per site per generation is 10(-8).     

Persistence of Mhc heterozygosity in homozygous clonal killifish,Rivulus marmoratus: implications for the origin of hermaphroditism

The mangrove killifish Rivulus marmoratus, a neotropical fish in the order Cyprinodontiformes, is the only known obligatorily selfing, synchronous hermaphroditic vertebrate. To shed light on its population structure and the origin of hermaphroditism, major histocompatibility complex (Mhc) class I genes of the killifish from seven different localities in Florida, Belize, and the Bahamas were cloned and sequenced. Thirteen loci and their alleles were identified and classified into eight groups. The loci apparently arose approximately 20 million years ago (MYA) by gene duplications from a single common progenitor in the ancestors of R. marmoratus and its closest relatives. Distinct loci were found to be restricted to different populations and different individuals in the same population. Up to 44% of the fish were heterozygotes at Mhc loci, as compared to near homozygosity at non-Mhc loci. Large genetic distances between some of the Mhc alleles revealed the presence of ancestral allelic lineages. Computer simulation designed to explain these findings indicated that selfing is incomplete in R. marmoratus populations, that Mhc allelic lineages must have diverged before the onset of selfing, and that the hermaphroditism arose in a population containing multiple ancestral Mhc lineages. A model is proposed in which hermaphroditism arose stage-wise by mutations, each of which spread through the entire population and was fixed independently in the emerging clones.     

Effects of habitat fragmentation and changes of dispersal behaviour after a recent population decline on the genetic variability of noncoding and coding DNA of a monogamous Malagasy rodent

While interactions among demography, behaviour and genetic structure are well-documented for neutral genetic markers, the role of these parameters and the effects of genetic drift and selection are considerably less well understood in functional genes, such as the major histocompatibility complex (MHC). In this study, the consequences of habitat fragmentation and the effects of a current population decline on noncoding (mitochondrial DNA) and two coding MHC loci (DQA, DRB) with different functional importance were investigated in the small remnant subdivided population of the endangered Malagasy giant jumping rat (Hypogeomys antimena). Both neutral and selective markers revealed a significant genetic differentiation between the two remnant subpopulations. The FST values were much lower in the MHC DQA and DRB genes than in the mitochondrial data. The MHC DRB loci display the effects of both balancing selection (high sequence diversity, four times higher nonsynonymous than synonymous substitutions in the functionally important antigen-binding site positions, twice the average heterozygosity of individual amino acids at the positions identified as part of the antigen-binding site (ABS) than those outside the ABS and nonselective forces including genetic drift. Simultaneously with a current population decline offspring reduced their dispersal distances. No substantial effects were detected within the first 6 years of reduced gene flow in either mitochondrial or MHC markers.     

Molecular evolution of animal antimicrobial peptides: widespread moderate positive selection

An increasing number of studies in both vertebrates and invertebrates show that the evolution of antimicrobial peptides is driven by positive selection. Because these diverse molecules show potential for therapeutic applications, they are currently the targets of much structural and functional research, providing extensive background data for evolutionary studies. In this paper, patterns of molecular evolution in antimicrobial peptide genes are reviewed. Evidence for positive selection on antimicrobial peptides includes an excess of nonsynonymous nucleotide substitutions, an excess of charge-changing amino acid substitutions, nonneutral patterns of allelic variation, and functional assays in vivo and in vitro that show improved antimicrobial effects for derived sequence variants. Positive selection on antimicrobial peptides may be as common as, but perhaps weaker than, selection on the best-known example of adaptively evolving immunity genes, the major histocompatibility complex. Thus, antimicrobial peptides present a useful and underutilized model for the study of adaptive molecular evolution.     

Lineage-Specific Differences in the Amino Acid Substitution Process

In Darwinian evolution, mutations occur approximately at random in a gene, turned into amino acid mutations by the genetic code. Some mutations are fixed to become substitutions and some are eliminated from the population. Partitioning pairs of closely related species with complete genome sequences by average population size of each pair, we looked at the substitution matrices generated for these partitions and compared the substitution patterns between species. We estimated a population genetic model that relates the relative fixation probabilities of different types of mutations to the selective pressure and population size. Parameterizations of the average and distribution of selective pressures for different amino acid substitution types in different population size comparisons were generated with a Bayesian framework. We found that partitions in population size as well as in substitution type are required to explain the substitution data. Selection coefficients were found to decrease with increasingly radical amino acid substitution and with increasing effective population size.To further explore the role of underlying processes in amino acid substitution, we analyzed embryophyte (plant) gene families from TAED (The Adaptive Evolution Database), where solved structures for at least one member exist in the Protein Data Bank. Using PAML, we assigned branches to three categories: strong negative selection, moderate negative selection/neutrality, and positive diversifying selection. Focusing on the first and third categories, we identified sites changing along gene family lineages and observed the spatial patterns of substitution. Selective sweeps were expected to create primary sequence clustering under positive diversifying selection. Co-evolution through direct physical interaction was expected to cause tertiary structural clustering. Under both positive and negative selection, the substitution patterns were found to be nonrandom. Under positive diversifying selection, significant independent signals were found for primary and tertiary sequence clustering, suggesting roles for both selective sweeps and direct physical interaction. Under strong negative selection, the signals were not found to be independent. All together, a complex interplay of population genetic and protein thermodynamics forces is suggested.     

Genetic drift vs. natural selection in a long-term small isolated population: major histocompatibility complex class II variation in the Gulf of California endemic porpoise (Phocoena sinus)

Although many studies confirm long-term small isolated populations (e.g. island endemics) commonly sustain low neutral genetic variation as a result of genetic drift, it is less clear how selection on adaptive or detrimental genes interplay with random forces. We investigated sequence variation at two major histocompatibility complex (Mhc) class II loci on a porpoise endemic to the upper Gulf of California, México (Phocoena sinus, or vaquita). Its unique declining population is estimated around 500 individuals. Single-strand conformation polymorphism analysis revealed one putative functional allele fixed at the locus DQB (n = 25). At the DRB locus, we found two presumed functional alleles (n = 29), differing by a single nonsynonymous nucleotide substitution that could increase the stability at the dimer interface of alphabeta-heterodimers on heterozygous individuals. Identical trans-specific DQB1 and DRB1 alleles were identified between P. sinus and its closest relative, the Burmeister's porpoise (Phocoena spinipinnis). Comparison with studies on four island endemic mammals suggests fixation of one allele, due to genetic drift, commonly occurs at the DQA or DQB loci (effectively neutral). Similarly, deleterious alleles of small effect are also effectively neutral and can become fixed; a high frequency of anatomical malformations on vaquita gave empirical support to this prediction. In contrast, retention of low but functional polymorphism at the DRB locus was consistent with higher selection intensity. These observations indicated natural selection could maintain (and likely also purge) some crucial alleles even in the face of strong and prolonged genetic drift and inbreeding, suggesting long-term small populations should display low inbreeding depression. Low levels of Mhc variation warn about a high susceptibility to novel pathogens and diseases in vaquita.     

Phylogeography, population dynamics, and molecular evolution of European bat lyssaviruses

European bat lyssaviruses types 1 and 2 (EBLV-1 and EBLV-2) are widespread in Europe, although little is known of their evolutionary history. We undertook a comprehensive sequence analysis to infer the selection pressures, rates of nucleotide substitution, age of genetic diversity, geographical origin, and population growth rates of EBLV-1. Our study encompassed data from 12 countries collected over a time span of 35 years and focused on the glycoprotein (G) and nucleoprotein (N) genes. We show that although the two subtypes of EBLV-1--EBLV-1a and EBLV-1b--have both grown at a low exponential rate since their introduction into Europe, they have differing population structures and dispersal patterns. Furthermore, there were strong constraints against amino acid change in both EBLV-1 and EBLV-2, as reflected in a low ratio of nonsynonymous to synonymous substitutions per site, particularly in EBLV-1b. Our inferred rate of nucleotide substitution in EBLV-1, approximately 5 x 10(-5) substitutions per site per year, was also one of the lowest recorded for RNA viruses and implied that the current genetic diversity in the virus arose 500 to 750 years ago. We propose that the slow evolution of EBLVs reflects their distinctive epidemiology in bats, where they occupy a relatively stable fitness peak.     

Variation on islands: major histocompatibility complex (Mhc) polymorphism in populations of the Australian bush rat

Loss of genetic variation in small, isolated populations is commonly observed at neutral or nearly neutral loci. In this study, the loss of genetic variation was assessed in island populations for a locus of major histocompatibility complex (Mhc), a locus shown to be under the influence of balancing selection. A total of 36 alleles was found at the second exon of RT1.Ba in 14 island and two mainland populations of Rattus fuscipes greyii. Despite this high overall diversity, a substantial lack of variation was observed in the small island populations, with 13 islands supporting only one to two alleles. Two populations, Waldegrave and Williams Islands, showed moderately high levels of heterozygosity (52-56%) which were greater than expected under neutrality, suggesting the action of balancing selection. However, congruence between the level of variation at this Mhc locus and in previous allozyme electrophoresis and mitochondrial DNA studies highlights the dominant influence of genetic drift and population factors, such as bottlenecks and structuring in the founding population, in the loss of genetic variation in these small, isolated populations.     

Sequence diversity of Mhc genes in lake whitefish

The sequence variation of three exons of the major histocompatibility complex ( Mhc ) was examined in a lake whitefish Coregonus sp., population from the Swiss lake of Hallwil. DNA sequences from the Mhc class I A1 , A2 and class II B1 exons, corresponding to the α1, α2 and β1 domains of the Mhc glycoproteins, were obtained by the polymerase chain reaction followed by cloning and sequencing. The numbers of variable sequences detected for each exon were 15 ( A1 ), 11 ( A2 ) and 20 ( B1 ). Levels of nucleotide similarity ranged from 82 to 99% for the A1 exon, 58–96% for the A2 and 88–99% for the B1 exon. At the A1 and B1 exons, the nonsynonymous substitution rates ( dn ) exceeded synonymous substitution rates ( ds ) greatly within the peptide binding regions, indicating the effect of balancing selection. Sequence diversity at the A2 exon did not seem to be maintained by balancing selection ( ds > dn ). Phylogenetic comparison of whitefish Mhc sequences with sequences from other salmonid species and more distantly related teleosts indicated shared ancestral (trans-species) polymorphism.     

Evolution of MHC class IIB in the genome of wild and ornamental guppies, Poecilia reticulata

This is the first study to quantify genomic sequence variation of the major histocompatibility complex (MHC) in wild and ornamental guppies, Poecilia reticulata. We sequenced 196-219 bp of exon 2 MHC class IIB (DAB) in 56 wild Trinidadian guppies and 14 ornamental strain guppies. Each of two natural populations possessed high allelic richness (15-16 alleles), whereas only three or fewer DAB alleles were amplified from ornamental guppies. The disparity in allelic richness between wild and ornamental fish cannot be fully explained by fixation of alleles by inbreeding, nor by the presence of non-amplified sequences (ie null alleles). Rather, we suggest that the same allele is fixed at duplicated MHC DAB loci owing to gene conversion. Alternatively, the number of loci in the ornamental strains has contracted during >100 generations in captivity, a hypothesis consistent with the accordion model of MHC evolution. We furthermore analysed the substitution patterns by making pairwise comparisons of sequence variation at the putative peptide binding region (PBR). The rate of non-synonymous substitutions (dN) only marginally exceeded synonymous substitutions (dS) in PBR codons. Highly diverged sequences showed no evidence for diversifying selection, possibly because synonymous substitutions have accumulated since their divergence. Also, the substitution pattern of similar alleles did not show evidence for diversifying selection, plausibly because advantageous non-synonymous substitutions have not yet accumulated. Intermediately diverged sequences showed the highest relative rate of non-synonymous substitutions, with dN/dS>14 in some pairwise comparisons. Consequently, a curvilinear relationship was observed between the dN/dS ratio and the level of sequence divergence.     

Major histocompatibility alleles associated with local resistance to malaria in a passerine

Malaria parasites are a major cause of human mortality in tropical countries and a potential threat for wildlife, as witnessed by the malaria-induced extinction of naive Hawaiian avifauna. Identifying resistance mechanisms is therefore crucial both for human health and wildlife conservation. Patterns of malaria resistance are known to be highly polygenic in both humans and mice, with marked contributions attributed to major histocompatibility (Mhc) genes. Here we show that specific Mhc variants are linked to both increased resistance and susceptibility to malaria infection in a wild passerine species, the house sparrow (Passer domesticus). In addition, links between host immunogenetics and resistance to malaria involved population-specific alleles, suggesting local adaptation in this host-parasite interaction. This is the first evidence for a population-specific genetic control of resistance to malaria in a wild species.     

Genetic variation in MHC class II expression and interactions with MHC sequence polymorphism in three-spined sticklebacks

Genes of the major histocompatibility complex (MHC) have been studied for several decades because of their pronounced allelic polymorphism. Structural allelic polymorphism is, however, not the only source of variability subjected to natural selection. Genetic variation may also exist in gene expression patterns. Here, we show that in a natural population of three-spined sticklebacks (Gasterosteus aculeatus) the expression of MHC class IIB genes was positively correlated with parasite load, which indicates increased immune activation of the MHC when infections are frequent. To experimentally study MHC expression, we used laboratory-bred sticklebacks that were exposed to three naturally occurring species of parasite. We found strong differences in MHC class IIB expression patterns among fish families, which were consistent over two generations, thus demonstrating a genetic component. The average number of MHC class IIB sequence variants within families was negatively correlated to the MHC expression level suggesting compensatory up-regulation in fish with a low (i.e. suboptimal) MHC sequence variability. The observed differences among families and the negative correlation with individual sequence diversity imply that MHC expression is evolutionary relevant for the onset and control of the immune response in natural populations.     

MHC class II variation in the endangered European mink <Emphasis Type="Italic">Mustela lutreola</Emphasis> (L. 1761)—consequences for species conservation

The polymorphic major histocompatibility complex (MHC) has gained a specific relevance in pathogen resistance and mate choice. Particularly the antigen-binding site (ABS), encoded by exon 2 of the DRB class II gene, exhibits numerous alleles and extensive sequence variations between alleles. A lack of MHC variability has attributed to instances such as bottleneck effects or relaxed selection pressure and has a certain impact on the long-term viability of the species concerned. As a result of seriously decreased population density during the last century, the current population of the endangered European mink (Mustela lutreola, L. 1761) has suffered from geographic isolation. In this study, we amplified a partial sequence of the MHC class II DRB exon 2 (229 bp), assessed the degree of genetic variation and compared the variability with those of other Mustelidae. As a result, nine alleles were detected in 20 investigated individuals, which differ from each other by four to 25 nucleotide substitutions (two to 11 amino acid substitutions). Whilst an equal ratio for synonymous and non-synonymous substitutions was found inside the ABS, synonymous substitutions were significantly higher than non-synonymous substitutions in the non-ABS region. Results might indicate that no positive selection exists within the ex situ population of M. lutreola, at least in the analysed fragment. In addition, phylogenetic analyses support the trans-species model of evolution. Becker and Nieberg have contributed equally to this work.     

Within-host evolution of CD8+-TL epitopes encoded by overlapping and non-overlapping reading frames of simian immunodeficiency virus

In order to understand the impact of overlapping reading frames on natural selection by host CD8+ T lymphocytes (CD8(+)-TL), we analyzed the pattern of nucleotide substitution in simian immunodeficiency virus (SIV) genomes sampled from populations at time of death in 35 rhesus monkeys. Both the mean number of nonsynonymous nucleotide substitutions per nonsynonymous site (d(N)) and the mean number of synonymous nucleotide substitutions per synonymous site (d(S)) were elevated in overlap regions in comparison to non-overlap regions. Mean d(N) exceeded mean d(S) in CD8(+)-TL epitopes restricted by the host's class I major histocompatibility complex molecules. This pattern, which is indicative of positive Darwinian selection favoring amino acid changes in these epitopes, was seen in both overlap and non-overlap regions; but mean d(N) was particularly elevated in restricted CD8(+)-TL epitopes encoded in overlap regions. Amino acid changes from the inoculum were defined as parallel if the same amino acid change occurred at the same site independently in two or more monkeys, and a surprisingly high proportion (71.9%) of observed amino acid changes throughout the SIV genome occurred in parallel in different monkeys. The proportion of parallel changes in restricted epitopes encoded by overlapping reading frames was still higher (80%), supporting the hypothesis that the interaction of positive selection and overlapping reading frames enhances the probability of convergent or parallel amino acid change.     

Identification of Biased Amino Acid Substitution Patterns in Human Immunodeficiency Virus Type 1 Isolates from Patients Treated with Protease Inhibitors

Human immunodeficiency virus type 1 (HIV-1) amino acid substitutions observed during antiretroviral drug therapy may be caused by drug selection, non-drug-related evolution, or sampling error introduced by the sequencing process. We analyzed HIV-1 sequences from 371 untreated patients and from 178 patients receiving a single protease inhibitor. Amino acid substitution patterns during treatment were compared with inferred substitution patterns arising evolutionarily without treatment. Our results suggest that most treatment-associated amino acid substitutions are caused by selective drug pressure, including substitutions not previously associated with drug resistance.     

A 39-kb sequence around a blackbird Mhc class II gene: ghost of selection past and songbird genome architecture

To gain an understanding of the evolution and genomic context of avian major histocompatibility complex (Mhc) genes, we sequenced a 38.8-kb Mhc-bearing cosmid insert from a red-winged blackbird (Agelaius phoeniceus). The DNA sequence, the longest yet retrieved from a bird other than a chicken, provides a detailed view of the process of gene duplication, divergence, and degeneration ("birth and death") in the avian Mhc, as well as a glimpse into major noncoding features of a songbird genome. The peptide-binding region (PBR) of the single Mhc class II B gene in this region, Agph-DAB2, is almost devoid of polymorphism, and a still-segregating single-base-pair deletion and other features suggest that it is nonfunctional. Agph-DAB2 is estimated to have diverged about 40 MYA from a previously characterized and highly polymorphic blackbird Mhc gene, Aph-DAB1, and is therefore younger than most mammalian Mhc paralogs and arose relatively late in avian evolution. Despite its nonfunctionality, Agph-DAB2 shows very high levels of nonsynonymous divergence from Agph-DAB1 and from reconstructed ancestral sequences in antigen-binding PBR codons-a strong indication of a period of adaptive divergence preceding loss of function. We also found that the region sequenced contains very few other unambiguous genes, a partial Mhc- class II gene fragment, and a paucity of simple-sequence and other repeats. Thus, this sequence exhibits some of the genomic streamlining expected for avian as compared with mammalian genomes, but is not as densely packed with functional genes as is the chicken Mhc.     

Rapid evolution of the MH class I locus results in different allelic compositions in recently diverged populations of Atlantic salmon

We compared major histocompatibility class I allelic diversity in two currently reproductively isolated Atlantic salmon (Salmo salar) populations (Irish and Norwegian) with a common postglacial origin in order to test for among-population differences in allelic composition and patterns of recombination and point mutation. We also examined the evidence for adaptive molecular divergence at this locus by analyzing the rate of amino acid replacement in relation to a neutral expectation. Contrary to our prediction, and in contrast to the situation for other genetic markers, the two populations have almost nonoverlapping sets of major histocompatibility class I alleles. Although there is a strong signal of point mutation that predates population divergence, recent recombination, acting in similar, but not identical, ways in both populations appears to be a significant force in creating new alleles. Moreover, selection acting on peptide-binding residues seems to favor new recombinant alleles and is likely to be responsible for the rapid divergence between populations.