Next Article in Journal
The Highly Efficient Expression System of Recombinant Human Prolidase and the Effect of N-Terminal His-Tag on the Enzyme Activity
Previous Article in Journal
The SOCE Machinery: An Unbalanced Knowledge between Left and Right Ventricular Pathophysiology
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 11
Issue 20
10.3390/cells11203283
cells-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

RNA Sequencing of Intestinal Enterocytes Pre- and Post-Roux-en-Y Gastric Bypass Reveals Alteration in Gene Expression Related to Enterocyte Differentiation, Restitution, and Obesity with Regulation by Schlafen 12

by
Emilie E. Vomhof-DeKrey
1,2,
Sonalika Singhal
3,
Sandeep K. Singhal
3,
Allie D. Stover
2,
Odele Rajpathy
2,
Elizabeth Preszler
2,
Luis Garcia
4 and
Marc D. Basson
1,2,3,*
1
Department of Surgery, School of Medicine and the Health Sciences, University of North Dakota, Grand Forks, ND 58202, USA
2
Department of Biomedical Sciences, School of Medicine and the Health Sciences, University of North Dakota, Grand Forks, ND 58202, USA
3
Department of Pathology, School of Medicine and the Health Sciences, University of North Dakota, Grand Forks, ND 58202, USA
4
Sanford Health Clinic, Sioux Falls, ND 57117, USA
*
Author to whom correspondence should be addressed.
Cells 2022, 11(20), 3283; https://doi.org/10.3390/cells11203283
Submission received: 23 September 2022 / Revised: 10 October 2022 / Accepted: 14 October 2022 / Published: 18 October 2022
(This article belongs to the Section Tissues and Organs)
Browse Figures
Versions Notes

Abstract

:
Background: The intestinal lining renews itself in a programmed fashion that can be affected by adaptation to surgical procedures such as gastric bypass. Methods: To assess adaptive mechanisms in the human intestine after Roux-en-Y gastric bypass (RYGB), we biopsied proximal jejunum at the anastomotic site during surgery to establish a baseline and endoscopically re-biopsied the same area 6–9 months after bypass for comparison. Laser microdissection was performed on pre- and post-RYGB biopsies to isolate enterocytes for RNA sequencing. Results: RNA sequencing suggested significant decreases in gene expression associated with G2/M DNA damage checkpoint regulation of the cell cycle pathway, and significant increases in gene expression associated with the CDP-diacylglycerol biosynthesis pathway TCA cycle II pathway, and pyrimidine ribonucleotide salvage pathway after RYGB. Since Schlafen 12 (SLFN12) is reported to influence enterocytic differentiation, we stained mucosa for SLFN12 and observed increased SLFN12 immunoreactivity. We investigated SLFN12 overexpression in HIEC-6 and FHs 74 Int intestinal epithelial cells and observed similar increased expression of the following genes that were also increased after RYGB: HES2, CARD9, SLC19A2, FBXW7, STXBP4, SPARCL1, and UTS. Conclusions: Our data suggest that RYGB promotes SLFN12 protein expression, cellular mechanism and replication pathways, and genes associated with differentiation and restitution (HES2, CARD9, SLC19A2), as well as obesity-related genes (FBXW7, STXBP4, SPARCL1, UTS).

1. Introduction

Patients undergoing Roux-en-Y gastric bypass (RYGB) surgery lose weight and exhibit improved glycemic and lipidemic control and other obesity-related co-morbidities [1,2,3]. Interestingly, such metabolic improvements are not entirely dependent on weight loss, as these improvements are often clinically observed before significant weight loss is even achieved [4]. It is not fully understood how bariatric surgery induces these metabolic improvements, especially since obesity-related metabolic disturbances are generated as a consequence of an interaction between lifestyle, aging, environment, and genomic and epigenetic factors [1]. Additionally, the RYGB method alters gastrointestinal continuity [5]. Within the jejunal segment that is positioned at the proximal end of the Roux-limb, the intraluminal milieu is substantially altered by its more proximal location in the chain of oral intake. This results in an augmentation in microbiota and mechanical load as well as the chemical effects of undigested food [6]. However, little is known about how enterocyte gene expression is affected by RYGB.
We previously studied the role of SLFN12 (and its mouse ortholog, Slfn3) in the regulation of enterocyte differentiation [7,8,9,10,11,12,13]. SLFN12 is one of a family of proteins that also regulate hematopoiesis, the immune response, and bone biology [14,15,16,17]. SLFN proteins are classified into three groups based on size. Group I, long SLFN proteins, have a helicase domain with a nuclear targeting sequence that allows for nuclear DNA binding and then modulation of transcription. However, intermediate Group II SLFN proteins such as SLFN12 and the short Group III SLFN proteins lack this nuclear targeting sequence [16]. Instead, SLFN12 influences enterocyte differentiation through a pathway involving SLFN12, serpinB12, and the deubiquitylases UCHL5 and USP14 that in turn alter protein levels of the transcription factor CDX2 [8,9]. In this study, we collected operative RYGB proximal jejunum tissue at the anastomotic site and endoscopically biopsied the same area 6–9 months after RYGB for comparison. After sectioning and staining, laser capture dissection isolated enterocytes along the villi and RNASeq were performed on RNA isolated from the laser capture dissected enterocytes. We stained for SLFN12 to evaluate the effects of RYGB on enterocytic SLFN12 and then investigated whether SLFN12 might influence some of the RNASeq genes that were affected by RYGB.

2. Materials and Methods

2.1. Human Samples

Human samples and non-identifying patient information were obtained under a protocol approved under IRB-201606-406 from the University of North Dakota in Grand Forks, ND, USA and by the IRB of the Sanford Heath System in Fargo, ND, USA. Patients who were already scheduled for a Roux-en-Y gastric bypass surgery for clinical indications at Sanford Health were recruited for participation in this research study. The study population thus consisted of 11 consecutive patients scheduled for RYGB who had consented to participation in the research. During the scheduled RYGB surgery, the most proximal 2 cm portion of the Roux limb (proximal jejunum) was collected and divided into 2 mL tubes for either direct freezing in liquid nitrogen, or storage in DNA/RNA shield solution (Zymo Research, Irvine, CA, USA) or 10% formalin (Sigma-Aldrich, St. Louis, MO, USA). Six to nine months after RYGB, each patient was scheduled for an upper endoscopy. Four to six endoscopic mucosal biopsies were obtained from the jejunum. Each biopsy consisted of a piece of mucosa that was approximately 0.4 cm × 0.2 cm × 0.1 cm. Biopsies were collected into 2 mL tubes and preserved in a similar way to the RYGB samples above. After all the samples had been collected and deidentified with a subject ID number, they were transported to the research laboratory at the University of North Dakota together with corresponding deidentified health information for each research subject including age, sex, race, medical conditions, current medications, weight, and BMI prior to surgery and at the time of endoscopy.

2.2. Tissue Sectioning, Staining, and Laser Capture Microscopy

Frozen tissue was embedded in an optimal cutting temperature compound and 8 µm sections were frozen sectioned on a microtome-cryostat. Sections were mounted onto PET membrane steel frame slides (1.4 µm, RNase-free, Leica LMD slides; NCI, Inc. Brooklyn Park, MN, USA) and sectioned slides were stored at −80 °C. Applied 50 µL of RNase Inhibitor (ProtectRNA 500x R7397, Sigma-Aldrich, St. Louis, MO, USA) was added to all staining jars (per 25 mL of the following ethanol and isopropanol jars, not in a xylene jar). Frozen slides were immediately submerged into 95% ethanol for 30 s. with frequent dipping up and down. Hematoxylin stain was directly applied to the tissue sections on the slides for 60 s. These were then de-stained in 100% isopropanol, via 10 dips for 15 s each. Eosin stain was applied directly to tissue sections on the slide for 2 s and they were then de-stained in 100% isopropanol for 10 dips for 30 s each. The slides were subsequently de-stained in 75% ethanol for five dips for 15 s, then 95% ethanol for 10 dips, then 100% ethanol for 10 dips for 10 s. Finally, the slides were submerged in anhydrous 100% ethanol for 30 s. Then the slides were repeatedly dipped in xylene for 15 s with a final submersion for 10 min. Each slide was air-dried for 2 min. Laser capture microscopy immediately followed in which two sections of about 10–30 intestinal epithelial cells along with the villi were cut and collected into the collection cap (Leica LMD 6, Deer Park, IL, USA). Sectioned cells were lysed with 5 µL of lysis buffer using the Tecan LCMSeq kit and tubes were stored at −80 °C. The RNA library was prepared using the Tecan Ovation SoLo RNA-Seq library preparation kit as per the manufacturer’s instructions. Libraries were quantified using the Quant-iT PicoGreen dsDNA assay kit and read using a BioTek FLx800 fluorescence microplate reader. Libraries were pooled at equal molarity and the sample pool was quantified using Biorad ddPCR. Libraries were sequenced using an Illumina HiSeq X Ten at Psomagen (Rockville, MD, USA).

