Amphibian phylogeography in the Antipodes: Refugia and postglacial colonization explain mitochondrial haplotype distribution in the Patagonian frog Eupsophus calcaratus (Cycloramphidae)

Local adaptation of Dromiciops marsupials (Microbiotheriidae) from southern South America: Implications for species management facing climate change

The two species of the microbiotheriid marsupial genus Dromiciops (Dromiciops bozinovici: "Panchos's monito del monte" and Dromiciops gliroides: "monito del monte") exhibit a marked latitudinal genetic differentiation. Nevertheless, it is unclear whether this differentiation results from neutral processes or can be explained, to some extent, by local adaptation to different environmental conditions. Here, we used an SNP panel gathered by Rad-seq and searched for footprints of local adaptation (putative loci under selection) by exploring genetic associations with environmental variables in the two species of Dromiciops in Chilean and Argentinean populations. We applied three methods for detecting outlier SNPs and two genotype-environment associations approaches to quantify associations between allelic frequencies and environmental variables. Both species display strong genetic structure. D. bozinovici exhibited three distinct genetic groups, marking the first report of such structuring in this species using SNPs. In contrast, D. gliroides displayed four genetic clusters, consistent with previous studies. Both species exhibited an association of their genetic structure with environmental variables. D. bozinovici exhibited significant associations of allelic frequencies with elevation, precipitation during the warmest periods, and seasonality in the thermal regime. For D. gliroides, genetic variation appeared to be associated with more variables than D. bozinovici, including precipitation and temperature-related variables, isothermality, and elevation. All the outlier SNPs were mapped to the D. gliroides reference genome to explore if they fell within functionally known genes. These results represent a necessary first step toward identifying the genome regions that harbor genes associated with climate adaptations in Dromiciops. Notably, we identified genes involved in various functions, including carbohydrate synthesis (ALG8), muscle and neuronal regulation (MEF2D), and stress responses (PTGES3). Ultimately, this study contributes valuable insights that can inform targeted conservation strategies aimed at preserving the genetic diversity of Dromiciops in the face of environmental challenges.

Genomic insights into the Montseny brook newt (Calotriton arnoldi), a Critically Endangered glacial relict

The Montseny brook newt (Calotriton arnoldi), considered the most endangered amphibian in Europe, is a relict salamandrid species endemic to a small massif located in northeastern Spain. Although conservation efforts should always be guided by genomic studies, those are yet scarce among urodeles, hampered by the extreme sizes of their genomes. Here, we present the third available genome assembly for the order Caudata, and the first genomic study of the species and its sister taxon, the Pyrenean brook newt (Calotriton asper), combining whole-genome and ddRADseq data. Our results reveal significant demographic oscillations which accurately mirrored Europe's climatic history. Although severe bottlenecks have led to depauperate genomic diversity and long runs of homozygosity along a gigantic genome, inbreeding might have been avoided by assortative mating strategies. Other life history traits, however, seem to have been less advantageous, and the lack of land dispersal has driven to exceptional levels of population fragmentation.

Easternmost distribution of Bufo bufo (Linnaeus, 1758) in Türkiye: implications for the putative contact zone between B. bufo and B. verrucosissimus

The geographic range of a species is crucial for obtaining information on the exact distribution of the species. The geographic data are important for delimiting distinct species or exploring the degree of differentiation among different populations of a species. The local details of species boundaries facilitate the study of the importance of phylogeographic background, secondary contacts, and hybrid zones, along with the relations between the species and its extrinsic environmental factors. In the present study, the range boundaries of Bufo bufo and Bufo verrucosissimus in the north-eastern region of Türkiye were delineated using an integrative taxonomic approach that utilized a combination of molecular and morphological data. According to the mtDNA results of the present study, B. bufo inhabits a single distribution from İyidere town to Çayeli town in Rize, while B. verrucosissimus is distributed from Şavşat town of Artvin to Ardeşen town in Rize. In addition, the two species coexist in Pazar, Hemşin, and Çamlıhemşin towns in Rize. The demographic analyses indicated a distinct population expansion for the B. verrucosissimus species after the Last Glacial Maximum, while the same did not occur for B. bufo. The univariate and multivariate statistical analyses conducted for the morphological data of the two species corroborated the presence of a putative contact zone between B. bufo and B. verrucosissimus. In summary, the present study resolved the non-distinct geographic boundaries between B. bufo and B. verrucosissimus species and also revealed the easternmost distribution of B. bufo in Türkiye. In addition, important evidence on the putative contact zone between the two species was indicated using an integrative taxonomic approach.

ddRAD-Seq reveals evolutionary insights into population differentiation and the cryptic phylogeography of Hyporhamphus intermedius in Mainland China

Species differentiation and local adaptation in heterogeneous environments have attracted much attention, although little is known about the mechanisms involved. Hyporhamphus intermedius is an anadromous, brackish-water halfbeak that is widely distributed in coastal areas and hyperdiverse freshwater systems in China, making it an interesting model for research on phylogeography and local adaptation. Here, 156 individuals were sampled at eight sites from heterogeneous aquatic habitats to examine environmental and genetic contributions to phenotypic divergence. Using double-digest restriction-site-associated DNA sequencing (ddRAD-Seq) in the specimens from the different watersheds, 5498 single nucleotide polymorphisms (SNPs) were found among populations, with obvious population differentiation. We find that present-day Mainland China populations are structured into distinct genetic clusters stretching from southern and northern ancestries, mirroring geography. Following a transplant event in Plateau Lakes, there were virtually no variations of genetic diversity occurred in two populations, despite the fact two main splits were unveiled in the demographic history. Additionally, dorsal, and anal fin traits varied widely between the southern group and the others, which highlighted previously unrecognized lineages. We then explore genotype-phenotype-environment associations and predict candidate loci. Subgroup ranges appeared to correspond to geographic regions with heterogeneous hydrological factors, indicating that these features are likely important drivers of diversification. Accordingly, we conclude that genetic and phenotypic polymorphism and a moderate amount of genetic differentiation occurred, which might be ascribed to population subdivision, and the impact of abiotic factors.

