Latitudinal Patterns in Immunogenetic and Microbiome Diversity in two anuran species: Bufo bufo and Bufo spinosus

Evolutionary and demographic processes such as selection, drift and migration shape the genetic variation of populations. Genetic diversity is often lower in populations toward higher latitudes. This decrease potentially threatens their survival as several factors are putting more pressure on the populations, including the spread of infectious diseases. In this study, we combined whole-genome re-sequencing with MHC class II genotyping and skin microbiome profiling in Bufo bufo and B. spinosus, two closely related European toad species. We investigated the underlying immunogenetic and microbial variation resulting from different demographic histories and environmental conditions to identify their potential impact on infection outcomes in these two species. We found lower immunogenetic diversity in B. bufo compared to B. spinosus, with highly significant differences in genes related to adaptive and innate immunity. We found lower overall MHC class II diversity and skin microbiome diversity at the species level in B. bufo, compared with B. spinosus. In contrast, at the individual level, B. bufo showed higher MHC allelic diversity and greater diversity in the core skin microbiota than B. spinosus. Together, our findings suggest that divergence in immunogenetic background and host-associated microbial communities may underlie differences in susceptibility to emerging infectious diseases. This integrative framework provides new insight into how host genetics and microbial communities jointly influence disease outcomes across environmental gradients.