Ancestral polymorphism in exon 2 of bluethroat (Luscinia svecica) MHC class II B genes

The genes of the major histocompatibility complex (MHC) are important model genes for understanding selective forces in evolution. Here, we document, using a cloning and sequencing approach, high polymorphism at the exon 2 of the MHC class II B (MHCIIB) genes in the bluethroat (Luscinia svecica); a minimum of 61 unique alleles were detected in 20 individuals, and at least 11 functional loci. In addition, several pseudogenes were revealed. The specimens originated from three different bluethroat subspecies (azuricollis, cyanecula and svecica), and we also analysed four specimens of the closely related thrush nightingale (L. luscinia) for comparison. Phylogenetic analyses of the functional alleles revealed 258 equally parsimonious trees with poor statistical support for the majority of nodes. The distribution of the sequences in the trees point to an ancestral origin of the polymorphism in MHC class II B genes, a portion of which predated the phylogenetic split between the bluethroat and the thrush nightingale. Strong signatures of balancing selection were uncovered for the codons coding for the peptide‐binding residues of the functional MHCIIB exon 2 alleles. Our results highlight the importance of duplication and recombination events for shaping passerine MHC and give insights in the evolutionary dynamics of MHC variation among closely related taxa.     

Major histocompatibility complex class II variation in the giant panda (Ailuropoda melanoleuca)

Habitat destruction and human activity have greatly impacted the natural history of the giant panda (Ailuropoda melanoleuca). Although the genetic diversity of neutral markers has been examined in this endangered species, no previous work has examined adaptive molecular polymorphisms in the giant panda. Here, the major histocompatibility complex (MHC) class II DRB locus was investigated in the giant panda, using single-strand conformation polymorphism (SSCP) and sequence analysis. Comparisons of DNA samples extracted from faecal and blood samples from the same individual revealed that the two materials yielded similar quantities and qualities of DNA, as well as identical SSCP patterns and allelic sequences, demonstrating the reliability of DNA isolation from panda faeces. Analysis of faecal samples from 60 giant pandas revealed relatively low number of alleles: seven alleles. However, the alleles were quite divergent, varying from each other by a range of 7-47 nucleotide substitutions (4-25 amino acid substitutions). Construction of a neighbour-joining tree and comparisons among DRB alleles from other species revealed that both similar and highly divergent alleles survived in the bottlenecked panda populations. Despite species-specific primers used and excellent faecal DNA isolated, a lower level of heterozygosity than expected was still observed due to inbreeding. There were three types of evidence supporting the presence of balancing selection in the giant panda: (i) an obvious excess of nonsynonymous substitutions over synonymous at the antigen-binding positions; (ii) trans-species evolution of two alleles between the giant panda and other felids; and (iii) a more even distribution of alleles than expected from neutrality.     

Spatial pattern of MHC class II variation in the great snipe (Gallinago media)

The genes of the major histocompatibility complex (MHC) code for proteins involved in antigen recognition and triggering of the adaptive immune response, and are therefore likely to be under selection from parasites. These selection regimes may vary in space and time. Here we report a strong geographical structure in MHC class II B genes of a migrating bird, the great snipe (Gallinago media). Genetic differentiation in the MHC between two ecologically distinct distributional regions (Scandinavian mountain populations vs. East European lowland populations) was still present after statistically controlling for the effect of selectively neutral variation (microsatellites) using partial Mantel tests. This suggests a role for selection in generating this spatial structure and that it represents local adaptation to different environments. Differentiation between populations within the two regions was negligible. Overall, we found a high number of MHC alleles (50, from 175 individuals). This, together with a tendency for a higher rate of nonsynonymous than synonymous substitutions in the peptide binding sites, and high Tajima's D in certain regions of the gene, suggests a history of balancing selection. MHC variation is often thought to be maintained by some form of balancing selection, but the nature of this selection remains unclear. Our results support the hypothesis that spatial variation in selection regimes contributes to the high polymorphism.     

On the relative roles of selection and genetic drift in shaping MHC variation

