WGCNA-Based Identification of Hub Genes and Key Pathways Involved in Nonalcoholic Fatty Liver Disease

NDUFA9 and its crotonylation modification promote browning of white adipocytes by activating mitochondrial function in mice

Protein crotonylation plays a role in regulating cellular metabolism, gene expression, and other biological processes. NDUFA9 (NADH dehydrogenase [ubiquinone] 1 alpha subcomplex subunit 9) is closely associated with the activity and function of mitochondrial respiratory chain complex I. Mitochondrial function and respiratory chain are closely related to browning of white adipocytes, it's speculated that NDUFA9 and its crotonylation are associated with browning of white adipocytes. Firstly, the effect of NDUFA9 on white adipose tissue was verified in white fat browning model mice, and it was found that NDUFA9 promoted mitochondrial respiration, thermogenesis, and browning of white adipose tissue. Secondly, in cellular studies, it was discovered that NDUFA9 facilitated browning of white adipocytes by enhancing mitochondrial function, mitochondrial complex I activity, ATP synthesis, and mitochondrial respiration. Again, the level of NDUFA9 crotonylation was increased by treating cells with vorinostat (SAHA)+sodium crotonate (NaCr) and overexpressing NDUFA9, it was found that NDUFA9 crotonylation promoted browning of white adipocytes. Meanwhile, the acetylation level of NDUFA9 was increased by treating cells with SAHA+sodium acetate (NaAc) and overexpressing NDUFA9, the assay revealed that NDUFA9 acetylation inhibited white adipocytes browning. Finally, combined with the competitive relationship between acetylation and crotonylation, it was also demonstrated that NDUFA9 crotonylation promoted browning of white adipocytes. Above results indicate that NDUFA9 and its crotonylation modification promote mitochondrial function, which in turn promotes browning of white adipocytes. This study establishes a theoretical foundation for the management and intervention of obesity, which is crucial in addressing obesity and related medical conditions in the future.

Curculigoside, a traditional Chinese medicine monomer, ameliorates oxidative stress in Alzheimer's disease mouse model via suppressing ferroptosis

Alzheimer's disease (AD) is a neurodegenerative disorder where oxidative stress, induced by ferroptosis, has been linked to neuronal damage and cognitive deficits. The objective of this study is to investigate if the potential therapeutic agent, Curculigoside (CUR), could ameliorate AD by inhibiting ferroptosis. The potential therapeutic targets, such as GPX4 and SLC7A11, were identified using weighted gene co-expression network analysis (WGCNA). Concurrently, CUR was also screened against these potential targets using various analytical methods. For the in vivo studies, intragastric administration of CUR significantly ameliorated cognitive impairment in AD model mice induced by scopolamine and okadaic acid (OA). In vitro, CUR protected neuronal cells by altering the levels of ferroptosis-related specific markers in OA and scopolamine-induced neurotoxicity. The administration of CUR through intragastric route significantly reduced the levels of AD-promoting factors (such as Aβ1-42, p-tau) and ferroptosis-promoting factors in the hippocampus and cortex of AD mice. Furthermore, CUR up-regulated the expression of GPX4 and decreased the expression of SLC7A11 in the ferroptosis signaling pathway, thereby increasing the ratio of glutathione (GSH)/oxidized glutathione (GSSG) in vivo and vitro. In conclusion, the cumulative results suggest that the natural compound CUR may serve as a promising therapeutic agent to ameliorate AD by inhibiting ferroptosis.

Altered splicing factor and alternative splicing events in a mouse model of diet- and polychlorinated biphenyl-induced liver disease

