Evolution and gene expression of matrix metalloproteinase gene family during gonadal development in Scatophagus argus.

BACKGROUND: During the reproductive cycle of Scatophagus argus (S. argus), their gonads undergo degeneration and re-maturation including the degradation of proteins in the extracellular matrix (ECM). Matrix metalloproteinases (MMPs), a family of zinc proteases, play a crucial role in ECM degradation. OBJECTIVE: Our aim was to identify the MMP gene family of S. argus and determine their gene expression levels across various stages of gonadal development. METHODS: The MMP gene family of S. argus in the genome was identified by using basic Local Alignment Search Tool (BLAST) and HMMER. Phylogenetic tree and synteny analysis were performed to investigate the evolutionary past of the MMP gene family. The gonads of 18 S. argus (9 males and 9 females) were dissected and a quantitative reverse transcription polymerase chain reaction (RT-qPCR) was conducted to investigate the expression levels of MMP genes across different stages of gonadal development. RESULTS: Twenty-three MMP family genes were identified in the genome of S. argus. We divided the MMP gene family into 4 categories and found that teleosts exhibit a higher MMP gene copy number relative to other vertebrates. By sampling S. argus at different stages of gonadal development, we observed an upregulation in relative expression levels of 11 MMP genes in the testis or ovary. Ten MMP genes (mmp2, 9, 14a, 15a, 15b, 16a, 17a, 23b and 24) showed higher mRNA expression in the testis compared to the ovary and mmp28 had higher expression during ovarian development. The tissue distribution results demonstrated that the gills exhibited the lowest relative expression level among all tissues examined. However, 6 genes (mmp2, 9, 14a, 15a, 15b, and 16a) had relatively high expression in all tissues. CONCLUSION: The result suggested that teleost-special whole genome duplication was mainly responsible for the formation of the MMP gene family in teleosts. Expression patterns of MMP genes indicated that mmp2, 9, 14a, 15a, 15b, 16a, 17a, 23b and 24 played a vital role in testicular development while mmp28 was more important for ovarian development. Limitaion: Further studies are needed to determine their protein expressions in gonadal development and precise mechanism in gonadal differentiation. The study enhances our understanding of the MMP gene family in evolution of teleost and provides valuable insights for further research on MMPs in S. argus.

A novel study of brain microvascular endothelial cells induced by astrocyte conditioned medium for constructing blood brain barrier model in vitro: A promising tool for meningitis of teleost.

The blood-brain barrier (BBB) is mainly composed of specialized endothelial cells, which can resist harmful substances, transport nutrients, and maintain the stability of the brain environment. In this study, an endothelial cell line from tilapia (Oreochromis niloticus) named TVEC-01 was successfully established. During the earlier establishment phase of the cell line, the TVEC-01 cells were persistently exposed to an astrocyte-conditioned medium (ACM). TVEC-01 cells were identified as an endothelial cell line. TVEC-01 cells retained the multiple functions of endothelial cells and were capable of performing various experiments in vitro. Furthermore, TVEC-01 cells efficiently expressed BBB-related tight junctions and key efflux transporters. From the results of the qRT-PCR, we found that the TVEC-01 cell line did not gradually lose BBB characteristics after persistent and repetitive passages, which was different from the vast majority of immortalized endothelial cells. The results showed that ACM induced up-regulation of the expression levels of multiple BBB-related genes in TVEC-01 cells. We confirmed that Streptococcus agalactiae was capable of invading the TVEC-01 cells and initiating a series of immune responses, which provided a theoretical basis for S. agalactiae to break through the BBB of teleost through the transcellular traversal pathway. In summary, we have successfully constructed an endothelial cell line of teleost, named TVEC-01, which can be used in many experiments in vitro and even for constructing BBB in vitro. Moreover, it was confirmed that S. agalactiae broke through the BBB of teleost through the transcellular traversal pathway and caused meningitis.