Molecular data reveal multiple lineages of Scinax nebulosus (Spix, 1824) (Anura: Hylidae) with Plio-Pleistocene diversification in different Brazilian regions

To understand the organism's history, we can start assessing the complexity of the biome where they occur. In this study, we used a region of the mitochondrial genome, the rRNA 16S, to evaluate the genetic differentiation in Scinax nebulosus along with its geographical range highlighting important Brazilian biomes as Restinga, Cerrado, Amazon, and Atlantic Forest. Geographically structured genetic divergence was observed within the species S. nebulosus. The values of the fixation index (Фst) and the pairwise Fst index were high and significant regarding this structuring. Besides, the haplotype network corroborates these results with the haplotypes arrangement found by separating the S. nebulosus populations in two major groups: North and Northeast. The lineage delimitation analyses indicate the occurrence of several lineages with divergence mainly between the samples from the Northeast group. Thus, we can suggest that S. nebulosus may present itself as a group of cryptic species due to the genetic characteristics found. The existence of a mosaic of heterogeneous habitats may explain the genetic divergence found, which justifies the existence of cryptic species in this group. However, this hypothesis needs more detail in molecular studies, including large sample sizes and other population and demographic analyses.

Genomic diversity and Demographic History of the Dromiciops genus (Marsupialia: Microbiotheriidae)

Three orders represent the South American fauna of marsupials. Of these, Microbiotheria was until recently known as a monotypic genus with the only surviving species Dromiciops gliroides (monito del monte). The recent proposal of a new Dromiciops species (Dromiciops bicinovici), together with new information on the origin and diversification of living microbioterians has changed the prevailing paradigm around the evolutionary history of these emblematic marsupials. Here, we used a RADseq approach to test for evidence of admixture and past or current gene flow among both species of Dromiciops and evaluate the genetic structure within D. gliroides. We analyzed 127 samples of Dromiciops distributed across the known distribution range of both species. We also inferred the joint demographic history of these lineages, thus corroborating the status of D. bozinovici as a distinct species. Demographic history reconstruction indicated that D. bozinovici diverged from D. gliroides around 4my ago and has remained isolated and demographically stable ever since. In contrast, D. gliroides is subdivided into three subclades that experienced recent expansions and moderate gene flow among them (mostly from north to south). Furthermore, genetic distances among populations within D. gliroides were significantly correlated with geographic distances. These results suggest that some of the D. gliroides populations would have survived in glacial refuges, with posterior expansions after ice retreat. Our results have important implications for the systematics of the genus and have profound conservation consequences for the new species, especially considering the fragmentation level of the temperate rainforest.

Assessing signals of selection and historical demography to develop conservation strategies in the Chilean emblematic Araucaria araucana

Loss of genetic diversity reduces the ability of species to evolve and respond to environmental change. Araucaria araucana is an emblematic conifer species from southern South America, with important ethnic value for the Mapuche people (Pehuenche); the Chilean Government has catalogued its conservation status as vulnerable. Climatic fluctuations were potentially a major impact in the genetic variation within many tree species. In this context, the restricted geographic distribution of A. araucana in Chile appears to be a consequence of the Last Glacial Maximum (LGM). During the past two centuries, strong human intervention has also affected the geographical distribution and population sizes of A. araucana. Reduction of population size may cause loss of genetic diversity, which could affect frequency of adaptive loci. The aims of this study were to know the existence of potential loci under selection and populations with genetic, demographic disequilibrium in the Chilean distribution of A. araucana. Based on 268 polymorphic AFLP loci, we have investigated potential loci under selection and genetic, demographic disequilibrium within seven Chilean populations of Araucaria araucana. Correlation of 41 outlier loci with the environmental variables of precipitation and temperature reveals signatures of selection, whereas 227 neutral loci provide estimates of demographic equilibrium and genetic population structure. Three populations are recommended as priorities for conservation.

Phylogeography and molecular species delimitation reveal cryptic diversity in Potamolithus (Caenogastropoda: Tateidae) of the southwest basin of the Andes

The species of the genus Potamolithus inhabiting the southwestern basin of the Andes are difficult to distinguish due to small size and similar shell morphology. Only Potamolithus australis and Potamolithus santiagensis have been traditionally recognized in this region, but the occurrence of several morphologically similar undescribed populations could increase the regional richness. Here we delimit described and potentially undescribed cryptic species of the genus using partial sequences of the mitochondrial cytochrome c oxidase subunit I (COI) gene. Network analysis and diversity indices inferred six highly differentiated haplogroups, many of them sympatric and widespread in the study area. Phylogeographic analyses suggest a scenario of recent diversification and the occurrence of multiple refuges during the successive Pleistocene glaciations. Phylogenetic analysis also recovered six major clades that showed no relationship with physiography. Species delimitation analyses consistently recognized three or four candidate species apart from P. australis and P. santiagensis. Divergence times indicate that speciation of Chilean Potamolithus began at the end of the Pliocene, probably driven by climatic rather than geographic events. Considering the high inter- and intra-basin genetic diversity, conservation efforts should be focused on protecting sympatric taxa in the basins with the highest species richness.

