A cyclical switch of gametogenic pathways in hybrids depends on the ploidy level.

The cellular and molecular mechanisms governing sexual reproduction are conserved across eukaryotes. Nevertheless, hybridization can disrupt these mechanisms, leading to asexual reproduction, often accompanied by polyploidy. In this study, we investigate how ploidy level and ratio of parental genomes in hybrids affect their reproductive mode. We analyze the gametogenesis of sexual species and their diploid and triploid hybrids from the freshwater fish family Cobitidae, using newly developed cytogenetic markers. We find that diploid hybrid females possess oogonia and oocytes with original (diploid) and duplicated (tetraploid) ploidy. Diploid oocytes cannot progress beyond pachytene due to aberrant pairing. However, tetraploid oocytes, which emerge after premeiotic genome endoreplication, exhibit normal pairing and result in diploid gametes. Triploid hybrid females possess diploid, triploid, and haploid oogonia and oocytes. Triploid and haploid oocytes cannot progress beyond pachytene checkpoint due to aberrant chromosome pairing, while diploid oocytes have normal pairing in meiosis, resulting in haploid gametes. Diploid oocytes emerge after premeiotic elimination of a single-copied genome. Triploid hybrid males are sterile due to aberrant pairing and the failure of chromosomal segregation during meiotic divisions. Thus, changes in ploidy and genome dosage may lead to cyclical alteration of gametogenic pathways in hybrids.

Genome assembly of the Korean intertidal mud-creeper Batillaria attramentaria.

Batillaridae is a common gastropod family that occurs abundantly in the shallow coastal zone of the intertidal mudflats of the northwest Pacific Ocean, Australasia, and North America. In this family, Batillaria attramentaria is known for its biological invasion and colonization in estuarine and intertidal zones. It can endure and adapt the harsh intertidal conditions such as frequent temperature alteration, salinity, and air exposure. Therefore, we sequenced and assembled this Korean batillariid genome to get insight into its intertidal adaptive features. Approximately 53 Gb of DNA sequences were generated, and 863 scaffolds were assembled into a draft genome of 0.715 Gb with 97.1% BUSCO completeness value. A total of 40,596 genes were predicted. We estimated that B. attramentaria and Conus consors diverged about 230 million years ago (MYA) based on the phylogenetic analysis of closely related gastropod species. This genome study sets the footstep for genomics studies among native and introduced Batillaria populations and the Batillaridae family members.

A regulatory hydrogenase gene cluster observed in the thioautotrophic symbiont of Bathymodiolus mussel in the East Pacific Rise.

The mytilid mussel Bathymodiolus thermophilus lives in the deep-sea hydrothermal vent regions due to its relationship with chemosynthetic symbiotic bacteria. It is well established that symbionts reside in the gill bacteriocytes of the mussel and can utilize hydrogen sulfide, methane, and hydrogen from the surrounding environment. However, it is observed that some mussel symbionts either possess or lack genes for hydrogen metabolism within the single-ribotype population and host mussel species level. Here, we found a hydrogenase cluster consisting of additional H2-sensing hydrogenase subunits in a complete genome of B. thermophilus symbiont sampled from an individual mussel from the East Pacific Rise (EPR9N). Also, we found methylated regions sparsely distributed throughout the EPR9N genome, mainly in the transposase regions and densely present in the rRNA gene regions. CRISPR diversity analysis confirmed that this genome originated from a single symbiont strain. Furthermore, from the comparative analysis, we observed variation in genome size, gene content, and genome re-arrangements across individual hosts suggesting multiple symbiont strains can associate with B. thermophilus. The ability to acquire locally adaptive various symbiotic strains may serve as an effective mechanism for successfully colonizing different chemosynthetic environments across the global oceans by host mussels.

Divergent paths in the evolutionary history of maternally transmitted clam symbionts.

Vertical transmission of bacterial endosymbionts is accompanied by virtually irreversible gene loss that results in a progressive reduction in genome size. While the evolutionary processes of genome reduction have been well described in some terrestrial symbioses, they are less understood in marine systems where vertical transmission is rarely observed. The association between deep-sea vesicomyid clams and chemosynthetic Gammaproteobacteria is one example of maternally inherited symbioses in the ocean. Here, we assessed the contributions of drift, recombination and selection to genome evolution in two extant vesicomyid symbiont clades by comparing 15 representative symbiont genomes (1.017-1.586 Mb) to those of closely related bacteria and the hosts' mitochondria. Our analyses suggest that drift is a significant force driving genome evolution in vesicomyid symbionts, though selection and interspecific recombination appear to be critical for maintaining symbiont functional integrity and creating divergent patterns of gene conservation. Notably, the two symbiont clades possess putative functional differences in sulfide physiology, anaerobic respiration and dependency on environmental vitamin B12, which probably reflect adaptations to different ecological habitats available to each symbiont group. Overall, these results contribute to our understanding of the eco-evolutionary processes shaping reductive genome evolution in vertically transmitted symbioses.

Experimental data supporting adaptive locomotor responses to salt stress in the mud-tidal gastropod populations (Batillaria).

This article describes the experimental locomotor data used to study the general and adaptive responses to salt stress of the northern Pacific intertidal gastropod Batillaria attramentaria. The data were obtained from a series of 30-day experiments on snails acclimated to different salinity regimes. Snails were collected from coastal areas on the eastern and western sides of the North Pacific Ocean. The data consist of three parts: 1) raw videos recording the locomotion of the snails when exposed to novel artificial salinity regimes in laboratory settings, 2) Spectral Time-Lapse results of movement distance of the snails extracted from the recorded videos, and 3) CO1-gene sequences isolated from individuals collected from four sampling sites. A Linear Mixed-effect Model inference procedure was applied in an attempt to assess the impacts of geographic distribution and genetic composition on the locomotor response to salt stress in the snail B. attramentaria. The locomotor dataset we present are the first reports of locomotor response to salt stress of the snail B. attramentaria, that is valuable for further exploration and understanding of the impacts of environmental changes on the physiology and adaptive capacity of living marine molluscs.

Locomotor responses to salt stress in native and invasive mud-tidal gastropod populations (Batillaria).

Plasticity in salt tolerance can be crucial for successful biological invasions of novel habitats by marine gastropods. The intertidal snail Batillaria attramentaria, which is native to East Asia but invaded the western shores of North America from Japan 80 years ago, provides an opportunity to examine how environmental salinity may shape behavioral and morphological traits. In this study, we compared the movement distance of four B. attramentaria populations from native (Korea and Japan) and introduced (United States) habitats under various salinity levels (13, 23, 33, and 43 PSU) during 30 days of exposure in the lab. We sequenced a partial mitochondrial CO1 gene to infer phylogenetic relationships among populations and confirmed two divergent mitochondrial lineages constituting our sample sets. Using a statistical model-selection approach, we investigated the effects of geographic distribution and genetic composition on locomotor performance in response to salt stress. Snails exposed to acute low salinity (13 PSU) reduced their locomotion and were unable to perform at their normal level (the moving pace of snails exposed to 33 PSU). We did not detect any meaningful differences in locomotor response to salt stress between the two genetic lineages or between the native snails (Japan vs. Korea populations), but we found significant locomotor differences between the native and introduced groups (Japan or Korea vs. the United States). We suggest that the greater magnitude of tidal salinity fluctuation at the US location may have influenced locomotor responses to salt stress in introduced snails.

Dataset supporting description of the new mussel species of genus Gigantidas (Bivalvia: Mytilidae) and metagenomic data of bacterial community in the host mussel gill tissue.

