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

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

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.

Diversity, Taxonomic Novelty, and Encoded Functions of Salar de Ascotán Microbiota, as Revealed by Metagenome-Assembled Genomes

Salar de Ascotán is a high-altitude arsenic-rich salt flat exposed to high ultraviolet radiation in the Atacama Desert, Chile. It hosts unique endemic flora and fauna and is an essential habitat for migratory birds, making it an important site for conservation and protection. However, there is limited information on the resident microbiota's diversity, genomic features, metabolic potential, and molecular mechanisms that enable it to thrive in this extreme environment. We used long- and short-read metagenomics to investigate the microbial communities in Ascotán's water, sediment, and soil. Bacteria predominated, mainly Pseudomonadota, Acidobacteriota, and Bacteroidota, with a remarkable diversity of archaea in the soil. Following hybrid assembly, we recovered high-quality bacterial (101) and archaeal (6) metagenome-assembled genomes (MAGs), including representatives of two putative novel families of Patescibacteria and Pseudomonadota and two novel orders from the archaeal classes Halobacteriota and Thermoplasmata. We found different metabolic capabilities across distinct lineages and a widespread presence of genes related to stress response, DNA repair, and resistance to arsenic and other metals. These results highlight the remarkable diversity and taxonomic novelty of the Salar de Ascotán microbiota and its rich functional repertoire, making it able to resist different harsh conditions. The highly complete MAGs described here could serve future studies and bioprospection efforts focused on salt flat extremophiles, and contribute to enriching databases with microbial genome data from underrepresented regions of our planet.

Genetic Diversity, Morphometric Characterization, and Conservation Reassessment of the Critically Endangered Freshwater Snail, Heleobia atacamensis, in the Atacama Saltpan, Northern Chile

Evaporitic ecosystems of the Atacama Desert contain a rich endemic fauna, including mollusk species. A recent study performed in the freshwater snail Heleobia atacamensis, endemic to the Atacama Saltpan, revealed a strong interdependence of genetic patterns with climatic fluctuations and landscape physiography. The species is currently listed as Critically Endangered at regional scale and as Data Deficient on the International Union for Conservation of Nature (IUCN) Red List. Here, we studied genetic diversity and demographic history of several populations of the species occurring on a connectivity gradient, including snails from new peripherical localities (Peine and Tilomonte), which were compared with topotype specimens. In addition, we reassessed the conservation status using the IUCN Red List categories and criteria considering species-specific idiosyncrasy. Phylogenetic and phylogeographical analyses indicated that snails from Peine and Tilomonte belong to H. atacamensis. We discovered significant differentiation in shell morphology, which was generally greater in geographically isolated populations. We also inferred six genetic clusters and a demographic expansion congruent with the wet periods that occurred at the end of the Pleistocene. Considering the highest risk category obtained, H. atacamensis was reassessed as Endangered at regional scale. Future conservation plans should consider the genetic assemblages as conservation units.

Phylogeny of Telmatobius marmoratus complex (Anura, Telmatobiidae) reveals high cryptic diversity in the Andean Altiplano

Telmatobius is the most diverse group of anurans in the Andean Altiplano (highlands) Morphologically, these amphibians have a generally conserved morphology but in turn present large intraspecific variation, which has led to a complex taxonomy and systematics. T. marmoratus has the widest distribution of the genus and forms a complex composed of at least two Telmatobius species. Partial systematic studies based on molecular evidence reveal the existence of three lineages with a complex spatial distribution. However, these studies did not include the entire distribution of T. marmoratus. Our study aims to reassess the current systematic scenario including the complete distribution of the complex. For this, we used a multilocus approach based on mitochondrial (16S, Cytb) and nuclear (RAG1-1, BFIB) DNA sequences to build a phylogenetic hypothesis based on Bayesian inference, maximum likelihood and maximum parsimony. Subsequently, we performed single-locus (ABGD and PTP) and multilocus (STACEY) species delimitation analyses to verify the diversity of nominal species within the complex. The analyses suggest seven non-sibling lineages and 6-10 candidate species within the marmoratus complex. Only one of the two lineages restricted to the central northern plateau correspond to T. marmoratus sensu stricto. South-central marbled water frogs belong to completely new lineages closer to T. gigas and T. culeus, evidencing the polyphyletic condition of the marmoratus complex. The findings of several sympatric lineages in some localities reveal a complex history of ancient water connections in south-central Altiplano.

