Sexual transmission may drive pair similarity of the cloacal microbiome in a polyandrous species.

All animals host a microbial community within and on their reproductive organs, known as the reproductive microbiome. In free-living birds, studies on the sexual transmission of bacteria have typically focused on a few pathogens instead of the bacterial community as a whole, despite a potential link to reproductive function. Theory predicts higher sexual transmission of the reproductive microbiome in females via the males' ejaculates and higher rates of transmission in promiscuous systems. We studied the cloacal microbiome of breeding individuals of a socially polyandrous, sex-role-reversed shorebird, the red phalarope (Phalaropus fulicarius). We expected (i) higher microbial diversity in females compared to males; (ii) low compositional differentiation between sexes; (iii) lower variation in composition between individuals (i.e. microbiome dispersion) in females than in males; (iv) convergence in composition as the breeding season progresses as a consequence of sexual transmission and/or shared habitat use; and (v) higher similarity in microbial composition between social pair members than between two random opposite-sex individuals. We found no or small between-sex differences in cloacal microbiome diversity/richness and composition. Dispersion of predicted functional pathways was lower in females than in males. As predicted, microbiome dispersion decreased with sampling date relative to clutch initiation of the social pair. Microbiome composition was significantly more similar among social pair members than among two random opposite-sex individuals. Pair membership explained 21.5% of the variation in taxonomic composition and 10.1% of functional profiles, whereas temporal and sex effects explained only 0.6%-1.6%. Consistent with evidence of functional convergence of reproductive microbiomes within pairs, some select taxa and predicted functional pathways were less variable between social pair members than between random opposite-sex individuals. As predicted if sexual transmission of the reproductive microbiome is high, sex differences in microbiome composition were weak in a socially polyandrous system with frequent copulations. Moreover, high within-pair similarity in microbiome composition, particularly for a few taxa spanning the spectrum of the beneficial-pathogenic axis, demonstrates the link between mating behaviour and the reproductive microbiome. Our study is consistent with the hypothesis that sexual transmission plays an important role in driving reproductive microbiome ecology and evolution.

Transcending sea turtles: First report of hatching failure in eggs of an Amazonian freshwater turtle with symptoms of the fungal emerging disease fusariosis.

During the last few decades, fungal pathogens have caused devastating population declines across a broad range of taxa. A newly emerging fungal disease, sea turtle egg fusariosis, caused by members of the Fusarium solani species complex (FSSC), has been reported to be responsible for hatching failure in sea turtles worldwide. However, this has not been detected in fresh water turtle species. Here, using relocated clutches and artificial incubation, we report high hatching failure in eggs symptomatic of fusariosis in the yellow-spotted Amazon River turtle (Podocnemis unifilis) inhabiting a pristine environment in the Ecuadorian Amazon. In 2020, we screened 680 eggs of the yellow-spotted Amazon River turtle, relocated from wild nesting areas to artificial nests, for visual symptoms of fusariosis and to estimate hatchability despite infection. We selected 68 eggs sampled in 2019 to confirm Fusarium infection by PCR amplification of the TEF-1α gene and sequenced seven of those amplicons on an Illumina Miseq to assess FSSC membership. We observed fusariosis symptoms in 42% of the 680 eggs. The proportion of symptomatic eggs within nests was negatively linked to the proportion of eggs that hatched. Hatchability was 8% for symptomatic eggs compared with 72% of asymptomatic eggs. Through PCR testing, 58% of symptomatic and 8% of asymptomatic eggs sampled in 2019 tested positive for Fusarium spp., and sequencing revealed that nine sequence variants from three asymptomatic and four symptomatic eggs corresponded to F. keratoplasticum, F. solani and F. falciforme, the three major FSSC pathogens reported in sea turtle egg fusariosis. Our study suggests that hatching failure in eggs linked to symptoms of fusariosis appears to be partially caused by Fusarium pathogens within FSSC in a freshwater turtle. Thus, fusariosis is more widespread among the Testudines than previously reported and is not limited to sea environments, findings of particular conservation concern.

Egg microbiota is the starting point of hatchling gut microbiota in the endangered yellow-spotted Amazon river turtle.

