Rabbitfish (Siganidae) are coral reef fish that are distributed across diverse habitats that include estuaries, mangroves, reefs, and even seaweed mats. Given their ecological diversity and natural widespread distributions across the Indo-Pacific region, we were interested to investigate the evolutionary history of this group and patterns of divergence that have contributed to their present-day distributions. In the present study, samples were collected from the South China Sea to study taxonomic and phylogenetic relationships, and divergence times. We investigated the taxonomic relationships among modern rabbitfish species, reconstructed their molecular phylogeny, and estimated divergence times among selected lineages based on a fragment of the mtDNA cytochrome oxidase I (COI) and sequences of the nuclear rhodopsin retrogene (RHO). Our results indicate that modern rabbitfish likely originated in the Indo-West Pacific during the late Eocene [37.4 million years ago (mya)], following which they diverged into three major clades during the Pliocene/Pleistocene. Subsequent diversification and origins of the majority of siganids may likely be associated with episodes of paleo-oceanographic events, including greenhouse and glaciation events (Eocene-Miocene) as well as major plate tectonic events (Pliocene-Pleistocene). Some modern siganid species may naturally hybridize with congeneric species where their geographical ranges overlap. A comprehensive taxonomic analysis revealed that the phylogeny of Siganidae (cladogenesis of Clades I, II, and III) is characterized by divergence in several external morphological characters and morphometric parameters. Our study demonstrates that morphological characteristics, geographical heterogeneity, and environmental change have contributed to siganids' historical diversification.

The migration of anadromous fish in heterogenic environments unceasingly imposes a selective pressure that results in genetic variation for local adaptation. However, discrimination of anadromous fish populations by fine-scale local adaptation is challenging because of their high rate of gene flow, highly connected divergent population, and large population size. Recent advances in next-generation sequencing (NGS) have expanded the prospects of defining the weakly structured population of anadromous fish. Therefore, we used NGS-based restriction site-associated DNA (NextRAD) techniques on 300 individuals of an anadromous Hilsa shad (Tenualosa ilisha) species, collected from nine strategic habitats, across their diverse migratory habitats, which include sea, estuary, and different freshwater rivers. The NextRAD technique successfully identified 15,453 single nucleotide polymorphism (SNP) loci. Outlier tests using the FST OutFLANK and pcadapt approaches identified 74 and 449 SNPs (49 SNPs being common), respectively, as putative adaptive loci under a divergent selection process. Our results, based on the different cluster analyses of these putatively adaptive loci, suggested that local adaptation has divided the Hilsa shad population into two genetically structured clusters, in which marine and estuarine collection sites were dominated by individuals of one genetic cluster and different riverine collection sites were dominated by individuals of another genetic cluster. The phylogenetic analysis revealed that all the riverine populations of Hilsa shad were further subdivided into the north-western riverine (turbid freshwater) and the north-eastern riverine (clear freshwater) ecotypes. Among all of the putatively adaptive loci, only 36 loci were observed to be in the coding region, and the encoded genes might be associated with important biological functions related to the local adaptation of Hilsa shad. In summary, our study provides both neutral and adaptive contexts for the observed genetic divergence of Hilsa shad and, consequently, resolves the previous inconclusive findings on their population genetic structure across their diverse migratory habitats. Moreover, the study has clearly demonstrated that NextRAD sequencing is an innovative approach to explore how dispersal and local adaptation can shape genetic divergence of non-model anadromous fish that intersect diverse migratory habitats during their life-history stages.

We present chromosome-level genome assemblies from representative species of three independently evolved seagrass lineages: Posidonia oceanica, Cymodocea nodosa, Thalassia testudinum and Zostera marina. We also include a draft genome of Potamogeton acutifolius, belonging to a freshwater sister lineage to Zosteraceae. All seagrass species share an ancient whole-genome triplication, while additional whole-genome duplications were uncovered for C. nodosa, Z. marina and P. acutifolius. Comparative analysis of selected gene families suggests that the transition from submerged-freshwater to submerged-marine environments mainly involved fine-tuning of multiple processes (such as osmoregulation, salinity, light capture, carbon acquisition and temperature) that all had to happen in parallel, probably explaining why adaptation to a marine lifestyle has been exceedingly rare. Major gene losses related to stomata, volatiles, defence and lignification are probably a consequence of the return to the sea rather than the cause of it. These new genomes will accelerate functional studies and solutions, as continuing losses of the 'savannahs of the sea' are of major concern in times of climate change and loss of biodiversity.