Population transcriptomics reveals weak parallel genetic basis in repeated marine and freshwater divergence in nine‐spined sticklebacks

The role of plastic and evolved DNA methylation in parallel adaptation of threespine stickleback (Gasterosteus aculeatus)

Repeated phenotypic patterns among populations undergoing parallel evolution in similar environments provide support for the deterministic role of natural selection. Epigenetic modifications can mediate plastic and evolved phenotypic responses to environmental change and might make important contributions to parallel adaptation. While many studies have explored the genetic basis of repeated phenotypic divergence, the role of epigenetic processes during parallel adaptation remains unclear. The parallel evolution of freshwater ecotypes of threespine stickleback fish (Gasterosteus aculeatus) following colonization of thousands of lakes and streams from the ocean is a classic example of parallel phenotypic and genotypic adaptation. To investigate epigenetic modifications during parallel adaptation of threespine stickleback, we reanalysed three independent data sets that investigated DNA methylation variation between marine and freshwater ecotypes. Although we found widespread methylation differentiation between ecotypes, there was no significant tendency for CpG sites associated with repeated methylation differentiation across studies to be parallel versus nonparallel. To next investigate the role of plastic versus evolved changes in methylation during freshwater adaptation, we explored if CpG sites exhibiting methylation plasticity during salinity change were more likely to also show evolutionary divergence in methylation between ecotypes. The directions of divergence between ecotypes were generally in the opposite direction to those observed for plasticity when ecotypes were challenged with non-native salinity conditions, suggesting that most plastic responses are likely to be maladaptive during colonization of new environments. Finally, we found a greater number of CpG sites showing evolved changes when ancestral marine ecotypes are acclimated to freshwater environments, whereas plastic changes predominate when derived freshwater ecotypes transition back to their ancestral marine environments. These findings provide evidence for an epigenetic contribution to parallel adaptation and demonstrate the contrasting roles of plastic and evolved methylation differences during adaptation to new environments.

Complex population history affects admixture analyses in nine-spined sticklebacks

Introgressive hybridization is an important process in evolution but challenging to identify, undermining the efforts to understand its role and significance. On the contrary, many analytical methods assume direct descent from a single common ancestor, and admixture among populations can violate their assumptions and lead to seriously biased results. A detailed analysis of 888 whole-genome sequences of nine-spined sticklebacks (Pungitius pungitius) revealed a complex pattern of population ancestry involving multiple waves of gene flow and introgression across northern Europe. The two recognized lineages were found to have drastically different histories, and their secondary contact zone was wider than anticipated, displaying a smooth gradient of foreign ancestry with some curious deviations from the expected pattern. Interestingly, the freshwater isolates provided peeks into the past and helped to understand the intermediate states of evolutionary processes. Our analyses and findings paint a detailed picture of the complex colonization history of northern Europe and provide backdrop against which introgression and its role in evolution can be investigated. However, they also expose the challenges in analyses of admixed populations and demonstrate how hidden admixture and colonization history misleads the estimation of admixture proportions and population split times.

Contrasting population differentiation in two sympatric Triplophysa loaches on the Qinghai–Tibet Plateau

Genetic differentiation in aquatic organisms is usually shaped by drainage connectivity. Sympatric aquatic species are thus expected to show similar population differentiation patterns and similar genetic responses to their habitats. Water bodies on the Qinghai–Tibet Plateau (QTP) have recently experienced dramatic physicochemical changes, threatening the biodiversity of aquatic organisms on the “roof of the world.” To uncover ecological genetics in Tibetan loaches (Triplophysa)—the largest component of the QTP ichthyofauna—we characterized population differentiation patterns and adaptive mechanisms to salinity change in two sympatric and phylogenetically closely related Tibetan loaches, T. stewarti and T. stenura, by integrating population genomic, transcriptomic, and electron probe microanalysis approaches. Based on millions of genome-wide SNPs, the two Tibetan loach species show contrasting population differentiation patterns, with highly geographically structured and clear genetic differentiation among T. stewarti populations, whereas there is no such observation in T. stenura, which is also supported by otolith microchemistry mapping. While limited genetic signals of parallel adaption to salinity changes between the two species are found from either genetic or gene expression variation perspective, a catalog of genes involved in ion transport, energy metabolism, structural reorganization, immune response, detoxification, and signal transduction is identified to be related to adaptation to salinity change in Triplophysa loaches. Together, our findings broaden our understanding of the population characteristics and adaptive mechanisms in sympatric Tibetan loach species and would contribute to biodiversity conservation and management of aquatic organisms on the QTP.

