ORIGINAL RESEARCH article

Front. Med., 14 December 2021

Sec. Precision Medicine

Volume 8 - 2021 | https://doi.org/10.3389/fmed.2021.797195

YTHDF2 Gene rs3738067 A>G Polymorphism Decreases Neuroblastoma Risk in Chinese Children: Evidence From an Eight-Center Case-Control Study

  • 1. Department of Pediatric Surgery, Guangzhou Institute of Pediatrics, Guangdong Provincial Key Laboratory of Research in Structural Birth Defect Disease, Guangzhou Women and Children's Medical Center, Guangzhou Medical University, Guangzhou, China

  • 2. Department of Clinical Laboratory, Biobank, Harbin Medical University Cancer Hospital, Harbin, China

  • 3. Department of Pediatric Surgery, The Second Affiliated Hospital of Xi'an Jiaotong University, Xi'an, China

  • 4. Department of Pediatric Surgery, Hunan Children's Hospital, Changsha, China

  • 5. Department of Pediatric Surgery, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China

  • 6. Department of Pediatric Surgery, Shengjing Hospital of China Medical University, Shenyang, China

  • 7. Kunming Key Laboratory of Children Infection and Immunity, Yunnan Key Laboratory of Children's Major Disease Research, Yunnan Institute of Pediatrics Research, Yunnan Medical Center for Pediatric Diseases, Kunming Children's Hospital, Kunming, China

  • 8. Department of Hematology, The Second Affiliated Hospital and Yuying Children's Hospital of Wenzhou Medical University, Wenzhou, China

  • 9. Department of Pathology, Children Hospital and Women Health Center of Shanxi, Taiyuan, China

Abstract

Neuroblastoma is a primary malignancy mainly occurring in children. We have reported that polymorphisms of several N6-methyladenosine (m6A) RNA modification-related genes contributed to neuroblastoma risk in previous studies. YTHDF2, a “reader” of RNA m6A modification, is involved in cancer progression. Here, we estimated the association between a YTHDF2 gene rs3738067 A>G polymorphism and neuroblastoma susceptibility in 898 neuroblastoma patients and 1,734 healthy individuals from China. We found that the rs3738067 A>G could decrease neuroblastoma risk [AG vs. AA: adjusted odds ratio (OR) = 0.76, 95% confidence interval (CI) = 0.64–0.90, P = 0.002; AG/GG vs. AA: adjusted OR = 0.81, 95% CI = 0.69–0.95, P = 0.011). Besides, the rs3738067 AG/GG genotype was related to reduced neuroblastoma risk in the following subgroups: children aged 18 months and under, boys, patients with tumors originating from retroperitoneal, patients at clinical stage IV, and cases at clinical stages III plus IV. Importantly, false-positive report probability analysis proved our significant results worthy of close attention of. The expression quantitative trait locus analysis results revealed that the rs3738067 was associated with the expression of YTHDF2.

Introduction

Neuroblastoma is a prevalent malignancy originating from precursor cells of the sympathetic nervous system, and it mainly affects infants and children under 5 years of age (). Neuroblastoma with high aggressiveness often progresses quickly, leading to a disappointing prognosis and high recurrence rate. Although some patients experience mild or no treatment exhibit spontaneous regression (), more than half of patients with high-risk neuroblastoma die even with multimodality treatment (). Due to the complex nature of the disease, the pathogenesis of neuroblastoma is still far from clear. Increasing evidence suggests that the gradual accumulation of adverse genetic alterations leads to the transformation of normal cells to cancer cells (). Therefore, it is essential to uncover the detrimental genetic changes in neuroblastoma to screen for high-risk individuals and explore potentially effective treatment.

