Transcriptomic analysis between Normal and high-intake feeding geese provides insight into adipose deposition and susceptibility to fatty liver in migratory birds

Exploring the dynamic three-dimensional chromatin architecture and transcriptional landscape in goose liver tissues underlying metabolic adaptations induced by a high-fat diet

BACKGROUND

Goose, descendants of migratory ancestors, have undergone extensive selective breeding, resulting in their remarkable ability to accumulate fat in the liver and exhibit a high tolerance for significant energy intake. As a result, goose offers an excellent model for studying obesity, metabolic disorders, and liver diseases in mammals. Although the impact of the three-dimensional arrangement of chromatin within the cell nucleus on gene expression and transcriptional regulation is widely acknowledged, the precise functions of chromatin architecture reorganization during fat deposition in goose liver tissues still need to be fully comprehended.

RESULTS

In this study, geese exhibited more pronounced changes in the liver index and triglyceride (TG) content following the consumption of the high-fat diet (HFD) than mice without significant signs of inflammation. Additionally, we performed comprehensive analyses on 10 goose liver tissues (5 HFD, 5 normal), including generating high-resolution maps of chromatin architecture, conducting whole-genome gene expression profiling, and identifying H3K27ac peaks in the livers of geese and mice subjected to the HFD. Our results unveiled a multiscale restructuring of chromatin architecture, encompassing Compartment A/B, topologically associated domains, and interactions between promoters and enhancers. The dynamism of the three-dimensional genome architecture, prompted by the HFD, assumed a pivotal role in the transcriptional regulation of crucial genes. Furthermore, we identified genes that regulate chromatin conformation changes, contributing to the metabolic adaptation process of lipid deposition and hepatic fat changes in geese in response to excessive energy intake. Moreover, we conducted a cross-species analysis comparing geese and mice exposed to the HFD, revealing unique characteristics specific to the goose liver compared to a mouse. These chromatin conformation changes help elucidate the observed characteristics of fat deposition and hepatic fat regulation in geese under conditions of excessive energy intake.

CONCLUSIONS

We examined the dynamic modifications in three-dimensional chromatin architecture and gene expression induced by an HFD in goose liver tissues. We conducted a cross-species analysis comparing that of mice. Our results contribute significant insights into the chromatin architecture of goose liver tissues, offering a novel perspective for investigating mammal liver diseases.

Multi-omics reveals goose fatty liver formation from metabolic reprogramming

To comprehensively provide insight into goose fatty liver formation, we performed an integrative analysis of the liver transcriptome, lipidome, and amino acid metabolome, as well as peripheral adipose tissue transcriptome analysis using samples collected from the overfed geese and normally fed geese. Transcriptome analysis showed that liver metabolism pathways were mainly enriched in glucolipid metabolism, amino acid metabolism, inflammation response, and cell cycle; peripheral adipose tissue and the liver cooperatively regulated liver lipid accumulation during overfeeding. Liver lipidome patterns obviously changed after overfeeding, and 157 different lipids were yielded. In the liver amino acid metabolome, the level of Lys increased after overfeeding. In summary, this is the first study describing goose fatty liver formation from an integrative analysis of transcriptome, lipidome, and amino acid metabolome, which will provide a whole new dimension to understanding the mechanism of goose fatty liver formation.

Insights into genetic diversity and phenotypic variations in domestic geese through comprehensive population and pan-genome analysis

BACKGROUND

Domestic goose breeds are descended from either the Swan goose (Anser cygnoides) or the Greylag goose (Anser anser), exhibiting variations in body size, reproductive performance, egg production, feather color, and other phenotypic traits. Constructing a pan-genome facilitates a thorough identification of genetic variations, thereby deepening our comprehension of the molecular mechanisms underlying genetic diversity and phenotypic variability.

RESULTS

To comprehensively facilitate population genomic and pan-genomic analyses in geese, we embarked on the task of 659 geese whole genome resequencing data and compiling a database of 155 RNA-seq samples. By constructing the pan-genome for geese, we generated non-reference contigs totaling 612 Mb, unveiling a collection of 2,813 novel genes and pinpointing 15,567 core genes, 1,324 softcore genes, 2,734 shell genes, and 878 cloud genes in goose genomes. Furthermore, we detected an 81.97 Mb genomic region showing signs of genome selection, encompassing the TGFBR2 gene correlated with variations in body weight among geese. Genome-wide association studies utilizing single nucleotide polymorphisms (SNPs) and presence-absence variation revealed significant genomic associations with various goose meat quality, reproductive, and body composition traits. For instance, a gene encoding the SVEP1 protein was linked to carcass oblique length, and a distinct gene-CDS haplotype of the SVEP1 gene exhibited an association with carcass oblique length. Notably, the pan-genome analysis revealed enrichment of variable genes in the "hair follicle maturation" Gene Ontology term, potentially linked to the selection of feather-related traits in geese. A gene presence-absence variation analysis suggested a reduced frequency of genes associated with "regulation of heart contraction" in domesticated geese compared to their wild counterparts. Our study provided novel insights into gene expression features and functions by integrating gene expression patterns across multiple organs and tissues in geese and analyzing population variation.

