Skip to main content

Advertisement

Genetic variation at MHC class II loci influences both olfactory signals and scent discrimination in ring-tailed lemurs

Article metrics

Abstract

Background

Diversity at the Major Histocompatibility Complex (MHC) is critical to health and fitness, such that MHC genotype may predict an individual’s quality or compatibility as a competitor, ally, or mate. Moreover, because MHC products can influence the components of bodily secretions, an individual’s body odors may signal its MHC composition and influence partner identification or mate choice. Here, we investigated MHC-based signaling and recipient sensitivity by testing for odor-gene covariance and behavioral discrimination of MHC diversity and pairwise dissimilarity in a strepsirrhine primate, the ring-tailed lemur (Lemur catta).

Methods

First, we coupled genotyping of the MHC class II gene, DRB, with gas chromatography-mass spectrometry of genital gland secretions to investigate if functional genetic diversity is signaled by the chemical diversity of lemur scent secretions. We also assessed if the chemical similarity between individuals correlated with their MHC-DRB similarity. Next, we assessed if lemurs discriminated this chemically encoded, genetic information in opposite-sex conspecifics.

Results

We found that both sexes signaled overall MHC-DRB diversity and pairwise MHC-DRB similarity via genital secretions, but in a sex- and season-dependent manner. Additionally, the sexes discriminated absolute and relative MHC-DRB diversity in the genital odors of opposite-sex conspecifics, suggesting that lemur genital odors function to advertise genetic quality.

Conclusions

In summary, genital odors of ring-tailed lemurs provide honest information about an individual’s absolute and relative MHC quality. Complementing evidence in humans and Old World monkeys, we suggest that reliance on scent signals to communicate MHC quality may be important across the primate lineage.

Background

The Major Histocompatibility Complex (MHC) is an extremely polymorphic group of genes within the adaptive immune system of vertebrates that plays a critical role in disease resistance [93]. Because genetic diversity at the MHC is fundamentally linked to parasite resistance, survivorship, and reproductive success [93, 113], an individual’s MHC genotype is hypothesized to be an important predictor of its quality as a mate. If MHC-based information is recognizable to others, animals could increase their reproductive success by selecting mates that possess particular MHC genotypes, such as diverse alleles or specific alleles which convey disease resistance [74, 91]. Although researchers have found evidence that MHC genotype influences mate choice or its proxies in many species (reviewed in [58]), the mechanism by which animals assess the MHC of conspecifics is still under investigation (reviewed in [98]). Given that the protein products of the MHC can influence body odor, scientists have implicated an olfactory-based mechanism (reviewed in [8, 126]); however, researchers rarely combine chemical and behavioral approaches within the same study to test the purported mechanism of information transfer [70, 79]. Here, using the ring-tailed lemur (Lemur catta) – a strepsirrhine primate for which there is strong evidence of condition-dependent olfactory signaling [10, 16, 23, 44] – we test for both olfactory-based MHC advertisement and recognition.

Because MHC diversity is critical to an individual’s current health and the health of its offspring, potential mates or social partners might be chosen for their MHC diversity (i.e., quality), for their possession of a particular disease-resistant allele or for their MHC dissimilarity relative to the chooser [78, 113]. For mating and social behavior to be influenced by the MHC, however, individuals must both indicate their respective MHC genotype and be able to evaluate the MHC information in the signals of conspecifics [8, 47]. Previously, researchers have shown that condition-dependent signals of quality can be used by both sexes to assess potential partners [18, 54, 89, 90]. Although evidence of correlation with MHC genotype has derived primarily from visual signals, such as antler size [25] or bright coloration [107], chemical signals could prove more reliable for advertising MHC genotype [7, 74, 88, 126]: Notably, because degraded MHC molecules are shed from the cell surface and found in body fluids (e.g. serum, saliva, sweat, urine, and glandular secretions), they may function directly as olfactory cues [8, 79, 111]. MHC molecules may also bind relevant volatile compounds, forming a ligand-MHC molecule complex that may stimulate the olfactory senses ([3, 70, 122], but see [62]). Lastly, the MHC may influence the composition of the host’s microbiota [5, 64, 127], including those dwelling within scent glands that contribute to volatile chemical production [37, 65, 67, 118]. Among taxa that display MHC-associated mate choice, researchers have implicated the operation of an olfactory mechanism in fish [1, 79, 95], reptiles [82], birds [31, 66, 68], and mammals [94, 124], including humans (reviewed in [45, 123]).

The ring-tailed lemur is a fitting model for an odor-based test of MHC advertisement (e.g. [61]) and discrimination of conspecific quality. Endemic to Madagascar, ring-tailed lemurs are an endangered species [2] that, owing to population decline and habitat fragmentation [22, 101], faces the threat of inbreeding and inbreeding depression, whether in the wild or in captivity [17, 42, 85]. They live in female-dominated, multi-male, multi-female societies characterized by strictly seasonal breeding and an elaborate system of olfactory reproductive advertisement [55]. Beyond the male’s specialized antebrachial and brachial scent glands [80, 102] and associated wrist-marking behavior [59], both sexes possess genital scent glands, the secretions of which are unusually chemically complex [106]. We focus on these labial and scrotal secretions because both sexes deposit these secretions through genital marking and these labial and scrotal secretions share ~ 170 volatile compounds [9]. The diversity and relative abundance of these chemicals in these genital secretions contain information about the signaler’s sex, breeding condition, injury status, individual identity, and genome-wide microsatellite diversity (or neutral heterozygosity), as well as its relatedness to other individuals [9, 10, 16, 23, 24, 44, 106]. Moreover, this chemically encoded information is salient and distinguishable to conspecifics [18, 23, 44, 105]. Thus, lemur genital odors honestly advertise at least one measure of genetic quality and relatedness in both sexes.

We combined MHC genotyping with chemical analyses of genital secretions and behavioral tests of scent discrimination to ask 1) if lemurs also advertise their MHC-DRB quality and dissimilarity via chemical cues and 2) if opposite-sex conspecifics can detect this olfactory information. We genotyped captive ring-tailed lemurs (N = 62) at the most diverse class II MHC gene, DRB [41], and analyzed the volatile chemical composition of their genital secretions. We used next generation sequencing to genotype the MHC-DRB gene [41], investigating both allelic MHC-DRB diversity and functional MHC-DRB diversity by collapsing alleles into ‘supertypes,’ or groups of alleles with similar immunogenetic binding properties despite different nucleotide sequences [42]. Using gas chromatography and mass spectrometry, we analyzed both the overall, volatile chemical composition of genital secretions, as well as a subset of compounds, including fatty acids (FAs) and fatty acid esters (FAEs), identified a priori based on their putative linkage to fertility in some female primates [26, 73, 77] and their relation to microsatellite diversity in female ring-tailed lemurs [10]. Lastly, we used behavioral testing to determine if conspecifics can discriminate between absolute diversity and relative dissimilarity in MHC-DRB genotypes based on scent alone.

Results

Signaling of individual MHC quality via odor-gene covariance

We found that both sexes of ring-tailed lemurs signaled their individual MHC-DRB quality via the chemical compounds expressed in their genital secretions. Male MHC-DRB diversity was significantly and positively correlated with chemical diversity (N = 23, Z = 2.17, P = 0.03), regardless of season (Fig. 1; Table 1; Additional file 1: Tables S3A and S4). By contrast, female MHC-DRB diversity was unrelated to overall chemical diversity in either season (N = 20, Z = 0.24, P = 0.81; Fig. 2; Table 1; Additional file 1: Table S3B). Nevertheless, female MHC-DRB diversity was significantly and negatively correlated to the diversity of an important subset of chemicals, FAs, but only during the nonbreeding season (Z = − 3.75, P = 0.001; Fig. 2; Table 1; Additional file 1: Table S5A). Female MHC-DRB diversity was not related to the diversity of FAEs, regardless of season (Z = − 0.34, P = 0.740; Fig. 2; Table 1; Additional file 1: Table S5B). For both sexes, removal of the most MHC-diverse individual did not alter the results (Additional file 1: Tables S3, S4 and S5).

Fig. 1
figure1

Linear regression (black line) showing the seasonal relationships between chemical diversity of all compounds in scrotal secretions and MHC-DRB supertype diversity in male ring-tailed lemurs in the a breeding season (closed circles) and b nonbreeding season (open circles)

Table 1 Relationships between the Simpson index of chemical diversity and MHC-DRB diversity in ring-tailed lemurs across seasons, with significant relationships indicated in bold
Fig. 2
figure2

Linear regression (black line) showing seasonal relationships between different measures of chemical diversity (Simpson index of all compounds: a, b; Simpson index of FAs: c, d; Simpson index of FAEs: e, f) in labial secretions and MHC-DRB supertype diversity in female ring-tailed lemurs in the a, c, e breeding season (closed circles) and b, d, f nonbreeding season (open circles)

Signaling of relatedness via dyadic, odor-gene covariance