Genes of the major histocompatibility complex (MHC) have provided some of the clearest examples of how natural selection generates discordances between adaptive and neutral variation in natural populations. The type and intensity of selection as well as the strength of genetic drift are believed to be important in shaping the resulting pattern of MHC diversity. However, evaluating the relative contribution of multiple microevolutionary forces is challenging, and empirical studies have reported contrasting results. For instance, balancing selection has been invoked to explain high levels of MHC diversity and low population differentiation in comparison with other nuclear markers. Other studies have shown that genetic drift can sometimes overcome selection and then patterns of genetic variation at adaptive loci cannot be discerned from those occurring at neutral markers. Both empirical and simulated data also indicate that loss of genetic diversity at adaptive loci can occur faster than at neutral loci when selection and population bottlenecks act simultaneously. Diversifying selection, on the other hand, explains accelerated MHC divergence as the result of spatial variation in pathogen‐mediated selective regimes. Because of all these possible scenarios and outcomes, collecting information from as many study systems as possible, is crucial to enhance our understanding about the evolutionary forces driving MHC polymorphism. In this issue, Miller and co‐workers present an illuminating contribution by combining neutral markers (microsatellites) and adaptive MHC class I loci during the investigation of genetic differentiation across island populations of tuatara Sphenodon punctatus. Their study of geographical variation reveals a major role of genetic drift in shaping MHC variation, yet they also discuss some support for diversifying selection.     

Disentangling the roles of natural selection and genetic drift in shaping variation at MHC immunity genes

The major histocompatibility complex (MHC) forms an integral component of the vertebrate immune response and, due to strong selection pressures, is one of the most polymorphic regions of the entire genome. Despite over 15 years of research, empirical studies offer highly contradictory explanations of the relative roles of different evolutionary forces, selection and genetic drift, acting on MHC genes during population bottlenecks. Here, we take a meta-analytical approach to quantify the results of studies into the effects of bottlenecks on MHC polymorphism. We show that the consequences of selection acting on MHC loci prior to a bottleneck event, combined with drift during the bottleneck, will result in overall loss of MHC polymorphism that is ～15% greater than loss of neutral genetic diversity. These results are counter to general expectations that selection should maintain MHC polymorphism, but do agree with the results of recent simulation models and at least two empirical studies. Notably, our results suggest that negative frequency-dependent selection could be more important than overdominance for maintaining high MHC polymorphism in pre-bottlenecked populations.     

Major histocompatibility complex variation at class II DQA locus in the brown hare (Lepus europaeus)

The major histocompatability complex (MHC) is a multigene family of receptors that bind and present antigenic peptides to T-cells. Genes of the MHC are characterized by an outstanding genetic polymorphism, which is considered to be maintained by positive selection. Sites involved in peptide binding form binding pockets (P) that are collectively termed the peptide-binding region (PBR). In this study, we examined the level of MHC genetic diversity within and among natural populations of brown hare ( Lepus europaeus ) from Europe and Anatolia choosing for analysis of the second exon of the DQA locus, one of the most polymorphic class II loci. We aimed at an integrated population genetic analysis of L. europeaus by (i) correlating MHC polymorphism to genetic variability and phylogenetic status estimated previously from maternally (mtDNA) and biparentally (allozymes, microsatellites) inherited loci; and (ii) comparing full-length exon amino acid polymorphism with functional polymorphism in the PBR and the binding pockets P1, P6 and P9. A substantial level of DQA exon 2 polymorphism was detected with two completely different set of alleles between the Anatolian and European populations. However, the phylogeny of full-length exon 2 Leeu-DQA alleles did not show a strong phylogeographic signal. The presence of balancing selection was supported by a statistically significant excess of nonsynonymous substitutions over synonymous in the PBR and a trans-species pattern of evolution detected after phylogenetic reconstruction. The differentiating patterns detected between genetic and functional polymorphism, i.e. the number and the distribution of pocket variants within and among populations, indicated a hierarchical action of selection pressures.     

Disentangling the effects of recombination,selection, and demography on the genetic variation at a major histocompatibility complex class II gene in the alpine chamois