Genetic Diversity and Demographic History of the Shaggy Soft-Haired Mouse Abrothrix hirta (Cricetidae; Abrotrichini)

Genetic information on species can inform decision making regarding conservation of biodiversity since the response of organisms to changing environments depend, in part, on their genetic makeup. Territories of central-southern Chile and Argentina have undergone a varying degree of impact during the Quaternary, where the response of local fauna and flora was rather species-specific. Here, we focus on the sigmodontine rodent Abrothrix hirta, distributed from 35° S in Chile and Argentina to northern Tierra del Fuego. Based on 119,226 transcriptome-derived SNP loci from 46 individuals of A. hirta, we described the geographic distribution of the genetic diversity of this species using a maximum likelihood tree, principal component and admixture analyses. We also addressed the demographic history of the main intraspecific lineages of A. hirta using GADMA. We found that A. hirta exhibited four allopatric intraspecific lineages. Three main genetic groups were identified by a Principal Component Analysis and by Ancestry analysis. The demographic history of A. hirta was characterized by recent population stability for populations at the northernmost part of the range, while southern populations experienced a recent population expansion.

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.

Contrasting evolutionary responses in two co-distributed species of Galaxias (Pisces, Galaxiidae) in a river from the glaciated range in Southern Chile

Life-history traits are among the most important factors affecting population abundance and genetic diversity of species. Here, we analysed the genetic patterns of two Galaxias species with different life-history traits to investigate how these biological differences impacted their evolution in the Valdivia River basin, Southern Chile. We analysed mitochondrial DNA (mtDNA) sequences from 225 individuals of Galaxias maculatus and 136 of G. platei to compare patterns of genetic diversity, structure and demographic growth across the basin. Galaxias maculatus presented higher genetic diversity and higher genetic structure than G. platei. Demographic analyses showed G. maculatus kept a higher population size over time, with a signal of demographic expansion in the last 250 kyr. Whereas Galaxias platei, exhibited lower, but constant population size over time. Furthermore, haplotype networks revealed higher lineage diversity in G. maculatus with a tendency to occupy different areas of the basin. Coalescent simulations ruled out that genetic differences between species could be explained by stochastic processes (genetic drift), suggesting species-specific biological differences as responsible for the observed genetic differences. We discuss how differences in life-history traits and past glaciations interact to shape the evolutionary history of the two Galaxias species.

Sufficient versus optimal climatic stability during the Late Quaternary: using environmental quality to guide phylogeographic inferences in a Neotropical montane system

Quaternary climatic oscillations affected species distributions worldwide, creating cycles of connectivity and isolation that impacted population demography and promoted lineage divergence. These effects have been well studied in temperate regions. Taxa inhabiting mesic montane habitats in tropical ecosystems show high levels of endemism and diversification in the distinct mountain ranges they inhabit; such a pattern has commonly been ascribed to past climatic oscillations, but few phylogeographic studies have tested this hypothesis. Here, we combine ecological niche models of species distributions with molecular data to study phylogeographic patterns in two rodents endemic to the highlands of Costa Rica and western Panama (Reithrodontomys creper and Nephelomys devius). In so doing, we apply a novel approach that incorporates a basic ecological principle: the expected positive relationship between environmental suitability and population abundance. Specifically, we use niche models to predict potential patterns of population connectivity and stability of different suitability levels during climatic extremes of the last glacial–interglacial cycle; we then test these predictions with population genetic analyses of a mitochondrial and a nuclear marker. The detailed predictions arising from the different levels of suitability were moderately to highly congruent with the molecular data depending on the species. Overall, results suggest that in these tropical montane ecosystems, cycles of population connectivity and isolation followed a pattern opposite to that typically described for temperate or lowland tropical ecosystems: namely, higher connectivity during the colder glacials, with isolation in montane refugia during the interglacials, including today. Nevertheless, the individualistic patterns for each species indicate a potentially wide gamut of phylogeographic histories reflecting particularities of their niches. Taken together, this study illustrates how phylogeographic inferences may benefit from niche model outputs that provide more detailed predictions of connectivity and finer characterizations of potential refugia through time.

Taxonomy, systematics and geographic distribution of ground frogs (Alsodidae, Eupsophus): a comprehensive synthesis of the last six decades of research

The genus Eupsophus (ground frogs) inhabits exclusively the temperate forests of southern South America (Chile and Argentina). The current delimitation of the genus was reached in the late 1970s, when only two species were recognized, but since then the number of described species steadily increased, reaching a maximum of 11 by 2012. Subsequent studies that applied explicit species delimitation approaches decreased the number of species to six in 2017 and raised it again to 11 the following year, including an undescribed putative species. Despite these taxonomic changes, the two species groups traditionally recognized, roseus and vertebralis, have been maintained. Another recent contribution to the taxonomy of the genus was the explicit recognition of the extremely high level of external phenotypic variation exhibited by species of the roseus group, which undermines the utility of some diagnostic characters. Here we provide a critical review of the extensive taxonomic and systematic literature on the genus over the last six decades, to examine the evidence behind the recurrent taxonomic changes and advances in its systematics. We also update and complete a 2017 review of geographic information, provide additional qualitative observations of external characters commonly used in the diagnoses of species of the roseus group, and reassess the phylogenetic position of a putative new species from Tolhuaca (Chile), which was not included in the last species delimitation study. The present review shows that: 1) there is no congruence between the patterns of phenotypic and genetic/phylogenetic differentiation among species of both groups; 2) in the roseus group, the intraspecific variation in some external characters is as high as the differences described among species; 3) there is little morphological and bioacoustic differentiation within species groups, and inconsistencies in the chromosomal evidence at the genus level; 4) under the latest taxonomic proposal (2018), species of the roseus group still lack consistent and reliable diagnoses and their distribution limits are poorly defined; and 5) the population from Tolhuaca represents an additional undescribed species under the most recent taxonomic framework. Finally, we discuss the implications of these findings for the taxonomy and biogeography of the genus, pointing out some areas that require further research to understand their patterns and processes of diversification.

