DOI: https://doi.org/10.21203/rs.3.rs-1344846/v1
Background: Loss of adipose tissue is one of the most important characteristics of cancer cachexia. The loss degree of subcutaneous adipose tissue (SAT) and visceral adipose tissue (VAT) is different in the same cachectic patient. However, it is unclear which genes lead to the differential fat loss in cancer cachexia. The purpose of this study is to screen the differentially expressed genes (DEGs) profile of SAT and VAT and clarify its function in cancer cachexia.
Methods: We screened the DEGs of SAT and VAT in cachectic patients with gastric cancer by RNA sequencing. Gene enrichment analyses were preformed to analyze the biological function of DEGs. Interaction network was constructed to illustrate the relationship among key genes. Clinical significance of IRX1 related to cancer cachexia was illustrated by Spearman’s correlation analysis. SAT and VAT samples were collected and performed histomorphology analyses. Mouse model of cancer cachexia was established and differences between SAT and VAT were compared. The function of IRX1 on lipid metabolism was clarified by Oil Red O staining, qPCR, and western blotting in adipocytes.
Results: A total of 455 DEGs were screened between SAT and VAT in cachectic patients with gastric cancer. Gene enrichment analysis indicated that the immune response might contribute to the pathological difference of SAT and VAT in caner cachexia. Patients with higher expression of IRX1 in SAT than VAT revealed significantly higher weight loss, IL-6 and TNFa concentration, as well as lower BMI, apoE, SAT and VAT area. A greater loss in SAT was observed in cachectic mice comparing to non-cachectic mice while the expression of IRX1 in SAT of cachectic mice was significantly higher. Oil Red O staining showed that IRX1 inhibited lipogenesis in adipocytes. Significant changes in the genes related to adipocyte differentiation and adipogenesis were observed in both the cases of IRX1 overexpression and knockdown.
Conclusion: Higher expression of IRX1 was observed in SAT than VAT in gastric cancer patients with cachexia and thus IRX1 promotes loss of SAT by inhibiting adipocyte differentiation and adipogenesis. IRX1 is expected to be a potential target for the prevention and treatment of SAT loss in cancer cachexia.
Cancer cachexia is a clinical syndrome characterized by muscle atrophy and fat loss[1]. As an important feature of patients with cancer cachexia, fat loss can reduce chemotherapy tolerance and is an independent prognostic factor for poor prognosis[2]. Our previous studies have also shown that fat loss is related to poor nutritional status and high inflammatory factors in patients with cancer cachexia[3, 4]. It is also reported that fat loss occurs earlier than muscle atrophy in cachexic patients[5]. Therefore, the role and mechanism of fat loss in cancer cachexia have attracted much attention.
The anatomical location of the fat depot influences adipose tissue metabolism and function. SAT comprises more than 80% of the total body fat, whereas VAT comprises up to 10% and 20% of the total body fat, respectively[6]. CT image analysis enables precise quantification of adipose depots and has emerged as the golden standard to segment SAT and VAT[7].Not only do SAT and VAT differ in anatomic location, but also in endocrine function, adipokine secretion and lipolytic activity [8, 9]. For example, the accumulation of VATcontributes to obesity related diseases, such as metabolic syndrome[10]. Accumulation of VAT is also a high-risk factor for colorectal cancer and other tumors[11, 12]. On the contrary, SAT accumulation,if not brings a protective effect, is unlikely todo harminmetabolism of patients.These metabolic differences demonstrate the necessity to evaluate the effect based on SAT and VAT respectively rather than total adipose tissue. Our previous research has already found that the loss of SAT was greater than the loss of VAT in a cohort study of 411 gastric cancer patients with cachexia[13].And the SAT area index was an independent risk factor for poor prognosis in patients with gastric cancer cachexia. Hence, it is rather essential to evaluate the different gene profile between SAT and VAT, and clarify the key genes that prompt the greater loss of SAT than VAT. Although genes were reported to be differentially expressed in SAT and VAT in obese patients[14],theDEGsbetween SAT and VAT in gastric cancer patients with cachexia has not been fully explored before.
It has beenconfirmed that the main mechanisms of adipose tissue lossin cancer cachexia include activation of lipolysis, white-to-brown trans-differentiation of WAT (WAT browning), and dysfunctional adipogenesis[15]. Our previous research has proved these points from multiple perspectives. Interleukin-6 (IL-6)was proved to induce adipose tissue loss in patients with cancer cachexia by promoting lipolysis and WAT browning[4]. LncRNA-CAAlnc1 couldsuppress adipogenesis by blocking the binding of HuR in adipose tissue of cancer cachexia patients[16]. MiR-410-3p could down-regulate IRS-1 and PPAR-γ to inhibitadipogenesis and lipid accumulation[17]. CircPTK2 was upregulated in adipose tissues of cachectic patients and induced adipose tissue loss by promoting lipolysis and inhibiting adipogenesis[18]. However, these studies have mainly focused on the DEGs of SAT between cachectic patients and non-cachectic patients. The DEGs betweenSAT and VAT in cachectic patients and the key genes among them remain unclear.
In this study, RNA sequencing was performed to screen DEGsbetween SAT and VAT of gastric cancer patients with cachexia. Clinicopathological significance of DEGs was analyzed in a large number of samples of gastric cancer patients with cachexia. In addition, DEGs were further verified in cancer cachexia mouse models. The function of the DEGswas analyzedby gain and loss of function in preadipocytes.
Patients
Cachectic patients with gastric cancerwho underwentsurgery in Zhongshan Hospital of Fudan University from 2018.06 to 2020.06 were includedin this study. Inclusive criteria for the study were: (1) patients diagnosed with gastric adenocarcinomacancer (not gastric stromal tumor orlymphoma); (2) patients received surgical treatmentwithout preoperative radiotherapy or chemotherapy; (3) patients with abdominal CT examinations and complete clinical data; (4) patients with weight loss >5% in recent 6 months before surgery. This study was approved by the Ethics Committee of Zhongshan Hospital of Fudan University (B2019-193R). Written informed consents were obtained from all patients.
Clinical data collection
Height, weight, gender, and age were extracted from the preoperative medical records. Body mass index (BMI) was calculated as body weight (kg)/height2 (m2).Cachexia-related indicators were extracted from preoperative blood biochemical examinations. Areas of SAT and VAT from CT scans at the third lumbar vertebra were measuredas described before[13].
Human tissue specimens
Adipose tissues were obtained from enrolled patients. At the beginning of the operation, about 500 mg of subcutaneous adipose tissue near the median abdominal incision was obtained. About 500 mg of omental adipose tissue was taken as visceral adipose tissue within 30 minutes after gastric cancer specimens were isolated. The adipose tissue was immediately stored in liquid nitrogen at -80 ℃ or transferred into tissue fixative for further analysis.
Mouse model of cancer cachexia
The mouse model of cancer cachexia refers to the previous methods[4]. In brief, Cachectic mice were induced by subcutaneous injection of colon-26 adenocarcinoma cells into the right flank of the mice. The littermate control mice received PBS injection only. Micewere euthanized at day 21 post-injection and were dissected toharvestinguinal white adipose tissue (SAT) and epididymal white adipose tissue (VAT). Weight of SAT and VAT was recorded at day 0 (control mice) and day 21 (cachectic mice). The proportion of adipose tissues loss was calculated by weight change (weight at day 0subtract weight at day 21) divided to initial weight. Allanimal studies were performed in accordance with the guidelinesprovided by Animal Care Committee of Fudan University.
RNA sequencing
The total RNA was extracted from 3 paired SAT and VAT from gastric cancer patients with cachexia. After quantification and qualification,a total amount of 1 μg RNA per sample was used as input material for the RNAsample preparations. Sequencing libraries were generated using NEBNextUltraTM RNA Library Prep Kit for Illumina (NEB, USA) following manufacturer’s recommendations and index codes were added to attribute sequences to each sample. The clustering of the index-coded samples was performed on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3-cBot-HS (Illumina, USA) according to the manufacturer’s instructions. After cluster generation, the library preparations were sequenced on an Illumina Novaseq platform and 150bp paired-end reads were generated. FeatureCounts v1.5.0-p3 was used to count the reads numbers mapped to each gene. Fragments per kilobase million (FPKM) of each gene was calculated based on the length of the gene and reads count mapped to this gene. The raw sequencing dataset that supported the resultsof this study was deposited in the NCBI GEO database. The data areaccessible through GEO: GSE186466.
RNA sequencing data analysis
Differential expression analysis of twogroupswas performed using theDESeq2 R package. Differentially expressed transcripts between the two groups were identified when |logFoldChange| >0and the p value <0.05. Gene Ontology (GO) enrichment analysis of DEGs wasimplemented by the ClusterProfiler R package, in which gene length biaswascorrected. We also used ClusterProfiler R package to test thestatistical enrichment of differential expression genes in KEGG pathways. GO terms and KEGG pathways with corrected p value <0.05 were consideredsignificantly enriched by differential expressed genes.
Conventional enrichment analysis based on hypergeometric distribution depends on significantly up-regulated or down-regulated genes, and it is easy to omit some genes with insignificant differential expression but important biological significance. Gene set enrichment analysis (GSEA) does not need to specify a clear differential gene threshold. All genes are sorted according to the degree of differential expression in the two groups of samples, and then statistical methods are used to test whether the preset gene set is enriched at the top or low section of the sorting table. GSEA mainly includes three steps: calculation of enrichment score; estimation of the significance level of enrichment score; multiple hypothesis tests.
The PPI network of DEGs was predicted using the Search Tool for the Retrieval of Interacting Genes (STRING) database. The interaction score threshold of 0.4 was set as the cut-off criterion. The PPI network was constructed using Cytoscape. Comprehensive experimentally validated miRNA-gene interaction data were collected from TargetBase. Transcription factor and gene target data derived from the ENCODE ChIP-seq data. Only peak intensity signal <500 and the predicted regulatory potential score <1 is used (using BETA Minus algorithm).
Cell culture and differentiation
The mouse immortalized white preadipocytes werekindly provided by Professor Qiurong Ding from the Shanghai Institute of Nutrition and Health, Chinese Academy of Sciences. Thiscell line has been described previously and used in several studies toassess the effects of different factors on adipose differentiation andfunction[19]. The culture and differentiation methods of preadipocytes cell lines were as previously reported[18].
Plasmid and siRNA construction
The IRX1 expression plasmid pcDNA3.1þ/IRX1 and empty plasmid pcDNA3.1þ were designed and synthesized by GenePharma (Shanghai, China). Small interfering RNAs (siRNAs) targeting IRX1 were also designed by GenePharma. Cell transfection was conducted using Lipofectamine RNAiMAX Transfection Reagent kit (Invitrogen,USA) or Lipofectamine 2000 Transfection Reagent (Invitrogen, USA) according to the manufacturer’s instruction
Hematoxylin-eosin and Immunohistochemical staining
All the samples were transferred totissue fixativeafter being harvested. The protocols were previously described[18]. In brief, histological sections of adipose tissue were stained with hematoxylin–eosin to evaluate morphological changes and the adipocyte cross-sectional area (CSA). A total of 10 randomly selected fields for each sectionwere captured and analyzed to evaluate adipocyte CSA with a computerized imaging software (ImageJ, USA).For immunohistochemistry, the positive cellsin 10 randomly selected fields per section were counted and evaluated by two independent researchers. And the mean number of positive cells per field was calculated.
Oil Red O (ORO) staining
Mature adipocytes were fixed with 4% formaldehyde for 30 min, then washed twice with PBS. They were stained with 0.3% ORO solution and washed three times with distilled water. To assess lipid accumulation, the dye retained in the cells was dissolved in isopropanol and the absorbance of the resulting solution at 520 nm was examined.
RNA isolation and qRT-PCR
Total RNA was isolated from adipose tissues and adipocytes using TRIzol Reagent (Invitrogen, USA) according to the manufacturer’s recommendations.cDNA was synthesized from 1 μg total RNA using FastKing RT Kit (Tiangen, China).Gene expression analysis was performed using Prime-Script RT master mix (Takara, Japan) in StepOnePlusReal-Time system (Applied Biosystems, USA). Expression levels of targeted geneswerenormalized to the expression of GAPDH. qRT-PCR was performedaccording to the manufacturer’s instructions and the relativefold change was calculated by the 2-ΔΔCt method. Primers were designed and synthesized by Sangon Biotech (Shanghai, China) and are listed in Additional file: Table S2. All experiments were repeated at least three times.
Western blot analysis
Preparation of total protein lysates and western blot analysis were performed as previously described[16]. Primary antibodiesagainst IRX1 (Immunoway,YT2412), CEBPα (Cell signaling technology, #2295),AdipoQ(Cell signaling technology, #2789), and FABP4(Cell signaling technology, #2120)were used. Tubulin expression wasused as an endogenous control.
Statistical analyses
Statistical analyses were performed using GraphPad Prism software. Data calculated from independent experiments were presented as the mean ±standard deviation and a student’s t-test was performed to compare the differences between two groups. To analyze the correlation between IRX1 mRNA levels in SAT and the clinicopathological factors in gastric cancer patients with cachexia, we divided 61 patients into two groups according to the IRX1expression in SATcompared to VAT. Comparisons between these two groups were made using thet-test for continuous data and χ2 test for categorical data.p<0.05 was considered statistically significant.
Identification of DEGs between SAT and VAT
To investigate DEGs between SAT and VAT in patients with cachexia, we performed whole transcriptome RNA sequencing of SAT and VAT from 3 patients. Their clinical features and typical CT image at the third lumbar vertebra were shown in Table S1 and Figure S1.Gene expression distributions of each sample were shown in Figure S2.As illustrated by volcano plot in Figure 1A, a total of 455 DEGs were detected in SAT compared to VAT (p < 0.05 and |logFC| > 0), including 342 downregulated mRNAs and 113 upregulated mRNAs. Cluster analysis of DEGs was shown in Figure 1B by heatmap.Thetop 10downregulated DEGs and top 10 upregulated DEGs were shown in Table 1.To better understand the functional implications of 455 DEGs betweenSAT and VAT in cachectic patients, we performed GO analysis and KEGG analysis, and GSEA analysis. The GO terms were shown in Figure 1C. We found that the leukocyte migration, humoral immune response, complement activation, adaptive immune response , and humoral immune response mediated by circulating immunoglobulin were among the top 10 GO terms, which indicated that the immune response might play an essential role in different functions and mechanisms between SAT and VAT.Interestingly, cytokine-cytokine receptor interaction, complement and coagulation cascades, and chemokine signaling pathway weretop 3 of the KEGG pathways(Figure 1D), which also indicated thesignificance of immune response between SAT and VAT.
GSEA analysis and interaction network of DEGs
GSEA analysis was conducted to achieve further insight intothe biologic pathways involved in different adipose tissue distribution of cachectic patients. As shown in Figure 2A-2C, we found that lipid catabolic process, regulation of fatty acid beta-oxidation, and positive regulation of triglyceride metabolic process were enriched in SAT group. On the contrary, we found that negative regulation of chondrocyte differentiation, positive regulation of calcium-mediated signaling, and negative regulation of muscle cell apoptotic process were enriched in VATgroup (Figure 2D-2F).The network of DEGsbetween SAT and VAT wasconstructed by STRING, with a threshold of 0.4. The PPI network with 205 nodes and 720 edges wasdisplayed by Cytoscape-MCODE (Figure 3A). Then we constructed an interaction network linking top 10 DEGs to their relevant miRNAs/translation factors.As shown in Figure 3B-3C,HOXC8, HOXC10, and IRX5 were key genes of the regulating network.
Expression of IRX1 was higher in SAT than VAT in cachectic patients
Among the most differentially expressed genes, we found that IRX1, which encodes a member of the Iroquois homeobox protein family, is highly expressed in SAT compared to VAT.To further clarify the clinical significance of different expression of IRX1 between SAT and VAT, we performed qRT-PCR assay and showed that IRX1 was significantly upregulated in SATcompared with VAT (p< 0.05, Figure 4A).Among the 61 pairs of adipose tissues, we found that IRX1 was higher expressed in SAT than VAT in 37 pairs and lower expressed in SAT than VAT in 24 pairs. Then, we divided them into two groups and analyzed the clinicopathological correlation between groups. The analysis showed that weight loss, IL6, and TNFα were significantly increased in cachectic patients with higherIRX1 expression in SAT than VAT. While BMI, apoE, SAT area and VAT area were significantly decreased in cachectic patients with higherIRX1 expression in SAT than VAT (Table 2).In addition, correlationanalysis was performed between IRX1 expression and adipose tissue area. The results of the Spearman’s correlation analysis showed that the expression level of IRX1in WAT was negatively correlated with WAT area (Figure 4B), and the expression level of IRX1in VAT revealed no correlation with VAT area(Figure 4C). To further validate our hypothesis, HE staining and IHC were performed in WAT and VAT from cachectic patients.The representative images of HE staining revealed that there is no significant difference between adipocytes of SAT and VAT in terms of morphology(Figure 4D). The representative imagesof IHC confirmed that there is a significant increase in IRX1 expression of adipocytes in SAT with respect to VAT (Figure 4E, p<0.05). These results suggested that IRX1is involved inthe lipid metabolism ofthe pathological process in cachexia.
Expression of IRX1 was higher in SAT than VAT in mouse model of cachexia
To further validate our findings in animal models, we construct a mouse modelof cancer cachexiaas described before. As shown in Figure 5A-5B, it became clear that adipocytesmorphologydid not reveal significant difference between different fat depots of the same cachectic mice. Yet the IRX1 expression of adipocytes in SAT is significantly higher than that in VAT (Figure 5C-5D). To calculate the loss ratio of different fat depots, we recordweight of SAT and VAT at day 0 (control mice) and day 21 (cachectic mice).The result showed that loss ratio of SAT was significantly higher than that of VAT in cachectic mice (Figure 5E). Then we compared the mRNA and protein levels of IRX1 in SAT and VAT of cancercachexia mice. The results showed that the mRNA and protein expression levels of IRX1 inSAT were higher than that in VAT in cancer cachexia mice (Figure 5F-5G), which was consistent with the resultsabove.
IRX1 inhibited adipogenesis in adipocytes
To investigate the role of IRX1 in adipocyte metabolism, IRX1 was upregulated by transfecting preadipocyte with IRX1 overexpression vectors, followed by induction and differentiation. The qPCRresult showed that IRX1 expression had increased about ten-fold bydifferentiation day 6 in the IRX1 overexpression group compared tothe control group (Figure 6A).We detected the changes of adipogenesis-related genes involved in the lipid metabolismafterIRX1 overexpression. The results of qPCR and western blottingshowed that AdipoQ, FABP4, and CEBPα were significantly down-regulated when IRX1 was overexpressed (Figure 6B-6C). IRX1 overexpression also resulted in decreasing in lipid dropletaccording to the Oil Red O staining (Figure 6D).To further demonstrate whether IRX1 was required for adipocyte metabolism, adipocytes weretransfected with three independent siRNAs.Effective knockdown was achieved by two of the siRNAs (Figure 6E).The results of qPCR and western blottingshowed that AdipoQ, FABP4, and CEBPαwere also significantly up-regulatedwhen IRX1 was knockdown (Figure 6F-6G). The results of Oil Red O stainingwere consistent with the findings above (Figure 6H).
The role of fat loss in patients with cancer cachexia is becoming a heat topic. SAT and VAT play different roles in the development of cancer cachexia. However, there is no systematic study on the difference of gene expression between SAT and VAT in gastric cancer patients with cachexia. In this study, the mRNA expression profiles of SAT and VAT in patients with cancer cachexia were screened by RNA sequencing for the first time.342 downregulated and 113 upregulated mRNAs were screened compared SAT to VAT in patients with gastric cancer cachexia. IRX1 was one of the most highly expressed mRNAs in SATboth in cachexic patients and mice. It is observed that weight loss, IL6, and TNFα were significantly increased while BMI, apoE, SAT area and VAT area were significantly decreased in cachectic patients with higher IRX1 expression in SAT than VAT.In vitro experiments showed that IRX1 promotesfat loss by inhibiting lipid differentiation and adipogenesis.
The metabolic differences between SAT and VAT were mostly studied in obese patients. The adipose tissues in both fat depotsare similar in their morphologyand lipid storage functions. As shown in present research, adipocytes from SAT an VAT in both cachectic patients and mice revealed evident atrophy, yet no significant difference was observed regarding the morphology of adipocytes. Nonetheless, emerging evidence has suggested that SAT and VAT have different characteristicsand functional roles in metabolic regulation[20, 21]. Thevisceral depot is often considered more deleterious than thesubcutaneous depot. Two independent research demonstrated that VAT hypertrophy is predominantly associated with dyslipidemia whereas SAT hypertrophy is mainly associated with insulin resistance in humans[22]. However, metabolic differences between SAT and VAT in cancer cachexia have not been reported. In this study, we detected hundreds of DEGs between SAT and VAT in gastriccancer patients with cachexia by mRNA sequencing. GO and KEGG analysis showed that immune response might be involved in the difference between SAT and VAT, which was consistent with the differential expression of immune cells in SAT and VAT by single cell sequencing[23]. However, the functions of specific types of immune cells and their mechanisms during the development of cancer cachexiaremain further clarified.
IRX1, a member of the Iroquois homeobox family of transcription factors, played a crucial role inembryonic development[24]. Recently, IRX1 was identified as apotential tumorsuppressor gene in head and neck squamous cellcarcinoma and gastric cancer[25-27]. However,IRX1 also acted as a metastaticoncogene which was activated by hypomethylation[28]. As far as we know, the role ofIRX1 in the progression of fat loss during cancer cachexia has not been reported. In this study, we found IRX1 was significantly upregulated in SAT compared with VAT in both human and mice. IRX1 expression in SAT was found to benegatively correlated with SAT area andwas involved in the lipid metabolism of the pathological process in cachexia. These results indicated IRX1 gene plays an important role in promoting the loss of SAT during cancer cachexia.
We firstly reported that overexpression of IRX1 inhibited adipogenesisin preadipocytes. Activation of lipolysis and WAT browning or dysfunctional adipose differentiation and adipogenesis were reported to be the main mechanism of fat loss in cancer cachexia[29]. By overexpression and knockdown of IRX1, we clarifiedthat IRX1 exerted its function by inhibiting adipogenesis in the lipid metabolism. These results suggested that adipose differentiation and adipogenesis related genes might be important in the regulation of lipid metabolism during cancer cachexia, which was consistent with our previous foundation.In future research, we will further explore the specific molecular mechanism of IRX1 regulating adipose differentiation and adipogenesis related genes.
Inconclusion, we screened mRNA expression profiles of SAT and VAT in patients with cancer cachexia by RNA sequencing. Hundreds of DEGs were detected inSAT and VAT and IRX1 wasone of thehighly expressed mRNAs in SAT. Clinical studies showed thatIRX1 expression in SAT was negatively correlated with SAT area.In vitro experiments showed that IRX1 inhibited lipogenesis by modulating differentiation and adipogenesis related genes. Our results demonstrated that IRX1 was expected to be a potential target for the prevention and treatment of SAT loss in cancer cachexia.
WAT: White adipose tissue; SAT: Subcutaneous adipose tissue; VAT: Visceral adipose tissue; FPKM: Fragments per kilobase million; GO: Gene Ontology; KEGG: Kyoto Encyclopedia and Genes and Genomes; GSEA: Gene set enrichment analysis; ORO: Oil Red O; BMI: Body mass index; lncRNA: Long non-coding RNA; circRNA: Circular RNA; siRNA: Small interfering RNAs; ceRNA: Competing endogenouse RNA; IL-6: Interleukin- 6; TNF-α: Tumor necrosis factor-α; FABP4: Fatty acid binding protein 4; CEBPα: CCAAT/enhancer binding proteins α; IRX1: Iroquois homeobox 1.
Acknowledgements
We thank Prof. Qiurong Ding from Shanghai Institute of Nutrition and Health, Shanghai Institutes for Biological Sciences,Chinese Academy of Sciences for kindly providing the immortalized subcutaneous white preadipocytes.
Authors’ Contribution
J.H. and G.W. conceived, designed, and drafted the manuscript. Q.Z., L.S., F.Y. and S.S. collected the clinical data. J.H and Z.D performed the experiments. J.H., Z.D. and Q.Z. conducted the statistical analyses. J.H. drafted the manuscript. G.W. and J.H. revised the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by the Youth of National Natural Science Foundation of China (81803091).
Availability of Data and Materials
RNA sequencing data were deposited at the GEO database with the accession number GSE1186466.All data generated or analyzed during this study are included in this published article and its supplementary information files.
Ethics Approval and Consent to Participate
This study was conducted according to the standards of the Declaration of Helsinki. Approval by the Ethics Committee of Zhongshan Hospital, Fudan University was obtained. Written informed consents were obtained from all participants. All procedures described involving animals were approved by the Animal Care Committee of Fudan University and conducted in accordance with the guidelines of Fudan University on the experimental use of laboratory animals.
Competing Interests
The authors declare that they have no competing interests.
Consent for Publication
Not applicable.
Table 1. Top 10 downregulated mRNAs and 10 upregulated mRNAs in SAT compared to VAT
Gene id |
Gene name |
log2 Fold Change |
p value |
Downregulated mRNAs |
|||
ENSG00000184937 |
WT1 |
-10.42020286 |
7.24E-09 |
ENSG00000119919 |
NKX2-3 |
-9.340609993 |
3.06E-12 |
ENSG00000183242 |
WT1-AS |
-8.433724032 |
1.59E-10 |
ENSG00000211892 |
IGHG4 |
-8.083413442 |
1.00E-09 |
ENSG00000211668 |
IGLV2-11 |
-7.842925465 |
4.15E-08 |
ENSG00000054803 |
CBLN4 |
-7.209658278 |
2.30E-06 |
ENSG00000205038 |
PKHD1L1 |
-7.107755456 |
5.64E-06 |
ENSG00000187848 |
P2RX2 |
-7.045386839 |
2.30E-06 |
ENSG00000242221 |
PSG2 |
-6.846739335 |
6.57E-06 |
ENSG00000211663 |
IGLV3-19 |
-6.843669107 |
9.34E-06 |
Upregulated mRNAs |
|||
ENSG00000112246 |
SIM1 |
10.22187292 |
1.84E-15 |
ENSG00000176842 |
IRX5 |
7.352859596 |
2.08E-08 |
ENSG00000089225 |
TBX5 |
6.447165153 |
2.06E-05 |
ENSG00000170549 |
IRX1 |
6.369599248 |
2.82E-22 |
ENSG00000180318 |
ALX1 |
6.250090981 |
0.000155648 |
ENSG00000180818 |
HOXC10 |
6.123438559 |
7.65E-24 |
ENSG00000251151 |
HOXC-AS3 |
5.976364812 |
0.001044655 |
ENSG00000154162 |
CDH12 |
5.92279196 |
0.00024887 |
ENSG00000255399 |
TBX5-AS1 |
5.831484757 |
3.40E-06 |
ENSG00000250451 |
HOXC-AS1 |
5.825341418 |
3.79E-22 |
Table 2. ClinicalCharacteristics of 61 Participants
Clinical characteristics |
IRX1 |
p Value |
|
High (n=37) |
Low (n=24) |
||
Age |
64.71±9.06 |
65.79±8.38 |
0.633 |
BMI |
20.36±2.12 |
22.82±1.61 |
0.006* |
Weight Loss |
6.53±2.19 |
5.12±1.93 |
0.015* |
Disease stage(III/IV) |
23/37 |
14/24 |
0.772 |
IL6(mmol/L) |
10.67±5.87 |
7.52±4.38 |
0.027* |
TNFa(mmol/L) |
14.24±6.36 |
11.80±4.11 |
0.019* |
Alb (g/L) |
37.89±3.43 |
39.04±4.45 |
0.261 |
PAb(mg/L) |
207.59±52.98 |
220.17±58.42 |
0.093 |
FAA(mmol/L) |
0.58±0.13 |
0.53±0.17 |
0.098 |
TC(mmol/L) |
4.01±0.95 |
4.26±1.05 |
0.351 |
TG(mmol/L) |
1.07±0.41 |
1.10±0.37 |
0.829 |
LDL(mmol/L) |
2.37±0.84 |
2.60±0.91 |
0.317 |
HDL(mmol/L) |
1.18±0.36 |
1.25±0.57 |
0.599 |
ApoA(g/L) |
1.16±0.28 |
1.13±0.26 |
0.445 |
ApoB(g/L) |
0.79±0.25 |
0.85±0.21 |
0.268 |
ApoE(mg/L) |
35.78±14.30 |
44.13±17.70 |
0.047* |
SAT area(cm2) |
108.54±19.18 |
127.09±24.13 |
<0.001* |
VAT area(cm2) |
87.83±17.78 |
101.09±19.85 |
0.015* |
BMI Body mass index, ALB Albumin, PAb Prealbumin, TC Totalcholesterol, TG Tri-glyceride, LDL Low-density lipoprotein, HDL High-density lipoprotein, ApoA Apolipoprotein A, ApoB Apolipoprotein B, ApoE Apolipoprotein E, FFA Free fatty acid, IL-6 Interleukin 6, TNF-α Tumor Necrosis Factor-α