Playing Peekaboo with a Master Manipulator: Metagenetic Detection and Phylogenetic Analysis of Wolbachia Supergroups in Freshwater Invertebrates.

The infamous "master manipulators"-intracellular bacteria of the genus Wolbachia-infect a broad range of phylogenetically diverse invertebrate hosts in terrestrial ecosystems. Wolbachia has an important impact on the ecology and evolution of their host with documented effects including induced parthenogenesis, male killing, feminization, and cytoplasmic incompatibility. Nonetheless, data on Wolbachia infections in non-terrestrial invertebrates are scarce. Sampling bias and methodological limitations are some of the reasons limiting the detection of these bacteria in aquatic organisms. In this study, we present a new metagenetic method for detecting the co-occurrence of different Wolbachia strains in freshwater invertebrates host species, i.e., freshwater Arthropoda (Crustacea), Mollusca (Bivalvia), and water bears (Tardigrada) by applying NGS primers designed by us and a Python script that allows the identification of Wolbachia target sequences from the microbiome communities. We also compare the results obtained using the commonly applied NGS primers and the Sanger sequencing approach. Finally, we describe three supergroups of Wolbachia: (i) a new supergroup V identified in Crustacea and Bivalvia hosts; (ii) supergroup A identified in Crustacea, Bivalvia, and Eutardigrada hosts, and (iii) supergroup E infection in the Crustacea host microbiome community.

Taxonomic classification of the bacterial endosymbiont Wolbachia based on next-generation sequencing: is there molecular evidence for its presence in tardigrades?

We used high-throughput sequencing of 16S rRNA to test whether tardigrade species are infected with Wolbachia parasites. We applied SILVA and Greengenes databases that allowed taxonomic classification of bacterial sequences to OTUs. The results obtained from both databases differed considerably in the number of OTUs, and only the Greengenes database allowed identification of Wolbachia (infection was also supported by comparison of sequences to NCBI database). The putative bacterial endosymbiont Wolbachia was discovered only in adult eutardigrades, while bacteria identified down to the order Rickettsiales were detected in both eutardigrade eggs and adult specimens. Nevertheless, the frequency of Wolbachia in the bacterial communities of the studied eutardigrades was low. Similarly, in our positive control, i.e., a fairy shrimp Streptocephalus cafer, which was found to be infected with Wolbachia in our previous study using Sanger sequencing, only the Rickettsiales were detected. We also carried out phylogenetic reconstruction using Wolbachia sequences from the SILVA and Greengenes databases, Alphaproteobacteria putative endosymbionts and Rickettsiales OTUs obtained in previous studies on the microbial community of tardigrades, and Rickettsiales and Wolbachia OTUs obtained in the current study. Our discovery of Wolbachia in tardigrades can fuel new research to uncover the specifics of this interaction.

Integrative description of bisexual Paramacrobiotus experimentalis sp. nov. (Macrobiotidae) from republic of Madagascar (Africa) with microbiome analysis.

In a moss samples collected on Madagascar two populations of Paramacrobiotus experimentalis sp. nov. were found. Paramacrobiotus experimentalis sp. nov. with the presence of a microplacoid and areolatus type of eggs is similar to Pam. danielae, Pam. garynahi, Pam. hapukuensis, Pam. peteri, Pam. rioplatensis and Pam. savai, but it differs from them by some morphological and morphometric characters of the eggs. The p-distance between two COI haplotypes of Pam. experimentalis sp. nov. was 0.17%. In turn, the ranges of uncorrected genetic p-distances of all Paramacrobiotus species available in GenBank was from 18.27% (for Pam. lachowskae) to 25.26% (for Pam. arduus) with an average distance of 20.67%. We also found that Pam. experimentalis sp. nov. is bisexual. This observation was congruent on three levels: (i) morphological - specimen size dimorphism; (ii) structural (primary sexual characteristics) - females have an unpaired ovary while males have an unpaired testis and (iii) molecular - heterozygous and homozygous strains of the ITS-2 marker. Although symbiotic associations of hosts with bacteria (including endosymbiotic bacteria) are common in nature and these interactions exert various effects on the evolution, biology and reproductive ecology of hosts, there is still very little information on the bacterial community associated with tardigrades. To fill this gap and characterise the bacterial community of Pam. experimentalis sp. nov. populations and microbiome of its microhabitat, high throughput sequencing of the V3-V4 hypervariable regions in the bacterial 16S rRNA gene fragment was performed. The obtained 16S rRNA gene sequences ranged from 92,665 to 131,163. In total, 135 operational taxonomic units (OTUs) were identified across the rarefied dataset. Overall, both Pam. experimentalis sp. nov. populations were dominated by OTUs ascribed to the phylum Proteobacteria (89-92%) and Firmicutes (6-7%). In the case of samples from tardigrades' laboratory habitat, the most abundant bacterial phylum was Proteobacteria (51-90%) and Bacteroides (9-48%). In all compared microbiome profiles, only 16 of 137 OTUs were shared. We found also significant differences in beta diversity between the partly species-specific microbiome of Pam. experimentalis sp. nov. and its culturing environment. Two OTUs belonging to a putative bacterial endosymbiont were identified - Rickettsiales and Polynucleobacter. We also demonstrated that each bacterial community was rich in genes involved in membrane transport, amino acid metabolism, and carbohydrate metabolism.

