Upregulation of the PPAR signaling pathway and accumulation of lipids are related to the morphological and structural transformation of the dragon-eye goldfish eye

A Comparative Transcriptome and Proteome Analysis of the Molecular Mechanism Underlying Anterior to Dorsal Eye Rotation in the Celestial-Eye Goldfish (Carassius auratus)

Goldfish (Carassius auratus), subjected to millennia of artificial selection and breeding, have diversified into numerous ornamental varieties, such as the celestial-eye (CE) goldfish, noted for its unique dorsal eye rotation. Previous studies have primarily focused on anatomical modifications in CE goldfish eyes, yet the molecular underpinnings of their distinctive eye orientation remain poorly understood. This study employed high-throughput transcriptome and proteome sequencing on 110-day-old full-sibling CE goldfish, which displayed either anterior or upward eye rotations. Verification of these findings was conducted using quantitative PCR (qPCR) for transcriptomic data and parallel reaction monitoring (PRM) for proteomic analysis. Our research identified 73,685 genes and 7717 proteins, pinpointing 8 common differentially expressed genes (DEGs) and proteins (DEPs) implicated in cytoskeleton remodeling, cell adhesion, apoptosis, and optic nerve regeneration. Enrichment analyses further delineated pathways associated with apoptosis, necroptosis, and cell adhesion molecules. The results indicated a significant role for genes involved in cytoskeletal dynamics, nervous system function, and apoptotic processes in the dorsal eye rotation of CE goldfish. Analyses of abnormalities in ocular membrane structures, along with disturbances in lipid and protein synthesis metabolism and energy metabolism during developmental stages, provided compelling evidence for the potential use of CE goldfish as a model organism in studying human eye-related disorders. This investigation provided the first comprehensive transcriptomic and proteomic overview of eye rotation in CE goldfish, offering insights crucial for the genetic breeding of new ornamental fish varieties.

A novel goldfish orthotopic xenograft model of hepatocellular carcinoma developed to evaluate antitumor drug efficacy

Tumor xenograft animal models play a crucial role in hepatocellular carcinoma (HCC) research. Mice xenograft models are time consuming, laborious and expensive while zebrafish tumor xenograft models are cost-effective and effortless. However, the development of orthotopic xenograft models for HCC in zebrafish embryos has been challenging due to the small size of zebrafish livers. In this study, we utilized 7-day-old goldfish embryos as hosts and successfully established an orthotopic xenograft model of HCC in goldfish livers. Through injecting fluorescence labeled HCC cells into the liver of goldfish, we could visualize the proliferation and migration of tumor cells in vivo. In addition, we found that the temperature of 36 °C was better for tumor cell survival in goldfish larvae compared to 28 °C, assessed by EdU and TUNEL assays. Moreover, macrophage infiltration in the goldfish liver could be evaluated by neutral red staining. Finally, we evaluated the efficacy of the targeted therapy drug Sorafenib and the traditional Chinese medicine, Huaier granules, alone or in combination in the goldfish HCC orthotopic xenograft model. We found that the combination therapy showed the best efficacy against HCC cells in terms of macrophage infiltration, polarization as well as tumor cells proliferation, metastasis and apoptosis. In conclusion, the proposed goldfish HCC orthotopic xenograft model opens new avenues for HCC related research, including evaluation of tumor progression, cell interactions in the immune microenvironment, drug efficacy, and screening of anti-tumor drugs.

A Network of Circular RNA and MicroRNA Sequencing Provides Insights into Pigment Deposition of Changshun Blue Eggshell Chickens

Eggshell color plays important biological roles and attracts the attention of both egg retailers and researchers. However, whether non-coding RNAs are involved in pigment deposition among different eggshell colors remains unknown. In this study, RNA sequencing was used to analyse the uterine gland transcriptome (CircRNA and miRNA) of Changshun chicken blue-shell hens producing four different eggshell color eggs including dark blue PK(DB) and light blue (LB), dark brown and greenish (between blue and pink, DP) and pink (p). We found that miR-192-x, targeting SLC16a7, was expressed in DB, DP, and LB groups compared with the PK group, which indicates that miR-192-x may play a role in the blue eggshell color. KEGG and GO analyses showed that the "metabolic pathways" with targeted genes such BLVRA and HMOX1 were detected in dark and light blue color eggshell chickens, which confirms the different ratios of biliverdin and HO-1 involved in the deposition of blue color. As annotated by connectivity analysis, RASGRF1 and RASGRF2, belonging to the RASGRF family, are involved in the Ras signaling pathway, which plays an important role in cell growth, differentiation, metastasis and apoptosis. Our findings enrich the database of circRNA, miRNAs and genes for chicken uterine tissue, which will be useful in accelerating molecular selection for blue eggshell color layers.

