Genotyping strategy matters when analyzing hypervariable major histocompatibility complex‐Experience from a passerine bird

Genotyping of classical major histocompatibility complex (MHC) genes is challenging when they are hypervariable and occur in multiple copies. In this study, we used several different approaches to genotype the moderately variable MHC class I exon 3 (MHCIe3) and the highly polymorphic MHC class II exon 2 (MHCIIβe2) in the bluethroat (Luscinia svecica). Two family groups (eight individuals) were sequenced in replicates at both markers using Ion Torrent technology with both a single‐ and a dual‐indexed primer structure. Additionally, MHCIIβe2 was sequenced on Illumina MiSeq. Allele calling was conducted by modifications of the pipeline developed by Sommer et al. (BMC Genomics, 14, 2013, 542) and the software AmpliSAS. While the different genotyping strategies gave largely consistent results for MHCIe3, with a maximum of eight alleles per individual, MHCIIβe2 was remarkably complex with a maximum of 56 MHCIIβe2 alleles called for one individual. Each genotyping strategy detected on average 50%–82% of all MHCIIβe2 alleles per individual, but dropouts were largely allele‐specific and consistent within families for each strategy. The discrepancies among approaches indicate PCR biases caused by the platform‐specific primer tails. Further, AmpliSAS called fewer alleles than the modified Sommer pipeline. Our results demonstrate that allelic dropout is a significant problem when genotyping the hypervariable MHCIIβe2. As these genotyping errors are largely nonrandom and method‐specific, we caution against comparing genotypes across different genotyping strategies. Nevertheless, we conclude that high‐throughput approaches provide a major advance in the challenging task of genotyping hypervariable MHC loci, even though they may not reveal the complete allelic repertoire.     

MHC Class IIB Exon 2 Polymorphism in the Grey Partridge (Perdix perdix) Is Shaped by Selection,Recombination and Gene Conversion

Among bird species, the most studied major histocompatibility complex (MHC) is the chicken MHC. Although the number of studies on MHC in free-ranging species is increasing, the knowledge on MHC variation in species closely related to chicken is required to understand the peculiarities of bird MHC evolution. Here we describe the variation of MHC class IIB (MHCIIB) exon 2 in a population of the Grey partridge (Perdix perdix), a species of high conservation concern throughout Europe and an emerging galliform model in studies of sexual selection. We found 12 alleles in 108 individuals, but in comparison to other birds surprisingly many sites show signatures of historical positive selection. Individuals displayed between two to four alleles both on genomic and complementary DNA, suggesting the presence of two functional MHCIIB loci. Recombination and gene conversion appear to be involved in generating MHCIIB diversity in the Grey partridge; two recombination breakpoints and several gene conversion events were detected. In phylogenetic analysis of galliform MHCIIB, the Grey partridge alleles do not cluster together, but are scattered through the tree instead. Thus, our results indicate that the Grey partridge MHCIIB is comparable to most other galliforms in terms of copy number and population polymorphism.     

Genetic Variation of the Major Histocompatibility Complex (MHC Class II B Gene) in the Threatened Hume’s Pheasant,Syrmaticus humiae

Major histocompatibility complex (MHC) genes are the most polymorphic genes in vertebrates and encode molecules that play a crucial role in pathogen resistance. As a result of their diversity, they have received much attention in the fields of evolutionary and conservation biology. Here, we described the genetic variation of MHC class II B (MHCIIB) exon 2 in a wild population of Hume’s pheasant (Syrmaticus humiae), which has suffered a dramatic decline in population over the last three decades across its ranges in the face of heavy exploitation and habitat loss. Twenty-four distinct alleles were found in 73 S. humiae specimens. We found seven shared alleles among four geographical groups as well as six rare MHCIIB alleles. Most individuals displayed between one to five alleles, suggesting that there are at least three MHCIIB loci of the Hume’s pheasant. The d _N ⁄ d _S ratio at putative antigen-binding sites (ABS) was significantly greater than one, indicating balancing selection is acting on MHCIIB exon 2. Additionally, recombination and gene conversion contributed to generating MHCIIB diversity in the Hume’s pheasant. One to three recombination events and seventy-five significant gene conversion events were observed within the Hume’s pheasant MHCIIB loci. The phylogenetic tree and network analysis revealed that the Hume’s pheasant alleles do not cluster together, but are scattered through the tree or network indicating a trans-species evolutionary mode. These findings revealed the evolution of the Hume’s pheasant MHC after suffering extreme habitat fragmentation.     

