A Chromosome-Level Genome Assembly of the Spotted Scat (Scatophagus argus)

Sequence identification and expression characterization of leptin in the spotted scat, Scatophagus argus

Scatophagus argus is an important marine culture fish in South and South-East Asia, including Southeast coastal areas of China. Artificial propagation technology for S. argus is not optimum; thus further studies on its reproduction biology are required. Although previous studies have shown that leptin (Lep) can regulate fish reproduction, the role of lep genes in S. argus is unknown. Herein, in silico analysis showed that S. argus has two lep genes (lepa and lepb). Protein 3D-structure prediction showed that Lepa has four α-helices (similar to mammals), while Lepb only has three. Tissue distribution analysis showed that lepa is highly expressed in the liver, whereas lepb was not detected in any tissue. Notably, lepr was expressed in all tissues. Lepa mRNA expression levels in the liver and serum Lep, estradiol (E2) and vitellogenin (Vtg) levels of female fish were significantly higher in ovaries at stage IV than in ovaries at stage II. Serum E2 levels were significantly positively correlated with Vtg levels in female fish at different development stages, while serum E2 was not correlated with Lep levels. Consistently, in vitro incubation of the liver with E2 significantly up-regulated vtga, while it did not affect lepa expression. Recombinant Lep (10 nM) significantly up-regulated chicken gonadotropin-releasing hormone (cGnRH/GnRH-II) in the hypothalamus and GnRH receptor (GnRHR) and luteinizing hormone beta (Lhb) in the pituitary. These results suggest that lepa regulates female reproduction in S. argus.

Chromosome-level genome assembly of Nibea coibor using PacBio HiFi reads and Hi-C technologies

Nibea coibor belongs to Sciaenidae and is distributed in the South China Sea, East China Sea, India and the Philippines. In this study, we sequenced the DNA of a male Nibea coibor using PacBio long-read sequencing and generated chromatin interaction data. The genome size of Nibea coibor was estimated to be 611.85~633.88 Mb based on k-mer counts generated with Jellyfish. PacBio sequencing produced 29.26 Gb of HiFi reads, and Hifiasm was used to assemble a 627.60 Mb genome with a contig N50 of 10.66 Mb. We further found the canonical telomeric repeats “TTAGGG” to be present at the telomeres of all 24 chromosomes. The completeness of the assembly was estimated to be 98.9% and 97.8% using BUSCO and Merqury, respectively. Using the combination of ab initio prediction, protein homology and RNAseq annotation, we identified a total of 21,433 protein-coding genes. Phylogenetic analyses showed that Nibea coibor and Nibea albiflora are closely related. The results provide an important basis for research on the genetic breeding and genome evolution of Nibea coibor.

Measurement(s)	Whole-Genome Shotgun Sequencing • transcription profiling assay
Technology Type(s)	single molecular realtime sequencing • RNA sequencing
Sample Characteristic - Organism	Nibea coibor
Sample Characteristic - Location	China

Chromosome-Level Genome Assembly and Transcriptome Comparison Analysis of Cephalopholis sonnerati and Its Related Grouper Species

The tomato hind, Cephalopholis sonnerati, is a bottom-dwelling coral reef fish, which is widely distributed in the Indo-Pacific and Red Sea. C. sonnerati also features complex social structures and behaviour mechanisms. Here, we present a high-quality, chromosome-level genome assembly for C. sonnerati that was derived using PacBio sequencing and Hi-C technologies. A 1043.66 Mb genome with an N50 length of 2.49 Mb was assembled, produced containing 795 contigs assembled into 24 chromosomes. Overall, 97.2% of the complete BUSCOs were identified in the genome. A total of 26,130 protein-coding genes were predicted, of which 94.26% were functionally annotated. Evolutionary analysis revealed that C. sonnerati diverged from its common ancestor with E. lanceolatus and E. akaara approximately 41.7 million years ago. In addition, comparative genome analyses indicated that the expanded gene families were highly enriched in the sensory system. Finally, we found the tissue-specific expression of 8108 genes. We found that these tissue-specific genes were highly enriched in the brain. In brief, the high-quality, chromosome-level reference genome will provide a valuable genome resource for studies of the genetic conservation, resistance breeding, and evolution of C. sonnerati.

