Whole-genome sequencing of leopard coral grouper ( Plectropomus leopardus) and exploration of regulation mechanism of skin color and adaptive evolution

DNA Methylation and Transcriptome Profiling Reveal the Role of the Antioxidant Pathway and Lipid Metabolism in Plectropomus leopardus Skin Color Formation

Leopard coral grouper (Plectropomus leopardus), possessing a distinct red body color, is an important species in commercial markets; however, the high ratio of black individuals under intensive cultivation has limited the commercial value of the species. To dissect the regulatory mechanisms underlying the red skin trait in P. leopardus, gene expression and DNA methylation modifications were compared between red and black skin tissues after astaxanthin treatment. Astaxanthin effectively increased the redness value a* and body weight. Multi-omics analyses revealed the crucial roles of pathways related to antioxidants and lipid metabolism, particularly "Tyrosine metabolism", "Melanogenesis", "Fatty acid metabolism", "Fatty acid elongation", and "Biosynthesis of unsaturated acids", in red skin coloration. A molecular network for the regulation of red skin coloration in P. leopardus was constructed, and pmel, tyr, tyrp1a, tyrp1b, dct, slc24a5, wnt1, acsl4, elovl1, elovl6l.1, elovl6l.2, and elovl7 were identified as key genes. Notably, pmel, acsl4, and elovl7 were negatively regulated by differential DNA methylation. Our results provide new insight into the molecular and epigenetic mechanisms of body color variation, representing a significant step towards breeding for the red skin trait in P. leopardus.

Analysis of early efficacy and immune reconstitution after autologous hematopoietic stem cell transplantation in multiple myeloma

This retrospective study aimed to stress the advantages of autologous hematopoietic stem cell transplantation (auto-HSCT) in treating primary MM. Ninety-four MM patients who underwent initial parallel sequential auto-HSCT were selected. Data on efficacy (efficacy evaluation, renal function and hemoglobin recovery), immune reconstitution (B-cell subsets, immunoglobulin levels, T-cell subsets, NK cells, neutrophil-to-lymphocyte ratio (NLR), lymphocyte-to-monocyte ratio (LMR)) and hematopoietic reconstitution times were collected and analyzed. Whether in all selected patients or in groups R-ISS II-III, there was a notable increase in the proportion of patients achieving in a very good partial response (VGPR) or better (P < 0.001, P = 0.02) and a complete response (CR) or better (P = 0.007, P = 0.014) after transplantation compared to the pre-transplant status. Post-Transplant Immune Reconstitution Analysis (Baseline vs. Pre-Transplant and Pre-Transplant vs. Post-Transplant): The level of CD19 + B cells, CD20 + B cells, CD22 + B cells, CD3 + T cells, IgG and LMR showed the same change trend, that is, it decreased before transplantation (P < 0.001, P < 0.001, P < 0.001, P < 0.001, P<0.007, P < 0.001) and then increased significantly after transplantation(P < 0.001, P < 0.001, P < 0.001, P < 0.001, P < 0.001, P < 0.001). CD3 + CD4 + T cells from 545.97 (342.11,708.60)/µL to 342.93 (168.38, 475.52)/µL (P < 0.001) and then to 251.48 (188.52, 406.98)/µL (P = 0.348); CD3 + CD8 + T cells from 391.36 (242.19, 563.37)/µL to 337.23 (192.54, 505.96)/µL (P = 0.065) and then to 797.96 (514.49, 1198.03)/µL (P < 0.001), so the CD3 + CD4+/CD3 + CD8 + T cell ratio still remained inverted post-transplant. NK cells changed from 309.86 (206.33, 460.96)/µL to 258.31 (160.75, 436.68)/ µL (P = 0.229) and then to 151.08 (108.17, 240.84)/µL (P = 0.007). Auto-HSCT can promote prolonged remission in patients with MM and also overcome some high-risk factors to achieve superior efficacy in group R-ISS II-III. Patients were immunodeficient before transplantation and auto-HSCT facilitated immune reconstitution.

Gap-free telomere-to-telomere haplotype assembly of the tomato hind (Cephalopholis sonnerati)

The tomato hind (Cephalopholis sonnerati) is an emerging economically important grouper in recent years. With the increasing maturity of sequencing technologies and assembly methodologies, a higher quality reference genome has become both accessible and necessary. In this study, we present two telomere-to-telomere (T2T) gap-free haplotype assemblies of the tomato hind with lengths of 1039.53 Mb (YSFRI_Csonn_HA_1.0, N50 43.83 Mb) and 1039.91 Mb (YSFRI_Csonn_HB_1.0, N50 44.09 Mb). Reads from next-generation sequencing, ONT ultra-long sequencing, and PacBio HiFi sequencing exhibited mapping rates exceeding 99.8% when aligned to these two assemblies. Evaluation using Merqury indicated high accuracy for both assemblies, with average quality values of 51.80 and 51.83, respectively. Percentages of 97.9% and 97.8% of complete BUSCOs were achieved, and a total of 23,270 and 23,184 protein-code genes were inferred in each assembly. Moreover, telomere identification, centromere prediction, and repetitive sequence annotation were also successfully performed. These two assemblies provide robust foundation for the genetic analysis and development of molecular genetic breeding technologies in C. sonnerati.

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.

Effect of Dietary Astragalus polysaccharides (APS) on the Growth Performance, Antioxidant Responses, Immunological Parameters, and Intestinal Microbiota of Coral Trout (Plectropomus leopardus)

