Evolution and ecology of MHC molecules: from genomics to sexual selection

Paternity data reveal high MHC diversity among sires in a polygynandrous, egalitarian primate

Evidence from human and nonhuman primates suggests that females avoid breeding with close kin and may choose mates based on MHC diversity, which can improve offspring survival. In despotic societies, female mate choice may be hindered by male sexual coercion, but in egalitarian societies, females may be less constrained. Among northern muriquis-an egalitarian, polygynandrous primate with male philopatry-analyses of new data on paternity and variation at microsatellite and MHC loci, combined with behavioural and life-history data, revealed that sires showed higher MHC diversity than expected by chance and were never close kin of dams, consistent with predictions of female mate choice and close inbreeding avoidance. However, females did not differentially reproduce with males who were more distantly related to them or more dissimilar at the MHC than expected by chance, nor with those who had more MHC alleles distinct from their own. The lack of male dominance may permit females to identify and reproduce preferentially with non-offspring males and with males who are more diverse at the MHC. Nonetheless, the absence of disassortative mating at the MHC and neutral loci suggests that female mate choice may be limited by other factors impacting male fertilization success.

Immunological tolerance in the evolution of male pregnancy

The unique male pregnancy in pipefishes and seahorses ranges from basic attachment (pouch-less species: Nerophinae) of maternal eggs to specialized internal gestation in pouched species (e.g. Syngnathus and Hippocampus) with many transitions in between. Due to this diversity, male pregnancy offers a unique platform for assessing physiological and molecular adaptations in pregnancy evolution. These insights will contribute to answering long-standing questions of why and how pregnancy evolved convergently in so many vertebrate systems. To understand the molecular congruencies and disparities in male pregnancy evolution, we compared transcriptome-wide differentially expressed genes in four syngnathid species, at four pregnancy stages (nonpregnant, early, late and parturition). Across all species and pregnancy forms, metabolic processes and immune dynamics defined pregnancy stages, especially pouched species shared expression features akin to female pregnancy. The observed downregulation of adaptive immune genes in early-stage pregnancy and its reversed upregulation during late/parturition in pouched species, most notably in Hippocampus, combined with directionless expression in the pouch-less species, suggests immune modulation to be restricted to pouched species that evolved placenta-like systems. We propose that increased foeto-paternal intimacy in pouched syngnathids commands immune suppression processes in early gestation, and that the elevated immune response during parturition coincides with pouch opening and reduced progeny reliance. Immune response regulation in pouched species supports the recently described functional MHC II pathway loss as critical in male pregnancy evolution. The independent co-option of similar genes and pathways both in male and female pregnancy highlights immune modulation as crucial for the evolutionary establishment of pregnancy.

Both selection and drift drive the spatial pattern of adaptive genetic variation in a wild mammal

The major histocompatibility complex (MHC) has been intensively studied for the relative effects of different evolutionary forces in recent decades. Pathogen-mediated balancing selection is generally thought to explain the high polymorphism observed in MHC genes, but it is still unclear to what extent MHC diversity is shaped by selection relative to neutral drift. In this study, we genotyped MHC class II DRB genes and 15 neutral microsatellite loci across 26 geographic populations of European badgers (Meles meles) covering most of their geographic range. By comparing variation of microsatellite and diversity of MHC at different levels, we demonstrate that both balancing selection and drift have shaped the evolution of MHC genes. When only MHC allelic identity was investigated, the spatial pattern of MHC variation was similar to that of microsatellites. By contrast, when functional aspects of the MHC diversity (e.g., immunological supertypes) were considered, balancing selection appears to decrease genetic structuring across populations. Our comprehensive sampling and analytical approach enable us to conclude that the likely mechanisms of selection are heterozygote advantage and/or rare-allele advantage. This study is a clear demonstration of how both balancing selection and genetic drift simultaneously affect the evolution of MHC genes in a widely distributed wild mammal.

Decisive points for pregnancy losses in beef cattle

Beef cattle producers rely on each of their cows to produce a marketable calf each year to maintain a sustainable operation. Within the first month of gestation, pregnancy failures have been recorded to be upwards of 40-50%. From fertilisation to birth, there are numerous factors contributing to pregnancy failure. From the beginning of gestation oocyte competence is often a large factor impacting fertility as the dam contributes all mRNA for initial embryo development. Other factors contributing to early embryonic infertility include hormonal concentration and heat stress. After the embryo enters the uterus, it becomes critical for the uterus to be receptive to the developing conceptus. The embryo then begins to elongate and secrete interferon-tau to initiate maternal recognition of pregnancy; a requirement to establish and maintain bovine pregnancies. After a pregnancy completes these steps, placentation actively begins around day 22 of pregnancy and lasts until organogenesis. The fetal phase follows the embryonic phase where disease and/or toxins are often the cause of pregnancy failure at this period. However, fetal mortality has been reported to occur in less than 10% of pregnancies. Understanding of the many factors influencing infertility needs to be further investigated to increase pregnancy success in beef cattle.

Evaluation of Genetic Diversity and Parasite-Mediated Selection of MHC Class I Genes in Emberiza godlewskii (Passeriformes: Emberizidae)

The major histocompatibility complex (MHC) is a multi-copy immune gene family in vertebrates. Its genes are highly variable and code for antigen-presenting molecules. Characterization of MHC genes in different species and investigating the mechanisms that shape MHC diversity is an important goal in understanding the evolution of biological diversity. Here we developed a next-generation sequencing (NGS) protocol to genotype the MHC class I genes of 326 Godlewski’s buntings (Emberiza godlewskii) sampled in the Western mountain area of Beijing from 2014 to 2016. A total of 184 functional alleles were identified, including both non-classical and classical alleles, clustering into nine supertypes. Compared with other passerine birds, the number of MHC class I alleles per individual in Godlewski’s buntings is high (mean 16.1 ± 3.3, median 16). In addition, we demonstrated signatures of historical and contemporary selection on MHC genes. Reflecting historical selection, ten amino acid sites in the antigen-binding domain showed signatures of balancing selection, eight of which exhibit high amino acid polymorphism. In terms of contemporary selection, we found that specific MHC supertypes were nominally associated with the infection of two malaria parasite lineages. These findings indicate the action of historical and possibly also contemporary balancing selection and suggest negative frequency-dependent or fluctuating selection as possible selection mechanisms.