2.3. Immunohistochemical Staining for Slfn12 Protein

A Leica Bondmax Stainer in modified protocol F was used to stain for Slfn12 protein in the gastric bypass and EGD tissue sections. Paraffin sections were dewaxed and antigen retrieval was performed with Leica Bond Epitope Retrieval Solution 1 for 20 min. For immunostaining, slides were sequentially washed with bond wash buffer (BWB), incubated with SLFN12 antibody (Rabbit polyclonal, LS-B4757, 1:200 dilution; LSBiosystems, Seattle, WA, USA) for 30 min, bonded after the primary antibody for 20 min, washed with BWB, bonded with polymer for 20 min, washed with BWB, bonded with peroxidase block for 10 min, washed with BWB, bonded with DAB for 10 min, washed with deionized water, bonded with hematoxylin for 5 min, and washed with BWB. The slide was dehydrated, and a coverslip was applied. Slides were then deidentified as to the peri-operative or post-operative nature of the specimens. Two blinded observers were asked to score the intensity of epithelial SLFN12 staining on a 0 to 3 scale in each image based upon a predetermined scale. These scores were then aggregated and subjected to statistical analysis by a Student t-test in a GraphPad Prism 8.

2.4. Cell Culture and RNA Expression Analysis

Human HIEC-6 cells (ATCC, Manassas, VA, USA) were cultured in a base medium of OptiMEM 1 Reduced Serum (Gibco, Thermofisher, Waltham, MS, USA), 20 mM HEPES (Thermofisher, Waltham, MS, USA), 10 mM GlutaMAX (Thermofisher, Waltham, MS, USA), 10 ng/mL epidermal growth factor (EGF, Biolegend, San Diego, CA, USA) and 4% fetal bovine serum (FBS, Genesee Scientific, El Cajon, CA, USA). FHS 74 Int cells (ATCC, Manassas, VA, USA) were cultured in a base medium of Hybri-Care ATCC 46-X, 30 ng/mL EGF, and 10% FBS. Both cell lines were incubated at 37 °C with 5% CO2. Cells in 6-well plates with 300,000 cells/well were incubated with 4000 vp/cell Slfn12 adenovirus (Adv-Slfn12) or Adv-CMV for 2 h. Media was replaced and the cells harvested 72 h later with 600 µL of RLT buffer. RNA was isolated using the RNeasy Mini Kit and Qiacube as per the manufacturer’s protocols (Qiagen, Valencia, CA, USA). A SmartScribe RT kit was used for cDNA synthesis (Takara Clontech, Mountain View, CA, USA). The BioRad CFX96 Touch Real-Time PCR Detection System and the PrimeTime Gene Expression Master Mix (IDT, Coralville, IA, USA) were utilized for qPCR analysis. Expression levels were ascertained from the threshold cycle (Ct) values using the 2−∆∆Ct method, where RPLP0 was used as a housekeeping gene. The following primer/probe sets were used from BioRad (Hercules, CA, USA) and are proprietary: RPLP0 (Assay ID: qHsaCEP0041375, HEX), SPARCL1 (Assay ID: qHsaCEP0057695, TEX615), SLC19A2 (Assay ID: qHsaCIP0031109, Cy5.5), MBOAT2 (Assay ID: qHsaCEP0052422, Cy5.5), and FBXW7 (Assay ID: qHsaCIP0031037). Primer/probe sets from IDT were as follows: CARD9 forward 5′-GA CTG GAA GAT GGC TCA C-3′, reverse 5′-AGC TGC AAA GGG CTG TT-3′, probe 5′-/56-FAM/TTT GTC TGA/ZEN/GAG CTG GGT GTC CTC/3IABkFQ/-3′; STXBP4 forward 5′-CCT AAA TCC CTC TGT TCG CTT-3′, reverse 5′-AGT CTG CTT GTA CTT GCT GTC-3′, probe 5′-/56-FAM/TTG GTG TTG/ZEN/TTC CTT TGT GGG CTG/3IABkFQ/-3′; UTS2 forward 5′-CCT TCT ACA GAT ACT GCC AGA-3′, reverse 5′-AAC TTT CTC AAA TTT CCT CTT GGG-3′, probe 5′-/56-FAM/ATA TCC CCT/ZEN/CTT TCT GCA CCC AGC/3IABkFQ/-3′; HES2 forward 5′-CTG CTC GAA GCT AGA GAA GG-3′, reverse 5′-GTA GCT GTC GCA AGG CA-3′, probe 5′-/5Cy5/CGC ACG GTC ATT TCC AGG ACG T/3IAbRQSp/-3′; TRPC5 forward 5′-TCG ATG AGC CTA ACA ACT GC-3′, reverse 5-CCA AAT ACA GAC CAG AAG AGT GA-3′, probe 5′-/5Cy5/TCA AAG AGC GTG GAG AAG GCA TTG T/3IAbRQSp/-3′. qPCR cycle conditions were 1 cycle of 2 min at 95 °C, 50 cycles of 10 s at 95 °C and 45 s at the annealing temperature of 55 °C.

2.5. Processing of Sequence Data

RNA-Seq raw reads were stored in FASTQ format files containing the sequence data for this cohort study (11 patients). The quality of each raw FASTQ file was explored using fastQC [18]. As per the quality control report on the raw data, files were subjected to pre-processing steps including trimming, quality filtering, per-read quality pruning using a fastp [19] tool. After the data had been preprocessed, the FASTQ file reads were then aligned to a reference genome (Homo sapiens GRCh38, available on ENSEMBL [20]) using a STAR aligner v.2.5.2b to generate the read counts. The normalized distribution of the counts was checked to identify the outliers.

2.6. Statistical Analysis

Two sub-selections of the samples were made to classify pre-surgery and post-surgery samples and to compare the effects of gastric bypass on the gene profiles of patients. Principal component analysis (PCA) was performed on the RNA-Seq to describe the relationship between the samples by visualizing the samples and groups in two-dimensional space. The statistically significant genes that were differentially expressed in the sub-selections of the samples were identified using a two-tailed Student’s t-test. The significance level of the p-value < 0.05 was employed as a standard to filter genes. Further pathway and function enrichment analyses of the significant genes were performed using the Ingenuity Pathway Analysis (IPA, Ingenuity Systems; Qiagen, Valencia, CA, USA; www.qiagen.com/ingenuity, accessed on 4 June 2021) to identify the biological functions, canonical pathways, and regulatory networks of the functional mRNA target interactions. Comma delimited text files containing gene IDs, expression data (fold change), and p-values were uploaded into the IPA for core analysis. The top functions, pathways and gene-signaling networks were calculated using IPA-generated negative logarithm p-values i.e., −log10(p-value) and associated Z- and network scores. The investigation into the outcome of the surgery was observed by looking at the change in BMI (body mass index) of the patients. To differentiate between the gene expression among patients who were successful in weight loss and those who did not lose substantial weight, the change in BMI as a linear predictor was taken for post-surgery samples.

3. Results

3.1. Comparison between Pre- and Post-Surgery Gene Expression