The complete maternal and paternal mitochondrial genomes of Unio crassus: Mitochondrial molecular clock and the overconfidence of molecular dating.

The availability of a rapidly growing number of complete mitochondrial genome sequences provokes high confidence dating approaches. However, even if the congruence between mitochondrial and nuclear markers is reasonable, the resulting topologies are frequently questionable. The unique opportunity to study the evolutionary history of two independent mitochondrial genomes in one phylogenetic context exists in the freshwater mussels family Unionidae. The two lineages function under doubly uniparental inheritance since well before the emergence of the family. Despite the relatively high number of available complete sequences of maternally inherited genomes, comparative analyses are limited by the small number of sequences of counterpart paternally inherited genomes. We have sequenced for the first time the representative set of five sequences (two maternal and three paternal) from the species Unio crassus. Comparative analysis of the phylogenies reconstructed using relevant mitogenomic data available in GenBank (13 species in total) reveal that single - genome inferences are congruent only if the relaxed clock is assumed.

Morphological and genetic variability in cosmopolitan tardigrade species - Paramacrobiotus fairbanksi Schill, Förster, Dandekar & Wolf, 2010.

Paramacrobiotus fairbanksi was described from Alaska (USA) based on integrative taxonomy and later reported from various geographical locations making it a true cosmopolitan species. The 'Everything is Everywhere' (EiE) hypothesis assumes that microscopic organisms have unique features that help them to inhabit many different environments, meaning they can be considered cosmopolitan. In the present work we report four new populations of Pam. fairbanksi from the Northern Hemisphere which suggests that the 'EiE' hypothesis is true, at least for some tardigrade species. We also compared all known populations of Pam. fairbanksi at the genetic and morphological levels. The p-distances between COI haplotypes of all sequenced Pam. fairbanksi populations from Albania, Antarctica, Canada, Italy, Madeira, Mongolia, Spain, USA and Poland ranged from 0.002% to 0.005%. In total, twelve haplotypes (H1-H12) of COI gene fragments were identified. We also report statistically significant morphometrical differences of species even though they were cultured and bred in the same laboratory conditions, and propose epigenetic factor as a main cause rather than temperature, predation risk and food availability. Furthermore, we also discuss differences in the potential distribution of two Paramacrobiotus species.

Contrasting Marine Microbial Communities of the Fram Strait with the First Confirmed Record of Cyanobacteria Prochlorococcus marinus in the Arctic Region.

The seawater microbiome is crucial in marine ecosystems because of its role in food chains and biogeochemical cycles; thus, we studied the composition of the pelagic marine microbiome collected in the upper 50 m on the opposite sides of Fram Strait: Spitsbergen and Greenland shelves. We found out that it differed significantly, with salinity being the main environmental variable responsible for these differences. The Spitsbergen shelf was dominated by Atlantic Waters, with a rather homogenous water column in terms of salinity and temperature down to 300 m; hence, the marine microbial community was also homogenous at all sampled depths (0, 25, 50 m). On the contrary, stations on the Greenland shelf were exposed to different water masses of both Arctic and Atlantic origin, which resulted in a more diverse microbial community there. Unexpectedly, for the very first time, we identified cyanobacterium Prochlorococcus marinus in Arctic waters (Spitsbergen shelf, 75-77° N). Till now, the distribution of this cyanobacteria in oceans has been described only between 40° N and 40° S. Considering the accelerated rate of climate warming in the Arctic, our results indicated that the seawater microbiome can be viewed as an amplifier of global change and that the Atlantification is in progress.

