Next Article in Journal
Molecular Signature of Neuroinflammation Induced in Cytokine-Stimulated Human Cortical Spheroids
Next Article in Special Issue
Mitochondrial Phenotype as a Driver of the Racial Dichotomy in Obesity and Insulin Resistance
Previous Article in Journal
Targeting Src-Hic-5 Signal Cascade for Preventing Migration of Cholangiocarcinoma Cell HuCCT1
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomedicines
Volume 10
Issue 5
10.3390/biomedicines10051024
biomedicines-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Morbid Obesity in Women Is Associated with an Altered Intestinal Expression of Genes Related to Cancer Risk and Immune, Defensive, and Antimicrobial Response

by
Ailec Ho-Plágaro
1,2,3,
Cristina Rodríguez-Díaz
1,2,
Concepción Santiago-Fernández
1,2,
Carlos López-Gómez
1,2,
Sara García-Serrano
1,4,5,
Flores Martín-Reyes
1,2,
Francisca Rodríguez-Pacheco
1,2,5,
Alberto Rodríguez-Cañete
6,
Guillermo Alcaín-Martínez
1,2,
Luis Vázquez-Pedreño
7,
Sergio Valdés
1,4,5,*,
Lourdes Garrido-Sánchez
1,8,† and
Eduardo García-Fuentes
1,2,9,*,†
1
Instituto de Investigación Biomédica de Málaga-IBIMA, 29010 Malaga, Spain
2
Unidad de Gestión Clínica de Aparato Digestivo, Hospital Universitario Virgen de la Victoria, 29010 Malaga, Spain
3
Departamento de Biología Celular, Genética y Fisiología, Universidad de Málaga, 29010 Malaga, Spain
4
Unidad de Gestión Clínica de Endocrinología y Nutrición, Hospital Regional Universitario de Málaga, 29009 Malaga, Spain
5
CIBER de Diabetes y Enfermedades Metabólicas Asociadas (CIBERDEM), 29010 Malaga, Spain
6
Unidad de Gestión Clínica de Cirugía General, Digestiva y Trasplantes, Hospital Regional Universitario, 29010 Malaga, Spain
7
Unidad de Gestión Clínica de Aparato Digestivo, Hospital Regional Universitario, 29010 Malaga, Spain
8
Unidad de Gestión Clínica de Endocrinología y Nutrición, Hospital Universitario Virgen de la Victoria, 29010 Malaga, Spain
9
CIBER de Enfermedades Hepáticas y Digestivas (CIBEREHD), 29010 Malaga, Spain
*
Authors to whom correspondence should be addressed.
†
These authors contribute equally to this work.
Biomedicines 2022, 10(5), 1024; https://doi.org/10.3390/biomedicines10051024
Submission received: 6 April 2022 / Revised: 22 April 2022 / Accepted: 26 April 2022 / Published: 29 April 2022
(This article belongs to the Special Issue Emerging Paradigms in Insulin Resistance 2.0)
Browse Figures
Versions Notes

Abstract

:
Background: Little is known about the relation between morbid obesity and duodenal transcriptomic changes. We aimed to identify intestinal genes that may be associated with the development of obesity regardless of the degree of insulin resistance (IR) of patients. Material and Methods: Duodenal samples were assessed by microarray in three groups of women: non-obese women and women with morbid obesity with low and high IR. Results: We identified differentially expressed genes (DEGs) associated with morbid obesity, regardless of IR degree, related to digestion and lipid metabolism, defense response and inflammatory processes, maintenance of the gastrointestinal epithelium, wound healing and homeostasis, and the development of gastrointestinal cancer. However, other DEGs depended on the IR degree. We mainly found an upregulation of genes involved in the response to external organisms, hypoxia, and wound healing functions in women with morbid obesity and low IR. Conclusions: Regardless of the degree of IR, morbid obesity is associated with an altered expression of genes related to intestinal defenses, antimicrobial and immune responses, and gastrointestinal cancer. Our data also suggest a deficient duodenal immune and antimicrobial response in women with high IR.

1. Introduction

Obesity has reached epidemic proportions in recent decades. It is associated with numerous disorders such as a low-grade systemic inflammation, insulin resistance, dyslipidemia, and hypertension, which in turn increase the risk of heart disease, stroke, and type 2 diabetes mellitus (T2DM) [1,2]. Overfeeding, changes in the composition of the diet (higher caloric density), and a sedentary lifestyle are some of the main underlying causes for overweight and obesity [3]. The gastrointestinal (GI) tract is the first point of interaction between the host, the microbiota, and antigens coming from the diet; therefore, unhealthy food consumption habits disturb, in the first instance, the composition of the intestinal microbiota and the homeostasis of the GI tract. In this sense, it is known that obesity is associated with dysbiosis, altered intestinal mobility, permeability, and inflammation, as well as gut–brain disorders [4,5]. These imbalances can not only modulate the metabolism, promoting and aggravating the development of obesity, but can also be the starting point that triggers a series of events involved in the development of comorbidities associated with obesity, such as insulin resistance and T2DM [6].
Most studies on obesity and T2DM are focused on the study of metabolism modulating tissues such as adipose, skeletal muscle, or liver [7]. The remission of T2DM, as well as the improvement of other pathologies in patients with morbid obesity undergoing bariatric surgery, has revealed the important involvement of the small intestine in metabolism and body homeostasis [8]. However, the GI tract is not only involved in digestive functions. The intestine, through the secretion of specific molecules, and the immune system modulate the composition of the microflora, which is considered an important virtual endocrine organ [9]. On the other hand, the GI has a direct role in the regulation of the body’s energy balance through the coordinated action of hormonal, neural, and immunological signals [10]. Although studies have been conducted on the matter, little is known about the intestinal transcriptomic changes in patients with morbid obesity and how these changes may be different depending on the degree of insulin resistance. To date, there are few studies in this regard, with the majority being carried out in animal models [11,12,13,14]. In previous studies we have observed significant changes in the protein expression of various cytokines and chemokines related to pro/anti-inflammatory processes in the duodenum of patients with morbid obesity and different grades of insulin resistance, as well as in intestinal permeability [15]. Moreover, we have observed changes in the duodenal transcriptome of women according to their insulin resistance level, independently of the patient’s body mass index (BMI) [16].
To provide more information, in this study we aimed to identify intestinal genes that may be associated with the development of obesity regardless of the degree of insulin resistance of patients. To this effect, we identified genes with the same expression profile in the duodenum of women with morbid obesity alongside both high and low insulin resistance, which are the same genes that are associated with BMI regardless of insulin resistance.

2. Material and Methods

2.1. Subjects

The study included 15 patients with morbid obesity (MO) (BMI > 40 kg/m2) and 6 healthy, non-obese (NO) women (BMI < 30 kg/m2) (Table 1) (see Supplementary Table S1 for more detailed information about the women). The MO women included in this study underwent sleeve gastrectomy (SG) at the Regional University Hospital of Malaga (Spain). One month before SG, they underwent programmed gastroscopy to discard alterations in the stomach, at which time duodenal biopsies were obtained. Subjects were excluded if they had T2DM with anti-diabetic oral or insulin treatment, acute inflammatory disease, or infectious disease or if the patient did not consent. The MO women were classified into two groups according to the homeostasis model assessment of IR (HOMA-IR) level (low HOMA-IR value (<4.7) (MO-low-IR, n = 7) or high HOMA-IR value (>4.7) (MO-high-IR, n = 8)) (both groups without treatment for T2DM) [15,16]. The cut-off point for the HOMA-IR was taken from previous studies carried out in our population of patients with morbid obesity [15,16,17]. The non-obese women were selected among those who underwent programmed gastroscopy at the Virgen de la Victoria University Hospital (Málaga, Spain) with non-pathological results, at which time duodenal biopsies were taken. They had a similar average age to the MO group and reported that their body weight had been stable for at least 3 months prior to the study. We included only one group of NO patients who had low insulin resistance (NO-low-IR). This cut-off point for the HOMA-IR was obtained from the 75th percentile of the HOMA-IR value for non-obese subjects in our area with normal glucose metabolism according to the 1998 American Diabetes Association classification [18]. These subjects (NO and MO women) had been included in previous studies [15,16]. All individuals included were of Caucasian origin. The samples were processed and frozen immediately after arrival at the Regional University Hospital Biobank (Andalusian Public Health System Biobank). Participants gave their written informed consent. The study was carried out in accordance with the Code of Ethics of the World Medical Association (Declaration of Helsinki) and approved by the Malaga Provincial Research Ethics Committee, Spain (PI12/00338).

2.2. Analytical Procedures

Serum glucose, cholesterol, and triglycerides (Randox Laboratories Ltd., Antrium, UK) were measured in fasting state by standard enzymatic methods. Insulin was analyzed using an immunoradiometric kit (DIAsource ImmunoAssays SA, Louvain-la-Neuve, Belgium). The homeostasis model assessment of insulin resistance (HOMA-IR) was calculated with the following equation: HOMA-IR = fasting insulin (µIU/mL) × fasting glucose (mmol/L)/22.5 [17].

2.3. Duodenal Samples

Duodenal biopsies from NO and MO women from the 2nd part of the duodenum were obtained in a fasting state during a gastroscopy that the patients underwent [15,16]. The mucosa was washed with physiological saline solution, immediately frozen in liquid nitrogen, and maintained at −80 °C until analysis. Total RNA isolation from frozen duodenal biopsy samples was performed using an RNeasy Mini Kit (Qiagen GmbH, Hilden, Germany).

2.4. Microarray Procedure

The study of the differential expression profiling, with a SurePrint-G3 Human GE 8 × 60 K microarray kit (ID 028004; GPL13607) (Agilent Technologies, Madrid, Spain), was carried out following an experimental design in which a two-pairwise comparison was performed. Samples for each experimental condition were labeled, hybridized, washed, and scanned according to the two-color microarray-based gene expression analysis v6.5 protocol of the Genetic Diagnostic Bioarray facilities (Bioarray, Alicante, Spain). This microarray contained 62,976 probes, of which 58,717 had no controls. These probes correspond to 21,414 genes (filtered by gene identification). From these 58,717 probes, 45,283 showed mRNA expression. The data discussed in this publication have been deposited in the NCBI’s Gene Expression Omnibus (GEO) (GEO Series access number GSE147562).

2.5. Microarray Data Analysis

Data analysis was performed with Agilent Feature Extraction Software v.10.7 (Agilent Technologies), using the latest gene annotations available. The expression of each gene was reported as the ratio of the value obtained after each condition relative to the control condition after normalization of the data against the median of the control samples. A filter was applied to select the genes with significant differential expression (DEGs) that displayed an adjusted FDR of less than 0.05 by a nonparametric analysis (Rank Product). This analysis has the advantage of being less sensitive to the variability of the samples; thus, it is suitable in systems where high variability is expected. Prior to any further analysis, the array points were filtered in order to discard replicated genes, filtering first by probe name and later by systematic name. Due to the direction of the hybridizations, the DEGs in this study were those that were up- or downregulated in the MO groups when compared with the NO group. We obtained two lists of DEGs (FDR < 0.05) between (a) the MO-low-IR group with respect to the NO-low-IR group, and between (b) the MO-high-IR group with respect to the NO-low-IR group. These two gene lists were overlapped by a Venn diagram obtaining the upregulated and downregulated DEGs shared by MO-low-IR and MO-high-IR women with respect to the NO-low-IR group (genes in intersection), as well as the exclusive DEGs of MO-low-IR and MO-high-IR.

2.6. Functional Enrichment

