Mhc supertypes confer both qualitative and quantitative resistance to avian malaria infections in a wild bird population

MHC Class II Supertypes Affect Survival and Lifetime Reproductive Success in a Migratory Songbird

The major histocompatibility complex (MHC) plays a critical role in the immune response against pathogens. Its high polymorphism is thought to be mainly the consequence of host-pathogen co-evolution, but elucidating the mechanism(s) driving MHC evolution remains challenging for natural populations. We investigated the diversity of MHC class II genes in a wild population of pied flycatchers Ficedula hypoleuca and tested its associations with two key components of individual fitness: lifetime reproductive success and survival. Among 180 breeding adults in our study population, we found 182 unique MHC class II exon 2 alleles. The alleles showed a strong signal of positive selection and grouped into nine functional supertypes based on physicochemical properties at the inferred antigen-binding sites. Three supertypes were found in > 98% of the sampled individuals, indicating that they are nearly fixed in the population. We found no rare supertypes in the population, as all supertypes were present in > 70% of individuals. Three supertypes were related to different components of individual fitness: two were associated with lower offspring production over time, while the third was positively associated with survival. Overall, the substantial allelic and functional diversity and the relationship between specific supertypes and fitness are in accordance with the notion that balancing selection maintains MHC class II diversity in the study population, possibly with fluctuating selection as the underlying mechanism. The absence of rare supertypes in the population suggests that the balancing selection is not driven by rare-allele advantage.

Genotype-environment associations reveal genes potentially linked to avian malaria infection in populations of an endemic island bird

Patterns of pathogen prevalence are, at least partially, the result of coevolutionary host-pathogen interactions. Thus, exploring the distribution of host genetic variation in relation to infection by a pathogen within and across populations can provide important insights into mechanisms of host defence and adaptation. Here, we use a landscape genomics approach (Bayenv) in conjunction with genome-wide data (ddRADseq) to test for associations between avian malaria (Plasmodium) prevalence and host genetic variation across 13 populations of the island endemic Berthelot's pipit (Anthus berthelotii). Considerable and consistent spatial heterogeneity in malaria prevalence was observed among populations over a period of 15 years. The prevalence of malaria infection was also strongly positively correlated with pox (Avipoxvirus) prevalence. Multiple host loci showed significant associations with malaria prevalence after controlling for genome-wide neutral genetic structure. These sites were located near to or within genes linked to metabolism, stress response, transcriptional regulation, complement activity and the inflammatory response, many previously implicated in vertebrate responses to malarial infection. Our findings identify diverse genes - not just limited to the immune system - that may be involved in host protection against malaria and suggest that spatially variable pathogen pressure may be an important evolutionary driver of genetic divergence among wild animal populations, such as Berthelot's pipit. Furthermore, our data indicate that spatio-temporal variation in multiple different pathogens (e.g. malaria and pox in this case) may have to be studied together to develop a more holistic understanding of host pathogen-mediated evolution.

Limited associations between MHC diversity and reproductive success in a bird species with biparental care

The selective pressure from pathogens on individuals can have direct consequences on reproduction. Genes from the major histocompatibility complex (MHC) are central to the vertebrate adaptive immune system and pathogen resistance. In species with biparental care, each sex has distinct reproductive roles and levels of investment, and due to a trade-off with immunity, one can expect different selective regimes acting upon the MHC of each parent. Here, we addressed whether couples combine each other's variation at MHC loci to increase their breeding success. Specifically, we used a 23-year dataset from a barn owl population (Tyto alba) to understand how MHC class Iα and IIβ functional divergence and supertypes of each parent were associated with clutch size and fledging success. We did not detect associations between MHC diversity and supertypes with the clutch size or with the fledging success. In addition, to understand the relative contribution from the MHC of the genetic parents and the social parents, we analyzed the fledging success using only a cross-fostered dataset. We found several associations of weak-to-moderate effect sizes between the father's MHC and fledging success: (i) lower MHC-Iα divergence in the genetic father increases fledging success, which might improve paternal care during incubation, and (ii) one and two MHC-IIβ DAB2 supertypes in the social father decrease and increase, respectively, fledging success, which may affect the paternal care after hatching. Furthermore, fledging success increased when both parents did not carry MHC-IIβ DAB1 supertype 2, which could suggest conditional effects of this supertype. Although our study relied on a substantial dataset, we showed that the associations between MHC diversity and reproductive success remain scarce and of complex interpretation in the barn owl. Moreover, our results highlighted the need to incorporate more than one proxy of reproductive success and several MHC classes to capture more complex associations.

Immunogenetics, sylvatic plague and its vectors: insights from the pathogen reservoir Mastomys natalensis in Tanzania

Yersinia pestis is a historically important vector-borne pathogen causing plague in humans and other mammals. Contemporary zoonotic infections with Y. pestis still occur in sub-Saharan Africa, including Tanzania and Madagascar, but receive relatively little attention. Thus, the role of wildlife reservoirs in maintaining sylvatic plague and spillover risks to humans is largely unknown. The multimammate rodent Mastomys natalensis is the most abundant and widespread rodent in peri-domestic areas in Tanzania, where it plays a major role as a Y. pestis reservoir in endemic foci. Yet, how M. natalensis' immunogenetics contributes to the maintenance of plague has not been investigated to date. Here, we surveyed wild M. natalensis for Y. pestis vectors, i.e., fleas, and tested for the presence of antibodies against Y. pestis using enzyme-linked immunosorbent assays (ELISA) in areas known to be endemic or without previous records of Y. pestis in Tanzania. We characterized the allelic and functional (i.e., supertype) diversity of the major histocompatibility complex (MHC class II) of M. natalensis and investigated links to Y. pestis vectors and infections. We detected antibodies against Y. pestis in rodents inhabiting both endemic areas and areas considered non-endemic. Of the 111 nucleotide MHC alleles, only DRB*016 was associated with an increased infestation with the flea Xenopsylla. Surprisingly, we found no link between MHC alleles or supertypes and antibodies of Y. pestis. Our findings hint, however, at local adaptations towards Y. pestis vectors, an observation that more exhaustive sampling could unwind in the future.

Physiological impacts of chronic and experimental Plasmodium infection on breeding-condition male songbirds

While Plasmodium parasitism is common in songbirds, its impact on avian reproduction is unclear owing to conflicting reports in the existing literature. Particularly understudied is the impact of phase of infection on variation in host reproductive physiology in wild, breeding-condition birds. However, assessing the full impact of Plasmodium on reproductive success in the wild can be difficult because individuals experiencing severe effects of parasitism may not enter the breeding population and may be less likely to be captured during field studies. To address these factors, we quantified metrics of health and reproductive physiology in wild-caught, breeding-condition male dark-eyed juncos (Junco hyemalis hyemalis) before and after experimental Plasmodium inoculation in a captive setting. Metrics of health and reproductive physiology included activity rate, hematocrit, scaled body mass, testosterone and sperm production. Individuals already infected at capture (i.e., chronically infected) had higher levels of hematocrit than males without chronic infections. Experimentally infected males showed a larger reduction in hematocrit and activity rate as compared to controls. However, chronic infection status did not influence the extent of metric decline. Testosterone production did not vary by treatment and most birds produced sperm following inoculation. Broadly, our results suggest that male juncos exposed to Plasmodium during the breeding season likely experience declines in general health, but Plasmodium infections do not negatively impact reproductive physiology. We conclude that physiological tradeoffs in males may favor maintenance of reproductive function despite infection.