Morphological and genetic variability in cosmopolitan tardigrade species-Paramacrobiotus fairbanksi Schill, Förster, Dandekar & Wolf, 2010.

Paramacrobiotus fairbanksi was described from Alaska (USA) based on integrative taxonomy and later reported from various geographical localities making it a true cosmopolitan species. The 'Everything is Everywhere' (EiE) hypothesis assumes that the geographic distribution of microscopic organisms is not limited by dispersal but by local environmental conditions, making them potentially cosmopolitan. In the present work we report four new populations of P. fairbanksi from the Northern Hemisphere which suggests that the 'EiE' hypothesis is true, at least for some tardigrade species. We also compared all known populations of P. fairbanksi at the genetic and morphological levels. The p-distances between COI haplotypes of all sequenced P. fairbanksi populations from Albania, Antarctica, Canada, Italy, Madeira, Mongolia, Spain, USA and Poland ranged from 0.002 to 0.005%. In total, twelve haplotypes (H1-H12) of COI gene fragments were identified. We also report statistically significant morphometrical differences of species even though they were cultured and bred in the same laboratory conditions. Furthermore, we also discuss differences in the potential distribution of two Paramacrobiotus species.

Integrative taxonomy reveals new, widely distributed tardigrade species of the genus Paramacrobiotus (Eutardigrada: Macrobiotidae).

In a moss sample collected in Ribeiro Frio, Madeira, Paramacrobiotus gadabouti sp. nov. was found and described using the integrative taxonomy approach. The new species is described based on morphological and morphometric data from both phase-contrast light microscopy (PCM), as well as scanning electron microscopy (SEM). Moreover, four DNA markers, three nuclear (18S rRNA, 28S rRNA, ITS-2) and one mitochondrial (COI) markers, were used to elucidate the phylogenetic position of the new species within the family Macrobiotidae. The new species has a microplacoid that placed it within Parmacrobiotus richtersi group and exhibit richtersi-type eggs having processes terminated with cap-like structures. Paramacrobiotus gadabouti sp. nov. is most similar to Pam. alekseevi, Pam. filipi and Pam. garynahi, but differs from them mainly in details of egg morphology and morphometrics. Unlike other species from this group, which were confirmed as bisexual and showed limited distribution, Paramacrobiotus gadabouti sp. nov. is yet another parthenogenetic species with a wide distribution, demonstrating that at least some tardigrades confirm to the hypothesis of 'everything is everywhere'.

An integrative description of Diploechiniscus oihonnae/ (Richters, 1903) population from near the original type locality in Merok (Norway).

In the present study, we used integrative taxonomy to describe a population of Diploechiniscus oihonnae (Richters, 1903) from the neotype locality in Merok (Norway). We found no differences in the chaetotaxy formula between the life stages and sex of Dpl. oihonnae. The presence of filaments at Bd in some specimens of Dpl. oihonnae appears to be random and most likely represents the morphological variability of this species. We also obtained DNA sequences of 18S rRNA, 28S rRNA, ITS-2, and COI of Dpl. oihonnae from the neotype locality for comparison with the sequences available in GenBank, which showed low genetic differences between the neotypic population and specimens from other localities. Therefore, we decided to establish our specimens from Merok as neotype and neoparatypes of Dpl. oihonnae. Additionally, based on morphological characters, Dpl. horningi (Schuster Grigarick, 1971) was synonymised with Dpl. oihonnae.

Pleistocene allopatric differentiation followed by recent range expansion explains the distribution and molecular diversity of two congeneric crustacean species in the Palaearctic.