Non-alcoholic fatty liver disease (NAFLD) is associated with human environmental exposure to polychlorinated biphenyls (PCBs). Alternative splicing (AS) is dysregulated in steatotic liver disease and is regulated by splicing factors (SFs) and N-6 methyladenosine (m6A) modification. Here integrated analysis of hepatic mRNA-sequencing data was used to identify differentially expressed SFs and differential AS events (ASEs) in the livers of high fat diet-fed C57BL/6 J male mice exposed to Aroclor1260, PCB126, Aroclor1260 + PCB126, or vehicle control. Aroclor1260 + PCB126 co-exposure altered 100 SFs and replicate multivariate analysis of transcript splicing (rMATS) identified 449 ASEs in 366 genes associated with NAFLD pathways. These ASEs were similar to those resulting from experimental perturbations in m6A writers, readers, and erasers. These results demonstrate specific hepatic SF and AS regulatory mechanisms are disrupted by HFD and PCB exposures, contributing to the expression of altered isoforms that may play a role in NAFLD progression to NASH.

Identification and validation of disulfidptosis-associated molecular clusters in non-alcoholic fatty liver disease

Objective: Non-alcoholic fatty liver disease (NAFLD) is the most prevalent liver disease in the world, and its pathogenesis is not fully understood. Disulfidptosis is the most recently reported form of cell death and may be associated with NAFLD progression. Our study aimed to explore the molecular clusters associated with disulfidptosis in NAFLD and to construct a predictive model. Methods: First, we analyzed the expression profile of the disulfidptosis regulators and immune characteristics in NAFLD. Using 104 NAFLD samples, we investigated molecular clusters based on differentially expressed disulfidptosis-related genes, along with the related immune cell infiltration. Cluster-specific differentially expressed genes were then identified by using the WGCNA method. We also evaluated the performance of four machine learning models before choosing the optimal machine model for diagnosis. Nomogram, calibration curves, decision curve analysis, and external datasets were used to confirm the prediction effectiveness. Finally, the expression levels of the biomarkers were assessed in a mouse model of a high-fat diet. Results: Two differentially expressed DRGs were identified between healthy and NAFLD patients. We revealed the expression profile of DRGs in NAFLD and the correlation with 22 immune cells. In NAFLD, two clusters of molecules connected to disulfidptosis were defined. Significant immunological heterogeneity was shown by immune infiltration analysis among the various clusters. A significant amount of immunological infiltration was seen in Cluster 1. Functional analysis revealed that Cluster 1 differentially expressed genes were strongly linked to energy metabolism and immune control. The highest discriminatory performance was demonstrated by the SVM model, which had a higher area under the curve, relatively small residual and root mean square errors. Nomograms, calibration curves, and decision curve analyses were used to show how accurate the prediction of NAFLD was. Further analysis revealed that the expression of three model-related genes was significantly associated with the level of multiple immune cells. In animal experiments, the expression trends of DDO, FRK and TMEM19 were consistent with the results of bioinformatics analysis. Conclusion: This study systematically elucidated the complex relationship between disulfidptosis and NAFLD and developed a promising predictive model to assess the risk of disease in patients with disulfidptosis subtypes and NAFLD.

Effects of Methyl Sulfonyl Methane and Selenium Yeast on Fatty Liver Syndrome in Laying Hens and Their Biological Mechanisms

The purpose of this study was to explore the effects of MSM and Se-Y on FLS in laying hens during the late peak laying period and the underlying biological mechanisms. Therefore 240 55-week-old Jing-fen No. 6 laying hens were randomly divided into five groups, with eight replicates in each group and six laying hens in each replicate. The hens were fed a basal diet (Control) and diets supplemented with 350 and 700 mg/kg MSM and 25 and 50 mg/kg Se-Y, respectively, for four weeks. The results showed that MSM and Se-Y had no significant effects on the performance of laying hens. With the increasing dosage of MSM and Se-Y, the symptoms of liver steatosis in laying hens were reduced, and MSM and Se-Y could significantly reduce the content of malondialdehyde (MDA) in serum and liver (p < 0.05) and increase the contents of total superoxide dismutase (T-SOD) and glutathione peroxidase (GPX) in serum and liver (p < 0.05). The RNA-seq results showed that 700 mg/kg MSM significantly downregulated the expression levels of the ATP5I, ATP5G1, CYCS, and UQCRQ genes in the liver, and 50 mg/kg Se-Y significantly downregulated the expression levels of MAPK10, SRC, BMP2, and FGF9 genes in the liver. In conclusion, dietary supplementation with MSM and Se-Y can effectively reduce the FLS of laying hens in the late peak laying period and increase their antioxidant capacity. The underlying biological mechanism may be related to the downregulation of genes involved in liver oxidative phosphorylation and inflammation-related pathways.

