High-throughput identification of regulatory elements and functional assays to uncover susceptibility genes for nasopharyngeal carcinoma.

Large-scale genetic association studies have identified multiple susceptibility loci for nasopharyngeal carcinoma (NPC), but the underlying biological mechanisms remain to be explored. To gain insights into the genetic etiology of NPC, we conducted a follow-up study encompassing 6,907 cases and 10,472 controls and identified two additional NPC susceptibility loci, 9q22.33 (rs1867277; OR = 0.74, 95% CI = 0.68-0.81, p = 3.08 × 10-11) and 17q12 (rs226241; OR = 1.42, 95% CI = 1.26-1.60, p = 1.62 × 10-8). The two additional loci, together with two previously reported genome-wide significant loci, 5p15.33 and 9p21.3, were investigated by high-throughput sequencing for chromatin accessibility, histone modification, and promoter capture Hi-C (PCHi-C) profiling. Using luciferase reporter assays and CRISPR interference (CRISPRi) to validate the functional profiling, we identified PHF2 at locus 9q22.33 as a susceptibility gene. PHF2 encodes a histone demethylase and acts as a tumor suppressor. The risk alleles of the functional SNPs reduced the expression of the target gene PHF2 by inhibiting the enhancer activity of its long-range (4.3 Mb) cis-regulatory element, which promoted proliferation of NPC cells. In addition, we identified CDKN2B-AS1 as a susceptibility gene at locus 9p21.3, and the NPC risk allele of the functional SNP rs2069418 promoted the expression of CDKN2B-AS1 by increasing its enhancer activity. The overexpression of CDKN2B-AS1 facilitated proliferation of NPC cells. In summary, we identified functional SNPs and NPC susceptibility genes, which provides additional explanations for the genetic association signals and helps to uncover the underlying genetic etiology of NPC development.

Functional TET2 gene polymorphisms increase the risk of neuroblastoma in Chinese children.

The 5-methylcytosine (m5C) is the key chemical modification in RNAs. As one of the demethylases in m5C, TET2 has been shown as a tumor suppressor. However, the impact of TET2 gene polymorphisms on neuroblastoma has not been elucidated. 402 neuroblastoma patients and 473 controls were genotyped for TET2 gene polymorphisms using the TaqMan method. The impact of TET2 gene polymorphisms on neuroblastoma susceptibility was determined using multivariate logistic regression analysis. We also adopted genotype-tissue expression database to explore the impact of TET2 gene polymorphisms on the expression of host and nearby genes. We used the R2 platform and Sangerbox tool to analyze the relationship between gene expression and neuroblastoma risk and prognosis through non-parametric testing and Kaplan-Meier analysis, respectively. We found the TET2 gene polymorphisms (rs10007915 G > C and rs7670522 A > C) and the combined 2-5 risk genotypes can significantly increase neuroblastoma risk. Stratification analysis showed that these significant associations were more prominent in certain subgroups. TET2 rs10007915 G > C and rs7670522 A > C are significantly associated with reduced expression of TET2 mRNA. Moreover, lower expression of TET2 gene is associated with high risk, MYCN amplification, and poor prognosis of neuroblastoma. The rs10007915 G > C and rs7670522 A > C are significantly related to the increased expression of inorganic pyrophosphatase 2 mRNA, and higher expression of PPA2 gene is associated with high risk, MYCN amplification, and poor prognosis of neuroblastomas. In summary, TET2 rs10007915 G > C and rs7670522 A > C significantly confer neuroblastoma susceptibility, and further research is needed to investigate the underlying mechanisms.

LMO family gene polymorphisms and Wilms tumor susceptibility in Chinese children: a five-center case-control study.

