Dating and quantifying past gene flow due to hybridization events in Iberian chubs

Hybridization - the successful reproductive cross between individuals from different species - is increasingly recognized as a relatively common phenomenon in nature, mostly due to the increasing availability of genomic data for many organisms and the development of sophisticated methods to detect past gene flow. As a result, evidence of past gene flow has been detected in the genomes of many species, allowing us to investigate the evolutionary forces that act upon hybrid genomes. Biological systems where multiple instances of hybridization have been detected at different time scales are especially useful for this, as they allow investigating the effects of selection and recombination at different time scales post-hybridization. The Iberian chubs ( Squalius spp. ) present an ideal study system for this. In this group of primary fish inhabiting the rivers of the Iberian Peninsula, multiple instances of past hybridization at different geographical and time scales were reported between different species. However, before we investigate the effects of selection and recombination, it is crucial that we date the hybridization events and quantify the contribution from each parental species. In this study, we infer the demographic history of two Iberian chub hybrid lineages, one resulting from a proposed ancient hybridization event ( S. pyrenaicus ) and another potentially much more recent (São Martinho). We generated whole genome resequencing data for 52 individuals from the two hybrid lineages and their putative parental populations. Our demographic modelling results confirm the hypothesis of ancient hybridization between S. carolitertii and S. tartessicus around 600,000-800,000 years ago, giving rise to the S. pyrenaicus lineage. We also confirm that S. carolitertii is the major parental, with S. tartessicus only contributing <10% to the S. pyrenaicus genome. Furthermore, we confirm the more recent hybridization between S. caetobrigu s and S. pyrenaicus (around 4500 years ago) in the São Martinho population. We observe similar contributions from both parental species (around 50%) to the genome of the individuals of this population. Overall, our results give us a deeper understanding of how and when these past gene flow events happened and the contributions of each parental species to the genomes of the individuals of the hybrid lineages.