Identification and Evaluation of Hub Long Noncoding RNAs and mRNAs in High Fat Diet Induced Liver Steatosis

Nonalcoholic fatty liver disease (NAFLD) is considered the most prevalent chronic liver disease, but the understanding of the mechanism of NAFLD is still limited. The aim of our study was to explore hub lncRNAs and mRNAs and pathological processes in high-fat diet (HFD)-induced and lycopene-intervened liver steatosis. We analyzed the gene profiles in the GSE146627 dataset from the Gene Expression Omnibus (GEO) database to identify differentially expressed lncRNAs and mRNAs, and we constructed coexpression networks based on weighted gene coexpression network analysis (WGCNA). The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases were utilized for functional enrichment analysis. We found that the turquoise, blue, brown, yellow, green, and black modules were significantly correlated with NAFLD. Functional enrichment analysis revealed that some hub lncRNAs (Smarca2, Tacc1, Flywch1, and Mef2c) might be involved in the regulation of the inflammatory and metabolic pathways (such as TNF signaling, metabolic, mTOR signaling, MAPK signaling, and p53 signaling pathways) in NAFLD. The establishment of an NAFLD mouse model confirmed that lycopene supply attenuated hepatic steatosis in HFD-induced NAFLD. Our analysis revealed that the inflammatory and metabolic pathways may be crucially involved in the pathogenesis of NAFLD, and hub lncRNAs provide novel biomarkers, therapeutic ideas, and targets for NAFLD. Moreover, lycopene has the potential to be a phytochemical for the prevention of HFD-induced liver steatosis.

Identification and Verification of Potential Biomarkers in Renal Ischemia-Reperfusion Injury by Integrated Bioinformatic Analysis

Background

Renal ischemia-reperfusion injury (RIRI) plays an important role in the poor prognosis of patients with renal transplants. However, the potential targets and mechanism of IRI are still unclear.

Method

Differential gene expression (DEG) analysis and weighted correlation network analysis (WGCNA) were performed on the GSE27274 dataset. Pathway enrichment analysis on the DEGs was performed. To identify the hub DEGs, we constructed a protein-protein interaction (PPI) network. Finally, the hub genes were verified, and candidate drugs were screened from the DsigDB database.

Results

A hundred DEGs and four hub genes (Atf3, Psmb6, Psmb8, and Psmb10) were screened out. Pathway enrichment analysis revealed that 100 DEGs were mainly enriched in apoptosis and the TNF signaling pathway. The four hub genes were verified in animal models and another dataset (GSE148420). Thereafter, a PPI network was used to identify the four hub genes (Atf3, Psmb6, Psmb8, and Psmb10). Finally, eight candidate drugs were identified as potential drugs.

Conclusion

Three hub genes (Psmb6, Psmb8, and Psmb10) were associated with RIRI and could be potential novel biomarkers for RIRI.

Integrative Multiomics and Regulatory Network Analyses Uncovers the Role of OAS3, TRAFD1, miR-222-3p, and miR-125b-5p in Hepatitis E Virus Infection