Evolution of MHC IIB Diversity Across Cichlid Fish Radiations

The genes of the major histocompatibility complex (MHC) are among the most polymorphic genes in vertebrates and crucial for their adaptive immune response. These genes frequently show inconsistencies between allelic genealogies and species phylogenies. This phenomenon is thought to be the result of parasite-mediated balancing selection maintaining ancient alleles through speciation events (trans-species polymorphism [TSP]). However, allele similarities may also arise from postspeciation mechanisms, such as convergence or introgression. Here, we investigated the evolution of MHC class IIB diversity in the cichlid fish radiations across Africa and the Neotropics by a comprehensive review of available MHC IIB DNA sequence information. We explored what mechanism explains the MHC allele similarities found among cichlid radiations. Our results showed extensive allele similarity among cichlid fish across continents, likely due to TSP. Functionality at MHC was also shared among species of the different continents. The maintenance of MHC alleles for long evolutionary times and their shared functionality may imply that certain MHC variants are essential in immune adaptation, even in species that diverged millions of years ago and occupy different environments.

Do female bluethroats without extra-pair offspring have more MHC-compatible social mates?

Abstract

Genes of the major histocompatibility complex (MHC) are crucial for adaptive immunity in jawed vertebrates, and theory predicts that there should be mate choice for optimizing MHC constitution in the offspring. In a previous study, we demonstrated a non-random female choice of extra-pair males in the bluethroat (Luscinia svecica), yielding offspring that was closer to an intermediate MHC class II (MHCII) allele count than their within-pair halfsiblings. The present study tests whether social pairs with only within-pair young (WPY) in their brood, in the same study population, had a combined MHC-constitution closer to a presumed intermediate optimum, than social pairs with extra-pair young (EPY), with a corresponding pattern in their offspring. As expected, we found that WPY from pure WPY-broods were more MHC-optimal than WPY from mixed broods, but only in broods of young (second year) males. Correspondingly, there was a tendency for social pairs with only WPY in their brood to be more MHC-compatible than social pairs with EPY in their brood, when the male was young. Older bluethroat males have considerably larger testes than young males, and their higher sperm competitiveness could help them secure paternity in their own brood, also when they are not MHC-compatible. In other words, in the sexual conflict over paternity, females may be more likely to realise their preference for a MHC-compatible mate when paired to a young male. As a possible fitness indicator, immune responsiveness to an injected antigen (PHA) was elevated for offspring closer to “the golden mean” in MHCII allele count.

Significance statement

This study contributes to our understanding of MHC-based mate choice in extra-pair mating systems, by showing that female bluethroats (Luscinia svecica) with an MHCII-compatible social mate tend to have no extra-pair young in their brood, but only when the social male is young. This elucidates a possible sexual conflict, in which older social males are able to override female preferences and prevent other males from gaining paternity in their brood through higher sperm production. Studying systems in which extra-pair paternity occurs offers an insight into the genetic benefits of mate choice, as extra-pair males, in contrast to social males, generally contribute only sperm. Further, the strict and thorough genotyping scheme applied in this study enabled us to demonstrate a preference for “the golden mean” in MHC-diversity in a species with one of the highest MHC class II-diversity known to date.

Immune gene expression in the mosquito vector Culex quinquefasciatus during an avian malaria infection

Plasmodium relictum is the most widespread avian malaria parasite in the world. It is listed as one of the 100 most dangerous invasive species, having been responsible for the extinction of several endemic bird species, and the near-demise of several others. Here we present the first transcriptomic study focused on the effect of P. relictum on the immune system of its vector (the mosquito Culex quinquefasciatus) at different times post-infection. We show that over 50% of immune genes identified as being part of the Toll pathway and 30%-40% of the immune genes identified within the Imd pathway are overexpressed during the critical period spanning the parasite's oocyst and sporozoite formation (8-12 days), revealing the crucial role played by both these pathways in this natural mosquito-Plasmodium combination. Comparison of infected mosquitoes with their uninfected counterparts also revealed some unexpected immune RNA expression patterns earlier and later in the infection: significant differences in expression of several immune effectors were observed as early as 30 min after ingestion of the infected blood meal. In addition, in the later stages of the infection (towards the end of the mosquito lifespan), we observed an unexpected increase in immune investment in uninfected, but not in infected, mosquitoes. In conclusion, our work extends the comparative transcriptomic analyses of malaria-infected mosquitoes beyond human and rodent parasites and provides insights into the degree of conservation of immune pathways and into the selective pressures exerted by Plasmodium parasites on their vectors.

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

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

Diversifying selection and climatic effects on major histocompatibility complex class II gene diversity in the greater horseshoe bat

Heterogeneous pathogenic stress can shape major histocompatibility complex (MHC) diversity by influencing the functional plasticity of the immune response. Therefore, MHC diversity could reflect environmental stress, demonstrating its importance in uncovering the mechanisms of adaptive genetic variation. In this study, we combined neutral microsatellite loci, an immune-related MHC II-DRB locus, and climatic factors to unravel the mechanisms affecting the diversity and genetic differentiation of MHC genes in the greater horseshoe bat (Rhinolophus ferrumequinum), a species with a wide geographical distribution that has three distinct genetic lineages in China. First, increased genetic differentiation at the MHC locus among populations compared using microsatellites indicated diversifying selection. Second, the genetic differentiation of MHC and microsatellites were significantly correlated, suggesting that demographic processes exist. However, MHC genetic differentiation was significantly correlated with geographical distance among populations, even after controlling for the neutral markers, suggesting a major effect of selection. Third, although the MHC genetic differentiation was larger than that for microsatellites, there was no significant difference in the genetic differentiation between the two markers among genetic lineages, indicating the effect of balancing selection. Fourth, combined with climatic factors, MHC diversity and supertypes showed significant correlations with temperature and precipitation, but not with the phylogeographic structure of R. ferrumequinum, suggesting an effect of local adaptation driven by climate on MHC diversity. Moreover, the number of MHC supertypes varied between populations and lineages, suggesting regional characteristics and support for local adaptation. Taken together, the results of our study provide insights into the adaptive evolutionary driving forces at different geographic scales in R. ferrumequinum. In addition, climate factors may have played a vital role in driving adaptive evolution in this species.

Meta-analysis of major histocompatibility complex (MHC) class IIA reveals polymorphism and positive selection in many vertebrate species

Pathogen-mediated selection and sexual selection are important drivers of evolution. Both processes are known to target genes of the major histocompatibility complex (MHC), a gene family encoding cell-surface proteins that display pathogen peptides to the immune system. The MHC is also a model for understanding processes such as gene duplication and trans-species allele sharing. The class II MHC protein is a heterodimer whose peptide-binding groove is encoded by an MHC-IIA gene and an MHC-IIB gene. However, our literature review found that class II MHC papers on infectious disease or sexual selection included IIA data only 18% and 9% of the time, respectively. To assess whether greater emphasis on MHC-IIA is warranted, we analysed MHC-IIA sequence data from 50 species of vertebrates (fish, amphibians, birds, mammals) to test for polymorphism and positive selection. We found that the number of MHC-IIA alleles within a species was often high, and covaried with sample size and number of MHC-IIA genes assayed. While MHC-IIA variability tended to be lower than that of MHC-IIB, the difference was only ~25%, with ~3 fewer IIA alleles than IIB. Furthermore, the unexpectedly high MHC-IIA variability showed clear signatures of positive selection in most species, and positive selection on MHC-IIA was stronger in fish than in other surveyed vertebrate groups. Our findings underscore that MHC-IIA can be an important target of selection. Future studies should therefore expand the characterization of MHC-IIA at both allelic and genomic scales, and incorporate MHC-IIA into models of fitness consequences of MHC variation.

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

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

