Genomic data and multi-species demographic modelling uncover past hybridization between currently allopatric freshwater species.

Evidence for ancient interspecific gene flow through hybridization has been reported in many animal and plant taxa based on genetic markers. The study of genomic patterns of closely related species with allopatric distributions allows the assessment of the relative importance of vicariant isolating events and past gene flow. Here, we investigated the role of gene flow in the evolutionary history of four closely related freshwater fish species with currently allopatric distributions in western Iberian rivers-Squalius carolitertii, S. pyrenaicus, S. torgalensis and S. aradensis-using a population genomics dataset of 23,562 SNPs from 48 individuals, obtained through genotyping by sequencing (GBS). We uncovered a species tree with two well-differentiated clades: (i) S. carolitertii and S. pyrenaicus; and (ii) S. torgalensis and S. aradensis. By using D-statistics and demographic modelling based on the site frequency spectrum, comparing alternative demographic scenarios of hybrid origin, secondary contact and isolation, we found that the S. pyrenaicus North lineage is likely the result of an ancient hybridization event between S. carolitertii (contributing ~84%) and S. pyrenaicus South lineage (contributing ~16%), consistent with a hybrid speciation scenario. Furthermore, in the hybrid lineage, we identify outlier loci potentially affected by selection favouring genes from each parental lineage at different genomic regions. Our results suggest that ancient hybridization can affect speciation and that freshwater fish species currently in allopatry are useful to study these processes.

Adaptation and convergence in circadian-related genes in Iberian freshwater fish.

BACKGROUND: The circadian clock is a biological timing system that improves the ability of organisms to deal with environmental fluctuations. At the molecular level it consists of a network of transcription-translation feedback loops, involving genes that activate (bmal and clock - positive loop) and repress expression (cryptochrome (cry) and period (per) - negative loop). This is regulated by daily alternations of light but can also be affected by temperature. Fish, as ectothermic, depend on the environmental temperature and thus are good models to study its integration within the circadian system. Here, we studied the molecular evolution of circadian genes in four Squalius freshwater fish species, distributed across Western Iberian rivers affected by two climatic types with different environmental conditions (e.g., light and temperature). S. carolitertii and S. pyrenaicus inhabit the colder northern region under Atlantic climate type, while S. torgalensis, S. aradensis and some populations of S. pyrenaicus inhabit the warmer southern region affected by summer droughts, under Mediterranean climate type. RESULTS: We identified 16 circadian-core genes in the Squalius species using a comparative transcriptomics approach. We detected evidence of positive selection in 12 of these genes using methods based on dN/dS. Positive selection was mainly found in cry and per genes of the negative loop, with 55 putatively adaptive substitutions, 16 located on protein domains. Evidence for positive selection is predominant in southern populations affected by the Mediterranean climate type. By predicting protein features we found that changes at sites under positive selection can impact protein thermostability by changing their aliphatic index and isoelectric point. Additionally, in nine genes, the phylogenetic clustering of species that belong to different clades but inhabit southern basins with similar environmental conditions indicated evolutionary convergence. We found evidence for increased nonsynonymous substitution rate in convergent lineages, likely due to positive selection at 27 sites, mostly in cry genes. CONCLUSIONS: Our results support that temperature may be a selective pressure driving the evolution of genes involved in the circadian system. By integrating sequence-based functional protein prediction with dN/dS-based methods to detect selection we uncovered adaptive convergence in the southern populations, probably related to their similar thermal conditions.

Different ecophysiological responses of freshwater fish to warming and acidification.

Future climate change scenarios predict threatening outcomes to biodiversity. Available empirical data concerning biological response of freshwater fish to climate change remains scarce. In this study, we investigated the physiological and biochemical responses of two Iberian freshwater fish species (Squalius carolitertii and the endangered S. torgalensis), inhabiting different climatic conditions, to projected future scenarios of warming (+3°C) and acidification (ΔpH=-0.4). Herein, metabolic enzyme activities of glycolytic (citrate synthase - CS, lactate dehydrogenase - LDH) and antioxidant (glutathione S-transferase, catalase and superoxide dismutase) pathways, as well as the heat shock response (HSR) and lipid peroxidation were determined. Our results show that, under current water pH, warming causes differential interspecific changes on LDH activity, increasing and decreasing its activity in S. carolitertii and in S. torgalensis, respectively. Furthermore, the synergistic effect of warming and acidification caused an increase in LDH activity of S. torgalensis, comparing with the warming condition. As for CS activity, acidification significantly decreased its activity in S. carolitertii whereas in S. torgalensis no significant effect was observed. These results suggest that S. carolitertii is more vulnerable to climate change, possibly as the result of its evolutionary acclimatization to milder climatic condition, while S. torgalensis evolved in the warmer Mediterranean climate. However, significant changes in HSR were observed under the combined warming and acidification (S. carolitertii) or under acidification (S. torgalensis). Our results underlie the importance of conducting experimental studies and address species endpoint responses under projected climate change scenarios to improve conservation strategies, and to safeguard endangered freshwater fish.