The hepatitis E virus (HEV) is a long-ignored virus that has spread globally with time. It ranked 6th among the top risk-ranking viruses with high zoonotic spillover potential; thus, considering its viral threats is a pressing priority. The molecular pathophysiology of HEV infection or the underlying cause is limited. Therefore, we incorporated an unbiased, systematic methodology to get insights into the biological heterogeneity associated with the HEV. Our study fetched 93 and 2016 differentially expressed genes (DEGs) from chronic HEV (CHEV) infection in kidney-transplant patients, followed by hub module selection from a weighted gene co-expression network (WGCN). Most of the hub genes identified in this study were associated with interferon (IFN) signaling pathways. Amongst the genes induced by IFNs, the 2'-5'-oligoadenylate synthase 3 (OAS3) protein was upregulated. Protein-protein interaction (PPI) modular, functional enrichment, and feed-forward loop (FFL) analyses led to the identification of two key miRNAs, i.e., miR-222-3p and miR-125b-5p, which showed a strong association with the OAS3 gene and TRAF-type zinc finger domain containing 1 (TRAFD1) transcription factor (TF) based on essential centrality measures. Further experimental studies are required to substantiate the significance of these FFL-associated genes and miRNAs with their respective functions in CHEV. To our knowledge, it is the first time that miR-222-3p has been described as a reference miRNA for use in CHEV sample analyses. In conclusion, our study has enlightened a few budding targets of HEV, which might help us understand the cellular and molecular pathways dysregulated in HEV through various factors. Thus, providing a novel insight into its pathophysiology and progression dynamics.

Identification of polyunsaturated fatty acids related key modules and genes in metabolic dysfunction-associated fatty liver disease using WGCNA analysis

Polyunsaturated fatty acids (PUFAs) play important roles in the aetiology and pathogenesis of metabolic dysfunction-associated fatty liver disease (MAFLD). However, the underlying molecular mechanisms are not understood. We analysed a public GEO dataset, GSE89632, to identify differentially expressed genes (DEGs) in MAFLD. Weighted gene coexpression network analysis (WGCNA) was used to reveal the core gene regulation network and to explore the PUFA-related hub genes in MAFLD. We experimentally verified these genes by quantitative reverse transcription PCR in high-fat diet (HFD)-fed mice. A total of 286 common DEGs (89 upregulated; 197 downregulated), mostly related to inflammatory and immune responses, were identified. Six modules were constructed using WGCNA, and 2 modules showed significant correlations with PUFAs. After combining these 2 modules with DEGs, the top 10 hub genes were identified. We further established a MAFLD mouse model with liver steatosis, as proved by HE and Oil Red O staining. Of the hub genes, ADAM metallopeptidase with thrombospondin type 1 motif 1 (adamts1) (p = 0.005) and transforming growth factor β3 (tgfβ3) (p < 0.001) showed significantly lower mRNA expression in MAFLD in vivo. adamts1 and tgfβ3 bridged PUFAs and MAFLD, which might be potential causative genes and therapeutic targets of MAFLD.

Identification of DNA methylation-regulated differentially expressed genes in RA by integrated analysis of DNA methylation and RNA-Seq data

OBJECTIVE

To identify novel DNA methylation-regulated differentially expressed genes (MeDEGs) in RA by integrated analysis of DNA methylation and RNA-Seq data.

METHODS

The transcription and DNA methylation profiles of 9 RA and 15 OA synovial tissue were generated by RNA-Seq and Illumina 850K DNA methylation BeadChip. Gene set enrichment analysis (GSEA) and Weighted gene co-expression network analysis (WGCNA) were used to analyze methylation-regulated expressed genes by R software. The differentially expressed genes (DEGs), differentially methylated probes (DMPs), differentially methylated genes (DMGs) were analyzed by DESeq and ChAMP R package. The functional correlation of MeDEGs was analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). The protein-protein interaction (PPI) network of MeDEGs was constructed by STRING and Reactome FI Cytoscape Plugin. Correlation analysis between methylation level and mRNA expression was conducted with R software.

RESULTS