Sperm divergence in a passerine contact zone: Indication of reinforcement at the gametic level

Postcopulatory sexual selection may promote evolutionary diversification in sperm form, but the contribution of between‐species divergence in sperm morphology to the origin of reproductive isolation and speciation remains little understood. To assess the possible role of sperm diversification in reproductive isolation, we studied sperm morphology in two closely related bird species, the common nightingale (Luscinia megarhynchos) and the thrush nightingale (Luscinia luscinia), that hybridize in a secondary contact zone spanning Central and Eastern Europe. We found: (1) striking divergence between the species in total sperm length, accompanied by a difference in the length of the mitochondrial sperm component; (2) greater divergence between species in sperm morphology in sympatry than in allopatry, with evidence for character displacement in sperm head length detected in L. megarhynchos; (3) interspecific hybrids showing sperm with a length intermediate between the parental species, but no evidence for decreased sperm quality (the proportion of abnormal spermatozoa in ejaculates). Our results demonstrate that divergence in sperm morphology between the two nightingale species does not result in intrinsic postzygotic isolation, but may contribute to postcopulatory prezygotic isolation. This isolation could be strengthened in sympatry by reinforcement.     

Combining molecular evolution and environmental genomics to unravel adaptive processes of MHC class IIB diversity in European minnows (Phoxinus phoxinus)

Host–pathogen interactions are a major evolutionary force promoting local adaptation. Genes of the major histocompatibility complex (MHC) represent unique candidates to investigate evolutionary processes driving local adaptation to parasite communities. The present study aimed at identifying the relative roles of neutral and adaptive processes driving the evolution of MHC class IIB (MHCIIB) genes in natural populations of European minnows (Phoxinus phoxinus). To this end, we isolated and genotyped exon 2 of two MHCIIB gene duplicates (DAB1 and DAB3) and 1′665 amplified fragment length polymorphism (AFLP) markers in nine populations, and characterized local bacterial communities by 16S rDNA barcoding using 454 amplicon sequencing. Both MHCIIB loci exhibited signs of historical balancing selection. Whereas genetic differentiation exceeded that of neutral markers at both loci, the populations' genetic diversities were positively correlated with local pathogen diversities only at DAB3. Overall, our results suggest pathogen‐mediated local adaptation in European minnows at both MHCIIB loci. While at DAB1 selection appears to favor different alleles among populations, this is only partially the case in DAB3, which appears to be locally adapted to pathogen communities in terms of genetic diversity. These results provide new insights into the importance of host–pathogen interactions in driving local adaptation in the European minnow, and highlight that the importance of adaptive processes driving MHCIIB gene evolution may differ among duplicates within species, presumably as a consequence of alternative selective regimes or different genomic context.     

Sequential organization in the song of thrush nightingale (Luscinia luscinia): clustering and sequential order of the song types

Sequential organization (or syntax) is widely recognized as one of the most intriguing traits of avian song. We examined the singing of thrush nightingale (Luscinia luscinia) with emphasis on the clustering and sequencing of different song types. The individual repertoire size averaged 12.2 ± 3.7 song types (range 7–23 song types). The song types most common in the population were mutually associated within individual repertoires in non-random combinations. The vast majority of transitions between these song types were predominantly unidirectional. During singing, different song types were mutually associated in the fixed vocal programmes. Each programme was composed of 2–6 song types, which were usually delivered in a relatively fixed order. The compositions of some programmes were found to be identical across the whole of Moscow and even far outside it. We suggest that, during song learning, thrush nightingales memorize the programme not as a set of separate song types but as a unitary vocal pattern. Most frequently, the performance of the programs started from the very beginning. The number of complete and incomplete renditions that began from initial song types exceeded considerably the number of renditions that began from any intermediate song types. In the conclusion, we discuss some differences and similarities between song systems of the thrush and common nightingale (Luscinia megarhynchos).     

Characterization of MHC class II B polymorphism in bottlenecked New Zealand saddlebacks reveals low levels of genetic diversity

