Comparative population genetics of two dominant plant species of high Andean wetlands reveals complex evolutionary histories and conservation perspectives in Chile’s Norte Chico

Spatial autocorrelation signatures of ecological determinants on plant community characteristics in high Andean wetlands

Understanding how biological communities are shaped is a central tenet of community ecology. Recent evidence highlights the potential of decoupling diversity spatial autocorrelation into its positive and negative components to reveal community assembly processes that would otherwise remain undetected, as well as to improve understanding of their impacts on different facets of diversity. Yet, such approaches have only been implemented to investigate the effects of a few assembly drivers on a small number of diversity components. Here, we used high Andean wetland plant communities over a strong latitudinal gradient to investigate the effects of various ecological factors on spatial autocorrelation patterns of nine community metrics with different informative values, including measures of richness, dominance, evenness and beta-diversity. By combining Moran's Eigenvector Maps, partial least squares structural equation modeling, and regression analyses, we revealed two groups of community parameters presenting contrasting spatial patterns due to specific sensitivities to ecological factors. While environmental variation and wetland connectivity increased positive spatial autocorrelation in richness and dominance-related parameters, species co-occurrence promoted negative spatial autocorrelation in evenness-related parameters. These results offer new insights regarding both how ecological processes affect species assembly, as well as the information captured by classical taxonomic parameters.

Genome-wide genetic diversity yields insights into genomic responses of candidate climate-selected loci in an Andean wetland plant

Assessing population evolutionary potential has become a central tenet of conservation biology. Since adaptive responses require allelic variation at functional genes, consensus has grown that genetic variation at genes under selection is a better surrogate for adaptive evolutionary potential than neutral genetic diversity. Although consistent with prevailing theory, this argument lacks empirical support and ignores recent theoretical advances questioning the very concept of neutral genetic diversity. In this study, we quantified genome-wide responses of single nucleotide polymorphism loci linked to climatic factors over a strong latitudinal gradient in natural populations of the high Andean wetland plant, Carex gayana, and then assessed whether genetic variation of candidate climate-selected loci better predicted their genome-wide responses than genetic variation of non-candidate loci. Contrary to this expectation, genomic responses of climate-linked loci only related significantly to environmental variables and genetic diversity of non-candidate loci. The effects of genome-wide genetic diversity detected in this study may be a result of either the combined influence of small effect variants or neutral and demographic factors altering the adaptive evolutionary potential of C. gayana populations. Regardless of the processes involved, our results redeem genome-wide genetic diversity as a potentially useful indicator of population adaptive evolutionary potential.

Partitioning genetic and species diversity refines our understanding of species-genetic diversity relationships

Disentangling the origin of species-genetic diversity correlations (SGDCs) is a challenging task that provides insight into the way that neutral and adaptive processes influence diversity at multiple levels. Genetic and species diversity are comprised by components that respond differently to the same ecological processes. Thus, it can be useful to partition species and genetic diversity into their different components to infer the mechanisms behind SGDCs. In this study, we applied such an approach using a high-elevation Andean wetland system, where previous evidence identified neutral processes as major determinants of the strong and positive covariation between plant species richness and AFLP genetic diversity of the common sedge Carex gayana. To tease apart putative neutral and non-neutral genetic variation of C. gayana, we identified loci putatively under selection from a dataset of 1,709 SNPs produced using restriction site-associated DNA sequencing (RAD-seq). Significant and positive relationships between local estimates of genetic and species diversities (α-SGDCs) were only found with the putatively neutral loci datasets and with species richness, confirming that neutral processes were primarily driving the correlations and that the involved processes differentially influenced local species diversity components (i.e., richness and evenness). In contrast, SGDCs based on genetic and community dissimilarities (β-SGDCs) were only significant with the putative non-neutral datasets. This suggests that selective processes influencing C. gayana genetic diversity were involved in the detected correlations. Together, our results demonstrate that analyzing distinct components of genetic and species diversity simultaneously is useful to determine the mechanisms behind species-genetic diversity relationships.

Landscape connectivity among remnant populations of guanaco (Lama guanicoe Müller, 1776) in an arid region of Chile impacted by global change

Connectivity between populations plays a key role in the long-term persistence of species in fragmented habitats. This is of particular concern for biodiversity preservation in drylands, since water limited landscapes are typically characterized by little suitable habitat cover, high habitat fragmentation, harsh matrices, and are being rapidly degraded at a global scale. In this study, we modelled landscape connectivity between 11 guanaco Lama guanicoe populations in Chile’s arid Norte Chico, a region that supports the last remnant coastal populations of this emblematic herbivore indigenous to South America. We produced a habitat suitability model to derive a regional surface resistance map, and used circuit theory to map functional connectivity, investigate the relative isolation between populations, and identify those that contribute most to the patch connectivity network. Predicted suitable habitat for L. guanicoe represented about 25% of the study region (i.e., 29,173 km²) and was heterogeneously distributed along a continuous stretch along the Andes, and discontinuous patches along the coast. As a result, we found that high connectivity current flows in the mid and high Andes formed a wide, continuous connectivity corridor, enabling connectivity between all high Andean populations. Coastal populations, in contrast, were more isolated. These groups demonstrate no inter-population connectivity between themselves, only with higher altitude populations, and for two of them, animal movement was linked to the effectiveness of wildlife crossings along the Pan-American highway. Our results indicate that functional connectivity is an issue of concern for L. guanicoe in Chile’s Norte Chico, implying that future conservation and management plans should emphasize strategies aimed at conserving functional connectivity between coastal and Andean populations, as well as the protection of habitat patches likely to act as stepping stones within the connectivity network.