Management of avian malaria in populations of high conservation concern

Avian malaria is a vector-borne disease that is caused by Plasmodium parasites. These parasites are transmitted via mosquito bites and can cause sickness or death in a wide variety of birds, including many threatened and endangered species. This Primer first provides contextual background for the avian malaria system including the life cycle, geographic distribution and spread. Then, we focus on recent advances in understanding avian malaria ecology, including how avian malaria can lead to large ecosystem changes and variation in host immune responses to Plasmodium infection. Finally, we review advances in avian malaria management in vulnerable bird populations including genetic modification methods suitable for limiting the effects of this disease in wild populations and the use of sterile insect techniques to reduce vector abundance.

MHCtools - an R package for MHC high-throughput sequencing data: genotyping, haplotype and supertype inference, and downstream genetic analyses in non-model organisms

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.

High MHC diversity confers no advantage for phenotypic quality and reproductive performance in a wild bird

1. Genes of the major histocompatibility complex (MHC) encode antigen binding molecules and are an integral part of the acquired immune response of vertebrates. In general, high individual MHC diversity is expected to increase fitness by broadening the spectrum of pathogens recognized by the immune system, in accordance with the heterozygote advantage mechanism. On the other hand, the optimality hypothesis assumes that individuals with optimal (intermediate), rather than maximum diversity of the MHC will achieve the highest fitness because of inherent costs associated with expressing diverse MHC alleles. 2. Here, we tested for associations between individual diversity of the MHC class I and class II genes (binding antigens of intra- and extra-cellular pathogens, respectively) and a range of fitness-related traits (condition, ornament expression and reproduction) in an urban population of the Eurasian coot Fulica atra. 3. Contrary to our expectation, we found that high within-individual allelic diversity of MHC genes (both class I and II) was associated with poorer condition (lower blood haemoglobin concentrations), weaker expression of the putative ornament (smaller frontal shield), later onset of breeding and smaller clutches. An analysis of functional MHC allele clusters (supertypes) provided further support for negative associations of MHC diversity with phenotypic quality and reproductive performance, but most of these relationships could not be explained by the presence of specific maladaptive supertypes. Finally, we found little empirical support for the optimality hypothesis in the Eurasian coot. 4. Our results suggest that the costs of high MHC diversity outweighed any benefits associated with broad MHC repertoire, which could be driven by depauperate pathogen diversity in an urban landscape. To the best of our knowledge, this is one of the first studies providing consistent evidence for negative associations of MHC diversity with a range of fitness-related traits in a natural avian population.

Interaction between MHC diversity and constitution, gut microbiota and Astrovirus infections in a neotropical bat

Astroviruses (AstVs) infect numerous mammalian species including reservoirs such as bats. Peptides encoded by the genes of the highly polymorphic Major Histocompatibility Complex (MHC) form the first line of host defence against pathogens. Aside from direct involvement in mounting adaptive immune responses, MHC class II genes are hypothesized to regulate gut commensal diversity and shape the production of immune-modulatory substances by microbes, indirectly affecting host susceptibility. Despite initial empirical evidence for the link between host MHC and the microbiota, associations among these factors remain largely unknown. To fill this gap, we examined MHC allelic diversity and constitution, the gut bacterial community and abundance pattern of a wild population of the neotropical bat (Artibeus jamaicensis) challenged by AstV infections. First, we show an age-dependent relationship between the host MHC class II diversity and constitution and the gut microbiota in AstV uninfected bats. Crucially, these associations changed in AstV infected bats. Additionally, we identify changes in abundance of specific bacterial taxa linked to the presence of certain MHC supertypes and AstV infection. We suggest changes in the microbiota to be either a result of AstV infection or the MHC-mediated modulation of microbial communities. The latter could subsequently affect microbe-mediated immunity and resistance against AstV infection. Our results underscore that the reciprocal nature of host immune genetics, gut microbial diversity and pathogen infection requires attention, which is particularly important given its repercussions for disease susceptibility and severity in wild animal populations with a history of zoonotic spillover and frequent human contact.

Divergent and non-parallel evolution of MHC IIB in the Neotropical Midas cichlid species complex

Background

Ecological diversification is the result of divergent natural selection by contrasting habitat characteristics that favours the evolution of distinct phenotypes. This process can happen in sympatry and in allopatry. Habitat-specific parasite communities have the potential to drive diversification among host populations by imposing selective pressures on their host's immune system. In particular, the hyperdiverse genes of the major histocompatibility complex (MHC) are implicated in parasite-mediated host divergence. Here, we studied the extent of divergence at MHC, and discuss how it may have contributed to the Nicaraguan Midas cichlid species complex diversification, one of the most convincing examples of rapid sympatric parallel speciation.

Results

We genotyped the MHC IIB for individuals from six sympatric Midas cichlid assemblages, each containing species that have adapted to exploit similar habitats. We recovered large allelic and functional diversity within the species complex. While most alleles were rare, functional groups of alleles (supertypes) were common, suggesting that they are key to survival and that they were maintained during colonization and subsequent radiations. We identified lake-specific and habitat-specific signatures for both allelic and functional diversity, but no clear pattern of parallel divergence among ecomorphologically similar phenotypes.

Conclusions

Colonization and demographic effects of the fish could have contributed to MHC evolution in the Midas cichlid in conjunction with habitat-specific selective pressures, such as parasites associated to alternative preys or environmental features. Additional ecological data will help evaluating the role of host–parasite interactions in the Midas cichlid radiations and aid in elucidating the potential role of non-parallel features differentiating crater lake species assemblages.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01997-9.

Immunogenetic variation shapes the gut microbiome in a natural vertebrate population

BACKGROUND

The gut microbiome (GM) can influence many biological processes in the host, impacting its health and survival, but the GM can also be influenced by the host's traits. In vertebrates, Major Histocompatibility Complex (MHC) genes play a pivotal role in combatting pathogens and are thought to shape the host's GM. Despite this-and the documented importance of both GM and MHC variation to individual fitness-few studies have investigated the association between the GM and MHC in the wild.

RESULTS

We characterised MHC class I (MHC-I), MHC class II (MHC-II) and GM variation in individuals within a natural population of the Seychelles warbler (Acrocephalus sechellensis). We determined how the diversity and composition of the GM varied with MHC characteristics, in addition to environmental factors and other host traits. Our results show that the presence of specific MHC alleles, but not MHC diversity, influences both the diversity and composition of the GM in this population. MHC-I alleles, rather than MHC-II alleles, had the greatest impact on the GM. GM diversity was negatively associated with the presence of three MHC-I alleles (Ase-ua3, Ase-ua4, Ase-ua5), and one MHC-II allele (Ase-dab4), while changes in GM composition were associated with the presence of four different MHC-I alleles (Ase-ua1, Ase-ua7, Ase-ua10, Ase-ua11). There were no associations between GM diversity and TLR3 genotype, but GM diversity was positively correlated with genome-wide heterozygosity and varied with host age and field period.

CONCLUSIONS

These results suggest that components of the host's immune system play a role in shaping the GM of wild animals. Host genotype-specifically MHC-I and to a lesser degree MHC-II variation-can modulate the GM, although whether this occurs directly, or indirectly through effects on host health, is unclear. Importantly, if immune genes can regulate host health through modulation of the microbiome, then it is plausible that the microbiome could also influence selection on immune genes. As such, host-microbiome coevolution may play a role in maintaining functional immunogenetic variation within natural vertebrate populations. Video abstract.

