Population Genomics of an Anadromous Hilsa Shad Tenualosa ilisha Species across Its Diverse Migratory Habitats: Discrimination by Fine-Scale Local Adaptation

Seascape genomics reveals limited dispersal and suggests spatially varying selection among European populations of sea lamprey (Petromyzon marinus)

Sea lamprey Petromyzon marinus is an anadromous and semelparous fish without homing behaviors. Despite being a freshwater, free-living organism for a large part of their life cycle, its adulthood is spent as a parasite of marine vertebrates. In their native European range, while it is well-established that sea lampreys comprise a single nearly-panmictic population, few studies have further explored the evolutionary history of natural populations. Here, we performed the first genome-wide characterization of sea lamprey's genetic diversity in their European natural range. The objectives were to investigate the connectivity among river basins and explore evolutionary processes mediating dispersal during the marine phase, with the sequencing of 186 individuals from 8 locations spanning the North Eastern Atlantic coast and the North Sea with double-digest RAD-sequencing, obtaining a total of 30,910 bi-allelic SNPs. Population genetic analyses reinforced the existence of a single metapopulation encompassing freshwater spawning sites within the North Eastern Atlantic and the North Sea, though the prevalence of private alleles at northern latitudes suggested some limits to the species' dispersal. Seascape genomics suggested a scenario where oxygen concentration and river runoffs impose spatially varying selection across their distribution range. Exploring associations with the abundance of potential hosts further suggested that hake and cod could also impose selective pressures, although the nature of such putative biotic interactions was unresolved. Overall, the identification of adaptive seascapes in a panmictic anadromous species could contribute to conservation practices by providing information for restoration activities to mitigate local extinctions on freshwater sites.

Cross-species gene enrichment revealed a single population of Hilsa shad (Tenualosa ilisha) with low genetic variation in Bangladesh waters

Tenualosa ilisha is a popular anadromous and significant trans-boundary fish. For sustainable management and conservation of this fish, drawing an appropriate picture reflecting population status of this species is very essential based on their all-strategic habitats in Bangladesh. In this study, 139 samples from 18 sites were collected and cross-species gene enrichment method was applied. Like most of the Clupeiforms, nucleotide diversity of this shad was very low (0.001245-0.006612). Population differences between most of the locations were low and not significant (P > 0.05). However, P values of a few locations were significant (P < 0.05) but their pairwise FST values were very poor (0.0042-0.0993), which is inadequate to recognize any local populations. Our study revealed that the presence of a single population in the Bangladesh waters with some admixtured individuals, which may contain partial genes from other populations. Most of the individuals were admixed without showing any precise grouping in the ML IQtree and Network, which might due to their highly migratory nature. Fishes from haors and small coastal rivers were not unique and no genetic differences between migratory cohorts. The hilsa shad fishery should be managed considering it as a single panmictic population in Bangladesh with low genetic diversity.

Skim-Sequencing Based Genotyping Reveals Genetic Divergence of the Wild and Domesticated Population of Black Tiger Shrimp (Penaeus monodon) in the Indo-Pacific Region

The domestication of a wild-caught aquatic animal is an evolutionary process, which results in genetic discrimination at the genomic level in response to strong artificial selection. Although black tiger shrimp (Penaeus monodon) is one of the most commercially important aquaculture species, a systematic assessment of genetic divergence and structure of wild-caught and domesticated broodstock populations of the species is yet to be documented. Therefore, we used skim sequencing (SkimSeq) based genotyping approach to investigate the genetic structure of 50 broodstock individuals of P. monodon species, collected from five sampling sites (n = 10 in each site) across their distribution in Indo-Pacific regions. The wild-caught P. monodon broodstock population were collected from Malaysia (MS) and Japan (MJ), while domesticated broodstock populations were collected from Madagascar (MMD), Hawaii, HI, USA (MMO), and Thailand (MT). After various filtering process, a total of 194,259 single nucleotide polymorphism (SNP) loci were identified, in which 4983 SNP loci were identified as putatively adaptive by the pcadapt approach. In both datasets, pairwise F_ST estimates high genetic divergence between wild and domesticated broodstock populations. Consistently, different spatial clustering analyses in both datasets categorized divergent genetic structure into two clusters: (1) wild-caught populations (MS and MJ), and (2) domesticated populations (MMD, MMO and MT). Among 4983 putatively adaptive SNP loci, only 50 loci were observed to be in the coding region. The gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses suggested that non-synonymous mutated genes might be associated with the energy production, metabolic functions, respiration regulation and developmental rates, which likely act to promote adaptation to the strong artificial selection during the domestication process. This study has demonstrated the applicability of SkimSeq in a highly duplicated genome of P. monodon specifically, across a range of genetic backgrounds and geographical distributions, and would be useful for future genetic improvement program of this species in aquaculture.