Association of CMYC polymorphisms with hepatoblastoma risk
Introduction
Hepatoblastoma is the most common hepatic tumor of childhood (1,2). The incidence of hepatoblastoma is about 0.5–1.5 cases per million, and the mortality rate can be as high as 35–50% for high-risk patients (3). Over the past decades, efforts have been made to improve the outcome of hepatoblastoma. However, treatment has not changed significantly in the past 20 years (4). In recent years, several unique genetic features have been identified to be associated with hepatoblastoma, providing new insights into the understanding of hepatoblastoma (5). The elucidation of the genetic features of hepatoblastoma is thus of critical importance.
Single-nucleotide polymorphisms (SNPs) are the most common sources of genetic variation in the genome and are frequently associated with potential cancer risk (6). Some SNPs contributing to the progression of hepatoblastoma have been identified. Arai et al. revealed that MDM4 polymorphisms are significantly correlated with the outcomes of hepatoblastoma (7). Based on high-density SNP genotyping microarrays, Suzuki et al. demonstrated that expression levels of IGF2 and H19 were significantly correlated with hepatoblastoma (8). c-Myc is a well-known human transcription factor involved in cell cycle, growth, metabolism, and apoptosis (9). A previous study showed that the CMYC rs6883267 polymorphism is significantly associated with CMYC transcription efficiency and poor prognosis in colorectal cancer (10). However, the association between CMYC polymorphisms and hepatoblastoma remains unclear. This study therefore aimed to investigate the association of CMYC polymorphisms with hepatoblastoma susceptibility.
Methods
Patients
Patients less than 18 years old with a pathologic diagnosis of hepatoblastoma were enrolled. Cancer-free control subjects matched for age and sex were recruited from the same area. All patients and control subjects were genetically unrelated members of the Chinese Han population. Written informed consent was acquired from all participants’ legal guardians or parents. The institutional review board of Guangzhou Women and Children’s Medical Center approved this study. All patient data were anonymous or de-identified prior to analysis.
CMYC genotyping
Allelic discrimination of the rs4645943 and rs2070583 polymorphisms of CMYC was performed using TaqMan reagents (Applied Biosystems, Foster City, CA, USA) according to the manufacturer’s protocol, as reported previously (11-14). Control samples of known genotype were also included in each test, including blank, homozygous wild-type, homozygous mutant, and heterozygous samples. Quality control was performed with eight negative and positive control samples on each of the 384-well plates; 10% of the samples were also randomly selected for a second round of genotyping, and the concordance rate was 100%.
Statistical analysis
All statistical analyses were performed with SAS software (version 9.1; SAS Institute, Cary, NC, USA). Continuous variables were analyzed using Student’s t-test or one-way analysis of variance. Categorical variables were analyzed by χ2 test. Differences in allele or genotype frequencies between patients and controls were determined by χ2 test. Hardy-Weinberg equilibrium (HWE) was calculated using a goodness-of-fit χ2 test for biallelic markers. Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated for evaluation of the strength of the association of interest (15-17). Adjusted ORs were calculated using multivariate analysis after adjusting for age, sex, and clinical stage. Differences were considered significant at P<0.05.
Results
Characteristics of participants enrolled in this study
A total of 213 hepatoblastoma patients and 958 control subjects were recruited from Guangdong, Henan, Shaanxi, and Shanxi provinces in China. Males made up the majority of both the hepatoblastoma and control groups, accounting for 60.56% and 60.44% of individuals, respectively. Most of the patients had stage II disease (n=55), followed by stage I (n=42), stage III (n=40), and stage IV (n=15); stage information was lacking for 61 patients (Table S1). There were no significant differences between cases and controls regarding the distributions of age and sex (P>0.05, Table 1).
Table 1
Variables | Cases (n=213), N (%) | Controls (n=958), N (%) | P† |
---|---|---|---|
Age range, months | 0.23–149.97 | 0.004–156.000 | 0.105 |
Mean ± SD | 23.62±24.36 | 23.75±18.30 | |
<17 | 114 (53.52) | 454 (47.39) | |
≥17 | 99 (46.48) | 504 (52.61) | |
Sex | 0.973 | ||
Female | 84 (39.44) | 379 (39.56) | |
Male | 129 (60.56) | 579 (60.44) | |
Clinical stages | – | ||
I | 42 (19.72) | – | |
II | 55 (25.82) | – | |
III | 40 (18.78) | – | |
IV | 15 (7.04) | – | |
NA‡ | 61 (28.64) | – |
†, Two-sided χ2 test for distributions between hepatoblastoma patients and cancer-free controls; ‡, stage information was absent. SD, standard deviation; NA, not applicable.
Association between CMYC polymorphisms and hepatoblastoma risk
Genotype distributions and associations between CMYC gene polymorphisms and hepatoblastoma risk are summarized in Table 2. For rs4645943, compared with carriers of the CC genotype, carriers of the CT (OR, 1.10; 95% CI, 0.81–1.51; P=0.532) or TT (OR, 1.10; 95% CI, 0.63–1.92; P=0.726) genotypes showed no significant associations with hepatoblastoma risk. Moreover, there was no significant association between rs4645943 and hepatoblastoma risk under the additive (OR, 1.07; 95% CI, 0.85–1.35; P=0.550), dominant (OR, 1.10; 95% CI, 0.82–1.49; P=0.512), or recessive models (OR, 1.06; 95% CI, 0.62–1.81; P=0.842). For rs2070583, compared with carriers of the AA genotype, carriers of the AG (OR, 1.12; 95% CI, 0.80–1.55; P=0.516) and GG (OR, 0.84; 95% CI, 0.35–2.04; P=0.699) genotypes exhibited no significant associations with hepatoblastoma risk. Similarly, there was no significant association between rs2070583 and hepatoblastoma risk under the additive (OR, 1.04; 95% CI, 0.79–1.36; P=0.783), dominant (OR, 1.08; 95% CI, 0.79–1.49; P=0.618), or recessive models (OR, 0.81; 95% CI, 0.34–1.97; P=0.645).
Table 2
Genotype | Cases (N=213) | Controls (N=958) | P† | Crude OR (95% CI) | P | AOR (95% CI)‡ | P‡ |
---|---|---|---|---|---|---|---|
rs4645943 C>T (HWE, 0.850) | |||||||
CC | 105 (49.30) | 496 (51.77) | 1.00 | 1.00 | |||
CT | 90 (42.25) | 385 (40.19) | – | 1.10 (0.81–1.51) | 0.533 | 1.10 (0.81–1.51) | 0.532 |
TT | 18 (8.45) | 77 (8.04) | – | 1.10 (0.63–1.92) | 0.726 | 1.10 (0.63–1.92) | 0.726 |
Additive | – | – | 0.807 | 1.07 (0.85–1.35) | 0.550 | 1.07 (0.85–1.35) | 0.550 |
Dominant | 108 (50.70) | 462 (48.23) | 0.513 | 1.10 (0.82–1.49) | 0.513 | 1.10 (0.82–1.49) | 0.512 |
Recessive | 195 (91.55) | 881 (91.96) | 0.842 | 1.06 (0.62–1.81) | 0.842 | 1.06 (0.62–1.81) | 0.842 |
rs2070583 A>G (HWE, 0.319) | |||||||
AA | 143 (67.14) | 660 (68.89) | 1.00 | 1.00 | |||
AG | 64 (30.05) | 265 (27.66) | – | 1.12 (0.80–1.55) | 0.516 | 1.12 (0.80–1.55) | 0.516 |
GG | 6 (2.82) | 33 (3.44) | – | 0.84 (0.35–2.04) | 0.699 | 0.84 (0.35–2.04) | 0.699 |
Additive | – | – | 0.727 | 1.04 (0.79–1.36) | 0.783 | 1.04 (0.79–1.36) | 0.783 |
Dominant | 70 (32.86) | 298 (31.11) | 0.617 | 1.08 (0.79–1.49) | 0.617 | 1.08 (0.79–1.49) | 0.618 |
Recessive | 207 (97.18) | 925 (96.56) | 0.644 | 0.81 (0.34–1.96) | 0.645 | 0.81 (0.34–1.97) | 0.645 |
Combined effect of risk genotypes§ | |||||||
0–1 | 111 (52.11) | 529 (55.22) | 1.00 | 1.00 | |||
2 | 102 (47.89) | 429 (44.78) | 0.410 | 1.13 (0.84–1.53) | 0.410 | 1.13 (0.84–1.53) | 0.410 |
†, χ2 test for genotype distributions between hepatoblastoma patients and cancer-free controls; ‡, adjusted for age and sex; §, risk genotypes were carriers with rs4645943 CT/TT and rs2070583 AA/AG genotypes. AOR, adjusted odds ratio; CI, confidence interval; HWE, Hardy-Weinberg equilibrium.
In addition, we found no significant association between hepatoblastoma risk and the combination of the rs4645943 CT/TT genotype with the rs2070583 AA/AG genotype (OR, 1.13; 95% CI, 0.84–1.53; P=0.410).
Stratification analysis of CMYC genotypes and hepatoblastoma risk
Further analysis showed that neither CMYC polymorphism was significantly associated with hepatoblastoma risk in any of the subgroups of hepatoblastoma patients (Table 3), which were stratified according to age, sex, and clinical tumor stage (all P>0.05). In addition, the combination of the rs4645943 CT/TT and rs2070583 AA/AG genotypes was not significantly associated with hepatoblastoma risk in any subgroups stratified by age, sex, or clinical tumor stage (all P>0.05). These findings suggest that CMYC polymorphisms are not significantly associated with hepatoblastoma susceptibility.
Table 3
Variables | rs4645943 (case/control) | AOR (95% CI)† | P† | rs2070583 (case/control) | AOR (95% CI)† | P† | Combine genotypes (case/control) | AOR (95% CI)† | P† | |||
---|---|---|---|---|---|---|---|---|---|---|---|---|
CC | CT/TT | AA | AG/GG | 0–1 | 2 | |||||||
Age, months | ||||||||||||
<17 | 53/235 | 61/219 | 1.24 (0.82–1.86) | 0.315 | 74/308 | 40/146 | 1.14 (0.74–1.75) | 0.565 | 56/250 | 58/204 | 1.27 (0.84–1.91) | 0.257 |
≥17 | 52/261 | 47/243 | 0.97 (0.63–1.49) | 0.892 | 69/352 | 30/152 | 1.01 (0.63–1.61) | 0.979 | 55/279 | 44/225 | 0.99 (0.64–1.53) | 0.971 |
Sex | ||||||||||||
Female | 38/196 | 46/183 | 1.31 (0.81–2.10) | 0.273 | 54/266 | 30/113 | 1.32 (0.80–2.17) | 0.280 | 42/211 | 42/168 | 1.27 (0.79–2.03) | 0.332 |
Male | 67/300 | 62/279 | 1.00 (0.68–1.46) | 0.986 | 89/394 | 40/185 | 0.96 (0.64–1.45) | 0.846 | 69/318 | 60/261 | 1.06 (0.72–1.56) | 0.763 |
Clinical stage | ||||||||||||
I + II | 43/496 | 54/462 | 1.35 (0.88–2.05) | 0.166 | 61/660 | 36/298 | 1.31 (0.85–2.02) | 0.225 | 48/529 | 49/429 | 1.26 (0.83–1.91) | 0.285 |
III + IV | 27/496 | 28/462 | 1.11 (0.65–1.91) | 0.703 | 40/660 | 15/298 | 0.83 (0.45–1.53) | 0.556 | 28/529 | 27/429 | 1.19 (0.69–2.05) | 0.536 |
†, Adjusted for age and sex, omitting the corresponding stratification factor. AOR, adjusted odds ratio; CI, confidence interval.
Discussion
Our results showed that the CMYC rs4645943 and rs2070583 polymorphisms were not associated with hepatoblastoma susceptibility. Further stratification analysis based on age, sex, and clinical stage found similar results.
CMYC, encoding the c-Myc protein, is an important oncogene involved in many steps of tumorigenesis, such as proliferation, survival, apoptosis, migration, and invasion (18). A previous study revealed that the expression of c-Myc and cyclin-D1 was significantly elevated in pretreated hepatoblastoma samples but decreased after chemotherapy (19). Myc-expressing mice can present with hepatocellular carcinoma and hepatoblastoma-like tumors, but tumor regression can be induced by inhibiting the expression of Myc (20). Hartwell et al. demonstrated that prolactin suppresses hepatocellular carcinoma by inhibiting the innate immune activation of c-Myc in a mouse model (21). Han et al. found that miR-148a-5p and miR-363-3p negatively regulate the expression of c-Myc to modulate hepatocarcinogenesis (22). These findings suggest that the abnormal expression of CMYC may play a critical role in the development of liver cancer.
SNPs may be associated with gene transcriptional activity (20,23). For example, CMYC polymorphisms are cis-regulated in the immortalized lymphocytes of HapMap individuals (23). Lee et al. revealed that the CMYC rs4645943 polymorphism was associated with the warfarin dose requirement in patients undergoing cardiac valve replacement (24). Moreover, the CMYC rs2070583 polymorphism is significantly associated with coronary heart disease in African Americans (25). However, in the current study, no significant associations were found between the CMYC rs4645943 and rs2070583 polymorphisms and hepatoblastoma susceptibility in a Han Chinese population. Therefore, we speculate that abnormal expression of CMYC in hepatoblastoma may not be attributed to CMYC gene polymorphisms. Wang et al. demonstrated that the role of CMYC in hepatoblastoma is to impose mutually dependent alterations in gene expression and metabolic re-programming that are not obtained in non-transformed cells and that cooperate to promote tumor growth (26). The activation of β-catenin is one of the hallmarks of hepatoblastoma, inducing its translocation to the nucleus and activating target genes, including CMYC, MMP genes, and VEGF to regulate cell proliferation, invasion, and angiogenesis (27). As a member of the Wnt signaling pathway, Wnt ligand binding suppresses the phosphorylation of β-catenin to inhibit its downstream target genes, such as CMYC, repressing cell proliferation (28). Taken together, this evidence suggests that the abnormal expression of CMYC (or the protein c-Myc) in hepatoblastoma may largely depend on the regulation of upstream effectors, rather than its genetically encoded information.
Some limitations of this study should be mentioned. First, although we tried to recruit a large number of hepatoblastoma patients, the sample size in this study was still relatively small, and more patients are required to further validate our findings. Second, due to a lack of detailed information on the patients, associations between CMYC polymorphisms and clinical characteristics, such as tumor size and lymph node metastasis, were not analyzed in this study. Lastly, the study population does not represent the complete Chinese population.
Conclusions
In summary, the CMYC rs4645943 and rs2070583 polymorphisms may not be associated with hepatoblastoma risk. The abnormal regulation of CMYC in hepatoblastoma may therefore require further investigations and explanation.
Table S1
Variables | Guangdong province | Henan province | Shaanxi province | Shanxi province | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Cases (n=146), n (%) | Controls (n=438), n (%) | P† | Cases (n=42), n (%) | Controls (n=176), n (%) | P† | Cases (n=15), n (%) | Controls (n=186), n (%) | P† | Cases (n=10), n (%) | Controls (n=158), n (%) | P† | ||||
Age range, months | 0.63–149.97 | 0.07–156.00 | 0.214 | 0.83–108.00 | 0.10–108.00 | 0.285 | 3.60–72.00 | 0.03–60.00 | 0.286 | 0.23–72.00 | 0.004–60.00 | 0.785 | |||
Mean ± SD | 23.16±24.59 | 23.11±18.62 | 26.73±24.96 | 27.28±18.87 | 21.50±24.20 | 23.66±16.66 | 20.32±20.38 | 21.70±18.28 | |||||||
<17 | 79 (54.11) | 211 (48.17) | 21 (50.00) | 72 (40.91) | 9 (60.00) | 85 (45.70) | 5 (50.00) | 86 (54.43) | |||||||
≥17 | 67 (45.89) | 227 (51.83) | 21 (50.00) | 104 (59.09) | 6 (40.00) | 101 (54.30) | 5 (50.00) | 72 (45.57) | |||||||
Gender | 0.961 | 0.830 | 0.544 | 0.912 | |||||||||||
Female | 58 (39.73) | 175 (39.95) | 15 (35.71) | 66 (37.50) | 7 (46.67) | 72 (38.71) | 4 (40.00) | 66 (41.77) | |||||||
Male | 88 (60.27) | 263 (60.05) | 27 (64.29) | 110 (62.50) | 8 (53.33) | 114 (61.29) | 6 (60.00) | 92 (58.23) | |||||||
Clinical stages | |||||||||||||||
I | 6 (4.11) | – | 19 (45.24) | – | 15 (100.00) | – | 2 (20.00) | – | |||||||
II | 46 (31.51) | – | 3 (7.14) | – | – | – | 6 (60.00) | – | |||||||
III | 37 (25.34) | – | 3 (7.14) | – | – | – | 0 (0.00) | – | |||||||
IV | 12 (8.22) | – | 1 (2.38) | – | – | – | 2 (20.00) | – | |||||||
NAb | 45 (30.82) | – | 16 (38.10) | – | – | – | – | – |
†, Two-sided 2 test for distributions between hepatoblastoma patients and cancer-free controls. SD, standard deviation; NA, not applicable.
Acknowledgments
Funding: This research was funded by
Footnote
Conflicts of Interest: The authors have completed the ICMJE uniform disclosure form (available at http://dx.doi.org/10.21037/tcr.2019.12.19). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). Written informed consent was acquired from all participants’ legal guardians or parents. The institutional review board of Guangzhou Women and Children’s Medical Center approved this study (No. 2017120101). All patient data were anonymous or de-identified prior to analysis.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Lim IIP, Bondoc AJ, Geller JI, et al. Hepatoblastoma-the evolution of biology, surgery, and transplantation. Children (Basel) 2018; [Crossref] [PubMed]
- Bennett J, Kirby-Allen M, Ng V, et al. Hepatoblastoma in a child with early-onset cirrhosis. J Pediatr Hematol Oncol 2019;41:e30-3. [Crossref] [PubMed]
- Liu B, Cui X, Zheng S, et al. Aberrant KLK4 gene promoter hypomethylation in pediatric hepatoblastomas. Oncol Lett 2017;13:1360-4. [Crossref] [PubMed]
- Stafman LL, Williams AP, Garner EF, et al. Targeting PIM kinases affects maintenance of CD133 tumor cell population in hepatoblastoma. Transl Oncol 2019;12:200-8. [Crossref] [PubMed]
- Tomlinson GE, Kappler R. Genetics and epigenetics of hepatoblastoma. Pediatr Blood Cancer 2012;59:785-92. [Crossref] [PubMed]
- Qi P, Wang L, Zhou B, et al. Associations of miRNA polymorphisms and expression levels with breast cancer risk in the Chinese population. Genet Mol Res 2015;14:6289-96. [Crossref] [PubMed]
- Arai Y, Honda S, Haruta M, et al. Genome-wide analysis of allelic imbalances reveals 4q deletions as a poor prognostic factor and MDM4 amplification at 1q32.1 in hepatoblastoma. Genes Chromosomes Cancer 2010;49:596-609. [PubMed]
- Suzuki M, Kato M, Yuyan C, et al. Whole-genome profiling of chromosomal aberrations in hepatoblastoma using high-density single-nucleotide polymorphism genotyping microarrays. Cancer Sci 2008;99:564-70. [Crossref] [PubMed]
- Sun L, Song L, Wan Q, et al. cMyc-mediated activation of serine biosynthesis pathway is critical for cancer progression under nutrient deprivation conditions. Cell Res 2015;25:429-44. [Crossref] [PubMed]
- Takatsuno Y, Mimori K, Yamamoto K, et al. The rs6983267 SNP is associated with MYC transcription efficiency, which promotes progression and worsens prognosis of colorectal cancer. Ann Surg Oncol 2013;20:1395-402. [Crossref] [PubMed]
- Chang J, Tian J, Yang Y, et al. A rare missense variant in TCF7L2 associates with colorectal cancer risk by interacting with a GWAS-identified regulatory variant in the MYC enhancer. Cancer Res 2018;78:5164-72. [Crossref] [PubMed]
- Chang J, Tian J, Zhu Y, et al. Exome-wide analysis identifies three low-frequency missense variants associated with pancreatic cancer risk in Chinese populations. Nat Commun 2018;9:3688. [Crossref] [PubMed]
- Chang J, Zhong R, Tian J, et al. Exome-wide analyses identify low-frequency variant in CYP26B1 and additional coding variants associated with esophageal squamous cell carcinoma. Nat Genet 2018;50:338-43. [Crossref] [PubMed]
- Li J, Chang J, Tian J, et al. A rare variant P507L in TPP1 interrupts TPP1-TIN2 interaction, influences telomere length, and confers colorectal cancer risk in Chinese population. Cancer Epidemiol Biomarkers Prev 2018;27:1029-35. [Crossref] [PubMed]
- He J, Wang F, Zhu J, et al. Association of potentially functional variants in the XPG gene with neuroblastoma risk in a Chinese population. J Cell Mol Med 2016;20:1481-90. [Crossref] [PubMed]
- He J, Wang F, Zhu J, et al. The TP53 gene rs1042522 C>G polymorphism and neuroblastoma risk in Chinese children. Aging (Albany NY) 2017;9:852-9. [Crossref] [PubMed]
- He J, Zou Y, Wang T, et al. Genetic variations of GWAS-identified genes and neuroblastoma susceptibility: a replication study in southern Chinese children. Transl Oncol 2017;10:936-41. [Crossref] [PubMed]
- Stine ZE, Walton ZE, Altman BJ, et al. MYC, metabolism, and cancer. Cancer Discov 2015;5:1024-39. [Crossref] [PubMed]
- Armengol C, Cairo S, Fabre M, et al. Wnt signaling and hepatocarcinogenesis: the hepatoblastoma model. Int J Biochem Cell Biol 2011;43:265-70. [Crossref] [PubMed]
- Kress TR, Pellanda P, Pellegrinet L, et al. Identification of MYC-dependent transcriptional programs in oncogene-addicted liver tumors. Cancer Res 2016;76:3463-72. [Crossref] [PubMed]
- Hartwell HJ, Petrosky KY, Fox JG, et al. Prolactin prevents hepatocellular carcinoma by restricting innate immune activation of c-Myc in mice. Proc Natl Acad Sci U S A 2014;111:11455-60. [Crossref] [PubMed]
- Han H, Sun D, Li W, et al. A c-Myc-microRNA functional feedback loop affects hepatocarcinogenesis. Hepatology 2013;57:2378-89. [Crossref] [PubMed]
- Solé X, Hernández P, de Heredia ML, et al. Genetic and genomic analysis modeling of germline c-MYC overexpression and cancer susceptibility. BMC Genomics 2008;9:12. [Crossref] [PubMed]
- Lee KE, Chang BC, Park S, et al. Effects of single nucleotide polymorphisms in c-Myc on stable warfarin doses in patients with cardiac valve replacements. Pharmacogenomics 2015;16:1101-8. [Crossref] [PubMed]
- Franceschini N, Hu Y, Reiner AP, et al. Prospective associations of coronary heart disease loci in African Americans using the MetaboChip: the PAGE study. PLoS One 2014;9:e113203. [Crossref] [PubMed]
- Wang H, Lu J, Edmunds LR, et al. Coordinated activities of multiple Myc-dependent and Myc-independent biosynthetic pathways in hepatoblastoma. J Biol Chem 2016;291:26241-51. [Crossref] [PubMed]
- Purcell R, Childs M, Maibach R, et al. HGF/c-Met related activation of β-catenin in hepatoblastoma. J Exp Clin Cancer Res 2011;30:96. [Crossref] [PubMed]
- Gammons M, Bienz M. Multiprotein complexes governing Wnt signal transduction. Curr Opin Cell Biol 2018;51:42-9. [Crossref] [PubMed]