The distribution of data for the gene expression levels before and after surgery was observed through the application of principal component analysis (PCA). Figure 1A shows the variation between pre- and post-RYGB samples as a group and individuals using the first two principal components. The overall gene expression of pre- and post- surgery samples did not form distinct clusters, suggesting that not all patients conformed to a similar gene expression pattern. The total variation observed between the groups was 23.4% to 6.4% for principal component-1 (PC1) and principal component-2 (PC2), respectively. The paired samples ranged in variations from very low (such as GBP11-EGD11, GBP15-EGD15) to very high (such as GBP02-EGD02, GBP03--EGD03, seen as outliers) between pre- and post-surgery on the PC1 and PC2 scales.
Differential expression analysis between Post-RYGB and Pre-RYGB samples identified 1116 (916 upregulated and 200 downregulated) genes that changed significantly to meet the threshold p < 0.05 (Table S1). To gain further insights into the functional association of these genes, all significant genes were investigated by IPA core analysis. This analysis identified eight pathways (Figure 1B) that attained the desired level of significance (p < 0.05). The interaction between these pathways suggests two clusters in which 14 pathways are functionally interlinked. The human pathways significantly enriched were Cell Cycle: G2/M DNA damage Checkpoint Regulation (7.66 × 10−3), pyridoxal 5′ phosphate salvage (7.66 × 10−3), salvage pathways of pyrimidine ribonucleotides (1.47 × 10−2), and interferon signaling (1.42 × 10−2) based on their significance value (p-value ≤ 0.05) (Figure 1C).

3.2. Changes in the Gene Profiles for Prominent Responders of Gastric Bypass

