Phylogenetic Analyses of Lizards from the Chilean Humboldt Archipelago Reveal a New Species for the Chañaral Island (Squamata: Liolaemidae)

The Humboldt Archipelago, situated on Chile's north-central coast, is renowned for its exceptional biodiversity. However, lizards of the Liolaemus genus are a particularly understudied group in this archipelago. Liolaemus genus is divided into two clades: chiliensis and nigromaculatus. Within the nigromaculatus clade the zapallarensis group is restricted to the semi-arid and arid coastal habitats of the Atacama Desert in north-central Chile. While it has been reported that lizards from the zapallarensis group inhabit various islands within the Humboldt Archipelago, there has been limited knowledge regarding their specific species identification. To identify the lizard species inhabiting these islands, we conducted phylogenetic analyses using a mitochondrial gene and examined morphological characteristics. Our findings reveal that lizards from the Damas, Choros, and Gaviota islands belong to Liolaemus silvai. In contrast, the lizards on Chañaral Island form a distinct and previously unrecognised group, clearly distinguishable from Liolaemus silvai. In conclusion, our study not only confirms the presence of L. silvai on the Damas, Choros, and Gaviota islands but also describes a new lizard species on Chañaral Island named Liolaemus carezzae sp. nov. These findings contribute valuable insights into the biodiversity of these islands and introduce a newly discovered endemic taxon to the region, enriching our understanding of Chile's unique island ecosystems.

Life Stages and Phylogenetic Position of the New Scale-Mite of the Genus Neopterygosoma (Acariformes: Pterygosomatidae) from Robert's Tree Iguana

A new pterygosomatid mite species, Neopterygosoma robertmertensi sp. n. (Acariformes: Pterygosomatidae) was collected from two specimens of Liolaemus robertmertensi (Liolaemidae) from Argentina. This new species is described based on active stages: adults (female and male) and juveniles (deutonymphs, protonymphs, and larvae) and quiescent stages (nymphchrysalis, deutochrysalis and imagochrysalis). The changes in morphological characters that occur during the ontogeny of N. robertmertensi have been analyzed in detail. A difference in larval sex morphology was observed for the first time in the family Pterygosomatidae (female larvae differ from male larvae in terms of the shape and size of the idiosoma and the position of the genital area). This new mite species is most similar to N. cyanogasteri but can be distinguished by the presence of different leg chaetotaxy patterns of genua IV and femora IV, four to six genital setae, three to five dorsomedial setae, and two to three ventromedial setae. Phylogenetic analysis was conducted based on 120 morphological characters of all Neopterygosoma spp. and four outgroup species using the maximum parsimony approach. The results indicated that this species is nested within mites of the chilensis group of Neopterygosoma associated with host species of the section chiliensis of Liolaemus s. str. An updated diagnosis of the chilensis group of Neopterygosoma and an identification key for all species of this genus has been provided.

Integrating genetics, physiology and morphology to study desert adaptation in a lizard species

Integration of multiple approaches is key to understand the evolutionary processes of local adaptation and speciation. Reptiles have successfully colonized desert environments, that is, extreme and arid conditions that constitute a strong selective pressure on organisms. Here, we studied genomic, physiological and morphological variations of the lizard Liolaemus fuscus to detect adaptations to the Atacama Desert. By comparing populations of L. fuscus inhabiting the Atacama Desert with populations from the Mediterranean forests from central Chile, we aimed at characterizing features related to desert adaptation. We combined ddRAD sequencing with physiological (evaporative water loss, metabolic rate and selected temperature) and morphological (linear and geometric morphometrics) measurements. We integrated the genomic and phenotypic data using redundancy analyses. Results showed strong genetic divergence, along with a high number of fixed loci between desert and forest populations. Analyses detected 110 fixed and 30 outlier loci located within genes, from which 43 were in coding regions, and 12 presented non-synonymous mutations. The candidate genes were associated with cellular membrane and development. Desert lizards presented lower evaporative water loss than those from the forest. Morphological data showed that desert lizards had smaller body size, different allometry, larger eyeballs and more dorsoventrally compressed heads. Our results suggest incipient speciation between desert and forest populations. The adaptive signal must be cautiously interpreted since genetic drift could also contribute to the divergence pattern. Nonetheless, we propose water and resource availability, and changes in habitat structure, as the most relevant challenges for desert reptiles. This study provides insights of the mechanisms that allow speciation as well as desert adaptation in reptiles at multiple levels, and highlights the benefit of integrating independent evidence.

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.