Pigment Identification and Gene Expression Analysis during Erythrophore Development in Spotted Scat (Scatophagus argus) Larvae.

Red coloration is considered an economically important trait in some fish species, including spotted scat, a marine aquaculture fish. Erythrophores are gradually covered by melanophores from the embryonic stage. Despite studies of black spot formation and melanophore coloration in the species, little is known about erythrophore development, which is responsible for red coloration. 1-phenyl 2-thiourea (PTU) is a tyrosinase inhibitor commonly used to inhibit melanogenesis and contribute to the visualization of embryonic development. In this study, spotted scat embryos were treated with 0.003% PTU from 0 to 72 h post fertilization (hpf) to inhibit melanin. Erythrophores were clearly observed during the embryonic stage from 14 to 72 hpf, showing an initial increase (14 to 36 hpf), followed by a gradual decrease (36 to 72 hpf). The number and size of erythrophores at 36 hpf were larger than those at 24 and 72 hpf. At 36 hpf, LC-MS and absorbance spectrophotometry revealed that the carotenoid content was eight times higher than the pteridine content, and ß-carotene and lutein were the main pigments related to red coloration in spotted scat larvae. Compared with their expression in the normal hatching group, rlbp1b, rbp1.1, and rpe65a related to retinol metabolism and soat2 and apoa1 related to steroid hormone biosynthesis and steroid biosynthesis were significantly up-regulated in the PTU group, and rh2 associated with phototransduction was significantly down-regulated. By qRT-PCR, the expression levels of genes involved in carotenoid metabolism (scarb1, plin6, plin2, apoda, bco1, and rep65a), pteridine synthesis (gch2), and chromatophore differentiation (slc2a15b and csf1ra) were significantly higher at 36 hpf than at 24 hpf and 72 hpf, except for bco1. These gene expression profiles were consistent with the developmental changes of erythrophores. These findings provide insights into pigment cell differentiation and gene function in the regulation of red coloration and contribute to selective breeding programs for ornamental aquatic animals.

RNA Sequencing (RNA-Seq) Analysis Reveals Liver Lipid Metabolism Divergent Adaptive Response to Low- and High-Salinity Stress in Spotted Scat (Scatophagus argus).

Spotted scat (Scatophagus argus) can tolerate a wide range of salinity fluctuations. It is a good model for studying environmental salinity adaptation. Lipid metabolism plays an important role in salinity adaptation in fish. To elucidate the mechanism of lipid metabolism in the osmoregulation, the liver transcriptome was analyzed after 22 d culture with a salinity of 5 ppt (Low-salinity group: LS), 25 ppt (Control group: Ctrl), and 35 ppt (High-salinity group: HS) water by using RNA sequencing (RNA-seq) in spotted scat. RNA-seq analysis showed that 1276 and 2768 differentially expressed genes (DEGs) were identified in the LS vs. Ctrl and HS vs. Ctrl, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the pathways of steroid hormone biosynthesis, steroid biosynthesis, glycerophospholipid metabolism, glycerolipid metabolism, and lipid metabolism were significantly enriched in the LS vs. Ctrl. The genes of steroid biosynthesis (sqle, dhcr7, and cyp51a1), steroid hormone biosynthesis (ugt2a1, ugt2a2, ugt2b20, and ugt2b31), and glycerophospholipid metabolism (cept1, pla2g4a, and ptdss2) were significantly down-regulated in the LS vs. Ctrl. The pathways related to lipid metabolisms, such as fatty acid metabolism, fatty acid biosynthesis, peroxisome proliferator-activated receptor (PPAR) signaling pathway, adipocytokine signaling pathway, fatty acid degradation, and unsaturated fatty acid biosynthesis, were significantly enriched in the HS vs. Ctrl. The genes of unsaturated fatty acid biosynthesis (scd1, hacd3, fads2, pecr, and elovl1) and adipocytokine signaling pathway (g6pc1, socs1, socs3, adipor2, pck1, and pparα) were significantly up-regulated in the HS vs. Ctrl. These results suggest that the difference in liver lipid metabolism is important to adapt to low- and high-salinity stress in spotted scat, which clarifies the molecular regulatory mechanisms of salinity adaptation in euryhaline fish.