CONCLUSION

This accomplishment originates from the discernment of a multitude of selection signals and candidate genes associated with a wide array of traits, thereby markedly enhancing our understanding of the processes underlying domestication and breeding in geese. Moreover, assembling the pan-genome for geese has yielded a comprehensive apprehension of the goose genome, establishing it as an indispensable asset poised to offer innovative viewpoints and make substantial contributions to future geese breeding initiatives.

Liver Transcriptome Profiling Identifies Key Genes Related to Lipid Metabolism in Yili Geese

The Yili goose is the only indigenous goose breed that originates from Anser anser in China, known for its adaptability, strong flying ability, and tender meat with a low body lipid content. The liver plays a crucial role in lipid and glucose metabolism, including the intake, secretion, transportation, and storage of fatty acids (FAs). In this study, RNA-sequencing (RNA-seq) technology was performed to analyze the liver differentially expressed genes of Yili geese and their hybrid geese to investigate differences in liver lipid and glucose metabolism. A total of 452 differentially expressed genes (Q-value < 0.05) were identified. Notably, in KEGG enrichment analysis, four pathways (Q-value < 0.05) were enriched to be associated with lipid and glucose metabolism, including the metabolic pathway, PI3K-Akt signaling pathway, glycolysis/gluconeogenesis, and steroid biosynthesis. This study provides insights into potential candidate genes and metabolic pathways that affect the liver lipid metabolism of Yili goose. These findings provide a better understanding of animal liver lipid deposition and metabolism.

Edwardsiella tarda Isolated from a Kidney Mass in a Common Loon (Gavia immer)

A Common Loon (Gavia immer) was found recumbent at Island Beach State Park, New Jersey, US, and euthanized. Necropsy revealed a caseous mass in the kidney, from which bacteria were isolated and phenotypically and molecularly identified as Edwardsiella tarda.

Comparative transcriptome analysis reveals mechanisms of restriction feeding on lipid metabolism in ducks

Presently, excessive fat deposition is the main reason to limit the development of duck industry. In the production, the methods of restricted feeding (RF) were widely used to reduce the lipid deposition of ducks. The liver (L), abdominal adipose (AA), and subcutaneous adipose (SA) were the main tissues of lipid metabolism and deposition of ducks. However, the mechanisms of lipid metabolism and deposition of ducks under RF have not been fully clarified. In this study, in order to better understand the mechanisms of lipid metabolism and deposition in ducks under RF, a total of 120 male Nonghua ducks were randomly divided into a free feeding group (FF, n = 60) and RF group (RF, n = 60), then comparative transcriptomic analysis of L, AA, and SA between FF (n = 3) and RF (n = 3) ducks was performed at 56 d of age. Phenotypically, L, AA, and SA index of FF group was higher than that in RF group. There were 279, 390, and 557 differentially expressed genes (DEGs) in L, AA, and SA. Functional enrichment analysis revealed that ECM-receptor interaction and metabolic pathways were significantly enriched in L, AA, and SA. Lipid metabolism-related pathways including fatty acid metabolism, unsaturated fatty acid synthesis, and steroidogenesis were significantly enriched in AA and SA. Moreover, through integrated analysis weighted gene coexpression network (WGCNA) and protein-protein interaction network, 10 potential candidate genes involved in the ECM-receptor interaction and lipid metabolism pathways were identified, including 3-hydroxy-3-methylglutaryl-CoA synthase 2 (HMGCS2), aldolase B (ALDOB), formimidoyltransferase cyclodeaminase(FTCD), phosphoenolpyruvate carboxykinase 1 (PCK1), tyrosine aminotransferase (TAT), stearoyl-CoA desaturase (SCD), squalene epoxidase (SQLE), phosphodiesterase 4B (PDE4B), choline kinase A (CHKA), and elongation of very-long-chain fatty acids-like 2 (ELOVL2), which could play a key role in lipid metabolism and deposition of ducks under RF. Our study reveals that the liver might regulate the lipid metabolism of abdominal adipose and subcutaneous adipose through ECM-receptor interaction and metabolic pathways (fatty acid metabolism, unsaturated fatty acid synthesis, and steroid synthesis), thus to reduce the lipid deposition of ducks under RF. These results provide novel insights into the avian lipid metabolism and will help better understand the underlying molecular mechanisms.

