Genomic, Transcriptomic, and Epigenomic Features Differentiate Genes That Are Relevant for Muscular Polyunsaturated Fatty Acids in the Common Carp

Contribution of elovl5a to Docosahexaenoic Acid (DHA) Synthesis at the Transcriptional Regulation Level in Common Carp, Cyprinus carpio

Docosahexaenoic acid (DHA) is an essential nutrient for humans and plays a critical role in human development and health. Freshwater fish, such as the common carp (Cyprinus carpio), have a certain degree of DHA biosynthesis ability and could be a supplemental source of human DHA needs. The elongase of very-long-chain fatty acid 5 (Elovl5) is an important enzyme affecting polyunsaturated fatty acid (PUFA) biosynthesis. However, the function and regulatory mechanism of the elovl5 gene related to DHA synthesis in freshwater fish is not clear yet. Previous studies have found that there are two copies of the elovl5 gene, elovl5a and elovl5b, which have different functions. Our research group found significant DHA content differences among individuals in Yellow River carp (Cyprinus carpio var.), and four candidate genes were found to be related to DHA synthesis through screening. In this study, the expression level of elovl5a is decreased in the high-DHA group compared to the low-DHA group, which indicated the down-regulation of elovl5a in the DHA synthesis pathways of Yellow River carp. In addition, using a dual-luciferase reporter gene assay, we found that by targeting the 3'UTR region of elovl5a, miR-26a-5p could regulate DHA synthesis in common carp. After CRISPR/Cas9 disruption of elovl5a, the DHA content in the disrupted group was significantly higher than in the wildtype group; meanwhile, the expression level of elovl5a in the disrupted group was significantly reduced compared with the wildtype group. These results suggest that elovl5a may be down-regulating DHA synthesis in Yellow River carp. This study could provide useful information for future research on the genes and pathways that affect DHA synthesis.

Regulation of omega-3 fatty acids production by different genes in freshwater fish species: a review

The present study aims to compare the gene expression of three different fish species (common carp, tilapia, and trout) with varying levels of fatty acids (FA). Based on transcriptome analysis and RNA sequencing, various genes and their associated metabolic pathways are identified. Pathways are categorized based on the genes they encode. Genes that were differentially expressed and their promoter's methylation patterns were revealed by RNA-seq analysis in common carp. Furthermore, fatty acid-enriched pathways, such as ARA4 and adipocytokine signaling, were also identified. Many genes and pathways may influence tilapia's growth and omega-3 content. Using the mTOR pathway, trout with differential expression were discovered to be involved in producing omega-3 fatty acids. This study revealed major pathways in fish species to produce omega-3 fatty acids.

Immunity of turbot Induced by inactivated vaccine of Aeromonas salmonicida from the perspective of DNA methylation

Introduction

DNA methylation was one of the most important modification in epigenetics and played an important role in immune response. Since the introduction of Scophthalmus maximus, the scale of breeding has continued to expand, during which diseases caused by various bacteria, viruses and parasites have become increasingly serious. Therefore, the inactivated vaccines have been widely researched and used in the field of aquatic products with its unique advantages. However, the immune mechanism that occurred in turbot after immunization with inactivated vaccine of Aeromonas salmonicida was not clear.

Methods

In this study, differentially methylated regions (DMRs) were screened by Whole Genome Bisulfite Sequencing (WGBS) and significantly differentially expressed genes (DEGs) were screened by Transcriptome sequencing. Double luciferase report assay and DNA pull-down assay were further verified the DNA methylation state of the gene promoter region affected genes transcriptional activity after immunization with inactivated vaccine of Aeromonas salmonicida.

Results

