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

Consistent Pattern of Local Adaptation during an Experimental Heat Wave in a Pipefish-Trematode Host-Parasite System

Extreme climate events such as heat waves are expected to increase in frequency under global change. As one indirect effect, they can alter magnitude and direction of species interactions, for example those between hosts and parasites. We simulated a summer heat wave to investigate how a changing environment affects the interaction between the broad-nosed pipefish (Syngnathus typhle) as a host and its digenean trematode parasite (Cryptocotyle lingua). In a fully reciprocal laboratory infection experiment, pipefish from three different coastal locations were exposed to sympatric and allopatric trematode cercariae. In order to examine whether an extreme climatic event disrupts patterns of locally adapted host-parasite combinations we measured the parasite's transmission success as well as the host's adaptive and innate immune defence under control and heat wave conditions. Independent of temperature, sympatric cercariae were always more successful than allopatric ones, indicating that parasites are locally adapted to their hosts. Hosts suffered from heat stress as suggested by fewer cells of the adaptive immune system (lymphocytes) compared to the same groups that were kept at 18°C. However, the proportion of the innate immune cells (monocytes) was higher in the 18°C water. Contrary to our expectations, no interaction between host immune defence, parasite infectivity and temperature stress were found, nor did the pattern of local adaptation change due to increased water temperature. Thus, in this host-parasite interaction, the sympatric parasite keeps ahead of the coevolutionary dynamics across sites, even under increasing temperatures as expected under marine global warming

Molecular signatures of the rediae, cercariae and adult stages in the complex life cycles of parasitic flatworms (Digenea: Psilostomatidae)

BACKGROUND: Parasitic flatworms (Trematoda: Digenea) represent one of the most remarkable examples of drastic morphological diversity among the stages within a life cycle. Which genes are responsible for extreme differences in anatomy, physiology, behavior, and ecology among the stages? Here we report a comparative transcriptomic analysis of parthenogenetic and amphimictic generations in two evolutionary informative species of Digenea belonging to the family Psilostomatidae. METHODS: In this study the transcriptomes of rediae, cercariae and adult worm stages of Psilotrema simillimum and Sphaeridiotrema pseudoglobulus, were sequenced and analyzed. High-quality transcriptomes were generated, and the reference sets of protein-coding genes were used for differential expression analysis in order to identify stage-specific genes. Comparative analysis of gene sets, their expression dynamics and Gene Ontology enrichment analysis were performed for three life stages within each species and between the two species.RESULTS: Reference transcriptomes for P. simillimum and S. pseudoglobulus include 21,433 and 46,424 sequences, respectively. Among 14,051 orthologous groups (OGs), 1354 are common and specific for two analyzed psilostomatid species, whereas 13 and 43 OGs were unique for P. simillimum and S. pseudoglobulus, respectively. In contrast to P. simillimum, where more than 60% of analyzed genes were active in the redia, cercaria and adult worm stages, in S. pseudoglobulus less than 40% of genes had such a ubiquitous expression pattern. In general, 7805 (36.41%) and 30,622 (65.96%) of genes were preferentially expressed in one of the analyzed stages of P. simillimum and S. pseudoglobulus, respectively. In both species 12 clusters of co-expressed genes were identified, and more than a half of the genes belonging to the reference sets were included into these clusters. Functional specialization of the life cycle stages was clearly supported by Gene Ontology enrichment analysis.CONCLUSIONS: During the life cycles of the two species studied, most of the genes change their expression levels considerably, consequently the molecular signature of a stage is not only a unique set of expressed genes, but also the specific levels of their expression. Our results indicate unexpectedly high level of plasticity in gene regulation between closely related species. Transcriptomes of P. simillimum and S. pseudoglobulus provide high quality reference resource for future evolutionary studies and comparative analyses

New Data on Spermatogenic Cyst Formation and Cellular Composition of the Testis in a Marine Gastropod, Littorina saxatilis

Knowledge of the testis structure is important for gastropod taxonomy and phylogeny, particularly for the comparative analysis of sympatric Littorina species. Observing fresh tissue and squashing fixed tissue with gradually increasing pressure, we have recently described a peculiar type of cystic spermatogenesis, rare in mollusks. It has not been documented in most mollusks until now. The testis of adult males consists of numerous lobules filled with multicellular cysts containing germline cells at different stages of differentiation. Each cyst is formed by one cyst cell of somatic origin. Here, we provide evidence for the existence of two ways of cyst formation in Littorina saxatilis. One of them begins with a goniablast cyst formation; it somewhat resembles cyst formation in Drosophila testes. The second way begins with capture of a free spermatogonium by the polyploid cyst cell which is capable to move along the gonad tissues. This way of cyst formation has not been described previously. Our data expand the understanding of the diversity of spermatogenesis types in invertebrates

Data on RNA-seq analysis of the oviducts of five closely related species genus <i>Littorina</i> (Mollusca, Caenogastropoda): <i>L. saxatilis, L. arcana, L. compressa, L. obtusata, L. fabalis</i>.

In the evolution of invertebrates, the transition from egg-layers to brooders occurred many times. However, the molecular mechanisms underlying this transition are still not well understood. Recently diverged species genus Littorina (Mollusca, Gastropoda, Caenogastropoda, Littorinimorpha): Littorina saxatilis, L. arcana, L. compressa, L. obtusata and L. fabalis might be a fruitful model for elucidation of these mechanisms. All five species sympatrically inhabit an intertidal zone. Only L. saxatilis is ovoviviparous while the other four species form clutches. Although in L. saxatilis jelly gland of the pallial oviduct function as a brood pouch, it is not deeply modified at the morphological level in comparison to egg-laying relatives. Comparative analysis of transcriptomic profiles of the pallial oviducts of these closely related species might help to uncover the molecular mechanisms of the egg-laying to brooding transition. Unraveling of the mechanisms underlying this transition in L. saxatilis is important not only in aspects of reproduction biology and strategy, but also in a broader view as an example of relatively fast evolutionary transformations. We generated an RNA-seq dataset (224 104 446 clean reads) for oviducts of five species genus Littorina. Libraries of all five species were sequenced using Illumina HiSeq 2500; additional reads for L. arcana were obtained using Illumina NovaSeq 6000. Transcriptomic profiles were analyzed in pooled samples (of three individuals) with two biological replicates for each species (each biological replicate was prepared and sequenced as a separate library). The transcriptome was assembled de novo and annotated with five assembles corresponding to each species. The raw data were uploaded to the SRA database, the BioProject IDs are PRJNA662103 ("obtusata" group) and PRJNA707549 ("saxatilis" group)