Genetic Predisposition to Hepatocarcinogenesis in Inbred and Outbred Mouse Lines Selected for High or Low Inflammatory Response

de Carvalho, Lilian Rego; Borrego, Andrea; Jensen, José Ricardo; Cabrera, Wafa Hanna Koury; Santos, Aline Marques; Ribeiro, Orlando Garcia; Starobinas, Nancy; De Franco, Marcelo; Dragani, Tommaso A.; Manenti, Giacomo; Ibañez, Olga Célia Martinez

doi:https://doi.org/10.1155/2019/5298792

Journal of Immunology Research

On this page

Abstract Introduction Methods Results Discussion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 5298792 | https://doi.org/10.1155/2019/5298792

Genetic Predisposition to Hepatocarcinogenesis in Inbred and Outbred Mouse Lines Selected for High or Low Inflammatory Response

Lilian Rego de Carvalho,¹Andrea Borrego,¹José Ricardo Jensen,¹Wafa Hanna Koury Cabrera,¹Aline Marques Santos,²Orlando Garcia Ribeiro,¹Nancy Starobinas,¹Marcelo De Franco,^1,3Tommaso A. Dragani,⁴Giacomo Manenti,⁴and Olga Célia Martinez Ibañez¹

Academic Editor: Theresa Hautz

Received31 Oct 2018

Accepted24 Jan 2019

Published31 Mar 2019

Abstract

AIRmax and AIRmin mouse strains phenotypically selected for high and low acute inflammatory responsiveness (AIR) are, respectively, susceptible or resistant to developing hepatocellular carcinoma (HCC) induced by the chemical carcinogens urethane and diethylnitrosamine (DEN). Early production of TNF-α, IL-1β, and IL-6 in the liver after DEN treatment correlated with tumor development in AIRmax mice. Transcriptome analysis of livers from untreated AIRmax and AIRmin mice showed specific gene expression profiles in each line, which might play a role in their differential susceptibility to HCC. Linkage analysis with SNP markers in F2 (AIRmax×AIRmin) intercross mice revealed two quantitative trait loci (QTL) in chromosomes 2 and 9, which are significantly associated with the number and progression of urethane-induced liver tumors. An independent linkage analysis with an intercross population from A/J and C57BL/6J inbred mice mapped regions in chromosomes 1 and 7 associated with the progression of urethane-induced liver tumors, evidencing the heterogeneity of HCC genetic control.

1. Introduction

Many human cancers such as those affecting the liver are etiologically related to processes of inflammation and/or chronic infection. Human hepatocellular carcinoma (HCC) is the third most common cause of cancer mortality and the seventh in terms of cancer incidence worldwide [1]. The incidence of HCC varies widely according to racial and ethnic groups and to geographic location, as a result of the influence of genetic factors and of regional variations in exposure to risk factors [2]. Among the most important risk factors of liver cancer are hepatitis B, chronic hepatitis C virus infection, hereditary hemochromatosis, and cirrhosis of almost any cause. Tobacco and alcohol abuse, environmental toxins and dietary factors, diabetes mellitus, nonalcoholic fatty liver disease, alpha-1 antitrypsin deficiency, and autoimmune hepatitis are all associated with increased risk for HCC [3]. The incidence of HCC is approximately three times higher in men than in women, and most of the experimental models confirm this sex difference [4].

We developed two mouse lines that differ widely in inflammatory responsiveness to investigate the genetic control of inflammatory processes. These lines were produced by bidirectional phenotypic selection on the basis of the intensity of the local acute inflammatory response (AIR) to a nonimmunogenic inflammatory stimulus (Bio-Gel P-100 polyacrylamide beads) and were called AIRmax (for high response) and AIRmin (for low response) [5]. The acute inflammatory response was measured by the number of infiltrated cells and protein concentration in the 24-hour inflammatory exudate induced by subcutaneous Bio-Gel injection. The foundation population was constituted by mixing 8 inbred mouse laboratory lines, and after about 25 generations of selective breeding, the phenotypic divergence between AIRmax and AIRmin lines reached 20 to 30 times for the number of infiltrating leukocytes (85% neutrophils) and 2.5-fold for exudate protein concentration [6].

At the end of the selection process, these mouse lines also differed in susceptibility to several chemically induced tumors. AIRmax are more susceptible than AIRmin to colon carcinogenesis [7] while AIRmin are more susceptible to skin, lung, and kidney carcinogenesis [8–10]. The two strains have distinct gene expression profiles in the bone marrow and in the normal lung [11, 12] which correlated with their differential inflammation capacity and tumor susceptibility, respectively.

To better investigate the genetic relationship between inflammation and cancer susceptibility, we carried out linkage analyses using a large pedigree of about 700 F2 (AIRmax×AIRmin) intercross mice. In this experiment, newborn mice were treated with urethane, the inflammatory reaction to subcutaneous Bio-Gel injection was measured 3 months later, and tumors in internal organs were analyzed at 7 months of age. Genetic linkage analysis results of lung and kidney tumor susceptibility were published [10, 13], and here, we show the results of liver carcinogenesis.

Urethane treatment induced large and multiple liver tumors in all male AIRmax mice, whereas few AIRmin male mice developed small hepatic lesions. The liver tumor load in the whole pedigree of F2 (AIRmax×AIRmin) mice was used in the present study for linkage analysis with SNPs. Similar analysis was carried out in the fourth generation of an advanced intercross population named ABF4 between A/J and C57BL/6J inbred mice.

The transcriptome of normal livers as well as the response of AIRmax and AIRmin mice to another genotoxic chemical carcinogen, N-nitrosodiethylamine (DEN), was also analyzed in order to confirm the divergent susceptibility of these two mouse lines to liver carcinogenesis and to assess the role of tissue cells in the control of immune responses against tumors.