Arctic introgression and chromatin regulation facilitated rapid Qinghai-Tibet Plateau colonization by an avian predator

The Qinghai-Tibet Plateau (QTP), possesses a climate as cold as that of the Arctic, and also presents uniquely low oxygen concentrations and intense ultraviolet (UV) radiation. QTP animals have adapted to these extreme conditions, but whether they obtained genetic variations from the Arctic during cold adaptation, and how genomic mutations in non-coding regions regulate gene expression under hypoxia and intense UV environment, remain largely unknown. Here, we assemble a high-quality saker falcon genome and resequence populations across Eurasia. We identify female-biased hybridization with Arctic gyrfalcons in the last glacial maximum, that endowed eastern sakers with alleles conveying larger body size and changes in fat metabolism, predisposing their QTP cold adaptation. We discover that QTP hypoxia and UV adaptations mainly involve independent changes in non-coding genomic variants. Our study highlights key roles of gene flow from Arctic relatives during QTP hypothermia adaptation, and cis-regulatory elements during hypoxic response and UV protection.

Regulation of the cecal microbiota community and the fatty liver deposition by the addition of brewers’ spent grain to feed of Landes geese

The effects of brewers’ spent grain (BSG) diets on the fatty liver deposition and the cecal microbial community were investigated in a total of 320 healthy 5-day-old Landes geese. These geese were randomly and evenly divided into 4 groups each containing 8 replicates and 10 geese per replicate. These four groups of geese were fed from the rearing stage (days 5–60) to the overfeeding stage (days 61–90). The Landes geese in group C (control) were fed with basal diet (days 5–90); group B fed first with basal diet in the rearing stage and then basal diet + 4% BSG in the overfeeding stage; group F first with basal diet + 4% BSG during the rearing stage and then basal diet in the overfeeding stage; and group W with basal diet + 4% BSG (days 5–90). The results showed that during the rearing stage, the body weight (BW) and the average daily gain (ADG) of Landes geese were significantly increased in groups F and W, while during the overfeeding stage, the liver weights of groups W and B were significantly higher than that of group C. The taxonomic structure of the intestinal microbiota revealed that during the overfeeding period, the relative abundance of Bacteroides in group W was increased compared to group C, while the relative abundances of Escherichia–Shigella and prevotellaceae_Ga6A1_group were decreased. Results of the transcriptomics analysis showed that addition of BSG to Landes geese diets altered the expression of genes involved in PI3K-Akt signaling pathway and sphingolipid metabolism in the liver. Our study provided novel experimental evidence based on the cecal microbiota to support the application of BSG in the regulation of fatty liver deposition by modulating the gut microbiota in Landes geese.

Comparative Transcriptome Analysis Provides Novel Insights into the Effect of Lipid Metabolism on Laying of Geese

Simple Summary

The importance of lipid metabolism in the egg production of poultry has been widely reported. Meanwhile, geese have lower egg production and unique lipid metabolism patterns compared with chicken and duck. It is of great significance to further improve egg laying performance to explore the differences of fat metabolism and the molecular mechanisms in geese with different egg laying performance. This study compared the phenotypic differences of liver and abdominal fat, as well as the transcriptome level differences of liver, abdominal fat, and ovarian stroma among high-, low-, and no-egg production groups. The results reveal that lipid metabolism regulated by the circadian rhythm of the liver may directly or indirectly affect ovarian function through the inflammation and hormone secretion of abdominal fat.

Abstract

The lower egg production of geese (20~60 eggs per year) compared with chicken and duck limits the development of the industry, while the yolk weight and fatty liver susceptibility of geese was higher than that of other poultry. Therefore, the relationship between lipid metabolism and the laying performance of geese remains to be explored. Phenotypically, we observed that the liver fat content of the high-, low-, and no-egg production groups decreased in turn, while the abdominal fat weight increased in turn. For transcriptional regulation, the KEGG pathways related to lipid metabolism were enriched in all pairwise comparisons of abdominal fat and liver through functional analysis. However, some KEGG pathways related to inflammation and the circadian rhythm pathway were enriched by DEGs only in abdominal fat and the liver, respectively. The DEGs in ovarian stroma among different groups enriched some KEGG pathways related to ovarian steroidogenesis and cell adhesion. Our research reveals that lipid metabolism regulated by the circadian rhythm of the liver may directly or indirectly affect ovarian function through the inflammation and hormone secretion of abdominal fat. These results offer new insights into the regulation mechanisms of goose reproductive traits.

