Indication of ongoing amphipod speciation in Lake Baikal by genetic structures within endemic species

Museum Genomics Reveals the Hybrid Origin of an Extinct Crater Lake Endemic

Crater lake fishes are common evolutionary model systems, with recent studies suggesting a key role for gene flow in promoting rapid adaptation and speciation. However, the study of these young lakes can be complicated by human-mediated extinctions. Museum genomics approaches integrating genetic data from recently extinct species are, therefore, critical to understanding the complex evolutionary histories of these fragile systems. Here, we examine the evolutionary history of an extinct Southern Hemisphere crater lake endemic, the rainbowfish Melanotaenia eachamensis. We undertook a comprehensive sampling of extant rainbowfish populations of the Atherton Tablelands of Australia alongside historical museum material to understand the evolutionary origins of the extinct crater lake population and the dynamics of gene flow across the ecoregion. The extinct crater lake species is genetically distinct from all other nearby populations due to historic introgression between 2 proximate riverine lineages, similar to other prominent crater lake speciation systems, but this historic gene flow has not been sufficient to induce a species flock. Our results suggest that museum genomics approaches can be successfully combined with extant sampling to unravel complex speciation dynamics involving recently extinct species.

Lake Baikal amphipods and their genomes, great and small

Endemic amphipods (Crustacea: Amphipoda) of Lake Baikal represent an outstanding example of large species flocks occupying a wide range of ecological niches and originating from a handful of ancestor species. Their development took place at a restricted territory and is thus open for comprehensive research. Such examples provide unique opportunities for studying behavioral, anatomic, or physiological adaptations in multiple combinations of environmental conditions and thus attract considerable attention. The existing taxonomies of this group list over 350 species and subspecies, which, according to the molecular phylogenetic studies of marker genes, full transcriptomes and mitochondrial genomes, originated from at least two introductions into the lake. The studies of allozymes and marker genes have revealed a significant cryptic diversity in Baikal amphipods, as well as a large variance in genetic diversity within some morphological species. Crossing experiments conducted so far for two morphological species suggest that the differences in the mitochondrial marker (cytochrome c oxidase subunit I gene) can potentially be applied for making predictions about reproductive isolation. For about one-tenth of the Baikal amphipod species, nuclear genome sizes and chromosome numbers are known. While genome sizes vary within one order of magnitude, the karyotypes are relatively stable (2n = 52 for most species studied). Moreover, analysis of the diversity of repeated sequences in nuclear genomes showed significant between-species differences. Studies of mitochondrial genomes revealed some unusual features, such as variation in length and gene order, as well as duplications of tRNA genes, some of which also underwent remolding (change in anticodon specificity due to point mutations). The next important steps should be (i) the assembly of whole genomes for different species of Baikal amphipods, which is at the moment hampered by complicated genome structures with high repeat content, and (ii) updating species taxonomy taking into account all the data.

Leeches Baicalobdella torquata feed on hemolymph but have a low effect on the cellular immune response of amphipod Eulimnogammarus verrucosus from Lake Baikal

Lake Baikal is one of the largest and oldest freshwater reservoirs on the planet with a huge endemic diversity of amphipods (Amphipoda, Crustacea). These crustaceans have various symbiotic relationships, including the rarely described phenomenon of leech parasitism on amphipods. It is known that leeches feeding on hemolymph of crustacean hosts can influence their physiology, especially under stressful conditions. Here we show that leeches Baicalobdella torquata (Grube, 1871) found on gills of Eulimnogammarus verrucosus (Gerstfeldt, 1858), one of the most abundant amphipods in the Baikal littoral zone, indeed feed on the hemolymph of their host. However, the leech infection had no effect on immune parameters such as hemocyte concentration or phenoloxidase activity and also did not affect glycogen content. The intensity of hemocyte reaction to foreign bodies in a primary culture was identical between leech-free and leech-infected animals. Artificial infection with leeches also had only a subtle effect on the course of a model microbial infection in terms of hemocyte concentration and composition. Despite we cannot fully exclude deleterious effects of the parasites, our study indicates a low influence of a few leeches on E. verrucosus and shows that leech-infected amphipods can be used at least for some types of ecophysiological experiments.

Comparison of the evolutionary patterns of DNA repeats in ancient and young invertebrate species flocks of Lake Baikal

