Whole-genome assembly and annotation for the little yellow croaker (Larimichthys polyactis) provide insights into the evolution of hermaphroditism and gonochorism

The evolutionary direction of gonochorism and hermaphroditism is an intriguing mystery to be solved. The special transient hermaphroditic stage makes the little yellow croaker (Larimichthys polyactis) an appealing model for studying hermaphrodite formation. However, the origin and evolutionary relationship between of L. polyactis and Larimichthys crocea, the most famous commercial fish species in East Asia, remain unclear. Here, we report the sequence of the L. polyactis genome, which we found is ~706 Mb long (contig N50 = 1.21 Mb and scaffold N50 = 4.52 Mb) and contains 25,233 protein-coding genes. Phylogenomic analysis suggested that L. polyactis diverged from the common ancestor, L. crocea, approximately 25.4 million years ago. Our high-quality genome assembly enabled comparative genomic analysis, which revealed several within-chromosome rearrangements and translocations, without major chromosome fission or fusion events between the two species. The dmrt1 gene was identified as the male-specific gene in L. polyactis. Transcriptome analysis showed that the expression of dmrt1 and its upstream regulatory gene (rnf183) were both sexually dimorphic. Rnf183, unlike its two paralogues rnf223 and rnf225, is only present in Larimichthys and Lates but not in other teleost species, suggesting that it originated from lineage-specific duplication or was lost in other teleosts. Phylogenetic analysis shows that the hermaphrodite stage in male L. polyactis may be explained by the sequence evolution of dmrt1. Decoding the L. polyactis genome not only provides insight into the genetic underpinnings of hermaphrodite evolution, but also provides valuable information for enhancing fish aquaculture.