A genetically linked pair of NLR immune receptors shows contrasting patterns of evolution

Throughout their evolution, plant nucleotide-binding leucine-rich-repeat receptors (NLRs) have acquired widely divergent unconventional integrated domains that enhance their ability to detect pathogen effectors. However, the functional dynamics that drive the evolution of NLRs with integrated domains (NLR-IDs) remain poorly understood. Here, we reconstructed the evolutionary history of an NLR locus prone to unconventional domain integration and experimentally tested hypotheses about the evolution of NLR-IDs. We show that the rice (Oryza sativa) NLR Pias recognizes the effector AVR-Pias of the blast fungal pathogen Magnaporthe oryzae. Pias consists of a functionally specialized NLR pair, the helper Pias-1 and the sensor Pias-2, that is allelic to the previously characterized Pia pair of NLRs: the helper RGA4 and the sensor RGA5. Remarkably, Pias-2 carries a C-terminal DUF761 domain at a similar position to the heavy metal-associated (HMA) domain of RGA5. Phylogenomic analysis showed that Pias-2/RGA5 sensor NLRs have undergone recurrent genomic recombination within the genus Oryza, resulting in up to six sequence-divergent domain integrations. Allelic NLRs with divergent functions have been maintained transspecies in different Oryza lineages to detect sequence-divergent pathogen effectors. By contrast, Pias-1 has retained its NLR helper activity throughout evolution and is capable of functioning together with the divergent sensor-NLR RGA5 to respond to AVR-Pia. These results suggest that opposite selective forces have driven the evolution of paired NLRs: highly dynamic domain integration events maintained by balancing selection for sensor NLRs, in sharp contrast to purifying selection and functional conservation of immune signaling for helper NLRs.

Characterization, Selection, and Trans-Species Polymorphism in the MHC Class II of Heermann’s Gull (Charadriiformes)

The major histocompatibility complex (MHC) enables vertebrates to cope with pathogens and maintain healthy populations, thus making it a unique set of loci for addressing ecology and evolutionary biology questions. The aim of our study was to examine the variability of Heermann’s Gull MHC class II (MHCIIB) and compare these loci with other Charadriiformes. Fifty-nine MHCIIB haplotypes were recovered from sixty-eight Heermann’s Gulls by cloning, of them, twelve were identified as putative true alleles, forty-five as unique alleles, and two as pseudogenes. Intra and interspecific relationships indicated at least two loci in Heermann’s Gull MHCIIB and trans-species polymorphism among Charadriiformes (coinciding with the documented evidence of two ancient avian MHCIIB lineages, except in the Charadriidae family). Additionally, sites under diversifying selection revealed a better match with peptide-binding sites inferred in birds than those described in humans. Despite the negative anthropogenic activity reported on Isla Rasa, Heermann’s Gull showed MHCIIB variability consistent with population expansion, possibly due to a sudden growth following conservation efforts. Duplication must play an essential role in shaping Charadriiformes MHCIIB variability, buffering selective pressures through balancing selection. These findings suggest that MHC copy number and protected islands can contribute to seabird conservation.

MHC class I-α population differentiation in a commercial fish, the European sea bass (Dicentrarchus labrax)

Identifying population structuring in highly fecund marine species with high dispersal rates is challenging, but critical for conservation and stock delimitation for fisheries management. European sea bass (Dicentrarchus labrax) is a commercial species of fisheries and aquaculture relevance whose stocks are declining in the North Atlantic, despite management measures to protect them and identifying their fine population structure is needed for managing their exploitation. As for other marine fishes, neutral genetic markers indicate that eastern Atlantic sea bass form a panmictic population and is currently managed as arbitrarily divided stocks. The genes of the major histocompatibility complex (MHC) are key components of the adaptive immune system and ideal candidates to assess fine structuring arising from local selective pressures. We used Illumina sequencing to characterise allelic composition and signatures of selection at the MHC class I-α region of six D. labrax populations across the Atlantic range. We found high allelic diversity driven by positive selection, corresponding to moderate supertype diversity, with 131 alleles clustering into four to eight supertypes, depending on the Bayesian information criterion threshold applied, and a mean number of 13 alleles per individual. Alleles could not be assigned to particular loci, but private alleles allowed us to detect regional genetic structuring not found previously using neutral markers. Our results suggest that MHC markers can be used to detect cryptic population structuring in marine species where neutral markers fail to identify differentiation. This is particularly critical for fisheries management, and of potential use for selective breeding or identifying escapes from sea farms.

Comparing Class II MHC DRB3 Diversity in Colombian Simmental and Simbrah Cattle Across Worldwide Bovine Populations

The major histocompatibility complex (MHC) exerts great influence on responses to infectious diseases and vaccination due to its fundamental role in the adaptive immune system. Knowledge about MHC polymorphism distribution among breeds can provide insights into cattle evolution and diversification as well as population-based immune response variability, thus guiding further studies. Colombian Simmental and Simbrah cattle’s BoLA-DRB3 genetic diversity was compared to that of taurine and zebuine breeds worldwide to estimate functional diversity. High allele richness was observed for Simmental and Simbrah cattle; nevertheless, high homozygosity was associated with individual low sequence variability in both the β1 domain and the peptide binding region (PBR), thereby implying reduced MHC-presented peptide repertoire size. There were strong signals of positive selection acting on BoLA-DRB3 in all populations, some of which were poorly structured and displayed common alleles accounting for their high genetic similarity. PBR sequence correlation analysis suggested that, except for a few populations exhibiting some divergence at PBR, global diversity regarding potential MHC-presented peptide repertoire could be similar for the cattle populations analyzed here, which points to the retention of functional diversity in spite of the selective pressures imposed by breeding.