Genome-wide survey and genetic characteristics of Ophichthus evermanni based on Illumina sequencing platform

Ophichthidae fishes limit to continental shelf of all tropical and subtropical oceans and contain more than 350 species, representing the greatest specialization diversity in the order Anguiliformes. In the present study, we conducted a genome survey sequencing (GSS) analysis of Ophichthus evermanni by Illumina sequencing platform to briefly reveal its genomic characteristics and phylogenetic relationship. The first de novo assembled 1.97 Gb draft genome of O. evermanni was predicted based on K-mer analysis without obvious nucleotide bias. The heterozygosity ratio was 0.70%, and the sequence repeat ratio was calculated to be 43.30%. A total of 9016 putative coding genes were successfully predicted, in which 3587 unigenes were identified by gene ontology (GO) analysis and 4375 unigenes were classified into cluster of orthologous groups for enkaryotic complete genomes (KOG) functional categories. About 2,812,813 microsatellite motifs including mono-, di-, tri-, tetra-, penta- and hexanucleotide motifs were identified, with an occurrence frequency of 23.32%. The most abundant type was dinucleotide repeat motifs, accounting for 49.19% of the total repeat types. The mitochondrial genome, as a byproduct of GSS, was assembled to investigate the evolutionary relationships between O. evermanni and its relatives. Bayesian inference (BI) phylogenetic tree inferring from concatenated 12 protein-coding genes (PCGs) showed complicated relationships among Ophichthidae species, indicating a polyphyletic origin of the family. The results would achieve more thorough genetic information of snake eels and provide a theoretical basis and reference for further genome-wide analysis of O. evermanni.

High Polymorphism in the Dmrt2a Gene Is Incompletely Sex-Linked in Spotted Scat, Scatophagus argus

Unlike mammals and birds, many fishes have young sex chromosomes, providing excellent models to study sex chromosome differentiation at early stages. Previous studies showed that spotted scat possesses an XX-XY sex determination system. The X has a complete Dmrt3 copy (termed normal) and a truncated copy of Dmrt1 (called Dmrt1b), while the Y has the opposite (normal Dmrt1, which is male-specific, and a truncated Dmrt3 called Dmrt3△-Y). Dmrt1 is the candidate sex determination gene, while the differentiation of other sex-linked genes remains unknown. The spotted scat has proven to be a good model to study the evolution of sex chromosomes in vertebrates. Herein, we sequenced a neighbor gene of this family, Dmrt2, positioned farther from Dmrt1 and closer to Dmrt3 in the spotted scat, and analyzed its sequence variation and expression profiles. The physical locations of the three genes span across an estimated size of >40 kb. The open reading frames of Dmrt2a and its paralog Dmrt2b are 1578 bp and 1311 bp, encoding peptides of 525 and 436 amino acid residues, respectively. Dmrt2a is positioned close to Dmrt3 but farther from Dmrt1 on the same chromosome, while Dmrt2b is not. Sequence analysis revealed several mutations; insertions, and deletions (indels) on Dmrt2a non-coding regions and single-nucleotide polymorphisms (SNPs) on the Dmrt2a transcript. These indels and SNPs are sex-linked and showed high male heterogeneity but do not affect gene translation. The markers designed to span the mutation sites tested on four different populations showed varied concordance with the genetic sexes. Dmrt2a is transcribed solely in the gonads and gills, while Dmrt2b exists in the gonads, hypothalamus, gills, heart, and spleen. The Dmrt2a and Dmrt2b transcripts are profoundly expressed in the male gonads. Analyses of the transcriptome data from five other fish species (Hainan medaka (Oryzias curvinotus), silver sillago (Sillago sihama), Nile tilapia (Oreochromis niloticus), Hong Kong catfish (Clarias fuscus), and spot-fin porcupine fish (Diodon hystrix)) revealed testes-biased expression of Dmrt1 in all, similar to spotted scat. Additionally, the expression of Dmrt2a is higher in the testes than the ovaries in spotted scat and Hainan medaka. The Dmrt2a transcript was not altered in the coding regions as found in Dmrt1 and Dmrt3 in spotted scat. This could be due to the functional importance of Dmrt2a in development. Another possibility is that because Dmrt2a is positioned farther from Dmrt1 and the chromosome is still young, meaning it is only a matter of time before it differentiates. This study undeniably will aid in understanding the functional divergence of the sex-linked genes in fish.

