Complex Patterns of Diversification in the Gray Zone of Speciation: Model-Based Approaches Applied to Patagonian Liolaemid Lizards (Squamata: Liolaemus kingii clade).

In this study we detangled the evolutionary history of the Patagonian lizard clade Liolaemus kingii, coupling dense geographic sampling and novel computational analytical approaches. We analyzed nuclear and mitochondrial data (restriction site-associated DNA sequencing and cytochrome b) to hypothesize and evaluate species limits, phylogenetic relationships, and demographic histories. We complemented these analyses with posterior predictive simulations to assess the fit of the genomic data to the multispecies coalescent model. We also employed a novel approach to time-calibrate a phylogenetic network. Our results show several instances of mito-nuclear discordance and consistent support for a reticulated history, supporting the view that the complex evolutionary history of the kingii clade is characterized by extensive gene flow and rapid diversification events. We discuss our findings in the contexts of the "gray zone" of speciation, phylogeographic patterns in the Patagonian region, and taxonomic outcomes. [Model adequacy; multispecies coalescent; multispecies network coalescent; phylogenomics; species delimitation.].

ExplorATE: a new pipeline to explore active transposable elements from RNA-seq data.

MOTIVATION: Transposable elements (TEs) are ubiquitous in genomes and many remain active. TEs comprise an important fraction of the transcriptomes with potential effects on the host genome, either by generating deleterious mutations or promoting evolutionary novelties. However, their functional study is limited by the difficulty in their identification and quantification, particularly in non-model organisms. RESULTS: We developed a new pipeline [explore active transposable elements (ExplorATE)] implemented in R and bash that allows the quantification of active TEs in both model and non-model organisms. ExplorATE creates TE-specific indexes and uses the Selective Alignment (SA) to filter out co-transcribed transposons within genes based on alignment scores. Moreover, our software incorporates a Wicker-like criteria to refine a set of target TEs and avoid spurious mapping. Based on simulated and real data, we show that the SA strategy adopted by ExplorATE achieved better estimates of non-co-transcribed elements than other available alignment-based or mapping-based software. ExplorATE results showed high congruence with alignment-based tools with and without a reference genome, yet ExplorATE required less execution time. Likewise, ExplorATE expands and complements most previous TE analyses by incorporating the co-transcription and multi-mapping effects during quantification, and provides a seamless integration with other downstream tools within the R environment. AVAILABILITY AND IMPLEMENTATION: Source code is available at https://github.com/FemeniasM/ExplorATEproject and https://github.com/FemeniasM/ExplorATE_shell_script. Data available on request. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Rapid Radiation and Rampant Reticulation: Phylogenomics of South American Liolaemus Lizards.

Understanding the factors that cause heterogeneity among gene trees can increase the accuracy of species trees. Discordant signals across the genome are commonly produced by incomplete lineage sorting (ILS) and introgression, which in turn can result in reticulate evolution. Species tree inference using the multispecies coalescent is designed to deal with ILS and is robust to low levels of introgression, but extensive introgression violates the fundamental assumption that relationships are strictly bifurcating. In this study, we explore the phylogenomics of the iconic Liolaemus subgenus of South American lizards, a group of over 100 species mostly distributed in and around the Andes mountains. Using mitochondrial DNA (mtDNA) and genome-wide restriction site-associated DNA sequencing (RADseq; nDNA hereafter), we inferred a time-calibrated mtDNA gene tree, nDNA species trees, and phylogenetic networks. We found high levels of discordance between mtDNA and nDNA, which we attribute in part to extensive ILS resulting from rapid diversification. These data also reveal extensive and deep introgression, which combined with rapid diversification, explain the high level of phylogenetic discordance. We discuss these findings in the context of Andean orogeny and glacial cycles that fragmented, expanded, and contracted species distributions. Finally, we use the new phylogeny to resolve long-standing taxonomic issues in one of the most studied lizard groups in the New World.[Andes; ddRADSeq; introgression; lizards; mtDNA; reptiles; SNPs.].

An integrative approach to address species limits in the southernmost members of the Liolaemus kingii group (Squamata: Liolaemini).

