Charybdis japonica genome provides insights into desiccation adaptation and sex-determining region

DEAR EDITOR,

Shore swimming crabs, such asCharybdisjaponica,predominantly inhabit the intertidal zones and exhibit high tolerance to desiccation. Here, we present the first chromosome-levelC.japonicagenome, containing 51 pseudochromosomes, with a revised genome length of 1 431.02 Mb and scaffold N50 size of 29.67 Mb. In addition,824.02 Mb of repetitive elements, 30 900 protein-coding genes, 474 microRNAs (miRNAs), 15 570 transfer RNAs(tRNAs), 309 ribosomal RNAs (rRNAs), and 157 small nuclear RNAs (snRNAs) were identified in the genome. Wholegenome resequencing data were used to identify sex-related single-nucleotide polymorphisms (SNPs) and insertiondeletion (indel) mutations. The 0-10 120 000 bp region of chromosome 37 was identified as the sex-determining region ofC.japonica. Comparative genomic analysis identified 832C.japonica-specific gene families. Phylogenetic analysis indicated thatC.japonicawas closely related toPortunus trituberculatus, which also belongs to the family Portunidae.Compared with other species, metabolic rate, oxygen supply,oxidative stress, and various transporter-related genes were expanded or underwent positive selection inC.japonica,which may contribute to its ability to overcome multiple stresses in dry environments. Decoding this genome provides valuable information for revealing the mechanisms underlying desiccation adaptation and sex determination inC.japonicaand enriches current genetic information to explore the evolutionary history and environmental adaptation strategies of other Portunidae crabs.

In this study, we characterized a high-quality chromosomeanchored reference genome ofCharybdisjaponica(H. Milne Edwards, 1861)(Figure 1A) using Illumina short reads, PacBio long reads, and Hi-C reads (see Supplementary Materials and Methods). We correlated the reported genomic data with the biology, evolutionary history, and desiccation tolerance mechanisms ofC.japonica. In addition, we predicted the genomic regions associated with sex determination based on whole-genome resequencing data (see Supplementary Materials and Methods). These data provide insights into the genetic mechanisms associated with sex differentiation and desiccation tolerance inC.japonicaand other crustaceans.

The chromosome-levelC.japonicagenome was 1 431.02 Mb long, with a scaffold N50 of 29.67 Mb (Figure 1B). The genome size was confirmed by survey sequencing (1 398 Mb)after Illumina short read correction and PacBio long read correction. In total, 51 pseudochromosomes were assembled in the genome based on Hi-C reads. Some chromosomes were extremely small and rich in repetitive elements, and most chromosomes in the mitosis metaphase were spot-shaped and difficult to distinguish; therefore, identifying allC.japonicachromosomes through Hi-C interactions was difficult. Many sequencing reads (Illumina short reads and PacBio long reads) were successfully mapped to the assembled genome sequence (Supplementary Table S1), indicating high read consistency. In addition, 86.40% complete BUSCOs(Benchmarking Universal Single-Copy Orthologs) were found in theC.japonicagenome, suggesting that the assembled genome was relatively complete (Supplementary Table S2). A total of 824.02 Mb of repetitive elements were identified,accounting for 57.65% of the genome (Supplementary Table S3). The high proportion of repetitive elements may be a key reason for the large-sized genome ofC.japonica. Previous analysis of whole-genome repetitive sequences and genome size in 44 plants and 68 vertebrates confirmed that the proportion of repeats is positively correlated with genome size(Gao et al., 2018). We also predicted 824.02 Mb of repetitive elements, 30 900 protein-coding genes (Supplementary Table S4), 474 miRNAs, 15 570 tRNAs, 309 rRNAs, and 157 snRNAs (Supplementary Table S5). In conclusion, the assembled genome provides a useful resource for exploring the biological processes ofC.japonicaat the genomic level.

Figure 1 Genome analyses of Charybdis japonica

Based on whole-genome resequencing of 20C.japonicaindividuals (both male and female), we identified nine SNP variants (number: 1 811 614; Supplementary Table S6) and 12 indel variants (number: 1 671 995; Supplementary Table S7).We also predicted the sex-determining region (approximately 10 Mb) on chromosome 37 (Supplementary Figure S1A, B).Notably, several genes in this region may be involved in sex differentiation, including theF-boxand WD repeat domain(WDR) genes. Previous research has indicated that the F-box protein encoded by the recombinantFBXWgene (containingF-boxandWDRgenes) is specific for substrate recognition in ubiquitin-mediated proteolysis (Roberts et al., 2020).Furthermore, the protein can bind with S-phase kinaseassociated protein 1 (SKP1) and cullin to form an evolutionarily conserved Skp-Cullin-F-box (SCF) ubiquitin ligase complex (Supplementary Figure S2), which can selectively bind to degraded proteins and participate in the ubiquitin-proteasome pathway (UPP) (Roberts et al., 2020). In plants, the UPP system causes functional loss of male germ cells by regulating the conversion of functional proteins(Kipreos & Pagano, 2000). The mechanism of sexdetermination in crustaceans is still in the early stage of evolution, and the strength of genetic factors is lower than that of environmental factors. Therefore, the recombinantFBXWregulated UPP may potentially limit testis development,reduce sperm production, and even promote vitellus formation in maleC.japonica. However, this mechanism is speculative and requires further verification.

Our assembled genome can be a useful resource for better understanding the evolutionary history ofC.japonica.Differentiation of conserved single-copy orthologous genes can lead to species divergence. Such genes can therefore be used to explore evolutionary history. In the present study, 340 single-copy orthologous genes were obtained (Figure 1C) and used to determine the phylogenetics ofC.japonica(Figure 1D). The constructed phylogeny was congruent with the prevailing morphological and molecular view that crabs from the family Portunidae (e.g.,P.trituberculatusandC.japonica) cluster together on a single branch. We also identified 832 unique gene families in the genome(Supplementary Table S8), which likely contribute toC.japonica-specific adaptation.