Sublethal acetamiprid affects reproduction, development and disrupts gene expression in Binodoxys communis

At present, understanding of neonicotinoid toxicity in arthropods remains limited. We here evaluated the lethal and sublethal effects of acetamiprid in F0 and F1 generations of Binodoxys communis using a range of sublethal concentrations. The 10% lethal concentration (LC10) and half lethal concentration (LC25) of ACE had negative effects on the B. communis survival rate, adult longevity, parasitism rate, and emergence rate, and significantly prolonged the duration of the developmental cycle. ACE also had intergenerational effects, with some biological indices affected in the F1 generation after pesticide exposure. Transcriptomic analysis demonstrated that differentially expressed genes were enriched in specific pathways including the amino acid metabolism, carbohydrate metabolism, energy metabolism, exogenous metabolism, signal transduction, and glutathione metabolism pathways. These results indicated strong contact toxicity of ACE to B. communis, which may inhibit their biological control capacity. These results improve our understanding of the toxicological mechanisms of parasitic natural enemies in response to insecticide exposure.

Multi-omics analysis revealed the differences in lipid metabolism of the gut between adult and juvenile yellowfin tuna (Thunnus albacares)

Introduction

Tuna has a cost-effective energy supply to support the regional endothermic and high-speed swimming performance. The gut symbiotic microbiotas and their metabolites play essential roles in tuna's diet digestion, absorption, and energy acquirement, which are often highly related to the ontogenetic development of tuna.

Methods

We compared gut microbial compositions and metabolites, as well as mRNA expression of the intestine between juvenile and adult yellowfin tuna using 16S rRNA sequencing, metabolomic and transcriptomic, respectively.

Results and discussion

The results revealed that adults had a significantly higher microbial diversity and abundance of Acinetobacter than juveniles. Regarding the gut microbiota-derived metabolites, fatty acids, especially glycerophospholipid and sphingolipid, were significantly enriched in adults than in juveniles. Moreover, the short-chain fatty acid (butyrate and isobutyrate) contents were significantly higher in adults than in juveniles. To find the relationship between gut microbiotas and host physiology, intestinal transcriptome analysis demonstrated that the enriched pathways of differential expression genes (DEGs) in adult tuna were the lipid metabolism pathway, including "fat digestion and absorption," "cholesterol metabolism," "steroid hormone biosynthesis," "glycerolipid metabolism," and "glycerophospholipid metabolism." However, protein digestion and absorption and pancreatic secretion pathways were significantly enriched in the juveniles. The conjoint analysis indicated that the enriched pathways of both differential metabolites (DMs) and DEGs were remarkably related to the regulation of glycerophospholipids metabolism in adult tunas. This study highlights the role of gut microbiotas in fish nutrition metabolism. These findings provide new insights into the view of ontogenetic shifts of gut microbiotas and their metabolites on host health and gut function in endothermic and high-speed swimming marine fish species.

Translocation of gut microbes to epididymal white adipose tissue drives lipid metabolism disorder under heat stress

Heat stress induces multi-organ damage and serious physiological dysfunction in mammals, and gut bacteria may translocate to extra-intestinal tissues under heat stress pathology. However, whether gut bacteria translocate to the key metabolic organs and impair function as a result of heat stress remains unknown. Using a heat stress-induced mouse model, heat stress inhibited epididymal white adipose tissue (eWAT) expansion and induced lipid metabolic disorder but did not damage other organs, such as the heart, liver, spleen, or muscle. Microbial profiling analysis revealed that heat stress shifted the bacterial community in the cecum and eWAT but not in the inguinal white adipose tissue, blood, heart, liver, spleen, or muscle. Notably, gut-vascular barrier function was impaired, and the levels of some bacteria, particularly Lactobacillus, were higher in the eWAT, as confirmed by catalyzed reporter deposition fluorescence in situ hybridization (CARD-FISH) staining when mice were under heat stress. Moreover, integrated multi-omics analysis showed that the eWAT microbiota was associated with host lipid metabolism, and the expression of genes involved in the lipid metabolism in eWAT was upregulated under heat stress. A follow-up microbial supplementation study after introducing Lactobacillus plantarum to heat-stressed mice revealed that the probiotic ameliorated heat stress-induced loss of eWAT and dyslipidemia and reduced gut bacterial translocation to the eWAT by improving gut barrier function. Overall, our findings suggest that gut bacteria, particularly Lactobacillus spp., play a crucial role in heat stress-induced lipid metabolism disorder and that there is therapeutic potential for using probiotics, such as Lactobacillus plantarum.

ACSM5 inhibits ligamentum flavum hypertrophy by regulating lipid accumulation mediated by FABP4/PPAR signaling pathway