Expression and regulation of 42Sp50 in spotted scat (Scatophagus argus).

42Sp50 is an isoform of the eukaryotic translation elongation factor 1 A (eEF1A) and is vital for fish ovarian development. Spotted scat (Scatophagus argus) is a popular marine cultured fish species in Southern Asia and China, and its artificial reproduction is complicated, with a relatively low success ratio in practice. In this study, the 42Sp50 gene was cloned from spotted scat. Tissue distribution analysis showed that 42Sp50 was mainly expressed in the ovary. qRT-PCR showed that 42Sp50 expression levels gradually decreased insignificantly in the ovaries from phase II to IV. Western blot analysis showed that 42Sp50 was highly expressed in the ovary, while it was almost undetectable in the testis. Immunohistochemistry analysis stained 42Sp50 mainly in the cytoplasm of the previtellogenic oocytes in ovaries of normal XX-female and sex-reversed XY-female. Aside from fish and amphibians, 42Sp50 was also identified in some reptile species using genomic database searching. Analyses of the transcriptome data from four different fish species (Hainan medaka (Oryzias curvinotus), silver sillago (Sillago sihama), Nile tilapia (Oreochromis niloticus), and Hong Kong catfish (Clarias fuscus)) revealed ovaries biased expression of 42Sp50 in all, similar to spotted scat. While the neighbor genes of 42Sp50 did not show ovary biased expression in the fish species analyzed. Bisulfite Sequencing PCR (BSP) results showed that the DNA methylation level of 42Sp50 promoter was low in ovaries, testes, and muscles. The luciferase reporter assay demonstrated that Dmrt4 activated 42Sp50 expression in the presence of Sf1 or Foxh1. These results suggest that 42Sp50 may be involved in regulating the early phase oocytes development of spotted scat.

First account of a transient intersex in spotted scat, Scatophagus argus: a marine gonochoristic fish.

This study presents the first incidence of intersex associated with testis-ova in spotted scat (Scatophagus argus) reared in a controlled environment. The testis-ova is associated with the abnormal occurrence of primary oocytes (POs) in some male testis and is referred to as ectopic primary oocytes (Ecto-PO), whiles individuals with Ecto-PO are called "Ecto-PO gonad/individuals." We investigated gonads of 129 male spotted scat aged 4-12 and 18 months after hatch (mah) by histological studies for the presence of female sexual characteristics. A total of 20 out of 88 gonads representing 22.7% of male fish aged 6-12, or 15.5% of all male fish sampled, were found to have visible Ecto-PO. At least, the Ecto-PO had an average of 7 oocytes per gonadal section, indicating high severity. The Ecto-PO appears after sex differentiation and degenerates during sexual maturation. The Ecto-PO did not significantly influence the expression pattern of male and female sex-related genes performed using qPCR. Immunofluorescence of 42sp50 specifically stained the Ecto-PO without influence from the surrounding testicular tissues. In addition, temperature did not correlate with the proliferation of the Ecto-PO, but rather gonad developmental strategy. The results show that the naturally occurring Ecto-PO might not be detrimental to testis development and could be considered a frequent-high-level incidence of natural aberration. This study highlights the intricacy of fish sex differentiation and provides a new research chapter to ascertain the mystery behind the occurrence of Ecto-PO.

Establishment of an astrocyte-like cell line from the brain of tilapia (Oreochromis niloticus) for virus pathogenesis and a vitro model of the blood-brain barrier.

