-
随着现代社会生活方式和饮食结构的改变，肥胖（obesity）的发生率不断升高，目前肥胖及其相关并发症已经成为世界性的公共卫生问题。2020年发布的《中国居民营养与慢性病状况报告》显示，我国超重人口已达3.25亿，占据成年总人口的30.1%，其中肥胖人群数量接近9 000万[1]。当摄入能量大于消耗能量时，机体会将过剩能量转变为脂肪进行贮存，进而导致肥胖的发生发展。进一步根据肥胖症患者的相关代谢指标水平（如血糖、血压和血脂水平等）进行分类，可将其分为代谢健康的肥胖患者（metaboli⁃ cally healthy obesity，MHO）[2-3] 和代谢不健康的肥胖患者（metabolically unhealthy obesity，MUO）[4]。除生活方式和药物干预外，减重代谢手术已成为目前治疗重度肥胖及其并发症的有效手段[5-6]。
-
脂肪组织作为一个复杂的内分泌器官，除了脂肪细胞外，还包含了免疫细胞、血管基质细胞、结缔组织和神经组织等多种成分。既往研究表明，增生肥大的脂肪组织中通常伴随各类免疫细胞的异常浸润，导致肥胖症患者体内持续的低水平慢性炎症状态，并影响机体的糖脂代谢能力，造成相关临床指标的异常改变，进而引起相关并发症，例如胰岛素抵抗（insulin resistance，IR）[7]、心血管疾病[8-10] 以及非酒精性脂肪肝（nonalcoholic fatty liver disease， NAFLD）[11-12] 等。
-
根据部位的不同，脂肪组织主要分为皮下脂肪组织（subcutaneous adipose tissue，SAT）和内脏脂肪组织（visceral adipose tissue，VAT）两大类。VAT 中的固有免疫细胞群与肥胖状态下脂肪组织稳态及低水平炎症的调节有关，并伴有数量和功能的改变[13]。针对肥胖人群的研究证实，VAT中巨噬细胞浸润最为明显，相关炎症因子水平与VAT中巨噬细胞水平呈正相关，并且糖代谢受损亚组人群脂肪组织表现为巨噬细胞优先向VAT浸润的表型[14]；另有研究显示，VAT分布对于代谢综合征（metabolic syn⁃drome，MS）的发生具有预测作用[15-16]，均提示 VAT 的异常扩增伴随炎症细胞浸润可能是影响不同肥胖个体代谢功能的关键因素。近期相关研究发现，特殊类型的B细胞所产生的IgG可能是介导肥胖炎症发展的关键因子，再次强调了免疫细胞在肥胖相关炎症发生过程中的重要作用[17]。
-
本课题组前期的脂肪组织RNA⁃seq研究结果表明，与体重正常（lean）组相比，MHO 组及MUO 组间 VAT差异表达基因富集于炎症免疫通路[18]，从而引起思考：不同体重人群的VAT内免疫细胞浸润情况是否存在差异？其与代谢相关差异表达基因是否存在关联？故本研究针对以上3组间VAT中各类免疫细胞占比及变化情况进行探究，并通过MHO组及 MUO 组间差异表达的免疫细胞与 RNA⁃seq 中差异表达的炎症及免疫反应通路基因进行关联性分析，旨在评估肥胖患者代谢能力受损与各类免疫细胞占比改变的相关性。
-
1 对象和方法
-
1.1 对象
-
收集于南京医科大学第一附属医院因肥胖行代谢手术治疗的48例患者及体重和代谢指标正常但因其他非炎症性疾病（疝气、胆结石及肾囊肿等） 行腹腔镜手术的11例患者的内脏脂肪组织标本，同时汇总了患者相关病史及临床检查数据。该项目获得南京医科大学第一附属医院伦理委员会批准 （2016⁃SR⁃220），且所有标本的收集均在术前获得患者书面的知情同意。
-
MUO组的定义[4]：患者体重指数（body mass in⁃ dex，BMI）≥28 kg/m2，且至少满足以下标准中的 3 项：空腹血糖≥5.6 mmol/L，收缩压≥130 mmHg，舒张压≥85 mmHg，甘油三酯（triacylglycerol，TAG）＞1.7 mmol/L，高密度脂蛋白胆固醇（high⁃density lipopro⁃ tein cholesterol，HDL ⁃ c）＜1.3 mmol/L（女性）或＜1.04 mmol/L（男性）。MHO 组的定义[2，19]：患者 BMI≥28 kg/m2，空腹血糖 ＜5.6 mmol/L，收缩压 ＜130 mmHg，舒张压＜85 mmHg，TAG＜1.7 mmol/L，且既往无高脂血症、高血压、糖尿病等病史。Lean 组的定义：患者BMI＜28 kg/m2。
-
1.2 方法
-
1.2.1 RNA提取及RNA⁃seq测序
-
收集的内脏脂肪组织标本均通过TRIzol法提取 RNA，使用Qiagen RNeasy 试剂盒纯化后，通过琼脂糖凝胶电泳、Nanodrop及Agilent2100确认每个样品的 RNA 质量。基于 Illumina 平台的 NEBNext Ultra RNA Library Prep Kit 制备每个样品的 RNA⁃seq 库，并在Illumina HiSeq 4000平台上进行RNA⁃seq。
-
1.2.2 免疫细胞构成比计算
-
利用在线数据库 CIBERSORT（https：//cibersortx.stanford.edu/）[20]，对 RNA⁃seq 数据进行基于线性支持向量回归原理的反卷积分析，可计算获得22 种不同类型和功能的免疫细胞在内脏脂肪组织样本中的占比分数。
-
1.3 统计学方法
-
应用 SPSS 25 软件进行相关数据的统计分析，并对各项指标进行正态分布和方差齐性检验。3组间比较应用单因素方差分析，事后比较应用邓尼特t 检验；相关性分析应用 Pearson 相关分析。P＜0.05 为差异有统计学意义。图片由GraphPad Prism 8软件绘制。
-
2 结果
-
2.1 3组人群分组情况及相关临床指标
-
首先对相关患者的代谢特征数据（如年龄、性别比、BMI及空腹血糖等）按照分类标准进行统计，最终符合标准的Lean组、MHO组及MUO组入组人数分别为11例、22例及26例，其具体临床数据见表1。
-
2.2 3组人群内脏脂肪组织中的免疫细胞占比情况
-
将 3 组人群的内脏脂肪标本提取 RNA 并进行 RNA⁃seq测序，利用在线数据库CIBERSORT进行数据处理分析，可以推算得到共22种免疫细胞在VAT 中的占比情况。
-
HbA1c：糖化血红蛋白；TC：总胆固醇；TAG：甘油三酯；HDL⁃c：高密度脂蛋白胆固醇；LDL⁃c：低密度脂蛋白胆固醇。
-
与Lean组比较，MHO及MUO组人群VAT中的巨噬细胞占比变化明显，其中以 M2 型巨噬细胞表达增加为主（图1）。
-
通过分析各类免疫细胞在组间的表达差异，发现调节性 T 细胞（regulatory T cell，Treg）在 3 组间的比例呈依次递减的趋势[Lean 组（7.08±4.04）％， MHO 组（5.69±3.45）％，MUO 组（4.08±2.59）％]，且 Lean组与MUO组的Treg比例存在统计学差异（Lean 组与MUO组比较，P=0.031，图2A）。而被激活的自然杀伤（natural killer，NK）细胞则在 Lean 组中占比较高，为（4.53±2.80）%，在肥胖的两组人群间比较， MUO 组被激活的 NK 细胞占比有高于 MHO 组的趋势，但两组间无统计学差异[MHO组（1.78±1.73）％， MUO 组（2.32±1.87）％，MHO 组与 MUO 组比较，P= 0.361，图2B]。免疫细胞中的巨噬细胞亚群中，M1 型巨噬细胞的占比在3组间没有明显改变[Lean组（0.04 ± 0.14）％，MHO 组（0.11 ± 0.42）％，MUO 组 （0.21±0.38）％，图2C]。M0 型巨噬细胞的占比在 MUO组明显升高，而在代谢正常的两组间没有明显差异[Lean 组（2.69 ± 4.45）％，MHO 组（3.13 ± 7.46）％，MUO 组（12.00±12.79）％，Lean 组与 MUO 组比较，P=0.011；MHO组与MUO组比较，P=0.003，图2D]。发挥抗炎作用的M2型巨噬细胞占比的改变最为明显，无论是MHO组还是MUO组，对比Lean组，M2 型巨噬细胞的占比都明显增加[（Lean 组 （20.63±7.72）％，MHO 组（28.54±9.76）％，MUO 组 （28.51±7.19）％，Lean 组与 MHO 组比较，P=0.013； Lean组与MUO组比较，P=0.011，图2E]。
-
图1 3组人群内脏脂肪组织RNA⁃seq分析所得相关免疫细胞占比情况
-
Figure1 Proportion of immune cells in VAT among three groups from RNA⁃seq analysis
-
通过分析各亚型巨噬细胞在整群巨噬细胞中的占比（图3）发现，相较MHO组，MUO组VAT中的 M2 型巨噬细胞在巨噬细胞整群中占比明显下降 [Lean 组 ｛88.99 ± 15.01）％，MHO 组（90.75 ± 18.79）％，MUO组（73.82±23.55）％，MHO组与MUO 组比较，P=0.012]，而促炎的M1型巨噬细胞在巨噬细胞整群中占比则有相对升高的趋势[Lean 组 （0.14 ± 0.46）％，MHO 组（0.27 ± 0.92）％，MUO 组 （0.46±0.86）%]，提示MUO组内VAT中存在较高的炎症水平，且这种持续慢性低炎症水平状态与肥胖患者代谢紊乱的发生有关。
-
2.3 MHO与MUO组间GPR183和C5AR1与巨噬细胞亚群水平存在相反关联性
-
在既往对相关RNA⁃seq数据的分析统计中，已经明确了内脏脂肪组织中存在30余种基因在Lean 组、MHO组及MUO组间皆存在差异性表达，尤其是部分炎症通路和免疫反应相关基因，如 C5AR1、 GPR183、IL1RN、ITGAX、CD300LB及MSR1等，与空腹血糖、收缩压等代谢指标关联[18]，提示机体内部炎症反应与代谢健康状态具有一定相关性。利用 The Human Protein Atlas 数据库查询 C5AR1 等差异表达基因在脂肪组织中的单细胞测序结果，可以发现这些基因多数仅在以巨噬细胞为主的免疫细胞中存在高表达（图4）。
-
图2 3组人群内脏脂肪中差异表达的免疫细胞亚群
-
Figure2 Differentially expressed immune cell subclasses in VAT among three groups
-
图3 3组人群内脏脂肪中M1（A）和M2型（B）巨噬细胞在总巨噬细胞中的占比
-
Figure3 The proportion of M1（A）and M2（B）macro⁃ phages in total macrophages in VAT among three groups
-
接着，又将已明确的31种在3组间差异表达的基因[18] 与有差异的免疫细胞（Treg细胞，激活的NK 细胞，M0、M1及M2型巨噬细胞）占比进行相关性分析（图5）。在MHO与MUO组中，与代谢指标无相关性的差异表达基因与巨噬细胞亚类的关联性相对较弱，与代谢指标具有相关性的基因多与巨噬细胞亚类存在较强关联。进一步地对与代谢指标相关的基因进行分析，发现数种基因的表达水平，如 SPP1、CHIT1、ITGAX、CD300LB、IL1RN 及 DC⁃ STAMP 等，与 M0、M1 型巨噬细胞的占比呈正相关关系，而与 M2 型巨噬细胞的占比大多数呈负相关关系。但MHO与MUO两组间的不同之处在于C5a 受体 1（Complement fragment 5a receptor 1，C5AR1） 与 G 蛋白偶联受体 183（G ⁃ protein coupled receptor 183，GPR183）基因与巨噬细胞亚型的关联性，这两种基因在MHO组中表现为与M0及M1型巨噬细胞占比正相关，与M2型巨噬细胞占比负相关；在MUO 组中则呈现相反的关联性。
-
3 讨论
-
结果证明，相较体重正常组，两组肥胖患者的巨噬细胞整体表达水平都明显提高；同时在两组肥胖患者间比较巨噬细胞亚群的整体占比，抗炎的M2 型巨噬细胞的比例在MUO组中较MHO组下降，提示代谢异常的肥胖患者体内的炎症水平相对更高。造成这种改变的具体机制为肥胖时增生的脂肪组织通过释放白介素⁃6（interleukin⁃6，IL⁃6）、白介素⁃1β（interleukin⁃1β，IL⁃1β）、肿瘤坏死因子α（tu⁃ mor necrosis factor⁃α，TNF⁃α）及单核细胞趋化蛋白⁃1 （monocyte chemotactic protein⁃1，MCP⁃1）等炎症因子和细胞因子，激活并募集以巨噬细胞为主的免疫细胞[21-22]，成为病理性免疫反应的激活部位，表现为肥胖机体内的慢性持续性低水平炎症状态，最终导致机体代谢水平改变。
-
图4 脂肪组织单细胞测序中C5AR1、GPR183、IL1RN、ITGAX、CD300LB和MSR1的表达水平
-
Figure4 Expression levels of C5AR1、GPR183、IL1RN、ITGAX、CD300LB and MSR1 in single⁃cell sequencing of adipose tissue
-
图5 3组人群中免疫反应及炎症通路基因表达与免疫细胞比例的相关性分析
-
Figure5 The correlation of gene expression of immune response and inflammation and proportion of immune cells among three groups
-
除了巨噬细胞亚群外，Treg表达水平在MUO组中下降明显。Treg主要在免疫反应中调节效应性T 细胞和B细胞的作用以维持免疫稳态，同时也通过分泌细胞因子白介素⁃10（interleukin⁃10，IL⁃10）等参与抗炎调节[23]。Treg被证明与促炎型巨噬细胞存在拮抗作用并可改善胰岛素抵抗[24-26]，故肥胖状态下 Treg数量与功能下调与代谢紊乱的发生较为相关。
-
在激活的NK细胞分组中，MUO组相较MHO组的细胞占比有上升趋势。NK细胞与巨噬细胞在功能上存在双向调节作用[27]，并且Lee等[28] 的动物研究证明，在肥胖状态下NK细胞能够调节巨噬细胞向促炎方向极化，从而引起胰岛素抵抗等代谢状态紊乱的发生。但本研究中，激活的NK细胞水平在体重正常的对照组较肥胖组别更高，可能是由于部分体重正常的对照者术中取材部位源自于肿块周边的脂肪组织，进而引起了激活的NK细胞水平表达差异。
-
在前期研究中发现，3 组间的差异表达基因主要集中于免疫反应与炎症通路，其中部分基因与代谢参数相关[18]，且多数仅在脂肪组织的巨噬细胞等免疫细胞中高表达。继续对肥胖状态下上调的炎症通路基因与 VAT 中差异表达的免疫细胞占比进行关联分析验证，结果提示肥胖组中SPP1、CHIT1、 ITGAX、CD300LB、IL1RN 及 DCSTAMP 等差异表达基因与VAT中发挥促炎作用的M1型巨噬细胞占比正相关，与发挥抗炎作用的 M2 型巨噬细胞占比负相关，再次印证了肥胖个体VAT中巨噬细胞浸润增加，进而导致组织内高炎症水平状态。
-
与此同时，C5AR1 及 GPR183 作为同样在巨噬细胞中高表达并且与代谢指标相关的差异表达基因，在MHO与MUO组间，对二者与VAT中巨噬细胞占比进行统计却表现出了相反的相关性，这种相反的趋势提示我们思考：这两种基因是否存在作为预测肥胖患者代谢状态改变的相关因子的可能性？
-
C5AR1为补体C5a的对应受体，可以激活多条重要信号通路途径以正向调控炎症反应[29-31]。既往动物实验证明，饮食诱导肥胖小鼠模型的脂肪组织中C5AR1的表达水平明显升高；全身敲除C5AR1的小鼠在高脂喂养的条件下与对照组相比胰岛素敏感性得到改善，这可能与敲除C5AR1后脂肪组织中 M1型巨噬细胞数量减少及抗炎因子IL⁃10的表达增加有关[32]。在高脂饮食诱导下，全身敲除C5AR1的小鼠性腺脂肪及腹股沟脂肪质量较野生型小鼠明显减少，伴随更高的餐后脂质清除率和较低的血浆甘油三酯及游离脂肪酸水平，该过程可能与C5AR1 及C5AR样受体（C5aR⁃like receptor 2，C5L2）具有互作关系并形成异二聚体有关[33-34]。
-
GPR183作为在多种免疫细胞表面表达的一种 G蛋白偶联受体，与免疫反应的激活与免疫细胞的迁移定位有直接关系[35]。Copperi 等[36] 最新研究发现，GPR183 在棕色脂肪组织（brown adipose tissue， BAT）中高表达，通过应用其内源性配体7α，25⁃二羟基胆固醇（7α，25⁃dihydroxycholestero）可以激活BAT 中的 GPR183 从而抑制 BAT 产热功能，降低能量消耗，但GPR183在白色脂肪组织（white adipose tissue， WAT）中发挥何种作用尚未明确。故 C5AR1 和 GPR183 确可影响脂肪组织功能，但这二者在脂肪组织扩增过程中是如何引起机体代谢能力改变的，其具体机制仍需要更深入的基础实验进行验证。
-
综上，本文通过RNA⁃seq技术分析了正常体重者与不同代谢状态的肥胖症患者的 VAT 中各型免疫细胞占比变化与部分免疫反应通路基因的关联性，再次验证了代谢正常及异常的肥胖患者内脏脂肪组织中免疫细胞亚类的构成变化情况，并且分析得到C5AR1及GPR183基因与肥胖患者巨噬细胞亚群占比及代谢状态变化有关，具有作为预测肥胖代谢状态生物学因子的潜能。
-
参考文献
-
[1] 中国居民营养与慢性病状况报告（2020年）［J］.营养学报，2020，42（06）：521
-
[2] SMITH G I，MITTENDORFER B，KLEIN S.Metabolical⁃ ly healthy obesity：facts and fantasies［J］.J Clin Invest，2019，129（10）：3978-3989
-
[3] BLÜHER M.Metabolically healthy obesity［J］.Endocr Rev，2020，41（3）：bnaa004
-
[4] PHILLIPS C M.Metabolically healthy obesity：definitions，determinants and clinical implications［J］.Rev Endocr Metab Disord，2013，14（3）：219-227
-
[5] CUMMINGS D E，RUBINO F.Metabolic surgery for the treatment of type 2 diabetes in obese individuals［J］.Dia⁃ betologia，2018，61（2）：257-264
-
[6] 韩涛，冯文焕.儿童青少年肥胖病代谢手术应用进展 ［J］.南京医科大学学报（自然科学版），2022，42（10）：1482-1487
-
[7] BURHANS M S，HAGMAN D K，KUZMA J N，et al.Con⁃ tribution of Adipose Tissue Inflammation to the Develop⁃ ment of Type 2 Diabetes Mellitus［J］.Compr Physiol，2018，9（1）：1-58
-
[8] NEELAND I J，ROSS R，DESPRÉS J P，et al.Visceral and ectopic fat，atherosclerosis，and cardiometabolic dis⁃ ease：a position statement［J］.Lancet Diabetes Endocri⁃ nol，2019，7（9）：715-725
-
[9] PÉREZ⁃PEVIDA B，DÍAZ⁃GUTIÉRREZ J，MIRAS A D，et al.High body adiposity drives glucose intolerance and increases cardiovascular risk in normoglycemic subjects ［J］.Obesity（Silver Spring，Md），2018，26（4）：672-682
-
[10] GAO M，LV J，YU C，et al.Metabolically healthy obesity，transition to unhealthy metabolic status，and vascular dis⁃ ease in Chinese adults：a cohort study［J］.PLoS Med，2020，17（10）：e1003351
-
[11] ZHANG X，WU M，LIU Z，et al.Increasing prevalence of NAFLD/NASH among children，adolescents and young adults from 1990 to 2017：a population⁃based observation⁃ al study［J］.BMJ Open，2021，11（5）：e042843
-
[12] LU F，HU E，XU L，et al.The relationship between obesi⁃ ty and the severity of non⁃alcoholic fatty liver disease：sys⁃ tematic review and meta⁃analysis［J］.Expert Rev Gastro⁃ enterol Hepatol，2018，12（5）：491-502
-
[13] MAGNUSON A M，REGAN D P，FOUTS J K，et al.Diet⁃ induced obesity causes visceral，but not subcutaneous，lymph node hyperplasia via increases in specific immune cell populations［J］.Cell Prolif，2017，50（5）：e12365
-
[14] HARMAN⁃BOEHM I，BLÜHER M，REDEL H，et al.Mac⁃ rophage infiltration into omental versus subcutaneous fat across different populations：effect of regional adiposity and the comorbidities of obesity［J］.J Clin Endocrinol Metab，2007，92（6）：2240-2247
-
[15] KANG Y M，JUNG C H，CHO Y K，et al.Visceral adiposi⁃ ty index predicts the conversion of metabolically healthy obesity to an unhealthy phenotype［J］.PLoS One，2017，12（6）：e0179635
-
[16] LV F，CAI X，LI Y，et al.Association between indices of body composition and metabolically unhealthy phenotype in China：a cross ⁃ sectional study［J］.Front Endocrinol（Lausanne），2022，13：891327
-
[17] HÄGGLÖF T，VANZ C，KUMAGAI A，et al.T⁃bet B cells accumulate in adipose tissue and exacerbate metabolic disorder during obesity［J］.Cell Metab，2022，34（8）：1121-36.e6
-
[18] ZHOU Q，FU Z，GONG Y，et al.Metabolic health status contributes to transcriptome alternation in human visceral adipose tissue during obesity［J］.Obesity（Silver Spring，Md），2020，28（11）：2153-2162
-
[19] STEFAN N，HäRING H ⁃U，HU F B，et al.Metabolically healthy obesity：epidemiology，mechanisms，and clinical implications［J］.Lancet Diabetes Endocrinol，2013，1（2）：152-162
-
[20] NEWMAN A M，LIU C L，GREEN M R，et al.Robust enu⁃ meration of cell subsets from tissue expression profiles ［J］.Nat Methods，2015，12（5）：453-457
-
[21] WEISBERG S P，MCCANN D，DESAI M，et al.Obesity is associated with macrophage accumulation in adipose tis⁃ sue［J］.J Clin Invest，2003，112（12）：1796-808
-
[22] WANG Y，SMITH W，HAO D，et al.M1 and M2 macro⁃ phage polarization and potentially therapeutic naturally occurring compounds［J］.Int Immunopharmacol，2019，70：459-66
-
[23] LI C，SPALLANZANI R G，MATHIS D.Visceral adipose tissue Tregs and the cells that nurture them［J］.Immunol Rev，2020，295（1）：114-125
-
[24] MARTINEZ ⁃ SANCHEZ M E，HIRIART M，ALVAREZ ⁃ BUYLLA E R.The CD4+ T cell regulatory network medi⁃ ates inflammatory responses during acute hyperinsu⁃ linemia：a simulation study［J］.BMC Syst Biol，2017，11（1）：64
-
[25] ZHAO X Y，ZHOU L，CHEN Z，et al.The obesity ⁃ in⁃ duced adipokine sST2 exacerbates adipose T and ILC2 depletion and promotes insulin resistance［J］.Sci Adv，2020，6（20）：eaay6191
-
[26] TAO L，LIU H，GONG Y.Role and mechanism of the Th17/Treg cell balance in the development and progres⁃ sion of insulin resistance［J］.Mol Cell Biochem，2019，459（1⁃2）：183-188
-
[27] CORTESI F，DELFANTI G，CASORATI G，et al.The pathophysiological relevance of the iNKT Cell/Mononucle⁃ ar phagocyte crosstalk in tissues［J］.Front Immunol，2018，9：2375
-
[28] LEE B⁃C，KIM M⁃S，PAE M，et al.Adipose Natural Killer Cells Regulate Adipose Tissue Macrophages to Promote Insulin Resistance in Obesity［J］.Cell Metab，2016，23（4）：685-698
-
[29] LAUMONNIER Y，KARSTEN C M，KöHL J.Novel in⁃ sights into the expression pattern of anaphylatoxin recep⁃ tors in mice and men［J］.Mol Immunol，2017，89：44-58
-
[30] ZAAL A，DIEKER M，OUDENAMPSEN M，et al.Ana⁃ phylatoxin C5a regulates 6 ⁃ Sulfo ⁃ LacNAc dendritic cell function in human through crosstalk with toll⁃Like recep⁃ tor⁃induced CREB signaling［J］.Front Immunol，2017，8：818
-
[31] ZAAL A，VAN HAM S M，TEN BRINKE A.Differential effects of anaphylatoxin C5a on antigen presenting cells，roles for C5aR1 and C5aR2［J］.Immunol Lett，2019，209：45-52
-
[32] PHIELER J，CHUNG K J，CHATZIGEORGIOU A，et al.The complement anaphylatoxin C5a receptor contributes to obese adipose tissue inflammation and insulin resis⁃ tance［J］.J Immunol，2013，191（8）：4367-4374
-
[33] POURSHARIFI P，LAPOINTE M，FISETTE A，et al.C5aR and C5L2 act in concert to balance immunometabo⁃ lism in adipose tissue［J］.Mol Cell Endocrinol，2014，382（1）：325-333
-
[34] ROY C，GUPTA A，FISETTE A，et al.C5a receptor defi⁃ ciency alters energy utilization and fat storage［J］.PLoS One，2013，8（5）：e62531
-
[35] MISSELWITZ B，WYSS A，RASELLI T，et al.The oxyster⁃ ol receptor GPR183 in inflammatory bowel diseases［J］.Br J Pharmacol，2021，178（16）：3140-3156
-
[36] COPPERI F，SCHLEIS I，ROUMAIN M，et al.EBI2 is a negative modulator of brown adipose tissue energy expen⁃ diture in mice and human brown adipocytes［J］.Commun Biol，2022，5（1）：280
-
摘要
目的：探究体重正常（lean）、肥胖代谢正常（metabolically healthy obesity，MHO）与肥胖代谢异常（metabolically un- healthy obesity，MUO）个体内脏脂肪组织（visceral adipose tissue，VAT）中免疫细胞占比差异。方法：内脏脂肪组织标本来源于南京医科大学第一附属医院因肥胖行代谢手术治疗的48例患者及体重和代谢指标正常但因其他非炎症性疾病（疝气、胆结石及肾囊肿等）行腹腔镜手术的11例患者。提取脂肪组织RNA，进行RNA-seq检测，根据RNA-seq结果统计分析3组间内脏脂肪免疫细胞构成比情况和差异基因表达情况，通过单因素方差分析及Pearson相关性分析等分析方法进行数据统计。结果：各组间调节性T细胞，激活的NK细胞及M0、M2型巨噬细胞占比存在差异，P值分别为0.031、0.002、0.004和0.024；与其他两组相比，MUO组M2型细胞占比下降（P=0.014）；两组肥胖患者的内脏脂肪中SPP1、CHIT1、ITGAX、CD300LB、IL1RN及DCSTAMP等基因与M0、M1型巨噬细胞占比呈正相关关系，而与M2型巨噬细胞占比多呈负相关关系。C5a受体1（complement fragment 5a receptor 1，C5AR1）及G蛋白偶联受体183（G-protein coupled receptor 183，GPR183）基因在MHO组中表现为与M0及M1型巨噬细胞占比正相关，与M2型巨噬细胞占比负相关；在MUO组中则呈现相反的关联性。结论：与体重正常组相比较，肥胖组人群的内脏脂肪组织存在显著巨噬细胞浸润；肥胖两组间比较，代谢异常者巨噬细胞亚类中M1型巨噬细胞占比升高，提示机体内高炎症水平与代谢能力受损有关；C5AR1、GPR183可能作为预测肥胖患者代谢状态的生物学指标。
Abstract
Objective：The study aims to investigate the different proportions of immune cells in visceral adipose tissue（VAT）among lean group，metabolically healthy obesity（MHO）group and metabolically unhealthy obesity（MUO）group. Methods：Samples of VAT came from 48 obese patients who received metabolic surgeries and 11 other patients who had non-inflammatory disease（such as hernia， cholecystolithiasis and renal cyst）and received laparoscopic surgery at the First Affiliated Hospital of Nanjing Medical University. After extracting RNA from VAT，we applied RNA-seq technology to count and analyze proportions of immune cells and differentially expressed genes in VAT among three groups. One - way ANOVA and Pearson correlation analysis were used for statistical analysis. Results：The proportion of regulatory T cells，natural killer cells，M0 macrophages，and M2 macrophages varied among three groups（P= 0.031；P=0.002；P=0.004；P=0.024）. In MUO group，the proportion of M2 macrophages decreased compared with other two groups（P= 0.014）. SPP1，CHIT1，ITGAX，CD300LB，IL1RN，DCSTAMP and some other genes had positive correlation with the proportion of M0 and M1 macrophages in two obese groups，but had negative correlation with the proportion of M2 macrophages. C5AR1 and GPR183 had positive correlation with the proportion of M0 and M1 macrophages in MHO group and had negative correlations with the proportion of M2 macrophages. However，the correlation between these two genes and macrophages is opposite in MUO group compared with MHO group. Conclusion：Compared with lean group，the infiltration status of immune cells in VAT is more severe in MHO and MUO group. Comparing between MHO and MUO group，the proportion of M1 macrophages in VAT lifted，indicating that higher level of inflammation in vivo connects with impaired metabolic ability. The differently expressed gene C5AR1 and GPR183 between two obese groups may play a role in predicting alterations in metabolic ability.
Keywords
obesity ； visceral adipose tissue ； RNA-seq ； metabolic status ； immune cell