Swimming behaviour in two ecologically similar three-spined (Gasterosteus aculeatus L.) and nine-spined sticklebacks (Pungitius pungitius L.): a comparative approach for modelling the toxicity of metal mixtures

Sticklebacks (Gasterosteiformes) are increasingly used in ecological and evolutionary research and have become well established as role model species for biologists. However, ecotoxicology studies concerning behavioural effects in sticklebacks regarding stress responses, mainly induced by chemical mixtures, have hardly been addressed. For this purpose, we investigated the swimming behaviour (including mortality rate based on 96-h LC₅₀ values) of two ecologically similar three-spined (Gasterosteus aculeatus) and nine-spined sticklebacks (Pungitius pungitius) to short-term (up to 24 h) metal mixture (MIX) exposure. We evaluated the relevance and efficacy of behavioural responses of test species in the early toxicity assessment of chemical mixtures. Fish exposed to six (Zn, Pb, Cd, Cu, Ni, and Cr) metals in the mixture were either singled out by the Water Framework Directive as priority or as relevant substances in surface water, which was prepared according to the environmental quality standards (EQSs) of these metals set for inland waters in the European Union (EU) (Directive 2013/39/EU). The performed behavioural analysis showed the main effect on the interaction between time, species, and treatment variables. Although both species exposed to MIX revealed a decreasing tendency in swimming activity, these species' responsiveness to MIX was somewhat different. Substantial changes in the activity of G. aculeatus were established after a 3-h exposure to MIX solutions, which was 1.43-fold lower, while in the case of P. pungitius, 1.96-fold higher than established 96-h LC₅₀ values for each species. This study demonstrated species-specific differences in response sensitivity to metal-based water pollution, indicating behavioural insensitivity of P. pungitius as model species for aquatic biomonitoring and environmental risk assessments.

Comparing genome scans among species of the stickleback order reveals three different patterns of genetic diversity

Comparing genome scans among species is a powerful approach for investigating the patterns left by evolutionary processes. In particular, this offers a way to detect candidate genes that drive convergent evolution. We compared genome scan results to investigate if patterns of genetic diversity and divergence are shared among divergent species within the stickleback order (Gasterosteiformes): the threespine stickleback (Gasterosteus aculeatus), ninespine stickleback (Pungitius pungitus), and tubesnout (Aulorhynchus flavidus). Populations were sampled from the southern and northern edges of each species' range, to identify patterns associated with latitudinal changes in genetic diversity. Weak correlations in genetic diversity (F ST and expected heterozygosity) and three different patterns in the genomic landscape were found among these species. Additionally, no candidate genes for convergent evolution were detected. This is a counterexample to the growing number of studies that have shown overlapping genetic patterns, demonstrating that genome scan comparisons can be noisy due to the effects of several interacting evolutionary forces.

Integrated miRNA and transcriptome profiling to explore the molecular determinism of convergent adaptation to corn in two lepidopteran pests of agriculture

Background

The degree to which adaptation to same environment is determined by similar molecular mechanisms, is a topic of broad interest in evolutionary biology, as an indicator of evolutionary predictability. We wished to address if adaptation to the same host plant in phytophagous insects involved related gene expression patterns. We compared sRNA-Seq and RNA-Seq data between two pairs of taxa of Ostrinia and Spodoptera frugiperda sharing maize as host-plant. For the latter, we had previously carried out a reciprocal transplant experiment by feeding of the larvae of the Corn strain (Sf-C) and the Rice strain (Sf-R) on corn versus rice and characterized the mRNA and miRNA responses.

Results

