Population Genetic Structure and Evidence of Demographic Expansion of the Ayu (Plecoglossus altivelis) in East Asia

Population structure and demographic history for year cohort dynamics of landlocked ayu Plecoglossus altivelis altivelis in dam reservoir of Japan

Ayu Plecoglossus altivelis altivelis is a valuable osmeroid species for inland fishery in Japan. It is classified into two ecological forms of amphidromous migrating between rivers and sea and landlocked migrating between rivers and lakes or dam reservoirs. The number of dams and their reservoirs has remarkably increased in the twenty-first century under climate change, because of their respective roles in hydropower generation with negligible carbon emissions and in flood control. Dam reservoirs therefore become increasingly important as inland nursery grounds of ayu. In this study, we investigated the reproduction status of landlocked ayu migrating in the Haidzuka Dam reservoir and the Tabusa River in western Japan by molecular phylogenetic analysis based on population structure and demographic history for year cohort dynamics. A total of 849 individuals were collected monthly from October 2018 to September 2021 according to an annual life cycle of ayu. Nucleotide sequences of the partial mitochondrial DNA control region yielded 31 haplotypes, consisting of 4 shared haplotypes among the 2019, 2020 and 2021 cohorts and 27 unique haplotypes. The overall haplotype diversity and nucleotide diversity were calculated to be relatively low at 0.3503 ± 0.0206 and 0.0077 ± 0.0045, respectively, suggesting a founder event by dominant haplotypes. Star-shaped radiational haplotypes from dominant shared haplotypes on the median-joining network likely support a founder event. Although pairwise ФST values were determined to be very low among the year cohorts, only the 2019 cohort was found to have a significant difference from the 2020 and 2021 cohorts, for both of which Tajima's D values were also statistically significant. For the overall population, multimodal mismatch distribution and negative Tajima's D and Fu's Fs values in the neutrality test suggested population expansion or population subdivision. The native riverine population in the Tabusa River suffered habitat fragmentation and population bottleneck from dam construction, and therefore severe founder effect remained behind the artificially landlocked population with a low level of genetic diversity in the Haidzuka Dam reservoir.

Improvement of the ayu (Plecoglossus altivelis) draft genome using Hi-C sequencing

OBJECTIVE

The ayu or sweetfish Plecoglossus altivelis is ray-finned fish that is widely distributed in East Asia. The genome size of ayu was estimated at approximately 420 Mb. Previously, we reported on ayu draft genome assembly by whole-genome shotgun using Illumina short reads and PacBio long reads; however, the assembly was not to chromosome level. Therefore, to improve the draft genome sequence of ayu to chromosome level, we performed in situ Hi-C sequencing as a source of linkage information.

RESULTS

The ayu genome assembly yielded 28 large scaffolds that corresponded to the karyotype of ayu (n = 28). The resulting ayu genome assembly has a N50 scaffold length of 17.0 Mb, improved from 4.3 Mb. The high-quality reference genome will be helpful for phylogenetic research on bony fishes and for breeding programs in ayu aquaculture.

A Y-linked anti-Müllerian hormone type-II receptor is the sex-determining gene in ayu, Plecoglossus altivelis

Whole-genome duplication and genome compaction are thought to have played important roles in teleost fish evolution. Ayu (or sweetfish), Plecoglossus altivelis, belongs to the superorder Stomiati, order Osmeriformes. Stomiati is phylogenetically classified as sister taxa of Neoteleostei. Thus, ayu holds an important position in the fish tree of life. Although ayu is economically important for the food industry and recreational fishing in Japan, few genomic resources are available for this species. To address this problem, we produced a draft genome sequence of ayu by whole-genome shotgun sequencing and constructed linkage maps using a genotyping-by-sequencing approach. Syntenic analyses of ayu and other teleost fish provided information about chromosomal rearrangements during the divergence of Stomiati, Protacanthopterygii and Neoteleostei. The size of the ayu genome indicates that genome compaction occurred after the divergence of the family Osmeridae. Ayu has an XX/XY sex-determination system for which we identified sex-associated loci by a genome-wide association study by genotyping-by-sequencing and whole-genome resequencing using wild populations. Genome-wide association mapping using wild ayu populations revealed three sex-linked scaffolds (total, 2.03 Mb). Comparison of whole-genome resequencing mapping coverage between males and females identified male-specific regions in sex-linked scaffolds. A duplicate copy of the anti-Müllerian hormone type-II receptor gene (amhr2bY) was found within these male-specific regions, distinct from the autosomal copy of amhr2. Expression of the Y-linked amhr2 gene was male-specific in sox9b-positive somatic cells surrounding germ cells in undifferentiated gonads, whereas autosomal amhr2 transcripts were detected in somatic cells in sexually undifferentiated gonads of both genetic males and females. Loss-of-function mutation for amhr2bY induced male to female sex reversal. Taken together with the known role of Amh and Amhr2 in sex differentiation, these results indicate that the paralog of amhr2 on the ayu Y chromosome determines genetic sex, and the male-specific amh-amhr2 pathway is critical for testicular differentiation in ayu.

Genetic characteristics of the amphidromous fish Ayu Plecoglossus altivelis altivelis (Osmeriformes: Plecoglossidae) on Yaku-shima Island in Japan, the southernmost population of the subspecies

The Ayu (Plecoglossus altivelis altivelis) population on Yaku-shima Island, Kagoshima Prefecture, Japan represents the southernmost population of the subspecies and is considered to be facing extinction. We investigated the genetic characteristics of Ayu in the Miyanoura River (MYU) in the northeast of the island and the Kurio River (KRO) in the southwest of the island, using partial sequences of mitochondrial (mt) control region and polymorphic microsatellite (simple sequence repeat: SSR) loci. The mtDNA analysis revealed that the Yaku-shima Island population was significantly differentiated from other Ayu populations; additionally, the MYU and KRO populations were significantly different from each other in terms of mtDNA (pairwise Φ_ST = 0.5826, P < 0.01) and SSR (pairwise Φ_ST = 0.1598, P < 0.01) analyses, indicating very little or no gene flow between them. The mtDNA haplotype diversity values were minimal for KRO and somewhat lower for MYU (h = 0.8176) than for each population from the mainland of Japan (Honshu-Kyushu) and the Korean Peninsula (h = 0.9905-1.0000). The mean values of expected heterozygosity (He) of SSRs were also lower in KRO (mean He = 0.555) than in MYU (mean He = 0.649). A considerably small effective population size (N_e = 100.1 for MYU, 151.2 for KRO) and a bottleneck effect for Yaku-shima Island population were suggested by SSR analysis. These findings underscore the importance of monitoring the genetic diversity of Ayu on Yaku-shima Island and the necessity of designing conservation policies for each river's population.

Study on population genetics of Sillago aeolus (Perciformes: Sillaginidae) in the Coast of China

Sillago aeolus is a species from Sillaginidae, which is a widely distributed species with important scientific and economic value in the coast of China. Its population genetics have not been studied. This study investigated the population genetics of S. aeolus in the eastern and southern coast of China based on the mitochondrial control region markers obtained from 248 individuals of 9 locations. The population was characterized by a high population diversity with a low nucleotide diversity. There were no branches corresponding to the sampling sites according to the haplotype network and NJ tree. Recent asymmetric gene flow exchanges and significant genetic differences were detected between the Haikou population and the other populations. AMOVA result indicated slight genetic structures with homogeneity suggested. The neutral test and the mismatch distribution revealed a recent population expansion event. Historical geographic events may be the reason for the homogeneity within the population.

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.