The potential effects of Astragalus polysaccharides (APS) were evaluated in coral trout (Plectropomus leopardus). Five APS levels (0%, 0.05%, 0.10%, 0.15%, and 0.20%) were added to the diet of coral trout, and a 56-day growth trial (initial weight 18.62 ± 0.05 g) was conducted. Dietary APS enhanced growth performance, with the highest improvement observed in fish fed the 0.15% APS diet. This concentration also enhanced the antioxidant capacity and immunomodulation of the fish by regulating the expression of genes associated with antioxidant enzymes and immune responses. Intestinal microbiota analysis revealed that APS supplementation significantly increased the Chao1 index and relative abundance of beneficial bacteria (Firmicutes and Bacillus). A high level of APS (0.20%) did not provide additional benefits for growth and health compared to a moderate level (0.15%). These findings indicate that an optimal APS dose promotes growth, enhances antioxidant activity, supports immune function, and improves intestinal microbiota in coral trout. Based on a cubic regression analysis of the specific growth rate, the optimal APS level for the maximal growth of coral trout was determined to be 0.1455%.

SGIV VP82 inhibits the interferon response by degradation of IRF3 and IRF7

During virus-host co-evolution, viruses have developed multiple strategies to dampen IFN response and prevent its antiviral activity in host cells. To date, the interactions between host IFN response and the immune evasion strategies exploited by fish iridoviruses still remain largely uncertain. Here, a potential immune evasion protein candidate of Singapore grouper iridovirus (SGIV), VP82 (encoded by SGIV ORF82) was screened and its roles during viral replication were investigated in detail. Firstly, VP82 overexpression dramatically decreased IFN or ISRE promoter activity and the transcription levels of IFN stimulated genes (ISGs) stimulated by grouper cyclic GMP-AMP synthase (EccGAS)/stimulator of interferon genes (EcSTING), TANK-binding kinase 1 (EcTBK1), IFN regulatory factor 3 (EcIRF3)and EcIRF7. Secondly, Co-IP assays indicated that VP82 interacted with EcIRF3 and EcIRF7, but not EcSTING and EcTBK1, which was consistent with the co-localization between VP82 and EcIRF3 or EcIRF7. Furthermore, VP82 promoted the degradation of EcIRF3 and EcIRF7 in a dose-dependent manner via the autophagy pathway. Finally, VP82 overexpression accelerated SGIV replication, evidenced by the increased transcriptions of viral core genes and viral production. Moreover, the antiviral action of EcIRF3 or EcIRF7 was significantly depressed in VP82 overexpressed cells. Together, VP82 was speculated to exert crucial roles for SGIV replication by inhibiting the IFN response via the degradation of IRF3 and IRF7. Our findings provided new insights into understanding the immune evasion strategies utilized by fish iridovirus through IFN regulation.

Karyotypic changes and diversification time in Epinephelidae groupers (Perciformes). Implications on reproductive isolation

Groupers (Epinephelidae and Serranidae) have attracted special attention to fish farming, and their species offer good opportunities for successful hybridizations. Cytogenetic data allow a better understanding of the role of karyotypic diversification in the acquisition of post-zygotic reproductive isolation (RI). Thus, chromosomal analyses were performed on E. striatus (Caribbean Sea), E. coioides and E. tauvina (Indo-Pacific Region), using standard procedures and mapping of six repetitive DNA classes by the in situ hybridization. The three species have 2n=48 chromosomes. The karyotypes of E. coioides and E. striatus are composed only of acrocentric chromosomes (FN=48), while E. tauvina has 8 submetacentric chromosomes (FN=56). Heterochromatin has a preferential centromeric distribution, and the microsatellite repeats are dispersed throughout the chromosomes of all species. The 18S and 5S rDNA sites are unique but show a colocalization arrangement in E. tauvina and E. striatus. The chromosomal organization suggests that the three species still maintain a significant amount of syntenic regions. The range of the karyotype divergence and the RI levels showed low, but goes turn proportionally greater in relation to the divergence time between the parental species. The slow acquisition of postzygotic RI is consistent with the high karyotype homogeneity presented by Epinephelidae family.

Synchronously sexual maturity in hermaphrodite fish as revealed by transcriptome analysis in Plectropomus leopardus

Leopard coral grouper (Plectropomus leopardus) is a type of hermaphrodite fish, but the mechanisms of gonadal development and gametogenesis remain unclear. In the present study, we performed histological observation and transcriptomic analysis during the process of sexual differentiation in P. leopardus. According to the histological results, sexual differentiation was completed at 15 months old, developed synchronously in male and female individuals at 2 years old, and matured synchronously at 3 years old. Comparative transcriptomic analyses showed that the gonadal had differentiated by 15 months old, with enrichment of pathways associated with cell proliferation, transcriptional metabolism, and germline stem cell differentiation. Furthermore, cilium movement and fatty acid anabolism, which are associated with spermatogenesis and oocyte growth, were significantly enriched at 3 years old. In addition, key genes associated with male and female sex differentiation, such as amh, dmrt1, dmrt2a, zp4, sox3, gdf9, and gsdf, were identified by weighted gene co-expression network analysis (WGCNA). Finally, the localization and expression of the key genes amh and sox3 were observed in different cell types within the testes and ovaries, reflecting the development of the testes and ovaries, respectively. All the evidence indicates that P. leopardus is a hermaphrodite and synchronously sexually mature fish. Our study complements the gonadal development patterns of hermaphroditic fish by providing new insights into the molecular mechanisms underlying sexual differentiation and sex change in hermaphroditic groupers.

Supplementation of Extender with Melatonin Improves the Motility, Mitochondrial Membrane Potential, and Fertilization Ability of Cryopreserved Brown-Marbled Grouper Sperm