This article contains supplementary data from the research paper entitled "A newly discovered Gigantidas bivalve mussel from the Onnuri Vent Field on the northern Central Indian Ridge" [1], describes a new mussel species within the subfamily Bathymodiolinae named Gigantidas vrijenhoeki. Data are comprised of two parts: 1) shell image and molecular analyses of G. vrijenhoeki and 2) metagenomic community analyses of gill-associated symbiotic bacteria on G. vrijenhoeki. G. vrijenhoeki data were obtained from type specimens described in Jang et al. 2020 [1]. The molecular analysis was conducted by calculating genetic distance at intra- and inter-specific level within genus Gigantidas based on the sequence data of two mitochondrial genes (COI and ND4). The metagenomic dataset of gill-associated symbionts were generated by Illumina Miseq sequencing of the V3-V4 region of 16S rRNA from 12 specimens of G. vrijenhoeki collected from the same vent site, Onnuri Vent Field.

Impacts of Salt Stress on Locomotor and Transcriptomic Responses in the Intertidal Gastropod Batillaria attramentaria.

Salinity is one of the most crucial environmental factors that structures biogeographic boundaries of aquatic organisms, affecting distribution, abundance, and behavior. However, the association between behavior and gene regulation underlying acclimation to changes in salinity remains poorly understood. In this study, we investigated the effects of salinity stress on behavior (movement distance) and patterns of gene expression (using RNA sequencing) of the intertidal gastropod Batillaria attramentaria. We examined responses to short-term (1-hour) and long-term (30-day) acclimation to a range of salinities (43, 33 [control], 23, 13, and 3 psu). We found that the intertidal B. attramentaria is able to tolerate a broad range of salinity from 13 to 43 psu but not the acute low salinity of 3 psu. Behavioral experiments showed that salt stress significantly influenced snails' movement, with lower salinity resulting in shorter movement distance. Transcriptomic analyses revealed critical metabolic pathways and genes potentially involved in acclimation to salinity stress, including ionic and osmotic regulation, signal and hormonal transduction pathways, water exchange, cell protection, and gene regulation or epigenetic modification. In general, our study presents a robust, integrative laboratory-based approach to investigate the effects of salt stress on a nonmodel gastropod facing detrimental consequences of environmental change. The current genetic results provide a wealth of reference data for further research on mechanisms of ionic and osmotic regulation and adaptive evolution of this coastal gastropod.

Bidirectional mitochondrial introgression between Korean cobitid fish mediated by hybridogenetic hybrids.

Genomic introgression through interspecific hybridization has been observed in some species of the freshwater fish family Cobitidae. Within this family, a Cobitis hankugensis-Iksookimia longicorpa diploid-triploid hybrid species complex on the Korean peninsula is unique in displaying hybridogenesis, a unisexual reproduction mode that allows hybrids to mediate the transfer of mitochondrial DNA (but not nuclear DNA) between the two parent species. However, populations of the parental species in the wild have never been examined for the potential effect of introgression on their genomes. To address the genetic consequences of unisexual hybridization on the parental species, we examined genetic structure of the two parental species, C. hankugensis and I. longicorpa, in three independent natural habitats where they coexist with their hybrid complex using DNA sequence data of one mitochondrial gene and three nuclear genes. We found that mitochondrial introgression between the two species was extensive in all the examined localities, while there was no evidence of nuclear introgression across the species boundary. This result indicates that the hybridogenetic individuals mediate mitochondrial introgression from one species to the other, producing mito-nuclear mosaic genomes such as C. hankugensis nuclear genomes associated with I. longicorpa mitochondrial DNA and the reverse. The direction and degree of introgression varied among the three localities, but the underlying mechanisms for this observation proved elusive. Introgression might depend on which species serves as the predominant sperm or ovum donor or the environmental conditions of the localities. The present study suggests that introgressive hybridization between pure C. hankugensis and I. longicorpa species is highly likely where the two species co-occur with hybridogenetic individuals, but the consequence of introgression could be variable due to the history and environmental characteristics of particular populations across the parental species' ranges.

A new yeti crab phylogeny: Vent origins with indications of regional extinction in the East Pacific.