The major histocompatibility complex (MHC) is integral to the vertebrate adaptive immune system. Characterizing diversity at functional MHC genes is invaluable for elucidating patterns of adaptive variation in wild populations, and is particularly interesting in species of conservation concern, which may suffer from reduced genetic diversity and compromised disease resilience. Here, we use next generation sequencing to investigate MHC class II B (MHCIIB) diversity in two sister taxa of New Zealand birds: South Island saddleback (SIS), Philesturnus carunculatus, and North Island saddleback (NIS), Philesturnus rufusater. These two species represent a passerine family outside the more extensively studied Passerida infraorder, and both have experienced historic bottlenecks. We examined exon 2 sequence data from populations that represent the majority of genetic diversity remaining in each species. A high level of locus co-amplification was detected, with from 1 to 4 and 3 to 12 putative alleles per individual for South and North Island birds, respectively. We found strong evidence for historic balancing selection in peptide-binding regions of putative alleles, and we identified a cluster combining non-classical loci and pseudogene sequences from both species, although no sequences were shared between the species. Fewer total alleles and fewer alleles per bird in SIS may be a consequence of their more severe bottleneck history; however, overall nucleotide diversity was similar between the species. Our characterization of MHCIIB diversity in two closely related species of New Zealand saddlebacks provides an important step in understanding the mechanisms shaping MHC diversity in wild, bottlenecked populations.     

Polymorphisms in major histocompatibility complex class IIα genes are associated with resistance to infectious hematopoietic necrosis in rainbow trout,Oncorhynchus mykiss (Walbaum, 1792)

Infectious hematopoietic necrosis is a serious viral disease of salmonids, including rainbow trout Oncorhynchus mykiss, and causes tremendous economic losses to the rainbow trout farming industry. Major histocompatibility complex (MHC) genes are crucial elements of adaptive immunity in vertebrate organisms and have been linked with the resistance to numerous pathogenic diseases. In this study, polymerase chain reaction‐single strand conformation polymorphism (PCR‐SSCP) followed by cloning and sequencing were used to examine polymorphisms in the DAA genes (specifically DAA exon 2 of MHC class IIα) of rainbow trout and investigate their association with the infectious hematopoietic necrosis virus (IHNV) resistance in rainbow trout. Seventeen alleles were resolved, including 13 novel alleles. Individuals possessed between two and five alleles, indicating that the genome harbours at least three closely‐related DAA exon 2 loci. The ratio of non‐synonymous to synonymous nucleotide substitutions suggested that DAA exon 2 is under positive selection. A greater variability of amino acids and non‐synonymous nucleotide substitution rate was evident in the peptide‐binding region (PBR) than in the non‐PBR (27.75%). Importantly, the analyses revealed that certain MHC class IIα alleles appear to confer resistance to IHNV in rainbow trout, while others confer susceptibility. The most common alleles in the resistant populations of rainbow trout, Onmy‐DAA*1301 and Onmy‐DAA*0304, confer resistance to IHNV and were not present in the susceptible population. Hence, these alleles may be ideal molecular markers that can assist the breeding of IHNV resistance in rainbow trout.     

Extra‐pair mating in a passerine bird with highly duplicated major histocompatibility complex class II: Preference for the golden mean

Genes of the major histocompatibility complex (MHC) are essential in vertebrate adaptive immunity, and they are highly diverse and duplicated in many lineages. While it is widely established that pathogen‐mediated selection maintains MHC diversity through balancing selection, the role of mate choice in shaping MHC diversity is debated. Here, we investigate female mating preferences for MHC class II (MHCII) in the bluethroat (Luscinia svecica), a passerine bird with high levels of extra‐pair paternity and extremely duplicated MHCII. We genotyped family samples with mixed brood paternity and categorized their MHCII alleles according to their functional properties in peptide binding. Our results strongly indicate that females select extra‐pair males in a nonrandom, self‐matching manner that provides offspring with an allelic repertoire size closer to the population mean, as compared to offspring sired by the social male. This is consistent with a compatible genes model for extra‐pair mate choice where the optimal allelic diversity is intermediate, not maximal. This golden mean presumably reflects a trade‐off between maximizing pathogen recognition benefits and minimizing autoimmunity costs. Our study exemplifies how mate choice can reduce the population variance in individual MHC diversity and exert strong stabilizing selection on the trait. It also supports the hypothesis that extra‐pair mating is adaptive through altered genetic constitution in offspring.     