In all same-sex lemur dyads, genital olfactory cues encoded information about MHC-DRB distance, but in a season-dependent fashion (Fig. 3; Table 2; Additional file 1: Table S6). After controlling for covariates, chemical distances between MM dyads positively correlated with unique MHC-DRB supertypes during the breeding season (N = 22 males as 231 MM dyads, r = 0.408, P < 0.001, Fig. 3a), but not during the nonbreeding season (N = 20 males as 190 MM dyads, r = − 0.079, P = 0.270, Fig. 3b). Similarly, for FF dyads, we observed a significant, positive correlation during the breeding season between the number of unique MHC-DRB supertypes and chemical distance (N = 17 females as 136 FF dyads, r = 0.313, P < 0.001, Fig. 3c), that was not apparent during the nonbreeding season (N = 18 females as 153 FF dyads, r = 0.027, P = 0.729, Fig. 3d).

Fig. 3
figure3

Linear relationships (black lines) between the chemical distance (relative Euclidean) and the genetic distance (number of unique MHC-DRB supertypes, i.e., MHCsupertype diff) between a-b male-male and c-d female-female ring-tailed lemur dyads during the breeding season (closed circles; a, c) and the nonbreeding seasons (open circles; b, d). The numbers of dyads are provided below each data point

Table 2 Partial Mantel tests showing the seasonal relationships between the relative Euclidian chemical distance (for genital odorants) and MHC-based genetic distance for same-sex (MM and FF) dyads of ring-tailed lemurs

We could not detect any relationship between chemical distance and MHC distance between MF dyads in the breeding season (N = 39 subjects of both sexes combined as 374 MF dyads, r = 0.0014, P = 0.8280; Fig. 4a), but there was a trending negative relationship between mixed-sex dyads during the nonbreeding season (N = 38 subjects of both sexes combined as 360 MF dyads, r = − 0.0099, P = 0.0647; Fig. 4b).

Fig. 4
figure4

Linear relationships, indicated by black line, between the chemical distance (relative Euclidean) and the number of unique MHC-DRB supertypes (i.e., MHCsupertype diff) for male-female dyads of ring-tailed lemurs during the a breeding season (closed circles) and b) nonbreeding season (open circles). The numbers of dyads are provided below each data point

Olfactory discrimination of MHC genotype between mixed-sex conspecifics

Although we could only detect the chemical signaling of MHC-DRB diversity in males in the breeding season, behaviorally both male (Table 3; Fig. 5) and female (Fig. 6) recipients showed significant discrimination between the genital secretions of opposite-sex, conspecific donors based on their possession of different MHC-DRB genotypes. The pattern of response to conspecific secretions, however, differed between the sexes.

Table 3 Relationship between the MHC-DRB genotype of female odorant donors and the behavior that male recipients directed toward the female’s scent mark, with significant relationships indicated in bold. Explanatory variables with superscripts indicate the quadratic variable, whereas those without superscripts are linear
Fig. 5
figure5

Behavioral response (a: sniff mark duration; b: lick mark duration; c: shoulder rub frequency) of male ring-tailed lemurs to the odorants of female conspecifics. The line shows the regression and points are jittered to avoid overlap of the data

Fig. 6
figure6

Behavioral response of female ring-tailed lemurs, i.e., time spent sniffing the area on the dowel adjacent to the odorant of the male conspecifics. The line shows the regression and points are jittered to avoid overlap of the data

Male recipients (N = 18) investigated female secretions more (i.e., spent more time sniffing and licking them) if the donors (N = 29) were of intermediate diversity at the MHC-DRB than if they were at either end of the MHC-DRB diversity spectrum (Table 3; Fig. 5). Additionally, as the relative MHC-DRB dissimilarity of female donors increased, male recipients performed fewer shoulder rubs (Table 3).

Female recipients (N = 9) did not investigate (sniff: Z = − 0.86, P = 0.39; lick: Z = − 1.2, P = 0.23) male scent ‘marks’ according to the MHC-DRB diversity of the donors (N = 17); nevertheless, they differentiated their responses towards the areas adjacent to the male’s mark. Specifically, as the MHC-DRB supertype dissimilarity of the male donor increased, female recipients spent more time sniffing areas adjacent to the mark (slope = 0.226, Z = 2.37, P = 0.018; Fig. 6).

Discussion

Owing to its role in survival and reproductive success, immunogenetic diversity is an important predictor of individual quality and may be signaled via visual or chemical means. Our study provides support for the socially salient, chemical signaling of genetic quality in a strepsirrhine primate. Despite sex differences in the chemical ‘indicators’ of quality and their seasonal emergence, ring-tailed lemurs of both sexes signaled their MHC-DRB diversity and dissimilarity to conspecifics via the volatile component of their genital secretions. Moreover, both sexes were able to use these and potentially other olfactory cues to discriminate relevant information about the MHC genotypes of opposite-sex conspecifics. These results confirm the functional significance of our previous work showing detectable relationships between chemical diversity and microsatellite diversity in both sexes [10, 16, 18]. Our results also provide a foundation from which to explore if, using data on reproductive success from wild populations, ring-tailed lemurs actually choose mates according to diversity or dissimilarity of their MHC-DRB genotype.

Male ring-tailed lemurs appear to advertise their genetic quality both via MHC-DRB diversity and via microsatellite diversity. Moreover, they appear to do so in a similar fashion, in that both measures of genetic diversity were positively correlated with the overall chemical diversity of scrotal secretions. Although the relationship with microsatellite diversity only emerged in the breeding season [16], the MHC odor-gene covariance in males emerged regardless of season (albeit more strongly in the breeding season). Female ring-tailed lemurs, on the other hand, signaled their genetic diversity via certain chemicals, specifically FAs. Previously, we had shown that females signaled their increased microsatellite diversity via a negative relationship with the diversity of FAs, a relationship that was only evident during the breeding season [10]. Here, we show that females also advertise MHC-DRB diversity via a negative relationship with FA diversity, but this time the relationship was evident only during the nonbreeding season. It could be that we lacked the power to detect these relationships reliably in both seasons. Alternatively, it could be that contrasting demands in these seasons influence the differential expression of genetic quality in these odor-gene relationships.

In this female-dominant primate, in which female choice is likely to influence male mating success, males may benefit from advertising their genetic quality to females during the breeding season. Males may also benefit from advertising their quality to male competitors year-round. Females, however, may be relaying different information depending on the season. During the breeding season, signaling genome-wide microsatellite diversity and relatedness may be critical to avoid inbreeding [9, 10]. In contrast, signaling MHC-specific diversity and health during the nonbreeding season might convey competitive ability during periods of intense female-female competition (e.g. [44, 71]) and energetically expensive lactation [83]. Additionally, intragroup female competition for access to resources increases during the nonbreeding season [39, 100]. During these social disputes, the killing of vulnerable infants, committed by both sexes, is a significant risk [19, 48, 53, 56, 60]. Signaling one’s health and vitality may reduce the likelihood of aggressive encounters that could lead to infanticide by competing females (reviewed in [116]).

Our results contrast the lack of odor-gene covariance found in mandrills (Mandrillus sphinx), the only other primate in which a relationship between chemical secretions and MHC diversity has been investigated. In both male and female mandrills, MHC-DRB diversity was unrelated to the chemical diversity of secretions obtained from the surface of the sternal gland [109]. MHC information, however, may be signaled through other aspects of the animals’ olfactory signatures that were not analyzed by these authors. For instance, just as female ring-tailed lemurs signal MHC-DRB and microsatellite diversity through a subset of chemicals (e.g. FAs; [10]), so too might MHC-DRB information be contained in the ratios or relative abundances of specific odoriferous compounds. Alternatively, socially relevant information may be encoded in the non-volatile portion of secretions [6, 11, 52, 63] or be signaled through the composition of the microbiota present in the scent glands [5, 8, 65, 86, 127] and the odorants they produce [37, 67, 118]. Further exploration of individual compounds, specific subsets of chemicals, or the non-volatile fraction of secretions might yield a signaling pattern that conveys information about MHC genotype. Such evidence would support findings that male mandrills appear to use the MHC genotype of a potential mate for mate-guarding decisions [110] and that MHC diversity is correlated with male reproductive success [108].

The chemical composition of lemur genital secretions also signals MHC-DRB dissimilarity between male-male, female-female, and male-female dyads, echoing previous results demonstrating the same pattern for microsatellite diversity [9, 10, 16]. Signaling relatedness to any potential social ‘partner’ is likely to be relevant throughout the year, to avoid related competitors or to beneficially direct nepotism [16, 18]. Signaling relatedness or compatibility to opposite-sex conspecifics would be particularly important during the breeding season, to avoid inbreeding and maximize offspring diversity [13, 81, 119]. Evidence now exists that odorants signal MHC dissimilarity within same-sex and opposite-sex dyads in two taxa formerly thought to be primarily visually oriented, namely birds (black-legged kittiwake: [66]; song sparrows: [112]) and anthropoid primates (mandrills: [109]), suggesting greater relevance of olfactory cues than previously suspected.