Associations between MHC class II variation and phenotypic traits in a free-living sheep population

Pathogen-mediated selection (PMS) is thought to maintain the high level of allelic diversity observed in the major histocompatibility complex (MHC) class II genes. A comprehensive way to demonstrate contemporary selection is to examine associations between MHC variation and individual fitness. As individual fitness is hard to measure, many studies examine associations between MHC variation and phenotypic traits, including direct or indirect measures of adaptive immunity thought to contribute to fitness. Here, we tested associations between MHC class II variation and five phenotypic traits measured in free-living sheep captured in August: weight, strongyle faecal egg count, and plasma IgA, IgE and IgG immunoglobulin titres against the gastrointestinal nematode parasite Teladorsagia circumcincta. We found no association between MHC class II variation and weight or strongyle faecal egg count. We did, however, find associations between MHC class II variation and immunoglobulin levels which varied with isotype, age and sex. Our results suggest associations between MHC and phenotypic traits are more likely to be found for traits more closely associated with pathogen defence than integrative traits such as bodyweight and highlight the association between MHC variation and antibodies in wild populations.

Major histocompatibility complex class II DR and DQ evolution and variation in wild capuchin monkey species (Cebinae)

The major histocompatibility complex (MHC) is an important gene complex contributing to adaptive immunity. Studies of platyrrhine MHC have focused on identifying experimental models of immune system function in the equivalent Human Leukocyte Antigen (HLA). These genes have thus been explored primarily in captive platyrrhine individuals from research colonies. However, investigations of standing MHC variation and evolution in wild populations are essential to understanding its role in immunity, sociality and ecology. Capuchins are a promising model group exhibiting the greatest habitat diversity, widest diet breadth and arguably the most social complexity among platyrrhines, together likely resulting in varied immunological challenges. We use high-throughput sequencing to characterize polymorphism in four Class II DR and DQ exons for the first time in seven capuchin species. We find evidence for at least three copies for DQ genes and at least five for DRB, with possible additional unrecovered diversity. Our data also reveal common genotypes that are inherited across our most widely sampled population, Cebus imitator in Sector Santa Rosa, Costa Rica. Notably, phylogenetic analyses reveal that platyrrhine DQA sequences form a monophyletic group to the exclusion of all Catarrhini sequences examined. This result is inconsistent with the trans-species hypothesis for MHC evolution across infraorders in Primates and provides further evidence for the independent origin of current MHC genetic diversity in Platyrrhini. Identical allele sharing across cebid species, and more rarely genera, however, does underscore the complexity of MHC gene evolution and the need for more comprehensive assessments of allelic diversity and genome structure.

Avian blood parasite richness decreases with major histocompatibility complex class I loci number

Major histocompatibility complex (MHC) genes are among the most polymorphic in the vertebrate genome. The high allele diversity is believed to be maintained primarily by sexual and pathogen-mediated balancing selection. The number of MHC loci also varies greatly across vertebrates, most notably across birds. MHC proteins play key roles in presenting antigens on the cell surface for recognition by T cells, with class I proteins specifically targeting intracellular pathogens. Here, we explore the hypothesis that MHC class I diversity (measured as loci number) coevolves with haemosporidian parasite burden of the host. Using data on 54 bird species, we demonstrate that high-MHC class I diversity is associated with significantly lower richness of Plasmodium, Haemoproteus as well as overall haemosporidian parasite lineages, the former thus indicating more efficient protection against intracellular pathogens. Nonetheless, the latter associations were only detected when MHC diversity was assessed using cloning and not 454 pyrosequencing-based studies, nor across all genotyping methods combined. Our results indicate that high-MHC class I diversity might play a key role in providing qualitative resistance against diverse haemosporidian parasites in birds, but further clarification is needed for the origin of contrasting results when using different genotyping methods for MHC loci quantification.

Nucleotide sequence variation of the major histocompatibility complex class II DQA1 gene in different cattle breeds from Nigeria and New Zealand

The major histocompatibility complex (MHC) plays a role in immune response. Among other activities, the bovine MHC genes (BoLA) trigger immune responses, including the activation of antibody-producing B-cells. In this study, White Fulani (n = 24), Red Bororo (n = 5) and Holstein-White × Fulani-cross (n = 11) cattle from Nigeria, and New Zealand Holstein-Friesian × Jersey-cross (n = 40) cattle were used to investigate variability in exon 2 of BoLA-DQA1. Ten alleles were identified using a PCR-Single Strand Conformation Polymorphism (SSCP) approach and their nucleotide sequences confirmed by DNA sequencing. A total of 12.60 % of all nucleotide positions analysed were revealed to be variable and two novel BoLA-DQA1 alleles are reported here for the first time.

Genotyping and De Novo Discovery of Allelic Variants at the Brassicaceae Self-Incompatibility Locus from Short-Read Sequencing Data

Plant self-incompatibility (SI) is a genetic system that prevents selfing and enforces outcrossing. Because of strong balancing selection, the genes encoding SI are predicted to maintain extraordinarily high levels of polymorphism, both in terms of the number of functionally distinct S-alleles that segregate in SI species and in terms of their nucleotide sequence divergence. However, because of these two combined features, documenting polymorphism of these genes also presents important methodological challenges that have so far largely prevented the comprehensive analysis of complete allelic series in natural populations, and also precluded the obtention of complete genic sequences for many S-alleles. Here, we develop a powerful methodological approach based on a computationally optimized comparison of short Illumina sequencing reads from genomic DNA to a database of known nucleotide sequences of the extracellular domain of SRK (eSRK). By examining mapping patterns along the reference sequences, we obtain highly reliable predictions of S-genotypes from individuals collected from natural populations of Arabidopsis halleri. Furthermore, using a de novo assembly approach of the filtered short reads, we obtain full-length sequences of eSRK even when the initial sequence in the database was only partial, and we discover putative new SRK alleles that were not initially present in the database. When including those new alleles in the reference database, we were able to resolve the complete diploid SI genotypes of all individuals. Beyond the specific case of Brassicaceae S-alleles, our approach can be readily applied to other polymorphic loci, given reference allelic sequences are available.