Recent conceptual and methodological advances have enabled an increasing number of studies to address the problem of species delimitation in a comprehensive manner. This is of particular interest in cases of species whose divergence times are recent and/or effective population sizes are large, where the conclusions obtained from a single source of evidence may lead to erroneous estimations of true species numbers or incorrect assignment of individuals to species. Iguanian lizards of the Liolaemus kingii group (13 species) comprise an important component of the endemic fauna of Patagonia. The southernmost species of this group (namely L. baguali, L. escarchadosi, L. sarmientoi, and L. tari) show widely overlapping distributions across southern Patagonia, also, their phylogenetic relationships are ambiguous and species boundaries have not been explicitly tested. Here we use a comprehensive approach to assess species limits through the use of molecular and morphological information (mitochondrial cytb, nuclear sequences collected by ddRADseq, and linear, meristic and landmark-based morphometrics). We found support for the current taxonomy given that the different analyses recognized the nominal species (4 entities), also a candidate species was supported by mitochondrial and morphological data. In addition, we detected signs of admixture between some of the species. Our results indicate that the L. kingii group can serve as a model system in studies of diversification accompanied by hybridization in nature, which in turn might have been promoted by past climatic oscillations and generalist morphologies. We emphasize the importance of using multiple lines of evidence in order to solve evolutionary stories, and minimizing potential erroneous results that may arise when relying on a single source of information.

The enigmatic Leiosaurae clade: Phylogeography, species delimitation, phylogeny and historical biogeography of its southernmost species.

The clade Leiosaurae is composed of poorly-known species endemic to the southern region of South America. The difficulties of finding these lizards in the field, and their highly conserved morphology, have limited our taxonomic knowledge and understanding of their evolutionary histories. Here, we use data collected over 9 years to study the phylogenetic history, genetic diversity, and biogeographic history of almost all the southernmost species of Leiosaurae (except P. nigroigulus), including: Leiosaurus bellii, Diplolaemus darwinii, D. bibronii, D. sexcinctus and D. leopardinus. We use a fragment of the mitochondrial cytochrome-b gene to resolve general phylogeographic patterns, and add another mitochondrial gene and eight nuclear genes to perform species delimitation and phylogenetic analyses associated with divergence times. We found evidence for three putative new species-level taxa within L. bellii and five within Diplolaemus species, indicating high levels of geographic structure. We use a time-calibrated phylogeny to estimate ranges of ancestral distributions and to generate new hypotheses about their historical biogeography.

Phylogenomic data resolve higher-level relationships within South American Liolaemus lizards.

Phylogenomic approaches now generate hundreds of loci representative of the whole genome that can be used for phylogenetic analyses. The South American lizard genus Liolaemus is the most species-rich vertebrate radiation from temperate zones (more than 265 described species), yet most higher-level phylogenetic relationships within Liolaemus remain poorly resolved. In this study, we used 584 nuclear loci collected using targeted sequenced capture to estimate the phylogenetic relationships among 26 species representing the two subgenera within Liolaemus (Eulaemus + Liolaemus), and all major groups within Eulaemus. Previous molecular and morphological-based phylogenetic analyses of Eulaemus based on a limited number of characters resolved few higher-level relationships, although one point of agreement is that the early divergence within Eulaemus corresponds to the lineomaculatus section, followed by the diversification of eight main clades that are strongly supported and recognized. Liolaemus probably experienced relatively rapid divergences during parts of its evolutionary history, and a phylogenomic approach was used to resolve the relationships among the major groups. The new analyses presented here support the division of Liolaemus into two subgenera, and resolve relationships among many of the major clades of Eulaemus with strong support. A Bayesian divergence dating analysis using 44 protein-coding genes provides an estimation of the split of the two Liolaemus subgenera of approximately 19,7 ma (95% HPD = 16,94-23,04), while diversification within Eulaemus started at 15,05 ma (95% HPD = 12,94 - 17,59) among the L. lineomaculatus and the L. montanus series by Mid Miocene. A novel phylogenetic network analyses for SNP data identified two hybridizing edges among different groups of Eulaemus at different points in time. Having a solid phylogenetic hypothesis of the main Eulaemus clades opens new opportunities to test a variety of macroevolutionary questions for this unique radiation.

Coalescent-based species delimitation in the sand lizards of the Liolaemus wiegmannii complex (Squamata: Liolaemidae).

