MHCtools 1.5: Analysis of MHC Sequencing Data in R.

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.

Only rare classical MHC-I alleles are highly expressed in the European house sparrow.

The exceptional polymorphism observed within genes of the major histocompatibility complex (MHC), a core component of the vertebrate immune system, has long fascinated biologists. The highly polymorphic classical MHC class-I (MHC-I) genes are maintained by pathogen-mediated balancing selection (PMBS), as shown by many sites subject to positive selection, while the more monomorphic non-classical MHC-I genes show signatures of purifying selection. In line with PMBS, at any point in time, rare classical MHC alleles are more likely than common classical MHC alleles to confer a selective advantage in host-pathogen interactions. Combining genomic and expression data from the blood of wild house sparrows Passer domesticus, we found that only rare classical MHC-I alleles were highly expressed, while common classical MHC-I alleles were lowly expressed or not expressed. Moreover, highly expressed rare classical MHC-I alleles had more positively selected sites, indicating exposure to stronger PMBS, compared with lowly expressed classical alleles. As predicted, the level of expression was unrelated to allele frequency in the monomorphic non-classical MHC-I alleles. Going beyond previous studies, we offer a fine-scale view of selection on classical MHC-I genes in a wild population by revealing differences in the strength of PMBS according to allele frequency and expression level.

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 belong 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.

Impact of Chronic Exposure to Two Neonicotinoids on Honey Bee Antennal Responses to Flower Volatiles and Pheromonal Compounds.

Volatile compounds provide important olfactory cues for honey bees (Apis mellifera L.), which are essential for their ecology, behavior, and social communication. In the external environment bees locate food sources by the use of floral scents, while inside the hive, pheromones such as the queen mandibular pheromone (QMP) and alarm pheromones serve important functions in regulating colony life and inducing aggressive responses against intruders and parasites. Widely reported alterations of various behaviors in- and outside the hive following exposure to pesticides could therefore be associated with a disturbance of odor sensitivity. In the present study, we tested the effects of neonicotinoid pesticides at field concentrations on the ability of honey bees to perceive volatiles at the very periphery of the olfactory system. Bee colonies were subjected to treatments during the summer with either Imidacloprid or Thiacloprid at sublethal concentrations. Antennal responses to apple (Malus domestica L.) flower volatiles were studied by GC-coupled electro-antennographic detection (GC-EAD), and a range of volatiles, a substitute of the QMP, and the alarm pheromone 2-heptanone were tested by electroantennography (EAG). Short-term and long-term effects of the neonicotinoid treatments were investigated on bees collected in the autumn and again in the following spring. Treatment with Thiacloprid induced changes in antennal responses to specific flower VOCs, with differing short- and long-term effects. In the short term, increased antennal responses were observed for benzyl-alcohol and 1-hexanol, which are common flower volatiles but also constituents of the honey bee sting gland secretions. The treatment with Thiacloprid also affected antennal responses to the QMP and the mandibular alarm pheromone 2-heptanone. In the short term, a faster signal degeneration of the response signal to the positive control citral was recorded in the antennae of bees exposed to Thiacloprid or Imidacloprid. Finally, we observed season-related differences in the antennal responses to multiple VOCs. Altogether, our results suggest that volatile-specific alterations of antennal responses may contribute to explaining several behavioral changes previously observed in neonicotinoid-exposed bees. Treatment effects were generally more prominent in the short term, suggesting that adverse effects of neonicotinoid exposure may not persist across generations.

Genomics of host-pathogen interactions: challenges and opportunities across ecological and spatiotemporal scales.