Phylogeographic analysis and species distribution modelling of the wood frog Batrachyla leptopus (Batrachylidae) reveal interglacial diversification in south western Patagonia

Background

The evolutionary history of southern South American organisms has been strongly influenced by Pleistocene climate oscillations. Amphibians are good models to evaluate hypotheses about the influence of these climate cycles on population structure and diversification of the biota, because they are sensitive to environmental changes and have restricted dispersal capabilities. We test hypotheses regarding putative forest refugia and expansion events associated with past climatic changes in the wood frog Batrachyla leptopus distributed along ∼1,000 km of length including glaciated and non-glaciated areas in southwestern Patagonia.

Methods

Using three mitochondrial regions (D-loop, cyt b, and coI) and two nuclear loci (pomc and crybA1), we conducted multilocus phylogeographic analyses and species distribution modelling to gain insights of the evolutionary history of this species. Intraspecific genealogy was explored with maximum likelihood, Bayesian, and phylogenetic network approaches. Diversification time was assessed using molecular clock models in a Bayesian framework, and demographic scenarios were evaluated using approximate Bayesian computation (ABC) and extended Bayesian skyline plot (EBSP). Species distribution models (SDM) were reconstructed using climatic and geographic data.

Results

Population structure and genealogical analyses support the existence of four lineages distributed north to south, with moderate to high phylogenetic support (Bootstrap > 70%; BPP > 0.92). The diversification time of B. leptopus’ populations began at ∼0.107 mya. The divergence between A and B lineages would have occurred by the late Pleistocene, approximately 0.068 mya, and divergence between C and D lineages was approximately 0.065 mya. The ABC simulations indicate that lineages coalesced at two different time periods, suggesting the presence of at least two glacial refugia and a postglacial colonization route that may have generated two southern lineages (p = 0.93, type I error: <0.094, type II error: 0.134). EBSP, mismatch distribution and neutrality indexes suggest sudden population expansion at ∼0.02 mya for all lineages. SDM infers fragmented distributions of B. leptopus associated with Pleistocene glaciations. Although the present populations of B. leptopus are found in zones affected by the last glacial maximum (∼0.023 mya), our analyses recover an older history of interglacial diversification (0.107–0.019 mya). In addition, we hypothesize two glacial refugia and three interglacial colonization routes, one of which gave rise to two expanding lineages in the south.

Connection, isolation and reconnection: Quaternary climatic oscillations and the Andes shaped the phylogeographical patterns of the Patagonian bee Centris cineraria (Apidae)

The joint effect of the Andes as a geographical barrier and the Quaternary glaciations as promoters of genetic divergence remains virtually unexplored in southern South America. To help fill this knowledge gap, in this study we investigated the demographic history of Centris cineraria, a solitary bee mainly distributed in Patagonia. We used mitochondrial and nuclear markers and performed phylogeographical and dating analyses, adjusted spatio-temporal diffusion and species distribution models, and used Approximate Bayesian Computation to identify likely historical demographic scenarios. Our results revealed that during glacial periods the Andes represented a barrier due to the extent of the ice-sheets and the occurrence of unsuitable habitats, while interglacials allowed for gene flow across the Andes. Secondary contact between previously isolated lineages was evident across at least two low-altitude Andean areas, the northern one being a putative glacial refugium. Our findings also suggest that C. cineraria has persisted in situ in four periglacial refugia located along a north–south transect, congruent with the maximum extent of the ice sheet during the Greatest Patagonian Glaciation. As the first phylogeographical study of Patagonian insects, our work reveals that the interaction between Quaternary climatic oscillations and the Andes as a barrier was the main driver of the spatial and demographic history of C. cineraria.

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.

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.

Description and phylogeny of a new species of Liolaemus (Iguania: Liolaemidae) endemic to the south of the Plurinational State of Bolivia