Coalescent-based algorithms coupled with the access to genome-wide data have become powerful tools for assessing questions on recent or rapid diversification, as well as delineating species boundaries in the absence of reciprocal monophyly. In southern South America, the diversification of Liolaemus lizards during the Pleistocene is well documented and has been attributed to the climatic changes that characterized this recent period of time. Past climatic changes had harsh effects at extreme latitudes, including Patagonia, but habitat changes at intermediate latitudes of South America have also been recorded, including expansion of sand fields over northern Patagonia and Pampas). In this work, we apply a coalescent-based approach to study the diversification of the Liolaemus wiegmannii species complex, a morphologically conservative clade that inhabits sandy soils across northwest and south-central Argentina, and the south shores of Uruguay. Using four standard sequence markers (mitochondrial DNA and three nuclear loci) along with ddRADseq data we inferred species limits and a time-calibrated species tree for the L. wiegmannii complex in order to evaluate the influence of Quaternary sand expansion/retraction cycles on diversification. We also evaluated the evolutionary independence of the recently described L. gardeli and inferred its phylogenetic position relative to L. wiegmannii. We find strong evidence for six allopatric candidate species within L. wiegmannii, which diversified during the Pleistocene. The Great Patagonian Glaciation (∼1 million years before present) likely split the species complex into two main groups: one composed of lineages associated with sub-Andean sedimentary formations, and the other mostly related to sand fields in the Pampas and northern Patagonia. We hypothesize that early speciation within L. wiegmannii was influenced by the expansion of sand dunes throughout central Argentina and Pampas. Finally, L. gardeli is supported as a distinct lineage nested within the L. wiegmannii complex.

Evidence of body size and shape stasis driven by selection in Patagonian lizards of the Phymaturus patagonicus clade (Squamata: Liolaemini).

During the speciation process sibling lineages accumulate differences in time (e.g. genetic, morphological, and/or ecological). Phenotypic traits such as size or shape, however, could experience rapid changes or show stasis depending on their role in survival and reproduction. The clade Phymaturus patagonicus includes 26 species characterized by a conservative morphology, and all inhabit rock crevice microhabitats in arid environments. In this study we quantify levels of morphological divergence (size and shape) among the multiple species relative to interspecific molecular divergence, and show that most species have not diverged significantly in size and/or shape to permit unambiguous species diagnosis with morphological data alone. The influence of stabilizing selection for an adaptive optimum in body size and head shape was detected for 13 of the 16 variables analyzed in an Ornstein-Uhlenbeck model. The strict dependence of these species to rock-crevice microenvironments likely explains the observed morphological stasis across the many species of the Phymaturus patagonicus group.

Phylogenomic evidence for a recent and rapid radiation of lizards in the Patagonian Liolaemus fitzingerii species group.

Rapid evolutionary radiations are difficult to resolve because divergence events are nearly synchronous and gene flow among nascent species can be high, resulting in a phylogenetic "bush". Large datasets composed of sequence loci from across the genome can potentially help resolve some of these difficult phylogenetic problems. A suitable test case is the Liolaemus fitzingerii species group of lizards, which includes twelve species that are broadly distributed in Argentinean Patagonia. The species in the group have had a complex evolutionary history that has led to high morphological variation and unstable taxonomy. We generated a sequence capture dataset for 28 ingroup individuals of 580 nuclear loci, alongside a mitogenomic dataset, to infer phylogenetic relationships among species in this group. Relationships among species were generally weakly supported with the nuclear data, and along with an inferred age of ∼2.6 million years old, indicate either rapid evolution, hybridization, incomplete lineage sorting, non-informative data, or a combination thereof. We inferred a signal of mito-nuclear discordance, indicating potential hybridization between L. melanops and L. martorii, and phylogenetic network analyses provided support for 5 reticulation events among species. Phasing the nuclear loci did not provide additional insight into relationships or suspected patterns of hybridization. Only one clade, composed of L. camarones, L. fitzingerii, and L. xanthoviridis was recovered across all analyses. Genomic datasets provide molecular systematists with new opportunities to resolve difficult phylogenetic problems, yet the lack of phylogenetic resolution in Patagonian Liolaemus is biologically meaningful and indicative of a recent and rapid evolutionary radiation. The phylogenetic relationships of the Liolaemus fitzingerii group may be best modeled as a reticulated network instead of a bifurcating phylogeny.

Hybridization could be a common phenomenon within the highly diverse lizard genus Liolaemus.

