Deep ancestral structure, admixture and peripheral persistence shape diversity in the Papilio bianor species complex

Rapid diversification can generate regional biodiversity through ancestral polymorphism, admixture, and ecological differentiation, while leaving complex genomic signals. We analysed whole-genome variation in 325 butterflies, including 166 individuals of the Papilio bianor species complex and six outgroup species, to reconstruct lineage relationships and demographic history. Nuclear genomic analyses recovered three major clades corresponding to the recognized species and placed the island lineage as sister to continental P. bianor. Analyses of 18,635 autosomal window trees revealed extensive genealogical discordance, with incomplete lineage sorting predominating among continental lineages and the strongest introgression signals involving island related ancestry. Demographic modelling supported early island continental divergence followed by continental diversification and ancient admixture. The best fitting fastsimcoal2 model inferred an admixed origin of the Eastern China lineage, while BPP analyses supported ancient island-related introgression into the Eastern China ancestor, Western Himalaya and southwestern Thailand. Peripheral populations exhibited pronounced genetic differentiation and distinct demographic histories; Western Himalaya additionally showed mitochondrial affinity with Yaeyama despite its placement within the Himalayan nuclear lineage. Together, these findings suggest that introgression among deeply differentiated ancestral populations and incomplete lineage sorting during rapid divergence jointly shaped lineage diversity in the P. bianor complex. This study provides an empirical framework for investigating how lineage diversity emerges rapidly during the early stages of speciation, particularly in groups with deep ancestral population structure.