DNA repeat composition of low coverage (0.1-0.5) genomic libraries of four amphipods species endemic to Lake Baikal (East Siberia) and four endemic gastropod species of the fam. Baicaliidae have been compared to each other. In order to do so, a neighbor joining tree was inferred for each quartet of species (amphipods and mollusks) based on the ratio of repeat classes shared in each pair of species. The topology of this tree was compared to the phylogenies inferred for the same species from the concatenated protein-coding mitochondrial nucleotide sequences. In all species analyzed, the fraction of DNA repeats involved circa half of the genome. In relatively more ancient amphipods (most recent common ancestor, MRCA, existed approximately sixty millions years ago), the most abundant were species-specific repeats, while in much younger Baicaliidae (MRCA equal to ca. three millions years) most of the DNA repeats were shared among all four species. If the presence/absence of a repeat is regarded as a separate independent trait, and the ratio of shared to total numbers of repeats in a species pair is used as the measure of distance, the topology of the NJ tree is the same as the quartet phylogeny inferred for the mitogenomes protein coding nucleotide sequences. Meanwhile, in each group of species, a substantial number of repeats were detected pointing to the possibility of non-neutral evolution or a horizontal transfer between species occupying the same biotope. These repeats were shared by non-sister groups while being absent in the sister genomes. On the other hand, in such cases some traits of ecological significance were also shared.

Cold-adapted amphipod species upon heat stress: Proteomic responses and their correlation with transcriptomic responses

The cellular heat shock response (HSR) comprises transcriptomic and proteomic reactions to thermal stress. It was here addressed, how the proteomic, together with the transcriptomic HSR, relate to the thermal sensitivities of three cold-adapted but differently thermo-sensitive freshwater amphipod species. The proteomes of thermosensitive Eulimnogammarus verrucosus and thermotolerant Eulimnogammarus cyaneus, both endemic to Lake Baikal, and of thermotolerant Holarctic Gammarus lacustris were investigated upon 24 h exposure to the species-specific 10 % lethal temperatures (LT10). Furthermore, correlations of heat stress induced changes in proteomes (this study) and transcriptomes (previous study with identical experimental design) were examined. Proteomes indicated that the HSR activated processes encompassed (i) proteostasis maintenance, (ii) maintenance of cell adhesion, (iii) oxygen transport, (iv) antioxidant response, and (v) regulation of protein synthesis. Thermo-sensitive E. verrucosus showed the most pronounced proteomic HSR and the lowest correlation of transcriptomic and proteomic HSRs. For proteins related to translation (ribosomal proteins, elongation factors), transcriptomic and proteomic changes were inconsistent: transcripts were downregulated in many cases, with levels of corresponding proteins remaining unchanged. In the Eulimnogammarus species, levels of hemocyanin protein but not transcript were increased upon heat stress, suggesting a HSR also directed to enhance oxygen transport. Thermosensitive E. verrucosus showed the most pronounced relocation of transcription/translation activity to proteostasis maintenance, which may indicate that the general species-specific stability of protein structure could be a fundamental determinant of thermotolerance. By combining transcriptomic and proteomic response data, this study provides a comprehensive picture of the cellular HSR components in the studied amphipods.

Experimental Crossing Confirms Reproductive Isolation between Cryptic Species within Eulimnogammarus verrucosus (Crustacea: Amphipoda) from Lake Baikal

Ancient lakes are known speciation hotspots. One of the most speciose groups in the ancient Lake Baikal are gammaroid amphipods (Crustacea: Amphipoda: Gammaroidea). There are over 350 morphological species and subspecies of amphipods in Baikal, but the extent of cryptic variation is still unclear. One of the most common species in the littoral zone of the lake, Eulimnogammarus verrucosus (Gerstfeldt, 1858), was recently found to comprise at least three (pseudo)cryptic species based on molecular data. Here, we further explored these species by analyzing their mitogenome-based phylogeny, genome sizes with flow cytometry, and their reproductive compatibility. We found divergent times of millions of years and different genome sizes in the three species (6.1, 6.9 and 8 pg), further confirming their genetic separation. Experimental crossing of the western and southern species, which are morphologically indistinguishable and have adjacent ranges, showed their separation with a post-zygotic reproductive barrier, as hybrid embryos stopped developing roughly at the onset of gastrulation. Thus, the previously applied barcoding approach effectively indicated the separate biological species within E. verrucosus. These results provide new data for investigating genome evolution and highlight the need for precise tracking of the sample origin in any studies in this morphospecies.

Cellular Immune Response of an Endemic Lake Baikal Amphipod to Indigenous Pseudomonas sp