Hybridization is likely to occur more often between closely related taxa that have had insufficient time to diverge to the point of reproductive incompatibility; hybridization between deeply divergent lineages is rare. In squamate reptiles, hybridization has been proposed as a possible explanation for the extensive paraphyly observed in mitochondrial gene trees in several species complexes of the South American lizard genus Liolaemus. One of the best-documented cases is within the L. boulengeri and L. rothi complexes, which diverged ~5.5 million years ago. Here, we describe a comprehensive study for approaching the hybridization hypothesis between these lizard species complexes. We explored the level of gene tree discordance using the novel 'extra lineage contribution' statistics (XLC, presented in this study) that quantifies the level of gene tree discordance contribution per individual within a species. We included molecular data (12 nuclear and two mitochondrial genes) from 127 individuals, and results of a coalescent model-based analysis show that the most likely explanation for the gene tree-species tree discordance is interspecific hybridization. Our best-supported hypothesis suggests current and past hybridization between L. rothi (rothi complex) and L. tehuelche (boulengeri complex), and independently between L. rothi and L. boulengeri and L. telsen (boulengeri complex). The hybrid descendants are characterized by intermediate phenotypes between the parental species, but are more similar to L. rothi in body size. We discuss the possible role of hybridization in Liolaemus evolution.

Hidden diversity within the lizard genus Liolaemus: Genetic vs morphological divergence in the L. rothi complex (Squamata:Liolaeminae).

Currently, Liolaemus is the second most species-rich reptile genus in the world (257 species), and predictions of its real diversity suggest that it may be the most diverse genus. Originally, Liolaemus species were described as widely distributed and morphologically variable taxa, but extensive sampling in previously unexplored geographic areas, coupled with molecular and more extensive morphological studies, have discovered an unexpectedly high number of previously undetected species. Here, we study the level of molecular vs. morphological divergence within the L. rothi complex, combining a total of 14 loci (2 mitochondrial and 12 nuclear loci) for 97 individuals, as well as morphological data (nine morphometric and 15 color pattern variables), that represent all six described species of the L. rothi complex, plus two candidate species. We use the multi-coalescent species delimitation program iBPP and resolve strong differences in molecular divergence; and each species is inferred as an independent lineage supported by high posterior probabilities. However, morphological differences are not that clear, and our modeling of morphological characters suggests differential selection pressures implying some level of morphological stasis. We discuss the role of natural selection on phenotypic traits, which may be an important factor in "hiding" the real diversity of the genus.

Tracing evolutionary trajectories in the presence of gene flow in South American temperate lizards (Squamata: Liolaemus kingii group).

Evolutionary processes behind lineage divergence often involve multidimensional differentiation. However, in the context of recent divergences, the signals exhibited by each dimension may not converge. In such scenarios, incomplete lineage sorting, gene flow, and scarce phenotypic differentiation are pervasive. Here, we integrated genomic (RAD loci of 90 individuals), phenotypic (linear and geometric traits of 823 and 411 individuals, respectively), spatial, and climatic data to reconstruct the evolutionary history of a speciation continuum of liolaemid lizards (Liolaemus kingii group). Specifically, we (a) inferred the population structure of the group and contrasted it with the phenotypic variability; (b) assessed the role of postdivergence gene flow in shaping phylogeographic and phenotypic patterns; and (c) explored ecogeographic drivers of diversification across time and space. We inferred eight genomic clusters exhibiting leaky genetic borders coincident with geographic transitions. We also found evidence of postdivergence gene flow resulting in transgressive phenotypic evolution in one species. Predicted ancestral niches unveiled suitable areas in southern and eastern Patagonia during glacial and interglacial periods. Our study underscores integrating different data and model-based approaches to determine the underlying causes of diversification, a challenge faced in the study of recently diverged groups. We also highlight Liolaemus as a model system for phylogeographic and broader evolutionary studies.

Quaternary range and demographic expansion of Liolaemus darwinii (Squamata: Liolaemidae) in the Monte Desert of Central Argentina using Bayesian phylogeography and ecological niche modelling.