The Liolaemus montanus group is a diverse group of lizards that ranges from central Peru to southwestern Mendoza, Argentina, including much of the Plurinational State of Bolivia ("Bolivia") and Chile. The species of this group mainly inhabit high elevation areas with cold temperatures. In the last years, several species of this group have been described, mostly in Argentina and Chile. In Bolivia, there are at least thirteen valid species belonging to the L. montanus group. In this study, we describe a new species of the L. montanus group with a marked endemism in the Cordillera de Sama of the Tarija Department, Bolivia, and a combination of unique character states that allows its formal description as a new species. The phylogenetic relationships based on analysis of 159 morphological characters suggest that it belongs to the L. montanus group, and that it is closest to Liolaemus pulcherrimus, which is found allopatrically in a small area of the Jujuy Province, Argentina. The multivariate analyses of 66 morphological characters support the phylogenetic relationships. Statistical analyses of inter-species comparisons of morphological characters are not considered the only methods due to the non-independence of some characters states among species; thus, a phylogenetic analysis is recommended. The detailed revision of specimens of the L. montanus group held in the collections of Bolivia is filling major geographic gaps and improving our understanding of the phylogenetic and biogeographic relationships of this widely distributed group of South American lizards.

Unravelling interspecific relationships among highland lizards: first phylogenetic hypothesis using total evidence of the Liolaemus montanus group (Iguania: Liolaemidae)

The South American lizard genus Liolaemus comprises &gt; 260 species, of which &gt; 60 are recognized as members of the Liolaemus montanus group, distributed throughout the Andes in central Peru, Bolivia, Chile and central Argentina. Despite its great morphological diversity and complex taxonomic history, a robust phylogenetic estimate is still lacking for this group. Here, we study the morphological and molecular diversity of the L. montanus group and present the most complete quantitative phylogenetic hypothesis for the group to date. Our phylogeny includes 103 terminal taxa, of which 91 are members of the L. montanus group (58 are assigned to available species and 33 are of uncertain taxonomic status). Our matrix includes 306 morphological and ecological characters and 3057 molecular characters. Morphological characters include 48 continuous and 258 discrete characters, of which 70% (216) are new to the literature. The molecular characters represent five mitochondrial markers. We performed three analyses: a morphology-only matrix, a molecular-only matrix and a matrix including both morphological and molecular characters (total evidence hypothesis). Our total evidence hypothesis recovered the L. montanus group as monophyletic and included ≥ 12 major clades, revealing an unexpectedly complex phylogeny.

But they move! Vicariance and dispersal in southern South America: Using two methods to reconstruct the biogeography of a clade of lizards endemic to South America

This study aims to identify events that modeled the historical biogeography of Phymaturus, using three methodologies: Spatial Analysis of Vicariance (VIP), Statistical Dispersal-Vicariance Analysis (S-DIVA), and Bayesian Binary Method MCMC (BBM). In order to assign areas for the Dispersal-Vicariance and the BBM analyses, we preferred not to use predefined areas, but to identify areas defined via an endemism analysis of Phymaturus species. The analyses were conducted using the same basic topology, which we obtained by constructing a metatree with two recent phylogenies, both morphology and molecular-based. This topology was also used to obtain time divergence estimates in BEAST, using more outgroups than for the metatree in order to get more accurate estimates. The S-DIVA analysis based on the metatree found 25 vicariance events, 20 dispersals and two extinctions; the S-DIVA analysis based on the BEAST tree yielded 30 vicariance events, 42 dispersal events and five extinctions, and the BBM analysis yielded 63 dispersal events, 28 vicariance events and 1 extinction event. According to the metatree analysis, the ancestral area for Phymaturus covers northern Payunia and southern Central Monte. A vicariant event fragmented the ancestral distribution of the genus, resulting in northern Payunia and southern Central Monte as ancestral area for the P. palluma group, and southern Payunia for the P. patagonicus group. The analysis based on the BEAST tree showed a more complex reconstruction, with several dispersal and extinction events in the ancestral node. The Spatial Analysis of Vicariance identified 41 disjunct sister nodes and removed 10 nodes. The barrier that separates the P. palluma group from the P. patagonicus group is roughly congruent with the southern limit of the P. palluma group. The ancestral range for the genus occupies a central position relative to the distribution of the group, which implies that the species must have migrated to the north (P. palluma group) and to the south (P. patagonicus group). To answer questions related to the specific timing of the events, a molecular clock for Phymaturus was obtained, using a Liolaemus fossil for calibration. The present contribution provides a hypothetical framework for the events that modeled the distribution of Phymaturus.

