Cold seeps and hydrothermal vents are deep-sea reducing environments that are characterized by lacking oxygen and photosynthesis-derived nutrients. Most animals acquire nutrition in cold seeps or hydrothermal vents by maintaining epi- or endosymbiotic relationship with chemoautotrophic microorganisms. Although several seep- and vent-dwelling animals hosting symbiotic microbes have been well-studied, the genomic basis of adaptation to deep-sea reducing environment in nonsymbiotic animals is still lacking. Here, we report a high-quality genome of Chiridota heheva Pawson & Vance, 2004, which thrives by extracting organic components from sediment detritus and suspended material, as a reference for nonsymbiotic animal’s adaptation to deep-sea reducing environments. The expansion of the aerolysin-like protein family in C. heheva compared with other echinoderms might be involved in the disintegration of microbes during digestion. Moreover, several hypoxia-related genes (Pyruvate Kinase M2, PKM2; Phospholysine Phosphohistidine Inorganic Pyrophosphate Phosphatase, LHPP; Poly(A)-specific Ribonuclease Subunit PAN2, PAN2; and Ribosomal RNA Processing 9, RRP9) were subject to positive selection in the genome of C. heheva, which contributes to their adaptation to hypoxic environments.

Background: Dissecting the genetic basis of variation in the regulation of gene expression is essential for understanding phenotypic evolution. Structural variants intersecting the cis-regulatory elements are found to cause gene expression variation in several developmental genes, resulting in morphological divergence between species. Due to the difficulty of identifying structural variants accurately across the genome, a comprehensive study of impacts of structural variants in cis-regulatory divergence of closely related species, especially fish species, is still scarce. Recently identified broad H3K4me3 domains are essential for the regulation of genes involved in several biological processes. However, the role of broad H3K4me3 domains in phenotypic divergence remain poorly understood. Siniperca chuatsi and S. scherzeri are two closely related fish species diverge in several phenotypic traits, making them an ideal model to study cis-regulatory evolution in closely related species. Results: We generated chromosome-level genomes of S. chuatsi and S. scherzeri. The evolutionary histories of S. chuatsi and S. scherzeri were studied by inferring the dynamic changes in the ancestral population sizes. The genetic basis of adaptation in S. chuatsi and S. scherzeri was dissected by performing gene family expansion and contraction analysis and identifying positively selected genes (PSGs). To investigate the role of SVs in cis-regulatory divergence of closely related fish species, we identified high-quality SVs between S. chuatsi and S. scherzeri, as well as H3K27ac and H3K4me3 domains. Integrated analysis revealed that cis-regulatory divergence caused by SVs played an essential role in the differentiation of metabolism, skin pigmentation, and immunity between S. chuatsi and S. scherzeri. Additionally, divergent broad H3K4me3 domains were found to mostly associate with cancer-related genes in S. chuatsi and S. scherzeri and contribute to their phenotypic divergence. Conclusions: Our analysis reveals SVs play an essential role in cis-regulatory variation between the two sister fish species, which in turn contributes to their phenotypic divergence. The divergence of broad H3K4me3 domains contributes to phenotypic divergence between closely related species. Additionally, the association of broad H3K4me3 domains and cancer-related genes has an ancient origin.

Due to the difficulty in accurately identifying structural variants (SVs) across genomes, their impact on cis -regulatory divergence of closely related species, especially fish, remains to be explored. Recently identified broad H3K4me3 domains are essential for the regulation of genes involved in several biological processes. However, the role of broad H3K4me3 domains in phenotypic divergence remains poorly understood. Siniperca chuatsi and S. scherzeri are closely related but divergent in several phenotypic traits, making them an ideal model to study cis -regulatory evolution in sister species. Here, we generated chromosome-level genomes of S. chuatsi and S. scherzeri , with assembled genome sizes of 716.35 and 740.54 Mb, respectively. The evolutionary histories of S. chuatsi and S. scherzeri were studied by inferring dynamic changes in ancestral population sizes. To explore the genetic basis of adaptation in S. chuatsi and S. scherzeri , we performed gene family expansion and contraction analysis and identified positively selected genes (PSGs). To investigate the role of SVs in cis -regulatory divergence of closely related fish species, we identified high-quality SVs as well as divergent H3K27ac and H3K4me3 domains in the genomes of S. chuatsi and S. scherzeri . Integrated analysis revealed that cis -regulatory divergence caused by SVs played an essential role in phenotypic divergence between S. chuatsi and S. scherzeri . Additionally, divergent broad H3K4me3 domains were mostly associated with cancer-related genes in S. chuatsi and S. scherzeri and contributed to their phenotypic divergence.