The recent discovery of two new species of kiwaid squat lobsters on hydrothermal vents in the Pacific Ocean and in the Pacific sector of the Southern Ocean has prompted a re-analysis of Kiwaid biogeographical history. Using a larger alignment with more fossil calibrated nodes than previously, we consider the precise relationship between Kiwaidae, Chirostylidae and Eumunididae within Chirostyloidea (Decapoda: Anomura) to be still unresolved at present. Additionally, the placement of both new species within a new "Bristly" clade along with the seep-associated Kiwa puravida is most parsimoniously interpreted as supporting a vent origin for the family, rather than a seep-to-vent progression. Fossil-calibrated divergence analysis indicates an origin for the clade around the Eocene-Oligocene boundary in the eastern Pacific ~33-38 Ma, coincident with a lowering of bottom temperatures and increased ventilation in the Pacific deep sea. Likewise, the mid-Miocene (~10-16 Ma) rapid radiation of the new Bristly clade also coincides with a similar cooling event in the tropical East Pacific. The distribution, diversity, tree topology and divergence timing of Kiwaidae in the East Pacific is most consistent with a pattern of extinctions, recolonisations and radiations along fast-spreading ridges in this region and may have been punctuated by large-scale fluctuations in deep-water ventilation and temperature during the Cenozoic; further affecting the viability of Kiwaidae populations along portions of mid-ocean ridge.

Multiple Modes of Positive Selection Shaping the Patterns of Incomplete Selective Sweeps over African Populations of Drosophila melanogaster.

It remains a challenge in evolutionary genetics to elucidate how beneficial mutations arise and propagate in a population and how selective pressures on mutant alleles are structured over space and time. By identifying "sweeping haplotypes (SHs)" that putatively carry beneficial alleles and are increasing (or have increased) rapidly in frequency, and surveying the geographic distribution of SH frequencies, we can indirectly infer how selective sweeps unfold in time and thus which modes of positive selection underlie those sweeps. Using population genomic data from African Drosophila melanogaster, we identified SHs from 37 candidate loci under selection. At more than half of loci, we identify single SHs. However, many other loci harbor multiple independent SHs, namely soft selective sweeps, either due to parallel evolution across space or a high beneficial mutation rate. At about a quarter of the loci, intermediate SH frequencies are found across multiple populations, which cannot be explained unless a certain form of frequency-dependent positive selection, such as heterozygote advantage, is invoked given the reasonable range of migration rates between African populations. At one locus, many independent SHs are observed over multiple populations but always together with ancestral haplotypes. This complex pattern is compatible with a large number of mutational targets in a gene and frequency-dependent selection on new variants. We conclude that very diverse modes of positive selection are operating at different sets of loci in D. melanogaster populations.

Geographical structure of endosymbiotic bacteria hosted by Bathymodiolus mussels at eastern Pacific hydrothermal vents.

BACKGROUND: Chemolithoautotrophic primary production sustains dense invertebrate communities at deep-sea hydrothermal vents and hydrocarbon seeps. Symbiotic bacteria that oxidize dissolved sulfur, methane, and hydrogen gases nourish bathymodiolin mussels that thrive in these environments worldwide. The mussel symbionts are newly acquired in each generation via infection by free-living forms. This study examined geographical subdivision of the thiotrophic endosymbionts hosted by Bathymodiolus mussels living along the eastern Pacific hydrothermal vents. High-throughput sequencing data of 16S ribosomal RNA encoding gene and fragments of six protein-coding genes of symbionts were examined in the samples collected from nine vent localities at the East Pacific Rise, Galápagos Rift, and Pacific-Antarctic Ridge. RESULTS: Both of the parapatric sister-species, B. thermophilus and B. antarcticus, hosted the same numerically dominant phylotype of thiotrophic Gammaproteobacteria. However, sequences from six protein-coding genes revealed highly divergent symbiont lineages living north and south of the Easter Microplate and hosted by these two Bathymodiolus mussel species. High heterogeneity of symbiont haplotypes among host individuals sampled from the same location suggested that stochasticity associated with initial infections was amplified as symbionts proliferated within the host individuals. The mussel species presently contact one another and hybridize along the Easter Microplate, but the northern and southern symbionts appear to be completely isolated. Vicariance associated with orogeny of the Easter Microplate region, 2.5-5.3 million years ago, may have initiated isolation of the symbiont and host populations. Estimates of synonymous substitution rates for the protein-coding bacterial genes examined in this study were 0.77-1.62%/nucleotide/million years. CONCLUSIONS: Our present study reports the most comprehensive population genetic analyses of the chemosynthetic endosymbiotic bacteria based on high-throughput genetic data and extensive geographical sampling to date, and demonstrates the role of the geographical features, the Easter Microplate and geographical distance, in the intraspecific divergence of this bacterial species along the mid-ocean ridge axes in the eastern Pacific. Altogether, our results provide insights into extrinsic and intrinsic factors affecting the dispersal and evolution of chemosynthetic symbiotic partners in the hydrothermal vents along the eastern Pacific Ocean.