Sperm cryopreservation is a valuable tool for breeding, conservation, and genetic improvement in aquatic resources, while oxidative damage will cause a decline in sperm quality during this progress. Melatonin (MT), a natural antioxidant hormone, is used as an additive in sperm cryopreservation to reduce cellular damage from oxidative stress. Here, we aimed to investigate the effect of adding MT to the freezing medium in sperm cryopreservation of brown-marbled grouper (Epinephelus fuscoguttatus). Different concentrations of MT (0, 0.1, 0.25, and 0.5 mg/mL) were tested. We evaluated sperm motility, viability, apoptosis, mitochondrial membrane potential (MMP), and fertilization ability to assess the effects of MT supplementation. Our results demonstrated that the addition of MT to the extender improved the post-thaw motility, MMP, and fertilization ability of brown-marbled grouper sperm. The total motility, curvilinear velocity, straight linear velocity, and average path velocity in MT-treated groups (0.1 and 0.25 mg/mL) exhibited significantly higher values than that of the control group. A higher MMP (p < 0.05) was observed in the group treated with 0.25 mg/mL MT, suggesting that supplementation of MT in the extender might be able to protect mitochondrial membrane integrity effectively. Regarding fertilizing ability, 0.25 mg/mL MT yielded a significantly higher hatching rate than the control. An adverse effect was found with the concentration of MT up to 0.5 mg/mL, suggesting the possible toxicity of a high-dose addition. In this study, we optimized the sperm cryopreservation protocol of brown-marbled grouper, which might be valuable for sperm cryopreservation and sample commercialization of groupers and other fish.

Deciphering the population structure and genetic basis of growth traits from whole-genome resequencing of the leopard coral grouper ( Plectropomus leopardus)

The leopard coral grouper ( Plectropomus leopardus) is a species of significant economic importance. Although artificial cultivation of P. leopardus has thrived in recent decades, the advancement of selective breeding has been hindered by the lack of comprehensive population genomic data. In this study, we identified over 8.73 million single nucleotide polymorphisms (SNPs) through whole-genome resequencing of 326 individuals spanning six distinct groups. Furthermore, we categorized 226 individuals with high-coverage sequencing depth (≥14×) into eight clusters based on their genetic profiles and phylogenetic relationships. Notably, four of these clusters exhibited pronounced genetic differentiation compared with the other populations. To identify potentially advantageous loci for P. leopardus, we examined genomic regions exhibiting selective sweeps by analyzing the nucleotide diversity ( θπ) and fixation index ( F ST) in these four clusters. Using these high-coverage resequencing data, we successfully constructed the first haplotype reference panel specific to P. leopardus. This achievement holds promise for enabling high-quality, cost-effective imputation methods. Additionally, we combined low-coverage sequencing data with imputation techniques for a genome-wide association study, aiming to identify candidate SNP loci and genes associated with growth traits. A significant concentration of these genes was observed on chromosome 17, which is primarily involved in skeletal muscle and embryonic development and cell proliferation. Notably, our detailed investigation of growth-related SNPs across the eight clusters revealed that cluster 5 harbored the most promising candidate SNPs, showing potential for genetic selective breeding efforts. These findings provide a robust toolkit and valuable insights into the management of germplasm resources and genome-driven breeding initiatives targeting P. leopardus.

Identification and effective regulation of scarb1 gene involved in pigmentation change in autotetraploid Carassius auratus

The autotetraploid Carassius auratus (4nRR, 4 n=200, RRRR) is derived from whole-genome duplication of Carassius auratus red var. (RCC, 2 n=100, RR). In the current study, we demonstrated that chromatophores and pigment changes directly caused the coloration and variation of 4nRR skin (red in RCC, brownish-yellow in 4nRR). To further explore the molecular mechanisms underlying coloration formation and variation in 4nRR, we performed transcriptome profiling and molecular functional verification in RCC and 4nRR. Results revealed that scarb1, associated with carotenoid metabolism, underwent significant down-regulation in 4nRR. Efficient editing of this candidate pigment gene provided clear evidence of its significant role in RCC coloration. Subsequently, we identified four divergent scarb1 homeologs in 4nRR: two original scarb1 homeologs from RCC and two duplicated ones. Notably, three of these homeologs possessed two highly conserved alleles, exhibiting biased and allele-specific expression in the skin. Remarkably, after precise editing of both the original and duplicated scarb1 homeologs and/or alleles, 4nRR individuals, whether singly or multiply mutated, displayed a transition from brownish-yellow skin to a cyan-gray phenotype. Concurrently, the proportional areas of the cyan-gray regions displayed a gene-dose correlation. These findings illustrate the subfunctionalization of duplicated scarb1, with all scarb1 genes synergistically and equally contributing to the pigmentation of 4nRR. This is the first report concerning the functional differentiation of duplicated homeologs in an autopolyploid fish, substantially enriching our understanding of coloration formation and change within this group of organisms.

Genetic mechanism of body size variation in groupers: Insights from phylotranscriptomics

Animal body size variation is of particular interest in evolutionary biology, but the genetic basis remains largely unknown. Previous studies have shown the presence of two parallel evolutionary genetic clusters within the fish genus Epinephelus with evident divergence in body size, providing an excellent opportunity to investigate the genetic basis of body size variation in vertebrates. Herein, we performed phylotranscriptomic analysis and reconstructed the phylogeny of 13 epinephelids originating from the South China Sea. Two genetic clades with an estimated divergence time of approximately 15.4 million years ago were correlated with large and small body size, respectively. A total of 180 rapidly evolving genes and two positively selected genes were identified between the two groups. Functional enrichment analyses of these candidate genes revealed distinct enrichment categories between the two groups. These pathways and genes may play important roles in body size variation in groupers through complex regulatory networks. Based on our results, we speculate that the ancestors of the two divergent groups of groupers may have adapted to different environments through habitat selection, leading to genetic variations in metabolic patterns, organ development, and lifespan, resulting in body size divergence between the two locally adapted populations. These findings provide important insights into the genetic mechanisms underlying body size variation in groupers and species differentiation.

