Loss of genetic diversity means loss of geological information: the endangered Japanese crayfish exhibits remarkable historical footprints

Environmental DNA unveils deep phylogeographic structure of a freshwater fish

Phylogeography bears an important part in ecology and evolution. However, current phylogeographic studies are largely constrained by limited numbers of individual samples. Using an environmental DNA (eDNA) assay for phylogeographic analyses, this study provides detailed information regarding the history of Siberian stone loach Barbatula toni, a primary freshwater fish across the whole range of Hokkaido, Japan. Based on an eDNA metabarcoding on 293 river water samples, we detected eDNA from B. toni in 189 rivers. A total of 51 samples, representing the entire island, were then selected from the B. toni eDNA-positive sample set for the subsequent analyses. To elucidate the phylogeographic structure of B. toni, newly developed eDNA metabarcoding primers (Barba-cytb-F/R) were applied to these samples, specifically targeting their haplotypic variation in cytochrome b. After a bioinformatic processing to mitigate haplotypic false positives, a total of 50 eDNA haplotypes were identified. Two regionally restricted, genetically distinct lineages of the species were revealed as a result of phylogeographic analyses on the haplotypes and tissue-derived DNA from B. toni. According to a molecular clock analysis, they have been genetically isolated for at least 1.5 million years, suggesting their ancient origin and colonisation of Hokkaido, presumably in the glacial periods. These results demonstrate how freshwater fishes can alter their distributions over evolutionary timescales and how eDNA assay can deepen our understanding of phylogeography.

Evolutionary history of a cold-adapted limnephilid caddisfly: Effects of climate change and topography on genetic structure

The distribution of organisms is influenced by complex factors such as the phylogenetic evolutionary histories of species, the physiological and ecological characteristics of organisms, climate, and geographical and geohistorical features. In this study, we focused on a caddisfly, Asynarchus sachalinensis (Trichoptera: Limnephilidae), which has adapted to cold habitats. From phylogeographic analyses based on the mitochondrial DNA (mtDNA) cytochrome c oxidase subunit I (COI) and 16S rRNA regions and the nuclear DNA (nDNA) 18S rRNA, 28S rRNA, carbamoyl-phosphate synthetase (CAD), elongation factor-1 alpha (EF1-α), and RNA polymerase II (POLII) regions, two distinct genetic clades were detected. Clade I was shown to be widely distributed from Sakhalin to Honshu, whereas Clade II was only distributed within Honshu. The distributions of these clades overlapped in Honshu. The habitats were located at relatively lower altitudes for Clade I and higher altitudes for Clade II. The divergence time of these clades was estimated to be during the Pleistocene, indicating that repeated climatic changes facilitated distributional shifts. Haplotype network and demographic analyses based on the mtDNA COI region showed contrasting genetic structures in the two clades. It was indicated that the population sizes of Clade I had expanded rapidly in a recent period, whereas Clade II had maintained stable population sizes. The habitats of Clade II were typically isolated and scattered at high altitudes, resulting in restricted migration and dispersal because of their discontinuous "Sky Island" habitats. The habitats of Clade I were located at relatively low altitudes, and it was assumed that the populations were continuous, which resulted in a higher frequency of migration and dispersal between populations. Thus, differences in the spatial scale of the adapted habitats of each clade may have resulted in different patterns of population connectivity and fragmentation associated with repeated climatic changes during the Pleistocene. Our study provided new insight into the distributional patterns of cold-adapted aquatic insects in the Japanese Archipelago. Furthermore, the distributional shifts predicted by ecological niche modeling under future climatic change conditions were different for each clade. Therefore, different principles are required in the assessment of each clade to predict temporal changes in their distributions.

Phylogeography of the true freshwater crab, Geothelphusa dehaani: Detected dual dispersal routes via land and sea

Dispersal is an important factor that determines the potential for colonization to pioneer sites. Although most decapods employ seaward migration for reproduction with a planktonic larval phase, true freshwater crabs spend their entire life cycle in freshwater. Therefore, it is expected that genetic regionality can be easily detected. In this study, we focused on true freshwater crabs, Geothelphusa Stimpson, 1858. Herein, we reveal the evolutionary history and dispersal patterns of freshwater crustaceans. We collected and genetically analyzed 283 specimens at 138 localities across the Japanese Islands. Phylogenetic analyses were conducted on the combined dataset (mtDNA COI, 16S, and nDNA ITS1, histone H3 regions) and the data set based on the mtDNA COI region. The phylogenetic relationships detected 10 clades that were highly monophyletic. The highlights of this study were the discovery of several cryptic species or undescribed species, and the completely different heterogeneous dual dispersal pathways within a single species; i.e., both land and ocean routes. Although it was concluded that Japanese crabs are basically genetically divided by straits, strong evidence for dispersion via ocean currents was also detected (i.e., a "sweepstake"). It was also confirmed that Geothelphusa dehaani (White, 1847) could survive in seawater.

Origin of scarlet gynogenetic triploid Carassius fish: Implications for conservation of the sexual-gynogenetic complex

Conservation of sperm-dependent asexual (gynogenetic) species is challenging due to their complicated ecological dynamics, which requires the stable coexistence with their sperm-providing sexual relatives, who often share similar niches. A symbolic but vulnerable gynogenetic animal is the scarlet Carassius fish, or Hibuna, which is mainly found in Lake Harutori on Hokkaido, Japan. Although Hibuna in Lake Harutori has been protected as a symbol of the Natural Monument of Japan, it has recently suffered population decline. To establish an effective conservation strategy for Hibuna, we investigated its origin, reproductive mode, and genetic diversity, with reference to the surrounding wild populations, using nuclear microsatellites and mitochondrial gene sequences. Our genetic analyses revealed that the main ploidy of Hibuna was triploid or tetraploid, and it reproduces gynogenetically. However, no co-existing sexual diploid Carassius was detected among our samples, suggesting that the sexual diploids and the gynogenetic population including Hibuna would be at risk of co-extirpation. In addition, Hibuna showed high genetic/clonal diversity and most Hibuna had nonindigenous mitochondrial haplotypes that are mostly identical to those reported from goldfish. These results indicate that Hibuna most probably originated from hybridization between indigenous gynogenetic triploids and goldfish introduced about 100 years ago, involving rare sexual reproduction. This spontaneous long-term field experiment exemplifies the recently documented diversification process of gynogenetic Carassius via complex interploidy gene flow. Although the priority to be placed on the conservation of Hibuna is controversial, the maintenance of gynogenetic Carassius, including Hibuna, requires strategic conservation of sexual populations.

Phylogeography of Kisaura Ross (Trichoptera: Philopotamidae) of the Japanese Archipelago and the character displacement evolution observed in a secondary contact area between genetically differentiated intra-specific lineages

In this study we present evidence of an ongoing speciation event in the Japanese Archipelago. The Kisaura nozakii species complex (Philopotamidae) is a group of a small trichopteran insects, and it has been reported that this species has large differences between populations in the copulatory organs (genitalia) of males. Therefore, we conducted molecular phylogenetic analyses based on the mtDNA COI region (658 bp) and genome-wide nDNA single nucleotide polymorphisms (SNPs) (16 254 loci) using the K. nozakii species complex and a sister-species, K. borealis. We also conducted principal component analyses (PCA) based on the forewing length and seven morphological characteristics of the male genitalia. It was revealed that the K. nozakii species complex is composed of multiple allopatrically differentiated genetic lineages. However, two genetic lineages were distributed sympatrically/parapatrically in the western area of the Japanese Archipelago, and the occurrence of ‘character displacement’ in the male genitalia was detected in this area. We consider that the evolutionary event in this species complex has high potential to become a model case for elucidating the processes of speciation and species diversification.

Genetic Diversity of Ligidium Isopods in Hokkaido and Niigata, Northern Japan, Based on Mitochondrial DNA Analysis

The genetic diversity of the genus Ligidium in Hokkaido and Niigata, northern Japan, was investigated by analyzing the cytochrome c oxidase subunit 1 (CO1) region in the mitochondrial DNA (mtDNA). The genetic diversity in Hokkaido was much lower than that in Niigata. Nine different operational taxonomic units (OTUs) were identified. Only a single OTU, most likely Ligidium japonicum, was found in Hokkaido, whereas all nine OTUs were found in Niigata. Using the mtDNA evolutionary rate determined for the marine invertebrate Haptosquilla pulchella (Miers, 1880), population expansion for OTU1 in Hokkaido was estimated to have occurred at 12,600 years BP, suggesting that Ligidium underwent a bottleneck due to glacial cooling, and the population then expanded after postglacial warming. Assuming that the expansion of the OTU1 population occurred at 9600 years BP, when the sea surface temperature rose offshore of Tokachi in the Northwestern Pacific, the evolutionary rate (µ) of the mtDNA CO1 region in Ligidium is calculated as: 0.087 (95% confidence intervals: min: 0.042-max: 0.12) (substitutions/site/million years). The presence of a haplotype common to Hokkaido and Niigata implies that the haplotype migrated across the Tsugaru Strait. Considering that geological evidence indicates that the Tsugaru Strait was continuously present even during the last glacial maximum when the sea level was at its lowest, accidental transport by human beings or animals might have been critical to the migration of Ligidium.

The future of endangered crayfish in light of protected areas and habitat fragmentation

The long-term survival of a species requires, among other things, gene flow between populations. Approaches for the evaluation of fragmentation in the frame of freshwater habitats consider only a small amount of the information that combined demography and geography are currently able to provide. This study addresses two species of Austropotamobius crayfish in the light of population genetics, spatial ecology and protected areas of the Carpathians. Advancing the classical approaches, we defined ecological distances upon the rasterised river network as a surrogate of habitat resistance to migration, quantifying the deviations from the species´ suitability range for a set of relevant geospatial variables in each cell of the network. Molecular analyses revealed the populations of the two Austropotamobius crayfish species are clearly distinct, lacking hybridisation. Comparing pairs of populations, we found, in some cases, a strong disagreement regarding genetic and ecological distances, potentially due to human-mediated translocations or the geophysical phenomena of regressive erosion, which may have led to unexpected colonisation routes. Protected areas were found to offer appropriate local habitat conditions but failed to ensure connectivity. The methodology applied in this study allowed us to quantify the contribution of each geospatial (environmental) variable to the overall effect of fragmentation, and we found that water quality was the most important variable. A multilevel approach proved to reveal a better understanding of drivers behind the distribution patterns, which can lead to more adequate conservation measures.

Next generation sequencing-aided comprehensive geographic coverage sheds light on the status of rare and extinct populations of Aporia butterflies (Lepidoptera: Pieridae)

The Black-veined White Aporia crataegi (Linnaeus, 1758), a common and widespread butterfly ranging from northwestern Africa to Europe and Asia, has been extinct in Britain since the 1920s and is on a steady decline in several other parts of its range. In order to investigate genetic diversity within A. crataegi and its correspondence with current subspecies-level taxonomy, we barcoded 173 specimens from across its range including, for the first time, extinct populations from Britain and Korea. Using next generation sequencing we also obtained a sequence for Aporia joubini, a peculiar taxon from China known only by its type specimen collected in the early twentieth century. Our phylogenetic analysis placed A. joubini sister to A. oberthuri, although further taxon sampling may reveal a different scheme. Within A. crataegi, we observed a shallow and weak mitogenomic structure with only a few distinct lineages in North Africa, Sicily, Iran, and Japan. Eurasian populations, including those extinct in Britain and Korea, clustered into a large set of closely allied lineages, consistent with a recent expansion during the Late Pleistocene glacial period. This study highlights the importance of museum collections and the unique opportunities they provide in documenting species diversity and helping conservation efforts.

Genealogical Position of Japanese Populations of the Globally Distributed Mayfly Cloeon dipterum and Related Species (Ephemeroptera, Baetidae): A Molecular Phylogeographic Analysis

In the present study, we add genetic data of the mayfly Cloeon dipterum collected from the Japanese Islands to the established molecular phylogenetic knowledge in the mitochondrial COI gene of Cloeon mayflies. Cloeon dipterum is a typical cosmopolitan species that includes six intraspecific haplotype groups. The present phylogenetic analysis revealed that haplotypes of the Japanese C. dipterum constitute a seventh group together with a haplotype from Korea. This East Asian group forms a sister group with previously known European and North American haplotype groups (i.e., the clade CT1 to CT3). The present phylogenetic analysis further revealed the occurrence of two described species (C. dipterum and Cloeon ryogokuensis) and possibly three species (Cloeon sp. 1 to 3) in Japan. Consideration is given to the degree of genetic differentiation, divergence time, and differentiation process among these seven genetic groups.

Phylogeography of freshwater planorbid snails reveals diversification patterns in Eurasian continental islands

Background

Islands have traditionally been the centre of evolutionary biological research, but the dynamics of immigration and differentiation at continental islands have not been well studied. Therefore, we focused on the Japanese archipelago, the continental islands located at the eastern end of the Eurasian continent. While the Japanese archipelago is characterised by high biodiversity and rich freshwater habitats, the origin and formation mechanisms of its freshwater organisms are not clear. In order to clarify the history of the planorbid gastropod fauna, we conducted phylogenetic analysis, divergence time estimation, ancestral state reconstruction, and lineage diversity estimations.

Results

Our analyses revealed the formation process of the planorbid fauna in the Japanese archipelago. Most lineages in the Japanese archipelago have closely related lineages on the continent, and the divergence within the Japanese lineages presumably occurred after the late Pliocene. In addition, each lineage is characterised by different phylogeographical patterns, suggesting that immigration routes from the continent to the Japanese archipelago differ among lineages. Furthermore, a regional lineage diversity plot showed that the present diversity in the Japanese archipelago potentially reflects the differentiation of lineages within the islands after the development of the Japanese archipelago.

Conclusions

Although additional taxon sampling and genetic analysis focused on each lineage are needed, our results suggest that immigration from multiple routes just prior to the development of the Japanese archipelago and subsequent diversification within the islands are major causes of the present-day diversity of the Japanese planorbid fauna.

Electronic supplementary material

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

ddRAD-seq based phylogeographic study of Sargassum thunbergii (Phaeophyceae, Heterokonta) around Japanese coast

Using genome-wide SNP data obtained from high-throughput techniques based on double digest restriction site-associated DNA sequencing (ddRAD-seq), we elucidated the migration history and genetic diversity of the Japanese population of the ecologically important brown seaweed Sargassum thunbergii (Mertens ex Roth) Kuntze. STRUCTURE and NeighborNet analyses showed a clear genetic differentiation among populations of four geographic regions: Kyushu (POP1); Sea of Japan (POP2); Hokkaido and Tohoku (POP3); and Pacific coast from Kyushu to Kanto (POP4). Approximate Bayesian Computation (ABC) analysis indicated that POP4 diverged first, followed by the separation between POP2 (the largest effective population size) and POP3; POP1 was the last to form, shaped by the mixture of POP2 (73%) and POP4 (27%). High genetic diversity was detected in POP1 and POP2, whereas low genetic diversity was detected in POP3 and POP4. These results indicated that S. thunbergii populations of Kyushu and the Sea of Japan might have been maintained as large and stable populations gathered different lineages from China, Korea and/or Japan.

Molecular systematics reveals multiple lineages and cryptic speciation in the freshwater crayfish Parastacus brasiliensis (von Martens, 1869) (Crustacea : Decapoda : Parastacidae)

The characterisation of intraspecific genetic diversity in representatives of the South American crayfish genus Parastacus Huxley, 1879 is here carried out for the first time by comparing populations of Parastacus brasiliensis (von Martens, 1869) as currently defined. Phylogenetic reconstructions based on mitochondrial and nuclear markers indicate the existence of multiple lineages, of which only one can be considered as P. brasiliensis sensu stricto. In addition, there are seven other lineages, one of which is the subspecies Parastacus brasiliensis promatensis Fontoura & Conter, 2008, which is here elevated to species level. We thereby increase to 14 the number of recognised species within Parastacus in South America. Genetic distances among P. brasiliensis sensu stricto and the lineages ‘A’, ‘B’ and ‘C’ increase with geographical distances, suggesting isolation by distance as an important driver of diversification, and eventually speciation, in these burrowing crayfishes. Parastacus brasiliensis occurs mainly in the Guaíba Lake basin and studied populations show limited connectivity and gene flow, probably due to habitat fragmentation. On the basis of genetic and distribution data, the conservation status of P. brasiliensis sensu stricto is confirmed as Near Threatened (NT). These findings lead us to encourage the establishment of preservation areas for isolated populations. The importance of preserving the newly recognised distinct gene pools in order to maintain overall genetic diversity is emphasised.

Unrecognized diversity in New Guinean crayfish species (Decapoda, Parastacidae): The evidence from molecular data

The phylogenetic relationships among imported ornamental crayfish belonging to the genus Cherax were inferred from a combined dataset of 3 mitochondrial genes (COI, 16S and 12S) and by comparison with available GenBank sequences of 14 Cherax species. Furthermore, the concordance of previously described species obtained from a wholesaler (Cherax boesemani, C. holthuisi and C. peknyi) with available GenBank sequences was verified based on COI with special respect to comparison with sequences assigned as Cherax species. Recently described species C. gherardiae, C. pulcher and C. subterigneus belong to the northern group of Cherax species. Comparison and analysis with other GenBank COI sequences show previously unreported diversity of New Guinean species, suggesting 5 putative new species. Surprisingly, species assigned to the subgenus Astaconephrops do not form a monophyletic clade; this subgenus should be reappraised relative to the purported typical morphological characteristic of the uncalcified patch on male chelae. Increasing importation of crayfish underscores the importance of accurate species identification. Use of basic molecular methods is a necessary requisite for documenting occurrence, abundance and population trends of target species. Consequently, it helps to support eventual conservation decision-making by stakeholders.

Comparative riverscape genetics reveals reservoirs of genetic diversity for conservation and restoration of Great Plains fishes

We used comparative landscape genetics to examine the relative roles of historical events, intrinsic traits and landscape factors in determining the distribution of genetic diversity of river fishes across the North American Great Plains. Spatial patterns of diversity were overlaid on a patch-based graphical model and then compared within and among three species that co-occurred across five Great Plains watersheds. Species differing in reproductive strategy (benthic vs. pelagic-spawning) were hypothesized to have different patterns of genetic diversity, but the overriding factor shaping contemporary patterns of diversity was the signature of past climates and geological history. Allelic diversity was significantly higher at southern latitudes for Cyprinella lutrensis and Hybognathus placitus, consistent with northward expansion from southern Pleistocene refugia. Within the historical context, all species exhibited lowered occupancy and abundance in heavily fragmented and drier upstream reaches, particularly H. placitus; a pelagic-spawning species, suggesting rates of extirpation have outpaced losses of genetic diversity in this species. Within most tributary basins, genetically diverse populations of each species persisted. Hence, reconnecting genetically diverse populations with those characterized by reduced diversity (regardless of their position within the riverine network) would provide populations with greater genetic and demographic resilience. We discuss cases where cross-basin transfer may be appropriate to enhance genetic diversity and mitigate negative effects of climate change. Overall, striking similarities in genetic patterns and in response to fragmentation and dewatering suggest a common strategy for genetic resource management in this unique riverine fish assemblage.