Local persistence of Mann's soft-haired mouse Abrothrix manni (Rodentia, Sigmodontinae) during Quaternary glaciations in southern Chile

Quaternary climatic oscillations have impacted Patagonian sigmodontine fauna, leaving traceable genetic footprints. In southern Chile, changes in the landscape included transitions to different vegetation formations as well as the extension of ice sheets. In this study, we focus on the Valdivian forest endemic and recently described sigmodontine species Abrothrix manni. We aim to assess the genetic structure of this species, testing for the existence of intraspecific lineages, and inferring the recent demographic history of the species. Analyses were based on the first 801 bp of the mitochondrial gene Cytocrhome-b from 49 individuals of A. manni collected at 10 localities that covers most part of its geographic distribution. Genealogical analyses recovered two main intraspecific lineages that are geographically segregated and present an intermediate site of secondary contact. Historical demography shows signal of recent population decrease. Based on these results, we proposed that current genetic diversity of A. manni differentiated in at least two distinct refugial areas in southern Chile. This scenario, in addition to be unique among those uncovered for the so far studied Valdivian forest rodents, is noteworthy because of the reduced geographic scale inhabited by the species.

Speciation in a biodiversity hotspot: Phylogenetic relationships, species delimitation, and divergence times of Patagonian ground frogs from the Eupsophus roseus group (Alsodidae)

The alsodid ground frogs of the Eupsophus genus are divided into two groups, the roseus (2n = 30) and vertebralis (2n = 28), which are distributed throughout the temperate Nothofagus forests of South America. Currently, the roseus group is composed by four species, while the vertebralis group consists of two. Phylogenetic relationships and species delimitation within each group are controversial. In fact, previous analyses considered that the roseus group was composed of between four to nine species. In this work, we evaluated phylogenetic relationships, diversification times, and species delimitation within the roseus group using a multi-locus dataset. For this purpose, mitochondrial (D-loop, Cyt b, and COI) and nuclear (POMC and CRYBA1) partial sequences from 164 individuals were amplified, representing all species. Maximum Likelihood (ML) and Bayesian approaches were used to reconstruct phylogenetic relationships. Species tree was estimated using BEAST and singular value decomposition scores for species quartets (SVDquartets). Species limits were evaluated with six coalescent approaches. Diversification times were estimated using mitochondrial and nuclear rates with LogNormal relaxed clock in BEAST. Nine well-supported monophyletic lineages were recovered in Bayesian, ML, and SVDquartets, including eight named species and a lineage composed by specimens from the Villarrica population (Bootstrap:>70, PP:> 0.99). Single-locus species delimitation analyses overestimated the species number in E. migueli, E. calcaratus, and E. roseus lineages, while multi-locus analyses recovered as species the nine lineages observed in phylogenetic analyses (Ctax = 0.69). It is hypothesized that Eupsophus diversification occurred during Mid-Pleistocene (0.42-0.14 Mya), with most species having originated after the Last Southern Patagonian Glaciation (0.18 Mya). Our results revitalize the hypothesis that the E. roseus group is composed of eight species and support the Villarrica lineage as a new putative species.

Phylogeographic analysis reveals high genetic structure with uniform phenotypes in the paper wasp Protonectarina sylveirae (Hymenoptera: Vespidae)

Swarm-founding wasps are endemic and common representatives of neotropical fauna and compose an interesting social tribe of vespids, presenting both complex social characteristics and uncommon traits for a eusocial group, such as the absence of castes with distinct morphology. The paper wasp Protonectarina sylveirae (Saussure) presents a broad distribution from Brazil, Argentina and Paraguay, occurring widespread in the Atlantic rainforest and arboreal Caatinga, being absent in the Amazon region. Given the peculiar distribution among swarm-founding wasps, an integrative approach to reconstruct the evolutionary history of P. sylveirae in a spatial-temporal framework was performed to investigate: the presence of genetic structure and its relationship with the geography, the evolution of distinct morphologic lineages and the possible historical event(s) in Neotropical region, which could explain the observed phylogeographic pattern. Individuals of P. sylveirae were obtained from populations of 16 areas throughout its distribution for DNA extraction and amplification of mitochondrial genes 12S, 16S and COI. Analysis of genetic diversity, construction of haplotype net, analysis of population structure and dating analysis of divergence time were performed. A morphometric analysis was also performed using 8 measures of the body of the adult (workers) to test if there are morphological distinction among populations. Thirty-five haplotypes were identified, most of them exclusively of a group and a high population structure was found. The possibility of genetic divergence because of isolation by distance was rejected. Morphological analysis pointed to a great uniformity in phenotypes, with only a small degree of differentiation between populations of south and the remaining. Divergence time analysis showed a Middle/Late Miocene origin, a period where an extensive marine ingression occurred in South America. Divergence of haplogroups began from the Plio/Pleistocene boundary and the last glacial maximum most likely modeled the current distribution of species, even though it was not the cause of genetic breaks.