A new cell line was established from the brain of a cultured fish, tilapia (Oreochromis niloticus), designated as TA-02 (Tilapia Astrocyte clone 02 cell line). The TA-02 cells are grown for 300 days in an L-15 medium supplemented with 10% fetal bovine serum (FBS). This cell line showed excellent proliferative capacity and expressed various neuroglial cell markers, including SOX2, SOX10, Hes1, Notch1, Occludin, E-cadherin, and GFAP. In addition, TA-02 cells were susceptible to Tilapia Lake Virus (TiLV) as demonstrated by the presence of a severe cytopathic effect (CPE), virus particle in a transmission electron microscope (TEM), and PCR positive signal. Bacterial cytotoxicity studies showed that Streptococcus agalactiae was toxic to TA-02 cells. When co-culture with trans-well, TA-02 exhibited prominent barrier properties, manifested by tight intercellular junctions and increased trans-endothelial electrical resistance (TEER). In addition, the barrier is effective against Escherichia coli (non-meningitis pathogenic bacteria). In contrast, S. agalactiae (meningitis pathogenic bacteria) can pass through the membrane comprising the cells in the trans-well insert. The newly established TA-02 cell line provided a valuable tool for virus pathogenesis and a vitro model of the fish blood-brain barrier.

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.

Homozygous Mutation of gsdf Causes Infertility in Female Nile Tilapia (Oreochromis niloticus).

Gonadal somatic cell-derived factor (Gsdf) is a member of the TGF-ß superfamily, which exists mainly in fishes. Homozygous gsdf mutations in Japanese medaka and zebrafish resulted in infertile females, and the reasons for their infertility remain unknown. This study presents functional studies of Gsdf in ovary development using CRISPR/Cas9 in Nile tilapia (Oreochromis niloticus). The XX wild type (WT) female fish regularly reproduced from 12 months after hatching (mah), while the XX gsdf-/- female fish never reproduced and were infertile. Histological observation showed that at 24 mah, number of phase IV oocyte in the XX gsdf-/- female fish was significantly lower than that of the WT fish, although their gonadosomatic index (GSI) was similar. However, the GSI of the XX gsdf-/- female at 6 mah was higher than that of the WT. The mutated ovaries were hyperplastic with more phase I oocytes. Transcriptome analysis identified 344 and 51 up- and down-regulated genes in mutants compared with the WT ovaries at 6 mah. Some TGF-ß signaling genes that are critical for ovary development in fish were differentially expressed. Genes such as amh and amhr2 were up-regulated, while inhbb and acvr2a were down-regulated in mutant ovaries. The cyp19a1a, the key gene for estrogen synthesis, was not differentially expressed. Moreover, the serum 17ß-estradiol (E2) concentrations between XX gsdf-/- and WT were similar at 6 and 24 mah. Results from real-time PCR and immunofluorescence experiments were similar and validated the transcriptome data. Furthermore, Yeast-two-hybrid assays showed that Gsdf interacts with TGF-ß type II receptors (Amhr2 and Bmpr2a). Altogether, these results suggest that Gsdf functions together with TGF-ß signaling pathway to control ovary development and fertility. This study contributes to knowledge on the function of Gsdf in fish oogenesis.

RNA Sequencing Analysis Reveals Divergent Adaptive Response to Hypo- and Hyper-Salinity in Greater Amberjack (Seriola dumerili) Juveniles.

Salinity significantly affects physiological and metabolic activities, breeding, development, survival, and growth of marine fish. The greater amberjack (Seriola dumerili) is a fast-growing species that has immensely contributed to global aquaculture diversification. However, the tolerance, adaptation, and molecular responses of greater amberjack to salinity are unclear. This study reared greater amberjack juveniles under different salinity stresses (40, 30, 20, and 10 ppt) for 30 days to assess their tolerance, adaptation, and molecular responses to salinity. RNA sequencing analysis of gill tissue was used to identify genes and biological processes involved in greater amberjack response to salinity stress at 40, 30, and 20 ppt. Eighteen differentially expressed genes (DEGs) (nine upregulated and nine downregulated) were identified in the 40 vs. 30 ppt group. Moreover, 417 DEGs (205 up-regulated and 212 down-regulated) were identified in the 20 vs. 30 ppt group. qPCR and transcriptomic analysis indicated that salinity stress affected the expression of genes involved in steroid biosynthesis (ebp, sqle, lss, dhcr7, dhcr24, and cyp51a1), lipid metabolism (msmo1, nsdhl, ogdh, and edar), ion transporters (slc25a48, slc37a4, slc44a4, and apq4), and immune response (wnt4 and tlr5). Furthermore, KEGG pathway enrichment analysis showed that the DEGs were enriched in steroid biosynthesis, lipids metabolism, cytokine-cytokine receptor interaction, tryptophan metabolism, and insulin signaling pathway. Therefore, this study provides insights into the molecular mechanisms of marine fish adaptation to salinity.