Pleistocene glaciations had a tremendous impact on the biota across the Palaearctic, resulting in strong phylogeographic signals of range contraction and rapid postglacial recolonization of the deglaciated areas. Here, we explore the diversity patterns and history of two sibling species of passively dispersing taxa typical of temporary ponds, fairy shrimps (Anostraca). We combine mitochondrial (COI) and nuclear (ITS2 and 18S) markers to conduct a range-wide phylogeographic study including 56 populations of Branchinecta ferox and Branchinecta orientalis in the Palaearctic. Specifically, we investigate whether their largely overlapping ranges in Europe resulted from allopatric differentiation in separate glacial refugia followed by a secondary contact and reconstruct their postglacial recolonization from the inhabited refugia. Our results suggest the existence of distinct refugia for the two species, with genetic divergence among intraspecific lineages consistent with late Pleistocene glacial cycles. While B. ferox lineages originated from Mediterranean refugia, the origin of B. orientalis lineages was possibly located on the Pannonian Plain. We showed that most dispersal events predominantly happened within 100 km, coupled with several recent long-distance events (> 1000 km). Hence the regional habitat density of suitable habitats in Central Europe is possibly a key to the co-existence of the two species. Overall, our study illustrates how isolation in combination with stochastic effects linked to glacial periods are important drivers of the allopatric differentiation of Palaearctic taxa.

The Patterns and Puzzles of Genetic Diversity of Endangered Freshwater Mussel Unio crassus Philipsson, 1788 Populations from Vistula and Neman Drainages (Eastern Central Europe).

Mussels of the family Unionidae are important components of freshwater ecosystems. Alarmingly, the International Union for Conservation of Nature and Natural Resources Red List of Threatened Species identifies almost 200 unionid species as extinct, endangered, or threatened. Their decline is the result of human impact on freshwater habitats, and the decrease of host fish populations. The Thick Shelled River Mussel Unio crassus Philipsson, 1788 is one of the examples that has been reported to show a dramatic decline of populations. Hierarchical organization of riverine systems is supposed to reflect the genetic structure of populations inhabiting them. The main goal of this study was an assessment of the U. crassus genetic diversity in river ecosystems using hierarchical analysis. Different molecular markers, the nuclear ribosomal internal transcribed spacer ITS region, and mitochondrial DNA genes (cox1 and ndh1), were used to examine the distribution of U. crassus among-population genetic variation at multiple spatial scales (within rivers, among rivers within drainages, and between drainages of the Neman and Vistula rivers). We found high genetic structure between both drainages suggesting that in the case of the analyzed U. crassus populations we were dealing with at least two different genetic units. Only about 4% of the mtDNA variation was due to differences among populations within drainages. However, comparison of population differentiation within drainages for mtDNA also showed some genetic structure among populations within the Vistula drainage. Only one haplotype was shared among all Polish populations whereas the remainder were unique for each population despite the hydrological connection. Interestingly, some haplotypes were present in both drainages. In the case of U. crassus populations under study, the Mantel test revealed a relatively strong relationship between genetic and geographical distances. However, in detail, the pattern of genetic diversity seems to be much more complicated. Therefore, we suggest that the observed pattern of U. crassus genetic diversity distribution is shaped by both historical and current factors i.e. different routes of post glacial colonization and history of drainage systems, historical gene flow, and more recent habitat fragmentation due to anthropogenic factors.

Instances of erroneous DNA barcoding of metazoan invertebrates: Are universal cox1 gene primers too "universal"?

The cytochrome c oxidase subunit I (cox1) gene is the main mitochondrial molecular marker playing a pivotal role in phylogenetic research and is a crucial barcode sequence. Folmer's "universal" primers designed to amplify this gene in metazoan invertebrates allowed quick and easy barcode and phylogenetic analysis. On the other hand, the increase in the number of studies on barcoding leads to more frequent publishing of incorrect sequences, due to amplification of non-target taxa, and insufficient analysis of the obtained sequences. Consequently, some sequences deposited in genetic databases are incorrectly described as obtained from invertebrates, while being in fact bacterial sequences. In our study, in which we used Folmer's primers to amplify COI sequences of the crustacean fairy shrimp Branchipus schaefferi (Fischer 1834), we also obtained COI sequences of microbial contaminants from Aeromonas sp. However, when we searched the GenBank database for sequences closely matching these contaminations we found entries described as representatives of Gastrotricha and Mollusca. When these entries were compared with other sequences bearing the same names in the database, the genetic distance between the incorrect and correct sequences amplified from the same species was c.a. 65%. Although the responsibility for the correct molecular identification of species rests on researchers, the errors found in already published sequences data have not been re-evaluated so far. On the basis of the standard sampling technique we have estimated with 95% probability that the chances of finding incorrectly described metazoan sequences in the GenBank depend on the systematic group, and variety from less than 1% (Mollusca and Arthropoda) up to 6.9% (Gastrotricha). Consequently, the increasing popularity of DNA barcoding and metabarcoding analysis may lead to overestimation of species diversity. Finally, the study also discusses the sources of the problems with amplification of non-target sequences.