BACKGROUND

Ligamentum flavum (LF) hypertrophy is the main cause of lumbar spinal canal stenosis (LSCS). Previous studies have shown that LF hypertrophy tissue exhibits abnormal lipid accumulation, but the regulatory mechanism remains unclear. The objective of this study was to explore the function and potential mechanism of ACSM5 in LF lipid accumulation.

METHODS

To assess the ACSM5 expression levels, lipid accumulation and triglyceride (TG) level in LF hypertrophy and normal tissue, we utilized RT-qPCR, western blot, oil red O staining, and TG assay kit. The pearson correlation coefficient assay was used to analyze the correlation between ACSM5 levels and lipid accumulation or TG levels in LF hypertrophy tissue. The role of ACSM5 in free fatty acids (FFA)-induced lipid accumulation in LF cells was assessed in vitro, and the role of ACSM5 in LF hypertrophy in mice was verified in vivo. To investigate the underlying mechanisms of ACSM5 regulating lipid accumulation in LF, we conducted the mRNA sequencing, bioinformatics analysis, and rescue experiments.

RESULTS

In this study, we found that ACSM5, which was significantly down-regulated in LF tissues, correlated with lipid accumulation. In vitro cell experiments demonstrated that overexpression of ACSM5 significantly inhibited FFA-induced lipid accumulation and fibrosis in LF cells. In vivo animal experiments further confirmed that overexpression of ACSM5 inhibited LF thickening, lipid accumulation, and fibrosis. Mechanistically, ACSM5 inhibited lipid accumulation of LF cells by inhibiting FABP4-mediated PPARγ signaling pathway, thereby improving hypertrophy and fibrosis of LF.

CONCLUSIONS

our findings elucidated the important role of ACSM5 in the regulation of LF lipid accumulation and provide insight into potential therapeutic interventions for the treatment of LF hypertrophy. This study further suggested that therapeutic strategies targeting lipid deposition may be an effective potential approach to treat LF hypertrophy-induced LSCS.

Long-Term Tetrabromobisphenol A Exposure Induces Gut Microbiota Imbalance and Metabolic Disorders via the Peroxisome Proliferator-Activated Receptor Signaling Pathway in the Regenerated Gut of Apostichopus japonicus

Tetrabromobisphenol A (TBBPA), a commonly utilized brominated flame retardant, is found in many types of abiotic and biotic matrices. TBBPA can increase oxidative stress, disrupt the endocrine system, cause neurodevelopmental disorders and activate peroxisome proliferator-activated receptors to modulate lipid deposits in aquatic animals. However, the toxic mechanism of TBBPA on the gut microbiota and intestinal health remains unclear. Apostichopus japonicus is an ideal model for studying the relationship between environmental contaminants and intestinal health due to its unique capacity for evisceration and quickly regenerated intestine. In the present study, we investigated the toxic mechanism of TBBPA on the gut microbiota and intestinal health in the regenerated intestine of A. japonicus. The results show that TBBPA exposure decreased the health of the regenerated intestine and the enzymatic activities, alpha diversity indices, and the relative abundance of the gut microbiota. Transcriptome analysis shows that TBBPA exposure affected lipid metabolism via the PPAR signaling pathway during the process of intestinal regeneration in A. japonicus, suggesting that TBBPA exposure can affect the composition and function of the gut microbiota and intestinal health in the regenerated intestine of A. japonicus. These results provide a basis for further research on the potential toxicity of TBBPA to the intestinal health in animals.

PANoptosis-based molecular subtyping and HPAN-index predicts therapeutic response and survival in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is a highly prevalent and fatal cancer. The role of PANoptosis, a novel form of programmed cell death, in HCC is yet to be fully understood. This study focuses on identifying and analyzing PANoptosis-associated differentially expressed genes in HCC (HPAN_DEGs), aiming to enhance our understanding of HCC pathogenesis and potential treatment strategies.

Methods

We analyzed HCC differentially expressed genes from TCGA and IGCG databases and mapped them to the PANoptosis gene set, identifying 69 HPAN_DEGs. These genes underwent enrichment analyses, and consensus clustering analysis was used to determine three distinct HCC subgroups based on their expression profiles. The immune characteristics and mutation landscape of these subgroups were evaluated, and drug sensitivity was predicted using the HPAN-index and relevant databases.

Results

The HPAN_DEGs were mainly enriched in pathways associated with the cell cycle, DNA damage, Drug metabolism, Cytokines, and Immune receptors. We identified three HCC subtypes (Cluster_1, SFN+PDK4-; Cluster_2, SFN-PDK4+; Cluster_3, SFN/PDK4 intermediate expression) based on the expression profiles of the 69 HPAN_DEGs. These subtypes exhibited distinct clinical outcomes, immune characteristics, and mutation landscapes. The HPAN-index, generated by machine learning using the expression levels of 69 HPAN_DEGs, was identified as an independent prognostic factor for HCC. Moreover, the high HPAN-index group exhibited a high response to immunotherapy, while the low HPAN-index group showed sensitivity to small molecule targeted drugs. Notably, we observed that the YWHAB gene plays a significant role in Sorafenib resistance.

Conclusion

This study identified 69 HPAN_DEGs crucial to tumor growth, immune infiltration, and drug resistance in HCC. Additionally, we discovered three distinct HCC subtypes and constructed an HPAN-index to predict immunotherapeutic response and drug sensitivity. Our findings underscore the role of YWHAB in Sorafenib resistance, presenting valuable insights for personalized therapeutic strategy development in HCC.

Genomic and Transcriptional Profiles of Kelch-like (klhl) Gene Family in Polyploid Carassius Complex

Genome duplication supplies raw genetic materials and has been thought to be essential for evolutionary innovation and ecological adaptation. Here, we select Kelch-like (klhl) genes to study the evolution of the duplicated genes in the polyploid Carassius complex, including amphidiploid C. auratus and amphitriploid C. gibelio. Phylogenetic, chromosomal location and read coverage analyses indicate that most of Carassius klhl genes exhibit a 2:1 relationship with zebrafish orthologs and confirm two rounds of polyploidy, an allotetraploidy followed by an autotriploidy, occurred during Carassius evolution. The lineage-specific expansion and biased retention/loss of klhl genes are also found in Carassius. Transcriptome analyses across eight adult tissues and seven embryogenesis stages reveal varied expression dominance and divergence between the two species. The expression of klhls in response to Carassius herpesvirus 2 infection shows different expression changes corresponding to distinct herpesvirus resistances in three C. gibelio gynogenetic clones. Finally, we find that most C. gibelio klhl genes possess three alleles except eight genes that have lost one or two alleles due to genome rearrangement. The allele expression bias is prosperous for Cgklhl genes and varies during embryogenesis owning to the sequential expression manner of the alleles. The current study provides global insights into the genomic and transcriptional evolution of duplicated genes in a given superfamily resulting from multiple rounds of polyploidization.

Gonadal transcriptomes reveal sex-biased expression genes associated with sex determination and differentiation in red-tail catfish (Hemibagrus wyckioides)

BACKGROUND

Red-tail catfish (Hemibagrus wyckioides) is an important commercially farmed catfish in southern China. Males of red-tail catfish grow faster than females, suggesting that all-male catfish will produce more significant economic benefits in aquaculture practice. However, little research has been reported on sex determination and gonadal development in red-tail catfish.

RESULTS

In this study, we performed the first transcriptomic analysis of male and female gonads at four developmental stages at 10, 18, 30, and 48 days post hatching (dph) using RNA-seq technology. A total of 23,588 genes were screened in 24 sequenced samples, of which 28, 213, 636, and 1381 differentially expressed genes (DEGs) were detected at four developmental stages, respectively. Seven candidate genes of sex determination and differentiation were further identified. Real-time quantitative PCR (RT-qPCR) further confirmed that anti-Mullerian hormone (amh), growth differentiation factor 6a (gdf6a), testis-specific gene antigen 10 (tsga10), and cytochrome P450 family 17 subfamily A (cyp17a) were highly expressed mainly in the male, while cytochrome P450 family 19 subfamily A polypeptide 1b (cyp19a1b), forkhead box L2 (foxl2), and hydroxysteroid 17-beta dehydrogenase 1 (hsd17b1) were highly expressed in the female. The KEGG pathway enrichment data showed that these identified DEGs were mainly involved in steroid hormone biosynthesis and TGF-β signaling pathways.

CONCLUSIONS

Based on RNA-seq data of gonads at the early developmental stages, seven DEGs shared by the four developmental stages were identified, among which amh and gdf6a may be the male-biased expression genes, while foxl2, cyp19a1b and hsd17b1 may be the female-biased expression genes in red-tail catfish. Our study will provide crucial genetic information for the research on sex control in red-tail catfish, as well as for exploring the evolutionary processes of sex determination mechanisms in fish.

Barbel regeneration and function divergence in red-tail catfish (Hemibagrus wyckioides) based on the chromosome-level genomes and comparative transcriptomes