BACKGROUND: Wilms tumor is the most prevalent embryonal kidney malignancy in children worldwide. Previous genome-wide association study (GWAS) identified that LIM domain only 1 (LMO1) gene polymorphisms affected the susceptibility to develop certain tumor types. Apart from LMO1, the LMO gene family members also include LMO2-4, each of which has oncogenic potential. METHODS: We conducted this five-center case‒control study to assess the correlations between single nucleotide polymorphisms in LMO family genes and Wilms tumor susceptibility. Odds ratios and 95% confidence intervals were calculated to evaluate the strength of the association. RESULTS: We found LMO1 rs2168101 G > T and rs11603024 C > T as well as LMO2 rs7933499 G > A were significantly associated with Wilms tumor risk. Stratified analysis demonstrated a protective role of rs2168101 GT/TT genotypes against Wilms tumor in the subgroups of age ≤ 18 months, males and clinical stages I/II compared to the rs2168101 GG genotype. Nevertheless, carriers with the rs11603024 TT genotype were more likely to have an increased risk of Wilms tumor than those with rs11603024 CC/CT genotypes in age > 18 months. And the rs11603024 was identified as a protective polymorphism for reducing the risk of Wilms tumor in the sex- and gender- subgroup. Likewise, carriers with the rs7933499 GA/AA genotypes were at significantly elevated risk of Wilms tumor in age ≤ 18 months and clinical stages I/II. CONCLUSION: Overall, our study identified the importance of LMO family gene polymorphisms on Wilms tumor susceptibility in Chinese children. Further investigations are needed to validate our conclusions.

Quantitative detection of Epstein-Barr virus DNA methylation in the diagnosis of nasopharyngeal carcinoma by blind brush sampling.

Detecting EBV DNA load in nasopharyngeal (NP) brushing samples for the diagnosis of nasopharyngeal carcinoma (NPC) has attracted widespread attentions. Currently, NP brush sampling mostly relies on endoscopic guidance, and there are few reports on diagnostic markers suitable for nonguided conditions (blind brush sampling), which is of great significance for extending its application. One hundred seventy nasopharyngeal brushing samples were taken from 98 NPC patients and 72 non-NPC controls under the guidance of endoscope, and 305 blind brushing samples were taken without endoscopic guidance from 164 NPC patients and 141 non-NPC controls (divided into discovery and validation sets). Among these, 38 cases of NPC underwent both endoscopy-guided NP brushing and blind brushing. EBV DNA load targeting BamHI-W region and EBV DNA methylation targeting 11029 bp CpG site located at Cp-promoter region were detected by quantitative polymerase chain reaction (q-PCR). EBV DNA load showed good classification accuracy for NPC in endoscopy-guided brushing samples (AUC = 0.984). However, in blind bushing samples, the diagnostic performance was greatly reduced (AUC = 0.865). Unlike EBV DNA load, the accuracy of EBV DNA methylation was less affected by brush sampling methods, whether in endoscopy-guided brushing (AUC = 0.923) or blind brushing (AUC = 0.928 in discovery set and AUC = 0.902 in validation set). Importantly, EBV DNA methylation achieved a better diagnostic accuracy than EBV DNA load in blind brushing samples. Overall, detection of EBV DNA methylation with blind brush sampling shows great potential in the diagnosis of NPC and may facilitate its use in nonclinical screening of NPC.

Saliva biopsy: Detecting the difference of EBV DNA methylation in the diagnosis of nasopharyngeal carcinoma.

Saliva sampling is a non-invasive method, and could be performed by donors themselves. However, there are few studies reporting biomarkers in saliva in the diagnosis of NPC. A total of 987 salivary samples were used in this study. First, EBV DNA methylation was profiled by capture sequencing in the discovery cohort (n = 36). Second, a q-PCR based method was developed and five representative EBV DNA CpG sites (11 029 bp, 45 849 bp, 57 945 bp, 66 226 bp and 128 102 bp) were selected and quantified to obtain the methylated density in the validation cohort1 (n = 801). Third, a validation cohort2 (n = 108) was used to further verify the differences of EBV methylation in saliva. A significant increase of EBV methylation was found in NPC patients compared with controls. The methylated score of EBV genome obtained by capture sequencing could distinguish NPC from controls (sensitivity 90%, specificity 100%). Further, the methylated density of EBV DNA CpG sites revealed by q-PCR showed a good diagnostic performance. The sensitivity and specificity of detecting a single CpG site (11 029 bp) could reach 75.4% and 99.7% in the validation cohort1, and 78.2% and 100% in the validation cohort2. Besides, the methylated density of the CpG site was found to decrease below the COV in NPC patients after therapy, and increase above the COV after recurrence. Our study provides an appealing alternative for the non-invasive detection of NPC without clinical setting. It paves the way for conducting a home-based large-scale screening in the future.

Clinical and genome-wide association analysis of chemoradiation-induced hearing loss in nasopharyngeal carcinoma.