Detection of bacterial endosymbionts in freshwater crustaceans: the applicability of non-degenerate primers to amplify the bacterial 16S rRNA gene.

Bacterial endosymbionts of aquatic invertebrates remain poorly studied. This is at least partly due to a lack of suitable techniques and primers for their identification. We designed a pair of non-degenerate primers which enabled us to amplify a fragment of ca. 500 bp of the 16S rRNA gene from various known bacterial endosymbiont species. By using this approach, we identified four bacterial endosymbionts, two endoparasites and one uncultured bacterium in seven, taxonomically diverse, freshwater crustacean hosts from temporary waters across a wide geographical area. The overall efficiency of our new WOLBSL and WOLBSR primers for amplification of the bacterial 16S rRNA gene was 100%. However, if different bacterial species from one sample were amplified simultaneously, sequences were illegible, despite a good quality of PCR products. Therefore, we suggest using our primers at the first stage of bacterial endosymbiont identification. Subsequently, genus specific primers are recommended. Overall, in the era of next-generation sequencing our method can be used as a first simple and low-cost approach to identify potential microbial symbionts associated with freshwater crustaceans using simple Sanger sequencing. The potential to detected bacterial symbionts in various invertebrate hosts in such a way will facilitate studies on host-symbiont interactions and coevolution.

Gender-Associated Mitochondrial DNA Heteroplasmy in Somatic Tissues of the Endangered Freshwater Mussel Unio crassus (Bivalvia: Unionidae): Implications for Sex Identification and Phylogeographical Studies.

Some bivalve species possess two independent mitochondrial DNA lineages: maternally (F-type) and paternally (M-type) inherited. This phenomenon is called doubly uniparental inheritance. It is generally agreed that F-type mtDNA is typically present in female somatic and gonadal tissues as well as in male somatic tissues, whereas the M-type mtDNA occurs only in male germ line and gonadal tissue. In the present study, the mtDNA heteroplasmy (for both F and M genomes) in male somatic tissues of Unio crassus (Philipsson, 1788), species threatened with extinction, has been confirmed. Taking advantage from the presence of Mcox1 marker only in male somatic tissues, we developed a new method of sex identification in this endangered species, using nondestructive tissue sampling. Probability of correct sex identification was estimated at 97.5%. The present study is the first report on gender-associated mitochondrial DNA heteroplasmy in male somatic tissues of thick-shelled river mussel and first approach to U. crassus sex identification at molecular level. Our study also confirmed the utility of paternally inherited Mcox1 gene fragment as a complementary molecular tool for resolving phylogeographical relationships among populations of thick-shelled river mussel.

Identification and characterization of the first microsatellite loci for the thick-shelled river mussel Unio crassus (Bivalvia: Unionidae).

The thick-shelled river mussel Unio crassus is critically endangered throughout its range as a result of increasing human activity and habitat loss. Next generation DNA sequencing was used to develop a set of microsatellite markers that can be used for future ecological and population genetics studies of U. crassus. A total of 11 polymorphic loci were identified and characterized using 57 individuals from two Polish populations. Numbers of alleles ranged from 2 to 15 and expected heterozygosity levels ranged from 0.069 to 0.899. Deficiency of heterozygous genotypes was observed in four loci. Marker independence was confirmed with tests for linkage disequilibrium, however, analyses indicated evidence of null alleles in four loci. The microsatellite markers developed specifically for U. crassus provide a valuable tool for future ecological, population genetic assessments, and conservation management of this species.