Adaptive evolution of major histocompatibility complex class I immune genes and disease associations in coastal juvenile sea turtles

Characterizing polymorphism at the major histocompatibility complex (MHC) genes is key to understanding the vertebrate immune response to disease. Despite being globally afflicted by the infectious tumour disease fibropapillomatosis (FP), immunogenetic variation in sea turtles is minimally explored. We sequenced the α 1 peptide-binding region of MHC class I genes (162 bp) from 268 juvenile green (Chelonia mydas) and 88 loggerhead (Caretta caretta) sea turtles in Florida, USA. We recovered extensive variation (116 alleles) and trans-species polymorphism. Supertyping analysis uncovered three functional MHC supertypes corresponding to the three well-supported clades in the phylogeny. We found significant evidence of positive selection at seven amino acid sites in the class I exon. Random forest modelling and risk ratio analysis of Ch. mydas alleles uncovered one allele weakly associated with smooth FP tumour texture, which may be associated with disease outcome. Our study represents the first characterization of MHC class I diversity in Ch. mydas and the largest sample of sea turtles used to date in any study of adaptive genetic variation, revealing tremendous genetic variation and high adaptive potential to viral pathogen threats. The novel associations we identified between MHC diversity and FP outcomes in sea turtles further highlight the importance of evaluating genetic predictors of disease, including MHC and other functional markers.

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

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

Functional immunogenetic variation, rather than local adaptation, predicts ectoparasite infection intensity in a model fish species

Natural host populations differ in their susceptibility to infection by parasites, and these intrapopulation differences are still an incompletely understood component of host‐parasite dynamics. In this study, we used controlled infection experiments with wild‐caught guppies (Poecilia reticulata) and their ectoparasite Gyrodactylus turnbulli to investigate the roles of local adaptation and host genetic composition (immunogenetic and neutral) in explaining differences in susceptibility to infection. We found differences between our four study host populations that were consistent between two parasite source populations, with no indication of local adaptation by either host or parasite at two tested spatial scales. Greater values of host population genetic variability metrics broadly aligned with lower population mean infection intensity, with the best alignments associated with major histocompatibility complex (MHC) “supertypes”. Controlling for intrapopulation differences and potential inbreeding variance, we found a significant negative relationship between individual‐level functional MHC variability and infection: fish carrying more MHC supertypes experienced infections of lower severity, with limited evidence for supertype‐specific effects. We conclude that population‐level differences in host infection susceptibility probably reflect variation in parasite selective pressure and/or host evolutionary potential, underpinned by functional immunogenetic variation.

Spatio-temporal analyses reveal infectious disease-driven selection in a free-ranging ungulate

Infectious diseases play an important role in wildlife population dynamics by altering individual fitness, but detecting disease-driven natural selection in free-ranging populations is difficult due to complex disease-host relationships. Chronic wasting disease (CWD) is a fatal infectious prion disease in cervids for which mutations in a single gene have been mechanistically linked to disease outcomes, providing a rare opportunity to study disease-driven selection in wildlife. In Wyoming, USA, CWD has gradually spread across mule deer (Odocoileus hemionus) populations, producing natural variation in disease history to evaluate selection pressure. We used spatial variation and a novel temporal comparison to investigate the relationship between CWD and a mutation at codon 225 of the mule deer prion protein gene that slows disease progression. We found that individuals with the 'slow' 225F allele were less likely to test positive for CWD, and the 225F allele was more common in herds exposed to CWD longer. We also found that in the past 2 decades, the 225F allele frequency increased more in herds with higher CWD prevalence. This study expanded on previous research by analysing spatio-temporal patterns of individual and herd-based disease data to present multiple lines of evidence for disease-driven selection in free-ranging wildlife.

Balancing selection, genetic drift, and human-mediated introgression interplay to shape MHC (functional) diversity in Mediterranean brown trout

The extraordinary polymorphism of major histocompatibility complex (MHC) genes is considered a paradigm of pathogen-mediated balancing selection, although empirical evidence is still scarce. Furthermore, the relative contribution of balancing selection to shape MHC population structure and diversity, compared to that of neutral forces, as well as its interaction with other evolutionary processes such as hybridization, remains largely unclear. To investigate these issues, we analyzed adaptive (MHC-DAB gene) and neutral (11 microsatellite loci) variation in 156 brown trout (Salmo trutta complex) from six wild populations in central Italy exposed to introgression from domestic hatchery lineages (assessed with the LDH gene). MHC diversity and structuring correlated with those at microsatellites, indicating the substantial role of neutral forces. However, individuals carrying locally rare MHC alleles/supertypes were in better body condition (a proxy of individual fitness/parasite load) regardless of the zygosity status and degree of sequence dissimilarity of MHC, hence supporting balancing selection under rare allele advantage, but not heterozygote advantage or divergent allele advantage. The association between specific MHC supertypes and body condition confirmed in part this finding. Across populations, MHC allelic richness increased with increasing admixture between native and domestic lineages, indicating introgression as a source of MHC variation. Furthermore, introgression across populations appeared more pronounced for MHC than microsatellites, possibly because initially rare MHC variants are expected to introgress more readily under rare allele advantage. Providing evidence for the complex interplay among neutral evolutionary forces, balancing selection, and human-mediated introgression in shaping the pattern of MHC (functional) variation, our findings contribute to a deeper understanding of the evolution of MHC genes in wild populations exposed to anthropogenic disturbance.

Towards a more healthy conservation paradigm: integrating disease and molecular ecology to aid biological conservation†

Parasites, and the diseases they cause, are important from an ecological and evolutionary perspective because they can negatively affect host fitness and can regulate host populations. Consequently, conservation biology has long recognized the vital role that parasites can play in the process of species endangerment and recovery. However, we are only beginning to understand how deeply parasites are embedded in ecological systems, and there is a growing recognition of the important ways in which parasites affect ecosystem structure and function. Thus, there is an urgent need to revisit how parasites are viewed from a conservation perspective and broaden the role that disease ecology plays in conservation-related research and outcomes. This review broadly focusses on the role that disease ecology can play in biological conservation. Our review specifically emphasizes on how the integration of tools and analytical approaches associated with both disease and molecular ecology can be leveraged to aid conservation biology. Our review first concentrates on disease-mediated extinctions and wildlife epidemics. We then focus on elucidating how host–parasite interactions has improved our understanding of the eco-evolutionary dynamics affecting hosts at the individual, population, community and ecosystem scales. We believe that the role of parasites as drivers and indicators of ecosystem health is especially an exciting area of research that has the potential to fundamentally alter our view of parasites and their role in biological conservation. The review concludes with a broad overview of the current and potential applications of modern genomic tools in disease ecology to aid biological conservation.

The strength of selection is consistent across both domains of the MHC class I peptide-binding groove in birds

Background

The Major Histocompatibility Complex (MHC) codes for the key vertebrate immune receptors responsible for pathogen recognition. Foreign antigens are recognized via their compatibility to hyper-variable region of the peptide-binding groove (PBR), which consists of two separate protein domains. Specifically, the PBR of the MHC class I receptors, which recognize intra-cellular pathogens, has two α domains encoded by exon 2 (α₁) and exon 3 (α₂) of the same gene. Most research on avian MHC class I polymorphism has traditionally focused exclusively on exon 3 and comparisons of selection between the two domains have been hampered by the scarcity of molecular data for exon 2. Thus, it is not clear whether the two domains vary in their specificity towards different antigens and whether they are subject to different selective pressure.