Establishment and development of gut microbiota during vertebrates' early life are likely to be important predictors of health and fitness. Host-parental and host-environment interactions are essential to these processes. In oviparous reptiles whose nests represent a source of the parent's microbial inocula, the relative role of host-selection and stochastic environmental factors during gut microbial assemblage remains unknown. We sampled eggs incubated in artificial nests as well as hatchlings and juveniles (up to 30 days old) of the yellow-spotted Amazon river turtle (Podocnemis unifilis) developing in tubs filled with river water. We examined the relative role of the internal egg microbiota and the abiotic environment on hatchling and juvenile turtle's cloacal microbiota assemblages during the first 30 days of development. A mean of 71% of ASVs in hatched eggs could be traced to the nest environmental microbiota and in turn a mean of 77% of hatchlings' cloacal ASVs were traced to hatched eggs. Between day 5 and 20 of juvenile turtle's development, the river water environment plays a key role in the establishment of the gut microbiota (accounting for a mean of 13%-34.6% of cloacal ASVs) and strongly influences shifts in microbial diversity and abundance. After day 20, shifts in gut microbiota composition were mainly driven by host-selection processes. Therefore, colonization by environmental microbiota is key in the initial stages of establishing the host's gut microbiota which is subsequently shaped by host-selection processes. Our study provides a novel quantitative understanding of the host-environment interactions during gut microbial assemblage of oviparous reptiles.

Evidence of MHC class I and II influencing viral and helminth infection via the microbiome in a non-human primate.

Until recently, the study of major histocompability complex (MHC) mediated immunity has focused on the direct link between MHC diversity and susceptibility to parasite infection. However, MHC genes can also influence host health indirectly through the sculpting of the bacterial community that in turn shape immune responses. We investigated the links between MHC class I and II gene diversity gut microbiome diversity and micro- (adenovirus, AdV) and macro- (helminth) parasite infection probabilities in a wild population of non-human primates, mouse lemurs of Madagascar. This setup encompasses a plethora of underlying interactions between parasites, microbes and adaptive immunity in natural populations. Both MHC classes explained shifts in microbiome composition and the effect was driven by a few select microbial taxa. Among them were three taxa (Odoribacter, Campylobacter and Prevotellaceae-UCG-001) which were in turn linked to AdV and helminth infection status, correlative evidence of the indirect effect of the MHC via the microbiome. Our study provides support for the coupled role of MHC diversity and microbial flora as contributing factors of parasite infection.

Human encroachment into wildlife gut microbiomes.

In the Anthropocene, humans, domesticated animals, wildlife, and their environments are interconnected, especially as humans advance further into wildlife habitats. Wildlife gut microbiomes play a vital role in host health. Changes to wildlife gut microbiomes due to anthropogenic disturbances, such as habitat fragmentation, can disrupt natural gut microbiota homeostasis and make animals vulnerable to infections that may become zoonotic. However, it remains unclear whether the disruption to wildlife gut microbiomes is caused by habitat fragmentation per se or the combination of habitat fragmentation with additional anthropogenic disturbances, such as contact with humans, domesticated animals, invasive species, and their pathogens. Here, we show that habitat fragmentation per se does not impact the gut microbiome of a generalist rodent species native to Central America, Tome's spiny rat Proechimys semispinosus, but additional anthropogenic disturbances do. Indeed, compared to protected continuous and fragmented forest landscapes that are largely untouched by other human activities, the gut microbiomes of spiny rats inhabiting human-disturbed fragmented landscapes revealed a reduced alpha diversity and a shifted and more dispersed beta diversity. Their microbiomes contained more taxa associated with domesticated animals and their potential pathogens, suggesting a shift in potential metagenome functions. On the one hand, the compositional shift could indicate a degree of gut microbial adaption known as metagenomic plasticity. On the other hand, the greater variation in community structure and reduced alpha diversity may signal a decline in beneficial microbial functions and illustrate that gut adaption may not catch up with anthropogenic disturbances, even in a generalist species with large phenotypic plasticity, with potentially harmful consequences to both wildlife and human health.

Phylogeny- and Abundance-Based Metrics Allow for the Consistent Comparison of Core Gut Microbiome Diversity Indices Across Host Species.