Until recently, most phylogeographic approaches have been unable to distinguish between demographic and range expansion processes, making it difficult to test for the possibility of range expansion without population growth and vice versa. In this study, we applied a Bayesian phylogeographic approach to reconstruct both demographic and range expansion in the lizard Liolaemus darwinii of the Monte Desert in Central Argentina, during the Late Quaternary. Based on analysis of 14 anonymous nuclear loci and the cytochrome b mitochondrial DNA gene, we detected signals of demographic expansion starting at ~55 ka based on Bayesian Skyline and Skyride Plots. In contrast, Bayesian relaxed models of spatial diffusion suggested that range expansion occurred only between ~95 and 55 ka, and more recently, diffusion rates were very low during demographic expansion. The possibility of population growth without substantial range expansion could account for the shared patterns of demographic expansion during the Last Glacial Maxima (OIS 2 and 4) in fish, small mammals and other lizards of the Monte Desert. We found substantial variation in diffusion rates over time, and very high rates during the range expansion phase, consistent with a rapidly advancing expansion front towards the southeast shown by palaeo-distribution models. Furthermore, the estimated diffusion rates are congruent with observed dispersal rates of lizards in field conditions and therefore provide additional confidence to the temporal scale of inferred phylogeographic patterns. Our study highlights how the integration of phylogeography with palaeo-distribution models can shed light on both demographic and range expansion processes and their potential causes.

A molecular phylogeny of the lizard genus Phymaturus (Squamata, Liolaemini): implications for species diversity and historical biogeography of southern South America.

The lizard genus Phymaturus is widely distributed in Argentina and along the eastern edge of Chile between 25° and 45° south. We sampled 27 of the 38 currently recognized species plus 22 candidate species using two mitochondrial genes (cytb and 12S), four protein coding nuclear genes and seven anonymous nuclear loci, and present the first comprehensive molecular phylogenetic hypothesis for the clade. We recovered two large clades (the palluma or northern group and patagonicus or southern group) previously recognized on the basis of morphological and mitochondrial sequence evidence, and compared results obtained from concatenated-gene analyses with results of a coalescent-based species-tree approach (BEST). With both methods we identified four main clades within the palluma group (mallimaccii, roigorum, verdugo, and vociferator) and five main clades within the patagonicus group (calcogaster, indistinctus, payuniae, somuncurensis, and spurcus). We found several instances of non-monophyly with cytb and cases of incongruence between mitochondrial vs nuclear data for which we discuss alternative hypotheses. Although with lower support values, combined BEST results are more congruent with concatenated nuclear data than with combined concatenated analyses, suggesting that BEST is less influenced by demographic processes than combined concatenated analyses. We discuss the taxonomic, biogeographic and conservation implications of these results and how the future integration of phylogeographic and morphological approaches will allow the further testing of demographic and biogeographic hypotheses.

Accuracy and precision of species trees: effects of locus, individual, and base pair sampling on inference of species trees in lizards of the Liolaemus darwinii group (Squamata, Liolaemidae).

Molecular phylogenetics has entered a new era in which species trees are estimated from a collection of gene trees using methods that accommodate their heterogeneity and discordance with the species tree. Empirical evaluation of species trees is necessary to assess the performance (i.e., accuracy and precision) of these methods with real data, which consists of gene genealogies likely shaped by different historical and demographic processes. We analyzed 20 loci for 16 species of the South American lizards of the Liolaemus darwinii species group and reconstructed a species tree with *BEAST, then compared the performance of this method under different sampling strategies of loci, individuals, and sequence lengths. We found an increase in the accuracy and precision of species trees with the number of loci, but for any number of loci, accuracy substantially decreased only when using only one individual per species or 25% of the full sequence length (∼ 147 bp). In addition, locus "informativeness" was an important factor in the accuracy/precision of species trees when using a few loci, but it became increasingly irrelevant with additional loci. Our empirical results combined with the previous simulation studies suggest that there is an optimal range of sampling effort of loci, individuals, and sequence lengths for a given speciation history and information content of the data. Future studies should be directed toward further assessment of other factors that can impact performance of species trees, including gene flow, locus "informativeness," tree shape, missing data, and errors in species delimitation.

A new lizard species of the Liolaemus kingii group (Squamata: Liolaemidae) from northwestern Chubut province (Argentina).

We describe Liolaemus attenboroughi sp. nov., a lizard distributed in the northwestern Patagonian Steppe of Chubut province (Argentina) previously confused with L. kingii (Bell 1843). Recent studies based on molecular evidence supports its evolutionary independence. Here we provide a morphological diagnosis of this lineage, comparisons between three molecular species delimitation methods, and an updated phylogeny of the L. kingii group. Based on current knowledge of its distribution, this new species is allopatric with geographically close species of the L. kingii group.