A total of 17,736 genes, 25,578 methylated genes and 755,852 methylation probes were detected. A total of 16,421 methylation-regulated expressed genes were obtained. The GSEA showed that these genes are associated with activation of immune response, adaptive immune response, Inflammatory response in C5 (ontology gene sets). For KEGG analysis, these genes are associated with rheumatoid arthritis, NF-kappa B signaling pathway, T cell receptor signaling pathway. The WGCNA showed that the turquoise module exhibited the strongest correlation with RA (R = 0.78, P = 1.27 × 10^- 05), 660 genes were screened in the turquoise module. A total of 707 MeDEGs were obtained. GO analysis showed that MeDEGs were enriched in signal transduction, cell adhesion for BP, enriched in plasma membrane, integral component of membrane for CC, and enriched in identical protein binding, calcium ion binding for MF. The KEGG pathway analysis showed that the MeDEGs were enriched in calcium signaling pathway, T cell receptor signaling pathway, NF-kappa B signaling pathway, Rheumatoid arthritis. The PPI network containing 706 nodes and 882 edges, and the enrichment p value < 1.0 × 10^- 16. With Cytoscape, based on the range of more than 10 genes, a total of 8 modules were screened out. Spearman correlation analysis showed RGS1(cg10718027), RGS1(cg02586212), RGS1(cg10861751) were significantly correlated with RA.

CONCLUSIONS

RGS1 can be used as novel methylated biomarkers for RA.

Fucose as a potential therapeutic molecule against the immune-mediated inflammation in IgA nepharopathy: An unrevealed link

Background

IgA nephropathy (IgAN) is an autoimmune disease that affects people of any age and is an important cause of end-stage renal disease. However, the pathogenesis and pathophysiology of IgAN is not clear. This article aimed to explore the immune-mediated inflammation and genetic mechanisms in IgAN.

Methods

The transcriptome sequencing data of IgAN glomeruli in the Gene Expression Omnibus database were downloaded. Single-sample gene set enrichment analysis was used to estimate the immune microenvironment of the merged microarray data and GSE141295. IgAN samples were divided into two clusters by cluster analysis. “limma” and “DEseq2” package in R were used to identify differentially expressed genes (DEGs). The weighted gene co-expression network analysis (WGCNA) was used to identify the co-expression modules related to inflammation in IgAN. R software package “clusterProfiler” was used for enrichment analysis, whereas Short Time-Series Expression Miner (STEM) analysis was used to identify the trend of gene expression. Machine-learn (ML) was performed using the shiny app. Finally, Drug Signatures Database (DSigDB) was used to identify potential molecules for treating IgAN.

Results

The infiltration of macrophages in IgAN glomeruli was increased, whereas CD4+ T cells, especially inducedregulatory T cells (iTregs) were decreased. A total of 1,104 common DEGs were identified from the merged data and GSE141295. Brown module was identified to have the highest inflammatory correlation with IgAN using WGCNA, and 15 hub genes were screened from this module. Among these 15 hub genes, 14 increased with the severity of IgAN inflammation based on STEM analysis. Neural network (nnet) is considered as the best model to predict the severity of IgAN. Fucose identified from DSigDB has a potential biological activity to treat IgAN.

Conclusion

The increase of macrophages and the decrease of iTregs in glomeruli represent the immune-mediated inflammation of IgAN, and fucose may be a potential therapeutic molecule against IgAN because it affects genes involved in the severe inflammation of IgAN.

Genome-Wide Identification of Brassicaceae Hormone-Related Transcription Factors and Their Roles in Stress Adaptation and Plant Height Regulation in Allotetraploid Rapeseed

Phytohormone-related transcription factors (TFs) are involved in regulating stress responses and plant growth. However, systematic analysis of these TFs in Brassicaceae is limited, and their functions in stress adaptation and plant height (PH) regulation remain unclear. In this study, 2115 hormone-related TFs were identified in nine Brassicaceae species. Specific domains were found in several Brassicaceae hormone-related TFs, which may be associated with diverse functions. Syntenic analysis indicated that expansion of these genes was mainly caused by segmental duplication, with whole-genome duplication occurring in some species. Differential expression analysis and gene co-expression network analysis identified seven phytohormone-related TFs (BnaWRKY7, 21, 32, 38, 52, BnaGL3-4, and BnaAREB2-5) as possible key genes for cadmium (Cd) toxicity, salinity stress, and potassium (K) and nitrogen (N) deficiencies. Furthermore, BnaWRKY42 and BnaARR21 may play essential roles in plant height. Weighted gene co-expression network analysis (WGCNA) identified 15 phytohormone-related TFs and their potential target genes regulating stress adaptation and plant height. Among the above genes, BnaWRKY56 and BnaWRKY60 responded to four different stresses simultaneously, and BnaWRKY42 was identified in two dwarf rapeseeds. In summary, several candidate genes for stress resistance (BnaWRKY56 and BnaWRKY60) and plant height (BnaWRKY42) were identified. These findings should help elucidate the biological roles of Brassicaceae hormone-related TFs, and the identified candidate genes should provide a genetic resource for the potential development of stress-tolerant and dwarf oilseed plants.

