The draft genome of the large yellow croaker reveals well-developed innate immunity

Molecular and functional identification of a LC3C homolog in large yellow croaker (Larimichthys crocea)

Autophagy is a conserved cellular process in response to stress that sustains normal cell growth by degrading and recycling unnecessary cytosolic components in autophagosomes and autolysosomes. Microtubule-associated protein 1-light chain 3, MAP1LC3 or LC3, is a major autophagosome component and reliable autophagy marker. In humans, there exist three LC3 orthologues, LC3A, LC3B, and LC3C. To date, only two types of LC3 orthologues, LC3A and LC3B, have been recognized in fish species. This study identified a LC3C (LcLC3C) gene in the large yellow croaker (Larimichthys crocea). The deduced LcLC3C protein possesses conserved characteristic domains, including a GABARAP domain, a C-terminal glycine residue, Atg7 binding site, tubulin binding site, and a lipid site. The LcLC3C gene was constitutively expressed in all examined tissues, and in primary head kidney macrophages (PKM), primary head kidney lymphocytes (PKL), primary head kidney granulocytes (PKG), and large yellow croaker head kidney (LYCK) cell line. After stimulation by poly (I:C), a viral dsRNA mimic, and Vibrio alginolyticus, the mRNA levels of the LcLC3C were downregulated in the head kidney and spleen. Subcellular location indicated that LcLC3C was evenly distributed in the cytoplasm and nucleus of both LYCK cells and epithelioma papulosum cyprini (EPC) cells. Overexpression of LcLC3C in EPC cells promoted the replication of spring viremia of carp virus (SVCV), and its overexpression in both LYCK cells and EPC cells downregulated type I IFN response, suggesting that LcLC3C may promote the replication of virus by negatively regulating type I IFN response. Consequently, our findings enhance the comprehension of the role played by the LC3C homolog in fish.

Target miRNA identification for the LPL gene in large yellow croaker (Larimichthys crocea)