Phylogeny, time divergence, and historical biogeography of the South American Liolaemus alticolor-bibronii group (Iguania: Liolaemidae)

The genus Liolaemus comprises more than 260 species and can be divided in two subgenera: Eulaemus and Liolaemus sensu stricto. In this paper, we present a phylogenetic analysis, divergence times, and ancestral distribution ranges of the Liolaemus alticolor-bibronii group (Liolaemus sensu stricto subgenus). We inferred a total evidence phylogeny combining molecular (Cytb and 12S genes) and morphological characters using Maximum Parsimony and Bayesian Inference. Divergence times were calculated using Bayesian MCMC with an uncorrelated lognormal distributed relaxed clock, calibrated with a fossil record. Ancestral ranges were estimated using the Dispersal-Extinction-Cladogenesis (DEC-Lagrange). Effects of some a priori parameters of DEC were also tested. Distribution ranged from central Perú to southern Argentina, including areas at sea level up to the high Andes. The L. alticolor-bibronii group was recovered as monophyletic, formed by two clades: L. walkeri and L. gracilis, the latter can be split in two groups. Additionally, many species candidates were recognized. We estimate that the L. alticolor-bibronii group diversified 14.5 Myr ago, during the Middle Miocene. Our results suggest that the ancestor of the Liolaemus alticolor-bibronii group was distributed in a wide area including Patagonia and Puna highlands. The speciation pattern follows the South-North Diversification Hypothesis, following the Andean uplift.

Macroevolutionary diversification with limited niche disparity in a species-rich lineage of cold-climate lizards

BACKGROUND

Life diversifies via adaptive radiation when natural selection drives the evolution of ecologically distinct species mediated by their access to novel niche space, or via non-adaptive radiation when new species diversify while retaining ancestral niches. However, while cases of adaptive radiation are widely documented, examples of non-adaptively radiating lineages remain rarely observed. A prolific cold-climate lizard radiation from South America (Phymaturus), sister to a hyper-diverse adaptive radiation (Liolaemus), has extensively diversified phylogenetically and geographically, but with exceptionally minimal ecological and life-history diversification. This lineage, therefore, may offer unique opportunities to investigate the non-adaptive basis of diversification, and in combination with Liolaemus, to cover the whole spectrum of modes of diversification predicted by theory, from adaptive to non-adaptive. Using phylogenetic macroevolutionary modelling performed on a newly created 58-species molecular tree, we establish the tempo and mode of diversification in the Phymaturus radiation.

RESULTS

Lineage accumulation in Phymaturus opposes a density-dependent (or 'niche-filling') process of diversification. Concurrently, we found that body size diversification is better described by an Ornstein-Uhlenbeck evolutionary model, suggesting stabilizing selection as the mechanism underlying niche conservatism (i.e., maintaining two fundamental size peaks), and which has predominantly evolved around two major adaptive peaks on a 'Simpsonian' adaptive landscape.

CONCLUSIONS

Lineage diversification of the Phymaturus genus does not conform to an adaptive radiation, as it is characterised by a constant rate of species accumulation during the clade's history. Their strict habitat requirements (rocky outcrops), predominantly invariant herbivory, and especially the constant viviparous reproduction across species have likely limited their opportunities for adaptive diversifications throughout novel environments. This mode of diversification contrasts dramatically with its sister lineage Liolaemus, which geographically overlaps with Phymaturus, but exploits all possible microhabitats in these and other bioclimatic areas. Our study contributes importantly to consolidate these lizards (liolaemids) as promising model systems to investigate the entire spectrum of modes of species formations, from the adaptive to the non-adaptive extremes of the continuum.

Phylogenomic analysis of the Chilean clade of Liolaemus lizards (Squamata: Liolaemidae) based on sequence capture data