Protein analysis and gene expression indicate differential vulnerability of Iberian fish species under a climate change scenario.

Current knowledge on the biological responses of freshwater fish under projected scenarios of climate change remains limited. Here, we examine differences in the protein configuration of two endemic Iberian freshwater fish species, Squalius carolitertii and the critically endangered S. torgalensis that inhabit in the Atlantic-type northern and in the Mediterranean-type southwestern regions, respectively. We performed protein structure modeling of fourteen genes linked to protein folding, energy metabolism, circadian rhythms and immune responses. Structural differences in proteins between the two species were found for HSC70, FKBP52, HIF1α and GPB1. For S. torgalensis, besides structural differences, we found higher thermostability for two proteins (HSP90 and GBP1), which can be advantageous in a warmer environment. Additionally, we investigated how these species might respond to projected scenarios of 3° climate change warming, acidification (ΔpH = -0.4), and their combined effects. Significant changes in gene expression were observed in response to all treatments, particularly under the combined warming and acidification. While S. carolitertii presented changes in gene expression for multiple proteins related to folding (hsp90aa1, hsc70, fkbp4 and stip1), only one such gene was altered in S. torgalensis (stip1). However, S. torgalensis showed a greater capacity for energy production under both the acidification and combined scenarios by increasing cs gene expression and maintaining ldha gene expression in muscle. Overall, these findings suggest that S. torgalensis is better prepared to cope with projected climate change. Worryingly, under the simulated scenarios, disturbances to circadian rhythm and immune system genes (cry1aa, per1a and gbp1) raise concerns for the persistence of both species, highlighting the need to consider multi-stressor effects when evaluating climate change impacts upon fish. This work also highlights that assessments of the potential of endangered freshwater species to cope with environmental change are crucial to help decision-makers adopt future conservation strategies.

Non-canonical expression patterns and evolutionary rates of sex-biased genes in a seasonal fish.

Sex determination is a highly variable process that utilizes many different mechanisms to initiate the cascade of differentiation processes. The molecular pathways controlling sexual development are less conserved than previously assumed, and appear to require active maintenance in some species; indeed, the developmental decision of gonad phenotype in gonochoristic species is not fixed at an early developmental stage. Much of the knowledge about sex determination mechanisms was derived from research on gonochoristic, non-seasonal breeders. In this study, the transcriptome of resting adult gonads of a seasonal breeder, the endangered Iberian cyprinid fish Squalius pyrenaicus, was analyzed to assess the expression patterns and evolutionary rates of sex-biased genes that could be involved in maintenance of gonad identity as well as in sex determination. Remarkably, some crucial female genes-such as aromatase cyp19a1a, estrogen receptor esr1a, and foxl2-were expressed more abundantly in S. pyrenaicus testis than in ovaries. Moreover, contrary to the higher evolutionary rate changes observed in male-biased genes, higher dN /dS ratios were observed for female-biased genes than for male-biased genes in S. pyrenaicus. These results help unravel the impact of seasonality in sex determination mechanisms and the evolution of genes, and highlight the need to study fish at different gonadal maturation states to understand the function of sex-biased genes. Mol. Reprod. Dev. 83: 1102-1115, 2016. © 2016 Wiley Periodicals, Inc.

Gene copy silencing and DNA methylation in natural and artificially produced allopolyploid fish.

Allelic silencing is an important mechanism for coping with gene dosage changes in polyploid organisms that is well known in allopolyploid plants. Only recently, it was shown in the allotriploid fish Squalius alburnoides that this process also occurs in vertebrates. However, it is still unknown whether this silencing mechanism is common to other allopolyploid fish, and which mechanisms might be responsible for allelic silencing. We addressed these questions in a comparative study between Squalius alburnoides and another allopolyploid complex, the Amazon molly (Poecilia formosa). We examined the allelic expression patterns for three target genes in four somatic tissues of natural allo-anorthoploids and laboratory-produced tri-genomic hybrids of S. alburnoides and P. formosa. Also, for both complexes, we evaluated the correlation between total DNA methylation level and the ploidy status and genomic composition of the individuals. We found that allelic silencing also occurs in other allopolyploid organisms besides the single one that was previously known. We found and discuss disparities within and between the two considered complexes concerning the pattern of allele-specific expression and DNA methylation levels. Disparities might be due to intrinsic characteristics of each genome involved in the hybridization process. Our findings also support the idea that long-term evolutionary processes have an effect on the allele expression patterns and possibly also on DNA methylation levels.

Genomic Resources Notes Accepted 1 June 2015-31 July 2015.

