Time-course microarray analysis for identifying candidate genes involved in obesity-associated pathological changes in the mouse colon

Background Obesity is known to increase the risk of colorectal cancer. However, mechanisms underlying the pathogenesis of obesity-induced colorectal cancer are not completely understood. The purposes of this study were to identify differentially expressed genes in the colon of mice with diet-induced obesity and to select candidate genes as early markers of obesity-associated abnormal cell growth in the colon. Methods C57BL/6N mice were fed normal diet (11% fat energy) or high-fat diet (40% fat energy) and were euthanized at different time points. Genome-wide expression profiles of the colon were determined at 2, 4, 8, and 12 weeks. Cluster analysis was performed using expression data of genes showing log2 fold change of ≥1 or ≤−1 (twofold change), based on time-dependent expression patterns, followed by virtual network analysis. Results High-fat diet-fed mice showed significant increase in body weight and total visceral fat weight over 12 weeks. Time-course microarray analysis showed that 50, 47, 36, and 411 genes were differentially expressed at 2, 4, 8, and 12 weeks, respectively. Ten cluster profiles representing distinguishable patterns of genes differentially expressed over time were determined. Cluster 4, which consisted of genes showing the most significant alterations in expression in response to high-fat diet over 12 weeks, included Apoa4 (apolipoprotein A-IV), Ppap2b (phosphatidic acid phosphatase type 2B), Cel (carboxyl ester lipase), and Clps (colipase, pancreatic), which interacted strongly with surrounding genes associated with colorectal cancer or obesity. Conclusions Our data indicate that Apoa4, Ppap2b, Cel, and Clps are candidate early marker genes associated with obesity-related pathological changes in the colon. Genome-wide analyses performed in the present study provide new insights on selecting novel genes that may be associated with the development of diseases of the colon. Electronic supplementary material The online version of this article (doi:10.1186/s12263-016-0547-x) contains supplementary material, which is available to authorized users.


Background
Obesity is a major global health problem that is closely associated with non-communicable diseases with rapidly increasing incidence, including type 2 diabetes, hypertension, cardiovascular diseases, and some cancers [1]. Excess energy intake contributes to abnormal intermediate conditions such as hyperinsulinemia, hyperglycemia, and dyslipidemia, leading to the development of obesityrelated metabolic complications [2].
Epidemiological evidence indicates that excess body fat is associated with an increased risk of colorectal cancer (CRC) [3]. The risk of CRC increases by 7% with an increase in body mass index (BMI) by 2% [4]. Experimental studies also indicate that diet-induced obesity causes pathological changes in the colon. The number of polyps is significantly higher, and the areas of hyperplasia in the colonic mucosa and inflammatory foci throughout the gastrointestinal tract are broader in high-fat diet (HFD)fed mice than in control mice [5]. Mice fed HFD for two third of their life span and not treated with carcinogenic chemicals show substantially higher incidence and multiplicity of colon tumor than mice fed a control diet [6]. Increased circulating concentrations of insulin and leptin are linked to abnormal hyperproliferation of colorectal tissue and inflammation possibly by controlling transcription factors involved in the expression of cell growth-regulating molecules [7][8][9][10][11][12]. Whole-colon proteomic analyses of wild-type and leptin-deficient ob/ob mice suggest that 40 differently expressed proteins are associated with obesity-related pathological changes in the colon [13]. However, to our knowledge, no study has identified candidate molecules involved in obesityassociated pathological changes in the colon of HFD-fed mice. Moreover, limited information is available on mechanisms underlying the pathophysiological changes in the colon tissue of obese animals.
Interactions between nutritional factors and cellular events in the biological system are extremely complicated. Traditional nutrition research design involving one or two molecular targets often cannot explain phenotypic changes induced by missing responses of other important targets to nutritional stimuli. Recent developments in genome-wide analyses have been used to identify biomarkers that respond to nutritional intervention such as HFD. Several studies indicate that dietinduced obesity changes gene expression patterns in various tissues. Expression of key adipose transcription factors that regulate adipogenesis and insulin sensitivity, including leptin, resistin, uncoupling protein-2, tumor necrosis factor-alpha (TNF-α), CCAAT/enhancer-binding protein α, peroxisome proliferator-activated receptor, sterol regulatory element-binding transcription factor 1, and hydroxysteroid 11-beta dehydrogenase 1, is changed in the gonadal fat tissue of HFD-fed animals [14,15].
HFD also alters the expression of interferon-gamma, interleukin-4, interleukin-10, interleukin-12, and TNF-α in the liver tissue [16]. Despite a strong association between obesity and pathophysiological changes in the small intestine and colon that lead to the development of ulcerative colitis, irritable bowel syndrome (IBD), and CRC, only few studies have examined the association between dietinduced obesity and gene expression pattern of the intestinal tissue [17,18]. A recent study reported substantial changes in lipid metabolism-related gene expression in the small intestine of animals fed long-chain fatty acids of marine origin [19]. Our present study is the first to report global transcriptional changes at different time points during the development of diet-induced obesity in the colon of HFD-fed animals. In addition, we performed bioinformatics analyses to identify candidate early marker genes that might be involved in obesity-related pathological events such as CRC and IBD.

Animals
This study was performed in accordance with the Guide for the Care and Use of Laboratory Animals developed by the Institute of Laboratory Animal Resources of the National Research Council [20] and was approved by the Institutional Animal Care and Use Committee of Yonsei University in Seoul, Republic of Korea (Permit No.: 2010-0039). Eighty 5-week-old male C57BL/6N mice (Orient, Gyeonggi-do, Korea) were housed in a temperature (21 ± 2°C)-and humidity (50 ± 5%)-controlled room with a 12-h light/12-h dark cycle. The mice were fed a commercial diet (Purina, St. Louis, MO, USA) for 1 week and were randomly assigned to receive normal diet (ND, n = 40) and HFD (n = 40). HFD contained 200 g fat/kg (170 g lard plus 30 g corn oil) and 1% cholesterol by weight. It was formulated to provide 40% of the total energy from fat by replacing carbohydrates with lard and corn oil; however, it contained the same amount of vitamins and minerals per kilocalorie as those in the ND. Compositions of the experimental diets are presented in Additional file 1: Table S1. The mice were fed the experimental diets and water ad libitum. Food intake of the mice was recorded daily, and their body weights were measured every 3 days. Ten mice per group were sacrificed at 2, 4, 8, and 12 weeks of feeding the experimental diets by fasting them overnight and by anesthetizing them with diethyl ether. Their colons were laid flat on a glass plate, and the colonic mucosa was scraped using a glass slide. The colon samples were stored at −80°C until their use.

Time-course microarray analysis
Total RNA was isolated from the colon tissue of each mouse, using TRIzol (Invitrogen Life Technologies, Carlsbad, CA, USA), and was purified using RNeasy column (Qiagen, Valencia, CA, USA), according to the manufacturer's protocols. RNA purity and integrity were evaluated by denaturing gel electrophoresis, OD 260 / OD 280 ratio, and analyzed on the Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). The RNA Integrity Number (RIN) score was generated on the Agilent software, and the average RIN score of all samples used for microarray analysis was 8.5 ± 0.9 (mean ± SD). To reduce individual variability in gene expression, identical amounts of total colonic RNA were pooled from 10 mice in each experimental group and a pooled RNA sample representing the ND and HFD group at 2, 4, 8, and 12 weeks was subjected to microarray experiment as described previously [21]. Total RNA was amplified and purified using the Illumina® TotalPrepTM-96 RNA Amplification Kit (Ambion, Austin, TX, USA) to produce biotinylated complementary RNA (cRNA), according to the manufacturer's instructions. The cRNA obtained was quantified using an ND-1000 Spectrophotometer (NanoDrop, Wilmington, DE, USA). The biotinylated cRNA was hybridized onto the Illumina Mouse WG-6 v2.0 Expression BeadChip (Illumina, Inc., San Diego, CA, USA) containing 45,281 probes representing 30,584 genes. After washing and staining, the BeadChip was scanned with the Illumina Bead Array Reader Confocal Scanner according to the manufacturer's instructions. Raw data were exported and analyzed using BeadStudio v3.1.3 (Gene Expression Module v3.3.8; Illumina). All the data analyses and visualization of differentially expressed genes were conducted using ArrayAssist® (Stratagene, La Jolla, CA, USA). Values are expressed as log 2 fold change and were obtained by comparing the gene expression profiles of HFD-fed mice with those of ND-fed mice. Genes showing log 2 fold change of ≥1 or ≤−1 (fold change of ≥2 or ≤−2) were selected, and functional analysis was performed using PANTHER database system (www.patherdb.org). Clustering analysis was performed using genes showing similar expression trends over time. MultiExperiment Viewer program was used to evaluate Kmeans algorism [22]. A gene cluster showing the highest fluctuation over time was selected, and biological processes associated with these HFD-responsive genes over time were analyzed using Database for Annotation, Visualization and Integrated Discovery (DAVID, https:// david.ncifcrf.gov/) [23]. Virtual interaction networktargeted genes in the selected cluster were determined using Michigan Molecular Interactions software [24,25]. In this network, genes that interacted with genes in the protein interaction data consolidated from seven public databases ( [26][27][28] (Fig. 1).

Real-time quantitative polymerase chain reaction analysis
Real-time quantitative polymerase chain reaction (qPCR) was conducted to validate microarray data of several differentially expressed genes that were selected based on the clustering and network analyses and that were associated with the biological function of interest, including CRC and obesity. Template RNA isolated from the colon tissue was reverse transcribed using Superscript™ II RT-PCR System (Invitrogen, Karlsruhe, Germany), according to the manufacturer's instructions, for performing dT 20-primed complementary DNA (cDNA) synthesis. Next, real-time qPCR was performed using an ABI PRISM 7900HT Sequence Detection System (Applied Biosystems, Foster City, CA, USA) in 384-well microtiter plates containing a final reaction volume of 10 μl. Four primer/TaqMan probe combinations were designed based on the following sequences obtained from an NCBI public database: Apoa4, Mm00431814_m1; Cel, Mm00486975_m1; Clps, Mm00517960_m1; and Ppap2b, Mm00504516_m1. Amplifications were performed using the following protocol: initial template denaturation at 95°C for 10 min, followed by 40 cycles at 95°C for 15 s and 60°C for 1 min. All the samples were amplified in triplicate, and data were analyzed using Sequence Detector software (Applied Biosystems).

Statistical analysis
Differences among mice in the two dietary groups were analyzed by Student's t test, with SAS 9.4 (SAS Institute, Inc., Cary, NC, USA). Results were considered statistically significant if two-tailed P values were <0.05.

Results
Time-course of changes in body weight, visceral fat pad weight, and food efficiency ratio during the development of diet-induced obesity C57BL/6N mice fed HFD for 2 weeks gained significantly more weight than mice fed ND (P < 0.001; Fig. 2a). At the end of 12 weeks, HFD-fed mice gained 22.3 g weight compared with ND-fed mice that gained 15.3 g weight (P < 0.001). Total visceral fat weight of HFD-fed mice was higher than that of ND-fed mice at as early as 2 weeks of the experiment (P < 0.001; Fig. 2b). Food efficiency ratio also increased significantly for HFD-fed mice at all the time points compared with that for NDfed mice (P < 0.001; Additional file 2: Table S2).

Time-course of transcriptional changes in the colon tissue during the development of diet-induced obesity
Two-dimensional hierarchical clustering showed different gene expression patterns at different time points in HFDand ND-fed C57BL/6N mice (Additional file 3: Figure S1). The number of HFD-responsive genes at different time points in the colon of C57BL/6N mice is presented in Table 1. We found that 41, 35, 1, and 33 genes were upregulated and 9, 12, 35, and 378 genes were downregulated at 2, 4, 8, and 12 weeks, respectively, in response to HFD (Table 1). Genes affected by HFD at different time points in the colon tissue of C57BL/6N mice are listed in Additional file 4: Table S3.
Functional analysis was performed using PANTHER classification system to identify biological processes associated with HFD-responsive genes in the colon tissue of mice. The biological processes associated with HFDresponsive genes in the colon tissue of mice are presented in Table 2. At week 12, HFD affected several biological processes, including immunity and defense; nucleoside, nucleotide, and nucleic acid metabolism; signal transduction; and cell cycle ( Table 2). Biological processes associated with HFD-responsive genes at different time points in the colon tissue of C57BL/6N mice are listed in Additional file 5: Table S4.
We also identified HFD-responsive genes showing log 2 fold change of ≥1 or ≤−1 (corresponding to a fold change of ≥2 or ≤−2) at multiple time points (>3 times) over 12 weeks (Table 3). Most HFD-responsive genes were associated with digestive enzymes such as trypsin, carboxypeptidase, and amylase. Overall, these genes were upregulated at weeks 2 and 4 and were downregulated at week 12 in HFD-fed mice compared with those ND-fed mice ( Table 3). Cfd, complement factor D (adipsin), was downregulated at weeks 4, 8, and 12 in HFD-fed mice. Adipsin is suggested to activate an alternative complement pathway for inducing natural defense against infectious agents and red cell lysis and to regulate systemic energy balance [29,30]. A previous study reported that adipsin expression in the small intestine is a potential marker of changes in normal intestinal epithelial differentiation [31]. Pla2g1b, pancreatic phospholipase A2, was upregulated at weeks 2 and 4 and was downregulated at week 12 in HFD-fed mice. Pancreatic phospholipase A2 catalyzes the release of fatty acids from dietary phospholipids. Diet is the ultimate source of arachidonic acid present in cellular phospholipids, which serve as precursors of eicosanoid signaling molecules and are involved in inflammation, cell proliferation, and colorectal carcinogenesis. Arachidonic acid is metabolized by PTGS (COX)/ LOX pathway to prostaglandins and leukotrienes, which are associated with carcinogenesis, specifically of colonic carcinogenesis [32,33].
Cluster and network analyses for identifying candidate early marker genes associated with diet-induced obesity Week 2 Week 4 Week 8 Week 12 Normal diet High-fat diet Total visceral fat weight (mg/g body weight) Fig. 2 Time-course changes in body weight and total visceral fat weight during the development of diet-induced obesity. a Body weight. b Total visceral fat weight. Data are presented as mean ± SEM (n = 10/diet group). *P < 0.01 and **P < 0.001 indicate significant difference compared with mice in the normal diet group (Fig. 3). The number of genes in each cluster was as follows: cluster 1, 45 genes; cluster 2, 32 genes; cluster 3, 17 genes; cluster 4, 44 genes; cluster 5, 35 genes; cluster 6, 24 genes; cluster 7, 8 genes; cluster 8, 78 genes; cluster 9, 103 genes; and cluster 10, 76 genes. Virtual network analysis was performed for genes in cluster 4 that showed the most significant alterations in response to HFD over 12 weeks. The genes in cluster 4 are listed in Table 4. Gene ontology (GO) biological pathway analysis showed that genes in cluster 4 were involved in proteolysis, lipid catabolic process, digestion, defense response, and acute-phase response ( Table 5). Results of the virtual network analysis showed that Apoa4 (apolipoprotein A-IV), Ppap2b (phosphatidic acid phosphatase type 2B), Cel (carboxyl ester lipase), and Clps (colipase, pancreatic) strongly interacted with surrounding genes (Fig. 4). Previous studies have reported that these core genes are involved in pathological changes associated with CRC or obesity [34][35][36]. Results of microarray-based analysis of the expression of these genes were confirmed by performing real-time qPCR at each time point. Overall, the changes in the transcription profiles of Apoa4, Ppap2b, Cel, and Clps determined by real-time qPCR were consistent with the results of microarray analysis (Fig. 5). The direction of change between the two analyses was consistent for the significantly regulated genes except Ppap2b at week 4 (log 2 fold change −0.13) and Cel at week 8 (log 2 fold change 0.04).

