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.

Web building and silk properties functionally covary among species of wolf spider.

Although phylogenetic studies have shown covariation between the properties of spider major ampullate (MA) silk and web building, both spider webs and silks are highly plastic so we cannot be sure whether these traits functionally covary or just vary across environments that the spiders occupy. As MaSp2-like proteins provide MA silk with greater extensibility, their presence is considered necessary for spider webs to effectively capture prey. Wolf spiders (Lycosidae) are predominantly non-web building, but a select few species build webs. We accordingly collected MA silk from two web-building and six non-web-building species found in semirural ecosystems in Uruguay to test whether the presence of MaSp2-like proteins (indicated by amino acid composition, silk mechanical properties and silk nanostructures) was associated with web building across the group. The web-building and non-web-building species were from disparate subfamilies so we estimated a genetic phylogeny to perform appropriate comparisons. For all of the properties measured, we found differences between web-building and non-web-building species. A phylogenetic regression model confirmed that web building and not phylogenetic inertia influences silk properties. Our study definitively showed an ecological influence over spider silk properties. We expect that the presence of the MaSp2-like proteins and the subsequent nanostructures improves the mechanical performance of silks within the webs. Our study furthers our understanding of spider web and silk co-evolution and the ecological implications of spider silk properties.

High genetic diversity but low population structure in the frog Pseudopaludicola falcipes (Hensel, 1867) (Amphibia, Anura) from the Pampas of South America.

Relative to South America's ecoregions, the temperate grasslands of the Pampas have been poorly studied from a phylogeographic perspective. Based on an intermediate biogeographic setting between subtropical forest (Atlantic Forest) and arid ecosystems (Chaco and Patagonia), Pampean species are expected to show unstable demographic histories due to the Quaternary climatic oscillations. Herein, we investigate the phylogenetic relatedness and phylogeographic history of Pseudopaludicola falcipes, a small and common frog that is widely distributed across the Pampean grasslands. First, we use molecular data to assess if P. falcipes represents a single or multiple, separately evolving cryptic lineages. Because P. falcipes is a small-size species (<20mm) with extensive coloration and morphological variation, we suspected that it might represent a complex of cryptic species. In addition, we expected strong genetic and geographic structuring within Pseudopaludicola falcipes due to its large geographic distribution, potentially short dispersal distances, and multiple riverine barriers. We found that P. falcipes is a single evolutionary lineage with poor geographic structuring. Furthermore, current populations of P. falcipes have a large effective population size, maintain ancestral polymorphisms, and have a complex network of gene flow. We conclude that the demographic history of P. falcipes, combined with its ecological attributes and the landscape features of the Pampas, favored a unique combination among anurans of small body size, large population size, high genetic variability, but high cohesiveness of populations over a wide geographic distribution.

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.

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.

PCAtest: testing the statistical significance of Principal Component Analysis in R.

Principal Component Analysis (PCA) is one of the most broadly used statistical methods for the ordination and dimensionality-reduction of multivariate datasets across many scientific disciplines. Trivial PCs can be estimated from data sets without any correlational structure among the original variables, and traditional criteria for selecting non-trivial PC axes are difficult to implement, partially subjective or based on ad hoc thresholds. PCAtest is an R package that implements permutation-based statistical tests to evaluate the overall significance of a PCA, the significance of each PC axis, and of contributions of each observed variable to the significant axes. Based on simulation and empirical results, I encourage R users to routinely apply PCAtest to test the significance of their PCA before proceeding with the direct interpretation of PC axes and/or the utilization of PC scores in subsequent evolutionary and ecological analyses.

Lizards as model organisms for linking phylogeographic and speciation studies.