The major histocompatibility complex (MHC) harbours some of the most polymorphic loci in vertebrate genomes. MHC genes are thought to be subject to some form of balancing selection, most likely pathogen‐mediated selection. Hence, MHC genes are excellent candidates for exploring adaptive processes. In this study, we investigated the genetic variation at exon 2 of the DRB class II MHC locus in 191 alpine chamois (Rupicapra rupicapra) from 10 populations in the eastern Alps of Italy. In particular, we were interested in distinguishing and estimating the relative impact of selective and demographic factors, while taking into account the confounding effect of recombination. The extremely high d_n/d_s ratio and the presence of trans‐species polymorphisms suggest that a strong long‐term balancing selection effect has been operating at this locus throughout the evolutionary history of this species. We analysed patterns of genetic variation within and between populations, and the mitochondrial D‐loop polymorphism patterns were analysed to provide a baseline indicator of the effects of demographic processes. These analyses showed that (i) the chamois experienced a demographic decline in the last 5000–30 000 years, most likely related to the postglacial elevation in temperature; (ii) this demographic process can explain the results of neutrality tests applied to MHC variation within populations, but cannot justify the much weaker divergence between populations implied by MHC as opposed to mitochondrial DNA; (iii) similar sets of divergent alleles are probably maintained with similar frequencies by balancing selection in different populations, and this mechanism is also operating in small isolated populations, which are strongly affected by drift.     

Limited polymorphism of the functional MHC class II B gene in the black‐spotted frog (Pelophylax nigromaculatus) identified by locus‐specific genotyping

Amphibians can be more vulnerable to environmental changes and diseases than other species because of their complex life cycle and physiological requirements. Therefore, understanding the adaptation of amphibians to environmental changes is crucial for their conservation. Major histocompatibility complex (MHC) presents an excellent tool for the investigation of adaptive variations and the assessment of adaptive potential because it can be under strong diversifying selection. Here, we isolated the MHC class II B (MHCIIB) gene from cDNA sequences of the black‐spotted frog (Pelophylax nigromaculatus), a widespread amphibian species in China, and designed locus‐specific primers to characterize adaptive variability of this amphibian. Ten alleles were identified from 67 individual frogs of three populations and no more than two alleles were present in each individual animal. Furthermore, none of the sequences had indels or/and stop codons, which is in good agreement with locus‐specific amplification of a functional gene. However, we found low polymorphism at both nucleotide and amino acid levels, even in the antigen‐binding region. Purifying selection acting at this locus was supported by the findings that the d_N/d_S ratio across all alleles was below 1 and that negatively selected sites were detected by different methods. Allele frequency distributions were significantly different among geographic populations, indicating that physiographic factors may have strong effect on the genetic structure of the black‐spotted frog. This study revealed limited polymorphism of three adjacent black‐spotted frog populations at the functional MHCIIB locus, which may be attributed to region‐specific differences. The locus‐specific genotyping technique developed in this study would provide a foundation for future studies on adaptive divergence among different frog populations.     

Stabilizing selection maintains exuberant colour polymorphism in the spider Theridion californicum (Araneae, Theridiidae)

Genetically controlled colour polymorphisms provide a physical manifestation of the operation of selection and how this can vary according to the spatial or temporal arrangement of phenotypes, or their frequency in a population. Here, we examine the role of selection in shaping the exuberant colour polymorphism exhibited by the spider Theridion californicum. This species is part of a system in which several distantly related spiders in the same lineage, but living in very different geographical areas, exhibit remarkably convergent polymorphisms. These polymorphisms are characterized by allelic inheritance and the presence of a single common cryptic morph and, in the case of T. californicum and its congener the Hawaiian happy-face spider Theridion grallator, numerous rare patterned morphs. We compare population differentiation estimated from colour phenotypic data to differentiation at neutral amplified fragment length polymorphisms (AFLP) loci and demonstrate that the colour polymorphism appears to be maintained by balancing selection. We also examine the patterns of selection in the genome-wide sample of AFLP loci and compare approaches to detecting signatures of selection in this context. Our results have important implications regarding balancing selection, suggesting that selective agents can act in a similar manner across disparate taxa in globally disjunct locales resulting in parallel evolution of exuberant polymorphism.     