The definitive rediscovery of Telmatobiushalli (Anura, Telmatobiidae) at its historic type locality and its synonymy with T.dankoi and T.vilamensis

Telmatobiushalli was the first representative of its genus to be described exclusively for Chile, yet for 85 years no new individuals could be located due to the vagueness with which its type locality was described. The type series was collected by one of the members of the International High Altitude Expedition to Chile (IHAEC) of 1935. Recently, three studies successively claimed to have located the type locality in different places. The third study proved, according to the chronicles of the IHAEC, that the actual locality is Miño, at the origin of the Loa River, where currently there are no published records of Telmatobius. In this study, additional documentary antecedents and graphic material are provided that corroborate that Miño is indeed the type locality of T.halli. Additionally, the recently rediscovered Telmatobius population from Miño and the environment it inhabits are described. The external characteristics of the frogs are consistent with the description of T.halli. Furthermore, molecular phylogenetic analyses were performed that showed that T.halli, T.dankoi, and T.vilamensis, all known only from their type localities in Chile, comprise a clade without internal resolution. A detailed comparison of the diagnoses of the three species revealed that the few phenotypic differences between these taxa were based on characteristics that vary widely within and between populations of the genus, hence their conspecificity is proposed. The implications of this synonymy for the taxonomy, biogeography, and conservation of the Telmatobius from the extreme south of its distribution in Chile are discussed.

Genome sequencing and transcriptomic analysis of the Andean killifish Orestias ascotanensis reveals adaptation to high-altitude aquatic life

Orestias ascotanensis (Cyprinodontidae) is a teleost pupfish endemic to springs feeding into the Ascotan saltpan in the Chilean Altiplano (3,700 m.a.s.l.) and represents an opportunity to study adaptations to high-altitude aquatic environments. We have de novo assembled the genome of O. ascotanensis at high coverage. Comparative analysis of the O. ascotanensis genome showed an overall process of contraction, including loss of genes related to G-protein signaling, chemotaxis and signal transduction, while there was expansion of gene families associated with microtubule-based movement and protein ubiquitination. We identified 818 genes under positive selection, many of which are involved in DNA repair. Additionally, we identified novel and conserved microRNAs expressed in O. ascotanensis and its closely-related species, Orestias gloriae. Our analysis suggests that positive selection and expansion of genes that preserve genome stability are a potential adaptive mechanism to cope with the increased solar UV radiation to which high-altitude animals are exposed to.

Phylogeography of high Andean killifishes Orestias (Teleostei: Cyprinodontidae) in Caquena and Lauca sub-basins of the Altiplano (Chile): mitochondrial and nuclear analysis of an endangered fish

From the early Miocene, the uplift of the Andes Mountains, intense volcanic activity and the occurrence of successive periods of dryness and humidity would have differentially influenced the modification of Altiplano watersheds, and consequently the evolutionary history of the taxa that live there. We analyzed Orestias populations from the Caquena and Lauca Altiplanic sub-basins of northern Chile to determine their genetic differentiation and relationship to their geographical distribution using mitochondrial (D-loop) and nuclear (microsatellite) molecular markers and to reconstruct its biogeographic history on these sub-basins. The results allowed reconstructing and reevaluating the evolutionary history of the genus in the area; genic diversity and differentiation together with different founding genetic groups suggest that Orestias have been spread homogeneously in the study area and would have experienced local disturbances that promoted isolation and diversification in restricted zones of their distribution.

Influenced but not determined by historical events: genetic, demographic and morphological differentiation in Heleobia ascotanensis from the Chilean Altiplano

In the present study, we focus on the phylogeographic pattern, demographic history and morphological differentiation of Heleobia ascotanensis, a freshwater gastropod restricted to the Ascotán saltpan in the Chilean Altiplano. The current distribution of the species is limited to twelve isolated or partially isolated springs that were affected by transitions between humid and arid periods during last glaciations. The genetic analysis of 322 specimens showed that H. ascotanensis is subdivided into three genetically divergent populations, with low and moderate degrees of historical gene flow among them and incipient morphological differentiation as a consequence of genetic and geographical isolation. Molecular analyses revealed different demographic histories among populations which seem to respond independently to climatic events, probably due to an environmental imposition and idiosyncratic strategies developed to cope with water availability. The results of this study and co-distributed taxa support the hypothesis that contemporary and historical events have influenced microevolutionary differentiation of these snails, although there is a need to complement further information to predict genetic or morphological divergence at microgeographic scale.