Impact of human activities on changes of ichthyofauna in Dongjin River of Korea in the past 30 years.

Ichthyofauna and fish community were investigated at 17 representative stations of the Dongjin River drainage system from spring to fall in 2014. The survey resulted in a list of 53 species belonging to 14 families structured into 4 distinctive parts along the river: uppermost-stream, upper-stream, mid-stream, and lower-stream. Comparison of species lists with 30-year interval exhibited significant decreases in peripheral freshwater fishes, Acheilognathinae, endemic, and indigeneity species, but increases in exotic, epipelagic, and lentic species. Moreover, in the estuary of the Dongjin River drainage system, peripheral freshwater fish species were replaced by pure freshwater fish species due to the Saemangeum sea-wall project. In the upper region of the river, introduced eight alien species from Seomjin River via water diversion tunnels. In the mid-lower region, the construction of floodgates and numerous small weirs caused expansion of lentic water areas, facilitating the spread of problematic exotic species such as Micropterus salmoides, Lepomis macrochirus, and Carassius cuvieri. Also, water deterioration in this region resulted in an increase of tolerant species and a decrease of sensitive and endemic species. Our results suggest that a recovery strategy for a healthy ecosystem in the Dongjin River drainage system should reflect this compartmentalized cause and effect on the changes of icthyofauna.

Small fishes crossed a large mountain range: Quaternary stream capture events and freshwater fishes on both sides of the Taebaek Mountains.

The Taebaek Mountains in Korea serve as the most apparent biogeographic barrier for Korean freshwater fishes, resulting in 2 distinct ichthyofaunal assemblages on the eastern (East/Japan Sea slope) and western (Yellow Sea and Korea Strait slopes) sides of the mountain range. Of nearly 100 species of native primary freshwater fishes in Korea, only 18 species occur naturally on both sides of the mountain range. Interestingly, there are 5 rheophilic species (Phoxinus phoxinus, Coreoleuciscus splendidus, Ladislavia taczanowskii, Iksookimia koreensis and Koreocobitis rotundicaudata) found on both sides of the Taebaek Mountains that are geographically restricted to the Osip River (and several neighboring rivers, for L. taczanowskii and I. koreensis) on the eastern side of the mountain range. The Osip River and its neighboring rivers also shared a rheophilic freshwater fish, Liobagrus mediadiposalis, with the Nakdong River on the western side of the mountain range. We assessed historical biogeographic hypotheses on the presence of these rheophilic fishes, utilizing DNA sequence data from the mitochondrial cytochrome b gene. Results of our divergence time estimation indicate that ichthyofaunal transfers into the Osip River (and several neighboring rivers in East Sea slope) have occurred from the Han (Yellow Sea slope) and Nakdong (Korea Strait slope) Rivers since the Late Pleistocene. The inferred divergence times for the ichthyofaunal transfer across the Taebaek Mountains were consistent with the timing of hypothesized multiple reactivations of the Osip River Fault (Late Pleistocene), suggesting that the Osip River Fault reactivations may have caused stream capture events, followed by ichthyofaunal transfer, not only between the Osip and Nakdong Rivers, but also between the Osip and Han Rivers.