Glucose participates in the formation of goose fatty liver by regulating the expression of miRNA-33/CROT

Glucose oversupply promotes formation of fatty liver, and fatty liver is usually accompanied with hyperglycemia. However, the mechanism by which glucose promotes formation of fatty liver is not very clear. In this study, fatty liver was successfully induced in Landes goose by 19 days of overfeeding with corn-based feed, the overfed geese had a significantly higher level of blood glucose than the normally fed geese (control group). In goose primary liver cells, high level of glucose promoted fat deposition and induced the expression of SREBF2(or SREBP2), a key regulator of lipid metabolism, and its intronic gene, miR-33. Moreover, overexpression of miRNA-33(miR-33) promotes lipid accumulation in goose primary liver cells. Consistently, miR-33 inhibitor suppressed glucose induced lipid accumulation in liver cells. Interestingly, the relative abundance of miR-33 in goose fatty liver was significantly higher than that in normal liver, while the relative mRNA and protein abundances of CROT, the target gene of miR-33, in goose fatty liver were significantly lower than those in goose normal liver. Taken together, these findings suggest that miR-33 mediates glucose promotion of lipid accumulation in goose primary liver cells, and that glucose participates in formation of goose fatty liver by regulating the expression of miR-33/CROT.

A Prognostic Model for Brain Glioma Patients Based on 9 Signature Glycolytic Genes

Objective

To screen glycolytic genes linked to the glioma prognosis and construct the prognostic model.

Methods

The relevant data of glioma were downloaded from TCGA and GTEx databases. GSEA of glycolysis-related pathways was carried out, and enriched differential genes were extracted. Screening out prognostic-related genes with conspicuous significance and construction of the prognostic model were conducted by multivariate Cox regression analysis and Lasso regression analysis. The model was evaluated, and cBioPortal was used to analyze the mutation of the model gene. The expression of the model gene in tumor and normal colon tissue was analyzed. The model was used to evaluate the prognosis of patients in different groups to verify the applicability of the model.

Results

339 differentially glycolytic-related genes were enriched in REACTOME_GLYCOLYSIS, GLYCOLYTIC_PROCESS, HALLMARK_GLYCOLYSIS, and other pathways. We obtained 9 key prognostic genes and constructed the prognostic evaluation model. The 3-year AUC values of the ROC curve display model are greater than 0.75, which indicates that the accuracy of the model is good. The relation of age and risk score to prognosis is shown by univariate and multivariate Cox analysis. The expression of SRD5A3, MDH2, and B3GAT3 genes was significantly upregulated in the tumor tissues, while the HDAC4 and G6PC2 genes were downregulated. The mutation rate of MDH2 and HDAC4 genes was the highest. This model could effectively distinguish the risk of poor prognosis of patients in any age stage.

Conclusion

The prognostic assessment models based on glycolysis-related nine-gene signature could accurately predict the prognosis of patients with GBM.

Identification and characterization of the adipogenesis in intramuscular and subcutaneous adipocytes of the goose (Anser cygnoides)

Chinese geese are domesticated from wild swan (Anser cygnoides), which have maintained a strong capacity of fat deposit. Fat mainly distributes subcutaneous, abdominal, intermuscular or intramuscular in poultry, and they display some special physiological and biochemical characteristics in different parts. This study aimed to characterize the adipogenesis in intramuscular (IM) and subcutaneous (SC) adipocytes of the goose. Here, IM and SC preadipocytes were isolated from the 26-day-old Yangzhou goose embryos, and subsequently induced them to differentiate into mature adipocytes. The results showed that SC preadipocytes grew a little faster than IM preadipocytes during the logarithmic multiplication phase (p < 0.05). Meanwhile, SC adipocytes accumulated more lipid than IM adipocytes during the differentiation process in vitro (p < 0.01). In addition, the expression level of key genes involved in adipogenesis, including peroxisome proliferator activated receptor γ/α (PPARγ/α), CCAAT/enhancer binding protein-α/β (C/EBPα/β), adipocyte fatty acid binding protein 4 (FABP4), and lipoprotein lipase (LPL) were detected. PPARγ, C/EBPα, FABP4, and LPL, were predominantly expressed in SC adipocytes, whereas C/EBPβ was highly expressed in IM adipocytes. Taken together, these results demonstrated that SC preadipocytes tended to grow faster and accumulate more lipid than IM adipocytes, and show greater potential for adipogenesis.