This article documents the public availability of (i) microbiomes in diet and gut of larvae from the dipteran Dilophus febrilis using massive parallel sequencing, (ii) SNP and SSR discovery and characterization in the transcriptome of the Atlantic mackerel (Scomber scombrus, L) and (iii) assembled transcriptome for an endangered, endemic Iberian cyprinid fish (Squalius pyrenaicus).

Hybrid speciation in a marine mammal: the clymene dolphin (Stenella clymene).

Natural hybridization may result in the exchange of genetic material between divergent lineages and even the formation of new taxa. Many of the Neo-Darwinian architects argued that, particularly for animal clades, natural hybridization was maladaptive. Recent evidence, however, has falsified this hypothesis, instead indicating that this process may lead to increased biodiversity through the formation of new species. Although such cases of hybrid speciation have been described in plants, fish and insects, they are considered exceptionally rare in mammals. Here we present evidence for a marine mammal, Stenella clymene, arising through natural hybridization. We found phylogenetic discordance between mitochondrial and nuclear markers, which, coupled with a pattern of transgressive segregation seen in the morphometric variation of some characters, support a case of hybrid speciation. S. clymene is currently genetically differentiated from its putative parental species, Stenella coerueloalba and Stenella longisrostris, although low levels of introgressive hybridization may be occurring. Although non-reticulate forms of evolution, such as incomplete lineage sorting, could explain our genetic results, we consider that the genetic and morphological evidence taken together argue more convincingly towards a case of hybrid speciation. We anticipate that our study will bring attention to this important aspect of reticulate evolution in non-model mammal species. The study of speciation through hybridization is an excellent opportunity to understand the mechanisms leading to speciation in the context of gene flow.

Different levels of hsp70 and hsc70 mRNA expression in Iberian fish exposed to distinct river conditions.

Comprehension of the mechanisms by which ectotherms, such as fish, respond to thermal stress is paramount for understanding the threats that environmental changes may pose to wild populations. Heat shock proteins are molecular chaperones with an important role in several stress conditions such as high temperatures. In the Iberian Peninsula, particularly in Portugal, freshwater fish of the genus Squalius are subject to daily and seasonal temperature variations. To examine the extent to which different thermal regimes influence the expression patterns of hsp70 and hsc70 transcripts we exposed two species of Squalius (S. torgalensis and S. carolitertii) to different temperatures (20, 25, 30 and 35 °C). At 35 °C, there was a significant increase in the expression of hsp70 and hsc70 in the southern species, S. torgalensis, while the northern species, S. carolitertii, showed no increase in the expression of these genes; however, some individuals of the latter species died when exposed to 35 °C. These results suggest that S. torgalensis may cope better with harsher temperatures that are characteristic of this species natural environment; S. carolitertii, on the other hand, may be unable to deal with the extreme temperatures faced by the southern species.

Ploidy mosaicism and allele-specific gene expression differences in the allopolyploid Squalius alburnoides.

BACKGROUND: Squalius alburnoides is an Iberian cyprinid fish resulting from an interspecific hybridisation between Squalius pyrenaicus females (P genome) and males of an unknown Anaecypris hispanica-like species (A genome). S. alburnoides is an allopolyploid hybridogenetic complex, which makes it a likely candidate for ploidy mosaicism occurrence, and is also an interesting model to address questions about gene expression regulation and genomic interactions. Indeed, it was previously suggested that in S. alburnoides triploids (PAA composition) silencing of one of the three alleles (mainly of the P allele) occurs. However, not a whole haplome is inactivated but a more or less random inactivation of alleles varying between individuals and even between organs of the same fish was seen.In this work we intended to correlate expression differences between individuals and/or between organs to the occurrence of mosaicism, evaluating if mosaics could explain previous observations and its impact on the assessment of gene expression patterns. RESULTS: To achieve our goal, we developed flow cytometry and cell sorting protocols for this system generating more homogenous cellular and transcriptional samples. With this set-up we detected 10% ploidy mosaicism within the S. alburnoides complex, and determined the allelic expression profiles of ubiquitously expressed genes (rpl8; gapdh and ß-actin) in cells from liver and kidney of mosaic and non-mosaic individuals coming from different rivers over a wide geographic range. CONCLUSIONS: Ploidy mosaicism occurs sporadically within the S. alburnoides complex, but in a frequency significantly higher than reported for other organisms. Moreover, we could exclude the influence of this phenomenon on the detection of variable allelic expression profiles of ubiquitously expressed genes (rpl8; gapdh and ß-actin) in cells from liver and kidney of triploid individuals. Finally, we determined that the expression patterns previously detected only in a narrow geographic range is not a local restricted phenomenon but is pervasive in rivers where S. pyrenaicus is sympatric with S. alburnoides.We discuss mechanisms that could lead to the formation of mosaic S. alburnoides and hypothesise about a relaxation of the mechanisms that impose a tight control over mitosis and ploidy control in mixoploids.