The genus Liolaemus is one of the most ecologically diverse and species-rich genera of lizards worldwide. It currently includes more than 250 recognized species, which have been subject to many ecological and evolutionary studies. Nevertheless, Liolaemus lizards have a complex taxonomic history, mainly due to the incongruence between morphological and genetic data, incomplete taxon sampling, incomplete lineage sorting and hybridization. In addition, as many species have restricted and remote distributions, this has hampered their examination and inclusion in molecular systematic studies. The aims of this study are to infer a robust phylogeny for a subsample of lizards representing the Chilean clade (subgenus Liolaemus sensu stricto), and to test the monophyly of several of the major species groups. We use a phylogenomic approach, targeting 541 ultra-conserved elements (UCEs) and 44 protein-coding genes for 16 taxa. We conduct a comparison of phylogenetic analyses using maximum-likelihood and several species tree inference methods. The UCEs provide stronger support for phylogenetic relationships compared to the protein-coding genes; however, the UCEs outnumber the protein-coding genes by 10-fold. On average, the protein-coding genes contain over twice the number of informative sites. Based on our phylogenomic analyses, all the groups sampled are polyphyletic. Liolaemus tenuis tenuis is difficult to place in the phylogeny, because only a few loci (nine) were recovered for this species. Topologies or support values did not change dramatically upon exclusion of L. t. tenuis from analyses, suggesting that missing data did not had a significant impact on phylogenetic inference in this data set. The phylogenomic analyses provide strong support for sister group relationships between L. fuscus, L. monticola, L. nigroviridis and L. nitidus, and L. platei and L. velosoi. Despite our limited taxon sampling, we have provided a reliable starting hypothesis for the relationships among many major groups of the Chilean clade of Liolaemus that will help future work aimed at resolving the Liolaemus phylogeny.

What defines an adaptive radiation? Macroevolutionary diversification dynamics of an exceptionally species-rich continental lizard radiation

BACKGROUND

Adaptive radiation theory posits that ecological opportunity promotes rapid proliferation of phylogenetic and ecological diversity. Given that adaptive radiation proceeds via occupation of available niche space in newly accessed ecological zones, theory predicts that: (i) evolutionary diversification follows an 'early-burst' process, i.e., it accelerates early in the history of a clade (when available niche space facilitates speciation), and subsequently slows down as niche space becomes saturated by new species; and (ii) phylogenetic branching is accompanied by diversification of ecologically relevant phenotypic traits among newly evolving species. Here, we employ macroevolutionary phylogenetic model-selection analyses to address these two predictions about evolutionary diversification using one of the most exceptionally species-rich and ecologically diverse lineages of living vertebrates, the South American lizard genus Liolaemus.

RESULTS

Our phylogenetic analyses lend support to a density-dependent lineage diversification model. However, the lineage through-time diversification curve does not provide strong support for an early burst. In contrast, the evolution of phenotypic (body size) relative disparity is high, significantly different from a Brownian model during approximately the last 5 million years of Liolaemus evolution. Model-fitting analyses also reject the 'early-burst' model of phenotypic evolution, and instead favour stabilizing selection (Ornstein-Uhlenbeck, with three peaks identified) as the best model for body size diversification. Finally, diversification rates tend to increase with smaller body size.

CONCLUSIONS

Liolaemus have diversified under a density-dependent process with slightly pronounced apparent episodic pulses of lineage accumulation, which are compatible with the expected episodic ecological opportunity created by gradual uplifts of the Andes over the last ~25My. We argue that ecological opportunity can be strong and a crucial driver of adaptive radiations in continents, but may emerge less frequently (compared to islands) when major events (e.g., climatic, geographic) significantly modify environments. In contrast, body size diversification conforms to an Ornstein-Uhlenbeck model with multiple trait optima. Despite this asymmetric diversification between both lineages and phenotype, links are expected to exist between the two processes, as shown by our trait-dependent analyses of diversification. We finally suggest that the definition of adaptive radiation should not be conditioned by the existence of early-bursts of diversification, and should instead be generalized to lineages in which species and ecological diversity have evolved from a single ancestor.

Two new species of the Liolaemuselongatus-kriegi complex (Iguania, Liolaemidae) from Andean highlands of southern Chile

