The nature of selection on the major histocompatibility complex

Diversity begets diversity: do parasites promote variation in protective symbionts?

Insects commonly possess heritable microbial symbionts that increase their resistance to particular parasites. A diverse community of defensive symbionts may thus provide hosts with effective and specific protection against multiple parasites, although costs might constrain the accumulation of many symbionts. In parallel to the allelic diversity in the MHC complex of the vertebrate immune system, parasite diversity could be the driving force behind symbiont diversity. There is indeed evidence that parasites have the ability to drive frequencies of defensive symbionts in their hosts, and that these symbionts influence parasite communities, but direct evidence that parasite diversity can promote symbiont diversity is still lacking. We provide suggestions to investigate this potential link.

Behavioural and neuronal basis of olfactory imprinting and kin recognition in larval fish

Imprinting is a specific form of long-term memory of a cue acquired during a sensitive phase of development. To ensure that organisms memorize the right cue, the learning process must happen during a specific short time period, mostly soon after hatching, which should end before irrelevant or misleading signals are encountered. A well-known case of olfactory imprinting in the aquatic environment is that of the anadromous Atlantic and Pacific salmon, which prefer the olfactory cues of natal rivers to which they return after migrating several years in the open ocean. Recent research has shown that olfactory imprinting and olfactory guided navigation in the marine realm are far more common than previously assumed. Here, we present evidence for the involvement of olfactory imprinting in the navigation behaviour of coral reef fish, which prefer their home reef odour over that of other reefs. Two main olfactory imprinting processes can be differentiated: (1) imprinting on environmental cues and (2) imprinting on chemical compounds released by kin, which is based on genetic relatedness among conspecifics. While the first process allows for plasticity, so that organisms can imprint on a variety of chemical signals, the latter seems to be restricted to specific genetically determined kin signals. We focus on the second, elucidating the behavioural and neuronal basis of the imprinting process on kin cues using larval zebrafish (Danio rerio) as a model. Our data suggest that the process of imprinting is not confined to the central nervous system but also triggers some changes in the olfactory epithelium.

MHC class IIA polymorphisms and their association with resistance-susceptibility to Streptococcus agalactiae in Nile tilapia, Oreochromis niloticus

The association between major histocompatibility complex (MHC) class IIA polymorphisms and the severity of infection by Streptococcus agalactiae was investigated using 40 susceptible and 40 resistant individuals of Nile tilapia Oreochromis niloticus. Twenty-five alleles were identified from 80 individuals, which belong to 22 major allele types. High polymorphism of mhcIIa gene and at least two loci were discovered in O. niloticus. In peptide-binding region (PBR) and non-PBR, the ratio of nonsynonymous substitution (dN) to synonymous substitution (dS) was 1.294 (>1) and 1.240 (>1), suggesting that the loci are evolving under positive balancing selection. Association analysis showed that the allele, orni-daa*0501, was significantly associated with resistance to S. agalactiae, while the alleles, orni-daa*1101, orni-daa*1301, orni-daa*1401 and orni-daa*1201, were associated with susceptibility to S. agalactiae. To confirm these correlations, another independent challenge experiment was performed in the Huizhou population of the O. niloticus. The frequency distribution showed that the orni-daa*1101 allele was significantly more frequent in the Huizhou-Susceptible group (HZ-SG) than in the Huizhou-Resistant group (HZ-RG) (P < 0.05), which was consistent with the first challenge. However, orni-daa*0501 did not present in HZ-SG and HZ-RG and the distribution frequencies of the orni-daa*1201, orni-daa*1301 and orni-daa*1401 alleles were not significantly more frequent in HZ-SG than in HZ-RG. These results indicate that the orni-daa*1101 allele confers susceptibility to S. agalactia infection. These results suggest that the diversity of exon 2 of mcaIIa alleles could be used to explore the association between disease susceptibility or resistance and the multiformity of mcaIIa and to achieve the molecular-assisted selection of O. niloticus with enhanced disease resistance.

Improving the sustainability of ex situ populations with mate choice

Many breeding programs managed by the Association of Zoos and Aquariums' (AZA) Species Survival Plans® (SSPs) are not meeting goals for population size and genetic diversity due to failure of recommended pairs to breed successfully. According to AZA Population Management Center analyses, as many as 80% of recommended breeding pairs fail to produce young before the next breeding and transfer plan is issued. Determining reasons for failure and ensuring that a specific pairing produces offspring can be challenging. Mate incompatibility, one possible reason for failure, might be addressed by allowing mate choice. Although many SSP® coordinators and breeding managers, who implement breeding recommendations at their institutions, recognize the potential benefits of mate choice, examples and models for presenting and assessing choice are lacking. Here we review examples from birds, rodents, lagomorphs, marsupials, carnivores, fishes, and insects where mate choice has been incorporated. These examples provide strong evidence that free mate choice and mating with preferred partners increase a variety of reproductive success measurements when compared to assigned mate pairings. We suggest innovative housing and breeding arrangements for better incorporating mate choice into the management strategies for species held ex situ. Further, we discuss the fitness consequences and welfare implications of allowing choice. We advocate for a more systematic use of behavioral research in cooperative breeding programs. Behavioral management for mating can yield more successful programs, thus ensuring SSP® genetic and demographic goals are met, while simultaneously improving welfare.

Positive selection in coding regions and motif duplication in regulatory regions of bottlenose dolphin MHC class II genes