What happened in the South American Gran Chaco? Diversification of the endemic frog genus Lepidobatrachus Budgett, 1899 (Anura: Ceratophryidae)

The Chaco is one the most neglected and least studied regions of the world. This highly-seasonal semiarid biome is an extensive continuous plain without any geographic barrier, and in spite of its high species diversity, the events and processes responsible have never been assessed. Miocene marine introgressions and Pleistocene glaciations have been mentioned as putative drivers of diversification for some groups of vertebrates in adjacent biomes of southern South America. Here we used multilocus data (one mitochondrial and six nuclear loci) from the three species of the endemic frog genus Lepidobatrachus (Lepidobatrachus asper, Lepidobatrachus laevis, and Lepidobatrachus llanensis) to determine if any of the historical events suggested as drivers of vertebrate diversification in southern South America are related to the diversification of the genus and if the Chaco is indeed a biome without barriers. Using fossil calibration in a coalescent framework we estimated that the genus diversified in the second half of the Miocene, coinciding with marine introgressions. Genetic patterns and historical demography suggest an important role of old archs and cratons as refuges during floods. In one species of the genus, L. llanensis, genetic structure reveals some breaks along the landscape, the main one of which corresponds to an area of the central Chaco that may act as a climatic barrier. Additionally, we found differential effects of the main Chacoan rivers on species of Lepidobatrachus that could be related to the time of persistence of populations in the areas influenced by these rivers.

Species delimitation in frogs from South American temperate forests: The case of Eupsophus, a taxonomically complex genus with high phenotypic variation

One of the most characteristic and abundant amphibian taxa of South American temperate forests is Eupsophus. The ten currently recognized species of the genus have been divided in two species groups, roseus and vertebralis, but most of them, eight, belong to the roseus group. Recent phylogeographic and phylogenetic studies have suggested that species diversity of the roseus group could be underestimated. An examination of the literature shows that species of the roseus group exhibit high levels of variation in their external characteristics, particularly those used as diagnostic characters, which compromises their taxonomy and hinders their field recognition. High levels of variation were also observed in several new populations of the roseus group discovered in southern Chile (36°-40°S), which could not be identified to the species level by their external characteristics. On the other hand, the literature reveals a scarse karyotype differentiation and a high bioacoustic uniformity among the species of the roseus group. We performed a Bayesian phylogenetic analysis using mitochondrial and nuclear genes to reevaluate the species diversity of the roseus group, including all the nominal species of Eupsophus and new populations. This analysis was complemented with three species delimitation approaches, General Mixed Yule Coalescent, multi-rate Poisson Tree Process and Automatic Barcode Gap Discovery. We favored a conservative delimitation of only four species for the roseus group, a result more consistent with the distribution of pairwise genetic distances, and the available chromosome and bioacoustic evidence. The four recognized lineages, which have nearly completely allopatric distributions, are named after the earliest nominal species that they include, but because high levels of phenotypic variation, they are not diagnosable by consistent differences in external morphology. We discuss the implications of this new proposal for the taxonomy and conservation of the genus, and the possible causes of the difficulty to estimate its species diversity.

Tabula rasa in the Patagonian Channels? The phylogeography of Oreobolus obtusangulus (Cyperaceae)

The extent of the Pleistocene glaciations in the Patagonian Channel region (southwesternmost South America) and their impact on the vegetation there are largely unknown. Whether the regional flora was wiped out completely (tabula rasa) or survived in ice-free pockets (in situ survival) is still an open question. The molecular imprint of either scenario should still be visible in extant populations. Therefore, DNA sequence data of Oreobolus obtusangulus Gaudich. (Cyperaceae) were analysed. This species is an abundant constituent of Patagonian cushion peat bogs, one of the Patagonian Channel region's major vegetation types. Three hundred and eighty-four individuals from 48 populations were sequenced for two chloroplast (ycf3-psaA and trnQ^UUG -psbK intergenic spacers) and 14 nuclear loci containing simple sequence repeats (SSRs; microsatellites). Phylogenetic reconstructions and the geographic distribution of genetic diversity revealed that the species was split into three main lineages whose general distributions comprise three separate major regions, that is, south-central Chile, Fuego-Patagonia and the East Patagonian Andes, which probably constitute glacial refugia. Postglacial migration fronts formed a suture zone with high levels of genetic diversity in the Northwest Patagonian Andes, where remnants of a supposedly ancestral lineage were also found to be locally restricted to a single population (Huinay). The heavily glaciated Patagonian Channels were likely recolonized from the northwest, and partly from the south. Although the westernmost Patagonian Channel population (Estero Bachem) harboured private SSR alleles (singletons) and showed slightly elevated genetic diversity, it remained unclear whether this population actually survived in situ. This study helps fill a major gap in reconstructing the Pleistocene vegetation history of West and Andean Patagonia.

Genetics, Gene Flow, and Glaciation: The Case of the South American Limpet Nacella mytilina