In recent years, researchers have achieved dramatic advancements in the genetic etiology of neuroblastoma (). Of note, genome-wide association studies (GWASs) have emerged as a powerful tool for exploring the causal genetic mechanisms of human diseases, including tumors (). Chromosome instability was considered as one of the major causes in neuroblastoma oncogenesis (). Two studies demonstrated that neuroblastoma shares common DNA variants with malignant cutaneous melanoma () and congenital heart disease (). Currently, some single-nucleotide polymorphisms (SNPs) related to neuroblastoma susceptibility have been identified by GWASs and studies with candidate gene strategy, including LMO1 (, ), METTL14 (), PARP1 (), MTHFR and VDR (). Nevertheless, the genetic variations known presently are not sufficient to fulfill the genetic landscape in neuroblastoma.

N6-methyladenosine (m6A) is the most popular post-transcriptional modification of RNAs in eukaryotes, particularly in messenger RNAs (mRNAs) (). RNA m6A modification is a dynamic and reversible process regulated by methyltransferases (known as writers) and demethylases (known as erasers). RNAs with m6A modifications can be recognized by some RNA binding proteins (named readers), which decide the different destinies of the modified RNA (). As a member of the YTH domain family, YTHDF2 functions as an m6A reader to modulate the translation, location, and stability of targeted mRNA (). Emerging evidence has suggested that dysregulated m6A modifications are tightly implicated in various diseases, especially cancers (). Many studies have demonstrated the involvement of YTHDF2 in the regulation of m6A modified targets in cancer development (). However, there are few reports about SNPs in the YTHDF2 gene and tumor risk.

We carried out a multi-center epidemiology study among Chinese children to analyze the association between the SNPs in the key m6A modification modulator gene YTHDF2 and neuroblastoma susceptibility.

Materials and Methods

Sample Selection

This work was conducted with the approval of the Institutional Review Board of Guangzhou Women and Children's Medical Center, with 898 neuroblastoma patients registered in eight hospitals (Guangzhou, Zhengzhou, Wenzhou, Xi'an, Taiyuan, Kunming, Changsha, Shenyang) in China and 1,734 age- and gender-matched healthy controls involved in previous studies (Supplementary Table S1) (, ). All participants have signed informed consent.

Polymorphism Selection and Genotyping