Evolution of MHC-<Emphasis Type="Italic">DRB</Emphasis> class II polymorphism in the genus <Emphasis Type="Italic">Apodemus</Emphasis> and a comparison of <Emphasis Type="Italic">DRB</Emphasis> sequences within the family Muridae (Mammalia: Rodentia)

Allelic diversity at major histocompatibility complex (MHC) genes is thought to be maintained by balancing selection over long periods of time, even across multiple speciation events. Trans-species sharing of MHC alleles among genera has been supported by many studies on mammals and fish, but in rodents, the results are ambiguous. We investigated natural levels of MHC-DRB variability and evolutionary processes in the wood mouse (Apodemus sylvaticus) and the yellow-necked mouse (Apodemus flavicollis), which are common, sympatric murid rodents in European forests. Using single-strand conformation polymorphism analysis and DNA sequencing, 38 DRB exon 2 alleles were detected among 162 A. sylvaticus from nine different locations in Germany and Switzerland, and 15 DRB exon 2 alleles were detected among 60 A. flavicollis from three different locations in northern Germany. There was evidence for balancing selection in both species. Phylogenetic analysis, including additional murid taxa, showed that the DRB exon 2 sequences did not separate according to species, consistent with trans-species evolution of the MHC in these taxa.     

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.     

Avian Clock gene polymorphism: evidence for a latitudinal cline in allele frequencies

In comparison with most animal behaviours, circadian rhythms have a well-characterized molecular genetic basis. Detailed studies of circadian clock genes in 'model' organisms provide a foundation for interpreting the functional and evolutionary significance of polymorphic circadian clock genes found within free-living animal populations. Here, we describe allelic variation in a region of the avian Clock orthologue which encodes a functionally significant polyglutamine repeat (ClkpolyQcds), within free-living populations of two passerine birds, the migratory bluethroat (Luscinia svecica) and the predominantly nonmigratory blue tit (Cyanistes caeruleus). Multiple ClkpolyQcds alleles were found within populations of both species (bluethroat: 12 populations, 7 alleles; blue tit: 14 populations, 9 alleles). Some populations of both species were differentiated at the ClkpolyQcds locus as measured by F(ST) and R(ST) values. Among the blue tit, but not bluethroat populations, we found evidence of latitudinal clines in (i) mean ClkpolyQcds repeat length, and (ii) the proportions of three ClkpolyQcds genotype groupings. Parallel analyses of microsatellite allele frequencies, which are considered to reflect selectively neutral processes, indicate that interpopulation allele frequency variation at the ClkpolyQcds and microsatellite loci does not reflect the same underlying demographic processes. The possibility that the observed interpopulation ClkpolyQcds allele frequency variation is, at least in part, maintained by selection for microevolutionary adaptation to photoperiodic parameters correlated with latitude warrants further study.     

Evolution of a new nonclassical MHC class I locus in two Old World primate species

 HLA-G is a nonclassical major histocompatibility complex (MHC) class I molecule that is expressed only in the human placenta, suggesting that it plays an important role at the fetal-maternal interface. In rhesus monkeys, which have similar placentation to humans, the HLA-G orthologue is a pseudogene. However, rhesus monkeys express a novel placental MHC class I molecule, Mamu-AG, which has HLA-G-like characteristics. Phylogenetic analysis of AG alleles in two Old World primate species, the baboon and the rhesus macaque, revealed limited diversity characteristic of a nonclassical MHC class I locus. Gene trees constructed using classical and nonclassical primate MHC class I alleles demonstrated that the AG locus was most closely related to the classical A locus. Interestingly, gene tree analyses suggested that the AG alleles were most closely related to a subset of A alleles which are the products of an ancestral interlocus recombination event between the A and B loci. Calculation of the rates of synonymous and nonsynonymous substitution at the AG locus revealed that positive selection was not acting on the codons encoding the peptide binding region. In exon 4, however, the rate of nonsynonymous substitution was significantly lower than the rate of synonymous substitution, suggesting that negative selection was acting on these codons. Received: 22 April 1998 / Revised: 15 July 1998     

Extended sequence of the turkey MHC <Emphasis Type="Italic">B</Emphasis>-locus and sequence variation in the highly polymorphic B-G loci