Cobitis nalbanti, a new species of spined loach from South Korea, and redescription of Cobitis lutheri (Teleostei: Cobitidae).

Cobitis nalbanti, new species, is described from the South Korean Han and Geum Rivers, draining to the Yellow Sea. It differs from its congeners by having a wide, ovoid lamina circularis; rounded scales with a large, slightly-displaced focal zone; a relatively long, protruded snout; the suborbital spine not reaching the center of the eye; the mandibular barbel not reaching the anterior edge of the eye; usually 12-14 dark brown blotches in the fourth Gambetta's zone; a single elongated black spot on the upper part of the caudal-fin base; the only streak on the head running from the tip of the snout to the nape, across the eye; 2n=48-51 and NF=66. Cobitis nalbanti was previously identified as C. lutheri. Morphological analysis of C. lutheri shows that it is distinct in having a short obtuse snout, the suborbital spine reaching the center of the eye, the mandibular barbel reaching the anterior edge of the eye, usually less than 12 dark brown blotches in the fourth Gambetta's zone, two spots at the caudal-fin base, several dark streaks on the head, and 2n=50, NF=70. A lectotype is selected for C. lutheri. Cobitis lutheri has a wide distribution in the Amur River basin, several rivers flowing into the Sea of Japan (East Sea), and in the northern part of Sakhalin Island.

Population subdivision of hydrothermal vent polychaete Alvinella pompejana across equatorial and Easter Microplate boundaries.

BACKGROUND: The Equator and Easter Microplate regions of the eastern Pacific Ocean exhibit geomorphological and hydrological features that create barriers to dispersal for a number of animals associated with deep-sea hydrothermal vent habitats. This study examined effects of these boundaries on geographical subdivision of the vent polychaete Alvinella pompejana. DNA sequences from one mitochondrial and eleven nuclear genes were examined in samples collected from ten vent localities that comprise the species' known range from 23°N latitude on the East Pacific Rise to 38°S latitude on the Pacific Antarctic Ridge. RESULTS: Multi-locus genotypes inferred from these sequences clustered the individual worms into three metapopulation segments - the northern East Pacific Rise (NEPR), southern East Pacific Rise (SEPR), and northeastern Pacific Antarctic Ridge (PAR) - separated by the Equator and Easter Microplate boundaries. Genetic diversity estimators were negatively correlated with tectonic spreading rates. Application of the isolation-with-migration (IMa2) model provided information about divergence times and demographic parameters. The PAR and NEPR metapopulation segments were estimated to have split roughly 4.20 million years ago (Mya) (2.42-33.42 Mya, 95 % highest posterior density, (HPD)), followed by splitting of the SEPR and NEPR segments about 0.79 Mya (0.07-6.67 Mya, 95 % HPD). Estimates of gene flow between the neighboring regions were mostly low (2 Nm < 1). Estimates of effective population size decreased with southern latitudes: NEPR > SEPR > PAR. CONCLUSIONS: Highly effective dispersal capabilities allow A. pompejana to overcome the temporal instability and intermittent distribution of active hydrothermal vents in the eastern Pacific Ocean. Consequently, the species exhibits very high levels of genetic diversity compared with many co-distributed vent annelids and mollusks. Nonetheless, its levels of genetic diversity in partially isolated populations are inversely correlated with tectonic spreading rates. As for many other vent taxa, this pioneering colonizer is similarly affected by local rates of habitat turnover and by major dispersal filters associated with the Equator and the Easter Microplate region.

Complete mitochondrial genome of the hydrothermal vent tubeworm, Ridgeia piscesae (Polychaeta, Siboglinidae).