First, we predicted the genes encoding miRNA in Ostrinia nubilalis (On) and O. scapulalis (Os). Respectively 67 and 65 known miRNA genes, as well as 196 and 190 novel ones were predicted with Os genome using sncRNAs extracted from whole larvae feeding on corn or mugwort. In On, a read counts analysis showed that 37 (55.22%) known miRNAs and 19 (9.84%) novel miRNAs were differentially expressed (DE) on mugwort compared to corn (in Os, 25 known miRs (38.46%) and 8 novel ones (4.34%)). Between species on corn, 8 (12.5%) known miRNAs and 8 (6.83%) novel ones were DE while only one novel miRNA showed expression variation between species on mugwort. Gene target prediction led to the identification of 2953 unique target genes in On and 2719 in Os, among which 11.6% (344) were DE when comparing species on corn. 1.8% (54) of On miR targets showed expression variation upon a change of host-plant.

We found molecular changes matching convergent phenotype, i.e., a set of nine miRNAs that are regulated either according to the host-plant both in On and Sf-C or between them on the same plant, corn. Among DE miR target genes between taxa, 13.7% shared exactly the same annotation between the two pairs of taxa and had function related to insect host-plant interaction.

Conclusion

There is some similarity in underlying genetic mechanisms of convergent evolution of two distant Lepidopteran species having adopted corn in their host range, highlighting possible adaptation genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07905-7.

Invasion history shapes host transcriptomic response to a body-snatching parasite

By shuffling biogeographical distributions, biological invasions can both disrupt long-standing associations between hosts and parasites and establish new ones. This creates natural experiments with which to study the ecology and evolution of host-parasite interactions. In estuaries of the Gulf of Mexico, the white-fingered mud crab (Rhithropanopeus harrisii) is infected by a native parasitic barnacle, Loxothylacus panopaei (Rhizocephala), which manipulates host physiology and behaviour. In the 1960s, L. panopaei was introduced to the Chesapeake Bay and has since expanded along the southeastern Atlantic coast, while host populations in the northeast have so far been spared. We use this system to test the host's transcriptomic response to parasitic infection and investigate how this response varies with the parasite's invasion history, comparing populations representing (i) long-term sympatry between host and parasite, (ii) new associations where the parasite has invaded during the last 60 years and (iii) naïve hosts without prior exposure. A comparison of parasitized and control crabs revealed a core response, with widespread downregulation of transcripts involved in immunity and moulting. The transcriptional response differed between hosts from the parasite's native range and where it is absent, consistent with previous observations of increased susceptibility in populations lacking exposure to the parasite. Crabs from the parasite's introduced range, where prevalence is highest, displayed the most dissimilar response, possibly reflecting immune priming. These results provide molecular evidence for parasitic manipulation of host phenotype and the role of gene regulation in mediating host-parasite interactions.

Adaptive Divergence under Gene Flow along an Environmental Gradient in Two Coexisting Stickleback Species

There is a general and solid theoretical framework to explain how the interplay between natural selection and gene flow affects local adaptation. Yet, to what extent coexisting closely related species evolve collectively or show distinctive evolutionary responses remains a fundamental question. To address this, we studied the population genetic structure and morphological differentiation of sympatric three-spined and nine-spined stickleback. We conducted genotyping-by-sequencing and morphological trait characterisation using 24 individuals of each species from four lowland brackish water (LBW), four lowland freshwater (LFW) and three upland freshwater (UFW) sites in Belgium and the Netherlands. This combination of sites allowed us to contrast populations from isolated but environmentally similar locations (LFW vs. UFW), isolated but environmentally heterogeneous locations (LBW vs. UFW), and well-connected but environmentally heterogenous locations (LBW vs. LFW). Overall, both species showed comparable levels of genetic diversity and neutral genetic differentiation. However, for all three spatial scales, signatures of morphological and genomic adaptive divergence were substantially stronger among populations of the three-spined stickleback than among populations of the nine-spined stickleback. Furthermore, most outlier SNPs in the two species were associated with local freshwater sites. The few outlier SNPs that were associated with the split between brackish water and freshwater populations were located on one linkage group in three-spined stickleback and two linkage groups in nine-spined stickleback. We conclude that while both species show congruent evolutionary and genomic patterns of divergent selection, both species differ in the magnitude of their response to selection regardless of the geographical and environmental context.