Regarding behavior, our male recipients responded most to the scent of females that had intermediate MHC-DRB diversity, and they responded least to the odorants of females that were at the extremes of MHC-DRB diversity. It may be that increased investigation reflects a preference, whereas decreased investigation reflects an aversion. For example, the reduced responsiveness of males could indicate avoidance of extreme inbreeding and outbreeding depression [33, 34, 114]. Increased male investigation could reflect that more processing time was required to decipher the female’s potential as a mate, reflecting a trade-off between speed and discrimination accuracy (reviewed in [21]). For example, rats increase the accuracy of their ability to discriminate between odors the longer they sniff the odor, and, for more difficult discrimination tasks, the rate of increase in accuracy is slower [96]. Accordingly, it may have been more challenging for males to identify the potential quality or compatibility of females that had mid-range MHC-DRB diversity. Previously, in a study of microsatellite diversity, we had found that male ring-tailed lemurs spent more time sniffing the secretions of less-related females [18], a pattern that has since been replicated in chimpanzees [46], and which could be explained as a preference for unrelated females and/or as a greater processing demand. Regardless of the direction of the behavioral responses, both sets of findings indicate that male ring-tailed lemurs are minimally able to discriminate conspecifics according to both overall genetic relatedness and MHC-DRB diversity/dissimilarity. Male choice maybe be important in this species due to several factors, particularly limited availability of fertile females and differences in female quality [84]: In this species, females are strictly seasonal and generally fertile only 1–3 times per year for a period of less than 24 h [32, 121], and often cycle somewhat synchronously with other females in the social group [92](Pereira 1991). Thus, both sexes should be choosy about the competitive effort directed towards their potential partners.

Lastly, our finding that female ring-tailed lemurs spent the most time sniffing the vicinity of secretions from MHC-DRB dissimilar males complements previous work showing that females of other species show greater responsiveness to the scents of more MHC dissimilar males than of more MHC-diverse males (e.g. [1, 14]).

We have confirmed an honest olfactory mechanism of ornamentation and potential mate choice, namely via genital odor-MHC gene covariance and discrimination, in both sexes of ring-tailed lemurs. Olfactory information about immunogenetic quality and similarity may also influence general social behavior, specifically for prioritizing agonistic or nepotistic interactions. Female lemurs are expected to be choosy under the traditional paradigm of sexual selection [120]; however, mate choice may be equally important for male ring-tailed lemurs [84]. Our data extend the potential for olfactory-based MHC discrimination across the primate order and add to a growing body of literature suggesting that choice of social partner or mate may depend on both MHC dissimilarity and diversity [58].

Methods

Subjects

Our subjects (N = 62) derived from three captive populations of ring-tailed lemurs, located at the Duke Lemur Center (DLC, N = 24 males, 24 females) in Durham, NC, USA, the Indianapolis Zoo (N = 4 males, 8 females) in Indianapolis, IN, USA, and the Cincinnati Zoo (N = 2 females) in Cincinnati, OH, USA. All of the animals were healthy adults that were reproductively intact (i.e., neither gonadectomized nor hormonally contracepted) at the time of the study. They were housed in mixed-sex pairs or groups, with similar living conditions and provisioning routines across all three institutions (for more details about DLC housing, see [106]). Notably, all of the animals at all three facilities were fed Purina Monkey Chow with assorted fruits and vegetables and had free access to water. Animal care met with institutional guidelines and was in accordance with regulations of the United States Department of Agriculture. The research protocols were approved by the Institutional Animal Care and Use Committee of Duke University (protocol numbers A245–03-07 & A143–12-05) and by the research directors of each zoo.

MHC genotyping

Using DNA extracted from whole blood or tissue, we genotyped all of the subjects at the MHC-DRB loci using parallel tagged next-generation sequencing [41]. Briefly, blood samples were obtained by staff veterinarians from the femoral vessels of gently hand-restrained subjects or tissue samples were acquired banked from deceased subjects. These samples were stored at − 20 °C until processing. DNA was extracted using either DNA miniprep kits (Sigma, St. Louis, MO, USA) or DNeasy® Blood and Tissue kits (Qiagen, Valencia, CA, USA). We amplified a 171-bp fragment, excluding primers, of the 270-bp second exon of the MHC-DRB gene. This fragment is the most frequently genotyped MHC loci in non-model primate species, especially in lemur species for which genomic data to design primers are scarce (e.g. [50, 57, 87, 103, 115]). Because this fragment excludes several variable amino acids within the MHC-DRB gene, the total MHC-DRB variability may be underestimated. Nonetheless, because the genotyped fragment represents the most variable part of exon 2, we can use this 171-bp fragment as a proxy of diversity across the 6 exons of MHC-DRB. In previous work on ring-tailed lemurs, we have shown that diversity at this MHC-DRB fragment is representative of diversity across other class II MHC genes [41].

To generate MHC-DRB genotypes, we sequenced pooled amplicons using parallel tagged sequencing on two platforms: Ion Torrent PGM® 314v2 chips (Life Technologies, Grand Island, NY, USA) and 454 Titanium® 1/8th lanes (Roche, Nutley, NJ, USA). True MHC-DRB alleles were distinguished from artefacts using a published workflow [41]. Each ring-tailed lemur possessed a mean ± S.D. of 2.22 ± 0.92 MHC-DRB alleles (range = 1–4; see Additional file 1: Table S1, adapted from [42]).

Because of the degeneracy of the genetic code and similarity in the physiochemical properties of some amino acids, researchers can quantify both nucleotide sequence diversity and ‘functional’ diversity, the latter reflecting the diversity of pathogen proteins that an individual’s MHC proteins can bind. We thus organized the MHC-DRB alleles (n = 20) into MHC-DRB ‘supertypes’ (n = 13; [42]). Supertypes are groups of MHC alleles that, despite having different nucleotide sequences, have similar antigen binding properties [27, 103], and, thus, are likely to bind the same subset of pathogen peptides. Owing to the functional overlap in their peptide binding properties, alleles within a supertype are also likely to be subject to identical selection pressures [49].

We classified supertypes using a protocol [27] widely used in primate MHC-DRB supertype classification [49, 103, 107, 108]. We first determined the allelic reading frame by aligning the MHC-DRB sequences with the human HLA-DRB sequence [12] to identify antigen binding sites. Then, we identified any amino acid sites under positive selection using the CODEML analysis in PAML (Version 4.7; [125]). For amino acids identified as being under putative positive selection, we imported their physiochemical properties, including hydrophobicity, steric bulk, polarity, and electronic effects [99], into a matrix in Genesis 1.7.6 [117]. Lastly, using hierarchical clustering via Cosine, Euclidean, and Pearson correlation distance methods, we identified supertypes based on antigen binding similarity. A single supertype was defined as the terminal group with no further branching points. All three distance methods clustered all 64 ring-tailed lemur alleles into 27 identical clusters or supertypes. The range in the number of alleles that were collapsed into each supertype grouping was 1–8, with a mean ± S.D. of 2.01 ± 1.54 alleles.

Genital secretion sample collection

We obtained genital gland secretions from a subset (N = 57) of the subjects, hereafter scent ‘donors’ at two of the facilities (see Additional file 1 for a description of the factors limiting sample collection from all subjects). We collected samples at the DLC over a period of 10 years (2003–2013), including during the breeding and nonbreeding seasons (N = 24 males, 24 females). We also collected samples at the Indianapolis Zoo during the breeding season of 2011 (N = 1 male, 9 females). No secretions were collected from subjects at the Cincinnati Zoo. Because our subjects were in the Northern Hemisphere, we considered samples collected from November to March to be ‘breeding season’ samples and those collected from May to August to be ‘nonbreeding season’ samples [28, 106].

At the DLC, trained handlers carefully caught and gently restrained the animals, which were awake and habituated to these procedures. At the Indianapolis Zoo, collections occurred during the annual physical examinations, performed by Zoo staff members, while the animals were under anesthesia (see Additional file 1 for a discussion of the null effects of handling method on genital secretions). Following published methods [106], we used cotton swabs and forceps, pre-cleaned with methanol and pentane, to collect triplicate samples of genital secretions, per subject, at each collection. We gently rubbed the cotton swab against the scrotal or labial glandular field for 5–10 s, placed the scented swabs in pre-cleaned chromatography vials, and stored the vials at − 80 °C. We have previously shown that individual-specific scent signatures are stable across both years and storage time [23, 30, 44, 106]. Each odorant sample was used only once, for either chemical analyses or bioassay presentation, based upon the season of collection, the number of odorant samples available per individual, and the number of possible recipients to which the odorant could be presented. To maximize the possible bioassay presentations, we prioritized achieving an appropriate sample size for chemical analyses to detect statistical differences rather than analyzing the chemistry of every individual.

Gas chromatography mass spectrometry (GCMS) and chemical diversity indices

All of the chemical analyses were performed on a subset of the genital secretions collected from subjects (N = 43) at the DLC. We used previously published GCMS methods and resulting chemical data to quantify the volatile chemical composition of these secretions (collected from N = 23 males, 20 females; [10, 16, 44, 106]). Briefly, we extracted the volatile components of the secretions into 1.5 ml of methyl-tert-butyl ether, concentrated the extraction, and analyzed the components on a Shimadzu GCMS-QP2010 instrument (Shimadzu Scientific Instruments) equipped with a Shimadzu AOC-20 series autosampler. The compounds were detected using the automatic peak detector (SOLUTION WORKSTATION software, Shimadzu Scientific Instruments) and the peaks individually verified via consultation with the National Institute of Standards and Technology library (for further details, see [30]).

