Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs

Abstract

Subdivided Pleistocene glacial refugia, best known as “refugia within refugia”, provided opportunities for diverging populations to evolve into incipient species and/or to hybridize and merge following range shifts tracking the climatic fluctuations, potentially promoting extensive cytonuclear discordances and “ghost” mtDNA lineages. Here, we tested which of these opposing evolutionary outcomes prevails in northern Iberian areas hosting multiple historical refugia of common frogs (Rana cf. temporaria), based on a genomic phylogeography approach (mtDNA barcoding and RAD‐sequencing). We found evidence for both incipient speciation events and massive cytonuclear discordances. On the one hand, populations from northwestern Spain (Galicia and Asturias, assigned to the regional endemic R. parvipalmata), are deeply‐diverged at mitochondrial and nuclear genomes (~4 My of independent evolution), and barely admix with northeastern populations (assigned to R. temporaria sensu stricto) across a narrow hybrid zone (~25 km) located in the Cantabrian Mountains, suggesting that they represent distinct species. On the other hand, the most divergent mtDNA clade, widespread in Cantabria and the Basque country, shares its nuclear genome with other R. temporaria s. s. lineages. Patterns of population expansions and isolation‐by‐distance among these populations are consistent with past mitochondrial capture and/or drift in generating and maintaining this ghost mitochondrial lineage. This remarkable case study emphasizes the complex evolutionary history that shaped the present genetic diversity of refugial populations, and stresses the need to revisit their phylogeography by genomic approaches, in order to make informed taxonomic inferences.

Italian odonates in the Pandora's box: A comprehensive DNA barcoding inventory shows taxonomic warnings at the Holarctic scale

“…For example, by integrating morphology, multilocus genetics, ecology or other information, the molecularly defined species could be either confirmed or rejected (Dufresnes et al., 2019; Kajtoch, Montagna, & Wanat, 2018). On the other hand, morphological species that are inconsistently delimited at the molecular level could be designated for further analyses, such as geometric morphometrics (Solano et al., 2018) or population genomics (Dufresnes et al., 2020).…”

The toad’s warts: Discordance creates bumpy expectations of mitochondrial‐nuclear evolution between species

“…A number of possible evolutionary explanations for mito‐nuclear discordance have been posited, with mechanisms comprising selective causes, neutral demographic processes or both. The most frequent explanations of discordance (Dufresnes et al., 2020) include sex‐biased dispersal (e.g. Dai, Wang, & Lei, 2013), introgression (e.g.…”

Phylogeographic analysis of shrubby beardtongues reveals range expansions during the Last Glacial Maximum and implicates the Klamath Mountains as a hotspot for hybridization

“…The onset of hybridization often does not have adaptive advantages, however. In cases where hybrids are less fit in their environment than either parent, hybridization may facilitate speciation via reinforcement (Butlin, 1987; Dufresnes et al . 2020).…”