DOI: https://doi.org/10.21203/rs.3.rs-2241031/v1
Objective:Ferroptosis-related genes may have a critical regulatory role in the pathogenetic process of Crohn’ disease(CD).The purpose of this study was to identify genes expressed in CD that are associated with ferroptosis and provide direction in the diagnosis and therapy of Crohn's disease.
Methods: The data for CD mRNA expression were first gathered from the Gene Expression Omnibus (GEO) database, and two gene sets were selected as major targets (GSE75214 and GSE102133) and analyzed differentially expressed genes.Next,R software (version 4.1.2) was used to analyze the common genes in CD differential expressed genes and ferroptosis-related genes.GO enrichment analysis,KEGG pathway analysis were used to identify differential related pathways and functions.Protein-protein interaction(PPI) analysis was performed to identify target genes.DSigDB website was used to predict potential target drugs for hub genes.Finally, qRT-PCR method were used to detect the expression of these ferroptosis related genes in clinical samples obtained from healthy control and CD patients.
Result: According to two GEO datasets, we finally identified 13 ferroptosis DEGs(10 upregulated genes and 2 downregulated genes) in crohn disease with the threshold of p-value < 0.05 and |log2 FC| > 1 and selected for continued analysis.Go enrichment analysis and KEGG pathways results were shown in the following figures.PPI analysis indicate the mutual effect between these genes and filtered out 5 hub genes.Top 10 potential targeted drugs were selected. Finally, the result of qRT-PCR shown that the expression of three genes IL-6,PTGS2 and DUOX2 were different between CD samples and healthy samples.This result was consistent with the results we obtained in the biological information analysis.
Conclusion: Bioinformatics analysis identified a total of 13 iron death-associated genes in CD. Three differential genes IL-6,PTG32 and DUOX2 were detected in tissue experiments.Our findings might provide new biomarkers and promising treatment targets in CD.
Inflammatory bowel diseases(IBD) are currently known chronic inflammatory diseases of the gastrointestinal tract, including ulcerative colitis(UC) and Crohn's disease(CD)1 .The incidence of CD has been increasing in recent years around the global scope23.CD is usually found at the end of the ileum and can invade and spread to the surrounding intestinal tissues4 it also has a tendency to develop dysplasia and colorectal cancer(CRC)5.Patients with long-term CD have about 22% chance of developing cancer6. Ferroptosis is a newly discovered form of programmed cell death that is associated with many benign and malignant diseases of the digestive system7.However, the mechanisms involved in CD caused by ferroptosis-related genes are not fully understood.
Ferroptosis is a form of cell death that differs from apoptosis and pyroptosis, The usual causes are lipid peroxidation, iron accumulation, and cell membrane breakdown.8.Some other research had found that ferroptosis can be participant in the regulation of inflammatory bowel disease9. Huang et.al indicated that ferroptosis-related hub gene STAT3(signal transducer and activator of transcription 3) mediated the ferroptosis process and promoted ulcerative colitis occured10. Ferroptosis also has a significant part in the regulation of CRC.Glutathione peroxidase 4 (GPX4),a ferroptosis inhibitor,was an important central regulatory molecule that was involved in the progress of the CRC development11.
Although ferroptosis is strongly associated with these digestive tract-related diseases, the mechanisms involved in ferroptosis in CD were not understood.So we search for the ferroptosis-related genes in CD to find out the correlation between CD and ferroptosis.GSE102133 and GSE75214 were found as target datasets to find different express genes in CD and such genes made cross-comparison with ferroptosis-related genes to finally determine the object of study. Continuous The result will provide new biomarkers for further study in CD research,to validate the potential effect of clinical diagnosis and treatment.
We searched the CD-related microarray sequencing datasets on the GEO database (https://www.ncbi.nlm.nih.gov/geo/) to analysis in this study,the keyword is “crohn disease”.There are the following screening criteria:(a)the samples from human tissues.(b)the dataset was from Microarray expression data.(c)the dataset was required from CD ileal mucosa and healthy ileal mucosa.(d)the totally number of samples was more than 10.(e)the number of different expression genes(DEGs) was greater than 100.By this way, we finally determine two GEO datasets GSE75214 and GSE102133 to enroll in our study.Both of these datasets were present on GPL6244 (Affymetrix Human Gene 1.0 ST Array),GSE102133 included 65 CD patients and 12 controls.GSE75214 contained 67 CD patients and 11 controls,two GSE datasets specific informations was shown in Table 1.Ferroptosis-related genes get from FerrDb website (http://www.zhounan.org/ferrdb/index.html), there are 259 genes were found.
GSE dataset Name | GSE75214 | GSE102133 |
---|---|---|
Platform | GPL6244 | GPL6244 |
Organism | Homo sapiens | Homo sapiens |
Contributor(s) | Arijs I et.al | Verstockt S et.al |
Organization name | TARGID - IBD Leuven,Clinical and Experimental Medicine of KU Leuven | KU Leuven |
Sample site | Terminal ileum | Ileal mucosal biopsies |
Sample Number | 78/194 | 77/77 |
Last update date | Jul 26, 2018 | Jul 25, 2021 |
All microarray data were from the GEO database,the relative extracted data were normalized by log-2 transformation.Probes were converted to gene symbols according to the annotation information of the normalized data in the platform.The principal component analysis (PCA) plot was performed between two datasets.Data was standardized using the “Limma” toolkit in the R software (version 4.1.2). A | log2 (fold change)| > 1 and a P-value < 0.05 were used as the criteria for difffferential gene expression. The heatmap and volcano plot were drawn using the “Heatmap” and “ggplot2” packages in the R software (version 4.1.2). Then the DEGs made intersected with ferroptosis-related genes from FerrDb and generated a venn diagram by using online Bioinformatics tool (http://bioinformatics.psb.ugent.be/webtools/Venn/).
To further investigate the biological function of the genes involved,We used the Gene Ontology (GO) enrichment analyses which included biological process (BP)、molecular function (MF) and cellular component (CC) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses for evaluating relating functions.R package clusterProfifiler toolkit in R software version(4.1.2) was used for such calculation.A p-value < 0.05 result was defined as significant. Gene Correlation Analysis Using Pearson correlation test to Assess Reciprocal Relationships at the mRNA Level in CD.
To further investigate the interactions between DEGs promoting proteins, PPI network analysis was performed by using the Search Tool for the Retrieval of Interacting Genes (STRING) database.The PPI network analysis map was visualized by Cytoscape software(version 3.9.1).Firstly, STRING database was used to establish the PPI network according to DEGs,after that the analysis results were then exported to Cytoscape software(version 3.9.1) for mapping.Nodes represent proteins, lines represent interactions between proteins, and different colors represent different functions.Hub gene was obtained from the cytoHubba plugin in Cytoscape (version 3.9.1) through maximal clique centrality (MCC) methods. The GeneMANIA website (https://genemania.org/) was used for co-expression, physical interaction, co-localization, pathways and other analyses of the 5 hub gene.
DSigDBv1.0 (http://dsigdb.tanlab.org/DSigDBv1.0/) was used to predict the five hub genes of potential drug targets, and select the top 10 scores as candidate drugs.
From October 2021 to September 2022,we totally collected 3 pairs of tissues(3 CDs and 3 controls) from patients and healthy individuals at The Qingdao Municipal Hospital(Qingdao, China). Written informed consent was obtained from all patients and approved by the medical ethics committee of the hospital (2022-066).
The R software (version 4.1.2) is used to perform statistical analysis of biological information data. Results were analyzed by Student t test for comparison. In the qRT-PCR experiment, each sample was repeated three times.A p-value < 0.05 was considered statistically significant.
As shown in the flow chart (Fig. 1).We chose two CD dataset files (GSE102133 and GSE75214,Table 1) and downloaded them using the same platform files. A total of 155 samples were included in the study (132 CD samples and 23 control samples).Principal component analysis (PCA) was performed between two datasets(Supplementary Fig. 1A).The differential genes between the two groups of samples are shown in the volcano (Fig. 2A). 13 differential genes in total(11 upregulated genes and 2 downregulated genes)were obtained after crossover with 259 ferroptosis-related genes showing in Venn diagram (Fig. 2) and made heatmap plot between CD samples and normal samples.Expression of the 13 differential genes in CD's samples and healthy samples in box plots(Fig. 3).The specific information with the differential genes is shown in Table2.
Gene symbol | Description | logFC | UP/DOWN | P-value |
---|---|---|---|---|
DUOX2 | Dual oxidase 2 | 4.643206134 | Up | 4.23153E-25 |
MUC1 | Mucin 1, cell surface associated | 2.986183308 | Up | 1.65188E-23 |
NOS2 | Nitric oxide synthase 2 | 2.605279822 | Up | 5.37157E-20 |
HIF1A | Hypoxia inducible factor 1 alpha | 1.055152853 | Up | 8.37052E-15 |
ACSL4 | Acyl-CoA synthetase long chain family member 4 | 1.553297113 | Up | 6.5226E-13 |
SLC7A11 | Solute carrier family 7 member 11 | 1.72699847 | Up | 1.37298E-11 |
CXCL2 | C-X-C motif chemokine ligand 2 | 1.109787323 | Up | 1.3473E-10 |
PTGS2 | Prostaglandin-endoperoxide synthase 2 | 2.034361857 | Up | 1.76391E-08 |
NCF2 | Neutrophil cytosolic factor 2 | 1.213650296 | Up | 2.63299E-07 |
SLC2A12 | Solute carrier family 7 member 12 | 1.23521218 | Up | 3.29996E-05 |
IL6 | Interleukin 6 | 1.432941749 | Up | 7.12667E-05 |
MT1G | Metallothionein 1G | -1.077752314 | Down | 1.2887E-11 |
ACSF2 | Acyl-CoA synthetase family member 2 | -1.171871909 | Down | 9.68187E-09 |
GO enrichment and KEGG pathway and Gene interrelationships analysis of differential genes related to ferroptosis in CD's
According to the results of the GO enrichment analysis(Fig. 4A,B), the biological process (BP) function was predominantly enriched in the response to oxidative stress, reactive oxygen metabolic processes, cellular response to oxidative stress, and vascular endothelial growth factor production.Changes in molecular function (MF) and cellular component (CC) were mainly in oxygen reductase activity and NADPH oxidase complex.KEGG pathway(Fig. 4A) showed that genes were mainly enriched in TNF signaling pathway and IL-17 signaling pathway. Pearson correlation test showed that some of the 13 ferroptosis-related genes had significant positive and negative correlations in CD(Figure 4C).
PPI analysis was used to show the interaction relationships between the genes of interest(Supplementary Fig. 1B,Fig. 5A).3 of the 13 genes were excluded which did not interact with the other. The other 10 genes were mapped for the PPI network.5 hub genes IL6,PTGS2,HIF1A,DUOX2, and NOS2 were discovered using the Cytoscape software's cytohubba plugin(Figure 5B).Gene co-expression, functional analysis, and physical binding properties of the five hub genes were predicted using the GeneMANIA website in the Fig. 5C.
The DSigDB database was used to predict 5 hub gene-related target drugs that may treat CD by controlling the occurrence of ferroptosis. A total of 1378 possible drugs were predicted, and the top 10 target drugs were selected after ranking based on a composite score(Table 3).
Drug Name | P-value | Odds Ratio | Combined Score |
---|---|---|---|
Zinc protoporphyrin CTD 00000915 | 2.74E-08 | 1152.057692 | 20062.49288 |
Chitosamine CTD 00006030 | 1.66E-09 | 959.6144578 | 19397.65623 |
Benzoyl Peroxide CTD 00005495 | 2.75E-06 | 1480.444444 | 18956.85145 |
Ferric Citrate CTD 00001186 | 3.30E-06 | 1332.333333 | 16817.52991 |
Heme TTD 00008407 | 3.30E-06 | 1332.333333 | 16817.52991 |
Flufenamic acid CTD 00005976 | 4.48E-08 | 966 | 16346.25293 |
Roxithromycin CTD 00007073 | 3.90E-06 | 1211.151515 | 15085.68949 |
Kaempferol CTD 00000297 | 3.85E-09 | 772.5048544 | 14966.43303 |
NSC267099 CTD 00001568 | 4.54E-06 | 1110.166667 | 13656.83282 |
Deferoxamine CTD 00005759 | 7.00E-09 | 641 | 11953.71327 |
To further validate our results, we verified the expression of 5 hub genes in 3 pairs of CD's clinical samples(Fig. 6). Among them, IL-6, PTGS2 and DUOX2 expressions were elevated and statistically significant in CD. Although the other two genes, HIF1A and NOS2, also showed some trends, the results were not statistically significant(Fig. 7).The primer information were listed in Table 4.
Primers | Sequence(5’→3’) |
---|---|
IL6-F | ACTCACCTCTTCAGAACGAATTG |
IL6-R | CCATCTTTGGAAGGTTCAGGTTG |
PTGS2-F | CTGGCGCTCAGCCATACAG |
PTGS2-R | CGCACTTATACTGGTCAAATCCC |
DUOX2-F | CTGGGTCCATCGGGCAATC |
DUOX2-R | GTCGGCGTAATTGGCTGGTA |
HIF1A | GAACGTCGAAAAGAAAAGTCTCG |
HIF1A | CCTTATCAAGATGCGAACTCACA |
NOS2-F | TTCAGTATCACAACCTCAGCAAG |
NOS2-R | TGGACCTGCAAGTTAAAATCCC |
CD is a long-term, chronic, recurring intestinal disease, collectively known with ulcerative colitis as inflammatory bowel disease(IBD)12.The incidence of CD was more pronounced in the Americas and Europe before the 20th century, however, the prevalence of CD has been steadily rising in Asian countries in recent years13.It mostly develops at the end of the ileum and can cause lesions from the mouth to the anus, as well as extraintestinal complications. CD were usually clinically insidious and are not easily detected and diagnosed, patients are often overlooked and missed when having a colonoscopy14. Patients' symptoms are dominated by diarrhea, abdominal pain, fever and anemia.Extraintestinal symptoms include inflammatory arthropathy, osteoporosis, scleritis, IgA nephropathy and others1516.Review of the Pubmed database revealed that the current findings on risk factors for CD include environmental factors, genetic factors, nutritional status,intestinal barrier disorders, immune response, microbial dysbiosis, gut viral and fungal317.However, the specific mechanisms or involving pathways are not clear yet.
Ferroptosis is a newly discovered iron-dependent, a non-apoptotic form of cell death, usually caused by intracellular lipid peroxidation and regulated by the oxidative and antioxidant systems in vivo, including redox homeostasis, mitochondrial activation, amino acid lipid and sugar metabolism, and other signaling pathways7. It has been observed to have a modulating impact in a number of benign and malignant disorders, including tumours, neurological disorders, kidney damage and Cardiovascular disease181920.
Ferroptosis has an important relationship with intestinal barrier and immune response, which are also the key mechanism causing CD.Liu et.al found that ferroptosis reverses intestinal mucosal barrier integrity maintained by Sestrin221.Ma et.al announced that CD36-directed CD8+ T cells ferroptosis suppresses the anti-tumour effects in vivo22.Xu et.al suggested that ferroptosis could exert anti-tumour effects by reversing the immune microenvironment23. For this reason we further speculate that there is some connection between ferroptosis and CD.
Based on the GEO databases and the ferrdb database, we obtained a total of 13 differential genes associated with ferroptosis and further explored the related functions of these genes by GO and KEGG enrichment analysis.As the figure shows in GO biological processes(BP) are most related to oxidative stress including response to oxidative stress,ROSmetabolic process and cellular response to oxidative stress.This result is more in line with our speculation, as oxidative stress processes are a key part of the occurrence of ferroptosis2425. In the course of the metabolic pathway of the cell due to the imbalance of the redox balance eventually leads to ferroptosis of the cell.Oxidative stress is also an important factor in inflammation, which activates a variety of proinflammatory factors such as NF-κB and p5326. The increased expression of P53 in Crohn promotes apoptosis of intestinal epithelial cells and aggravates inflammation2728. We hypothesized that oxidative stress-related ferroptosis may also be involved in the development of CD.
At the same time, we also obtained another BP enrichment result-vascular endothelial growth factor (VEGF) production.VEGF is a biomolecule that It is crucial for mitosis and preventing apoptosis in vascular endothelial cells29, including several members such as VEGF-A, VEGF-B, and VEGF-C.Inflammation and tumorigenesis are often accompanied by pathological vascular proliferation, which also implies an imbalance of oxygen homeostasis in the microenvironment.The concentration of ROS dependent on NADPH oxidase production is also significantly increased, which seems to form a circular pathway30. High production of ROS increases the likelihood of ferroptosis in cells31.This is also consistent with the enrichment results we obtained in Molecular function (MF) and cellular component (CC).
In other studies, the biomodulatory role of VEGF in IBD has been partially reported in other research.Scaldaferri et.al found that overexpression of VEGF-A in mice with DSS-induced IBD inflammation can exacerbate disease severity32. According to Eder et al,the degree of VEGF expression affected how severe CD was.33. However, the existing mechanism between VEGF and ferroptosis is still unclear. It is initially speculated that VEGF factors may lead to ferroptosis of intestinal epithelial cells and promote the occurrence of Crohn.
KEGG pathway analysis further revealed that the pertinent genes were primarily abundant in the TNF pathway, IL-17 pathway and Rheumatoid arthritis.Tumour necrosis factor (TNF) is an important cell regulatory molecule involved in the regulation of cell survival, apoptosis, as well as inflammation, and immune response34. Interleukin 17 (IL-7) is a inflammatory factor secreted by T helper 17 (Th17) cells35. Bishu et.al found elevated expression of both TNF and IL-17 in CD, which may be associated with a type CD 4+ T cell36.However, the regulatory relationship between ferroptosis with TNF factors and IL-17 in the currently available studies of CD remains unclear.
Next,We mapped the PPI network in Cytoscape and screened 5 hub genes using CytoHubba software IL6、PTGS2、HIF1A、DUOX2、NOS2. Three pairs of samples were performed qRC-PCR experiments to analyze the expression of the genes according to these five genes, in which the expression of IL6, DUOX2 and PTGS2 were different and p < 0.05.
Interleukin-6 (IL-6) is a pro-inflammatory factor that belongs to the interleukin family and crucial in the host's body for anti-inflammatory and emergency processes37.Shi et.al suggested that Toll-like receptor 4 ( TLR4) aggravates intestinal injury by regulating the expression of IL-638.In a study by V Gross et.al, who found significantly elevated IL-6 in CD in both active and inactive phases39.IL-6 was not only play a role in inflammation.Han et.al found that in asthma disease IL-6 disrupts iron homeostasis in BEAS-2B cells by promoting lipid peroxidation, which in turn promotes ferroptosis in bronchial epithelial cells40. Li et.al recently demonstrated that the expression of an amino acid countertransporter (xCT) in head and neck squamous cell carcinoma is associated with levels of IL-6, and xCT induces cellular ferroptosis41.Despite various results suggesting a link between IL-6 and inflammation and ferroptosis, there is no conclusive evidence for an association between IL-6 and ferroptosis in CD.
Prostaglandin endoperoxide synthase 2 (PTGS2), also called cyclooxygenase(COX-2), an enzyme that can be suppressed by nonsteroidal anti-inflammatory drugs (NSAID) such as aspirin and ibuprofen played a regulatory role in many types of inflammatory or neoplastic diseases42.Roberts et al. found through densitometric analysis that the expression of COX-2 in IBD tissues increased by 6–8 times43. Singer et al. Et al. found that COX-2 was predominantly present in apical epithelial cells and in lamina propria mononuclear cells in CD. Interestingly, both studies showed no expression of COX-2 in normal tissues44.This demonstrates that COX-2 has an important regulatory role in CD. However, no relevant studies have found whether the overexpression of COX-2 in Crohn is associated with ferroptosis.
DUOX2,Dual oxidase 2,is an enzyme first identified in the mammalian thyroid and a member of NADPH oxidase involved in ROS production.In a study of pediatric CD.Haberman et.al found that within the villous enterocytes of the CD ileum, DUOX2 was strongly expressed and may have exacerbated intestinal mucosal damage45.Lipinski et.al's study found that DUOX2 is involved in the production of ROS46, ROS is an important factor mediating the occurrence of ferroptosis.We, therefore, speculate that DUOX2 may promote ferroptosis in CD. However, the specific mechanism is not yet clear.
In the result of the 5 hub genes drug predicting targets,Kaempferol had found to inhibit the onset of cellular iron death, while other studies have shown that it can also be used to improve intestinal inflammation and intestinal barrier disorders4748.
Our research has some drawbacks.Firstly,the data set in this study was obtained from the GEO database and the same microarray was used for both data sets, and there was no further validation in tissue samples or cells or animals. Therefore, it is hoped that the association between iron death and CD will be further explored in future studies.
Overall, 13 ferroptosis-associated genes were identified in CD by our bioinformatics study. Among them IL-6,DUOX2 and PTGS2 are most likely to regulate ferroptosis in CD involved in CD development and progression. But specific mechanistic studies need to be further explored. Our study provides theoretical support for CD and ferroptosis, increasing our understanding of CD and treatment strategies.
Funding
The authors declare that no funds, grants, or other support were received during the preparation of this manuscript.
Competing Interests
The authors have no relevant financial or non-financial interests to disclose.
Author Contributions
All authors contributed to the study conception and design.Material preparation, data collection and analysis were performed by WenquanZhang and Zhaoshui Li. The first draft of the manuscript was written by Wenquan Zhang and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.
Hongbo Li conceived the concept, designed the manuscript, coordinated, and critically revised manuscript, and was responsible for its financial supports and the corresponding works.
Ethics Approval
This study was performed in line with the principles of the Declaration of Helsinki. Approval was granted by the Ethics Committee of Qingdao municipal hospital.(2022-066)
Consent to participate
Informed consent was obtained from all individual participants included in the study.
Availability of data and materials
The datasets presented in this study can be found in online repositories. Further inquiries can be directed to the corresponding authors.