ABSTRACT
Background Ulcerative colitis (UC) and irritable bowel syndrome (IBS) are both chronic bowel diseases involving stress. To identify genes differentially expressed in UC and IBS, and to determine whether psychological stress can influence those gene expressions, we conducted this pilot study.
Methods Patients of UC, irritable bowel syndrome (IBS) and normal controls (NC) (n=5 in each group) were recruited. Perceived stress scale (PSS) was adopted to assess psychological stress level. Sigmoid mucosa samples were collected during colonoscopy. Differentially expressed genes (DEGs) of both mRNA and microRNA (DEG-mRNA and DEG-miRNA) were identified by RNA-Seq and microarray, respectively. Weighted gene co-expression network analysis (WGCNA), gene ontology, and microRNA target analysis were performed to identify regulatory relationships and pathways involved. Pearson correlation was performed to identify the relationship between DEGs and PSS score.
Results In total, 1,770 and 938 DEG-mRNAs, 107 and 3 DEG-miRNAs were identified in UC and IBS (nominal P<0.05), respectively. The transcriptome changes in UC and IBS were highly correlated but transcriptome severity was larger in UC. There were 268 overlapped DEG-mRNAs and 3 overlapped DEG-miRNAs between UC and IBS. Median PSS score was 27.5 (24.7, 40.0) in UC, 27.0 (17.5, 32.5) in IBS and 22.0 (12.0, 28.5) in NC. Five of the UC’s DEG-mRNAs were significantly correlated to PSS scores, including SLC4A7, PLCB1, SPATA33, DECR2, RP11-792A8.4. There were no DEGs correlate with PSS scores in IBS. One module enriched for immunological pathways was upregulated in UC. Metabolic pathways were enriched in IBS. Within the UC-related co-expression module, we identified 18 DEG-mRNAs that were also targets of DEG-miRNA.
Conclusions Integration of microRNA and mRNA expression profiles identified changes of gene expression in both UC and IBS. Although the two diseases have significant overlapped transcriptomic changes, they have distinct signature. Immunological pathways are involved prominently in UC while metabolic pathways in IBS. Psychological stress is uniquely associated with the transcriptome changes of UC.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
This work was supported by Scientific research projects of public welfare industry (Grant No. 201002020, 201502005) and Medical science and health technology innovation projects of the Chinese academy of medical science (Grant No. 2016-I2M-3-001, 2019-I2M-2-007).
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
The study protocols were submitted to, and approved by the ethics committee of Peking Union Medical College Hospital (PUMCH). The approval number was S-723. Informed consent was obtained from all subjects or, if subjects are under 18, from a parent. All methods were carried out in accordance with relevant guidelines and regulations along with the approval.
All necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.
Yes
Data Availability
The sequencing data generated in this study has been deposited in the NCBI SRA database (https://www.ncbi.nlm.nih.gov/sra/PRJNA674301).
The microarray data can be downloaded from the NCBI GEO database at http://www.ncbi.nlm.nih.gov/geo/ using accession number GSE160900.
LIST OF ABBREVIATIONS
- UC
- Ulcerative colitis
- IBS
- Irritable bowel syndrome
- IBS-C
- Constipation-predominant IBS
- IBS-D
- Diarrhea-predominant IBS
- PSS
- Perceived stress scale
- DEG
- Differentially expressed gene
- mRNA
- Message RNA
- miRNA
- microRNA
- WGCNA
- Weighted gene co-expression network analysis
- GO
- Gene ontology
- GWAS
- Genome-wide-associations studies
- IBD
- Inflammatory bowel disease
- HLA-DRB1
- Human leukocyte antigen-DR beta 1
- TAP2
- Transporter 2, ATP binding cassette subfamily B member
- INADL
- InaD-like gene
- PUMCH
- Peking Union Medical College Hospital
- NC
- Normal control
- RSEM
- Robust Structural Equation Modeling
- TPM
- Transcripts per million
- PCA
- Principle component analysis
- Log2FC
- Log2 fold-change
- FDR
- False discovery rate
- ME
- Module eigengene
- WebGastalt
- WEB-based Gene SeT AnaLysis Toolkit
- IQR
- Interquartile range
- SLC4A7
- Solute carrier family 4 member 7
- PLCB1
- Phospholipase C Beta 1
- SPATA33
- Spermatogenesis Associated 33
- DECR2
- 2,4-Dienoyl-CoA Reductase 2
- RP11-792A8.4 AHR
- Aryl Hydrocarbon Receptor
- GPR18
- G Protein-Coupled Receptor 18
- CD28
- CD28 Molecule
- FAM26F
- Family With Sequence Similarity 26, Member F
- AC019330.1 LYRM9
- LYR Motif Containing 9
- VCAM1
- Vascular Cell Adhesion Molecule 1
- TNFSF8
- Tumor Necrosis Factor (Ligand) Superfamily, Member 8
- TNFSF14
- Tumor Necrosis Factor (Ligand) Superfamily, Member 14
- STAT2
- Signal Transducer And Activator Of Transcription 2
- SP140L
- SP140 Nuclear Body Protein Like
- PPP1R16B
- Protein Phosphatase 1 Regulatory Subunit 16B
- LPXN
- Leupaxin
- LILRB2
- Leukocyte Immunoglobulin Like Receptor B2
- IKZF3
- IKAROS Family Zinc Finger 3
- GNB4
- G Protein Subunit Beta 4
- DUSP2
- Dual Specificity Phosphatase 2
- CYTIP
- Cytohesin 1 Interacting Protein
- CLEC7A
- C-Type Lectin Domain Containing 7A
- CIITA
- Class II Major Histocompatibility Complex Transactivator
- CD274
- CD274 Molecule(同时也是 PD-L1)
- CD180
- CD180 Molecule
- CARD8
- Caspase Recruitment Domain Family Member 8
- APOL2
- Apolipoprotein L2
- CFTR
- Cystic fibrosis transmembrane conductance regulator