The vertebrate immune response is mediated through highly adaptive, quickly evolving cell surface receptors, the major histocompatibility complex (MHC). MHC molecules bind and present a diverse array of pathogenic molecules and trigger a cascade of defenses. Use of MHC variation as a marker for population health has also evolved quickly following advances in sequencing methods. We applied a combination of traditional and next generation sequencing methodology to characterize coding (peptide binding region) and regulatory (proximal promoter) sequence variation in MHC Class II DQA and DQB genes between estuarine and coastal populations of the bottlenose dolphin, Tursiops truncatus, an apex predator whose health status is indicative of anthropogenic impacts on the ecosystem. The coding regions had 10 alleles each at DQA and DQB; the promoters had 6 and 7 alleles at DQA and DQB, respectively with variation within key regulatory motifs. Positive selection was observed for the coding regions of both genes while both coding and promoter regions exhibited geographic differences in allele composition that likely indicates diversifying selection across habitats. Most notable was the discovery of a complete duplication of a 14-bp T-box motif in the DQA promoter. Four class II promoter regions (DQA, DQB, DRA, DRB) were characterized in species from four cetacean families (Delphinidae, Monodontidae, Lipotidae, and Physeteridae) and revealed substantial promoter structural diversity across this order. Peptide binding regions may not be the only source of adaptive potential within cetacean MHC for responding to pathogenic threats. These findings are the first analysis of cetacean MHC regulatory motifs, which may divulge unique immunogenetic strategies among cetaceans and reveal how MHC transcriptional control continues to evolve. The combined MHC regulatory and coding data provide new genetic context for distinct vulnerability profiles between coastal and estuarine populations, which are key concerns for health and risk management.

Signatures of balancing selection in toll-like receptor (TLRs) genes - novel insights from a free-living rodent

Selective pressure from pathogens is considered a key selective force driving the evolution of components of the immune system. Since single components of the immune system may interact with many pathogens, and single pathogens may be recognized by multiple components of the immune system, gaining a better understanding of the mechanisms of parasite-driven selection requires the study of multiple genes and pathogens. Toll-like receptors (TLRs) are a large gene family that code for antigen-presenting components of the innate immune response. In the present paper we characterize polymorphism and signatures of selection in seven TLRs in free-living bank voles Myodes glareolus. We report the first evidence of balancing selection in several TLR genes, supported by positive values of Fu and Li's D* in TLR2 and TLR5, and positive values of Tajima's D in LRR regions within TLR1 and TLR2. We further found significant associations between amino-acid alleles of TLR1 and TLR5 and susceptibility to infection with the blood pathogen Bartonella. Interestingly, selection patterns in TLRs presenting virus-derived motifs (TLR7 and TLR9) differed considerably from those interacting with bacterial PAMPs. In contrast to the highly variable TLRs presenting bacterial motifs, TLR7 and TLR9 had low polymorphism and displayed signatures of directional selection. These findings suggest different functional responses across the TLR gene family and highlight the complexity of parasite-driven selection.

Polymorphic sites preferentially avoid co-evolving residues in MHC class I proteins

Major histocompatibility complex class I (MHC-I) molecules are critical to adaptive immune defence mechanisms in vertebrate species and are encoded by highly polymorphic genes. Polymorphic sites are located close to the ligand-binding groove and entail MHC-I alleles with distinct binding specificities. Some efforts have been made to investigate the relationship between polymorphism and protein stability. However, less is known about the relationship between polymorphism and MHC-I co-evolutionary constraints. Using Direct Coupling Analysis (DCA) we found that co-evolution analysis accurately pinpoints structural contacts, although the protein family is restricted to vertebrates and comprises less than five hundred species, and that the co-evolutionary signal is mainly driven by inter-species changes, and not intra-species polymorphism. Moreover, we show that polymorphic sites in human preferentially avoid co-evolving residues, as well as residues involved in protein stability. These results suggest that sites displaying high polymorphism may have been selected during vertebrates’ evolution to avoid co-evolutionary constraints and thereby maximize their mutability.

Amino acid co-evolution represents cases of simultaneous substitution of amino acids at distinct positions in protein sequences. In the MHC-I protein family, such co-evolution could result from either amino acid changes across species or changes within species due to the high polymorphism of MHC-I molecules. Here we show that signals captured by global methods such as Direct Coupling Analysis (DCA) to estimate co-evolution primarily result from changes across species. Moreover, our results indicate that polymorphic sites in MHC-I molecules tend to be decoupled from co-evolving ones. This could suggest that they have been selected to maximize their mutability, which is known to be functionally important to entail MHC-I molecules with a wide repertoire of binding specificities for antigen presentation.

Meta-populational demes constitute a reservoir for large MHC allele diversity in wild house mice (Mus musculus)

Background

The MHC class I and II loci mediate the adaptive immune response and belong to the most polymorphic loci in vertebrate genomes. In fact, the number of different alleles in a given species is often so large that it remains a challenge to provide an evolutionary model that can fully account for this.

Results

We provide here a general survey of MHC allele numbers in house mouse populations and two sub-species (M. m. domesticus and M. m. musculus) for H2 class I D and K, as well as class II A and E loci. Between 50 and 90% of the detected different sequences constitute new alleles, confirming that the discovery of new alleles is indeed far from complete. House mice live in separate demes with small effective population sizes, factors that were proposed to reduce, rather than enhance the possibility for the maintenance of many different alleles. To specifically investigate the occurrence of alleles within demes, we focused on the class II H2-Aa and H2-Eb exon 2 alleles in nine demes of M. m. domesticus from two different geographic regions. We find on the one hand a group of alleles that occur in different sampling regions and three quarters of these are also found in both sub-species. On the other hand, the larger group of different alleles (56%) occurs only in one of the regions and most of these (89%) only in single demes. We show that most of these region-specific alleles have apparently arisen through recombination and/or partial gene conversion from already existing alleles.

Conclusions

Demes can act as sources of alleles that outnumber the set of alleles that are shared across the species range. These findings support the reservoir model proposed for human MHC diversity, which states that large pools of rare MHC allele variants are continuously generated by neutral mutational mechanisms. Given that these can become important in the defense against newly emerging pathogens, the reservoir model complements the selection based models for MHC diversity and explains why the exceptional diversity exists.

Electronic supplementary material

The online version of this article (10.1186/s12983-018-0266-9) contains supplementary material, which is available to authorized users.

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

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

Clonal polymorphism and high heterozygosity in the celibate genome of the Amazon molly

