Skip to main content

Population structure and microscale morphological differentiation in a freshwater snail from the Chilean Altiplano

Abstract

Background

The diversity and population genetic structure of many species have been shaped by historical and contemporary climatic changes. For the species of the South American Altiplano, the historical climatic changes are mainly related to the wet events of great magnitude and regional influence that occurred during the Pleistocene climatic oscillations (PCOs). In contrast, contemporary climate changes are associated with events of lesser magnitude and local influence related to intensifications of the South American Summer Monsoon (SASM). Although multiple studies have analyzed the effect of PCOs on the genetic patterns of highland aquatic species, little is known about the impact of contemporary climate changes in recent evolutionary history. Therefore, in this study, we investigated the change in population structure and connectivity using nuclear and mitochondrial markers throughout the distribution range of Heleobia ascotanensis, a freshwater Cochliopidae endemic to the Ascotán Saltpan. In addition, using geometric morphometric analyses, we evaluated the concomitance of genetic divergence and morphological differentiation.

Results

The mitochondrial sequence analysis results revealed the presence of highly divergent co-distributed and geographically nested haplotypes. This pattern reflects an extension in the distribution of groups that previously would have differentiated allopatrically. These changes in distribution would have covered the entire saltpan and would be associated with the large-scale wet events of the PCOs. On the other hand, the microsatellite results defined five spatially isolated populations, separated primarily by geographic barriers. Contemporary gene flow analyses suggest that post-PCO, climatic events that would have connected all populations did not occur. The morphometric analyses results indicate that there is significant morphological differentiation in the populations that are more isolated and that present the greatest genetic divergence.

Conclusions

The contemporary population structure and morphological variation of H. ascotanensis mainly reflect the post-PCO climatic influence. Although both markers exhibit high genetic structuring, the microsatellite and morphology results show the preponderant influence of fragmentation in recent evolutionary history. The contemporary genetic pattern shows that in species that have limited dispersal capabilities, genetic discontinuities can appear rapidly, erasing signs of historical connectivity.

Peer Review reports

Background

The Altiplano landscape has experienced a complex environmental history influenced mainly by Andean orogenesis, glaciations and intense volcanic activity. During the Pleistocene climatic oscillations (PCOs), numerous drastic changes occurred in the temperature and rainfall of the region [1,2,3]. As a consequence of these changes, a series of humid and arid periods occurred in the Altiplano. The humid periods led to the formation of extensive paleolakes that connected currently isolated hydrological systems [4,5,6]. During these lake phases, multiple basins coalesced into paleolakes that would have connected, in extreme events, practically the entire Altiplano basin (e.g., the Titicaca-Desaguadero-Tauca system) during the Last Glacial Maximum [7, 8]. Later, during periods of aridity, fragmentation of the paleolakes would have occurred, giving rise to isolated hydrological systems we see today.

Various studies suggest that there is an association between the climatic changes that occurred during the Quaternary Period and the differentiation of the Altiplano taxa. Among the strictly aquatic animal models that support this hypothesis, research carried out on gastropods of the genus Heleobia Stimpson, 1865, and fish of the genus Orestias Valenciennes, 1839, stand out. In both cases, the results indicate that the species of these taxa would have diverged allopatrically during the transitions of wet and arid periods that occurred in the region [9,10,11,12]. On the other hand, the distribution patterns of co-distributed taxa, including gastropods of the genus Biomphalaria Preston, 1910, and amphipods of the genus Hyalella Smith, 1874, suggest that the evolutionary scenario of divergence of the highland biota could have involved multiple speciation processes (e.g., allopatric, secondary contact and dispersion) [13,14,15,16,17]. It is important to mention that while all the co-distributed taxa with Heleobia showed an influence of climatic changes in their evolutionary history, the differences in the patterns and processes are likely related to their life history traits. For example, Biomphalaria is a pulmonated gastropod, phylogenetically distant from Heleobia, and it exhibits physiological and phenotypic differences that are probably linked to the differences in the phylogeographic pattern observed between these species.

Although multiple studies have analyzed the effect of the PCOs on the diversification of the aquatic fauna of the Altiplano, few studies have focused on evaluating the effect of recent fragmentation (i.e., contemporary or post-PCO) on genetic divergence. This line of inquiry is particularly rare in macroinvertebrate groups that have historically received little attention and that in most cases lack taxonomic descriptions or systematic studies. On a recent timescale, the main environmental changes that have modified the landscape in the Altiplano are due to intensification events of the South American Summer Monsoon (SASM) [18,19,20,21,22,23,24,25]. These events have local effects and correspond to periods of increased in rainfall [26, 27] and, in extreme cases, produce catastrophic floods, such as those reported in 2015 in the Atacama Desert in Chile [28]. In this sense, although these events are not of the magnitude of the PCOs, it is likely that they have generated connections between nearby hydrological systems.

The Ascotán Saltpan (Antofagasta Region, Chile) is an endorheic basin of the southern Altiplano. The saltpan has a surface covered mainly by evaporite deposits, and the surface water is provided by a series of springs with different levels of connectivity (Fig. 1). For the strictly aquatic organisms of the saltpan, fragmentation of the habitat would have favored the differentiation of the populations, as terrestrial barriers prevented their dispersal and decreased the genetic flow between populations [14, 29,30,31]. Furthermore, in desert hydrological systems similar to the Ascotán Saltpan, it has been suggested that for freshwater gastropods, the isolation of the resident basins would induce considerable genetic differentiation [32]. Among the species that are distributed in the Ascotán Saltpan is Heleobia ascotanensis (Courty, 1907), and five subspecies that have been described based on conchological characters [33, 34]. Heleobia ascotanensis is a strictly aquatic snail; therefore, its distribution within the saltpan is restricted to the springs that make up the only bodies of freshwater within the system [30, 35]. It is important to mention that this restricted distribution to the springs of the saltpan could be modified, as passive dispersal mediated by birds may occur, as reported in similar freshwater snails from desert hydrological systems in Australia [36,37,38] or even over greater distances in the case of migratory birds [39]. However, there is no evidence that this bird-mediated dispersal occurs in Heleobia, and it has not been reported for other species from the Ascotán Saltpan or the Chilean Altiplano.

Fig. 1
figure 1