2. Methods

2.1. Mice and Tumor Induction

The formal stock designations of the mouse lines used are Ibut:AIRH and Ibut:AIRL at the ILAR (Institute for Laboratory Animal Research, National Academy of Sciences, USA), but they are referred to as AIRmax and AIRmin in this paper as well as in previous publications. These lines and crosses for production of the F2 (AIRmax×AIRmin) population were developed and maintained at the animal facilities of the Laboratory of Immunogenetics of the Butantan Institute (São Paulo, Brazil). Male and female 15-day-old mice were used in experiments for long-term studies of hepatic cancer development and also for analysis of cytokine production in the liver after injection of the carcinogen diethylnitrosamine (N-nitrosodiethylamine (DEN), Sigma-Aldrich, St. Louis, MO) (25 mg/kg bw ip). Urethane (carbamic acid ethyl ester, Sigma-Aldrich, St. Louis, MO) (300 mg/kg bw ip) was injected in groups of male and female mice seven days after birth. Mice were euthanized at 34 weeks after treatment for analysis of tumor formation in internal organs. Liver tumor incidence, multiplicity, and volumes (calculated by measuring the mean diameters of the lesions) were recorded. FFPE liver fragments were cut into 4 μm sections and stained with hematoxylin and eosin.

We used, for linkage analysis, the liver carcinogenesis results of 693 F2 (AIRmax and AIRmin) intercross mice and of 183 F4 intercross mice (ABF4), generated by mating C57BL/6 (B) with A/J (A) mice (details in [14]). F2 (AIRmax×AIRmin) mice were treated with one dose of urethane (300 mg/kg bw ip) 7 days after birth, and ABF4 mice were treated with a single intraperitoneal injection of urethane (1 g/kg body weight) at 4 weeks of age; the mice were killed at 34 and 40 weeks after treatment, respectively.

All procedures were approved by the Institutional Animal Care and Use Committees of Butantan Institute and Istituto Tumori, and all animals received humane care, according to the criteria outlined in the “Guide for the Care and Use of Laboratory Animals” prepared by the National Academy of Sciences and published by the National Institutes of Health.

2.2. Cytokine Analysis

Groups of 6 AIRmax and 6 AIRmin mice were euthanized at 2, 8, and 24 h after DEN injection to analyze the acute cytokine production in the liver. Total liver homogenates were prepared by homogenizing the whole liver from each mouse in 1500 μl of sterile PBS, and 100 μl of supernatant from each liver homogenate was analyzed with the OptEIA Mouse ELISA Sets (BD Biosciences, San Jose, CA, USA) for IL-6, TNF-α, and IL-1β according to the manufacturer’s instructions. Cytokine concentration was calculated as pg/ml for each sample.

2.3. Genome-Wide SNP Genotyping

Genomic DNA was extracted from tail tips using the E.Z.N.A.® Tissue DNA Kit (Omega Bio-Tek, Norcross, USA) and quantified using the Quant-iT™ PicoGreen® dsDNA Assay Kit (Invitrogen, Carlsbad, CA). SNP genotyping was carried out in all mice of both the F2 (AIRmax×AIRmin) intercross pedigree and intercrossed ABF4 mice using the 1449-SNP loci mouse medium density linkage panel (Illumina Inc., San Diego, CA) at a density of approximately three SNPs per 5 Mb across the whole genome as described in [15]. QTLs affecting the liver tumor phenotypes were mapped through genome-wide linkage analyses by interval mapping using GridQTL version 3.1.0 [16] that uses a linear model to fit phenotype data according to genotypes. Additive and dominant effects at the QTL were included with other explanatory variables such as sex and family. The significance thresholds of phenotype-genotype associations were estimated by genome-wide permutation analysis ( or 10,000 permutations).

2.4. Affymetrix Microarray Analysis

Total RNA from the livers of untreated adult mice were isolated using the RNAspin Mini Kit (GE Healthcare, Buckinghamshire, UK). The concentrations and quality of the purified RNA were determined in the NanoDrop (Thermo Fisher Scientific) and in the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA) equipment, respectively. The Mouse Gene 1.0 ST Array (Affymetrix, Santa Clara, CA) with 28,853 well-annotated genes was used for transcriptome analysis. The reactions of the array were performed following the manufacturers’ protocols at the AFIP-UNIFESP Molecular Facility, Federal University of São Paulo, Brazil. Four biological replicates for each group were run individually. Differentially expressed genes (DEGs) were identified by unpaired one-way ANOVA within the Transcriptome Analysis Console 3.0 (TAC) software (Affymetrix) with (false discovery rate).

2.5. Real-Time Quantitative RT-PCR

Microarray data were validated by quantitative real-time PCR. Amplification mixtures containing 1 μl template cDNA, 12.5 μl SYBR Green PCR Master Mix (Invitrogen, Carlsbad, CA, USA), and 0.3 μM specific PCR primers were run on a StepOne (Applied Biosystems, Foster City, CA) machine. The mouse Hprt (hypoxanthine-guanine phosphoribosyltransferase) and Ppia (cyclophilin A) genes were used as housekeeping controls for possible differences in cDNA amounts. Target and housekeeping genes were amplified with the following primers: Hprt F:5-CGTCGTGATTAGCGATGATGA-3 and R:5-CCAAATCCTCGGCATAATGATT-3; Ppia F:5-AGCGTTTTGGGTCCAGGAAT-3 and R:5-AAATGCCCGCAAGTCAAAAG-3; H2Ea prime time predesigned primers set-Mm.PT.56a.42692423.g; Defb1: prime time predesigned primers set-Mm.PT.56a.5676756; Vnn3 prime time predesigned primers Mm.PT.56a.9417908; TNF-α F:5-TCTCATCAGTTCTATGGCCC-3 and R:5-GGGAGTAGACAAGGTACAAC-3; Il6 F:5-GTTCTCTGGGAAATCGTGGA-3 and R:5-TGTACTCCAGGTAGCTATGG-3; Cyp1a1: F:5-TGGAGACCTTCCGGCATTC-3 and R:5-GCCATTCAGACTTGTATCTCTTGTG-3; and Cyp1b1 F:5TGGCTGCTCATCCTCTTTAC-3 and R:5-AGGTTGGGCTGGTCACTCAT-3.