Discussion
In the present study, we determined global transcriptional changes at different time points during the development of diet-induced obesity in the colon of mice. We also performed bioinformatics analyses to identify candidate genes that could be used as early markers of obesity-related pathological events. Diet-induced obesity is associated with many chronic diseases, including CRC and IBD. Epidemiological studies have reported a significant association between BMI and colon cancer (HR, 1.05; 95% CI, 1.02-1.09) [37]. Obese C57BL/6 mice develop colonic epithelial hyperplasia, and the risk of colon cancer increases by 42% after long-term (18 months) western-style diet feeding [38,39]. A recent study indicated that HFD increased the number of polyps in the colon and the area of hyperplasia in the mucous membrane tissue of the colon [5]. We previously observed that HFD-fed mice (45% total calories from fat) developed two-times more number of colonic tumors than ND-fed mice possibly because of adipokine-mediated signaling of phosphatidylinositol 3-kinase/Akt pathway [40]. However, limited information is available on mechanisms underlying the associations between obesity and pathophysiological changes in the colon. In the present study, genes showing differential expression in response to HFD were subjected to clustering and networking analyses. Clustering algorithms are frequently used to group genes with similar expression profiles [41]. This facilitates the visualization of coexpressed genes and allows the identification of genes that concurrently respond to stimuli. We clustered genes that were expressed differentially over time into 10 patterns. Of the 10 clusters, cluster 4 included genes that were the most responsive to HFD. Many of these genes were upregulated after the initiation of HFD and were downregulated gradually as the mice became obese. We postulated that these genes could be used as early markers of the initiation of metabolic changes in the colon.
We found that genes in cluster 4 were involved in proteolysis, lipid catabolic process, digestion, defense response, and acute-phase response. These results indicate that HFD Table 2 The biological processes associated with high-fat diet-responsive genes at week 12 in the colon tissue of C57BL/6N mice     upregulated the expression of genes involved in lipid catabolism and that these genes were downregulated over time possibly due to interactions with other compensatory and/ or adaptive mechanisms. Extracellular proteolysis is critical for tumor growth [42]. Trypsin activates protease-activated receptor-2 (PAR-2) and increases COX-2 expression through PAR-2 in Caco-2 cells [43]. These proteolytic activities may promote tumor cell growth and invasion, suggesting that HFD increases the risk of tumor development by facilitating proteolytic activity. Oxidative stress and geneenvironment interactions play a significant role in the development of colon cancer [44]. Oxidative stress results from an imbalance in the production of reactive oxygen species and cellular antioxidant defense system. In the present study, genes associated with defense response tended to be downregulated over time during HFD administration and before colon cancer initiation. This result suggests that continuous HFD administration affects defense mechanisms, which in turn may increase the risk of CRC. For further analysis, genes in cluster 4 were subjected to network analysis by using BIND, DIP, IntAct, Mint, Reactome, CCSB, and HPRD protein-protein interaction databases. Among the genes in cluster 4, four genes showing the most significant relationship with surrounding genes were selected and their expression was verified. Previous studies indicate that these four genes are associated with pathological changes in the colon or with obesity. APOA4 is an intestinally and cerebrally synthesized antiatherogenic plasma apolipoprotein that functions as a satiety factor and anti-inflammatory protein. Intestinal APOA4 synthesis is stimulated by fat intake and is attenuated by intravenous leptin infusion, indicating a close association between fat and energy intake [45]. Apoa4 expression is altered along with that of other genes involved in epithelial junctional integrity in the intestinal mucosa of patients with IBD [46]. APOA4 stabilizes adherent junctions by interacting with α-catenin and may be involved in the maintenance of junctional integrity. Epithelial tight junctions form a barrier to prevent the movement of pathogens, toxins, and allergens from the intestinal lumen into the tissue, and disruption of these tight junctions may play an important role in the pathogenesis of gastrointestinal diseases [47,48].
Lipid phosphate phosphatase 3 (LPP3) encoded by Ppap2b is an integral membrane glycoprotein that catalyzes the dephosphorylation of several bioactive lipid mediators, including lysophosphatidic acid, sphingosine 1-phosphate, and phosphatidic acid. Moreover, LPP3 functions as a cell-associated integrin ligand [49,50]. A recent study reported that LPP3 does not promote tumor formation but amplifies β-catenin signaling and cyclin-D1 activity to potentiate the growth of SW480 colon carcinoma [51]. Aberrant activation of PI3K/Akt/ mTOR and MAPK/ERK pathways may induce colon tumor growth and progression by increasing β-catenin and cyclin-D1 expression [52,53].  Carboxyl ester lipase (CEL) encoded by Cel is a 74-kDa lipolytic enzyme that hydrolyzes cholesteryl esters, triacylglycerol, phospholipids, and lysophospholipids [54,55]. This enzyme is synthesized in acinar cells of the pancreas and is stored in zymogen granules. Upon food ingestion, CEL is released into the intestinal lumen where it constitutes 1-5% of total proteins in the pancreatic juice [56]. CEL plays a significant role in catalyzing the absorption of cholesteryl esters from the intestinal lumen and in promoting the formation of large chylomicron [57,58]. A recent study reported that Celknockout mice developed a mild diabetic phenotype after the administration of 60% HFD [59]. Since insulin resistance is a risk factor of colon cancer, differential expression of Cel in obese animals may be responsible for the association of obesity with the pathophysiological changes in the colon.
Clps encodes colipase that is secreted from the exocrine pancreas into the gastrointestinal tract [60]. Colipase may interact with pancreatic triglyceride lipase to facilitate the digestion of dietary fats. HFD-fed Clps −/− mice develop hyperphagia, and procolipase performs essential functions by regulating body weight set point [61]. Also, Clps genetic variability is associated with insulin secretory function in non-diabetic humans, suggesting that Clps is a novel candidate gene associated with the development of type 2 diabetes [36]. Regulation of insulin secretion is important for metabolic homeostasis in various tissues, including the liver, adipose tissue, and colon [62]. Therefore, Clps expression would be a potential early marker of the development of obesity, insulin resistance, and/or colon cancer.

Conclusions
In conclusion, our data indicate that a few genes primarily involved in lipid metabolism play a functional role in diet-induced pathological changes in the colon. Genome-wide analyses performed in the present study provide new insights on selecting novel genes that may Ethics approval This study was approved by the Institutional Animal Care and Use Committee of Yonsei University in Seoul, Republic of Korea (Permit #: 2011-0061). All institutional and national guidelines for the care and use of laboratory animals were followed.