Results

Here, we took advantage of rapidly accumulating genomic resources to test for the differences in selection patterns between both MHC class I domains of the peptide-binding groove in birds. For this purpose, we compiled a dataset of MHC class I exon 2 and 3 sequences for 120 avian species from 46 families. Our phylogenetically-robust approach provided strong evidence for highly consistent levels of selection on the α₁ and α₂ domains. There were strong correlations in all selection measures (number of positively/negatively selected residues and dN/dS ratios) between both PBR exons. Similar positive associations were found for the level of amino acid polymorphism across the two domains.

Conclusions

We conclude that the strength of selection and the level of polymorphism are highly consistent between both peptide-binding domains (α₁ and α₂) of the avian MHC class I.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01812-x.

Influence of land use and host species on parasite richness, prevalence and co-infection patterns

Tropical forests are experiencing increasing impacts from a multitude of anthropogenic activities such as logging and conversion to agricultural use. These perturbations are expected to have strong impacts on ecological interactions and on the transmission dynamics of infectious diseases. To date, no clear picture of the effects of deforestation on vector-borne disease transmission has emerged. This is associated with the challenge of studying complex systems where many vertebrate hosts and vectors co-exist. To overcome this problem, we focused on an innately simplified system - a small oceanic island (São Tomé, Gulf of Guinea). We analyzed the impacts of human land-use on host-parasite interactions by sampling the bird community (1735 samples from 30 species) in natural and anthropogenic land use at different elevations, and screened individuals for haemosporidian parasites from three genera (Plasmodium, Haemoproteus, Leucocytozoon). Overall, Plasmodium had the highest richness but the lowest prevalence, while Leucocytozoon diversity was the lowest despite having the highest prevalence. Interestingly, co-infections (i.e. intra-host diversity) involved primarily Leucocytozoon lineages (95%). We also found marked differences between bird species and habitats. Some bird species showed low prevalence but harbored high diversity of parasites, while others showed high prevalence but were infected with fewer lineages. These infection dynamics are most likely driven by host specificity of parasites and intrinsic characteristics of hosts. In addition, Plasmodium was more abundant in disturbed habitats and at lower elevations, while Leucocytozoon was more prevalent in forest areas and at higher elevations. These results likely reflect the ecological requirements of their vectors: mosquitoes and black flies, respectively.

Sex and hatching order modulate the association between MHC-II diversity and fitness in early-life stages of a wild seabird

Genes of the major histocompatibility complex (MHC) play a pivotal role in parasite resistance, and their allelic diversity has been associated with fitness variations in several taxa. However, studies report inconsistencies in the direction of this association, with either positive, quadratic or no association being described. These discrepancies may arise because the fitness costs and benefits of MHC diversity differ among individuals depending on their exposure and immune responses to parasites. Here, we investigated in black-legged kittiwake (Rissa tridactyla) chicks whether associations between MHC class-II diversity and fitness vary with sex and hatching order. MHC-II diversity was positively associated with growth and tick clearance in female chicks, but not in male chicks. Our data also revealed a positive association between MHC-II diversity and survival in second-hatched female chicks (two eggs being the typical clutch size). These findings may result from condition-dependent parasite infections differentially impacting sexes in relation to hatching order. We thus suggest that it may be important to account for individual heterogeneities in traits that potentially exert selective pressures on MHC diversity in order to properly predict MHC-fitness associations.

Association between louse abundance and MHC II supertypes in Galápagos mockingbirds

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.

Towards a more healthy conservation paradigm: integrating disease and molecular ecology to aid biological conservation†

Parasites, and the diseases they cause, are important from an ecological and evolutionary perspective because they can negatively affect host fitness and can regulate host populations. Consequently, conservation biology has long recognized the vital role that parasites can play in the process of species endangerment and recovery. However, we are only beginning to understand how deeply parasites are embedded in ecological systems, and there is a growing recognition of the important ways in which parasites affect ecosystem structure and function. Thus, there is an urgent need to revisit how parasites are viewed from a conservation perspective and broaden the role that disease ecology plays in conservation-related research and outcomes. This review broadly focusses on the role that disease ecology can play in biological conservation. Our review specifically emphasizes on how the integration of tools and analytical approaches associated with both disease and molecular ecology can be leveraged to aid conservation biology. Our review first concentrates on disease mediated extinctions and wildlife epidemics. We then focus on elucidating how host-parasite interactions has improved our understanding of the eco-evolutionary dynamics affecting hosts at the individual, population, community and ecosystem scales. We believe that the role of parasites as drivers and indicators of ecosystem health is especially an exciting area of research that has the potential to fundamentally alter our view of parasites and their role in biological conservation. The review concludes with a broad overview of the current and potential applications of modern genomic tools in disease ecology to aid biological conservation.

Episodic positive diversifying selection on key immune system genes in major avian lineages

The major histocompatibility complex (MHC) of the adaptive immune system and the toll-like receptor (TLR) family of the innate immune system are involved in the detection of foreign invaders, and thus are subject to parasite-driven molecular evolution. Herein, we tested for macroevolutionary signatures of selection in these gene families within and among all three major clades of birds (Paleognathae, Galloanserae, and Neoaves). We characterized evolutionary relationships of representative immune genes (Mhc1 and Tlr2b) and a control gene (ubiquitin, Ubb), using a relatively large and phylogenetically diverse set of species with complete coding sequences (34 orthologous loci for Mhc1, 29 for Tlr2b, and 37 for Ubb). Episodic positive diversifying selection was found in the gene-wide phylogenies of the two immune genes, as well as at specific sites within each gene (8.5% of codon sites in Mhc1 and 2.7% in Tlr2b), but not in the control gene (Ubb). We found 20% of lineages under episodic diversifying selection in Mhc1 versus 9.1% in Tlr2b. For Mhc1, selection was relaxed in the Galloanserae and intensified in the Neoaves relative to the other clades, but no differences were detected among clades in the Tlr2b gene. In summary, we provide evidence of episodic positive diversifying selection in key immune genes and demonstrate differential strengths of selection within Class Aves, with the adaptive gene showing an increased divergence and evolutionary rate over the innate gene, contributing to the growing understanding of vertebrate immune gene evolution.

Lack of statistical power as a major limitation in understanding MHC-mediated immunocompetence in wild vertebrate populations

Disentangling the sources of variation in developing an effective immune response against pathogens is of major interest to immunoecology and evolutionary biology. To date, the link between immunocompetence and genetic variation at the major histocompatibility complex (MHC) has received little attention in wild animals, despite the key role of MHC genes in activating the adaptive immune system. Although several studies point to a link between MHC and immunocompetence, negative findings have also been reported. Such disparate findings suggest that limited statistical power might be affecting studies on this topic, owing to insufficient sample sizes and/or a generally small effect of MHC on the immunocompetence of wild vertebrates. To clarify this issue, we investigated the link between MHC variation and seven immunocompetence proxies in a large sample of barn owls and estimated the effect sizes and statistical power of this and published studies on this topic. We found that MHC poorly explained variation in immunocompetence of barn owls, with small-to-moderate associations between MHC and immunocompetence in owls (effect size: .1 ≥ r ≤ .3) similar to other vertebrates studied to date. Such small-to-moderate effects were largely associated with insufficient power, which was only sufficient (>0.8) to detect moderate-to-large effect sizes (r ≥ .3). Thus, studies linking MHC variation with immunocompetence in wild populations are underpowered to detect MHC effects, which are likely to be of generally small magnitude. Larger sample sizes (>200) will be required to achieve sufficient power in future studies aiming to robustly test for a link between MHC variation and immunocompetence.