Chemoradiation-induced hearing loss (CRIHL) is one of the most devasting side effects for nasopharyngeal carcinoma (NPC) patients, which seriously affects survivors' long-term quality of life. However, few studies have comprehensively characterized the risk factors for CRIHL. In this study, we found that age at diagnosis, tumor stage, and concurrent cisplatin dose were positively associated with chemoradiation-induced hearing loss. We performed a genome-wide association study (GWAS) in 777 NPC patients and identified rs1050851 (within the exon 2 of NFKBIA), a variant with a high deleteriousness score, to be significantly associated with hearing loss risk (HR = 5.46, 95% CI 2.93-10.18, P = 9.51 × 10-08). The risk genotype of rs1050851 was associated with higher NFKBIA expression, which was correlated with lower cellular tolerance to cisplatin. According to permutation-based enrichment analysis, the variants mapping to 149 hereditary deafness genes were significantly enriched among GWAS top signals, which indicated the genetic similarity between hereditary deafness and CRIHL. Pathway analysis suggested that synaptic signaling was involved in the development of CRIHL. Additionally, the risk score integrating genetic and clinical factors can predict the risk of hearing loss with a relatively good performance in the test set. Collectively, this study shed new light on the etiology of chemoradiation-induced hearing loss, which facilitates high-risk individuals' identification for personalized prevention and treatment.

Genomic Landscapes of Epstein-Barr Virus in Pulmonary Lymphoepithelioma-Like Carcinoma.

Epstein-Barr virus (EBV) infection is associated with multiple malignancies, including pulmonary lymphoepithelioma-like carcinoma (pLELC), a particular subtype of primary lung cancer. However, the genomic characteristics of EBV related to pLELC remain unclear. Here, we obtained the whole-genome data set of EBV isolated from 78 pLELC patients and 37 healthy controls using EBV-captured sequencing. Compared with the reference genome (NC_007605), a total of 3,995 variations were detected across pLELC-derived EBV sequences, with the mutational hot spots located in latent genes. Combined with 180 published EBV sequences derived from healthy people in Southern China, we performed a genome-wide association study and identified 32 variations significantly related to pLELC (P < 2.56 × 10-05, Bonferroni correction), with the top signal of single nucleotide polymorphism (SNP) coordinate T7327C (OR = 1.22, P = 2.39 × 10-15) locating in the origin of plasmid replication (OriP). The results of population structure analysis of EBV isolates in East Asian showed the EBV strains derived from pLELC were more similar to those from nasopharyngeal carcinoma (NPC) than other EBV-associated diseases. In addition, typical latency type-II infection were recognized for EBV of pLELC at both transcription and methylation levels. Taken together, we defined the global view of EBV genomic profiles in pLELC patients for the first time, providing new insights to deepening our understanding of this rare EBV-associated primary lung carcinoma. IMPORTANCE Pulmonary lymphoepithelioma-like carcinoma (pLELC) is a rare, distinctive subtype of primary lung cancer closely associated with Epstein-Barr virus (EBV) infection. Here, we gave the first overview of pLELC-derived EBV at the level of genome, methylation and transcription. We obtained the EBV sequences data set from 78 primary pLELC patients, and revealed the sequences diversity across EBV genome and detected variability in known immune epitopes. Genome-wide association analysis combining 217 healthy controls identifies significant variations related to the risk of pLELC. Meanwhile, we characterized the integration landscapes of EBV at the genome-wide level. These results provided new insight for understanding EBV's role in pLELC tumorigenesis.

Peptidome-wide association analysis of Epstein-Barr virus identifies epitope repertoires associated with nasopharyngeal carcinoma.

Human leukocyte antigen (HLA) molecules are essential for presenting Epstein-Barr virus (EBV) antigens and are closely related to nasopharyngeal carcinoma (NPC). This study aims to systematically investigate the association between HLA-bound EBV peptides and NPC risk through in silico HLA-peptide binding prediction. A total of 455 NPC patients and 463 healthy individuals in NPC endemic areas were recruited, and HLA-target sequencing was performed. HLA-peptide binding prediction for EBV, followed by peptidome-wide logistic regression and motif analysis, was applied. Binding affinity changes for EBV peptides carrying high-risk mutations were analyzed. We found that NPC-associated EBV peptides were significantly enriched in immunogenic proteins and core linkage disequilibrium (LD) proteins related to evolution, especially those binding HLA-A alleles (p = 3.10 × 10-4 for immunogenic proteins and p = 8.10 × 10-5 for core LD proteins related to evolution). These peptides were clustered and showed binding motifs of HLA supertypes, among which supertype A02 presented an NPC-risk effect (padj = 3.77 × 10-4 ) and supertype A03 presented an NPC-protective effect (padj = 4.89 × 10-4 ). Moreover, a decreased binding affinity toward risk HLA supertype A02 was observed for the peptide carrying the NPC-risk mutation BNRF1 V1222I (p = 0.0078), and an increased binding affinity toward protective HLA supertype A03 was observed for the peptide carrying the NPC-risk mutation BALF2 I613V (p = 0.022). This study revealed the distinct preference of EBV peptides for binding HLA supertypes, which may contribute to shaping EBV population structure and be involved in NPC development.