Catfish (Siluriformes) are one of the most diverse vertebrate orders and are characterized by whisker-like barbels, which are important sensory organs in most of teleosts. However, their specific biological functions are still unclear. Red-tail catfish (Hemibagrus wyckioides) is well-known catfish species with four pairs of barbels, of which the maxillary barbels reach two-thirds of the body length. In this study, a 776.58 Mb high-quality chromosome-level genome was assembled into 29 chromosomes. Comparative genome data indicated that the barbeled regeneration gene ccl33 has expanded into 11 tandemly duplicated copies. Transcriptome data revealed the functional differentiation of different barbels and suggested that the maxillary barbel might be necessary for water temperature perception. Taste receptor genes were also characterized in teleosts with different food habits. Selection pressures were revealed to affect the sugar-based solute transport domain of the sweet taste receptor gene t1r2 in carnivorous fishes. In addition, the bitter taste receptor gene t2r200 was found to be lost from the genomes of four catfish species. Therefore, our study provides a genomic foundation for understanding the regeneration and functional differentiation of barbels in red-tail catfish and also reveals novel insights into the feeding evolution of fish species with different feeding habits.

Reproductive Toxicity and Teratogenicity of Fluorene-9-bisphenol on Chinese Medaka (Oryzias sinensis): A Study from Laboratory to Field

Fluorene-9-bisphenol (BHPF), a bisphenol A (BPA) substitute, has been increasingly used as a material in syntheses of polymers that are widely used in road markings, artificial tracks, coating floors, building paints, etc., increasing the likelihood of BHPF contamination in the aquatic environment due to its release from the products. However, to date, it is unknown whether it may have actual impacts on fish in real environments. In this study, a 105-day exposure experiment of BHPF at various concentrations (0.01, 0.1, 1, and 10 μg/L) on Chinese medaka (Oryzias sinensis) was performed under laboratory conditions and found decreased fecundity, such as lower egg qualities and quantities, retarded oogenesis, and atretic follicles in the fish and deformed eyes and bodies in its F1 generation. Toxico-transcriptome analyses showed that estrogen-responsive genes were significantly suppressed by BHPF, indicating that antagonist properties of BHPF on estrogen receptors might be causes for the decreased fecundity. Field investigations (Beijing) demonstrated that BHPF was detectable in 60% surface waters, with a mean concentration of 10.49 ± 6.33 ng/L, by gas chromatography-mass spectrometry, and similar effects in wild Chinese medaka were also observed, some of which the parameters were found to be obviously correlated with the BHPF levels in corresponding waters.

Comprehensive Analysis of the lncRNA-miRNA-mRNA Regulatory Network for Intramuscular Fat in Pigs

Intramuscular fat (IMF) is an essential trait closely related to meat quality. The IMF trait is a complex quantitative trait that is regulated by multiple genes. In order to better understand the process of IMF and explore the key factors affecting IMF deposition, we identified differentially expressed mRNA, miRNA, and lncRNA in the longissimus dorsi muscle (LD) between Songliao Black (SL) pigs and Landrace pigs. We obtained 606 differentially expressed genes (DEGs), 55 differentially expressed miRNAs (DEMs), and 30 differentially expressed lncRNAs (DELs) between the SL pig and Landrace pig. Enrichment results from GO and KEGG indicate that DEGs are involved in fatty acid metabolism and some pathways related to glycogen synthesis. We constructed an lncRNA-miRNA-mRNA interaction network with 18 DELs, 11 DEMs, and 42 DEGs. Finally, the research suggests that ARID5B, CPT1B, ACSL1, LPIN1, HSP90AA1, IRS1, IRS2, PIK3CA, PIK3CB, and PLIN2 may be the key genes affecting IMF deposition. The LncRNAs MSTRG.19948.1, MSTRG.13120.1, MSTRG.20210.1, and MSTRG.10023.1, and the miRNAs ssc-miRNA-429 and ssc-miRNA-7-1, may play a regulatory role in IMF deposition through their respective target genes. Our research provides a reference for further understanding the regulatory mechanism of IMF.

Comparative genome anatomy reveals evolutionary insights into a unique amphitriploid fish

Triploids are rare in nature because of difficulties in meiotic and gametogenic processes, especially in vertebrates. The Carassius complex of cyprinid teleosts contains sexual tetraploid crucian carp/goldfish (C. auratus) and unisexual hexaploid gibel carp/Prussian carp (C. gibelio) lineages, providing a valuable model for studying the evolution and maintenance mechanism of unisexual polyploids in vertebrates. Here we sequence the genomes of the two species and assemble their haplotypes, which contain two subgenomes (A and B), to the chromosome level. Sequencing coverage analysis reveals that C. gibelio is an amphitriploid (AAABBB) with two triploid sets of chromosomes; each set is derived from a different ancestor. Resequencing data from different strains of C. gibelio show that unisexual reproduction has been maintained for over 0.82 million years. Comparative genomics show intensive expansion and alterations of meiotic cell cycle-related genes and an oocyte-specific histone variant. Cytological assays indicate that C. gibelio produces unreduced oocytes by an alternative ameiotic pathway; however, sporadic homologous recombination and a high rate of gene conversion also exist in C. gibelio. These genomic changes might have facilitated purging deleterious mutations and maintaining genome stability in this unisexual amphitriploid fish. Overall, the current results provide novel insights into the evolutionary mechanisms of the reproductive success in unisexual polyploid vertebrates.