Extensive variation at MHC DRB in the New Zealand sea lion (Phocarctos hookeri) provides evidence for balancing selection

Marine mammals are often reported to possess reduced variation of major histocompatibility complex (MHC) genes compared with their terrestrial counterparts. We evaluated diversity at two MHC class II B genes, DQB and DRB, in the New Zealand sea lion (Phocarctos hookeri, NZSL) a species that has suffered high mortality owing to bacterial epizootics, using Sanger sequencing and haplotype reconstruction, together with next-generation sequencing. Despite this species'' prolonged history of small population size and highly restricted distribution, we demonstrate extensive diversity at MHC DRB with 26 alleles, whereas MHC DQB is dimorphic. We identify four DRB codons, predicted to be involved in antigen binding, that are evolving under adaptive evolution. Our data suggest diversity at DRB may be maintained by balancing selection, consistent with the role of this locus as an antigen-binding region and the species'' recent history of mass mortality during a series of bacterial epizootics. Phylogenetic analyses of DQB and DRB sequences from pinnipeds and other carnivores revealed significant allelic diversity, but little phylogenetic depth or structure among pinniped alleles; thus, we could neither confirm nor refute the possibility of trans-species polymorphism in this group. The phylogenetic pattern observed however, suggests some significant evolutionary constraint on these loci in the recent past, with the pattern consistent with that expected following an epizootic event. These data may help further elucidate some of the genetic factors underlying the unusually high susceptibility to bacterial infection of the threatened NZSL, and help us to better understand the extent and pattern of MHC diversity in pinnipeds.     

A population genetic study of the endangered plant species Limonium dufourii (Plumbaginaceae) based on amplified fragment length polymorphism (AFLP)

Limonium dufourii ( Plumbaginaceae ) is a triploid species with obligate apomictic reproduction and is endemic to the East Mediterranean coast of Spain, where it is present in only six populations, most of which have a very low number of individuals. Genetic variation and population structure in this species was studied using amplified fragment length polymorphisms (AFLPs) as markers, using the same individuals as in a previous study with random amplified polymorphic DNA (RAPD). Three different primers provided 252 bands of which 51 were polymorphic among the 152 individuals analysed. Those polymorphic bands were able to define 65 different phenotypes, of which all but two were present in only one population. The comparative analyses of data from AFLPs with those from RAPDs show a high degree of concordance. Additionally, and given the nature of these markers, we propose the estimation of nucleotide divergences from AFLP patterns. Relationships among the different AFLP patterns and the estimates of population genetic parameters obtained with this evolutionary distance are in good agreement with previous results. These analyses show that substantial genetic variability and differentiation exist within and among populations of L. dufourii . Their higher reproducibility and the possibility of obtaining estimates of nucleotide divergence make AFLPs a much better DNA fingerprinting technique.     

Extreme copy number variation at a tRNA ligase gene affecting phenology and fitness in yellow monkeyflowers

Copy number variation (CNV) is a major part of the genetic diversity segregating within populations, but remains poorly understood relative to single nucleotide variation. Here, we report on a tRNA ligase gene (Migut.N02091; RLG1a) exhibiting unprecedented, and fitness‐relevant, CNV within an annual population of the yellow monkeyflower Mimulus guttatus. RLG1a variation was associated with multiple traits in pooled population sequencing (PoolSeq) scans of phenotypic and phenological cohorts. Resequencing of inbred lines revealed intermediate‐frequency three‐copy variants of RLG1a (trip+; 5/35 = 14%), and trip+ lines exhibited elevated RLG1a expression under multiple conditions. trip+ carriers, in addition to being over‐represented in late‐flowering and large‐flowered PoolSeq populations, flowered later under stressful conditions in a greenhouse experiment (p < 0.05). In wild population samples, we discovered an additional rare RLG1a variant (high+) that carries 250–300 copies of RLG1a totalling ~5.7 Mb (20–40% of a chromosome). In the progeny of a high+ carrier, Mendelian segregation of diagnostic alleles and qPCR‐based copy counts indicate that high+ is a single tandem array unlinked to the single‐copy RLG1a locus. In the wild, high+ carriers had highest fitness in two particularly dry and/or hot years (2015 and 2017; both p < 0.01), while single‐copy individuals were twice as fecund as either CNV type in a lush year (2016: p < 0.005). Our results demonstrate fluctuating selection on CNVs affecting phenological traits in a wild population, suggest that plant tRNA ligases mediate stress‐responsive life‐history traits, and introduce a novel system for investigating the molecular mechanisms of gene amplification.     