Evolutionary genomics has recently entered a new era in the study of host-pathogen interactions. A variety of novel genomic techniques has transformed the identification, detection and classification of both hosts and pathogens, allowing a greater resolution that helps decipher their underlying dynamics and provides novel insights into their environmental context. Nevertheless, many challenges to a general understanding of host-pathogen interactions remain, in particular in the synthesis and integration of concepts and findings across a variety of systems and different spatiotemporal and ecological scales. In this perspective we aim to highlight some of the commonalities and complexities across diverse studies of host-pathogen interactions, with a focus on ecological, spatiotemporal variation, and the choice of genomic methods used. We performed a quantitative review of recent literature to investigate links, patterns and potential tradeoffs between the complexity of genomic, ecological and spatiotemporal scales undertaken in individual host-pathogen studies. We found that the majority of studies used whole genome resolution to address their research objectives across a broad range of ecological scales, especially when focusing on the pathogen side of the interaction. Nevertheless, genomic studies conducted in a complex spatiotemporal context are currently rare in the literature. Because processes of host-pathogen interactions can be understood at multiple scales, from molecular-, cellular-, and physiological-scales to the levels of populations and ecosystems, we conclude that a major obstacle for synthesis across diverse host-pathogen systems is that data are collected on widely diverging scales with different degrees of resolution. This disparity not only hampers effective infrastructural organization of the data but also data granularity and accessibility. Comprehensive metadata deposited in association with genomic data in easily accessible databases will allow greater inference across systems in the future, especially when combined with open data standards and practices. The standardization and comparability of such data will facilitate early detection of emerging infectious diseases as well as studies of the impact of anthropogenic stressors, such as climate change, on disease dynamics in humans and wildlife.

Evidence for sexual conflict over major histocompatibility complex diversity in a wild songbird.

Sex differences in parasite load and immune responses are found across a wide range of animals, with females generally having lower parasite loads and stronger immune responses than males. Intrigued by these general patterns, we investigated if there was any sign of sex-specific selection on an essential component of adaptive immunity that is known to affect fitness, the major histocompatibility complex class I (MHC-I) genes, in a 20-year study of great reed warblers. Our analyses on fitness related to MHC-I diversity showed a highly significant interaction between MHC-I diversity and sex, where males with higher, and females with lower, MHC-I diversity were more successful in recruiting offspring. Importantly, mean MHC-I diversity did not differ between males and females, and consequently neither sex reached its MHC-I fitness optimum. Thus, there is an unresolved genetic sexual conflict over MHC-I diversity in great reed warblers. Selection from pathogens is known to maintain MHC diversity, but previous theory ignores that the immune environments are considerably different in males and females. Our results suggest that sexually antagonistic selection is an important, previously neglected, force in the evolution of vertebrate adaptive immunity, and have implications for evolutionary understanding of costs of immune responses and autoimmune diseases.

Sex differences in immune responses: Hormonal effects, antagonistic selection, and evolutionary consequences.

Males and females differ in both parasite load and the strength of immune responses and these effects have been verified in humans and other vertebrates. Sex hormones act as important modulators of immune responses; the male sex hormone testosterone is generally immunosuppressive while the female sex hormone estrogen tends to be immunoenhancing. Different sets of T-helper cells (Th) have important roles in adaptive immunity, e.g. Th1 cells trigger type 1 responses which are primarily cell-mediated, and Th2 cells trigger type 2 responses which are primarily humoral responses. In our review of the literature, we find that estrogen and progesterone enhance type 2 and suppress type 1 responses in females, whereas testosterone suppresses type 2 responses and shows an inconsistent pattern for type 1 responses in males. When we combine these patterns of generally immunosuppressive and immunoenhancing effects of the sex hormones, our results imply that the sex differences in immune responses should be particularly strong in immune functions associated with type 2 responses, and less pronounced with type 1 responses. In general the hormone-mediated sex differences in immune responses may lead to genetic sexual conflicts on immunity. Thus, we propose the novel hypothesis that sexually antagonistic selection may act on immune genes shared by the sexes, and that the strength of this sexually antagonistic selection should be stronger for type 2- as compared with type 1-associated immune genes. Finally, we put the consequences of sex hormone-induced effects on immune responses into behavioral and ecological contexts, considering social mating system, sexual selection, geographical distribution of hosts, and parasite abundance.