All patients lost weight after surgery, as evidenced by a substantially decreased body mass index (BMI (Figure 2). We compared the gene expression level with the BMI to find the genes most significantly associated with weight loss. Three different approaches were applied to find the most significant genes associated with the pre- and post-RYGB and BMI changes together. First, we identified the genes changing in the samples after RYGB with respect to the change in BMI (pre- vs. post-RYGB) using a regression analysis. A total of 1750 genes changed significantly with p < 0.05 (Table S2). We then compared these 1750 genes (Table S2) with the 1116 genes identified between the pre- and the post-RYGB in Table S1. There were 74 genes in common between the two analyses (Figure 3). Second, to test whether there was any significant statistical relationship between the BMI and pre/post-RYGB together with the changes in gene expression, we introduced an interaction term to the regression model. We found 1200 genes that had changed significantly with the combined effect of RYGB and BMI (Table S3). Out of these 1200 genes, 540 genes were upregulated and 660 were downregulated. Third, in another approach, we identified the genes whose expression changed significantly (in terms of fold change) with respect to changes in percentages of BMI when measured pre- and post-the RYGB condition. A total of 867 genes were highly correlated (Spearman correlation coefficient > 0.70) with 789 upregulated and 78 downregulated (Table S4). Figure 4 shows the changes in 19 genes from this list which significantly changed their expression with respect to the change in BMI (Table S5).

3.3. SLFN12 Protein Expression Is Increased after Gastric Bypass

Human enterocytic differentiation is regulated by SLFN12 via a pathway involving SerpinB12, the deubiquitylases UCHL5 and USP14, and the transcription factor, Cdx2 [9]. Indeed, SLFN12 protein immuno-reactivity was clearly greater after surgery as compared to the preoperative baseline (typical images in Figure 5a, right panel) and statistical analysis of scoring by blinded observers confirmed this difference (Figure 5a, left panel). The SLFN12 protein immunoreactivity average intensity values were subsequently compared with the pre- and post-operative BMIs (Figure 5B). Furthermore, the change in SLFN12 protein intensity after surgery in individual patients had a trending correlation that could indicate an increase in SLFN12 expression, with a decrease in BMI after surgery (Figure 5C).

3.4. Similar RNA Expression Change in Gastric Bypass-Specific Genes Compared to HIEC-6 or FHs 74 Int Cells Overexpressing SLFN12

To investigate whether this higher SLFN12 expression after RYGB could contribute to some of the post-operative gene changes that we observed, we overexpressed SLFN12 with an adenovirus (AdvSLFN12) in two human intestinal cell lines, HIEC-6 and FHs 74 Int (Figure 6A). The following genes were selected for study from among those genes identified by the RNAseq analysis (Table S1 and Figure 3) because the literature suggests that they are involved in the differentiation and restitution of intestinal mucosa: HES2, CARD9, TRPC5, and SLC19A2 [21,22,23,24,25]. We also analyzed the expression of FBXW7, MBOAT2, STXBP4, SPARCL1, and UTS after SLFN12 overexpression because these are obesity-related genes [26,27,28,29,30,31,32]. HES2 and SPARCL1 were increased in RNA expression in the HIEC-6 + AdvSLFN12 cells compared to the CMV control. However, the FHs 74 Int + AdvSLFN12 treated cells only had a modest but not significant increase (Figure 6B,I). CARD9 was only increased in the FHs 74 Int + AdvSLFN12 treated cells compared to the CMV control (Figure 6C). SLC19A2, FBXW7, STXBP4, and UTS2 all increased in RNA expression in both the HIEC-6 and FHs 74 Int + AdvSLFN12 treated cells in comparison to controls (Figure 6D,F,H,J), while MBOAT2 RNA expression decreased (Figure 6G). Interestingly, TRPC5 decreased in the HIEC-6 + AdvSLFN12 treated cells, while in the FHs 74 Int + AdvSLFN12 treated cells the TRPC5 RNA expression increased in comparison to the CMV control (Figure 6E). The changes observed in the genes above were compared to the RNA Seq expression data in Table 1. Most of the genes that we studied had similar increases or decreases in expression between the RNA Seq data and the AdvSLFN12 overexpressed intestinal cell lines, except for HES2, TRPC5, and SPARCL1 in the FHs 74 Int AdvSLFN12 treated cells, CARD9 in the HIEC-6 AdvSLFN12 treated cells, and MBOAT2, which decreased in the intestinal AdvSLFN12-treated cells compared to the RNA Seq data indicating an increase in MBOAT2 expression (Table 1).

4. Discussion

Instead of collecting RNA from whole intestinal mucosal tissue, we performed laser capture microscopy to collect intestinal epithelial cells only from the villi of the proximal jejunum biopsies from the Roux-en-Y gastric bypass surgery and then the EGD biopsies 6–9 months post-operatively. Patients served as their own controls for a paired analysis when examining the genes in laser-dissected enterocytes that had changed after RYGB. RNA sequencing identified 1116 statistically significant genes and pathways that changed expression in the enterocytes after RYGB (Table S1, Figure 1B). The most statistically significant canonical pathways that were affected were the cell cycle G2/M DNA damage checkpoint regulation, the CDP-diacylglycerol biosynthesis, the phosphatidylglycerol biosynthesis, the TCA cycle, and the pyrimidine ribonucleotide interconversion (derived from Figure 1 and summarized in Table 2). Jorsal et al. performed mRNA sequencing on 20 patients’ intestinal mucosal biopsies that had been collected 3 months prior to RYGB surgery and 3 months after RYGB surgery [3]. They identified 2437 differentially expressed genes after RYGB and then focused on the specific expression changes of 16 prohormones of enteroendocrine or neuroendocrine peptides and 39 G-protein coupled receptor genes [3]. They identified 16 robustly upregulated or downregulated prohormone genes. Notably, we did not observe a change in gene expression for these same genes in our study [3]. Additionally, Jorsal stated that 39 G-protein coupled receptors GPCRs had differential expression (up- or down-regulation was not noted) after RYGB [3]. Two GPCRs that were differentially expressed in both their study and ours were ADORA1 and S1PR3 [3]. We also observed expression changes in the following related but different GPCRs: GPR158, ADGRA3, ADGRG3, and ADRM1. Another study by Ben-Zvi et al. analyzed jejunal mucosa biopsies, from at the time of RYGB surgery and 1 month after RYGB [4]. They observed expression changes related to glucose, fatty acid, and cholesterol metabolism [4]. Although we did not find the same genes to change expression in our study similar to those reported by Ben-Zyi and colleagues, Ben-Zyi did observe an increase in the Paneth cell marker gene DEFA5 and the tuft cell marker gene POU2F3, while we also observed an increase in the related genes DEFB112 and POU4F2. Notably, these previous RNA sequencing studies were performed on heterogeneous mucosal populations of intestinal epithelial, immune and other cells within the mucosa, whereas our laser capture microscopy focused our analysis on the effects of RYGB, specifically on intestinal enterocytes, avoiding potential obfuscation by changes in gene expression in other cell types in the mucosa. In addition, each of these studies analyzed intestinal mucosal gene expression relatively early after surgery, while our study focused on a longer time interval by which a more steady state would have been expected to be achieved. This could have additionally contributed to differences between our observations and those of these other two studies.
The most statistically significant changing canonical pathway (Table 2), the cell cycle G2/M DNA damage checkpoint regulation, exhibited decreases in the AURKA, CKS1B, CKS2, and PLK1 genes after RYGB. The AURKA gene codes for the Aurora A protein which is involved in the disassembly of primary cilia necessary for the cell to enter into the cell cycle [33]. Dysregulation of cilium assembly/disassembly has been linked to a risk of obesity [34,35,36]. Cyclin-dependent kinases (CKs) (i.e., CKS1B, CKS2, Table S1) are also involved with cell cycle regulation, and inhibitors of CKs have drawn attention to their use in controlling diet-induced obesity [37]. Finally, Plk1 mRNA expression was upregulated in the visceral adipose tissue of obese non-diabetic patients compared to both control and obese diabetic patients [38].
The CDP-diacylglycerol biosynthesis pathway displayed increases in the SDHB and SDHC genes after RYGB surgery. CDP-diacylglycerol synthases are involved in regulating the growth in lipid droplets, fat storage, and adipocyte development [39,40]. The loss of SDHC in TH+ glomus cells leads to increased obesity in a mouse model. [41] However, SDHB-deficient mice were protected from obesity induced by a high fat diet [42]. These data and ours indicate that CDP-diacylglycerol synthases have a complex role in obesity that may be affected by diet.
After RYGB, the TCA cycle II pathway included increases in FAM20B and MAP3K9. FAM20B encodes for glycosaminoglycan xylosylkinase that is involved in the TCA cycle II pathway and is also involved in proteoglycan synthesis [43]. FAM20B was linked to obesity in a weight-discordant monozygotic twin study, where the adipose tissue matrisome had a positive association to the percent of liver fat, circulating triglycerides, and homeostatic model assessment of insulin resistance, each of which may contribute to metabolic dysfunction [43]. Single-nucleotide polymorphisms (SNPs) in MAP3K9 have been shown to interact with pre-menopausal obesity and diabetes [44]. Additionally, MAP3K9 is known to have a key role in regulating the JNK pathway, and is involved in the development of obesity [45].
Finally, the pyrimidine ribonucleotides interconversion pathway was characterized by increases in IFI35, IFNAR2, MED14, and TYK2. An interesting study involving human maternal-diet effect on offspring obesity examined the SNP interaction effect on the expression of various genes using DNA isolated from the blood buffy coat of offspring born before the mothers underwent gastrointestinal bypass surgery (BMS) or those born after bypass surgery (AMS) [46]. SNP rs35447805 had the most significant interaction with treatment for expression level of IFI35 where the IFI35 expression levels were a 0.94-fold decrease in expression for BMS heterozygous carriers and a 1.42-fold increase in expression for AMS heterozygous carriers [46]. Another study found IFI35 to have a decrease in expression from the subcutaneous adipose tissue of obese pigs [47]. MED14 has been implicated in the regulation of lipid homeostasis and adipogenesis via its interactions with the glucocorticoid receptor and PPARγ [48,49,50]. Finally, TYK2 has a role in the differentiation of Myf5+ brown adipocytes and when Tyk2 is lost in Tyk2−/− mice they become obese with age [51]. Furthermore, diet-induced obese mice and obese humans have decreased Tyk2 expression; this can be reversed by constitutively active Stat3, which restores BAT differentiation and energy balance [51].
Next, we used weight loss as a measure of how well the patients responded to the gastric bypass treatment. We compared gene expression differences to weight loss based on changes in patient BMI (Figure 3 and Tables S2–S5). Additionally, we observed that the human enterocytic differentiation regulator, SLFN12, increased in protein expression in the post-gastric bypass endoscopic biopsies. We previously showed that the expression of SLFN12 immunostaining expression mirrors the mRNA expression in control vs. fasting duodenal human tissue [9] and that immunostaining for rat Slfn3 (orthologous to human SLFN12) mirrors Western blot-assayed protein and mRNA expression in rat intestinal tissue [8]. A study involving an in vitro human ulcerative colitis colon organoid model found that after 60 weeks of inflammatory stimulation the SLFN11 mRNA and protein expression was increased [52]. Additionally, the investigators described higher SLFN11 mRNA and protein expression in human endoscopic biopsy samples from patients with ulcerative colitis in comparison to endoscopic biopsies from healthy patients [52]. Given that little is known about the mechanisms that control SLFN protein expression, it is possible that similar mechanisms may influence GI mucosal changes in SLFN12 and SLFN11 in various disease states. Further analysis was performed to determine whether SLFN12 may regulate specific genes involved with differentiation and restitution of intestinal mucosa and obesity-related genes from the significant 74 genes observed above (Figure 3). The overexpression of SLFN12 in human HIEC-6 and/or FHs 74 Int cells led to a similar augmentation in the expression of the following differentiation and restitution genes: HES2, CARD9, and SLC19A2 and the following obesity-related genes: FBXW7, STXBP4, SPARCL1, and UTS in comparison with the increased expression observed in the RNA Seq data. The few genes that did not change similarly (HES2, TRPC5, SPARCL1 in FHs 74 Int, CARD9 in HIEC-6, and MBOAT5) between the RYGB RNA Seq and the HIEC-6/FHs 74 Int cell lines overexpressing SLFN12 could reflect differences in the biology of the individual cell lines or differences in the interaction between SLFN12 intracellular effects and those of other possible factors found in vivo but not in vitro (such as neuroendocrine [53], hormonal [54,55], growth factors [56,57,58,59], nutrients [60,61,62], extracellular matrix proteins [56,57,58], and physical forces [61,63,64,65,66]). Furthermore, since little is known about the factors that regulate SLFN12 expression itself, it is possible that some of the differences that we observed in vivo but not in vitro could have been upstream of SLFN12 and thus changed in association with, but not because of, SLFN12. However, most of the genes in Table 1 did change similarly between the RYGB RNA Seq and cell culture data, suggesting that SLFN12 may contribute to these gene changes that we observed after RYGB.
While little is known about how SLFN12 itself is regulated, we have previously shown that the rodent ortholog Slfn3 is predominantly cytosolic and modestly nuclear in localization with exogenously expressed Slfn3 in IEC-6 cells. We then transfected Slfn3-null Caco2 cells with a Slfn3 plasmid containing a nuclear exclusion sequence and found that Slfn3 was confined to the cytoplasm and was still able to induce villin promoter activity [7]. SLFN12 itself is also predominantly cytosolic and only modestly nuclear in localization when exogenously adenovirally expressed in Caco-2 cells. [9] Indeed, and unlike long SLFNs such as SLFN11, SLFN12 modulates enterocytic differentiation in the cytosol by affecting the deubiquitylation and proteasomal degradation of the transcription factor CDX2, thus altering levels of newly synthesized CDX2 levels in the cytosol before that transcription factor moves back into the nucleus to bind to the chromatin [9].
Finally, with regard to genes which changed significantly in expression with respect to BMI (Figure 4, Table S5), these genes might reflect or contribute to the efficacy of the RYGB surgery. Ten of these genes are yet to be characterized and the other nine genes do not appear to be related to obesity or intestinal differentiation and restitution based on current literature review. Therefore, future studies would be warranted to determine the role of these 19 genes in adaptive mechanisms after RYGB and weight loss.

5. Conclusions

This study investigated the genes that specifically change only in intestinal enterocytes after human RYGB. RNA sequencing yielded canonical pathways that were affected favorably for obesity and its comorbidities. Correlating gene expression changes to the patient BMI, we were able to define gene sets associated with RYGB effectiveness. Finally, we demonstrated an increase in SLFN12 protein after RYGB and that overexpression of SLFN12 in intestinal cells leads to similar gene expression changes that were altered after RYGB.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cells11203283/s1. Table S1: Genes that changed significantly with differential expression analysis between Post-RYGB and Pre-RYGP samples. Table S2: Genes that changed significantly for post-RYGB samples with respect to the change in BMI using a regression analysis. Table S3: Genes that changed significantly with the combined effect of RYGB and BMI. Table S4: Genes that highly correlated with respect to changes in percentages of BMI when measured pre- and post- the RYGB condition. Table S5: 19 genes that changed their expression significantly from Table S4 with respect to the change in BMI.

Author Contributions

Conceptualization, E.E.V.-D. and M.D.B.; methodology, E.E.V.-D., L.G. and M.D.B.; software, S.S. and S.K.S.; formal analysis, E.E.V.-D., S.S. and S.K.S.; investigation, E.E.V.-D., A.D.S., O.R., E.P. and L.G.; resources, M.D.B.; data curation, E.E.V.-D., S.S. and S.K.S.; writing—original draft preparation, E.E.V.-D., S.S. and S.K.S.; writing—review and editing, E.E.V.-D., S.S., S.K.S., L.G. and M.D.B.; visualization, E.E.V.-D., S.S. and S.K.S.; supervision, M.D.B.; funding acquisition, M.D.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Institutes of Health (NIH) Grant RO1-DK-096137, M.D.B. and UNDSMHS funds. Genomics Core Research reported in this publication was supported by the National Institute of General Medical Sciences of the National Institutes of Health under Award Number U54GM128729 and Award number 2P20GM104360.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board of the University of North Dakota and Sanford Health System (protocol code IRB-201606-406 and date of approval 04-07-2017) for studies involving humans.

Informed Consent Statement

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

Data Availability Statement

RNA sequencing data can be found at the following GEO link: https://www.ncbi.nlm.nih.gov/geo/info/seq.html (accessed on 1 November 2022).

Acknowledgments

The authors would like to thank the Pathology department for their assistance in the SLFN12 histology staining.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fraszczyk, E.; Luijten, M.; Spijkerman, A.M.W.; Snieder, H.; Wackers, P.F.K.; Bloks, V.W.; Nicoletti, C.F.; Nonino, C.B.; Crujeiras, A.B.; Buurman, W.A.; et al. The effects of bariatric surgery on clinical profile, DNA methylation, and ageing in severely obese patients. Clin. Epigenet. 2020, 12, 14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Jakobsen, G.S.; Smastuen, M.C.; Sandbu, R.; Nordstrand, N.; Hofso, D.; Lindberg, M.; Hertel, J.K.; Hjelmesaeth, J. Association of Bariatric Surgery vs Medical Obesity Treatment With Long-term Medical Complications and Obesity-Related Comorbidities. JAMA 2018, 319, 291–301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Jorsal, T.; Christensen, M.M.; Mortensen, B.; Nygaard, E.B.; Zhang, C.; Rigbolt, K.; Wandall, E.; Langholz, E.; Friis, S.; Worm, D.; et al. Gut Mucosal Gene Expression and Metabolic Changes After Roux-en-Y Gastric Bypass Surgery. Obesity (Silver Spring) 2020, 28, 2163–2174. [Google Scholar] [CrossRef] [PubMed]
  4. Ben-Zvi, D.; Meoli, L.; Abidi, W.M.; Nestoridi, E.; Panciotti, C.; Castillo, E.; Pizarro, P.; Shirley, E.; Gourash, W.F.; Thompson, C.C.; et al. Time-Dependent Molecular Responses Differ between Gastric Bypass and Dieting but Are Conserved Across Species. Cell Metab. 2018, 28, 310–323.e316. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Spak, E.; Bjorklund, P.; Helander, H.F.; Vieth, M.; Olbers, T.; Casselbrant, A.; Lonroth, H.; Fandriks, L. Changes in the mucosa of the Roux-limb after gastric bypass surgery. Histopathology 2010, 57, 680–688. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Sundbom, M.; Elphick, D.A.; Mahida, Y.R.; Cunliffe, R.N.; Midtvedt, T.; Engstrand, L.; Rubio, C.; Axelsson, L.G. Alteration in human defensin-5 expression following gastric bypass surgery. J. Clin. Pathol. 2007, 60, 1029–1034. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Chaturvedi, L.; Sun, K.; Walsh, M.F.; Kuhn, L.A.; Basson, M.D. The P-loop region of Schlafen 3 acts within the cytosol to induce differentiation of human Caco-2 intestinal epithelial cells. Biochim. Biophys. Acta 2014, 1843, 3029–3037. [Google Scholar] [CrossRef] [Green Version]
  8. Kovalenko, P.L.; Yuan, L.; Sun, K.; Kunovska, L.; Seregin, S.; Amalfitano, A.; Basson, M.D. Regulation of epithelial differentiation in rat intestine by intraluminal delivery of an adenoviral vector or silencing RNA coding for Schlafen 3. PLoS ONE 2013, 8, e79745. [Google Scholar] [CrossRef] [Green Version]
  9. Basson, M.D.; Wang, Q.; Chaturvedi, L.S.; More, S.; Vomhof-DeKrey, E.E.; Al-Marsoummi, S.; Sun, K.; Kuhn, L.A.; Kovalenko, P.; Kiupel, M. Schlafen 12 Interaction with SerpinB12 and Deubiquitylases Drives Human Enterocyte Differentiation. Cell Physiol. Biochem. 2018, 48, 1274–1290. [Google Scholar] [CrossRef]
  10. Vomhof-DeKrey, E.E.; Umthun, J.; Basson, M.D. Loss of Schlafen3 influences the expression levels of Schlafen family members in ileum, thymus, and spleen tissue. PeerJ 2020, 8, e8461. [Google Scholar] [CrossRef]
  11. Vomhof-DeKrey, E.E.; Stover, A.D.; Labuhn, M.; Osman, M.R.; Basson, M.D. Vil-Cre specific Schlafen 3 knockout mice exhibit sex-specific differences in intestinal differentiation markers and Schlafen family members expression levels. PLoS ONE 2021, 16, e0259195. [Google Scholar] [CrossRef] [PubMed]
  12. Vomhof-DeKrey, E.E.; Lee, J.; Lansing, J.; Brown, C.; Darland, D.; Basson, M.D. Schlafen 3 knockout mice display gender-specific differences in weight gain, food efficiency, and expression of markers of intestinal epithelial differentiation, metabolism, and immune cell function. PLoS ONE 2019, 14, e0219267. [Google Scholar] [CrossRef] [PubMed]
  13. Vomhof-DeKrey, E.E.; Lansing, J.T.; Darland, D.C.; Umthun, J.; Stover, A.D.; Brown, C.; Basson, M.D. Loss of Slfn3 induces a sex-dependent repair vulnerability after 50% bowel resection. Am. J. Physiol. Gastrointest. Liver Physiol. 2021, 320, G136–G152. [Google Scholar] [CrossRef] [PubMed]
  14. Fletcher, S.J.; Johnson, B.; Lowe, G.C.; Bem, D.; Drake, S.; Lordkipanidze, M.; Guiu, I.S.; Dawood, B.; Rivera, J.; Simpson, M.A.; et al. SLFN14 mutations underlie thrombocytopenia with excessive bleeding and platelet secretion defects. J. Clin. Investig. 2015, 125, 3600–3605. [Google Scholar] [CrossRef] [Green Version]
  15. Lee, N.K.; Choi, H.K.; Yoo, H.J.; Shin, J.; Lee, S.Y. RANKL-induced schlafen2 is a positive regulator of osteoclastogenesis. Cell Signal. 2008, 20, 2302–2308. [Google Scholar] [CrossRef]
  16. Mavrommatis, E.; Fish, E.N.; Platanias, L.C. The schlafen family of proteins and their regulation by interferons. J. Interferon. Cytokine Res. 2013, 33, 206–210. [Google Scholar] [CrossRef] [Green Version]
  17. Puck, A.; Aigner, R.; Modak, M.; Cejka, P.; Blaas, D.; Stockl, J. Expression and regulation of Schlafen (SLFN) family members in primary human monocytes, monocyte-derived dendritic cells and T cells. Results Immunol. 2015, 5, 23–32. [Google Scholar] [CrossRef] [Green Version]
  18. Andrews, S.; Krueger, F.; Seconds-Pichon, A.; Biggins, F.; Wingett, S. FastQC: A quality control tool for high throughput sequence data. In Babraham Bioinformatics; Babraham Institute: Cambridge, UK, 2015. [Google Scholar]
  19. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. Fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef]
  20. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  21. Nabokina, S.M.; Reidling, J.C.; Said, H.M. Differentiation-dependent up-regulation of intestinal thiamin uptake: Cellular and molecular mechanisms. J. Biol. Chem. 2005, 280, 32676–32682. [Google Scholar] [CrossRef]
  22. Rao, J.N.; Platoshyn, O.; Li, L.; Guo, X.; Golovina, V.A.; Yuan, J.X.; Wang, J.Y. Activation of K(+) channels and increased migration of differentiated intestinal epithelial cells after wounding. Am. J. Physiol. Cell Physiol. 2002, 282, C885–C898. [Google Scholar] [CrossRef] [Green Version]
  23. Sokol, H.; Conway, K.L.; Zhang, M.; Choi, M.; Morin, B.; Cao, Z.; Villablanca, E.J.; Li, C.; Wijmenga, C.; Yun, S.H.; et al. Card9 mediates intestinal epithelial cell restitution, T-helper 17 responses, and control of bacterial infection in mice. Gastroenterology 2013, 145, 591–601.e593. [Google Scholar] [CrossRef] [Green Version]
  24. Cao, L.; Qin, X.; Peterson, M.R.; Haller, S.E.; Wilson, K.A.; Hu, N.; Lin, X.; Nair, S.; Ren, J.; He, G. CARD9 knockout ameliorates myocardial dysfunction associated with high fat diet-induced obesity. J. Mol. Cell Cardiol. 2016, 92, 185–195. [Google Scholar] [CrossRef] [Green Version]
  25. Oittinen, M.; Popp, A.; Kurppa, K.; Lindfors, K.; Maki, M.; Kaikkonen, M.U.; Viiri, K. Polycomb Repressive Complex 2 Enacts Wnt Signaling in Intestinal Homeostasis and Contributes to the Instigation of Stemness in Diseases Entailing Epithelial Hyperplasia or Neoplasia. Stem Cells 2017, 35, 445–457. [Google Scholar] [CrossRef] [Green Version]
  26. Jiang, Z.; Michal, J.J.; Tobey, D.J.; Wang, Z.; Macneil, M.D.; Magnuson, N.S. Comparative understanding of UTS2 and UTS2R genes for their involvement in type 2 diabetes mellitus. Int. J. Biol. Sci. 2008, 4, 96–102. [Google Scholar] [CrossRef] [Green Version]
  27. Meissburger, B.; Perdikari, A.; Moest, H.; Muller, S.; Geiger, M.; Wolfrum, C. Regulation of adipogenesis by paracrine factors from adipose stromal-vascular fraction—A link to fat depot-specific differences. Biochim. Biophys. Acta 2016, 1861, 1121–1131. [Google Scholar] [CrossRef] [Green Version]
  28. Hubal, M.J.; Nadler, E.P.; Ferrante, S.C.; Barberio, M.D.; Suh, J.H.; Wang, J.; Dohm, G.L.; Pories, W.J.; Mietus-Snyder, M.; Freishtat, R.J. Circulating adipocyte-derived exosomal MicroRNAs associated with decreased insulin resistance after gastric bypass. Obesity (Silver Spring) 2017, 25, 102–110. [Google Scholar] [CrossRef] [Green Version]
  29. Liu, K.D.; Acharjee, A.; Hinz, C.; Liggi, S.; Murgia, A.; Denes, J.; Gulston, M.K.; Wang, X.; Chu, Y.; West, J.A.; et al. Consequences of Lipid Remodeling of Adipocyte Membranes Being Functionally Distinct from Lipid Storage in Obesity. J. Proteome Res. 2020, 19, 3919–3935. [Google Scholar] [CrossRef]
  30. Wang, C.; Chao, Y.; Xu, W.; Liu, Z.; Wang, H.; Huang, K. Myeloid FBW7 deficiency disrupts redox homeostasis and aggravates dietary-induced insulin resistance. Redox Biol. 2020, 37, 101688. [Google Scholar] [CrossRef]
  31. Bengoechea-Alonso, M.T.; Ericsson, J. The ubiquitin ligase Fbxw7 controls adipocyte differentiation by targeting C/EBPalpha for degradation. Proc. Natl. Acad. Sci. USA 2010, 107, 11817–11822. [Google Scholar] [CrossRef]
  32. Zhao, J.; Xiong, X.; Li, Y.; Liu, X.; Wang, T.; Zhang, H.; Jiao, Y.; Jiang, J.; Zhang, H.; Tang, Q.; et al. Hepatic F-Box Protein FBXW7 Maintains Glucose Homeostasis Through Degradation of Fetuin-A. Diabetes 2018, 67, 818–830. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Seeley, E.S.; Nachury, M.V. The perennial organelle: Assembly and disassembly of the primary cilium. J. Cell Sci. 2010, 123, 511–518. [Google Scholar] [CrossRef] [Green Version]
  34. Nishimura, Y.; Yamakawa, D.; Shiromizu, T.; Inagaki, M. Aurora A and AKT Kinase Signaling Associated with Primary Cilia. Cells 2021, 10, 3602. [Google Scholar] [CrossRef]
  35. Plotnikova, O.V.; Seo, S.; Cottle, D.L.; Conduit, S.; Hakim, S.; Dyson, J.M.; Mitchell, C.A.; Smyth, I.M. INPP5E interacts with AURKA, linking phosphoinositide Signaling to primary cilium stability. J. Cell Sci. 2015, 128, 364–372. [Google Scholar] [CrossRef] [Green Version]
  36. Mariman, E.C.; Vink, R.G.; Roumans, N.J.; Bouwman, F.G.; Stumpel, C.T.; Aller, E.E.; van Baak, M.A.; Wang, P. The cilium: A cellular antenna with an influence on obesity risk. Br. J. Nutr. 2016, 116, 576–592. [Google Scholar] [CrossRef] [Green Version]
  37. Benot-Dominguez, R.; Cimini, A.; Barone, D.; Giordano, A.; Pentimalli, F. The Emerging Role of Cyclin-Dependent Kinase Inhibitors in Treating Diet-Induced Obesity: New Opportunities for Breast and Ovarian Cancers? Cancers 2022, 14, 2709. [Google Scholar] [CrossRef]
  38. Zineldeen, D.H.T.; Mohamed, N.; El-Sayed, W.R.; El-Khadrawy, O.H. Polo Like Kinase 1 (Plk1) expression in visceral adipose tissue of morbidly obese non diabetic women. Am. J. Res. Commun. 2013, 1, 177–208. [Google Scholar]
  39. Liu, Y.; Wang, W.; Shui, G.; Huang, X. CDP-diacylglycerol synthetase coordinates cell growth and fat storage through phosphatidylinositol metabolism and the insulin pathway. PLoS Genet. 2014, 10, e1004172. [Google Scholar] [CrossRef]
  40. Qi, Y.; Kapterian, T.S.; Du, X.; Ma, Q.; Fei, W.; Zhang, Y.; Huang, X.; Dawes, I.W.; Yang, H. CDP-diacylglycerol synthases regulate the growth of lipid droplets and adipocyte development. J. Lipid Res. 2016, 57, 767–780. [Google Scholar] [CrossRef] [Green Version]
  41. Al Khazal, F.; Kang, S.; Nelson Holte, M.; Choi, D.S.; Singh, R.; Ortega-Saenz, P.; Lopez-Barneo, J.; Maher, L.J., 3rd. Unexpected obesity, rather than tumorigenesis, in a conditional mouse model of mitochondrial complex II deficiency. FASEB J. 2021, 35, e21227. [Google Scholar] [CrossRef]
  42. Mu, C.; Ma, B.; Zhang, C.; Geng, G.; Zhang, X.; Chen, L.; Wang, M.; Li, J.; Zhao, T.; Cheng, H.; et al. Loss of SDHB reprograms energy metabolisms and inhibits high fat diet induced metabolic syndromes. bioRxiv 2018. bioRxiv:10.1101/259226. [Google Scholar] [CrossRef]
  43. Kaartinen, M.T.; Hang, A.; Barry, A.; Arora, M.; Heinonen, S.; Lundbom, J.; Hakkarainen, A.; Lundholm, N.; Rissanen, A.; Kaprio, J.; et al. Matrisome alterations in obesity—Adipose tissue transcriptome study on monozygotic weight-discordant twins. Matrix Biol. 2022, 108, 1–19. [Google Scholar] [CrossRef] [PubMed]
  44. Slattery, M.L.; Lundgreen, A.; John, E.M.; Torres-Mejia, G.; Hines, L.; Giuliano, A.R.; Baumgartner, K.B.; Stern, M.C.; Wolff, R.K. MAPK genes interact with diet and lifestyle factors to alter risk of breast cancer: The Breast Cancer Health Disparities Study. Nutr. Cancer 2015, 67, 292–304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Hirosumi, J.; Tuncman, G.; Chang, L.; Gorgun, C.Z.; Uysal, K.T.; Maeda, K.; Karin, M.; Hotamisligil, G.S. A central role for JNK in obesity and insulin resistance. Nature 2002, 420, 333–336. [Google Scholar] [CrossRef]
  46. Guenard, F.; Lamontagne, M.; Bosse, Y.; Deshaies, Y.; Cianflone, K.; Kral, J.G.; Marceau, P.; Vohl, M.C. Influences of gestational obesity on associations between genotypes and gene expression levels in offspring following maternal Gastrointest.inal bypass surgery for obesity. PLoS ONE 2015, 10, e0117011. [Google Scholar] [CrossRef]
  47. Kogelman, L.J.; Cirera, S.; Zhernakova, D.V.; Fredholm, M.; Franke, L.; Kadarmideen, H.N. Identification of co-expression gene networks, regulatory genes and pathways for obesity based on adipose tissue RNA Sequencing in a porcine model. BMC Med. Genomics. 2014, 7, 57. [Google Scholar] [CrossRef] [Green Version]
  48. Youn, D.Y.; Xiaoli, A.M.; Pessin, J.E.; Yang, F. Regulation of metabolism by the Mediator complex. Bioph. Rep. 2016, 2, 69–77. [Google Scholar] [CrossRef] [Green Version]
  49. Hittelman, A.B.; Burakov, D.; Iniguez-Lluhi, J.A.; Freedman, L.P.; Garabedian, M.J. Differential regulation of glucocorticoid receptor transcriptional activation via AF-1-associated proteins. EMBO J. 1999, 18, 5380–5388. [Google Scholar] [CrossRef] [Green Version]
  50. Grontved, L.; Madsen, M.S.; Boergesen, M.; Roeder, R.G.; Mandrup, S. MED14 tethers mediator to the N-terminal domain of peroxisome proliferator-activated receptor gamma and is required for full transcriptional activity and adipogenesis. Mol. Cell Biol. 2010, 30, 2155–2169. [Google Scholar] [CrossRef] [Green Version]
  51. Derecka, M.; Gornicka, A.; Koralov, S.B.; Szczepanek, K.; Morgan, M.; Raje, V.; Sisler, J.; Zhang, Q.; Otero, D.; Cichy, J.; et al. Tyk2 and Stat3 regulate brown adipose tissue differentiation and obesity. Cell Metab. 2012, 16, 814–824. [Google Scholar] [CrossRef] [Green Version]
  52. Watanabe, S.; Nishimura, R.; Shirasaki, T.; Katsukura, N.; Hibiya, S.; Kirimura, S.; Negi, M.; Okamoto, R.; Matsumoto, Y.; Nakamura, T.; et al. Schlafen 11 Is a Novel Target for Mucosal Regeneration in Ulcerative Colitis. J. Crohns. Colitis. 2021, 15, 1558–1572. [Google Scholar] [CrossRef]
  53. Gershon, M.D. Review article: Serotonin receptors and transporters—Roles in normal and abnormal Gastrointest.inal motility. Aliment Pharmacol. Ther. 2004, 20 (Suppl. 7), 3–14. [Google Scholar] [CrossRef]
  54. Litwack, G. Obesity, 1st ed.; Elsevier Academic Press: Amsterdam, The Netherlands; Boston, MA, USA, 2013; p. xviii. [Google Scholar]
  55. Ye, L.; Liddle, R.A. Gastrointest.inal hormones and the gut connectome. Curr. Opin. Endocrinol. Diabetes Obes. 2017, 24, 9–14. [Google Scholar] [CrossRef]
  56. Basson, M.D.; Modlin, I.M.; Flynn, S.D.; Jena, B.P.; Madri, J.A. Independent modulation of enterocyte migration and proliferation by growth factors, matrix proteins, and pharmacologic agents in an in vitro model of mucosal healing. Surgery 1992, 112, 299–307, discussion 307–298. [Google Scholar]
  57. Basson, M.D.; Beidler, D.R.; Turowski, G.; Zarif, A.; Modlin, I.M.; Jena, B.P.; Madri, J.A. Effect of tyrosine kinase inhibition on basal and epidermal growth factor-stimulated human Caco-2 enterocyte sheet migration and proliferation. J. Cell Physiol. 1994, 160, 491–501. [Google Scholar] [CrossRef]
  58. Basson, M.D.; Modlin, I.M.; Madri, J.A. Human enterocyte (Caco-2) migration is modulated in vitro by extracellular matrix composition and epidermal growth factor. J. Clin. Investig. 1992, 90, 15–23. [Google Scholar] [CrossRef]
  59. Owen, C.R.; Yuan, L.; Basson, M.D. Smad3 knockout mice exhibit impaired intestinal mucosal healing. Lab. Investig. 2008, 88, 1101–1109. [Google Scholar] [CrossRef] [Green Version]
  60. Basson, M.D.; Liu, Y.W.; Hanly, A.M.; Emenaker, N.J.; Shenoy, S.G.; Gould Rothberg, B.E. Identification and comparative analysis of human colonocyte short-chain fatty acid response genes. J. Gastrointest. Surg. 2000, 4, 501–512. [Google Scholar] [CrossRef]
  61. Murnin, M.; Kumar, A.; Li, G.D.; Brown, M.; Sumpio, B.E.; Basson, M.D. Effects of glutamine isomers on human (Caco-2) intestinal epithelial proliferation, strain-responsiveness, and differentiation. J. Gastrointest. Surg. 2000, 4, 435–442. [Google Scholar] [CrossRef]
  62. Perdikis, D.A.; Basson, M.D. Basal nutrition promotes human intestinal epithelial (Caco-2) proliferation, brush border enzyme activity, and motility. Crit. Care Med. 1997, 25, 159–165. [Google Scholar] [CrossRef]
  63. Basson, M.D. Paradigms for mechanical Signal. transduction in the intestinal epithelium. Category: Molecular, cell, and developmental biology. Digestion 2003, 68, 217–225. [Google Scholar] [CrossRef] [PubMed]
  64. Chaturvedi, L.S.; Wang, Q.; More, S.K.; Vomhof-DeKrey, E.E.; Basson, M.D. Schlafen 12 mediates the effects of butyrate and repetitive mechanical deformation on intestinal epithelial differentiation in human Caco-2 intestinal epithelial cells. Hum. Cell 2019, 32, 240–250. [Google Scholar] [CrossRef] [PubMed]
  65. Basson, M.D. Effects of repetitive deformation on intestinal epithelial cells. Inflammopharmacology 2007, 15, 109–114. [Google Scholar] [CrossRef] [PubMed]
  66. Gayer, C.P.; Basson, M.D. The effects of mechanical forces on intestinal physiology and pathology. Cell Signal. 2009, 21, 1237–1244. [Google Scholar] [CrossRef]