Transcriptomic analysis of pituitary in female and male spotted scat (Scatophagus argus) after 17ß-estradiol injection.

Spotted scat (Scatophagus argus) is a popular species of marine fish cultured in China. It shows normal sexual growth dimorphism. Female spotted scat grows quicker and bigger than males. Growth and reproduction are the most important traits in aquaculture. In vertebrates, the pituitary gland occupies an important position in the growth and reproduction axis. Estrogen is involved in regulating growth and reproduction in the pituitary gland in an endocrine fashion. Transcriptome sequencing of the pituitary was performed in female and male fish at 6 h after 17ß-estradiol injection (4.0 µg E2/g body weight, BW). Compared with the pituitary of female and male groups, 144 and 64 genes [|log2(fold change)| ≥ 1.0 and false discovery rate (FDR) < 0.05] were significantly differentially expressed in E2-injected females and males, respectively (p < 0.05). Of these, 59 and 48 were up-regulated, and 85 and 16 were down-regulated. According to the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) pathway analyses, DEGs were involved in signal pathways, such as growth, reproduction, oocyte meiosis and steroid biosynthesis. Of these, estrogen affected the expression of some sex steroid synthesis and receptor genes in the pituitary gland through feedback, such as hsd17b7, pgr and cyp19a1b, regulating the reproductive activities. Besides, some growth-related genes, such as gap43, junbb, mstn2 and insm1a responded to estrogen. E2 might affect the expression level of gh mRNA by regulating the expression levels of growth-related genes. Our results provide a theoretical basis for studying the molecular mechanism of growth and reproduction regulation at the pituitary level of spotted scat responded to E2.

Comparative Physiological and Transcriptomic Profiling Offers Insight into the Sexual Dimorphism of Hepatic Metabolism in Size-Dimorphic Spotted Scat (Scatophagus argus).

The spotted scat (Scatophagus argus) is an economically important cultured marine fish that exhibits a typical sexual size dimorphism (SSD). SSD has captivated considerable curiosity for farmed fish production; however, up till now the exact underlying mechanism remains largely unclear. As an important digestive and metabolic organ, the liver plays key roles in the regulation of fish growth. It is necessary to elucidate its significance as a downstream component of the hypothalamic-pituitary-liver axis in the formation of SSD. In this study, the liver physiological differences between the sexes were evaluated in S. argus, and the activity of several digestive and metabolic enzymes were affected by sex. Females had higher amylase, protease, and glucose-6-phosphate dehydrogenase activities, while males exhibited markedly higher hepatic lipase and antioxidant enzymes activities. A comparative transcriptomics was then performed to characterize the responsive genes. Illumina sequencing generated 272.6 million clean reads, which were assembled into 79,115 unigenes. A total of 259 differentially expressed genes were identified and a few growth-controlling genes such as igf1 and igfbp1 exhibited female-biased expression. Further analyses showed that several GO terms and pathways associated with metabolic process, particularly lipid and energy metabolisms, were significantly enriched. The male liver showed a more active mitochondrial energy metabolism, implicating an increased energy expenditure associated with reproduction. Collectively, the female-biased growth dimorphism of S. argus may be partially attributed to sexually dimorphic metabolism in the liver. These findings would facilitate further understanding of the nature of SSD in teleost fish.

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.

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

