The pleistocene species pump past its prime:Evidence from European butterfly sister species

This is the final version. Available on open access from Wiley via the DOI in this recordData availability statement: Read data are available from the ENA at PRJEB43082. Sequence alignments for the COI barcode locus were obtained from the dataset DS-EUGENMAP (dx.doi.org/10.5883/DS-EUGENMAP) on BOLD at www.boldsystems.org and were originally produced by Dincӑ et al., (2021). The script used for calculating diversity and divergence is available at https://github.com/samebdon/orthodiver/blob/master/orthodiver.py.The Pleistocene glacial cycles had a profound impact on the ranges and genetic make-up of organisms. Whilst it is clear that the contact zones that have been described for many sister taxa are secondary and have formed during the last interglacial, it is unclear when the taxa involved began to diverge. Previous estimates based on small numbers of loci are unreliable given the stochasticity of genetic drift and the contrasting effects of incomplete lineage sorting and gene flow on gene divergence. Here we use genome-wide transcriptome data to estimate divergence for 18 sister species pairs of European butterflies showing either sympatric or contact zone distributions. We find that in most cases species divergence predates the mid-Pleistocene transition or even the entire Pleistocene period. We also show that although post divergence gene flow is restricted to contact zone pairs, they are not systematically younger than sympatric pairs. This suggests that contact zones are not limited to the initial stages of the speciation process, but can involve notably old taxa. Finally, we show that mitochondrial and nuclear divergence are only weakly correlated and mitochondrial divergence is higher for contact-zone pairs.Biotechnology & Biological Sciences Research Council (BBSRC)Natural Environment Research Council (NERC)European Union Horizon 202