The elongatus-kriegi complex is one of the most diverse clades of the Liolaemus (sensu stricto) subgenus of lizards. There are currently 29 species recognized in this group distributed between Chile and Argentina. Based on molecular evidence, there seem to be five main clades nested within this complex: the elongatus, leopardinus, kriegi, petrophilus and punmahuida clades. Liolaemusbuergeri and Liolaemuskriegi, both of the kriegi clade, were believed to inhabit the surroundings of the Laja Lagoon, in the Biobío Region of Chile. Moreover, this Chilean population of Liolaemuskriegi was recently recognized as an undescribed taxon called “Liolaemus sp. A” based on molecular phylogenetics. In this work, we studied these two populations of the Laja Lagoon and provided the morphological diagnosis to describe them as two new species: Liolaemusscorialissp. n. and Liolaemuszabalaisp. n., previously considered Liolaemusbuergeri and “Liolaemuskriegi/Liolaemus sp. A” respectively. Additionally, we identified another population of Liolaemusscorialis in the vicinity of La Mula Lagoon in the Araucanía Region of Chile. Liolaemusscorialis differs from almost all of the species of the elongatus-kriegi complex by its considerably smaller size. Nevertheless, without molecular data we cannot assign it to any particular subclade. Liolaemuszabalai belongs to the kriegi clade based on published molecular phylogenies. Finally, we provide some natural history data on both species and we document for the first time the presence of Liolaemusneuquensis in Chile from a museum specimen from La Mula Lagoon.

Influence of gene flow on divergence dating - implications for the speciation history of Takydromus grass lizards

Dating the time of divergence and understanding speciation processes are central to the study of the evolutionary history of organisms but are notoriously difficult. The difficulty is largely rooted in variations in the ancestral population size or in the genealogy variation across loci. To depict the speciation processes and divergence histories of three monophyletic Takydromus species endemic to Taiwan, we sequenced 20 nuclear loci and combined with one mitochondrial locus published in GenBank. They were analysed by a multispecies coalescent approach within a Bayesian framework. Divergence dating based on the gene tree approach showed high variation among loci, and the divergence was estimated at an earlier date than when derived by the species-tree approach. To test whether variations in the ancestral population size accounted for the majority of this variation, we conducted computer inferences using isolation-with-migration (IM) and approximate Bayesian computation (ABC) frameworks. The results revealed that gene flow during the early stage of speciation was strongly favoured over the isolation model, and the initiation of the speciation process was far earlier than the dates estimated by gene- and species-based divergence dating. Due to their limited dispersal ability, it is suggested that geographical isolation may have played a major role in the divergence of these Takydromus species. Nevertheless, this study reveals a more complex situation and demonstrates that gene flow during the speciation process cannot be overlooked and may have a great impact on divergence dating. By using multilocus data and incorporating Bayesian coalescence approaches, we provide a more biologically realistic framework for delineating the divergence history of Takydromus.

Integrative taxonomy and preliminary assessment of species limits in the Liolaemus walkeri complex (Squamata, Liolaemidae) with descriptions of three new species from Peru

Species delimitation studies based on integrative taxonomic approaches have received considerable attention in the last few years, and have provided the strongest hypotheses of species boundaries. We used three lines of evidence (molecular, morphological, and niche envelopes) to test for species boundaries in Peruvian populations of the Liolaemus walkeri complex. Our results show that different lines of evidence and analyses are congruent in different combinations, for unambiguous delimitation of three lineages that were "hidden" within known species, and now deserve species status. Our phylogenetic analysis shows that L. walkeri, L. tacnae and the three new species are strongly separated from other species assigned to the alticolor-bibronii group. Few conventional morphological characters distinguish the new species from closely related taxa and this highlights the need to integrate other sources of data to erect strong hypothesis of species limits. A taxonomic key for known Peruvian species of the subgenus Lioalemus is provided.

Undersampling taxa will underestimate molecular divergence dates: an example from the South american lizard clade liolaemini