Genetic variation and selection in the major histocompatibility complex Class II gene in the Guizhou pony

The Guizhou pony (GZP) is an indigenous species of equid found in the mountains of the Guizhou province in southwest China. We selected four regions of the equine leukocyte antigen (ELA), including DQA, DRA, DQB, and DRB, and used them to assess the diversity of the major histocompatibility complex (MHC) class II gene using direct sequencing technology. DRA had the lowest d_N/d_S ratio (0.560) compared with the other three loci, indicating that DRA was conserved and could be conserved after undergoing selective processes. Nine DQA, five DQB, nine DRA, and seven DRB codons were under significant positive selection at the antigen binding sites (ABS), suggesting that the selected residues in ABS may play a significant role in the innate immune system of the GZP. Two GZP alleles were shared with Przewalski’s horse, and six older GZP haplotypes had a better relationship with other horse species by one or two mutational steps, indicating that the GZP may be a natural ancient variety of equid. The specific diversity of ABS and the numbers of unique haplotypes in the evolutionary process affords this species a better genetic fitness and ability to adapt to the native environment.

Allelic diversity and patterns of selection at the major histocompatibility complex class I and II loci in a threatened shorebird, the Snowy Plover (Charadrius nivosus)

Background

Understanding the structure and variability of adaptive loci such as the major histocompatibility complex (MHC) genes is a primary research goal for evolutionary and conservation genetics. Typically, classical MHC genes show high polymorphism and are under strong balancing selection, as their products trigger the adaptive immune response in vertebrates. Here, we assess the allelic diversity and patterns of selection for MHC class I and class II loci in a threatened shorebird with highly flexible mating and parental care behaviour, the Snowy Plover (Charadrius nivosus) across its broad geographic range.

Results

We determined the allelic and nucleotide diversity for MHC class I and class II genes using samples of 250 individuals from eight breeding population of Snowy Plovers. We found 40 alleles at MHC class I and six alleles at MHC class II, with individuals carrying two to seven different alleles (mean 3.70) at MHC class I and up to two alleles (mean 1.45) at MHC class II. Diversity was higher in the peptide-binding region, which suggests balancing selection. The MHC class I locus showed stronger signatures of both positive and negative selection than the MHC class II locus. Most alleles were present in more than one population. If present, private alleles generally occurred at very low frequencies in each population, except for the private alleles of MHC class I in one island population (Puerto Rico, lineage tenuirostris).

Conclusion

Snowy Plovers exhibited an intermediate level of diversity at the MHC, similar to that reported in other Charadriiformes. The differences found in the patterns of selection between the class I and II loci are consistent with the hypothesis that different mechanisms shape the sequence evolution of MHC class I and class II genes. The rarity of private alleles across populations is consistent with high natal and breeding dispersal and the low genetic structure previously observed at neutral genetic markers in this species.

Diversity analysis at MHC class II DQA locus in buffalo (Bubalus bubalis) indicates extensive duplication and trans-species evolution

Variation at MHC Class II-DQA locus in riverine and swamp buffaloes (Bubu) has been explored in this study. Through sequencing of buffalo DQA, 48 nucleotide variants identified from 17 individuals, reporting 42 novel alleles, including one pseudogene. Individual animal displayed two to seven variants, suggesting the presence of more than two Bubu-DQA loci, as an evidence of extensive duplication. dN values were found to be higher than dS values at peptide binding sites, separately for riverine and swamp buffaloes, indicating locus being under positive selection. Evolutionary analysis revealed numerous trans-species polymorphism with alleles from water buffalo assigned to at least three different loci (Bubu-DQA1, DQA2, DQA3). Alleles of both the sub-species intermixed within the cluster, showing convergent evolution of MHC alleles in bovines. The results thus suggest that both riverine and swamp buffaloes share con-current arrangement of DQA region, comparable to cattle in terms of copy number and population polymorphism.

No evidence for a relationship between MHC heterozygosity and life history strategy in a sample of North American undergraduates

Although allelic diversity at the major histocompatibility complex (MHC) has implications for adaptive immunity, mate choice, and social signalling, how diversity at the MHC influences the calibration of life history strategies remains largely uninvestigated. The current study investigated whether greater MHC heterozygosity was associated with markers of slower life history strategies in a sample of 789 North American undergraduates. Contrary to preregistered predictions and to previously published findings, MHC heterozygosity was not related to any of the psychological life history-relevant variables measured (including short- vs. long-term sexual strategy, temporal discounting, the Arizona life history battery, past and current health, disgust sensitivity, and Big Five personality traits). Further, no meaningful effects emerged when analysing women and men separately. Possible reasons for why the current results are inconsistent with previous work are discussed.

Latitudinal diversity gradient and cetaceans from the perspective of MHC genes

Pathogen diversity is a key source of selective pressure on immune system genes, shaping molecular evolution mainly on widely distributed or migratory organisms such as cetaceans. Here, we investigated the effects of latitudinal span migration, different biomes occupation, and pathogen-mediated selection on MHC DQB locus divergence on cetaceans. We applied some evolutionary genetics methods using a dataset of 15 species and 121 sequences, and we found a trend on greater MHC divergence on tropical species when compared with either temperate or migratory species. In addition, oceanic cetaceans exhibit greater MHC divergence. Here, we show that, despite there was a correlation between the diversity of MHC DQB alleles with the distribution of organisms, the pattern of diversity found is not completely explained by pathogenic pressure, suggesting that other factors must be investigated for a better understanding of the processes related to the diversity of MHC in cetaceans.