Avian MHC Evolution in the Era of Genomics: Phase 1.0

Birds are a wonderfully diverse and accessible clade with an exceptional range of ecologies and behaviors, making the study of the avian major histocompatibility complex (MHC) of great interest. In the last 20 years, particularly with the advent of high-throughput sequencing, the avian MHC has been explored in great depth in several dimensions: its ability to explain ecological patterns in nature, such as mating preferences; its correlation with parasite resistance; and its structural evolution across the avian tree of life. Here, we review the latest pulse of avian MHC studies spurred by high-throughput sequencing. Despite high-throughput approaches to MHC studies, substantial areas remain in need of improvement with regard to our understanding of MHC structure, diversity, and evolution. Recent studies of the avian MHC have nonetheless revealed intriguing connections between MHC structure and life history traits, and highlight the advantages of long-term ecological studies for understanding the patterns of MHC variation in the wild. Given the exceptional diversity of birds, their accessibility, and the ease of sequencing their genomes, studies of avian MHC promise to improve our understanding of the many dimensions and consequences of MHC variation in nature. However, significant improvements in assembling complete MHC regions with long-read sequencing will be required for truly transformative studies.

Elevated Heterozygosity in Adults Relative to Juveniles Provides Evidence of Viability Selection on Eagles and Falcons

Viability selection yields adult populations that are more genetically variable than those of juveniles, producing a positive correlation between heterozygosity and survival. Viability selection could be the result of decreased heterozygosity across many loci in inbred individuals and a subsequent decrease in survivorship resulting from the expression of the deleterious alleles. Alternatively, locus-specific differences in genetic variability between adults and juveniles may be driven by forms of balancing selection, including heterozygote advantage, frequency-dependent selection, or selection across temporal and spatial scales. We use a pooled-sequencing approach to compare genome-wide and locus-specific genetic variability between 74 golden eagle (Aquila chrysaetos), 62 imperial eagle (Aquila heliaca), and 69 prairie falcon (Falco mexicanus) juveniles and adults. Although genome-wide genetic variability is comparable between juvenile and adult golden eagles and prairie falcons, imperial eagle adults are significantly more heterozygous than juveniles. This evidence of viability selection may stem from a relatively smaller imperial eagle effective population size and potentially greater genetic load. We additionally identify ~2000 single-nucleotide polymorphisms across the 3 species with extreme differences in heterozygosity between juveniles and adults. Many of these markers are associated with genes implicated in immune function or olfaction. These loci represent potential targets for studies of how heterozygote advantage, frequency-dependent selection, and selection over spatial and temporal scales influence survivorship in avian species. Overall, our genome-wide data extend previous studies that used allozyme or microsatellite markers and indicate that viability selection may be a more common evolutionary phenomenon than often appreciated.

Is MHC diversity a better marker for conservation than neutral genetic diversity? A case study of two contrasting dolphin populations

Genetic diversity is essential for populations to adapt to changing environments. Measures of genetic diversity are often based on selectively neutral markers, such as microsatellites. Genetic diversity to guide conservation management, however, is better reflected by adaptive markers, including genes of the major histocompatibility complex (MHC). Our aim was to assess MHC and neutral genetic diversity in two contrasting bottlenose dolphin (Tursiops aduncus) populations in Western Australia-one apparently viable population with high reproductive output (Shark Bay) and one with lower reproductive output that was forecast to decline (Bunbury). We assessed genetic variation in the two populations by sequencing the MHC class II DQB, which encompasses the functionally important peptide binding regions (PBR). Neutral genetic diversity was assessed by genotyping twenty-three microsatellite loci. We confirmed that MHC is an adaptive marker in both populations. Overall, the Shark Bay population exhibited greater MHC diversity than the Bunbury population-for example, it displayed greater MHC nucleotide diversity. In contrast, the difference in microsatellite diversity between the two populations was comparatively low. Our findings are consistent with the hypothesis that viable populations typically display greater genetic diversity than less viable populations. The results also suggest that MHC variation is more closely associated with population viability than neutral genetic variation. Although the inferences from our findings are limited, because we only compared two populations, our results add to a growing number of studies that highlight the usefulness of MHC as a potentially suitable genetic marker for animal conservation. The Shark Bay population, which carries greater adaptive genetic diversity than the Bunbury population, is thus likely more robust to natural or human-induced changes to the coastal ecosystem it inhabits.

Immunogenetic response of the bananaquit in the face of malarial parasites

BACKGROUND

In the arms race between hosts and parasites, genes involved in the immune response are targets for natural selection. Toll-Like Receptor (TLR) genes play a role in parasite detection as part of the innate immune system whereas Major Histocompatibility Complex (MHC) genes encode proteins that display antigens as part of the vertebrate adaptive immune system. Thus, both gene families are under selection pressure from pathogens. The bananaquit (Coereba flaveola) is a passerine bird that is a common host of avian malarial parasites (Plasmodium sp. and Haemoproteus sp.). We assessed molecular variation of TLR and MHC genes in a wild population of bananaquits and identified allelic associations with resistance/susceptibility to parasitic infection to address hypotheses of avian immune response to haemosporidian parasites.

RESULTS

We found that allele frequencies are associated with infection status at the immune loci studied. A consistent general trend showed the infected groups possessed more alleles at lower frequencies, and exhibited unique alleles, compared to the uninfected group.

CONCLUSIONS

Our results support the theory of natural selection favoring particular alleles for resistance while maintaining overall genetic diversity in the population, a mechanism which has been demonstrated in some systems in MHC previously but understudied in TLRs.

Major histocompatibility complex class I diversity limits the repertoire of T cell receptors

Major histocompatibility complex (MHC) genes encode proteins that initiate adaptive immune responses through the presentation of foreign antigens to T cells. The high polymorphism found at these genes, thought to be promoted and maintained by pathogen-mediated selection, contrasts with the limited number of MHC loci found in most vertebrates. Although expressing many diverse MHC genes should broaden the range of detectable pathogens, it has been hypothesized to also cause deletion of larger fractions of self-reactive T cells, leading to a detrimental reduction of the T cell receptor (TCR) repertoire. However, a key prediction of this TCR depletion hypothesis, that the TCR repertoire should be inversely related to the individual MHC diversity, has never been tested. Here, using high-throughput sequencing and advanced sequencing error correction, we provide evidence of such an association in a rodent species with high interindividual variation in the number of expressed MHC molecules, the bank vole (Myodes glareolus). Higher individual diversity of MHC class I, but not class II, was associated with smaller TCR repertoires. Our results thus provide partial support for the TCR depletion model, while also highlighting the complex, potentially MHC class-specific mechanisms by which autoreactivity may trade off against evolutionary expansion of the MHC gene family.

Specific MHC class I supertype associated with parasite infection and color morph in a wild lizard population