2.6. Statistical Analysis

Differences between means were determined by Student’s -test or analysis of variance (ANOVA); differences between proportions of affected and nonaffected animals were analyzed by Fisher’s exact test; differences between parameters of tumor growth were analyzed by nonparametric Mann-Whitney test. Pearson correlation analysis was carried out to compare the microarray with qPCR results. was established as the minimum level of significance.

Gene expression of each array was log-transformed to approximate Gaussian distributions and standardized over the array to adjust for systematic differences. The Significance Analysis of Microarray (SAM) software (two-class unpaired, %) was used to detect differentially expressed genes (DEGs), in four biological replicates from each line [17].

3. Results

3.1. Urethane and DEN Treatments Induce Liver Tumors in AIRmax and AIRmin Mice

Mice were injected with urethane or DEN, and the tumor burden in internal organs was determined after 34 weeks. A reliable parameter to quantify tumor development consists in the multiplicity of lesions with a mm. Figure 1 shows the mean number of ≥2 mm liver tumors after DEN or urethane injection in AIRmax and AIRmin mice. Irrespective of the carcinogen used, AIRmax mice are significantly more susceptible than AIRmin mice to developing liver tumors. Spontaneous development of liver tumors, up to the age of 1.5 years, was not observed in these mice.

Histopathological analysis revealed the occurrence of hepatocarcinomas with both treatments (Figure 2).

(a)

(b)

(c)

(d)

Figure 2

Macroscopic aspect of the AIRmax male mouse liver at 34 weeks after DEN treatment (a). HE-stained histological section showing high-grade hepatocarcinoma with atypical hepatocytes arranged in macrotrabecules forming pseudoacini (b, 200x). AIRmax liver after urethane treatment (c); HE-stained section showing high-grade hepatocarcinoma, presence of steatosis, cellular pleomorphism, alteration of lobular architecture of the hepatic parenchyma, and intratumoral inflammatory infiltrate. Arrow indicates an aberrant mitosis (d, 1000x).

3.2. Inflammatory Cytokine Levels Increase in the Liver after DEN Treatment

The levels of IL-6, TNF-α, and IL-1β in liver macerate supernatants were measured in order to evaluate the early DEN-induced liver inflammation. As shown in Figure 3, DEN induced cytokine production as soon as 2 h after injection. Coincidently, the cytokine levels were higher in the livers of susceptible AIRmax mice, suggesting a role for carcinogen-induced inflammatory signaling in HCC development.

(a)

(b)

(c)

3.3. Linkage Analysis of Hepatocarcinogenesis Modifier Loci

We aimed to map the genetic factors that contribute to the divergent susceptibility to liver carcinogenesis of AIRmax and AIRmin mice. Then, we carried out linkage analysis in an intercross F2 (AIRmax×AIRmin) population treated with urethane. Details of urethane-induced liver tumors in the groups of male mice representing the parental and F1 hybrids used for the production of the F2 intercross population are summarized in Table 1. Females were much more resistant, with only 2% incidence in the F2 population.

The logarithm of odds (LOD) thresholds for linkage of the phenotype N_2 (number of liver tumors with a mm), determined by permutation analysis ( permutations and sex interaction), were 4.5 and 5.4 for and 0.01 genome-wide statistical probabilities, respectively. Because N_2 did not follow a normal distribution, the N_2 values were square root-transformed for interval mapping.

Two loci with significant LOD scores were detected in chromosomes 2 and 9. The LOD score in chromosome 2 was localized at SNP CEL_2 63143553 (70 Mb), and for the QTL in chromosome 9, the LOD score was localized at rs13480103 (17 Mb). An additional region with suggestive linkage (; LOD ) was detected at 143 Mb in chromosome 5 (Figure 4).

F2 mice carrying the AIRmax-derived genotype GG (78/334) at CEL_2 63143553 in chromosome 2 had a 2-fold higher N_2 than mice carrying AIRmin-derived genotype AA (86/334; ). Heterozygous mice had intermediate N_2 values. In chromosome 9, inheritance of the AA genotype (81/334) at rs13480103 from AIRmax was associated with 2.8-fold increased N_2 when compared to the group carrying the AIRmin-derived GG genotype (87/334; ). In mice carrying the combined genotypes for both SNPs (21 AIRmin-like, 87 heterozygous, and 25 AIRmax-like genotypes), the difference in N_2 values increased to 4.4-fold (), evidencing an additive effect of the QTLs (Figure 5).

We carried out a similar genetic linkage study in an advanced (A/J×C57BL/6) F4 intercross population (ABF4) with 183 urethane-treated male mice. In this population, 548 SNPs were informative (polymorphic) and nonredundant. Simple interval mapping (SIM) was performed to detect QTLs for the square root-transformed total liver tumor volume calculated by the sum of every visible liver tumor volumes. LOD thresholds for 10,000 permutations were 3.83 and 3.48 for significance at and , respectively. SIM of urethane-induced liver tumors mapped one significant locus at 44 Mb in chromosome 7 (LOD ) and one suggestive locus at 114 Mb in chromosome 1 (LOD , Figure 6). Similar results were obtained for liver tumor multiplicity.