A total of 8149 differentially methylated regions (DMRs) were screened, in which there were many immune-related genes with altered DNA methylation status. Meanwhile, 386 significantly differentially expressed genes (DEGs) were identified, many of which were significantly enriched in Toll-like receptor signaling pathway, NOD-like receptor signaling pathway and C-type lectin receptor signaling pathway. Combined analysis of WGBS results and RNA-seq results, a total of 9 DMRs of negatively regulated genes are located in the promoter region, including 2 hypermethylated genes with lower expression and 7 hypomethylated genes with higher expression. Then, two immune-related genes C5a anaphylatoxin chemotactic receptor 1-like (C5ar1-Like) and Eosinophil peroxidase-like (EPX-Like), were screened to explore the regulation mechanism of DNA methylation modification on their expression level. Moreover, the DNA methylation state of the gene promoter region affected genes transcriptional activity by inhibiting the binding of transcription factors, which lead to changes in the expression level of the gene.

Discussion

We jointly analyzed WGBS and RNA-seq results and revealed the immune mechanism that occurred in turbot after immunized with inactivated vaccine of A. salmonicida from the perspective of DNA methylation.

Characterization of fatty acid compositions in longissimus thoracis muscle and identification of candidate gene and SNPs related to polyunsaturated fatty acid in Hu sheep

Fatty acid (FA) composition contributes greatly to the quality and nutritional value of lamb meat. In the present study, FA was measured in longissimus thoracis (LT) muscles of 1,085 Hu sheep using gas chromatography. Comparative transcriptomic analysis was conducted in LT muscles to identify differentially expressed genes (DEGs) between six individuals with high polyunsaturated fatty acids (H-PUFA, 15.27% ± 0.42%) and six with low PUFA (L-PUFA, 5.22% ± 0.25%). Subsequently, the single nucleotide polymorphisms (SNPs) in a candidate gene PLIN2 were correlated with FA traits. The results showed a total of 29 FA compositions and 8 FA groups were identified, with the highest content of monounsaturated fatty acids (MUFA, 46.54%, mainly C18:1n9c), followed by saturated fatty acids (SFA, 44.32%, mainly C16:0), and PUFA (8.72%, mainly C18:2n6c), and significant correlations were observed among the most of FA traits. Transcriptomic analyses identified 110 upregulated and 302 downregulated DEGs between H-PUFA and L-PUFA groups. The functional enrichment analysis revealed three significant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and 17 gene ontology (GO) terms, in which regulation of lipolysis in adipocytes, the AMPK signaling pathway, and the PPAR signaling pathway may play important roles in FA metabolism and biosynthesis. In addition, weighted gene co-expression network analysis (WGCNA) identified 37 module genes associated with PUFA-related traits. In general, PLIN1, LIPE, FABP4, LEP, ACACA, ADIPOQ, SCD, PCK2, FASN, PLIN2, LPL, FABP3, THRSP, and ACADVL may have a great impact on PUFA metabolism and lipid deposition. Four SNPs within PLIN2 were significantly associated with FA. Of those, SNP1 (g.287 G>A) was significantly associated with C18:1n9c and MUFA, and SNP4 (g.7807 T>C) was significantly correlated with PUFA (C18:3n3). In addition, the combined genotype of SNP1 (g.287 G>A), SNP3 (g.7664 T>C), and SNP4 (g.7807 T>C) were significantly correlated with C16:1, C17:0, C18:1C6, PUFA (C18:3n3, C22:6n3), and n-6/n-3 PUFA. These results contribute to the knowledge of the biological mechanisms and genetic markers involved in the composition of FA in Hu sheep.

The Biological Hierarchy, Time, and Temporal 'Omics in Evolutionary Biology: A Perspective