The filtering of gut microbial datasets to retain high prevalence taxa is often performed to identify a common core gut microbiome that may be important for host biological functions. However, prevalence thresholds used to identify a common core are highly variable, and it remains unclear how they affect diversity estimates and whether insights stemming from core microbiomes are comparable across studies. We hypothesized that if macroecological patterns in gut microbiome prevalence and abundance are similar across host species, then we would expect that increasing prevalence thresholds would yield similar changes to alpha diversity and beta dissimilarity scores across host species datasets. We analyzed eight gut microbiome datasets based on 16S rRNA gene amplicon sequencing and collected from different host species to (1) compare macroecological patterns across datasets, including amplicon sequence variant (ASV) detection rate with sequencing depth and sample size, occupancy-abundance curves, and rank-abundance curves; (2) test whether increasing prevalence thresholds generate universal or host-species specific effects on alpha and beta diversity scores; and (3) test whether diversity scores from prevalence-filtered core communities correlate with unfiltered data. We found that gut microbiomes collected from diverse hosts demonstrated similar ASV detection rates with sequencing depth, yet required different sample sizes to sufficiently capture rare ASVs across the host population. This suggests that sample size rather than sequencing depth tends to limit the ability of studies to detect rare ASVs across the host population. Despite differences in the distribution and detection of rare ASVs, microbiomes exhibited similar occupancy-abundance and rank-abundance curves. Consequently, increasing prevalence thresholds generated remarkably similar trends in standardized alpha diversity and beta dissimilarity across species datasets until high thresholds above 70%. At this point, diversity scores tended to become unpredictable for some diversity measures. Moreover, high prevalence thresholds tended to generate diversity scores that correlated poorly with the original unfiltered data. Overall, we recommend that high prevalence thresholds over 70% are avoided, and promote the use of diversity measures that account for phylogeny and abundance (Balance-weighted phylogenetic diversity and Weighted Unifrac for alpha and beta diversity, respectively), because we show that these measures are insensitive to prevalence filtering and therefore allow for the consistent comparison of core gut microbiomes across studies without the need for prevalence filtering.

Bioaccumulation of trace elements affects chick body condition and gut microbiome in greater flamingos.

Elevated concentrations of trace elements represent a major concern to wetland ecosystems, since river estuaries are geochemical endpoints that accumulate pollution. Although the negative impact of environmental exposure of highly toxic elements such as Pb and Hg has received substantial attention, we still lack a comprehensive understanding of the effects that these and other common trace elements have on natural populations. We used greater flamingos as a study system within three sites that represent a gradient of pollution. Controlling for environmental sediment exposure, we assessed if signatures of bioaccumulation in feathers for ten trace elements (As, Cd, Cr, Cu, Hg, Ni, Pb, Se, Sn and Zn) are associated with two known proxies of health: body condition and the gut bacterial microbiome. We found evidence of an adverse effect of Se, Hg, and Pb bioaccumulation on body condition. Furthermore, bioaccumulation of the elements As, Cu, Se, Pb and Zn influenced different aspects of the gut microbiome. Bioaccumulation of Se led to a shift in the microbiome composition, largely driven by an enrichment of Bacteroides plebeius, which is linked to the breakdown of sulphated polysaccharides of algae. Bacteroides plebeius was negatively associated with chick body condition, suggesting an adverse effect of a microalgae diet rich in Se. Pb bioaccumulation was linked with a decrease in microbial diversity (adjusted-R2 = 10.4%) and an increase in heterogeneity of the microbial community (adjusted-R2 = 10.5%), an indication of impaired gut homeostasis. As, Cu and Zn had more nuanced effects on gut microbiome heterogeneity according to breeding site and bioaccumulation concentration. Our results therefore suggest that in addition to well-studied elements, bioaccumulation of poorly studied elements also adversely affect health of natural populations.

A novel workflow to improve genotyping of multigene families in wildlife species: An experimental set-up with a known model system.