First steps towards assessing the evolutionary history and phylogeography of a widely distributed Neotropical grassland bird (Motacillidae: Anthus correndera)

Grasslands in southern South America are extensive ecosystems which harbor a unique biodiversity; however, studies on the evolution of their taxa are scarce. Here we studied the phylogeography and population history of the Correndera Pipit (Anthus correndera), a grassland specialist bird with a large breeding distribution in southern South America, with the goals of investigating its phylogeographic history and relate it to the historical development of South American grasslands. The mitochondrial NADH dehydrogenase subunit II gene (ND2) was sequenced in 66 individuals from 19 localities and the intron 9 of the sex-linked gene for aconitase (ACOI9) was sequenced from a subset of those individuals, including all five subspecies of A. correndera, as well as the closely related A. antarcticus. Phylogenetic analysis revealed two distinct lineages within the complex: the first (A) corresponding to Andean subspecies A. c. calcaratus and A. c. catamarcae and the second (B) including birds traditionally assigned to A. c. correndera, A. c. chilensis, A. c. grayi and some individuals of A. c. catamarcae. A. antarcticus is nested within this second lineage. These results were also supported by evidence of niche divergence for variables associated with precipitation. The oldest split between clade A and B was estimated at c. 0.37 Mya, during the middle Pleistocene. Species distribution models for the present and the Last Glacial Maximum (LGM) suggest that grassland areas in southern South America remained relatively stable, in contrast to the general view of a reduction in grassland cover in South America since the LGM. Recent divergences and low phylogeographic structure (for lowland vs. highland geographic groups, intra-population genetic variance was greater than inter-groups; e.g., for ACOI9: 95.47% and ND2: 51.51% respectively), suggest widespread gene flow between lowland populations.

Morphological and taxonomic descriptions of a new genus and species of killifishes (Teleostei: Cyprinodontiformes) from the high Andes of northern Chile

A new genus and species, Pseudorestias lirimensis, is described from the southern part of the Chilean Altiplano. While sharing several characters that clearly align the new species with Orestias, this new fish is characterized by numerous autapomorphies: the Meckel cartilage is a continuous cartilage that broadly expands posteriorly (in large specimens, it keeps its anterior part and is resorbed posteriorly), the basibranchials are fused into one long element, the second pharyngobranchial is not displaced dorsally over pharyngobranchial tooth plate 3+4, but they are aligned, the anterior and posterior ceratohyals are closely articulated keeping a scarce amount of cartilage between both bones and ventral to them, ossified middle and distal dorsal radials are present in females as well as ossified middle and distal anal radials. Pseudorestias lirimensis presents strong sexual dimorphism associated to size. Females are almost twice as large and long than males, neuromast lines are absent in males, a mesethmoid is present in males, squamation on head is reduced in males, and ossified middle and distal radial of dorsal fin are cartilaginous in males. Pseudorestias and Orestias are suggested as the sole members of the tribe Orestiini. A list of characters diagnosing the tribe is provided. The presence of the new genus is interpreted as a possible result of the ecosystem isolation where the fish is living from surrounding basins-as early as possibly from the Miocene-Pliocene times-and its physical and chemical characteristics. Small populations, living conditions, small habitat, and reduced distribution make this species a strong candidate to be considered critically endangered, a situation already established for all other Chilean species living in the Altiplano. There is high probability it will become extinct due to water demands and climate change in the region.

Pattern of genetic differentiation of an incipient speciation process: The case of the high Andean killifish Orestias

