1
|
Lenz TL. HLA Genes: A Hallmark of Functional Genetic Variation and Complex Evolution. Methods Mol Biol 2024; 2809:1-18. [PMID: 38907887 DOI: 10.1007/978-1-0716-3874-3_1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/24/2024]
Abstract
The major histocompatibility complex (MHC) with its highly polymorphic HLA genes represents one of the most intensely studied genomic regions in the genome. MHC proteins play a key role in antigen-specific immunity and are associated with a wide range of complex diseases. Despite decades of research and many advances in the field, the characterization and interpretation of its genetic and genomic variability remain challenging. Here an overview is provided of the MHC, the nature of its exceptional variability, and the complex evolutionary processes assumed to drive this variability. Highlighted are also recent advances in the field that promise to improve our understanding of the variability in the MHC and in antigen-specific immunity more generally.
Collapse
Affiliation(s)
- Tobias L Lenz
- Research Unit for Evolutionary Immunogenomics, Department of Biology, University of Hamburg, Hamburg, Germany.
| |
Collapse
|
2
|
Roved J. MHCtools 1.5: Analysis of MHC Sequencing Data in R. Methods Mol Biol 2024; 2809:275-295. [PMID: 38907904 DOI: 10.1007/978-1-0716-3874-3_18] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/24/2024]
Abstract
The genes of the major histocompatibility complex (MHC) play a vital role in the vertebrate immune system and have attracted considerable interest in evolutionary biology. While the MHC has been characterized in detail in humans (human leukocyte antigen, HLA) and in model organisms such as the mouse, studies in non-model organisms often lack prior knowledge about structure, genetic variability, and evolutionary properties of this locus. MHC genotyping in non-model species commonly relies on PCR-based amplicon sequencing, and while several published protocols facilitate generation of MHC sequence data, there is a lack of transparent and standardized tools for downstream data analysis.Here, I present the R package MHCtools version 1.5, which contains 15 tools that (i) assist accurate MHC genotyping from high-throughput amplicon sequencing data, and provide standardized methods to analyze (ii) MHC diversity, (iii) MHC supertypes, and (iv) MHC haplotypes.I hope that MHCtools will be helpful in future studies of the MHC in non-model species and that it may help to advance our understanding of the important roles of the MHC in ecology and evolution.
Collapse
Affiliation(s)
- Jacob Roved
- Section for Molecular Ecology and Evolution, GLOBE Institute, University of Copenhagen, Copenhagen, Denmark.
| |
Collapse
|
3
|
Radwan J, Kohi C, Ejsmond M, Paganini J, Pontarotti P. Integration of the immune memory into the pathogen-driven MHC polymorphism hypothesis. HLA 2023; 102:653-659. [PMID: 37688391 DOI: 10.1111/tan.15216] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2023] [Revised: 06/01/2023] [Accepted: 08/26/2023] [Indexed: 09/10/2023]
Abstract
Major histocompatibility complex (MHC) genes (referred to as human leukocyte antigen or HLA in humans) are a key component of vertebrate immune systems, coding for proteins which present antigens to T-cells. These genes are outstanding in their degree of polymorphism, with important consequences for human and animal health. The polymorphism is thought to arise from selection pressures imposed by pathogens on MHC allomorphs, which differ in their antigen-binding capacity. However, the existing theory has not considered MHC selection in relation to the formation of immune memory. In this paper, we argue that this omission limits our understanding of the evolution of MHC polymorphism and its role in disease. We review recent evidence that has emerged from the massive research effort related to the SARS-CoV-2 pandemics, and which provides new evidence for the role of MHC in shaping immune memory. We then discuss why the inclusion of immune memory within the existing theory may have non-trivial consequence for our understanding of the evolution of MHC polymorphism. Finally, we will argue that neglecting immune memory hinders our interpretation of empirical findings, and postulate that future studies focusing on pathogen-driven MHC selection would benefit from stratifying the available data according to the history of infection (and vaccination, if relevant).
Collapse
Affiliation(s)
- Jacek Radwan
- Evolutionary Biology Group, Faculty of Biology, Adam Mickiewicz University, Poznań, Poland
| | - Chirine Kohi
- MEPHI, Aix Marseille Université, Marseille, France
| | - Maciej Ejsmond
- Institute of Environmental Sciences, Jagiellonian University, Krakow, Poland
| | | | - Pierre Pontarotti
- MEPHI, Aix Marseille Université, Marseille, France
- Institute of Environmental Sciences, Jagiellonian University, Krakow, Poland
- SNC 5039 CNRS, Marseille, France
| |
Collapse
|
4
|
Peters MAE, Mideo N, MacPherson A. The maintenance of genetic diversity under host-parasite coevolution in finite, structured populations. J Evol Biol 2023; 36:1328-1341. [PMID: 37610056 DOI: 10.1111/jeb.14207] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Revised: 06/12/2023] [Accepted: 06/27/2023] [Indexed: 08/24/2023]
Abstract
As a corollary to the Red Queen hypothesis, host-parasite coevolution has been hypothesized to maintain genetic variation in both species. Recent theoretical work, however, suggests that reciprocal natural selection alone is insufficient to maintain variation at individual loci. As highlighted by our brief review of the theoretical literature, models of host-parasite coevolution often vary along multiple axes (e.g. inclusion of ecological feedbacks or abiotic selection mosaics), complicating a comprehensive understanding of the effects of interacting evolutionary processes on diversity. Here we develop a series of comparable models to explore the effect of interactions between spatial structures and antagonistic coevolution on genetic diversity. Using a matching alleles model in finite populations connected by migration, we find that, in contrast to panmictic populations, coevolution in a spatially structured environment can maintain genetic variation relative to neutral expectations with migration alone. These results demonstrate that geographic structure is essential for understanding the effect of coevolution on biological diversity.
Collapse
Affiliation(s)
- Madeline A E Peters
- Department of Microbiology and Molecular Genetics, Michigan State University, East Lansing, Michigan, USA
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Ontario, Canada
| | - Nicole Mideo
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Ontario, Canada
| | - Ailene MacPherson
- Department of Mathematics, Simon Fraser University, Burnaby, British Columbia, Canada
| |
Collapse
|
5
|
Lam DK, Frantz AC, Burke T, Geffen E, Sin SYW. Both selection and drift drive the spatial pattern of adaptive genetic variation in a wild mammal. Evolution 2023; 77:221-238. [PMID: 36626810 DOI: 10.1093/evolut/qpac014] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2021] [Revised: 10/03/2022] [Accepted: 11/04/2022] [Indexed: 01/12/2023]
Abstract
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.
Collapse
Affiliation(s)
- Derek Kong Lam
- School of Biological Sciences, The University of Hong Kong, Hong Kong, China
| | - Alain C Frantz
- Musée National d'Histoire Naturelle, Luxembourg, Luxembourg
| | - Terry Burke
- Ecology and Evolutionary Biology, School of Biosciences, University of Sheffield, Sheffield, United Kingdom
| | - Eli Geffen
- School of Zoology, Tel Aviv University, Tel Aviv, Israel
| | - Simon Yung Wa Sin
- School of Biological Sciences, The University of Hong Kong, Hong Kong, China
| |
Collapse
|
6
|
Cortazar-Chinarro M, Meurling S, Schroyens L, Siljestam M, Richter-Boix A, Laurila A, Höglund J. Major Histocompatibility Complex Variation and Haplotype Associated Survival in Response to Experimental Infection of Two Bd-GPL Strains Along a Latitudinal Gradient. Front Ecol Evol 2022. [DOI: 10.3389/fevo.2022.915271] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
While both innate and adaptive immune system mechanisms have been implicated in resistance against the chytrid fungus Batrachochytrium dendrobatidis (Bd), studies on the role of specific MHC haplotypes on Bd infection are rare. Here, we studied variation in MHC Class IIB loci in the common toad Bufo bufo along a latitudinal gradient across Sweden. In general, Swedish toad populations had few MHC Class IIB haplotypes and MHC diversity declined from south (13 haplotypes) to the north (four haplotypes). The low diversity may compromise the ability of northern populations to fight emerging disease, such as Bd. In a laboratory experiment, we infected newly metamorphosed toads with two strains of the Global Pandemic Lineage of the fungus (Bd-GPL) and compared survival with sham controls. Bd-infected toads had lower survival compared to controls. Moreover, survival was dependent on the Bd-strain and northern toads had lower Bd-mediated survival than southern individuals. MHC diversity was lower in northern toads. All northern experimental animals were monomorphic for a single MHC haplotype, whereas we found seven different haplotypes in southern experimental animals. In southern toads, survival was dependent on both Bd-strain and MHC haplotype suggesting differential infection dynamics depending on both Bd-strain and host immune system characteristics.
Collapse
|
7
|
Roved J, Hansson B, Stervander M, Hasselquist D, Westerdahl H. MHCtools - an R package for MHC high-throughput sequencing data: genotyping, haplotype and supertype inference, and downstream genetic analyses in non-model organisms. Mol Ecol Resour 2022; 22:2775-2792. [PMID: 35587892 PMCID: PMC9543685 DOI: 10.1111/1755-0998.13645] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2020] [Revised: 04/11/2022] [Accepted: 05/11/2022] [Indexed: 11/30/2022]
Abstract
The major histocompatibility complex (MHC) plays a central role in the vertebrate adaptive immune system and has been of long-term interest in evolutionary biology. While several protocols have been developed for MHC genotyping, there is a lack of transparent and standardized tools for downstream analysis of MHC data. Here, we present the R package MHCtools and demonstrate the use of its functions to (i) assist accurate MHC genotyping from high-throughput amplicon-sequencing data, (ii) infer functional MHC supertypes using bootstrapped clustering analysis, (iii) identify segregating MHC haplotypes from family data, and (iv) analyse functional and genetic distances between MHC sequences. We employed MHCtools to analyse MHC class I (MHC-I) amplicon data of 559 great reed warblers (Acrocephalus arundinaceus). We identified 390 MHC-I alleles which clustered into 14 functional supertypes. A phylogenetic analysis and analyses of positive selection suggested that the MHC-I alleles belonged to several distinct functional groups. We furthermore identified 107 segregating haplotypes among 116 families, and found substantial variation in diversity with 4-21 MHC-I alleles and 3-13 MHC-I supertypes per haplotype. Finally, we show that the great reed warbler haplotypes harboured combinations of MHC-I supertypes with greater functional divergence than observed in simulated populations of possible haplotypes, a result that is in accordance with the divergent allele advantage hypothesis. Our study demonstrates the power of MHCtools to support genotyping and analysis of MHC in non-model species, which we hope will encourage broad implementation among researchers in MHC genetics and evolution.
Collapse
Affiliation(s)
- Jacob Roved
- GLOBE Institute, Section for Evolutionary Genomics, University of Copenhagen, 1350, Copenhagen K, Denmark
| | - Bengt Hansson
- Department of Biology, Molecular Ecology and Evolution Lab, Lund University, Ecology Building, 223 62, Lund, Sweden
| | - Martin Stervander
- Department of Biology, Molecular Ecology and Evolution Lab, Lund University, Ecology Building, 223 62, Lund, Sweden.,Department of Biology and Environmental Science, Faculty of Health and Life Sciences, Linnaeus University, 391 82, Kalmar, Sweden.,Bird Group, Department of Life Sciences, Natural History Museum, Akeman Street, Hertfordshire, HP23 6AP, UK
| | - Dennis Hasselquist
- Department of Biology, Molecular Ecology and Evolution Lab, Lund University, Ecology Building, 223 62, Lund, Sweden.,Bird Group, Department of Life Sciences, Natural History Museum, Akeman Street, Hertfordshire, HP23 6AP, UK
| | - Helena Westerdahl
- Department of Biology, Molecular Ecology and Evolution Lab, Lund University, Ecology Building, 223 62, Lund, Sweden
| |
Collapse
|
8
|
Migalska M, Przesmycka K, Alsarraf M, Bajer A, Behnke-Borowczyk J, Grzybek M, Behnke JM, Radwan J. Long term patterns of association between MHC and helminth burdens in the bank vole support Red Queen dynamics. Mol Ecol 2022; 31:3400-3415. [PMID: 35510766 PMCID: PMC9325469 DOI: 10.1111/mec.16486] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2021] [Revised: 03/28/2022] [Accepted: 04/04/2022] [Indexed: 11/28/2022]
Abstract
Major histocompatibility complex (MHC) genes encode proteins crucial for adaptive immunity of vertebrates. Negative frequency-dependent selection (NFDS), resulting from adaptation of parasites to common MHC types, has been hypothesized to maintain high, functionally relevant polymorphism of MHC, but demonstration of this relationship has remained elusive. In particular, differentiation of NFDS from fluctuating selection, resulting from changes in parasite communities in time and space (FS), has proved difficult in short-term studies. Here, we used temporal data, accumulated through long-term monitoring of helminths infecting bank voles (Myodes glareolus), to test specific predictions of NFDS on MHC class II. Data were collected in three, moderately genetically differentiated subpopulations in Poland, which were characterized by some stable spatiotemporal helminth communities but also events indicating introduction of new species and loss of others. We found a complex association between individual MHC diversity and species richness, where intermediate numbers of DRB supertypes correlated with lowest species richness, but the opposite was true for DQB supertypes - arguing against universal selection for immunogenetic optimality. We also showed that particular MHC supertypes explain a portion of the variance in prevalence and abundance of helminths, but this effect was subpopulation-specific, which is consistent with both NFDS and FS. Finally, in line with NFDS, we found that certain helminths that have recently colonized or spread in a given subpopulation, more frequently or intensely infected voles with MHC supertypes that have been common in the recent past. Overall, our results highlight complex spatial and temporal patterns of MHC-parasite associations, the latter being consistent with Red Queen coevolutionary dynamics.
Collapse
Affiliation(s)
- Magdalena Migalska
- Institute of Environmental Sciences, Faculty of Biology, Jagiellonian University, Gronostajowa 7, 30-387, Krakow, Poland
| | - Karolina Przesmycka
- Institute of Environmental Biology, Faculty of Biology, Adam Mickiewicz University in Poznań, Uniwersytetu Poznańskiego 6, 61-614, Poznań, Poland
| | - Mohammed Alsarraf
- Department of Eco-epidemiology of Parasitic Diseases, Institute of Developmental Biology and Biomedical Sciences, Faculty of Biology, University of Warsaw, Miecznikowa 1, 02-096, Warsaw, Poland
| | - Anna Bajer
- Department of Eco-epidemiology of Parasitic Diseases, Institute of Developmental Biology and Biomedical Sciences, Faculty of Biology, University of Warsaw, Miecznikowa 1, 02-096, Warsaw, Poland
| | - Jolanta Behnke-Borowczyk
- Department of Forest Entomology and Pathology, Faculty of Forestry and Wood Technology, Poznań University of Life Sciences, Wojska Polskiego 71c, 60-625, Poznań, Poland
| | - Maciej Grzybek
- Institute of Maritime and Tropical Medicine, Medical University of Gdansk, Powstania Styczniowego 9B, 81-429, Gdynia, Poland
| | - Jerzy M Behnke
- School of Life Science, University of Nottingham, University Park, Nottingham, NG7 2RD, UK
| | - Jacek Radwan
- Institute of Environmental Biology, Faculty of Biology, Adam Mickiewicz University in Poznań, Uniwersytetu Poznańskiego 6, 61-614, Poznań, Poland
| |
Collapse
|
9
|
Slatkin M. Joint estimation of selection intensity and mutation rate under balancing selection with applications to HLA. Genetics 2022; 221:6569836. [PMID: 35435218 PMCID: PMC9157114 DOI: 10.1093/genetics/iyac058] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2022] [Accepted: 04/07/2022] [Indexed: 11/13/2022] Open
Abstract
A composite likelihood method is introduced for jointly estimating the intensity of selection and the rate of mutation, both scaled by the effective population size, when there is balancing selection at a single multi-allelic locus in an isolated population at demographic equilibrium. The performance of the method is tested using simulated data. Average estimated mutation rates and selection intensities are close to the true values but there is considerable variation about the averages. Allowing for both population growth and population subdivision does not result in qualitative differences but the estimated mutation rates and selection intensities do not in general reflect the current effective population size. The method is applied to three class I (HLA-A, HLA-B and HLA-C) and two class II loci (HLA-DRB1 and HLA-DQA1) in the 1000 Genomes populations. Allowing for asymmetric balancing selection has only a slight effect on the results from the symmetric model. Mutations that restore symmetry of the selection model are preferentially retained because of the tendency of natural selection to maximize average fitness. However, slight differences in selective effects result in much longer persistence time of some alleles. Trans-species polymorphism (TSP), which is characteristic of major-histocompatibility loci in vertebrates, is more likely when there are small differences in allelic fitness than when complete symmetry is assumed. Therefore, variation in allelic fitness expands the range of parameter values consistent with observations of TSP.
Collapse
Affiliation(s)
- Montgomery Slatkin
- Department of Integrative Biology, University of California, Berkeley, Berkeley, CA 94720-3140, USA
| |
Collapse
|
10
|
Huang W, Dicks KL, Hadfield JD, Johnston SE, Ballingall KT, Pemberton JM. Contemporary selection on MHC genes in a free-living ruminant population. Ecol Lett 2022; 25:828-838. [PMID: 35050541 PMCID: PMC9306867 DOI: 10.1111/ele.13957] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2021] [Revised: 10/21/2021] [Accepted: 12/08/2021] [Indexed: 11/27/2022]
Abstract
Genes within the major histocompatibility complex (MHC) are the most variable identified in vertebrates. Pathogen-mediated selection is believed to be the main force maintaining MHC diversity. However, relatively few studies have demonstrated contemporary selection on MHC genes. Here, we examine associations between MHC variation and several fitness measurements including total fitness and five fitness components, in 3400 wild Soay sheep (Ovis aries) monitored between 1989 and 2012. In terms of total fitness, measured as lifetime breeding success of all individuals born, we found haplotypes named C and D were associated with decreased and increased male total fitness respectively. In terms of fitness components, juvenile survival was associated with haplotype divergence while individual haplotypes (C, D and F) were associated with adult fitness components. Consistent with the increased male total fitness, the rarest haplotype D has increased in frequency throughout the study period more than expected under neutral expectations. Our results demonstrate contemporary natural selection is acting on MHC class II genes in Soay sheep and the mode of selection on specific fitness components can be different mode from selection on total fitness.
Collapse
Affiliation(s)
- Wei Huang
- Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh, UK
| | - Kara L Dicks
- Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh, UK.,Royal Zoological Society of Scotland, Edinburgh, UK
| | - Jarrod D Hadfield
- Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh, UK
| | - Susan E Johnston
- Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh, UK
| | | | - Josephine M Pemberton
- Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh, UK
| |
Collapse
|
11
|
Ametrano CG, Lumbsch HT, Di Stefano I, Sangvichien E, Muggia L, Grewe F. Should we hail the Red King? Evolutionary consequences of a mutualistic lifestyle in genomes of lichenized ascomycetes. Ecol Evol 2022; 12:e8471. [PMID: 35136549 PMCID: PMC8809443 DOI: 10.1002/ece3.8471] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2021] [Revised: 11/18/2021] [Accepted: 11/19/2021] [Indexed: 11/17/2022] Open
Abstract
The Red Queen dynamic is often brought into play for antagonistic relationships. However, the coevolutionary effects of mutualistic interactions, which predict slower evolution for interacting organisms (Red King), have been investigated to a lesser extent. Lichens are a stable, mutualistic relationship of fungi and cyanobacteria and/or algae, which originated several times independently during the evolution of fungi. Therefore, they represent a suitable system to investigate the coevolutionary effect of mutualism on the fungal genome. We measured substitution rates and selective pressure of about 2000 protein-coding genes (plus the rDNA region) in two different classes of Ascomycota, each consisting of closely related lineages of lichenized and non-lichenized fungi. Our results show that independent lichenized clades are characterized by significantly slower rates for both synonymous and non-synonymous substitutions. We hypothesize that this evolutionary pattern is connected to the lichen life cycle (longer generation time of lichenized fungi) rather than a result of different selection strengths, which is described as the main driver for the Red Kind dynamic. This first empirical evidence of slower evolution in lichens provides an important insight on how biotic cooperative interactions are able to shape the evolution of symbiotic organisms.
Collapse
Affiliation(s)
- Claudio G. Ametrano
- Grainger Bioinformatics Center and Negaunee Integrative Research Center, Science and EducationField Museum of Natural HistoryChicagoIllinoisUSA
| | - H. Thorsten Lumbsch
- Grainger Bioinformatics Center and Negaunee Integrative Research Center, Science and EducationField Museum of Natural HistoryChicagoIllinoisUSA
| | - Isabel Di Stefano
- Grainger Bioinformatics Center and Negaunee Integrative Research Center, Science and EducationField Museum of Natural HistoryChicagoIllinoisUSA
| | - Ek Sangvichien
- Department of BiologyFaculty of ScienceRamkhamhaeng UniversityBangkokThailand
| | | | - Felix Grewe
- Grainger Bioinformatics Center and Negaunee Integrative Research Center, Science and EducationField Museum of Natural HistoryChicagoIllinoisUSA
| |
Collapse
|
12
|
Diversity of MHC IIB genes and parasitism in hybrids of evolutionarily divergent cyprinoid species indicate heterosis advantage. Sci Rep 2021; 11:16860. [PMID: 34413384 PMCID: PMC8376869 DOI: 10.1038/s41598-021-96205-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2020] [Accepted: 08/02/2021] [Indexed: 02/07/2023] Open
Abstract
The genes of the major histocompatibility complex (MHC) are an essential component of the vertebrate immune system and MHC genotypes may determine individual susceptibility to parasite infection. In the wild, selection that favors MHC variability can create situations in which interspecies hybrids experience a survival advantage. In a wild system of two naturally hybridizing leuciscid fish, we assessed MHC IIB genetic variability and its potential relationships to hosts' ectoparasite communities. High proportions of MHC alleles and parasites were species-specific. Strong positive selection at specific MHC codons was detected in both species and hybrids. MHC allele expression in hybrids was slightly biased towards the maternal species. Controlling for a strong seasonal effect on parasite communities, we found no clear associations between host-specific parasites and MHC alleles or MHC supertypes. Hybrids shared more MHC alleles with the more MHC-diverse parental species, but expressed intermediate numbers of MHC alleles and positively selected sites. Hybrids carried significantly fewer ectoparasites than either parent species, suggesting a hybrid advantage via potential heterosis.
Collapse
|
13
|
Özer O, Lenz TL. Unique pathogen peptidomes facilitate pathogen-specific selection and specialization of MHC alleles. Mol Biol Evol 2021; 38:4376-4387. [PMID: 34110412 PMCID: PMC8476153 DOI: 10.1093/molbev/msab176] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
A key component of pathogen-specific adaptive immunity in vertebrates is the presentation of pathogen-derived antigenic peptides by major histocompatibility complex (MHC) molecules. The excessive polymorphism observed at MHC genes is widely presumed to result from the need to recognize diverse pathogens, a process called pathogen-driven balancing selection. This process assumes that pathogens differ in their peptidomes—the pool of short peptides derived from the pathogen’s proteome—so that different pathogens select for different MHC variants with distinct peptide-binding properties. Here, we tested this assumption in a comprehensive data set of 51.9 Mio peptides, derived from the peptidomes of 36 representative human pathogens. Strikingly, we found that 39.7% of the 630 pairwise comparisons among pathogens yielded not a single shared peptide and only 1.8% of pathogen pairs shared more than 1% of their peptides. Indeed, 98.8% of all peptides were unique to a single pathogen species. Using computational binding prediction to characterize the binding specificities of 321 common human MHC class-I variants, we investigated quantitative differences among MHC variants with regard to binding peptides from distinct pathogens. Our analysis showed signatures of specialization toward specific pathogens especially by MHC variants with narrow peptide-binding repertoires. This supports the hypothesis that such fastidious MHC variants might be maintained in the population because they provide an advantage against particular pathogens. Overall, our results establish a key selection factor for the excessive allelic diversity at MHC genes observed in natural populations and illuminate the evolution of variable peptide-binding repertoires among MHC variants.
Collapse
Affiliation(s)
- Onur Özer
- Research Group for Evolutionary Immunogenomics, Max Planck Institute for Evolutionary Biology, 24306 Plön, Germany.,Research Unit for Evolutionary Immunogenomics, Department of Biology, Universität Hamburg, 20146 Hamburg, Germany
| | - Tobias L Lenz
- Research Group for Evolutionary Immunogenomics, Max Planck Institute for Evolutionary Biology, 24306 Plön, Germany.,Research Unit for Evolutionary Immunogenomics, Department of Biology, Universität Hamburg, 20146 Hamburg, Germany
| |
Collapse
|
14
|
Peng F, Ballare KM, Hollis Woodard S, den Haan S, Bolnick DI. What evolutionary processes maintain MHC IIꞵ diversity within and among populations of stickleback? Mol Ecol 2021; 30:1659-1671. [PMID: 33576071 PMCID: PMC8049082 DOI: 10.1111/mec.15840] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2020] [Revised: 01/29/2021] [Accepted: 02/08/2021] [Indexed: 12/01/2022]
Abstract
Major Histocompatibility Complex (MHC) genes code for proteins that recognize foreign protein antigens to initiate T-cell-mediated adaptive immune responses. They are often the most polymorphic genes in vertebrate genomes. How evolution maintains this diversity remains of debate. Three main hypotheses seek to explain the maintenance of MHC diversity by invoking pathogen-mediated selection: heterozygote advantage, frequency-dependent selection, and fluctuating selection across landscapes or through time. Here, we use a large-scale field parasite survey in a stickleback metapopulation to test predictions derived from each of these hypotheses. We identify over 1000 MHC IIβ variants (alleles spanning paralogous genes) and find that many of them covary positively or negatively with parasite load, suggesting that these genes contribute to resistance or susceptibility. However, despite our large sample-size, we find no evidence for the widely cited stabilizing selection on MHC heterozygosity, in which individuals with an intermediate number of MHC variants have the lowest parasite burden. Nor do we observe a rare-variant advantage, or widespread fluctuating selection across populations. In contrast, we find that MHC diversity is best predicted by neutral genome-wide heterozygosity and between-population genomic divergence, suggesting neutral processes are important in shaping the pattern of metapopulation MHC diversity. Thus, although MHC IIβ is highly diverse and relevant to the type and intensity of macroparasite infection in these populations of stickleback, the main models of MHC evolution still provide little explanatory power in this system.
Collapse
Affiliation(s)
- Foen Peng
- Department of Ecology and Evolutionary BiologyUniversity of ConnecticutStorrsCTUSA
| | - Kimberly M. Ballare
- Department of Ecology and Evolutionary BiologyUniversity of CaliforniaSanta CruzCAUSA
| | | | | | - Daniel I. Bolnick
- Department of Ecology and Evolutionary BiologyUniversity of ConnecticutStorrsCTUSA
| |
Collapse
|
15
|
MacPherson A, Keeling MJ, Otto SP. Coevolution fails to maintain genetic variation in a host-parasite model with constant finite population size. Theor Popul Biol 2020; 137:10-21. [PMID: 33340528 DOI: 10.1016/j.tpb.2020.12.001] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2020] [Revised: 11/30/2020] [Accepted: 12/02/2020] [Indexed: 01/17/2023]
Abstract
Coevolutionary negative frequency-dependent selection has been hypothesized to maintain genetic variation in host and parasites. Despite the extensive literature pertaining to host-parasite coevolution, the temporal dynamics of genetic variation have not been examined in a matching-alleles model (MAM) with a finite population size relative to the expectation under neutral genetic drift alone. The dynamics of the MA coevolution in an infinite population, in fact, suggests that genetic variation in these coevolving populations behaves neutrally. By comparing host heterozygosity to the expectation in a single-species model of neutral genetic drift we find that while this is also largely true in finite populations two additional phenomena arise. First, reciprocal natural selection acting on stochastic perturbations in host and pathogen allele frequencies results in a slight increase or decrease in genetic variation depending on the parameter conditions. Second, following the fixation of an allele in the parasite, selection in the MAM becomes directional, which then rapidly erodes genetic variation in the host. Hence, rather than maintain it, we find that, on average, matching-alleles coevolution depletes genetic variation.
Collapse
Affiliation(s)
- Ailene MacPherson
- Department of Zoology and Biodiversity Research Centre, University of British Columbia, Vancouver, Canada.
| | - Matthew J Keeling
- Zeeman Institute of Systems Biology and Infectious Disease Research (SBIDER), University of Warwick, Coventry, United Kingdom
| | - Sarah P Otto
- Department of Zoology and Biodiversity Research Centre, University of British Columbia, Vancouver, Canada
| |
Collapse
|
16
|
Nandakumar M, Ishtiaq F. Genetic drift and bottleneck do not influence diversity in Toll-like receptor genes at a small spatial scale in a Himalayan passerine. Ecol Evol 2020; 10:12246-12263. [PMID: 33209285 PMCID: PMC7663051 DOI: 10.1002/ece3.6855] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2020] [Revised: 07/14/2020] [Accepted: 07/29/2020] [Indexed: 12/15/2022] Open
Abstract
Genetic diversity is important for long-term viability of a population. Low genetic diversity reduces persistence and survival of populations and increases susceptibility to diseases. Comparisons of the neutral markers with functional loci such as immune genes [Toll-like receptors; TLR] can provide useful insights into evolutionary potential of a species and how the diversity of pathogens and selection pressures on their hosts are directly linked to their environment. In this study, we compare genetic diversity in neutral (eleven microsatellite loci) and adaptive (seven TLR loci) loci to determine genetic variation in a nonmigratory western Himalayan passerine, the black-throated tit (Aegithalos concinnus), distributed across an elevation gradient with varying degree of pathogen-mediated selection pressure. We further compare the diversity in TLR loci with a high-elevation sister species, the white-throated tit (Aegithalos niveogularis). Our results indicate a lack of population genetic structure in the black-throated tit and signatures of a past bottleneck. In contrast, we found high diversity in TLR loci and locus-specific (TLR7) signatures of pathogen-mediated selection, which was comparable to diversity in the white-throated tit. Levels of diversity at TLR5 locus corresponded very closely with neutral microsatellite variation. We found evidence of positive selection in TLR1LA, TLR5, and TLR7 loci highlighting the importance in pathogen recognition. Our finding demonstrates that reduction in neutral variation does not necessarily lead to reduction in functional genetic diversity and probably helps in revival of population in a widespread species.
Collapse
Affiliation(s)
- Mridula Nandakumar
- Centre for Ecological SciencesIndian Institute of ScienceBangaloreIndia
- Present address:
Department of BiologyLund UniversityLundSweden
| | - Farah Ishtiaq
- Centre for Ecological SciencesIndian Institute of ScienceBangaloreIndia
| |
Collapse
|
17
|
Fair BJ, Blake LE, Sarkar A, Pavlovic BJ, Cuevas C, Gilad Y. Gene expression variability in human and chimpanzee populations share common determinants. eLife 2020; 9:59929. [PMID: 33084571 PMCID: PMC7644215 DOI: 10.7554/elife.59929] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2020] [Accepted: 10/20/2020] [Indexed: 12/20/2022] Open
Abstract
Inter-individual variation in gene expression has been shown to be heritable and is often associated with differences in disease susceptibility between individuals. Many studies focused on mapping associations between genetic and gene regulatory variation, yet much less attention has been paid to the evolutionary processes that shape the observed differences in gene regulation between individuals in humans or any other primate. To begin addressing this gap, we performed a comparative analysis of gene expression variability and expression quantitative trait loci (eQTLs) in humans and chimpanzees, using gene expression data from primary heart samples. We found that expression variability in both species is often determined by non-genetic sources, such as cell-type heterogeneity. However, we also provide evidence that inter-individual variation in gene regulation can be genetically controlled, and that the degree of such variability is generally conserved in humans and chimpanzees. In particular, we found a significant overlap of orthologous genes associated with eQTLs in both species. We conclude that gene expression variability in humans and chimpanzees often evolves under similar evolutionary pressures.
Collapse
Affiliation(s)
| | - Lauren E Blake
- Department of Human Genetics, University of Chicago, Chicago, United States
| | - Abhishek Sarkar
- Department of Human Genetics, University of Chicago, Chicago, United States
| | - Bryan J Pavlovic
- Department of Neurology, University of California, San Francisco (UCSF), San Francisco, United States
| | - Claudia Cuevas
- Department of Human Genetics, University of Chicago, Chicago, United States
| | - Yoav Gilad
- Department of Medicine, University of Chicago, Chicago, United States.,Department of Human Genetics, University of Chicago, Chicago, United States
| |
Collapse
|
18
|
Balancing selection versus allele and supertype turnover in MHC class II genes in guppies. Heredity (Edinb) 2020; 126:548-560. [PMID: 32985616 DOI: 10.1038/s41437-020-00369-7] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2020] [Revised: 09/03/2020] [Accepted: 09/07/2020] [Indexed: 01/10/2023] Open
Abstract
Selection pressure from parasites is thought to be a major force shaping the extreme polymorphism of the major histocompatibility complex (MHC) genes, but the modes and consequences of selection remain unclear. Here, we analyse MHC class II and microsatellite diversity in 16 guppy (Poecilia reticulata) populations from two islands (Trinidad and Tobago) that have been separated for at least 10 ky. Within-population MHC diversity was high, but allele sharing was limited within islands and even lower between islands, suggesting relatively fast turnover of alleles. Allelic lineages strongly supported in phylogenetic analyses tended to be island-specific, suggesting rapid lineage sorting, and an expansion of an allelic lineage private to Tobago was observed. New alleles appear to be generated locally at a detectably high frequency. We did not detect a consistent signature of local adaptation, but FST outlier analysis suggested that balancing selection may be the more general process behind spatial variation in MHC allele frequencies in this system, particularly within Trinidad. We found no evidence for divergent allele advantage within populations, or for decreased genetic structuring of MHC supertypes compared to MHC alleles. The dynamic and complex nature of MHC evolution we observed in guppies, coupled with some evidence for balancing selection shaping MHC allele frequencies, are consistent with Red Queen processes of host-parasite coevolution.
Collapse
|
19
|
An alternative route of bacterial infection associated with a novel resistance locus in the Daphnia-Pasteuria host-parasite system. Heredity (Edinb) 2020; 125:173-183. [PMID: 32561843 PMCID: PMC7490384 DOI: 10.1038/s41437-020-0332-x] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2019] [Revised: 06/07/2020] [Accepted: 06/07/2020] [Indexed: 11/23/2022] Open
Abstract
To understand the mechanisms of antagonistic coevolution, it is crucial to identify the genetics of parasite resistance. In the Daphnia magna–Pasteuria ramosa host–parasite system, the most important step of the infection process is the one in which P. ramosa spores attach to the host’s foregut. A matching-allele model (MAM) describes the host–parasite genetic interactions underlying attachment success. Here we describe a new P. ramosa genotype, P15, which, unlike previously studied genotypes, attaches to the host’s hindgut, not to its foregut. Host resistance to P15 attachment shows great diversity across natural populations. In contrast to P. ramosa genotypes that use foregut attachment, P15 shows some quantitative variation in attachment success and does not always lead to successful infections, suggesting that hindgut attachment represents a less-efficient infection mechanism than foregut attachment. Using a Quantitative Trait Locus (QTL) approach, we detect two significant QTLs in the host genome: one that co-localizes with the previously described D. magna PR locus of resistance to foregut attachment, and a second, major QTL located in an unlinked genomic region. We find no evidence of epistasis. Fine mapping reveals a genomic region, the D locus, of ~13 kb. The discovery of a second P. ramosa attachment site and of a novel host-resistance locus increases the complexity of this system, with implications for both for the coevolutionary dynamics (e.g., Red Queen and the role of recombination), and for the evolution and epidemiology of the infection process.
Collapse
|
20
|
Radwan J, Babik W, Kaufman J, Lenz TL, Winternitz J. Advances in the Evolutionary Understanding of MHC Polymorphism. Trends Genet 2020; 36:298-311. [DOI: 10.1016/j.tig.2020.01.008] [Citation(s) in RCA: 106] [Impact Index Per Article: 26.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2019] [Revised: 01/13/2020] [Accepted: 01/14/2020] [Indexed: 12/26/2022]
|
21
|
Märkle H, Tellier A. Inference of coevolutionary dynamics and parameters from host and parasite polymorphism data of repeated experiments. PLoS Comput Biol 2020; 16:e1007668. [PMID: 32203545 PMCID: PMC7156111 DOI: 10.1371/journal.pcbi.1007668] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2019] [Revised: 04/14/2020] [Accepted: 01/19/2020] [Indexed: 01/27/2023] Open
Abstract
There is a long-standing interest in understanding host-parasite coevolutionary dynamics and associated fitness effects. Increasing amounts of genomic data for both interacting species offer a promising source to identify candidate loci and to infer the main parameters of the past coevolutionary history. However, so far no method exists to perform the latter. By coupling a gene-for-gene model with coalescent simulations, we first show that three types of biological costs, namely, resistance, infectivity and infection, define the allele frequencies at the internal equilibrium point of the coevolution model. These in return determine the strength of selective signatures at the coevolving host and parasite loci. We apply an Approximate Bayesian Computation (ABC) approach on simulated datasets to infer these costs by jointly integrating host and parasite polymorphism data at the coevolving loci. To control for the effect of genetic drift on coevolutionary dynamics, we assume that 10 or 30 repetitions are available from controlled experiments or several natural populations. We study two scenarios: 1) the cost of infection and population sizes (host and parasite) are unknown while costs of infectivity and resistance are known, and 2) all three costs are unknown while populations sizes are known. Using the ABC model choice procedure, we show that for both scenarios, we can distinguish with high accuracy pairs of coevolving host and parasite loci from pairs of neutrally evolving loci, though the statistical power decreases with higher cost of infection. The accuracy of parameter inference is high under both scenarios especially when using both host and parasite data because parasite polymorphism data do inform on costs applying to the host and vice-versa. As the false positive rate to detect pairs of genes under coevolution is small, we suggest that our method complements recently developed methods to identify host and parasite candidate loci for functional studies.
Collapse
Affiliation(s)
- Hanna Märkle
- Section of Population Genetics, TUM School of Life Sciences Weihenstephan, Technical University of Munich, Freising, Germany
| | - Aurélien Tellier
- Section of Population Genetics, TUM School of Life Sciences Weihenstephan, Technical University of Munich, Freising, Germany
| |
Collapse
|
22
|
Arora J, Pierini F, McLaren PJ, Carrington M, Fellay J, Lenz TL. HLA Heterozygote Advantage against HIV-1 Is Driven by Quantitative and Qualitative Differences in HLA Allele-Specific Peptide Presentation. Mol Biol Evol 2020; 37:639-650. [PMID: 31651980 PMCID: PMC7038656 DOI: 10.1093/molbev/msz249] [Citation(s) in RCA: 45] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Pathogen-mediated balancing selection is regarded as a key driver of host immunogenetic diversity. A hallmark for balancing selection in humans is the heterozygote advantage at genes of the human leukocyte antigen (HLA), resulting in improved HIV-1 control. However, the actual mechanism of the observed heterozygote advantage is still elusive. HLA heterozygotes may present a broader array of antigenic viral peptides to immune cells, possibly resulting in a more efficient cytotoxic T-cell response. Alternatively, heterozygosity may simply increase the chance to carry the most protective HLA alleles, as individual HLA alleles are known to differ substantially in their association with HIV-1 control. Here, we used data from 6,311 HIV-1-infected individuals to explore the relative contribution of quantitative and qualitative aspects of peptide presentation in HLA heterozygote advantage against HIV. Screening the entire HIV-1 proteome, we observed that heterozygous individuals exhibited a broader array of HIV-1 peptides presented by their HLA class I alleles. In addition, viral load was negatively correlated with the breadth of the HIV-1 peptide repertoire bound by an individual's HLA variants, particularly at HLA-B. This suggests that heterozygote advantage at HLA-B is at least in part mediated by quantitative peptide presentation. We also observed higher HIV-1 sequence diversity among HLA-B heterozygous individuals, suggesting stronger evolutionary pressure from HLA heterozygosity. However, HLA heterozygotes were also more likely to carry certain HLA alleles, including the highly protective HLA-B*57:01 variant, indicating that HLA heterozygote advantage ultimately results from a combination of quantitative and qualitative effects in antigen presentation.
Collapse
Affiliation(s)
- Jatin Arora
- Research Group for Evolutionary Immunogenomics, Max Planck Institute for Evolutionary Biology, Plön, Germany
| | - Federica Pierini
- Research Group for Evolutionary Immunogenomics, Max Planck Institute for Evolutionary Biology, Plön, Germany
| | - Paul J McLaren
- JC Wilt Infectious Diseases Research Center, National HIV and Retrovirology Laboratory, Public Health Agency of Canada, Winnipeg, MB, Canada
- Department of Medical Microbiology and Infectious Diseases, University of Manitoba, Winnipeg, MB, Canada
| | - Mary Carrington
- Basic Science Program, Frederick National Laboratory for Cancer Research, Frederick, MD
- Ragon Institute of Massachusetts General Hospital, Massachusetts Institute of Technology and Harvard University, Cambridge, MA
| | - Jacques Fellay
- Global Health Institute, School of Life Sciences, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
- Precision Medicine Unit, Lausanne University Hospital and University of Lausanne, Lausanne, Switzerland
| | - Tobias L Lenz
- Research Group for Evolutionary Immunogenomics, Max Planck Institute for Evolutionary Biology, Plön, Germany
| |
Collapse
|
23
|
Bentkowski P, Radwan J. Mating preferences can drive expansion or contraction of major histocompatibility complex gene family. Proc Biol Sci 2020; 287:20192706. [PMID: 32097586 DOI: 10.1098/rspb.2019.2706] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
Major histocompatibility complex (MHC)-based mating rules can evolve as a way to avoid inbreeding or to increase offspring immune competence. While the role of mating preference in shaping the MHC diversity in vertebrates has been acknowledged, its impact on individual MHC diversity has not been considered. Here, we use computer simulations to investigate how simple mating rules favouring MHC-dissimilar partners affect the evolution of the number of MHC variants in individual genomes, accompanying selection for resistance to parasites. We showed that the effect of such preferences could sometimes be dramatic. If preferences are aimed at avoiding identical alleles, the equilibrium number of MHC alleles is much smaller than under random mating. However, if the mating rule minimizes the ratio of shared to different alleles in partners, MHC number is higher than under random mating. Additionally, our simulations revealed that a negative correlation between the numbers of MHC variants in mated individuals can arise from simple rules of MHC-disassortative mating. Our results reveal unexpected potential of MHC-based mating preferences to drive MHC gene family expansions or contractions and highlight the need to study the mechanistic basis of such preferences, which is currently poorly understood.
Collapse
Affiliation(s)
- Piotr Bentkowski
- Evolutionary Biology Group, Faculty of Biology, Adam Mickiewicz University in Poznań, Poznań 61-614, Poland
| | - Jacek Radwan
- Evolutionary Biology Group, Faculty of Biology, Adam Mickiewicz University in Poznań, Poznań 61-614, Poland
| |
Collapse
|
24
|
Jaworska J, Ropka-Molik K, Wocławek-Potocka I, Siemieniuch M. Inter- and intrabreed diversity of the major histocompatibility complex (MHC) in primitive and draft horse breeds. PLoS One 2020; 15:e0228658. [PMID: 32012208 PMCID: PMC6996847 DOI: 10.1371/journal.pone.0228658] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2019] [Accepted: 01/21/2020] [Indexed: 12/14/2022] Open
Abstract
Background Polymorphism of major histocompatibility complex (MHC) genes ensures effective immune responses against a wide array of pathogens. However, artificial selection, as performed in the case of domestic animals, may influence MHC diversity. Here, we investigate and compare the MHC diversity of three populations of horses, for which different breeding policies were applied, to evaluate the impact of artificial selection and the environment on MHC polymorphism. Methods Samples of DNA were taken from 100 Polish draft horses, 38 stabled Konik Polski horses and 32 semiferal Konik Polski horses. MHC alleles and haplotype diversity within and between these populations of horses was estimated from 11 MHC microsatellite loci. Results MHC diversity measured based on allelic richness, observed heterozygosity, expected heterozygosity and polymorphism content was similar across the MHC microsatellite loci in all three populations. The highest expected heterozygosity was detected in semiferal primitive horses (He = 0.74), while the lowest was calculated for draft horses (He = 0.65). In total, 203 haplotypes were determined (111 in Polish draft horses, 43 in semiferal Konik Polski horses and 49 in stabled Konik Polski horses), and four haplotypes were shared between the two populations of Koniks. None of these haplotypes were present in any of the previously investigated horse breeds. Intra-MHC recombination events were detected in all three populations. However, the population of semiferal Konik horses showed the highest recombination frequency among the three horse populations. In addition, three recombination events were detected. Conclusions These results showed that despite the different breeding policies, the MHC allele and haplotype diversity was similarly high in all three horse populations. Nevertheless, the proportion of new haplotypes in the offspring was the highest in semiferal Konik Polski horses, which indicates the influence of the environment on MHC diversity in horses. Thus, we speculate that the genetic makeup of the domestic horse MHC might be more strongly influenced by the environment than by artificial selection. Moreover, intra-MHC conversion, insertion, and deletion and intra-MHC recombination may be proposed as mechanisms underlying the generation of new MHC haplotypes in horses.
Collapse
Affiliation(s)
- Joanna Jaworska
- Department of Gamete and Embryo Biology, Institute of Animal Reproduction and Food Research, Polish Academy of Sciences, Olsztyn, Poland
- * E-mail:
| | - Katarzyna Ropka-Molik
- Department of Animal Molecular Biology, National Research Institute of Animal Production, Balice, Poland
| | - Izabela Wocławek-Potocka
- Department of Gamete and Embryo Biology, Institute of Animal Reproduction and Food Research, Polish Academy of Sciences, Olsztyn, Poland
| | - Marta Siemieniuch
- Research Station of the Institute of Reproduction and Food Research, Polish Academy of Sciences in Popielno, Ruciane-Nida, Poland
| |
Collapse
|
25
|
Vlček J, Štefka J. Association between louse abundance and MHC II supertypes in Galápagos mockingbirds. Parasitol Res 2020; 119:1597-1605. [PMID: 32006226 DOI: 10.1007/s00436-020-06617-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2019] [Accepted: 01/22/2020] [Indexed: 01/30/2023]
Abstract
Major histocompatibility complex class II (MHC II) is an essential molecule triggering the adaptive immune response by the presentation of pathogens to helper T cells. The association between individual MHC II variants and various parasites has become a frequent finding in studies of vertebrate populations. However, although bird ectoparasites have a significant effect on their host's fitness, and the host's immune system can regulate ectoparasitic infections, no study has yet investigated the association between MHC II polymorphism and ectoparasite infection in the populations of free-living birds. Here, we test whether an association exists between the abundance of a chewing louse (Myrsidea nesomimi) and MHC II polymorphism of its hosts, the Galápagos mockingbirds (Mimus). We have found that the presence of two MHC II supertypes (functionally differentiated clusters) was significantly associated with louse abundance. This pattern supports the theory that a co-evolutionary interaction stands behind the maintenance of MHC polymorphism. Moreover, we have found a positive correlation between louse abundance and heterophil/lymphocyte ratio (an indicator of immunological stress) that serves as an additional piece of evidence that ectoparasite burden is affected by immunological state of Galápagos mockingbirds.
Collapse
Affiliation(s)
- Jakub Vlček
- Institute of Parasitology, Biology Centre CAS, Branišovská 1160/31, 37005, České Budějovice, Czech Republic. .,Faculty of Science, University of South Bohemia, Branišovská, 1760, České Budějovice, Czech Republic.
| | - Jan Štefka
- Institute of Parasitology, Biology Centre CAS, Branišovská 1160/31, 37005, České Budějovice, Czech Republic.,Faculty of Science, University of South Bohemia, Branišovská, 1760, České Budějovice, Czech Republic
| |
Collapse
|
26
|
Talarico L, Babik W, Marta S, Pietrocini V, Mattoccia M. MHC structuring and divergent allele advantage in a urodele amphibian: a hierarchical multi-scale approach. Heredity (Edinb) 2019; 123:593-607. [PMID: 31036951 PMCID: PMC6972932 DOI: 10.1038/s41437-019-0221-3] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2018] [Revised: 03/31/2019] [Accepted: 03/31/2019] [Indexed: 12/27/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Lorenzo Talarico
- Department of Biology, University of Rome Tor Vergata, Via della Ricerca Scientifica 1, Rome, 00133, Italy.
| | - Wiesław Babik
- Institute of Environmental Sciences, Jagiellonian University, Gronostajowa 7, Kraków, 30-387, Poland
| | - Silvio Marta
- Department of Environmental Sciences and Policy, University of Milan, Via G. Celoria 26, Milan, 20133, Italy
| | - Venusta Pietrocini
- Department of Biology, University of Rome Tor Vergata, Via della Ricerca Scientifica 1, Rome, 00133, Italy
| | - Marco Mattoccia
- Department of Biology, University of Rome Tor Vergata, Via della Ricerca Scientifica 1, Rome, 00133, Italy
| |
Collapse
|
27
|
Dudek K, Gaczorek TS, Zieliński P, Babik W. Massive introgression of major histocompatibility complex (MHC) genes in newt hybrid zones. Mol Ecol 2019; 28:4798-4810. [PMID: 31574568 DOI: 10.1111/mec.15254] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2019] [Revised: 09/25/2019] [Accepted: 09/26/2019] [Indexed: 12/22/2022]
Abstract
Variation in the vertebrate major histocompatibility complex (MHC) genes is crucial for fighting pathogen assault. Because new alleles confer a selective advantage, MHC should readily introgress between species, even under limited hybridization. Using replicated transects through two hybrid zones between strongly reproductively isolated European newts, Lissotriton montandoni and L. vulgaris, we demonstrated recent and ongoing MHC class I and II introgression in the Carpathian region. The extent of introgression correlated with the age of contact. In the older zone, MHC similarity between species within transects exceeded similarity between transects within species, implying pervasive introgression - a massive exchange of MHC genes, not limited to specific variants. In simulations, the observed pattern emerged under the combined action of balancing selection and hybridization, but not when these processes acted separately. Thus, massive introgression at advanced stages of divergence can introduce novel and restore previously lost MHC variation, boosting the adaptive potential of hybridizing taxa. In consequence, MHC genes may be the last to stop introgressing between incipient species.
Collapse
Affiliation(s)
- Katarzyna Dudek
- Institute of Environmental Sciences, Faculty of Biology, Jagiellonian University, Kraków, Poland
| | - Tomasz S Gaczorek
- Institute of Environmental Sciences, Faculty of Biology, Jagiellonian University, Kraków, Poland
| | - Piotr Zieliński
- Institute of Environmental Sciences, Faculty of Biology, Jagiellonian University, Kraków, Poland
| | - Wiesław Babik
- Institute of Environmental Sciences, Faculty of Biology, Jagiellonian University, Kraków, Poland
| |
Collapse
|
28
|
Sella G, Barton NH. Thinking About the Evolution of Complex Traits in the Era of Genome-Wide Association Studies. Annu Rev Genomics Hum Genet 2019; 20:461-493. [DOI: 10.1146/annurev-genom-083115-022316] [Citation(s) in RCA: 123] [Impact Index Per Article: 24.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Many traits of interest are highly heritable and genetically complex, meaning that much of the variation they exhibit arises from differences at numerous loci in the genome. Complex traits and their evolution have been studied for more than a century, but only in the last decade have genome-wide association studies (GWASs) in humans begun to reveal their genetic basis. Here, we bring these threads of research together to ask how findings from GWASs can further our understanding of the processes that give rise to heritable variation in complex traits and of the genetic basis of complex trait evolution in response to changing selection pressures (i.e., of polygenic adaptation). Conversely, we ask how evolutionary thinking helps us to interpret findings from GWASs and informs related efforts of practical importance.
Collapse
Affiliation(s)
- Guy Sella
- Department of Biological Sciences, Columbia University, New York, NY 10027, USA
- Department of Systems Biology, Columbia University, New York, NY 10032, USA
- Program for Mathematical Genomics, Columbia University, New York, NY 10032, USA
| | - Nicholas H. Barton
- Institute of Science and Technology Austria, 3400 Klosterneuburg, Austria
| |
Collapse
|
29
|
Strotz LC, Simões M, Girard MG, Breitkreuz L, Kimmig J, Lieberman BS. Getting somewhere with the Red Queen: chasing a biologically modern definition of the hypothesis. Biol Lett 2019; 14:rsbl.2017.0734. [PMID: 29720444 DOI: 10.1098/rsbl.2017.0734] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2017] [Accepted: 04/13/2018] [Indexed: 01/24/2023] Open
Abstract
The Red Queen hypothesis (RQH) is both familiar and murky, with a scope and range that has broadened beyond its original focus. Although originally developed in the palaeontological arena, it now encompasses many evolutionary theories that champion biotic interactions as significant mechanisms for evolutionary change. As such it de-emphasizes the important role of abiotic drivers in evolution, even though such a role is frequently posited to be pivotal. Concomitant with this shift in focus, several studies challenged the validity of the RQH and downplayed its propriety. Herein, we examine in detail the assumptions that underpin the RQH in the hopes of furthering conceptual understanding and promoting appropriate application of the hypothesis. We identify issues and inconsistencies with the assumptions of the RQH, and propose a redefinition where the Red Queen's reign is restricted to certain types of biotic interactions and evolutionary patterns occurring at the population level.
Collapse
Affiliation(s)
- Luke C Strotz
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA .,Biodiversity Institute, University of Kansas, Lawrence, KS 66045, USA
| | - Marianna Simões
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA.,Biodiversity Institute, University of Kansas, Lawrence, KS 66045, USA
| | - Matthew G Girard
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA.,Biodiversity Institute, University of Kansas, Lawrence, KS 66045, USA
| | - Laura Breitkreuz
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA.,Biodiversity Institute, University of Kansas, Lawrence, KS 66045, USA
| | - Julien Kimmig
- Biodiversity Institute, University of Kansas, Lawrence, KS 66045, USA
| | - Bruce S Lieberman
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA.,Biodiversity Institute, University of Kansas, Lawrence, KS 66045, USA
| |
Collapse
|
30
|
Evolution of major histocompatibility complex gene copy number. PLoS Comput Biol 2019; 15:e1007015. [PMID: 31095555 PMCID: PMC6541284 DOI: 10.1371/journal.pcbi.1007015] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2018] [Revised: 05/29/2019] [Accepted: 04/09/2019] [Indexed: 01/05/2023] Open
Abstract
MHC genes, which code for proteins responsible for presenting pathogen-derived antigens to the host immune system, show remarkable copy-number variation both between and within species. However, the evolutionary forces driving this variation are poorly understood. Here, we use computer simulations to investigate whether evolution of the number of MHC variants in the genome can be shaped by the number of pathogen species the host population encounters (pathogen richness). Our model assumed that while increasing a range of pathogens recognised, expressing additional MHC variants also incurs costs such as an increased risk of autoimmunity. We found that pathogen richness selected for high MHC copy number only when the costs were low. Furthermore, the shape of the association was modified by the rate of pathogen evolution, with faster pathogen mutation rates selecting for increased host MHC copy number, but only when pathogen richness was low to moderate. Thus, taking into account factors other than pathogen richness may help explain wide variation between vertebrate species in the number of MHC genes. Within population, variation in the number of unique MHC variants carried by individuals (INV) was observed under most parameter combinations, except at low pathogen richness. This variance gave rise to positive correlations between INV and host immunocompetence (proportion of pathogens recognised). However, within-population variation in host immunocompetence declined with pathogen richness. Thus, counterintuitively, pathogens can contribute more to genetic variance for host fitness in species exposed to fewer pathogen species, with consequences to predictions from “Hamilton-Zuk” theory of sexual selection. Highly polymorphic genes of the Major Histocompatibility Complex (MHC) code for proteins responsible for presenting antigens to lymphocytes, thus initiating adaptive immune response. The polymorphism is driven by coevolution with parasites which are selected to evade recognition by MHC proteins. Expressing many MHC molecules could ensure that an individual could present antigens of most pathogen species encountered, but this comes at a cost, such as enhanced negative selection on lymphocytes leading to holes in T-cell receptor repertoire. Our simulations showed that evolution of the number of MHC genes in the genome is driven by a complex interaction between three factors we explored: pathogen richness, the intrinsic cost of expressing additional MHC variants, and pathogen mutation rate. In contrast to verbal arguments, our results indicate that pathogen richness does not always selects for MHC gene family expansion. Taking into account factors other than pathogen richness, in particular costs of expressing additional MHC variants which are still poorly understood, may help explain striking interspecific variation in the number of MHC genes. Counterintuitively, our results also demonstrated that opportunity for selection on immunocompetence should decrease with MHC gene family expansion.
Collapse
|
31
|
Manczinger M, Boross G, Kemény L, Müller V, Lenz TL, Papp B, Pál C. Pathogen diversity drives the evolution of generalist MHC-II alleles in human populations. PLoS Biol 2019; 17:e3000131. [PMID: 30703088 PMCID: PMC6372212 DOI: 10.1371/journal.pbio.3000131] [Citation(s) in RCA: 48] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2018] [Revised: 02/12/2019] [Accepted: 01/15/2019] [Indexed: 02/03/2023] Open
Abstract
Central players of the adaptive immune system are the groups of proteins encoded in the major histocompatibility complex (MHC), which shape the immune response against pathogens and tolerance to self-peptides. The corresponding genomic region is of particular interest, as it harbors more disease associations than any other region in the human genome, including associations with infectious diseases, autoimmune disorders, cancers, and neuropsychiatric diseases. Certain MHC molecules can bind to a much wider range of epitopes than others, but the functional implication of such an elevated epitope-binding repertoire has remained largely unclear. It has been suggested that by recognizing more peptide segments, such promiscuous MHC molecules promote immune response against a broader range of pathogens. If so, the geographical distribution of MHC promiscuity level should be shaped by pathogen diversity. Three lines of evidence support the hypothesis. First, we found that in pathogen-rich geographical regions, humans are more likely to carry highly promiscuous MHC class II DRB1 alleles. Second, the switch between specialist and generalist antigen presentation has occurred repeatedly and in a rapid manner during human evolution. Third, molecular positions that define promiscuity level of MHC class II molecules are especially diverse and are under positive selection in human populations. Taken together, our work indicates that pathogen load maintains generalist adaptive immune recognition, with implications for medical genetics and epidemiology. Whereas specialist major histocompatibility complex (MHC) molecules initiate immune response against only relatively few pathogens, generalists provide protection against a broad range. Accordingly, this study shows that the geographical distribution of generalist MHC alleles in human populations reflects exposure to diverse infectious diseases. Variation in the human genome influences our susceptibility to infectious diseases, but the causal link between disease and underlying mutation often remains enigmatic. Major histocompatibility complex II (MHC class II) molecules shape both our immune response against pathogens and our tolerance of self-peptides. The genomic region that encodes MHC molecules is of particular interest, as it is home to more genetic disease associations than any other region in the human genome, including associations with infectious diseases, autoimmune disorders, cancers, and neuropsychiatric diseases. Here, we propose that MHC class II molecules can be categorized into two major types; specialists initiate effective immune response against only relatively few pathogens, while generalists provide protection against a broad range of pathogens. As support, we demonstrate that generalist MHC class II variants are more prevalent in human populations residing in pathogen-rich areas.
Collapse
Affiliation(s)
- Máté Manczinger
- Synthetic and Systems Biology Unit, Institute of Biochemistry, Biological Research Centre, Hungarian Academy of Sciences, Szeged, Hungary
- Department of Dermatology and Allergology, University of Szeged, Szeged, Hungary
- MTA-SZTE Dermatological Research Group, University of Szeged, Szeged, Hungary
| | - Gábor Boross
- Synthetic and Systems Biology Unit, Institute of Biochemistry, Biological Research Centre, Hungarian Academy of Sciences, Szeged, Hungary
| | - Lajos Kemény
- Department of Dermatology and Allergology, University of Szeged, Szeged, Hungary
- MTA-SZTE Dermatological Research Group, University of Szeged, Szeged, Hungary
| | - Viktor Müller
- Institute of Biology, Eötvös Loránd University, Budapest, Hungary
| | - Tobias L. Lenz
- Research Group for Evolutionary Immunogenomics, Max Planck Institute for Evolutionary Biology, Plön, Germany
| | - Balázs Papp
- Synthetic and Systems Biology Unit, Institute of Biochemistry, Biological Research Centre, Hungarian Academy of Sciences, Szeged, Hungary
- * E-mail: (CP); (BP)
| | - Csaba Pál
- Synthetic and Systems Biology Unit, Institute of Biochemistry, Biological Research Centre, Hungarian Academy of Sciences, Szeged, Hungary
- * E-mail: (CP); (BP)
| |
Collapse
|
32
|
Ejsmond MJ, Phillips KP, Babik W, Radwan J. The role of MHC supertypes in promoting trans-species polymorphism remains an open question. Nat Commun 2018; 9:4362. [PMID: 30341302 PMCID: PMC6195607 DOI: 10.1038/s41467-018-06821-x] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2018] [Accepted: 09/27/2018] [Indexed: 11/09/2022] Open
Affiliation(s)
- Maciej J Ejsmond
- Institute of Environmental Sciences, Jagiellonian University, ul. Gronostajowa 7, 30-387, Kraków, Poland.,Centre for Ecology and Evolution in Microbial model Systems - EEMiS, Linnaeus University, 39182, Kalmar, Sweden
| | - Karl P Phillips
- School of Biological, Earth & Environmental Sciences, University College Cork, Cork, T23 N73K, Ireland.,Marine Institute, Furnace, Newport, Co. Mayo, F28 PF65, Ireland
| | - Wiesław Babik
- Institute of Environmental Sciences, Jagiellonian University, ul. Gronostajowa 7, 30-387, Kraków, Poland
| | - Jacek Radwan
- Evolutionary Biology Group, Institute of Environmental Biology, Adam Mickiewicz University, ul. Umultowska 89, Poznan, Poland.
| |
Collapse
|
33
|
Fijarczyk A, Dudek K, Niedzicka M, Babik W. Balancing selection and introgression of newt immune-response genes. Proc Biol Sci 2018; 285:20180819. [PMID: 30111606 PMCID: PMC6111169 DOI: 10.1098/rspb.2018.0819] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2018] [Accepted: 07/18/2018] [Indexed: 12/20/2022] Open
Abstract
The importance of interspecific introgression as a source of adaptive variation is increasingly recognized. Theory predicts that beneficial genetic variants cross species boundaries easily even when interspecific hybridization is rare and gene flow is strongly constrained throughout the genome. However, it remains unclear whether certain classes of genes are particularly prone to adaptive introgression. Genes affected by balancing selection (BS) may constitute such a class, because forms of BS that favour novel, initially rare alleles, should facilitate introgression. We tested this hypothesis in hybridizing newts by comparing 13 genes with signatures of BS, in particular an excess of common non-synonymous polymorphisms, to the genomic background (154 genes). Parapatric hybridizing taxa were less differentiated in BS candidate genes than more closely related allopatric lineages, while the opposite was observed in the control genes. Coalescent and forward simulations that explored neutral and BS scenarios under isolation and migration showed that processes other than differential gene flow are unlikely to account for this pattern. We conclude that BS, probably involving a form of novel allele advantage, promotes introgression. This mechanism may be a source of adaptively relevant variation in hybridizing species over prolonged periods.
Collapse
Affiliation(s)
- Anna Fijarczyk
- Institute of Environmental Sciences, Jagiellonian University, Gronostajowa 7, 30-387 Kraków, Poland
- Institut de Biologie Intégrative et des Systèmes, Département de Biologie, Université Laval, 1030, Avenue de la Médecine, Québec, Canada G1V 0A6
| | - Katarzyna Dudek
- Institute of Environmental Sciences, Jagiellonian University, Gronostajowa 7, 30-387 Kraków, Poland
| | - Marta Niedzicka
- Institute of Environmental Sciences, Jagiellonian University, Gronostajowa 7, 30-387 Kraków, Poland
| | - Wiesław Babik
- Institute of Environmental Sciences, Jagiellonian University, Gronostajowa 7, 30-387 Kraków, Poland
| |
Collapse
|
34
|
Phillips KP, Cable J, Mohammed RS, Herdegen-Radwan M, Raubic J, Przesmycka KJ, van Oosterhout C, Radwan J. Immunogenetic novelty confers a selective advantage in host-pathogen coevolution. Proc Natl Acad Sci U S A 2018; 115:1552-1557. [PMID: 29339521 PMCID: PMC5816137 DOI: 10.1073/pnas.1708597115] [Citation(s) in RCA: 70] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Abstract
The major histocompatibility complex (MHC) is crucial to the adaptive immune response of vertebrates and is among the most polymorphic gene families known. Its high diversity is usually attributed to selection imposed by fast-evolving pathogens. Pathogens are thought to evolve to escape recognition by common immune alleles, and, hence, novel MHC alleles, introduced through mutation, recombination, or gene flow, are predicted to give hosts superior resistance. Although this theoretical prediction underpins host-pathogen "Red Queen" coevolution, it has not been demonstrated in the context of natural MHC diversity. Here, we experimentally tested whether novel MHC variants (both alleles and functional "supertypes") increased resistance of guppies (Poecilia reticulata) to a common ectoparasite (Gyrodactylus turnbulli). We used exposure-controlled infection trials with wild-sourced parasites, and Gyrodactylus-naïve host fish that were F2 descendants of crossed wild populations. Hosts carrying MHC variants (alleles or supertypes) that were new to a given parasite population experienced a 35-37% reduction in infection intensity, but the number of MHC variants carried by an individual, analogous to heterozygosity in single-locus systems, was not a significant predictor. Our results provide direct evidence of novel MHC variant advantage, confirming a fundamental mechanism underpinning the exceptional polymorphism of this gene family and highlighting the role of immunogenetic novelty in host-pathogen coevolution.
Collapse
Affiliation(s)
- Karl P Phillips
- Evolutionary Biology Group, Faculty of Biology, Adam Mickiewicz University, 60-614 Poznań, Poland
- School of Biological, Earth & Environmental Sciences, University College Cork, Cork, Ireland
| | - Joanne Cable
- School of Biosciences, Cardiff University, CF10 3AX Cardiff, United Kingdom
| | - Ryan S Mohammed
- Department of Life Sciences, The University of the West Indies, St. Augustine, Trinidad and Tobago
| | - Magdalena Herdegen-Radwan
- Department of Behavioural Ecology, Institute of Environmental Biology, Faculty of Biology, Adam Mickiewicz University, 61-614 Poznań, Poland
| | - Jarosław Raubic
- Evolutionary Biology Group, Faculty of Biology, Adam Mickiewicz University, 60-614 Poznań, Poland
| | - Karolina J Przesmycka
- Evolutionary Biology Group, Faculty of Biology, Adam Mickiewicz University, 60-614 Poznań, Poland
| | - Cock van Oosterhout
- School of Environmental Sciences, University of East Anglia, Norwich Research Park, NR4 7TJ Norwich, United Kingdom
| | - Jacek Radwan
- Evolutionary Biology Group, Faculty of Biology, Adam Mickiewicz University, 60-614 Poznań, Poland;
| |
Collapse
|
35
|
|
36
|
Bento G, Routtu J, Fields PD, Bourgeois Y, Du Pasquier L, Ebert D. The genetic basis of resistance and matching-allele interactions of a host-parasite system: The Daphnia magna-Pasteuria ramosa model. PLoS Genet 2017; 13:e1006596. [PMID: 28222092 PMCID: PMC5340410 DOI: 10.1371/journal.pgen.1006596] [Citation(s) in RCA: 39] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2016] [Revised: 03/07/2017] [Accepted: 01/20/2017] [Indexed: 11/18/2022] Open
Abstract
Negative frequency-dependent selection (NFDS) is an evolutionary mechanism suggested to govern host-parasite coevolution and the maintenance of genetic diversity at host resistance loci, such as the vertebrate MHC and R-genes in plants. Matching-allele interactions of hosts and parasites that prevent the emergence of host and parasite genotypes that are universally resistant and infective are a genetic mechanism predicted to underpin NFDS. The underlying genetics of matching-allele interactions are unknown even in host-parasite systems with empirical support for coevolution by NFDS, as is the case for the planktonic crustacean Daphnia magna and the bacterial pathogen Pasteuria ramosa. We fine-map one locus associated with D. magna resistance to P. ramosa and genetically characterize two haplotypes of the Pasteuria resistance (PR-) locus using de novo genome and transcriptome sequencing. Sequence comparison of PR-locus haplotypes finds dramatic structural polymorphisms between PR-locus haplotypes including a large portion of each haplotype being composed of non-homologous sequences resulting in haplotypes differing in size by 66 kb. The high divergence of PR-locus haplotypes suggest a history of multiple, diverse and repeated instances of structural mutation events and restricted recombination. Annotation of the haplotypes reveals striking differences in gene content. In particular, a group of glycosyltransferase genes that is present in the susceptible but absent in the resistant haplotype. Moreover, in natural populations, we find that the PR-locus polymorphism is associated with variation in resistance to different P. ramosa genotypes, pointing to the PR-locus polymorphism as being responsible for the matching-allele interactions that have been previously described for this system. Our results conclusively identify a genetic basis for the matching-allele interaction observed in a coevolving host-parasite system and provide a first insight into its molecular basis.
Collapse
Affiliation(s)
- Gilberto Bento
- Basel University, Zoological Institute, Vesalgasse 1, Basel, Switzerland
- * E-mail:
| | - Jarkko Routtu
- Basel University, Zoological Institute, Vesalgasse 1, Basel, Switzerland
| | - Peter D. Fields
- Basel University, Zoological Institute, Vesalgasse 1, Basel, Switzerland
| | - Yann Bourgeois
- Basel University, Zoological Institute, Vesalgasse 1, Basel, Switzerland
| | - Louis Du Pasquier
- Basel University, Zoological Institute, Vesalgasse 1, Basel, Switzerland
| | - Dieter Ebert
- Basel University, Zoological Institute, Vesalgasse 1, Basel, Switzerland
| |
Collapse
|
37
|
Positive selection in cathelicidin host defense peptides: adaptation to exogenous pathogens or endogenous receptors? Heredity (Edinb) 2016; 118:453-465. [PMID: 27925615 PMCID: PMC5564380 DOI: 10.1038/hdy.2016.117] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2016] [Revised: 09/28/2016] [Accepted: 10/03/2016] [Indexed: 11/24/2022] Open
Abstract
The cause of adaptive protein evolution includes internal (for example, co-evolution of ligand-receptor pairs) and external (for example, adaptation to different ecological niches) mechanisms. Host defense peptides (HDPs) are a class of vertebrate-specific cationic antimicrobial peptides evolving under positive selection. Besides their antibiotic activity, HDPs also exert an effect on multiple host immune cells, thus providing an ideal model to study selective agents driving their evolution. On the basis of a combination of computational and experimental approaches, we studied the evolution of LL-37-type HDPs in mammals, the mature peptide of cathelicidin CAP18 (herein termed CAP18-MP) and investigated the driving force behind the evolution. Using codon-substitution maximum likelihood models, we analyzed CAP18-MPs in 40 species belonging to nine mammalian Orders and identified four positively selected sites (PSSs) that all are located on two terminal unordered regions of CAP18-MPs. Grafting the two positively selected regions of human or whale CAP18-MP to the α-helical scaffold of a rabbit homolog (substituting its corresponding parts) led to no alterations in antibacterial activity, spectrum and action mode. Likewise, further deletion of the two terminal regions did not alter its functional features. Evolutionary conservation analysis of mammalian FPR2, a receptor known to interact with the C-terminal positively selected region of LL-37, revealed high evolutionary variability in its ligand-binding extracellular loop domains, matching sequence diversity of the unordered regions in CAP18-MPs. This is the first report describing that the signature of positive selection of cathelicidins is not associated with their direct bactericidal activity, but rather with the evolutionary variability of their endogenous receptors.
Collapse
|
38
|
Schwensow N, Mazzoni CJ, Marmesat E, Fickel J, Peacock D, Kovaliski J, Sinclair R, Cassey P, Cooke B, Sommer S. High adaptive variability and virus-driven selection on major histocompatibility complex (MHC) genes in invasive wild rabbits in Australia. Biol Invasions 2016. [DOI: 10.1007/s10530-016-1329-5] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
39
|
Migalska M, Sebastian A, Konczal M, Kotlík P, Radwan J. De novo transcriptome assembly facilitates characterisation of fast-evolving gene families, MHC class I in the bank vole (Myodes glareolus). Heredity (Edinb) 2016; 118:348-357. [PMID: 27782121 DOI: 10.1038/hdy.2016.105] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2016] [Accepted: 09/20/2016] [Indexed: 02/07/2023] Open
Abstract
The major histocompatibility complex (MHC) plays a central role in the adaptive immune response and is the most polymorphic gene family in vertebrates. Although high-throughput sequencing has increasingly been used for genotyping families of co-amplifying MHC genes, its potential to facilitate early steps in the characterisation of MHC variation in nonmodel organism has not been fully explored. In this study we evaluated the usefulness of de novo transcriptome assembly in characterisation of MHC sequence diversity. We found that although de novo transcriptome assembly of MHC I genes does not reconstruct sequences of individual alleles, it does allow the identification of conserved regions for PCR primer design. Using the newly designed primers, we characterised MHC I sequences in the bank vole. Phylogenetic analysis of the partial MHC I coding sequence (2-4 exons) of the bank vole revealed a lack of orthology to MHC I of other Cricetidae, consistent with the high gene turnover of this region. The diversity of expressed alleles was characterised using ultra-deep sequencing of the third exon that codes for the peptide-binding region of the MHC molecule. High allelic diversity was demonstrated, with 72 alleles found in 29 individuals. Interindividual variation in the number of expressed loci was found, with the number of alleles per individual ranging from 5 to 14. Strong signatures of positive selection were found for 8 amino acid sites, most of which are inferred to bind antigens in human MHC, indicating conservation of structure despite rapid sequence evolution.
Collapse
Affiliation(s)
- M Migalska
- Evolutionary Biology Group, Faculty of Biology, Adam Mickiewicz University, Poznan, Poland
| | - A Sebastian
- Evolutionary Biology Group, Faculty of Biology, Adam Mickiewicz University, Poznan, Poland
| | - M Konczal
- Evolutionary Biology Group, Faculty of Biology, Adam Mickiewicz University, Poznan, Poland
| | - P Kotlík
- Laboratory of Molecular Ecology, Institute of Animal Physiology and Genetics, The Czech Academy of Sciences, Liběchov, Czech Republic
| | - J Radwan
- Evolutionary Biology Group, Faculty of Biology, Adam Mickiewicz University, Poznan, Poland
| |
Collapse
|
40
|
Bateson ZW, Hammerly SC, Johnson JA, Morrow ME, Whittingham LA, Dunn PO. Specific alleles at immune genes, rather than genome‐wide heterozygosity, are related to immunity and survival in the critically endangered Attwater's prairie‐chicken. Mol Ecol 2016; 25:4730-44. [DOI: 10.1111/mec.13793] [Citation(s) in RCA: 50] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2016] [Revised: 07/17/2016] [Accepted: 07/19/2016] [Indexed: 12/20/2022]
Affiliation(s)
- Zachary W. Bateson
- Behavioral and Molecular Ecology Group Department of Biological Sciences University of Wisconsin‐Milwaukee PO Box 413 Milwaukee WI 53201 USA
| | - Susan C. Hammerly
- Department of Biological Sciences Institute of Applied Sciences University of North Texas 1155 Union Circle Denton TX 76203 USA
| | - Jeff A. Johnson
- Department of Biological Sciences Institute of Applied Sciences University of North Texas 1155 Union Circle Denton TX 76203 USA
| | - Michael E. Morrow
- Attwater Prairie Chicken National Wildlife Refuge PO Box 519 Eagle Lake TX 77434 USA
| | - Linda A. Whittingham
- Behavioral and Molecular Ecology Group Department of Biological Sciences University of Wisconsin‐Milwaukee PO Box 413 Milwaukee WI 53201 USA
| | - Peter O. Dunn
- Behavioral and Molecular Ecology Group Department of Biological Sciences University of Wisconsin‐Milwaukee PO Box 413 Milwaukee WI 53201 USA
| |
Collapse
|
41
|
Rabajante JF, Tubay JM, Ito H, Uehara T, Kakishima S, Morita S, Yoshimura J, Ebert D. Host-parasite Red Queen dynamics with phase-locked rare genotypes. SCIENCE ADVANCES 2016; 2:e1501548. [PMID: 26973878 PMCID: PMC4783124 DOI: 10.1126/sciadv.1501548] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/31/2015] [Accepted: 01/08/2016] [Indexed: 06/05/2023]
Abstract
Interactions between hosts and parasites have been hypothesized to cause winnerless coevolution, called Red Queen dynamics. The canonical Red Queen dynamics assume that all interacting genotypes of hosts and parasites undergo cyclic changes in abundance through negative frequency-dependent selection, which means that any genotype could become frequent at some stage. However, this prediction cannot explain why many rare genotypes stay rare in natural host-parasite systems. To investigate this, we build a mathematical model involving multihost and multiparasite genotypes. In a deterministic and controlled environment, Red Queen dynamics occur between two genotypes undergoing cyclic dominance changes, whereas the rest of the genotypes remain subordinate for long periods of time in phase-locked synchronized dynamics with low amplitude. However, introduction of stochastic noise in the model might allow the subordinate cyclic host and parasite types to replace dominant cyclic types as new players in the Red Queen dynamics. The factors that influence such evolutionary switching are interhost competition, specificity of parasitism, and degree of stochastic noise. Our model can explain, for the first time, the persistence of rare, hardly cycling genotypes in populations (for example, marine microbial communities) undergoing host-parasite coevolution.
Collapse
Affiliation(s)
- Jomar F. Rabajante
- Graduate School of Science and Technology, Shizuoka University, 3-5-1 Johoku, Naka-ku, Hamamatsu 432-8561, Japan
- Mathematics Division, Institute of Mathematical Sciences and Physics, University of the Philippines Los Baños, College, Laguna 4031, Philippines
| | - Jerrold M. Tubay
- Graduate School of Science and Technology, Shizuoka University, 3-5-1 Johoku, Naka-ku, Hamamatsu 432-8561, Japan
- Mathematics Division, Institute of Mathematical Sciences and Physics, University of the Philippines Los Baños, College, Laguna 4031, Philippines
| | - Hiromu Ito
- Graduate School of Science and Technology, Shizuoka University, 3-5-1 Johoku, Naka-ku, Hamamatsu 432-8561, Japan
| | - Takashi Uehara
- Graduate School of Science and Technology, Shizuoka University, 3-5-1 Johoku, Naka-ku, Hamamatsu 432-8561, Japan
- Department of Preschool Education, Nagoya College, Toyoake, Aichi 470-1193, Japan
| | - Satoshi Kakishima
- Graduate School of Science and Technology, Shizuoka University, 3-5-1 Johoku, Naka-ku, Hamamatsu 432-8561, Japan
| | - Satoru Morita
- Department of Mathematical and Systems Engineering, Shizuoka University, 3-5-1 Johoku, Naka-ku, Hamamatsu 432-8561, Japan
| | - Jin Yoshimura
- Graduate School of Science and Technology, Shizuoka University, 3-5-1 Johoku, Naka-ku, Hamamatsu 432-8561, Japan
- Department of Mathematical and Systems Engineering, Shizuoka University, 3-5-1 Johoku, Naka-ku, Hamamatsu 432-8561, Japan
- Marine Biosystems Research Center, Chiba University, Uchiura, Kamogawa, Chiba 299-5502, Japan
- Department of Environmental and Forest Biology, State University of New York College of Environmental Science and Forestry, Syracuse, NY 13210, USA
| | - Dieter Ebert
- Zoological Institute, University of Basel, Vesalgasse 1, Basel 4051, Switzerland
| |
Collapse
|