Investigating the change in gene expression profile of blood mononuclear cells post-laparoscopic sleeve gastrectomy in Chinese obese patients
Authors: Liu Na , Chen Xiaolei , Ran Jianghua , Yin Jianhui , Zhang Lijun , Yang Yuelin , Cen Jianchang , Dai Hongmei , Zhou Jiali , Gao Kui , Zhang Jihong , Liu Liyin , Chen Zhiyuan , Wang Haibin 
Published: March 2023

Background

Laparoscopic sleeve gastrectomy (LSG) is a sustainable technique that effectively treats morbid obesity. However, the molecular mechanisms underlying the improvement of metabolic health following this process warrants more investigation. This study investigates LSG-related molecules and uses bulk RNA-sequencing high-throughput analysis to unravel their regulatory mechanisms.

Methods

Peripheral blood mononuclear cells (PBMC) were collected from ten obese patients with BMI ≥ 32.5 kg/m2 in the Department of General Surgery of Kunming First People’s Hospital. After LSG, patients were followed up for one month, and blood samples were retaken. Blood samples from ten patients before and after LSG and bulk RNA-Seq data were analyzed in this study. LSG-associated gene expression was detected by weighted gene coexpression network analysis (WGCNA) and differential analysis. Subsequently, essential signature genes were identified using logistic least absolute shrinkage and selection operator (LASSO) and support vector machine-recursive feature elimination (SVM-RFE) algorithms. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and single-sample gene set enrichment analysis (ssGSEA) were utilized to reveal the potential functions of the target genes. Furthermore, the Pearson correlation of signature genes with leptin and lipocalin was also explored. Finally, we constructed a robust endogenous RNA (ceRNA) network based on miRWalk and starBase databases.

Results

We identified 18 overlapping genes from 91 hub genes, and 165 differentially expressed mRNAs (DE-mRNA), which were revealed to be significantly associated with immune cells, immune response, inflammatory response, lipid storage, and localization upon functional enrichment analysis. Three signature genes, IRF1, NFKBIA, and YRDC, were identified from the 18 overlapping genes by LASSO and SVM-REF algorithms. The logistic regression model based on the three signature genes highlighted how robustly they discriminated between samples. ssGSEA indicated these genes to be involved in lipid metabolism and degradation pathways. Moreover, leptin levels were significantly reduced in patients undergoing LSG, and NFKBIA significantly negatively correlated with leptin. Finally, we identified how the long non-coding RNA (lncRNA) ATP2B1-AS1 regulated the expression of the signature genes by competitively binding to six microRNAs (miRNAs), which were hsa-miR-6509-5p, hsa-miR-330-5P, hsa-miR-154-5P, hsa-miR-145-5P, hsa-miR4726-5P and hsa-miR-134-5P.

Conclusion

This study identified three critical regulatory genes significantly differentiated between patients before and after LSG treatment and highlighted their potentially crucial role after bariatric surgery. This provides novel insights to increase our understanding of the underlying mechanisms of weight loss and associated metabolic improvement after bariatric surgery.