The spotted scat, Scatophagus argus is a member of the family Scatophagidae found in Indo-Pacific coastal waters. It is an emerging commercial aquaculture species, particularly in East and Southeast Asia. In this study, the first chromosome-level genome of S. argus was constructed using PacBio and Hi-C sequencing technologies. The genome is 572.42 Mb, with a scaffold N50 of 24.67 Mb. Using Hi-C data, 563.28 Mb (98.67% of the genome) sequences were anchored and oriented in 24 chromosomes, ranging from 12.57 Mb to 30.38 Mb. The assembly is of high integrity, containing 94.26% conserved single-copy orthologues, based on BUSCO analysis. A total of 24,256 protein-coding genes were predicted in the genome, and 96.30% of the predicted genes were functionally annotated. Evolutionary analysis showed that S. argus diverged from the common ancestor of Japanese puffer (Takifugu rubripes) approximately 114.8 Ma. The chromosomes of S. argus showed significant correlation to T. rubripes chromosomes. A comparative genomic analysis identified 49 unique and 90 expanded gene families. These genomic resources provide a solid foundation for functional genomics studies to decipher the economic traits of this species.

Chromosomal-Level Genome Assembly of Silver Sillago (Sillago sihama).

Silver sillago, Sillago sihama is a member of the family Sillaginidae and found in all Chinese inshore waters. It is an emerging commercial marine aquaculture species in China. In this study, high-quality chromosome-level reference genome of S. sihama was first constructed using PacBio Sequel sequencing and high-throughput chromosome conformation capture (Hi-C) technique. A total of 66.16 Gb clean reads were generated by PacBio sequencing platforms. The genome-scale was 521.63 Mb with 556 contigs, and 13.54 Mb of contig N50 length. Additionally, Hi-C scaffolding of the genome resulted in 24 chromosomes containing 96.93% of the total assembled sequences. A total of 23,959 protein-coding genes were predicted in the genome, and 96.51% of the genes were functionally annotated in public databases. A total of 71.86 Mb repetitive elements were detected, accounting for 13.78% of the genome. The phylogenetic relationships of silver sillago with other teleosts showed that silver sillago was separated from the common ancestor of Sillago sinica ∼7.92 Ma. Comparative genomic analysis of silver sillago with other teleosts showed that 45 unique and 100 expansion gene families were identified in silver sillago. In this study, the genomic resources provide valuable reference genomes for functional genomics research of silver sillago.

The reproductive regulation of LPXRFa and its receptor in the hypothalamo-pituitary-gonadal axis of the spotted scat (Scatophagus argus).

Gonadotropin-inhibitory hormone (GnIH) plays a critical role in regulating gonadotropin-releasing hormone (GnRH), gonadotropin hormone (GtH), and steroidogenesis. The Lpxrfa (the piscine ortholog of GnIH) system has been found to regulate fish reproduction. To gain insight into the role of Lpxrfa in the regulation of spotted scat (Scatophagus argus) reproduction, spotted scat Lpxrfa (ssLpxrfa), and its receptor (ssLpxrfa-r) were cloned and analyzed. Tissue distribution and expression patterns at the hypothalamo-pituitary-gonadal axis (HPG axis) of sslpxrfa and sslpxrfa-r mRNA were also investigated during gonadal development of spotted scat. The open reading frame (ORF) of the sslpxrfa was 606 bp encoding 201 amino acids and includes a putative signal peptide and two mature ssLpxrfa peptides with LPXRFamide motif at their C-terminus. The sslpxrfa-r ORF was 1449 bp encoding 482 amino acids and contracted a seven-hydrophobic transmembrane (TM) domain structure. The tissue distribution showe d that the sslpxrfa was highly expressed in hypothalami, gill, and the gonads. In addition, sslpxrfa-r was highly expressed in hypothalami, pituitaries, and the gonads. Quantitative real-time polymerase chain reaction (qPCR) revealed that sslpxrfa had the highest expression in the hypothalami and pituitaries, and the lowest expression in the gonads in stage V. During gonadal development, the expression of sslpxrfa-r was gradually increased in the hypothalami but reduced in the gonads. However, no obvious trend was observed in the pituitaries. The expression of sslpxrfa and sslpxrfa-r decreased significantly after injection with 17ß-estradiol (E2). However, the expression of both sslpxrfa and sslpxrfa-r was not changed after injection with 17α-methyltestosterone(17α-MT) in the hypothalami. In addition, no changes were observed in the expression of fshß and lhß in the pituitaries after injecting ssLpxrfa-1. However, ssLpxrfa-2 could downregulate the expression of sbgnrh and fshß in the hypothalami and pituitaries, respectively. Taken together, these findings suggested that ssLpxrfa may participate in E2 feedback in reproduction and regulate the reproductive axis of spotted scat.