Cranial shape differentiation in three closely related delphinid cetacean species: insights into evolutionary history.

The present study investigates the pattern of differentiation of cranial shape in three closely related delphinid cetacean species of the complex Delphinus-Stenella-Tursiops: Delphinus delphis, Stenella coeruleoalba and Tursiops truncatus. Dorsal and ventral aspects of the cranium were analysed using landmark-based geometric morphometric methods. While there was no evidence of sexual dimorphism for shape or size, multivariate statistical analyses showed that there were interspecific differences in skull morphology. Skull shape differences between the three studied species were related with cranial width and differences in the length of the rostrum relative to the cranial portion of the skull. D. delphis and S. coeruleoalba showed high cranial shape similarity, which is indicative of their evolutionary proximity when compared with T. truncatus. Phenetic clusters based on cranial shape similarities were found to be concordant with the molecular phylogenetic clades obtained from mitochondrial DNA genes. Geometric morphometric methods can thus be an exceptionally useful tool for the study of differentiation of delphinid cetacean species and therefore provide some insights into their evolutionary history.

Dosage compensation by gene-copy silencing in a triploid hybrid fish.

In mammals, the increase in gene dosage, in the form of polyploidy or involving chromosomal fragments, has deleterious effects [1]. Regulation of appropriate gene product amounts has to be warranted by complex dosage-compensation mechanisms. Lower vertebrates, on the other hand, cope very well with ploidy increase [2-4], implying either effective compensation or a lack of necessity for such mechanisms. Unfortunately, nothing is known about the genetic and molecular mechanisms underlying this phenomenon. For an experimental approach, we have studied gene expression in the allotriploid form of Squalius alburnoides. In these organisms, different genomes are joined through hybridization; thus, sequence differences can be used to follow expression of different alleles [5, 6]. We found that a compensation mechanism exists, reducing transcript levels to the diploid state. Our data suggest a silencing of one of the three alleles. Unexpectedly, it is not a whole haplome that is inactivated. The allelic expression patterns differ between genes and between different tissues for one and the same gene. Our data provide the first evidence of a regulation mechanism involving gene-copy silencing in a triploid vertebrate.

Speciation towards tetraploidization after intermediate processes of non-sexual reproduction.

Polyploidy, hybridization and variation in mating systems are central issues for a deeper understanding of animal evolution. The Iberian species Squalius alburnoides represents an example combining all three phenomena. Previous studies showed that S. alburnoides populations are mainly composed of triploid and diploid hybrid forms (mainly females), and that the tetraploid forms are rare or absent. Both populations from the Douro drainage reveal a distinct scenario: tetraploid individuals represent 85.6-97.5% of the population, with no sex ratio bias observed. Based on the flow cytometry measurements of blood and spermatozoa cells, microsatellite loci and experimental crosses, we describe here, for the first time, two symmetric allotetraploid populations (CCAA) that resumed normal meiosis after undergoing intermediate processes of non-sexual reproduction to give rise to a new sexually reproducing polyploid species. Prezygotic (habitat selection and assortative mating) and postzygotic mechanisms (nonviable embryos) are responsible for the reproductive isolation from other forms of the S. alburnoides complex (e.g. CA, CAA). This example illustrates how hybrid polyploid complexes may lead to speciation.

Modes of reproduction of the hybridogenetic fish Squalius alburnoides in the Tejo and Guadiana rivers: an approach with microsatellites.

The Squalius alburnoides complex was produced by hybridization between female S. pyrenaicus (PP genome) and an hypothetical paternal ancestor related with Anaecypris hispanica (AA genome). This study examined a diversity of mating types and found that there is the potential for considerable gene exchange among diploid, triploid and tetraploid hybrids. Using microsatellites, genomes were attributed to Squalius pyrenaicus (P) or reconstituted "nuclear non-hybrid"S. alburnoides (A), and subsequently confirmed in hybrids. Recombination of AA genomes in the "nuclear non-hybrid males" and recombination of the homogametic genomes (AA or PP) after exclusion of the heterogametic genome in triploid females (PAA) were observed by analysing parents and progeny of breeding experiments. Reproduction of tetraploids, generating a symmetric tetraploid genotype (PPAA) in the progeny, suggests a process that could potentially lead to the formation of a new bisexual species. Present results also support: (i) previously hypothesized pathways, in which PPA S. alburnoides females exclude the A genome, exhibit meiotic recombination between the P genomes and generate haploid eggs; (ii) reconstitution of the diploid maternal ancestor genome (PP) as well as of the unknown paternal ancestor (AA); (iii) the occurrence of the same genomic reproductive mechanisms when Anaecypris hispanica is involved; and (iv) the existence of an A. hispanica-like ancestor as the paternal ancestor of S. alburnoides.