Genotyping complex multigene families in novel systems is particularly challenging. Target primers frequently amplify simultaneously multiple loci leading to high PCR and sequencing artefacts such as chimeras and allele amplification bias. Most genotyping pipelines have been validated in nonmodel systems whereby the real genotype is unknown and the generation of artefacts may be highly repeatable. Further hindering accurate genotyping, the relationship between artefacts and genotype complexity (i.e. number of alleles per genotype) within a PCR remains poorly described. Here, we investigated the latter by experimentally combining multiple known major histocompatibility complex (MHC) haplotypes of a model organism (chicken, Gallus gallus, 43 artificial genotypes with 2-13 alleles per amplicon). In addition to well-defined 'optimal' primers, we simulated a nonmodel species situation by designing 'cross-species' primers based on sequence data from closely related Galliform species. We applied a novel open-source genotyping pipeline (ACACIA; https://gitlab.com/psc_santos/ACACIA), and compared its performance with another, previously published pipeline (AmpliSAS). Allele calling accuracy was higher when using ACACIA (98.5% versus 97% and 77.8% versus 75% for the 'optimal' and 'cross-species' data sets, respectively). Systematic allele dropout of three alleles owing to primer mismatch in the 'cross-species' data set explained high allele calling repeatability (100% when using ACACIA) despite low accuracy, demonstrating that repeatability can be misleading when evaluating genotyping workflows. Genotype complexity was positively associated with nonchimeric artefacts, chimeric artefacts (nonlinearly by levelling when amplifying more than 4-6 alleles) and allele amplification bias. Our study exemplifies and demonstrates pitfalls researchers should avoid to reliably genotype complex multigene families.

Vessel diameter is related to amount and spatial arrangement of axial parenchyma in woody angiosperms.

Parenchyma represents a critically important living tissue in the sapwood of the secondary xylem of woody angiosperms. Considering various interactions between parenchyma and water transporting vessels, we hypothesize a structure-function relationship between both cell types. Through a generalized additive mixed model approach based on 2,332 woody angiosperm species derived from the literature, we explored the relationship between the proportion and spatial distribution of ray and axial parenchyma and vessel size, while controlling for maximum plant height and a range of climatic factors. When factoring in maximum plant height, we found that with increasing mean annual temperatures, mean vessel diameter showed a positive correlation with axial parenchyma proportion and arrangement, but not for ray parenchyma. Species with a high axial parenchyma tissue fraction tend to have wide vessels, with most of the parenchyma packed around vessels, whereas species with small diameter vessels show a reduced amount of axial parenchyma that is not directly connected to vessels. This finding provides evidence for independent functions of axial parenchyma and ray parenchyma in large vesselled species and further supports a strong role for axial parenchyma in long-distance xylem water transport.

Very high MHC Class IIB diversity without spatial differentiation in the mediterranean population of greater Flamingos.

BACKGROUND: Selective pressure from pathogens is thought to shape the allelic diversity of major histocompatibility complex (MHC) genes in vertebrates. In particular, both local adaptation to pathogens and gene flow are thought to explain a large part of the intraspecific variation observed in MHC allelic diversity. To date, however, evidence that adaptation to locally prevalent pathogens maintains MHC variation is limited to species with limited dispersal and, hence, reduced gene flow. On the one hand high gene flow can disrupt local adaptation in species with high dispersal rates, on the other hand such species are much more likely to experience spatial variation in pathogen pressure, suggesting that there may be intense pathogen mediated selection pressure operating across breeding sites in panmictic species. Such pathogen mediated selection pressure operating across breeding sites should therefore be sufficient to maintain high MHC diversity in high dispersing species in the absence of local adaptation mechanisms. We used the Greater Flamingo, Phoenicopterus roseus, a long-lived colonial bird showing a homogeneous genetic structure of neutral markers at the scale of the Mediterranean region, to test the prediction that higher MHC allelic diversity with no population structure should occur in large panmictic populations of long-distance dispersing birds than in other resident species. RESULTS: We assessed the level of allelic diversity at the MHC Class IIB exon 2 from 116 individuals born in four different breeding colonies of Greater Flamingo in the Mediterranean region. We found one of the highest allelic diversity (109 alleles, 2 loci) of any non-passerine avian species investigated so far relative to the number of individuals and loci genotyped. There was no evidence of population structure between the four major Mediterranean breeding colonies. CONCLUSION: Our results suggest that local adaptation at MHC Class IIB in Greater Flamingos is constrained by high gene flow and high MHC diversity appears to be maintained by population wide pathogen-mediated selection rather than local pathogen-mediated selection. Further understanding of how pathogens vary across space and time will be crucial to further elucidate the mechanisms maintaining MHC diversity in species with large panmictic populations and high dispersal rates.

Home-Made Cost Effective Preservation Buffer Is a Better Alternative to Commercial Preservation Methods for Microbiome Research.

The investigation of wildlife gastrointestinal microbiomes by next-generation sequencing approaches is a growing field in microbial ecology and conservation. Such studies often face difficulties in sample preservation if neither freezing facilities nor liquid nitrogen (LQN) are readily available. Thus, in order to prevent microbial community changes because of bacterial growth after sampling, preservation buffers need to be applied to samples. However, the amount of microbial community variation attributable to the different preservation treatments and potentially affecting biological interpretation is hardly known. Here, we sampled feces of 11 sheep (Ovis aries sp.) by using swabs and analyzed the effect of air-drying, an inexpensive self-made nucleic acid preservation buffer (NAP), DNA/RNA Shield™, and RNAlater®, each together with freezing (for 10 days) or storing at room temperature (for 10 days) prior to 16S rRNA gene high-throughput sequencing to determine bacterial communities. Results revealed that the proportions of operational taxonomic units (OTUs) belonging to a bacterial phylum were affected by the preservation treatments, and that alpha diversities [observed OTUs, Shannon index, and phylogenetic diversity (PD)] were lower in all preservation treatments than in samples taken by forensic swabs and immediately frozen which is considered as the favored preservation treatment in the absence of any logistic constraints. Overall, NAP had better preservation qualities than RNAlater® and DNA/RNA Shield™ making this self-made buffer a valuable solution in wildlife microbiome studies.

A global analysis of parenchyma tissue fractions in secondary xylem of seed plants.

Parenchyma is an important tissue in secondary xylem of seed plants, with functions ranging from storage to defence and with effects on the physical and mechanical properties of wood. Currently, we lack a large-scale quantitative analysis of ray parenchyma (RP) and axial parenchyma (AP) tissue fractions. Here, we use data from the literature on AP and RP fractions to investigate the potential relationships of climate and growth form with total ray and axial parenchyma fractions (RAP). We found a 29-fold variation in RAP fraction, which was more strongly related to temperature than with precipitation. Stem succulents had the highest RAP values (mean ± SD: 70.2 ± 22.0%), followed by lianas (50.1 ± 16.3%), angiosperm trees and shrubs (26.3 ± 12.4%), and conifers (7.6 ± 2.6%). Differences in RAP fraction between temperate and tropical angiosperm trees (21.1 ± 7.9% vs 36.2 ± 13.4%, respectively) are due to differences in the AP fraction, which is typically three times higher in tropical than in temperate trees, but not in RP fraction. Our results illustrate that both temperature and growth form are important drivers of RAP fractions. These findings should help pave the way to better understand the various functions of RAP in plants.

A multi-locus inference of the evolutionary diversification of extant flamingos (Phoenicopteridae).

BACKGROUND: Modern flamingos (Phoenicopteridae) occupy a highly specialized ecology unique among birds and represent a potentially powerful model system for informing the mechanisms by which a lineage of birds adapts and radiates. However, despite a rich fossil record and well-studied feeding morphology, molecular investigations of the evolutionary progression among modern flamingos have been limited. Here, using three mitochondrial (mtDNA) markers, we present the first DNA sequence-based study of population genetic variation in the widely distributed Chilean Flamingo and, using two mtDNA and 10 nuclear (nDNA) markers, recover the species tree and divergence time estimates for the six extant species of flamingos. Phylogenetic analyses include likelihood and Bayesian frameworks and account for potential gene tree discordance. Analyses of divergence times are fossil calibrated at the divergence of Mirandornithes (flamingos + grebes) and the divergence of crown grebes. RESULTS: mtDNA sequences confirmed the presence of a single metapopulation represented by two minimally varying mtDNA barcodes in Chilean flamingos. Likelihood and Bayesian methods recovered identical phylogenies with flamingos falling into shallow-keeled (comprising the Greater, American and Chilean Flamingos) and deep-keeled (comprising the Lesser, Andean and James's Flamingos) sub-clades. The initial divergence among flamingos occurred at or shortly after the Mio-Pliocene boundary (6-3 Ma) followed by quick consecutive divergences throughout the Plio-Pleistocene. There is significant incongruence between the ages recovered by the mtDNA and nDNA datasets, likely due to mutational saturation occurring in the mtDNA loci. CONCLUSION: The finding of a single metapopulation in the widespread Chilean Flamingo confirms similar findings in other widespread flamingo species. The robust species phylogeny is congruent with previous classifications of flamingos based on feeding morphology. Modern phoenicopterids likely originated in the New World with each sub-clade dispersing across the Atlantic at least once. Our divergence time estimates place flamingos among the youngest families of birds, counter to the classical notion of flamingos as among the oldest based on biogeography and the fossil record. Finally, we designate 'Phoeniconaias' as a junior synonym of 'Phoenicoparrus' and redefine the latter genus as containing all flamingos more closely related to Phoenicoparrus andinus than Phoenicopterus roseus.

Evidence for an association between post-fledging dispersal and microsatellite multilocus heterozygosity in a large population of greater flamingos.

Dispersal can be divided into three stages: departure, transience and settlement. Despite the fact that theoretical studies have emphasized the importance of heterozygosity on dispersal strategies, empirical evidence of its effect on different stages of dispersal is lacking. Here, using multi-event capture-mark-recapture models, we show a negative association between microsatellite multilocus heterozygosity (MLH; 10 loci; n = 1023) and post-fledging dispersal propensity for greater flamingos, Phoenicopterus roseus, born in southern France. We propose that the negative effects of inbreeding depression affects competitive ability and therefore more homozygous individuals are more likely to disperse because they are less able to compete within the highly saturated natal site. Finally, a model with the effect of MLH on propensity of post-fledgling dispersers to disperse to the long-distance sites of Africa was equivalent to the null model, suggesting that MLH had low to no effect on dispersal distance. Variations in individual genetic quality thus result in context-dependent heterogeneity in dispersal strategies at each stage of dispersal. Our results have important implications on fitness since sites visited early in life are known to influence site selection later on in life and future survival.

Cryptic female choice favours sperm from major histocompatibility complex-dissimilar males.

Cryptic female choice may enable polyandrous females to avoid inbreeding or bias offspring variability at key loci after mating. However, the role of these genetic benefits in cryptic female choice remains poorly understood. Female red junglefowl, Gallus gallus, bias sperm use in favour of unrelated males. Here, we experimentally investigate whether this bias is driven by relatedness per se, or by similarity at the major histocompatibility complex (MHC), genes central to vertebrate acquired immunity, where polymorphism is critical to an individual's ability to combat pathogens. Through experimentally controlled natural matings, we confirm that selection against related males' sperm occurs within the female reproductive tract but demonstrate that this is more accurately predicted by MHC similarity: controlling for relatedness per se, more sperm reached the eggs when partners were MHC--dissimilar. Importantly, this effect appeared largely owing to similarity at a single MHC locus (class I minor). Further, the effect of MHC similarity was lost following artificial insemination, suggesting that male phenotypic cues might be required for females to select sperm differentially. These results indicate that postmating mechanisms that reduce inbreeding may do so as a consequence of more specific strategies of cryptic female choice promoting MHC diversity in offspring.

Genetic polymorphism in dopamine receptor D4 is associated with early body condition in a large population of greater flamingos, Phoenicopterus roseus.

Body condition is an important determinant of fitness in many natural populations. However, as for many fitness traits, the underlying genes that regulate body condition remain elusive. The dopamine receptor D4 gene (DRD4) is a promising candidate as dopamine is known to play an important role in the regulation of food intake and the metabolism of both glucose and lipids in vertebrates. In this study, we take advantage of a large data set of greater flamingos, Phoenicopterus roseus, to test whether DRD4 polymorphism predicts early body condition (EBC) while controlling for whole-genome effects of inbreeding and outbreeding using microsatellite multilocus heterozygosity (MLH). We typed 670 of these individuals for exon 3 of the homologue of the human DRD4 gene and 10 microsatellite markers. When controlling for the effects of yearly environmental variations and differences between sexes, we found strong evidence of an association between exon 3 DRD4 polymorphisms and EBC, with 2.2-2.3% of the variation being explained by DRD4 polymorphism, whereas there was only weak evidence that MLH predicts EBC. Because EBC is most likely a polygenic trait, this is a considerable amount of variation explained by a single gene. This is to our knowledge, the first study to show an association between exon 3 DRD4 polymorphism and body condition in non-human animals. We anticipate that the DRD4 gene as well as other genes coding for neurotransmitters and their receptors may play an important role in explaining variation in traits that affect fitness.

Cryptic preference for MHC-dissimilar females in male red junglefowl, Gallus gallus.

An increasing number of studies test the idea that females increase offspring fitness by biasing fertilization in favour of genetically compatible partners; however, few have investigated or controlled for corresponding preferences in males. Here, we experimentally test whether male red junglefowl, Gallus gallus, prefer genetically compatible females, measured by similarity at the major histocompatibility complex (MHC), a key gene complex in vertebrate immune function. Theory predicts that because some degree of MHC heterozygosity favours viability, individuals should prefer partners that carry MHC alleles different from their own. While male fowl showed no preference when simultaneously presented with an MHC-similar and an MHC-dissimilar female, they showed a 'cryptic' preference, by allocating more sperm to the most MHC-dissimilar of two sequentially presented females. These results provide the first experimental evidence that males might respond to the MHC similarity of a female through differential ejaculate expenditure. By revealing that cryptic male behaviours may bias fertilization success in favour of genetically compatible partners, this study demonstrates the need to experimentally disentangle male and female effects when studying preferences for genetically compatible partners.