Genome sequencing and haplotype assembly of two cyprinid teleosts, a sexual tetraploid and an unisexual hexaploid, reveal insights into the evolutionary mechanisms underpinning the reproductive success of unisexual polyploid vertebrates.

Two duplicated gsdf homeologs cooperatively regulate male differentiation by inhibiting cyp19a1a transcription in a hexaploid fish

Although evolutionary fates and expression patterns of duplicated genes have been extensively investigated, how duplicated genes co-regulate a biological process in polyploids remains largely unknown. Here, we identified two gsdf (gonadal somatic cell-derived factor) homeologous genes (gsdf-A and gsdf-B) in hexaploid gibel carp (Carassius gibelio), wherein each homeolog contained three highly conserved alleles. Interestingly, gsdf-A and gsdf-B transcription were mainly activated by dmrt1-A (dsx- and mab-3-related transcription factor 1) and dmrt1-B, respectively. Loss of either gsdf-A or gsdf-B alone resulted in partial male-to-female sex reversal and loss of both caused complete sex reversal, which could be rescued by a nonsteroidal aromatase inhibitor. Compensatory expression of gsdf-A and gsdf-B was observed in gsdf-B and gsdf-A mutants, respectively. Subsequently, we determined that in tissue culture cells, Gsdf-A and Gsdf-B both interacted with Ncoa5 (nuclear receptor coactivator 5) and blocked Ncoa5 interaction with Rora (retinoic acid-related orphan receptor-alpha) to repress Rora/Ncoa5-induced activation of cyp19a1a (cytochrome P450, family 19, subfamily A, polypeptide 1a). These findings illustrate that Gsdf-A and Gsdf-B can regulate male differentiation by inhibiting cyp19a1a transcription in hexaploid gibel carp and also reveal that Gsdf-A and Gsdf-B can interact with Ncoa5 to suppress cyp19a1a transcription in vitro. This study provides a typical case of cooperative mechanism of duplicated genes in polyploids and also sheds light on the conserved evolution of sex differentiation.

Two Duplicated Ptpn6 Homeologs Cooperatively and Negatively Regulate RLR-Mediated IFN Response in Hexaploid Gibel Carp

Src homology region 2 domain-containing phosphatase 1 (SHP1), encoded by the protein tyrosine phosphatase nonreceptor type 6 (ptpn6) gene, belongs to the family of protein tyrosine phosphatases (PTPs) and participates in multiple signaling pathways of immune cells. However, the mechanism of SHP1 in regulating fish immunity is largely unknown. In this study, we first identified two gibel carp (Carassius gibelio) ptpn6 homeologs (Cgptpn6-A and Cgptpn6-B), each of which had three alleles with high identities. Then, relative to Cgptpn6-B, dominant expression in adult tissues and higher upregulated expression of Cgptpn6-A induced by polyinosinic-polycytidylic acid (poly I:C), poly deoxyadenylic-deoxythymidylic (dA:dT) acid and spring viremia of carp virus (SVCV) were uncovered. Finally, we demonstrated that CgSHP1-A (encoded by the Cgptpn6-A gene) and CgSHP1-B (encoded by the Cgptpn6-B gene) act as negative regulators of the RIG-I-like receptor (RLR)-mediated interferon (IFN) response via two mechanisms: the inhibition of CaTBK1-induced phosphorylation of CaMITA shared by CgSHP1-A and CgSHP1-B, and the autophagic degradation of CaMITA exclusively by CgSHP1-A. Meanwhile, the data support that CgSHP1-A and CgSHP1-B have sub-functionalized and that CgSHP1-A overwhelmingly dominates CgSHP1-B in the process of RLR-mediated IFN response. The current study not only sheds light on the regulative mechanism of SHP1 in fish immunity, but also provides a typical case of duplicated gene evolutionary fates.

Global Analysis of Transcriptome and Translatome Revealed That Coordinated WNT and FGF Regulate the Carapacial Ridge Development of Chinese Soft-Shell Turtle