The lists of up- and downregulated genes offered by the Venn diagram were functionally analyzed using Gene Set Enrichment Analysis (GSEA; GSEA/MSigDB web site v6.4 version; MSigDB database v7.2 updated September 2020; https://www.gsea-msigdb.org/gsea/index.jsp, accessed on 27 March 2021) [19,20], which is a computational method that determines whether an a priori-defined set of genes shows statistically significant, concordant differences between two biological states. We computed overlaps between our gene sets and gene sets in the MSigDB database. Using this platform, we identified gene sets significantly overrepresented (FDR q-value < 0.05) in different selected collections and sub-collections. The FDR q-value is a false discovery rate analog of hypergeometric p-value after correction for multiple hypothesis testing according to Benjamini and Hochberg. Overrepresentation analysis is a technique for determining whether a set of terms is present more than it would be expected. The collections and sub-collections used in the analysis were the following: gene ontology gene sets (GO terms), chemical and genetic perturbations (CGP), hallmark (H), human phenotype ontology (HPO), KEGG subset of canonical pathways (CP:KEGG), REACTOME subset of canonical pathways (CP:REACTOME), and cancer modules (CM). GO defines the function used to describe gene function with respect to three aspects: molecular function (MF) (molecular activities of gene products), cellular component (CC) (where gene products are active), and biological process (BP) (pathways and larger processes made up of the activities of multiple gene products).

2.7. Technical Validation of Microarray Data by Real-Time–Quantitative PCR (RT–qPCR)

The technical validation of the data obtained in the microarray was carried out in our previous study performed in these same patients [16].

2.8. Statistical Methods

We designed the experiment taking into account the anticipated number of undifferentially expressed genes in the microarray (20,000), 1 as the number of false positives, 0.80 as the desired power, 2 as the mean difference in log-expression between two groups, and 1 as the anticipated standard deviation of the difference in log-expression between two groups. With this design, the sample size for each group is 6 [21,22]. The statistical analysis was performed with R statistical software, version 2.8.1 (Department of Statistics, University of Auckland, Auckland, NZ; http://www.r-project.org/). Differences between groups were established using the Mann–Whitney test. Values were considered to be statistically significant when p ≤ 0.05. The results are given as the median (interquartile range). The statistical significance of the microarray gene expression is described above.

3. Results

3.1. Differentially Expressed Genes (DEGs) in the Microarray Hybridizations

We obtained two lists of DEGs (FDR < 0.05): between the MO-low-IR and NO-low-IR groups and between the MO-high-IR and NO-low-IR groups. Excluding the unnamed genes, we found a total of 175 DEGs among the MO-low-IR and NO-low-IR groups (80 upregulated and 95 downregulated) and 138 DEGs between the MO-high-IR and NO-low-IR groups (39 upregulated and 99 downregulated). In the intersection of these two groups of DEGs, we found 73 DEGs that were associated with BMI, regardless of the degree of insulin resistance (Figure 1). It is worth noting the high degree of similarity in the up/downregulation patterns of these shared genes, where we found 20 upregulated and 48 downregulated genes in both groups of obese patients. We only found five genes shared between both groups (MUCL3, PGC, TCN1, TFF2, BPIFB1) that were upregulated in MO-low-IR but downregulated in MO-high-IR when compared with women without obesity (Figure 1).
However, other DEGs were differentially expressed between MO and non-obese groups that were dependent on insulin resistance degree; there were 102 DEGs exclusively in the MO-low-IR group and 65 DEGs exclusively in the MO-high-IR group with respect to the NO-low-IR group. In this case, we observed a greater presence of exclusive upregulated DEGs in MO-low-IR and a predominance of downregulated genes in MO-high-IR (Figure 1).

3.2. Functional Enrichment

With the aim of delving into the processes, routes, or biological aspects related to the set of differentially expressed genes, we selected a wide range of gene sets collected in The Molecular Signatures Database (MSigDB) for functional analysis with GSEA. The complete analysis of the overrepresented gene sets (FDR < 0.05) and the associated p-value and FDR q-value are represented in Supplementary Table S2. To simplify the information, we have grouped the gene sets into the most relevant biological functions, processes, or responses in which gene expression is upregulated and downregulated in the duodenum of women with morbid obesity (Table 2); selected gene sets for each function are indicated in Supplementary Table S2.
Our results indicated an alteration in important processes such as inflammation, immune response, hypoxia, digestion, maintenance of the gastrointestinal epithelium, wound healing, homeostasis, and cell proliferation; as well as in the response to biotic stimuli, drugs, xenobiotics, and toxic substances, among others (Table 2). On the other hand, we found the dysregulation of a group of genes that are associated with the development of gastrointestinal cancer. A more detailed analysis of the most relevant results obtained on the processes and biological functions is shown in the following sections.

3.2.1. Digestion and Lipid Metabolism

In the intersection of the two groups of DEGs, the overrepresented gene set (FDR < 0.05) showed an altered expression of genes in the duodenum that are related to digestion and lipid metabolism (Table 2 and Figure 2). There was a group of genes up- or downregulated in both MO groups regardless of the degree of insulin resistance. Lipases such as LPL, LIPF, and PNLIPRP2 were downregulated, while the inhibitors of lipase activity ANGPLT4 and APOC3 were upregulated in the duodenum of both MO groups. In addition, the mucin MUC6 was downregulated in both MO groups. Other genes, such as TFF2 and PGC were upregulated in MO-low-IR but downregulated in MO-high-IR (Table 2 and Figure 2).
However, there were genes involved in digestion that were expressed exclusively in each group of patients (Figure 2 and Table 2). For example, SCT and RBP4 were upregulated and LCT was downregulated in the MO-low-IR group.

3.2.2. Defense Response and Inflammatory Processes

We found relevant results in the intersection of the two groups of DEGs and through the analysis of the genes involved in defensive functions and inflammation (Figure 3 and Table 2). According to the analysis with GSEA, many of these genes were related to gene sets of innate and adaptive immune responses (Supplementary Table S2).
We found shared upregulated and downregulated genes involved in both types of defensive responses regardless of the degree of insulin resistance. However, our results point to a trend towards increased expression of genes involved in adaptive and innate immune responses in MO-low-IR and a decrease in MO-high-IR (Table 2 and Figure 3).
Interestingly, according to the overrepresented gene set (FDR < 0.05) related to inflammatory processes and bibliographic analysis, both pro- and anti-inflammatory processes were altered in MO (Figure 4). Morbid obesity, regardless of the degree of insulin resistance, was associated with changes in the expression of genes involved in pro- (such as CD86 and S100A8) and anti-inflammatory (such as CCL18 and ZBTB16) responses (Figure 4). These genes are involved in different functions (Table 3). However, we found a different gene expression pattern depending on the insulin resistance level, mainly in pro-inflammatory responses. In the MO-low-IR group, there was a significant upregulation of genes associated with pro-inflammatory immune processes (such as XCL1, ENPP3, NOS2, CLC, CEBPE, RBP4, MMP9) (Figure 4 and Table 3). Additionally, we found a downregulation of certain genes related with pro-inflammatory processes (CXCL9, SERPINA3) in the MO-high-IR group (Figure 4).
On the other hand, there was an overrepresentation (FDR < 0.05) of a gene set related to antimicrobial processes and responses to another organism (Supplementary Table S2). Our results indicate a defective response to external organisms in MO-high-IR, while in MO-low-IR it appears to be increased (Figure 5).

3.2.3. Hypoxia

In the MO-low-IR group, there was an upregulation of genes involved in overrepresented (FDR < 0.05) gene sets related to hypoxia, although the hypoxia process does not seem to play a special role in MO-high-IR and in the intersection between the MO-low-IR and MO-high-IR groups (Table 2).

3.2.4. Epithelial Maintenance, Wound Healing, and Homeostasis

In the intersection of the two groups of DEGs, we found changes in the expression of genes (FDR < 0.05) related to the maintenance of the gastrointestinal epithelium and epithelial structure (Table 2 and Supplementary Table S2). We found an upregulation of several genes in both MO groups regardless of the degree of insulin resistance (Table 2).
However, the functional analysis revealed mainly an upregulation of genes involved in wound healing functions and tissue homeostasis in the MO-low-IR group but not in the MO-high-IR group (Table 2).

3.2.5. Cancer

In the analysis with GSEA, we found a significant overrepresentation (FDR < 0.05) of gene sets involved in the development of colorectal adenoma and gastric cancer in both MO groups (Table 2 and Supplementary Table S2). As shown in Supplementary Table S3, the changes in gene expression found in this study seem to be associated with the promotion of the migration, invasion, or proliferation of cancer cells.

4. Discussion

In this study, we analyzed changes in gene expression in the duodenum of two groups of women with morbid obesity, one with high and one with low insulin resistance, and compared them with a group of metabolically healthy women. First, we found a group of genes with a similar expression profile in both groups of women with morbid obesity compared to the control group; that is, genes upregulated and downregulated in morbid obesity regardless of the degree of insulin resistance. It is known that obesity is associated with a higher degree of systemic and local inflammation, where hypoxia and the release of cytokines and pro-inflammatory hormones are some of the factors that promote this process [39]. Thus, we have observed changes both in the expression of certain pro-inflammatory and anti-inflammatory genes, such as S100A8, GREM1, GREM2, and CCL18. S100A8 stimulates and promotes the migration of neutrophils and monocytes [40], with a critical role in intestinal pro-inflammatory responses [24,41]. In the same line, the increase in GREM1 and GREM2 expression could be associated with a hyperplasia of the duodenal crypts, which is characteristic of epithelial damage [42]. Nevertheless, we have also found an increase in the expression of the relevant anti-inflammatory chemokine CCL18, which is inducible by Th2 cytokines [43] and involved in M2 macrophage maturation [23], and a subexpression of CD86, which is associated with a pro-inflammatory state [44,45]. Thus, these results suggest the activation not only of pro-inflammatory responses in the duodenum in morbid obesity but also of synergistic anti-inflammatory mechanisms that could negatively regulate the possible collateral damage that occurs during inflammation [46].
We observed other DEGs that were related to lipid metabolism and digestion irrespective of the degree of insulin resistance. The downregulation of LPL, accompanied by the upregulation of two potent inhibitors of LPL enzyme, ANGPTL4 and APOC3 [47,48], was found. Increasing evidence suggests that LPL is regulated in a tissue-specific manner [49]. However, little is known at the intestinal level. In line with these results, we observed the downregulation of other lipases such as LIPF and PNLIPRP2. Overall, our data suggest deficient metabolism in the lipid degradation at the intestinal level in women with morbid obesity, which could constitute a protection mechanism in response to excessive fat intake but at the same time may contribute to the development of dyslipidemia, which could have a relevant role in the etiology of insulin resistance [50].
In morbid obesity, the microbiota could play a role as a regulator of metabolism [51]. Regardless of insulin resistance, the two groups of women with morbid obesity showed upregulation of LYPD8 and downregulation of a greater number of genes, such as REG1B, GKN2, LPL, PNLIPRP2, FKBP5, and CD86, that are related to responses to bacterial stimuli and antimicrobial activity. The microbiota could be involved in the change of CD86 expression since lipopolysaccharides (LPS) can alter its expression in both human peritoneal and M2 macrophages [52]. LYPD8 plays an important role in inhibiting the attachment of flagellated microbiota to colonic epithelia [53] and REG1B has antimicrobial effects [54]. On the other hand, the downregulation of LPL and PNLIPRP2 is not only involved in the regulation of lipid metabolism, but both also have antimicrobial effects. PNLIPRP2 displays its highest phospholipase activities on phosphatidylglycerol and phosphatidylethanolamine, two major constituents of bacterial membranes [55] with a supposed protective role against pathogenic microbiota [56,57]. In this sense, dietary lipids can exert pro- or anti-inflammatory functions on cells of the innate immune system and influence antigen presentation in cells of the adaptive immune system [58]. These results suggest a greater vulnerability of the intestine against microorganisms in patients with morbid obesity.
Nonetheless, our study has tried to delve deeper into the effects of high insulin resistance on the transcriptomic profile of the duodenum. Our results agree with a previous study where we observed an increase in the number of cytokines and chemokines in the duodenum of morbidly obese patients with low insulin resistance [15]. In the current study, we found a large number of upregulated genes involved in inflammatory processes related to the presence of various types of immune cells, such as M1 macrophages, T lymphocytes, NKT, eosinophils, basophils, and neutrophils, in women with morbid obesity and low insulin resistance. There was an increase in RBP4, which is expressed during the differentiation of monocytes into primary macrophages [59] and plays an important role in the development of insulin resistance in adipose tissue [60]. Moreover, an increase in NOS2 expression, which is induced in inflammatory processes by LPS [61] and hypoxia [62], was found. Its expression is characteristic of classically activated M1 macrophages in adipose tissue [31]. At the same time, we observed an increase in XCL1 expression, which is produced by T, NK, and NKT cells during infectious and inflammatory responses [63,64] and has bactericidal activity [65]. The expression of other genes is also increased, such as BTNL8, CLC, and PHD3 (Egln3). BTNL8 is highly expressed by neutrophils [34], which supports the idea that it has a regulatory role in inflammation [66]; along the same lines, other studies have shown the involvement of PHD3 in the hypoxic regulation of neutrophilic inflammation in humans and mice [67]. Moreover, the inhibition of PHD3 improves insulin sensitivity and ameliorates diabetes by specifically stabilizing HIF-2α [68]. On the other hand, CLC is typically expressed in eosinophils, and its increase is typically associated with altered epithelial barrier functions, including food-allergic enteropathies and inflammatory bowel diseases [69]. Although we found a large number of genes involved in pro-inflammatory responses, we have also observed the upregulation of S100G, a molecule that has been linked with anti-inflammatory effects [70]. However, the epithelial barrier function could be damaged by the observed upregulation of CEACAM6 and downregulation of LGALS2 in the MO-low-IR group. The increase in CEACAM6 in the intestinal epithelium has been linked with the massive colonization of adherent invasive Escherichia coli [71]. Additionally, LGALS2 confers maximum protection against exposure to pathogens [72].
It is known that obesity is associated with a greater degree of hypoxia in adipose and intestinal tissue, resulting in adverse metabolic effects like insulin resistance [73]. Hypoxia could contribute significantly to the changes found in the expression of genes in the duodenum of women with morbid obesity. In this sense, we have found an upregulation of different genes related to hypoxia regardless of the degree of insulin resistance, such as ANGPTL4, LPL, S100A8, SLC11A2, HBB, HBA2, and HBD. An increased expression and secretion of ANGPTL4 under hypoxic conditions has been observed in human adipocytes [74], and chronic intermittent hypoxia inhibits LPL by upregulating ANGPTL4 [75]. The activation of hypoxia signaling could induce the expression of the transporter SLC11A2, thus increasing the uptake of iron [76]. Additionally, an increase in the expression of S100A8 mediated by hypoxia has been observed in prostate cancer [77]. In addition, the upregulation of the three hemoglobins HBB, HBA2 and HBD may be indicative of a hypoxia situation [78]. However, there is another group of upregulated genes related to hypoxia in women with morbid obesity and low insulin resistance. This is in contrast to the literature that shows an association between hypoxia and insulin resistance [79]. However, part of these genes, in addition to participating in the regulation of hypoxia, could be involved in other different pathways, playing a beneficial role and supporting the maintenance of immune responses or tissue integrity and thereby influencing tissue recovery. For example, MMP7, which is induced by microbial products, can regulate tissue repair and plays an important role in the maintenance of innate immunity in the intestine, where it activates anti-bacterial peptides such as pro-defensins [80]. On the other hand, Egln3 (PHD3) not only functions as a tumor suppressor but may also promote fibrosis and anti-inflammatory responses and prevent neutrophil apoptosis under hypoxic conditions [81,82]. In addition, NOS2 produces NO, which is involved in the immune response as a defense mechanism; NO is responsible for inhibiting the production of IL-12 and macrophages. The possible functions of the rest of the genes related to hypoxic conditions are more unknown.
The downregulated genes found in the MO-high-IR group are mainly associated with defensive responses against bacteria (SPINK5, CXCL9, GBP3, SLPI, and MUC1). We have found a relation between both changes in microbiota and duodenal immune response in the presence of high insulin resistance in morbid obesity in our previous studies [15,83]. Moreover, we found that neurotensin expression was increased in women with high insulin resistance. The secretion of this molecule is stimulated by glucose [84] and inhibits the pro-inflammatory status of macrophages when administered under hyperglycemic conditions [85]. Together, our results suggest that there could be a relationship between the decrease in the duodenal immune response and the downregulation of genes involved in the maintenance of the gastrointestinal epithelium in the presence of high insulin resistance. This could constitute an underlying mechanism involved in the increase in the translocation of bacterial components, such as LPS, through the intestinal barrier in individuals with high insulin resistance [51].
Overweight and obesity are factors associated with an increased risk of cancer [86]. In this study, we found in the duodenum of women with morbid obesity a downregulation of a large number of genes that may be involved in antitumor processes, as well as the upregulation of genes with protumor activity in various types of cancer, mainly gastrointestinal cancers. For example, a lower expression of genes that may play an important role in suppressing gastric cancer (GKN1, LIPF, ANXA10, MUC6, PSCA, and SNHG5) was found. On the other hand, we have observed the upregulation of genes that are involved in the proliferation, migration, and invasion of cancer cells in colorectal and gastric cancer (ANGPTL4, S100A8, LINC00668). In the same way, we have obtained similar results in the exclusive DEGs of women with morbid obesity and high or low insulin resistance. Our results show the differential expression of many genes related with different types of cancer (see Supplementary Table S3) in the duodenum of women with morbid obesity, in whom there are still no signs of cancer development. However, we do not know if this altered gene expression will lead to the development of cancer in the future or if other factors are the main drivers of cancer development in this type of subject. Although small intestine cancer is not the most common form of gastrointestinal cancer, some studies point to an increased risk of developing it in patients with obesity [87,88]. However, due to the characteristics of this study, we do not know if the alteration in the expression of these genes may be related to an increased risk of developing small intestine cancer.
This study has several limitations. Firstly, this study was only carried out in women, so we cannot extrapolate the results to men. Additionally, although our results have shown an association between morbid obesity and insulin resistance with the expression of genes related to different pathways, e.g., cancer, more longitudinal studies are needed to confirm our findings.

5. Conclusions

In this study we have shown that morbid obesity, regardless of the degree of insulin resistance, is associated with alterations in intestinal defensive processes and antimicrobial responses and with the activation of pro- and anti-inflammatory genes, which could be associated with possible greater exposure and perhaps vulnerability of the intestine to products and metabolites derived from microorganisms and the diet. Similarly, we also found an alteration in pathways related to the lipid degradation at the intestinal level. More interestingly, we found that morbid obesity is associated with changes in the expression of certain genes that have been related to the development of cancer. Our data also suggest a deficient duodenal immune and antimicrobial response in women with high insulin resistance.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biomedicines10051024/s1, Supplementary Table S1. Anthropometric and biochemical variables for each patient of the three groups of women; Supplementary Table S2. Complete analysis of the overrepresented gene sets (FDR < 0.05), the associated p-value and FDR q-value, and the differentially expressed genes (DEGs); Supplementary Table S3. Bibliographical analysis of differentially expressed genes (DEGs) (FDR < 0.05) involved in the development of cancer [89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122,123,124,125,126,127,128,129,130,131,132,133].

Author Contributions

Conceptualization, A.H.-P. and E.G.-F.; formal analysis, A.H.-P. and E.G.-F.; funding acquisition, E.G.-F.; methodology, A.H.-P., C.R.-D., C.S.-F., C.L.-G., S.G.-S., F.M.-R., F.R.-P., A.R.-C., G.A.-M., L.V.-P. and S.V.; writing—original draft, A.H.-P.; writing—review and editing, A.H.-P., S.V., L.G.-S. and E.G.-F. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by a grant from the Instituto de Salud Carlos III, Spain (PI12/00338). This study was co-funded by FEDER funds (“A way to make Europe”). CIBEREHD and CIBERDEM are funded by ISCIII. AHP is supported by a grant from the Ministerio de Educación y Formación Professional (Spain) (FPU14/01972). FMR is supported by a grant from the ISCIII (Spain) (FI19/00189). EGF is supported by the Nicolas Monardes program from the Consejería de Salud de Andalucía (Spain) (C-0031-2016). CLG is supported by a grant from the Ministerio de Educación y Formación Professional (Spain) (IJCI-2017-31466). CRD is supported by a grant from the ISCIII (Spain) (CD18/00188). FRP is supported by a grant from the ISCIII (Spain) (CP18/00042). LGS is supported by the Miguel Servet program from the ISCIII (Spain) (“Miguel Servet II” program, CPII18/00030) and Nicolas Monardes program from the Consejería de Salud de Andalucía (Spain) (C-0028-2018).

Institutional Review Board Statement

This study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Malaga Provincial Research Ethics Committee of the Biomedical Research Institute of Malaga (IBIMA) (protocol code PI12/00338).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data discussed in this publication have been deposited in the NCBI’s Gene Expression Omnibus (GEO) (GEO Series access number GSE147562).

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Saklayen, M.G. The Global Epidemic of the Metabolic Syndrome. Curr. Hypertens. Rep. 2018, 20, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Samson, S.L.; Garber, A.J. Metabolic syndrome. Endocrinol. Metab. Clin. N. Am. 2014, 43, 1–23. [Google Scholar] [CrossRef] [PubMed]
  3. González-Muniesa, P.; Mártinez-González, M.A.; Hu, F.B.; Després, J.P.; Matsuzawa, Y.; Loos, R.J.F.; Moreno, L.A.; Bray, G.A.; Martinez, J.A. Obesity. Nat. Rev. Dis. Primers 2017, 3, 17034. [Google Scholar] [CrossRef] [PubMed]
  4. Emerenziani, S.; Guarino, M.P.L.; Trillo Asensio, L.M.; Altomare, A.; Ribolsi, M.; Balestrieri, P.; Cicala, M. Role of Overweight and Obesity in Gastrointestinal Disease. Nutrients 2019, 12, 111. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Agustí, A.; García-Pardo, M.P.; López-Almela, I.; Campillo, I.; Maes, M.; Romaní-Pérez, M.; Sanz, Y. Interplay Between the Gut-Brain Axis, Obesity and Cognitive Function. Front. Neurosci. 2018, 12, 155. [Google Scholar] [CrossRef]
  6. Tran, C.D.; Grice, D.M.; Wade, B.; Kerr, C.A.; Bauer, D.C.; Li, D.; Hannan, G.N. Gut permeability, its interaction with gut microflora and effects on metabolic health are mediated by the lymphatics system, liver and bile acid. Future Microbiol. 2015, 10, 1339–1353. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Petersen, M.C.; Shulman, G.I. Mechanisms of Insulin Action and Insulin Resistance. Physiol. Rev. 2018, 98, 2133–2223. [Google Scholar] [CrossRef] [Green Version]
  8. Nguyen, N.T.; Varela, J.E. Bariatric surgery for obesity and metabolic disorders: State of the art. Nat. Rev. Gastroenterol. Hepatol. 2017, 14, 160–169. [Google Scholar] [CrossRef] [PubMed]
  9. Evans, J.M.; Morris, L.S.; Marchesi, J.R. The gut microbiome: The role of a virtual organ in the endocrinology of the host. J. Endocrinol. 2013, 218, R37–R47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Bauer, P.V.; Hamr, S.C.; Duca, F.A. Regulation of energy balance by a gut-brain axis and involvement of the gut microbiota. Cell Mol. Life Sci. 2016, 73, 737–755. [Google Scholar] [CrossRef] [PubMed]
  11. Mucunguzi, O.; Melouane, A.; Ghanemi, A.; Yoshioka, M.; Boivin, A.; Calvo, E.L.; St-Amand, J. Identification of the principal transcriptional regulators for low-fat and high-fat meal responsive genes in small intestine. Nutr. Metab. 2017, 14, 66. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Primeaux, S.D.; Braymer, H.D.; Bray, G.A. High fat diet differentially regulates the expression of olfactory receptors in the duodenum of obesity-prone and obesity-resistant rats. Dig. Dis. Sci. 2013, 58, 72–76. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. de Wit, N.J.; Bosch-Vermeulen, H.; de Groot, P.J.; Hooiveld, G.J.; Bromhaar, M.M.; Jansen, J.; Müller, M.; van der Meer, R. The role of the small intestine in the development of dietary fat-induced obesity and insulin resistance in C57BL/6J mice. BMC Med. Genom. 2008, 1, 14. [Google Scholar] [CrossRef] [PubMed]
  14. Yoshioka, M.; Bolduc, C.; Raymond, V.; St-Amand, J. High-fat meal-induced changes in the duodenum mucosa transcriptome. Obesity 2008, 16, 2302–2307. [Google Scholar] [CrossRef] [PubMed]
  15. Ho-Plagaro, A.; Santiago-Fernandez, C.; García-Serrano, S.; Rodriguez, C.; Garrido-Sanchez, L.; Escamilla, A.; Gonzalo, M.; Montiel-Casado, C.; Alcaín-Martínez, G.; Garcia-Muñoz, B.; et al. A lower duodenal immune response is associated with an increase of insulin resistance in patients with morbid obesity. Int. J. Obes. 2020, 44, 340–352. [Google Scholar] [CrossRef] [PubMed]
  16. Ho-Plagaro, A.; Santiago-Fernandez, C.; Rodríguez-Díaz, C.; Lopez-Gómez, C.; Garcia-Serrano, S.; Rodríguez-Pacheco, F.; Valdes, S.; Rodríguez-Cañete, A.; Alcaín-Martínez, G.; Ruiz-Santana, N.; et al. Different Expression of Duodenal Genes Related to Insulin Resistance Between Nonobese Women and Those with Severe Obesity. Obesity 2020, 28, 1708–1717. [Google Scholar] [CrossRef] [PubMed]
  17. García-Fuentes, E.; García-Almeida, J.M.; García-Arnés, J.; Rivas-Marín, J.; Gallego-Perales, J.L.; González-Jiménez, B.; Cardona, I.; García-Serrano, S.; Garriga, M.J.; Gonzalo, M.; et al. Morbidly obese individuals with impaired fasting glucose have a specific pattern of insulin secretion and sensitivity: Effect of weight loss after bariatric surgery. Obes. Surg. 2006, 16, 1179–1188. [Google Scholar] [CrossRef] [PubMed]
  18. Rojo-Martínez, G.; Esteva, I.; de Adana, S.R.; Catalá, M.; Merelo, M.J.; Tinahones, F.; Gómez-Zumaquero, J.M.; Cuesta, A.L.; Cardona, F.; Soriguer, F. Patterns of insulin resistance in the general population of southeast Spain. Diabetes Res. Clin. Pract. 2004, 65, 247–256. [Google Scholar] [CrossRef] [PubMed]
  19. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [Green Version]
  20. Mootha, V.K.; Lindgren, C.M.; Eriksson, K.F.; Subramanian, A.; Sihag, S.; Lehar, J.; Puigserver, P.; Carlsson, E.; Ridderstråle, M.; Laurila, E.; et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat. Genet. 2003, 34, 267–273. [Google Scholar] [CrossRef] [PubMed]
  21. Lee, M.L.T. Analysis of Microarray Gene Expression Data; Kluwer Academic Publishers: New York, NY, USA, 2004; ISBN 0-7923-7087-2. Available online: https://brac.umd.edu/sample-size-and-power-calculations-microarray-studies.html (accessed on 20 April 2022).
  22. Tsai, C.A.; Wang, S.J.; Chen, D.T.; Chen, J.J. Sample size for gene expression microarray experiments. Bioinformatics 2005, 21, 1502–1508. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Schraufstatter, I.U.; Zhao, M.; Khaldoyanidi, S.K.; Discipio, R.G. The chemokine CCL18 causes maturation of cultured monocytes to macrophages in the M2 spectrum. Immunology 2012, 135, 287–298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Fujita, Y.; Khateb, A.; Li, Y.; Tinoco, R.; Zhang, T.; Bar-Yoseph, H.; Tam, M.A.; Chowers, Y.; Sabo, E.; Gerassy-Vainberg, S.; et al. Regulation of S100A8 Stability by RNF5 in Intestinal Epithelial Cells Determines Intestinal Inflammation and Severity of Colitis. Cell Rep. 2018, 24, 3296–3311.e6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Bealer, J.F.; Colgin, M. S100A8/A9: A potential new diagnostic aid for acute appendicitis. Acad. Emerg. Med. 2010, 17, 333–336. [Google Scholar] [CrossRef] [PubMed]
  26. Kovalovsky, D.; Uche, O.U.; Eladad, S.; Hobbs, R.M.; Yi, W.; Alonzo, E.; Chua, K.; Eidson, M.; Kim, H.J.; Im, J.S.; et al. The BTB-zinc finger transcriptional regulator PLZF controls the development of invariant natural killer T cell effector functions. Nat. Immunol. 2008, 9, 1055–1064. [Google Scholar] [CrossRef]
  27. Van Coillie, S.; Wiernicki, B.; Xu, J. Molecular and Cellular Functions of CTLA-4. Adv. Exp. Med. Biol. 2020, 1248, 7–32. [Google Scholar] [CrossRef] [PubMed]
  28. Ohta, T.; Sugiyama, M.; Hemmi, H.; Yamazaki, C.; Okura, S.; Sasaki, I.; Fukuda, Y.; Orimo, T.; Ishii, K.J.; Hoshino, K.; et al. Crucial roles of XCR1-expressing dendritic cells and the XCR1-XCL1 chemokine axis in intestinal immune homeostasis. Sci. Rep. 2016, 6, 23505. [Google Scholar] [CrossRef] [Green Version]
  29. Furuta, Y.; Tsai, S.H.; Kinoshita, M.; Fujimoto, K.; Okumura, R.; Umemoto, E.; Kurashima, Y.; Kiyono, H.; Kayama, H.; Takeda, K. E-NPP3 controls plasmacytoid dendritic cell numbers in the small intestine. PLoS ONE 2017, 12, e0172509. [Google Scholar] [CrossRef]
  30. Döhler, C.; Zebisch, M.; Sträter, N. Crystal structure and substrate binding mode of ectonucleotide phosphodiesterase/pyrophosphatase-3 (NPP3). Sci. Rep. 2018, 8, 10874. [Google Scholar] [CrossRef]
  31. Lumeng, C.N.; Bodzin, J.L.; Saltiel, A.R. Obesity induces a phenotypic switch in adipose tissue macrophage polarization. J. Clin. Investig. 2007, 117, 175–184. [Google Scholar] [CrossRef] [Green Version]
  32. Calafat, J.; Janssen, H.; Knol, E.F.; Weller, P.F.; Egesten, A. Ultrastructural localization of Charcot-Leyden crystal protein in human eosinophils and basophils. Eur. J. Haematol. 1997, 58, 56–66. [Google Scholar] [CrossRef] [PubMed]
  33. Shyamsunder, P.; Shanmugasundaram, M.; Mayakonda, A.; Dakle, P.; Teoh, W.W.; Han, L.; Kanojia, D.; Lim, M.C.; Fullwood, M.; An, O.; et al. Identification of a novel enhancer of CEBPE essential for granulocytic differentiation. Blood 2019, 133, 2507–2517. [Google Scholar] [CrossRef] [PubMed]
  34. Guo, Y.; Wang, A.Y. Novel Immune Check-Point Regulators in Tolerance Maintenance. Front. Immunol. 2015, 6, 421, Erratum in Front. Immunol. 2016, 7, 38. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Steinhoff, J.S.; Lass, A.; Schupp, M. Biological Functions of RBP4 and Its Relevance for Human Diseases. Front. Physiol. 2021, 12, 659977. [Google Scholar] [CrossRef] [PubMed]
  36. Heissig, B.; Nishida, C.; Tashiro, Y.; Sato, Y.; Ishihara, M.; Ohki, M.; Gritli, I.; Rosenkvist, J.; Hattori, K. Role of neutrophil-derived matrix metalloproteinase-9 in tissue regeneration. Histol. Histopathol. 2010, 25, 765–770. [Google Scholar] [CrossRef] [PubMed]
  37. Couture, A.; Garnier, A.; Docagne, F.; Boyer, O.; Vivien, D.; Le-Mauff, B.; Latouche, J.B.; Toutirais, O. HLA-Class II Artificial Antigen Presenting Cells in CD4+ T Cell-Based Immunotherapy. Front. Immunol. 2019, 10, 1081. [Google Scholar] [CrossRef] [PubMed]
  38. Tokunaga, R.; Zhang, W.; Naseem, M.; Puccini, A.; Berger, M.D.; Soni, S.; McSkane, M.; Baba, H.; Lenz, H.J. CXCL9, CXCL10, CXCL11/CXCR3 axis for immune activation—A target for novel cancer therapy. Cancer Treat. Rev. 2018, 63, 40–47. [Google Scholar] [CrossRef]
  39. Karczewski, J.; Śledzińska, E.; Baturo, A.; Jończyk, I.; Maleszko, A.; Samborski, P.; Begier-Krasińska, B.; Dobrowolska, A. Obesity and inflammation. Eur. Cytokine Netw. 2018, 29, 83–94. [Google Scholar] [CrossRef] [PubMed]
  40. Ryckman, C.; Vandal, K.; Rouleau, P.; Talbot, M.; Tessier, P.A. Proinflammatory activities of S100: Proteins S100A8, S100A9, and S100A8/A9 induce neutrophil chemotaxis and adhesion. J. Immunol. 2003, 170, 3233–3242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Leach, S.T.; Yang, Z.; Messina, I.; Song, C.; Geczy, C.L.; Cunningham, A.M.; Day, A.S. Serum and mucosal S100 proteins, calprotectin (S100A8/S100A9) and S100A12, are elevated at diagnosis in children with inflammatory bowel disease. Scand. J. Gastroenterol. 2007, 42, 1321–1331. [Google Scholar] [CrossRef] [PubMed]
  42. Senger, S.; Sapone, A.; Fiorentino, M.R.; Mazzarella, G.; Lauwers, G.Y.; Fasano, A. Celiac Disease Histopathology Recapitulates Hedgehog Downregulation, Consistent with Wound Healing Processes Activation. PLoS ONE 2015, 10, e0144634. [Google Scholar] [CrossRef] [Green Version]
  43. Allavena, P.; Garlanda, C.; Borrello, M.G.; Sica, A.; Mantovani, A. Pathways connecting inflammation and cancer. Curr. Opin. Genet. Dev. 2008, 18, 3–10. [Google Scholar] [CrossRef] [PubMed]
  44. Vuckovic, S.; Florin, T.H.; Khalil, D.; Zhang, M.F.; Patel, K.; Hamilton, I.; Hart, D.N. CD40 and CD86 upregulation with divergent CMRF44 expression on blood dendritic cells in inflammatory bowel diseases. Am. J. Gastroenterol. 2001, 96, 2946–2956. [Google Scholar] [CrossRef] [PubMed]
  45. Rogler, G.; Hausmann, M.; Spöttl, T.; Vogl, D.; Aschenbrenner, E.; Andus, T.; Falk, W.; Schölmerich, J.; Gross, V. T-cell co-stimulatory molecules are upregulated on intestinal macrophages from inflammatory bowel disease mucosa. Eur. J. Gastroenterol. Hepatol. 1999, 11, 1105–1111. [Google Scholar] [CrossRef] [PubMed]
  46. Sirisinha, S. Insight into the mechanisms regulating immune homeostasis in health and disease. Asian Pac. J. Allergy Immunol. 2011, 29, 1–14. [Google Scholar] [PubMed]
  47. Basu, D.; Goldberg, I.J. Regulation of lipoprotein lipase-mediated lipolysis of triglycerides. Curr. Opin. Lipidol. 2020, 31, 154–160. [Google Scholar] [CrossRef]
  48. Gaudet, D.; Brisson, D.; Tremblay, K.; Alexander, V.J.; Singleton, W.; Hughes, S.G.; Geary, R.S.; Baker, B.F.; Graham, M.J.; Crooke, R.M.; et al. Targeting APOC3 in the familial chylomicronemia syndrome. N. Engl. J. Med. 2014, 371, 2200–2206. [Google Scholar] [CrossRef] [PubMed]
  49. Wang, H.; Eckel, R.H. Lipoprotein lipase: From gene to obesity. Am. J. Physiol. Endocrinol. Metab. 2009, 297, E271–E288. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Li, N.; Fu, J.; Koonen, D.P.; Kuivenhoven, J.A.; Snieder, H.; Hofker, M.H. Are hypertriglyceridemia and low HDL causal factors in the development of insulin resistance? Atherosclerosis 2014, 233, 130–138. [Google Scholar] [CrossRef] [PubMed]
  51. Saad, M.J.; Santos, A.; Prada, P.O. Linking Gut Microbiota and Inflammation to Obesity and Insulin Resistance. Physiology 2016, 31, 283–293. [Google Scholar] [CrossRef] [PubMed]
  52. Xu, W.; Schlagwein, N.; Roos, A.; van den Berg, T.K.; Daha, M.R.; van Kooten, C. Human peritoneal macrophages show functional characteristics of M-CSF-driven anti-inflammatory type 2 macrophages. Eur. J. Immunol. 2007, 37, 1594–1599. [Google Scholar] [CrossRef] [PubMed]
  53. Okumura, R.; Kodama, T.; Hsu, C.C.; Sahlgren, B.H.; Hamano, S.; Kurakawa, T.; Iida, T.; Takeda, K. Lypd8 inhibits attachment of pathogenic bacteria to colonic epithelia. Mucosal. Immunol. 2020, 13, 75–85. [Google Scholar] [CrossRef] [PubMed]
  54. van Beelen Granlund, A.; Østvik, A.E.; Brenna, Ø.; Torp, S.H.; Gustafsson, B.I.; Sandvik, A.K. REG gene expression in inflamed and healthy colon mucosa explored by in situ hybridisation. Cell Tissue Res. 2013, 352, 639–646. [Google Scholar] [CrossRef] [Green Version]
  55. Sias, B.; Ferrato, F.; Pellicer-Rubio, M.T.; Forgerit, Y.; Guillouet, P.; Leboeuf, B.; Carrière, F. Cloning and seasonal secretion of the pancreatic lipase-related protein 2 present in goat seminal plasma. Biochim. Biophys. Acta 2005, 1686, 169–180. [Google Scholar] [CrossRef] [PubMed]
  56. Alves, B.N.; Leong, J.; Tamang, D.L.; Elliott, V.; Edelnant, J.; Redelman, D.; Singer, C.A.; Kuhn, A.R.; Miller, R.; Lowe, M.E.; et al. Pancreatic lipase-related protein 2 (PLRP2) induction by IL-4 in cytotoxic T lymphocytes (CTLs) and reevaluation of the negative effects of its gene ablation on cytotoxicity. J. Leukoc. Biol. 2009, 86, 701–712. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Gilleron, M.; Lepore, M.; Layre, E.; Cala-De Paepe, D.; Mebarek, N.; Shayman, J.A.; Canaan, S.; Mori, L.; Carrière, F.; Puzo, G.; et al. Lysosomal Lipases PLRP2 and LPLA2 Process Mycobacterial Multi-acylated Lipids and Generate T Cell Stimulatory Antigens. Cell Chem. Biol. 2016, 23, 1147–1156. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. López-Fandiño, R. Role of dietary lipids in food allergy. Crit. Rev. Food Sci. Nutr. 2020, 60, 1797–1814. [Google Scholar] [CrossRef] [PubMed]
  59. Broch, M.; Ramírez, R.; Auguet, M.T.; Alcaide, M.J.; Aguilar, C.; Garcia-España, A.; Richart, C. Macrophages are novel sites of expression and regulation of retinol binding protein-4 (RBP4). Physiol. Res. 2010, 59, 299–303. [Google Scholar] [CrossRef] [PubMed]
  60. Moraes-Vieira, P.M.; Castoldi, A.; Aryal, P.; Wellenstein, K.; Peroni, O.D.; Kahn, B.B. Antigen Presentation and T-Cell Activation Are Critical for RBP4-Induced Insulin Resistance. Diabetes 2016, 65, 1317–1327. [Google Scholar] [CrossRef] [Green Version]
  61. Cullen, J.J.; Mercer, D.; Hinkhouse, M.; Ephgrave, K.S.; Conklin, J.L. Effects of endotoxin on regulation of intestinal smooth muscle nitric oxide synthase and intestinal transit. Surgery 1999, 125, 339–344. [Google Scholar] [CrossRef] [PubMed]
  62. Xi, L.; Tekin, D.; Gursoy, E.; Salloum, F.; Levasseur, J.E.; Kukreja, R.C. Evidence that NOS2 acts as a trigger and mediator of late preconditioning induced by acute systemic hypoxia. Am. J. Physiol. Heart Circ. Physiol. 2002, 283, H5–H12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Lei, Y.; Takahama, Y. XCL1 and XCR1 in the immune system. Microbes Infect. 2012, 14, 262–267. [Google Scholar] [CrossRef] [PubMed]
  64. Jubin, V.; Ventre, E.; Leverrier, Y.; Djebali, S.; Mayol, K.; Tomkowiak, M.; Mafille, J.; Teixeira, M.; Teoh, D.Y.; Lina, B.; et al. T inflammatory memory CD8 T cells participate to antiviral response and generate secondary memory cells with an advantage in XCL1 production. Immunol. Res. 2012, 52, 284–293. [Google Scholar] [CrossRef] [PubMed]
  65. Nevins, A.M.; Subramanian, A.; Tapia, J.L.; Delgado, D.P.; Tyler, R.C.; Jensen, D.R.; Ouellette, A.J.; Volkman, B.F. A Requirement for Metamorphic Interconversion in the Antimicrobial Activity of Chemokine XCL1. Biochemistry 2016, 55, 3784–3793. [Google Scholar] [CrossRef] [PubMed]
  66. Chapoval, A.I.; Smithson, G.; Brunick, L.; Mesri, M.; Boldog, F.L.; Andrew, D.; Khramtsov, N.V.; Feshchenko, E.A.; Starling, G.C.; Mezes, P.S. BTNL8, a butyrophilin-like molecule that costimulates the primary immune response. Mol. Immunol. 2013, 56, 819–828. [Google Scholar] [CrossRef] [PubMed]
  67. Walmsley, S.R.; Chilvers, E.R.; Thompson, A.A.; Vaughan, K.; Marriott, H.M.; Parker, L.C.; Shaw, G.; Parmar, S.; Schneider, M.; Sabroe, I.; et al. Prolyl hydroxylase 3 (PHD3) is essential for hypoxic regulation of neutrophilic inflammation in humans and mice. J. Clin. Investig. 2011, 121, 1053–1063. [Google Scholar] [CrossRef] [PubMed]
  68. Taniguchi, C.M.; Finger, E.C.; Krieg, A.J.; Wu, C.; Diep, A.N.; LaGory, E.L.; Wei, K.; McGinnis, L.M.; Yuan, J.; Kuo, C.J.; et al. Cross-talk between hypoxia and insulin signaling through Phd3 regulates hepatic glucose and lipid metabolism and ameliorates diabetes. Nat. Med. 2013, 19, 1325–1330. [Google Scholar] [CrossRef] [Green Version]
  69. De Re, V.; Simula, M.P.; Caggiari, L.; Orzes, N.; Spina, M.; Da Ponte, A.; De Appollonia, L.; Dolcetti, R.; Canzonieri, V.; Cannizzaro, R. Proteins specifically hyperexpressed in a coeliac disease patient with aberrant T cells. Clin. Exp. Immunol. 2007, 148, 402–409. [Google Scholar] [CrossRef]
  70. Ishiguro, K.; Watanabe, O.; Nakamura, M.; Yamamura, T.; Ando, T.; Goto, H.; Hirooka, Y. S100G expression and function in fibroblasts on colitis induction. Int. Immunopharmacol. 2016, 39, 92–96. [Google Scholar] [CrossRef] [PubMed]
  71. Strober, W. Impact of the gut microbiome on mucosal inflammation. Trends Immunol. 2013, 34, 423–430. [Google Scholar] [CrossRef] [Green Version]
  72. Thomsen, M.K.; Hansen, G.H.; Danielsen, E.M. Galectin-2 at the enterocyte brush border of the small intestine. Mol. Membr. Biol. 2009, 26, 347–355. [Google Scholar] [CrossRef] [PubMed]
  73. Gonzalez, F.J.; Xie, C.; Jiang, C. The role of hypoxia-inducible factors in metabolic diseases. Nat. Rev. Endocrinol. 2018, 15, 21–32. [Google Scholar] [CrossRef] [PubMed]
  74. González-Muniesa, P.; de Oliveira, C.; Pérez de Heredia, F.; Thompson, M.P.; Trayhurn, P. Fatty acids and hypoxia stimulate the expression and secretion of the adipokine ANGPTL4 (angiopoietin-like protein 4/ fasting-induced adipose factor) by human adipocytes. J. Nutrigenet. Nutr. 2011, 4, 146–153. [Google Scholar] [CrossRef] [PubMed]
  75. Yao, Q.; Shin, M.K.; Jun, J.C.; Hernandez, K.L.; Aggarwal, N.R.; Mock, J.R.; Gay, J.; Drager, L.F.; Polotsky, V.Y. Effect of chronic intermittent hypoxia on triglyceride uptake in different tissues. J. Lipid Res. 2013, 54, 1058–1065. [Google Scholar] [CrossRef] [Green Version]
  76. Shah, Y.M.; Matsubara, T.; Ito, S.; Yim, S.H.; Gonzalez, F.J. Intestinal hypoxia-inducible transcription factors are essential for iron absorption following iron deficiency. Cell Metab. 2009, 9, 152–164. [Google Scholar] [CrossRef] [Green Version]
  77. Grebhardt, S.; Veltkamp, C.; Ströbel, P.; Mayer, D. Hypoxia and HIF-1 increase S100A8 and S100A9 expression in prostate cancer. Int. J. Cancer. 2012, 131, 2785–2794. [Google Scholar] [CrossRef] [PubMed]
  78. Grek, C.L.; Newton, D.A.; Spyropoulos, D.D.; Baatz, J.E. Hypoxia up-regulates expression of hemoglobin in alveolar epithelial cells. Am. J. Respir. Cell Mol. Biol. 2011, 44, 439–447. [Google Scholar] [CrossRef]
  79. Gabryelska, A.; Karuga, F.F.; Szmyd, B.; Białasiewicz, P. HIF-1α as a Mediator of Insulin Resistance, T2DM, and Its Complications: Potential Links With Obstructive Sleep Apnea. Front. Physiol. 2020, 11, 1035. [Google Scholar] [CrossRef]
  80. Burke, B. The role of matrix metalloproteinase 7 in innate immunity. Immunobiology 2004, 209, 51–56, Erratum in Immunobiology 2005, 209, 755. [Google Scholar] [CrossRef] [PubMed]
  81. Xu, Y.; Gao, Q.; Xue, Y.; Li, X.; Xu, L.; Li, C.; Qin, Y.; Fang, J. Prolyl hydroxylase 3 stabilizes the p53 tumor suppressor by inhibiting the p53-MDM2 interaction in a hydroxylase-independent manner. J. Biol. Chem. 2019, 294, 9949–9958. [Google Scholar] [CrossRef]
  82. Beneke, A.; Guentsch, A.; Hillemann, A.; Zieseniss, A.; Swain, L.; Katschinski, D.M. Loss of PHD3 in myeloid cells dampens the inflammatory response and fibrosis after hind-limb ischemia. Cell Death Dis. 2017, 8, e2976. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Gutiérrez-Repiso, C.; Moreno-Indias, I.; Martín-Núñez, G.M.; Ho-Plágaro, A.; Rodríguez-Cañete, A.; Gonzalo, M.; García-Fuentes, E.; Tinahones, F.J. Mucosa-associated microbiota in the jejunum of patients with morbid obesity: Alterations in states of insulin resistance and metformin treatment. Surg. Obes. Relat. Dis. 2020, 16, 1575–1585. [Google Scholar] [CrossRef] [PubMed]
  84. Kuhre, R.E.; Bechmann, L.E.; Wewer Albrechtsen, N.J.; Hartmann, B.; Holst, J.J. Glucose stimulates neurotensin secretion from the rat small intestine by mechanisms involving SGLT1 and GLUT2, leading to cell depolarization and calcium influx. Am. J. Physiol. Endocrinol. Metab. 2015, 308, E1123–E1130. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Moura, L.I.; Silva, L.; Leal, E.C.; Tellechea, A.; Cruz, M.T.; Carvalho, E. Neurotensin modulates the migratory and inflammatory response of macrophages under hyperglycemic conditions. Biomed. Res. Int. 2013, 2013, 941764. [Google Scholar] [CrossRef] [Green Version]
  86. Avgerinos, K.I.; Spyrou, N.; Mantzoros, C.S.; Dalamaga, M. Obesity and cancer risk: Emerging biological mechanisms and perspectives. Metabolism 2019, 92, 121–135. [Google Scholar] [CrossRef] [PubMed]
  87. Barsouk, A.; Rawla, P.; Barsouk, A.; Thandra, K.C. Epidemiology of Cancers of the Small Intestine: Trends, Risk Factors, and Prevention. Med. Sci. 2019, 7, 46. [Google Scholar] [CrossRef] [Green Version]
  88. Pan, S.Y.; Morrison, H. Epidemiology of cancer of the small intestine. World J. Gastrointest. Oncol. 2011, 3, 33–42. [Google Scholar] [CrossRef]
  89. Kim, S.-H.; Park, Y.-Y.; Kim, S.-W.; Lee, J.-S.; Wang, D.; Dubois, R.N. ANGPTL4 Induction by Prostaglandin E2 under Hypoxic Conditions Promotes Colorectal Cancer Progression. Cancer Res. 2011, 71, 7010–7020. [Google Scholar] [CrossRef] [Green Version]
  90. Kwon, C.H.; Moon, H.J.; Park, H.J.; Choi, J.H.; Park, D.Y. S100A8 and S100A9 promotes invasion and migration through p38 mitogen-activated protein kinase-dependent NF-κB activation in gastric cancer cells. Mol. Cells 2013, 35, 226–234. [Google Scholar] [CrossRef] [Green Version]
  91. Dianatpour, A.; Ghafouri-Fard, S. Long Non Coding RNA Expression Intersecting Cancer and Spermatogenesis: A Systematic Review. Asian Pac. J. Cancer Prev. 2017, 18, 2601–2610. [Google Scholar] [CrossRef]
  92. Hui, A.W.-H.; Lau, H.-W.; Cao, C.Y.; Zhou, J.-W.; Lai, P.B.-S.; Tsui, S.K.-W. Downregulation of PLZF in human hepatocellular carcinoma and its clinical significance. Oncol. Rep. 2014, 33, 397–402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. Peng, Y.; Luo, G.; Zhou, J.; Wang, X.; Hu, J.; Cui, Y.; Li, X.C.; Tan, J.; Yang, S.; Zhan, R.; et al. CD86 Is an Activation Receptor for NK Cell Cytotoxicity against Tumor Cells. PLoS ONE 2013, 8, e83913. [Google Scholar] [CrossRef] [PubMed]
  94. Mao, W.; Chen, J.; Peng, T.-L.; Yin, X.-F.; Chen, L.-Z.; Chen, M.-H. Downregulation of gastrokine-1 in gastric cancer tissues and restoration of its expression induced gastric cancer cells to apoptosis. J. Exp. Clin. Cancer Res. 2012, 31, 49. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Yoon, J.H.; Choi, W.S.; Kim, O.; Choi, B.J.; Nam, S.W.; Lee, J.Y.; Park, W.S. Gastrokine 1 inhibits gastric cancer cell migration and invasion by downregulating RhoA expression. Gastric Cancer 2016, 20, 274–285. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Menheniott, T.R.; O’Connor, L.; Chionh, Y.T.; Däbritz, J.; Scurr, M.; Rollo, B.N.; Ng, G.; Jacobs, S.; Catubig, A.; Kurklu, B.; et al. Loss of gastrokine-2 drives premalignant gastric inflammation and tumor progression. J. Clin. Investig. 2016, 126, 1383–1400. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. Kong, Y.; Zheng, Y.; Jia, Y.; Li, P.; Wang, Y. Decreased LIPF expression is correlated with DGKA and predicts poor outcome of gastric cancer. Oncol. Rep. 2016, 36, 1852–1860. [Google Scholar] [CrossRef] [Green Version]
  98. Kim, J.K.; Kim, P.J.; Jung, K.H.; Noh, J.H.; Eun, J.W.; Bae, H.J.; Xie, H.J.; Shan, J.M.; Ping, W.Y.; Park, W.S.; et al. Decreased ex-pression of annexin A10 in gastric cancer and its overexpression in tumor cell growth suppression. Oncol. Rep. 2010, 24, 607–612. [Google Scholar]
  99. Leir, S.-H.; Harris, A. MUC6 mucin expression inhibits tumor cell invasion. Exp. Cell Res. 2011, 317, 2408–2419. [Google Scholar] [CrossRef]
  100. Heinrich, M.-C.; Göbel, C.; Kluth, M.; Bernreuther, C.; Sauer, C.; Schroeder, C.; Möller-Koop, C.; Hube-Magg, C.; Lebok, P.; Burandt, E.; et al. PSCA expression is associated with favorable tumor features and reduced PSA recurrence in operated prostate cancer. BMC Cancer 2018, 18, 612. [Google Scholar] [CrossRef] [Green Version]
  101. Zhao, L.; Guo, H.; Zhou, B.; Feng, J.; Li, Y.; Han, T.; Liu, L.; Li, L.; Zhang, S.; Liu, Y.; et al. Long non-coding RNA SNHG5 suppresses gastric cancer progression by trapping MTA2 in the cytosol. Oncogene 2016, 35, 5770–5780. [Google Scholar] [CrossRef]
  102. Bergström, S.H.; Järemo, H.; Nilsson, M.; Adamo, H.H.; Bergh, A. Prostate tumors downregulate microseminoprotein-beta (MSMB) in the surrounding benign prostate epithelium and this response is associated with tumor aggressiveness. Prostate 2017, 78, 257–265. [Google Scholar] [CrossRef] [PubMed]
  103. Dahlman, A.; Rexhepaj, E.; Brennan, D.J.; Gallagher, W.M.; Gaber, A.; Lindgren, A.; Jirström, K.; Bjartell, A. Evaluation of the prognostic significance of MSMB and CRISP3 in prostate cancer using automated image analysis. Mod. Pathol. 2011, 24, 708–719. [Google Scholar] [CrossRef] [Green Version]
  104. Whitaker, H.C.; Kote-Jarai, Z.; Ross-Adams, H.; Warren, A.Y.; Burge, J.; George, A.; Bancroft, E.; Jhavar, S.; Leongamornlert, D.; Tymrakiewicz, M.; et al. The rs10993994 Risk Allele for Prostate Cancer Results in Clinically Relevant Changes in Microseminoprotein-Beta Expression in Tissue and Urine. PLoS ONE 2010, 5, e13363. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  105. Han, Y.; Yang, Y.-N.; Yuan, H.-H.; Zhang, T.-T.; Sui, H.; Wei, X.-L.; Liu, L.; Huang, P.; Zhang, W.-J.; Bai, Y.-X. UCA1, a long non-coding RNA up-regulated in colorectal cancer influences cell proliferation, apoptosis and cell cycle distribution. Pathology 2014, 46, 396–401. [Google Scholar] [CrossRef] [PubMed]
  106. Xia, Y.-J.; Jiang, X.-T.; Jiang, S.-B.; He, X.-J.; Luo, J.-G.; Liu, Z.-C.; Wang, L.; Tao, H.-Q.; Chen, J.-Z. PHD3 affects gastric cancer progression by negatively regulating HIF1A. Mol. Med. Rep. 2017, 16, 6882–6889. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  107. Yerushalmi, H.F.; Besselsen, D.G.; Ignatenko, N.A.; Blohm-Mangone, K.A.; Padilla-Torres, J.L.; Stringer, D.E.; Cui, H.; Holubec, H.; Payne, C.M.; Gerner, E.W. The role of NO synthases in arginine-dependent small intestinal and colonic carcinogenesis. Mol. Carcinog. 2005, 45, 93–105. [Google Scholar] [CrossRef] [PubMed]
  108. Ji, Y.; Li, J.; Li, P.; Wang, L.; Yang, H.; Jiang, G. C/EBPβ Promotion of MMP3-Dependent Tumor Cell Invasion and Association with Metastasis in Colorectal Cancer. Genet. Test. Mol. Biomark. 2018, 22, 5–10. [Google Scholar] [CrossRef]
  109. Hua, Y.; Liang, C.; Zhu, J.; Miao, C.; Yu, Y.; Xu, A.; Zhang, J.; Li, P.; Li, S.; Bao, M.; et al. Expression of lactate dehydrogenase C correlates with poor prognosis in renal cell carcinoma. Tumor Biol. 2017, 39. [Google Scholar] [CrossRef] [Green Version]
  110. Garg, M.; Braunstein, G.; Koeffler, H.P. LAMC2 as a therapeutic target for cancers. Expert Opin. Ther. Targets 2014, 18, 979–982. [Google Scholar] [CrossRef] [Green Version]
  111. Cho, N.H.; Shim, H.S.; Rha, S.Y.; Kang, S.H.; Hong, S.H.; Choi, Y.D.; Cho, S.H. Increased Expression of Matrix Metalloproteinase 9 Correlates with Poor Prognostic Variables in Renal Cell Carcinoma. Eur. Urol. 2003, 44, 560–566. [Google Scholar] [CrossRef]
  112. Noda, T.; Yamamoto, H.; Takemasa, I.; Yamada, D.; Uemura, M.; Wada, H.; Kobayashi, S.; Marubashi, S.; Eguchi, H.; Tanemura, M.; et al. PLOD2 induced under hypoxia is a novel prognostic factor for hepatocellular carcinoma after curative resection. Liver Int. 2011, 32, 110–118. [Google Scholar] [CrossRef] [PubMed]
  113. Du, W.; Liu, N.; Zhang, Y.; Liu, X.; Yang, Y.; Chen, W.; He, Y. PLOD2 promotes aerobic glycolysis and cell progression in colorectal cancer by upregulating HK2. Biochem. Cell Biol. 2020, 98, 386–395. [Google Scholar] [CrossRef] [PubMed]
  114. Yuan, K.; Gao, Z.J.; Yuan, W.D.; Yuan, J.Q.; Wang, Y. High expression of SLC6A10P contributes to poor prognosis in lung adeno-carcinoma. Int. J. Clin. Exp. Pathol. 2018, 11, 720–726. [Google Scholar] [PubMed]
  115. Lewis, C.; Murdoch, C. Macrophage Responses to Hypoxia: Implications for Tumor Progression and Anti-Cancer Therapies. Am. J. Pathol. 2005, 167, 627–635. [Google Scholar] [CrossRef]
  116. Wu, N.; Song, Y.; Pang, L.; Chen, Z. CRCT1 regulated by microRNA-520 g inhibits proliferation and induces apoptosis in esophageal squamous cell cancer. Tumor Biol. 2015, 37, 8271–8279. [Google Scholar] [CrossRef] [PubMed]
  117. Oh, S.-J.; Kim, S.-M.; Kim, Y.-O.; Chang, H.-K. Clinicopathologic Implications of PIWIL2 Expression in Colorectal Cancer. Korean J. Pathol. 2012, 46, 318–323. [Google Scholar] [CrossRef]
  118. Blumenthal, R.D.; Leon, E.; Hansen, H.J.; Goldenberg, D.M. Expression patterns of CEACAM5 and CEACAM6 in primary and metastatic cancers. BMC Cancer 2007, 7, 2. [Google Scholar] [CrossRef] [Green Version]
  119. Chen, D.; Sun, Y.; Wei, Y.; Zhang, P.; Rezaeian, A.H.; Teruya-Feldstein, J.; Gupta, S.; Liang, H.; Lin, H.-K.; Hung, M.-C.; et al. LIFR is a breast cancer metastasis suppressor upstream of the Hippo-YAP pathway and a prognostic marker. Nat. Med. 2012, 18, 1511–1517. [Google Scholar] [CrossRef] [Green Version]
  120. Zhao, J.-H.; Sun, J.-X.; Song, Y.-X.; Chen, X.-W.; Yang, Y.-C.; Ma, B.; Wang, J.; Gao, P.; Wang, Z.-N. A novel long noncoding RNA-LOWEG is low expressed in gastric cancer and acts as a tumor suppressor by inhibiting cell invasion. J. Cancer Res. Clin. Oncol. 2015, 142, 601–609. [Google Scholar] [CrossRef]
  121. Zhao, E.-H.; Shen, Z.-Y.; Liu, H.; Jin, X.; Cao, H. Clinical significance of human kallikrein 12 gene expression in gastric cancer. World J. Gastroenterol. 2012, 18, 6597–6604. [Google Scholar] [CrossRef]
  122. Sun, M.; Jin, F.-Y.; Xia, R.; Kong, R.; Li, J.-H.; Xu, T.-P.; Liu, Y.-W.; Zhang, E.-B.; Liu, X.-H.; De, W. Decreased expression of long noncoding RNA GAS5 indicates a poor prognosis and promotes cell proliferation in gastric cancer. BMC Cancer 2014, 14, 319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  123. Chen, Z.; Ju, H.; Yu, S.; Zhao, T.; Jing, X.; Li, P.; Jia, J.; Li, N.; Tan, B.; Li, Y. Prader–Willi region non-protein coding RNA 1 suppressed gastric cancer growth as a competing endogenous RNA of miR-425-5p. Clin. Sci. 2018, 132, 1003–1019. [Google Scholar] [CrossRef] [PubMed]
  124. Bao, W.; Fu, H.-J.; Jia, L.-T.; Zhang, Y.; Li, W.; Jin, B.-Q.; Yao, L.-B.; Chen, S.-Y.; Yang, A.-G. HER2-mediated upregulation of MMP-1 is involved in gastric cancer cell invasion. Arch. Biochem. Biophys. 2010, 499, 49–55. [Google Scholar] [CrossRef] [PubMed]
  125. Kim, C.J.; Cho, Y.G.; Jeong, S.W.; Kim, Y.S.; Kim, S.Y.; Nam, S.W.; Lee, S.H.; Yoo, N.J.; Lee, J.Y.; Park, W.S. Altered expression of KCNK9 in colorectal cancers. APMIS 2004, 112, 588–594. [Google Scholar] [CrossRef]
  126. Liot, S.; Aubert, A.; Hervieu, V.; El Kholti, N.; Schalkwijk, J.; Verrier, B.; Valcourt, U.; Lambert, E. Loss of Tenascin-X expression during tumor progression: A new pan-cancer marker. Matrix Biol. Plus 2020, 6-7, 100021. [Google Scholar] [CrossRef]
  127. Auslander, N.; E Cunningham, C.; Toosi, B.M.; McEwen, E.J.; Yizhak, K.; Vizeacoumar, F.S.; Parameswaran, S.; Gonen, N.; Freywald, T.; Bhanumathy, K.K.; et al. An integrated computational and experimental study uncovers FUT 9 as a metabolic driver of colorectal cancer. Mol. Syst. Biol. 2017, 13, 956. [Google Scholar] [CrossRef]
  128. Jing, Y.; Jia, D.; Wong, C.-M.; Ng, I.O.-L.; Zhang, Z.; Liu, L.; Wang, Q.; Zhao, F.; Li, J.; Yao, M.; et al. SERPINA5 inhibits tumor cell migration by modulating the fibronectin-integrin β1 signaling pathway in hepatocellular carcinoma. Mol. Oncol. 2013, 8, 366–377. [Google Scholar] [CrossRef]
  129. Bijsmans, I.T.G.W.; Smits, K.M.; De Graeff, P.; A Wisman, G.B.; Van Der Zee, A.G.J.; Slangen, B.F.; De Bruïne, A.P.; Van Engeland, M.; Sieben, N.L.; Van De Vijver, K.K. Loss of SerpinA5 protein expression is associated with advanced-stage serous ovarian tumors. Mod. Pathol. 2010, 24, 463–470. [Google Scholar] [CrossRef]
  130. Yang, X.; Wood, P.A.; Ansell, C.M.; Quiton, D.F.T.; Oh, E.-Y.; Du-Quiton, J.; Hrushesky, W.J.M. The Circadian Clock Gene PER1 Suppresses Cancer Cell Proliferation and Tumor Growth at Specific Times of Day. Chronobiol. Int. 2009, 26, 1323–1339. [Google Scholar] [CrossRef]
  131. Mostafaie, N.; Kállay, E.; Sauerzapf, E.; Bonner, E.; Kriwanek, S.; Cross, H.S.; Huber, K.R.; Krugluger, W. Correlated downregulation of estrogen receptor beta and the circadian clock gene Per1 in human colorectal cancer. Mol. Carcinog. 2009, 48, 642–647. [Google Scholar] [CrossRef]
  132. Bachelor, M.A.; Lu, Y.; Owens, D.M. l-3-Phosphoserine phosphatase (PSPH) regulates cutaneous squamous cell carcinoma proliferation independent of l-serine biosynthesis. J. Dermatol. Sci. 2011, 63, 164–172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  133. Li, X.; Wang, L.; Nunes, D.P.; Troxler, R.F.; Offner, G.D. Suppression of MUC1 synthesis downregulates expression of the epidermal growth factor receptor. Cancer Biol. Ther. 2005, 4, 968–973. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Biomedicines 10 01024 g001 550
Figure 1. Venn diagram of the differentially expressed genes (DEGs) (FDR < 0.05) in the microarray. We obtained two lists of DEGs between the MO-high-IR and NO-low-IR groups and between the MO-low-IR and NO-low-IR groups. NO: non-obese women. MO: women with morbid obesity.
Figure 1. Venn diagram of the differentially expressed genes (DEGs) (FDR < 0.05) in the microarray. We obtained two lists of DEGs between the MO-high-IR and NO-low-IR groups and between the MO-low-IR and NO-low-IR groups. NO: non-obese women. MO: women with morbid obesity.
Biomedicines 10 01024 g001
Biomedicines 10 01024 g002 550
Figure 2. Representation of the Log2 of fold change of the differentially expressed genes (DEGs) included in the overrepresented gene sets (FDR < 0.05) related to digestion and lipid metabolism processes.
Figure 2. Representation of the Log2 of fold change of the differentially expressed genes (DEGs) included in the overrepresented gene sets (FDR < 0.05) related to digestion and lipid metabolism processes.
Biomedicines 10 01024 g002
Biomedicines 10 01024 g003 550
Figure 3. Representation of the Log2 of fold change of the differentially expressed genes (DEGs) included in the overrepresented gene sets (FDR < 0.05) related to defense response, inflammation, and innate and adaptive immune responses.
Figure 3. Representation of the Log2 of fold change of the differentially expressed genes (DEGs) included in the overrepresented gene sets (FDR < 0.05) related to defense response, inflammation, and innate and adaptive immune responses.
Biomedicines 10 01024 g003
Biomedicines 10 01024 g004 550
Figure 4. Representation of the Log2 of fold change of the differentially expressed genes (DEGs) related to pro- and anti-inflammatory processes.
Figure 4. Representation of the Log2 of fold change of the differentially expressed genes (DEGs) related to pro- and anti-inflammatory processes.
Biomedicines 10 01024 g004
Biomedicines 10 01024 g005 550
Figure 5. Representation of the Log2 of fold change of the differentially expressed genes (DEGs) included in the overrepresented gene sets (FDR < 0.05) related to antimicrobial response and response to another organism.
Figure 5. Representation of the Log2 of fold change of the differentially expressed genes (DEGs) included in the overrepresented gene sets (FDR < 0.05) related to antimicrobial response and response to another organism.
Biomedicines 10 01024 g005
Table
Table 1. Anthropometric and biochemical variables of the three groups of women.
Table 1. Anthropometric and biochemical variables of the three groups of women.
NO-Lower-IRMO-Lower-IRMO-Higher-IR
N678
Age (years)43.6 (11.0)48.0 (10.0)45.0 (17.0)
Weight (kg)55.5 (14.0)110.6 (21.0) 3126.3 (34.4) 3
BMI (kg/m2)22.1 (6.9)46.02 (6.7) 350.9 (14.3) 3
Waist (cm)74.5 (14.3)116.0 (18.0) 3131.0 (10.5) 3
Hip (cm)103.5 (6.0)143 0 (31.0) 2151.5 (23.5) 2
Glucose (mg/dL)78.0 (13.0)84.0 (7.0)95.5 (19.0)
Insulin (µIU/mL)8.0 (2.5)10.5 (4.8)24.3 (8.9) 3,*
Cholesterol (mg/dL)194.5 (32.0)186.0 (53.0)188.0 (69.0)
Triglycerides (mg/dL)83.5 (39.0)110.0 (57.0)142.0 (84.0)
HOMA-IR1.47 (0.54)2.13 (1.10)5.2 (2.13) 3,#
Data given as median (interquartile range). Significant differences between MO-low-IR and MO-high-IR groups: * p < 0.05, # p < 0.01. Significant differences between NO-low-IR and MO-low-IR groups or between NO-low-IR and MO-high-IR groups: 2 p < 0.01, 3 p < 0.001. HOMA-IR, homeostatic model assessment of insulin resistance index; MO-high-IR, women with morbid obesity with high insulin resistance; MO-low-IR, women with morbid obesity with low insulin resistance; NO-low-IR, nonobese women with low insulin resistance.
Table
Table 2. Significantly overrepresented (FDR < 0.05) gene sets involved in different functions obtained from the analysis of differentially expressed genes (DEGs) found in the intersection, from exclusive DEGs found in MO-low-IR and from exclusive DEGs found in MO-high-IR.
Table 2. Significantly overrepresented (FDR < 0.05) gene sets involved in different functions obtained from the analysis of differentially expressed genes (DEGs) found in the intersection, from exclusive DEGs found in MO-low-IR and from exclusive DEGs found in MO-high-IR.
Function aExclusive DEGs in MO-Lower-IRDEGs in IntersectionExclusive DEGs in MO-Higher-IR
Defense response and inflammation [1]BTNL8, CDA, CEACAM6, CEBPE, CLC, CRCT1, CYP2S1, EGLN3, ENPP3, FBXO2, FCGR1B, HLA-DQA2, KBTBD13, LAIR2, LAMC2, MMP7, MMP9, NOS2, PLOD2, RBP4, SLC40A1, SORD, XCL1, ZDHHC11, BICC1, CYP4F11, EPB42, ERAP2, FAP, HLA-DQB1, KRT1, LIFR, PWRN1ANGPTL4, BPIFB1 *, CCL18, GREM1, GREM2, HBA2, HBB, LYPD8, ODAM, PGC *, PHEX, S100A8, TFF2 *, TCN1 *, CD86, ELF5, FKBP5, GKN2, LPL, MUC6, MYH7, PNLIPRP2, PSPHP1, REG1B, ZBTB16CR2, DNAJB13, MMP1, MMP3, SLCO4C1, XPNPEP3, APOD, CXCL9, ENTPD3, F5, FBXL13, FUT9, GBP3, KCNN4, LRRC17, MUC1, PCDHGC3, PER1, PSPH, SERPINA3, SERPINA5, SLPI, SPINK5, TNXB, ZG16B
Adaptative immune response [2]BTNL8, CLC, FBXO2, FCGR1B, HLA-DQA2, KBTBD13, LAIR2, XCL1, ERAP2, HLA-DQB1CD86, ZBTB16CR2, FBXL13
Innate immune response [3]CDA, CEACAM6, FCGR1B, HLA-DQA2, MMP9, NOS2, XCL1, ERAP2, HLA-DQB1, KRT1BPIFB1 *, CCL18, HBB, PGC *, PHEX, S100A8, TCN1 *, MUC6, REG1BCR2,SLCO4C1, CXCL9, F5, GBP3, MUC1, SERPINA3, SLPI, SPINK5
Inflammation [4]ENPP3,LAMC2, MMP9, NOS2, RBP4, XCL1, BICC1, CYP4F11, EPB42, KRT1, HLA-DQB1CCL18, HBA2, HBB, ODAM,S100A8, LPL, ZBTB16CR2, DNAJB13, MMP1, MMP3, XPNPEP3, APOD, CXCL9, F5, PER1, PSPH, SERPINA3, SPINK5
Mediators in the production, signaling and response to cytokines [5]FCGR1B, HLA-DQA2, MMP9, NOS2, XCL1, HLA-DQB1, LIFRCCL18, GREM2, TFF2 *, CD86MMP1, MMP3, CXCL9, GBP3, MUC1,
Antimicrobial response [6]CEBPE,FCGR1B, HLA-DQA2, NOS2, XCL1, ZDHHC11, HLA-DQB1, KRT1BPIFB1 *, CCL18, LYPD8, PGC *, S100A8, CD86, FKBP5, GKN2, LPL, MUC6, PNLIPRP2, REG1BCR2, CXCL9, GBP3, MUC1, SLPI, SPINK5
Hypoxia [7]ARSL, EGLN3, MMP7, NOS2, PLOD2, SLC6A8, SLC6A10P, CYP4F11, HLA-DQB1ANGPTL4APOD, MUC1
Epithelial maintenance and wound healing [8]CYP4F2, RBP4, CYP4F11, FAP, KRT1HBD, HBB, ODAM, S100A8, TFF2 *, MUC6HBG1, APOD, F5, SERPINA3, SERPINA5
Homeostasis [9]CYP4F2, ERN1, NOS2, PM20D1, RBP4, SCT, SLC30A10, SLC40A1, XCL1, EPB42, FTO, KRT1APOC3, ANGPTL4, S100A8, SLC11A2, TFF2 *, CCDC66, CLRN1, HMBOX1, JPH4, LPL, MUC6MYOC, SCN3B, CXCL9, KCNN4, SERPINA3, ZG16B
Fucosylation [10]FUT2 FUT9, FUT1
Cell proliferation [11]CEACAM6, CLC, EGLN3, ENPP3, ERN1, FBXO2, FUT2, LAMC2, MMP9, NOS2, RBP4, SPEG, XCL1, FTO, FAP, LIFRA4GNT, GREM1, ODAM, CD86, GKN2, GKN1, REG1B, ZBTB16APOD, CXCL9, FUT1, NCCRP1
Digestion and lipid metabolism [12]CLC, PLAAT2, RBP4, SCT, CAPN8, GUCA2A, LCTANGPTL4, APOC3, PGC *, TFF2 *, GKN1, MUC6, LPL, LIPF, PNLIPRP2UCN2, AQP5, KCNN4, SERPINA3
Response to biotic stimulus [13]CEBPE, FCGR1B, HLA-DQA2, NOS2, XCL1, ZDHHC11, FAP, HLA-DQB1, KRT1BPIFB1 *, CCL18, LYPD8, PGC *, S100A8, CD86, FKBP5, GKN2, LPL, MUC6, PNLIPRP2, REG1BCR2, CXCL9, GBP3, MUC1, SLPI, SPINK5
Metabolic process [14]CDA, CYP4F2, CYP2S1, FUT2, LDHC, MMP7, MMP9, NOS2, PLAAT2, PLOD2, PM20D1, RBP4, SCT, SORD, SULT1C2, SLC6A8, CYP4F11, FAP, UGT2B15APOC3, ANGPTL4, HBD, HBB, HBA2, LPL, LIPF, PNLIPRP2, PSPHP1HBG1, MMP3, MMP1, TYRP1, APOD, CYP2D6, ENTPD3, FADS6, FUT1, FUT9, PSPH, SERPINA3, TNXB
Transport [15]AQP12A, SCT, SLC30A10, CYP4F2, MMP9, NOS2, PM20D1, RBP4, SLC40A1, SLC6A8, XCL1, KCNJ13APOC3, GREM1, HBD, HBB, HBA2, S100A8, SLC11A2, TFF2 *, TCN1 *, JPH4, LPL, VPS18HBG1, KCNN2, KCNK9, SCN3B, SLCO4C1, AQP5, APOD, CXCL9, GABRB3, KCNN4, KCNE2, PER1
Proteolysis [16]EGLN3, FBXO2, MMP9, MMP7, PM20D1, XPNPEP2, ERAP2, FAP, KLK12, SPINK4PHEX, PGC *, S100A8, C17orf97MMP1, MMP3, XPNPEP3, CPO, FBXL13, KCNE2, NCCRP1, SLPI, SERPINA5, SPINK5, SERPINA3
Response to xenobiotics [17]CDA, CYP4F2, CYP2S1, RBP4,NOS2, SORD, UGT2B15LPL, NAT8CYP2D6, KCNE2
Response to drug [18]CYP2S1, NOS2, SORD, SPINK4LPL, NAT8APOD, CYP2D6, KCNE2
Response to toxic substance [19] HBA2, HBB, HBDHBG1, GSTT1
Gastrointestinal cancer [20]BTNL8, CDA, LAMC2, MMP7,PLAAT2, SULT1C2, SLC6A8, TMED6, XPNPEP2, ZDHHC11, GAS5, GUCA2A, HLA-DQB1, HOXC6, LGALS2, LIFR, UGT2B15APOC3, GREM2, LCN15, LYPD8, HBB, ODAM, PGC *, S100A8, TFF2 *, TCN1 *, C6orf58, FKBP5, GKN1, GKN2, LPL, LIPF, PNLIPRP2, REG1B, RERE, UCA1, ZBTB16CR2, FDCSP, MMP1, MMP3, MYOC, XPNPEP3, TYRP1, CLDN2, C16orf89, ENTPD3, FUT9, GABRB3, HLF, KCNE2, KCNN4, LRRC17, MUC1, SERPINA3, SLPI, SPINK5, VSIG2
Upregulated DEGs are marked in bold. * DEGs in the intersection that are upregulated in MO-low-IR and downregulated in MO-high-IR relative to the control group are in bold and marked with an asterisk. a The numbers of these functions are described in Supplementary Table S2 and they include different significantly overrepresented (FDR < 0.05) gene sets.
Table
Table 3. Bibliographic analysis of differentially expressed genes (DEGs) (FDR < 0.05) that are expressed or involved in the migration, proliferation, activation, or maturation of different immune cells.
Table 3. Bibliographic analysis of differentially expressed genes (DEGs) (FDR < 0.05) that are expressed or involved in the migration, proliferation, activation, or maturation of different immune cells.
Direction of the Gene Expression in MODEGExpression or Action in Immune CellsReference
DEGs from the intersectionUpregulatedCCL18CCL18 causes maturation of cultured monocytes to macrophages in the M2 spectrum.[23]
S100A8S100A8 induces mucosal CD4+ T cells with a Th1 pro-inflammatory response.[24]
S100A8/A9 is constitutively expressed in immune and epithelial cells of inflamed tissues.[25]
DownregulatedZBTB16ZBTB16 controls the development of invariant natural killer T cell effector functions.[26]
CD86CD86 is typically found on the surface of antigen-presenting cells and can either bind CD28 or CTLA-4, resulting in a costimulatory or a co-inhibitory response, respectively.[27]
Exclusive DEGs from MO-low-IRUpregulatedXCL1XCL1 is produced mainly by NK and activated CD8+ T cells and facilitates the activation and migration of intestinal dendritic cells.[28]
ENPP3ENPP3 prevents a decrease in plasmacytoid dendritic cell numbers in the small intestine.[29]
ENPP3 is highly expressed in activated basophils and mast cells and is rapidly induced by IgE.[30]
NOS2An increase in Nos2 expression is characteristic of classically activated M1 macrophages.[31]
CLCCLC has lysophospholipase activity and is a characteristic constituent of eosinophils and basophils.[32]
CEBPECEBPE is an essential transcription factor for granulocytic differentiation.[33]
BTNL8BTNL8 is highly expressed on neutrophils.[34]
RBP4RBP4 induces antigen-presenting cells as the drivers of an inflammatory response.[35]
MMP9MMP9 recruits neutrophils to sites of inflammation.[36]
DownregulatedHLA-DQB1HLA-DQB1 is expressed in antigen presenting cells.[37]
Exclusive DEGs from MO-high-IRDownregulatedCXCL9CXCL9 is secreted in response to IFN-γ. Mainly secreted by monocytes, endothelial cells, fibroblasts, and cancer cells in response to IFN-γ[38]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ho-Plágaro, A.; Rodríguez-Díaz, C.; Santiago-Fernández, C.; López-Gómez, C.; García-Serrano, S.; Martín-Reyes, F.; Rodríguez-Pacheco, F.; Rodríguez-Cañete, A.; Alcaín-Martínez, G.; Vázquez-Pedreño, L.; et al. Morbid Obesity in Women Is Associated with an Altered Intestinal Expression of Genes Related to Cancer Risk and Immune, Defensive, and Antimicrobial Response. Biomedicines 2022, 10, 1024. https://doi.org/10.3390/biomedicines10051024

AMA Style

Ho-Plágaro A, Rodríguez-Díaz C, Santiago-Fernández C, López-Gómez C, García-Serrano S, Martín-Reyes F, Rodríguez-Pacheco F, Rodríguez-Cañete A, Alcaín-Martínez G, Vázquez-Pedreño L, et al. Morbid Obesity in Women Is Associated with an Altered Intestinal Expression of Genes Related to Cancer Risk and Immune, Defensive, and Antimicrobial Response. Biomedicines. 2022; 10(5):1024. https://doi.org/10.3390/biomedicines10051024

Chicago/Turabian Style

Ho-Plágaro, Ailec, Cristina Rodríguez-Díaz, Concepción Santiago-Fernández, Carlos López-Gómez, Sara García-Serrano, Flores Martín-Reyes, Francisca Rodríguez-Pacheco, Alberto Rodríguez-Cañete, Guillermo Alcaín-Martínez, Luis Vázquez-Pedreño, and et al. 2022. "Morbid Obesity in Women Is Associated with an Altered Intestinal Expression of Genes Related to Cancer Risk and Immune, Defensive, and Antimicrobial Response" Biomedicines 10, no. 5: 1024. https://doi.org/10.3390/biomedicines10051024

APA Style

Ho-Plágaro, A., Rodríguez-Díaz, C., Santiago-Fernández, C., López-Gómez, C., García-Serrano, S., Martín-Reyes, F., Rodríguez-Pacheco, F., Rodríguez-Cañete, A., Alcaín-Martínez, G., Vázquez-Pedreño, L., Valdés, S., Garrido-Sánchez, L., & García-Fuentes, E. (2022). Morbid Obesity in Women Is Associated with an Altered Intestinal Expression of Genes Related to Cancer Risk and Immune, Defensive, and Antimicrobial Response. Biomedicines, 10(5), 1024. https://doi.org/10.3390/biomedicines10051024

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop