Differential introgression among loci across a hybrid zone of the intermediate horseshoe bat (Rhinolophus affinis)

Background: Hybrid zones formed by the secondary contact of divergent lineages represent natural laboratories for studying the genetic basis of speciation. Here we tested for patterns of differential introgression among three X-linked and 11 autosomal regions to identify candidate loci related to either reproductive isolation or adaptive introgression across a hybrid zone between two Chinese mainland subspecies of the intermediate horseshoe bat Rhinolophus affinis: R. a. himalayanus and R. a. macrurus. Results: Our results support the previous suggestion that macrurus formed when a third subspecies (R. a. hainanus) recolonized the mainland from Hainan Island, and that himalayanus is the ancestral taxon. However, this overall evolutionary history was not reflected in all loci examined, with considerable locus-wise heterogeneity seen in gene tree topologies, levels of polymorphism, genetic differentiation and rates of introgression. Coalescent simulations suggested levels of lineage mixing seen at some nuclear loci might result from incomplete lineage sorting. Isolation with migration models supported evidence of gene flow across the hybrid zone at one intronic marker of the hearing gene Prestin. Conclusions: We suggest that phylogenetic discordance with respect to the species tree seen here is likely to arise via a combination of incomplete lineage sorting and a low incidence of introgression although we cannot rule out other explanations such as selection and recombination. Two X-linked loci and one autosomal locus were identified as candidate regions related to reproductive isolation across the hybrid zone. Our work highlights the importance of including multiple genomic regions in characterizing patterns of divergence and gene flow across a hybrid zone.