Sampling localities of Heleobia ascotanensis used in the present study. A Map showing the Ascotán Saltpan and the 13 springs analyzed. B Detail of springs 2 to 6. C and D correspond to photographs of springs 2 and 9, respectively. Each photograph shows representative individuals of each population (the black circle represents 5 mm). All photographs by Moisés A. Valladares. The map was made using QGIS Geographic Information System v3.4.9 with ESRI world imagery (ESRI, DigitalGlobe, GeoEye, Earthstar Geographics, CNES/Airbus DS, USDA, USGS, AeroGRID, IGN, and the GIS User Community) (http://www.qgis.org, accessed on April 10, 2023)

According to the environmental regulations in Chile, H. ascotanensis has been classified by the Ministry of the Environment as an endangered species (RCE DS38/2015), a classification based mainly on the decrease in the species’ habitat quality and since the species is estimated to occupy a range of less than 18 km2 [40]. According to the phylogenetic analysis carried out by Collado et al. [12] the populations of H. ascotanensis are genetically related to the highland populations of El Tatio (northern Chile) and Laguna Blanca (southern Bolivia); however, these relationships have not been evaluated in detail. Recently, a phylogeographic study of H. ascotanensis using the COI mitochondrial genetic marker suggested that three populations coexist in the system, which present low and moderate degrees of historical connectivity, in addition to exhibiting incipient morphological differentiation as a consequence of genetic and geographic isolation [30]. However, to date, no study has evaluated the structure and contemporary genetic flow of H. ascotanensis populations, for which more variable genetic markers are needed to describe the population structure and identify recent demographic processes (i.e., on a generational scale).

Based on the aforementioned background, the present study examined the genetic diversity, population structure and connectivity of H. ascotanensis populations using mitochondrial and nuclear DNA sequences (microsatellites). Considering the properties of each marker, this approach allowed us to contrast the historical population history obtained from the mitochondrial information with the contemporary history inferred from the microsatellite information. Additionally, the congruence between the genetic divergence and the degree of morphological differentiation among populations was evaluated using geometric morphometric analyses. By integrating multiple lines of evidence, we generated a detailed description of the evolutionary and demographic history of H. ascotanensis, illustrating how historical and contemporary habitat fragmentation processes have modulated population genetic diversity and promoted phenotypic differentiation in this species.

Results

Mitochondrial genetic structure and diversity

A total of 310 individuals from the 13 springs of the Ascotán Saltpan were sequenced for the COI gene. The alignment had a total length of 658 bp and presented 50 variable sites (24 were parsimony informative) that defined 56 haplotypes. Table 1 shows the values of the genetic diversity indices calculated from the COI sequences. The lowest genetic variability was found in Spring 9, where a single polymorphic site and two haplotypes were detected. In contrast, the greatest variability was found in Spring 13 (western margin of the saltpan), which presented 14 polymorphic sites and 16 haplotypes. On the other hand, Springs 2 to 6 presented the greatest intra-spring differentiation since, for all these systems, the average number of differences per pair of sequences was higher than 4.7, with a maximum of 5.437 in Spring 3.

Table 1 Genetic diversity indices obtained from COI gene sequences for each spring of the Ascotán Saltpan. The geographic coordinates for each spring, sample size (N), number of polymorphic sites (S), number of haplotypes (K), haplotype diversity (H) and average of differences between pairs of sequences (Π) are indicated

In the haplotype network (Fig. 2), it was observed that within the saltpan, there were no clearly differentiated haplogroups, and only Spring 11 presented solely exclusive haplotypes. On the other hand, a divergent haplotype was detected in Spring 1 (separated by five mutational steps from the rest of the network, A in the figure), and a high-frequency haplotype was present only between Springs 2 to 6 (separated by six mutational steps from the rest of the network, B in the figure). In general, the network showed that many of the haplotypes are shared between springs, with the exception of those that occur in Spring 11. The grouping of the haplotypes of Springs 11 and 13 shows a dispersed and extended network. On the other hand, for the rest of the springs, the grouping pattern exhibits star-shaped branching. The network shows a central haplotype (C in the figure) that encompasses all springs, except for Spring 11. Spring 13 included multiple exclusive haplotypes (D and E in the figure) and representatives of the main haplotype (C in the figure).

Fig. 2
figure 2

Haplotype network of Heleobia ascotanensis from the Ascotán Saltpan generated using COI gene sequences. The segments on the lines represent the number of mutational steps. The size of the circumferences represents the relative frequency of each haplotype. A, B, D and E indicate divergent haplotypes, and C indicates the higher-frequency haplotype (details in main text)

The genetic differences between pairs of springs of the Ascotán Saltpan obtained from the COI gene are shown in Table 2 (below the diagonal). The results indicate that Spring 11 is the most differentiated system since it presents the highest FST values, and all comparisons were significant. Spring 13 (western margin of the saltpan) also shows high differentiation values. In general, the results of the FST test suggest that genetic divergence is associated with the spatial positions of the springs, since more distant springs presented higher differentiation values.

Table 2 FST values between pairs of springs of the Ascotán Saltpan. The distances calculated using the mitochondrial genetic sequences of the COI marker are shown below the diagonal, and the microsatellite results are shown above the diagonal. For the COI marker, the significant results (p < 0.05) are highlighted in bold. For microsatellites, the analysis was performed using the ENA correction for null alleles and 10,000 bootstrap replicates

The Mantel test was used to analyze the significance of the correlation between the matrix of geographic distances and genetic distances (FST) obtained using COI sequences; a Pearson’s correlation coefficient r = 0.33 (p < 0.05) was obtained. This result suggests that there is a significant relationship between geographic distance and genetic differentiation. Therefore, the pattern of genetic divergence of the Heleobia populations of the Ascotán Saltpan is congruent with a scenario of isolation by distance.

Nuclear microsatellite markers

Genotypes were obtained from 390 individuals from the 13 springs of the Ascotán Saltpan (10 microsatellite loci). This set of data contained all the individuals of the mitochondrial analyses; in addition, the sample size was increased in all springs. All the loci (except the locus Hel45) showed deviations from Hardy-Weinberg equilibrium in at least one of the springs (Table S1). The diversity indices showed that eight of the 13 springs presented private or exclusive alleles (Table 3); in this sense, Springs 8 and 13 showed the highest number of private alleles (AP = 5) in the saltpan. In 10 of the 13 springs, a deficit of heterozygotes (HO < HE) was detected, which was confirmed by high values of the fixation index (F). Similarly, in 11 springs, positive values of internal relatedness (i.e., inbreeding) were detected; in particular, the values were higher for individuals of Springs 4, 6 and 8. Only individuals of Springs 10 and 12 presented IR values close to zero (i.e., unrelated parents).

Table 3 Genetic diversity indices of Heleobia ascotanensis populations determined using 10 microsatellite loci. The number of samples (N), number of alleles (A), number of deprived alleles (AP), average allelic richness (AR), percentage of polymorphic loci (% PL), observed heterozygosity (HO), expected heterozygosity (HE), fixation index (F) and internal relatedness (IR) are shown. For % PL, HO, HE, F and IR, the population mean of the parameter and its standard error are shown

The paired FST(SSR) test results (Table 2, above the diagonal) show that in general, the geographically close springs presented low differentiation values. This phenomenon was detected, for example, among Springs 1 to 7 (except for some comparisons of Spring 2) and between Springs 8 and 12. Spring 13 presented the highest differentiation values with respect to the rest of the systems, except for the comparison with Spring 10. The Mantel test between the matrix of geographic distances and microsatellite genetic distances (FST(SSR)) yielded a Pearson’s correlation coefficient r = 0.54 (p < 0.05). The correlation coefficient was higher than that obtained using mitochondrial data, confirming the importance of isolation and geographic distance in the genetic structure of H. ascotanensis.

Population genetic structure

The most likely number of populations in the initial Structure round, determined using Evanno’s method, was KSSR = 2 (Fig. 3). This partitioning of the dataset clearly segregated the northern and southern springs of the Ascotán Saltpan. A second round of Structure analysis was performed for each of the groups separately. In the second round of the northern group, no obvious structuring was detected; therefore, Springs 1 to 6 constituted a population (hereinafter Population 1). On the other hand, the second round of the southern group defined three populations (KSSR = 3), and this new partitioning of the data also reflected the spatial arrangement of the springs in the saltpan. In the third round of the southern group, the population formed by Springs 7 and 8 (Population 2) and the population composed of Springs 9, 10 and 12 (Population 3) were defined. Finally, in the fourth round of Structure, the most isolated springs of the saltpan were analyzed, and both were defined as independent populations (Spring 11 constituted Population 4, and Spring 13 Population 5).

Fig. 3
figure 3

Population structure of Heleobia ascotanensis from 10 microsatellite loci using a hierarchical Structure analysis

The migration rates estimated using BayesAss showed high self-recruitment values for Populations 1, 2 and 3 (98.7, 95.7 and 96.5%, respectively) (Fig. 4). In addition, none of these populations was a source of migration to another population with values greater than 1% (Table S2). With respect to the most isolated populations, Population 4 (Spring 11) presented the lowest self-recruitment value in the saltpan, at 67.6%, and contributed 24.4% to Population 5 (Spring 13) and 4.9% to Population 2 (Springs 7 and 8). Finally, Population 5 presented a moderate rate of self-recruitment (93.5%) and contributed 3.0% of migrants to Population 3 (Springs 9, 10 and 12).

Fig. 4
figure 4

Estimated migration rates for Heleobia ascotanensis populations from 10 microsatellite loci in the Ascotán Saltpan. The groupings correspond to the populations defined by Structure. Rates greater than 2% are indicated, and the self-recruitment rate is highlighted in bold

No significant signs of genetic bottlenecks were detected in any of the populations defined by Structure using both the two-phase model (TPM) and the stepwise mutational model (SMM) in Bottleneck (all p values > 0.1 in the Wilcoxon tests). Similarly, the five populations presented normal L-shaped distributions, indicating the absence of abrupt demographic declines in the populations.

Geometric morphometric analyses

Significant differences in the size distributions (response variable) considering the populations defined by Structure were evaluated by means of a Procrustes ANOVA. In all analyses, the logarithm of the centroid size [Log (Centroid Size)] was used as a proxy for the size of the specimens. The results showed significant differences in the size distribution (F = 13.81; Z = 5.60; p < 0.05) between populations. Post hoc comparisons indicated that Populations 3 and 4 (Springs 9, 10 and 12 and Spring 11, respectively) presented significant differences in the size of the centroid compared to the rest of the populations (Fig. 5).

Fig. 5
figure 5

Morphological differentiation associated with shell size among Heleobia ascotanensis populations. Differences of the logarithm of the size of the centroid classifying the individuals according to the populations defined by Structure. The population mean of the centroid size is shown, and significant differences are indicated (Bonferroni post hoc)

The principal component analysis of the individuals classified by population showed an overlap in the morphological space (Fig. 6), and only Population 4 (Spring 11) encompassed a bounded region in the center of the morphospace. The deformation grids for extreme values of Principal Component 1 (PC1) showed that this axis was mainly associated with the size and elongation of the shell. On the other hand, PC2 was associated with the position and shape of the peristome and showed the differentiation in the relative width of the shell of the populations of H. ascotanensis. The MANOVA results showed significant differences between populations (F = 2.08; Z = 12.56; p < 0.05). The penalized discriminant analysis (PDA) for the individuals classified by population had an overall precision (for the entire dataset) of 90% (p < 0.05). The percentages of correct classification by population were Population 1 = 95%; Population 2 = 87%; Population 3 = 74%; Population 4 = 89%; and Population 5 = 96%.

Fig. 6
figure 6

Morphological space of the populations of Heleobia ascotanensis. Principal component analysis using the GPA matrix of the populations defined by Structure. The deformation grids for the extreme shapes of each axis of variation in relation to the global consensus are shown in the corners of the graph

Discussion

In the present study, genetic diversity was estimated from mitochondrial and nuclear data, and the demographic history and population structure of the gastropod Heleobia ascotanensis from the Ascotán Saltpan were inferred. The results indicate that the analyzed samples are structured into five populations that present different degrees of diversity, differentiation and genetic flow. This pattern is expected considering that hydrological systems composed of springs with a certain degree of isolation are habitats where endemic species have evolved in a variety of groups, which apparently have species in relatively stable and isolated conditions [35, 41,42,43,44]. On the other hand, the morphometric analyses results indicate that there is incipient morphological differentiation in the populations that are more isolated. Furthermore, although there is considerable overlap in the multivariate space, the analyses revealed significant differences in all comparisons. The results obtained using mitochondrial sequences are consistent with those reported by Valladares et al. [30]. In both cases, high genetic structuring was detected among the Heleobia populations. However, differences were also detected in the composition of the populations, and this situation is probably related to the incorporation of a new locality (Spring 13).

The genetic pattern obtained from microsatellite loci was partially discordant with that obtained from mitochondrial data. In this sense, although both sources of evidence suggest a marked phylogeographic structure influenced by the geographic isolation between the springs, the microsatellite data reflect a more prevalent modulation of the landscape. A result that supports this proposition is that for both sets of data, a significant and positive relationship was detected between the geographic distance of the springs and the genetic distance of the populations. This correlation is consistent with signals of a pattern associated with isolation by distance. However, the Pearson correlation index of the Mantel tests was higher for microsatellites (rmicrosatellites = 0.54) than for mitochondrial sequences (rCOI = 0.33). It is likely that this discordance is due to the rapid rate of variation of microsatellite markers compared to that of mitochondrial sequences [45, 46]. Therefore, these markers are much more sensitive to the contemporary geographic and environmental influence imposed by the system [47].

The mitochondrial sequence analysis results showed great genetic variability and high intra-spring differentiation in Springs 2 to 6 (north of the saltpan). In the haplotype network, two haplotypes separated by multiple mutational steps (B and C in Fig. 2) were detected in these springs, a finding accompanied by high values of differences between pairs of sequences (Π). According to Valladares et al. [30], this pattern reflects an allopatric divergence between groups (i.e., they come from different more recent common ancestors, MRCA), followed by an extension in the distribution range of the C haplotype from the south of the saltpan to Springs 2–6. Considering the great geographic extension of haplotype C, it is probable that the extension of the range of this species occurred during the climatic events of great magnitude that characterized the PCO and that connected the more isolated springs (Springs 1 and 13). This hypothesis is congruent with our results and is also consistent with the results of phylogeographic studies carried out in the endemic fish Orestias ascotanensis Parenti, 1984, of the Ascotán Saltpan [29, 31]. In contrast, the microsatellite results did not show evident genetic segregation among individuals from the northern springs of the saltpan (Population 1, Springs 1–6). However, the inferred ancestry ratios in Structure suggest that the individuals came from two ancestral populations. On the one hand, this result supports the hypothesis that two genetically divergent groups co-occurred in the northern springs of the saltpan and later intermixed. On the other hand, unlike the mitochondrial pattern, the microsatellite results suggest that recently (i.e., post-PCO), climatic events of great magnitude that connected all the springs have not occurred. In contrast, habitat fragmentation would be the predominant condition, and population dynamics would mainly reflect the influence of geographic isolation.

Gene flow analysis using microsatellite markers suggests that migration between populations occurs asymmetrically and in a unidirectional direction from south to north. In this sense, Population 1 (north of the saltpan, Springs 1–6) would be the most isolated since it showed the highest rate of self-recruitment (98.7%) and would not receive migrants from other populations. Populations 2 and 3 (center of the saltpan) also showed high rates of self-recruitment and would also receive migrants from Populations 4 and 5, respectively. Although Population 5 (Spring 13) is located on the western margin of the saltpan and a great distance from the rest of the systems, the results show that it would receive a large proportion of migrants from Population 4 (Spring 11). Finally, Population 4 is the southernmost in the saltpan and presented the lowest self-recruitment rate (67.6%). This unidirectional flow from south to north can be associated with the altitudinal profile of the saltpan since the southern springs are higher than those located in the northern part of the system [29, 48]. Therefore, in events of flooding of the saltpan, a product of anomalies in the periods of precipitation [27], runoff would occur mainly from the south to the north. This hypothesis is compatible with the connectivity models that have been proposed for the Ascotán Saltpan and that have been evaluated in different aquatic organisms [14, 29,30,31].

The morphological analysis results showed that there are significant differences in the size and shape of the shell between the populations of H. ascotanensis. The geometric morphometrics approach was used to detect subtle differences in the shell and is a valuable tool that could be used in systematic studies. This approach is important considering that in the shell of Heleobia, there are no meristic features such as spines, ribs or other morphological characters that can be used as diagnostic features in the identification or delimitation of species. Although all populations showed significant differences in shell shape, Population 4 (Spring 11) also showed differences in size. This result is important because molecular analyses consistently show genetic differentiation between individuals in this population and the rest of the system. Considering both sources of evidence (morphology and genetics), the divergence process is more pronounced in Population 4 with respect to the rest of the populations of the species. According to Valladares et al. [30], the population of Spring 11 would have been isolated from the rest of the system during the last three lake cycles that occurred in the Altiplano (Sajsi [c. 21–25 kya], Tauca [c. 14–18 kya] and Coipasa [c. 11–13 kya]). Therefore, it is likely that the temporal and geographic isolation in which Population 4 has evolved has led to morphological differentiation and genetic divergence. A similar pattern was detected in Population 5 (Spring 13; western margin of the saltpan); however, the morphological differentiation was less than that detected in Population 4 and presented haplotypes shared with other springs.

The results of the present study show that it is necessary to resort to different sources of evidence to obtain a robust understanding of the pattern of population structuring. The latter is of particular importance when the species of interest is distributed in a complex hydrological system of springs with different degrees of connectivity. Our results and studies carried out in similar desert systems [29,30,31, 44, 47, 49,50,51,52] show that defining populations a priori based on geographical aspects is not appropriate. Recognizing this precaution is crucial in the case of threatened species such as H. ascotanensis since the correct description of the population structure and evolutionary history are necessary to undertake adequate conservation strategies [53,54,55]. In this sense, although genetic bottlenecks were not observed in any of the populations, the elevated levels of inbreeding in the majority of springs raise concerns. This pattern is probably modulated by geographic isolation, as shown by the high values of internal relatedness of the springs that make up the most isolated populations. However, this pattern is likely due to the intrinsic components related to the low vagility and reproductive strategy (direct development) reported for the species of the genus Heleobia [47, 56]. Our results contribute to the knowledge of the species and fauna of macroinvertebrates of the South American Altiplano. Future studies should integrate more variable markers (e.g., SNPs) to decipher the finer genetic pattern but should also undertake monitoring that assesses temporary changes in the population structure.

Conclusions

Multiple complementary methods were used to evaluate the genetic diversity, population structure, demographic history and morphological differentiation of Heleobia ascotanensis. The results of the analysis of the mitochondrial sequences suggest that the historical genetic structure would have been influenced mainly by the wet events associated with the PCO that would have connected currently isolated springs. In contrast, the contemporary population structure inferred using microsatellite markers reflects the fragmentation of the system due to the arid conditions that have prevailed in the post-PCO Altiplano. Five populations were defined that were consistent with the geography of the saltpan and the degree of isolation of the springs. The results of the gene flow analysis showed that migration occurs only from south to north, responding to the altitudinal profile of the saltpan. Finally, the morphometrics study showed that genetic differentiation was accompanied by significant differences in the shape and size of the shell.

Methods

Study and sampling area

Individuals of Heleobia ascotanensis were collected from 13 springs of the Ascotán Saltpan (Antofagasta Region, Chile) in 2019. The sampling was authorized by the Undersecretary of Fisheries and Aquaculture (SUBPESCA) in accordance with Exempt Resolution No. 2604/19. The sampling sites corresponded to 12 springs of the eastern margin of the saltpan (Springs 1 to 12, S1–S12) and a spring located in the western margin of the system (Spring 13, S13) (Fig. 1). The individuals used for the molecular analyses were previously photographed to be included in the subsequent morphological analyses. To avoid problems associated with allometric changes during ontogeny, the morphological study used adult individuals that were in the upper third of the size distribution of each sample collected.

Extraction and amplification of DNA sequences

The DNA samples were extracted and amplified by PCR following the protocol of Valladares et al. [30, 47]. For the COI marker, both strands of the amplification product were sequenced by Macrogen Inc. (Seoul, Korea). The nucleotide sequences were edited in the CodonCode Aligner v.8.0.1 program (CodonCode Corporation, Dedham, MA, USA) and were then aligned using the MAFFT v7.505 algorithm using the web server [57, 58].

For microsatellite markers, seven specific primers were developed for H. ascotanensis following the protocol of Fabres et al. [59], and three other primers previously reported for the congeneric species Heleobia atacamensis (Philippi, 1860) were used [59]. The amplification products were analyzed by the sequencing service of the Pontificia Universidad Católica de Chile (Santiago, Chile). Then, fragment analysis was performed using GeneMapper v4.0 software (Applied Biosystems, Waltham, MA). Subsequently, the excess of homozygotes or heterozygotes was checked with Micro-Checker v2.2.3 [60], and for each locus, the presence of stutters and null alleles were searched. The details of the mitochondrial and microsatellites are indicated in Table S3.

Mitochondrial DNA sequence analysis

The genetic diversity of each spring was described by the following indices calculated in DnaSP v6.0.1 [61]: number of polymorphic sites (S), number of haplotypes (K), haplotype diversity (H) and average number of differences between pairs of sequences (Π). The relationships between the haplotypes recovered for the COI gene were visualized by constructing a haplotype network using the median-joining algorithm [62] in PopART v1.7 [63]. To examine the differentiation between each pair of systems, a paired FST test was performed in the Arlequin v3.5 program [64]. The significance of each paired FST index was obtained based on 10,000 permutations. The existence of a pattern of genetic divergence associated with a phenomenon of isolation by distance (IBD) was tested by means of a Mantel test using matrices of genetic distances (dependent variable) and geographic distances (independent variable). The spatial distance matrix corresponded to the geographic distance in meters standardized between the systems studied. The genetic distance matrix was linearized according to the formula [FST/(1- FST)] proposed by Rousset [65]. The analysis was performed through 10,000 permutations using the vegan package v2.5–6 [66] in R software v3.6.1 [67].

Analysis of microsatellite sequences

Possible deviations from Hardy-Weinberg equilibrium were evaluated by loci and populations using an exact test in Genepop v4.7.5 [68]. Indices of genetic diversity considering the number of alleles (A), private number of alleles (AP), percentage of polymorphic loci (% PL), observed heterozygosity (HO) and expected heterozygosity (HE) and fixation index (F) were calculated in GenAlEx v6.5 [69]. Allelic richness (AR) was calculated in the PopGenReport v3.0.4 package [70] using rarefaction considering the population with the smallest sample size. To estimate the degree of inbreeding, the population mean internal relatedness (IR) was calculated in the Rhh v1.0.2 package [71]. The IR index assigns positive or negative values, where positive values suggest inbred individuals and negative values suggest exogamous individuals, while values close to zero suggest that the individuals come from unrelated parents [72]. Finally, a paired matrix of FST(SSR) (where SSR denotes microsatellite results) was calculated in FreeNA v1.0 [73] (using the null allele exclusion correction, ENA), which was subsequently used to evaluate the existence of a pattern of genetic divergence associated with an IBD phenomenon following the same procedure used for mitochondrial sequences.

The number of genetic groups (KSSR) was evaluated in Structure v2.3.4 [74,75,76] by performing multiple rounds of analysis using a hierarchical approach [77]. The decision to apply a hierarchical approach was based on the evolutionary history of H. ascotanensis. Previous studies have identified two genetically differentiated groups in the saltpan and proposed the occurrence of secondary contact between these groups [12, 30, 35]. This background suggests the presence of two groups with substantial genetic differentiation that could potentially mask the detection of additional populations. Therefore, given that Structure “accurately detects the uppermost hierarchical level of structure” [78], we opted to employ a hierarchical approach. In the first round, the complete dataset was analyzed to estimate the number of groups considering all the springs of the saltpan. Subsequently, the process was repeated recursively for each group detected in the previous round. The analysis was repeated until no population structure was detected in the estimated groups. In all rounds, a Bayesian analysis was performed using the MCMC algorithm that considered 600,000 iterations and a burn-in of 300,000 parameters. In addition, the results were obtained based on 15 runs. The structuring parameters were calculated for each K value using the admixture model and allele frequencies correlated between groups [74, 75]. The most likely number of populations (KSSR) was evaluated considering the value of the log likelihood of the observed data (LnP [D]) and the rate of second-order change in the logarithm of the data in different runs of K (ΔK) described in Evanno et al. [78]. ΔK was calculated using the Structure Harvester v0.6.94 online platform [79]. The results were summarized and compared in Clumpak v1.1 [80].

To estimate the contemporary gene flow between populations defined by Structure, a Bayesian analysis was performed using BayesAss v3.0.8 software [81]. This analysis was performed using 50 million iterations and a 30% burn-in period, with sampling every 5000 iterations. The mixing parameters for allelic frequencies, migration rates and inbreeding coefficients were defined as 0.30, 0.10 and 0.30, respectively. Ten independent runs were performed to examine the consistency in the results, and the convergence of each run was visualized with an MCMC trace plot. The Bayesian deviation measure was used to determine the run that showed the best fit of the model [82, 83].

Recent bottlenecks for the H. ascotanensis populations defined by Structure were estimated using Bottleneck v1.2 [84], based on the graphical analysis of Luikart et al. [85] and by Wilcoxon signed-rank tests (WSRs). Considering the nature of the dataset, the two-phase model (TPM) and the stepwise mutational model (SMM) were used to simulate the mutation-drift equilibria. For TPM, a two-phase mutation model was defined with 95% stepwise mutations, 5% multistep mutations and a variance (σ2g) of 12.

Geometric morphometric analyses

The external morphology of the individuals of H. ascotanensis was analyzed using geometric morphometrics, which allows an empirical comparison of the complex shapes presented by the shells. In Heleobia, a morphological change associated with sexual dimorphism has not been detected [12, 30, 56], so males and females were included in the analyses. The shell of each specimen was photographed from the ventral view using a digital camera (Moticam 5) coupled to a stereomicroscope (Motic SMZ-168). For the geometric morphometric analyses, the 310 individuals of the mitochondrial genetic study (COI) were analyzed, except for one individual from Spring 7 that presented damage to the peristome. The morphological study of the shell was carried out by classifying the individuals according to the populations defined by Structure.

Morphological analyses were carried out considering eight landmarks and 26 semi landmarks located in the extension of the entire shell. The landmarks were positioned at the apex of the shell and at the sutures present in the first, second and last turn of the shell. The semi landmarks were placed equidistant on the curves that connected the landmarks of the last turn of the shell and on the peristome (Fig. S1). The position of the landmarks was established to capture the main characteristics of the shape of the shell and functional traits associated with locomotion and feeding (opening of the shell). The location of the morphological landmarks was determined in the StereoMorph v1.6.7 package [86]. Using the geomorph v4.0.4 package [87], the landmark coordinates were aligned using a generalized Procrustes analysis (GPA). To determine the morphological differentiation associated with size, the GPA matrix (shell shape information) and the size of the log-transformed centroid (shell size information) were jointly analyzed. The shell morphospace was explored by principal component analysis (PCA) using the values of the morphological covariance matrix. After checking the assumption of data normality using the Shapiro-Wilk test in R, to quantify the degree of morphological variation associated with shell size between populations, a Procrustes ANOVA was performed on the morphological covariance matrix and the size of the centroid. The Bonferroni correction was applied to adjust the p-values of the multiple comparisons. To evaluate the variation in the shell shape between populations, a MANOVA was performed on the morphological covariance matrix. For the analyses, statistical significance was estimated through permutations (Residual Randomization in Permutation Procedures, RRPP; 10,000 iterations) [88,89,90] and a Pillai-Bartlett distribution was used. Finally, to determine whether individuals could be correctly assigned to each population using the shell shape (covariance matrix), a linear discriminant function analysis (LDA) was performed in R. The accuracy of the LDA classification was evaluated through a 10-fold cross-validation method, where the dataset is divided into 10 equally sized groups. Nine groups are utilized for training the model, while one group is set aside for validation.

Availability of data and materials

Sequences generated in this study have been deposited in GenBank under the Accession numbers OQ847095–OQ847404.

References

  1. Argollo J, Mourguiart P. Late Quaternary climate history of the Bolivian Altiplano. Quatern Int. 2000;72:37–51.

    Article  Google Scholar 

  2. Placzek CJ, Quade J, Patchett PJ. A 130ka reconstruction of rainfall on the Bolivian Altiplano. Earth Planet Sc Lett. 2013;363:97–108.

    Article  CAS  Google Scholar 

  3. Blard P-H, Lavé J, Farley KA, Fornari M, Jiménez N, Ramirez V. Late local glacial maximum in the central Altiplano triggered by cold and locally-wet conditions during the paleolake Tauca episode (17–15ka, Heinrich 1). Quaternary Sci Rev. 2009;28:3414–27.

    Article  Google Scholar 

  4. Lavenú A, Fornari M, Sebrier M. Existence de deux nouveaux épisodes lacustres quaternaires dans l’Altiplano Péruvo-Bolivien. Cahiers ORSTOM Série Géologie. 1984;14:103–14.

    Google Scholar 

  5. Fornari M, Risacher F, Féraud G. Dating of paleolakes in the central Altiplano of Bolivia. Palaeogeogr Palaeoclim Palaeoecol. 2001;172:269–82.

    Article  Google Scholar 

  6. Placzek C, Quade J, Patchett PJ. Geochronology and stratigraphy of late Pleistocene lake cycles on the southern Bolivian Altiplano: implications for causes of tropical climate change. Geol Soc Am Bull. 2006;118:515–32.

    Article  Google Scholar 

  7. Blard P-H, Sylvestre F, Tripati AK, Claude C, Causse C, Coudrain A, et al. Lake highstands on the Altiplano (Tropical Andes) contemporaneous with Heinrich 1 and the Younger Dryas: new insights from 14C, U–Th dating and δ18O of carbonates. Quaternary Sci Rev. 2011;30:3973–89.

    Article  Google Scholar 

  8. Martin LCP, Blard P-H, Lavé J, Condom T, Prémaillon M, Jomelli V, et al. Lake Tauca highstand (Heinrich Stadial 1a) driven by a southward shift of the Bolivian high. Sci Adv. 2018;4:eaar2514.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  9. Kroll O, Hershler R, Albrecht C, Terrazas EM, Apaza R, Fuentealba C, et al. The endemic gastropod fauna of Lake Titicaca: correlation between molecular evolution and hydrographic history. Ecol Evol. 2012;2:1517–30.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Wolff C, Albrecht C, Wilke T. Recovery from interglacial-related bottleneck likely triggered diversification of Lake Titicaca gastropod species flock. J Great Lakes Res. 2020;46:1199–206.

    Article  Google Scholar 

  11. Vila I, Morales P, Scott S, Poulin E, Véliz D, Harrod C, et al. Phylogenetic and phylogeographic analysis of the genus Orestias (Teleostei: Cyprinodontidae) in the southern Chilean Altiplano: the relevance of ancient and recent divergence processes in speciation. J Fish Biol. 2013;82:927–43.

    Article  PubMed  CAS  Google Scholar 

  12. Collado GA, Valladares MA, Méndez MA. Unravelling cryptic species of freshwater snails (Caenogastropoda, Truncatelloidea) in the Loa River basin. Atacama Desert Syst Biodivers. 2016;14:417–29.

    Article  Google Scholar 

  13. Collado GA, Vila I, Méndez MA. Monophyly, candidate species and vicariance in Biomphalaria snails (Mollusca: Planorbidae) from the southern Andean Altiplano. Zool Scr. 2011;40:613–22.

    Article  Google Scholar 

  14. Collado GA, Méndez MA. Microgeographic differentiation among closely related species of Biomphalaria (Gastropoda: Planorbidae) from the Andean Altiplano. Zool J Linnean Soc. 2013;169:640–52.

    Article  Google Scholar 

  15. Victoriano PF, Muñoz-Mendoza C, Sáez PA, Salinas HF, Muñoz-Ramírez C, Sallaberry M, et al. Evolution and conservation on top of the world: Phylogeography of the marbled water frog (Telmatobius marmoratus species complex; Anura, Telmatobiidae) in protected areas of Chile. J Hered. 2015;106:546–59.

    Article  PubMed  CAS  Google Scholar 

  16. Jurado-Rivera JA, Zapelloni F, Pons J, Juan C, Jaume D. Morphological and molecular species boundaries in the Hyalella species flock of Lake Titicaca (Crustacea: Amphipoda). Contrib Zool. 2020;89:353–72.

    Article  Google Scholar 

  17. Adamowicz SJ, Marinone MC, Menu-Marque S, Martin JW, Allen DC, Pyle MN, et al. The Hyalella (Crustacea: Amphipoda) species cloud of the ancient Lake Titicaca originated from multiple colonizations. Mol Phylogenet Evol. 2018;125:232–42.

    Article  PubMed  Google Scholar 

  18. Wade L. Chasing South America’s monsoon. Science. 2014;346:1042–3.

    Article  PubMed  CAS  Google Scholar 

  19. Betancourt JL, Latorre C, Rech JA, Quade J, Rylander KA. A 22,000-year record of monsoonal precipitation from northern Chile’s Atacama Desert. Science. 2000;289:1542–6.

    Article  PubMed  CAS  Google Scholar 

  20. Zhou J, Lau K-M. Does a monsoon climate exist over South America? J Clim. 1998;11:1020–40.

    Article  Google Scholar 

  21. Krishnamurthy V, Misra V. Observed ENSO teleconnections with the south American monsoon system. Atmos Sci Lett. 2010;11:7–12.

    Article  Google Scholar 

  22. Wang B, Liu J, Kim H-J, Webster PJ, Yim S-Y, Xiang B. Northern hemisphere summer monsoon intensified by mega-El Niño/southern oscillation and Atlantic multidecadal oscillation. Proc National Acad Sci. 2013;110:5347–52.

    Article  CAS  Google Scholar 

  23. Marengo JA, Liebmann B, Grimm AM, Misra V, Dias PLS, Cavalcanti IFA, et al. Recent developments on the south American monsoon system. Int J Climatol. 2012;32:1–21.

    Article  Google Scholar 

  24. Carvalho LMV, Cavalcanti IFA. The monsoons and climate change, observations and modeling. Springer Clim; 2015. p. 121–48.

    Google Scholar 

  25. Rojas M, Arias PA, Flores-Aqueveque V, Seth A, Vuille M. The South American Monsoon Variability over the Last Millennium in CMIP5/PMIP3 simulations. Clim Past Discuss. 2015;11:5651–81.

    Google Scholar 

  26. González-Pinilla FJ, Latorre C, Rojas M, Houston J, Rocuant MI, Maldonado A, et al. High- and low-latitude forcings drive Atacama Desert rainfall variations over the past 16,000 years. Sci Adv. 2021;7:eabg1333.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Morales MS, Christie DA, Villalba R, Argollo J, Pacajes J, Silva JS, et al. Precipitation changes in the south American Altiplano since 1300 AD reconstructed by tree-rings. Clim Past. 2012;8:653–66.

    Article  Google Scholar 

  28. Wilcox AC, Escauriaza C, Agredano R, Mignot E, Zuazo V, Otárola S, et al. An integrated analysis of the march 2015 Atacama floods. Geophys Res Lett. 2016;43:8035–43.

    Article  Google Scholar 

  29. Cruz-Jofré F, Morales P, Vila I, Esquer-Garrigos Y, Hugueny B, Gaubert P, et al. Geographical isolation and genetic differentiation: the case of Orestias ascotanensis (Teleostei: Cyprinodontidae), an Andean killifish inhabiting a highland salt pan. Biol J Linn Soc. 2016;117:747–59.

    Article  Google Scholar 

  30. Valladares MA, Méndez MA, Collado GA. Influenced but not determined by historical events: genetic, demographic and morphological differentiation in Heleobia ascotanensis from the Chilean Altiplano. Peerj. 2018;6:e5802.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Morales P, Vila I, Poulin E. Genetic structure in remnant populations of an endangered cyprinodontid fish, Orestias ascotanensis, endemic to the Ascotán salt pan of the Altiplano. Conserv Genet. 2011;12:1639–43.

    Article  Google Scholar 

  32. Johnson SG. Age, phylogeography and population structure of the microendemic banded spring snail, Mexipyrgus churinceanus. Mol Ecol. 2005;14:2299–311.

    Article  PubMed  CAS  Google Scholar 

  33. Collado GA, Méndez MA. Los taxa nominales de moluscos descritos por Courty del Salar de Ascotán, Altiplano chileno. Rev Chil Hist Nat. 2012;85:233–5.

    Article  Google Scholar 

  34. Courty G. Explorations géologiques dans l’Amerique du Sud. In: Mission Scientifique de G. de Crèqui Montefort et E. Sénéchal de la Grange. 1st ed. Impremiere Nationale; 1907.

    Google Scholar 

  35. Collado GA, Valladares MA, Méndez MA. Hidden diversity in spring snails from the Andean Altiplano, the second highest plateau on earth, and the Atacama Desert, the driest place in the world. Zoöl Stud. 2013;52:50.

    Article  Google Scholar 

  36. Van Leeuwen CH, Huig N, Van der Velde G, Van Alen TA, Wagemaker CA, Sherman CD, et al. How did this snail get here? Several dispersal vectors inferred for an aquatic invasive species. Freshw Biol. 2013;58:88–99.

    Article  Google Scholar 

  37. Van Leeuwen CHA, Tollenaar ML, Klaassen M. Vector activity and propagule size affect dispersal potential by vertebrates. Oecologia. 2012;170:101–9.

    Article  PubMed  PubMed Central  Google Scholar 

  38. van Leeuwen CHA, van der Velde G. Prerequisites for flying snails: external transport potential of aquatic snails by waterbirds. Freshw Sci. 2012;31:963–72.

    Article  Google Scholar 

  39. Saito T, Tatani M, Odaya Y, Chiba S. Direct evidence for intercontinental dispersal of a snail via a bird. Ecography. 2023;10:e06771.

    Article  Google Scholar 

  40. Collado GA. Heleobia ascotanensis (Courty, 1907). 11o Proceso de Clasificación de Especies Ministerio del Medio Ambiente, Gobierno de Chile. 2014. Available at: http://www.mma.gob.cl (accessed on April 10, 2023).

  41. Moline AB, Shuster SM, Hendrickson DA, Marks JC. Genetic variation in a desert aquatic snail (Nymphophilus minckleyi) from Cuatro Ciénegas, Coahuila, Mexico. Hydrobiologia. 2004;522:179–92.

    Article  Google Scholar 

  42. Sei M, Lang BK, Berg DJ. Genetic and community similarities are correlated in endemic-rich springs of the northern Chihuahuan Desert. Glob Ecol Biogeogr. 2009;18:192–201.

    Article  Google Scholar 

  43. Trajanovski S, Albrecht C, Schreiber K, Schultheiß R, Stadler T, Benke M, et al. Testing the spatial and temporal framework of speciation in an ancient lake species flock: the leech genus Dina (Hirudinea: Erpobdellidae) in Lake Ohrid. Biogeosciences. 2010;7:3387–402.

    Article  Google Scholar 

  44. Murphy NP, Adams M, Austin AD. Independent colonization and extensive cryptic speciation of freshwater amphipods in the isolated groundwater springs of Australia’s great Artesian Basin. Mol Ecol. 2009;18:109–22.

    Article  PubMed  CAS  Google Scholar 

  45. Bos DH, Gopurenko D, Williams RN, DeWoody AJ. Inferring population history and demography using microsatellites, mitochondrial DNA, and major histocompatibility complex (MHC) genes. Evolution. 2008;62:1458–68.

    Article  PubMed  CAS  Google Scholar 

  46. Castella V, Ruedi M, Excoffier L. Contrasted patterns of mitochondrial and nuclear structure among nursery colonies of the bat Myotis myotis. J Evol Biol. 2001;14:708–20.

    Article  Google Scholar 

  47. Valladares MA, Fabres AA, Collado GA, Sáez PA, Méndez MA. Coping with dynamism: Phylogenetics and phylogeographic analyses reveal cryptic diversity in Heleobia snails of Atacama saltpan. Chile Frontiers Ecol Evol. 2022;10:869626.

    Article  Google Scholar 

  48. Keller B, Soto D. Hydrogeologic influences on the preservation of Orestias ascotanensis (Teleostei Cyprinodontidae), in Salar de Ascotán, northern Chile. Rev Chil Hist Nat. 1998;71:147–56.

    Google Scholar 

  49. Collado GA, Fuentealba CG, Cazzaniga NJ, Valladares MA. Assessing biodiversity within the range of Heleobia chimbaensis (Caenogastropoda: Cochliopidae) on the Atacama Desert coast. Syst Biodivers. 2020;18:708–19.

    Article  Google Scholar 

  50. Murphy NP, Breed MF, Guzik MT, Cooper SJB, Austin AD. Trapped in desert springs: phylogeography of Australian desert spring snails. J Biogeogr. 2012;39:1573–82.

    Article  Google Scholar 

  51. Murphy NP, Guzik MT, Wilmer JW. The influence of landscape on population structure of four invertebrates in groundwater springs. Freshw Biol. 2010;55:2499–509.

    Article  Google Scholar 

  52. Liu H-P, Hershler R. Phylogeography of an endangered Western north American springsnail. Conserv Genet. 2012;13:299–305.

    Article  CAS  Google Scholar 

  53. Cook CN, Sgrò CM. Poor understanding of evolutionary theory is a barrier to effective conservation management. Conserv Lett. 2019;12.

  54. Buerki S, Callmander MW, Bachman S, Moat J, Labat J-N, Forest F. Incorporating evolutionary history into conservation planning in biodiversity hotspots. Philos Trans R Soc B Biol Sci. 2015;370:20140014.

    Article  Google Scholar 

  55. Collado GA, Torres-Díaz C, Vidal MA, Valladares MA. Genetic diversity, morphometric characterization, and conservation reassessment of the critically endangered freshwater snail, Heleobia atacamensis, in the Atacama saltpan, Northern Chile. Biology. 2023;12:791.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Collado GA, Chihuailaf E, Muñoz-Herrera N, Contreras M, Novoa F, Valladares MA. Reproductive aspects of the poorly known and critically endangered freshwater snail Heleobia atacamensis (Gastropoda: Truncatelloidea). Peerj. 2021;9:e11550.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  58. Katoh K, Rozewicki J, Yamada KD. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform. 2017;20:bbx108.

    Google Scholar 

  59. Fabres AA, Valladares MA, Sáez PA, Collado GA, Pastenes L, Méndez MA. Novel microsatellite markers for an endangered freshwater snail, Heleobia atacamensis (Caenogastropoda: Cochliopidae), from the Atacama saltpan. Molluscan Res. 2020;40:231–5.

    Article  Google Scholar 

  60. Oosterhout CV, Hutchinson WF, Wills DPM, Shipley P. micro-checker: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4:535–8.

    Article  Google Scholar 

  61. Rozas J, Ferrer-Mata A, Sánchez-DelBarrio JC, Guirao-Rico S, Librado P, Ramos-Onsins SE, et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol Biol Evol. 2017;34:3299–302.

    Article  PubMed  CAS  Google Scholar 

  62. Bandelt HJ, Forster P, Rohl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.

    Article  PubMed  CAS  Google Scholar 

  63. Leigh JW, Bryant D. Popart: full-feature software for haplotype network construction. Methods Ecol Evol. 2015;6:1110–6.

    Article  Google Scholar 

  64. Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and windows. Mol Ecol Resour. 2010;10:564–7.

    Article  PubMed  Google Scholar 

  65. Rousset F. Genetic differentiation and estimation of gene flow from F -statistics under isolation by distance. Genetics. 1997;145:1219–28.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  66. Oksanen J, Simpson G, Blanchet FG, Kindt R, Legendre P, Minchin P, et al. Vegan community ecology package version 2. 2022.

  67. R Development Core Team. R: A Language and Environment for Statistical Computing, Vienna, Austria. ISBN 3–900051–07-0. R Foundation for Statistical Computing. 2019.

  68. Rousset F. genepop’007: a complete re-implementation of the genepop software for windows and Linux. Mol Ecol Resour. 2008;8:103–6.

    Article  PubMed  Google Scholar 

  69. Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in excel. Population genetic software for teaching and research—an update. Bioinformatics. 2012;28:2537–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  70. Adamack AT, Gruber B. PopGenReport: simplifying basic population genetic analyses in R. Methods Ecol Evol. 2014;5:384–7.

    Article  Google Scholar 

  71. Alho JS, Välimäki K, Merilä J. Rhh: an R extension for estimating multilocus heterozygosity and heterozygosity–heterozygosity correlation. Mol Ecol Resour. 2010;10:720–2.

    Article  PubMed  Google Scholar 

  72. Amos W, Wilmer JW, Fullard K, Burg TM, Croxall JP, Bloch D, et al. The influence of parental relatedness on reproductive success. Proc Royal Soc Lond Ser B Biol Sci. 2001;268:2021–7.

    Article  CAS  Google Scholar 

  73. Chapuis M-P, Estoup A. Microsatellite null alleles and estimation of population differentiation. Mol Biol Evol. 2007;24:621–31.

    Article  PubMed  CAS  Google Scholar 

  74. Falush D, Stephens M, Pritchard JK. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003;164:1567–87.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  75. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  76. Falush D, Stephens M, Pritchard JK. Inference of population structure using multilocus genotype data: dominant markers and null alleles. Mol Ecol Notes. 2007;7:574–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  77. Vähä J-P, Erkinaro J, Niemelä E, Primmer CR. Life-history and habitat features influence the within-river genetic structure of Atlantic salmon. Mol Ecol. 2007;16:2638–54.

    Article  PubMed  Google Scholar 

  78. Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software structure: a simulation study. Mol Ecol. 2005;14:2611–20.

    Article  PubMed  CAS  Google Scholar 

  79. Earl DA, vonHoldt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4:359–61.

    Article  Google Scholar 

  80. Kopelman NM, Mayzel J, Jakobsson M, Rosenberg NA, Mayrose I. Clumpak: a program for identifying clustering modes and packaging population structure inferences across K. Mol Ecol Resour. 2015;15:1179–91.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  81. Wilson GA, Rannala B. Bayesian inference of recent migration rates using multilocus genotypes. Genetics. 2003;163:1177–91.

    Article  PubMed  PubMed Central  Google Scholar 

  82. Meirmans PG. Nonconvergence in Bayesian estimation of migration rates. Mol Ecol Resour. 2014;14:726–33.

    Article  PubMed  Google Scholar 

  83. Yamashiro T, Yamashiro A, Inoue M, Maki M. Genetic diversity and divergence in populations of the threatened grassland perennial Vincetoxicum atratum (Apocynaceae-Asclepiadoideae) in Japan. J Hered. 2016;107:455–62.

    Article  PubMed  Google Scholar 

  84. Piry S, Luikart G, Cornuet J-M. BOTTLENECK: a computer program for detecting recent reductions in the effective size using allele frequency data. J Hered. 1999;90:502–3.

    Article  Google Scholar 

  85. Luikart G, Allendorf F, Cornuet J-M, Sherwin W. Distortion of allele frequency distributions provides a test for recent population bottlenecks. J Hered. 1998;89:238–47.

    Article  PubMed  CAS  Google Scholar 

  86. Olsen AM, Westneat MW. StereoMorph: an R package for the collection of 3D landmarks and curves using a stereo camera set-up. Methods Ecol Evol. 2015;6:351–6.

    Article  Google Scholar 

  87. Adams DC, Otárola-Castillo E. Geomorph: an r package for the collection and analysis of geometric morphometric shape data. Methods Ecol Evol. 2013;4:393–9.

    Article  Google Scholar 

  88. Adams DC, Collyer ML. Permutation tests for phylogenetic comparative analyses of high-dimensional shape data: what you shuffle matters. Evolution. 2015;69:823–9.

    Article  PubMed  Google Scholar 

  89. Collyer ML, Sekora DJ, Adams DC. A method for analysis of phenotypic change for phenotypes described by high-dimensional data. Heredity. 2015;115:357–65.

    Article  PubMed  CAS  Google Scholar 

  90. Adams DC, Collyer ML. Phylogenetic ANOVA: group-clade aggregation, biological challenges, and a refined permutation procedure. Evolution. 2018;72:1204–15.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

MAV acknowledges to ANID-FONDECYT Postdoctorado 3220352. MAM acknowledges to ANID-FONDECYT 1200419, 1221214, ANID PIA/BASAL FB0002 and CHIC-ANID PIA/BASAL PFB210018. The authors thank Elizabeth Chihuailaf and Natalia C. Muñoz-Herrera for their collaboration in the field work. We also thank Paola Sáez and Maira Anderson for their assistance throughout the study. This research was partially supported by the supercomputing infrastructure of the NLHPC (ECM-02).

Funding

This work was supported by the Minera El Abra Ltda., and Centro de Ecología Aplicada Ltda., Chile.

Author information

Authors and Affiliations

Authors

Contributions

MAV and MAM conceptualized the study. MAV and AAF designed the experiments. MAV and GAC performed the sampling. AAF and FS-R performed the experiments. MAV, AAF and FS-R analyzed the data. MAM obtained funding for the original project idea and contributed in work proposal. MAV wrote the manuscript. All authors reviewed and approved the final manuscript.

Corresponding author

Correspondence to Marco A. Méndez.

Ethics declarations

Ethics approval and consent to participate

This research was carried out following approved standard procedures for working with invertebrate animals. The field work was carried out under the following permit from the Undersecretary of Fisheries and Aquaculture (SUBPESCA): Exempt Resolution No. 2604/19.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Characteristics of the 10 microsatellites used in the study and significant deviations of Hardy-Weinberg Equilibrium by populations after Bonferroni corrections (p < 0.05). Table S2. Estimated contemporary migration rate for Heleobia populations of the Ascotán Saltpan obtained using 10 microsatellite loci. Standard deviation is shown in parentheses. The self-recruitment rate is indicated along the diagonal. Table S3. List of primers sequences used in the study for the mitochondrial and microsatellites markers. Fig. S1. Geometric morphometrics in Heleobia ascotanensis. The landmarks and semilandmarks used in the morphological study are presented.

Additional file 2: Table S1.

Microsatellite allele data for 390 individuals genotyped at 10 loci.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Valladares, M.A., Fabres, A.A., Sánchez-Rodríguez, F. et al. Population structure and microscale morphological differentiation in a freshwater snail from the Chilean Altiplano. BMC Ecol Evo 24, 5 (2024). https://doi.org/10.1186/s12862-023-02196-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12862-023-02196-w

Keywords