The major histocompatibility complex (MHC) is a large gene family that plays a central role in the immune system of all jawed vertebrates. Nonavian reptiles are underrepresented within the MHC literature and little is understood regarding the mechanisms maintaining MHC diversity in this vertebrate group. Here, we examined the relative roles of parasite-mediated selection and sexual selection in maintaining MHC class I diversity of a color polymorphic lizard. We discovered evidence for parasite-mediated selection acting via rare-allele advantage or fluctuating selection as ectoparasite load was significantly lower in the presence of a specific MHC supertype (functional clustering of alleles): supertype four. Based on comparisons between ectoparasite prevalence and load, and assessment of the impact of ectoparasite load on host fitness, we suggest that supertype four confers quantitative resistance to ticks or an intracellular tickborne parasite. We found no evidence for MHC-associated mating in terms of pair genetic distance, number of alleles, or specific supertypes. An association was uncovered between supertype four and male throat color morph. However, it is unlikely that male throat coloration acts as a signal of MHC genotype to conspecifics because we found no evidence to suggest that male throat coloration predicts male mating status. Overall, our results suggest that parasite-mediated selection plays a role in maintaining MHC diversity in this population via rare-allele advantage and/or fluctuating selection. Further work is required to determine whether sexual selection also plays a role in maintaining MHC diversity in agamid lizards.

Quantitative host resistance drives the evolution of increased virulence in an emerging pathogen

Emergent infectious diseases can have a devastating impact on host populations. The high selective pressures on both the hosts and the pathogens frequently lead to rapid adaptations not only in pathogen virulence but also host resistance following an initial outbreak. However, it is often unclear whether hosts will evolve to avoid infection-associated fitness costs by preventing the establishment of infection (here referred to as qualitative resistance) or by limiting its deleterious effects through immune functioning (here referred to as quantitative resistance). Equally, the evolutionary repercussions these different resistance mechanisms have for the pathogen are often unknown. Here, we investigate the co-evolutionary dynamics of pathogen virulence and host resistance following the epizootic outbreak of the highly pathogenic bacterium Mycoplasma gallisepticum in North American house finches (Haemorhous mexicanus). Using an evolutionary modelling approach and with a specific emphasis on the evolved resistance trait, we demonstrate that the rapid increase in the frequency of resistant birds following the outbreak is indicative of strong selection pressure to reduce infection-associated mortality. This, in turn, created the ecological conditions that selected for increased bacterial virulence. Our results thus suggest that quantitative host resistance was the key factor underlying the evolutionary interactions in this natural host-pathogen system.

Blood parasites shape extreme major histocompatibility complex diversity in a migratory passerine

Pathogens are one of the main forces driving the evolution and maintenance of the highly polymorphic genes of the vertebrate major histocompatibility complex (MHC). Although MHC proteins are crucial in pathogen recognition, it is still poorly understood how pathogen-mediated selection promotes and maintains MHC diversity, and especially so in host species with highly duplicated MHC genes. Sedge warblers (Acrocephalus schoenobaenus) have highly duplicated MHC genes, and using data from high-throughput MHC genotyping, we were able to investigate to what extent avian malaria parasites explain temporal MHC class I supertype fluctuations in a long-term study population. We investigated infection status and infection intensities of two different strains of Haemoproteus, that is avian malaria parasites that are known to have significant fitness consequences in sedge warblers. We found that prevalence of avian malaria in carriers of specific MHC class I supertypes was a significant predictor of their frequency changes between years. This finding suggests that avian malaria infections partly drive the temporal fluctuations of the MHC class I supertypes. Furthermore, we found that individuals with a large number of different supertypes had higher resistance to avian malaria, but there was no evidence for an optimal MHC class I diversity. Thus, the two studied malaria parasite strains appear to select for a high MHC class I supertype diversity. Such selection may explain the maintenance of the extremely high number of MHC class I gene copies in sedge warblers and possibly also in other passerines where avian malaria is a common disease.

Characterization of Plasmodium relictum, a cosmopolitan agent of avian malaria

Background

Microscopic research has shown that Plasmodium relictum is the most common agent of avian malaria. Recent molecular studies confirmed this conclusion and identified several mtDNA lineages, suggesting the existence of significant intra-species genetic variation or cryptic speciation. Most identified lineages have a broad range of hosts and geographical distribution. Here, a rare new lineage of P. relictum was reported and information about biological characters of different lineages of this pathogen was reviewed, suggesting issues for future research.

Methods

The new lineage pPHCOL01 was detected in Common chiffchaff Phylloscopus collybita, and the parasite was passaged in domestic canaries Serinus canaria. Organs of infected birds were examined using histology and chromogenic in situ hybridization methods. Culex quinquefasciatus mosquitoes, Zebra finch Taeniopygia guttata, Budgerigar Melopsittacus undulatus and European goldfinch Carduelis carduelis were exposed experimentally. Both Bayesian and Maximum Likelihood analyses identified the same phylogenetic relationships among different, closely-related lineages pSGS1, pGRW4, pGRW11, pLZFUS01, pPHCOL01 of P. relictum. Morphology of their blood stages was compared using fixed and stained blood smears, and biological properties of these parasites were reviewed.

Results

Common canary and European goldfinch were susceptible to the parasite pPHCOL01, and had markedly variable individual prepatent periods and light transient parasitaemia. Exo-erythrocytic and sporogonic stages were not seen. The Zebra finch and Budgerigar were resistant. Neither blood stages nor vector stages of all examined P. relictum lineages can be distinguished morphologically.

Conclusion

Within the huge spectrum of vertebrate hosts, mosquito vectors, and ecological conditions, different lineages of P. relictum exhibit indistinguishable, markedly variable morphological forms. Parasites of same lineages often develop differently in different bird species. Even more, the variation of biological properties (parasitaemia dynamics, blood pathology, prepatent period) in different isolates of the same lineage might be greater than the variation in different lineages during development in the same species of birds, indicating negligible taxonomic value of such features. Available lineage information is excellent for parasite diagnostics, but is limited in predictions about relationships in certain host-parasite associations. A combination of experiments, field observations, microscopic and molecular diagnostics is essential for understanding the role of different P. relictum lineages in bird health.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2325-2) contains supplementary material, which is available to authorized users.

A high-density SNP chip for genotyping great tit (Parus major) populations and its application to studying the genetic architecture of exploration behaviour

High-density SNP microarrays ("SNP chips") are a rapid, accurate and efficient method for genotyping several hundred thousand polymorphisms in large numbers of individuals. While SNP chips are routinely used in human genetics and in animal and plant breeding, they are less widely used in evolutionary and ecological research. In this article, we describe the development and application of a high-density Affymetrix Axiom chip with around 500,000 SNPs, designed to perform genomics studies of great tit (Parus major) populations. We demonstrate that the per-SNP genotype error rate is well below 1% and that the chip can also be used to identify structural or copy number variation. The chip is used to explore the genetic architecture of exploration behaviour (EB), a personality trait that has been widely studied in great tits and other species. No SNPs reached genomewide significance, including at DRD4, a candidate gene. However, EB is heritable and appears to have a polygenic architecture. Researchers developing similar SNP chips may note: (i) SNPs previously typed on alternative platforms are more likely to be converted to working assays; (ii) detecting SNPs by more than one pipeline, and in independent data sets, ensures a high proportion of working assays; (iii) allele frequency ascertainment bias is minimized by performing SNP discovery in individuals from multiple populations; and (iv) samples with the lowest call rates tend to also have the greatest genotyping error rates.

Plasmodium Infections in Natural Populations of Anolis sagrei Reflect Tolerance Rather Than Susceptibility

