Mating system is correlated with immunogenetic diversity in sympatric species of Peromyscine mice

The number of reproductive partners per individual varies markedly across animal mating systems. This variation may be an important determinant of patterns of immunogenetic diversity, particularly at Major Histocompatibility Complex (MHC) Class I and II loci. To compare immunogenetic variation in taxa with markedly different mating systems, we used RNAseq-generated data to quantify genotypic diversity in three species of Peromyscine rodents: the monogamous California mouse (Peromyscus californicus) and the polygynandrous deer mouse (P. maniculatus) and brush mouse (P. boylii). By sampling populations of these species from multiple localities in California, we were able to conduct replicated analyses of the relationship between mating system and immunogenetic variation. Across the four localities sampled, diversity at MHC Class I and II genes was consistently higher in the two polygynandrous species. We found no evidence that sampling location (i.e., variation in habitat conditions) contributed to observed differences in MHC variation among populations or species. Collectively, our data indicate that immunogenetic variation in Peromyscine mice is associated with reproductive behavior, rather than geographic locality or habitat type. The consistently greater variability detected in the polygynandrous species examined suggests that balancing selection imposed by behaviorally-mediated pathogen exposure is important in maintaining variation at MHC genes in these animals.