Only one potential functional SNP in the YTHDF2 gene (rs3738067 A>G) was chosen and genotyped in this study. Selection criteria and genotyping by TaqMan methodology were described previously (, ). The YTHDF2 gene rs3738067 A>G is located in transcription factor binding sites (TFBS) and might affect transcription activity as predicted by SNPinfo (https://snpinfo.niehs.nih.gov/snpinfo/snpfunc.html).

Statistical Analysis

The Chi-square test was applied to measure the compliance of alleles at individual loci in controls with the Hardy-Weinberg equilibrium (HWE) and the differences of selected demographic variables between patients and controls. Logistic regression analyses determined crude or adjusted odds ratios (ORs, adjusted for age and gender) with respective 95% confidence intervals (CIs) to analyze the association of YTHDF2 gene polymorphism with neuroblastoma risk. False-positive report probability (FPRP) analysis was applied to estimate the deserving attention of YTHDF2 gene polymorphism in neuroblastoma as described before (). In brief, three parameters were used to determine FPRP values, including statistical power, P-value, and prior probability representing a real association between the SNP and disease. We set 0.2 as an FPRP threshold and assigned a prior probability of 0.1 to detect an OR of 1.5 (for risk effects) or 0.67 (for protective effects) for the association of genotypes with neuroblastoma susceptibility. The association of the rs3738067 A>G with YTHDF2 expression was determined in the GTEx portal (https://www.gtexportal.org/home/) via eQTLs analysis. P < 0.05 was taken as statistically significant. Analyses were processed with SAS 9.1 (SAS Institute).

Results

Association of YTHDF2 rs3738067 A>G With Neuroblastoma Risk

The genotyping of YTHDF2 was successfully screened in 896 neuroblastoma patients and 1,733 controls. The genotype distribution of YTHDF2 rs3738067 A>G polymorphism and its relation to neuroblastoma susceptibility was indicated in Table 1. The frequency of the YTHDF2 rs3738067 A>G genotype coincided with HWE among the controls (HWE=0.359). The minor allele frequency (MAF) of YTHDF2 rs3738067 A>G polymorphism, was 0.2591 for the controls and 0.2410 for the cases. Based on the results of the single-locus analysis, we found that the G carriers of the rs3738067 were associated with decreased neuroblastoma risk (AG vs. AA: adjusted OR = 0.76, 95% CI = 0.64–0.90, P = 0.002; AG/GG vs. AA: adjusted OR = 0.81, 95% CI = 0.69–0.95, P = 0.011).

Table 1

GenotypeCases (%) (N = 896)Controls (%) (N = 1,733)PaCrude OR (95% CI)PAdjusted OR (95% CI)bPb
rs3738067 (HWE=0.359)
AA535 (59.71)944 (54.47)1.001.00
AG292 (32.59)680 (39.24)0.76 (0.64–0.90)0.0020.76 (0.64–0.90)0.002
GG69 (7.70)109 (6.29)1.12 (0.81–1.54)0.4981.12 (0.82–1.55)0.478
Additive0.1340.91 (0.79–1.03)0.1340.91 (0.80–1.03)0.143
Dominant361 (40.29)789 (45.53)0.0100.81 (0.69–0.95)0.0100.81 (0.69–0.95)0.011
Recessive827 (92.30)1,624 (93.71)0.1721.24 (0.91–1.70)0.1731.25 (0.91–1.71)0.165

YTHDF2 rs3738067 A>G polymorphism and neuroblastoma susceptibility.

OR, odds ratio; CI, confidence interval; HWE, Hardy-Weinberg equilibrium.

Values were in bold if the P-values <0.05 or the 95% CIs excluded 1.

a

χ2 test for genotype distributions between neuroblastoma cases and cancer-free controls.

b

Adjusted for age and gender.

Stratification Analysis

After that, we assessed the relation between YTHDF2 gene polymorphism and neuroblastoma susceptibility in subgroups classified via age, gender, sites of origins as well as clinical stages. As presented in Table 2, we detected that the rs3738067 AG/GG genotype carriers were linked to reduced neuroblastoma risk in subgroups of children with the age of 18 months and under (adjusted OR = 0.76, 95% CI = 0.58–0.98, P = 0.036), males (adjusted OR = 0.75, 95% CI = 0.60–0.93, P = 0.009), patients with retroperitoneal tumors (adjusted OR = 0.69, 95% CI = 0.54–0.88, P = 0.003), patients at clinical stage IV (adjusted OR = 0.68, 95% CI = 0.51–0.91, P = 0.009) and those at clinical stages III+IV (adjusted OR = 0.74, 95% CI = 0.59–0.93, P = 0.009).

Table 2

Variablesrs3738067 (cases/controls)OR (95% CI)PAOR (95% CI)aPa
AAAG/GG
Age, month
≤ 18204/373140/3400.75 (0.58–0.98)0.0330.76 (0.58–0.98)0.036
>18331/571221/4490.85 (0.69–1.05)0.1280.86 (0.70–1.06)0.158
Gender
Females238/415168/3290.89 (0.70–1.14)0.3530.90 (0.71–1.15)0.411
Males297/529193/4600.75 (0.60–0.93)0.0100.75 (0.60–0.93)0.009
Sites of origin
Adrenal gland141/944107/7890.91 (0.69–1.19)0.4810.91 (0.69–1.19)0.477
Retroperitoneal202/944116/7890.69 (0.54–0.88)0.0030.69 (0.54–0.88)0.003
Mediastinum121/94492/7890.91 (0.68–1.21)0.5180.92 (0.69–1.22)0.552
Others67/94438/7890.68 (0.45–1.02)0.0630.68 (0.45–1.03)0.066
Clinical stages
I179/944130/7890.87 (0.68–1.11)0.2610.87 (0.68–1.12)0.275
II91/94469/7890.91 (0.65–1.26)0.5590.92 (0.66–1.27)0.601
III98/94465/7890.79 (0.57–1.10)0.1660.80 (0.58–1.11)0.186
IV147/94484/7890.68 (0.52–0.91)0.0090.68 (0.51–0.91)0.009
4 s14/9444/7890.34 (0.11–1.04)0.0590.37 (0.12–1.15)0.085
I+II+4 s270/944199/7890.88 (0.72–1.08)0.2320.88 (0.72–1.09)0.239
III+IV245/944149/7890.73 (0.58–0.91)0.0060.74 (0.59–0.93)0.009

Stratify analysis for YTHDF2 rs3738067 A>G polymorphism and neuroblastoma susceptibility.

OR, odds ratio; CI, confidence interval; AOR, adjusted odds ratio.

Values were in bold if the P-values <0.05 or the 95% CIs excluded 1.

a

Adjusted for age and gender, omitting the corresponding stratify factor.

FPRP Analysis

An FPRP analysis was implemented to verify whether our significant findings deserve attentions. As shown in Table 3, the significant association for rs3738067 A>G (AG vs. AA, AG/GG vs. AA, males, retroperitoneal, clinical stage IV, and III+IV) remained noteworthy at the prior probability level of 0.1.

Table 3

GenotypeCrude OR(95% CI)PaStatistical powerbPrior probability
0.250.10.010.0010.0001
AG vs. AA0.76 (0.64–0.90)0.0020.9290.0050.0160.1530.6460.948
AG/GG vs. AA0.81 (0.69–0.95)0.0100.9860.0300.0860.5080.9130.991
≤ 18 months0.75 (0.58–0.98)0.0330.8080.1080.2670.8000.9760.998
Males0.75 (0.60–0.93)0.0100.8330.0330.0930.5300.9190.991
Retroperitoneal0.69 (0.54–0.88)0.0030.5820.0150.0430.3300.8330.980
Stage IV0.68 (0.52–0.91)0.0090.5610.0450.1240.6080.9400.994
Stage III+IV0.73 (0.58–0.91)0.0060.7650.0210.0610.4160.8780.986

False-positive report probability analysis for significant findings for the association between YTHDF2 rs3738067 A>G polymorphism and neuroblastoma susceptibility.

OR, odds ratio; CI, confidence interval.

The results were in bold if the false-positive report probability <0.200.

a

Chi-square test was used to calculate the genotype frequency distributions.

b

Statistical power was calculated using the number of observations in the subgroup and the OR and P-values in this table.

Effect of rs3738067 A>G on the Expression of YTHDF2

To confirm the functional relevance of rs3738067 A>G to the mRNA expression of YTHDF2, Cis-expression quantitative trait loci (eQTLs) analysis of the rs3738067 A>G and YTHDF2 expression was estimated using GTEx data. Results manifested that the rs3738067 A allele was related to increased YTHDF2 expression levels in the whole blood [Figure 1, P=1.9*10−5, normalized effect size (NES)=0.084].

Figure 1

Discussion

Though many genetic variants linking to neuroblastoma susceptibility have been recognized, further efforts are needed to fully understand this disease's genetic landscape. The present work verified the YTHDF2 rs3738067 A>G could reduce neuroblastoma risk for the first time. YTHDF2 is an m6A modification “reader” which recognizes m6A-modified mRNAs to modulate the translation and stability of targeted mRNA (). The roles of YTHDF2 in tumors are critical but controversial. For instance, As shown by Zhong et al., YTHDF2 restrained tumor cell growth in hepatocellular carcinoma (). Shen et al. reported that YTHDF2 repressed cell growth in gastric cancer through modulating FOXC2 expression (). In contrast, Li and colleagues disclosed that YTHDF2 promoted cell proliferation and migration in ovarian cancer (). However, the role of YTHDF2 in neuroblastoma remains largely unknown.

Increasing evidence has indicated that genetic variations, including SNPs in m6A modification modulators, correlate closely with cancer progression (). Also, a report has pointed that YTHDF2 rs3738067 A>G polymorphism exhibits a significant inverse association with glioma risk (). Our study evaluated the association of the YTHDF2 gene SNP (rs3738067 A>G) with neuroblastoma susceptibility. The results showed that rs3738067 AG/GG genotype was related to reduced neuroblastoma risk in several subgroups, including children aged 18 months and under, males, patients with tumors originating from retroperitoneal, patients in clinical stages IV, and patients in clinical stages III+IV. We also performed FPRP tests to confirm if the obtained associations were noteworthy or not to provide further evidence of the reliability of our results. The association of the rs3738067 A>G with YTHDF2 expression was determined in the GTEx portal via eQTLs analysis. The integrative analyses of eQTL and SNP information may provide more understanding about the complex disease-modulating network (). However, further studies are needed to substantiate the association between the rs3738067 A>G polymorphism and mRNA expression levels of YTHDF2.

There are several limitations to this study. First, only one SNP (rs3738067 A>G) in the YTHDF2 gene was evaluated. More studies will be performed to find other potential functional SNPs in the YTHDF2 gene. Second, only Chinese children were involved in this study. Thus, the results of this study may not be applicable to other ethnic groups. Moreover, as neuroblastoma is a multifactorial tumor, only genetic analysis is not enough to estimate neuroblastoma risk, and this study failed to incorporate environmental and genetic-environmental factors.

Our work for the first time verified the significant correlation of YTHDF2 gene rs3738067 A>G polymorphism with neuroblastoma risk, and this polymorphism is an intriguing locus for in-depth researches. However, the underlying biological mechanisms remain to be explored.

Funding

This study was supported by grants from the National Natural Science Foundation of China (No. 82173593), Natural Science Foundation of Guangdong Province (No. 2019A1515010360), and the Major Science and Technology Special Project of Wenzhou (No. ZY2020021). The Natural science foundation of Guangdong Province (No. 2020A1515011569). Guangzhou Science and Technology Innovation Commission (No. 201607010395). Guangzhou Health science and Technology Project (No. 20201A010018). Guangzhou Health science and Technology Project (No. 20211A011033).

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Ethics statement

The studies involving human participants were reviewed and approved by Institutional Review Board of Guangzhou Women and Children's Medical Center. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

Author contributions

HZe, TY, and JH contributed to conception and design of the study. ML and JL organized the original data. JZhu, JC, YL, JZha, ZY, LL, HZh, and SL provided the clinical tissue and blood samples for the study. HX and YZ provided some technical guidance. All authors contributed to the article and approved the submitted version.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2021.797195/full#supplementary-material

References

Summary

Keywords

neuroblastoma, YTHDF2, rs3738067, polymorphism, susceptibility

Citation

Zeng H, Li M, Liu J, Zhu J, Cheng J, Li Y, Zhang J, Yang Z, Li L, Zhou H, Li S, Xia H, Zou Y, He J and Yang T (2021) YTHDF2 Gene rs3738067 A>G Polymorphism Decreases Neuroblastoma Risk in Chinese Children: Evidence From an Eight-Center Case-Control Study. Front. Med. 8:797195. doi: 10.3389/fmed.2021.797195

Received

18 October 2021

Accepted

26 November 2021

Published

14 December 2021

Volume

8 - 2021

Edited by

Mario Capasso, University of Naples Federico II, Italy

Reviewed by

Alessandro Testori, University of Milano-Bicocca, Italy; Ferdinando Bonfiglio, University of Naples Federico II, Italy

Updates

Copyright

*Correspondence: Jing He Tianyou Yang

†These authors have contributed equally to this work

This article was submitted to Precision Medicine, a section of the journal Frontiers in Medicine

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics