Enigmatic incongruence between mtDNA and nDNA revealed by multi-locus phylogenomic analyses in freshwater snails

How complex is the Naineris setosa species complex? First integrative study of a presumed cosmopolitan and invasive annelid (Sedentaria: Orbiniidae)

We performed a comparative study of the specimens from the Naineris setosa complex from the Pacific and the Atlantic Oceans and re-described the syntype of N. setosa, including the selection of the lectotype. Molecular phylogenetic and species delimitation analyses based on two mitochondrial (COI and 16S) and one nuclear (28S) marker revealed the presence of three species. One clade with wide Amphi-Atlantic distribution was attributed as Naineris setosa s. str. The second Atlantic clade restricted to Southern and Southeastern Brazil was described as a new species, Naineris lanai sp. n. The third clade, reported from the Northwestern Pacific, was identified as a new species but was not formally described due to the presence of only juvenile-sized worms in the studied material. Detailed morphological descriptions of several diagnostic characters in the Naineris setosa complex are provided.

Host-Shift Speciation Proceeded with Gene Flow in Algae Covering Shells

AbstractHost shifts represent the advancement of a novel niche and often lead to speciation in symbionts. However, its mechanisms are not well understood. Here, we focused on the alga Pseudocladophora conchopheria growing on the shells of intertidal snails. Previous surveys have shown that the alga has host specificity-only attaching to the shell of Lunella correensis-but we discovered that the alga attaches to the shells of multiple sympatric snails. A genome-wide single-nucleotide polymorphism analysis (MIG-seq) was performed to determine whether host-associated speciation occurred in the algae. As a result, there was no gene flow or limited gene flow among the algae from different hosts, and some algae were genetically differentiated among hosts. In addition, the demographic estimate revealed that speciation with gene flow occurred between the algae from different hosts. Therefore, these results support the idea that host-shift speciation gradually proceeded with gene flow in the algae, providing insight into the early evolution of host shifts.

Phylogenomic analyses reveal incongruences between divergence times and fossil records of freshwater snails in East Asia

Fossils provide important insight into our understanding of phylogenetic history by serving as calibration points for divergence time estimation. However, uncertainties in the fossil record due to parallel evolution and convergent evolution can critically affect estimates of node ages. Here, we compare and contrast estimates of phylogenetic divergence with geologic and fossil history for two freshwater snail genera of the family Viviparidae in East Asia (Cipangopaludina and Margarya). Cipangopaludina species are commonly widely distributed species in East Asia, but extant Margarya species are endemic to the ancient lakes in Yunnan, China. According to some previous studies, parallel evolution or convergent evolution of shell morphology has occurred in the family several times which may affect divergence time estimation using fossil records. In this study, we used SNP data derived from ddRAD-seq loci to investigate population demographic history of both genera. Our results show a common pattern of lake endemic lineages diversifying from widely distributed lineages in the Miocene, and multiple colonization to a single ancient lake occurred in the Pleistocene. Our results indicate substantial incongruence among estimated phylogenomic divergence times, some fossil records, and formation ages of ancient lakes. These findings suggest some fossil records may be misidentified in these groups and highlight the need to carefully evaluate geological evidence and fossil records when using these for divergence time estimation.

Taxonomic insights and evolutionary history in East Asian terrestrial slugs of the genus Meghimatium

East Asia, specifically the Japanese Archipelago, is a biodiversity hotspot of both vertebrates and invertebrates. Mollusks represent a burst of species diversity in this region due to the effects of biotic and abiotic factors on their morphological traits, such as shell shape and size. However, the evolutionary history of terrestrial slugs in East Asia remains unknown. In the present study, we investigated the molecular phylogeny of terrestrial slugs of the genus Meghimatium. This genus includes three described and eight undescribed species, and our study used all except for two. Based on phylogeny and the species delimitation tests, the genus Meghimatium was split into many putative species, suggesting higher species diversity than previously thought based on morphological and anatomical studies and that almost undescribed species may be inappropriate. Therefore, morphological traits, such as body size and colour, conventionally considered for classification may easily vary or be similar across geographic region. Moreover, the divergence time of this genus is almost concordant with the geographical time scale of the formation of the Japanese mainland. Our findings suggest that molecular phylogenetics helps classify Japanese Meghimatium slugs, but comprehensive taxonomic revisions using multi-locus analyses are needed.

Multiple introgressions shape mitochondrial evolutionary history in Drosophila paulistorum and the Drosophila willistoni group

Hybridization and the consequent introgression of genomic elements is an important source of genetic diversity for biological lineages. This is particularly evident in young clades in which hybrid incompatibilities are still incomplete and mixing between species is more likely to occur. Drosophila paulistorum, a representative of the Neotropical Drosophila willistoni subgroup, is a classic model of incipient speciation. The species is divided into six semispecies that show varying degrees of pre- and post-mating incompatibility with each other. In the present study, we investigate the mitochondrial evolutionary history of D. paulistorum and the willistoni subgroup. For that, we perform phylogenetic and comparative analyses of the complete mitochondrial genomes and draft nuclear assemblies of 25 Drosophila lines of the willistoni and saltans species groups. Our results show that the mitochondria of D. paulistorum are polyphyletic and form two non-sister clades that we name α and β. Identification and analyses of nuclear mitochondrial insertions further reveal that the willistoni subgroup has an α-like mitochondrial ancestor and strongly suggest that both the α and β mitochondria of D. paulistorum were acquired through introgression from unknown fly lineages of the willistoni subgroup. We also uncover multiple mitochondrial introgressions across D. paulistorum semispecies and generate novel insight into the evolution of the species.

Deep structure, long-distance migration and admixture in the colour polymorphic land snail Cepaea nemoralis

Although snails of the genus Cepaea have historically been important in studying colour polymorphism, an ongoing issue is that there is a lack of knowledge of the underlying genetics of the polymorphism, as well as an absence of genomic data to put findings in context. We, therefore, used phylogenomic methods to begin to investigate the post‐glacial history of Cepaea nemoralis, with a long‐term aim to understand the roles that selection and drift have in determining both European‐wide and local patterns of colour polymorphism. By combining prior and new mitochondrial DNA data from over 1500 individuals with ddRAD genomic data from representative individuals across Europe, we show that patterns of differentiation are primarily due to multiple deeply diverged populations of snails. Minimally, there is a widespread Central European population and additional diverged groups in Northern Spain, the Pyrenees, as well as likely Italy and South Eastern Europe. The genomic analysis showed that the present‐day snails in Ireland and possibly some other locations are likely descendants of admixture between snails from the Pyrenees and the Central European group, an observation that is consistent with prior inferences from mitochondrial DNA alone. The interpretation is that C. nemoralis may have arrived in Ireland via long‐distance migration from the Pyrenean region, subsequently admixing with arrivals from elsewhere. This work, therefore, provides a baseline expectation for future studies on the genetics of the colour polymorphism, as well as providing a comparator for similar species.

Resolving species-level diversity of Beringiana and Sinanodonta mussels (Bivalvia: Unionidae) in the Japanese archipelago using genome-wide data

Accurate species identification is of primary importance in ecology and evolutionary biology. For a long time, the unionid mussels Beringiana and Sinanodonta have puzzled researchers trying to unravel their diversity because of their poorly discernible morphologies. A recent study conducted species delineation of unionid mussels based on mitochondrial DNA variation, opening up a new avenue to grasp species diversity of the mussels. However, mtDNA-based classification may not align with species boundaries because mtDNA is prone to introgression and incomplete lineage sorting that cause discordance between species affiliation and gene phylogeny. In this study, we evaluated the validity of the mtDNA-based classification of unionid mussels Beringiana and Sinanodonta in Japan using mitochondrial sequence data, double digest restriction site-associated DNA library (ddRAD) sequencing, and morphological data. We found significant inconsistencies in the mitochondrial and nuclear DNA phylogenies, casting doubt on the reliability of the mtDNA-based classification in this group. In addition, nuclear DNA phylogeny revealed that there are at least two unionid lineages hidden in the mtDNA phylogeny. Although molecular dating technique indicates that Beringiana and Sinanodonta diverged >35 million years ago, their shell morphologies are often indistinguishable. Specifically, morphological analyses exhibited the parallel appearance of nearly identical ball-like shell forms in the two genera in Lake Biwa, which further complicates species identification and the morphological evolution of unionid mussels. Our study adds to a growing body of literature that accurate species identification of unionid mussels is difficult when using morphological characters alone. Although mtDNA-based classification is a simple and convenient way to classify unionid mussels, considerable caution is warranted for its application in ecological and evolutionary studies.

Species identification of introduced veronicellid slugs in Japan

Reliable identification of species is important for protecting native ecosystems against the invasion of non-native species. DNA barcoding using molecular markers, such as the mitochondrial cytochrome oxidase subunit 1 (COI) gene, helps researchers distinguish species. In this study, we focused on introduced veronicellid slugs in the Ryukyu Islands and some greenhouses on mainland Japan. Some veronicellids are medium-to-high risk pest species for humans. Identifying veronicellid species by their external morphology is difficult and unreliable because there is substantial overlap between intraspecific variation and interspecific differentiation. Therefore, internal morphologies such as male genitalia have been the primary traits used to distinguish veronicellids. To identify introduced veronicellid slugs in Japan to the species level, we used morphological assessment of male genitalia and DNA barcoding of the standard COI gene fragment. We also conducted species-delimitation analyses based on the genetic data. The results showed that five evolutionarily significant units, corresponding to four nominal species inhabit the Ryukyu Islands, of which two species were also found in the greenhouses of mainland Japan, including the first record of Sarasinula plebeia in Japan. The presence of non-native slug species could increase the transmission of parasites in Japan.

Integrated Taxonomy Revealed Genetic Differences in Morphologically Similar and Non-Sympatric Scoliodon macrorhynchos and S. laticaudus

Simple Summary

In this study, the species identities of similar-looking coastal spadenose sharks from different areas were clarified by adding new molecular markers and more individual body measurements, including animals from the Malaysian Peninsula that had not been examined previously. Collective evidence showed that there are two genetically distinct species that do not overlap in their spatial occurrence. The Malacca Strait acts as a boundary delineating the distribution range of the Pacific spadenose shark Scoliodon macrorhynchos to the east and, of the Northern Indian Ocean, S. laticaudus to the west. In addition, the need to determine the species status of Scoliodon animals from Indonesian waters was identified. The present study reinforced the need to rely on comprehensive genetic information in addition to external characteristics to assess the species identities and distribution range for small sharks and rays that have apparent contiguous coastal distribution with limited dispersal abilities.

Abstract

Previous examination of the mitochondrial NADH2 gene and morphological characteristics led to the resurrection of Scoliodon macrorhynchos as a second valid species in the genus, in addition to S. laticaudus. This study applied an integrated taxonomic approach to revisit the classification of the genus Scoliodon based on new materials from the Malaysian Peninsula, Malaysian Borneo and Eastern Bay of Bengal. Mitochondrial DNA data suggested the possibility of three species of Scoliodon in the Indo-West Pacific, while the nuclear DNA data showed partially concordant results with a monophyletic clade of S. macrorhynchos and paraphyletic clades of S. laticaudus and S. cf. laticaudus from the Malacca Strait. Morphological, meristic and dental characteristics overlapped between the three putative species. Collective molecular and morphological evidence suggested that the differences that exist among the non-sympatric species of Scoliodon are consistent with isolation by distance, and Scoliodon macrorhynchos remains as a valid species, while S. cf. laticaudus is assigned as S. laticaudus. The Malacca Strait acts as a spatial delineator in separating the Pacific S. macrorhynchos (including South China Sea) from the Northern Indian Ocean S. laticaudus. Future taxonomic work should focus on clarifying the taxonomic status of Scoliodon from the Indonesian waters.

High-throughput SNPs dataset reveal restricted population connectivity of marine gastropod within the narrow distribution range of peripheral oceanic islands

Molecular studies based on the high resolution genetic markers help us to grasp the factor shaping the genetic structure of marine organisms. Ecological factors linking to life history traits have often explained the process of genetic structuring in open and connectable oceanic environments. Besides, population genetic divergence can be affected by fragmented habitat, oceanic current, and past geographical events. In the present study, we demonstrated the genetic differentiation of marine gastropod Monodonta sp. within a narrow range of peripheral oceanic islands, the Ogasawara Islands. Genetic analyses were performed not only with a mitochondrial DNA marker but also with a high-throughput SNPs dataset obtained by ddRAD-seq. The results of the mtDNA analyses did not show genetic divergence among populations, while the SNPs dataset detected population genetic differentiation. Population demographic analyses and gene flow estimation suggested that the genetic structure was formed by sea level fluctuation associated with the past climatic change and regulated by temporal oceanographic conditions. These findings provide important insights into population genetic patterns in open and connectable environments.

Chronological Incongruences between Mitochondrial and Nuclear Phylogenies of Aedes Mosquitoes

One-third of all mosquitoes belong to the Aedini, a tribe comprising common vectors of viral zoonoses such as Aedes aegypti and Aedes albopictus. To improve our understanding of their evolution, we present an updated multigene estimate of Aedini phylogeny and divergence, focusing on the disentanglement between nuclear and mitochondrial phylogenetic signals. We first show that there are some phylogenetic discrepancies between nuclear and mitochondrial markers which may be caused by wrong taxa assignment in samples collections or by some stochastic effect due to small gene samples. We indeed show that the concatenated dataset is model and framework dependent, indicating a general paucity of signal. Our Bayesian calibrated divergence estimates point toward a mosquito radiation in the mid-Jurassic and an Aedes radiation from the mid-Cretaceous on. We observe, however a strong chronological incongruence between mitochondrial and nuclear data, the latter providing divergence times within the Aedini significantly younger than the former. We show that this incongruence is consistent over different datasets and taxon sampling and that may be explained by either peculiar evolutionary event such as different levels of saturation in certain lineages or a past history of hybridization throughout the genus. Overall, our updated picture of Aedini phylogeny, reveal a strong nuclear-mitochondrial incongruence which may be of help in setting the research agenda for future phylogenomic studies of Aedini mosquitoes.

Refugia within refugium of Geranium yesoense (Geraniaceae) in Japan were driven by recolonization into the southern interglacial refugium

Recent studies have found that geographical fragmentation and recurrent colonization result in complex genetic structures in refugial areas. This phenomenon, known as ‘refugia within refugium’, has been identified from many geographical locations. In Japan, the high-elevation mountains of central Honshu provided an interglacial refugium for alpine plants. Here we focused on the Geranium yesoense complex, which exhibits increased morphological variation in the refugial area, to determine whether this variation was shaped by recurrent colonization, range fragmentation or phenotypic changes independent of population history. We analysed single nucleotide polymorphism data and chloroplast genome sequences. Diversification in the G. yesoense species complex occurred in the mid-Pleistocene. The varieties are distinct entities and suggest the presence of a genetic cluster with highly disjunct distributions, occurring both in northern Japan and in southern refugial areas in central Honshu. Demographic analysis suggests that a single ancestral variety (var. nipponicum) evolved in the alpine region of central Honshu, and that subsequent migration from one of the two diverged northern varieties (var. pseudopratense) led to secondary contact with var. nipponicum during the last glacial period. Recolonization into refugial populations in central Honshu and hybridization between diverged populations have resulted in complex genetic structures among refugial populations.

Global Diversification Dynamics Since the Jurassic: Low Dispersal and Habitat-Dependent Evolution Explain Hotspots of Diversity and Shell Disparity in River Snails (Viviparidae)

The Viviparidae, commonly known as River Snails, is a dominant group of freshwater snails with a nearly worldwide distribution that reaches its highest taxonomic and morphological diversity in Southeast Asia. The rich fossil record is indicative of a probable Middle Jurassic origin on the Laurasian supercontinent where the group started to diversify during the Cretaceous. However, it remains uncertain when and how the biodiversity hotspot in Southeast Asia was formed. Here, we used a comprehensive genetic data set containing both mitochondrial and nuclear markers and comprising species representing 24 out of 28 genera from throughout the range of the family. To reconstruct the spatiotemporal evolution of viviparids on a global scale, we reconstructed a fossil-calibrated phylogeny. We further assessed the roles of cladogenetic and anagenetic events in range evolution. Finally, we reconstructed the evolution of shell features by estimating ancestral character states to assess whether the appearance of sculptured shell morphologies was driven by major habitat shifts. The molecular phylogeny supports the monophyly of the three subfamilies, the Bellamyinae, Lioplacinae, and Viviparinae, but challenges the currently accepted genus-level classification in several cases. The almost global distribution of River Snails has been influenced both by comparatively ancient vicariance and more recent founder events. In Southeast Asia, Miocene dispersal was a main factor in shaping the modern species distributions. A recurrent theme across different viviparid taxa is that many species living in lentic waters exhibit sculptured shells, whereas only one strongly sculptured species is known from lotic environments. We show that such shell sculpture is habitat-dependent and indeed evolved several times independently in lentic River Snails. Considerably high transition rates between shell types in lentic habitats probably caused the co-occurrence of morphologically distinct shell types in several lakes. In contrast, directional evolution toward smooth shells in lotic habitats, as identified in the present analyses, explains why sculptured shells are rarely found in these habitats. However, the specific factors that promoted changes in shell morphology require further work. [biogeographical analyses; fossil-calibrated phylogeny; fossil-constrained analyses; Southeast Asia; stochastic character mapping.].

New insights from museum specimens: a case of Viviparidae (Caenogastropoda: Mollusca) in Iwakawa's collection preserved in the National Museum of Nature and Science, Tokyo

Background

In this study, we clarify the classification of museum specimens of the family Viviparidae, which is composed of six species/subspecies in Japan, including three endangered species. We examined Viviparus sclateri specimens from the Tomotaro Iwakawa collection (1855-1933) in the National Museum of Nature and Science, Tokyo. The Iwakawa's collection was catalogued in 1919 and Viviparus sclateri, labelled with Naga-tanishi, the current Japanese name for H. longispira, which was, at the time, equivalent to Viviparus sclateri, was listed in this catalogue. The catalogue noted localities of Viviparus sclateri (Naga-tanishi) from outside Lake Biwa, including occurrences in Lake Kasumigaura and Lake Suwa. However, Heterogen longispira (Naga-tanishi) is currently considered to be endemic to Lake Biwa drainage. The actual status of Viviparus sclateri in Iwakawa (1919) has not been clarified until now.

New information

Our examination revealed that Viviparus sclateri from Iwakawa's catalogue included H. japonica, H. longispira and Sinotaia quadrata histrica, based on current taxonomy. Specimens assigned to H. longispira occurred only in Lake Biwa drainage. Heterogen japonica was confirmed to be present in all lots and some H. japonica from Lake Suwa had a distinctive morphology. Sinotaia quadrata histrica was only confirmed to occur in Lake Suwa. Furthermore, some specimens from southern Lake Biwa and the Seta River had intermediate characteristics between H. japonica and H. longispira and their populations are currently almost extinct.

Human-geographic effects on variations in the population genetics of Sinotaia quadrata (Gastropoda: Viviparidae) that historically migrated from continental East Asia to Japan

BACKGROUND

Anthropogenic factors potentially affect observed biogeographical patterns in population genetics, but the effects of ancient human activities on the original patterns created by natural processes are unknown. Sinotaia quadrata, a widely distributed freshwater snail species in East Asia, was used to investigate this issue. It is unclear whether S. quadrata in Japan was introduced from China and how different human uses and varying geographic patterns affect the contemporary population genetics between the two regions. Thus, we investigated the demography of S. quadrata and detected its genetic structure in Japan and continental East Asia.

RESULTS

Sinotaia quadrata populations first naturally migrated from continental East Asia to Japan, which is associated with the ancient period in Japanese geohistory (about 70,000 years ago). They were then artificially introduced in association with agriculture expansion by human movements in two recent periods (about 8,000 and 1,200 years ago). Populations in different parts of Japan have their own sources. Natural migration in the ancient period and artificial introduction in the recent period suggest that the population distribution is affected by both the geohistory of East Asia and the history of human expansion. In the background of the historical migration and introduction, contemporary populations in the two regions show different genetic patterns. Population divergence levels were significantly correlated with geographical patterns in Japan and significantly correlated with human interventions variables in continental East Asia, suggesting that long-term geographical isolation is likely the major factor that shaped the contemporary population genetics in Japan, while modern human uses are likely the major factor in continental East Asia.

CONCLUSIONS

Our preliminary results show a complex demography and unusual genetic patterns in the contemporary populations for a common freshwater snail and are of significance to determine the historical formation and contemporary patterns of biogeography in Japan and continental East Asia.

The efficiency of universal mitochondrial DNA barcodes for species discrimination of Pomacea canaliculata and Pomacea maculata

Invasive apple snails, Pomacea canaliculata and P. maculata, have a widespread distribution globally and are regarded as devastating pests of agricultural wetlands. The two species are morphologically similar, which hinders species identification via morphological approaches and species-specific management efforts. Advances in molecular genetics may contribute effective diagnostic tools to potentially resolve morphological ambiguity. DNA barcoding has revolutionized the field of taxonomy by providing an alternative, simple approach for species discrimination, where short sections of DNA, the cytochrome c oxidase subunit I (COI) gene in particular, are used as ‘barcodes’ to delineate species boundaries. In our study, we aimed to assess the effectiveness of two mitochondrial markers, the COI and 16S ribosomal deoxyribonucleic acid (16S rDNA) markers for DNA barcoding of P. canaliculata and P. maculata. The COI and 16S rDNA sequences of 40 Pomacea specimens collected from six localities in Peninsular Malaysia were analyzed to assess their barcoding performance using phylogenetic methods and distance-based assessments. The results confirmed both markers were suitable for barcoding P. canaliculata and P. maculata. The phylogenies of the COI and 16S rDNA markers demonstrated species-specific monophyly and were largely congruent with the exception of one individual. The COI marker exhibited a larger barcoding gap (6.06–6.58%) than the 16S rDNA marker (1.54%); however, the magnitude of barcoding gap generated within the barcoding region of the 16S rDNA marker (12-fold) was bigger than the COI counterpart (approximately 9-fold). Both markers were generally successful in identifying P. canaliculata and P. maculata in the similarity-based DNA identifications. The COI + 16S rDNA concatenated dataset successfully recovered monophylies of P. canaliculata and P. maculata but concatenation did not improve individual datasets in distance-based analyses. Overall, although both markers were successful for the identification of apple snails, the COI molecular marker is a better barcoding marker and could be utilized in various population genetic studies of P. canaliculata and P. maculata.