For analyses of the chemical data, we discarded compounds that had inconsistent retention times, or that did not comprise at least 0.05% of the overall area of the GCMS chromatogram. The remaining compounds (n = 203 compounds in scrotal secretions and n = 338 compounds in labial secretions) consisted of fatty acids, fatty acid esters, cholesterol derivatives, alkanes, and other unidentified compounds [10, 16, 44, 106]. To represent the overall chemical composition of lemur genital secretions [16], we used three measures of diversity: richness, the Shannon index, and the Simpson index [69, 75]. Richness reflects the absolute number of compounds present per chromatogram, regardless of relative abundance or rarity. By contrast, the Shannon and Simpson diversity indices reflect the relative abundances in different ways: The Shannon index is primarily influenced by common compounds of intermediate abundance, whereas the Simpson index gives more weight to compounds of the greatest relative abundance [16, 75]. We calculated these diversity indices for each individual’s overall chemical profile.

We also calculated these diversity indices for two subsets of chemicals, fatty acids (FAs) and fatty acid esters (FAEs), which are synthesized from FAs [20, 43]. Because FAs have been linked to fertility in certain female primates ([26, 73, 77], although see [36]), we had selected these types of compounds, a priori, for examining odor-gene covariance in previous studies [10]. We have shown that both chemical subsets, FAs and FAEs, are correlated with microsatellite diversity of female ring-tailed lemurs during the breeding season and, thus, might be indicators of individual quality in lemurs [10]. Here, we examined the three diversity indices for these compounds, specifically, in both sexes (n = 25 FAs in 203 total compounds in scrotal secretions and n = 33 FAs in 338 total compounds in labial secretions; n = 87 and 112 FAEs in male and female genital secretions, respectively; [10, 23]).

Behavioral bioassays

To test if ring-tailed lemurs, hereafter ‘recipients,’ can use the secretions of ‘donors’ to discriminate between the MHC genotypes of opposite-sex conspecifics, we conducted 300 behavioral trials or ‘bioassays’ [16, 40, 105]. We used recipients for whom the odorant donors were ‘unknown,’ defined as never having resided concurrently with the recipient in the same group and/or never having had their secretions presented to the recipient in prior bioassays (see Additional file 1). We thus used recipients (N = 27) from the multiple institutions, including at the DLC (N = 14 males, 5 females), Cincinnati Zoo (N = 2 females), and Indianapolis Zoo (N = 4 males, 2 females), and secretion samples from ‘unknown’ donors of the opposite sex at the DLC (N = 16 males, 20 females) and Indianapolis Zoo (N = 1 male, 9 females).

Following previously established protocols [16, 40, 105], we conducted bioassays during the breeding season of 2011 and 2012. Because the subjects lived socially in multi-chambered enclosures, focal animals were temporarily isolated for bioassays, a process to which they had been accustomed. We encouraged the focal animal into a room by itself, then closed the pass-through between this room and the rest of the enclosure. We allowed samples to thaw at ambient temperature, then secured a row of three fresh wooden dowels to the fence of the animal’s test enclosure (at a 45° angle to the ground and separated by 20 cm). Using pre-cleaned forceps, we removed the thawed swab and rubbed the donor’s secretions (for ~ 10–15 s) on a predetermined dowel. The center dowel served as an unscented control, whereas a ~ 2 cm area (at lemur nose level) of the outer dowels was rubbed with a scented swab. The outer dowels thus carried scent, each from different donors, simulating two naturally placed scent marks.

Each recipient underwent 1–3 trials per day over 4–6 days, with each trial lasting 10 min, ultimately participating in 8–12 trials in total. We presented the secretions to each recipient in a randomized order. We also maximized the number of donor dyads whose secretions could be presented across recipients, while minimizing the number of times we presented secretions from each donor to any recipient (average ± S.D. exposures = 1.85 ± 1.05, range = 0–6). Recipient-donor pairs were chosen blindly with respect to donor location or MHC-DRB genotype, and not all donors were presented to all recipients, owing to logistical constraints described in the Additional file 1. Upon completion of the day’s trials, the recipient was reunited with its group.

The bioassays were videotaped, and the videos were scored using an established ethogram [105], by three observers who were blind to the MHC genotypes of the bioassay donors and recipients. Prior to scoring experimental trials, we calculated inter-observer reliability [72] from five ‘practice’ trials. Differences in the labeling of an event or in the chronology or timing (> 1 s) were considered disagreements [105] and scoring of videos did not commence until inter-observer reliability scores exceeded 90%. The main behavior recorded included investigation (e.g. sniffing and licking) and scent marking behavior (e.g. genital marking and, for males only, shoulder rubbing and wrist marking; Additional file 1: Table S2, adapted from [105]). Sniffing allows the intake of volatile information via the nasal epithelium, whereas licking is generally thought to transport non-volatile chemicals directly from the scent source to the vomeronasal organ [29]. We also recorded where investigatory or scent-marking behavior occurred relative to each scent ‘mark’ (i.e., whether the behavior was directed at the mark itself, adjacent to the mark, but on the dowel, or within 15 cm of the dowel). The placement of countermarks could have particular significance: Overmarking or placing one’s mark directly on top of the original mark might mask the original mark, whereas adjacent-marking or placing one’s mark near the original mark leaves the original mark intact [29].

Statistical analyses

General analytical procedures

To examine the relationships between MHC-DRB genotype and olfactory ornamentation, as well as the ability of ring-tailed lemurs to discriminate MHC genotype via genital secretion, we analyzed the data in a series of generalized linear mixed models (GLMMs), using the package ‘glmmADMB’ (Version 0.7.7) in RStudio (Version 3.2.2; [97]). MHC diversity can be measured in various ways, including as the number of different MHC-DRB nucleotide sequences (or alleles) and as the number of MHC-DRB supertypes, which putatively predicts functional antigen binding capabilities [27, 103]. Because the number of alleles and the number of supertypes are positively correlated, we evaluated these explanatory genetic variables with independent GLMMs and used Akaike information criteria (AIC) values to determine the best-fit model [128]. We considered the model with the lowest AIC value (ΔAIC ≥2; [15]) to be the best and report only those models in the main text, although the AIC values for all models are reported in the Additional file 1: Tables S3, S4, S5 and S7). Because the sexes are dimorphic in their glands [106], and thus, their marking behavior, we treated the sexes separately in the analyses of both chemical diversity and behavioral discrimination.

To examine if the similarity in MHC between individuals was reflected in their chemical similarity, we used partial Mantel tests to compare the number of un-shared or unique MHC-DRB alleles and supertypes to the relative Euclidean distance matrices between male-male (MM), female-female (FF), and male-female (MF) dyads. For consistency, we report results for MHC-DRB supertypes in the main text and results for MHC-DRB alleles in the Additional file 1: Table S6.

Analyses of MHC-DRB diversity and chemical complexity in individual males and females

To examine the relationship between MHC-DRB diversity and chemical complexity of the labial or scrotal secretions, each chemical diversity index was evaluated in a separate series of GLMMs using Gaussian distribution with donor identity as a random variable. Explanatory variables included season (i.e., breeding and nonbreeding) and either the number of MHC-DRB alleles (MHCallele) or the number of MHC-DRB supertypes (MHCsupertype) possessed by an individual donor. Because of skew in the frequency of specific MHC-DRB supertypes (i.e., seven supertypes were found in fewer than five individuals, whereas one supertype was found in more than 85% of individuals), we were unable to examine if possession of a specific supertype could be signaled via the chemical complexity of genital secretions. For both sexes, we also analyzed genetic diversity in relation to the chemical diversity of FAs and FAEs [10]. Where patterns for chemical diversity of all the compounds in a genital secretion reflect the patterns for chemical diversity of FAs and FAEs, we report only the results for the chemical diversity of all compounds.

We were unable to control for neutral heterozygosity estimated via microsatellites (see [10, 16] for microsatellite methods) because microsatellite data were unavailable for > 20% of our subjects. Nevertheless, we assessed the correlation between microsatellite heterozygosity and MHC-DRB diversity to determine if microsatellite diversity might explain the pattern of results. Using a linear regression on the subset of subjects for which both genetic measures of diversity were available (N = 36), we found no correlation between microsatellite heterozygosity and the number of MHC-DRB alleles within an individual (slope = 1.056, F = 1.51, P = 0.227), and no correlation between microsatellite heterozygosity and the number of MHC-DRB supertypes (slope = − 0.8812, F = 1.461, P = 0.235; see also [42]). As we previously had detected no relationships between individual chemical diversity and adult age, the month of collection within season, or DLC housing condition [16], we did not include these co-variables in our analyses of the relationship between MHC-DRB diversity and chemical diversity.