How lizards survived blizzards: phylogeography of the Liolaemus lineomaculatus group (Liolaemidae) reveals multiple breaks and refugia in southern Patagonia and their concordance with other codistributed taxa.

Patagonia was shaped by a complex geological history, including the Miocene uplift of the Andes, followed by volcanism, marine introgressions, and extreme climatic oscillations during Pliocene-Pleistocene glaciation-deglaciation cycles. The distributional patterns and phylogenetic relationships of southern patagonian animals and plants were affected in different ways, and those imprints are reflected in the seven phylogeographic breaks and eight refugia that have been previously proposed. In this study, we estimated time-calibrated phylogenetic/phylogeographic patterns in lizards of the Liolaemus lineomaculatus group and relate them to historical Miocene-to-Pleistocene events of Patagonia and the previously proposed phylogeographic patterns. Individuals from 51 localities were sequenced for the mitochondrial marker (cyt-b) and a subsample of individuals from each mitochondrial lineage was sequenced for one nuclear (LDA12D) and one slow evolving mitochondrial gene (12S). Our analyses revealed strong phylogeographic structure among lineages and, in most cases, no signal of demographic changes through time. The lineomaculatus group is composed of three strongly supported clades (lineomaculatus, hatcheri and kolengh + silvanae), and divergence estimates suggested their origins associated with the oldest known Patagonian glaciation (7-5 Ma); subsequent diversification within the lineomaculatus clade coincided with the large Pliocene glaciations (~3.5 Ma). The lineomaculatus clade includes nine strongly genetically and geographically structured lineages, five of which are interpreted as candidate species. Our findings suggest that some Liolaemus lineages have persisted in situ, each of them in a different refugium, through several glaciation-deglaciation cycles without demographic fluctuations. We also summarize and update qualitative evidence of some shared phylogeographic breaks and refugia among plants, rodents and lizards.

Environmental correlates of phenotypic evolution in ecologically diverse Liolaemus lizards.

Evolutionary correlations between phenotypic and environmental traits characterize adaptive radiations. However, the lizard genus Liolaemus, one of the most ecologically diverse terrestrial vertebrate radiations on earth, has so far shown limited or mixed evidence of adaptive diversification in phenotype. Restricted use of comprehensive environmental data, incomplete taxonomic representation and not considering phylogenetic uncertainty may have led to contradictory evidence. We compiled a 26-taxon dataset for the Liolaemus gracilis species group, representing much of the ecological diversity represented within Liolaemus and used environmental data to characterize how environments occupied by species' relate to phenotypic evolution. Our analyses, explicitly accounting for phylogenetic uncertainty, suggest diversification in phenotypic traits toward the present, with body shape evolution rapidly evolving in this group. Body shape evolution correlates with the occupation of different structural habitats indicated by vegetation axes suggesting species have adapted for maximal locomotory performance in these habitats. Our results also imply that the effects of phylogenetic uncertainty and model misspecification may be more extensive on univariate, relative to multivariate analyses of evolutionary correlations, which is an important consideration in analyzing data from rapidly radiating adaptive radiations.

Lizards from the end of the world: phylogenetic relationships of the Liolaemus lineomaculatus section (Squamata: Iguania: Liolaemini).

The Liolaemus lineomaculatus section is a geographically widely distributed group of lizards from the Patagonian region of southern South America, and includes 18 described species representing the most southerly distributed Liolaemus taxa (the genus includes 228 species and extends from Tierra del Fuego north to south-central Peru). Despite high species diversity, the phylogenetic relationships of this section are unknown. In the present work we sampled all described species in the L. lineomaculatus section as well as currently undescribed candidate species to reconstruct the first complete phylogenetic hypothesis for the clade. Our data set included four anonymous nuclear loci, three nuclear protein-coding loci, and two mitochondrial genes. We compared results obtained with three different phylogenetic methods for the concatenated data set (Maximum Parsimony, Maximum Likelihood and Bayesian Inference) with a coalescent-based species tree approach (BEST), and recovered congruent, strongly-supported topological arrangements across all methods. We identified four main clades within the L. lineomaculatus section: the lineomaculatus, magellanicus, somuncurae, and kingii+archeforus groups, for which we estimated divergence times. We discuss the taxonomic implications of these results and how the future integration of phylogeographic, niche modeling and morphological approaches will allow testing biogeographical hypotheses in this clade.