Maternal killer-cell immunoglobulin-like receptors and paternal human leukocyte antigen ligands in recurrent pregnancy loss cases in Turkey

Objective

The survival of a semi-allogeneic fetus depends on several immunological mechanisms, and it has been suggested that recurrent pregnancy loss (RPL) could develop as a result of one or more immunological abnormalities.

Methods

Compatibility between partners for human leukocyte antigen (HLA) genotypes and the relationships between maternal killer-cell immunoglobulin-like receptor (KIR) and paternal HLA-Bw4/Bw6 and HLA-C1/C2 supra-groups were investigated in 25 couples with RPL in comparison to healthy couples with children. HLA and KIR genotyping was performed using polymerase chain reaction with sequence-specific primers and/or sequence-specific oligonucleotides.

Results

HLA class I incompatibility between partners, especially in HLA-B alleles, was more common in the RPL group (p= 0.01). HLA-C2 homozygosity was more frequent in the male partners of RPL couples than in other groups (p= 0.03). The KIR2DL5 gene frequency was significantly higher in both the female and male partners of RPL couples, whereas the KIR2DS3 gene frequency in male partners of RPL couples was significantly reduced (p= 0.03). The presence of KIR2DL3 in women with RPL was correlated with the presence of HLA-C2 alleles in their spouses (p= 0.03).

Conclusion

Our data from a Turkish population suggest that male HLA-C2 homozygosity may play an important role in RPL. Additionally, an incidental match between male HLA-C2 and female HLA-C1 ligand KIR receptors might perturb the balance between activatory and inhibitory KIR-ligand interactions during pregnancy in couples affected by RPL. The roles of orphan KIR2DL5 and orphan KIR2DS3 in RPL remain obscure.

Major Histocompatibility Complex Class II (DRB3) Genetic Diversity in Spanish Morucha and Colombian Normande Cattle Compared to Taurine and Zebu Populations

Bovine leukocyte antigens (BoLA) have been used as disease markers and immunological traits in cattle due to their primary role in pathogen recognition by the immune system. A higher MHC allele diversity in a population will allow presenting a broader peptide repertoire. However, loss of overall diversity due to domestication process can decrease a population's peptide repertoire. Within the context of zebu and taurine cattle populations, BoLA-DRB3 genetic diversity in Spanish Morucha and Colombian Normande cattle was analyzed and an approach to estimate functional diversity was performed. Sequence-based typing was used for identifying 29, 23, 27, and 28 alleles in Spanish Morucha, Nariño-, Boyacá-, and Cundinamarca-Normande cattle, respectively. These breeds had remarkably low heterozygosity levels and the Hardy-Weinberg principle revealed significant heterozygote deficiency. F_ST and D_A genetic distance showed that Colombian Normande populations had greater variability than other phenotypically homogeneous breeds, such as Holstein. It was also found that Spanish Morucha cattle were strongly differentiated from other cattle breeds. Spanish Morucha had greater divergence in the peptide-binding region regarding other cattle breeds. However, peptide-binding region covariation indicated that the potential peptide repertoire seemed equivalent among cattle breeds. Despite the genetic divergence observed, the extent of the potential peptide repertoire in the cattle populations studied appears to be similar and thus their pathogen recognition potential should be equivalent, suggesting that functional diversity might persist in the face of bottlenecks imposed by domestication and breeding.

Comparative analysis of host resistance to Sarcoptes scabiei var. cuniculi in two different rabbit breeds

BACKGROUND

Scabies, caused by infestation of the mite Sarcoptes scabiei, is one of the most severe ectoparasitic diseases in rabbits. Scabies seriously affects the commercial rabbit breeding, causing severe economic losses. Host resistance to S. scabiei is an important factor in further development of the rabbit industry. In the present study, we compared the host resistance to S. scabiei var. cuniculi of a new breed of domestic rabbit propagated by the Sichuan Animal Sciences Academy (QiXing rabbit, QX) compared with that of a traditional rabbit breed in the domestic rabbit industry (IRA rabbit, IRA).

METHODS

Both QX and IRA rabbits were experimentally infested with live S. scabiei var. cuniculi mites for 48 h. Then, during the course of four-week experimental infestation period, the body weight of rabbits was recorded every two weeks for calculating body-weight variations in comparison to the non-infested control rabbits. Skin lesions in the foot area were assessed on weekly basis and serum samples were tested weekly for the estimation of changes in the total antibody levels (IgG, IgE and IgM). Moreover, DNA extracted from the blood samples was amplified for analysis of the genetic diversity in the major histocompatibility complex, class II, DQ Alpha (MHC-DQA) gene.

RESULTS

Compared to the IRA rabbits, the QX rabbits showed a significantly higher (P < 0.05) relative body weight gain compared to the non-infested control rabbits and significantly lower (P < 0.05) scores for foot skin lesions and higher levels of IgG, IgE and IgM at weeks 1 to 4, week 2 and week 1 post-infestation, respectively. Furthermore, a polymorphism site at position 103 bp of exon two of MHC-DQA gene and a different gene frequency were found between two rabbit breeds, suggesting the genetic basis for the differential host resistance to the S. scabiei var. cuniculi between two rabbit breeds.

CONCLUSIONS

The QX rabbits showed higher host resistance to S. scabiei var. cuniculi compared to the IRA rabbits at the clinical, immunological and genetic levels. These results provide a reference for the breeding of rabbits with adequately improved and sustained host resistance to scabies in the domestic rabbit industry.

MHC structuring and divergent allele advantage in a urodele amphibian: a hierarchical multi-scale approach