The turtle carapace is composed of severely deformed fused dorsal vertebrae, ribs, and bone plates. In particular, the lateral growth in the superficial layer of turtle ribs in the dorsal trunk causes an encapsulation of the scapula and pelvis. The recent study suggested that the carapacial ridge (CR) is a new model of epithelial–mesenchymal transition which is essential for the arrangement of the ribs. Therefore, it is necessary to explore the regulatory mechanism of carapacial ridge development to analyze the formation of the turtle shell. However, the current understanding of the regulatory network underlying turtle carapacial ridge development is poor due to the lack of both systematic gene screening at different carapacial ridge development stages and gene function verification studies. In this study, we obtained genome-wide gene transcription and gene translation profiles using RNA sequencing and ribosome nascent-chain complex mRNA sequencing from carapacial ridge tissues of Chinese soft-shell turtle at different development stages. A correlation analysis of the transcriptome and translatome revealed that there were 129, 670, and 135 codifferentially expressed genes, including homodirection and opposite-direction differentially expressed genes, among three comparison groups, respectively. The pathway enrichment analysis of codifferentially expressed genes from the Kyoto Encyclopedia of Genes and Genomes showed dynamic changes in signaling pathways involved in carapacial ridge development. Especially, the results revealed that the Wnt signaling pathway and MAPK signaling pathway may play important roles in turtle carapacial ridge development. In addition, Wnt and Fgf were expressed during the carapacial ridge development. Furthermore, we discovered that Wnt5a regulated carapacial ridge development through the Wnt5a/JNK pathway. Therefore, our studies uncover that the morphogenesis of the turtle carapace might function through the co-operation between conserved WNT and FGF signaling pathways. Consequently, our findings revealed the dynamic signaling pathways acting on the carapacial ridge development of Chinese soft-shell turtle and provided new insights into uncover the molecular mechanism underlying turtle shell morphogenesis.

Divergent Antiviral Mechanisms of Two Viperin Homeologs in a Recurrent Polyploid Fish

Polyploidy and subsequent diploidization provide genomic opportunities for evolutionary innovations and adaptation. The researches on duplicated gene evolutionary fates in recurrent polyploids have seriously lagged behind that in paleopolyploids with diploidized genomes. Moreover, the antiviral mechanisms of Viperin remain largely unclear in fish. Here, we elaborate the distinct antiviral mechanisms of two viperin homeologs (Cgviperin-A and Cgviperin-B) in auto-allo-hexaploid gibel carp (Carassius gibelio). First, Cgviperin-A and Cgviperin-B showed differential and biased expression patterns in gibel carp adult tissues. Subsequently, using co-immunoprecipitation (Co-IP) screening analysis, both CgViperin-A and CgViperin-B were found to interact with crucian carp (C. auratus) herpesvirus (CaHV) open reading frame 46 right (ORF46R) protein, a negative herpesvirus regulator of host interferon (IFN) production, and to promote the proteasomal degradation of ORF46R via decreasing K63-linked ubiquitination. Additionally, CgViperin-B also mediated ORF46R degradation through autophagosome pathway, which was absent in CgViperin-A. Moreover, we found that the N-terminal α-helix domain was necessary for the localization of CgViperin-A and CgViperin-B at the endoplasmic reticulum (ER), and the C-terminal domain of CgViperin-A and CgViperin-B was indispensable for the interaction with degradation of ORF46R. Therefore, the current findings clarify the divergent antiviral mechanisms of the duplicated viperin homeologs in a recurrent polyploid fish, which will shed light on the evolution of teleost duplicated genes.

Genomic anatomy of male-specific microchromosomes in a gynogenetic fish

Unisexual taxa are commonly considered short-lived as the absence of meiotic recombination is supposed to accumulate deleterious mutations and hinder the creation of genetic diversity. However, the gynogenetic gibel carp (Carassius gibelio) with high genetic diversity and wide ecological distribution has outlived its predicted extinction time of a strict unisexual reproduction population. Unlike other unisexual vertebrates, males associated with supernumerary microchromosomes have been observed in gibel carp, which provides a unique system to explore the rationales underlying male occurrence in unisexual lineage and evolution of unisexual reproduction. Here, we identified a massively expanded satellite DNA cluster on microchromosomes of hexaploid gibel carp via comparing with the ancestral tetraploid crucian carp (Carassius auratus). Based on the satellite cluster, we developed a method for single chromosomal fluorescence microdissection and isolated three male-specific microchromosomes in a male metaphase cell. Genomic anatomy revealed that these male-specific microchromosomes contained homologous sequences of autosomes and abundant repetitive elements. Significantly, several potential male-specific genes with transcriptional activity were identified, among which four and five genes displayed male-specific and male-biased expression in gonads, respectively, during the developmental period of sex determination. Therefore, the male-specific microchromosomes resembling common features of sex chromosomes may be the main driving force for male occurrence in gynogenetic gibel carp, which sheds new light on the evolution of unisexual reproduction.

Genotypic Males Play an Important Role in the Creation of Genetic Diversity in Gynogenetic Gibel Carp