The extreme rarity of asexual vertebrates in nature is generally explained by genomic decay due to absence of meiotic recombination, thus leading to extinction of such lineages. We explore features of a vertebrate asexual genome, the Amazon molly, Poecilia formosa, and find few signs of genetic degeneration but unique genetic variability and ongoing evolution. We uncovered a substantial clonal polymorphism and, as a conserved feature from its interspecific hybrid origin, a 10-fold higher heterozygosity than in the sexual parental species. These characteristics seem to be a principal reason for the unpredicted fitness of this asexual vertebrate. Our data suggest that asexual vertebrate lineages are scarce not because they are at a disadvantage, but because the genomic combinations required to bypass meiosis and to make up a functioning hybrid genome are rarely met in nature.

Genetic variation and selection of MHC class I loci differ in two congeneric frogs

Major histocompatibility complex (MHC) genes encode proteins in the acquired immune response pathway that often show distinctive selection-driven patterns in wild vertebrate populations. We examined genetic variation and signatures of selection in the MHC class I alpha 1 (A1)- and alpha 2 (A2)-domain encoding exons of two frog congeners [Agalychnis callidryas (n = 20) and A. lemur (n = 20)] from a single locality in Panama. We also investigated how historical demographic processes may have impacted MHC genetic diversity by analyzing a neutral mitochondrial marker. We found that both MHC domains were highly variable in both species, with both species likely expressing three loci. Our analyses revealed different signatures of selection between the two species, most notably that the A. callidryas A2 domain had experienced positive selection while the A2 domain of A. lemur had not. Diversifying selection acted on the same number of A1 and A2 allelic lineages, but on a higher percentage of A1 sites compared to A2 sites. Neutrality tests of mitochondrial haplotypes predominately indicated that the two species were at genetic equilibrium when the samples were collected. In addition, two historical tests of demography indicated both species have had relatively stable population sizes over the past 100,000 years; thus large population size changes are unlikely to have greatly influenced MHC diversity in either species during this time period. In conclusion, our results suggest that the impact of selection on MHC diversity varied between these two closely related species, likely due to a combination of distinct ecological conditions and past pathogenic pressures.

Immunogenetic novelty confers a selective advantage in host-pathogen coevolution

The major histocompatibility complex (MHC) is crucial to the adaptive immune response of vertebrates and is among the most polymorphic gene families known. Its high diversity is usually attributed to selection imposed by fast-evolving pathogens. Pathogens are thought to evolve to escape recognition by common immune alleles, and, hence, novel MHC alleles, introduced through mutation, recombination, or gene flow, are predicted to give hosts superior resistance. Although this theoretical prediction underpins host-pathogen "Red Queen" coevolution, it has not been demonstrated in the context of natural MHC diversity. Here, we experimentally tested whether novel MHC variants (both alleles and functional "supertypes") increased resistance of guppies (Poecilia reticulata) to a common ectoparasite (Gyrodactylus turnbulli). We used exposure-controlled infection trials with wild-sourced parasites, and Gyrodactylus-naïve host fish that were F2 descendants of crossed wild populations. Hosts carrying MHC variants (alleles or supertypes) that were new to a given parasite population experienced a 35-37% reduction in infection intensity, but the number of MHC variants carried by an individual, analogous to heterozygosity in single-locus systems, was not a significant predictor. Our results provide direct evidence of novel MHC variant advantage, confirming a fundamental mechanism underpinning the exceptional polymorphism of this gene family and highlighting the role of immunogenetic novelty in host-pathogen coevolution.

AmpliSAS and AmpliHLA: Web Server Tools for MHC Typing of Non-Model Species and Human Using NGS Data

AmpliSAS and AmpliHLA are web server tools for automatic genotyping of MHC genes from high-throughput sequencing data. AmpliSAS is designed specifically to analyze amplicon sequencing data from non-model species and it is able to perform de-novo genotyping without any previous knowledge of the reference alleles. AmpliHLA is a human-specific version, it performs HLA typing by comparing sequenced variants against human reference alleles from the IMGT/HLA database. Here we describe four genotyping protocols: the first two use amplicon sequencing data to genotype the MHC genes of a passerine bird and human respectively; the third and fourth present the HLA typing of a human cell line starting from RNA and exome sequencing data respectively.

Social behaviour and gut microbiota in red-bellied lemurs (Eulemur rubriventer): In search of the role of immunity in the evolution of sociality

Vertebrate gut microbiota form a key component of immunity and a dynamic link between an individual and the ecosystem. Microbiota might play a role in social systems as well, because microbes are transmitted during social contact and can affect host behaviour. Combining methods from behavioural and molecular research, we describe the relationship between social dynamics and gut microbiota of a group-living cooperative species of primate, the red-bellied lemur (Eulemur rubriventer). Specifically, we ask whether patterns of social contact (group membership, group size, position in social network, individual sociality) are associated with patterns of gut microbial composition (diversity and similarity) between individuals and across time. Red-bellied lemurs were found to have gut microbiota with slight temporal fluctuations and strong social group-specific composition. Contrary to expectations, individual sociality was negatively associated with gut microbial diversity. However, position within the social network predicted gut microbial composition. These results emphasize the role of the social environment in determining the microbiota of adult animals. Since social transmission of gut microbiota has the potential to enhance immunity, microbiota might have played an escalating role in the evolution of sociality.

Shared evolutionary origin of major histocompatibility complex polymorphism in sympatric lemurs

Genes of the major histocompatibility complex (MHC) play a central role in adaptive immune responses of vertebrates. They exhibit remarkable polymorphism, often crossing species boundaries with similar alleles or allelic motifs shared across species. This pattern may reflect parallel parasite-mediated selective pressures, either favouring the long maintenance of ancestral MHC allelic lineages across successive speciation events by balancing selection ("trans-species polymorphism"), or alternatively favouring the independent emergence of functionally similar alleles post-speciation via convergent evolution. Here, we investigate the origins of MHC similarity across several species of dwarf and mouse lemurs (Cheirogaleidae). We examined MHC class II variation in two highly polymorphic loci (DRB, DQB) and evaluated the overlap of gut-parasite communities in four sympatric lemurs. We tested for parasite-MHC associations across species to determine whether similar parasite pressures may select for similar MHC alleles in different species. Next, we integrated our MHC data with those previously obtained from other Cheirogaleidae to investigate the relative contribution of convergent evolution and co-ancestry to shared MHC polymorphism by contrasting patterns of codon usage at functional vs. neutral sites. Our results indicate that parasites shared across species may select for functionally similar MHC alleles, implying that the dynamics of MHC-parasite co-evolution should be envisaged at the community level. We further show that balancing selection maintaining trans-species polymorphism, rather than convergent evolution, is the primary mechanism explaining shared MHC sequence motifs between species that diverged up to 30 million years ago.