The population genomics of repeated freshwater colonizations by Gulf pipefish

How organisms adapt to the novel challenges imposed by the colonization of a new habitat has long been a central question in evolutionary biology. When multiple populations of the same species independently adapt to similar environmental challenges, the question becomes whether the populations have arrived at their adaptations through the same genetic mechanisms. In recent years, genetic techniques have been used to tackle these questions by investigating the genome-level changes underlying local adaptation. Here, we present a genomic analysis of colonization of freshwater habitats by a primarily marine fish, the Gulf pipefish (Syngnathus scovelli). We sample pipefish from four geographically distinct freshwater locations and use double-digest restriction site associated DNA sequencing to compare them to 12 previously studied saltwater populations. The two most geographically distant and isolated freshwater populations are the most genetically distinct, although demographic analysis suggests that these populations are experiencing ongoing migration with their saltwater neighbours. Additionally, outlier regions were found genome-wide, showing parallelism across ecotype pairs. We conclude that these multiple freshwater colonizations involve similar genomic regions, despite the large geographical distances and different underlying mechanisms. These similar patterns are probably facilitated by the interacting effects of intrinsic barriers, gene flow among populations and ecological selection in the Gulf pipefish.

Genetic population structure constrains local adaptation in sticklebacks

Repeated and independent adaptation to specific environmental conditions from standing genetic variation is common. However, if genetic variation is limited, the evolution of similar locally adapted traits may be restricted to genetically different and potentially less optimal solutions or prevented from happening altogether. Using a quantitative trait locus (QTL) mapping approach, we identified the genomic regions responsible for the repeated pelvic reduction (PR) in three crosses between nine-spined stickleback populations expressing full and reduced pelvic structures. In one cross, PR mapped to linkage group 7 (LG7) containing the gene Pitx1, known to control pelvic reduction also in the three-spined stickleback. In the two other crosses, PR was polygenic and attributed to 10 novel QTL, of which 90% were unique to specific crosses. When screening the genomes from 27 different populations for deletions in the Pitx1 regulatory element, these were only found in the population in which PR mapped to LG7, even though the morphological data indicated large-effect QTL for PR in several other populations as well. Consistent with the available theory and simulations parameterized on empirical data, we hypothesize that the observed variability in genetic architecture of PR is due to heterogeneity in the spatial distribution of standing genetic variation caused by >2× stronger population structuring among freshwater populations and >10× stronger genetic isolation by distance in the sea in nine-spined sticklebacks as compared to three-spined sticklebacks.

Genomic basis of the loss of diadromy in Galaxias maculatus: Insights from reciprocal transplant experiments

Diadromy is known for having major effects on the distribution and richness of aquatic species, and so does its loss. The loss of diadromy has led to the diversification of many species, yet research focusing on understanding its molecular basis and consequences are limited. This is particularly true for amphidromous species despite being the most abundant group of diadromous species. Galaxias maculatus, an amphidromous species and one of the most widely distributed fishes in the Southern Hemisphere, exhibits many instances of nonmigratory or resident populations. The existence of naturally replicated resident populations in Patagonia can serve as an ideal system for the study of the mechanisms that lead to the loss of the diadromy and its ecological and evolutionary consequences. Here, we studied two adjacent river systems in which resident populations are genetically differentiated yet derived from the same diadromous population. By combining a reciprocal transplant experiment with genomic data, we showed that the two resident populations followed different evolutionary pathways by exhibiting a differential response in their capacity to survive in salt water. While one resident population was able to survive salt water, the other was not. Genomic analyses provided insights into the genes that distinguished (a) migratory from nonmigratory populations; (b) populations that can vs those that cannot survive a saltwater environment; and (c) between these resident populations. This study demonstrates that the loss of diadromy can be achieved by different pathways and that environmental (selection) and random (genetic drift) forces shape this dynamic evolutionary process.