MicroRNA (miRNA), a conservatively evolved single-stranded non-coding RNA, exerts pivotal control over the appearance of target genes and several biological processes. This study conducted a comprehensive screening of candidate microRNAs (miRNAs) associated with Lipoprotein Lipase (LPL) in the large yellow croaker (Larimichthys crocea), utilizing sophisticated bioinformatics techniques across the species' muscular and hepatic tissues. The bioinformatics analysis facilitated the compilation and examination of miRNA datasets specific to these tissues. The investigation culminated in the identification of miR-84a and miR-1231-5p as key miRNAs that modulate fat hydrolysis, highlighting their potential roles in lipid metabolism. Subsequent in-depth analysis further implicated these miRNAs, along with miR-891a, as prospective targets of LPL, suggesting their integral involvement in the regulation of this critical enzyme. Validation of these bioinformatics predictions was conducted through the construction of double luciferase reporters concealing the LPL 3' untranslated region (3'UTR), substantiating that miR-84a and miR-1231-5p can modulate LPL expression via the LPL 3'UTR. Conversely, miR-891a was not concerned with this regulatory mechanism. Site-directed mutagenesis experiments elucidated the specificity of the interaction sequences. Quantitative PCR assays suggested that miR-84a and miR-1231-5p might influence LPL expression during the starvation phase, intimating the regulatory role of miRNA in fatty acid metabolism within hepatic and muscular tissue under starvation. These findings offer a nuanced understanding of LPL's molecular functionality under stress conditions in fish, emphasizing the regulatory dynamics of miRNA during metabolic stress.

Efficacy of Larimichthys crocea TASOR protein-derived peptide FAM286 against Staphylococcus aureus

Staphylococcus aureus (S. aureus) is a major foodborne pathogen, could lead cause of intestinal infections in humans. Antimicrobial peptides, as emerging antimicrobial agents, are gradually replacing traditional agents due to their highly effective and safe antimicrobial activity. In this study, a novel peptide, designated as FAM286, was identified from TASOR protein of Larimichthys crocea, which had a strong antimicrobial activity against S. aureus with a minimum bactericidal concentration (MBC) of 3.9 μg/mL and completed bacteria killing by treatment for 1.5 h. The FAM286 could increase the permeability and disrupt the integrity of cell membranes. The cell showed aggregation phenomenon and entered the apoptosis stage. In addition, the non-bactericidal concentration of FAM286 could effectively inhibit the formation of biofilm and remove mature biofilm. Molecular docking experiments further verified the binding sites of FAM286 to S. aureus biofilm proteins SarA, AgrA, and Hld. FAM286 could also bind the bacterial DNA in an embedded manner, disrupting the structure of DNA and leading to the death of bacteria. This study comprehensively evaluated the antimicrobial mechanism of the FAM286 against S. aureus and provided a theoretical basis for the prevention and control of S. aureus.

Palmitoylation-mediated NLRP3 inflammasome activation in teleosts highlights evolutionary divergence in immune regulation

NLRP3 inflammasome activation is pivotal for cytokine secretion and pyroptosis in response to diverse stimuli, playing a crucial role in innate immunity. While extensively studied in mammals, the regulatory mechanisms governing NLRP3 activation in non-mammalian vertebrates remain largely unexplored. Teleosts, as basal vertebrates, represent an ideal model for exploring the evolutionary trajectory of inflammasome regulation. In this study, ABE assays, confocal microscopy, and biochemical analyses were applied to systematically characterize the mechanisms underlying NLRP3 inflammasome in teleosts, using large yellow croakers ( Larimichthys crocea, Lc) and zebrafish ( Danio rerio, Dr) as representative models. Our findings revealed a previously unrecognized palmitoylation-dependent regulatory mechanism essential for teleost NLRP3 activation. Specifically, zDHHC18-mediated palmitoylation at a teleost-specific cysteine residue (C946 in LcNLRP3, C1037 in DrNLRP3) was required for the translocation of NLRP3 to the dispersed trans-Golgi network, facilitating its subsequent recruitment to the microtubule-organizing center. This membrane trafficking was crucial for inflammasome assembly and downstream inflammatory responses. These findings provide new insights into the distinct regulatory mechanisms of NLRP3 activation in teleosts, highlighting an evolutionary divergence that contributes to innate immunity adaptation in early vertebrates.

Increased intake of marine fish contributed to a decreased odds of comorbid depressive symptoms and coronary heart disease in Chinese adults

BACKGROUNDS

Increased consumption of fish has beneficial impacts upon emotional health; however, this benefit for comorbid depressive symptoms and coronary heart disease (DCHD) is not fully clear. We aimed to investigate the relationship between consumption of marine fish and DCHD in Chinese adults.

METHODS

A cross-sectional study was conducted in 1,106 participants aged 25-95 years living in Taizhou, China. Fish intakes were assessed by using a validated food frequency questionnaire, with their tertiles as category levels of ≤1 time/week, 2-6 times/week, and ≥7 times/week. Coronary heart disease (CHD) was diagnosed using the coronary angiography, while the concurrent depressive symptoms was indicated using ≥8 scores from hospital anxiety and depression scales (HADS). Primary measurements were the prevalent DCHD, presented as multivariate-adjusted odds ratios (ORs) with 95% confidence intervals (CIs).

RESULTS

A total of 932 participants were included, 88 (9.44%) participants with depressive symptoms, 477 (51.18%) CHD, and 106 (11.37%) DCHD, respectively. Participants at the highest tertile of fish intake have a lower odds of DCHD compared with those at the lowest (OR: 0.34, 95% CI: 0.20, 0.58), with 42% reductions in odds of DCHD for per one-tertile (3 times/week) increase (OR: 0.58, 95% CI: 0.45, 0.76). The beneficial associations were pronounced with decreased odds of depressive symptom (OR: 0.31, 95% CI: 0.20, 0.47), but not with CHD (OR: 0.87, 95% CI: 0.59, 1.29).

CONCLUSIONS

Increased consumption of marine fish is associated with decreased severity of depressive symptoms, which might have great benefits toward comorbid depressive symptom and with coronary heart diseases.

Chromosome-level genome assembly for three geographical stocks of large yellow croaker (Larimichthys crocea)

Large yellow croaker (Larimichthys crocea) has been demonstrated to be divided into three geographical stocks from south to north along the coast of China, including Nanhai, Mindong, and Daiqu. Although multiple versions of L. crocea have been published, no high-quality Nanhai and Daiqu genomes have been assembled, hampering the assessment of the fine-scale genetic structure and adversely affecting wild stock conservation, fishery management, and germplasm exploitation of large yellow croaker. To fill the gap, we sequenced the genomes of three L. crocea stocks using a combination of PacBio and Hi-C technologies. We assembled each genome (~712 Mb) into 24 chromosomes with a contig N50 of 19.46-29.71 Mb and an integration efficiency of 88.13-92.80%. Furthermore, 26,851-28,133 protein-coding genes were predicted. The reference genomes of three geographical stocks of L. crocea provide vital resources for future research on the conservation and utilization of genetic diversity.

The whole chromosome-level genome provides resources and insights into the endangered fish Percocypris pingi evolution and conservation

BACKGROUND

Percocypris pingi (Tchang) was classified as Endangered on the Red List of China's Vertebrates in 2015 and is widely distributed in the Upper Yangtze River. Although breeding and release into wild habitats have been performed for this commercially important fish in recent years, low genetic diversity has been found in wild populations. Genomic resources are strongly recommended before formulating and carrying out conservation strategies for P. pingi. Thus, there is an urgent need to conserve germplasm resources and improve the population diversity of P. pingi. To date, the whole genome of P. pingi has not been reported.

RESULTS

In our study, we constructed the first chromosome-level genome of P. pingi by high-throughput chromosome conformation capture (Hi-C) technology and PacBio long-read sequencing. The assembled genome was 1.7 Gb in size, with an N50 of 17,692 bp and a GC content from circular consensus sequencing of 37.67%. The Hi-C results again demonstrated that P. pingi was tetraploid (n = 98), with the genome consisting of 24-type and 25-type chromosomes. Chr.19 of the 24-type chromosomes in P. pingi resulted from the fusion of chr.19 and chr.22 in zebrafish. The divergence times between 24-type and 25-type chromosomes was around 6.1 million years ago. A total of 25,198 and 25,291 protein-coding genes were obtained from the 24-type and 25-type chromosomes, respectively. The ploidy of P. pingi is an allotetraploid. A total of 8,741 genes of P. pingi were clustered into 4,378 gene families that were shared with 14 other species, and the P. pingi genome had 68 unique gene families. Phylogenetic analyses indicated that P. pingi was most closely related to Schizothorax oconnori, and the genes were clustered on one branch. We identified 166 significantly expanded gene families and 173 significantly contracted gene families in P. pingi. The most enriched positive protein-coding genes, such as Bmp-4, Etfdh, homeobox protein HB9, and ATG3, were screened.

CONCLUSION

Our study provides a high-quality chromosome-anchored reference genome for P. pingi and provides sufficient information on the chromosomes, which will lead to valuable resources for genetic, genomic, and biological studies of P. pingi and for improving the genetic diversity, population size, and scientific conservation of endangered fish and other key cyprinid species in aquaculture.

Single-cell transcriptomes reveal spermatogonial stem cells and the dynamic heterogeneity of spermatogenesis in a seasonal breeding teleost

Seasonal spermatogenesis in fish is driven by spermatogonial stem cells (SSCs), which undergo a complex cellular process to differentiate into mature sperm. In this study, we characterized spermatogenesis in the large yellow croaker (Larimichthys crocea), a marine fish of significant commercial value, based on a high-resolution single-cell RNA-sequencing atlas of testicular cells from three distinct developmental stages: juvenile, adult differentiating and regressed testes. We detailed a continuous developmental trajectory of spermatogenic cells, from spermatogonia to spermatids, elucidating the molecular events involved in spermatogenesis. We uncovered dynamic heterogeneity in cellular compositions throughout the annual reproductive cycle, accompanied by strong molecular signatures within specific testicular cells. Notably, we identified a distinct population of SSCs and observed a critical metabolic transition from glycolysis to oxidative phosphorylation, enhancing our understanding of the biochemical and molecular characteristics of SSCs. Additionally, we elucidated the interactions between somatic cells and spermatogonia, illuminating the mechanisms that regulate SSC development. Overall, this work enhances our understanding of spermatogenesis in seasonal breeding teleosts and provides essential insights for the further conservation and culture of SSCs.

Lipid fingerprinting and geographical origin identification of large yellow croaker (Larimichthys crocea) using UHPLC-QTOF-MS

The Large Yellow Croaker (LYC) with rich lipids offers numerous health benefits, yet its lipid profile remains underexplored. Therefore, the lipids of LYC were comprehensively profiled and compared based lipidomics. Higher nutritional value of LYC was identified based on lipid and fatty acid assessment, particularly ovary and brain. A total of 963 lipid species belonging to 47 lipid subclasses were identified, exhibiting higher levels of glycerophospholipids (GPs) and sphingolipids (SPs) in the brain, as well as high levels of glycolipids (GLs) in the muscle. Furthermore, unique lipid subclasses were detected in the brain (Acylcarnitine (CAR), Alpha-hydroxy-N-stearoyl phytosphingosine (Cer-AP)) and ovary (N-arachidonoyl glycine (NAGly)). Notably, 8 lipid subclasses were selected as potential contributors for four tissues differentiation. Additionally, identification of LYC from various origins was achieved through lipidomics for the first time. LYC from Zhejiang and Fujian provinces could be distinguished by 100 lipid biomarkers. Most different lipids significantly negatively correlated with seawater pH and dissolved oxygen, but positively correlated with netting density, depth, environmental temperature and salinity. This comprehensive analysis provided valuable insights into LYC's nutritional values and origin identified.

Advances in biotechnology and breeding innovations in China's marine aquaculture

Biotechnology is the key driving force behind the sustainable development of aquaculture, as biological innovation would significantly improve the capabilities of aquatic breeding and achieve independent and controllable seeding sources to ensure food safety. In this article, we have analyzed the current status and existing problems of marine aquaculture in China. Based on these data, we have summarized the recent (especially the last 10 years) biotechnological innovation and breeding progress of marine aquaculture in China, including whole genome sequencing, sex-related marker screening, genomic selection, and genome editing, as well as progress of improved marine fish varieties in China. Finally, the perspectives in this field have been discussed, and three future countermeasures have been proposed.

Decoding the fish genome opens a new era in important trait research and molecular breeding in China

Aquaculture represents the fastest-growing global food production sector, as it has become an essential component of the global food supply. China has the world's largest aquaculture industry in terms of production volume. However, the sustainable development of fish culture is hindered by several concerns, including germplasm degradation and disease outbreaks. The practice of genomic breeding, which relies heavily on genome information and genotypephenotype relationships, has significant potential for increasing the efficiency of aquaculture production. In 2014, the completion of the genome sequencing and annotation of the Chinese tongue sole signified the beginning of the fish genomics era in China. Since then, domestic researchers have made dramatic progress in functional genomic studies. To date, the genomes of more than 60 species of fish in China have been assembled and annotated. Based on these reference genomes, evolutionary, comparative, and functional genomic studies have revolutionized our understanding of a wide range of biologically and economically important traits of fishes, including growth and development, sex determination, disease resistance, metamorphosis, and pigmentation. Furthermore, genomic tools and breeding techniques such as SNP arrays, genomic selection, and genome editing have greatly accelerated genetic improvement through the incorporation of functional genomic information into breeding activities. This review aims to summarize the current status, advances, and perspectives of the genome resources, genomic study of important traits, and genomic breeding techniques of fish in China. The review will provide aquaculture researchers, fish breeders, and farmers with updated information concerning fish genomic research and breeding technology. The summary will help to promote the genetic improvement of production traits and thus will support the sustainable development of fish aquaculture.

Larimichthys crocea (large yellow croaker): A bibliometric study

Larimichthys crocea is an important economic fish of East Asia, and numerous studies have been conducted on its breeding, aquaculture, preservation and processing; however, there is no systematic review of the literature on the research of Larimichthys crocea. Derwent Data Analyzer (DDA) was used to analyze 1192 Larimichthys crocea research papers indexed by SCI-E, CSCD and KCI from 2001 to 2023. The number of research publications on Larimichthys crocea has rapidly increased, and institutions and scholars from China, the United States, South Korea, Japan, and Norway have conducted the majority of Larimichthys crocea research. The immune response, Pseudomonas plecoglossicida, gene expression, lipid immune response, transcriptomics and other areas have attracted the most attention. To increase the immunity and disease resistance of Larimichthys crocea and improve its survival, growth, storage and transport, researchers have carried out a large amount of research, which has promoted not only the culture of Larimichthys crocea but also the restoration of wild Larimichthys crocea and the rehabilitation of the ecological environment.

Comparative genomics of dusky kob (Argyrosomus japonicus, Sciaenidae) conspecifics: Evidence for speciation and the genetic mechanisms underlying traits

Dusky kob (Argyrosomus japonicus) is a commercially important finfish, indigenous to South Africa, Australia, and China. Previous studies highlighted differences in genetic composition, life history, and morphology of the species across geographic regions. A draft genome sequence of 0.742 Gb (N50 = 5.49 Mb; BUSCO completeness = 97.8%) and 22,438 predicted protein-coding genes was generated for the South African (SA) conspecific. A comparison with the Chinese (CN) conspecific revealed a core set of 32,068 orthologous protein clusters across both genomes. The SA genome exhibited 440 unique clusters compared to 1928 unique clusters in the CN genome. Transportation and immune response processes were overrepresented among the SA accessory genome, whereas the CN accessory genome was enriched for immune response, DNA transposition, and sensory detection (FDR-adjusted p < 0.01). These unique clusters may represent an adaptive component of the species' pangenome that could explain population divergence due to differential environmental specialisation. Furthermore, 700 single-copy orthologues (SCOs) displayed evidence of positive selection between the SA and CN genomes, and globally these genomes shared only 92% similarity, suggesting they might be distinct species. These genes primarily play roles in metabolism and digestion, illustrating the evolutionary pathways that differentiate the species. Understanding these genomic mechanisms underlying adaptation and evolution within and between species provides valuable insights into growth and maturation of kob, traits that are particularly relevant to commercial aquaculture.

Draft Genome of Akame (Lates Japonicus) Reveals Possible Genetic Mechanisms for Long-Term Persistence and Adaptive Evolution with Low Genetic Diversity

It is known that some endangered species have persisted for thousands of years despite their very small effective population sizes and low levels of genetic polymorphisms. To understand the genetic mechanisms of long-term persistence in threatened species, we determined the whole genome sequences of akame (Lates japonicus), which has survived for a long time with extremely low genetic variations. Genome-wide heterozygosity in akame was estimated to be 3.3 to 3.4 × 10-4/bp, one of the smallest values in teleost fishes. Analysis of demographic history revealed that the effective population size in akame was around 1,000 from 30,000 years ago to the recent past. The relatively high ratio of nonsynonymous to synonymous heterozygosity in akame indicated an increased genetic load. However, a detailed analysis of genetic diversity in the akame genome revealed that multiple genomic regions, including genes involved in immunity, synaptic development, and olfactory sensory systems, have retained relatively high nucleotide polymorphisms. This implies that the akame genome has preserved the functional genetic variations by balancing selection, to avoid a reduction in viability and loss of adaptive potential. Analysis of synonymous and nonsynonymous nucleotide substitution rates has detected signs of positive selection in many akame genes, suggesting adaptive evolution to temperate waters after the speciation of akame and its close relative, barramundi (Lates calcarifer). Our results indicate that the functional genetic diversity likely contributed to the long-term persistence of this species by avoiding the harmful effects of the population size reduction.

Stability of population genetic structure in large yellow croaker (Larimichthys crocea): Insights from temporal, geographical factors, and artificial restocking processes

Despite concerns about overfishing and the potential impact of release programs on wild populations, our study of 3116 individuals from 13 wild populations and 2787 individuals from two cultured populations in Zhejiang and Fujian provinces spanning 2008 to 2023 reveals a relatively stable genetic diversity in Larimichthys crocea. Surprisingly, the genetic diversity of wild large yellow croaker populations has remained consistent over the years, suggesting minimal influence from population declines due to overfishing. With the exception of populations in Sansha Bay and Luoyuan Bay, no significant genetic differences were observed among wild populations, indicating a single panmictic genetic population across the East and South China seas. Notably, significant genetic differentiation exists between cultured and wild populations, suggesting a possible limited genetic adaptation of cultured-released individuals to the wild environment. The genetic differences observed between the Sansha Bay, with its adjacent Luoyuan Bay populations, and other wild populations underscore the dual effects of habitat environment and farming activities on the genetic structure of large yellow croaker. Our findings suggest that, despite declines in population numbers due to overfishing and expands extensive cultured releases, the genetic diversity of L. crocea populations remains largely unaffected. Moreover, the L. crocea population along the Chinese coast appears to form a single panmictic population with considerable genetic diversity.

Effects of Gene Alternative Splicing Events on Resistance to Cryptocaryonosis of Large Yellow Croaker (Larimichthys crocea)

Large yellow croaker (L. crocea) is a productive species in marine aquaculture with great economic value in China. However, the sustainable development of large yellow croaker is hampered by various diseases including cryptocaryonosis caused by Cryptocaryon irritans. The genetic regulation processes for cryptocaryonosis in large yellow croaker are still unclear. In this present study, we analyzed differential alternative splicing events between a C. irritans resistance strain (RS) and a commercial strain (CS). We identified 678 differential alternative splicing (DAS) events from 453 genes in RS and 719 DAS events from 500 genes in CS. A set of genes that are specifically alternatively spliced in RS was identified including mfap5, emp1, and trim33. Further pathway analysis revealed that the specifically alternative spliced genes in RS were involved in innate immune responses through the PRR pathway and the Toll and Imd pathway, suggesting their important roles in the genetic regulation processes for cryptocaryonosis in large yellow croaker. This study would be helpful for the studies of the pathogenesis of cryptocaryonosis and dissection of C. irritans resistance for L. crocea.

Young duplicate genes show developmental stage- and cell type-specific expression and function in Caenorhabditis elegans

Gene duplication produces the material that fuels evolutionary innovation. The "out-of-testis" hypothesis suggests that sperm competition creates selective pressure encouraging the emergence of new genes in male germline, but the somatic expression and function of the newly evolved genes are not well understood. We systematically mapped the expression of young duplicate genes throughout development in Caenorhabditis elegans using both whole-organism and single-cell transcriptomic data. Based on the expression dynamics across developmental stages, young duplicate genes fall into three clusters that are preferentially expressed in early embryos, mid-stage embryos, and late-stage larvae. Early embryonic genes are involved in protein degradation and develop essentiality comparable to the genomic average. In mid-to-late embryos and L4-stage larvae, young genes are enriched in intestine, epidermal cells, coelomocytes, and amphid chemosensory neurons. Their molecular functions and inducible expression indicate potential roles in innate immune response and chemosensory perceptions, which may contribute to adaptation outside of the sperm.

Vitamin D3 can effectively and rapidly clear largemouth bass ranavirus by immunoregulation

Largemouth bass ranavirus (LMBV) is a highly destructive pathogen that causes significant mortality rates among largemouth bass populations. Unfortunately, there is a dearth of drug development efforts specifically aimed at treating LMBV. To address this, our study sought to investigate the potential effectiveness of incorporating varying doses of VD3 into the diet as a treatment for LMBV. Through qRT-PCR and semi-qPCR, we observed significant suppression and clearance of LMBV pathogens in largemouth bass fed with 15000 IU/Kg and 20000 IU/Kg of VD3 within 14 days. In addition, VD3 treatment significantly increased the expression levels of key immune-related genes such as IL-1β, IFN-γ, Mx, and IgM. Encouragingly, we observed that VD3 significantly increased antioxidant and immune activities such as TSOD, TAOC and C3 in serum and maintained total protein levels. Additionally, tissue pathology sections highlighted a dose-dependent relationship between VD3 supplementation and tissue damage, with the 15000 IU and 20000 IU groups exhibiting minimal damage. In conclusion, a reasonable concentration of VD3 effectively reduced LMBV replication and tissue damages, while improved immune-related genes expression and serum biochemical indices. These findings declare the considerable therapeutic potential of VD3 supplementation for combating LMBV disease and provide an alternative treatment option for fish farming.

Genome-Wide RADseq Reveals Genetic Differentiation of Wild and Cultured Populations of Large Yellow Croaker

Larimichthys crocea (also known as the large yellow croaker) is one of the most economically important marine fishes in China, and research on the ecology and genetics of this species is of immense significance. In this study, we performed restriction site-associated DNA sequencing (RAD-seq) of 54 individuals collected from four sites in China to analyze the genetic structure and diversity of large yellow croaker at the genome level. It revealed that the large yellow croaker populations in the Ningde and Zhoushan coastal waters can be clearly distinguished. Different genetic diversity indices were used to analyze the genetic diversity of the large yellow croaker, which showed that there was a differentiation trend between the wild and farmed populations in Ningde. Moreover, we identified genetically differentiated genomic regions between the populations. GO gene enrichment analysis identified genes that are related to fatty acid metabolism and growth. These findings enhance our understanding of genetic differentiation and adaptation to different living environments, providing a theoretical basis for the preservation and restoration of the genetic resources of the large yellow croaker.

Different cold tolerances among three strains of large yellow croaker: related to antioxidant defense and energy metabolism

The aim of this study was to compare low-temperature tolerances in different strains of large yellow croaker. Dai Qu (DQ), Min-Yue Dong (MY), and Quan Zhou (NZ) strains of large yellow croaker were subjected to cold stress (8.6 °C) for 12 h, 24 h, 48 h, and 96 h. Survival rate, histological observation, and antioxidant and energy metabolism indicators were determined. The results showed that compared with the DQ group and MY group, NZ group aggravated hepatic structure, enhanced ROS, lactate, and anaerobic metabolism (PK gene expression and activity), while inhibited ATP, GSH, antioxidant enzymes (mRNA levels and activities of SOD, GPx, and CAT), and aerobic metabolism enzymes (mRNA levels and activities of F-ATPase, SDH, and MDH), indicating the reduction of cold tolerance in the NZ group was closely correlated with the decrement of antioxidative capacity and energy metabolism efficiency. Nrf2 and AMPK gene expressions were correlated with antioxidant and energy metabolism mRNA levels, respectively, suggesting Nrf2 and AMPK might participate in the modulation of target genes during the cold-stress adaptation. In conclusion, low temperature tolerance of fish depended on the antioxidant defense and energy metabolism efficiency, which contributes to understanding the underlying mechanisms of cold adaptation in large yellow croaker.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Genome assembly and isoform analysis of a highly heterozygous New Zealand fisheries species, the tarakihi (Nemadactylus macropterus)

Although being some of the most valuable and heavily exploited wild organisms, few fisheries species have been studied at the whole-genome level. This is especially the case in New Zealand, where genomics resources are urgently needed to assist fisheries management. Here, we generated 55 Gb of short Illumina reads (92× coverage) and 73 Gb of long Nanopore reads (122×) to produce the first genome assembly of the marine teleost tarakihi [Nemadactylus macropterus (Forster, 1801)], a highly valuable fisheries species in New Zealand. An additional 300 Mb of Iso-Seq reads were obtained to assist in gene annotation. The final genome assembly was 568 Mb long with an N50 of 3.37 Mb. The genome completeness was high, with 97.8% of complete Actinopterygii Benchmarking Universal Single-Copy Orthologs. Heterozygosity values estimated through k-mer counting (1.00%) and bi-allelic SNPs (0.64%) were high compared with the same values reported for other fishes. Iso-Seq analysis recovered 91,313 unique transcripts from 15,515 genes (mean ratio of 5.89 transcripts per gene), and the most common alternative splicing event was intron retention. This highly contiguous genome assembly and the isoform-resolved transcriptome will provide a useful resource to assist the study of population genomics and comparative eco-evolutionary studies in teleosts and related organisms.

Linked selection, differential introgression and recombination rate variation promote heterogeneous divergence in a pair of yellow croakers

Understanding the mechanisms underlying heterogeneous genomic divergence is of particular interest in evolutionary biology. Highly differentiated genomic regions, known as genomic islands, often evolve between diverging lineages. These genomic islands may be related to selection promoting adaptation or reproductive isolation. Based on whole genome assembly and genome-wide RAD sequencing in a pair of yellow croakers (genus: Larimichthys), we investigated the evolutionary processes shaping genomic landscapes of divergence. Demographic modelling indicated that the two species diverged following a secondary contact scenario, where differential introgression and linked selection were suggested to be involved in heterogeneous genomic divergence. We identified reduced recombination rate in genomic islands and a relatively good conservation of both genetic diversity and recombination landscapes between species, which highlight the roles of linked selection and recombination rate variation in promoting heterogeneous divergence in the common ancestral lineage of the two species. In addition, we found a positive correlation between differentiation (F_ST ) and absolute sequence divergence (D_xy ), and elevated D_xy in genomic islands, indicating that the genomic landscape of divergence was not shaped by linked selection alone. Restricted gene flow in highly differentiated regions has probably remodelled the landscape of heterogeneous genomic divergence. This study highlights that highly differentiated genomic regions can also arise from a combination of linked selection and differential gene flow in interaction with varying recombination rates.

Integrative analysis of GWAS and transcriptome reveals p53 signaling pathway mediates resistance to visceral white-nodules disease in large yellow croaker

Visceral white-nodules disease (VWND), caused by Pseudomonas plecoglossicida, is one of the primary causes of morbidity and mortality in large yellow croaker aquaculture. Host disease resistance is a heritable trait that involves complex regulatory processes. However, the regulatory mechanism of bacterial resistance in large yellow croaker is still unclear. This study attempted to systematically evaluate the major genetic loci and transcriptional regulatory mechanisms associated with the resistance to VWND in large yellow croaker by crossover method studies. A large population of large yellow croaker was challenged with P. plecoglossicida, with survival time recorded and samples were taken for genotyping. Meanwhile, spleen samples that were used for RNA-seq to compare their transcriptomic profiles before and after infection were taken from resistant populations (RS) and susceptible control populations (CS) bred using the genomic selection (GS) technique. Genome-wide association analyses using 46 K imputed SNP genotypes highlighted that resistance is a polygenic trait. The integrative analysis results show the co-localization of the cd82a gene between disease resistance-related genetic loci and comparative transcriptional analysis. And functional enrichment analysis showed differential enrichment of the p53 signaling pathway in RS and CS groups, suggesting that there may be cd82a-mediated p53 signaling pathway activation for VWND resistance. This large-scale study provides further evidence for the heritability and transcriptional regulatory mechanisms of host inheritance of VWND resistance.

Transcriptome characterization and SSR discovery in the giant spiny frog Quasipaa spinosa

The giant spiny frog Quasipaa spinosa (Amphibia: Ranidae) is a large unique frog species found mainly in southern China with a low amount of fat and high protein, and it has become one of the most important aquaculture animal species in China. To better understand its genetic background and screen potential molecular markers for artificial breeding and species conservation, we constructed an expression profile of Q. spinosa with high-throughput RNA sequencing and acquired potential SSR markers. Approximately 81.7 Gb of data and 93,887 unigenes were generated. The transcriptome contains 2085 (80.7 %) complete BUSCOs, suggesting that our assembly methods were effective and accurate.These unigenes were functionally classified using 7 functional databases, yielding 17,482 Pfam-, 12,752 Sting-, 17,526 KEGG-, 24,341 Swiss-Prot-, 28,604 Nr-, 16,287 GO- and 12,752 COG-annotated unigenes. Among several amphibian species, Q. spinosa unigenes had the highest number of hits to Xenopus tropicalis (35.25 %), followed by Xenopus laevis (12.68 %). 1417 unigenes were assigned to the immune system. In addition, a total of 33,019 candidate SSR markers were identified from the constructed library. Further tests with 20 loci and 118 large-scale breeding specimens gathered from four culture farms in China showed that 15 (75 %) loci were polymorphic, with the number of alleles per locus varying from 3 to 9 (mean of 4.3). The PIC values for the SSR markers ranged from 0.19 to 0.82, with an average value of 0.43, indicating moderate polymorphism in Q. spinosa. The transcriptomic profile and SSR repertoire obtained in the present study will facilitate population genetic studies and the selective breeding of amphibian species.

Aquaculture omics: An update on the current status of research and data analysis

Aquaculture is the fast-growing agricultural sector and has the ability to meet the growing demand for protein nutritional security for future population. In future aquaculture is going to be the major source of fish proteins as capture fisheries reached at its maximum. However, several challenges need to overcome such as lack of genetically improved strains/varieties, lack of species-specific feed/functional feed, round the year availability of quality fish seed, pollution of ecosystems and increased frequencies of disease occurrence etc. In recent years, the continuous development of high throughput sequencing technology has revolutionized the biological sciences and provided necessary tools. Application of 'omics' in aquaculture research have been successfully used to resolve several productive and reproductive issues and thus ensure its sustainability and profitability. To date, high quality draft genomes of over fifty fish species have been generated and successfully used to develop large number of single nucleotide polymorphism markers (SNPs), marker panels and other genomic resources etc in several aquaculture species. Similarly, transcriptome profiling and miRNAs analysis have been used in aquaculture research to identify key transcripts and expression analysis of candidate genes/miRNAs involved in reproduction, immunity, growth, development, stress toxicology and disease. Metagenome analysis emerged as a promising scientific tool to analyze the complex genomes contained within microbial communities. Metagenomics has been successfully used in the aquaculture sector to identify novel and potential pathogens, antibiotic resistance genes, microbial roles in microcosms, microbial communities forming biofloc, probiotics etc. In the current review, we discussed application of high-throughput technologies (NGS) in the aquaculture sector.

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.

Molecular characterization, expression and functional analysis of large yellow croaker (Larimichthys crocea) peroxisome proliferator-activated receptor gamma

The peroxisome proliferator-activated receptor gamma (PPARγ) are nuclear receptors with distinct roles in energy metabolism and immunity. Although extensively studied in mammals, immunomodulatory roles of this molecule in teleost fish remain to be investigated. In this study, large yellow croaker (Larimichthys crocea) PPARγ (LcPPARγ) sequence was cloned, which encodes a polypeptide of 541 amino acids that include signature domains belonging to the nuclear receptor superfamily. Phylogenetically, LcPPARγ was most closely related to PPARγ derived from European sea bass (Dicentrarchus labrax). Quantitative analysis shown a ubiquitous expression of this molecule, with highest expression level detected in the intestine. The expression of LcPPARγ was decreased in the intestine, muscle, body kidney, spleen and head kidney-derived monocytes/macrophages (MO/MФs) over the course of Vibrio alginolyticus (V. alginolyticus) infection. In contrast, an up-regulation of LcPPARγ was observed in head kidney-derived MO/MФs following docosahexaenoic acid (DHA) treatment. This increase in LcPPARγ leads to an up-regulation of LcCD11b and LcCD18 and an enhancement of complement-mediated phagocytosis. Furthermore, cytokine secretions of V. alginolyticus-stimulated MO/MФs were modulated following LcPPARγ activations that up-regulated the expression of LcIL-10, while decreased the expression of LcIL-1β, LcTNF-α and LcTGF-β1. Overall, our results indicated that LcPPARγ plays a role in regulating functions of MO/MФs and likely contribute to MO/MФs polarization.

Linkage mapping, comparative genome analysis, and QTL detection for growth in a non-model teleost, the meagre Argyrosomus regius, using ddRAD sequencing

Meagre (Argyrosomus regius), is a benthopelagic species rapidly emerging in aquaculture, due to its low food to biomass conversion rate, good fillet yield and ease of production. Tracing a species genomic background along with describing the genetic basis of important traits can greatly influence both conservation strategies and production perspectives. In this study, we employed ddRAD sequencing of 266 fish from six F1 meagre families, to construct a high-density genetic map comprising 4529 polymorphic SNP markers. The QTL mapping analysis provided a genomic appreciation for the weight trait identifying a statistically significant QTL on linkage group 15 (LG15). The comparative genomics analysis with six teleost species revealed an evolutionarily conserved karyotype structure. The synteny observed, verified the already well-known fusion events of the three-spine stickleback genome, reinforced the evidence of reduced evolutionary distance of Sciaenids with the Sparidae family, reflected the evolutionary proximity with Dicentrarchus labrax, traced several putative chromosomal rearrangements and a prominent putative fusion event in meagre’s LG17. This study presents novel elements concerning the genome evolutionary history of a non-model teleost species recently adopted in aquaculture, starts to unravel the genetic basis of the species growth-related traits, and provides a high-density genetic map as a tool that can help to further establish meagre as a valuable resource for research and production.

Robust host source tracking building on the divergent and non-stochastic assembly of gut microbiomes in wild and farmed large yellow croaker

BACKGROUND

Given the lack of genetic background, the source tracking unknown individuals of fish species with both farmed and wild populations often cannot be robustly achieved. The gut microbiome, which is shaped by both deterministic and stochastic processes, can serve as a molecular marker of fish host source tracking, particularly as an alternative to the yet-to-be-established host genetic marker. A candidate for testing the feasibility is the large yellow croaker, Larimichthys crocea, which is carnivorous and ranks the top mariculture fish in China. Wild resource of this fish was depleted decades ago and might have potential problematic estimation because of escaping of farmed individuals.

RESULTS

The rectums of wild (n = 212) and farmed (n = 79) croakers from multiple batches were collected for the profiling of their gut bacterial communities. The farmed individuals had a higher alpha diversity and lower bacterial load than the wild individuals. The gut microbiota of the two sources exhibited divergence and high inter-batch variation, as featured by the dominance of Psychrobacter spp. in the wild group. Predicted functional capacity of the gut microbiome and representative isolates showed differences in terms of host source. This difference can be linked to the potential diet divergence between farmed and wild fishes. The non-stochastic distribution pattern of the core gut microbiota of the wild and farmed individuals supports the feasibility of microbiota-based host source tracking via the machine learning algorithm. A random forest classifier based on the divergence and non-stochastic assembly of the gut microbiome was robust in terms of host source tracking the individuals from all batches of croaker, including a newly introduced batch.

CONCLUSIONS

Our study revealed the divergence of gut microbiota and related functional profiles between wild and farmed croakers. For the first time, with representative datasets and non-stochastic patterns, we have verified that gut microbiota can be robustly applied to the tracking of host source even in carnivorous fish. Video abstract.

Genetic and Pathogenic Characterization of a New Iridovirus Isolated from Cage-Cultured Large Yellow Croaker (Larimichthys crocea) in China

Iridoviruses are an important pathogen of ectothermic vertebrates and are considered a significant threat to aquacultural fish production. Recently, one of the most economically important marine species in China, the large yellow croaker (Larimichthys crocea), has been increasingly reported to be the victim of iridovirus disease. In this study, we isolated and identified a novel iridovirus, LYCIV-ZS-2020, from cage-cultured large yellow croaker farms in Zhoushan island, China. Genome sequencing and subsequent phylogenetic analyses showed that LYCIV-ZS-2020 belongs to the genus Megalocytivirus and is closely related to the Pompano iridoviruses isolated in the Dominican Republic. LYCIV-ZS-2020 enriched from selected tissues of naturally infected large yellow croaker was used in an artificial infection trial and the results proved its pathogenicity in large yellow croaker. This is the first systematic research on the genetic and pathogenic characterization of iridovirus in large yellow croakers, which expanded our knowledge of the iridovirus.

Whole-Genome Sequencing and Genome-Wide Studies of Spiny Head Croaker (Collichthys lucidus) Reveals Potential Insights for Well-Developed Otoliths in the Family Sciaenidae

Spiny head croaker (Collichthys lucidus), belonging to the family Sciaenidae, is a small economic fish with a main distribution in the coastal waters of Northwestern Pacific. Here, we constructed a nonredundant chromosome-level genome assembly of spiny head croaker and also made genome-wide investigations on genome evolution and gene families related to otolith development. A primary genome assembly of 811.23 Mb, with a contig N50 of 74.92 kb, was generated by a combination of 49.12-Gb Illumina clean reads and 35.24 Gb of PacBio long reads. Contigs of this draft assembly were further anchored into chromosomes by integration with additional 185.33-Gb Hi-C data, resulting in a high-quality chromosome-level genome assembly of 817.24 Mb, with an improved scaffold N50 of 26.58 Mb. Based on our phylogenetic analysis, we observed that C. lucidus is much closer to Larimichthys crocea than Miichthys miiuy. We also predicted that many gene families were significantly expanded (p-value <0.05) in spiny head croaker; among them, some are associated with "calcium signaling pathway" and potential "inner ear functions." In addition, we identified some otolith-related genes (such as otol1a that encodes Otolin-1a) with critical deletions or mutations, suggesting possible molecular mechanisms for well-developed otoliths in the family Sciaenidae.

Chromosome-scale assembly and high-density genetic map of the yellow drum, Nibea albiflora

The yellow drum (Nibea albiflora) is an economically important sciaenid fish in East Asian countries. In this study, we sequenced and assembled a near-complete gynogenetic yellow drum genome. We generated 45.63 Gb of Illumina short-reads and 80.27 Gb of PacBio long-reads and assembled them into a 628.01-Mb genome with a contig N50 of 4.42 Mb. Twenty-four chromosomes with a scaffold N50 of 26.73 Mb were obtained using the Hi-C analysis. We predicted a set of 27,069 protein-coding genes, of which 1,581 and 2,583 were expanded and contracted gene families, respectively. The most expanded genes were categorised into the protein binding, zinc-ion binding and ATP binding functional pathways. We built a high-density genetic linkage map that spanned 4,300.2 cM with 24 linkage groups and a resolution of 0.69 cM. The high-quality reference genome and annotated profiles that we produced will not only increase our understanding of the genetic architecture of economic traits in the yellow drum, but also help us explore the evolution and unique biological characteristics of sciaenid fishes.

Genome-Wide Survey Reveals the Microsatellite Characteristics and Phylogenetic Relationships of Harpadon nehereus

Harpadon nehereus forms one of the most important commercial fisheries along the Bay of Bengal and the southeast coast of China. In this study, the genome-wide survey dataset first produced using next-generation sequencing (NGS) was used to provide general information on the genome size, heterozygosity and repeat sequence ratio of H. nehereus. About 68.74 GB of high-quality sequence data were obtained in total and the genome size was estimated to be 1315 Mb with the 17-mer frequency distribution. The sequence repeat ratio and heterozygosity were calculated to be 52.49% and 0.67%, respectively. A total of 1,027,651 microsatellite motifs were identified and dinucleotide repeat was the most dominant simple sequence repeat (SSR) motif with a frequency of 54.35%. As a by-product of whole genome sequencing, the mitochondrial genome is a powerful tool to investigate the evolutionary relationships between H. nehereus and its relatives. The maximum likelihood (ML) phylogenetic tree was constructed according to the concatenated matrix of amino acids translated from the 13 protein-coding genes (PCGs). Monophyly of two species of the genus Harpadon was revealed in the present study and they formed a monophyletic clade with Saurida with a high bootstrap value of 100%. The results would help to push back the frontiers of genomics and open the doors of molecular diversity as well as conservation genetics studies on this species.

Hypoxic acclimation improves mitochondrial bioenergetic function in large yellow croaker Larimichthys crocea under Cu stress

The purpose of this study was to investigate how pre-hypoxia exposure affected the mitochondrial structure and bioenergetic function of large yellow croaker in responding to Cu stress. Fish were acclimated to normoxia and 3.0 mg DO L-1 for 48 h, then subjected to 0 and 120 μg Cu L-1 for another 48 h. Hypoxic acclimation did not affect mitochondrial ultrastructure and reactive oxygen species (ROS), but reduced oxidative phosphorylation (OXPHOS) efficiency. Cu exposure impaired mitochondrial ultrastructure, increased ROS generation and inhibited OXPHOS efficiency. Compared with Cu exposure alone, hypoxic acclimation plus Cu exposure reduced ROS production and improved OXPHOS efficiency by enhancing mitochondrial respiratory control ratio, mitochondrial membrane potential, and activities and gene expressions of electron transport chain enzymes. In conclusion, hypoxic acclimation improved the mitochondrial energy metabolism of large yellow croaker under Cu stress, facilitating our understanding of the molecular mechanisms regarding adaptive responses of hypoxia-acclimated fish under Cu stress.

Whole-genome resequencing of large yellow croaker (Larimichthys crocea) reveals the population structure and signatures of environmental adaptation

Large yellow croaker is an economically important fish in China and East Asia. Despite its economic importance, genome-wide adaptions of domesticated large yellow croaker are largely unknown. Here, we performed whole-genome resequencing of 198 individuals of large yellow croaker obtained in the sea or from farmers in Zhoushan or Ningde. Population genomics analyses revealed the genetic population structure of our samples, reflecting the living environment. Each effective population size is estimated to be declining over generations. Moreover, we identified genetically differentiated genomic regions between the sea-captured population in the Zhoushan Sea area and that of the Ningde Sea area or between the sea-captured population and the farmed population in either area. Gene ontology analyses revealed the gene groups under selective sweep for the adaptation to the domesticated environment. All these results suggest that individuals of the large yellow croaker populations show genomic signatures of adaptation to different living environments.

Genomic characteristics and profile of microsatellite primers for Acanthogobius ommaturus by genome survey sequencing

Acanthogobius ommaturus is one of the suitable species to study the genetic mechanism of adaptive evolution, but there are few reports on its genetics. In the present study, the genomic survey sequencing method was used to analyze the genome characters of A. ommaturus. A total of 50.50 G high-quality sequence data were obtained in the present study. From the 19-mer distribution frequency, the estimated genome size was 928.01 Mb. The calculated sequence repeat rate was about 38.31%, the heterozygosity was approximately 0.17%, and the GC% content was approximately 40.88%. Moreover, 475,724 simple sequence repeats (SSRs) were identified. Among them, dinucleotide repeats were the most (53.70% of the total SSRs), followed by tri- (35.36%), hexa- (4.59%), tetra- (4.57%) and penta- (1.77%) nucleotide repeats type. This is the first genome-wide feature of this species to be reported.

Identification and functional characterization of the transcription factor coding Dp1 gene in large yellow croaker Pseudosciaena crocea

The transcription factor Dp1, as a binding partner, often forms a dimerization complex with typical E2F to play a central role in regulating gene expression during G1/S cell cycle progression. In this study, a full-length dp1 cDNA (Pcdp1) was successfully cloned and characterized from the large yellow croaker Pseudosciaena crocea. The nucleotidic sequence of Pcdp1 is 1,427 bp long with an open reading frame (ORF) of 1,239 bp encoding a putative protein of 412 amino acids, a 5′-untranslated region of 116 bp and a 3′-untranslated region of 70 bp. Prediction of protein domains showed that PcDp1 contains a DNA-binding domain (DBD) with a DEF box, a dimerization domain and an acidic region at C terminus with transcription activity. Homology comparisons indicated that PcDp1 shared the highest sequence identity of 98.55% with Oreochromis niloticus dp1, followed by 88.72% identity with Danio rerio dp1 and a relatively low identity of 78.91–80.55% with its mammalian and amphibian counterparts. The mRNA of Pcdp1 showed ubiquitously expression in all analyzed tissues, with the highest level of expression in the body kidney. Moderate expression levels of Pcdp1 was found in several immune-related tissues including the gills, head kidney and liver, indicating that PcDp1 might play an important role in osmotic pressure regulation and immune response of the large yellow croaker. The subcellular localization of PcDp1 revealed that it is mainly distributed in the cytoplasm both in COS-7 and parenchymal cells of the spleen, head kidney and kidney tissues. Furthermore, the recombinant PcDp1 exhibited DNA-binding activity to E2F site in vitro. In conclusion, these results indicated that PcDp1 may participate in immune regulation and provide a foundation for further study of the regulatory mechanism of Dp1 in teleosts.

Molecular Characterization and Evolution Analysis of Two Forms of TLR5 and TLR13 Genes Base on Larimichthys crocea Genome Data

TLRs (Toll-like receptors) are essential in host defense against pathogens. There are two types of TLR5, namely, membrane form of TLR5 (TLR5M) and soluble form of TLR5 (TLR5S), both of which perform a crucial role in flagellin response. TLR13 is a TLR that localizes to endosomes and recognizes nucleic acids released by internal microorganisms, including viruses, bacteria, and fungi. Here, the full-length coding sequence (CDS), protein structure, and immune response and subcellular localization of TLR5 (TLR5S) and TLR13 were characterized in large yellow croaker (Larimichthys crocea). These TLRs share high sequence homology with other ichthyic TLRs, while also having their own characters; qtPCR was determined and the results found that the three genes were constitutively expressed in all examined tissues: TLR5M was highly expressed in the spleen and liver; TLR13 expression was high in the kidney, liver, and spleen. And TLRs were upregulated following stimulation with Vibrio parahaemolyticus in the liver, spleen, and kidney. Immunofluorescence staining revealed that TLR5M were localized in the cytoplasm, while TLR5S and TLR13 were in the endosome. The evolutionary analysis has shown that TLR13 was clustered with TLR11, 19, 20, 21, and 22, while TLR5 and TLR3 were classified into a group; these results suggest that TLRs are vital in the defense of L. crocea against bacterial infection and further increase our understanding of TLR function in innate immunity in teleosts.

Molecular characterization of the interferon regulatory factor (IRF) family and functional analysis of IRF11 in the large yellow croaker (Larimichthys crocea)

Interferon regulatory factors (IRFs) are a family of transcription factors involved in regulating interferon (IFN) responses and immune cell development. A total of 11 IRFs have been identified in teleost fish. Here, a complete repertoire of 11 IRFs (LcIRFs) in the large yellow croaker (Larimichthys crocea) was characterized with the addition of five newly identified members, LcIRF2, LcIRF5, LcIRF6, LcIRF10, and LcIRF11. These five LcIRFs possess a DNA-binding domain (DBD) at the N-terminal that contains five to six conserved tryptophan residues and an IRF-association domain (IAD) or IAD2 at the C-terminal that is responsible for interaction with other IRFs or co-modulators. Phylogenetic analysis showed that the 11 LcIRFs were divided into four clades including the IRF1 subfamily, IRF3 subfamily, IRF4 subfamily, and IRF5 subfamily. These are evolutionarily related to their respective counterparts in other fish species. The 11 LcIRFs were constitutively expressed in all examined tissues, although at different expression levels. Upon polyinosinic: polycytidylic acid (poly (I:C)) stimulation, the expression of all 11 LcIRFs was significantly induced in the head kidney and reached the highest levels at 6 h post-stimulation (except LcIRF4). LcIRF1, LcIRF3, LcIRF7, LcIRF8, and LcIRF10 were more strongly induced by poly (I:C) than the other LcIRFs. Significant induction of all LcIRFs was observed in the spleen, with LcIRF2, LcIRF5, LcIRF6, LcIRF7, LcIRF9, and LcIRF11 reaching their highest levels at 48 h LcIRF3 and LcIRF11 showed a stronger response to poly (I:C) in the spleen than the other LcIRFs. In addition, LcIRF1, LcIRF3, LcIRF7, LcIRF9, LcIRF10, and LcIRF11 were significantly induced by Vibro alginolyticus in both the spleen and the head kidney, with LcIRF1 strongly induced. Thus, LcIRFs exhibited differential inducible expression patterns in response to different stimuli in different tissues, suggesting that LcIRFs have different functions in the regulation of immune responses. Furthermore, overexpression of LcIRF11 activated the promoters of LcIFNc, LcIFNd, and LcIFNh, and differentially induced the expression levels of LcIFNs and IFN-stimulated genes (ISGs). Overexpression of LcIRF11 in epithelioma papulosum cyprinid (EPC) cells inhibited the replication of viral genes after infection of spring viremia of carp virus (SVCV). These data suggested that LcIRF11 may function as a positive regulator in regulating the cellular antiviral response through induction of type I IFN expression. Taken together, the present study reported molecular characterization and expression analysis of 11 IRFs in the large yellow croaker, and investigated the role of LcIRF11 in the antiviral response, which laid a good foundation for further study on the evolution and functional characterization of fish IRFs.

Effect of replacement of dietary fish oil with four vegetable oils on prostaglandin E2 synthetic pathway and expression of inflammatory genes in marine fish Larimichthys crocea

As a lipid mediator with important immune function, prostaglandin E₂ (PGE₂) has been widely studied in mammals, whereas its synthetic pathway and immune function in fish have yet to be fully studied. To investigate the regulation of PGE₂ synthetic pathway and inflammatory genes expression by dietary different oils and the underlying relationship, a 10-week feeding experiment and an immune challenge were carried out in marine fish Larimichthys crocea. Replacement of dietary fish oil (FO) with four vegetable oils (VO), including soybean oil, linseed oil, palm oil, and olive oil, all reduced PGE₂ levels, and the decrease of arachidonic acid (ARA, substrate for PGE₂) could account for this decline. Meanwhile, the expression of PGE₂ synthesis related genes was basically upregulated, which seemed to be a feedback regulation, but it cannot compensate the deficiency of ARA. In addition, mRNA expression of inflammatory genes, including interleukin (IL)-1β, IL-6, tumor necrosis factor (TNF)α and interferon (IFN)γ was all upregulated in four VO groups compared with FO group, which was the opposite of PGE₂ levels. To verify the inflammatory regulation of PGE₂, an immune challenge was conducted, and PGE₂ alleviated LPS-induced expression of inflammatory genes, including IL-6, TNFα and IFNγ, and the similar downregulation of toll-like receptor (TLR) genes expression revealed that TLR signaling pathway participated in the anti-inflammatory regulation of PGE₂. In conclusion, replacement of dietary FO with four VO (lack of ARA) reduced the levels of PGE₂ that could alleviate LPS-induced inflammatory genes expression via TLR signaling pathway, which could be one of the reasons that VO induced inflammation in marine fish.

Development and Evaluation of a High-Throughput Single-Nucleotide Polymorphism Array for Large Yellow Croaker (Larimichthys crocea)

High-density single-nucleotide polymorphism (SNP) genotyping array is an essential tool for genetic analyses of animals and plants. Large yellow croaker (Larimichthys crocea) is one of the most commercially important marine fish species in China. Although plenty of SNPs have been identified in large yellow croaker, no high-throughput genotyping array is available. In this study, a high-throughput SNP array named NingXin-I with 600K SNPs was developed and evaluated. A set of 82 large yellow croakers were collected from different locations of China and re-sequenced. A total of 9.34M SNPs were identified by mapping sequence reads to the large yellow croaker reference genome. About 1.98M candidate SNPs were selected for further analyses by using criteria such as SNP quality score and conversion performance in the final array. Finally, 579.5K SNPs evenly distributed across the large yellow croaker genome with an average spacing of 1.19 kb were proceeded to array production. The performance of NingXin-I array was evaluated in 96 large yellow croaker individuals from five populations, and 83.38% SNPs on the array were polymorphic sites. A further test of the NingXin-I array in five closely related species in Sciaenidae identified 26.68–56.23% polymorphic SNP rate across species. A phylogenetic tree inferred by using the genotype data generated by NingXin-I confirmed the phylogenetic distance of the species in Sciaenidae. The performance of NingXin-I in large yellow croaker and the other species in Sciaenidae suggested high accuracy and broad application. The NingXin-I array should be valuable for quantitative genetic studies, such as genome-wide association studies (GWASs), high-density linkage map construction, haplotype analysis, and genome-based selection.

Genome-Wide Association Study of Body Shape-Related Traits in Large Yellow Croaker (Larimichthys crocea)

Large yellow croaker (Larimichthys crocea) is one of the most important cultured marine fish on the southeast coast of China. Its body shape is important for the aquaculture industry since it affects the behavior such as swimming, ingesting, and evading, as well as customer preference. Due to the greater consumer demand of small head, slender body large yellow croaker, selecting and breeding of slender individuals with the assistance of genetic markers will benefit the industry quickly. In this study, several traits were employed to represent body shape, including body depth/body length (BD/BL), body thickness/body length (BT/BL), caudal peduncle depth/caudal peduncle length (CPDLR), tail length/body length (TL/BL), and body area/head area (BA/HA). Genome-wide association study was conducted with a panmictic population of 280 individuals to identify SNP and genes potentially associated with body shape. A set of 20 SNPs on 12 chromosomes were identified to be significantly associated with body shape-related traits. Besides, 5 SNPs were identified to be suggestive associated with CPDLR and BT/BL. Surrounding these SNPs, we found some body shape-related candidate genes, including fabp1, acrv1, bcor, mstn, bambi, and neo1, which involved in lipid metabolism, TGF-β signaling, and BMP pathway and other important regulatory pathways. These results will be useful for the understanding of the genetic basis of body shape formation and helpful for body shape controlling of large yellow croaker by using marker-assisted selection.

Molecular characterization and functional analysis of tumor necrosis factor receptor-associated factor 2/7 and tumor necrosis factor receptor 1-associated death domain protein from Larimichthys crocea

In the present study, we characterized tumor necrosis factor receptor-associated factor 2/7 (lcTRAF2/7) and TNFR1-associated death domain protein (lcTRADD) in Larimichthys crocea (L. crocea) and examined their expression profiles in tissues of Vibrio-challenged and unchallenged fish. The coding sequences of lcTRAF2, lcTRAF7, and lcTRADD were 1488, 2454, and 744 nucleotides, and they encoded proteins of 495, 344, and 248 amino acids, respectively. The results of phylogenetic analysis revealed that lcTRAF2, lcTRAF7, and lcTRADD were closest to Oplegnathus fasciatus (85%), Xiphophorus maculatus (97%), and Acanthochromis polyacanthus (65%), respectively. Multiple sequence alignment showed that lcTRAF2 and lcTRAF7 were highly conserved with other vertebrate TRAFs in their functional domains; however, lcTRADD was poorly conserved. The results of quantitative real-time polymerase chain reaction analysis indicated that lcTRAF2, lcTRAF7, and lcTRADD were constitutively expressed in the spleen, liver, kidney, heart, brain, gill, bladder, skin, fin, eye, and muscle. After challenging fish with Vibrio parahaemolyticus, the mRNA expression levels of lcTRAF2, lcTRAF7, and lcTRADD were upregulated in liver, spleen, and kidney. Immunofluorescence staining revealed that lcTRAF2 and lcTRADD were cytoplasmic in localization, whereas lcTRAF7 targeted both the cytoplasm and nucleus. In addition, the NF-κB protein level was upregulated after lipopolysaccharide stimulation in lcTRAF2, lcTRAF7, or lcTRADD overexpressing cells. Taken collectively, these results have improved our understanding of the functions of TRAF2, TRAF7, and TRADD in pathogenic infections in teleosts.

Modulation of copper-induced antioxidant defense, Cu transport, and mitophagy by hypoxia in the large yellow croaker (Larimichthys crocea)

This study aimed to investigate the effects of hypoxia on Cu-induced antioxidant defense, Cu transport, and mitophagy in the liver of the large yellow croaker. Fish were exposed to hypoxia (3.0 mg L^-1), Cu (120 μg L^-1), and hypoxia (3.0 mg L^-1) plus Cu (120 μg L^-1) for 48 h. Hypoxia exposure increased antioxidant abilities to maintain cellular redox balance. Although Cu exposure alone improved antioxidant defense, Cu transport, and mitophagy, these stress responses could not completely neutralize Cu toxicity, as reflected by the elevated reactive oxygen species (ROS) and lipid peroxidation (LPO) and hepatic vacuoles. When compared with Cu stress alone, hypoxia increased Cu toxicity by inhibiting antioxidant defense, Cu transport, and mitophagy, leading to the increment of mortality, ROS, and LPO, and the deterioration of histological structure. The adverse effects of hypoxia on Cu-induced metal transport and mitophagy might be involved in metal-responsive element-binding transcription factor-1 (MTF-1) and Forkhead box O-3 (FoxO3) signaling pathways, respectively. Overall, hypoxia reduced antioxidant response, Cu transport, and mitophagy in fish exposed to Cu, which contributes to understanding the molecular mechanisms underlying negative effects of hypoxia on Cu toxicity in fish.

Effects of dietary exogenous xylanase supplementation on growth performance, intestinal health, and carbohydrate metabolism of juvenile large yellow croaker, Larimichthys crocea

An 8-week feeding trial was conducted to evaluate the effects of dietary xylanase supplementation on growth performance, digestive enzyme activity, intestinal morphology parameter, intestinal microbiome diversity, and carbohydrate metabolism for juvenile large yellow croaker (Larimichthys crocea). Four levels of xylanase were added to basal diets (0, 600, 1200, and 1800 U kg^-1). The results indicated that fish fed the 1200 U kg^-1 xylanase diet had higher weight gain than those fed the 0 and 600 U kg^-1 xylanase diet. The highest intestinal folds and microvillous height were observed at fish fed the 1200 U kg^-1 xylanase diet. High-throughput sequencing revealed that the majority of reads derived from the large yellow croaker digesta belonged to members of Proteobacteria followed by Chloroflex, Bacteroidetes, Spirochaetae, and Firmicute. Supplementation of xylanase in diets increased the relative abundance of Bacteroides and Gemmatimonadete. The higher hepatic glucokinase (GK) and glucose-6-phosphate dehydrogenase (G6PD) activities were observed in fish fed the xylanase supplementation diet. Accordingly, dietary xylanase supplementation upgraded the relative expressions of gk and g6pd genes in liver. In conclusion, optimum dietary xylanase supplementation (600-1200 U kg^-1) could improve the growth performance, optimize the intestinal morphology structure and microbiota constitution, and enhance the ability of carbohydrate utilization of juvenile large yellow croaker.