Sequencing data-genomics, transcriptomics, epigenomics, proteomics, and metabolomics-have revolutionized biological research, enabling a more detailed study of processes, ranging from subcellular to evolutionary, that drive biological organization. These processes, collectively, are responsible for generating patterns of phenotypic variation and can operate over dramatically different timescales (milliseconds to billions of years). While researchers often study phenotypic variation at specific levels of biological organization to isolate processes operating at that particular scale, the varying types of sequence data, or 'omics, can also provide complementary inferences to link molecular and phenotypic variation to produce an integrated view of evolutionary biology, ranging from molecular pathways to speciation. We briefly describe how 'omics has been used across biological levels and then demonstrate the utility of integrating different types of sequencing data across multiple biological levels within the same study to better understand biological phenomena. However, single-time-point studies cannot evaluate the temporal dynamics of these biological processes. Therefore, we put forward temporal 'omics as a framework that can better enable researchers to study the temporal dynamics of target processes. Temporal 'omics is not infallible, as the temporal sampling regime directly impacts inferential ability. Thus, we also discuss the role the temporal sampling regime plays in deriving inferences about the environmental conditions driving biological processes and provide examples that demonstrate the impact of the sampling regime on biological inference. Finally, we forecast the future of temporal 'omics by highlighting current methodological advancements that will enable temporal 'omics to be extended across species and timescales. We extend this discussion to using temporal multi-omics to integrate across the biological hierarchy to evaluate and link the temporal dynamics of processes that generate phenotypic variation.

Effects of Stocking Density on the Growth Performance, Physiological Parameters, Redox Status and Lipid Metabolism of Micropterus salmoides in Integrated Rice–Fish Farming Systems

Stocking density has been identified as one of the main factors affecting fish growth, welfare and behavior. However, few studies have focused on the effects of stocking density on fish health in integrated rice–fish farming systems. Thus, the aim of this study was to evaluate the effects of different stocking densities on the growth performance, physiological parameters, redox status and lipid metabolism of Micropterus salmoides in an integrated rice–fish farming system. The fish were reared at three densities: low density (LD, 40 g/m³), medium density (MD, 80 g/m³) and high density (HD, 120 g/m³) for 90 days. At the end of the experiment, fish reared in the MD and HD groups showed lower growth performance than those from the LD group. The HD treatment significantly altered the physiological parameters, including glucose and lactate. Meanwhile, the HD treatment induced oxidative stress and lipid peroxidation after 90 days of farming. Furthermore, transcriptomic analysis revealed that HD treatment led to abnormal lipid metabolism. Interestingly, we found the suppression of three key pathways related to lipid metabolism, including the PPAR, insulin and adipocytokine signaling pathways, in the HD group. Overall, our data indicated that the HD treatment inhibited growth and caused physiological responses, oxidative stress and abnormal hepatic lipid metabolism in M. salmoides in an integrated rice–fish farming system.

Genetic architecture and genomic selection of fatty acid composition predicted by Raman spectroscopy in rainbow trout

Background

In response to major challenges regarding the supply and sustainability of marine ingredients in aquafeeds, the aquaculture industry has made a large-scale shift toward plant-based substitutions for fish oil and fish meal. But, this also led to lower levels of healthful n−3 long-chain polyunsaturated fatty acids (PUFAs)—especially eicosapentaenoic (EPA) and docosahexaenoic (DHA) acids—in flesh. One potential solution is to select fish with better abilities to retain or synthesise PUFAs, to increase the efficiency of aquaculture and promote the production of healthier fish products. To this end, we aimed i) to estimate the genetic variability in fatty acid (FA) composition in visceral fat quantified by Raman spectroscopy, with respect to both individual FAs and groups under a feeding regime with limited n-3 PUFAs; ii) to study the genetic and phenotypic correlations between FAs and processing yields- and fat-related traits; iii) to detect QTLs associated with FA composition and identify candidate genes; and iv) to assess the efficiency of genomic selection compared to pedigree-based BLUP selection.

Results

Proportions of the various FAs in fish were indirectly estimated using Raman scattering spectroscopy. Fish were genotyped using the 57 K SNP Axiom™ Trout Genotyping Array. Following quality control, the final analysis contained 29,652 SNPs from 1382 fish. Heritability estimates for traits ranged from 0.03 ± 0.03 (n-3 PUFAs) to 0.24 ± 0.05 (n-6 PUFAs), confirming the potential for genomic selection. n-3 PUFAs are positively correlated to a decrease in fat deposition in the fillet and in the viscera but negatively correlated to body weight. This highlights the potential interest to combine selection on FA and against fat deposition to improve nutritional merit of aquaculture products. Several QTLs were identified for FA composition, containing multiple candidate genes with indirect links to FA metabolism. In particular, one region on Omy1 was associated with n-6 PUFAs, monounsaturated FAs, linoleic acid, and EPA, while a region on Omy7 had effects on n-6 PUFAs, EPA, and linoleic acid. When we compared the effectiveness of breeding programmes based on genomic selection (using a reference population of 1000 individuals related to selection candidates) or on pedigree-based selection, we found that the former yielded increases in selection accuracy of 12 to 120% depending on the FA trait.

