Genetic structure in remnant populations of an endangered cyprinodontid fish, Orestias ascotanensis, endemic to the Ascotán salt pan of the Altiplano

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.

Coping With Dynamism: Phylogenetics and Phylogeographic Analyses Reveal Cryptic Diversity in Heleobia Snails of Atacama Saltpan, Chile

The species that inhabit systems highly affected by anthropic activities usually exhibit this external influence in their gene pool. In this study, we investigated the genetic patterns of populations of Heleobia atacamensis, a freshwater microgastropod endemic to the Atacama Saltpan, a system historically exposed to environmental changes, and currently subjected to conditions associated with metallic and non-metallic mining and other anthropic activities. Molecular analyses based on nuclear and mitochondrial sequences indicate that the saltpan populations are highly fragmented and that the genetic structure is mainly explained by historical geographic isolation, with little influence of contemporary factors. The microsatellite results suggest a moderate genetic diversity and sharp differentiation mediated by isolation by distance. Additionally, despite the high environmental heterogeneity detected and the marked historical dynamism of the region, our data reveal no signs of demographic instability. The patterns of contemporary gene flow suggest a change in the current genetic structure, based on the geographic proximity and specific environmental conditions for each population. Our results, highlight the role of fragmentation as a modulator of genetic diversity, but also suggest that the historical persistence of isolated populations in naturally dynamic environments could explain the apparent demographic stability detected.

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.

Phylogenetic analysis of the complete mitochondrial genome of Orestias ascotanensis (Cyprinodontiformes, Cyprinodontidae)

The complete mitochondrial genome of this species was first determined in this study, which is 16,617 bp in length, containing 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, a putative control region, and 1 origin of replication on the light-strand. The overall base composition includes C(27.11%), A(26.68%), T(29.15%), G(17.04%) and three degenerate bases are R, R and S. Moreover, the 13 PCGs encode 3800 amino acids in total, 12 of which use the initiation codon ATG except COI that uses GTG. Most of them have TAA as the stop codon, whereas ND5 ends with AGA, four protein-coding genes (ND1, ND2, ND3 and Cytb) ended with TAG, and two protein-coding genes (COII and ND4) ended with the incomplete stop codon represented as a single T. The phylogenetic tree based on the Neighbor Joining method was constructed to provide relationship within Cyprinodontidae, which could be a useful basis for management of this species.

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.

Genetic divergence of isolated populations of the native micromoth Bucculatrix mirnae (Lepidoptera: Bucculatricidae) in the arid environments of Northern Chile

Analysis of maternally inherited genes is especially helpful in population studies of host-specialized insects, as female dispersal is key to find an adequate host plant to ensure larval survival. Bucculatrix mirnae (Lepidoptera: Bucculatricidae) is a little-known Neotropical micromoth native to the arid environments of northern Chile whose hypermetamorphic larvae are miners and skeletonizers on leaves of two species of Baccharis (Asteraceae) shrubs. This micromoth has been detected in three isolated locations embracing a narrow geographic range: two from the coastal valleys of the Atacama Desert near sea level and one from the western slopes of the Andes at about 3000 m elevation. As the dispersal of B. mirnae is mostly restricted to the small adult stage, the altitudinal gradient and desert areas among the three localities could be effective barriers, triggering genetic differentiation among populations. Sequences of the DNA barcode fragment of the cytocrome oxidase subunit I mitochondrial gene were analyzed to assess for the first time the patterns of genetic variation of B. mirnae. Fifteen haplotypes, each exclusive to one locality, were found in the 71 specimens analyzed. Genetic divergence (K2P) between haplotypes of different localities was at least 2.0%. A Bayesian analysis with sequences of congeneric species grouped all the B. mirnae haplotypes in a clade, in which three well-supported locality-specific haplogroups were found. In concordance with this pattern, an analysis of molecular variance showed that the highest genetic variation was found among populations. Furthermore, all the population pairwise comparisons (F_ST) were significant. These results suggest that female migration between isolated populations of B. mirnae is absent. This pattern must be considered in the current scenario of habitat destruction and modification in the arid environments of northern Chile.

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.

Chromosomal organization of four classes of repetitive DNA sequences in killifish Orestias ascotanensis Parenti, 1984 (Cyprinodontiformes, Cyprinodontidae)

Orestias Valenciennes, 1839 is a genus of freshwater fish endemic to the South American Altiplano. Cytogenetic studies of these species have focused on conventional karyotyping. The aim of this study was to use classical and molecular cytogenetic methods to identify the constitutive heterochromatin distribution and chromosome organization of four classes of repetitive DNA sequences (histone H3 DNA, U2 snRNA, 18S rDNA and 5S rDNA) in the chromosomes of O. ascotanensis Parenti, 1984, an endemic species restricted to the Salar de Ascotán in the Chilean Altiplano. All individuals analyzed had a diploid number of 48 chromosomes. C-banding identified constitutive heterochromatin mainly in the pericentromeric region of most chromosomes, especially a GC-rich heterochromatic block of the short arm of pair 3. FISH assay with an 18S probe confirmed the location of the NOR in pair 3 and revealed that the minor rDNA cluster occurs interstitially on the long arm of pair 2. Dual FISH identified a single block of U2 snDNA sequences in the pericentromeric regions of a subtelocentric chromosome pair, while histone H3 sites were observed as small signals scattered in throughout the all chromosomes. This work represents the first effort to document the physical organization of the repetitive fraction of the Orestias genome. These data will improve our understanding of the chromosomal evolution of a genus facing serious conservation problems.

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.

Historical dimensions of population structure in a continuously distributed marine species: The case of the endemic Chilean dolphin

The complementarity of historical and contemporary processes contributes to understanding the genetic structure of continuously distributed marine species with high dispersal capabilities. Cephalorhynchus eutropia, has a continuous coastal distribution with strong genetic differentiation identified by nuclear DNA markers. We explored the historical dimension of this genetic differentiation between northern and southern populations to evaluate phylogeographic structure. Additionally, we conducted mtDNA and microsatellite analyses to detect past and recent demographic changes. The southern population was characterized by lower genetic diversity with a signal of population expansion, likely associated with ice retreat and habitat extension after the Last Glacial Maximum (LGM). In contrast, structure within the northern population was more consistent with stable historical population size. Approximate Bayesian Computation analyses suggested that during the LGM, C. eutropia persisted in the northern area; while the south was colonized by dispersal ~11,000 years ago followed by population expansion. This study shows that Chilean dolphin population structure is consistent with predictions from the Expansion-Contraction biogeographic model, with a poleward post-glacial shift revealed in current genetic structure. The results also confirm the validity of the population units previously identified, demonstrating their historical origin and highlighting the utility of integrating genetic markers with different temporal scale resolutions.

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

This study presents phylogenetic molecular data of the Chilean species of Orestias to propose an allopatric divergence hypothesis and phylogeographic evidence that suggests the relevance of abiotic factors in promoting population divergence in this complex. The results reveal that diversification is still ongoing, e.g. in the Ascotán salt pan, where populations of Orestias ascotanensis restricted to individual freshwater springs exhibit strong genetic differentiation, reflecting putative independent evolutionary units. Diversification of Orestias in the southern Altiplano may be linked to historical vicariant events and contemporary variation in water level; these processes may have affected the populations from the Plio-Pleistocene until the present.