3.4. Transcriptome of Normal Liver Differs in AIRmax and AIRmin Mice

The gene expression profile of normal livers from adult untreated AIRmax and AIRmin mice, analyzed by unsupervised hierarchical clustering (considering all the 28,120 expressed probe sets), was distinct in the two mouse lines (Figure 6). We identified 145 transcripts (corresponding to 85 unique genes) discriminating the two mouse lines. Nineteen genes were upregulated and 66 were downregulated in AIRmax samples when compared to AIRmin samples (Figure 7).

Microarray results were validated by real-time PCR using the same individual mRNA samples. The fold changes between AIRmax and AIRmin for Vnn, H2ea-ps, Defb1, S1pr5, Cyp1a1, Cyp1b1, Tnfa, and Il6 genes detected by real-time PCR were significantly correlated with the gene expression levels detected by microarrays (Pearson correlation ).

4. Discussion

In the present study, we show that AIRmax male mice are highly susceptible and AIRmin mice are resistant to urethane-induced development of hepatocellular carcinoma whereas females from both strains are resistant. This phenotype was confirmed by the treatment of mice with the tobacco-related nitrosamine N-nitrosodiethylamine (NDEA or DEN), a typical genotoxic chemical carcinogen that forms DNA adducts after bioactivation by P450 isozymes, including cytochrome P450 2E1 in the liver [18, 19].

The DEN-treated animals also developed pulmonary nodules with a higher incidence in AIRmin, in agreement with the results of urethane-induced tumors [9]. Thus, the inverse susceptibility to liver versus lung carcinogenesis of AIRmax and AIRmin mice, following treatment with chemicals that induce tumors in both organs, indicates that tissue-specific genetic and environmental factors are responsible for these phenotypes [20].

In a short-term experiment, treatment with DEN induced TNF-α, IL-6, and IL-1β production in the liver. The proinflammatory cytokine IL-6, considered an acute phase protein, is increased in conditions of liver damage and hepatitis [21]. IL-6 concentrations were higher in the livers from AIRmax mice when compared to AIRmin mice. The difference was also observed for females, with variation during the time course of the experiment. Higher production of IL-1β and of TNF-α was also found in AIRmax liver macerates compared to AIRmin liver macerates, beginning 2 h after DEN injection and reaching peak concentrations at 14 and 24 h. TNF-α and IL-1β are major mediators of inflammation and display both pro- and antitumor activities. IL-1β has a role in tumor progression via induction of several tumor growth factors and of endothelial-mesenchymal transition [22, 23].

Several studies reported that the early DEN-elicited DNA damage contributes to necrotic liver cell death with the liberation of danger-associated molecular patterns (DAMPs). These in turn bind to toll like receptors (TLRs) that signal through the adaptor protein MyD88. MyD88 signaling induces translocation of NF-κB to the nucleus and the transcription of cytokines such as IL-6 and TNF-α that are essential for the initiation of DEN-induced liver cancer [24–29].

It has to be pointed out that IL-6, IL-1β, and TNF-α were not detected in the liver macerates of AIRmax or AIRmin females in the first two hours after DEN injection, whereas they were present in susceptible AIRmax male mice. The liver is a sexually dimorphic organ, and gender disparity seen in HCC progression involves sex hormones. Estrogen prevents necrosis of DEN-triggered hepatocytes in female mice; at concentrations present in females but not in males, estrogen modulates or suppresses early IL-6 cytokine production, a mechanism that could account for the inhibition of liver carcinogenesis observed in females [24]. On the other side, in studies using KO mice, IL-6 antagonists or estrogen analogs to prevent HCC indicate that other mechanisms must be induced in females for liver cancer protection [25–28].

The connection between chronic inflammation and carcinogenesis has long been described for several tumor types. In normal conditions, the initial inflammatory response to injury during the acute phase is followed by an intermediate downregulation, ending in the resolution phase of inflammation. Lack of resolution is followed by a chronic phase that can be protumorigenic. In this situation, local alterations start with early activated inflammatory cells interacting with nonimmune cells (epithelial, vascular, and neuronal) via secretion of acute phase mediators. Therefore, the acute inflammation triggered in the acute phase response to injury, if unresolved, promotes tumorigenesis. In our model, the possible failure in the ability to resolve the intense carcinogen-induced acute inflammatory response in the liver of AIRmax male mice could contribute to their susceptibility to developing HCC.

Next, we studied the genetic basis for the divergent susceptibility of AIRmax and AIRmin mice to chemically induced HCC. In humans, many genetic predisposing factors to HCC are related to known environmental risk, such as obesity, alcohol use, or viral hepatitis [21, 29–31]. Nevertheless, in the absence of these factors, genetic/epidemiological studies detected an association between germline mutations in genes that underlie other diseases and increased HCC risk [32].

In mice, the complex nature of HCC genetic control was demonstrated in animal models, and several susceptible QTLs were described: Hcs 1-8 (for HCC susceptibility locus) in chromosomes 1, 2, 5, 7, 8, 12, and 19 and Hcf 1-2 (for HCC female susceptibility locus) in chromosomes 1 and 17. Other QTLs were associated with resistance to HCC: Hcr2 in chromosome 10, Hpcr3 in central chromosome 15, and Hpcr4 in chromosome 1 [33, 34].