Conclusion

This study reveals the polygenic genetic architecture for FA composition in rainbow trout and confirms that genomic selection has potential to improve EPA and DHA proportions in aquaculture species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08062-7.

Deletion of SREBF1, a Functional Bone-Muscle Pleiotropic Gene, Alters Bone Density and Lipid Signaling in Zebrafish

Through a genome-wide analysis of bone mineral density (BMD) and muscle mass, identification of a signaling pattern on 17p11.2 recognized the presence of sterol regulatory element-binding factor 1 (SREBF1), a gene responsible for the regulation of lipid homeostasis. In conjunction with lipid-based metabolic functions, SREBF1 also codes for the protein, SREBP-1, a transcription factor known for its role in adipocyte differentiation. We conducted a quantitative correlational study. We established a zebrafish (ZF) SREBF1 knockout (KO) model and used a targeted customized lipidomics approach to analyze the extent of SREBF1 capabilities. For lipidomics profiling, we isolated the dorsal muscles of wild type (WT) and KO fishes, and we performed liquid chromatography-tandem mass spectrometry screening assays of these samples. In our analysis, we profiled 48 lipid mediators (LMs) derived from various essential polyunsaturated fatty acids to determine potential targets regulated by SREBF1, and we found that the levels of 11,12 epoxyeicosatrienoic acid (11,12-EET) were negatively associated with the number of SREBF1 alleles (P = 0.006 for a linear model). We also compared gene expression between KO and WT ZF by genome-wide RNA-sequencing. Significantly enriched pathways included fatty acid elongation, linoleic acid metabolism, arachidonic acid metabolism, adipocytokine signaling, and DNA replication. We discovered trends indicating that BMD in adult fish was significantly lower in the KO than in the WT population (P < 0.03). These studies reinforce the importance of lipidomics investigation by detailing how the KO of SREBF1 affects both BMD and lipid-signaling mediators, thus confirming the importance of SREBF1 for musculoskeletal homeostasis.

Differences in DNA Methylation Between Disease-Resistant and Disease-Susceptible Chinese Tongue Sole (Cynoglossus semilaevis) Families

DNA methylation, the most widely studied and most well-understood epigenetic modification, has been reported to play crucial roles in diverse processes. Although it has been found that DNA methylation can modulate the expression of immune-related genes in teleosts, a systemic analysis of epigenetic regulation on teleost immunity has rarely been performed. In this research, we employed whole-genome bisulfite sequencing to investigate the genome-wide DNA methylation profiles in select disease-resistant Cynoglossus semilaevis (DR-CS, family 14L006) and disease-susceptible C. semilaevis (DS-CS, family 14L104) against Vibrio harveyi infection. The results showed that following selective breeding, DR-CS had higher DNA methylation levels and different DNA methylation patterns, with 3,311 differentially methylated regions and 6,456 differentially methylated genes. Combining these data with the corresponding transcriptome data, we identified several immune-related genes that exhibited differential expression levels that were modulated by DNA methylation. Specifically, DNA methylation of tumor necrosis factor–like and lipopolysaccharide-binding protein-like was significantly correlated with their expression and significantly contributed to the disease resistance of the selected C. semilaevis family. In conclusion, we suggest that artificial selection for disease resistance in Chinese tongue sole causes changes in DNA methylation levels in important immune-related genes and that these epigenetic changes are potentially involved in multiple immune responses in Chinese tongue sole.