The interaction and mediation effects between the host genetic factors and Epstein-Barr virus VCA-IgA in the risk of nasopharyngeal carcinoma.

Previous studies have demonstrated strong associations between host genetic factors and Epstein-Barr virus (EBV) VCA-IgA with the risk of nasopharyngeal carcinoma (NPC). However, the specific interplay between host genetics and EBV VCA-IgA on NPC risk is not well understood. In this two-stage case-control study (N = 4804), we utilized interaction and mediation analysis to investigate the interplay between host genetics (genome-wide association study-derived polygenic risk score [PRS]) and EBV VCA-IgA antibody level in the NPC risk. We employed a four-way decomposition analysis to assess the extent to which the genetic effect on NPC risk is mediated by or interacts with EBV VCA-IgA. We consistently found a significant interaction between the PRS and EBV VCA-IgA on NPC risk (discovery population: synergy index [SI] = 2.39, 95% confidence interval [CI] = 1.85-3.10; replication population: SI = 3.10, 95% CI = 2.17-4.44; all pinteraction < 0.001). Moreover, the genetic variants included in the PRS demonstrated similar interactions with EBV VCA-IgA antibody. We also observed an obvious dose-response relationship between the PRS and EBV VCA-IgA antibody on NPC risk (all ptrend < 0.001). Furthermore, our decomposition analysis revealed that a substantial proportion (approximately 90%) of the genetic effects on NPC risk could be attributed to host genetic-EBV interaction, while the risk effects mediated by EBV VCA-IgA antibody were weak and statistically insignificant. Our study provides compelling evidence for an interaction between host genetics and EBV VCA-IgA antibody in the development of NPC. These findings emphasize the importance of implementing measures to control EBV infection as a crucial strategy for effectively preventing NPC, particularly in individuals at high genetic risk.

Genomic landscape of Epstein-Barr virus in familial nasopharyngeal carcinoma.

To better understand the genomic characteristics of Epstein-Barr virus (EBV) in familial nasopharyngeal carcinoma (NPC), we sequenced the EBV genomes by whole-genome capture in 38 unrelated patients with NPC family history in first-degree relatives and 47 healthy controls, including 13 with family history and 34 without. Compared with type 1 reference genome, mutation hotspots were observed in the latent gene regions of EBV in familial NPC cases. Population structure analysis showed that one cluster has a higher frequency in familial cases than in controls (OR=5.33, 95 % CI 2.50-11.33, P=1.42×10-5), and similar population structure composition was observed among familial and sporadic NPC cases in high-endemic areas. By genome-wide association analysis, four variants were found to be significantly associated with familial NPC. Consistent results were observed in the meta-analysis integrating two published case-control EBV sequencing studies in NPC high-endemic areas. High-risk haplotypes of EBV composed of 34 variants were associated with familial NPC risk (OR=13.85, 95 % CI 4.13-46.44, P=2.06×10-5), and higher frequency was observed in healthy blood-relative controls with NPC family history (9/13, 69.23 %) than those without family history (16/34, 47.06%). This study suggested the potential contribution of EBV high-risk subtypes to familial aggregation of NPC.

Association between HLA alleles and Epstein-Barr virus Zta-IgA serological status in healthy males from southern China.