Identification of Hub Genes Associated With Non-alcoholic Steatohepatitis Using Integrated Bioinformatics Analysis

Background and aims: As a major cause of liver disease worldwide, non-alcoholic fatty liver disease (NAFLD) comprises non-alcoholic fatty liver (NAFL) and non-alcoholic steatohepatitis (NASH). Due to the high prevalence and poor prognosis of NASH, it is critical to understand its mechanisms. However, the etiology and mechanisms remain largely unknown. In addition, the gold standard for the diagnosis of NASH is liver biopsy, which is an invasive procedure. Therefore, there is a pressing need to develop noninvasive tests for NASH diagnosis. The goal of the study is to discover key genes involved in NASH development and investigate their value as noninvasive biomarkers.

Methods: The Gene Expression Omnibus (GEO) database was used to obtain two datasets encompassing NASH patients and healthy controls. We used weighted gene co-expression network analysis (WGCNA) and differential expression analysis in order to investigate the association between gene sets and clinical features, as well as to discover co-expression modules. A protein-protein interaction (PPI) network was created to extract hub genes. The results were validated using another publicly available dataset and mice treated with a high-fat diet (HFD) and carbon tetrachloride (CCl4).

Results: A total of 24 differentially co-expressed genes were selected by WGCNA and differential expression analysis. KEGG analysis indicated most of them were enriched in the focal adhesion pathway. GO analysis showed these genes were mainly enriched in circadian rhythm, aging, angiogenesis and response to drug (biological process), endoplasmic reticulum lumen (cellular component), and protein binding (molecular function). As a result, eight genes (JUN, SERPINE1, GINS2, TYMS, HMMR, IGFBP2, BIRC3, TNFRSF12A) were identified as hub genes. Finally, three genes were found significantly changed in both the validation dataset and the mouse model.

Conclusion: Our research discovered genes that have the potential to mediate the process of NASH and might be useful diagnostic biomarkers for the disorder.

Integrative Analysis of Nanopore and Illumina Sequencing Reveals Alternative Splicing Complexity in Pig Longissimus Dorsi Muscle

Alternative splicing (AS) is a key step in the post-transcriptional regulation of gene expression that can affect intramuscular fat (IMF). In this study, longissimus dorsi muscles from 30 pigs in high- and low- IMF groups were used to perform Oxford Nanopore Technologies (ONT) full-length sequencing and Illumina strand-specific RNA-seq. A total of 43,688 full-length transcripts were identified, with 4,322 novel genes and 30,795 novel transcripts. Using AStalavista, a total of 14,728 AS events were detected in the longissimus dorsi muscle. About 17.79% of the genes produced splicing isoforms, in which exon skipping was the most frequent AS event. By analyzing the expression differences of mRNAs and splicing isoforms, we found that differentially expressed mRNAs with splicing isoforms could participate in skeletal muscle development and fatty acid metabolism, which might determine muscle-related traits. SERBP1, MYL1, TNNT3, and TNNT1 were identified with multiple splicing isoforms, with significant differences in expression. AS events occurring in IFI6 and GADD45G may cause significant differences in gene expression. Other AS events, such as ONT.15153.3, may regulate the function of ART1 by regulating the expression of different transcripts. Moreover, co-expression and protein-protein interaction (PPI) analysis indicated that several genes (MRPL27, AAR2, PYGM, PSMD4, SCNM1, and HNRNPDL) may be related to intramuscular fat. The splicing isoforms investigated in our research provide a reference for the study of alternative splicing regulation of intramuscular fat deposition.