Analysis of circRNA and miRNA expression profiles in IGF3-induced ovarian maturation in spotted scat (Scatophagus argus)

Insulin-like growth factor 3 (IGF3) induces ovarian maturation in teleosts; however, research on its molecular regulatory mechanism remains deficient. Circular RNAs (circRNAs) and microRNAs (miRNAs) are involved in various biological processes, including reproduction. In this study, circRNAs and miRNAs involved in IGF3-induced ovarian maturation were evaluated in spotted scat (Scatophagus argus). In ovarian tissues, we identified 176 differentially expressed (DE) circRNAs and 52 DE miRNAs between IGF3 treatment and control groups. Gene Ontology (GO) enrichment analyses showed that host genes of DE circRNAs and target genes of DE miRNAs were enriched for various processes with a high degree of overlap, including cellular process, reproduction, reproductive process, biological adhesion, growth, extracellular region, cell junction, catalytic activity, and transcription factor activity. Enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways included cell adhesion molecules, ECM-receptor interaction, regulation of actin cytoskeleton, focal adhesion, cell cycle, Hedgehog signaling pathway, phosphatidylinositol signaling system, PI3K-Akt signaling pathway, Apelin signaling pathway, Notch signaling pathway, insulin signaling pathway, and Rap1 signaling pathway. A circRNA-miRNA-mRNA regulatory network was constructed, including DE genes involved in reproduction (e.g., oocyte maturation, oocyte meiosis, and ECM remodeling), such as ccnd2, hecw2, dnm2, irs1, adam12, and cdh13. According to the regulatory network and tissue distribution, we identified one circRNA (Lachesis_group5:6245955|6270787) and three miRNAs (novel_miR_622, novel_miR_980, and novel_miR_64) that may exert regulatory effects in IGF3-induced ovarian maturation in S. argus. Taken together, this study provides a novel insight into the molecular mechanisms by which IGF3 functions in ovaries and highlights the effects of circRNAs and miRNAs in reproduction in S. argus.

Characterization, expression, and regulatory effects of nr0b1a and nr0b1b in spotted scat (Scatophagus argus)

Nuclear receptor subfamily 0 group B member 1 (Nr0b1) belongs to the nuclear receptor (NR) superfamily. It plays critical roles in sex determination, sex differentiation, and gonadal development in mammals. In this study, the duplicated genes nr0b1a and nr0b1b were identified in spotted scat (Scatophagus argus). Phylogenetic and synteny analyses revealed that, unlike nr0b1a, nr0b1b was retained in several species of teleosts after an nr0b1 gene duplication event but was secondarily lost in other fish species, amphibians, reptiles, birds, and mammals. In a sequence analysis, only 1.5 LXXLL-related repeat motifs were identified in spotted scat Nr0b1a, Nr0b1b, and non-mammalian Nr0b1a/Nr0b1, different from the 3.5 repeat motifs in mammalian Nr0b1. By qPCR, nr0b1a and nr0b1b were highly expressed in testes from stages IV to V and in ovaries from stages II to IV, respectively. Male-to-female sex reversal was induced in XY spotted scat by the administration of exogenous E2. A qPCR analysis showed that nr0b1b mRNA expression was higher in sex-reversed XY fish than in control XY fish, with no difference in nr0b1a. A luciferase assay showed that spotted scat Nr0b1a and Nr0b1b did not individually activate cyp19a1a gene transcription. As in mammals, spotted scat Nr0b1a suppressed Nr5a1-mediated cyp19a1a expression, despite containing only 1.5 LXXLL-related repeat motifs in its N-terminal region, while Nr0b1b stimulated Nr5a1-mediated cyp19a1a transcription. These results demonstrated that nr0b1a and nr0b1b in spotted scat have distinct expression patterns and regulatory effects and further indicate that nr0b1b might be involved in ovarian development by regulating Nr5a1-mediated cyp19a1a expression.