BACKGROUND: Nasopharyngeal carcinoma (NPC), an Epstein-Barr virus (EBV) associated cancer, exhibits an extremely high incidence in southern Chinese. Given that human leukocyte antigen (HLA) plays critical roles in antigen presentation and relates to NPC susceptibility, it is speculated that certain HLA variants may affect EBV reactivation, which is a key pathogenic factor of NPC. Therefore, we attempted to identify HLA alleles associated with the indicator of EBV reactivation, Zta-IgA, in healthy males from NPC endemic area. METHODS: HLA alleles of 1078 healthy males in southern China from the 21-RCCP study were imputed using genome-wide single nucleotide polymorphism data. EBV Zta-IgA in blood samples were measured using an enzyme-linked immunosorbent assay. Multiple logistic regression analysis was used to evaluate the effect of HLA allele on Zta-IgA serological status and its potential joint association with smoking. The binding affinity for Zta-peptide was predicted using NetMHCIIpan 4.0. RESULTS: HLA-DRB1*09:01 was found to be associated with a higher risk of Zta-IgA seropositivity (odds ratio = 1.80, 95% confidence interval = 1.32-2.45; p = 1.82 × 10-4 ). Compared with non-smokers without HLA-DRB1*09:01, the effect size increased to 2.19- and 3.70-fold for the light and heavy smokers carrying HLA-DRB1*09:01, respectively. Furthermore, HLA-DRB1*09:01 showed a stronger binding affinity to Zta peptide than other HLA-DRB1 alleles. CONCLUSIONS: Our study highlighted the pivotal role of genetic HLA variants in EBV reactivation and the etiology of NPC. Smokers with HLA-DRB1*09:01 have a significantly higher risk of being Zta-IgA seropositive, which indicates the necessity of smoking cessation in certain high-risk populations and also provide clues for further research on the etiology of NPC.

Nasopharyngeal Epstein-Barr virus DNA loads in high-risk nasopharyngeal carcinoma families: Familial aggregation and host heritability.

Nasopharyngeal carcinoma (NPC), the most common head and neck cancer, is characterized by distinct geographic distribution and familial aggregation. Multiple risk factors, including host genetics, environmental factor, and EBV infection, have been linked to the development of NPC, particularly in the familial clustering cases. However, the cause of NPC endemicity remains enigmatic due possibly to the complicated interplay between these risk factors. Recently, positive Epstein-Barr virus (EBV) DNA loads at nasopharyngeal (NP) cavity has been found to reflect NPC development and applied in NPC screening. To examine whether the increased NP EBV loads could aggregate in the families and be affected by host genetics and environmental factor, EBV loads were obtained by 510 NP brushing samples from eligible unaffected individuals, who have two or more relatives affected with NPC, in 116 high-risk NPC families. The correlation of relative pairs was estimated using S.A.G.E. (version 6.4, 2016), and host heritability of NP EBV loads was calculated with variance component models using SOLAR (version 8.4.2, 2019). In result, significant correlations of EBV loads were observed between parent-offspring pairs and sibling-sibling pairs (P < .001), but not in distant kin relationship pairs. Interestingly, after excluding the shared environmental factor within families, host genetics contributes significantly to NP EBV loads with a heritability of 56.41% (P = 1.00 × 10-7 ), and its effect was slightly elevated (68.86%, P = 3.40 × 10-6 ) in families with more NPC cases (≥3). These findings indicate that additional host-genetic variants involved in the EBV local NP mucosal behavior may be especially important for the development of NPC.

ALKBH1 rs2267755 C>T polymorphism decreases neuroblastoma risk in Chinese children.

Neuroblastoma is a highly malignant extracranial solid tumor in pediatrics. ALKBH1 as a recently discovered DNA N6-methyldeoxyadenosine (6mA) demethylase closely links to tumorigenesis. Whether the ALKBH1 polymorphism contributes to neuroblastoma risk remains unclear. In the present study, we genotyped the ALKBH1 single nucleotide polymorphisms (SNPs) in 402 neuroblastoma patients and 473 healthy controls by TaqMan assay. Odds ratios (ORs) and 95% confidence intervals (CIs) were also calculated to evaluate the strength of the association. Our result exhibited that the rs2267755 C>T (CT vs. CC, adjusted OR=0.69, 95% CI=0.50-0.94, P=0.019) is significantly associated with reduced neuroblastoma risk. And its protective effect is particularly significant in children with tumors originating from the retroperitoneal. Combined genotype analysis revealed that carriers with 1-2 protective genotypes are more susceptible to neuroblastoma than those with 3-4 protective genotypes (adjusted OR=0.71, 95% CI=0.53-0.97, P=0.028). Moreover, the rs2267755 C>T is significantly associated with messenger RNA (mRNA) expression of ALKBH1 and three of its surrounding genes, including SNWQ, ADCK1, and RPL21P10. These results suggest that the rs2267755 C>T may be a genetic variant to reduce neuroblastoma risk.