Selection on MHC in a Context of Historical Demographic Change in 2 Closely Distributed Species of Tuco-tucos (Ctenomys australis and C. talarum)

Selection necessarily acts within the same current and historical demographic framework as neutral evolutionary processes, and the outcome of the interplay between these forces may vary according to their relative strength. In this study, we compare the variation at a major histocompatibility complex (MHC) locus (DRB exon 2), typically subject to strong diversifying selection, and mitochondrial diversity (D-loop) across populations encompassing the entire distribution of 2 species of South American subterranean rodents: Ctenomys australis and C. talarum (tuco-tucos). Although these species are parapatric along most of their distribution, historically they have followed distinct demographic trajectories associated with sea level changes during the Quaternary. We surveyed 8 populations of C. australis and 15 of C. talarum, from which we analyzed 70 and 212 D-loop haplotypes and 91 and 346 DRB genotypes, respectively. Both species have gone through a recent demographic expansion; however, the signal of this process only encompasses the entire distribution of one of the species: C. australis. While balancing selection on MHC in C. talarum-enhanced DRB diversity at the local level compared to D-loop, although not promoting divergence among populations, in C. australis local diversifying selection may have driven higher population differentiation at DRB than at D-loop. Our findings reinforce the idea that the relative strength of selection acting on MHC genes varies spatially and temporally within and among species, even between species using the same macrohabitat and exposed to similar immune challenges.

Extreme MHC class I diversity in the sedge warbler (Acrocephalus schoenobaenus); selection patterns and allelic divergence suggest that different genes have different functions

Background

Recent work suggests that gene duplications may play an important role in the evolution of immunity genes. Passerine birds, and in particular Sylvioidea warblers, have highly duplicated major histocompatibility complex (MHC) genes, which are key in immunity, compared to other vertebrates. However, reasons for this high MHC gene copy number are yet unclear. High-throughput sequencing (HTS) allows MHC genotyping even in individuals with extremely duplicated genes. This HTS data can reveal evidence of selection, which may help to unravel the putative functions of different gene copies, i.e. neofunctionalization. We performed exhaustive genotyping of MHC class I in a Sylvioidea warbler, the sedge warbler, Acrocephalus schoenobaenus, using the Illumina MiSeq technique on individuals from a wild study population.

Results

The MHC diversity in 863 genotyped individuals by far exceeds that of any other bird species described to date. A single individual could carry up to 65 different alleles, a large proportion of which are expressed (transcribed). The MHC alleles were of three different lengths differing in evidence of selection, diversity and divergence within our study population. Alleles without any deletions and alleles containing a 6 bp deletion showed characteristics of classical MHC genes, with evidence of multiple sites subject to positive selection and high sequence divergence. In contrast, alleles containing a 3 bp deletion had no sites subject to positive selection and had low divergence.

Conclusions

Our results suggest that sedge warbler MHC alleles that either have no deletion, or contain a 6 bp deletion, encode classical antigen presenting MHC molecules. In contrast, MHC alleles containing a 3 bp deletion may encode molecules with a different function. This study demonstrates that highly duplicated MHC genes can be characterised with HTS and that selection patterns can be useful for revealing neofunctionalization. Importantly, our results highlight the need to consider the putative function of different MHC genes in future studies of MHC in relation to disease resistance and fitness.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0997-9) contains supplementary material, which is available to authorized users.

Rapid and reliable screening of HLA-B*15:02 in four Chinese populations using single-tube multiplex real-time PCR assay

AIM

 HLA-B*15:02 screening should be performed to prevent antiepileptic drug induced severe cutaneous adverse reactions in populations of Asian origin. This study aimed to develop fast and reliable HLA-B*15:02 genotyping method and to investigate the distribution of HLA-B*15:02 in different Chinese ethnicities.

MATERIALS & METHODS

A single-tube multiplex real-time PCR assay for HLA-B*15:02 genotyping was established by combining allele-specific primers with TaqMan probes.

RESULTS

The HLA-B*15:02 genotyping results in Bouyei (n = 100) by the established assay were completely consistent with the corresponding PCR sequence-based typing findings. The percentage of HLA-B*15:02 carrier in Bouyei (19%; n = 100) was significantly higher than those of Han (1%; n = 100), Tibetan (0%; n = 100) and Uyghur (0%; n = 50) populations (p < 0.001).

CONCLUSION

The novel method provides rapid, reliable and cost-effective detection of HLA-B*15:02 allele in clinical applications.

MHC VARIATION IN THE ENDANGERED GILA TOPMINNOW

Sequence variation at a major histocompatibility complex (MHC) gene, assumed to be involved in parasite and pathogen resistance, was examined in the endangered Gila topminnow (Poeciliopis o. occidentalis), from the four watersheds where they remain in the United States. This is the first estimate of variation in MHC genes in natural populations of an endangered species. The population that has experienced the most bottlenecks historically was monomorphic for MHC variation. Another population, which earlier had been found to be the only population polymorphic for allozymes, had five MHC alleles, four different from those found in the other populations. Overall, nine different alleles were found. The four populations were highly divergent at MHC with four of the six population pairs not sharing any alleles. However, the magnitude of differentiation between populations on the amino-acid level varied fivefold for the populations that shared no alleles. Using single-stranded conformational polymorphism (SSCP), these alleles segregated consistently with Mendelian expectations in families. Because of the high genetic differentiation between these populations for a potentially adaptive gene, we recommend that the four watersheds be examined further for separate conservation and management.

GENETIC DIVERSITY AND CAPILLARIA HEPATICA (NEMATODA) PREVALENCE IN MICHIGAN DEER MOUSE POPULATIONS

