Population transcriptomics reveals the effect of gene flow on the evolution of range limits

One of the most important questions in evolutionary biology is how the spatial distribution of species is limited. Asymmetric gene flow from core populations is suggested to increase the number of poorly adapted immigrants in the populations at the range edge. Genetic load due to migration, i.e., migration load, should prevent adaptation to the local habitat, leading to decreases in distribution range via local extinction or the limiting range expansion. However, few experimental studies have examined the effects of immigration on fitness and natural selection within recipient populations. To investigate the influence of migration load on the evolution of distribution range, we performed field and laboratory observations as well as population transcriptomics for the common river snail, Semisulcospira reiniana. This species meets the conditions that migration from source populations can prevent local adaptation in a sink population because they inhabit the broader range of environments, including middle/upper reaches of a river and estuaries within a single river and they may be more vulnerable to being swept away by water currents due to lowered spontaneous (upward) locomotion activity. We found that river steepness was related to the lower distribution limit of S. reiniana, with a narrower distribution range in the steeper river. Population transcriptomic analysis showed that gene flow was heavily asymmetric from the upstream populations to downstream ones in the steep river, suggesting a greater migration load in the steep river. The number of genes putatively involved in adaptation to the local habitat was lower in the steep river than in the gentle river. Gene expression profiles suggested that individuals achieve better local adaptation in the gentle river. Laboratory experiments suggested that evolutionary differences in salinity tolerance among local populations were only found in the gentle river. Our results consistent with the hypothesis that migration load owing to asymmetric gene flow disturbs local adaptation and restricts the distribution range of river snails.

Linking ecology, morphology, and metabolism: Niche differentiation in sympatric populations of closely related species of the genus Littorina (Neritrema)

Divergence of ecological niches in phylogenetically closely related species indicates the importance of ecology in speciation, especially for sympatric species are considered. Such ecological diversification provides an advantage of alleviating interspecies competition and promotes more efficient exploitation of environmental resources, thus being a basis for ecological speciation. We analyzed a group of closely related species from the subgenus Neritrema (genus Littorina, Caenogastropoda) from the gravel-bouldery shores. In two distant sites at the Barents and Norwegian Sea, we examined the patterns of snail distribution during low tide (quantitative sampling stratified by intertidal level, presence of macrophytes, macrophyte species, and position on them), shell shape and its variability (geometric morphometrics), and metabolic characteristics (metabolomic profiling). The studied species diversified microbiotopes, which imply an important role of ecological specification in the recent evolution of this group. The only exception to this trend was the species pair L. arcana / L. saxatilis, which is specifically discussed. The ecological divergence was accompanied by differences in shell shape and metabolomic characteristics. Significant differences were found between L. obtusata versus L. fabalis and L. saxatilis / L. arcana versus L. compressa both in shell morphology and in metabolomes. L. saxatilis demonstrated a clear variability depending on intertidal level which corresponds to a shift in conditions within the occupied microhabitat. Interestingly, the differences between L. arcana (inhabiting the upper intertidal level) and L. compressa (inhabiting the lower one) were analogous to those between the upper and lower fractions of L. saxatilis. No significant level-dependent changes were found between the upper and lower fractions of L. obtusata, most probably due to habitat amelioration by fucoid macroalgae. All these results are discussed in the contexts of the role of ecology in speciation, ecological niche dynamics and conservatism, and evolutionary history of the Neritrema species.

Multi-omics investigations within the Phylum Mollusca, Class Gastropoda: from ecological application to breakthrough phylogenomic studies

Gastropods are the largest and most diverse class of mollusc and include species that are well studied within the areas of taxonomy, aquaculture, biomineralization, ecology, microbiome and health. Gastropod research has been expanding since the mid-2000s, largely due to large-scale data integration from next-generation sequencing and mass spectrometry in which transcripts, proteins and metabolites can be readily explored systematically. Correspondingly, the huge data added a great deal of complexity for data organization, visualization and interpretation. Here, we reviewed the recent advances involving gastropod omics ('gastropodomics') research from hundreds of publications and online genomics databases. By summarizing the current publicly available data, we present an insight for the design of useful data integrating tools and strategies for comparative omics studies in the future. Additionally, we discuss the future of omics applications in aquaculture, natural pharmaceutical biodiscovery and pest management, as well as to monitor the impact of environmental stressors.

Proteomic similarity of the Littorinid snails in the evolutionary context

BACKGROUND

The introduction of DNA-based molecular markers made a revolution in biological systematics. However, in cases of very recent divergence events, the neutral divergence may be too slow, and the analysis of adaptive part of the genome is more informative to reconstruct the recent evolutionary history of young species. The advantage of proteomics is its ability to reflect the biochemical machinery of life. It may help both to identify rapidly evolving genes and to interpret their functions.

METHODS

Here we applied a comparative gel-based proteomic analysis to several species from the gastropod family Littorinidae. Proteomes were clustered to assess differences related to species, geographic location, sex and body part, using data on presence/absence of proteins in samples and data on protein occurrence frequency in samples of different species. Cluster support was assessed using multiscale bootstrap resampling and the stability of clustering-using cluster-wise index of cluster stability. Taxon-specific protein markers were derived using IndVal method. Proteomic trees were compared to consensus phylogenetic tree (based on neutral genetic markers) using estimates of the Robinson-Foulds distance, the Fowlkes-Mallows index and cophenetic correlation.

RESULTS