Genetic variants of m7G modification genes influence neuroblastoma susceptibility.

Objective: Neuroblastoma is a life-threatening pediatric solid tumor whose etiology remains unclear. N7-methylguanosine (m7G) is one of the most important epigenetic modifications of RNA, which plays a crucial role in tumorigenesis. The m7G-mediated genes METTL1 and WDR4 also have been reported to be dysregulated in various cancers. However, the implications of METTL1 and WDR4 in neuroblastoma have not been clarified. Methods: Given the oncogenic potential of m7G modification, we performed a case-control study to assess the association of METTL1 and WDR4 genes polymorphisms with neuroblastoma risk in a Chinese population consisting of 402 cases and 473 controls. Odds ratios (ORs) and 95 % confidence intervals (CIs) were applied to evaluate the associations between studied polymorphisms and neuroblastoma risk. The adjusted odds ratio (AOR) was adjusted for age and gender. Results: Overall, four polymorphisms were significantly associated with neuroblastoma risk, including METTL1 rs2291617 (recessive model: adjusted OR = 1.59, 95 % CI = 1.08-2.34, P = 0.019), WDR4 rs2156316 (dominant model: adjusted OR = 0.74, 95 % CI = 0.57-0.97, P = 0.028), WDR4 rs6586250 (dominant model: adjusted OR = 0.59, 95 % CI = 0.42-0.84, P = 0.004) and WDR4 rs15736 (dominant model: adjusted OR = 0.60, 95 % CI = 0.42-0.85, P = 0.004). Stratified analysis showed stronger correlations between significant polymorphisms and neuroblastoma risk among subgroups divided by age, gender, tumor origin, and clinical stage. Furthermore, expression quantitative trait loci (eQTL) analysis revealed that significant polymorphisms were associated with the expression of the adjacent genes. Conclusions: Our study indicated that four polymorphisms in m7G-mediated genes contribute to neuroblastoma susceptibility in the eastern Chinese population. However, our findings should be verified further by large-scale and well-designed studies.

N6-methyladenosine modification of circMARK2 enhances cytoplasmic export and stabilizes LIN28B, contributing to the progression of Wilms tumor.

BACKGROUND: The potential involvement of circular RNAs (circRNAs) and N6-methyladenosine (m6A) modification in the progression of Wilms tumor (WT) has not been fully elucidated. This study investigates the regulatory mechanisms and clinical significance of m6A-modified circMARK2 and its role in WT progression. METHODS: We identified dysregulated circRNAs through deep sequencing and validated their expression by qRT-PCR in WT tissues. The biological functions of circMARK2 were assessed using clone formation, transwell migration, and orthotopic animal models. To dissect the underlying mechanisms, we employed RNA immunoprecipitation, RNA pull-down, dual-luciferase reporter assays, Western blotting, and immunofluorescence and immunohistochemical staining. RESULTS: CircMARK2, upregulated in WT tissues, was found to be m6A-modified and promoted cytoplasmic export. It facilitated WT progression by stabilizing LIN28B mRNA through the circMARK2/IGF2BP2 interaction. In vitro and in vivo studies demonstrated that circMARK2 enhances the malignant behavior of WT cells. Clinically, higher circMARK2 levels in tumor tissues of WT patients were linked to increased tumor aggressiveness and reduced survival rates. CONCLUSIONS: Our study provides the first comprehensive evidence that m6A-modified circMARK2 contributes to WT progression by enhancing LIN28B mRNA stability, promoting cellular aggressiveness. CircMARK2 emerges as a potential biomarker for prognosis and a promising target for therapeutic intervention in WT, underscoring the clinical relevance of m6A modification in pediatric renal cancer.

EBV DNA methylation profiles and its application in distinguishing nasopharyngeal carcinoma and nasal NK/T-cell lymphoma.

