Effect of myeloid differentiation primary response gene 88 on expression profiles of genes during the development and progression of Helicobacter-induced gastric cancer

Background Gastric cancer is one of the most common and lethal type of cancer worldwide. Infection with Helicobacter pylori (H. pylori) is recognized as the major cause of gastric cancer. However, it remains unclear the mechanism by which Helicobacter infection leads to gastric cancer. Furthermore, the underlying molecular events involved during the progression of Helicobacter infection to gastric malignancy are not well understood. In previous studies, we demonstrated that that H. felis-infected Myd88 −/− mice exhibited dramatic pathology and an accelerated progression to gastric dysplasia; however, the MyD88 downstream gene targets responsible for this pathology have not been described. This study was designed to identify MyD88-dependent genes involved in the progression towards gastric cancer during the course of Helicobacter infection. Methods Wild type (WT) and Myd88 deficient mice (Myd88 −/−) were infected with H. felis for 25 and 47 weeks and global transcriptome analysis performed on gastric tissue using MouseWG-6 v2 expression BeadChips microarrays. Function and pathway enrichment analyses of statistically significant, differential expressed genes (p < 0.05) were performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) online tools. Results Helicobacter infection affected the transcriptional profile of more genes in Myd88 −/− mice compared to WT mice. Infection of Myd88 −/− mice resulted in the differential expression of 1,989 genes at 25 weeks (1031 up and 958 downregulated). At 47 weeks post-H.felis infection, 2,162 (1140 up and 1022 downregulated) were differentially expressed. The most significant differentially upregulated gene during Helicobacter infection in Myd88 −/− mice was chitinase-like 4 (chil4), which is involved in tissue remodeling and wound healing. Other highly upregulated genes in H. felis-infected Myd88 −/− mice included, Indoleamine 2,3-Dioxygenase 1 (Ido1), Guanylate binding protein 2 (Gbp2), ubiquitin D (Ubd), β2-Microglobulin (B2m), CD74 antigen (Cd74), which have been reported to promote cancer progression by enhancing angiogenesis, proliferation, migration, metastasis, invasion, and tumorigenecity. For downregulated genes, the highly expressed genes included, ATPase H+/K+ transporting, alpha subunit (Atp4a), Atp4b, Mucin 5 AC (Muc5ac), Apolipoprotein A-1 (Apoa1), and gastric intrinsic factor (Gif), whose optimal function is important in maintaining gastric hemostasis and lower expression has been associated with increased risk of gastric carcinogenesis. Conclusions These results provide a global transcriptional gene profile during the development and progression of Helicobacter-induced gastric cancer. The data show that our mouse model system is useful for identifying genes involved in gastric cancer progression. Electronic supplementary material The online version of this article (doi:10.1186/s12885-017-3114-y) contains supplementary material, which is available to authorized users.