There have been few field tests of the hypothesis that homozygous populations are prone to high levels of disease. I tested for a negative correlation between genetic diversity and parasitism by estimating the allozyme heterozygosity, population density, and proportion of individuals infected by Capillaria hepatica (Nematoda) in nine Michigan populations of deer mice (Peromyscus maniculatus). Parasite prevalence was correlated negatively with heterozygosity when the effects of density were held constant, but was not correlated with population density after controlling for the effects of genetic diversity. These data support the prediction that inbred populations will be more susceptible to parasite infestations.

EVOLUTION OF ALLORECOGNITION IN BOTRYLLID ASCIDIANS INFERRED FROM A MOLECULAR PHYLOGENY

Despite the functional and phyletic ubiquity of highly polymorphic genetic recognition systems, the evolution and maintenance of these remarkable loci remain an empirical and theoretical puzzle. Many clonal invertebrates use polymorphic genetic recognition systems to discriminate kin from unrelated individuals during behavioral interactions that mediate competition for space. Space competition may have been a selective force promoting the evolution of highly polymorphic recognition systems, or preexisting polymorphic loci may have been coopted for the purpose of mediating space competition. Ascidian species in the family Botryllidae have an allorecognition system in which fusion or rejection between neighboring colonies is controlled by allele-sharing at a single, highly polymorphic locus. The behavioral sequence involved in allorecognition varies in a species-specific fashion with some species requiring extensive intercolony tissue integration prior to the allorecognition response, while other species contact opposing colonies at only a few points on the outer surface before resolving space conflicts. Due to an apparent species-specific continuum of behavioral variation in the degree of intercolony tissue integration required for allorecognition, this system lends itself to a phylogenetic analysis of the evolution of an allorecognition system. We constructed a molecular phylogeny of the botryllids based on 18S rDNA sequence and mapped allorecognition behavioral variation onto the phylogeny. Our phylogeny shows the basal allorecognition condition for the group is the most internal form of the recognition reaction. More derived species show progressively more external allorecognition responses, and in some cases loss of some features of internal function. We suggest that external allorecognition appears to be a secondary function of a polymorphic discriminatory system that was already in place due to other selective pressures such as gamete, pathogen, or developmental cell lineage recognition.

Multiple major histocompatibility complex class I genes in Asian anurans: Ontogeny and phylogeny

Amphibians, as the first terrestrial vertebrates, offer a window into early major histocompatibility complex (MHC) evolution. We characterized the MHC class I of two Korean amphibians, the Asiatic toad (Bufo gargarizans) and the Japanese tree frog (Hyla japonica). We found at least four transcribed MHC class I (MHC I) loci, the highest number confirmed in any anuran to date. Furthermore, we identified MHC I transcripts in terrestrial adults, and possibly in aquatic larvae, of both species. We conducted a phylogenetic analysis based on MHC I sequence data and found that B. gargarizans and H. japonica cluster together in the superfamily Nobleobatrachia. We further identified three supertypes shared by the two species. Our results reveal substantial variation in the number of MHC I loci in anurans and suggest that certain supertypes have particular physiochemical properties that may confer pathogen resistance.

Cheetahs have a stronger constitutive innate immunity than leopards

As a textbook case for the importance of genetics in conservation, absence of genetic variability at the major histocompatibility complex (MHC) is thought to endanger species viability, since it is considered crucial for pathogen resistance. An alternative view of the immune system inspired by life history theory posits that a strong response should evolve in other components of the immune system if there is little variation in the MHC. In contrast to the leopard (Panthera pardus), the cheetah (Acinonyx jubatus) has a relatively low genetic variability at the MHC, yet free-ranging cheetahs are healthy. By comparing the functional competence of the humoral immune system of both species in sympatric populations in Namibia, we demonstrate that cheetahs have a higher constitutive innate but lower induced innate and adaptive immunity than leopards. We conclude (1) immunocompetence of cheetahs is higher than previously thought; (2) studying both innate and adaptive components of immune systems will enrich conservation science.

Molecular Mechanisms Regulating Cell Fusion and Heterokaryon Formation in Filamentous Fungi

For the majority of fungal species, the somatic body of an individual is a network of interconnected cells sharing a common cytoplasm and organelles. This syncytial organization contributes to an efficient distribution of resources, energy, and biochemical signals. Cell fusion is a fundamental process for fungal development, colony establishment, and habitat exploitation and can occur between hyphal cells of an individual colony or between colonies of genetically distinct individuals. One outcome of cell fusion is the establishment of a stable heterokaryon, culminating in benefits for each individual via shared resources or being of critical importance for the sexual or parasexual cycle of many fungal species. However, a second outcome of cell fusion between genetically distinct strains is formation of unstable heterokaryons and the induction of a programmed cell death reaction in the heterokaryotic cells. This reaction of nonself rejection, which is termed heterokaryon (or vegetative) incompatibility, is widespread in the fungal kingdom and acts as a defense mechanism against genome exploitation and mycoparasitism. Here, we review the currently identified molecular players involved in the process of somatic cell fusion and its regulation in filamentous fungi. Thereafter, we summarize the knowledge of the molecular determinants and mechanism of heterokaryon incompatibility and place this phenomenon in the broader context of biotropic interactions and immunity.

Temporal and spatial instability in neutral and adaptive (MHC) genetic variation in marginal salmon populations

The role of marginal populations for the long-term maintenance of species’ genetic diversity and evolutionary potential is particularly timely in view of the range shifts caused by climate change. The Centre-Periphery hypothesis predicts that marginal populations should bear reduced genetic diversity and have low evolutionary potential. We analysed temporal stability at neutral microsatellite and adaptive MHC genetic variation over five decades in four marginal Atlantic salmon populations located at the southern limit of the species’ distribution with a complicated demographic history, which includes stocking with foreign and native salmon for at least 2 decades. We found a temporal increase in neutral genetic variation, as well as temporal instability in population structuring, highlighting the importance of temporal analyses in studies that examine the genetic diversity of peripheral populations at the margins of the species’ range, particularly in face of climate change.

No evidence of prenatal diversifying selection at locus or supertype levels in the dog MHC class II loci

BACKGROUND