BACKGROUND: As an oncovirus, EBV is associated with multiple cancers, including solid tumors and hematological malignancies. EBV methylation plays an important role in regulating tumor occurrence. However, the EBV methylation profiles in EBV-associated tumor tissues are poorly understood. RESULTS: In this study, EBV methylation capture sequencing was conducted in several different tumor tissue samples, including NPC, EBVaGC, lung LELC and parotid LELC. Besides, EBV capture sequencing and following qMSP were performed on nasopharyngeal brushing samples from NPC and nasal NKTCL patients. Our results showed that the EBV genome among different types of tumors displayed specific methylation patterns. Among the four types of tumors from epithelial origin (NPC, EBVaGC, lung LELC and parotid LELC), the most significant differences were found between EBVaGC and the others. For example, in EBVaGC, all CpG sites within 1,44,189-1,45,136 bp of the EBV genome sequence on gene RPMS1 were hyper-methylated compared to the others. Differently, significant differences of EBV CpG sites, particularly those located on gene BILF2, were observed between NPC and nasal NKTCL patients in nasopharyngeal brushing samples. Further, the methylated level of BILF2 was further detected using qMSP, and a diagnostic model distinguishing NPC and nasal NKTCL was established. The AUC of the model was 0.9801 (95% CI 0.9524-1.0000), with the sensitivity and specificity of 98.81% (95% CI 93.63-99.94%) and 76.92% (95% CI 49.74-91.82%), respectively. CONCLUSIONS: Our study reveals more clues for further understanding the pathogenesis of EBV, and provides a possibility for distinguishing EBV-related tumor by detecting specific EBV CpG sites.

A comprehensive predictive model for radiation-induced brain injury in risk stratification and personalized radiotherapy of nasopharyngeal carcinoma.

BACKGROUND AND PURPOSE: Radiation-induced brain injury (RBI) is a severe radiotoxicity for nasopharyngeal carcinoma (NPC) patients, greatly affecting their long-term life quality and survival. We aim to establish a comprehensive predictive model including clinical factors and newly developed genetic variants to improve the precision of RBI risk stratification. MATERIALS AND METHODS: By performing a large registry-based retrospective study with magnetic resonance imaging follow-up on RBI development, we conducted a genome-wide association study and developed a polygenic risk score (PRS) for RBI in 1189 NPC patients who underwent intensity-modulated radiotherapy. We proposed a tolerance dose scheme for temporal lobe radiation based on the risk predicted by PRS. Additionally, we established a nomogram by combining PRS and clinical factors for RBI risk prediction. RESULTS: The 38-SNP PRS could effectively identify high-risk individuals of RBI (P = 1.42 × 10-34). Based on genetic risk calculation, the recommended tolerance doses of temporal lobes should be 57.6 Gy for individuals in the top 10 % PRS subgroup and 68.1 Gy for individuals in the bottom 50 % PRS. Notably, individuals with high genetic risk (PRS > P50) and receiving high radiation dose in the temporal lobes (D0.5CC > 65 Gy) had an approximate 50-fold risk over individuals with low PRS and receiving low radiation dose (HR = 50.09, 95 %CI = 24.27-103.35), showing an additive joint effect (Pinteraction < 0.001). By combining PRS with clinical factors including age, tumor stage, and radiation dose of temporal lobes, the predictive accuracy was significantly improved with C-index increased from 0.78 to 0.85 (P = 1.63 × 10-2). CONCLUSIONS: The PRS, together with clinical factors, could improve RBI risk stratification and implies personalized radiotherapy.

NSUN2 gene rs13181449 C>T polymorphism reduces neuroblastoma risk.

Neuroblastoma is the most common tumor in infants. RNA m5C modification regulates the survival, differentiation, and migration of cells affecting RNA function. However, the effects of the m5C modification methyltransferase gene NSUN2 polymorphism on neuroblastoma susceptibility have not been reported. TaqMan method was used to determine genotypes of four NSUN2 polymorphisms (rs4702373 C>T, rs13181449 C>T, rs166049 T>G, and rs8192120 A>C) in 402 patients with neuroblastoma and 473 cancer-free controls from Jiangsu province, China. Odds ratio (OR) and 95% confidence interval (CI) were used to evaluate the association of NSUN2 polymorphisms with neuroblastoma susceptibility. The association was also further assessed in subgroups stratified by age, sex, tumor origin, and stage. GTEx was used to analyze the effect of these polymorphisms on NSUN2 expression. We found the rs13181449 C>T was significantly associated with reduced neuroblastoma risk (CT vs. CC: adjusted OR = 0.68, 95% CI = 0.51-0.92, P = 0.012; CT/TT vs. CC: adjusted OR = 0.70, 95% CI = 0.53-0.92, P = 0.010). Compared with 0-2 protective genotypes, those with 3-4 protective genotypes could significantly reduce the neuroblastoma risk (adjusted OR = 0.68, 95% CI = 0.52 to 0.90, P = 0.006). Stratification analysis showed that the protective effect of rs13181449 polymorphism remained significant in children with age >18 months, boys, and those with early INSS stages. Moreover, children with more protective genotypes in the same subgroups also exhibited significantly reduced neuroblastoma risk. GTEx analysis showed that the rs13181449 T genotype was related with decreased NSUN2 gene expression. In conclusions, NSUN2 rs13181449 polymorphism is associated with decreased neuroblastoma risk, and the underlying mechanism in neuroblastoma needs further study.