Effects of 17ß-Estradiol on growth-related genes expression in female and male spotted scat (Scatophagus argus).

Growth hormone (GH) is the most important endocrine factor to regulate somatic growth. Spotted scat (Scatophagus argus) is a famous marine aquaculture species in China with a typical sexual growth dimorphism in which females grow faster and larger than males. In this study, gh messenger RNA (gh mRNA) and GH protein expression were examined in the pituitary glands of female and male spotted scat. Based on qPCR analysis, gh mRNA was mainly expressed in the pituitary gland, and weakly in the gonads and hypothalamus. Furthermore, gh mRNA expression in the pituitary gland was significantly higher in females at stages II-IV than in males at stages III-V. In addition, gh mRNA was highly expressed in the ovary and testis during mature development stages. In this study, spotted scat GH polyclonal antibody was produced. Western blot analysis showed that the molecular weight of spotted scat GH was about 21 KDa. Immunohistochemistry (IHC) in pituitary glands showed that GH was mainly expressed in the proximal pars distal (PPD) and a few cells were distributed in the rostral pairs distal (RPD). After injecting 17ß-Estradiol (E2) in vivo, gh mRNA expression was significantly up-regulated in the pituitary gland, whereas igf1 and ghr1 mRNA levels were down-regulated in the liver, which might regulate gh mRNA expression in the pituitary gland. These results provide valuable insight into the molecular mechanisms of E2 regulating gh expression in spotted scat.

Transcriptome analysis of liver provides insight into metabolic and translation changes under hypoxia and reoxygenation stress in silver sillago (Sillago sihama).

Hypoxia can lead to adverse effects on growth, reproduction, behavioral activities and survival in fish, and is one of the most critical factors in the aquatic environment. The liver is an important target organ for reducing toxin accumulation and hypoxia in fish. In this study, silver sillago (Sillago sihama) was exposed to normoxia (dissolved oxygen, DO = 8.0 mg/L), hypoxia for 1 h (hypoxia 1 h, DO = 1.5 mg/L), hypoxia for 4 h (hypoxia 4 h, DO = 1.5 mg/L) and reoxygenation for 4 h after hypoxia 4 h (reoxygenation 4 h, DO = 8.0 mg/L). Results showed that the expression of 506, 1721, and 1230 differentially expressed genes (DEGs) (|log2(fold change) > 1.0| and padj < 0.05) were identified at hypoxia 1 h, hypoxia 4 h, and reoxygenation 4 h in the liver, respectively. The enrichment analysis showed that the DEGs were significantly enriched in metabolic and translation changes pathways, including mapk signaling pathway, p53 signaling pathway, fatty acid metabolism, protein export, ribosome biogenesis in eukaryotes. The DEGs of 17 genes validated the RNA-seq results by quantitative real-time PCR (qRT-PCR). This study provides a comprehensive understanding of the transcriptional changes that occur in different hypoxia and insights into the mechanisms of hypoxia adaptation of the liver in S. sihama.