The complete mitochondrial genome of Ridgeia piscesae (Polychaeta, Siboglinidae), one of the dominant taxa in deep-sea hydrothermal vents, is reported here for the first time. The R. piscesae mitogenome is 15,002 bp in total length and includes 13 protein-coding gene sequences, small and large rRNA sequences and 22 tRNA sequences. All genes are encoded on the heavy strand. The mitochondrial genomes of R. piscesae and other six polychaete species have a conserved gene order.

Complete mitochondrial genome of the headwater livebearer, Poeciliopsis monacha: the mother of clones.

The sexually reproducing fish, Poeciliopsis monacha (Actinopterygii, Cyprinodontiformes, Poeciliidae), is the maternal ancestor of six hybrid biotypes that reproduce clonally. The gene content and order of its 16,818 bp mitochondrial genome is virtually identical with that of other sexually reproducing poeciliid fishes, providing no evidence for a mitochondrial involvement in the origins of all-female reproduction.

Postglacial range shift and demographic expansion of the marine intertidal snail Batillaria attramentaria.

To address the impacts of past climate changes, particularly since the last glacial period, on the history of the distribution and demography of marine species, we investigated the evolutionary and demographic responses of the intertidal batillariid gastropod, Batillaria attramentaria, to these changes, using the snail as a model species in the northwest Pacific. We applied phylogeographic and divergence population genetic approaches to mitochondrial COI sequences from B. attramentaria. To cover much of its distributional range, 197 individuals collected throughout Korea and 507 publically available sequences (mostly from Japan) were used. Finally, a Bayesian skyline plot (BSP) method was applied to reconstruct the demographic history of this species. We found four differentiated geographic groups around Korea, confirming the presence of two distinct, geographically subdivided haplogroups on the Japanese coastlines along the bifurcated routes of the warm Tsushima and Kuroshio Currents. These two haplogroups were estimated to have begun to split approximately 400,000 years ago. Population divergence analysis supported the hypothesis that the Yellow Sea was populated by a northward range expansion of a small fraction of founders that split from a southern ancestral population since the last glacial maximum (LGM: 26,000-19,000 years ago), when the southern area became re-submerged. BSP analyses on six geographically and genetically defined groups in Korea and Japan consistently demonstrated that each group has exponentially increased approximately since the LGM. This study resolved the phylogeography of B. attramentaria as a series of events connected over space and time; while paleoceanographic conditions determining the connectivity of neighboring seas in East Asia are responsible for the vicariance of this species, the postglacial sea-level rise and warming temperatures have played a crucial role in rapid range shifts and broad demographic expansions of its populations.

Genomic replacement of native Cobitis lutheri with introduced C. tetralineata through a hybrid swarm following the artificial connection of river systems.

River connections via artificial canals will bring about secondary contacts between previously isolated fish species. Here, we present a genetic consequence of such a secondary contact between Cobitis fish species, C. lutheri in the Dongjin River, and C. tetralineata in the Seomjin River in Korea. The construction of water canals about 80 years ago has unidirectionally introduced C. tetralineata into the native habitat of C. lutheri, and then these species have hybridized in the main stream section of the Dongjin River. According to the divergence population genetic analyses of DNA sequence data, the two species diverged about 3.3 million years ago, which is interestingly coincident with the unprecedented paleoceanographic change that caused isolations of the paleo-river systems in northeast Asia due to sea-level changes around the late Pliocene. Multilocus genotypic data of nine microsatellites and three nuclear loci revealed an extensively admixed structure in the hybrid zone with a high proportion of various post-F1 hybrids. Surprisingly, pure native C. lutheri was absent in the hybrid zone in contrast to the 7% of pure C. tetralineata. Such a biased proportion must have resulted from the dominant influence of continually introducing C. tetralineata on the native C. lutheri which has no supply of natives from other tributaries to the hybrid zone due to numerous low-head dams. In addition, mating experiments indicated that there is no discernible reproductive isolation between them. All the results suggest that the gene pool of native C. lutheri is being rapidly replaced by that of continually introducing C. tetralineata through a hybrid swarm for the last 80 years, which will ultimately lead to the genomic extinction of natives in this hybrid zone.