Proteins encoded by extraordinarily polymorphic major histocompatibility complex (MHC) genes are involved in the adaptive immune response. Balancing selection is believed to maintain MHC polymorphism in the long term, although neutral processes also play a role in shaping MHC diversity. However, the relative contribution of these processes is poorly understood. Here we characterized MHC class II variation of a low-dispersal, pond-breeding newt (Triturus carnifex) over a restricted, geographically structured area. We aimed to (1) evaluate the contribution of selection and neutral processes to shaping MHC diversity at two geographic scales, and (2) test for signatures of divergent allele advantage (DAA), which is a potentially important mechanism of balancing selection. The dominant role of selection in shaping MHC variation was suggested by the lack of correlation between MHC and neutral (microsatellite) variation. Although most variation occurred within populations for both types of markers, they differed in the extent of structuring at the two spatial scales. MHC structuring was more pronounced at local scales, suggesting the role of local selection, while structuring was not detectable at a larger scale, possibly due to the effect of balancing selection. Microsatellites showed the opposite pattern. As expected under DAA, the observed genotypes combined more sequence diversity than expected under a random association of alleles. Thus, DAA may contribute to maintaining MHC polymorphism, which is ancient, as supported by signatures of historical positive selection and trans-species polymorphism. Our results point to the importance of a multi-scale approach in studying MHC variation, especially in low-dispersal taxa, which are genetically structured at fine spatial scales.

A Developmental Switch Generating Phenotypic Plasticity Is Part of a Conserved Multi-gene Locus

Switching between alternative complex phenotypes is often regulated by "supergenes," polymorphic clusters of linked genes such as in butterfly mimicry. In contrast, phenotypic plasticity results in alternative complex phenotypes controlled by environmental influences rather than polymorphisms. Here, we show that the developmental switch gene regulating predatory versus non-predatory mouth-form plasticity in the nematode Pristionchus pacificus is part of a multi-gene locus containing two sulfatases and two α-N-acetylglucosaminidases (nag). We provide functional characterization of all four genes, using CRISPR-Cas9-based reverse genetics, and show that nag genes and the previously identified eud-1/sulfatase have opposing influences. Members of the multi-gene locus show non-overlapping neuronal expression and epistatic relationships. The locus architecture is conserved in the entire genus Pristionchus. Interestingly, divergence between paralogs is counteracted by gene conversion, as inferred from phylogenies and genotypes of CRISPR-Cas9-induced mutants. Thus, we found that physical linkage accompanies regulatory linkage between switch genes controlling plasticity in P. pacificus.

Balancing Selection Drives the Maintenance of Genetic Variation in Drosophila Antimicrobial Peptides

Genes involved in immune defense against pathogens provide some of the most well-known examples of both directional and balancing selection. Antimicrobial peptides (AMPs) are innate immune effector genes, playing a key role in pathogen clearance in many species, including Drosophila. Conflicting lines of evidence have suggested that AMPs may be under directional, balancing, or purifying selection. Here, we use both a linear model and control-gene-based approach to show that balancing selection is an important force shaping AMP diversity in Drosophila. In Drosophila melanogaster, this is most clearly observed in ancestral African populations. Furthermore, the signature of balancing selection is even more striking once background selection has been accounted for. Balancing selection also acts on AMPs in Drosophila mauritiana, an isolated island endemic separated from D. melanogaster by about 4 Myr of evolution. This suggests that balancing selection may be broadly acting to maintain adaptive diversity in Drosophila AMPs, as has been found in other taxa.

Sperm success and immunity

The moment of the fertilization of an egg by a spermatozoon-the point of "sperm success"-is a key milestone in the biology of sexually reproducing species and is a fundamental requirement for offspring production. Fertilization also represents the culmination of a suite of sexually selected processes in both sexes and is commonly used as a landmark to measure reproductive success. Sperm success is heavily dependent upon interactions with other key aspects of male and female biology, with the immune system among the most important. The immune system is vital to maintaining health in both sexes; however, immune reactions can also have antagonistic effects on sperm success. The effects of immunity on sperm success are diverse, and may include trade-offs in the male between investment in the production or protection of sperm, as well as more direct, hostile, immune responses to sperm within the female, and potentially the male, reproductive tract. Here, we review current understanding of where the biology of immunity and sperm meet, and identify the gaps in our knowledge.

Genomic Organization and Phylogeny of MHC Class II Loci in Cetaceans

Cetaceans are a suborder of secondarily adapted aquatic mammals with an enigmatic history involving a transition from land to sea approximately 55 Mya. During the transition period, cetaceans would have faced many new pathogen challenges, but limited information is available about the adaptive immune system of these mammals. The major histocompatibility complex (MHC) family plays a key role in antigen recognition and presentation in adaptive immunity, which is believed to have evolved in response to pathogens. In the present study, MHC class II loci were characterized in 7 published cetacean genome assemblies and the genomic organization of cetaceans was compared with that of their terrestrial relatives, the cow, sheep, and pig. A total of 9 MHC class II loci were identified in the cetacean genomes: DRA, DRB, DQA, DQB, DPB, DOA, DOB, DMA, and DMB. Sequences from 8 of the 9 genes included intact coding regions and were presumably functional. The organization of the MHC class II loci was conserved across the examined mammalian species, whereas the orientation and number of the alpha and beta genes varied among the species. The phylogenetic reconstruction of all MHC genes from Cetartiodactyla suggested that alpha and beta genes had different topologies. Additionally, based on a phylogenetic reconstruction of the multi-locus DRB, 2 (DRB1 and DRB2) of the 4 putative gene copies were hypothesized to have duplicated and evolved during the radiation of cetaceans.

Alternated selection mechanisms maintain adaptive diversity in different demographic scenarios of a large carnivore

Background

Different population trajectories are expected to impact the signature of neutral and adaptive processes at multiple levels, challenging the assessment of the relative roles of different microevolutionary forces. Here, we integrate adaptive and neutral variability patterns to disentangle how adaptive diversity is driven under different demographic scenarios within the Iberian wolf (Canis lupus) range. We studied the persistent, the expanding and a small, isolated group within the Iberian wolf population, using 3 MHC class II genes (DRB1, DQA1, and DQB1), which diversity was compared with 39 microsatellite loci.