Unisexual lineages are commonly considered to be short-lived in the evolutionary process as accumulation of deleterious mutations stated by Muller’s ratchet. However, the gynogenetic hexaploid gibel carp (Carassius gibelio) with existence over 0.5 million years has wider ecological distribution and higher genetic diversity than its sexual progenitors, which provides an ideal model to investigate the underlying mechanisms on countering Muller’s ratchet in unisexual taxa. Unlike other unisexual lineages, the wild populations of gibel carp contain rare and variable proportions of males (1–26%), which are determined via two strategies including genotypic sex determination and temperature-dependent sex determination. Here, we used a maternal gibel carp from strain F to be mated with a genotypic male from strain A⁺, a temperature-dependent male from strain A⁺, and a male from another species common carp (Cyprinus carpio), respectively. When the maternal individual was mated with the genotypic male, a variant of gynogenesis was initiated, along with male occurrence, accumulation of microchromosomes, and creation of genetic diversity in the offspring. When the maternal individual was mated with the temperature-dependent male and common carp, typical gynogenesis was initiated that all the offspring showed the same genetic information as the maternal individual. Subsequently, we found out that the genotypic male nucleus swelled and contacted with the female nucleus after fertilization although it was extruded from the female nucleus eventually, which might be associated with the genetic variation in the offspring. These results reveal that genotypic males play an important role in the creation of genetic diversity in gynogenetic gibel carp, which provides insights into the evolution of unisexual reproduction.

Comparative transcriptomic analysis reveals an association of gibel carp fatty liver with ferroptosis pathway

BACKGROUND

Fatty liver has become a main problem that causes huge economic losses in many aquaculture modes. It is a common physiological or pathological phenomenon in aquaculture, but the causes and occurring mechanism are remaining enigmatic.

METHODS

Each three liver samples from the control group of allogynogenetic gibel carp with normal liver and the overfeeding group with fatty liver were collected randomly for the detailed comparison of histological structure, lipid accumulation, transcriptomic profile, latent pathway identification analysis (LPIA), marker gene expression, and hepatocyte mitochondria analyses.

RESULTS

Compared to normal liver, larger hepatocytes and more lipid accumulation were observed in fatty liver. Transcriptomic analysis between fatty liver and normal liver showed a totally different transcriptional trajectory. GO terms and KEGG pathways analyses revealed several enriched pathways in fatty liver, such as lipid biosynthesis, degradation accumulation, peroxidation, or metabolism and redox balance activities. LPIA identified an activated ferroptosis pathway in the fatty liver. qPCR analysis confirmed that gpx4, a negative regulator of ferroptosis, was significantly downregulated while the other three positively regulated marker genes, such as acsl4, tfr1 and gcl, were upregulated in fatty liver. Moreover, the hepatocytes of fatty liver had more condensed mitochondria and some of their outer membranes were almost ruptured.

CONCLUSIONS

We reveal an association between ferroptosis and fish fatty liver for the first time, suggesting that ferroptosis might be activated in liver fatty. Therefore, the current study provides a clue for future studies on fish fatty liver problems.

Functional Divergence of Multiple Duplicated Foxl2 Homeologs and Alleles in a Recurrent Polyploid Fish

Evolutionary fates of duplicated genes have been widely investigated in many polyploid plants and animals, but research is scarce in recurrent polyploids. In this study, we focused on foxl2, a central player in ovary, and elaborated the functional divergence in gibel carp (Carassius gibelio), a recurrent auto-allo-hexaploid fish. First, we identified three divergent foxl2 homeologs (Cgfoxl2a-B, Cgfoxl2b-A, and Cgfoxl2b-B), each of them possessing three highly conserved alleles and revealed their biased retention/loss. Then, their abundant sexual dimorphism and biased expression were uncovered in hypothalamic-pituitary-gonadal axis. Significantly, granulosa cells and three subpopulations of thecal cells were distinguished by cellular localization of CgFoxl2a and CgFoxl2b, and the functional roles and the involved process were traced in folliculogenesis. Finally, we successfully edited multiple foxl2 homeologs and/or alleles by using CRISPR/Cas9. Cgfoxl2a-B deficiency led to ovary development arrest or complete sex reversal, whereas complete disruption of Cgfoxl2b-A and Cgfoxl2b-B resulted in the depletion of germ cells. Taken together, the detailed cellular localization and functional differences indicate that Cgfoxl2a and Cgfoxl2b have subfunctionalized and cooperated to regulate folliculogenesis and gonad differentiation, and Cgfoxl2b has evolved a new function in oogenesis. Therefore, the current study provides a typical case of homeolog/allele diversification, retention/loss, biased expression, and sub-/neofunctionalization in the evolution of duplicated genes driven by polyploidy and subsequent diploidization from the recurrent polyploid fish.