Here, we carried out two independent experiments for mapping urethane-induced HCC. We used 7-day-old F2 (AIRmax×AIRmin) mice and ABF4 mice in the adult phase (4 weeks after birth) aiming to compare the genetic control of urethane-induced HCC in both conditions and in different genetic backgrounds. Like many other carcinogens, urethane requires metabolic activation to become mutagenic. Compared to adults, newborn mice are more liable to the toxic and carcinogenic effects of urethane and other agents, because of lower metabolic capacity, influence of hormones, immunologic responses, and rapidly dividing cells in target tissues [35]. The different treatment protocols indeed pointed to different sets of loci in the two intercrossed populations, even if the liver tumor phenotypes used in linkage assays, scored at 34 or 40 weeks after urethane injection in F2 (AIRmax and AIRmin) and in ABF4 mice, respectively, were similar in both models.

In F2 (AIRmax×AIRmin) mice, we detected significant linkage in chromosomes 2 at 70 Mb (1-LOD score interval: 44-77 Mb) and 9 at 17 Mb (1-LOD score interval: 0-27 Mb). In the ABF4 population, one locus was detected in chromosome 7 at 52 Mb, spanning a region from 40 to 57 Mb. These three regions are novel and nonoverlapping with sets of loci reported in rat and mouse crosses. These results support the concept of genetic heterogeneity, where distinct genetic elements exerting pleiotropic effects control the same complex phenotype. Further support comes from one skin carcinogenesis study that revealed the role of genetic heterogeneity in the control of predisposition to inflammatory response and skin tumor susceptibility [36]. The association of diverse sets of loci, from different families, with the corresponding diseases is evidence that genetic heterogeneity also occurs in humans [37].

Furthermore, at these QTL regions, there are interesting coincidences with QTLs for other phenotypes described in independent experiments. In AIRmax and AIRmin mice, a linkage study including clustering by family identified a suggestive QTL for inflammatory response at a 53-56 Mb interval in chromosome 2 [38] suggesting that the selection process of AIRmax and AIRmin lines might have affected the HCC susceptibility QTL mapped at chromosome 2. The chromosome 9 QTL interval contains several genes involved in inflammation and/or in cancer such as Caspase 1 (Casp1) and pannexin 1 (Panx1), involved in the release of mature inflammatory IL-1β [39]. In ABF4 mice, possible candidate gene mapping in the chromosome 7 QTL is Fancf (Fanconi anemia complement group F). Complications of Fanconi anemia include the development of liver tumors and other cancers [40]; Fcgrt gene (Fc receptor, IgG, alpha chain transporter) codes for hepatic FcRn regulate albumin homeostasis and susceptibility to liver injury enhancing sensitivity to albumin bound hepatotoxins [41]; Fgf21 (fibroblast growth factor 21) is a novel biomarker for nonalcoholic fatty liver disease (NAFLD) in humans and limits ethanol-associated hepatotoxicity in mice [42]; and Saa1 (serum amyloid A1) is an evolutionary highly conserved proinflammatory acute phase protein. Saa is predominantly secreted by hepatocytes and modulates fibrinogenesis by inducing inflammation, proliferation, and cell death in hepatic stellate cells [43].

We observed also coincidences with other QTLs such as Lfibq14 and 21 for liver fibrosis involved in the regulatory network that determines the vulnerability to hepatic chronic injury and Ap7q, for alcohol preference [44, 45].

We determined the liver gene expression profile of AIRmax and AIRmin mice as an additional tool to identify candidate genes, because the transcriptome of normal tissue has been associated with genetic predisposition to tumorigenesis in different organs and experimental models. We observed that several transcripts with potential contribution to HCC susceptibility differ significantly in the normal liver from the two mouse lines. For example, genes grouped into functional categories such as defense response (lymphocyte antigens, major histocompatibility complex genes, integrins, and defensin) and cell cycle (cyclin) and serpina1d are upregulated in AIRmin livers. β-Defensin-1 behaves as a tumor suppressor protein in some human tumor types [46]. The vascular noninflammatory molecules Vanin pantetheinases (Vnn) are potential therapeutic targets in inflammatory diseases; Vanin genes influence the lipid profile and studies report a role for Vanin in inflammation, oxidative stress, cell migration, and numerous diseases [47]. H2Ea: mice bearing the H2b haplotype do not express I-E MHC class II molecules, due to deletion of the promoter region of the Ea gene. We verified before that the majority of AIRmax mice bear the H2b haplotype which might underlie the absence of H2Ea expression in the liver as well as in all tissues tested so far [11, 12]. AIRmax mice are susceptible to autoimmune diseases such as pristane-induced arthritis, and the E-alpha chain was described as protective to the development of systemic lupus erythematosus [48].

On the other hand, several CYP450 genes, implicated in xenobiotic metabolism, and genes involved in lipid metabolism are upregulated in AIRmax livers. The sphyngosine-1-receptor 5 (Sp1r5) gene maps at 21 Mb in chromosome 9, inside the QTL confidence interval. The bioactive lipid mediator sphingosine-1-phosphate (S1P), through activation of cell surface receptors S1PR1 to S1PR5, is involved in multiple cellular signaling systems and has a pivotal role in immune cell functions, inflammation, and cancer due to the control of immune cell trafficking. NK cells abundantly express S1PR5 and S1PR5 responses play a role in the trafficking of these cells from the bone marrow and lymph nodes into the tissue [49, 50]. It was previously demonstrated that AIRmax mice had higher relative numbers of NK cells in the spleen, as well as higher cytolytic activity against Yac.1 target cells than AIRmin mice [51]. NK cells are involved in the clearance of infectious agents and antitumor surveillance [52]; however, these cells can accelerate liver injury by producing proinflammatory cytokines and killing hepatocytes [53]. The S1P axis has been implicated in cancer and inflammatory diseases, and S1pr5 is a good candidate modulator of the differential inflammation and HCC susceptibility of AIRmax and AIRmin mice [54].