Glacial episodes of the Quaternary, and particularly the Last Glacial Maximum (LGM) drastically altered the distribution of the Southern-Hemisphere biota, principally at higher latitudes. The irregular coastline of Patagonia expanding for more than 84.000 km constitutes a remarkable area to evaluate the effect of Quaternary landscape and seascape shifts over the demography of near-shore marine benthic organisms. Few studies describing the biogeographic responses of marine species to the LGM have been conducted in Patagonia, but existing data from coastal marine species have demonstrated marked genetic signatures of post-LGM recolonization and expansion. The kelp-dweller limpet Nacella mytilina is broadly distributed along the southern tip of South America and at the Falkland/Malvinas Islands. Considering its distribution, abundance, and narrow bathymetry, N. mytilina represents an appropriate model to infer how historical and contemporary processes affected the distribution of intraspecific genetic diversity and structure along the southern tip of South America. At the same time, it will be possible to determine how life history traits and the ecology of the species are responsible for the current pattern of gene flow and connectivity across the study area. We conducted phylogeographic and demographic inference analyses in N. mytilina from 12 localities along Pacific Patagonia (PP) and one population from the Falkland/Malvinas Islands (FI). Analyses of the mitochondrial gene COI in 300 individuals of N. mytilina revealed low levels of genetic polymorphism and the absence of genetic differentiation along PP. In contrast, FI showed a strong and significant differentiation from Pacific Patagonian populations. Higher levels of genetic diversity were also recorded in the FI population, together with a more expanded genealogy supporting the hypothesis of glacial persistence of the species in these islands. Haplotype genealogy, and mismatch analyses in the FI population recognized an older and more complex demographic history than in PP. Demographic reconstructions along PP suggest a post-LGM expansion process (7.5 ka), also supported by neutrality tests, mismatch distribution and maximum parsimony haplotype genealogies. Migration rate estimations showed evidence of asymmetrical gene flow from PP to FI. The absence of genetic differentiation, the presence of a single dominant haplotype, high estimated migration rates, and marked signal of recent demographic growth, support the hypothesis of rapid post-glacial expansion in N. mytilina along PP. This expansion could have been sustained by larval and rafting-mediated dispersal of adults from northernmost populations following the Cape Horn Current System. Marked genetic differentiation between PP and FI could be explained through differences in their respective glacial histories. During the LGM, Pacific Patagonia (PP) was almost fully covered by the Patagonian Ice Sheet, while sheet coverage in the FI ice was restricted to small cirques and valleys. As previously recorded in the sister-species N. magellanica, the FI rather than represent a classical glacial refugium for N. mytilina, seems to represent a sink area and/or a secondary contact zone. Accordingly, historical and contemporary processes, contrasting glacial histories between the analyzed sectors, as well as life history traits constitute the main factors explaining the current biogeographical patterns of most shallow Patagonian marine benthic organisms.

Isolation Driven Divergence in Osmoregulation in Galaxias maculatus (Jenyns, 1848) (Actinopterygii: Osmeriformes)

Background

Marine species have colonized extreme environments during evolution such as freshwater habitats. The amphidromous teleost fish, Galaxias maculatus is found mainly migrating between estuaries and rivers, but some landlocked populations have been described in lakes formed during the last deglaciation process in the Andes. In the present study we use mtDNA sequences to reconstruct the historical scenario of colonization of such a lake and evaluated the osmoregulatory shift associated to changes in habitat and life cycle between amphidromous and landlocked populations.

Results

Standard diversity indices including the average number of nucleotide differences (Π) and the haplotype diversity index (H) indicated that both populations were, as expected, genetically distinctive, being the landlocked population less diverse than the diadromous one. Similarly, pairwise G_ST and N_ST comparison detected statistically significant differences between both populations, while genealogy of haplotypes evidenced a recent founder effect from the diadromous stock, followed by an expansion process in the lake. To test for physiological differences, individuals of both populations were challenged with a range of salinities from 0 to 30 ppt for 8 days following a period of progressive acclimation. The results showed that the landlocked population had a surprisingly wider tolerance to salinity, as landlocked fish survival was 100% from 0 to 20 ppt, whereas diadromous fish survival was 100% only from 10 to 15 ppt. The activity of ATPase enzymes, including Na⁺/K⁺-ATPase (NKA), and H⁺-ATPase (HA) was measured in gills and intestine. Activity differences were detected between the populations at the lowest salinities, including differences in ATPases other than NKA and HA. Population differences in mortality are not reflected in enzyme activity differences, suggesting divergence in other processes.

Conclusions

These results clearly demonstrate the striking adaptive changes of G. maculatus osmoregulatory system, especially at hyposmotic environments, associated to a drastic shift in habitat and life cycle at a scale of a few thousand years.

Late Glacial Demographic Expansion Motivates a Clock Overhaul for Population Genetics

The molecular clock hypothesis is fundamental in evolutionary biology as by assuming constancy of the molecular rate it provides a timeframe for evolution. However, increasing evidence shows time dependence of inferred molecular rates with inflated values obtained using recent calibrations. As recent demographic calibrations are virtually non-existent in most species, older phylogenetic calibration points (>1 Ma) are commonly used, which overestimate demographic parameters. To obtain more reliable rates of molecular evolution for population studies, I propose the calibration of demographic transition (CDT) method, which uses the timing of climatic changes over the late glacial warming period to calibrate expansions in various species. Simulation approaches and empirical data sets from a diversity of species (from mollusk to humans) confirm that, when compared with other genealogy-based calibration methods, the CDT provides a robust and broadly applicable clock for population genetics. The resulting CDT rates of molecular evolution also confirm rate heterogeneity over time and among taxa. Comparisons of expansion dates with ecological evidence confirm the inaccuracy of phylogenetically derived divergence rates when dating population-level events. The CDT method opens opportunities for addressing issues such as demographic responses to past climate change and the origin of rate heterogeneity related to taxa, genes, time, and genetic information content.

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.

Amphibian molecular ecology and how it has informed conservation

