Abstract

The ineffective short-term control of feeding behavior compromises energy homeostasis and can lead to obesity. The gastrointestinal tract secretes several regulatory peptides. However, little is known about the stomach peptide contribution to the acute regulation of intake. In an attempt to identify new gastric signals, the serial analysis of gene expression (SAGE) method was used for the transcription profiling of stomach mucosa in 7 groups of mice: fasting and sacrificed 30 minutes, 1 hour, 3 hours after a low-fat (LF) or high-fat (HF) ad libitum meal. In total, 35 genes were differentially modulated by LF and HF meals compared to fasting, including 15 mRNAs coding for digestive enzymes/secretory proteins, and 10 novel transcripts. Although the basic expression profile did not undergo substantial variations, both LF and HF meals influenced the transcription. This study represents the first global analysis of stomach transcriptome as induced by different nutritional stimuli. Further studies including the characterization of novel genes may help to identify new targets for the therapy and prevention of obesity.

1. Introduction

Obesity epidemic continues its worrying global progression although significant advances have been achieved in the knowledge of its causes and consequences. This condition, in concert with glucose intolerance/type 2 diabetes, dyslipidemia, and metabolic syndrome, widely contributes to what has been recently defined as “prosperity’s plague” [1]. A complex system of social, psychological, physiological, and biological factors has to be considered in order to successfully control this “plague” and prevent from its further spread [2, 3]. From a physiological point of view, it is fundamental to understand the specific relative importance of long-term and short-term mechanisms involved in the regulation of energy balance [4]. Food intake and daily overconsumption may have a predominant impact on body weight regulation, and this led large interest to focus on appetite/satiety balance as one of the key potential therapeutic targets [5]. Multiple sites in the gastrointestinal (GI) tract, including the stomach, proximal and distal small intestine, colon, and pancreas, are involved in the short-term regulation of energy homeostasis, which basically controls what, when, and how much we eat within a single day or a single meal [6]. In addition to mechanoreceptors and chemoreceptors, which are activated during a meal and signal to the brainstem through the vagal nerve [4], several gut-derived peptides and lipid mediators have a role in the regulation of food intake and energy homeostasis [7]. The only recognized “hunger” gut peptide so far discovered is ghrelin. This peptide is mainly produced by the stomach and directly influences the number of meals consumed per day though it has probably no direct effect on meal size [8]. However, the stomach plays an essential mechanical role in the regulation of satiety perception and meal termination. It is, in fact, the first organ to receive the bolus of food. Here, the food ingested is rapidly homogenized, partially digested, and finally delivered to the small intestine. Importantly, gastric distension and gastric emptying via the pylorus are finely regulated in order to match food delivery to the actual gut digestive and absorptive capacity [9]. As a consequence, gastric competence can be considered as the first limiting step of GI ingestive and digestive capacity and thus represents a relevant target for obesity prevention and treatment. Motility and physical mechanisms involved in gastric-mediated satiety have been extensively studied and successfully targeted in recent years [10, 11]. However, it is still unknown whether the gastric mucosa also releases any biochemical signal that may influence satiety in the very short term of meal consumption. In an attempt to explore this possibility, the serial analysis of gene expression (SAGE) method was used to identify the early transcriptional changes induced by a low-fat (LF) or high-fat (HF) meal in the gastric mucosa. The discovery potential of SAGE was determinant for choosing this technology instead of other comparable transcriptomic methods. SAGE, indeed, is a powerful and reliable sequencing-based technique [12] which allows to detect the regulation of novel transcripts as well as characterized genes, as we have already shown in a number of previous studies [13, 14].

2. Materials and Methods

2.1. Animals, Diet, and Sample Preparation