Cells 11 03283 g001 550
Figure 1. (A) PCA plot for pre- and post-surgery conditions. (B) IPA pathway enrichment showing top eight pathways for 1116 significant genes. The yellow line represents the threshold of p-value ≤ 0.05. (C) Their network of interaction of topmost pathways is shown. (C) Overlapping pathways: The figure shows overlapping pathways associated with differential genes. The nodes represent pathways and edges are labeled with the number of common genes connecting each node. The four most statistically significant pathways are marked with a * followed by the significant p-value.
Figure 1. (A) PCA plot for pre- and post-surgery conditions. (B) IPA pathway enrichment showing top eight pathways for 1116 significant genes. The yellow line represents the threshold of p-value ≤ 0.05. (C) Their network of interaction of topmost pathways is shown. (C) Overlapping pathways: The figure shows overlapping pathways associated with differential genes. The nodes represent pathways and edges are labeled with the number of common genes connecting each node. The four most statistically significant pathways are marked with a * followed by the significant p-value.
Cells 11 03283 g001
Cells 11 03283 g002 550
Figure 2. Comparison of weight loss in pre- and post-surgery conditions. For a better comparison, BMI was used to monitor changes in weight loss.
Figure 2. Comparison of weight loss in pre- and post-surgery conditions. For a better comparison, BMI was used to monitor changes in weight loss.
Cells 11 03283 g002
Cells 11 03283 g003 550
Figure 3. Comparison of genes changing in the post-RYGB samples with respect to the change in BMI (Table S2) vs. genes identified between the pre- vs. post-RYGB genes (Table S1).
Figure 3. Comparison of genes changing in the post-RYGB samples with respect to the change in BMI (Table S2) vs. genes identified between the pre- vs. post-RYGB genes (Table S1).
Cells 11 03283 g003
Cells 11 03283 g004 550
Figure 4. Linear relationship between gene expression fold change and percentage change in BMI. Each dot represents one patient. The blue line shows the linear relationship, and the red line shows the deviation from the blue line for each patient.
Figure 4. Linear relationship between gene expression fold change and percentage change in BMI. Each dot represents one patient. The blue line shows the linear relationship, and the red line shows the deviation from the blue line for each patient.
Cells 11 03283 g004
Cells 11 03283 g005 550
Figure 5. SLFN12 protein expression increases after gastric bypass. (A) Immunohistochemical staining in the intestinal mucosa for SLFN12 protein in pre- and post-RYGB samples and scoring by blinded observers. (Two observers each achieved similar results, * p ≤ 0.05.) (B) Comparison of patient BMI and SLFN12 protein intensity pre- and post-gastric bypass. (C) Correlation graph of the BMI difference vs. the SLFN12 protein intensity difference between pre- and post- gastric bypass samples. (n = 11).
Figure 5. SLFN12 protein expression increases after gastric bypass. (A) Immunohistochemical staining in the intestinal mucosa for SLFN12 protein in pre- and post-RYGB samples and scoring by blinded observers. (Two observers each achieved similar results, * p ≤ 0.05.) (B) Comparison of patient BMI and SLFN12 protein intensity pre- and post-gastric bypass. (C) Correlation graph of the BMI difference vs. the SLFN12 protein intensity difference between pre- and post- gastric bypass samples. (n = 11).
Cells 11 03283 g005
Cells 11 03283 g006 550
Figure 6. Overexpression of SLFN12 in human intestinal epithelial cell lines influences the expression of differentiation, restitution, and obesity-related genes. RNA expression of (A) SLFN12 was measured from either HIEC-6 or FHs 74 Int cells that were treated with either CMV adenovirus control or AdvSLFN12. The following genes are involved in differentiation and restitution of intestinal mucosa: (B) HES2, (C) CARD9, (D) SLC19A2, (E) TRPC5, while the following genes are related to obesity: (F) FBXW7, (G) MBOAT2, (H) STXBP4, (I) SPARCL1, and (J) UTS2. (n = 15–24, **** p < 0.0001, *** p < 0.001, ** p < 0.01).
Figure 6. Overexpression of SLFN12 in human intestinal epithelial cell lines influences the expression of differentiation, restitution, and obesity-related genes. RNA expression of (A) SLFN12 was measured from either HIEC-6 or FHs 74 Int cells that were treated with either CMV adenovirus control or AdvSLFN12. The following genes are involved in differentiation and restitution of intestinal mucosa: (B) HES2, (C) CARD9, (D) SLC19A2, (E) TRPC5, while the following genes are related to obesity: (F) FBXW7, (G) MBOAT2, (H) STXBP4, (I) SPARCL1, and (J) UTS2. (n = 15–24, **** p < 0.0001, *** p < 0.001, ** p < 0.01).
Cells 11 03283 g006
Table
Table 1. Comparison of expression of differentiation, restitution, and obesity-related genes between RNA Seq and HIEC-6 and FHs74 Int in the in vitro data.
Table 1. Comparison of expression of differentiation, restitution, and obesity-related genes between RNA Seq and HIEC-6 and FHs74 Int in the in vitro data.
RNA SeqHIEC-6 + AdvSLFN12FHs 74 Int + AdvSLFN12
Differentiation & Restitution genesHES2↑↑=
CARD9↑=↑
TRPC5↓↓↑
SLC19A2↑↑↑
Obesity–related genesFBXW7↑↑↑
MBOAT2↑↓↓
STXBP4↑↑↑
SPARCL1↑↑=
UTS↑↑↑
Table
Table 2. Major canonical pathways affected by RYGB.
Table 2. Major canonical pathways affected by RYGB.
Top Canonical Pathwaysp-ValueOverlap
Cell Cycle: G2/M DNA Damage Checkpoint Regulation1.21 × 10−36/49 (AURKA, CKS1B, CKS2, PKMYT1, PLK1, TOP2A)
CDP-diacylglycerol Biosynthesis I2.29 × 10−24/24 (ACO2, SDHB, SDHC, SUCLA2)
Phosphatidylglycerol Biosynthesis II (Nonplastidic)2.50 × 10−23/12 (POLB, POLG, XRCC1)
TCA Cycle II (Eukaryotic)2.5 × 10−26/66 (ADPGK, FAM20B, G6PC, GRK4, MAP3K9, PLK1)
Pyrimidine Ribonucleotides Interconversion3.42 × 10−24/36 (IFI35, IFNAR2, MED14, TYK2)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vomhof-DeKrey, E.E.; Singhal, S.; Singhal, S.K.; Stover, A.D.; Rajpathy, O.; Preszler, E.; Garcia, L.; Basson, M.D. RNA Sequencing of Intestinal Enterocytes Pre- and Post-Roux-en-Y Gastric Bypass Reveals Alteration in Gene Expression Related to Enterocyte Differentiation, Restitution, and Obesity with Regulation by Schlafen 12. Cells 2022, 11, 3283. https://doi.org/10.3390/cells11203283