The captive ring-tailed lemur population lacks genetic diversity compared to wild populations [42]. Notably, our captive subjects included only two (N = 1 male, 1 female) relatively MHC-diverse individual. Because these individuals were representative of the average MHC-DRB diversity present in wild populations [42], we did not consider them to be outliers; nevertheless, to verify that they were not driving the association between MHC-DRB diversity and chemical diversity, we re-ran the final GLMMs after removing the most diverse individuals the datasets (N = 1; Additional file 1: Tables S3, S4 and S5). Because all three measures of chemical diversity (i.e., richness, Shannon index, and Simpson index) showed the same patterns when compared to either measure of genetic diversity (i.e., allele or supertype number), we report only the GLMM with the lowest AIC values [15] in the main text.

Analysis of MHC relatedness and chemical similarity between all possible dyads

We used partial Mantel tests to investigate if the chemical similarity between dyads reflected the similarity in their MHC genotypes. First, we calculated matrices of genetic distances using the number of different MHC alleles and supertypes between each dyad. We then estimated the chemical distances between pairs of individuals, by analyzing all of the chemical compounds identified in secretion profiles of either sex (n = 203 compounds for males: [16]; n = 338 for females: [10]), or only those compounds shared by both sexes (n = 170 compounds: [9]). We calculated relative Euclidean distance matrices for same-sex (MM or FF) and mixed-sex (MF) dyads, respectively, using PC-ORD (version 7.0, [76]), and following published protocols [9, 16]. We calculated matrices separately as follows: breeding season (for N = 22 males, there were 231 MM dyads; for N = 17 females, there were 136 FF dyads; and for N = 39 males and females, there were 374 MF dyads); nonbreeding season (for N = 20 males, there were 190 MM dyads; for N = 18 females, there were 153 FF dyads; and for N = 38 males and females, there were 360 MF dyads). Because MM and FF matrices were square, we assessed linear relationships between chemical and MHC distances using partial Mantel tests in FSTAT (version 2.9.3.2, with 10,000 randomizations; [38]). As in previous studies [10, 16], we controlled for potentially confounding covariates, including the subject’s age, social housing condition, and the month of secretion sample collection. For the MF comparisons, we first generated full matrices using all possible MM, FF, and MF pairs (breeding season: n = 704 dyads; nonbreeding season: n = 741 dyads). We then extracted chemical, genetic, and covariate information for MF dyads only. Unlike MM and FF matrices, the MF matrix was not square. Therefore, we assessed relationships with 10,000 Spearman’s correlation permutation tests using the JMUOUTLIER package in R (Version 1.3; [35]), as in the study by Slade et al. [112].

Lastly, to confirm that our results were not being driven by the overall genetic similarity between dyads, rather than by allelic sharing at the MHC-DRB loci, we assessed the correlation between MHC similarity within dyads (i.e., the number of unique or unshared MHC-DRB alleles and supertypes between two individuals) with dyad relatedness, as measured by the Queller and Goodnight index (IDQG calculated in [10]). Although dyad relatedness was significantly and negatively correlated with MHC dissimilarity for both the number of MHC-DRB alleles (n = 629 dyads, slope = − 0.71, T-value = − 4.21, P = 0.000029) and the number of MHC-DRB supertypes (n = 629 dyads, slope = − 0.67, T-value = − 4.23, P = 0.000027), the negative relationships explained less than 3% of the variance in either correlation (R2 = 0.026 and R2 = 0.026, respectively). Because the partial Mantel tests for both the unique MHC alleles and the unique supertypes showed similar patterns, we report supertype results in the main text and allelic results in the Additional file 1.

Behavioral analyses of mixed-sex, recipient-donor combinations

We explored the relationship between the recipients’ behavioral responses to donor secretions and measures of absolute and relative MHC-DRB diversity between the mixed-sex, recipient-donor dyads using a separate series of GLMMs for each behavioral response, with a negative binomial distribution and log link function. In each GLMM, we controlled for the random variables of trial number on a given day (i.e., 1–3), the number of times that a recipient had been presented with the secretion from a given donor over the course of the study (i.e., 1–6), as well as the secretion donor ID nested under secretion recipient.

To test for odorant discrimination, we used two measures of dissimilarity and sequence divergence between each recipient-donor dyad, the number of unique alleles and the number of unique supertypes between dyads (Additional file 1: Table S8; [51, 104, 108]). We also used the donor’s number of MHC alleles to examine if conspecifics that had the greatest MHC diversity were distinguished by their scent alone, regardless of their genetic dissimilarity (Additional file 1: Table S8). Lastly, we examined non-linear relationships between MHC diversity and dissimilarity between dyads by including the quadratic forms of all genetic explanatory variables in our GLMMs. Quadratic terms were retained only if the AIC value was better than the GLMM that included only linear terms (Additional file 1: Table S7).

Although we recorded both the frequency and duration of all behavior, we analyzed only frequencies for events and durations for states [4, 72]. In our analyses, we excluded all recipient behavior that occurred in < 5% of trials and any behavior that was not directed significantly more toward the test dowels over the control dowel, as determined via Wilcoxon signed-rank tests [18]. Ultimately, for male recipients, we analyzed the duration of time spent in proximity to the dowels and sniffing and licking the mark and surrounding areas, as well as the frequency of shoulder rubs, and wrist marking the area adjacent to the mark. For females, we analyzed the duration of time spent sniffing and licking the mark and the adjacent area. For each behavioral response, we report the genetic explanatory variable (e.g. unshared MHC-DRB alleles, unshared MHC-DRB supertypes, or the number of donor supertypes) with the lowest AIC value. The AIC values for other models are reported in the Additional file 1: Table S7. We also verified that the behavioral responses to odorants were comparable for samples collected at different facilities (i.e., DLC or Indianapolis Zoo) and across trials, regardless of trial order (Additional file 1).

Availability of data and materials

Analyses reported in this article can be reproduced using the following datasets (i.e., CSV files of MHC, chemical, and behavioral data from bioassays) that have been deposited into Dryad at DOI: https://doi.org/10.5061/dryad.cb010bf.

Abbreviations

DLC:

Duke Lemur Center

FAEs:

Fatty acid esters

FAs:

Fatty acids

FF dyad:

Female-female dyad

GCMS:

Gas chromatography mass spectrometry

GLMM:

General linear mixed model

MF dyad:

Male-female dyad

MHC:

Major Histocompatibility Complex

MHC-DRB:

Major Histocompatibility Complex DRB gene

MM dyad:

Male-male dyad

