Speciation across mountains: Phylogenomics, species delimitation and taxonomy of the Liolaemus leopardinus clade (Squamata, Liolaemidae).

• We found important mito-nuclear discordance in the evolution of the Liolaemus leopardinus group. • Species delimitation using mtDNA under-estimates diversity. • SNP-based delimitation uncovers more accurate and complex diversity. • These patterns are likely due to introgression between taxa in the same mountain chain. • We sink one species and describe a new one for Liolaemus. Organisms rapidly diversifying across unstable environments such as mountain tops provide substantial challenges for resolving evolutionary histories and delimiting species. The Liolaemus leopardinus clade is a group of five species of lizards adapted to high altitudes in central Chile, with most species found in the Andes, but one species, L. frassinettii is found in the independent Costa Cordillera. Despite their allopatric distributions, they display shallow mitochondrial divergences, making phylogenetics and species delimitation of this clade hard to resolve. We use an integrative approach to delimit species by considering morphological data (linear and landmark-based), mitochondrial DNA (mtDNA), and nuclear DNA (Sequences and SNPs collected with ddRADseq). We find strong conflicting signals between phylogenetic analyses of the nuclear and mtDNA data. While mtDNA places L. frassinettii as sister to the rest of the clade, the SNPs support a south to north order of divergences, with southernmost species (new taxon described here) as sister to the rest of the clade. Moreover, species delimitation using mtDNA only supports two species (one in the Costa and one in the Andes), whereas combined analyses using the nuclear data and morphology support multiple Andean taxa, including a new one we describe here. Based on these results, population structure analyses and our knowledge of the geological and climatic history of the Andes, we argue that this mito-nuclear discordance is explained by past introgression among the Andean taxa, likely during glacial periods that forced these lizards to lower altitudes where they would hybridize. The complete isolation between the Costa and Andes cordilleras has prevented any further contact between taxa on either mountain chain. Our study highlights the importance of using multiple lines of evidence to resolve evolutionary histories, and the potential misleading results from relying solely on mtDNA. [ABSTRACT FROM AUTHOR]