Methods for estimating divergence times from molecular data have improved dramatically over the past decade, yet there are few studies examining alternative taxon sampling effects on node age estimates. Here, I investigate the effect of undersampling species diversity on node ages of the South American lizard clade Liolaemini using several alternative subsampling strategies for both time calibrations and taxa numbers. Penalized likelihood (PL) and Bayesian molecular dating analyses were conducted on a densely sampled (202 taxa) mtDNA-based phylogenetic hypothesis of Iguanidae, including 92 Liolaemini species. Using all calibrations and penalized likelihood, clades with very low taxon sampling had node age estimates younger than clades with more complete taxon sampling. The effect of Bayesian and PL methods differed when either one or two calibrations only were used with dense taxon sampling. Bayesian node ages were always older when fewer calibrations were used, whereas PL node ages were always younger. This work reinforces two important points: (1) whenever possible, authors should strongly consider adding as many taxa as possible, including numerous outgroups, prior to node age estimation to avoid considerable node age underestimation and (2) using more, critically assessed, and accurate fossil calibrations should yield improved divergence time estimates.

Global taxonomic diversity of living reptiles

Reptiles are one of the most ecologically and evolutionarily remarkable groups of living organisms, having successfully colonized most of the planet, including the oceans and some of the harshest and more environmentally unstable ecosystems on earth. Here, based on a complete dataset of all the world's diversity of living reptiles, we analyse lineage taxonomic richness both within and among clades, at different levels of the phylogenetic hierarchy. We also analyse the historical tendencies in the descriptions of new reptile species from Linnaeus to March 2012. Although (non-avian) reptiles are the second most species-rich group of amniotes after birds, most of their diversity (96.3%) is concentrated in squamates (59% lizards, 35% snakes, and 2% amphisbaenians). In strong contrast, turtles (3.4%), crocodilians (0.3%), and tuataras (0.01%) are far less diverse. In terms of species discoveries, most turtles and crocodilians were described early, while descriptions of lizards, snakes and amphisbaenians are multimodal with respect to time. Lizard descriptions, in particular, have reached unprecedented levels during the last decade. Finally, despite such remarkably asymmetric distributions of reptile taxonomic diversity among groups, we found that the distributions of lineage richness are consistently right-skewed, with most clades (monophyletic families and genera) containing few lineages (monophyletic genera and species, respectively), while only a few have radiated greatly (notably the families Colubridae and Scincidae, and the lizard genera Anolis and Liolaemus). Therefore, such consistency in the frequency distribution of richness among clades and among phylogenetic levels suggests that the nature of reptile biodiversity is fundamentally fractal (i.e., it is scale invariant). We then compared current reptile diversity with the global reptile diversity and taxonomy known in 1980. Despite substantial differences in the taxonomies (relative to 2012), the patterns of lineage richness remain qualitatively identical, hence reinforcing our conclusions about the fractal nature of reptile biodiversity.

Lizards on ice: evidence for multiple refugia in Liolaemus pictus (Liolaemidae) during the last glacial maximum in the Southern Andean beech forests

Historical climate changes and orogenesis are two important factors that have shaped intraspecific biodiversity patterns worldwide. Although southern South America has experienced such complex events, there is a paucity of studies examining the effects on intraspecific diversification in this part of the world. Liolaemus pictus is the southernmost distributed lizard in the Chilean temperate forest, whose genetic structure has likely been influenced by Pleistocene glaciations. We conducted a phylogeographic study of L. pictus in Chile and Argentina based on one mitochondrial and two nuclear genes recovering two strongly divergent groups, Northern and Southern clades. The first group is distributed from the northernmost limit of the species to the Araucanía region while the second group is distributed throughout the Andes and the Chiloé archipelago in Southern Chile. Our results suggest that L. pictus originated 751 Kya, with divergence between the two clades occurring in the late Pleistocene. Demographic reconstructions for the Northern and Southern clades indicate a decrease in effective population sizes likely associated with Pleistocene glaciations. Surprisingly, patterns of genetic variation, clades age and historical gene flow in populations distributed within the limits of the Last Glacial Maximum (LGM) are not explained by recent colonization. We propose an “intra-Andean multiple refuge” hypothesis, along with the classical refuge hypothesis previously proposed for the biota of the Chilean Coastal range and Eastern Andean Cordillera. Our hypothesis is supported by niche modelling analysis suggesting the persistence of fragments of suitable habitat for the species within the limits of the LGM ice shield. This type of refuge hypothesis is proposed for the first time for an ectothermic species.