Overall, the DNA-based phylogenetic tree and the proteomic similarity tree had consistent topologies. Further, we observed some interesting deviations of the proteomic littorinid tree from the neutral expectations. (1) There were signs of molecular parallelism in two Littoraria species that phylogenetically are quite distant, but live in similar habitats. (2) Proteome divergence was unexpectedly high between very closely related Littorina fabalis and L. obtusata, possibly reflecting their ecology-driven divergence. (3) Conservative house-keeping proteins were usually identified as markers for cryptic species groups ("saxatilis" and "obtusata" groups in the Littorina genus) and for genera (Littoraria and Echinolittorina species pairs), while metabolic enzymes and stress-related proteins (both potentially adaptively important) were often identified as markers supporting species branches. (4) In all five Littorina species British populations were separated from the European mainland populations, possibly reflecting their recent phylogeographic history. Altogether our study shows that proteomic data, when interpreted in the context of DNA-based phylogeny, can bring additional information on the evolutionary history of species.

Spermatogenesis and lobular cyst type of testes organization in marine gastropod Littorina saxatilis (Olivi 1792)

Although individual stages of spermatogenesis in Littorina saxatilis are well studied at the electron microscopic level, the gonad structure and the spatial localization of gametes at different stages of maturation remain unclear. Using differential-interference contrast (DIC) for observations of fresh tissue we show that the mature testis consists of numerous ovoid lobules forming larger lobes. The lobules of intact mature testes of L. saxatilis are filled with randomly arranged multicellular cysts containing gametes at different stages of maturation. Gametes within a cyst are highly synchronized in respect of the differentiation degree. At the same time, no spatial gradient in the arrangement of cysts according to the maturation degree of gametes in them was observed in any of the studied lobules. The male gonads contain cysts with early spermatids, mid, late spermatids, and spermatozoa. Using silver-staining, DAPI, and chromomycin A3 (CMA3) staining, we identify 20 main types of nucleus organization in differentiating sperm. Premature and mature male gonads contain cysts with a mosaic arrangement as well as rare solitary cyst cells, goniablast cysts, or separate spermatogonia in between them. Our data indicate that the testis structure in L. saxatilis cannot be attributed to the tubular type, as previously thought. It corresponds to the lobular cyst type but individual lobules contain cysts with gametes at the same stage of development. It is similar to the testis structure of several fishes, amphibians, and Drosophila melanogaster. This type of the gonad organization has never been described in gastropods before.

The evolution of ependymin-related proteins

Background

Ependymins were originally defined as fish-specific secreted glycoproteins involved in central nervous system plasticity and memory formation. Subsequent research revealed that these proteins represent a fish-specific lineage of a larger ependymin-related protein family (EPDRs). EPDRs have now been identified in a number of bilaterian animals and have been implicated in diverse non-neural functions. The recent discoveries of putative EPDRs in unicellular holozoans and an expanded EPDR family with potential roles in conspecific communication in crown-of-thorns starfish suggest that the distribution and diversity of EPDRs is significantly broader than currently understood.

Results

We undertook a systematic survey to determine the distribution and evolution of EPDRs in eukaryotes. In addition to Bilateria, EPDR genes were identified in Cnidaria, Placozoa, Porifera, Choanoflagellatea, Filasterea, Apusozoa, Amoebozoa, Charophyta and Percolozoa, and tentatively in Cercozoa and the orphan group Malawimonadidae. EPDRs appear to be absent from prokaryotes and many eukaryote groups including ecdysozoans, fungi, stramenopiles, alveolates, haptistans and cryptistans. The EPDR family can be divided into two major clades and has undergone lineage-specific expansions in a number of metazoan lineages, including in poriferans, molluscs and cephalochordates. Variation in a core set of conserved residues in EPDRs reveals the presence of three distinct protein types; however, 3D modelling predicts overall protein structures to be similar.

Conclusions

Our results reveal an early eukaryotic origin of the EPDR gene family and a dynamic pattern of gene duplication and gene loss in animals. This research provides a phylogenetic framework for the analysis of the functional evolution of this gene family.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1306-y) contains supplementary material, which is available to authorized users.

Impacts of the combined exposure to seawater acidification and arsenic on the proteome of Crassostrea angulata and Crassostrea gigas

Proteomic analysis was performed to compare the effects of Arsenic (As), seawater acidification (Low pH) and the combination of both stressors (Low pH + As) on Crassostrea angulata and Crassostrea gigas juveniles in the context of global environmental change. This study aimed to elucidate if two closely related Crassostrea species respond similarly to these environmental stressors, considering both single and combined exposures, to infer if the simultaneous exposure to both stressors induced a differentiated response. Identification of the most important differentially expressed proteins between conditions revealed marked differences in the response of each species towards single and combined exposures, evidencing species-related differences towards each experimental condition. Moreover, protein alterations observed in the combined exposure (Low pH + As) were substantially different from those observed in single exposures. Identified proteins and their putative biological functions revealed an array of modes of action in each condition. Among the most important, those involved in cellular structure (Actin, Atlastin, Severin, Gelsolin, Coronin) and extracellular matrix modulation (Ependymin, Tight junction ZO-1, Neprilysin) were strongly regulated, although in different exposure conditions and species. Data also revealed differences regarding metabolic modulation capacity (ATP β, Enolase, Aconitate hydratase) and oxidative stress response (Aldehyde dehydrogenase, Lactoylglutathione, Retinal dehydrogenase) of each species, which also depended on single or combined exposures, illustrating a different response capacity of both oyster species to the presence of multiple stressors. Interestingly, alterations of piRNA abundance in C. angulata suggested genome reconfiguration in response to multiple stressors, likely an important mode of action related to adaptive evolution mechanisms previously unknown to oyster species, which requires further investigation. The present findings provide a deeper insight into the complexity of C. angulata and C. gigas responses to environmental stress at the proteome level, evidencing different capacities to endure abiotic changes, with relevance regarding the ecophysiological fitness of each species and competitive advantages in a changing environment.