Key Genes Associated With Non-Alcoholic Fatty Liver Disease and Polycystic Ovary Syndrome
Yong Chen
Leikai Ma
Zhouling Ge3 
Lubin Xie- 1Department of Anesthesiology, The First Affiliated Hospital of Wenzhou Medical University, Wenzhou, China
- 2Department of Obstetrics and Gynecology, The First Affiliated Hospital of Wenzhou Medical University, Wenzhou, China
- 3Department of Respiratory Medicine, The Third Affiliated Hospital of Shanghai University (Wenzhou People’s Hospital), Wenzhou, China
Background: Polycystic ovary syndrome (PCOS) is the most common metabolic and endocrinopathies disorder in women of reproductive age and non-alcoholic fatty liver (NAFLD) is one of the most common liver diseases worldwide. Previous research has indicated potential associations between PCOS and NAFLD, but the underlying pathophysiology is still not clear. The present study aims to identify the differentially expressed genes (DEGs) between PCOS and NAFLD through the bioinformatics method, and explore the associated molecular mechanisms.
Methods: The microarray datasets GSE34526 and GSE63067 were downloaded from Gene Expression Omnibus (GEO) database and analyzed to obtain the DEGs between PCOS and NAFLD with the GEO2R online tool. Next, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis for the DEGs were performed. Then, the protein-protein interaction (PPI) network was constructed and the hub genes were identified using the STRING database and Cytoscape software. Finally, NetworkAnalyst was used to construct the network between the targeted microRNAs (miRNAs) and the hub genes.
Results: A total of 52 genes were identified as DEGs in the above two datasets. GO and KEGG enrichment analysis indicated that DEGs are mostly enriched in immunity and inflammation related pathways. In addition, nine hub genes, including TREM1, S100A9, FPR1, NCF2, FCER1G, CCR1, S100A12, MMP9, and IL1RN were selected from the PPI network by using the cytoHubba and MCODE plug-in. Then, four miRNAs, including miR-20a-5p, miR-129-2-3p, miR-124-3p, and miR-101-3p, were predicted as possibly the key miRNAs through the miRNA-gene network construction.
Conclusion: In summary, we firstly constructed a miRNA-gene regulatory network depicting interactions between the predicted miRNA and the hub genes in NAFLD and PCOS, which provides novel insights into the identification of potential biomarkers and valuable therapeutic leads for PCOS and NAFLD.
Introduction
PCOS is a heterogeneous disorder characterized by hyperandrogenemia, ovulatory dysfunction, and it is a common reproductive and endocrine disorder that mainly occurs in puberty and childbearing age, affecting 6%–20% of women (Rotterdam, 2004; Azziz et al., 2009; Rocha et al., 2019). Furthermore, PCOS is associated with multiple metabolic disorders, including insulin resistance (IR), obesity, dyslipidemia, type 2 diabetes (T2DM), and cardiovascular disease (Goyal and Dawood, 2017). PCOS is the most common cause of anovulatory infertility, and its high incidence has brought huge health and economic costs to the family and society (Homburg and obstetrics, 2004).
NAFLD as one of the most prevalent chronic liver diseases globally is a type of metabolic-related liver disease closely related to IR and genetic susceptibility that includes non-alcoholic simple hepatic steatosis (NAFL), non-alcoholic steatosis Hepatitis (NASH), cirrhosis and hepatocellular carcinoma (HCC) (Diehl and Day, 2017). NAFLD is commonly recognized as a hepatic manifestation of the metabolic syndrome (MS), and it is frequently accompanied by T2DM, obesity, and dyslipidemia (Wattacheril, 2020). According to statistics, the global prevalence of NAFLD was estimated to be 25% (Younossi et al., 2016). Therefore, NAFLD has become a health problem and an economic burden in the world.
In 2005, Brown et al. first found a potential relationship between PCOS and NAFLD (Brown et al., 2005). To date, there are more and more that studies have substantiated a significant link between NAFLD and PCOS (Kelley et al., 2014). Some reports have suggested that 68 the prevalence of NAFLD in women with PCOS is significantly higher than that in 69 non-PCOS women (Falzarano et al., 2022). PCOS patients can promote the occurrence and development of NAFLD through various factors such as hyperandrogenism, IR, obesity, dyslipidemia, chronic low-grade inflammation, and intestinal flora imbalance (Vassilatou, 2014; Boursier et al., 2016; Kim et al., 2017; Yurtda ş and Akdevelioğlu, 2020). Although the relationship between NAFLD and PCOS has attracted substantial interest recently, leading to a vast amount of research related to this topic, the related genetics research is still limited and needs further exploration. There is a need for further studies to investigate NAFLD in the context of PCOS, as well as potential therapeutic options.
In recent years, microarray technology and bioinformatic analysis have been widely used to screen genetic alterations at the genome level. In the present study, we for the first time analyzed two original microarray datasets that were downloaded from Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo) to obtain DEGs between NAFLD and PCOS. Subsequently, the identified DEGs were analyzed by Gene Ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and protein-protein interaction (PPI) analysis. Then, we predicted the potential target miRNAs of DEGs, and constructed a visualized miRNA-gene interaction network through Cytoscape software (http://www.cytoscape.org/). The present study aimed to identify hub genes and hot research topic.
Materials and Methods
Gene Expression Profile Data Collection
Based on the GPL570 platform, two microarray datasets (GSE63067, GSE34526) of PCOS, NAFLD and control samples were collected from the GEO database. The GSE34526 data set consisted of the gene expression profiles of seven PCOS patients and three normal controls. The GSE63067 data set included 18 samples, of which 11 were NAFLD patients, and seven were controls.
Differentially Expressed Genes Selection
DEGs were extracted and analyzed separately using the GEO2R online analysis tool (https://www.ncbi.nlm.nih.gov/geo/geo2r/), which is an R-based web application included in the GEO database. Adjusted p values were used to reduce the false-positive rate using Benjamini and Hochberg method. The threshold of DEGs screening was ∣log2 FC | ≥1 and p < 0.05. The DEGs obtained from the two datasets were visualized using the R packages “complexheatmap” and “ggplot2” to generate the heat maps and volcano maps, respectively. In addition, the overlapping DEGs between PCOS and NAFLD were delineated using the Venn diagrams with the Venn online platform (http://bioinformatics.psb.ugent.be/webtools/Venn/). These overlapping DEGs were retained for subsequent analysis.
Functional Classification and Pathway Enrichment of DEGs
The above overlapping DEGs were submitted to GO function enrichment analysis, which consisted of biological process (BP), cellular component (CC), and molecular function (MF), and the KEGG signaling pathway enrichment analysis using an R package “clusterProfiler”. The enriched GO terms and KEGG pathways with an adjusted p value < 0.05 was selected.
Protein-Protein Interaction Establishment and Identification of hub Genes
To further explore the interaction among the common genes obtained above, the Search Tool for the Retrieval of Interacting Genes (STRING) (http://string-db.org/) was used for PPI network construction. A minimum interaction score of above 0.4 was considered to be significant. Subsequently, Cytoscape software was used to visualize the PPI network. And then, we used the Cytoscape plug-in Minimal Common Oncology Data Elements (MCODE, http://apps.cytoscape.org/apps/mcode) to screen out key protein expression molecules. Then, the Maximal Clique Centrality (MCC) algorithm in the cytoHubba plug-in (http://hub.iis.sinica.edu.tw/cytohubba/) was applied to screen the hub genes with high connectivity in PPI networks.
MiRNAs Associated With hub Genes
NetworkAnalyst tool (version 3.0, https://www.networkanalyst.ca/) was used to construct the miRNA-gene interactions of the hub genes. Finally, these hub genes and miRNA were plotted through Cytoscape software.
Results
Identification of DEGs in PCOS and NAFLD
We downloaded the series GSE34526 dataset about PCOS and the series GSE63067 dataset about NAFLD from the NCBI GEO database. After screening with the threshold of an adjusted p value < 0.05 and ∣log2 FC | >1.0, 3003 DEGs (2176 upregulated and 827 downregulated) were identified in the GSE34526 dataset, 150 DEGs (133 upregulated and 17 downregulated) were identified in GSE63067 using the “limma” package in R software. The volcano plot and heatmap analyses were used to visualize the DEGs of the two data sets were shown in Figure 1 and Figure 2, respectively. In addition, a Venn diagram analysis was performed to evaluate the common DEGs between GSE34526 and GSE63067. As presented in Fig. 1C, 52 overlapping DEGs were identified.
FIGURE 1. The volcano plot and venn diagram of differentially expressed genes (DEGs). Panel (A), volcano plot of DEGs in GSE34526; Panel (B), volcano plot of DEGs in GSE63067; Panel (C), venn diagram of DEGs in GSE34526 and GSE63067 gene chips. There are also two vertical dashed lines in the figure, representing log2 FC at −1 and 1; The horizontal dashed line represents adjusted p value at 0.05. Abbreviations: DEGs, differentially expressed genes; PCOS, Polycystic ovary syndrome; NAFLD, non-alcoholic fatty liver.
FIGURE 2. The heatmaps of differentially expressed genes (DEGs). Panel (A), heatmap of DEGs in GSE34526; Panel (B), heatmap of DEGs in GSE63067. Red represents upregulated DEGs, blue represents downregulated DEGs, and white represents no significant changes. Abbreviations: DEGs, differentially expressed genes; PCOS, Polycystic ovary syndrome; NAFLD, non-alcoholic fatty liver.
GO and KEGG Enrichment Pathway Analysis of Overlapping DEGs
With the clusterProfiler package in R software, we performed GO and KEGG pathway enrichment analyses to better understand the biological functions of the identified DEGs. After screening with the threshold of adjusted p < 0.05, we select the top five significantly enriched GO terms and top five KEGG terms. The results showed that DEGs were enriched in biological processes, including neutrophil degranulation, neutrophil activation involved in immune response, neutrophil activation, neutrophil mediated immunity, and leukocyte migration (Figure 3, Supplementary Table S1). In terms of cell component, DEGs were principally associated with ficolin-1-rich granule, secretory granule membrane, tertiary granule, collagen-containing extracellular matrix, and external side of the plasma membrane (Figure 3, Supplementary Table S1). The analysis of molecular function indicated that DEGs significantly enriched in RAGE receptor binding, cytokine receptor activity, cytokine binding, G protein-coupled peptide receptor activity, and peptide receptor activity (Figure 3, Supplementary Table S1). The top five significant KEGG pathways of DEGs were enriched in Staphylococcus aureus infection, tumor necrosis factor (TNF) signaling pathway, Tuberculosis, Osteoclast differentiation, and Fluid shear stress and atherosclerosis (Figure 3, Supplementary Table S1).
FIGURE 3. GO and KEGG analysis results of overlapping DEGs. Panel (A), GO biological process enrichment results; Panel (B), GO cell component enrichment results. Panel (C), GO molecular function enrichment results. Panel (D), KEGG pathway enrichment results. An adjusted p < 0.05 was identified as significantly changed GOs and KEGG. The x-axis shows the gene ratio of each term. The y-axis represents the annotation items of GO-BP in A, the annotation of GO-CC in B, the annotation of GO-MF in C, and the annotation of KEGG in D. The bubble size represents the number of genes associated with each term. The color of each bubble represents the adjusted p value: the redder the color, the higher the enrichment. Abbreviations: GO, Gene Ontology; BP, biological process; CC, cellular component; MF, molecular function; KEGG, Kyoto Encyclopedia of Genes and Genomes.
PPI Network Construction and hub Genes Identification
The PPI network was first performed based on the STRING database to determine the interactions among the overlapping DEGs. Then, the obtained results were imported into Cytoscape software for visual analysis (Figure 4A). And the interaction number of each gene was also shown (Figure 4B). The PPI network was analyzed by Cytoscape plug-in cytoHubba to identify hub genes. Based on the MCC algorithm, the top 10 genes were identified as potential hub genes: Triggering receptor expressed in myeloid cells 1 (TREM1), S100 calcium binding protein A9 (S100A9), Formyl peptide receptor 1 (FPR1), Neutrophil cytosolic factor 2 (NCF2), Fc fragment Of IgE receptor Ig (FCER1G), C-C chemokine receptor type 1 (CCR1), S100A12 (S100 calcium binding protein A12), Matrix metalloproteinases 9 (MMP9), and IL-1 receptor antagonist (IL1RN) (Figure 5A). Next, the MCODE plug-in was used to identify significant gene clusters modules and obtain cluster scores (filter criteria: degree cut-off = 2; node score cut-off = 0.2; k-core = 2; max depth = 100), and three modules were obtained, which were shown in Figures 5B–D. These modules contained 13 potential hub genes: IL1RN, MMP9, S100A9, TREM1, S100A12, NCF2, FPR1, FCER1G, CCR1, VNN2, SOCS3, IRF1, and IFNGR1. We took the intersection of the hub gene obtained by cytoHubba and the hub gene obtained by MCODE to obtain nine hub genes: TREM1, S100A9, FPR1, NCF2, FCER1G, CCR1, S100A12, MMP9, and IL1RN.
FIGURE 4. Protein-protein interaction network. Panel (A), The PPI among overlapping DEGs. Panel (B), The interaction number of each DEG. Abbreviations: DEGs, differentially expressed genes; PPI, protein-protein interaction.
FIGURE 5. The PPI network analyzed by cytoHubba and MCODE. Panel (A), PPI network for the top 10 hub genes. Panel (B–D), three cluster modules extracted by MCODE. Abbreviations: PPI, protein-protein interaction; MCODE, Minimal Common Oncology Data Elements.
Target miRNAs Prediction and Integrated miRNAs-Targets Network Construction
The NetworkAnalyst databases were used to anticipate target miRNAs of hub genes. We used Cytoscape software to construct the miRNA-gene interaction network, which comprised 67 nodes and 96 edges. As shown in Figure 6, NCF2, S100A9, and FPR1 could be interacted with two common target miRNAs: miR-20a-5p and miR-101-3p. miR-129-2-3p was interacted with four hub genes, including TREM1, NCF2, IL1RN, CCR1, and FPR1. miR-124-3p was interacted with four hub genes, including NCF2, FCER1G, MMP9, and IL1RN. However, these findings need further validation.
FIGURE 6. Integrated miRNA-gene interaction networks for the nine hub genes. Red circles represent nine hub genes. Blue squares represent miRNA which has connectivity with hub genes. Abbreviations: miRNA, microRNA.
Discussion
Nowadays, more and more studies have confirmed the correlation between PCOS and NAFLD. A previous study showed that the prevalence of NAFLD was significantly higher in women with PCOS than in those without, with a prevalence of 51.56 and 29.64%, respectively (Falzarano et al., 2022). Meanwhile, a study revealed that PCOS was also more prevalent in female NAFLD patients with a prevalence of 50% up to a high of 70%, which is significantly higher than in those without NAFLD (Georgescu, 2022). Furthermore, studies have shown that both PCOS patients and NAFLD patients can increase the risk of cardiovascular disease (Dokras, 2019), and metabolic abnormalities were more pronounced in the PCOS patients with concurring NAFLD (Sarkar et al., 2020). In addition, cardiovascular disease is the most common cause of death in NAFLD patients (Söderberg et al., 2010). So far, the mechanism of the link between NAFLD and PCOS is not fully understood. Therefore, exploring the molecular mechanisms between PCOS and NAFLD, early identification and intervention, undoubtedly have important clinical significance.
In the present study, a series of bioinformatics analyses were performed on two independent gene chip databases of PCOS and NAFLD, and 52 common DEGs between PCOS and NAFLD were obtained based on the GEO database. The result of GO enrichment analysis indicated that the DEGs were mainly enriched in immune response, neutrophil activation, neutrophil mediated immunity, and leukocyte migration. This result suggests that the DEGs between PCOS and PCOS are related to inflammation and immune response, which is consistent with previous studies (Moulana, 2019; Oates et al., 2019). In addition, KEGG analysis showed that the DEGs were mainly enriched in Staphylococcus aureus infection, TNF signaling pathway, Tuberculosis, Osteoclast differentiation, and Fluid shear stress and atherosclerosis. We further classified these pathways according to the KEGG pathway Database. We found that they are mainly related to immunity and inflammation. Many studies have shown that immune inflammatory processes and atherosclerosis play a critical role in cardiovascular diseases (Zhang et al., 2020). Moreover, PCOS and NAFLD are strongly associated with cardiovascular diseases (Scicchitano et al., 2012; Niederseer et al., 2021). Our findings provide new ideas that patients with PCOS and NAFLD are more prone to cardiovascular disease.
Subsequently, according to the MCODE plug-in and cytoHubba plug-in of Cytoscape, we screened nine overlapping DEGs as hub genes in the PPI network, including TREM1, S100A9, FPR1, NCF2, FCER1G, CCR1, S100A12, MMP9, and IL1RN. These nine genes were all upregulated in both PCOS patients and NAFLD patients, suggesting these genes may play important role in the mechanism of PCOS and NAFLD.
S100A9 and S100A12 belong to a family of 25 homologous low‐molecular‐weight intracellular calcium‐binding proteins produced by cells of myeloid origin (Donato and biology, 2001). S100A9 is mainly released by activated granulocytes and associated with various types of inflammation-related pathways (Wang et al., 2018). A previous study showed that S100A9 could increase the production of inflammatory cytokines and disturb steroidogenesis via activating nuclear factor kappa B (NF‐κB) signaling pathway in PCOS (Li et al., 2020). Furthermore, S100A9 promotes inflammation and lipolysis in the liver during NAFLD progression (Wu et al., 2020). Another study demonstrated that S100A9 can differentiate hepatic and metabolic progression in NAFLD as a biomarker in NAFLD progression (Liu et al., 2015). Additionally, studies have shown that S100A12 interacts with S100A9 (Hatakeyama et al., 2004) and also has proinflammatory activity (Miranda et al., 2001). There are no reports of S100A12 in PCOS and NAFLD. We speculate that S100A9 and S100A12 synergistically promote inflammation in PCOS and NAFLD. Interestingly, S100A12 has also been shown to block the activity of S100A9 (Katte and Yu, 2018).
NCF2 is a rate-limiting cofactor of NADPH oxidase that is necessary for reactive oxygen species (ROS) production in phagocytes, which plays a critical role in innate immunity and phagocytic microbicidal activity (Muise et al., 2012). NCF2 is highly expressed in PCOS and NAFLD and contributes to oxidative stress in PCOS and NAFLD (Kaur et al., 2012; Garcia-Jaramillo et al., 2019).
IL1RN is originally discovered as a natural antagonist of IL-1 (Arend, 1991), which encodes the IL1 antagonist protein (IL1RA) (Gabay et al., 2010). And IL1RA has been reported as a critical mediator of inflammatory processes and plays a vital role in the pathogenesis of PCOS (Xia et al., 2013). IL1RN also plays an important role in the development of NAFLD (Wolfs et al., 2015).
MMP9 is one of the zinc-ion-dependent metalloproteinase family, playing a pathogenic role in chronic inflammatory autoimmune diseases (Akter et al., 2015). MMP9 was shown to be overexpressed in PCOS and involved in the pathogenesis of PCOS (Ranjbaran et al., 2016). MMP9 could increase fibrosis severity in NAFLD and be involved in the development and progression of human HCC metastasis (Shi et al., 2010; Goyale et al., 2021). TREM1 is expressed on myeloid lineage cells and played an important role in maintaining homeostasis and the normal immune responses. Dou et al. demonstrated that TREM1 in the liver was involved in the pathological process and therapeutic outcome of NASH (Dou et al., 2016). In addition, TREM1 promotes atherosclerosis by inducing monocytosis and foam cell formation (Zysset et al., 2016). FPR1 is a member of 7-transmembrane G protein-coupled receptors (GPCRs), which is mostly expressed by myeloid cells and is a key regulator of the inflammatory environment (Dahlgren et al., 2016; Mastromarino et al., 2019). FCER1G is engaged in a variety of immune responses as a constitutive component of the high-affinity immunoglobulin E (Ig E) receptor and interleukin-3 receptor complex (Liang et al., 2013). CCR1 is a member of the seven-transmembrane G-protein-coupled receptor family, which is widely expressed on myeloid cells and is involved in the activation and trafficking of immune cells (Schaller et al., 2008). Based on previous research, these nine hub genes have been shown to have key effects on immune response and inflammatory response. Moreover, it has been known that immune dysregulation and chronic inflammation are important factors in the development of PCOS and NAFLD (Anstee et al., 2013; Qin et al., 2016). Based on this, we speculate that: these genes may also contribute to the occurrence and development of PCOS and NAFLD, and may be potentially new treatment targets in the future.
MiRNAs are small non-coding RNA molecules (21–25 nucleotides long) that are partially or fully complementary to the 3′ URT of the target gene mRNA by inducing mRNA degradation or repressing mRNA translation (He and Hannon, 2004). In this study, we also constructed a miRNA-target gene network and selected the four miRNAs (miR-20a-5p, miR-129-2-3p, miR-124-3p, and miR-101-3p) that interacted the most with DEGs. miR-20a-5p was found to bind to pro-angiogenic genes and suppress their expression and thus exerting anti-angiogenic activity in PCOS (Patil et al., 2022). miR-20a-5p also played a protective role in lipid metabolism and IR of NAFLD (Wang et al., 2020; Zhang et al., 2021). Previous studies on miR‐129-2-3p mainly focused on cancer. miR-129-2-3p has not been reported in PCOS and NAFLD. A study has shown that miR-129-2-3p was involved in the functional regulation of neutrophils in chronic inflammatory processes. Thus, we speculate that miR-129-2-3 may be involved in the pathophysiology of PCOS through inflammatory processes. miR-124-3p was downregulated in PCOS patients (Ding et al., 2015) and was involved in lipid metabolism in NAFLD as a critical regulator in lipid homeostasis of the liver (Wang et al., 2020918). miR-101-3p was also downregulated in PCOS patients (Butler et al., 2020). Moreover, the downregulation of miR-101-3p was more strongly associated with hepatic injury in NAFLD (Meroni et al., 2019).
In conclusion, our study identified some hub genes and speculated that miRNA-gene regulatory networks contribute to the pathophysiological process of PCOS and NAFLD by bioinformatics analysis, which provided the potential diagnostic and therapeutic targets of PCOS and NAFLD. However, there are some limitations to our study. Firstly, only the top nine hub genes were involved in our current study. Secondly, a lack of research on detailed molecular mechanisms that the hub genes and miRNAs regulate in PCOS and NAFLD. Thirdly, the miRNA-gene interaction networks were only based on predictions from public databases. Thus, the molecular mechanism of these hub genes and miRNAs in the occurrence of PCOS and NAFLD need to be further studied.
Conclusion
Our study identified a miRNA-gene network potentially relevant for PCOS and NAFLD. The nine hub genes (including TREM1, S100A9, FPR1, NCF2, FCER1G, CCR1, S100A12, MMP9, and IL1RN) were significantly upregulated, which may have a critical influence on the pathophysiological mechanism of PCOS and NAFLD. Some potential target miRNAs (miR-20a-5p, miR-129-2-3p, miR-124-3p, and miR-101-3p) were also predicted and may participate in the pathophysiological process of PCOS and NAFLD through inflammation and immune response. These findings may contribute to the development of early diagnostic strategies, prognostic markers, and therapeutic targets.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
LX, YC, and ZG conceptualized the study design and analyze the data. YC and LM wrote the manuscript. LM and YP collected the data.
Funding
Funds were from Zhejiang Provincial natural scientific (LQ21H010002), Wenzhou Science and Technology Bureau (2021Y1398) and the Scientific Research Incubation Project of the First Affiliated Hospital of Wenzhou Medical University (Grant Number: FHY2019070).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2022.888194/full#supplementary-material
References
Akter, H., Park, M., Kwon, O.-S., Song, E. J., Park, W.-S., and Kang, M.-J. (2015). Activation of Matrix Metalloproteinase-9 (MMP-9) by Neurotensin Promotes Cell Invasion and Migration through ERK Pathway in Gastric Cancer. Tumor Biol. 36 (8), 6053–6062. doi:10.1007/s13277-015-3282-9
Anstee, Q. M., Targher, G., and Day, C. P. (2013). Progression of NAFLD to Diabetes Mellitus, Cardiovascular Disease or Cirrhosis. Nat. Rev. Gastroenterol. Hepatol. 10 (6), 330–344. doi:10.1038/nrgastro.2013.41
Arend, W. P. (1991). Interleukin 1 Receptor Antagonist. A New Member of the Interleukin 1 Family. J. Clin. Invest. 88 (5), 1445–1451. doi:10.1172/jci115453
Azziz, R., Carmina, E., Dewailly, D., Diamanti-Kandarakis, E., Escobar-Morreale, H. F., Futterweit, W., et al. (2009). The Androgen Excess and PCOS Society Criteria for the Polycystic Ovary Syndrome: the Complete Task Force Report. Fertil. Steril. 91 (2), 456–488. doi:10.1016/j.fertnstert.2008.06.035
Boursier, J., Mueller, O., Barret, M., Machado, M., Fizanne, L., Araujo-Perez, F., et al. (2016). The Severity of Nonalcoholic Fatty Liver Disease Is Associated with Gut Dysbiosis and Shift in the Metabolic Function of the Gut Microbiota. Hepatology 63 (3), 764–775. doi:10.1002/hep.28356
Brown, A. J., Tendler, D. A., McMurray, R. G., and Setji, T. L. (2005). Polycystic Ovary Syndrome and Severe Nonalcoholic Steatohepatitis: Beneficial Effect of Modest Weight Loss and Exercise on Liver Biopsy Findings. Endocr. Pract. 11 (5), 319–324. doi:10.4158/ep.11.5.319
Butler, A. E., Ramachandran, V., Sathyapalan, T., David, R., Gooderham, N. J., Benurwar, M., et al. (2020). microRNA Expression in Women with and without Polycystic Ovarian Syndrome Matched for Body Mass Index. Front. Endocrinol. 11, 206. doi:10.3389/fendo.2020.00206
Dahlgren, C., Gabl, M., Holdfeldt, A., Winther, M., and Forsman, H. (2016). Basic Characteristics of the Neutrophil Receptors that Recognize Formylated Peptides, a Danger-Associated Molecular Pattern Generated by Bacteria and Mitochondria. Biochem. Pharmacol. 114, 22–39. doi:10.1016/j.bcp.2016.04.014
Diehl, A. M., and Day, C. (2017). Cause, Pathogenesis, and Treatment of Nonalcoholic Steatohepatitis. N. Engl. J. Med. 377 (21), 2063–2072. doi:10.1056/nejmra1503519
Ding, C.-F., Chen, W.-Q., Zhu, Y.-T., Bo, Y.-L., Hu, H.-M., and Zheng, R.-H. (2015). Circulating microRNAs in Patients with Polycystic Ovary Syndrome. Hum. Fertil. 18 (1), 22–29. doi:10.3109/14647273.2014.956811
Dokras, A (2019). Does Body Weight Affect Cardiometabolic Risk in Women with Polycystic Ovary Syndrome? Fertil. Steril. 111 (1), 56–57. doi:10.1016/j.fertnstert.2018.11.001
Donato, R., and biology, c. (2001). S100: a Multigenic Family of Calcium-Modulated Proteins of the EF-Hand Type with Intracellular and Extracellular Functional Roles. Int. J. Biochem. Cell Biol. 33 (7), 637–668. doi:10.1016/s1357-2725(01)00046-2
Dou, Z. F., Guo, Y. R., Liu, J. C., Li, N., Chai, B., Li, X. Q., et al. (2016). Ameliorative Effects of glycine in an Experimental Nonalcoholic Steatohepatitis and its Correlation between TREM-1 and TREM-2. Am. J. Transl. Res. 8 (2), 284–297.
Falzarano, C., Lofton, T., Osei-Ntansah, A., Oliver, T., Southward, T., Stewart, S., et al. Nonalcoholic Fatty Liver Disease in Women and Girls with Polycystic Ovary Syndrome. J Clin Endocrinol Metab. 2022;107(1):258–272.doi:10.1210/clinem/dgab658
Gabay, C., Lamacchia, C., and Palmer, G. (2010). IL-1 Pathways in Inflammation and Human Diseases. Nat. Rev. Rheumatol. 6 (4), 232–241. doi:10.1038/nrrheum.2010.4
Garcia-Jaramillo, M., Spooner, M. H., Löhr, C. V., Wong, C. P., Zhang, W., and Jump, D. B. (2019). Lipidomic and Transcriptomic Analysis of Western Diet-Induced Nonalcoholic Steatohepatitis (NASH) in Female Ldlr -/- Mice. PLoS ONE 14 (4), e0214387. doi:10.1371/journal.pone.0214387
Georgescu, C. E. (2022). Polycystic Ovary Syndrome and Nonalcoholic Fatty Liver Disease. Elsevier, 187–216. doi:10.1016/b978-0-12-823045-9.00007-9
Goyal, M., and Dawood, A. S. (2017). Debates Regarding Lean Patients with Polycystic Ovary Syndrome: A Narrative Review. J. Hum. Reprod. Sci. 10 (3), 154–161. doi:10.4103/jhrs.JHRS_77_17
Goyale, A., Jain, A., Smith, C., Papatheodoridi, M., Misas, M. G., Roccarina, D., et al. (2021). Assessment of Non-alcoholic Fatty Liver Disease (NAFLD) Severity with Novel Serum-Based Markers: A Pilot Study. PLoS ONE 16 (11), e0260313. doi:10.1371/journal.pone.0260313
Hatakeyama, T., Okada, M., Shimamoto, S., Kubota, Y., and Kobayashi, R. Identification of Intracellular Target Proteins of the Calcium-Signaling Protein S100A12. Eur J Biochem. 2004;271(18):3765–3775.doi:10.1111/j.1432-1033.2004.04318.x
He, L., and Hannon, G. J. (2004). MicroRNAs: Small RNAs with a Big Role in Gene Regulation. Nat. Rev. Genet. 5 (7), 522–531. doi:10.1038/nrg1379
Homburg, R., and obstetrics, r. C. (2004). Management of Infertility and Prevention of Ovarian Hyperstimulation in Women with Polycystic Ovary Syndrome. Best. Pract. Res. Clin. Obstet. Gynaecol. 18 (5), 773–788. doi:10.1016/j.bpobgyn.2004.05.006
Katte, R., and Yu, C. (2018). Blocking the Interaction between S100A9 Protein and RAGE V Domain Using S100A12 Protein. PLoS ONE 13 (6), e0198767. doi:10.1371/journal.pone.0198767
Kaur, S., Archer, K. J., Devi, M. G., Kriplani, A., Strauss, J. F., Singh, R., et al. (2012). Differential Gene Expression in Granulosa Cells from Polycystic Ovary Syndrome Patients with and without Insulin Resistance: Identification of Susceptibility Gene Sets through Network Analysis. J. Clin. Endocrinol. Metabolism 97 (10), E2016–E2021. doi:10.1210/jc.2011-3441
Kelley, C. E., Brown, A. J., Diehl, A. M., Setji, T. L. J, and Wjog, W. (2014). Review of Nonalcoholic Fatty Liver Disease in Women with Polycystic Ovary Syndrome. Wjg 20 (39), 14172. doi:10.3748/wjg.v20.i39.14172
Kim, J. J., Kim, D., Yim, J. Y., Kang, J. H., Han, K. H., Kim, S. M., et al. (2017). Polycystic Ovary Syndrome with Hyperandrogenism as a Risk Factor for Non-obese Non-alcoholic Fatty Liver Disease. Aliment. Pharmacol. Ther. 45 (11), 1403–1412. doi:10.1111/apt.14058
Li, H., Huang, X., Chang, X., Yao, J., He, Q., Shen, Z., et al. (2020). S100‐A9 Protein in Exosomes Derived from Follicular Fluid Promotes Inflammation via Activation of NF‐κB Pathway in Polycystic Ovary Syndrome. J. Cell. Mol. Medi 24 (1), 114–125. doi:10.1111/jcmm.14642
Liang, Y., Zhao, M., Liang, G., Yin, H., and Lu, Q. (2013). Construction of Special Reporter to Detect DNA Methylation Regulatory Activity in FCER1G Gene Promoter through Patch-Methylation. Zhong Nan Da Xue Xue Bao Yi Xue Ban. 38 (2), 120–124. doi:10.3969/j.issn.1672-7347.2013.02.002
Liu, X., Wang, Y., Ming, Y., Song, Y., Zhang, J., Chen, X., et al. (2015). S100A9: a Potential Biomarker for the Progression of Non-alcoholic Fatty Liver Disease and the Diagnosis of Non-alcoholic Steatohepatitis. PLoS ONE 10 (5), e0127352. doi:10.1371/journal.pone.0127352
Mastromarino, M., Kirpotina, L. N., Schepetkin, I. A., Quinn, M. T., Lacivita, E., and Leopoldo, M. (2019). Aurantiamide-related Dipeptide Derivatives Are Formyl Peptide Receptor 1 Antagonists. Med. Chem. Commun. 10 (12), 2078–2088. doi:10.1039/c9md00336c
Meroni, M., Longo, M., Erconi, V., Valenti, L., Gatti, S., Fracanzani, A. L., et al. (2019). mir-101-3p Downregulation Promotes Fibrogenesis by Facilitating Hepatic Stellate Cell Transdifferentiation during Insulin Resistance. Nutrients 11 (11), 2597. doi:10.3390/nu11112597
Miranda, J., Tao, T., Jones, A., Chernushevich, I., Standing, K. G., Geczy, C. L., et al. 2001. Imaginaire et utopie dans la politique du Nordeste brésilien. 488(1-2):85–105.doi:10.4000/books.pul.11481
Moulana, M. (2019). Immunophenotypic Profile of Leukocytes in Hyperandrogenemic Female Rat an Animal Model of Polycystic Ovary Syndrome. Life Sci. 220, 44–49. doi:10.1016/j.lfs.2019.01.048
Muise, A. M., Xu, W., Guo, C.-H., Walters, T. D., Wolters, V. M., Fattouh, R., et al. (2012). NADPH Oxidase Complex and IBD Candidate Gene Studies: Identification of a Rare Variant inNCF2that Results in Reduced Binding to RAC2. Gut 61 (7), 1028–1035. doi:10.1136/gutjnl-2011-300078
Niederseer, D., Wernly, B., Aigner, E., Stickel, F., and Datz, C. NAFLD and Cardiovascular Diseases: Epidemiological, Mechanistic and Therapeutic Considerations. J Clin Med. 2021;10(3). 467. doi:10.3390/jcm10030467
Oates, J. R., McKell, M. C., Moreno-Fernandez, M. E., Damen, M. S. M. A., Deepe, G. S., Qualls, J. E., et al. (2019). Macrophage Function in the Pathogenesis of Non-alcoholic Fatty Liver Disease: The Mac Attack. Front. Immunol. 10, 2893. doi:10.3389/fimmu.2019.02893
Patil, K., Joseph, S., Shah, J., and Mukherjee, S. (2022). Genetics. An Integrated in Silico Analysis Highlighted Angiogenesis Regulating miRNA-mRNA Network in PCOS Pathophysiology, Springer, 1–14.
Qin, L., Xu, W., Li, X., Meng, W., Hu, L., Luo, Z., et al. (2016). Differential Expression Profile of Immunological Cytokines in Local Ovary in Patients with Polycystic Ovarian Syndrome: Analysis by Flow Cytometry. Eur. J. Obstetrics Gynecol. Reproductive Biol. 197, 136–141. doi:10.1016/j.ejogrb.2015.12.003
Ranjbaran, J., Farimani, M., Tavilani, H., Ghorbani, M., Karimi, J., Poormonsefi, F., et al. Matrix Metalloproteinases 2 and 9 and MMP9/NGAL Complex Activity in Women with PCOS. Reproduction. 2016;151(4):305–311.doi:10.1530/rep-15-0340
Rocha, A., Oliveira, F., Azevedo, R., Silva, V., Peres, T., Candido, A., et al. (2019). Recent Advances in the Understanding and Management of Polycystic Ovary Syndrome, F1000Res. 8. doi:10.12688/f1000research.15318.1
Rotterdam, E. A-S. P. C. W. G. (2004). Revised 2003 Consensus on Diagnostic Criteria and Long-Term Health Risks Related to Polycystic Ovary Syndrome. Fertil. Steril. 81 (1), 19–25. doi:10.1016/j.fertnstert.2003.10.004
Sarkar, M., Terrault, N., Chan, W., Cedars, M. I., Huddleston, H. G., Duwaerts, C. C., et al. (2020). Polycystic Ovary Syndrome (PCOS) Is Associated with NASH Severity and Advanced Fibrosis. Liver Int. 40 (2), 355–359. doi:10.1111/liv.14279
Schaller, M. A., Kallal, L. E., and Lukacs, N. W. (2008). A Key Role for CC Chemokine Receptor 1 in T-Cell-Mediated Respiratory Inflammation. Am. J. Pathology 172 (2), 386–394. doi:10.2353/ajpath.2008.070537
Scicchitano, P., Dentamaro, I., Carbonara, R., Bulzis, G., Dachille, A., Caputo, P., et al. (2012). Cardiovascular Risk in Women with PCOS. Int. J. Endocrinol. Metab. 10 (4), 611–618. doi:10.5812/ijem.4020
Shi, G.-M., Ke, A.-W., Zhou, J., Wang, X.-Y., Xu, Y., Ding, Z.-B., et al. (2010). CD151 Modulates Expression of Matrix Metalloproteinase 9 and Promotes Neoangiogenesis and Progression of Hepatocellular Carcinoma. Hepatology 52 (1), 183–196. doi:10.1002/hep.23661
Söderberg, C., Stål, P., Askling, J., Glaumann, H., Lindberg, G., Marmur, J., et al. (2010). Decreased Survival of Subjects with Elevated Liver Function Tests during a 28-year Follow-Up. Hepatology 51 (2), 595–602. doi:10.1002/hep.23314
Vassilatou, E. (2014). Nonalcoholic Fatty Liver Disease and Polycystic Ovary Syndrome. World J. Gastroenterol. 20 (26), 8351–8363. doi:10.3748/wjg.v20.i26.8351
Wang, G., Zou, H., Lai, C., Huang, X., Yao, Y., and Xiang, G. J. Fi. E. (2020). Repression of MicroRNA-124-3p Alleviates High-Fat Diet–Induced Hepatosteatosis by Targeting Pref-1. Frontiers in Endocrinology. 11. 589994. doi:10.3389/fendo.2020.589994
Wang, S., Song, R., Wang, Z., Jing, Z., Wang, S., and Ma, J. (2018). S100A8/A9 in Inflammation,Front Immunol. 9.1298. doi:10.3389/fimmu.2018.01298
Wang, X., Ma, Y., Yang, L-Y., Zhao, D. J, Fi, C., and Biology, D. (2020) MicroRNA-20a-5p Ameliorates Non-alcoholic Fatty Liver Disease via Inhibiting the Expression of CD36. Frontiers in Cell and Developmental Biology. 8 1413. doi:10.3389/fcell.2020.596329
Wattacheril, J. (2020). Extrahepatic Manifestations of Nonalcoholic Fatty Liver Disease. Gastroenterology Clin. N. Am. 49 (1), 141–149. doi:10.1016/j.gtc.2019.10.002
Wolfs, M. G. M., Gruben, N., Rensen, S. S., Verdam, F. J., Greve, J. W., Driessen, A., et al. (2015). Determining the Association between Adipokine Expression in Multiple Tissues and Phenotypic Features of Non-alcoholic Fatty Liver Disease in Obesity. Nutr Diabetes 5 (2), e146. doi:10.1038/nutd.2014.43
Wu, J., Zheng, L., Mo, J., Yao, X., Fan, C., and Bao, Y. (2020). Protective Effects of MitoTEMPO on Nonalcoholic Fatty Liver Disease via Regulating Myeloid-Derived Suppressor Cells and Inflammation in Mice. BioMed Res. Int. 2020. 1-9. doi:10.1155/2020/9329427
Xia, Y.-H., Yao, L., and Zhang, Z.-X. (2013). Correlation between IL-1β, IL-1Ra Gene Polymorphism and Occurrence of Polycystic Ovary Syndrome Infertility. Asian Pac. J. Trop. Med. 6 (3), 232–236. doi:10.1016/s1995-7645(13)60030-9
Younossi, Z. M., Koenig, A. B., Abdelatif, D., Fazel, Y., Henry, L., and Wymer, M. (2016). Global Epidemiology of Nonalcoholic Fatty Liver Disease-Meta-Analytic Assessment of Prevalence, Incidence, and Outcomes. Hepatology 64 (1), 73–84. doi:10.1002/hep.28431
Yurtdaş, G., and Akdevelioğlu, Y. A. (2020). A New Approach to Polycystic Ovary Syndrome: the Gut Microbiota. Journal of the American College of Nutrition. 39(4). 371–382. doi:10.1080/07315724.2019.1657515
Zhang, H., Niu, Q., Liang, K., Li, X., Jiang, J., Bian, C., et al. (2021). Effect of LncPVT1/miR-20a-5p on Lipid Metabolism and Insulin Resistance in NAFLD. Dmso 14, 4599–4608. doi:10.2147/dmso.s338097
Zhang, S., Cheng, M., Wang, Z., Liu, Y., Ren, Y., Rong, S., et al. (2020). Secoisolariciresinol Diglucoside Exerts Anti-inflammatory and Antiapoptotic Effects through Inhibiting the Akt/IκB/NF-Κb Pathway on Human Umbilical Vein Endothelial Cells, Mediators of Inflammation. 2020. 1-8.doi:10.1155/2020/3621261
Keywords: PCOS (polycystic ovarian syndrome, NAFLD (non alcoholic fatty liver disease), bioinformatics analysis, differentially expressed genes (DEG’s), MiRNA-mRNA regulatory network
Citation: Chen Y, Ma L, Ge Z, Pan Y and Xie L (2022) Key Genes Associated With Non-Alcoholic Fatty Liver Disease and Polycystic Ovary Syndrome. Front. Mol. Biosci. 9:888194. doi: 10.3389/fmolb.2022.888194
Received: 02 March 2022; Accepted: 25 April 2022;
Published: 25 May 2022.
Edited by:
Mario Bermudez De Leon, Mexican Social Security Institute (IMSS), MexicoReviewed by:
Stanley Andrisse, Howard University, United StatesIrene Martin-Estal, Escuela de Medicina y Ciencias de la Salud Tec Salud, Mexico
Copyright © 2022 Chen, Ma, Ge, Pan and Xie. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lubin Xie, xielubin1993@163.com
†These authors have contributed equally to this work and share first authorship