AMA Style

Vomhof-DeKrey EE, Singhal S, Singhal SK, Stover AD, Rajpathy O, Preszler E, Garcia L, Basson MD. RNA Sequencing of Intestinal Enterocytes Pre- and Post-Roux-en-Y Gastric Bypass Reveals Alteration in Gene Expression Related to Enterocyte Differentiation, Restitution, and Obesity with Regulation by Schlafen 12. Cells. 2022; 11(20):3283. https://doi.org/10.3390/cells11203283

Chicago/Turabian Style

Vomhof-DeKrey, Emilie E., Sonalika Singhal, Sandeep K. Singhal, Allie D. Stover, Odele Rajpathy, Elizabeth Preszler, Luis Garcia, and Marc D. Basson. 2022. "RNA Sequencing of Intestinal Enterocytes Pre- and Post-Roux-en-Y Gastric Bypass Reveals Alteration in Gene Expression Related to Enterocyte Differentiation, Restitution, and Obesity with Regulation by Schlafen 12" Cells 11, no. 20: 3283. https://doi.org/10.3390/cells11203283

APA Style

Vomhof-DeKrey, E. E., Singhal, S., Singhal, S. K., Stover, A. D., Rajpathy, O., Preszler, E., Garcia, L., & Basson, M. D. (2022). RNA Sequencing of Intestinal Enterocytes Pre- and Post-Roux-en-Y Gastric Bypass Reveals Alteration in Gene Expression Related to Enterocyte Differentiation, Restitution, and Obesity with Regulation by Schlafen 12. Cells, 11(20), 3283. https://doi.org/10.3390/cells11203283

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