A total of 140 male C57BL6 mice (12-week-old) were purchased from Charles River Canada Inc. and were acclimated for two weeks. Low-fat diet (LF, Research Diet no. 12450B: 10% calories from fat, 70% from carbohydrate, and 20% from protein; 3.85 kcal/g) and tap water were served ad libitum. In the last day of the acclimatization, the body weight paired mice were randomly distributed into seven groups and fasted for 12 hours during darkness of the light cycle. On the experimental day, one group of fasted mice was sacrificed before the meal. The other mice were fed ad libitum with high-fat (HF, Research Diet no. 12492: 60% calories from fat, 20% from carbohydrate, and 20% from protein; 5.24 kcal/g) or LF meal and sacrificed 30 min, 1 h, and 3 h after the beginning of the meal. The amount of macronutrients and energy ingested was recorded. In total, seven groups of mice under isoflurane anesthesia (fasting, HF 30 min/1 h/3 h, and LF 30 min/1 h/3 h) were alternatively exsanguinated by cardiac puncture after cervical dislocation. Stomach was opened vertically and flushed clean with saline, and the mucosa was removed by scrapping with a glass microscope slide. The samples were rapidly frozen in liquid nitrogen and stored at −80°C until RNA extraction. All animal experimentation was conducted in accordance with the requirements of the Canadian Council on Animal Care and approved by the Animal Protection Committee of Laval University.

2.2. Transcriptomic Analysis

The seven serial analyses of gene expression (SAGE) libraries were constructed as previously described [15]. Total RNA was isolated from pooled stomach mucosa for each group ( 𝑛 = 2 0 ) by Trizol (Invitrogen Canada Inc., Burlington, ON). The quality of total RNA was monitored by microcapillary electrophoresis (Bioanalizer 2100, Agilent Technologies, Mississauga, ON). Polyadenylated RNA was extracted (Oligotex mRNA Mini Kit, Qiagen Inc., Mississauga, ON), annealed with the biotin-5′-T18-3′ primer, and converted to cDNA (cDNA synthesis kit, Invitrogen Canada Inc.). The resulting cDNAs were digested with NlaIII (New England BioLabs Ltd., Pickering, ON), and the 3′ restriction fragments were isolated with streptavidin-coated magnetic beads (Dynal Biotech LLC, Brown Deer, WI) and separated into two populations. Each population was ligated to one of two annealed linkers and extensively washed to remove unligated linkers [15]. The tag beside the most 3′NlaIII restriction site (CATG) of each transcript was released by digestion with BsmFI (New England BioLabs Ltd.). The blunting kit from Takara Bio Inc. (Otsu, Japan) was used for the blunting and ligation of the two tag populations. The resulting ligation products containing the ditags were amplified by PCR and digested with NlaIII. The band containing the ditags was extracted from the 12% polyacrylamide gel. The purified ditags were self-ligated to form concatemers using T4 ligase (Invitrogen Canada Inc.). The isolated 500 bp to 1800 bp concatemers were isolated by agarose gel, and the resulting DNA fragments were ligated into the SphI site of pUC19 and cloned into OmniMAX 2T1 competent cells (Invitrogen Canada Inc.). White colonies were picked up, and the concatemer inserts were finally sequenced by the Applied Biosystems 3730 (Foster City, CA).

2.3. Bioinformatic Analysis