Parasites can represent formidable selection pressures for hosts, but the cost of infection is sometimes difficult to demonstrate in natural populations. While parasite exploitation strategies may, in some instances, actually inflict low costs on their hosts, the response of hosts to infection is also likely to determine whether or not these costs can be detected. Indeed, costs of infection may be obscured if infected individuals in the wild are those that are the most tolerant, rather than the most susceptible, to infection. Here we test this hypothesis in two natural populations of Anolis sagrei, one of the most common anole lizard of the Bahamas. Plasmodium parasites were detected in > 7% of individuals and belonged to two distinct clades: P. mexicanum and P. floriensis. Infected individuals displayed greater body condition than non-infected ones and we found no association between infection status, stamina, and survival to the end of the breeding season. Furthermore, we found no significant difference in the immuno-competence (measured as a response to phytohemagglutinin challenge) of infected versus non-infected individuals. Taken together, our results suggest that the infected individuals that are caught in the wild are those most able to withstand the cost of the infection and that susceptible, infected individuals have been removed from the population (i.e., through disease-induced mortality). This study highlights the need for caution when interpreting estimates of infection costs in natural populations, as costs may appear low either when parasites exploitation strategies truly inflict low costs on their hosts or when those costs are so high that susceptible hosts are removed from the population.

Characteristics of inflammatory reactions during development of liver abscess in hamsters inoculated with Entamoeba nuttalli

Background

Entamoeba nuttalli is an intestinal protozoan with pathogenic potential that can cause amebic liver abscess. It is highly prevalent in wild and captive macaques. Recently, cysts were detected in a caretaker of nonhuman primates in a zoo, indicating that E. nuttalli may be a zoonotic pathogen. Therefore, it is important to evaluate the pathogenicity of E. nuttalli in detail and in comparison with that of E. histolytica.

Methodology/Principal findings

Trophozoites of E. nuttalli GY4 and E. histolytica SAW755 strains were inoculated into liver of hamsters. Expression levels of proinflammatory factors of hamsters and virulence factors from E. histolytica and E. nuttalli were compared between the two parasites. Inoculations with trophozoites of E. nuttalli resulted in an average necrotic area of 24% in liver tissue in 7 days, whereas this area produced by E. histolytica was nearly 50%. Along with the mild liver tissue damage induced by E. nuttalli, expression levels of proinflammatory factors (TNF-α, IL-6 and IL-1β) and amebic virulence protein genes (lectins, cysteine proteases and amoeba pores) in local tissues were lower with E. nuttalli in comparison with E. histolytica. In addition, M2 type macrophages were increased in E. nuttalli-induced amebic liver abscesses in the late stage of disease progression and lysate of E. nuttalli trophozoites induced higher arginase expression than E. histolytica in vitro.

Conclusions/Significance

The results show that differential secretion of amebic virulence proteins during E. nuttalli infection triggered lower levels of secretion of various cytokines and had an impact on polarization of macrophages towards a M1/M2 balance. However, regardless of the degree of macrophage polarization, there is unambiguous evidence of an intense acute inflammatory reaction in liver of hamsters after infection by both Entamoeba species.

Entamoeba nuttalli is the phylogenetically closest protozoan to Entamoeba histolytica and is highly prevalent in macaques. Previous studies have indicated that E. nuttalli is virulent in a hamster model. In this study, we compared the immunopathological basis of formation of liver abscess in hamsters between E. nuttalli and E. histolytica. Mild liver tissue damage developed after intrahepatic injection of trophozoites of E. nuttalli, and lower expression levels of genes for host proinflammatory factors and amebic virulence proteins were detected at the edges of liver abscesses induced by E. nuttalli. In addition, alternatively activated macrophages were increased in E. nuttalli-induced liver abscesses in the late stage of disease progression. The lysate of E. nuttalli trophozoites also induced higher arginase expression than E. histolytica in vitro. Polarization of macrophages is likely to affect the degree of acute inflammatory reactions in liver in an animal model during E. nuttalli infection. Our data reveal new characteristics of abscess formation by E. nuttalli.

Major histocompatibility complex selection dynamics in pathogen-infected túngara frog (Physalaemus pustulosus) populations

Pathogen-driven selection can favour major histocompatibility complex (MHC) alleles that confer immunological resistance to specific diseases. However, strong directional selection should deplete genetic variation necessary for robust immune function in the absence of balancing selection or challenges presented by other pathogens. We examined selection dynamics at one MHC class II (MHC-II) locus across Panamanian populations of the túngara frog, Physalaemus pustulosus, infected by the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). We compared MHC-II diversity in highland túngara frog populations, where amphibian communities have experienced declines owing to Bd, with those in the lowland region that have shown no evidence of decline. Highland region frogs had MHC variants that confer resistance to Bd. Variant fixation appeared to occur by directional selection rather than inbreeding, as overall genetic variation persisted in populations. In Bd-infected lowland sites, however, selective advantage may accrue to individuals with only one Bd-resistance allele, which were more frequent. Environmental conditions in lowlands should be less favourable for Bd infection, which may reduce selection for specific Bd resistance in hosts. Our results suggest that MHC selection dynamics fluctuate in túngara frog populations as a function of the favourability of habitat to pathogen spread and the vulnerability of hosts to infection.

Evolutionary genetics of immunological supertypes reveals two faces of the Red Queen

Red Queen host-parasite co-evolution can drive adaptations of immune genes by positive selection that erodes genetic variation (Red Queen arms race) or results in a balanced polymorphism (Red Queen dynamics) and long-term preservation of genetic variation (trans-species polymorphism). These two Red Queen processes are opposite extremes of the co-evolutionary spectrum. Here we show that both Red Queen processes can operate simultaneously by analysing the major histocompatibility complex (MHC) in guppies (Poecilia reticulata and P. obscura) and swamp guppies (Micropoecilia picta). Sub-functionalisation of MHC alleles into 'supertypes' explains how polymorphisms persist during rapid host-parasite co-evolution. Simulations show the maintenance of supertypes as balanced polymorphisms, consistent with Red Queen dynamics, whereas alleles within supertypes are subject to positive selection in a Red Queen arms race. Building on the divergent allele advantage hypothesis, we show that functional aspects of allelic diversity help to elucidate the evolution of polymorphic genes involved in Red Queen co-evolution.

Genetic wealth, population health: Major histocompatibility complex variation in captive and wild ring-tailed lemurs (Lemur catta)

Across species, diversity at the major histocompatibility complex (MHC) is critical to individual disease resistance and, hence, to population health; however, MHC diversity can be reduced in small, fragmented, or isolated populations. Given the need for comparative studies of functional genetic diversity, we investigated whether MHC diversity differs between populations which are open, that is experiencing gene flow, versus populations which are closed, that is isolated from other populations. Using the endangered ring-tailed lemur (Lemur catta) as a model, we compared two populations under long-term study: a relatively "open," wild population (n = 180) derived from Bezà Mahafaly Special Reserve, Madagascar (2003-2013) and a "closed," captive population (n = 121) derived from the Duke Lemur Center (DLC, 1980-2013) and from the Indianapolis and Cincinnati Zoos (2012). For all animals, we assessed MHC-DRB diversity and, across populations, we compared the number of unique MHC-DRB alleles and their distributions. Wild individuals possessed more MHC-DRB alleles than did captive individuals, and overall, the wild population had more unique MHC-DRB alleles that were more evenly distributed than did the captive population. Despite management efforts to maintain or increase genetic diversity in the DLC population, MHC diversity remained static from 1980 to 2010. Since 2010, however, captive-breeding efforts resulted in the MHC diversity of offspring increasing to a level commensurate with that found in wild individuals. Therefore, loss of genetic diversity in lemurs, owing to small founder populations or reduced gene flow, can be mitigated by managed breeding efforts. Quantifying MHC diversity within individuals and between populations is the necessary first step to identifying potential improvements to captive management and conservation plans.