Molecular ecology has become one of the key tools in the modern conservationist's kit. Here we review three areas where molecular ecology has been applied to amphibian conservation: genes on landscapes, within-population processes, and genes that matter. We summarize relevant analytical methods, recent important studies from the amphibian literature, and conservation implications for each section. Finally, we include five in-depth examples of how molecular ecology has been successfully applied to specific amphibian systems.

Microevolution due to pollution in amphibians: A review on the genetic erosion hypothesis

The loss of genetic diversity, due to exposure to chemical contamination (genetic erosion), is a major threat to population viability. Genetic erosion is the loss of genetic variation: the loss of alleles determining the value of a specific trait or set of traits. Almost a third of the known amphibian species is considered to be endangered and a decrease of genetic variability can push them to the verge of extinction. This review indicates that loss of genetic variation due to chemical contamination has effects on: 1) fitness, 2) environmental plasticity, 3) co-tolerance mechanisms, 4) trade-off mechanisms, and 5) tolerance to pathogens in amphibian populations.

Phylogeography in Galaxias maculatus (Jenyns, 1848) along Two Biogeographical Provinces in the Chilean Coast

Major geologic and climatic changes during the Quaternary exerted a major role in shaping past and contemporary distribution of genetic diversity and structure of aquatic organisms in southern South America. In fact, the northern glacial limit along the Pacific coast, an area of major environmental changes in terms of topography, currents, and water salinity, represents a major biogeographic transition for marine and freshwater species. We used mitochondrial DNA sequences (D-loop) to investigate the consequences of Quaternary glacial cycles over the pattern of genetic diversity and structure of G. maculatus (Pisces: Galaxiidae) along two biogeographical provinces in the Chilean coast. Extreme levels of genetic diversity and strong phylogeographic structure characterize the species suggesting a low amount of influence of the last glacial cycle over its demography. However, we recognized contrasting patterns of genetic diversity and structure between main biogeographical areas here analyzed. Along the Intermediate Area (38°-41° S) each estuarine population constitutes a different unit. In contrast, Magellanic populations (43°-53° S) exhibited low levels of genetic differentiation. Contrasting patterns of genetic diversity and structure recorded in the species between the analyzed biogeographic areas are consistent with the marked differences in abiotic factors (i.e., different coastal configurations, Quaternary glacial histories, and oceanographic regimes) and to inherent characteristics of the species (i.e., salt-tolerance, physiology, and reproductive behavior).

Phylogeography of the ancient catfish family Diplomystidae: biogeographic, systematic, and conservation implications

The catfish family Diplomystidae is one of the earliest branching lineages within the diverse order Siluriformes and shows a deep phylogenetic split from all other extant and extinct major catfish groups. Despite its relevance in the evolution of siluriforms, phylogenetic relationships within the Diplomystidae are poorly understood, and prior to this study, no phylogenetic hypotheses using molecular data had been published. By conducting a phylogeographic study across the entire distribution of the family, that encompasses river systems from Central-South Chile and Argentina, we provide the first molecular phylogenetic hypothesis among all known species of Diplomystidae, and in addition, investigate how their evolutionary history relates to major historical events that took place in southern South America. Our phylogenetic analyses show four main lineages and nine sub-lineages strongly structured geographically. All Pacific basin populations, with one exception (those found in the Baker basin) clustered within three of the four main lineages (clades I-III), while all populations from Atlantic basins and those from the Baker basin clustered in a single main clade (clade IV). There was a tendency for genetic diversity to decrease from north to south for Pacific basins consistent with an increasing north-south ice coverage during the last glacial maximum. However, we did not find a statistically significant correlation between genetic diversity and latitude. Analysis of molecular variance (AMOVA) showed that river basins and the barrier created by the Andes Mountains explained a high percentage of the genetic variation. Interestingly, most of the genetic variation among drainages was explained among Pacific basins. Molecular phylogenetic analyses agree only partially with current systematics. The geographical distribution of main lineages did not match species distribution and suggests a new taxonomic hypothesis with support for four species of Diplomystes, three species distributed allopatrically from the Rapel to the Valdivia basin, and only one species distributed in Baker and Atlantic basins. High genetic differentiation among river basins suggests that conservation efforts should focus on protecting populations in each basin in order to preserve the genetic diversity of one of the oldest groups of catfishes on the earth today.

Phylogeographical patterns shed light on evolutionary process in South America

The South American continent is composed of several biogeographical regions harbouring the highest biodiversity on the globe, encompassing five of the world's biodiversity 'hot spots'. Nonetheless, the patterns and processes responsible for shaping its astonishing species diversity are largely unknown. Here, we present a review of current South American phylogeographical knowledge based on published articles on this topic. An appraisal of the literature reveals emerging phylogeographical patterns in the biota of South America. The striking phylogeographical divergence observed among organism lineages in South American studies is suggestive of high levels of undocumented species diversity. The interplay between Pleistocene climatic oscillations and Pliocene/Miocene orogenic events has contributed to shaping the current diversity and distribution of modern lineages in both the tropical and temperate regions of South America. Although older divergence times were observed for a range of species, most herpetofauna underwent an intraspecific lineage split much earlier than other organisms. The geographical ranges of species associated with forest habitats were reduced mainly during glacial cycles, whereas species associated with open vegetation domains have shown variable responses to climatic oscillations. The results suggest a highly complex mosaic of phylogeographical patterns in South America. We suggest future research directions to promote a better understanding of the origin and maintenance of the South American biota.

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.

Early stages of divergence: phylogeography, climate modeling, and morphological differentiation in the South American lizard Liolaemus petrophilus (Squamata: Liolaemidae)