Sequence files were analyzed using the SAGEana program, a modification of SAGEparser [12]. Identification of the transcripts was obtained by matching the 15 bp (sequence at the last CATG + 11 bp tags) with SAGEmap, UniGene, and GenBank databases. Classification of the transcripts was based upon the updated information of the genome directory [16] found at the TIGR website (http://www.tigr.org/), the SOURCE (http://genome-www5.stanford.edu/cgi-bin/source/sourceSearch), and the OMIM (http://www.ncbi.nlm.nih.gov/) as well as previously published literatures. We have previously shown that the SAGE method is very reproducible with 𝑟 2 = 0 . 9 6 between two SAGE libraries constructed from the same total RNA pool [12].

2.4. Validation by Quantitative Real-Time PCR (Q_RT-PCR)

First strand cDNA was synthesized using 5 μg of pooled RNA of each experimental group in a reaction containing 200 U of Superscript III Rnase H-RT (Invitrogen Canada Inc.), 50 ng of random hexamers, 300 ng of oligo-dT18, 50 mM Tris-HCl pH 8.3, 75 mM KCl, 3 mM MgCl2, 5 mM dithiothreitol, 0.5 mM deoxynucleotides triphosphate, and 40 U human RNase inhibitor (Roche) in a final volume of 50 μL. The resulting products were then treated with 1 μg of Rnase A for 30 minutes at 37°C and purified thereafter with Qiaquick PCR purification kits (Qiagen). The cDNA corresponding to 20 ng of total RNA was used to perform fluorescent-based real-time PCR quantification using the LightCycler real-time PCR apparatus (Roche Inc., Nutley, NJ) and the FastStart DNA Master SYBR green kit (Roche Diagnostics). Quantifications were conducted in triplicates. Reading of the fluorescence signal was taken at the end of the heating to avoid nonspecific signal. A melting curve was performed to assess nonspecific signals. Annealing temperature was selected based on contamination levels and melting curve results. Prior to mRNA quantification, RNA samples were also verified for genomic DNA contamination. Oligoprimer pairs that allow the amplification of approximately 250 bp were designed by GeneTools software (Biotools Inc., Edmonton, AB) and their specificity was verified by BLAST in GenBank database. Gene name, GenBank accession numbers, and regions used for the primer pairs were the following: chymotrypsin-like elastase family member 3B (elastase 3, Cela3b), NM_026419, 189–385; amylase 2a1, pancreatic (Amy2a1), NM_001042712, 848–1039; pancreatic lipase (Pnlip), NM_026925, 840–1067; major urinary protein 1 (Mup1), NM_031188, 127–405; protein disulfide isomerase associated 3 (Pdia3), BC033439, 994–1222; zymogen granule membrane protein 16 (Zg16), NM_026918, 348–522. The mRNA levels were calculated using a standard curve of crossing point (Cp) versus logarithm of the quantity and expressed as the number of copies per microgram of total RNA [17]. The LightCycler 3.5 program provided by the manufacturer (Roche Inc.) was used to calculate the Cp according to the second derivative and double-correction method previously described by Luu-The et al. [17]. The standard curve with efficiency coefficient 𝐸 = 2 was established using known cDNA amounts of 0, 102, 103, 104, 105, and 106 copies of ATP synthase O subunit.

2.5. Statistical Analysis

Food and macronutrient intakes were analyzed by the two-way ANOVA. When the ANOVA revealed a significant interaction between diet and time, the contrast analysis was performed to identify the significant difference between the HF and LF groups from the same time points ( 𝑃 < . 0 5 ). For the SAGE data, the comparative count display (CCD) test was used to identify the transcripts which were significantly differentially expressed ( 𝑃 . 0 5 ) between the groups with more than a two-fold change, as previously described by Lash et al. [18]. The data were normalized to 50000 tags for presentation. Q_RT-PCR data were analyzed by the two-tailed Student’s t-test ( 𝑃 < . 0 1 ) for the time points formerly determined by the SAGE method.

3. Results

3.1. Food Intake

Food and energy intakes are presented in Figure 1 as cumulative and 30 minutes average consumption. As expected, cumulative energy as well as protein and fat intakes were higher, whereas carbohydrate intake was lower in the HF groups compared to LF (Figure 2(a)). Interestingly, a distinct pattern of feeding behavior between LF- and HF-fed mice can be observed by comparing the average 30 minutes consumption of the two groups (Figure 1(c)). Mice assigned to LF meal ate a moderate amount of food during one hour, after which their consumption decreased to a minimum intake. In contrast, the HF group consumed a large amount of meal in the first 30 minutes although the intake dropped in the following 30 minutes, reaching the minimum level at 1 h. However, the HF-group mice ingestion increased again in the last two hours.

3.2. Transcriptomic Analysis

Seven SAGE libraries were generated to identify the transcripts differentially modulated in mouse gastric mucosa by the following experimental conditions: fasting; LF or HF meal at 30 min, 1 h and 3 h since the beginning of consumption (LF 30 m/1 h/3 h and HF 30 m/1 h/3 h, resp.). Among the 56382 SAGE tag species detected, a total of 35 transcripts were significantly regulated by LF and HF feeding compared to fasting, whereas 19 were specifically modulated by HF compared to LF.

3.2.1. Transcripts Coding for Digestive/Secretory Proteins

The most represented group includes transcripts coding for digestive enzymes and secretory pathway components (Table 1). Among the modulated mRNAs were amylase 2a1, pancreatic lipase, carboxyl ester lipase, elastase 1 and 3, and carboxypeptidase A1 and B1. Globally, three of them have lipolytic functions, nine code for proteolytic enzymes and one gene is involved in carbohydrate digestion. In addition, two genes involved in zymogen granule secretion were regulated, namely syncollin and zymogen granule membrane protein 16. Interestingly, most of these genes showed a common pattern of regulation since their expression was downregulated at LF 30 m and 1 h, as well as at HF 1 h and 3 h, compared to fasting. Moreover, the expression of these transcripts tended to increase at LF 3 h though the up-regulation was statistically significant only for six of them, such as amylase 2a1, chymotrypsinogen B1, and pancreatic lipase related protein 1. Remarkably, for 13 out of the 15 tags considered, the specific differential regulation between LF 3 h and HF 3 h achieved statistical significance. There is, therefore, a common downregulation of these transcripts following both LF and HF feeding although a temporal delay between the two groups can be observed. In addition, 3 hours after the beginning of the LF meal, their transcription tended to increase at higher levels than those observed at the fasting state.

3.2.2. Cellular Defence Related and Other Transcripts Differentially Regulated by Feeding

Feeding significantly regulated eight transcripts with protective roles (Table 2). Among these, the mRNAs coding for the heat shock protein 70 (Hsp70) 1A and expressed sequence tag (EST) Hsp70 1B were both upregulated by the LF meal, and the latter also by the HF meal, at 1 h following the beginning of ingestion. Gene expressions of cysteine-rich protein 1, metallothionein 2, and EST WD repeat domain 92 were reduced at LF 3 h compared to fasting, whereas onzin mRNA levels significantly decreased at HF 3 h. Moreover, HF specifically modulated the transcription of protein disulfide isomerase associated 3 and EST elastase 2A/neutrophil elastase, respectively, up- and downregulated at 3 h.

The results also showed the differential regulation of two transcripts coding for proprotein convertase proteins (furin and kallikrein) and an interesting new candidate as potential regulator of energy metabolism, namely Mup 1 (Table 2). In addition, LF and HF feeding significantly regulated 12 novel transcripts with no match in public databases (Table 3). In particular, the tags GGAGAACAGCG and CTGACTCAAAT were specifically modulated by HF feeding compared to LF and could represent potential targets for further characterization studies.

3.3. Q_RT-PCR Confirmation of SAGE Results

To validate the SAGE results, the Q_RT-PCR analysis was also performed for some of the genes differentially regulated by feeding. The chosen genes are representative of the functional groups discussed. As presented in Figure 3, the Q_RT-PCR results globally confirmed the changes in expression level as well as the significant modulation highlighted by the SAGE method.

4. Discussion

4.1. Fasting and Feeding States Modulated Pancreas-Related Genes in the Stomach

The main regulated functional group is represented by digestive enzyme-coding genes, many of which are mostly expressed by the pancreas [19]. In addition, two transcripts involved in the secretory pathway, syncollin, and zymogen granule membrane protein 16 were also modulated. Many of the physiological changes induced by food intake should arise within minutes. Therefore, the digestive enzymes and regulatory factors are normally synthesized and packed in secretory granules at rest, ready to be released when food ingestion and appropriate neurohumoral stimuli occur. Following the meal, another cycle of synthesis and packaging prepares the zymogenic cells to the next secretory events [9]. Concordantly, the present results showed a decreased transcription of digestive enzyme-coding and secretory genes following the start of ingestion, when the cells are more likely to invest energy in secretion. Moreover, in the LF-fed mice, the reinduction of these genes that occurred 3 hours after the beginning of intake seems a reasonable event, particularly if referred to their feeding pattern (Figure 2). Conversely, the digestion of a high-calorie and high-energy density meal, including its rate of emptying from the stomach, normally takes a longer time to be accomplished [20, 21]. The latter point may explain the delayed and prolonged downregulation of digestive transcripts in the HF group, which also suggests that the reinduction of transcription observed at LF 3 h may have started later in the HF-fed mice. Interestingly, though all the mice had ad libitum access to food, feeding behavior and total ingestion were dissimilar between LF and HF groups (Figures 1 and 2), possibly explaining the differences in the transcriptional regulation of digestive enzymes. However, the levels of digestive enzymes and secretory proteins might as well influence the intake. A physiological redundant excess of digestive capacity characterizes the GI system and guarantees the effectiveness of nutrition. Moreover, the excess capacity for nutrient uptake, including the excess of surface, specialized cells, digestive enzymes, and other secretory products, is proportionally related to body weight [9]. Hence, this may also contribute to compromise the efficient control of appetite/satiety balance in overweight and obese subjects. The stomach may represent a primary target for the control of meal size and satiation. Eventually, surgical options which physically affect gastric capacity and emptying mechanisms successfully modify the eating behavior and metabolic profile of obese patients [10] though they still are invasive and lifestyle-affecting methods [22]. Likewise, mechanisms other than mechanical may also affect gastric volume and emptying rates to regulate satiation during meal consumption. Presently, there still is a paucity of literature addressing the specific effects of HF intake on the neurohormonal control of gastric capacity and motility mechanisms. Therefore, the search for these pathways was the principal aim of the current study.

4.2. Plausible Reasons for the Expression of Pancreatic Genes in the Stomach

The observation that gastric mucosa could express pancreatic genes was unexpected. In particular, it would be useful to explain the specific role of these genes in the stomach and, most interestingly, their differential regulation in response to fasting and feeding. Considering the high level of transcription, such as for amylase 2a1 at fasting and LF 3 h, the expression of these genes should entail a needed role in gastric mucosa. In contrast, it is hard to explain how the digestive enzymes, normally active at a neutral pH, could conceivably work in such an acidic milieu. However, it should also be considered that these proteins are often released as precursors and eventually activated by specific signals or the appropriate pH.

This was not the first time that the expression and activity of pancreas-related genes were detected in the stomach and other nonpancreatic components of the GI system. Terada and colleagues [19] had already showed the expression of alpha-amylase, trypsin, chymotrypsin, and pancreatic lipase in normal and pathologic epithelial cells of gastric mucosa by immunohistochemistry and western blotting. They had also proved their enzymatic activity in stomach specimens, although to a far lesser extent than in pancreas. In that report, the authors explained the presence of pancreatic enzymes in nonpancreatic tissues as a result of the common embryonic origin (foregut) shared by the gastroenteric tissues. In addition, the mRNA expression levels of representative genes have been further confirmed by Q_RT-PCR in the present study. It can be hypothesized that the digestive enzymes expressed by the gastric mucosa would combine with the bolus of homogenized and partially digested food that finally reaches the small intestine. In normal conditions, the digestive enzymes would be in excess, but enough to guarantee the effectiveness of digestion in case of insufficient pancreatic secretion. The redundant production of these proteins along the digestive tract was also suggested by our previous study of the duodenum transcriptome, where analogue experimental conditions had been applied. In the intestine, the same pancreas-related transcripts were modulated ( 𝑛 = 9 ) [13] but showed the opposite trend, being upregulated at LF 1 h and HF 3 h. Interestingly, in both the duodenal and gastric mucosa, the regulation of digestive gene expression presented a temporal delay between the LF and HF groups.

4.3. LF and HF Modulated Genes Involved in the Protection and Defense of Gastric Mucosa

Protection is one of the main functions and challenges of gastric mucosa. The mucosal epithelium is a primary barrier, which defends the whole organism from external dangerous agents eventually ingested. Moreover, uncontrolled acid secretion, inflammation, oxidative stress, and the epithelial damages consequently engendered can compromise the physical and functional integrity of the mucosal barrier [23]. Among the transcripts differentially regulated in the gastric mucosa, seven were related to defense and protection. Metallothionein 2 [24, 25], Hsp70 (1A and 1B) [26, 27], and protein disulfide isomerase associated 3 genes [28, 29] code for multiple-task proteins, which can show chaperone activity and/or be involved in cell redox homeostasis control and apoptosis regulation. The latter is a very important role, since gastric epithelium is subject to constant renewal. The epithelial cells, in fact, rapidly turnover (1–3 days in humans), undergoing a cycle of division and differentiation before succumbing to apoptosis [9, 30]. In this study, EST Hsp70 1B was upregulated by both LF and HF at 1 h, and the same trend was observed for hsp70 1A. However, for the latter, the HF 1 h increase did not reach a statistical significance. Interestingly, a polymorphism of Hsp70 1B gene has been recently associated with obesity-related traits [31], thus stimulating further questions about its acute modulation by food intake.

4.4. Mup1, Already Found in Hypothalamus and Duodenum, was Expressed in the Gastric Mucosa and Regulated by Feeding

In the present study, the gene coding for Mup1 was specifically downregulated at HF 30 m compared to LF. In previous transcriptomic studies conducted with the SAGE method, Mup1 gene was also significantly regulated by feeding in the duodenum mucosa and hypothalamus of mice [13, 14]. The corresponding protein is a pheromone transporter normally expressed by the liver. However, recent studies in mice have surprisingly revealed that Mup1 is also involved in glucose and lipid metabolism [32], and that it might play an important role in the regulation of energy expenditure [33]. These findings raise the interest about the specific role of this molecule at the tissue level but also as a potential modulator of energy balance.

4.5. Other Interesting Genes Differentially Regulated by Feeding

Globally, the most regulated class of genes was the one coding for proteolytic enzymes, particularly the serine-protease type. This group of proteases is highly represented in nature and shows numerous and functionally diverse functions, ranging from digestion and coagulation to apoptosis and immunity [34]. In addition to the digestion-related genes described above, feeding regulated two other transcripts coding for highly important serine proteases, namely kallikrein and furin. EST kallikrein 1 was specifically regulated by HF 3 h compared to LF, whereas furin was downregulated at LF 1 h. These two molecules act as proprotein convertases in distinct regulatory pathways. The first cleaves kininogen to produce kinin peptide [35], whereas furin processes the precursors of a large variety of proteins, including growth factors and receptors [36]. Interestingly, the specific pathways involving kallikrein-kinin [37, 38] and other proteins of furin family [39, 40] are presently being studied for their potential contribution to obesity and cardiovascular disorders.

5. Conclusion

The present study was the first to analyze the global transcriptional changes acutely induced in mouse stomach mucosa by feeding and, in particular, by different nutritional stimuli. The principal aim was to identify new signals specifically induced by HF intake in the short term of meal consumption. Given the weakest satiation power of HF compared to LF foods, these signals may represent potential pharmacological targets for the early modulation of appetite/satiety balance. In this study, both LF and HF regulated gene expression in gastric mucosa, and 17 known genes were differentially modulated by HF compared to LF intake. In addition, a number of novel tags were significantly regulated, some of which may be good objects for future characterization studies. However, a lower number of genes was regulated in the stomach compared to duodenum, when the same experimental conditions have been applied [13]. This may suggest that gastric mucosa has a restricted role in the acute regulation of food intake and mainly centered on meal initiation than meal size/termination control. Another plausible hypothesis is that satiation signals eventually raising from the mucosa could be induced earlier than 30 min after the beginning of the meal, at least at the transcriptional level. Although it is still uncertain whether gastric mucosa releases an early molecular signal specifically involved in satiation control, the present study contributed to highlight some potential mediators of this process. In addition, the characterization of novel regulated genes could stimulate future investigations. Since signals secreted by gastric mucosa may be the optimal targets for appetite control and obesity therapeutic strategies, further research efforts are deserved.

Author’s Contributions and Conflict of Interests

The paper has been approved by all listed authors and there is no conflict of interest that would prejudice its impartiality. M. R. De Giorgio analyzed and interpreted the SAGE data, and drafted the paper. M. Yoshioka and J. St-Amand conceived the study, designed it and critically revised the paper. J. St-Amand gave the final approval of the version to be published.

Aknowledgments

The authors would like to thank Mrs. Jean-Philippe Dionne, René Paradis, and Philippe Rigault for the bioinformatic support. This work was supported by the Heart and Stroke Foundation of Canada (HSFC) and the Canadian Institutes of Health Research (CIHR). Dr. J. St-Amand is a senior investigator supported by Fonds de la Recherche en Santé du Québec.