Comparative genomics can provide insights into the unique adaptive plasticity of marine species (Lou et al., 2022).Previous studies have shown that crustaceans may be under constant stress from both dehydration and hypoxia in arid environments. Specifically, in crustaceans, desiccation can disrupt internal osmotic pressure and metabolic capacity and reduce the oxygen-binding ability of hemoglobin, resulting in hypoxic stress (Allen et al., 2012). Prolonged disruption of the respiratory mechanism may further impair immune function,eventually leading to death. Furthermore, antioxidant enzyme activity directly determines the desiccation tolerance of crustaceans. Here, based on annotations, the desiccationadaptive mechanisms ofC.japonicawere supported by the marked expansion of gene families or positive selection of genes related to metabolic rate, oxygen supply, oxidative stress, and transporters.

Carbohydrates provide a carbon skeleton and energy for growth, metabolic activity, and stress resistance of organisms.Glycometabolism connects the metabolism of proteins, lipids,nucleic acids, and secondary substances. Based on comparative genomic analysis, abundant glycometabolismrelated genes were expanded or positively selected inC.japonica, likely contributing to desiccation tolerance. The reasons why carbohydrates may improve desiccation tolerance inC.japonicainclude the following: (i)carbohydrates can accumulate and regulate osmotic pressure and enhance the water retention of cells; (ii) carbohydrates can protect biological substances, such as biofilms and biological macromolecules, to maintain cytoskeletal integrity under drought environments; (iii) carbohydrates can produce energy for resistance to desiccation; and (iv) carbohydrates can be used as a signaling molecule in certain desiccationresistant regulatory mechanisms.

Although hemocyanin (Figure 1E) has a low affinity for oxygen, it remains an essential oxidizing medium for crustacean tissues in the absence of oxygen (Figure 1F).Desiccation can inhibit oxygen supply and induce respiratory acidosis in crustaceans. Respiratory acidosis further reduces the ability of oxygen and hemoglobin to bind, thereby exacerbating oxygen deficiency and tissue hypoxia (Allen et al., 2012). However, desiccation-tolerant crustaceans can regulate ion (including Ca2+and Mg2+) concentrations to improve the oxygen-binding ability of hemocyanin and promote oxygen transfer between tissues. The expansion and positive selection ofIP3andPRPS1indicated thatC.japonicamay maintain Ca2+and Mg2+homeostasis to ensure the oxygen-binding capacity of hemocyanin to resist hypoxia in drought environments. Hypoxia-inducible factors (HIFs) also play important roles in hypoxia response by binding with hypoxia response elements near the downstream gene promoter for transcriptional activation (Martens et al., 2007).The expansion and positive selection ofSTKinC.japonicamay enhance HIF phosphorylation and increase the molecular weight of HIFs. TheGSRgene, which prevents the oxidative decomposition of hemocyanin, was also positively selected and expanded inC.japonica. Glutathione reductase encoded byGSRcan effectively maintain reduced glutathione in cells,which plays an important role in preventing the oxidative decomposition of hemocyanin (Kanzok et al., 2001). In conclusion, the maintenance of hemocyanin levels and oxygen-carrying capacity should helpC.japonicacope with hypoxic conditions in drought environments.

We also detected the expansion or positive selection of severalHSPgenes (includingHSP70andHSP90), which may encode prototypical chaperonin to prevent the accumulation of unfolded proteins and/or assist in the refolding or degradation of denatured proteins (Beck et al., 2009). The zinc finger(ZNF) gene family was also enriched in theC.japonicagenome. Zinc finger proteins are the most abundant proteins in eukaryotes and exhibit a wide range of structures and functions (Laity et al., 2001). Interestingly, positive selection/expansion of RING-type zinc finger proteins occurred inC.japonica. As a zinc finger protein widely found in plants and animals, the RING-type zinc finger protein confers strong drought tolerance inArabidopsisthaliana(Ryu et al., 2010), because it can regulate the UPP system to degrade misfolded proteins and thus improve organismal defense against environmental stress (Takai et al., 2002).

DATA AVAILABILITY

Raw sequencing data (including PacBio, Hi-C, Survey, and RNA-seq reads) forC. japonicagenome were deposited at the National Center for Biotechnology Information (NCBI)sequence read archive (SRA) (SRR16072374, SRR16072376,SRR16072375, and SRR16078881) under BioProjectID PRJNA766329. Raw whole-genome resequencing data of 20C. japonicawere also deposited in the SRA (SRR1689 4624-SRR16894643) under BioProjectID PRJNA779101.Genomic data were collected from the Genome Sequence Archive (GSA) database (https://ngdc.cncb.ac.cn/gwh/Assembly/25212/show) (Accession No. GWHBISL00000000).Genomic data were also archived in the Science Data Bank database (DOI: 10.57760/sciencedb.j00139.00026).

SUPPLEMENTARY DATA

Supplementary data to this article can be found online.

COMPETING INTERESTS

The authors declare that they have no competing interests.

AUTHORS’ CONTRIBUTIONS

Z.Q.H. and F.R.L. conceived and managed the project. Z.Q.H.and T.X.G. collected the sequencing samples. Q.L. extracted the DNA/RNA and performed genome sequencing. F.R.L.analyzed the data. Z.Q.H. and F.R.L. wrote the manuscript. All authors read and approved the final version of the manuscript.