References

  1. 1.

    Aeschlimann PB, Haberli MA, Reusch TBH, Boehm T, Milinski M. Female sticklebacks, Gasterosteus aculeatus, use self-reference to optimize MHC allele number during mate selection. Behav Ecol Sociobiol. 2003;54:119–26.

  2. 2.

    Andriaholinirina N, Baden A., Blanco M, Chikhi L, Cooke A, Davies N, Dolch R, Donati G, Ganzhorn J, Golden C, Groeneveld LF, Hapke A, Irwin M, Johnson S, Kappeler P, King T, Lewis R, Louis EE, Markolf M, Mass V, Mittermeier RA, Nichols R, Patel E, Rabarivola CJ, Raharivololona B, Rajaobelina S, Rakotoarisoa G, Rakotomanga B, Rakotonanahary J, Rakotondrainibe H, Rakotondratsimba G, Rakotondratsimba M, Rakotonirina L, Ralainasolo FB, Ralison J, Ramahaleo T, Ranaivoarisoa JF, Randrianahaleo SI, Randrianambinina B, Randrianarimanana L, Randrianasolo H, Randriatahina G, Rasamimananana H, Rasolofoharivelo T, Rasoloharijaona S, Ratelolahy F, Ratsimbazafy J, Ratsimbazafy N, Razafindraibe H, Razafindramanana J, Rowe N, Salmona J, Seiler M, Volampeno S, Wright P, Youssouf J, Zaonarivelo J, Zaramody A (2014) Lemur catta. The IUCN red list of threatened species 2014: e.T11496A62260437. Downloaded 05 February 2019.

  3. 3.

    Aksenov AA, Gojova A, Zhao W, Morgan JT, Sankaran S, Sandrock CE, Davis CE. Characterization of volatile organic compounds in human leukocyte antigen heterologous expression systems: a cell’s “chemical odor fingerprint”. Chembiochem. 2012;13:1053–9.

  4. 4.

    Altmann J. Observational study of behavior: sampling methods. Behaviour. 1974;49:227–67.

  5. 5.

    Archie E, Theis K. Animal behaviour meets microbial ecology. Anim Behav. 2011;82:425–36.

  6. 6.

    Beynon RJ, Veggerby C, Payne CE, Robertson DHL, Gaskell SJ, Humphries RE, Hurst JL. Polymorphism in major urinary proteins: molecular heterogeneity in a wild mouse population. J Chem Ecol. 2002;28(7):1429–46.

  7. 7.

    Blaustein AR. Sexual selection and mammalian olfaction. Am Nat. 1981;117:1006–10.

  8. 8.

    Boehm T, Zufall F. MHC peptides and the sensory evaluation of genotype. Trends Neurosci. 2006;29:100–7.

  9. 9.

    Boulet M, Charpentier MJE, Drea CM. Decoding an olfactory mechanism of kin recognition and inbreeding avoidance in a primate. BMC Evol Biol. 2009;9(1):281–92.

  10. 10.

    Boulet M, Crawford JC, Charpentier MJE, Drea CM. Honest olfactory ornamentation in a female-dominant primate. J Evol Biol. 2010;23:1558–63.

  11. 11.

    Brennan PA, Zufall F. Pheromonal communication in vertebrates. Nature. 2006;444:308–15.

  12. 12.

    Brown JH, Jardetzky TS, Gorga JC, Stern LJ, Urban RG, Strominger JL, Wiley DC. Three-dimensional structure of the human class II histocompatibility antigen HLA-DR1. Nature. 1993;364(6432):33–9.

  13. 13.

    Brown JL. The new heterozygosity theory of mate choice and the MHC. Genetica. 1999;104:215–21.

  14. 14.

    Burger D, Thomas S, Aepli H, Dreyer M, Fabre G, Marti E, Sieme H, Robinson MR, Wedekind C. Major histocompatibility complex-linked social signaling affects female fertility. Proc R Soc Lond B Biol Sci. 2017;284:20171824.

  15. 15.

    Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. New York: Springer; 2002.

  16. 16.

    Charpentier MJE, Boulet M, Drea CM. Smelling right: the scent of male lemurs advertises genetic quality and relatedness. Mol Ecol. 2008a;17:3225–33.

  17. 17.

    Charpentier MJE, Williams CV, Drea CM. Inbreeding depression in ring-tailed lemurs (Lemur catta): genetic diversity predicts parasitism, immunocompetence, and survivorship. Conserv Genet. 2008b;9:1605–15.

  18. 18.

    Charpentier MJE, Crawford JC, Boulet M, Drea CM. Message ‘scent’: lemurs detect the genetic relatedness and quality of conspecifics via olfactory cues. Anim Behav. 2010;80:101–8.

  19. 19.

    Charpentier MJE, Drea CM. Victims of infanticide and conspecific bite wounding in a female-dominant primate: a long-term study. PLoS One. 2013;8:e82830.

  20. 20.

    Cheng JB, Russell DW. Mammalian wax biosynthesis II. Expression cloning of wax synthase cDNAs encoding a member of the acyltransferase enzyme family. J Biol Chem. 2004;279:37798–807.

  21. 21.

    Chittka L, Skorupski P, Raine NE. Speed-accuracy tradeoffs in animal decision making. Trends Ecol Evol. 2009;27(7):400–7.

  22. 22.

    Clarke TA, Gray O, Gould L, Burrell AS. Genetic diversity of the ring-tailed lemur (Lemur catta) in south-Central Madagascar. Folia Primatol. 2015;86:76–84.

  23. 23.

    Crawford JC, Boulet M, Drea CM. Smelling wrong: hormonal contraception in lemurs alters critical female odour cues. Proc R Soc Lond B Biol Sci. 2011;278(1702):122–30.

  24. 24.

    Crawford J, Drea CM. Baby on board: olfactory cues indicate pregnancy and fetal sex in a non-human primate. Biol Lett. 2015;11(2):20140831.

  25. 25.

    Ditchkoff SS, Lochmiller RL, Masters RE, Hoofer SR, Van Den Bussche RA. Major histocompatibility complex associated variation in secondary sexual traits of white-tailed deer (Odocoileus virginianus): evidence for good genes advertisement. Evolution. 2001;55:616–25.

  26. 26.

    Doty RL, Ford M, Preti G, Huggins GR. Changes in intensity and pleasantness of human vaginal odors during menstrual-cycle. Science. 1975;190:1316–8.

  27. 27.

    Doytchinova IA, Flower D. In silico identification of supertypes for class II MHCs. J Immunol. 2005;174:7085–95.

  28. 28.

    Drea CM. Sex and seasonal differences in aggression and steroid secretion in Lemur catta: are socially dominant females hormonally ‘masculinized’? Horm Behav. 2007;51(4):555–67.

  29. 29.

    Drea CM. D’Scent of man: a comparative survey of primate chemosignaling in relation to sex. Horm Behav. 2015;68:117–33.

  30. 30.

    Drea CM, Boulet M, delBarco-Trillo J, Greene LK, Sacha CR, Goodwin TE, Dubay GR. The “secret” in secretions: methodological considerations in deciphering primate olfactory communication. Am J Primatol. 2013;75(7):621–42.

  31. 31.

    Ekblom R, Saether SA, Grahn M, Fiske P, Kålås JA, Höglund J. Major histocompatibility complex variation and mate choice in a lekking bird, the great snipe (Gallinago media). Mol Ecol. 2004;13:3821–8.

  32. 32.

    Evans CS, Goy RW. Social behaviour and reproductive cycles in captive ring-tailed lemurs (Lemur catta). J Zool. 1968;156:181–97.

  33. 33.

    Evans ML, Dionne M, Miller KM, Bernatchez L. Mate choice for major histocompatibility complex genetic divergence as a bet-hedging strategy in the Atlantic salmon (Salmo salar). Proc R Soc Lond B Biol Sci. 2012;279:379–86.

  34. 34.

    Forsberg LA, Dannewitz J, Petersson E, Grahn M. Influence of genetic dissimilarity in the reproductive success and mate choice of brown trout - females fishing for optimal MHC dissimilarity. J Evol Biol. 2007;20:1859–69.

  35. 35.

    Garren S. jmuOutlier: permutation tests for nonparametric statistics. In: R package version 1.1; 2016. See http://CRAN.R-project.org/package=jmuOutlier.

  36. 36.

    Goldfoot DA, Kravetz MA, Goy RW, Freeman SK. Lack of effect of vaginal lavages and aliphatic acids on ejaculatory responses in rhesus macaques: behavioral and chemical analyses. Horm Behav. 1976;7:1–27.

  37. 37.

    Gorman ML, Nedwell DB, Smith RM. An analysis of the contents of the anal scent pockets of Herpestes auropunctatus (Carnivora: Veverridae). J Zool. 1974;172:389–99.

  38. 38.

    Goudet J. 2001. FSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3). See http://www2.unil.ch/popgen/softwares/fstat.html.

  39. 39.

    Gould L, Sussman RW, Sauther ML. Demographic and life-history patterns in a population of ring-tailed lemurs (Lemur catta) at Beza Mahafaly reserve, Madagascar: a 15-year perspective. Am J Phys Anthropol. 2003;120(2):182–94.

  40. 40.

    Greene LK, Grogan KE, Smyth KN, Adams CA, Klager SA, Drea CM. Mix it and fix it: functions of composite olfactory signals in ring-tailed lemurs. R Soc Open Sci. 2016;3:160076.

  41. 41.

    Grogan KE, McGinnis GJ, Sauther ML, Cuozzo FP, Drea CM. Next-generation genotyping of hypervariable loci in many individuals of a non-model species: technical and theoretical implications. BMC Genomics. 2016;17:204.

  42. 42.

    Grogan KE, Sauther ML, Cuozzo FP, Drea CM. Genetic wealth, population health: major histocompatibility complex variation in captive and wild ring-tailed lemurs (Lemur catta). Ecol Evol. 2017;7(19):7638–49.

  43. 43.

    Hargrove JL, Greenspan P, Hartle DK. Nutritional significance and metabolism of very long chain fatty alcohols and acids from dietary waxes. Exp Biol Med. 2004;229(3):215–26.

  44. 44.

    Harris RL, Boulet M, Grogan KE, Drea CM. Costs of injury for scent signaling in a strepsirrhine primate. Sci Rep. 2018;8:9880.

  45. 45.

    Havlicek J, Roberts SC. MHC-correlated mate choice in humans: a review. Psychoneuroendocrinology. 2009;34(4):497–512.

  46. 46.

    Henkel S, Setchell JM. Group and kin recognition via olfactory cues in chimpanzees (Pan troglodytes). Proc R Soc B Biol Sci. 2018;285:20181527.

  47. 47.

    Hettyey A, Hegyi G, Puurtinen M, Hoi H, Török J, Penn DJ. Mate choice for genetic benefits: time to put the pieces together. Ethology. 2010;116:1–9.

  48. 48.

    Hood LC. Infanticide among ringtailed lemurs (Lemur catta) at Berenty reserve, Madagascar. Am J Primatol. 1994;33:65–9.

  49. 49.

    Huchard E, Weill M, Cowlishaw G, Raymond M, Knapp LA. Polymorphism, haplotype composition, and selection in the Mhc-DRB of wild baboons. Immunogenetics. 2008;60:585–98.

  50. 50.

    Huchard E, Albrecht C, Schliehe-Diecks S, Baniel A, Roos C, Kappeler PM, Brameier M. Large-scale MHC class II genotyping of a wild lemur population by next generation sequencing. Immunogenetics. 2012;64(12):895–913.

  51. 51.

    Huchard E, Baniel A, Schliehe-Diecks S, Kappeler PM. MHC-disassortative mate choice and inbreeding avoidance in a solitary primate. Mol Ecol. 2013;22:4071–86.

  52. 52.

    Hurst JL, Payne CE, Nevison CM, Marie AD, Humphries RE, Robertson DHL, Cavaggioni A, Beynon RJ. Individual recognition in mice mediated by major urinary proteins. Nature. 2001;414:631–4.

  53. 53.

    Ichino S. Attacks on a wild infant ring-tailed lemur (Lemur catta) by immigrant males at Berenty, Madagascar: interpreting infanticide by males. Am J Primatol. 2005;67:267–72.

  54. 54.

    Johansson BG, Jones TM. The role of chemical communication in mate choice. Biol Rev Camb Philos Soc. 2007;82:265–89.

  55. 55.

    Jolly A. Lemur behavior: a Madagascar field study. Chicago: University of Chicago Press; 1966.

  56. 56.

    Jolly A, Caless S, Cavigelli S, Gould L, Pereira ME, Pitts A, Pride RE, Rabenandrasana HD, Walker JD, Zafison T. Infant killing, wounding and predation in Eulemur and Lemur. Int J Primatol. 2000;21:21–40.

  57. 57.

    Kaesler E, Kappeler PM, Brameier M, Demeler J, Kraus C, Rakotoniaina JH, Hämäläinen AM, Huchard E. Shared evolutionary origin of MHC polymorphism in sympatric lemurs. Mol Ecol. 2017;26(20):5629–45.

  58. 58.

    Kamiya T, O’Dwyer K, Westerdahl H, Senior A, Nakagawa S. A quantitative review of MHC-based mating preference: the role of diversity and dissimilarity. Mol Ecol. 2014;23(21):5151–63.

  59. 59.

    Kappeler PM. To whom it may concern: the transmission and function of chemical signals in Lemur catta. Behav Ecol Sociobiol. 1998;42:411.421.

  60. 60.

    Kittler K, Dietzel S. Female infanticide and female-directed lethal targeted aggression in a group of ring-tailed lemurs (Lemur catta). Primate Biol. 2016;3(2):41–6.

  61. 61.

    Knapp LA, Robson J, Waterhouse JS. Olfactory signals and the MHC: a review and a case study in Lemur catta. Am J Primatol. 2006;68(6):568–84.

  62. 62.

    Kwak J, Willse A, Preti G, Yamazaki K, Beauchamp GK. In search of the chemical basis for MHC odourtypes. Proc R Soc Lond B Biol Sci. 2010;277:2417–25.

  63. 63.

    Kwatra SG & Drea CM. 2007. Proteomic analysis of ringtailed lemur scent gland secretions: glandular- and individual-specific protein profiles, poster presented at the 30th meeting of the American Society of Primatologists. Winston-Salem, NC. (Abstract: American Journal of Primatology 69 supplement 159: 108-109.

  64. 64.

    Lanyon CV, Rushton SP, O’Donnell AG, Goodfellow M, Ward AC, Petrie M, Jensen SP, Gosling LM, Penn DJ. Murine scent mark microbial communities are genetically determined. FEMS Microbiol Ecol. 2007;59:576–83.

  65. 65.

    Leclaire S, Nielson JF, Drea CM. Bacterial communities in meerkat anal scent secretions vary with host sex, age, and group membership. Behav Ecol. 2014a;25(4):996–1004.

  66. 66.

    Leclaire S, van Dongen WFD, Voccia S, Merkling T, Ducamp C, Hatch SA, Blanchard P, Danchin E, Wagner RH. Preen secretions encode information on MHC similarity in certain sex-dyads in a monogamous seabird. Sci Rep. 2014b;4:6920.

  67. 67.

    Leclaire S, Jacob S, Greene LK, Dubay GR, Drea CM. Social odours covary with bacterial community in the anal secretions of wild meerkats. Sci Rep. 2017a;7:3240.

  68. 68.

    Leclaire S, Strandh M, Mardon J, Westerdahl H, Bonadonna F. Odour-based discrimination of similarity at the major histocompatibility complex in birds. Proc R Soc Lond B Biol Sci. 2017b;284:20162466.

  69. 69.

    Legendre P, Legendre L. Numerical Ecology. Amsterdam: Elsevier Science BV; 1998.

  70. 70.

    Leinders-Zufall T, Brennan P, Widmayer P, Chandramani P, Maul-Pavicic A, Jager M, Li X-H, Breer H, Zufall F, Boehm T. MHC class 1 peptides as chemosensory signals in the vomeronasal organ. Science. 2004;306:1033–7.

  71. 71.

    López-Idiáquez D, Vergara P, Fargallo JA, Martinez-Padilla J. Female plumage coloration signals status to conspecifics. Anim Behav. 2016;121:101–6.

  72. 72.

    Martin P, Bateson P. Measuring behavior: an introductory guide. New York: Cambridge University Press; 2007.

  73. 73.

    Matsumoto-Oda A, Oda R, Hayashi Y, Murakami H, Maeda N, Kumazakid K, Shimizud K, Matsuzawa T. Vaginal fatty acids produced by chimpanzees during menstrual cycles. Int J Primatol. 2003;74:75–9.

  74. 74.

    Mays HL, Hill GE. Choosing mates: good genes versus genes that are a good fit. Trends Ecol Evol. 2004;19:554–9.

  75. 75.

    McCune B, Grace JB, Urban DL. Analysis of ecological communities. Gleneden Beach: MjM Software Design; 2002.

  76. 76.

    McCune B, Mefford MJ. PC-ORD. In: Multivariate analysis of ecological data version 7. Gleneden Beach: MjM Software Design; 2016.

  77. 77.

    Michael RP, Bonsall RW, Warner P. Human vaginal secretions: volatile fatty acid content. Science. 1974;186:1217–9.

  78. 78.

    Milinski M. Arms races, ornaments and fragrant genes: the dilemma of mate choice in fishes. Neurosci Biobehav Rev. 2014;46:567–72.

  79. 79.

    Milinski M, Griffiths S, Wegner KM, Reusch TBH, Haas-Assenbaum A, Boehm T. Mate choice decisions of stickleback females predictably modified by MHC peptide ligands. Proc Natl Acad Sci U S A. 2005;102:4414–8.

  80. 80.

    Montagna W, Yun JS. The skin of primates X. the skin of the ring-tailed lemur (Lemur catta). Am J Phys Anthropol. 1962;20:95–118.

  81. 81.

    Neff BD, Pitcher TE. Genetic quality and sexual selection: an integrated framework for good genes and compatible genes. Mol Ecol. 2005;14:19–38.

  82. 82.

    Olsson M, Madsen T, Ujvari B, Wapstra E. Fecundity and MHC affects ejaculation tactics and paternity bias in sand lizards. Evolution. 2004;58(4):906–9.

  83. 83.

    O’Mara MT, Hickey CM. The development of sex differences in ring-tailed lemur feeding ecology. Behav Ecol Sociobiol. 2014;68:1273–86.

  84. 84.

    Parga JA. Male mate choice in Lemur catta. Int J Primatol. 2006;27(1):107–31.

  85. 85.

    Pastorini J, Sauther ML, Sussman RW, Gould L, Cuozzo FP, Fernando P, Nievergelt CM, Mundy NI. Comparison of the genetic variation of captive ring-tailed lemurs with a wild population in Madagascar. Zoo Biol. 2015;34(5):463–72.

  86. 86.

    Pearce DS, Hoover BA, Jennings S, Nevitt GA, Docherty KM. Morphological and genetic factors shape the microbiome of a seabird species (Oceanodroma leucorhoa) more than environmental and social factors. Microbiome. 2017;5(1):146.

  87. 87.

    Pechouskova E, Dammhahn M, Brameier M, Fichtel C, Kappeler PM, Huchard E. MHC class II variation in a rare and ecological specialist mouse lemur reveals lower allelic richness and contrasting selection patterns compared to a generalist and widespread sympatric congener. Immunogenetics. 2015;67(4):229–45.

  88. 88.

    Penn DJ. The scent of genetic compatibility: sexual selection and the major histocompatibility complex. Ethology. 2002;108:1–21.

  89. 89.

    Penn D, Potts WK. Chemical signals and parasite-mediated sexual selection. Trends Ecol Evol. 1998a;13:391–5.

  90. 90.

    Penn D, Potts W. How do major histocompatibility complex genes influence odor and mating preferences? Adv Immunol. 1998b;69:411–36.

  91. 91.

    Penn DJ, Potts WK. The evolution of mating preferences and major histocompatibility complex genes. Am Nat. 1999;153:145–64.

  92. 92.

    Pereira ME, Weiss ML. Female mate choice, male migration, and the threat of infanticide in ringtailed lemurs. Behav Ecol Sociobiol. 1991;28(2).

  93. 93.

    Piertney SB, Oliver MK. The evolutionary ecology of the major histocompatibility complex. Heredity. 2006;96:7–21.

  94. 94.

    Radwan J, Tkacz A, Kloch A. MHC and preferences for male odour in the bank vole. Ethology. 2008;114:827–33.

  95. 95.

    Reusch TB, Haberli MA, Aeschlimann PB, Milinski M. Female sticklebacks count alleles in a strategy of sexual selection explaining MHC polymorphism. Nature. 2001;414:300–2.

  96. 96.

    Rinbert D, Koulakov A, Gelperin A. Speed-accuracy tradeoff in olfaction. Neuron. 2006;51:351–8.

  97. 97.

    RStudio Team. RStudio: integrated development for R. Boston: RStudio, Inc.; 2017. See http://www.rstudio.com/

  98. 98.

    Ruff JS, Nelson AC, Kubinak JL, Potts WK. MHC signaling during social communication. Pgs 290-313. In: López-Larrea C, editor. Self and Nonself. New York: Springer; 2012.

  99. 99.

    Sandberg M, Eriksson L, Jonsson J, Sjöström M, Wold S. New chemical descriptors relevant for the design of biologically active peptides. A multivariate characterization of 87 amino acids. J Med Chem. 1998;41(14):2481–91.

  100. 100.

    Sauther ML. The dynamics of feeding competition in wild populations of ringtailed lemurs (Lemur catta). In: Kappeler PM, Ganzhorn JU, editors. Lemur social systems and their ecological basis. New York: Plenum; 1993. p. 135–52.

  101. 101.

    Sauther ML, Gould L. Going, going, gone..Is the iconic ring-tailed lemur (Lemur catta) headed for imminent extinction? Primate Conserv. 2016;30:89–101.

  102. 102.

    Schilling A. A study of marking behaviour in Lemur catta. In: Martin RD, Doyle GA, Walker AC, editors. Prosimian Biology. Pittsburgh: University of Pittsburgh Press; 1974. p. 347–62.

  103. 103.

    Schwensow N, Fietz J, Dausmann K, Sommer S. MHC-associated mating strategies and the importance of overall genetic diversity in an obligate pair-living primate. Evol Ecol. 2007;22:617–36.

  104. 104.

    Schwensow N, Eberle M, Sommer S. Compatibility counts: MHC-associated mate choice in a wild promiscuous primate. Proc R Soc Lond B Biol Sci. 2008;275:555–64.

  105. 105.

    Scordato ES, Drea CM. Scents and sensibility: information content of olfactory signals in the ringtailed lemur (Lemur catta). Anim Behav. 2007;73(2):301–14.

  106. 106.

    Scordato ES, Dubay G, Drea CM. Chemical composition of scent marks in the ringtailed lemur (Lemur catta): glandular differences, seasonal variation and individual signatures. Chem Senses. 2007;32:493–504.

  107. 107.

    Setchell JM, Charpentier M, Abbott KA, Wickings EJ, Knapp LA. Is brightest best? Testing the Hamilton-Zuk hypothesis in mandrills. Int J Primatol. 2009;30:825–44.

  108. 108.

    Setchell JM, Charpentier MJE, Abbott KM, Wickings EJ, Knapp LA. Opposites attract: MHC-associated mate choice in a polygynous primate. J Evol Biol. 2010;23(1):136–48.

  109. 109.

    Setchell JM, Vaglio S, Abbott KM, Moggi-Cecchi J, Boscaro F, Pieraccini G, Knapp LA. Odour signals major histocompatibility complex genotype in an Old World monkey. Proc R Soc Lond B Biol Sci. 2011;278(1703):274–80.

  110. 110.

    Setchell JM, Richards SA, Abbott KM, Knapp LA. Mate-guarding by male mandrills (Mandrillus sphinx) is associated with female MHC genotype. Behav Ecol. 2016;27(6):1756–66.

  111. 111.

    Singh PB, Brown RE, Roser B. MHC antigens in urine as olfactory recognition cues. Nature. 1987;327:161–4.

  112. 112.

    Slade JWG, Watson MJ, Kelly TR, Gloor GB, Bernards MA. Chemical composition of preen wax reflects major histocompatibility complex similarity in songbirds. Proc R Soc B Biol Sci. 2016;283:20161966.

  113. 113.

    Sommer S. The importance of immune gene variability (MHC) in evolutionary ecology and conservation. Front Zool. 2005a;2:16.

  114. 114.

    Sommer S. Major histocompatibility complex and mate choice in a monogamous rodent. Behav Ecol Sociobiol. 2005b;58:181–9.

  115. 115.

    Sommer S, Courtiol A, Mazzoni CJ. MHC genotyping of non-model organisms using next-generation sequencing: a new methodology to deal with artefacts and allelic dropout. BMC Genomics. 2013;14(1):542.

  116. 116.

    Stockley P, Campbell A. Female competition and aggression: interdisciplinary perspectives. Philos Trans R Soc B. 2013;368:20130073.

  117. 117.

    Sturn A, Quackenbush J, Trajanoski Z. Genesis: cluster analysis of microarray data. Bioinformatics. 2002;18(11):207–8.

  118. 118.

    Theis KR, Venkataraman A, Dycus JA, Koonter KD, Schmitt-matzen EN. Symbiotic bacteria appear to mediate hyena social odors. Proc Natl Acad Sci U S A. 2013;110(49):19832–7.

  119. 119.

    Tregenza T, Wedell N. Genetic compatibility, mate choice and patterns of parentage. Mol Ecol. 2000;9:1013–27.

  120. 120.

    Trivers RL. Parental investment and sexual selection. In: Campbell B, editor. Sexual Selection and the Descent of Man. Chicago: Aldine; 1972. p. 136–79.

  121. 121.

    Van Horn RN, Resko JA. The reproductive cycle of the ringtailed lemur (Lemur catta): sex steroid levels and sexual receptivity under controlled photoperiods. Endocrinology. 1977;101:1579–86.

  122. 122.

    Willse A, Belcher AM, Preti G, Wahl JH, Thresher M, Yang P, Yamazaki K, Beauchamp GK. Identification of major histocompatibility complex-regulated body odorants by statistical analysis of a comparative gas chromatography/mass spectrometry experiment. Anal Chem. 2005;77:2348–61.

  123. 123.

    Winternitz J, Abbate JL, Huchard E, Havlicek J, Garamszegi LZ. Patterns of MHC-dependent mate selection in humans and non-human primates: a meta-analysis. Mol Ecol. 2016;26(2):668–88.

  124. 124.

    Yamazaki K, Boyse EA, Mike V, Thaler HT, Mathieson BJ, Abbott J, Boyse J, Zayas ZA, Thomas L. Control of mating preferences in mice by genes in the major histocompatibility complex. J Exp Med. 1976;144:1324–35.

  125. 125.

    Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24:1586–91.

  126. 126.

    Ziegler A, Kentenich H, Uchanska-Ziegler B. Female choice and the MHC. Trends Immunol. 2005;26:496–502.

  127. 127.

    Zomer S, Dixon SJ, Xu Y, Jensen SP, Wang H, Lanyon CV, O’Donnell AG, Clare AS, Gosling LM, Penn DJ, Brereton RG. Consensus multivariate methods in gas chromatography mass spectrometry and denaturing gradient gel electrophoresis: MHC-congenic and other strains of mice can be classified according to the profiles of volatiles and microflora in their scent-marks. Analyst. 2009;134:114–23.

  128. 128.

    Zuur AF, Ieno EN, Walker NJ, Saveliev AA, Smith GM. Mixed effects models and extensions in ecology with R. New York: Springer; 2009.

Download references

Acknowledgements

We would like to thank the staff and animal technicians of the Duke Lemur Center, the Indianapolis Zoo, and the Cincinnati Zoo for their assistance during sample collection and bioassay trials. We are especially grateful to Erin Ehmke, David Brewer, and Britt Keith at the Duke Lemur Center, to Lynn Villers, Robert Shumaker, and Holly Blaylock at the Indianapolis Zoo and to Terri Roth and Ronald Evans at the Cincinnati Zoo. We could not have accomplished this study without them. We would also like to thank Jillian Wisse for her help with conducting the bioassays, as well as Laura Damiani and Bernice Kwan for their help in scoring the videos and two anonymous reviewers for their suggestions. This manuscript is DLC publication #1433.

Funding

This work was financially supported by the National Science Foundation (BCS #0409367, IOS #0719003) to CMD, the Duke University Center for Science Education to KEG, Duke University to CMD and KEG, and NIH (F32 GM123634–01) to KEG. These funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

CMD and KEG conceived of the idea for this study with help from MB & RLH. KEG produced the MHC-DRB genotypes and conducted the behavioral bioassays, and MB & RLH produced the chemical data. KEG, MB, & RLH performed the analyses and KEG wrote the original draft of the manuscript. CMD critically revised the manuscript with assistance from KEG, RLH, and MB. All authors have approved the final manuscript for publication.

Correspondence to Kathleen E. Grogan.

Ethics declarations

Ethics approval and consent to participate

The research protocols were approved by the Institutional Animal Care and Use Committee of Duke University (protocol numbers A245–03-07 & A143–12-05) and by the research directors of Cincinnati Zoo & Botanical Garden and the Indianapolis Zoo. Written informed consent to have zoo animals participate in this study was given by the research directors of each zoo after approval of the protocols.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Supplementary methods and results. (DOCX 45 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Grogan, K.E., Harris, R.L., Boulet, M. et al. Genetic variation at MHC class II loci influences both olfactory signals and scent discrimination in ring-tailed lemurs. BMC Evol Biol 19, 171 (2019) doi:10.1186/s12862-019-1486-0

Download citation

Keywords

  • MHC-DRB
  • Genetic diversity
  • Strepsirrhine primate
  • Chemical signal
  • Scent mark
  • Sexual selection