Dietary antarctic krill improves antioxidant capacity, immunity and reduces lipid accumulation, insights from physiological and transcriptomic analysis of Plectropomus leopardus

BACKGROUND

Due to its enormous biomass, Antarctic krill (Euphausia superba) plays a crucial role in the Antarctic Ocean ecosystem. In recent years, Antarctic krill has found extensive application in aquaculture, emerging as a sustainable source of aquafeed with ideal nutritional profiles. However, a comprehensive study focused on the detailed effects of dietary Antarctic krill on aquaculture animals, especially farmed marine fishes, is yet to be demonstrated.

RESULTS

In this study, a comparative experiment was performed using juvenile P. leopardus, fed with diets supplemented with Antarctic krill (the krill group) or without Antarctic krill (the control group). Histological observation revealed that dietary Antarctic krill could reduce lipid accumulation in the liver while the intestine exhibited no obvious changes. Enzyme activity measurements demonstrated that dietary Antarctic krill had an inhibitory effect on oxidative stress in both the intestine and the liver. By comparative transcriptome analysis, a total of 1,597 and 1,161 differentially expressed genes (DEGs) were identified in the intestine and liver, respectively. Functional analysis of the DEGs showed multiple enriched terms significantly related to cholesterol metabolism, antioxidants, and immunity. Furthermore, the expression profiles of representative DEGs, such as dhcr7, apoa4, sc5d, and scarf1, were validated by qRT-PCR and fluorescence in situ hybridization. Finally, a comparative transcriptome analysis was performed to demonstrate the biased effects of dietary Antarctic krill and astaxanthin on the liver of P. leopardus.

CONCLUSIONS

Our study demonstrated that dietary Antarctic krill could reduce lipid accumulation in the liver of P. leopardus, enhance antioxidant capacities in both the intestine and liver, and exhibit molecular-level improvements in lipid metabolism, immunity, and antioxidants. It will contribute to understanding the protective effects of Antarctic krill in P. leopardus and provide insights into aquaculture nutritional strategies.

Chromosome-level genome assembly and annotation of the yellow grouper, Epinephelus awoara

Epinephelus awoara, as known as yellow grouper, is a significant economic marine fish that has been bred artificially in China. However, the genetic structure and evolutionary history of yellow grouper remains largely unknown. Here, this work presents the high-quality chromosome-level genome assembly of yellow grouper using PacBio single molecule sequencing technique (SMRT) and High-through chromosome conformation capture (Hi-C) technologies. The 984.48 Mb chromosome-level genome of yellow grouper was assembled, with a contig N50 length of 39.77 Mb and scaffold N50 length of 41.39 Mb. Approximately 99.76% of assembled sequences were anchored into 24 pseudo-chromosomes with the assistance of Hi-C reads. Furthermore, approximately 41.17% of the genome was composed of repetitive elements. In total, 24,541 protein-coding genes were predicted, of which 22,509 (91.72%) genes were functionally annotated. The highly accurate, chromosome-level reference genome assembly and annotation are crucial to the understanding of population genetic structure, adaptive evolution and speciation of the yellow grouper.

Integrated bioinformatics approaches to investigate alterations in transcriptomic profiles of monkeypox infected human cell line model

BACKGROUND

The recent re-emergence of the monkeypox (mpox) epidemic in nonendemic regions has raised concerns regarding a potential global outbreak. The mpox virus (MPV) is a smallpox-like virus belonging to the genus Orthopoxvirus (family: Poxviridae). Although studies suggest that MPV infection suppresses the Toll-like receptor-3- and tumor necrosis factor-α-related signaling pathways, whether MPV regulates other immune-related pathways remains unclear.

METHODS

In this study, two distinct temporal patterns were used for establishing an MPV-infected human immortal epithelial cancer cell line (HeLa). These two durations 2 and 12 h of incubation were selected to identify the coregulated genes and pathways affected by MPV infection.

RESULTS

The use of the Gene Ontology framework, Kyoto Encyclopedia of Genes and Genome database, and MetaCore software yielded valuable insights. Specifically, various pathways were found to be enriched in HeLa cells infected with MPV for 2 and 12 h. These pathways included Notch, CD40, CD95, hypoxia-inducible factor-1-α, interleukin (IL)- 1, IL-6, phosphoinositide 3-kinase, nuclear factor-κB, mitogen-activated protein kinase, and oxidative stress-induced signalling pathways. Clusters and pathways of metabolism and viral replication cycles were significantly associated with the 2-hour infection group. This association was identified based on the regulation of genes such as HSPG2, RHPN2, MYL1, ASPHD2, CA9, VIPR1, SNX12, MGC2752, SLC25A1, PEX19, and AREG. Furthermore, clusters and pathways related to immunity and cell movement were found to be associated with the 12-hour infection group. This association was identified based on the regulation of genes such as C1orf21, C19orf48, HRK, IL8, GULP1, SCAND2, ATP5C1, FEZ1, SGSH, TACC2, CYP4X1, MMP1, CPB1, P2RY13, WDR27, PRPF4, and ENDOD1.

CONCLUSIONS

This study can improve our understanding of the mechanisms underlying the pathophysiology and post-infection sequelae of mpox. Our findings provide valuable insights into the various modes of MPV infection.

Transcriptome analysis reveals candidate miRNAs involved in skin color differentiation of juvenile Plectropomus leopardus in response to different background colors