Despite decades of studying, the mechanisms maintaining high diversity in the genes of the Major Histocompatibility Complex (MHC) are still puzzling scientists. In addition to pathogen recognition and other functions, MHC molecules may act prenatally in mate choice and in maternal-foetal interactions. These interactions are potential selective mechanisms that increase genetic diversity in the MHC. During pregnancy, immune response has a dual role: the foetus represents foreign tissue compared to mother, but histo-incompatibility is required for successful pregnancy. We have studied the prenatal selection in MHC class II loci (DLA-DQA1, DLA-DQB1 and DLA-DRB1) in domestic dogs by comparing the observed and expected offspring genotype proportions in 110 dog families. Several potential selection targets were addressed, including the peptide-binding site, the MHC locus, three-locus haplotype and supertype levels. For the supertype analysis, the first canine supertype classification was created based on in silico analysis of peptide-binding amino-acid polymorphism.

RESULTS

In most loci and levels, no deviation from the expected genotype frequencies was observed. However, one peptide-binding site in DLA-DRB1 had an excess of heterozygotes among the offspring. In addition, if the father shared a DLA-DRB1 allele with the mother, that allele was inherited by the offspring more frequently than expected, suggesting the selective advantage of a histo-compatible foetus, in contrast to our expectations.

CONCLUSIONS

We conclude that there is some evidence of post-copulatory selection at nucleotide site level in the MHC loci of pet dogs. But due to no indication of selection at locus, three-locus, or supertype levels, we estimated that the prenatal selection coefficient is less than 0.3 in domestic dogs and very likely other factors are more important in maintaining the genetic diversity in MHC loci.

The investigation of the origin of Southern Tunisians using HLA genes

The south of Tunisia is characterized by marked ethnic diversity, highlighted by the coexistence of native Berbers with Blacks, Jews and Arab-speaking populations. Despite this heterogeneity, genetic anthropology studies investigating the origin of current Southern Tunisians were rarely reported. We examined human leukocyte antigen (HLA) class I (A, B) and class II (DRB1, DQB1) gene profiles of 250 unrelated Southern Tunisians, and compared them with those of Arab-speaking communities, along with Mediterranean and sub-Sahara African populations using genetic distances, neighbor-joining dendrograms, correspondence and haplotype analysis. In total, 137 HLA alleles were detected, which comprised 32 HLA-A, 52 HLA-B, 32 DRB1 and 21 DQB1 alleles. The most frequent alleles were HLA-A*02:01(18.02%), HLA-B*50:01 (9.11%), HLA-DRB1*07:01 (22.06%) and HLA-DQB1*02:01 (17.21%). All pairs of HLA loci show significant linkage disequilibrium. The four loci depict negative F_nd (the normalized deviate of the homozygosity) values indicating an overall trend to balancing selection. Southern Tunisians appear to be closely related to others Tunisian populations including Berbers, North Africans and Iberians. On the contrary, Southern Tunisians were distinct from Palestinian, Lebanese and Jordanian Middle Eastern Arab-speaking population, despite the deep Arab incursions and Arabization that affected Southern Tunisia. In addition, Southern Tunisians were distant from many sub-Saharan communities, evidenced by genetic distance analysis. Collectively, this indicates a limited genetic contribution of Arab invasion and Black caravans on the makeup of Southern Tunisian gene pool.

De novo transcriptome assembly facilitates characterisation of fast-evolving gene families, MHC class I in the bank vole (Myodes glareolus)

The major histocompatibility complex (MHC) plays a central role in the adaptive immune response and is the most polymorphic gene family in vertebrates. Although high-throughput sequencing has increasingly been used for genotyping families of co-amplifying MHC genes, its potential to facilitate early steps in the characterisation of MHC variation in nonmodel organism has not been fully explored. In this study we evaluated the usefulness of de novo transcriptome assembly in characterisation of MHC sequence diversity. We found that although de novo transcriptome assembly of MHC I genes does not reconstruct sequences of individual alleles, it does allow the identification of conserved regions for PCR primer design. Using the newly designed primers, we characterised MHC I sequences in the bank vole. Phylogenetic analysis of the partial MHC I coding sequence (2-4 exons) of the bank vole revealed a lack of orthology to MHC I of other Cricetidae, consistent with the high gene turnover of this region. The diversity of expressed alleles was characterised using ultra-deep sequencing of the third exon that codes for the peptide-binding region of the MHC molecule. High allelic diversity was demonstrated, with 72 alleles found in 29 individuals. Interindividual variation in the number of expressed loci was found, with the number of alleles per individual ranging from 5 to 14. Strong signatures of positive selection were found for 8 amino acid sites, most of which are inferred to bind antigens in human MHC, indicating conservation of structure despite rapid sequence evolution.

Sexual Selection and the differences between the sexes in Mandrills (Mandrillus sphinx)

Sexual selection has become a major focus in evolutionary and behavioral ecology. It is also a popular research topic in primatology. I use studies of mandrills (Mandrillus sphinx), a classic example of extravagant armaments and ornaments in animals, to exemplify how a long-term, multidisciplinary approach that integrates field observations with laboratory methods can contribute to on-going theoretical debates in the field of sexual selection. I begin with a brief summary of the main concepts of sexual selection theory and the differences between the sexes. I then introduce mandrills and the study population and review mandrill life history, the ontogeny of sex differences, and maternal effects. Next, I focus on male-male competition and female choice, followed by the less well-studied questions of female-female competition and male choice. This review shows how different reproductive priorities lead to very different life histories and divergent adaptations in males and females. It demonstrates how broadening traditional perspectives on sexual selection beyond the ostentatious results of intense sexual selection on males leads to an understanding of more subtle and cryptic forms of competition and choice in both sexes and opens many productive avenues in the study of primate reproductive strategies. These include the potential for studies of postcopulatory selection, female intrasexual competition, and male choice. These studies of mandrills provide comparison and, I hope, inspiration for studies of both other polygynandrous species and species with mating systems less traditionally associated with sexual selection.

MHC-correlated preferences in diestrous female horses (Equus caballus)