Studies of invertebrates have shown that the internal environment of crustaceans is not always sterile in normal conditions, and in many species, it can be populated by microorganisms even in the absence of any visible pathological processes in the body. This observation raises the question of whether genetically modified indigenous hemolymph microorganisms can be used for biotechnological purposes inside the crustacean either as local producers of some compounds or as sensors to physiological parameters. In this study, we tested the ability of the bacteria isolated from the hemolymph of the amphipod Eulimnogammarus verrucosus to hide from the cellular immune response of the host as the most important feature for their potential long-term application in vivo. 16S rDNA amplicon sequencing revealed five common bacterial genera in all analyzed samples of the amphipod hemolymph, among which Pseudomonas is most easily subjected to genome modification and, thus, the most prospective for biotechnological application. Cultivation of Pseudomonas gave us a number of strains undoubtedly derived from the amphipod hemolymph, and one of them (belonging to the Pseudomonas fluorescens group) was chosen for further tests. The primary culture of amphipod hemocytes was used to analyze the immunogenicity of the strain and showed a pronounced reaction of the immune cells to a high amount of the bacteria within six hours. This result indicates that modulation of cellular immune response to metabolically active bacterial cells is not mandatory for the survival and wide distribution of these microorganisms in the hemolymph of numerous amphipod individuals.

Changes of cellular stress response related hsp70 and abcb1 transcript and Hsp70 protein levels in Siberian freshwater amphipods upon exposure to cadmium chloride in the lethal concentration range

The induction of cellular stress response systems, heat shock protein hsp70/Hsp70 and multixenobiotic transporter abcb1, by cadmium chloride (CdCl₂) was explored in amphipod species with different stress adaptation strategies from the Lake Baikal area. Based on the lethal concentrations (LC) of CdCl₂, the sensitivities of the different species to CdCl₂ were ranked (24 hr LC50 in mg/L CdCl₂ (mean/95% confidence interval)): Gammarus lacustris (1.7/1.3–2.4) < Eulimnogammarus cyaneus (2.9/2.1–4.0) < Eulimnogammarus verrucosus (5.7/3.8–8.7) < Eulimnogammarus vittatus (18.1/12.4–26.6). Conjugated dienes, indicating lipid peroxidation, were significantly increased after 24 hr exposures to 5 mg/L CdCl₂ only in the more CdCl₂-sensitive species G. lacustris and E. cyaneus. Upon treatment with 0.54 to 5.8 mg/L CdCl₂ for 1, 6 and 24 hrs, hsp70 transcript levels were generally more increased after the longer exposure times and in the more CdCl₂-sensitive species. Relating the CdCl₂ exposure concentrations to LCx values revealed that across the species the increases of hsp70 transcript levels were comparatively low (up to 2.6-fold) at CdCl₂ concentrations ≤LC50. Relative hsp70 transcript levels were maximally increased in E. cyaneus by 5 mg/L CdCl₂ (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\hat {=}$\end{document}= ˆLC70) at 24 hrs (9.1-fold increase above the respective control). When G. lacustris was exposed to 5 mg/L CdCl₂ (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\hat {=}$\end{document}= ˆLC90) for 24 hrs, the increase in hsp70 was in comparison to E. cyaneus considerably less pronounced (3.0-fold increase in hsp70 levels relative to control). Upon exposure of amphipods to 5 mg/L CdCl₂, increases in Hsp70 protein levels compared to untreated controls were highest in E. cyaneus at 1 and 6 hrs (5 mg/L CdCl₂\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\hat {=}$\end{document}= ˆ LC70) and in E. verrucosus at 24 hrs (5 mg/L CdCl₂\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\hat {=}$\end{document}= ˆ LC45). Thus, when the fold increases in Hsp70 protein levels in the different amphipod species were related to the respective species-specific LCx values a similar bell-shaped trend as for hsp70 transcript levels was seen across the species. Transcript levels of abcb1 in CdCl₂exposed individuals of the different amphipod species varied up to 4.7-fold in relation to the respective controls. In contrast to hsp70/Hsp70, abcb1 transcripts in CdCl₂ exposed individuals of the different amphipod species did not indicate similar levels of induction of abcb1 at equal LCx levels across the species. Induction of hsp70 and abcb1 genes and Hsp70 proteins by CdCl₂ in the lethal concentration range shows that these cellular responses are rather insensitive to CdCl₂ stress in the examined amphipod species. Furthermore, the increase of expression of these cellular defense systems at such high stress levels suggests that induction of these genes is not related to the maintenance of normal metabolism but to mitigation of the effects of severe toxic stress.