Red skin color in Plectropomus leopardus is important to its ornamental and economic value. However, the color of P. leopardus can change during the rearing process, darkening and turning black due to the influence of environmental background color. The underlying molecular mechanisms that regulate this phenomenon remain unclear. MicroRNAs (miRNAs) are endogenous, small non-coding RNAs that play important roles in numerous biological processes, such as skin differentiation and color formation in many animals. Therefore, we performed miRNA sequencing of P. leopardus skin before (initial) and after rearing with three different background colors (white, black, and blue) using Illumina sequencing to identify candidate miRNAs that may contribute to skin color differentiation. In total, 154,271,376 clean reads were obtained, with over 92 % of them successfully mapped to the P. leopardus reference genome. The miRNA length distributions of all samples displayed peaks around a typical length of 22 nt. Within these sequences, 243 known and 287 novel miRNAs were identified. A total of 65 significantly differentially expressed miRNAs (DEMs) were identified (P < 0.05), including 40 known DEMs and 25 novel DEMs. These DEMs included novel_561, miR-141-3p, and miR-129-5p, whose target genes were primarily associated with pigmentation related processes, including tyrosine metabolism, melanogenesis, and the Wnt signaling pathway. These findings shed light on the potential roles of miRNAs in the darkening of skin color in P. leopardus, thus enhancing our understanding of the molecular mechanisms involved in skin pigmentation differentiation in this species.

Genome-wide association study (GWAS) analysis of black color trait in the leopard coral grouper (Plectropomus leopardus) using whole genome resequencing

The leopard coral grouper (Plectropomus leopardus) is a coral reef fish species that exhibits rapid and diverse color variation. However, the presence of melanoma and the high proportion of individuals displaying black color in artificial breeding have led to reduced economic and ornamental value. To pinpoint single nucleotide polymorphisms (SNPs) and potential genes linked to the black pigmentation characteristic in this particular species, This study gathered a cohort of 360 specimens from diverse origins and conducted a comprehensive genome-wide association analysis (GWAS) employing whole-genome resequencing. As a result, 57 SNPs related to the black skin trait were identified, and a grand total of 158 genes were annotated within 50 kb of these SNPs. Subsequently, GWAS was applied to three populations (LED, QHH, and QHL), and the corresponding results were compared with the analysis results of the total population. The results of the four GWAS models showed significant enrichment in Rap1 signaling pathway, melanin biosynthesis, metabolic pathways, tyrosine metabolism, cAMP signaling pathway, AMPK signaling pathway, PI3K-Akt signaling pathway, EGFR tyrosine kinase inhibitor resistance, HIF-1 signaling pathway, Ras signaling pathway, MAPK signaling pathway, etc. (p < 0.05), which were mainly associated with eleven genes (POL4, MET, E2F2, COMT, ZBED1, TYRP2, FOXP2, THIKA, LORF2, MYH16 and SOX2). Significant differences (p < 0.05) were observed in the expression of all 11 genes in the dorsal skin tissue, in 10 genes except COMT in the ventral skin tissue, and in all 11 genes in the caudal fin tissue. These findings imply that the control of body color in the P. leopardus is the result of the joint action of multiple genes and signaling pathways. These findings will contribute to a more profound comprehension of the genetic attributes that underlie the development of black skin in the vibrant P. leopardus, thus furnishing a theoretical foundation for genetic enhancement.

Potential use of garlic products in ruminant feeding: A review

The addition of antibiotics as growth promoters to ruminant feed can result in bacterial resistance and antibiotic residues in ruminant products. Correspondingly, there is serious public concern regarding the presence of antibiotic residue in ruminant products and the consequent threat to human health. As a result, the addition of plants and their products to ruminant feeds, as an alternative to antibiotics, has received much attention recently. Garlic and its products are rich in organosulphur compounds, which have a variety of biological activities and have been widely used as natural additives in animal production. This review presents recent knowledge on the addition of garlic products (powder, skin, oil, leaf and extracts) to the diets of ruminants. In this paper, garlic products are evaluated with respect to their chemical composition, bioactive compounds, and their impacts on the rumen ecosystem, antioxidant status, immune response, parasitic infection, growth performance and product quality of ruminants. This review provides valuable guidance and a theoretical basis for the development of garlic products as green, highly efficient and safe additives, with the aims of promoting ruminant growth and health, reducing methane emissions and improving ruminant product quality. Garlic extracts have the potential to control parasite infections by decreasing the faecal egg count. Garlic powder, oil and allicin are able to reduce the methane emissions of ruminants. Organosulphur compounds such as allicin, which is present in garlic products, have the potential to inhibit membrane lipid synthesis of the archaeal community, thus influencing the population of methanogenic archaea and resulting in a reduction in methane emissions. Some garlic products are also able to increase the average daily gain (garlic skin, water extract, and leaf) and the feed conversion ratio (garlic skin and leaf) of ruminants. Garlic stalk silage fed to sheep has the potential to improve the nutritional value of mutton by increasing the concentrations of linoleic and linolenic acids and essential amino acids. Sheep fed a diet containing garlic powder or oil are able to produce milk with higher concentrations of the conjugated linoleic acids and n-3 fatty acids, which has health benefits for consumers, due to the widely recognized positive impact of n-3 polyunsaturated fatty acids and conjugated linoleic acids on human heart health, improving platelet aggregation, vasodilation and thrombotic tendency. Overall, garlic products have the potential to enhance growth performance and product quality and reduce parasite infections, as well as methane emissions of ruminants.

Improved chromosomal-level genome assembly and re-annotation of leopard coral grouper