Lizards have been model organisms for ecological and evolutionary studies from individual to community levels at multiple spatial and temporal scales. Here we highlight lizards as models for phylogeographic studies, review the published population genetics/phylogeography literature to summarize general patterns and trends and describe some studies that have contributed to conceptual advances. Our review includes 426 references and 452 case studies: this literature reflects a general trend of exponential growth associated with the theoretical and empirical expansions of the discipline. We describe recent lizard studies that have contributed to advances in understanding of several aspects of phylogeography, emphasize some linkages between phylogeography and speciation and suggest ways to expand phylogeographic studies to test alternative pattern-based modes of speciation. Allopatric speciation patterns can be tested by phylogeographic approaches if these are designed to discriminate among four alternatives based on the role of selection in driving divergence between populations, including: (i) passive divergence by genetic drift; (ii) adaptive divergence by natural selection (niche conservatism or ecological speciation); and (iii) socially-mediated speciation. Here we propose an expanded approach to compare patterns of variation in phylogeographic data sets that, when coupled with morphological and environmental data, can be used to to discriminate among these alternative speciation patterns. [Correction made after online publication (28/07/2010): (minor deletion in the last line of the abstract)].

Phylogeography of the frog Leptodactylus validus (Amphibia: Anura): patterns and timing of colonization events in the Lesser Antilles.

The frog Leptodactylus validus occurs in northern South America, Trinidad and Tobago, and the southern Lesser Antilles (Grenada and St. Vincent). Mitochondrial DNA sequences were used to perform a nested clade phylogeographic analysis (NCPA), to date colonization events, and to analyze colonization patterns using on a relaxed molecular clock and coalescent simulations. L. validus originated on the mainland and first colonized Trinidad with subsequent independent colonizations of Tobago and the Lesser Antilles from Trinidad. The NCPA suggests a historical vicariant event between populations in Trinidad and Tobago from those in the Lesser Antilles. The colonization of Trinidad occurred approximately 1 million years ago (mya) and the colonization of the Lesser Antillean islands occurred approximately 0.4 mya. The coalescent approach supported the scenario where L. validus dispersed from Trinidad to St. Vincent and from there to Grenada, a dispersal event that could have been mediated by human introduction as recent as 1600 years ago.

Patterns, Mechanisms and Genetics of Speciation in Reptiles and Amphibians.

In this contribution, the aspects of reptile and amphibian speciation that emerged from research performed over the past decade are reviewed. First, this study assesses how patterns and processes of speciation depend on knowing the taxonomy of the group in question, and discuss how integrative taxonomy has contributed to speciation research in these groups. This study then reviews the research on different aspects of speciation in reptiles and amphibians, including biogeography and climatic niches, ecological speciation, the relationship between speciation rates and phenotypic traits, and genetics and genomics. Further, several case studies of speciation in reptiles and amphibians that exemplify many of these themes are discussed. These include studies of integrative taxonomy and biogeography in South American lizards, ecological speciation in European salamanders, speciation and phenotypic evolution in frogs and lizards. The final case study combines genomics and biogeography in tortoises. The field of amphibian and reptile speciation research has steadily moved forward from the assessment of geographic and ecological aspects, to incorporating other dimensions of speciation, such as genetic mechanisms and evolutionary forces. A higher degree of integration among all these dimensions emerges as a goal for future research.

Environmental and seasonal variation in the diet of Elachistocleis bicolor (Guérin-Méneville 1838) (Anura: Microhylidae) from northern Uruguay.

We studied the diet of Elachistocleis bicolor captured in pine, eucalypt, and unmodified environments in Uruguay. Comparisons between seasons (active, inactive) and among three environments (pine, eucalypt, unmodified), were made using rarefaction analyses, importance indices, and non-parametric tests. Elachistocleis bicolor has a specialized diet composed mainly of Pheidole and Solenopsis ants and termites. The diet of E. bicolor includes a high number of prey per individual, suggesting active search as a foraging strategy. The consumption of myrmicine ants (i.e., Solenopsis) represents a source for potential toxic skin secretions that in E. bicolor may be used to avoid being attacked in the ant nests used for shelter during aestivation, as occurs in the microhylid Phrynomantis microps. Diets in forested environments were richer in both periods, which may reflect the colonization of opportunistic ant species in these environments. Ants were more important in the eucalyptus plantations, particularly in the inactive period, whereas termites were more important in the pine plantations mainly in the inactive period. These environmental and seasonal differences in diet are consistent with the temperature and humidity tolerances of ants and termites, which are the main prey of E. bicolor.