In conclusion, our results show that the divergent susceptibility of AIRmax and AIRmin mice to chemically induced liver tumors correlates with specific hepatic gene expression profiles in untreated mice and the early production of inflammatory mediators in the liver following carcinogen injection. Linkage analysis in F2 (AIRmax×AIRmin) mice, as well as in ABF4 populations, revealed QTLs for HCC, in chromosomes 2 and 9, and 7, respectively. Genes and QTLs for phenotypes relevant for liver cancer, such as inflammatory response, liver fibrosis, alcoholism, obesity and lipid metabolism, and Fanconi anemia, also map to these regions. The experiments evidenced genetic heterogeneity in the control of predisposition to liver cancer and the role of inflammatory mediators secreted by resident or infiltrating cells in tumor development.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

The 2 first authors (Carvalho L. and Borrego A.) contributed equally to the study.

Acknowledgments

LRC had a fellowship from FAPESP. OCMI and MDF had grants from CNPq.

References

J. Ferlay, I. Soerjomataram, R. Dikshit et al., “Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012,” International Journal of Cancer, vol. 136, no. 5, pp. E359–E386, 2015.
View at: Publisher Site | Google Scholar
A. P. Venook, C. Papandreou, J. Furuse, and L. L. de Guevara, “The incidence and epidemiology of hepatocellular carcinoma: a global and regional perspective,” The Oncologist, vol. 15, Supplement 4, pp. 5–13, 2010.
View at: Publisher Site | Google Scholar
G. Montalto, M. Cervello, L. Giannitrapani, F. Dantona, A. Terranova, and L. A. Castagnetta, “Epidemiology, risk factors, and natural history of hepatocellular carcinoma,” Annals of the New York Academy of Sciences, vol. 963, pp. 13–20, 2002.
View at: Google Scholar
S. H. Yeh and P. J. Chen, “Gender disparity of hepatocellular carcinoma: the roles of sex hormones,” Oncology, vol. 78, no. 1, pp. 172–179, 2010.
View at: Publisher Site | Google Scholar
O. M. Ibanez, C. Stiffel, O. G. Ribeiro et al., “Genetics of nonspecific immunity: I. bidirectional selective breeding of lines of mice endowed with maximal or minimal inflammatory responsiveness,” European Journal of Immunology, vol. 22, no. 10, pp. 2555–2563, 1992.
View at: Publisher Site | Google Scholar
O. G. Ribeiro, D. A. Maria, S. Adriouch et al., “Convergent alteration of granulopoiesis, chemotactic activity, and neutrophil apoptosis during mouse selection for high acute inflammatory response,” Journal of Leukocyte Biology, vol. 74, no. 4, pp. 497–506, 2003.
View at: Publisher Site | Google Scholar
R. F. Di Pace, S. Massa, O. G. Ribeiro et al., “Inverse genetic predisposition to colon versus lung carcinogenesis in mouse lines selected based on acute inflammatory responsiveness,” Carcinogenesis, vol. 27, no. 8, pp. 1517–1525, 2006.
View at: Publisher Site | Google Scholar
G. Biozzi, O. G. Ribeiro, A. Saran et al., “Effect of genetic modification of acute inflammatory responsiveness on tumorigenesis in the mouse,” Carcinogenesis, vol. 19, no. 2, pp. 337–346, 1998.
View at: Publisher Site | Google Scholar
D. A. Maria, G. Manenti, F. Galbiati et al., “Pulmonary adenoma susceptibility 1 (Pas1) locus affects inflammatory response,” Oncogene, vol. 22, no. 3, pp. 426–432, 2003.
View at: Publisher Site | Google Scholar
J. R. Jensen, A. Galvan, A. Borrego et al., “Genetic control of renal tumorigenesis by the mouse Rtm1 locus,” BMC Genomics, vol. 14, no. 1, p. 724, 2013.
View at: Publisher Site | Google Scholar
P. S. dos S Carneiro, L. C. Peters, F. Vorraro et al., “Gene expression profiles of bone marrow cells from mice phenotype-selected for maximal or minimal acute inflammations: searching for genes in acute inflammation modifier loci,” Immunology, vol. 128, 1, part 2, pp. e562–e571, 2009.
View at: Publisher Site | Google Scholar
M. D. Franco, F. Colombo, A. Galvan et al., “Transcriptome of normal lung distinguishes mouse lines with different susceptibility to inflammation and to lung tumorigenesis,” Cancer Letters, vol. 294, no. 2, pp. 187–194, 2010.
View at: Publisher Site | Google Scholar
A. Galvan, W. Cabrera, F. Vorraro et al., “Genetic linkage analysis identifies Pas1 as the common locus modulating lung tumorigenesis and acute inflammatory response in mice,” Genes and Immunity, vol. 14, no. 8, pp. 512–517, 2013.
View at: Publisher Site | Google Scholar
A. Dassano, F. Colombo, G. Trincucci et al., “Mouse pulmonary adenoma susceptibility 1 locus is an expression QTL modulating Kras-4A,” PLoS Genetics, vol. 10, no. 4, article e1004307, 2014.
View at: Publisher Site | Google Scholar
F. Vorraro, A. Galvan, W. H. K. Cabrera et al., “Genetic control of IL-1β production and inflammatory response by the mouse Irm1 locus,” The Journal of Immunology, vol. 185, no. 3, pp. 1616–1621, 2010.
View at: Publisher Site | Google Scholar
G. Seaton, J. Hernandez, J. A. Grunchec et al., “GridQTL: a grid portal for QTL mapping of compute intensive datasets,” in Proceedings of the 8th world congress on genetics applied to livestock production, p. 27, Belo Horizonte Brazil, 2006.
View at: Google Scholar
V. G. Tusher, R. Tibshirani, and G. Chu, “Significance analysis of microarrays applied to the ionizing radiation response,” Proceedings of the National Academy of Sciences of the United States of America, vol. 98, no. 9, pp. 5116–5121, 2001.
View at: Publisher Site | Google Scholar
L. Verna, J. Whysner, and G. M. Williams, “N-Nitrosodiethylamine mechanistic data and risk assessment: bioactivation, DNA-adduct formation, mutagenicity, and tumor initiation,” Pharmacology Therapeutics, vol. 71, no. 1-2, pp. 57–81, 1996.
View at: Publisher Site | Google Scholar
J. S. Kang, H. Wanibuchi, K. Morimura, F. J. Gonzalez, and S. Fukushima, “Role of CYP2E1 in diethylnitrosamine-induced hepatocarcinogenesis in vivo,” Cancer Research, vol. 67, no. 23, pp. 11141–11146, 2007.
View at: Publisher Site | Google Scholar
H. H. H. Göring, “Tissue specificity of genetic regulation of gene expression,” Nature Genetics, vol. 44, no. 10, pp. 1077-1078, 2012.
View at: Publisher Site | Google Scholar
E. J. Park, J. H. Lee, G. Y. Yu et al., “Dietary and genetic obesity promote liver inflammation and tumorigenesis by enhancing IL-6 and TNF expression,” Cell, vol. 140, no. 2, pp. 197–208, 2010.
View at: Publisher Site | Google Scholar
A. M. Lewis, S. Varghese, H. Xu, and H. R. Alexander, “Interleukin-1 and cancer progression: the emerging role of interleukin-1 receptor antagonist as a novel therapeutic agent in cancer treatment,” Journal of Translational Medicine, vol. 4, no. 1, p. 48, 2006.
View at: Publisher Site | Google Scholar
C. A. Dinarello, “Immunological and inflammatory functions of the interleukin-1 family,” Annual Review of Immunology, vol. 27, no. 1, pp. 519–550, 2009.
View at: Publisher Site | Google Scholar
W. E. Naugler, T. Sakurai, S. Kim et al., “Gender disparity in liver cancer due to sex differences in MyD88-dependent IL-6 production,” Science, vol. 317, no. 5834, pp. 121–124, 2007.
View at: Publisher Site | Google Scholar
Z. Li, G. Tuteja, J. Schug, and K. H. Kaestner, “Foxa1 and Foxa2 are essential for sexual dimorphism in liver cancer,” Cell, vol. 148, no. 1-2, pp. 72–83, 2012.
View at: Publisher Site | Google Scholar
M. Di Maio, E. De Maio, A. Morabito et al., “Hormonal treatment of human hepatocellular carcinoma,” Annals of the New York Academy of Sciences, vol. 1089, no. 1, pp. 252–261, 2006.
View at: Publisher Site | Google Scholar
M. Kalra, J. Mayes, S. Assefa, A. K. Kaul, and R. Kaul, “Role of sex steroid receptors in pathobiology of hepatocellular carcinoma,” World Journal of Gastroenterology, vol. 14, no. 39, pp. 5945–5961, 2008, Review.
View at: Publisher Site | Google Scholar
T. Lawrence, T. Hageman, and F. Balkwill, “Cancer. Sex, cytokines, and cancer,” Science, vol. 317, no. 5834, pp. 51-52, 2007.
View at: Publisher Site | Google Scholar
S. Yoshimoto, T. M. Loo, K. Atarashi et al., “Obesity-induced gut microbial metabolite promotes liver cancer through senescence secretome,” Nature, vol. 499, no. 7456, pp. 97–101, 2013.
View at: Publisher Site | Google Scholar
J. A. Marrero, R. J. Fontana, S. Fu, H. S. Conjeevaram, G. L. Su, and A. S. Lok, “Alcohol, tobacco and obesity are synergistic risk factors for hepatocellular carcinoma,” Journal of Hepatology, vol. 42, no. 2, pp. 218–224, 2005.
View at: Publisher Site | Google Scholar
A. J. Walker, C. J. Peacock, V. Pedergnana, STOP-HCV Consortium, and W. L. Irving, “Host genetic factors associated with hepatocellular carcinoma in patients with hepatitis C virus infection: a systematic review,” Journal of Viral Hepatitis, vol. 25, no. 5, pp. 442–456, 2018.
View at: Publisher Site | Google Scholar
D. Dong, Y. Zou, P. Zhang, and Z. Wu, “Systematic analyses and comprehensive field synopsis of genetic association studies in hepatocellular carcinoma,” Oncotarget, vol. 7, no. 29, pp. 45757–45763, 2016, Review.
View at: Publisher Site | Google Scholar
P. Demant, “Cancer susceptibility in the mouse: genetics, biology and implications for human cancer,” Nature Reviews Genetics, vol. 4, no. 9, pp. 721–734, 2003.
View at: Publisher Site | Google Scholar
G. Manenti, A. Galvan, F. S. Falvella et al., “Genetic control of resistance to hepatocarcinogenesis by the mouse Hpcr3 locus,” Hepatology, vol. 48, no. 2, pp. 617–623, 2008.
View at: Publisher Site | Google Scholar
S. D. Vesselinovitch, M. Koka, N. Mihailovich, and K. V. N. Rao, “Carcinogenicity of diethylnitrosamine in newborn, infant, and adult mice,” Journal of Cancer Research and Clinical Oncology, vol. 108, no. 1, pp. 60–65, 1984.
View at: Publisher Site | Google Scholar
A. Galvan, F. Vorraro, W. H. K. Cabrera et al., “Genetic heterogeneity of inflammatory response and skin tumorigenesis in phenotypically selected mouse lines,” Cancer Letters, vol. 295, no. 1, pp. 54–58, 2010.
View at: Publisher Site | Google Scholar
T. A. Dragani, “Risk of HCC: genetic heterogeneity and complex genetics,” Journal of Hepatology, vol. 52, no. 2, pp. 252–257, 2010.
View at: Publisher Site | Google Scholar
A. Galvan, F. Vorraro, W. Cabrera et al., “Association study by genetic clustering detects multiple inflammatory response loci in non-inbred mice,” Genes and Immunity, vol. 12, no. 5, pp. 390–394, 2011.
View at: Publisher Site | Google Scholar
P. Pelegrin and A. Surprenant, “Pannexin-1 mediates large pore formation and interleukin-1beta release by the ATP-gated P2X7 receptor,” The EMBO Journal, vol. 25, no. 21, pp. 5071–5082, 2006.
View at: Publisher Site | Google Scholar
B. P. Alter, “Cancer in Fanconi anemia, 1927–2001,” Cancer, vol. 97, no. 2, pp. 425–440, 2003.
View at: Publisher Site | Google Scholar
M. Pyzik, T. Rath, T. T. Kuo et al., “Hepatic FcRn regulates albumin homeostasis and susceptibility to liver injury,” Proceedings of the National Academy of Sciences of the United States of America, vol. 114, no. 14, pp. E2862–E2871, 2017.
View at: Publisher Site | Google Scholar
B. N. Desai, G. Singhal, M. Watanabe et al., “Fibroblast growth factor 21 (FGF21) is robustly induced by ethanol and has a protective role in ethanol associated liver injury,” Molecular Metabolism, vol. 6, no. 11, pp. 1395–1406, 2017.
View at: Publisher Site | Google Scholar
S. V. Siegmund, M. Schlosser, F. A. Schildberg et al., “Serum amyloid A induces inflammation, proliferation and cell death in activated hepatic stellate cells,” PLoS One, vol. 11, no. 3, article e0150893, 2016.
View at: Publisher Site | Google Scholar
R. A. Hall, R. Liebe, K. Hochrath et al., “Systems genetics of liver fibrosis: identification of fibrogenic and expression quantitative trait loci in the BXD murine reference population,” PLoS One, vol. 9, no. 2, article e89279, 2014.
View at: Publisher Site | Google Scholar
A. A. Bachmanov, D. R. Reed, X. Li, S. Li, G. K. Beauchamp, and M. G. Tordoff, “Voluntary ethanol consumption by mice: genome-wide analysis of quantitative trait loci and their interactions in a C57BL/6ByJ x 129P3/J F2 intercross,” Genome Research, vol. 12, no. 8, pp. 1257–1268, 2002.
View at: Publisher Site | Google Scholar
N. Droin, J. B. Hendra, P. Ducoroy, and E. Solary, “Human defensins as cancer biomarkers and antitumour molecules,” Journal of Proteomics, vol. 72, no. 6, pp. 918–927, 2009.
View at: Publisher Site | Google Scholar
B. J. Kaskow, J. M. Proffit, J. Blangero, E. K. Moses, and L. J. Abraham, “Diverse biological activities of the vascular non-inflammatory molecules - the Vanin pantetheinases,” Biochemical and Biophysical Research Communications, vol. 417, no. 2, pp. 653–658, 2012.
View at: Publisher Site | Google Scholar
E. Martínez-Soria, M. L. Santiago-Raber, L. Ho, T. Moll, and S. Izui, “Protection of murine systemic lupus by the Ea transgene without expression of I-E heterodimers,” The Journal of Immunology, vol. 181, no. 5, pp. 3651–3657, 2008.
View at: Publisher Site | Google Scholar
K. Mayol, V. Biajoux, J. Marvel, K. Balabanian, and T. Walzer, “Sequential desensitization of CXCR4 and S1P5 controls natural killer cell trafficking,” Blood, vol. 118, no. 18, pp. 4863–4871, 2011.
View at: Publisher Site | Google Scholar
C. N. Jenne, A. Enders, R. Rivera et al., “T-bet–dependent S1P5 expression in NK cells promotes egress from lymph nodes and bone marrow,” Journal of Experimental Medicine, vol. 206, no. 11, pp. 2469–2481, 2009.
View at: Publisher Site | Google Scholar
L. Castoldi, M. A. Golim, O. G. R. Filho, G. G. Romagnoli, O. C. M. Ibañez, and R. Kaneno, “Enhanced natural killer activity and production of pro-inflammatory cytokines in mice selected for high acute inflammatory response (AIRmax),” Immunology, vol. 120, no. 3, pp. 372–379, 2007.
View at: Publisher Site | Google Scholar
P. Liu, L. Chen, and H. Zhang, “Natural killer cells in liver disease and hepatocellular carcinoma and the NK cell-based immunotherapy,” Journal of Immunology Research, vol. 2018, Article ID 1206737, 8 pages, 2018.
View at: Publisher Site | Google Scholar
B. Gao, S. Radaeva, and O. Park, “Liver natural killer and natural killer T cells: immunobiology and emerging roles in liver diseases,” Journal of Leukocyte Biology, vol. 86, no. 3, pp. 513–528, 2009.
View at: Publisher Site | Google Scholar
G. T. Kunkel, M. Maceyka, S. Milstien, and S. Spiegel, “Targeting the sphingosine-1-phosphate axis in cancer, inflammation and beyond,” Nature Reviews Drug Discovery, vol. 12, no. 9, pp. 688–702, 2013.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Lilian Rego de Carvalho et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1113

Downloads

959

Citations