Plectropomus leopardus, as known as leopard coral grouper, is a valuable marine fish that has gradually been bred artificially. To promote future conservation, molecular breeding, and comparative studies, we generated an improved high-quality chromosomal-level genome assembly of leopard coral grouper using Nanopore long-reads, Illumina short reads, and the Hi-C sequencing data. The draft genome is 849.74 Mb with 45 contigs and N50 of 35.59 Mb. Finally, a total of 846.49 Mb corresponding to 99.6% of the contig sequences was anchored to 24 pseudo-chromosomes using Hi-C technology. A final set of 25,965 genes is annotated after manual curation of the predicted gene models, and BUSCO analysis yielded a completeness score of 99.5%. This study significantly improves the utility of the grouper genome and provided a reference for the study of molecular breeding, genomics and biology in this species.

Molecular Mechanisms Underlying Vertebrate Adaptive Evolution: A Systematic Review

Adaptive evolution is a process in which variation that confers an evolutionary advantage in a specific environmental context arises and is propagated through a population. When investigating this process, researchers have mainly focused on describing advantageous phenotypes or putative advantageous genotypes. A recent increase in molecular data accessibility and technological advances has allowed researchers to go beyond description and to make inferences about the mechanisms underlying adaptive evolution. In this systematic review, we discuss articles from 2016 to 2022 that investigated or reviewed the molecular mechanisms underlying adaptive evolution in vertebrates in response to environmental variation. Regulatory elements within the genome and regulatory proteins involved in either gene expression or cellular pathways have been shown to play key roles in adaptive evolution in response to most of the discussed environmental factors. Gene losses were suggested to be associated with an adaptive response in some contexts. Future adaptive evolution research could benefit from more investigations focused on noncoding regions of the genome, gene regulation mechanisms, and gene losses potentially yielding advantageous phenotypes. Investigating how novel advantageous genotypes are conserved could also contribute to our knowledge of adaptive evolution.

Comparative transcriptome analysis of skin color-associated genes in leopard coral grouper (Plectropomus leopardus)

BACKGROUND

The leopard coral grouper (Plectropomus leopardus) is an important economic species in East Asia-Pacific countries. To meet the market demand, leopard coral grouper is facing overfishing and their population is rapidly declining. With the improvement of the artificial propagation technique, the leopard coral grouper has been successfully cultured by Fisheries Research Institute in Taiwan. However, the skin color of farmed individuals is often lacking bright redness. As such, the market price of farmed individuals is lower than wild-type.

RESULTS

To understand the genetic mechanisms of skin coloration in leopard coral grouper, we compared leopard coral grouper with different skin colors through transcriptome analysis. Six cDNA libraries generated from wild-caught leopard coral grouper with different skin colors were characterized by using the Illumina platform. Reference-guided de novo transcriptome data of leopard coral grouper obtained 24,700 transcripts, and 1,089 differentially expressed genes (DEGs) were found between red and brown skin color individuals. The results showed that nine candidate DEGs (epha2, sema6d, acsl4, slc7a5, hipk1, nol6, timp2, slc25a42, and kdf1) significantly associated with skin color were detected by using comparative transcriptome analysis and quantitative real-time polymerase chain reaction (qRT-PCR).

CONCLUSIONS

The findings may provide genetic information for further skin color research, and to boost the market price of farmed leopard coral grouper by selective breeding.

Chromosome-Level Genome Assembly of the Speckled Blue Grouper (Epinephelus cyanopodus) Provides Insight into Its Adaptive Evolution

Epinephelus cyanopodus is a coral reef-dwelling grouper with important economic and ecological value and is widely distributed in the western Pacific Ocean. The lack of genomic resources for E. cyanopodus hinders its adaptive evolution and phylogeny research. We constructed the first high-quality genome of E. cyanopodus based on DNBSEQ, PacBio, and Hic sequencing technologies, with a genome size of 998.82 Mb, contig N50 of 5.855 Mb, and scaffold N50 of 41.98 Mb. More than 99.7% of contigs were anchored to 24 pseudochromosomes, and 94.2% of BUSCO genes were found in the E. cyanopodus genome, indicating a high genome assembly completeness. A total of 26,337 protein-coding genes were predicted, of which 98.77% were functionally annotated. Phylogenetic analysis showed that E. cyanopodus separated from its closely related species Epinephelus akaara about 11.5-26.5 million years ago, and the uplift of the Indo-Australian archipelago may have provided an opportunity for its rapid radiation. Moreover, several gene families associated with innate and adaptive immunity were significantly expanded in speckled blue grouper compared to other teleost genomes. Additionally, we identified several genes associated with immunity, growth and reproduction that are under positive selection in E. cyanopodus compared to other groupers, suggesting that E. cyanopodus has evolved broad adaptability in response to complex survival environment, which may provide the genetic basis for its rapid radiation. In brief, the high-quality reference genome of the speckled blue grouper provides a foundation for research on its biological traits and adaptive evolution and will be an important genetic tool to guide aquaculture and resolve its taxonomic controversies in future studies.

The Insertion in the 3' UTR of Pmel17 Is the Causal Variant for Golden Skin Color in Tilapia

Understanding of the relationships between genotypes and phenotypes is a central problem in biology. Although teleosts have colorful phenotypes, not much is known about their underlying mechanisms. Our previous study showed that golden skin color in Mozambique tilapia was mapped in the major locus containing the Pmel gene, and an insertion in 3' UTR of Pmel17 was fully correlated with the golden color. However, the molecular mechanism of how Pmel17 determines the golden skin color is unknown. In this study, knockout of Pmel17 with CRISPR/Cas9 in blackish tilapias resulted in golden coloration, and rescue of Pmel17 in golden tilapias recovered the wild-type blackish color, indicating that Pmel17 is the gene determining the golden and blackish color. Functional analysis in vitro showed that the insertion in the 3' UTR of Pmel17 reduced the transcripts of Pmel17. Our data supplies more evidence to support that Pmel17 is the gene for blackish and golden colors, and highlights that the insertion in the 3' UTR of Pmel17 is the causative mutation for the golden coloration.

Analysis of opsin gene family of Crimson snapper (Lutjanus erythropterus)

Opsin is a fellow of the G protein-coupled receptors (GPCRs) superfamily. It can be divided into visual and non-visual opsin according to whether it is directly involved in visual imaging. Opsin plays an important role in visual image formation and the regulation of non-image forming functions such as circadian entrainment in the growth, development and evolution of fish. Crimson snapper belongs to Perciforme mainly found in the Indo-West Pacific and the South China Sea. It is one of the most influential economic fishes in the South China Sea. In order to study the existence and expression of opsin gene in Crimson snapper, we sequenced the genome and tissue sample transcriptome of Crimson snapper. In this study, 32 opsin genes were identified from the genome of Crimson snapper. The length of these genes ranged from 1061 bp to 86203 bp and were distributed on 15 different chromosomes. The analysis of opsin gene family of Crimson snapper showed that the sws2 had two extra copies as compared with that of Zebrafish. Domain and motif analysis revealed that all the 32 opsin genes have seven-(pass)-transmembrane domain receptors (7TM receptors) each, and the opsin family contained 10 common motifs. The expression level of opsin gene, confirmed by RT-qPCR, was analyzed by using nine tissues transcriptome databases of Crimson snapper. The results showed that almost all opsin genes were highly expressed in the retina and brain, except opn7a and opn7b which were expressed in intestine and red skin, and almost no expression in other tissues. Our results provide a comprehensive basic knowledge for the opsin gene family of Crimson snapper, which has significance for the study of the function of opsin in Lutjanidaes.

An alternative splicing hypothesis for neuropathology of schizophrenia: evidence from studies on historical candidate genes and multi-omics data

Alternative splicing of schizophrenia risk genes, such as DRD2, GRM3, and DISC1, has been extensively described. Nevertheless, the alternative splicing characteristics of the growing number of schizophrenia risk genes identified through genetic analyses remain relatively opaque. Recently, transcriptomic analyses in human brains based on short-read RNA-sequencing have discovered many "local splicing" events (e.g., exon skipping junctions) associated with genetic risk of schizophrenia, and further molecular characterizations have identified novel spliced isoforms, such as AS3MT^d2d3 and ZNF804A^E3E4. In addition, long-read sequencing analyses of schizophrenia risk genes (e.g., CACNA1C and NRXN1) have revealed multiple previously unannotated brain-abundant isoforms with therapeutic potentials, and functional analyses of KCNH2-3.1 and Ube3a1 have provided examples for investigating such spliced isoforms in vitro and in vivo. These findings suggest that alternative splicing may be an essential molecular mechanism underlying genetic risk of schizophrenia, however, the incomplete annotations of human brain transcriptomes might have limited our understanding of schizophrenia pathogenesis, and further efforts to elucidate these transcriptional characteristics are urgently needed to gain insights into the illness-correlated brain physiology and pathology as well as to translate genetic discoveries into novel therapeutic targets.

Chromosome Genome Assembly of Cromileptes altivelis Reveals Loss of Genome Fragment in Cromileptes Compared with Epinephelus Species

The humpback grouper (Cromileptes altivelis), an Epinephelidae species, is patchily distributed in the reef habitats of Western Pacific water. This grouper possesses a remarkably different body shape and notably low growth rate compared with closely related grouper species. For promoting further research of the grouper, in the present study, a high-quality chromosome-level genome of humpback grouper was assembled using PacBio sequencing and high-throughput chromatin conformation capture (Hi-C) technology. The assembled genome was 1.013 Gb in size with 283 contigs, of which, a total of 143 contigs with 1.011 Gb in size were correctly anchored into 24 chromosomes. Moreover, a total of 26,037 protein-coding genes were predicted, of them, 25,243 (96.95%) genes could be functionally annotated. The high-quality chromosome-level genome assembly will provide pivotal genomic information for future research of the speciation, evolution and molecular-assisted breeding in humpback groupers. In addition, phylogenetic analysis based on shared single-copy orthologues of the grouper species showed that the humpback grouper is included in the Epinephelus genus and clustered with the giant grouper in one clade with a divergence time of 9.86 Myr. In addition, based on the results of collinearity analysis, a gap in chromosome 6 of the humpback grouper was detected; the missed genes were mainly associated with immunity, substance metabolism and the MAPK signal pathway. The loss of the parts of genes involved in these biological processes might affect the disease resistance, stress tolerance and growth traits in humpback groupers. The present research will provide new insight into the evolution and origin of the humpback grouper.

Whole-genome sequencing of brown-marbled grouper (Epinephelus fuscoguttatus) provides insights into adaptive evolution and growth differences