During the Pleistocene and Holocene, the southwest Andean Altiplano (17°-22°S) was affected by repeated fluctuations in water levels, high volcanic activity and major tectonic movements. In the early Holocene the humid Tauca phase shifted to the arid conditions that have lasted until the present, producing endorheic rivers, lakes, lagoons and wetlands. The endemic fish Orestias (Cyprinodontidae) represents a good model to observe the genetic differentiation that characterizes an incipient speciation process in allopatry since the morphospecies described inhabit a restricted geographic area, with present habitat fragmentation. The genetic diversity and population structure of four endemic morphospecies of Orestias (Cyprinodontidae) found in the Lauca National Park (LNP) analyzed with mitochondrial markers (Control Region) and eight microsatellites, revealed the existence of genetic groups that matches the fragmentation of these systems. High values of genetic and phylogeographic differentiation indices were observed between Chungará Lake and Piacota lagoon. The group composed of the Lauca River, Copapujo and Chuviri wetlands sampling sites showed a clear signal of expansion, with a star-like haplotype network. Levels of genetic differentiation were lower than in Chungará and Piacota, suggesting that these localities would have differentiated after the bottlenecks linked to the collapse of Parinacota volcano. The Parinacota sample showed a population signal that differed from the other localities revealing greater genetic diversity and a disperse network, presenting haplotypes shared with other LNP localities. A mixing pattern of the different genetic groups was evident using the microsatellite markers. The chronology of the vicariance events in LNP may indicate that the partition process of the Orestias populations was gradual. Considering this, and in view of the genetic results, we may conclude that the morphospecies from LNP are populations in ongoing differentiation process.

Microbial Activity Response to Solar Radiation across Contrasting Environmental Conditions in Salar de Huasco, Northern Chilean Altiplano

In high altitude environments, extreme levels of solar radiation and important differences of ionic concentrations over narrow spatial scales may modulate microbial activity. In Salar de Huasco, a high-altitude wetland in the Andean mountains, the high diversity of microbial communities has been characterized and associated with strong environmental variability. Communities that differed in light history and environmental conditions, such as nutrient concentrations and salinity from different spatial locations, were assessed for bacterial secondary production (BSP, ³H-leucine incorporation) response from short-term exposures to solar radiation. We sampled during austral spring seven stations categorized as: (a) source stations, with recently emerged groundwater (no-previous solar exposure); (b) stream running water stations; (c) stations connected to source waters but far downstream from source points; and (d) isolated ponds disconnected from ground sources or streams with a longer isolation and solar exposure history. Very high values of 0.25 μE m^-2 s^-1, 72 W m^-2 and 12 W m^-2 were measured for PAR, UVA, and UVB incident solar radiation, respectively. The environmental factors measured formed two groups of stations reflected by principal component analyses (near to groundwater sources and isolated systems) where isolated ponds had the highest BSP and microbial abundance (35 microalgae taxa, picoeukaryotes, nanoflagellates, and bacteria) plus higher salinities and PO₄^3- concentrations. BSP short-term response (4 h) to solar radiation was measured by ³H-leucine incorporation under four different solar conditions: full sun, no UVB, PAR, and dark. Microbial communities established in waters with the longest surface exposure (e.g., isolated ponds) had the lowest BSP response to solar radiation treatments, and thus were likely best adapted to solar radiation exposure contrary to ground source waters. These results support our light history (solar exposure) hypothesis where the more isolated the community is from ground water sources, the better adapted it is to solar radiation. We suggest that factors other than solar radiation (e.g., salinity, PO₄^3-, NO₃^-) are also important in determining microbial productivity in heterogeneous environments such as the Salar de Huasco.

The cryptic origins of evolutionary novelty: 1000-fold faster trophic diversification rates without increased ecological opportunity or hybrid swarm

Ecological opportunity is frequently proposed as the sole ingredient for adaptive radiation into novel niches. An additional trigger may be genome-wide hybridization resulting from "hybrid swarm." However, these hypotheses have been difficult to test due to the rarity of comparable control environments lacking adaptive radiations. Here I exploit such a pattern in microendemic radiations of Caribbean pupfishes. I show that a sympatric three species radiation on San Salvador Island, Bahamas diversified 1445 times faster than neighboring islands in jaw length due to the evolution of a novel scale-eating adaptive zone from a generalist ancestral niche. I then sampled 22 generalist populations on seven neighboring islands and measured morphological diversity, stomach content diversity, dietary isotopic diversity, genetic diversity, lake/island areas, macroalgae richness, and Caribbean-wide patterns of gene flow. None of these standard metrics of ecological opportunity or gene flow were associated with adaptive radiation, except for slight increases in macroalgae richness. Thus, exceptional trophic diversification is highly localized despite myriad generalist populations in comparable environmental and genetic backgrounds. This study provides a strong counterexample to the ecological and hybrid swarm theories of adaptive radiation and suggests that diversification of novel specialists on a sparse fitness landscape is constrained by more than ecological opportunity and gene flow.

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