TRMT6 gene rs236110 C > A polymorphism increases the risk of Wilms tumor.

tRNA methyltransferase 6 (TRMT6)is an enzyme catalyzing N1-methyladenosine, a reversible modification in RNA, including tRNA, mRNA, rRNA, and lncRNA. Increasing evidence has shown the implications of this post-transcriptional modification and its regulators in carcinogenesis. However, its roles in Wilms tumor haven't been reported. In this study, four TRMT6 gene polymorphisms (rs236170 A > G, rs451571 T > C, rs236188 G > A, and rs236110 C > A) were tested for association with susceptibility to Wilms tumor, the most frequently diagnosed pediatric renal tumor. TaqMan method was adopted to analyze the genotypes of these polymorphisms in 414 cases and 1199 controls. Among the four TRMT6 gene polymorphisms, only the rs236110 C > A displayed a significant association with the risk of Wilms tumor [AA vs. CC, adjusted odds ratio (OR) = 1.93, 95 % confidence interval (CI) = 1.14-3.27, P = 0.015]. This association was confirmed under the recessive models (AA vs. CC/CA, OR = 1.92, 95 % CI = 1.14-3.23, P = 0.015). Furthermore, after stratifying by age, gender, and clinical stage, we mainly detected significant associations for the rs236110 C > A in children older than 18 months, boys, and those with stage IV or III + IV diseases. The rs236110 A allele was significantly associated with decreased expression of MCM8. In conclusion, we identified the rs236110 C > A in the TRMT6 gene as a Wilms tumor susceptibility locus, and this polymorphism warrants more validation studies to be translated into individualized risk prediction strategies for children.

WDR4 gene polymorphisms and Wilms tumor susceptibility in Chinese children: A five-center case-control study.

Wilms tumor is the most common embryonal renal malignancy in children. WDR4 is an indispensable noncatalytic subunit of the RNA N7-methylguanosine (m7G) methyltransferase complex and plays an essential role in tumorigenesis. However, the relationship between polymorphisms in the WDR4 gene and susceptibility to Wilms tumor remains to be fully investigated. We performed a large case-control study involving 414 patients and 1199 cancer-free controls to investigate whether single nucleotide polymorphisms (SNPs) in the WDR4 gene are associated with Wilms tumor susceptibility. WDR4 gene polymorphisms (rs2156315 C > T, rs2156316 C > G, rs6586250 C > T, rs15736 G > A, and rs2248490 C > G) were genotyped using the TaqMan assay. In addition, unconditioned logistic regression analysis was performed, odds ratios (ORs) and 95% confidence intervals (CIs) were used to assess the association between WDR4 gene SNPs and Wilms tumor susceptibility as well as the strength of the associations. We found that only the rs6586250 C>T polymorphism was significantly associated with an increased risk of Wilms tumor (adjusted OR=2.99, 95% CI = 1.28-6.97, P = 0.011 for the rs6586250 TT genotype; adjusted OR=3.08, 95% CI = 1.33-7.17, P = 0.009 for the rs6586250 CC/CT genotype). Furthermore, the stratification analysis revealed that patients with the rs6586250 TT genotype and carriers with 1-5 risk genotypes exhibited statistically significant associations with increased Wilms tumor risk in specific subgroups. However, the rs2156315 CT/TT genotype was identified as having a protective effect against Wilms tumor in the age >18 months subgroup compared with the rs2156315 CC genotype. In brief, our study demonstrated that the rs6586250 C > T polymorphism of the WDR4 gene was significantly associated with Wilms tumor. This finding may contribute to the understanding of the genetic mechanism of Wilms tumor.