Fasting and overfeeding affect the expression of the immunity- or inflammation-related genes in the liver of poultry via endogenous retrovirus

It is known that nutrition and immunity are connected, but the mechanism is not very clear. Endogenous retroviruses (ERV) account for 8 to 10% of the human and mouse genomes and play an important role in some biological processes of animals. Recent studies indicate that the activation of ERV can affect the expression of the immunity- or inflammation-related genes, and the activities of ERV are subjected to regulation of many factors including nutritional factors. Therefore, we hypothesize that nutritional status can affect the expression of the immunity- or inflammation-related genes via ERV. To verify this hypothesis, the nutritional status of animals was altered by fasting or overfeeding, and the expression of intact ERV (ERVK18P, ERVK25P) and immunity- or inflammation-related genes (DDX41, IFIH1, IFNG, IRF7, STAT3) in the liver was determined by quantitative PCR, followed by overexpressing ERVK25P in goose primary hepatocytes and determining the expression of the immunity- or inflammation-related genes. The data showed that compared with the control group (no fasting), the expression of ERV and the immunity- or inflammation-related genes was increased in the liver of the fasted chickens but decreased in the liver of the fasted geese. Moreover, compared with the control group (routinely fed), the expression of ERV and the immunity- or inflammation-related genes was increased in the liver of the overfed geese. In addition, overexpression of ERVK25P in goose primary hepatocytes can induce the expression of the immunity- or inflammation-related genes. In conclusion, these findings suggest that ERV mediate the effects of fasting and overfeeding on the expression of the immunity- or inflammation-related genes, the mediation varied with poultry species, and ERV and the immunity- or inflammation-related genes may be involved in the development of goose fatty liver. This study provides a potential mechanism for the connection between nutrition and immunity.

Pacific Biosciences assembly with Hi-C mapping generates an improved, chromosome-level goose genome

BACKGROUND

The domestic goose is an economically important and scientifically valuable waterfowl; however, a lack of high-quality genomic data has hindered research concerning its genome, genetics, and breeding. As domestic geese breeds derive from both the swan goose (Anser cygnoides) and the graylag goose (Anser anser), we selected a female Tianfu goose for genome sequencing. We generated a chromosome-level goose genome assembly by adopting a hybrid de novo assembly approach that combined Pacific Biosciences single-molecule real-time sequencing, high-throughput chromatin conformation capture mapping, and Illumina short-read sequencing.

FINDINGS

We generated a 1.11-Gb goose genome with contig and scaffold N50 values of 1.85 and 33.12 Mb, respectively. The assembly contains 39 pseudo-chromosomes (2n = 78) accounting for ∼88.36% of the goose genome. Compared with previous goose assemblies, our assembly has more continuity, completeness, and accuracy; the annotation of core eukaryotic genes and universal single-copy orthologs has also been improved. We have identified 17,568 protein-coding genes and a repeat content of 8.67% (96.57 Mb) in this genome assembly. We also explored the spatial organization of chromatin and gene expression in the goose liver tissues, in terms of inter-pseudo-chromosomal interaction patterns, compartments, topologically associating domains, and promoter-enhancer interactions.

CONCLUSIONS

We present the first chromosome-level assembly of the goose genome. This will be a valuable resource for future genetic and genomic studies on geese.

Exploring the genetic basis of fatty liver development in geese

Although geese possess an adaptive physiological capacity for lipid storage, few candidate genes contributing to this ability are characterised. By comparing the genomes of individuals with extremely high and low fatty liver weights (FLW), candidate genes were identified, including ARAP2, GABRE, and IL6. Single-nucleotide polymorphisms in or near these genes were significantly (p < 0.05) associated with carcass traits (FLW) and biochemical indexes (very-low-density lipoprotein and N-terminal procollagen III), suggesting contribution to trait variation. A common variant at the 5'-end of LCORL explained ~ 18% and ~ 26% of the phenotypic variance in body weight with/without overfeeding and had significant effects on FLW (p < 0.01). ZFF36L1, ARHGEF1 and IQCJ, involved in bile acid metabolism, blood pressure, and lipid concentration modulation, were also identified. The presence of highly divergent haplotypes within these genes suggested involvement in protection against negative effects from excessive lipids in the liver or circulatory system. Based on this and transcriptomic data, we concluded that geese hepatosteatosis results from severe imbalance between lipid accumulation and secretion, comparable to human non-alcohol fatty liver disease but involving other genes. Our results provided valuable insights into the genesis of geese fatty liver and detected potential target genes for treatment of lipid-related diseases.