Identification, functional characterization, and estrogen regulation on gonadotropin-releasing hormone in the spotted scat, Scatophagus argus.

Gonadotropin-releasing hormone (GnRH) is a key neuropeptide of the reproductive system. However, little is known about the role of GnRH in the spotted scat (Scatophagus argus). Here, three GnRH subtypes (cGnRH-II, sGnRH, and sbGnRH) were identified in the spotted scat. cGnRH-II and sGnRH were only expressed in the brains and gonads of both male and female fish, exhibiting a tissue-specific expression pattern, while sbGnRH was expressed at different transcription levels in all examined tissues. During ovarian maturation, hypothalamus-associated sbGnRH was upregulated, while the expression of sGnRH was variable and cGnRH-II first increased and then decreased. In vivo experiments showed that sbGnRH significantly promoted the expression of fsh and lh genes in a dose-dependent manner and exhibited a desensitization effect on lh expression at high concentrations. For sGnRH and cGnRH-II, only high concentrations could induce fsh and lh expression. Furthermore, treatment with highly concentrated sbGnRH peptide also induced fsh and lh expression, whereas the sGnRH and cGnRH-II peptides only induced fsh expression in vitro. 17ß-Estradiol (E2) significantly inhibited the expression of sbGnRH mRNA in a dose-dependent manner and did not impact sGnRH and cGnRH-II mRNA levels in vivo or in vitro. The inhibitory effect of E2 on sbGnRH expression was attenuated by the estrogen receptor (ER) broad-spectrum antagonist (fulvestrant) and the ERα-specific antagonist (methyl-piperidinopyrazole), respectively, implying that the feedback regulation on sbGnRH is mediated via ERα. This study provides a theoretical basis for the reproductive endocrinology of the spotted scat by studying GnRH.

ddRADseq-assisted construction of a high-density SNP genetic map and QTL fine mapping for growth-related traits in the spotted scat (Scatophagus argus).

BACKGROUND: Scatophagus argus is a popular farmed fish in several countries of Southeast Asia, including China. Although S. argus has a highly promising economic value, a significant lag of breeding research severely obstructs the sustainable development of aquaculture industry. As one of the most important economic traits, growth traits are controlled by multiple gene loci called quantitative trait loci (QTLs). It is urgently needed to launch a marker assisted selection (MAS) breeding program to improve growth and other pivotal traits. Thus a high-density genetic linkage map is necessary for the fine mapping of QTLs associated with target traits. RESULTS: Using restriction site-associated DNA sequencing, 6196 single nucleotide polymorphism (SNP) markers were developed from a full-sib mapping population for genetic map construction. A total of 6193 SNPs were grouped into 24 linkage groups (LGs), and the total length reached 2191.65 cM with an average marker interval of 0.35 cM. Comparative genome mapping revealed 23 one-to-one and 1 one-to-two syntenic relationships between S. argus LGs and Larimichthys crocea chromosomes. Based on the high-quality linkage map, a total of 44 QTLs associated with growth-related traits were identified on 11 LGs. Of which, 19 significant QTLs for body weight were detected on 9 LGs, explaining 8.8-19.6% of phenotypic variances. Within genomic regions flanking the SNP markers in QTL intervals, we predicted 15 candidate genes showing potential relationships with growth, such as Hbp1, Vgll4 and Pim3, which merit further functional exploration. CONCLUSIONS: The first SNP genetic map with a fine resolution of 0.35 cM for S. argus has been developed, which shows a high level of syntenic relationship with L. crocea genomes. This map can provide valuable information for future genetic, genomic and evolutionary studies. The QTLs and SNP markers significantly associated with growth-related traits will act as useful tools in gene mapping, map-based cloning and MAS breeding to speed up the genetic improvement in important traits of S. argus. The interesting candidate genes are promising for further investigations and have the potential to provide deeper insights into growth regulation in the future.