Evolutionarily significant units of the critically endangered leaf frog Pithecopus ayeaye (Anura, Phyllomedusidae) are not effectively preserved by the Brazilian protected areas network.

Protected areas (PAs) are essential for biodiversity conservation, but their coverage is considered inefficient for the preservation of all species. Many species are subdivided into evolutionarily significant units (ESUs) and the effectiveness of PAs in protecting them needs to be investigated. We evaluated the usefulness of the Brazilian PAs network in protecting ESUs of the critically endangered Pithecopus ayeaye through ongoing climate change. This species occurs in a threatened mountaintop ecosystem known as campos rupestres. We used multilocus DNA sequences to delimit geographic clusters, which were further validated as ESUs with a coalescent approach. Ecological niche modeling was used to estimate spatial changes in ESUs' potential distributions, and a gap analysis was carried out to evaluate the effectiveness of the Brazilian PAs network to protect P. ayeaye in the face of climate changes. We tested the niche overlap between ESUs to gain insights for potential management alternatives for the species. Pithecopus ayeaye contains at least three ESUs isolated in distinct mountain regions, and one of them is not protected by any PA. There are no climatic niche differences between the units, and only 4% of the suitable potential area of the species is protected in present and future projections. The current PAs are not effective in preserving the intraspecific diversity of P. ayeaye in its present and future range distributions. The genetic structure of P. ayeaye could represent a typical pattern in campos rupestres endemics, which should be considered for evaluating its conservation status.

Allocation trade-off under climate warming in experimental amphibian populations.

Climate change could either directly or indirectly cause population declines via altered temperature, rainfall regimes, food availability or phenological responses. However few studies have focused on allocation trade-offs between growth and reproduction under marginal resources, such as food scarce that may be caused by climate warming. Such critical changes may have an unpredicted impact on amphibian life-history parameters and even population dynamics. Here, we report an allocation strategy of adult anuran individuals involving a reproductive stage under experimental warming. Using outdoor mesocosm experiments we simulated a warming scenario likely to occur at the end of this century. We examined the effects of temperature (ambient vs. pre-/post-hibernation warming) and food availability (normal vs. low) on reproduction and growth parameters of pond frogs (Pelophylax nigromaculatus). We found that temperature was the major factor influencing reproductive time of female pond frogs, which showed a significant advancing under post-hibernation warming treatment. While feeding rate was the major factor influencing reproductive status of females, clutch size, and variation of body size for females, showed significant positive correlations between feeding rate and reproductive status, clutch size, or variation of body size. Our results suggested that reproduction and body size of amphibians might be modulated by climate warming or food availability variation. We believe this study provides some new evidence on allocation strategies suggesting that amphibians could adjust their reproductive output to cope with climate warming.

Species delimitation with ABC and other coalescent-based methods: a test of accuracy with simulations and an empirical example with lizards of the Liolaemus darwinii complex (Squamata: Liolaemidae).

Species delimitation is a major research focus in evolutionary biology because accurate species boundaries are a prerequisite for the study of speciation. New species delimitation methods (SDMs) can accommodate nonmonophyletic species and gene tree discordance as a result of incomplete lineage sorting via the coalescent model, but do not explicitly accommodate gene flow after divergence. Approximate Bayesian computation (ABC) can incorporate gene flow and estimate other relevant parameters of the speciation process while testing alternative species delimitation hypotheses. We evaluated the accuracy of BPP, SpeDeSTEM, and ABC for delimiting species using simulated data and applied these methods to empirical data from lizards of the Liolaemus darwinii complex. Overall, BPP was the most accurate, ABC showed an intermediate accuracy, and SpeDeSTEM was the least accurate under most simulated conditions. All three SDMs showed lower accuracy when speciation occurred despite gene flow, as found in previous studies, but ABC was the method with the smallest decrease in accuracy. All three SDMs consistently supported the distinctness of southern and northern lineages within L. darwinii. These SDMs based on genetic data should be complemented with novel SDMs based on morphological and ecological data to achieve truly integrative and statistically robust approaches to species discovery.