This study examines the phylogeographic structure within the Patagonian lizard Liolaemus petrophilus and tests for patterns of between-clade morphological divergence and sexual dimorphism, as well as demographic and niche changes associated with Pleistocene climate changes. We inferred intraspecific relationships, tested hypotheses for historical patterns of population expansion, and incorporated ecological niche modeling (ENM) with standard morphological and geometric morphometric analyses to examine between-clade divergence as indirect evidence for adaptation to different niches. The two inferred haploclades diverged during the early Pleistocene with the Southern clade depicting the genetic signature of a recent population increase associated with expanding niche envelope, whereas the Northern clade shows stable populations in a shrinking niche envelope. The combination of molecular evidence for postisolation demographic change and ENM, suggest that the two haploclades have responded differently to Pleistocene climatic events.

Towards a model of postglacial biogeography in shallow marine species along the Patagonian Province: lessons from the limpet Nacella magellanica (Gmelin, 1791)

BACKGROUND

Patagonia extends for more than 84,000 km of irregular coasts is an area especially apt to evaluate how historic and contemporary processes influence the distribution and connectivity of shallow marine benthic organisms. The true limpet Nacella magellanica has a wide distribution in this province and represents a suitable model to infer the Quaternary glacial legacy on marine benthic organisms. This species inhabits ice-free rocky ecosystems, has a narrow bathymetric range and consequently should have been severely affected by recurrent glacial cycles during the Quaternary. We performed phylogeographic and demographic analyses of N. magellanica from 14 localities along its distribution in Pacific Patagonia, Atlantic Patagonia, and the Falkland/Malvinas Islands.

RESULTS

Mitochondrial (COI) DNA analyses of 357 individuals of N. magellanica revealed an absence of genetic differentiation in the species with a single genetic unit along Pacific Patagonia. However, we detected significant genetic differences among three main groups named Pacific Patagonia, Atlantic Patagonia and Falkland/Malvinas Islands. Migration rate estimations indicated asymmetrical gene flow, primarily from Pacific Patagonia to Atlantic Patagonia (Nem=2.21) and the Falkland/Malvinas Islands (Nem=16.6). Demographic reconstruction in Pacific Patagonia suggests a recent recolonization process (< 10 ka) supported by neutrality tests, mismatch distribution and the median-joining haplotype genealogy.

CONCLUSIONS

Absence of genetic structure, a single dominant haplotype, lack of correlation between geographic and genetic distance, high estimated migration rates and the signal of recent demographic growth represent a large body of evidence supporting the hypothesis of rapid postglacial expansion in this species in Pacific Patagonia. This expansion could have been sustained by larval dispersal following the main current system in this area. Lower levels of genetic diversity in inland sea areas suggest that fjords and channels represent the areas most recently colonized by the species. Hence recolonization seems to follow a west to east direction to areas that were progressively deglaciated. Significant genetic differences among Pacific, Atlantic and Falkland/Malvinas Islands populations may be also explained through disparities in their respective glaciological and geological histories. The Falkland/Malvinas Islands, more than representing a glacial refugium for the species, seems to constitute a sink area considering the strong asymmetric gene flow detected from Pacific to Atlantic sectors. These results suggest that historical and contemporary processes represent the main factors shaping the modern biogeography of most shallow marine benthic invertebrates inhabiting the Patagonian Province.

Species replacement along a linear coastal habitat: phylogeography and speciation in the red alga Mazzaella laminarioides along the south east Pacific

Background

The Chilean shoreline, a nearly strait line of coast expanding across 35 latitudinal degrees, represents an interesting region to assess historical processes using phylogeographic analyses. Stretching along the temperate section of the East Pacific margin, the region is characterized by intense geologic activity and has experienced drastic geomorphological transformations linked to eustatic and isostatic changes during the Quaternary. In this study, we used two molecular markers to evaluate the existence of phylogeographic discontinuities and detect the genetic footprints of Pleistocene glaciations among Patagonian populations of Mazzaella laminarioides, a low-dispersal benthic intertidal red seaweed that inhabits along ~3,700 km of the Chilean coastal rocky shore.

Results

Three main genetic lineages were found within M. laminarioides. They are distributed along the Chilean coast in strict parapatry. The deep divergence among lineages suggests that they could be considered putative genetic sibling species. Unexpectedly, genetic breaks were not strictly concordant with the biogeographic breaks described in the region. A Northern lineage was restricted to a broad transition zone located between 30°S and 33°S and showed signals of a recent bottleneck. The reduction of population size could be related to warm events linked to El Niño Southern Oscillation, which is known to cause massive seaweed mortality in this region. To the south, we propose that transient habitat discontinuities driven by episodic tectonic uplifting of the shoreline around the Arauco region (37°S-38°S); one of the most active forearc-basins in the South East Pacific; could be at the origin of the Central/South genetic break. The large beaches, located around 38°S, are likely to contribute to the lineages’ integrity by limiting present gene flow. Finally, the Southern lineage, occupies an area affected by ice-cover during the last glaciations. Phylogeny suggested it is a derived clade and demographic analyses showed the lineage has a typical signature of postglacial recolonization from a northern glacial refugium area.

Conclusions

Even if environmental adaptation could have strengthened divergence among lineages in M. laminarioides, low dispersal capacity and small population size are sufficient to generate phylogeographic discontinuities determined by genetic drift alone. Interestingly, our results confirm that seaweed population connectivity over large geographic scales does not rely only on dispersal capacity but also seem to depend highly on substratum availability and population density of the receiving locality.