The Andean Altiplano has served as a complex setting throughout its history, driving dynamic processes of diversification in several taxa. We investigated phylogeographic processes in the Telmatobius marmoratus species complex occurring in this region by studying the geographic patterns of genetic variability, genealogies, and historical migration, using the cytochrome b (cyt-b) gene as a marker. DNA sequences from Telmatobius gigas and Telmatobius culeus, Bolivian species with an uncertain taxonomic status, were also included. Additionally, we evaluated the phylogenetic diversity (PD) represented within Chilean protected areas and the complementary contribution from unprotected populations. Phylogenetic reconstructions from 148 cyt-b sequences revealed 4 main clades, one of which corresponded to T. culeus. T. gigas was part of T. marmoratus clade indicating paraphyletic relationships. Haplotypes from Chilean and Bolivian sites were not reciprocally monophyletic. Geographic distribution of lineages, spatial Bayesian analysis, and migration patterns indicated that T. marmoratus displays a weaker geographic structure than expected based on habitat distribution and physiological requirements. Demographic and statistical phylogeography analyses pointed out to a scenario of recent population expansion and high connectivity events of a more recent age than the post Last Glacial Maximum, probably associated to more humid events in Altiplano. PD of T. marmoratus populations within protected areas represents 55.6% of the total estimated PD. The unprotected populations that would contribute the most to PD are Caquena and Quebe (21%). Recent evolutionary processes and paleoclimatic changes, potentially driving shifts in habitat connectivity levels and population sizes, could explain the phylogeographic patterns recovered herein.

The complete mitochondrial genome of the killifish Orestias sp. (Cyprinodontiformes, Cyprinodontidae) from the high Andean range

The killifish genus Orestias is endemic to freshwater ecosystems in the High Andes of Peru, Bolivia and Chile. Phylogenetic and phylogeographic studies with partial mitochondrial genes have failed to resolve relationship among species, needing more comprehensive approaches. In this study, we described the complete mitochondrial genome of Orestias sp. from Lirima, northern Chile, with the aim to provide useful data for phylogenetic purposes and species delimitation. The mitochondrial genome was assembled with 2.6 million of reads obtained through an Ion Torrent (chip 318) sequencer. The circular sequence of 16,617 bp showed the following nucleotide composition: A, 26.7%, C, 27.1%, G, 17.0%, and T, 29.2%. Gene composition and structure were similar to other fish sequences available, and comprised 13 protein-coding genes, 12S and 16S rRNA, 22 tRNA genes, and a control region.

Comparative phylogeography of co-distributed Phrygilus species (Aves, Thraupidae) from the Central Andes

The Neotropical ecoregion has been an important place of avian diversification where dispersal and allopatric events coupled with periods of active orogeny and climate change (Late Pliocene-Pleistocene) have shaped the biogeography of the region. In the Neotropics, avian population structure has been sculpted not only by geographical barriers, but also by non-allopatric factors such as natural selection and local adaptation. We analyzed the genetic variation of six co-distributed Phrygilus species from the Central Andes, based on mitochondrial and nuclear markers in conjunction with morphological differentiation. We examined if Phrygilus species share patterns of population structure and historical demography, and reviewed the intraspecific taxonomy in part of their geographic range. Our results showed different phylogeographic patterns between species, even among those belonging to the same phylogenetic clade. P. alaudinus, P. atriceps, and P. unicolor showed genetic differentiation mediated by allopatric mechanisms in response to specific geographic barriers; P. gayi showed sympatric lineages in northern Chile, while P. plebejus and P. fruticeti showed a single genetic group. We found no relationship between geographic range size and genetic structure. Additionally, a signature of expansion was found in three species related to the expansion of paleolakes in the Altiplano region and the drying phase of the Atacama Desert. Morphological analysis showed congruence with molecular data and intraspecific taxonomy in most species. While we detected genetic and phenotypic patterns that could be related to natural selection and local adaptation, our results indicate that allopatric events acted as a major factor in the population differentiation of Phrygilus species.