Genes of the major histocompatibility complex (MHC) have been shown to influence communication in many vertebrates, possibly with context-specific MHC-correlated reactions. Here we test for MHC-linked female preferences in the polygynous horse (Equus caballus) by repeatedly exposing 19 mares to a group of seven sexually experienced stallions. Each mare was tested four times during two consecutive reproductive cycles, twice during estrus and twice during diestrus. Male plasma testosterone concentrations were determined from weekly blood samples, and equine leukocyte antigen (ELA) class I and II alleles were determined serologically at the end of the experiments. Perception of male attractiveness was strongly dependent on estrous cycle: mean preference scores did not correlate for mares in diestrus and estrus and varied more during estrus than during diestrus. We found elevated female interests for MHC-dissimilar stallions, but only during diestrus, not during estrus. Female preferences were not significantly predicted by mean male testosterone plasma concentrations. However, testosterone concentrations changed during the 11 weeks of the experiment. By the end of the experiment, average testosterone concentration was significantly correlated to the average number of MHC alleles the stallions shared with the mares. We conclude that the MHC affects female preferences for stallions, but non-MHC linked male characteristics can overshadow effects of the MHC during estrus.

sNebula, a network-based algorithm to predict binding between human leukocyte antigens and peptides

Understanding the binding between human leukocyte antigens (HLAs) and peptides is important to understand the functioning of the immune system. Since it is time-consuming and costly to measure the binding between large numbers of HLAs and peptides, computational methods including machine learning models and network approaches have been developed to predict HLA-peptide binding. However, there are several limitations for the existing methods. We developed a network-based algorithm called sNebula to address these limitations. We curated qualitative Class I HLA-peptide binding data and demonstrated the prediction performance of sNebula on this dataset using leave-one-out cross-validation and five-fold cross-validations. This algorithm can predict not only peptides of different lengths and different types of HLAs, but also the peptides or HLAs that have no existing binding data. We believe sNebula is an effective method to predict HLA-peptide binding and thus improve our understanding of the immune system.

Egernia stokesii (gidgee skink) MHC I positively selected sites lack concordance with HLA peptide binding regions

Genes of the major histocompatibility complex (MHC) play an important role in vertebrate disease resistance, kin recognition and mate choice. Mammalian MHC is the most widely characterised of all vertebrates, and attention is often given to the peptide binding regions of the MHC because they are presumed to be under stronger selection than non-peptide binding regions. For vertebrates where the MHC is less well understood, researchers commonly use the amino acid positions of the peptide binding regions of the human leukocyte antigen (HLA) to infer the peptide binding regions within the MHC sequences of their taxon of interest. However, positively selected sites within MHC have been reported to lack correspondence with the HLA in fish, frogs, birds and reptiles including squamates. Despite squamate diversity, the MHC has been characterised in few snakes and lizards. The Egernia group of scincid lizards is appropriate for investigating mechanisms generating MHC variation, as their inclusion will add a new lineage (i.e. Scincidae) to studies of selection on the MHC. We aimed to identify positively selected sites within the MHC of Egernia stokesii and then determine if these sites corresponded with the peptide binding regions of the HLA. Six positively selected sites were identified within E. stokesii MHC I, only two were homologous with the HLA. E. stokesii positively selected sites corresponded more closely to non-lizard than other lizard taxa. The characterisation of the MHC of more intermediate taxa within the squamate order is necessary to understand the evolution of the MHC across all vertebrates.

Development of novel associations between MHC alleles and susceptibility to parasitic infections in an isolated population of an endangered mammal

The role of pathogens in dynamics of endangered species is not fully understood, and the effect of infection often interacts with other processes affecting those species, such as fragmentation and isolation or loss of genetic variation. Small, isolated populations are prone to losing functional alleles due to demographic processes and genetic drift, which may diminish their ability to resist infection if immune genes are affected. Demographic processes may also alter the selective pressure exerted by a parasite, as they influence the rate of parasite transmission between individuals. In the present paper we studied changes in parasite infection levels and genetic variability in an isolated population of spotted suslik (Spermophillus suslicus). Over a three-year period (approx. three generations), when the population size remained relatively stable, we observed a considerable increase in parasite prevalence and infection intensity, followed by the development of novel associations between MHC DRB alleles and parasite burden. Contrary to expectations, the change in MHC allele frequency over time was not consistent with the effect of the allele - for instance, Spsu-DRB*07, associated with higher intensity of infection with a nematode Capillaria sp., increased in frequency from 11.8 to 20.2%. Yet, we found no signatures of selection in the studied loci. Our results show that an isolated, stable population may experience a sudden increase in parasitic infections, resulting in a development of novel associations between MHC alleles and parasite susceptibility/resistance, even though no signatures of selection can be found.

Mate choice for neutral and MHC genetic characteristics in Alpine marmots: different targets in different contexts?

Sexual selection through female mate choice for genetic characteristics has been suggested to be an important evolutionary force maintaining genetic variation in animal populations. However, the genetic targets of female mate choice are not clearly identified and whether female mate choice is based on neutral genetic characteristics or on particular functional loci remains an open question. Here, we investigated the genetic targets of female mate choice in Alpine marmots (Marmota marmota), a socially monogamous mammal where extra‐pair paternity (EPP) occurs. We used 16 microsatellites to describe neutral genetic characteristics and two MHC loci belonging to MHC class I and II as functional genetic characteristics. Our results reveal that (1) neutral and MHC genetic characteristics convey different information in this species, (2) social pairs show a higher MHC class II dissimilarity than expected under random mate choice, and (3) the occurrence of EPP increases when social pairs present a high neutral genetic similarity or dissimilarity but also when they present low MHC class II dissimilarity. Thus, female mate choice is based on both neutral and MHC genetic characteristics, and the genetic characteristics targeted seem to be context dependent (i.e., the genes involved in social mate choice and genetic mate choice differ). We emphasize the need for empirical studies of mate choice in the wild using both neutral and MHC genetic characteristics because whether neutral and functional genetic characteristics convey similar information is not universal.

Fine-Scale Spatial Covariation between Infection Prevalence and Susceptibility in a Natural Population