Genetic variation in the major histocompatibility complex (MHC) is directly correlated to differences in disease resistance. Immunity is greatly dependent on highly polymorphic genes in the MHC, such as class I, class II, and class III complement genes. Preliminary studies of wild turkey populations show extreme polymorphisms in a family of genes exclusive to the avian MHC, the class IV or B-G genes. Significance of this variation is unclear as there are few and conflicting studies of the expression of these genes. Confounding understanding of B-G variation is the lack of a complete delineation of the number of loci in the turkey genome. Direct 454 sequencing of a clone from the CHORI-260 BAC library was used to extend the turkey MHC B-locus sequence, identifying five additional complete B-locus genes including two B-G loci. Sequences of the new B-G genes were compared with those of other turkey gene (BG1–3) and sequences available for other galliformes. Phylogenetic analysis shows species-specific gene evolution supporting a birth–death model of evolution for the B-G gene family. Analysis of variation within the signal peptide sequence (exon 1) found two clusters of polymorphism among the turkey B-G genes. Resequencing of exon 1 in a diverse sample including wild, heritage, and commercial turkeys confirmed multiple alleles at each B-G gene. Future studies aim to correlate B-G variation with group and individual immunological differences.     

Genetic diversity of MHC class II <Emphasis Type="Italic">DRB</Emphasis> alleles in the continental and Japanese populations of the sable <Emphasis Type="Italic">Martes zibellina</Emphasis> (Mustelidae,Carnivora, Mammalia)

The sable (Martes zibellina) is a medium-sized mustelid inhabiting forest environments in Siberia, northern China, the Korean Peninsula, and Hokkaido Island, Japan. To further understand the molecular evolution of the major histocompatibility complex (MHC), we sequenced part of exon 2 in MHC class II DRB genes, including codons encoding the antigen binding site, from 33 individuals from continental Eurasia and Japan. We identified 16 MHC class II DRB alleles (Mazi-DRBs), some of which were geographically restricted and others broadly distributed, and eight putative pseudogenes. A single-breakpoint recombination analysis detected a recombination site in the middle of exon 2. A mixed effects model of evolution analysis identified five amino acid sites presumably under positive selection. These sites were all located in the region 3′ to the recombination site, suggesting that positive selection and recombination could be committed to the diversity of the M. zibellina DRB gene. In a Bayesian phylogenetic tree, all Mazi-DRBs and the presumed pseudogenes grouped within a Mustelidae clade. The Mazi-DRBs showed trans-species polymorphism, with some alleles most closely related to alleles from other mustelid species. This result suggests that the sable DRBs have evolved under long-lasting balancing selection.     

Evolution of introns and exons of class II major histocompatibility complex genes of vertebrates

 The phylogenetic relationships and patterns of nucleotide substitution were compared for introns and exons of class II major histocompatibility complex (MHC) genes in three datasets: human DRB1, human DQA1, and cyprinid fish DAB1. In both human DRB1 and cyprinid DAB1, there was strong evidence that recombination events between alleles have occurred in such a way that intron and exon sequences of a given allele do not necessarily share the same evolutionary history. In the case of human DRB1, recombination was found to have homogenized intron 1 and intron 2 sequences relative to exon 2 sequences within lineages of alleles but not between lineages. As a result, mean divergence times of intron sequences are much more recent than those of exonic sequences. Thus, the divergence time of DRB1 introns cannot be used to date that of exons in the same alleles, and the hypothesis that most human DRB1 polymorphism is of very recent origin is not supported. Received: 5 September 1999 / Revised: 30 December 1999     

Sequence Polymorphism and Evolution of Three Cetacean MHC Genes

Sequence variability at three major histocompatibility complex (MHC) genes (DQB, DRA, and MHC-I) of cetaceans was investigated in order to get an overall understanding of cetacean MHC evolution. Little sequence variation was detected at the DRA locus, while extensive and considerable variability were found at the MHC-I and DQB loci. Phylogenetic reconstruction and sequence comparison revealed extensive sharing of identical MHC alleles among different species at the three MHC loci examined. Comparisons of phylogenetic trees for these MHC loci with the trees reconstructed only based on non-PBR sites revealed that allelic similarity/identity possibly reflected common ancestry and were not due to adaptive convergence. At the same time, trans-species evolution was also evidenced that the allelic diversity of the three MHC loci clearly pre-dated species divergence events according to the relaxed molecular clock. It may be the forces of balancing selection acting to maintain the high sequence variability and identical alleles in trans-specific manner at the MHC-I and DQB loci.