Results

Both the persistent and the expanding groups show evidence of balancing selection, revealed by a significant departure from neutrality at MHC loci, significant higher observed and expected heterozygosity and lower differentiation at MHC than at neutral loci, and signs of positive selection. However, despite exhibiting a significantly higher genetic diversity than the isolated group, the persistent group did not show significant excess of MHC heterozygotes. The expanding group, while showing a similar level of genetic diversity than the persistent group, displays by contrast a significant excess of MHC heterozygotes, which is compatible with the heterozygote advantage mechanism. Results are not clear regarding the role of drift and selection in the isolated group due to the small size of this population. Although diversity indices of MHC loci correspond to neutral expectations in the isolated group, accelerated MHC divergence, revealed by a higher differentiation at MHC than neutral loci, may indicate diversifying selection.

Conclusion

Different selective pressures were observed in the three different demographic scenarios, which are possibly driven by different selection mechanisms to maintain adaptive diversity.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1420-5) contains supplementary material, which is available to authorized users.

Challenges of next-generation sequencing in conservation management: Insights from long-term monitoring of corridor effects on the genetic diversity of mouse lemurs in a fragmented landscape

Long-term genetic monitoring of populations is essential for efforts aimed at preserving genetic diversity of endangered species. Here, we employ a framework of long-term genetic monitoring to evaluate the effects of fragmentation and the effectiveness of the establishment of corridors in restoring population connectivity and genetic diversity of mouse lemurs Microcebus ganzhorni. To this end, we supplement estimates of neutral genetic diversity with the assessment of adaptive genetic variability of the major histocompatibility complex (MHC). In addition, we address the challenges of long-term genetic monitoring of functional diversity by comparing the genotyping performance and estimates of MHC variability generated by single-stranded conformation polymorphism (SSCP)/Sanger sequencing with those obtained by high-throughput sequencing (next-generation sequencing [NGS], Illumina), an issue that is particularly relevant when previous work serves as a baseline for planning management strategies that aim to ensure the viability of a population. We report that SSCP greatly underestimates individual diversity and that discrepancies in estimates of MHC diversity attributable to the comparisons of traditional and NGS genotyping techniques can influence the conclusions drawn from conservation management scenarios. Evidence of migration among fragments in Mandena suggests that mouse lemurs are robust to the process of fragmentation and that the effect of corridors is masked by ongoing gene flow. Nonetheless, results based on a larger number of shared private alleles at neutral loci between fragment pairs found after the establishment of corridors in Mandena suggest that gene flow is augmented as a result of enhanced connectivity. Our data point out that despite low effective population size, M. ganzhorni maintains high individual heterozygosity at neutral loci and at MHC II DRB gene and that selection plays a predominant role in maintaining MHC diversity. These findings highlight the importance of long-term genetic monitoring in order to disentangle between the processes of drift and selection maintaining adaptive genetic diversity in small populations.

Cryptic haplotype-specific gamete selection yields offspring with optimal MHC immune genes

Females choose specific mates in order to produce fitter offspring. However, several factors interfere with females' control over fertilization of their eggs, including sneaker males and phenotypically unpredictable allele segregation during meiosis. Mate choice at the individual level thus provides only a poor approximation for obtaining the best genetic match. Consequently, postcopulatory sperm selection by female oocytes has been proposed as a mechanism to achieve complementary combinations of parental haplotypes. Here, using controlled in vitro fertilization of three‐spined stickleback eggs, we find haplotype‐specific fertilization bias toward gametes with complementary major histocompatibility complex (MHC) immunogenes. The resulting zygote (and thus offspring) genotypes exhibit an intermediate level of individual MHC diversity that was previously shown to confer highest pathogen resistance. Our finding of haplotype‐specific gamete selection thus represents an intriguing mechanism for fine‐tuned optimization of the offspring's immune gene composition and an evolutionary advantage in the Red Queen dynamics of host‐parasite coevolution.

Genomic analysis of MHC-based mate choice in the monogamous California mouse

Variation at Major Histocompatibility Complex (MHC) genes is thought to be an important mechanism underlying mate choice in vertebrates, with individuals typically predicted to prefer MHC-dissimilar reproductive partners. However, analyses based on individual MHC loci have generated contradictory results regarding the role of these genes in mate-choice decisions. To provide a more comprehensive assessment of relationships between MHC variation and mating behavior, we used an exome capture strategy to characterize variability at 13 MHC loci, 312 innate immune system genes, and 1044 nonimmune genes in 25 obligate monogamous pairs of California mice (Peromyscus californicus) from 2 free-living populations of this species in Monterey County, California. Pairwise genotypic comparisons and analyses of SNP-based allelic differences failed to detect disassortative mating based on MHC variability; reproductive partners were not more dissimilar than randomly generated male-female pairs at MHC, innate or nonimmune loci. Within populations, individuals tended to be more closely related at MHC genes than at innate or nonimmune genes. Consistent with the functional role of immunogenes, the 2 study populations were highly differentiated at MHC and innate genes but not at nonimmune loci. Collectively, our results suggest that MHC genetic variation in California mice reflects local differences in pathogen exposure rather than disassortative mating based on variability at MHC Class I and II genes.

MHC class IIα polymorphism and its association with resistance/susceptibility to Vibrio harveyi in golden pompano (Trachinotus ovatus)

