Simulated demographic histories inferred fromfastsimcoal2 . (A) Illustration of the best demographic model (Recentghost) revealed by fastsimcoal2 and the AIC values for TAI-IND-PLA. The model suggests the occurrence of two introgression events with different species that affected the genetic component of TAI during the last 8 kyr. We also included a ghost lineage representing individuals not covered by our sampling scheme. (B) Illustration of the best demographic model (hybrid speciation) revealed by fastsimcoal2 and the AIC values for HSI-IND-AUS. The effective population size (Ne) with 95% CI of each population is depicted on the plot. The number on the arrow represents the proportion of inherited gene flow. TAI, S. taiwanensis ; IND, S. indica ; PLA, S. playfairii ; AUS,S. austrotaiwanensis ; HSI, S. hsiehii ; TAS, S. tashiroi ; BAR, S. barbata .