Background
Gastric cancer is one of the most common causes of cancer-related death worldwide with an estimated 738,000 deaths each year [1]. Recently, H. pylori was recognized as the foremost cause of gastric cancer [2][3][4][5][6][7]. With an estimated half of the world's population being infected, Helicobacter infection contributes significantly to the worldwide gastric cancer burden [7,8]. Recognition of the factors leading up to the development and progression towards gastric cancer are critical in determination of cancer pathology. H. pylori-induced gastric carcinogenesis involves a multistep progression from normal gastric mucosa to superficial gastritis, chronic gastritis, atrophic gastritis, metaplasia, dysplasia, and finally gastric carcinoma [8,9]. Molecular events associated with disease progression to gastric malignancy have not been elucidated. Considerable amount of confirmatory evidence shows that host immune response to H. pylori is crucial in determining gastric cancer predisposition [10][11][12]. We have previously shown that a key signal transduction adaptor protein, myeloid differentiation primary response gene 88 (MyD88), regulates Helicobacter-induced gastric cancer progression in a mouse model of gastric cancer [13]. We demonstrated that H. felis-infected MyD88 deficient (Myd88 −/− ) mice exhibited severe gastric pathology and an accelerated progression to gastric dysplasia compared to wild type (WT) mice [13] However, the MyD88-dependent gene responsible for this pathology were not described.
MyD88 is a key adaptor molecule that is crucial in mediating innate immune signals from members of the toll-like receptor (TLR) and interleukin-1 (IL-1)/IL-18 families leading to downstream activation of nuclear factor (NF)-κB [14][15][16]. Consistent with involvement in these inflammatory pathways, MyD88 signaling has been associated with cancer progression, which stems from the understanding that inflammation is linked to cancer promotion [17,18]. Studies on the role of MyD88 cancer progression have been the subject of recent intense investigations. However, the data are contradictory, which indicate that the role of MyD88 in the development and progression of inflammation-associated cancers is complex [19]. Several studies using genetic or chemical carcinogenesis models involving Myd88 deficient mice have shown MyD88 to either promote [20][21][22][23][24][25][26][27] or suppress [13,[28][29][30][31][32][33][34] tumor development. The complex role of MyD88 in carcinogenesis is best typified by studies in colon cancer models [22,24,29,35] showing contradictory roles in the same tissue. The mechanistic basis for these opposing observation is still not fully understood and could be due to many factors including, the type of inflammation, the extent of tissue damage, and immune response elicited [35]. Further, the MyD88 dependent genes in this accelerated progression to dysplasia remain unknown. Therefore, this study was performed to identify potential genes involved in the accelerated progression of gastric cancer.

Gene expression and analysis
Prior to differential gene analysis, all data from 23,015 genes with a standard deviation of less than 0.1 were used for multiple dimensional scaling (MDS) analysis ( Fig. 1) to verify that Myd88 −/− and WT samples were differentiated according to gene expression in each sample with a relative p-value. Each sample is represented with distance between each one reflecting their approximate degree of correlation [36]. All genes included in the analysis had a minimum standard deviation of less than 0.1. The analysis showed that all uninfected mice were clustered together irrespective of genetic background or time point. For infected mice, WT and Myd88 −/− mice clustered distinctively separate indicating differential expression of their genes.
Statistical analysis of all 23,015 genes that went through the filtering process identified a total of 286 genes in WT and 4,151 in Myd88 −/− mice in response to H. felis infection with more genes differentially expressed at 47 than 25 weeks (Table 1). Comparing the number of upregulated genes between Myd88 −/− and WT at 47 weeks post infection, there were more upregulated genes (1140) in Myd88 −/− mice compared to WT mice (189 genes). A similar trend was observed for upregulated genes at 25 weeks and for downregulated genes at both time points, with more genes differentially expressed in Myd88 −/− than WT mice in response to H. felis infection. The number of differentiated genes at each time point in comparison to uninfected controls is illustrated in Fig. 2. Most genes overlapped between time points, however, there were a substantial number of genes that were unique to each set of time points that were differentially regulated (Fig. 2). Scatterplot depiction of differentially expressed genes shows significantly (P < 0.05) up-and downregulated genes in Myd88 −/− mice in response to H. felis infection (Fig. 3). A majority of genes were altered after 25 weeks of H. felis infection. Some of the new additional genes at 47 weeks post-H. felis infection included the ring finger protein 213 (Rnf213) and Furin, which have been reported to be involved in angiogenesis [37,38] and cancer progression [39], respectively indicating their role in advanced stages of cancer progression. Tables 2 and 3 show a list of the top 50 up-and downregulated genes in Myd88 −/− mice at 25 and 47 weeks post-H. felis infection compared to uninfected controls. The most highly upregulated gene during H. felis infection in Myd88 −/− mice included Chitinase-like (chil4), which is involved in tissue remodeling and wound healing [40][41][42]. Many of the upregulated genes in both 25 weeks and 47 weeks post-H. felis infection involved genes in the H2 Complex (murine major MHC), particularly the class I heavy chains, H2-K and H2-D. The light chain for this MHC complex consists of the β 2 -Microglobulin (B2m) [43]. MHC class II antigen presentation including the CD74 antigen (Cd74) was another gene that was upregulated in response to H. felis infection. High expression of Cd74 has been linked to chronic inflammation and carcinogenesis in the gastrointestinal tract [44]. Another highly expressed gene was Indoleamine 2,3-Dioxygenase 1 (Ido1), which is suggested to play a role in immune tolerance and high expression in colorectal cancer and is correlated with a poor clinical outcome (reviewed in [45]. The entire list of altered genes in response to infection with H. felis including those in WT mice have been submitted and can be uploaded as an excel file in Additional files 1: Table  S1 (inf vs uninfectd charts.xlsx).

Analysis of differentially expressed genes
STRING summary networks depicting protein-protein interactions among the top differentially expressed genes (DEGs) for both up-and downregulated genes in Myd88 −/− mice are shown in Figs. 4 and 5 for 25 and 47 weeks, respectively. Thicker lines connecting the genes indicate a stronger association between the genes. A confidence score of at least 0.70 (high) was used. One of the key central nodes in the top DEGs in Myd88 −/− mice at both 25 and 47 weeks post-infection was guanylate-binding protein 2 (Gbp2), which is considered a potential marker for esophageal squamous cell carcinoma [46].

Functional enrichment analysis of differentially expressed genes (DEGs)
To gain insights into the biological meaning and function of the differentially expressed genes, enrichment analysis was performed using the database for annotation, visualization and integrated discovery (DAVID) online analytical tools [47][48][49]. Annotation according to tissue expression, molecular function, cellular component and biological processing was done using Gene Ontology (GO) [49]. Enrichment analysis was performed to identify pathways, processes and gene categories that are over-represented in the list of DEGs compared to     the mouse genome. GO clustering analysis for biological processes showed that responses related to immune system processes were the most upregulated enriched process in Myd88 −/− mice at both 25 and 47 weeks in response to H. felis infection (Fig. 6a and c). For molecular functions, antigen and protein complex binding were the most enriched processes ( Fig. 6b and d). Downregulated enriched processes in Myd88 −/− mice in response to H. felis infection are presented as supplementary data (Additional files 2: Figure S1). A summary of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotation shown in Fig. 7 revealed that the most enriched pathway was antigen processing and presentation. This pathway shared the majority of its genes with the other major pathways and completely engulfed the other pathways by sharing more than 90% of the genes annotated. A breakdown of up-and downregulated KEGG pathway at 25 and 47 weeks in Myd88 −/− mice in response to H. felis infection is presented as supplementary data (Additional files 3: Figure S2).

Discussion
Gastric cancer develops and progresses through a stepwise sequence of events from inflammation to atrophy, metaplasia, dysplasia, and finally to gastric cancer [50]. We previously demonstrated using a mouse model of gastric that mice deficient in MyD88 signaling exhibited dramatic pathology and an accelerated progression to gastric neoplasia in response to H. felis infection [13]. In the present study, we used microarray gene expression analysis to identify the genes involved in this progression to gastric neoplasia. Although previous studies have investigated differential gene expression in mice stomachs in response to Helicobacter infection, most have focused on H. pylori [51][52][53], which does not result in neoplastic changes in mice [54,55]. The few studies that have examined gene expression profiles in mouse model of gastric cancer have used the insulin-gastrin (INS-GAS) transgenic gastric cancer mouse model [56,57]. These mice have been shown to spontaneously develop gastric cancer even in the absence of Helicobacter infection [58]. We have previously reported that Myd88 −/− mice do not exhibit abnormal pathology in the absence of Helicobacter infection [13]. The global transcriptional profiling of mouse gastric tissue identified a large number of significant differentially expressed genes in H. felis-infected Myd88 −/− mice compared to H. felis-infected WT mice. The most over expressed gene in Myd88 −/− mice during H. felis infection at 25 weeks was Chil4. Chitinase like proteins (CLPs) have been studied in relation to other cancers yet little has been investigated in relation to gastric cancer except for our present study and a couple other studies [57,59]. Upregulation of CLPs has been shown in a number of human cancers including brain, bone, breast, ovaries, lung, prostate, colon, thyroid, and liver [41,60]. For gastric cancer studies, chitinase protein 3 like 1 (Chil1) was upregulated in INS-GAS mice infected with H. felis, [57]. In our present study, Chil1 was not upregulated in response to H. felis infection. However, in addition to upregulation of Chil4, another CLP, Chil3 was also significantly over expressed in H. felis-infected mice at both 25 weeks (p = 0.03) and 47 weeks (p = 0.01) (gene not listed in Tables 2 and 3, only the top 50 are listed). An abundant over expression of Chil1, Chil4 as well as Chil3 has been reported in early preneoplastic stage in the epidermis [61]. Overall, CLPs have been implicated to play a role in chronic inflammation, tissue remodeling, and wound healing [40]. Up-regulation of genes involved in tissue remodeling is noteworthy because chronic inflammation and subsequent damage to the gastric epithelium has been suggested to play an important role in cancer development  and progression [62]. During chronic inflammation, the resulting prolonged tissue damage creates a loss of control over normal tissue repair mechanisms resulting in persistent hyper-tissue repair, which is accompanied with sustained proliferation [63] and ultimately advancing to precancerous lesions. Lost tissue is then replaced with stem and progenitor cells that are under a continuous stimulus of proliferation, leading to the accumulation of replacement cells with dysregulated and altered signaling pathways [63]. Further, studies investigating associations between chronic inflammation, tissue repair and carcinogenesis highlight the potential of these cellular changes in inducing both pro-oncogenic and tumor suppressor pathways [62,[64][65][66][67]. Our study, in addition to the work done by Li et al. [59] and Takaishi and Wang [57] show a need for further investigation into the role of CLPs in gastric carcinogenesis.
Other upregulated genes included Cd74, B2m, and interferon (IFN) induced genes such as GTPases (interferon gamma induced GTPase, lgtp, immune mediated GTPase family M member 2, lrgm2), Guanylate binding protein 2 (Gbp2), and transcription factor interferon regulatory factor 1 (Irf1). Cd74 or invariant chain (Ii) protein is a chaperone molecule responsible for regulating antigen presentation of MHC II molecules. It has been linked to chronic inflammation and carcinogenesis in the gastrointestinal tract [44]. Further, Cd74 was also shown to play a role as a receptor for migration inhibitory factor (MIF), a molecule reported to have procarcinogenic effects on gastric epithelial cells [68]. IFNs are known to activate signal transducer and activator of Genes with no known annotated name were excluded from the analysis. Genes were considered to be statistically significant when a threshold adjusted p value < 0.05 and Log FC > 2 were reached  [69,70] signaling leading to epithelial proliferation and inhibition of apoptosis [71,72].
Currently not much is known about IFNs in gastric cancer. Ubiquitin D (Ubd), which is associated with progression of colon cancer [73], was also significantly expressed genes in Myd88 −/− in response to H. felis infection.
For downregulated genes, the significantly expressed ones included, ATPase H+/K+ transporting, alpha subunit (Atp4a), Atp4b, Mucin 5 AC (Muc5ac), apolipoprotein A-1 (Apoa1), and gastric intrinsic factor (Gif ). The genes, Atp4a and Atp4b encode gastric H+/K + − ATPase alpha and beta subunits, respectively. Gastric H +/K + − ATPase alpha and beta subunits are expressed in parietal cells [74] and their loss has been associated with gastric dysplasia [58]. We observed downregulation of Atp4a and Atp4b in response to infection with H. felis, which may represent a loss of parietal cells that has been shown to precede gastric dysplasia. These results are in line with those observed in another fast progressing gastric cancer model involving the use of INS-GAS mice [58]. Muc5ac, which encodes gastric M1 mucin [75] has been reported to play a role in gastric carcinogenesis [76] was also downregulated in response to infection with H. felis in Myd88 −/− mice. Progression of gastric lesions has been reported to be associated with the gradual decrease in expression of Muc5ac [57,[77][78][79] followed by the transformation of the gastric epithelium [80] resulting in gastric dysplasia. Another highly downregulated gene we found in Myd88 −/− mice infected with H. felis was Apoa1, which was also reported to be downregulated in a fast progressing gastric cancer mouse model [57]. Proteomics approach in human gastric cancer also showed downregulation of Apoa1 [81], but its role in gastric carcinogenesis is unknown. In the lung, downregulation of Apoa1 was associated with an increased risk of lung cancer [82]. Data from Apoa1-deficient mice suggest antitumorigenic properties of Apoa1 via modulation of the immune system [83]. Gastric intrinsic factor (Gif), another downregulated gene is secreted by parietal cells and is required for Vitamin B12 absorption [84]. Concomitantly, the downregulation of Gif results in vitamin B12 deficiency (pernicious anemia). Gastric intrinsic factor was downregulated in H. felis-infected Myd88 −/− mice at both 25 and 47 weeks post-infection. A previous study using INS-GAS mice infected with H. felis also reported downregulation of Gif [57]. In human gastric cancer, Gif was one of the genes found by SAGE analysis to be downregulated [85]. Further, studies have shown an increased risk of gastric cancer in pernicious anemic patients [86]. We found a number of new genes in our fast progressing gastric mouse model, i.e., Myd88 −/− mice infected with H. felis. The genes included up-and downregulated genes, which had not been previously linked to Helicobacter-related gastric carcinogenesis including BPI fold containing family B member 1 (Bpifb1) and proteasome subunit beta 8 (Psmb8). These genes have been linked to cancer-related processes including, apoptosis and in some cases other cancers as well as prognosis indicators [87,88]. Psmb8 was shown to be significantly up regulated in cancers such as bladder, breast, kidney, lung, uterine, and head and neck [89]. A recent study by Kwon, et al. [90], which was published during the writing of our manuscript reported that Psmb8 may be a potential marker for prognosis in gastric cancer. Bpifb1 may be involved in the innate immune response particularly in response to bacterial exposure. The protein encoded for by Bpifb1 binds bacterial lipopolysaccharide (LPS) as well as modulates the cellular response to LPS [91]. Bpifb1 has been found to be overexpressed in mucous cells of salivary gland tumors of papillary cystadenocarcinoma [87]. Future studies using a gastric culture organoid system will validate these genes and some of the important novel genes we identified for their role in rapid progression of Helicobacter-induced gastric cancer.

Conclusions
In this study, we have identified genes that are involved in the rapid progression of Helicobacter-induced gastric cancer that are also potentially regulated by MyD88. The identification of these important genes could potentially serve as targets for disease prevention. In addition, we show that our model is a useful mouse model system to identify genes involved in gastric cancer progression.

Animals
Six-to ten-week-old wild type (WT) and MyD88 deficient (Myd88 −/− ) mice in the C57BL/6 background were Biological processes are depicted in figures A and C while molecular functions are depicted in B and D. For the biological processes, the top 20 processes are shown. All scores depicted are relative scores for number of genes in each function/ process relative to the number of total genes entered into the STRING database used in this study. WT mice were purchased from The Jackson Laboratory (Bar Harbor, ME). Myd88 −/− mice were from our breeding colony originally provided by Dr. Akira (Osaka University, Japan). All mice were housed together before infection with H. felis and for the duration of the study. The Institutional Animal Care and Use Committee at the University of California, San Diego, approved all animal procedures and performed using accepted veterinary standards.

Mouse infections
A well-characterized cancer mouse model, which involves infecting C57BL/6 mice with H. felis (strain CS1), a close relative of the human gastric pathogen H. pylori was used in this study. Mice were inoculated with 10 9organisms in 300 μL of BHI by oral gavage three times at 2-day intervals as previously described [13,92]. Control mice received BHI only. At 25 and 47 weeks postinfection, mice were euthanized and the stomachs removed under aseptic conditions and processed for assessment of gene expression.

RNA extraction and oligonucleotide microarray hybridization
Total RNA was extracted from gastric tissue obtained from H. felis-infected and uninfected WT and Myd88 −/− mice. Stomach tissue sections of 12 mice i.e., two uninfected controls and four infected WT (6) and MyD88 −/− (6) per time point (25 and 47 weeks) were analyzed. RNA was extracted from each section using the RNeasy miniprep kit (Qiagen) according to the manufacturer's instructions followed by digestion with DNase 1 to remove genomic DNA. RNA concentration was determined using a NanoDrop spectrophotometer (NanoDrop Technologies, Inc., Waltham, MA). Double stranded cDNA and biotin-labeled cRNA were synthesized following the recommended Illumina protocol. Integrity of purified cRNAs was assessed on an Agilent 2100 Bioanalyzer prior to hybridization. 1.5 μg of labeled cRNA was hybridized to MouseWG-6 v2 Expression BeadChips genome wide arrays, which analyzes 25,600 transcripts (Illumina, San Diego, CA) using recommended Illumina reagents and protocols.

Identification of differentially expressed genes
Probe profiles (each row corresponding to a given probe and different columns for each sample) were exported from Genome Studio v1.8 (Illumina). The resultant tabdelimited file (Additional file 4, Probe_Profile_f.txt) was used as input for the Bioconductor lumi v2.18 R package (http://bioconductor.org/packages/release/bioc/html/lumi. html) [93]. Sample information is provided in Additional file 5 (sampleInfotxt.txt). Following quality assessment [density plot of intensity, cumulative density (CD), MA (transformed data onto log ratios and mean average) and pairwise MDS plots], the data were transformed using the vst (variance-stabilizing transformation) algorithm [94] then normalized using the Robust Spline Normalization (rsn) algorithm [94]. A second round of quality control of the normalized data was done to ensure data quality. After normalization genes with minimal variance across samples and genes that were not expressed in any sample (determined by detection calls) were removed prior to differential gene expression using linear models as implemented in limma [95]. We compared gene expression profiles of uninfected and infected mice at 25 and 47 weeks post infection. All mice used were on the C57BL/6 background. In addition, after normalization, significantly expressed genes were identified through volcano plot generation on the R -platform. Genes exhibiting statistically significant differential mean expression values (p < 0.05) were subjected to hierarchical clustering.

Enrichment analysis and pathway generation
To analyze KEGG pathway enrichment, official gene symbols for differentially expressed genes were submitted to DAVID (Database for Annotation, Visualization and Integrated Discovery) (http://david.abcc.ncifcrf.gov). This tool allows for identification of over expression within a set of genes as compared to the entire genome of a specific animal [47]. In addition, it allows for the demonstration of Gene Ontology (GO) terms as well as pathway and functional processes (molecular functions, biological expression, and cellular expression) enrichment. Functionally enriched gene networks, KEGG Pathways and ontology terms were identified. The mouse, Mus musculus complete genome was used as background genes. All functionally annotated genes presented in gene networks or pathways had p values < 0.05.

Gene network analysis
Gene networks were built using data on protein-protein interactions from EMBL STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) database (http:// string-db.org). This search tool is used to identify interactions correlated to expression data and/or literature citations among other criteria [96]. All gene connections created using STRING had a combined confidence score higher than 0.7 as previously described [97], Garcia-Alonso, 2014 #613}.

Statistical analysis
Statistics were done on the R platform using the Limma package from Bioconductor. To control for multiple testing, the False Discovery Rate (FDR) method was used with a cutoff for statistical significance of P values of < 0.05 and a log fold expression of 2. Differentially expressed genes were determined at 25 and 47 weeks after removing background differences from both Myd88 −/− and WT mice by comparing infected to uninfected mice in the same background.