The prevalence of infection varies dramatically on a fine spatial scale. Many evolutionary hypotheses are founded on the assumption that this variation is due to host genetics, such that sites with a high frequency of alleles conferring susceptibility are associated with higher infection prevalence. This assumption is largely untested and may be compromised at finer spatial scales where gene flow between sites is high. We put this assumption to the test in a natural snail-trematode interaction in which host susceptibility is known to have a strong genetic basis. A decade of field sampling revealed substantial spatial variation in infection prevalence between 13 sites around a small lake. Laboratory assays replicated over 3 years demonstrate striking variation in host susceptibility among sites in spite of high levels of gene flow between sites. We find that mean susceptibility can explain more than one-third of the observed variation in mean infection prevalence among sites. We estimate that variation in susceptibility and exposure together can explain the majority of variation in prevalence. Overall, our findings in this natural host-parasite system argue that spatial variation in infection prevalence covaries strongly with variation in the distribution of genetically based susceptibility, even at a fine spatial scale.

Evaluating the adaptive potential of the European eel: is the immunogenetic status recovering?

The recent increased integration of evolutionary theory into conservation programs has greatly improved our ability to protect endangered species. A common application of such theory links population dynamics and indices of genetic diversity, usually estimated from neutrally evolving markers. However, some studies have suggested that highly polymorphic adaptive genes, such as the immune genes of the Major Histocompatibility Complex (MHC), might be more sensitive to fluctuations in population dynamics. As such, the combination of neutrally- and adaptively-evolving genes may be informative in populations where reductions in abundance have been documented. The European eel (Anguilla anguilla) underwent a drastic and well-reported decline in abundance in the late 20th century and still displays low recruitment. Here we compared genetic diversity indices estimated from neutral (mitochondrial DNA and microsatellites) and adaptive markers (MHC) between two distinct generations of European eels. Our results revealed a clear discrepancy between signatures obtained for each class of markers. Although mtDNA and microsatellites showed no changes in diversity between the older and the younger generations, MHC diversity revealed a contemporary drop followed by a recent increase. Our results suggest ongoing gain of MHC genetic diversity resulting from the interplay between drift and selection and ultimately increasing the adaptive potential of the species.

Does the parasite-mediated selection drive the MHC class IIB diversity in wild populations of European chub (Squalius cephalus)?

The genes of major histocompatibility complex (MHC) provide an excellent opportunity to study host-parasite relationships because they are expected to evolve in response to parasites and variation in parasite communities. In this study, we investigated the potential role of parasite-mediated selection acting on MHC class IIB (DAB) genes in European chub (Squalius cephalus) natural populations. We found significant differences between populations in metazoan parasites, neutral and adaptive genetic diversities. The analyses based on pairwise data revealed that populations with dissimilar MHC allelic profiles were geographically distant populations with significantly different diversity in microsatellites and a dissimilar composition of parasite communities. The results from the generalized estimating equations method (GEE) on the level of individuals revealed that metazoan parasite load in European chub was influenced by the diversity of DAB alleles as well as by the diversity of neutral genetic markers and host traits reflecting condition and immunocompetence. The multivariate co-inertia analysis showed specific associations between DAB alleles and parasite species. DAB1-like alleles were more involved in associations with ectoparasites, while DAB3-like alleles were positively associated with endoparasites which could suggest potential differences between DAB genes caused by different selection pressure. Our study revealed that parasite-mediated selection is not the only variable affecting MHC diversity in European chub; however, we strongly support the role of neutral processes as the main driver of DAB diversity across populations. In addition, our study contributes to the understanding of the evolution of MHC genes in wild living fish.

Surface expression, peptide repertoire, and thermostability of chicken class I molecules correlate with peptide transporter specificity

The chicken major histocompatibility complex (MHC) has strong genetic associations with resistance and susceptibility to certain infectious pathogens. The cell surface expression level of MHC class I molecules varies as much as 10-fold between chicken haplotypes and is inversely correlated with diversity of peptide repertoire and with resistance to Marek's disease caused by an oncogenic herpesvirus. Here we show that the average thermostability of class I molecules isolated from cells also varies, being higher for high-expressing MHC haplotypes. However, we find roughly the same amount of class I protein synthesized by high- and low-expressing MHC haplotypes, with movement to the cell surface responsible for the difference in expression. Previous data show that chicken TAP genes have high allelic polymorphism, with peptide translocation specific for each MHC haplotype. Here we use assembly assays with peptide libraries to show that high-expressing B15 class I molecules can bind a much wider variety of peptides than are found on the cell surface, with the B15 TAPs restricting the peptides available. In contrast, the translocation specificity of TAPs from the low-expressing B21 haplotype is even more permissive than the promiscuous binding shown by the dominantly expressed class I molecule. B15/B21 heterozygote cells show much greater expression of B15 class I molecules than B15/B15 homozygote cells, presumably as a result of receiving additional peptides from the B21 TAPs. Thus, chicken MHC haplotypes vary in several correlated attributes, with the most obvious candidate linking all these properties being molecular interactions within the peptide-loading complex (PLC).

Spatial pattern of adaptive and neutral genetic diversity across different biomes in the lesser anteater (Tamandua tetradactyla)

The genes of the major histocompatibility complex (MHC) code for proteins involved in antigen recognition and activation of the adaptive immune response and are thought to be regulated by natural selection, especially due to pathogen‐driven selective pressure. In this study, we investigated the spatial distribution of MHC class IIDRB exon 2 gene diversity of the lesser anteater (Tamandua tetradactyla) across five Brazilian biomes using next‐generation sequencing and compared the MHC pattern with that of neutral markers (microsatellites). We found a noticeable high level of diversity in DRB (60 amino acid alleles in 65 individuals) and clear signatures of historical positive selection acting on this gene. Higher allelic richness and proportion of private alleles were found in rain forest biomes, especially Amazon forest, a megadiverse biome, possibly harboring greater pathogen richness as well. Neutral markers, however, showed a similar pattern to DRB, demonstrating the strength of demography as an additional force to pathogen‐driven selection in shaping MHC diversity and structure. This is the first characterization and description of diversity of a MHC gene for any member of the magna‐order Xenarthra, one of the basal lineages of placental mammals.