Extraordinary MHC class II B diversity in a non‐passerine,wild bird: the Eurasian Coot Fulica atra (Aves: Rallidae)

The major histocompatibility complex (MHC) hosts the most polymorphic genes ever described in vertebrates. The MHC triggers the adaptive branch of the immune response, and its extraordinary variability is considered an evolutionary consequence of pathogen pressure. The last few years have witnessed the characterization of the MHC multigene family in a large diversity of bird species, unraveling important differences in its polymorphism, complexity, and evolution. Here, we characterize the first MHC class II B sequences isolated from a Rallidae species, the Eurasian Coot Fulica atra. A next‐generation sequencing approach revealed up to 265 alleles that translated into 251 different amino acid sequences (β chain, exon 2) in 902 individuals. Bayesian inference identified up to 19 codons within the presumptive peptide‐binding region showing pervasive evidence of positive, diversifying selection. Our analyses also detected a significant excess of high‐frequency segregating sites (average Tajima's D = 2.36, P < 0.05), indicative of balancing selection. We found one to six different alleles per individual, consistent with the occurrence of at least three MHC class II B gene duplicates. However, the genotypes comprised of three alleles were by far the most abundant in the population investigated (49.4%), followed by those with two (29.6%) and four (17.5%) alleles. We suggest that these proportions are in agreement with the segregation of MHC haplotypes differing in gene copy number. The most widespread segregating haplotypes, according to our findings, would contain one single gene or two genes. The MHC class II of the Eurasian Coot is a valuable system to investigate the evolutionary implications of gene copy variation and extensive variability, the greatest ever found, to the best of our knowledge, in a wild population of a non‐passerine bird.     

Selective pressures on MHC class II genes in the guppy (Poecilia reticulata) as inferred by hierarchical analysis of population structure

Genes of the major histocompatibility complex, which are the most polymorphic of all vertebrate genes, are a pre‐eminent system for the study of selective pressures that arise from host–pathogen interactions. Balancing selection capable of maintaining high polymorphism should lead to the homogenization of MHC allele frequencies among populations, but there is some evidence to suggest that diversifying selection also operates on the MHC. However, the pattern of population structure observed at MHC loci is likely to depend on the spatial and/or temporal scale examined. Here, we investigated selection acting on MHC genes at different geographic scales using Venezuelan guppy populations inhabiting four regions. We found a significant correlation between MHC and microsatellite allelic richness across populations, which suggests the role of genetic drift in shaping MHC diversity. However, compared to microsatellites, more MHC variation was explained by differences between populations within larger geographic regions and less by the differences between the regions. Furthermore, among proximate populations, variation in MHC allele frequencies was significantly higher compared to microsatellites, indicating that selection acting on MHC may increase population structure at small spatial scales. However, in populations that have significantly diverged at neutral markers, the population‐genetic signature of diversifying selection may be eradicated in the long term by that of balancing selection, which acts to preserve rare alleles and thus maintain a common pool of MHC alleles.     

Highly structured nucleotide variation within and among Arabidopsis lyrata populations at the FAH1 and DFR gene regions

Nucleotide variation at the FAH1 and DFR gene regions was surveyed in four populations of Arabidopsis lyrata (two European A. l. petraea and two North American A. l. lyrata populations). In contrast to previous results, levels of variation were not consistently lower in A. l. lyrata than in A. l. petraea, and similar degrees of genetic differentiation were detected between and within subspecies. These observations and the significant genetic differentiation detected among populations suggest population substructure and no real subdivision between subspecies. For each gene studied, genotypic data were obtained, which allowed comparing nucleotide diversity within individuals (between sequences from the same individual) and within populations (between sequences from the same population). The generally lower level of variation within than among individuals detected in each population yielded a significant deviation from panmixia within populations. In three of the four populations studied, two highly divergent alleles were detected within populations at the highly variable DFR locus. This pattern and the significant excess of derived variants detected in most populations suggest that most variation segregating within populations results from rare migration events between relatively small and isolated populations exhibiting reduced panmixia.     

Rapid loss of MHC class II variation in a bottlenecked population is explained by drift and loss of copy number variation

Population bottlenecks may reduce genetic variation and potentially increase the risk of extinction. Here, we present the first study to use historic samples to analyse loss of variation at the major histocompatibility complex (MHC), which plays a central role in vertebrate disease resistance. Balancing selection acts on the MHC and could moderate the loss of variation expected from drift; however, in a Wisconsin population of greater prairie-chickens (Tympanuchus cupido), the number of MHC class II B alleles per individual declined by 44% following a population bottleneck, compared to a loss of only 8% at microsatellites. Simulations indicate that drift likely reduced MHC variation at the population level, as well as within individuals by reducing the number of gene copies per individual or by fixing the same alleles across multiple loci. These multiple effects of genetic drift on MHC variation could have important implications for immunity and fitness.     

Genetic variation within and among fragmented populations of lesser prairie-chickens (Tympanuchus pallidicinctus)

As a result of recurrent droughts and anthropogenic factors, the range of the lesser prairie-chicken (Tympanuchus pallidicinctus) has contracted by 92% and the population has been reduced by approximately 97% in the past century, resulting in the smallest population size and most restricted geographical distribution of any North American grouse. We examined genetic variation through DNA sequence analysis of 478 base pairs of the mitochondrial genome and by assaying allelic variation at five microsatellite loci from lesser prairie-chickens collected on 20 leks in western Oklahoma and east-central New Mexico. Traditional population genetic analyses indicate that lesser prairie-chickens maintain high levels of genetic variation at both nuclear and mitochondrial loci. Although some genetic structuring among lesser prairie-chicken leks was detected within Oklahoma and New Mexico for both nuclear and mitochondrial loci, high levels of differentiation were detected between Oklahoma and New Mexico populations. Nested-clade analysis of mitochondrial haplotypes revealed that both historic and contemporary processes have influenced patterns of haplotype distributions and that historic processes have most likely led to the level of differentiation found between the Oklahoma and New Mexico populations.     

Genetic variation of early height growth traits at the xeric limits of Austrocedrus chilensis (Cupressaceae)

The study of the genetic variation of early height growth traits in seedlings helps to predict the possible outcomes of tree populations in the face of climate change. Second‐year height growth of 10 geographically marginal populations of Patagonian cypress (Austrocedrus chilensis (D. Don) Pic. Ser. et Bizzarri) (Cupressaceae) was characterized under greenhouse conditions. Variation among and within an average of 15 open‐pollinated families (comprising 21 seedlings per family) for each population was analysed for six size and timing traits obtained from fitted Boltzmann growth curves. The among‐family and among‐population variances were 4.03% and 2.74% of the total phenotypic variation, while the residual variance was 84.57% on average. Genetic differentiation among populations was low, except for the maximum growth rate (Q_ST = 0.35) and for growth initiation (Q_ST = 1). For most traits, genetic variation and heritability were variable across populations, except for growth initiation, which showed in general null intra‐population levels of genetic variance. Although no direct associations were found between the additive genetic variation and latitude or altitude, the north range of the distribution was more variable for the pool of the analysed traits. Although most extreme‐marginal populations of A. chilensis would be very limited in their ability to evolve if climate in north‐west Patagonia turns drier and warmer, their long‐term persistence could largely rely on a phenotypic diversification strategy.