The major histocompatibility complex (MHC) plays an important role in the vertebrate immune response to antigenic peptides, and it is essential for recognizing foreign pathogens in organisms. In this study, MHC class IIα (Trov-MHC IIα) from the golden pompano (Trachinotus ovatus) was first cloned and identified. The gene structure of Trov-MHC IIα was contained four exons and three introns. High levels of polymorphism were found in the exon 2 of Trov-MHC IIα. A total of 29 different MHC class IIα alleles with high polymorphism were identified from 80 individuals. The ratio of non-synonymous substitutions (dN) to synonymous substitutions (dS) was 3.157 (>1) in the peptide binding regions (PBRs) of Trov-MHC IIα, suggesting positive balancing selection. Six alleles were selected to analyze the association between alleles and resistance/susceptibility to Vibrio harveyi in golden pompano. The results showed that Trov-DAA*6401 and Trov-DAA*6702 alleles were associated with the resistance to V. harveyi in golden pompano, while alleles Trov-DAA*6304 and Trov-DAA*7301 were associated with the susceptibility to V. harveyi in golden pompano. This study confirmed the association between alleles of MHC class IIα and disease resistance, and also detected some alleles which might be correlated with high V. harveyi-resistance. These disease resistance-related MHC alleles could be used as potential genetic markers for molecular marker-assisted selective breeding in the golden pompano.

Ecology can inform genetics: Disassortative mating contributes to MHC polymorphism in Leach's storm-petrels (Oceanodroma leucorhoa)

Studies of MHC-based mate choice in wild populations often test hypotheses on species exhibiting female choice and male-male competition, which reflects the general prevalence of females as the choosy sex in natural systems. Here, we examined mutual mate-choice patterns in a small burrow-nesting seabird, the Leach's storm-petrel (Oceanodroma leucorhoa), using the major histocompatibility complex (MHC). The life history and ecology of this species are extreme: both partners work together to fledge a single chick during the breeding season, a task that requires regularly travelling hundreds of kilometres to and from foraging grounds over a 6- to 8-week provisioning period. Using a 5-year data set unprecedented for this species (n = 1078 adults and 925 chicks), we found a positive relationship between variation in the likelihood of female reproductive success and heterozygosity at Ocle-DAB2, a MHC class IIB locus. Contrary to previous reports rejecting disassortative mating as a mechanism for maintaining genetic polymorphism in this species, here we show that males make significant disassortative mate-choice decisions. Variability in female reproductive success suggests that the most common homozygous females (Ocle-DAB2*01/Ocle-DAB2*01) may be physiologically disadvantaged and, therefore, less preferred as lifelong partners for choosy males. The results from this study support the role of mate choice in maintaining high levels of MHC variability in a wild seabird species and highlight the need to incorporate a broader ecological framework and sufficient sample sizes into studies of MHC-based mating patterns in wild populations in general.

Characterisation of major histocompatibility complex class I transcripts in an Australian dragon lizard

Characterisation of squamate major histocompatibility complex (MHC) genes has lagged behind other taxonomic groups. MHC genes encode cell-surface glycoproteins that present self- and pathogen-derived peptides to T cells and play a critical role in pathogen recognition. Here we characterise MHC class I transcripts for an agamid lizard (Ctenophorus decresii) and investigate the evolution of MHC class I in Iguanian lizards. An iterative assembly strategy was used to identify six full-length C. decresii MHC class I transcripts, which were validated as likely to encode classical class I MHC molecules. Evidence for exon shuffling recombination was uncovered for C. decresii transcripts and Bayesian phylogenetic analysis of Iguanian MHC class I sequences revealed a pattern expected under a birth-and-death mode of evolution. This work provides a stepping stone towards further research on the agamid MHC class I region.

Selection outweighs drift at a fine scale: Lack of MHC differentiation within a family living lizard across geographically close but disconnected rocky outcrops

The highly polymorphic genes of the major histocompatibility complex (MHC) are involved in disease resistance, mate choice and kin recognition. Therefore, they are widely used markers for investigating adaptive variation. Although selection is the key driver, gene flow and genetic drift also influence adaptive genetic variation, sometimes in opposing ways and with consequences for adaptive potential. To further understand the processes that generate MHC variation, it is helpful to compare variation at the MHC with that at neutral genetic loci. Differences in MHC and neutral genetic variation are useful for inferring the relative influence of selection, gene flow and drift on MHC variation. To date, such investigations have usually been undertaken at a broad spatial scale. Yet, evolutionary and ecological processes can occur at a fine spatial scale, particularly in small or fragmented populations. We investigated spatial patterns of MHC variation among three geographically close, naturally discrete, sampling sites of Egernia stokesii, an Australian lizard. The MHC of E. stokesii has recently been characterized, and there is evidence for historical selection on the MHC. We found E. stokesii MHC weakly differentiated among sites compared to microsatellites, suggesting selection, acting similarly at each site, has outweighed any effects of low gene flow or of genetic drift on E. stokesii MHC variation. Our findings demonstrate the strength of selection in shaping patterns of MHC variation or consistency at a fine spatial scale.

A new hybrid approach for MHC genotyping: high-throughput NGS and long read MinION nanopore sequencing, with application to the non-model vertebrate Alpine chamois (Rupicapra rupicapra)

The major histocompatibility complex (MHC) acts as an interface between the immune system and infectious diseases. Accurate characterization and genotyping of the extremely variable MHC loci are challenging especially without a reference sequence. We designed a combination of long-range PCR, Illumina short-reads, and Oxford Nanopore MinION long-reads approaches to capture the genetic variation of the MHC II DRB locus in an Italian population of the Alpine chamois (Rupicapra rupicapra). We utilized long-range PCR to generate a 9 Kb fragment of the DRB locus. Amplicons from six different individuals were fragmented, tagged, and simultaneously sequenced with Illumina MiSeq. One of these amplicons was sequenced with the MinION device, which produced long reads covering the entire amplified fragment. A pipeline that combines short and long reads resolved several short tandem repeats and homopolymers and produced a de novo reference, which was then used to map and genotype the short reads from all individuals. The assembled DRB locus showed a high level of polymorphism and the presence of a recombination breakpoint. Our results suggest that an amplicon-based NGS approach coupled with single-molecule MinION nanopore sequencing can efficiently achieve both the assembly and the genotyping of complex genomic regions in multiple individuals in the absence of a reference sequence.