The brown-marbled grouper (Epinephelus fuscoguttatus) is an important species of fish in the coral reef ecosystem and marine aquaculture industry. In this study, a high-quality chromosome-level genome of brown-marbled grouper was assembled using Oxford Nanopore technology and Hi-C technology. The GC content and heterozygosity were approximately 42% and 0.35%, respectively. A total of 230 contigs with a total length of 1047 Mb and contig N50 of 13.8 Mb were assembled, and 228 contigs (99.13%) were anchored into 24 chromosomes. A total of 24,005 protein-coding genes were predicted, among which 23,862 (99.4%) predicted genes were annotated. Phylogenetic analysis showed that brown-marbled grouper and humpback grouper were clustered into one clade that separated approximately 11-23 million years ago. Collinearity analyses showed that there was no obvious duplication of large fragments between chromosomes in the brown-marbled grouper. Genomes of the humpback grouper and giant grouper showed a high collinearity with that of the brown-marbled grouper. A total of 305 expanded gene families were detected in the brown-marbled grouper genome, which is mainly involved in disease resistance. In addition, a genetic linkage map with 3061.88 cM was constructed. Based on the physical and genetic map, one growth-related quantitative trait loci was detected in 32,332,447 bp of chromosome 20, and meox1 and etv4 were considered candidate growth-related genes. This study provides pivotal genetic resources for further evolutionary analyses and artificial breeding of groupers.

Gene signatures of SARS-CoV/SARS-CoV-2-infected ferret lungs in short- and long-term models

Coronaviruses (CoVs) consist of six strains, and the severe acute respiratory syndrome coronavirus (SARS-CoV), newly found coronavirus (SARS-CoV-2) has rapidly spread leading to a global outbreak. The ferret (Mustela putorius furo) serves as a useful animal model for studying SARS-CoV/SARS-CoV-2 infection and developing therapeutic strategies. A holistic approach for distinguishing differences in gene signatures during disease progression is lacking. The present study discovered gene expression profiles of short-term (3 days) and long-term (14 days) ferret models after SARS-CoV/SARS-CoV-2 infection using a bioinformatics approach. Through Gene Ontology (GO) and MetaCore analyses, we found that the development of stemness signaling was related to short-term SARS-CoV/SARS-CoV-2 infection. In contrast, pathways involving extracellular matrix and immune responses were associated with long-term SARS-CoV/SARS-CoV-2 infection. Some highly expressed genes in both short- and long-term models played a crucial role in the progression of SARS-CoV/SARS-CoV-2 infection, including DPP4, BMP2, NFIA, AXIN2, DAAM1, ZNF608, ME1, MGLL, LGR4, ABHD6, and ACADM. Meanwhile, we revealed that metabolic, glucocorticoid, and reactive oxygen species-associated networks were enriched in both short- and long-term infection models. The present study showed alterations in gene expressions from short-term to long-term SARS-CoV/SARS-CoV-2 infection. The current result provides an explanation of the pathophysiology for post-infectious sequelae and potential targets for treatment.

Identification of a functional human-unique 351-bp Alu insertion polymorphism associated with major depressive disorder in the 1p31.1 GWAS risk loci

Genome-wide association studies (GWAS) have reported substantial single-nucleotide polymorphisms (SNPs) associated with major depressive disorder (MDD), but the underlying functional variations in the GWAS risk loci are unclear. Here we show that the European MDD genome-wide risk-associated allele of rs12129573 at 1p31.1 is associated with MDD in Han Chinese, and this SNP is in strong linkage disequilibrium (LD) with a human-unique Alu insertion polymorphism (rs70959274) in the 5' flanking region of a long non-coding RNA (lncRNA) LINC01360 (Long Intergenic Non-Protein Coding RNA 1360), which is preferably expressed in human testis in the currently available expression datasets. The risk allele at rs12129573 is almost completely linked with the absence of this Alu insertion. The Alu insertion polymorphism (rs70959274) is significantly associated with a lower RNA level of LINC01360 and acts as a transcription silencer likely through modulating the methylation of its internal CpG sites. Luciferase assays confirm that the presence of Alu insertion at rs70959274 suppresses transcriptional activities in human cells, and deletion of the Alu insertion through CRISPR/Cas9-directed genome editing increases RNA expression of LINC01360. Deletion of the Alu insertion in human cells also leads to dysregulation of gene expression, biological processes and pathways relevant to MDD, such as the alterations of mRNA levels of DRD2 and FLOT1, transcription of genes involved in synaptic transmission, neurogenesis, learning or memory, and the PI3K-Akt signaling pathway. In summary, we identify a human-unique DNA repetitive polymorphism in robust LD with the MDD risk-associated SNP at the prominent 1p31.1 GWAS loci, and offer insights into the molecular basis of the illness.

Further confirmation of netrin 1 receptor (DCC) as a depression risk gene via integrations of multi-omics data

Genome-wide association studies (GWAS) of major depression and its relevant biological phenotypes have been extensively conducted in large samples, and transcriptome-wide analyses in the tissues of brain regions relevant to pathogenesis of depression, e.g., dorsolateral prefrontal cortex (DLPFC), have also been widely performed recently. Integrating these multi-omics data will enable unveiling of depression risk genes and even underlying pathological mechanisms. Here, we employ summary data-based Mendelian randomization (SMR) and integrative risk gene selector (iRIGS) approaches to integrate multi-omics data from GWAS, DLPFC expression quantitative trait loci (eQTL) analyses and enhancer-promoter physical link studies to prioritize high-confidence risk genes for depression, followed by independent replications across distinct populations. These integrative analyses identify multiple high-confidence depression risk genes, and numerous lines of evidence supporting pivotal roles of the netrin 1 receptor (DCC) gene in this illness across different populations. Our subsequent explorative analyses further suggest that DCC significantly predicts neuroticism, well-being spectrum, cognitive function and putamen structure in general populations. Gene expression correlation and pathway analyses in DLPFC further show that DCC potentially participates in the biological processes and pathways underlying synaptic plasticity, axon guidance, circadian entrainment, as well as learning and long-term potentiation. These results are in agreement with the recent findings of this gene in neurodevelopment and psychiatric disorders, and we thus further confirm that DCC is an important susceptibility gene for depression, and might be a potential target for new antidepressants.