Exome-wide analyses identify low-frequency variant in CYP26B1 and additional coding variants associated with esophageal squamous cell carcinoma

PTPN23[Thr] variant reduces susceptibility and tumorigenesis in esophageal squamous cell carcinoma through dephosphorylation of EGFR

Post-translational modifications (PTMs) have emerged as pivotal regulators of the development of cancers, including esophageal squamous cell carcinoma (ESCC). Here, we conducted a comprehensive analysis of PTM-related genetic variants associated with ESCC risk using large-scale genome-wide and exome-wide association datasets. We observed significant enrichment of PTM-related variants in the ESCC risk loci and identified five variants that were significantly associated with ESCC risk. Among them, rs6780013 in PTPN23 exhibited the highest level of significance in ESCC susceptibility in 9,728 ESCC cases and 10,977 controls (odds ratio [OR] = 0.85, 95 % confidence interval [CI] = 0.81- 0.89, P = 9.77 × 10-14). Further functional investigations revealed that PTPN23[Thr] variant binds to EGFR and modulates its phosphorylation at Thr699. PTPN23[Thr] variant substantially inhibited ESCC cell proliferation both in vitro and in vivo. Our findings underscore the critical role of PTPN23[Thr]-EGFR interaction in ESCC development, providing more insights into the pathogenesis of this cancer.

Systematical identification of key genes and regulatory genetic variants associated with prognosis of esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) stands as a highly lethal malignancy characterized by pronounced recurrence and metastasis, resulting in a bleak 5-year survival rate. Despite extensive investigations, encompassing genome-wide association studies, the identification of robust prognostic markers has remained elusive. In this study, leveraging four independent data sets comprising 404 ESCC patients, we conducted a systematic analysis to unveil pivotal genes influencing overall survival. our meta-analysis identified 278 genes significantly associated with ESCC prognosis. Further exploration of the prognostic landscape involved an examination of expression quantitative trait loci for these genes, leading to the identification of six tag single nucleotide polymorphisms predictive of overall survival in a cohort of 904 ESCC patients. Notably, functional annotation spotlighted rs11227223, residing in the enhancer region of nuclear paraspeckle assembly transcript 1 (NEAT1), as a crucial variant likely exerting a substantive biological role. Through a series of biochemistry experiments, we conclusively demonstrated that the rs11227223-T allele, indicative of a poorer prognosis, augmented NEAT1 expression. Our results underscore the substantive role of NEAT1 and its regulatory variant in prognostic predictions for ESCC. This comprehensive analysis not only advances our comprehension of ESCC prognosis but also unveils a potential avenue for targeted interventions, offering promise for enhanced clinical outcomes.

Metabolites, Healthy Lifestyle, and Polygenic Risk Score Associated with Upper Gastrointestinal Cancer: Findings from the UK Biobank Study

Previous metabolomics studies have highlighted the predictive value of metabolites on upper gastrointestinal (UGI) cancer, while most of them ignored the potential effects of lifestyle and genetic risk on plasma metabolites. This study aimed to evaluate the role of lifestyle and genetic risk in the metabolic mechanism of UGI cancer. Differential metabolites of UGI cancer were identified using partial least-squares discriminant analysis and the Wilcoxon test. Then, we calculated the healthy lifestyle index (HLI) score and polygenic risk score (PRS) and divided them into three groups, respectively. A total of 15 metabolites were identified as UGI-cancer-related differential metabolites. The metabolite model (AUC = 0.699) exhibited superior discrimination ability compared to those of the HLI model (AUC = 0.615) and the PRS model (AUC = 0.593). Moreover, subgroup analysis revealed that the metabolite model showed higher discrimination ability for individuals with unhealthy lifestyles compared to that with healthy individuals (AUC = 0.783 vs 0.684). Furthermore, in the genetic risk subgroup analysis, individuals with a genetic predisposition to UGI cancer exhibited the best discriminative performance in the metabolite model (AUC = 0.770). These findings demonstrated the clinical significance of metabolic biomarkers in UGI cancer discrimination, especially in individuals with unhealthy lifestyles and a high genetic risk.

Single-cell analysis identified POSTN+ cells associated with the aggressive phenotype and risk of esophageal squamous cell carcinoma

Tumors are intricate and heterogeneous systems characterized by mosaic cancer cell populations with diverse expression profiles. Leveraging single-cell technologies, we employed the Scissor algorithm to delineate an epithelial subpopulation associated with the aggressive phenotype in esophageal squamous cell carcinoma (ESCC). This identified subpopulation exhibited elevated expression of genes involved in critical pathways, such as epithelial-mesenchymal transition and PI3K-Akt. Key signature genes within this subpopulation, namely CAV1, COL3A1, COL6A1, POSTN, and TAGLN, demonstrated significant upregulation concomitant with both tumorigenesis and tumor progression across independent single-cell datasets. Furthermore, we selected 1,450 expression quantitative trait loci of the top 62 signature genes of this cell subpopulation to investigate their potential in predicting ESCC risk. The results showed that the POSTN loci were predominantly associated with ESCC susceptibility. Through functional annotation and replication analyses, we identified that the rs1028728 in the POSTN promoter was significantly associated with increased ESCC risk in 7,049 ESCC cases and 8,063 controls (odds ratio = 1.29, 95% confidence interval: 1.18-1.42, p = 4.03 × 10-8). Subsequent biochemical experiments showed that the rs1028728[T] allele enhanced POSTN expression by affecting the binding of PRRX1 in the POSTN promoter. In summary, our meticulous single-cell analysis delineates an invasive epithelial subpopulation in ESCC, with POSTN emerging as an important marker for the aggressive phenotype. These findings offer more insights into potential strategies for the prevention and intervention of ESCC, enriching our understanding of this complex cancer landscape.

An atlas of causal and mechanistic drivers of interpatient heterogeneity in glioma

Grade IV glioma, formerly known as glioblastoma multiforme (GBM) is the most aggressive and lethal type of brain tumor, and its treatment remains challenging in part due to extensive interpatient heterogeneity in disease driving mechanisms and lack of prognostic and predictive biomarkers. Using mechanistic inference of node-edge relationship (MINER), we have analyzed multiomics profiles from 516 patients and constructed an atlas of causal and mechanistic drivers of interpatient heterogeneity in GBM (gbmMINER). The atlas has delineated how 30 driver mutations act in a combinatorial scheme to causally influence a network of regulators (306 transcription factors and 73 miRNAs) of 179 transcriptional "programs", influencing disease progression in patients across 23 disease states. Through extensive testing on independent patient cohorts, we share evidence that a machine learning model trained on activity profiles of programs within gbmMINER significantly augments risk stratification, identifying patients who are super-responders to standard of care and those that would benefit from 2 nd line treatments. In addition to providing mechanistic hypotheses regarding disease prognosis, the activity of programs containing targets of 2 nd line treatments accurately predicted efficacy of 28 drugs in killing glioma stem-like cells from 43 patients. Our findings demonstrate that interpatient heterogeneity manifests from differential activities of transcriptional programs, providing actionable strategies for mechanistically characterizing GBM from a systems perspective and developing better prognostic and predictive biomarkers for personalized medicine.

Genetic linkage analysis of head and neck cancer in a Spanish family

OBJECTIVES

To describe the genetic variants that may be associated with the development of head and neck cancer (HNC) and functionally validating the molecular implications.

MATERIALS AND METHODS

A prospective observational study was carried out on a family of 3 generations in which 3 members had developed HNC. Peripheral blood sample was taken in a routine procedure for exome sequencing in one relative and genotyping in the remaining twelve relatives. For the functional analysis all-trans retinoic acid (atRA) was extracted from saliva and serum and measured using ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). The presence of HPV-DNA.

RESULTS

None of the patients smoked or consumed alcohol. The presence of HPV DNA was not detected in any of the biopsied samples. A total amount of 6 members out of 13 (46.15%) carried out the same mutation of CYP26B1 (2p13.2; G>T). The mean plasma concentration of atRA was 3.3109 ± 1.4791 pg/mL for the study family and 4.7370 ± 1.5992 pg/mL for the controls (p = 0.042).

CONCLUSION

Lower levels of atRA were confirmed in the study family, which may open the way to the possible relationship between the polymorphism CYP26B1 (2p13.2; G>T) and HNC.

Segregation of nascent GPCRs in the ER-to-Golgi transport by CCHCR1 via direct interaction

G protein-coupled receptors (GPCRs) constitute the largest superfamily of cell surface signaling proteins that share a common structural topology. When compared with agonist-induced internalization, how GPCRs are sorted and delivered to functional destinations after synthesis in the endoplasmic reticulum (ER) is much less well understood. Here, we demonstrate that depletion of coiled-coil α-helical rod protein 1 (CCHCR1) by siRNA and CRISPR-Cas9 significantly inhibits surface expression and signaling of α2A-adrenergic receptor (α2A-AR; also known as ADRA2A), without affecting α2B-AR. Further studies show that CCHCR1 depletion specifically impedes α2A-AR export from the ER to the Golgi, but not from the Golgi to the surface. We also demonstrate that CCHCR1 selectively interacts with α2A-AR. The interaction is mediated through multiple domains of both proteins and is ionic in nature. Moreover, mutating CCHCR1-binding motifs significantly attenuates ER-to-Golgi export, surface expression and signaling of α2A-AR. Collectively, these data reveal a novel function for CCHCR1 in intracellular protein trafficking, indicate that closely related GPCRs can be sorted into distinct ER-to-Golgi transport routes by CCHCR1 via direct interaction, and provide important insights into segregation and anterograde delivery of nascent GPCR members.

Genome-wide association study of esophageal squamous cell cancer identifies shared and distinct risk variants in African and Chinese populations

Esophageal squamous cell carcinoma (ESCC) has a high disease burden in sub-Saharan Africa and has a very poor prognosis. Genome-wide association studies (GWASs) of ESCC in predominantly East Asian populations indicate a substantial genetic contribution to its etiology, but no genome-wide studies have been done in populations of African ancestry. Here, we report a GWAS in 1,686 African individuals with ESCC and 3,217 population-matched control individuals to investigate its genetic etiology. We identified a genome-wide-significant risk locus on chromosome 9 upstream of FAM120A (rs12379660, p = 4.58 × 10-8, odds ratio = 1.28, 95% confidence interval = 1.22-1.34), as well as a potential African-specific risk locus on chromosome 2 (rs142741123, p = 5.49 × 10-8) within MYO1B. FAM120A is a component of oxidative stress-induced survival signals, and the associated variants at the FAM120A locus co-localized with highly significant cis-eQTLs in FAM120AOS in both esophageal mucosa and esophageal muscularis tissue. A trans-ethnic meta-analysis was then performed with the African ESCC study and a Chinese ESCC study in a combined total of 3,699 ESCC-affected individuals and 5,918 control individuals, which identified three genome-wide-significant loci on chromosome 9 at FAM120A (rs12379660, pmeta = 9.36 × 10-10), chromosome 10 at PLCE1 (rs7099485, pmeta = 1.48 × 10-8), and chromosome 22 at CHEK2 (rs1033667, pmeta = 1.47 × 10-9). This indicates the existence of both shared and distinct genetic risk loci for ESCC in African and Asian populations. Our GWAS of ESCC conducted in a population of African ancestry indicates a substantial genetic contribution to ESCC risk in Africa.

Genome-wide enhancer-gene regulatory maps link causal variants to target genes underlying human cancer risk

Genome-wide association studies have identified numerous variants associated with human complex traits, most of which reside in the non-coding regions, but biological mechanisms remain unclear. However, assigning function to the non-coding elements is still challenging. Here we apply Activity-by-Contact (ABC) model to evaluate enhancer-gene regulation effect by integrating multi-omics data and identified 544,849 connections across 20 cancer types. ABC model outperforms previous approaches in linking regulatory variants to target genes. Furthermore, we identify over 30,000 enhancer-gene connections in colorectal cancer (CRC) tissues. By integrating large-scale population cohorts (23,813 cases and 29,973 controls) and multipronged functional assays, we demonstrate an ABC regulatory variant rs4810856 associated with CRC risk (Odds Ratio = 1.11, 95%CI = 1.05-1.16, P = 4.02 × 10-5) by acting as an allele-specific enhancer to distally facilitate PREX1, CSE1L and STAU1 expression, which synergistically activate p-AKT signaling. Our study provides comprehensive regulation maps and illuminates a single variant regulating multiple genes, providing insights into cancer etiology.

Wnt/β-Catenin Promotes the Osteoblastic Potential of BMP9 Through Down-Regulating Cyp26b1 in Mesenchymal Stem Cells

BACKGROUND

All-trans retinoic acid (ATRA) promotes the osteogenic differentiation induced by bone morphogenetic protein 9 (BMP9), but the intrinsic relationship between BMP9 and ATRA keeps unknown. Herein, we investigated the effect of Cyp26b1, a critical enzyme of ATRA degradation, on the BMP9-induced osteogenic differentiation in mesenchymal stem cells (MSCs), and unveiled possible mechanism through which BMP9 regulates the expression of Cyp26b1.

METHODS

ATRA content was detected with ELISA and HPLC-MS/MS. PCR, Western blot, and histochemical staining were used to assay the osteogenic markers. Fetal limbs culture, cranial defect repair model, and micro-computed tomographic were used to evaluate the quality of bone formation. IP and ChIP assay were used to explore possible mechanism.

RESULTS

We found that the protein level of Cyp26b1 was increased with age, whereas the ATRA content decreased. The osteogenic markers induced by BMP9 were increased by inhibiting or silencing Cyp26b1 but reduced by exogenous Cyp26b1. The BMP9-induced bone formation was enhanced by inhibiting Cyp26b1. The cranial defect repair was promoted by BMP9, which was strengthened by silencing Cyp26b1 and reduced by exogenous Cyp26b1. Mechanically, Cyp26b1 was reduced by BMP9, which was enhanced by activating Wnt/β-catenin, and reduced by inhibiting this pathway. β-catenin interacts with Smad1/5/9, and both were recruited at the promoter of Cyp26b1.

CONCLUSIONS

Our findings suggested the BMP9-induced osteoblastic differentiation was mediated by activating retinoic acid signalling, viadown-regulating Cyp26b1. Meanwhile, Cyp26b1 may be a novel potential therapeutic target for the treatment of bone-related diseases or accelerating bone-tissue engineering.

Helicobacter pylori plus N-Methyl-N’-nitro-N-nitrosoguanidine: DNA damage and repair; malignant transformation of human esophageal epithelial cells

N-Methyl-N'-nitro-N-nitrosoguanidine (MNNG), found in pickled foods and in chlorinated water, has been used to induce malignant transformation and gastrointestinal cancer in rats. Helicobacter pylori (HP) is implicated in human gastric cancer and possibly also in esophageal cancer. These two agents - one chemical and the other biological - might act together to induce esophageal cancer. In this study, human esophageal epithelial cells (HEECs) were divided into four groups: HP, MNNG, HP + MNNG, and control. The HP-to-HEEC ratio was 100:1. Cells were exposed for 6 h and then passaged until malignant transformation. HEEC at early, intermediate, and late stages of malignant transformation were used for proliferation, cell-cycle, and invasion assays. The alkaline comet assay was performed and expression of proteins, including γ-H2AX and PAXX, was studied by western blotting, to explore DNA damage and repair processes. Measurements of cell morphology, soft-agar clone formation, and invasiveness, and a nude mouse xenograft model, were used to examine malignancy. The effect of HP was stronger than that of MNNG. The combination HP + MNNG exerted a stronger malignant transformation effect than either HP or MNNG alone. Mechanisms of this combined carcinogenesis may include promotion of cell proliferation, perturbation of the cell cycle, promotion of invasiveness, DNA double-strand break induction, or PAXX inhibition.

Tenascin-X as a causal gene for classical-like Ehlers-Danlos syndrome

Tenascin-X (TNX) is an extracellular matrix glycoprotein for which a deficiency results in a recessive form of classical-like Ehlers-Danlos syndrome (clEDS), a heritable connective tissue disorder with hyperextensible skin without atrophic scarring, joint hypermobility, and easy bruising. Notably, patients with clEDS also suffer from not only chronic joint pain and chronic myalgia but also neurological abnormalities such as peripheral paresthesia and axonal polyneuropathy with high frequency. By using TNX-deficient (Tnxb−/−) mice, well-known as a model animal of clEDS, we recently showed that Tnxb−/− mice exhibit hypersensitivity to chemical stimuli and the development of mechanical allodynia due to the hypersensitization of myelinated A-fibers and activation of the spinal dorsal horn. Pain also occurs in other types of EDS. First, we review the underlying molecular mechanisms of pain in EDS, especially that in clEDS. In addition, the roles of TNX as a tumor suppressor protein in cancer progression have been reported. Recent in silico large-scale database analyses have shown that TNX is downregulated in various tumor tissues and that high expression of TNX in tumor cells has a good prognosis. We describe what is so far known about TNX as a tumor suppressor protein. Furthermore, some patients with clEDS show delayed wound healing. Tnxb−/− mice also exhibit impairment of epithelial wound healing in corneas. TNX is also involved in liver fibrosis. We address the molecular mechanism for the induction of COL1A1 by the expression of both a peptide derived from the fibrinogen-related domain of TNX and integrin α11.

Evolutionary Origins of Metabolic Reprogramming in Cancer

Metabolic changes that facilitate tumor growth are one of the hallmarks of cancer. These changes are not specific to tumors but also take place during the physiological growth of tissues. Indeed, the cellular and tissue mechanisms present in the tumor have their physiological counterpart in the repair of tissue lesions and wound healing. These molecular mechanisms have been acquired during metazoan evolution, first to eliminate the infection of the tissue injury, then to enter an effective regenerative phase. Cancer itself could be considered a phenomenon of antagonistic pleiotropy of the genes involved in effective tissue repair. Cancer and tissue repair are complex traits that share many intermediate phenotypes at the molecular, cellular, and tissue levels, and all of these are integrated within a Systems Biology structure. Complex traits are influenced by a multitude of common genes, each with a weak effect. This polygenic component of complex traits is mainly unknown and so makes up part of the missing heritability. Here, we try to integrate these different perspectives from the point of view of the metabolic changes observed in cancer.

Expression and function analysis of CRABP2 and FABP5, and their ratio in esophageal squamous cell carcinoma

Objective

The purpose of this study was to explore the effect of CRABP2 and FABP5, and their ratio on prognosis in esophageal squamous cell carcinoma.

Methods

The expression data of CRABP2 in esophageal cancer in TCGA and GEO were collected by the public database GEPIA. The expression levels of CRABP2 and FABP5 were examined using immunohistochemistry. The relationship between the two proteins and related clinicopathological parameters were analyzed by χ² test. Survival analysis was used to investigate the effect of CRABP2 and FABP5, and their ratio on prognosis.

Results

Compared with normal esophageal mucosal epithelium, there was lower CRABP2 gene mRNA in the esophageal cancer tissue, and the difference was statistically significant (p < 0.01). For the expression level, no significant difference was observed in patients with stages I–IV in esophageal cancer. Immunohistochemistry showed that CRABP2 and FABP5 were both highly expressed in normal esophageal squamous epithelial cells at 100 and 94.1%, while lower in ESCC (75.6 and 58.7%). There was a significant difference in the expression between cancer and adjacent tissues (p < 0.001). No inherent relationship was manifested between the CRABP2 expression and the clinical parameters of the ESCC. The expression of FABP5 was related to lymph node metastasis (p = 0.032), the depth of invasion (p = 0.041), and the AJCC stage (p = 0.013). The ratio of CRABP2 and FABP5 was related to ethnicity (p = 0.001), nerve invasion (p = 0.031), and postoperative treatment (p = 0.038). CRABP2 is positively associated with FABP5 (r = 0.156, p = 0.041) and the ratio (r = 0.334, p = 0.000), while there was a negative correlation between FABP5 and the ratio (r = −0.269, p = 0.000). Patients with CRABP2-positive expression had a significantly longer overall survival than patients with CRABP2-negative expression (p = 0.025).

Conclusion

CRABP2 as a suppressor factor is expected to be a potential prognosis marker for esophageal squamous cell carcinoma.

Research on esophageal cancer: With personal perspectives from studies in China and Kenya

The most common form of esophageal cancer (EC), esophageal squamous cell carcinoma (ESCC), is prevalent in many unindustrialized societies, among people with lower socioeconomic status and those who frequently use tobacco and alcohol. In some areas, ESCC mortality ranked top among all cancer. In this review, we begin with discussions of the extensive research on EC in Linxian in northern China that started 60 years ago and the recent studies in Kenya from our personal perspectives. Based on the results obtained from these studies and information from the literature, we summarize our current understanding about the risk factors for ESCC including lifestyle factors (smoking, alcohol, consumption of food and beverages at high temperature and other unhealthy habits), poor diet and nutritional insufficiencies and genetic susceptibility. Elimination or minimization of these environmental risk factors, as well as early detection and treatment of precancerous lesions, would be effective means for the prevention of ESCC. Current knowledge of molecular alterations in ESCC (gene mutations, hypermethylation and amplification or overexpression), as well as treatment of ESCC and the potential of targeted therapy, are also discussed. Finally, we propose effective approaches for the prevention of ESCC by adapting a healthy lifestyle, including a healthy diet that would also prevent other diseases. Community outreach, public education and international collaboration are important for achieving this public health goal.

Global spectrum of population-specific common missense variation in cytochrome P450 pharmacogenes

Next-generation sequencing technology has afforded the discovery of many novel variants that are of significance to inheritable pharmacogenomics (PGx) traits but a large proportion of them have unknown consequences. These include missense variants resulting in single amino acid substitutions in cytochrome P450 (CYP) proteins that can impair enzyme function, leading to altered drug efficacy and toxicity. While most unknown variants are rare, an overlooked minority are variants that are collectively rare but enriched in specific populations. Here, we analyzed sequence variation data in 141,456 individuals from across eight study populations in gnomAD for 38 CYP genes to identify such variants in addition to common variants. By further comparison with data from two PGx-specific databases (PharmVar and PharmGKB) and ClinVar, we identified 234 missense variants in 35 CYP genes, of which 107 were unknown to these databases. Most unknown variants (n = 83) were population-specific common variants and several (n = 7) were found in important CYP pharmacogenes (CYP2D6, CYP4F2, and CYP2C19). Overall, 29% (n = 31) of 107 unknown variants were predicted to affect CYP enzyme function although further biochemical characterization is necessary. These variants may elucidate part of the unexplained interpopulation differences observed in drug response.

Aberrant MCM10 SUMOylation induces genomic instability mediated by a genetic variant associated with survival of esophageal squamous cell carcinoma

BACKGROUND

Esophageal squamous cell carcinoma (ESCC) is one of the common gastrointestinal malignancy with an inferior prognosis outcome. DNA replication licensing aberration induced by dysregulation of minichromosome maintenance proteins (MCMs) causes genomic instability and cancer metastasis. SUMOylation modification plays a pivotal role in regulation of genomic integrity, while its dysregulation fueled by preexisting germline variants in cancers remains poorly understood.

METHODS

Firstly, we conducted two-stage survival analysis consisting of an exome-wide association study in 904 ESCC samples and another independent 503 ESCC samples. Then, multipronged functional experiments were performed to illuminate the potential biological mechanisms underlying the promising variants, and MCM10 influences the ESCC progression. Finally, we tested the effects of MCM10 inhibitors on ESCC cells.

RESULTS

A germline variant rs2274110 located at the exon 15 of MCM10 was identified to be significantly associated with the prognosis of ESCC patients. Individuals carrying rs2274110-AA genotypes confer a poor survival (hazard ratio = 1.61, 95% confidence interval = 1.35-1.93, p = 1.35 × 10-7 ), compared with subjects carrying rs2274110-AG/GG genotypes. Furthermore, we interestingly found that the variant can increase SUMOylation levels at K669 site (Lys[K]699Arg[R]) of MCM10 protein mediated by SUMO2/3 enzymes, which resulted in an aberrant overexpression of MCM10. Mechanistically, aberrant overexpression of MCM10 facilitated the proliferation and metastasis abilities of ESCC cells in vitro and in vivo by inducing DNA over-replication and genomic instability, providing functional evidence to support our population finding that high expression of MCM10 is extensively presented in tumor tissues of ESCC and correlated with inferior survival outcomes of multiple cancer types, including ESCC. Finally, MCM10 inhibitors Suramin and its analogues were revealed to effectively block the metastasis of ESCC cells.

CONCLUSIONS

These findings not only demonstrate a potential biological mechanism between aberrant SUMOylation, genomic instability and cancer metastasis, but also provide a promising biomarker aiding in stratifying ESCC individuals with different prognosis, as well as a potential therapeutic target MCM10.

Oncogenic SNORD12B activates the AKT-mTOR-4EBP1 signaling in esophageal squamous cell carcinoma via nucleus partitioning of PP-1α

Esophageal cancer is a complex malignancy and the sixth leading cause of cancer death worldwide. In Eastern Asia including China, about 90% of all incident cases have esophageal squamous cell carcinoma (ESCC). Mounting evidence elucidates that aberrant expression of various non-coding RNAs (ncRNAs) contributes to ESCC progression, but it remains unclear how small nucleolar RNAs (snoRNAs) are involved in ESCC development. We systemically screened clinically relevant snoRNAs in ESCC via integrative analyses of The Cancer Genome Atlas (TCGA) data and validation in ESCC tissues. We found that snoRNA SNORD12B was one of the most evidently upregulated snoRNAs in ESCC specimens and its high expression was significantly associated with poor prognosis of patients. SNORD12B profoundly promoted proliferation, migration, invasion, and metastasis of ESCC cells in vitro and in vivo, indicating its oncogene nature. In particular, SNORD12B could interact with PP-1α, one of the three catalytic subunits of serine/threonine protein phosphatase 1, which is a major phosphatase that directly dephosphorylates AKT to suppress its activation. Interestingly, high levels of SNORD12B in ESCC cells could break interactions between 14-3-3ζ and PP-1α, abolish the retention of PP-1α in the cytosol by 14-3-3ζ and relocate PP-1α from the cytosol to the nucleus. This led to sequestered PP-1α in the nucleus, enhanced phosphorylation of AKT in the cytosol, activated AKT-mTOR-4EBP1 signaling, and, thus, ESCC progression. These insights would improve our understanding of how snoRNAs contribute to tumorigenesis and highlight the potential of snoRNAs as future therapeutic targets against cancers.

Comprehensive Study of Germline Mutations and Double-Hit Events in Esophageal Squamous Cell Cancer

Esophageal squamous cell cancer (ESCC) is the eighth most common cancer around the world. Several reports have focused on somatic mutations and common germline mutations in ESCC. However, the contributions of pathogenic germline alterations in cancer susceptibility genes (CSGs), highly frequently mutated CSGs, and pathogenically mutated CSG-related pathways in ESCC remain unclear. We obtained data on 571 ESCC cases from public databases and East Asian from the 1000 Genomes Project database and the China Metabolic Analytics Project database to characterize pathogenic mutations. We detected 157 mutations in 75 CSGs, accounting for 25.0% (143/571) of ESCC cases. Six genes had more than five mutations: TP53 (n = 15 mutations), GJB2 (n = 8), BRCA2 (n = 6), RECQL4 (n = 6), MUTYH (n = 6), and PMS2 (n = 5). Our results identified significant differences in pathogenic germline mutations of TP53, BRCA2, and RECQL4 between the ESCC and control cohorts. Moreover, we identified 84 double-hit events (16 germline/somatic double-hit events and 68 somatic/somatic double-hit events) occurring in 18 tumor suppressor genes from 83 patients. Patients who had ESCC with germline/somatic double-hit events were diagnosed at younger ages than patients with the somatic/somatic double-hit events, though the correlation was not significant. Fanconi anemia was the most enriched pathway of pathogenically mutated CSGs, and it appeared to be a primary pathway for ESCC predisposition. The results of this study identified the underlying roles that pathogenic germline mutations in CSGs play in ESCC pathogenesis, increased our awareness about the genetic basis of ESCC, and provided suggestions for using highly mutated CSGs and double-hit features in the early discovery, prevention, and genetic counseling of ESCC.

Proteomic analysis of psoriatic skin lesions in a Chinese population

Psoriasis is a chronic skin disorder with undefined pathogenesis. Several biomarkers for this disease have been identified by proteomic analysis. We explored the whole-proteomic changes in 45 pairs of psoriatic and adjacent noninvolved skin tissues in a Chinese population. A total of 3686 proteins were identified, of which 3008 were quantified. A total of 102 and 124 proteins were upregulated and downregulated in lesional skin, respectively. SART1 (P = 3.55 × 10^-5) and GLTP (P = 1.54 × 10^-3) were the most significantly down- and upregulated proteins. Nearly 90% of these differentially regulated proteins exhibited the same expression trends as those in an online RNA sequencing dataset for psoriasis; 19 differentially regulated proteins exhibited a negative relationship with DNA methylation data for psoriatic lesions. The differentially expressed proteins were enriched in ribosomes, antigen processing and presentation, immune response, and IL-17 signalling pathways. This study identified multiple differentially regulated proteins in psoriatic lesions, which suggested that changes in the proteome play important regulatory roles in psoriasis-associated processes. SIGNIFICANCE: Proteomic analysis was performed in 45 pairs of psoriatic and adjacent noninvolved skin tissues in a Chinese population. More than 3000 proteins were quantified, of which 226 were differentially expressed in psoriatic skin tissues. These proteins were mainly enriched in the immune response, antigen processing and presentation and IL-17 signalling pathways, which have been reported to be associated with the pathogenesis of psoriasis.

Role of Genetic Variation in Cytochromes P450 in Breast Cancer Prognosis and Therapy Response

Breast cancer is the most frequent cancer in the female population worldwide. The role of germline genetic variability in cytochromes P450 (CYP) in breast cancer prognosis and individualized therapy awaits detailed elucidation. In the present study, we used the next-generation sequencing to assess associations of germline variants in the coding and regulatory sequences of all human CYP genes with response of the patients to the neoadjuvant cytotoxic chemotherapy and disease-free survival (n = 105). A total of 22 prioritized variants associating with a response or survival in the above evaluation phase were then analyzed by allelic discrimination in the large confirmation set (n = 802). Associations of variants in CYP1B1, CYP4F12, CYP4X1, and TBXAS1 with the response to the neoadjuvant cytotoxic chemotherapy were replicated by the confirmation phase. However, just association of variant rs17102977 in CYP4X1 passed the correction for multiple testing and can be considered clinically and statistically validated. Replicated associations for variants in CYP4X1, CYP24A1, and CYP26B1 with disease-free survival of all patients or patients stratified to subgroups according to therapy type have not passed a false discovery rate test. Although statistically not confirmed by the present study, the role of CYP genes in breast cancer prognosis should not be ruled out. In conclusion, the present study brings replicated association of variant rs17102977 in CYP4X1 with the response of patients to the neoadjuvant cytotoxic chemotherapy and warrants further research of genetic variation CYPs in breast cancer.

Identifying Novel Susceptibility Genes for Colorectal Cancer Risk From a Transcriptome-Wide Association Study of 125,478 Subjects

BACKGROUND AND AIMS

Susceptibility genes and the underlying mechanisms for the majority of risk loci identified by genome-wide association studies (GWAS) for colorectal cancer (CRC) risk remain largely unknown. We conducted a transcriptome-wide association study (TWAS) to identify putative susceptibility genes.

METHODS

Gene-expression prediction models were built using transcriptome and genetic data from the 284 normal transverse colon tissues of European descendants from the Genotype-Tissue Expression (GTEx), and model performance was evaluated using data from The Cancer Genome Atlas (n = 355). We applied the gene-expression prediction models and GWAS data to evaluate associations of genetically predicted gene-expression with CRC risk in 58,131 CRC cases and 67,347 controls of European ancestry. Dual-luciferase reporter assays and knockdown experiments in CRC cells and tumor xenografts were conducted.

RESULTS

We identified 25 genes associated with CRC risk at a Bonferroni-corrected threshold of P < 9.1 × 10-6, including genes in 4 novel loci, PYGL (14q22.1), RPL28 (19q13.42), CAPN12 (19q13.2), MYH7B (20q11.22), and MAP1L3CA (20q11.22). In 9 known GWAS-identified loci, we uncovered 9 genes that have not been reported previously, whereas 4 genes remained statistically significant after adjusting for the lead risk variant of the locus. Through colocalization analysis in GWAS loci, we additionally identified 12 putative susceptibility genes that were supported by TWAS analysis at P < .01. We showed that risk allele of the lead risk variant rs1741640 affected the promoter activity of CABLES2. Knockdown experiments confirmed that CABLES2 plays a vital role in colorectal carcinogenesis.

CONCLUSIONS

Our study reveals new putative susceptibility genes and provides new insight into the biological mechanisms underlying CRC development.

Unique roles of rare variants in the genetics of complex diseases in humans

Genome-wide association studies have identified >10,000 genetic variants associated with various phenotypes and diseases. Although the majority are common variants, rare variants with >0.1% of minor allele frequency have been investigated by imputation and using disease-specific custom SNP arrays. Rare variants sequencing analysis mainly revealed have played unique roles in the genetics of complex diseases in humans due to their distinctive features, in contrast to common variants. Unique roles are hypothesis-free evidence for gene causality, a precise target of functional analysis for understanding disease mechanisms, a new favorable target for drug development, and a genetic marker with high disease risk for personalized medicine. As whole-genome sequencing continues to identify more rare variants, the roles associated with rare variants will also increase. However, a better estimation of the functional impact of rare variants across whole genome is needed to enhance their contribution to improvements in human health.

HMGA2 gene polymorphisms and Wilms tumor susceptibility in Chinese children: a four-center case-control study

Wilms tumor is a kidney malignancy that typically occurs in children. Aberrant expression of HMGA2 gene is commonly seen in many malignant tumors. Yet, HMGA2 gene polymorphisms on Wilms tumor risk are not established. We carried out the first four-center case-control study with 355 patients and 1,070 controls to assess the association of HMGA2 polymorphisms (rs6581658 A>G, rs8756 A>C, and rs968697 T>C) with Wilms tumor risk. All of these three polymorphisms in single could not impact Wilms tumor risk. Stratified analysis revealed a contributing Wilms tumor risk role of rs968697 TC/CC in subgroup of male (TC/CC vs. TT: adjusted odds ratio [OR] = 1.46, 95% confidence interval [CI] = 1.03-2.08, P = 0.035). However, we found that presence of 1-3 protective genotypes were less likely to develop tumor in subgroup of female (adjusted OR = 0.69, 95% CI = 0.48-0.99, P = 0.045). Our findings suggest that HMGA2 gene polymorphisms might influence Wilms tumor predisposition in a weak manner, under certain circumstances.

Genetic variants in m6A modification genes are associated with esophageal squamous-cell carcinoma in the Chinese population

N 6-methyladenosine (m6A) is an abundant modification in RNAs that affects RNA metabolism, and it is reported to be closely related to cancer occurrence and metastasis. In this study, we focused on evaluating the associations between genetic variants in m6A modification genes and the risk of esophageal squamous-cell carcinoma (ESCC). By integrating data of our previous genome-wide association studies and the predictions of several annotation tools, we identified a single nucleotide polymorphism, rs2416282 in the promoter of YTHDC2, that was significantly associated with the susceptibility of ESCC (odds ratio = 0.84, 95% CI: 0.77-0.92, P = 2.81 × 10-4). Through further functional experiments in vitro, we demonstrated that rs2416282 regulated YTHDC2 expression. Knockdown of YTHDC2 substantially promoted the proliferation rate of ESCC cells by affecting several cancer-related signaling pathways. Our results suggested that rs2416282 contributed to ESCC risk by regulating YTHDC2 expression. This study provided us a valuable insight into the roles of genetic variants in m6A modification genes for ESCC susceptibility and may contribute to the prevention of this disease in the future.

Assessment of PPARGC1A, PPARGC1B, and PON1 Genetic Polymorphisms in Esophageal Squamous Cell Carcinoma Susceptibility in the Eastern Chinese Han Population: A Case-Control Study Involving 2351 Subjects

Previous studies suggested that alterations in the energy metabolism might be underlying cancer initiation and progression. Polymorphisms of genes involved in energy metabolism regulation, such as peroxisome proliferator-activated receptor gamma coactivator 1α (PPARGC1A), -β (PPARGC1B), and paraoxonase 1 (PON1), might confer susceptibility to esophageal squamous cell carcinoma (ESCC) and partially explain its pathogenesis. We investigated the effects of several single nucleotide polymorphisms (SNPs) in three metabolic-related genes (e.g., PPARGC1A, PPARGC1B, and PON1) on ESCC susceptibility. In total, 829 patients with sporadic ESCC and 1522 nontumor controls were enrolled in the study. SNPs were genotyped using PCR-ligase detection reaction. Our study revealed that the PPARGC1A rs3736265 G/A SNP significantly increased the risk for ESCC (GA vs. GG: adjusted odds ratio [OR] = 1.25, 95% confidence interval [95% CI] = 1.02-1.54, p = 0.034; GA+AA vs. GG: adjusted OR = 1.25, 95% CI = 1.03-1.52, p = 0.027]. In addition, a stratified analysis revealed that the PPARGC1A rs3736265 SNP was correlated with the development of ESCC in male and nondrinking subgroups. We also confirmed that the PPARGC1B rs17572019 G/A SNP promoted the risk of ESCC in subgroup with high alcohol intake. The PPARGC1A rs8192678 C/T polymorphism decreased the susceptibility of ESCC in men. These findings highlight that polymorphisms in PPARGC1A and PPARGC1B may contribute to ESCC susceptibility. In the future, further well-designed epidemiological studies are needed to confirm our findings.

Identification of Leukocyte telomere length-related genetic variants contributing to predisposition of Esophageal Squamous Cell Carcinoma

Background: Cancers may arise from cells with dysregulated telomeric functions due to shorten telomere length. We and others previously found that short leukocyte telomere length was associated with markedly evaluated risk of esophageal squamous cell carcinoma (ESCC). Hence, we hypothesized that single nucleotide polymorphisms (SNPs) associated with shorter telomere length may contribute to ESCC predisposition.

Methods: We systematically evaluated association between seven candidate seven SNPs (CXCR4 rs6430612, TERT rs13172201, TERT rs10069690, TERT rs2853676, TERT rs451360, OBFC1 rs4387287, and VPS34 rs2162440) and ESCC risk in two case-control sets consisting of 1588 ESCC cases and 1600 controls. Logistic regression models were utilized to estimate associations between SNPs and ESCC susceptibility and odds ratios (ORs) and their 95% confidence intervals (95% CIs) were computed.

Results: We firstly identified three SNPs (rs6430612, rs13172201 and rs4387287) which are significantly associated with telomere length in Chinese populations (all P<0.05). Importantly, CXCR4 rs6430612 and OBFC1 rs4387287 polymorphisms significantly confer reduced risk of ESCC (P=1.7×10^-7 and P=3.9×10^-5). On the contrary, we observed an evidently increased risk for ESCC development associated with TERT rs13172201 genetic variant (P=2.2×10^-4).

Conclusions: In summary, rs6430612, rs13172201 and rs4387287 might be key genetic components in complicated regulation of telomere length and contributing to ESCC predisposition. Our results elucidate the prevalent involvement of genetic variants in telomere biology and further provide pathogenic insights into the role of telomeres in cancer development.

Establishment and validation of an immune-based prognostic score model in glioblastoma

BACKGROUND

Immune escape is one of the landmark features of glioblastoma (GBM). Immunotherapy is undoubtedly a revolution in the field of tumor treatment, especially the application of immune checkpoint inhibitors and CAR-T cells, which have achieved amazing results in fighting against cancer. This study aimed to establish a TP53-related immune-based score model to improve the prognostic of GBM by investigating the gene mutations and the immune landscape of GBM.

METHODS

Data were obtained from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA) databases. Differentially expressed genes (DEGs) analysis between the TP53 mutated (TP53^MUT) and wild-type (TP53^WT) GBM patients was conducted. The CIBERSORT algorithm was applied to evaluate the proportion of immune cell types and RNA sequencing (RNA-seq) data from the TCGA and CGGA were used as discovery and validation cohorts, respectively, to build and validate an immune-related prognostic model (IPM). Genes in the IPM model were first screened by univariate Cox analysis, then filtered by the least absolute shrinkage and selection operator (LASSO) Cox regression method to eliminate collinearity among DEGs. A nomogram was finally established and evaluated by combining both the IPM and other clinical factors.

RESULTS

PTEN was the top most mutated gene in GBM patients (118/393), followed by TP53 (116/393). 332 immune-related genes were identified and the immune response in the TP53^WT group was remarkably greater than in the TP53^MUT group. The final IPM model composed three immune-related genes: IPM risk score = (0.392 × S100A8 expression) + (0.174 × CXCL1 expression) + (0.368 × IGLL5 expression), significantly correlated with the overall survival (OS) of GBM in the stratified TP53 status subgroups and was an independent prognostic variate for GBM. By integrating the IPM and clinical characteristics, a nomogram was generated to facilitate clinical utilization, with the results suggesting that it has better predictive performance for GBM prognosis than the IPM.

CONCLUSIONS

The IPM model can identify patients at high-risk and can be combined with other clinical factors to estimate the OS of GBM patients, demonstrating that it is a promising biomarker to optimize the prognosis of GBM.

AWESOME: a database of SNPs that affect protein post-translational modifications

Protein post-translational modifications (PTMs), including phosphorylation, ubiquitination, methylation, acetylation, glycosylation et al, are very important biological processes. PTM changes in some critical genes, which may be induced by base-pair substitution, are shown to affect the risk of diseases. Recently, large-scale exome-wide association studies found that missense single nucleotide polymorphisms (SNPs) play an important role in the susceptibility for complex diseases or traits. One of the functional mechanisms of missense SNPs is that they may affect PTMs and leads to a protein dysfunction and its downstream signaling pathway disorder. Here, we constructed a database named AWESOME (A Website Exhibits SNP On Modification Event, http://www.awesome-hust.com), which is an interactive web-based analysis tool that systematically evaluates the role of SNPs on nearly all kinds of PTMs based on 20 available tools. We also provided a well-designed scoring system to compare the performance of different PTM prediction tools and help users to get a better interpretation of results. Users can search SNPs, genes or position of interest, filter with specific modifications or prediction methods, to get a comprehensive PTM change induced by SNPs. In summary, our database provides a convenient way to detect PTM-related SNPs, which may potentially be pathogenic factors or therapeutic targets.

Functional characterization of a low-frequency V1937I variant in FASN associated with susceptibility to esophageal squamous cell carcinoma

Metabolic reprogramming has been regarded as one of the core hallmarks of cancer and increased de novo fatty acid synthesis has been documented in multiple tumors including esophageal squamous cell carcinoma (ESCC). Our previous exome-wide analyses found a Val1937Ile variant (rs17848945) in the 34th exon of fatty acid synthase (FASN) that showed a strong association with the risk of ESCC. In this study, we performed a series of functional assays to investigate the biological functions underlying this variant in the development of ESCC. We demonstrated that FASN was upregulated in ESCC and both knockdown and knockout of FASN significantly inhibited ESCC cell proliferation, suggesting a tumor promoter role for this gene in ESCC. Furthermore, the results showed that overexpression of FASN[I] in the ESCC cells substantially enhanced cell proliferation, compared with overexpression of FASN[V], or the control vector. Intriguingly, we found that the FASN[I] variant can enhance the enzyme activity of FASN, and, thus, increase the amount of the FASN end-product, palmitate in the ESCC cells. We also observed elevated palmitate levels in the plasma of the FASN[I] genotype carriers among a total of 632 healthy Chinese adults. In conclusion, our results suggested that the FASN V1937I variant influenced ESCC cell proliferation and susceptibility by altering the catabolic activity of FASN on palmitate. These findings may highlight an important role of palmitate metabolism in the development of ESCC and may contribute to the personalized medicine of this disease.

TERT-rs33963617 and CLPTM1L-rs77518573 reduce the risk of non-small cell lung cancer in Chinese population

BACKGROUND

Genome-wide association studies (GWAS) have identified 5p15.33 as a susceptible locus for lung cancer. However, for non-small cell lung cancer (NSCLC), low-frequency risk variants in this region have not been systematically studied. We intended to explore the associations between low-frequency variants on 5p15.33 and NSCLC using a next-generation sequencing based approach in this study.

METHODS

We have acquisited the variation spectrum of 400 NSCLC patients on 5p15.33 by sequencing the targeted region before. Candidate variants were primarily selected by restricting the minor allele frequency (MAF 1-5%) and then by comparing their frequency in 400 NSCLC patients with 1008 East Asians from The genome Aggregation Database (gnomAD). The associations between candidate variants and NSCLC were discovered and replicated in two case-control sets: discovery stage with 960 cases and 916 controls, and replication stage with 1596 cases and 1614 controls in total.

RESULTS

Five low-frequency variants were selected as our candidates and subsequent association analyses showed that 2 polymorphisms were significantly associated with risk of NSCLC, including rs33963617 (OR = 0.63, 95% CI: 0.53-0.76, P = 3.80 × 10^-7) in TERT and rs77518573 (OR = 0.73, 95% CI: 0.63-0.84, P = 2.00 × 10^-5) in upstream of CLPTM1L. When stratified by histologic subtype, a significant association was only investigated in adenocarcinoma for rs77518573. We also observed an obvious cumulative effect of the two significant variants.

CONCLUSIONS

We newly identified two NSCLC related variants on chromosome 5p15.33. Both TERT-rs33963617 and CLPTM1L-rs77518573 conferred reduced risk for NSCLC in Chinese Han population.

Regulating Retinoic Acid Availability during Development and Regeneration: The Role of the CYP26 Enzymes

This review focuses on the role of the Cytochrome p450 subfamily 26 (CYP26) retinoic acid (RA) degrading enzymes during development and regeneration. Cyp26 enzymes, along with retinoic acid synthesising enzymes, are absolutely required for RA homeostasis in these processes by regulating availability of RA for receptor binding and signalling. Cyp26 enzymes are necessary to generate RA gradients and to protect specific tissues from RA signalling. Disruption of RA homeostasis leads to a wide variety of embryonic defects affecting many tissues. Here, the function of CYP26 enzymes is discussed in the context of the RA signalling pathway, enzymatic structure and biochemistry, human genetic disease, and function in development and regeneration as elucidated from animal model studies.

The association of miR34b/c and TP53 gene polymorphisms with Wilms tumor risk in Chinese children

Wilms tumor is the most common pediatric malignancy in the kidney. The miR34b/c is a downstream target gene of the transcription factor p53. The important role of TP53 mutations, the methylation of miR34b/c, and the interaction between these two molecules in tumorigenesis have been well documented. Due to the biological connection between p53 and miR34b/c, in the present study, we investigated the association between polymorphisms in these two molecules and Wilms tumor susceptibility through genotyping two important functional polymorphisms (miR34b/c rs4938723 T>C and TP53 rs1042522 C>G) in 183 cases and 603 controls. Adjusted odds ratios (ORs) and 95% confidence intervals (CIs) derived from the logistic regression analysis were used to assess the correlation of miR34b/c rs4938723 and TP53 rs1042522 polymorphisms with Wilms tumor risk. Our results indicated that the association of miR34b/c rs4938723 and TP53 rs1042522 polymorphisms with Wilms tumor susceptibility was not statistically significant. Stratified analysis by age, gender, and clinical stage, as well as combined effect analysis were also performed, yet, no significant association was found. In conclusion, our study indicated a lack of association between the two selected polymorphisms and Wilms tumor susceptibility. Our findings need to be verified in studies with larger sample size in the future.

NRAS rs2273267 A>T polymorphism reduces neuroblastoma risk in Chinese children

Neuroblastoma is an extracranial solid tumor that mainly occurs in childhood. Mutations of NRAS gene have been described in several cancers. However, whether NRAS gene polymorphisms can predict the risk of neuroblastoma have not been investigated. We hypothesized that variations of NRAS gene contribute to neuroblastoma predisposition. Therefore, we conducted a multi-center case-control study using 263 cases and 715 controls to examine the association of NRAS gene rs2273267 A>T polymorphism and neuroblastoma risk. We calculated odds ratios (ORs) and corresponding 95% confidence intervals (CIs) to assess the strength of the associations. Relative to those with AA genotype, subjects with AT/TT genotype had reduced neuroblastoma risk (adjusted OR = 0.72, 95% CI = 0.54-0.96, P = 0.024). Stratified analysis revealed that rs2273267 AT/TT carriers were less likely to develop neuroblastoma for patients with tumor originating from the adrenal gland (adjusted OR = 0.67, 95% CI = 0.45-0.99, P = 0.047) and clinical stages III + IV (adjusted OR = 0.57, 95% CI = 0.36-0.90, P = 0.015). Our findings underline the likely importance of NRAS gene rs2273267 A>T in the risk of neuroblastoma. Further independent case-control studies with functional analysis are needed to verify the role of NRAS gene rs2273267 A>T polymorphism in the risk of neuroblastoma.

A functional variant in TNXB promoter associates with the risk of esophageal squamous-cell carcinoma

Our previous study identified a tag single-nucleotide polymorphism (SNP) rs204900 in TNXB associated with risk of esophageal squamous-cell carcinoma (ESCC) in the Chinese population. However, the functional role of TNXB and causal variants had not been interrogated in that study. In the present study, we explored the effects of TNXB expression in the development of ESCC and searched for functional variants in this gene. We found TNXB was downregulated in ESCC tumors. Using small interfering RNAs and CRISPR-Cas9 methods, we identified that both knockdown and knockout of TNXB significantly promoted ESCC cell growth in vitro, suggesting a tumor suppressor role of this gene in ESCC. Through further fine-mapping analysis, we identified that a noncoding variant in the promoter of TNXB, rs411337, predisposed to ESCC risk (odds ratio = 1.36, 95% confidence interval: 1.22-1.51, P = 9.10 × 10^-9 ). These findings revealed the functional mechanism of TNXB in the development of ESCC and may contribute to the prevention and treatment of this disease in the future.

GWAS in cancer: progress and challenges

The genome-wide association study (GWAS) is an effective method to detect single-nucleotide polymorphisms (SNPs) of multiple individual genes based on linkage disequilibrium (LD). GWAS examines genotypes and distinguishing gene characteristics that are exhibited in diseases. In the past few decades, more and more literature has reported the results of applying GWAS to study tumors. Although many pleiotropic loci associated with complex phenotypes have been identified by GWAS, the biological functions of many genetic variation loci remain unclear, and the genetic mechanisms of most complex phenotypes cannot be systematically explained. In this article, we will review the new findings of several tumor types, and categorize the new sites and mechanisms that have recently been discovered. We linked the mechanisms of action of various tumors and searched for links to related gene expression pathways. We found that susceptible sites can be divided into hub genes and peripheral genes; the two interact to link gene expression in a variety of diseases.

Inter- and intratumor DNA methylation heterogeneity associated with lymph node metastasis and prognosis of esophageal squamous cell carcinoma

Background: Esophageal squamous cell carcinoma (ESCC), one of the leading causes of cancer mortality worldwide, is a heterogeneous cancer with diverse clinical manifestations. However, little is known about the epigenetic heterogeneity and its clinical relevance for this prevalent cancer. Methods: We generated 7.56 Tb single-base resolution whole-genome bisulfite sequencing data for 84 ESCC and paired paraneoplastic tissues. The analysis identified inter- and intratumor DNA methylation (DNAm) heterogeneity, epigenome-wide DNAm alterations together with the functional regulators involved in the hyper- or hypomethylated regions, and their association with clinical features. We then validated the correlation between the methylation level of specific regions and clinical outcomes of 96 ESCC patients in an independent cohort. Results: ESCC manifested substantial inter- and intratumor DNAm heterogeneity. The high intratumor DNAm heterogeneity was associated with lymph node metastasis and worse overall survival. Interestingly, hypermethylated regions in ESCC were enriched in promoters of numerous transcription factors, and demethylated noncoding regions related to RXR transcription factor binding appeared to contribute to the development of ESCC. Furthermore, we identified numerous DNAm alterations associated with carcinogenesis and lymph node metastasis of ESCC. We also validated three novel prognostic markers for ESCC, including one each in the promoter of CLK1, the 3' untranslated region of ZEB2, and the intergenic locus surrounded by several lncRNAs. Conclusions: This study presents the first population-level resource for dissecting base-resolution DNAm variation in ESCC and provides novel insights into the ESCC pathogenesis and progression, which might facilitate diagnosis and prognosis for this prevalent malignancy.

Association of CMYC polymorphisms with hepatoblastoma risk

BACKGROUND

Single-nucleotide polymorphisms (SNPs) in genes may affect gene expression and contribute to cancer susceptibility. This study aimed to explore the association between CMYC gene polymorphisms and hepatoblastoma risk.

METHODS

Hepatoblastoma patients and cancer-free controls were recruited and matched by age and sex. Genotypes were determined by TaqMan, and the strength of the association of interest was determined by calculating odds ratios (ORs) and 95% confidence intervals (CIs). The distributions of various CMYC genotypes among subjects were recorded, followed by analyses of associations between CMYC polymorphisms and hepatoblastoma risk.

RESULTS

A total of 213 hepatoblastoma patients and 958 cancer-free controls were enrolled. No significant associations between the CMYC rs4645943 and rs2070583 polymorphisms and hepatoblastoma risk were found (all P>0.05). In stratification analysis based on age, sex, and clinical stage, the CMYC rs4645943 and rs2070583 polymorphisms were not associated with hepatoblastoma susceptibility (all P>0.05).

CONCLUSIONS

Thus, the CMYC rs4645943 and rs2070583 polymorphisms were not associated with hepatoblastoma risk in the study cohort.

Genetic Predisposition to Colon and Rectal Adenocarcinoma Is Mediated by a Super-enhancer Polymorphism Coactivating CD9 and PLEKHG6

BACKGROUND

Genome-wide association studies (GWAS) have identified dozens of loci associated with colon and rectal adenocarcinoma risk. As tissue-specific super-enhancers (SE) play important roles in tumorigenesis, we systematically investigate SEs and inner variants in established GWAS loci to decipher the underlying biological mechanisms.

METHODS

Through a comprehensive bioinformatics analysis on multi-omics data, we screen potential single-nucleotide polymorphisms (SNP) in cancer-specific SEs, and then subject them to a two-stage case-control study containing 4,929 cases and 7,083 controls from the Chinese population. A series of functional assays, including reporter gene assays, electrophoretic mobility shift assays (EMSA), CRISPR-Cas9 genome editing, chromosome conformation capture (3C) assays, and cell proliferation experiments, are performed to characterize the variant's molecular consequence and target genes.

RESULTS

The SNP rs11064124 in 12p13.31 is found significantly associated with the risk of colon and rectal adenocarcinoma with an odds ratio (OR) of 0.87 [95% confidence interval (CI), 0.82-0.92, P = 8.67E-06]. The protective rs11064124-G weakens the binding affinity with vitamin D receptor (VDR) and increases the enhancer's activity and interactions with two target genes' promoters, thus coactivating the transcription of CD9 and PLEKHG6, which are both putative tumor suppressor genes for colon and rectal adenocarcinoma.

CONCLUSIONS

Our integrative study highlights an SE polymorphism rs11064124 and two susceptibility genes CD9 and PLEKHG6 in 12p13.31 for colon and rectal adenocarcinoma.

IMPACT

These findings suggest a novel insight for genetic pathogenesis of colon and rectal adenocarcinoma, involving transcriptional coactivation of diverse susceptibility genes via the SE element as a gene regulation hub.

Cumulative Evidence for Associations between Genetic Variants and Risk of Esophageal Cancer

BACKGROUND

A large number of studies have been conducted to investigate associations between genetic variants and esophageal cancer risk in the past several decades. However, findings from these studies have been generally inconsistent. We aimed to provide a summary of the current understanding of the genetic architecture of esophageal cancer susceptibility.

METHODS

We performed a comprehensive field synopsis and meta-analysis to evaluate associations between 95 variants in 70 genes or loci and esophageal cancer risk using data from 304 eligible publications, including 104,904 cases and 159,797 controls, through screening a total of 21,328 citations. We graded levels of cumulative epidemiologic evidence of a significant association with esophageal cancer using the Venice criteria and false-positive report probability tests. We constructed functional annotations for these variants using data from the Encyclopedia of DNA Elements Project and other databases.

RESULTS

Thirty variants were nominally significantly associated with esophageal cancer risk. Cumulative epidemiologic evidence of a significant association with overall esophageal cancer, esophageal squamous cell carcinoma, or esophageal adenocarcinoma was strong for 13 variants in or near 13 genes (ADH1B, BARX1, CDKN1A, CHEK2, CLPTM1L, CRTC1, CYP1A1, EGF, LTA, MIR34BC, PLCE1, PTEN, and PTGS2). Bioinformatics analysis suggested that these variants and others correlated with them might fall in putative functional regions.

CONCLUSIONS

Our study summarizes the current literature on the genetic architecture of esophageal cancer susceptibility and identifies several potential polymorphisms that could be involved in esophageal cancer susceptibility.

IMPACT

These findings provide direction for future studies to identify new genetic factors for esophageal cancer.

ALDH3B2 Polymorphism Is Associated with Colorectal Cancer Susceptibility

Colorectal cancer (CRC) is the 5^th leading cancer in China. Alcohol consumption has been reported to be one of the risk factors of CRC. However, it remains unclear whether genetic variants of alcohol metabolic genes are associated with CRC risk. In this study, we tested the coding variants in the alcohol metabolic genes and the risk of CRC, by using 485 cases and 516 controls. A total of 16 germline coding variants in 10 alcohol metabolic genes were genotyped. We identified that rs3741178 in ALDH3B2 was significantly associated with CRC risk with odds ratio being 2.13 (95% CI: 1.24–3.68, P=0.0064). Further functional annotation suggested that this variant may damage the protein function of ALDH3B2. Our results suggested that ALDH3B2 in the alcohol metabolism pathway contributed to the development of CRC, which may contribute to the prevention of this disease in the future.

LIN28A gene polymorphisms modify neuroblastoma susceptibility: A four-centre case-control study

Neuroblastoma ranks the most common seen solid tumour in childhood. Overexpression of LIN28A gene has been linked to the development of multiple human malignancies, but the relationship between LIN28A single nucleotide polymorphisms (SNPs) and neuroblastoma susceptibility is still under debate. Herein, we evaluated the correlation of four potentially functional LIN28A SNPs (rs3811464 G>A, rs3811463 T>C, rs34787247 G>A, and rs11247957 G>A) and neuroblastoma susceptibility in 505 neuroblastoma patients and 1070 controls from four independent hospitals in China. The correlation strengths were determined by using odds ratios (ORs) and corresponding 95% confidence intervals (CIs). Among these SNPs, rs34787247 G>A exhibited a significant association with increased susceptibility in neuroblastoma (GA vs GG: adjusted OR = 1.30, 95% CI = 1.03-1.64; AA vs GG: adjusted OR = 2.51, 95% CI = 1.36-4.64, AA/GA vs GG: adjusted OR = 1.42, 95% CI = 1.12-1.80, AA vs GG/GA: adjusted OR = 2.39, 95% CI = 1.29-4.42). Furthermore, the combined analysis of risk genotypes revealed that subjects carrying three risk genotypes (adjusted OR = 1.64, 95% CI = 1.02-2.63) are more inclined to develop neuroblastoma than those without risk genotype, and so do carriers of 1-4 risk genotypes (adjusted OR = 1.26, 95% CI = 1.01-1.56). Stratification analysis further revealed risk effect of rs3811464 G>A, rs34787247 G>A and 1-4 risk genotypes in some subgroups. Haplotype analysis of these four SNPs yields two haplotypes significantly correlated with increased neuroblastoma susceptibility. Overall, our finding indicated that LIN28A SNPs, especially rs34787247 G>A, may increase neuroblastoma risk.

HMGA2 Gene rs8756 A>C Polymorphism Reduces Neuroblastoma Risk in Chinese Children: A Four-Center Case-Control Study

BACKGROUND

Neuroblastoma, mainly affecting children, is a lethal malignancy arising from the developing sympathetic nervous system. The genetic etiology of neuroblastoma remains mostly obscure. High mobility group AT-hook 2 (HMGA2), an oncogenic gene, is up-regulated in many tumors. Single nucleotide polymorphisms (SNPs) often modify cancer susceptibility. However, no studies are investigating the association between HMGA2 SNPs and neuroblastoma susceptibility.

METHODS

We conducted a four-center case-control study to evaluate the association between three HMGA2 polymorphisms (rs6581658 A>G, rs8756 A>C and rs968697 T>C) and neuroblastoma susceptibility in a Chinese population with 505 cases and 1070 controls. Logistic regression was performed to evaluate the strength of the association.

RESULTS

We found that the rs8756 AC/CC genotypes were associated with a reduced neuroblastoma risk when compared to rs8756 AA genotype [Adjusted odds ratio (OR)=0.74, 95% confidence interval (CI)=0.56-0.99, P=0.039]. Carriers with 3 protective genotypes have lower neuroblastoma susceptibility than those without or with 0-2 protective genotypes. The stratified analysis revealed that the protective effects of rs8756 AC/CC genotypes were more predominant among children of age > 18 months, males, and subgroups with the tumor in the mediastinum. Furthermore, haplotype analysis uncovered that haplotype ACC significantly reduced neuroblastoma risk.

CONCLUSION

Our study indicated HMGA2 rs8756 A>C polymorphism is significantly associated with decreased neuroblastoma risk.

Evaluation of polymorphisms in microRNA-binding sites and pancreatic cancer risk in Chinese population

As promising biomarkers and therapy targets, microRNAs (miRNAs) are involved in various physiological and tumorigenic processes. Genetic variants in miRNA‐binding sites can lead to dysfunction of miRNAs and contribute to disease. However, systematic investigation of the miRNA‐related single nucleotide polymorphisms (SNPs) for pancreatic cancer (PC) risk remains elusive. We performed integrative bioinformatics analyses to select 31 SNPs located in miRNA‐target binding sites using the miRNASNP v2.0, a solid database providing miRNA‐related SNPs for genetic research, and investigated their associations with risk of PC in two large case‐control studies totally including 1847 cases and 5713 controls. We observed that the SNP rs3802266 is significantly associated with increased risk of PC (odds ratio (OR) = 1.21, 95% confidence intervals (CI) = 1.11‐1.31, P = 1.29E‐05). Following luciferase reporter gene assays show that rs3802266‐G creates a stronger binding site for miR‐181a‐2‐3p in 3′ untranslated region (3′UTR) of the gene ZHX2. Expression quantitative trait loci (eQTL) analysis suggests that ZHX2 expression is lower in individuals carrying rs3802266‐G with increased PC risk. In conclusion, our findings highlight the involvement of miRNA‐binding SNPs in PC susceptibility and provide new clues for PC carcinogenesis.

Biochemical and physiological importance of the CYP26 retinoic acid hydroxylases

The Cytochrome P450 (CYP) family 26 enzymes contribute to retinoic acid (RA) metabolism and homeostasis in humans, mammals and other chordates. The three CYP26 family enzymes, CYP26A1, CYP26B1 and CYP26C1 have all been shown to metabolize all-trans-retinoic acid (atRA) it's 9-cisRA and 13-cisRA isomers and primary metabolites 4-OH-RA and 4-oxo-RA with high efficiency. While no crystal structures of CYP26 enzymes are available, the binding of various ligands has been extensively explored via homology modeling. All three CYP26 enzymes are inducible by treatment with atRA in various prenatal and postnatal tissues and cell types. However, current literature shows that in addition to regulation by atRA, CYP26 enzyme expression is also regulated by other endogenous processes and inflammatory cytokines. In humans and in animal models the expression patterns of CYP26 enzymes have been shown to be tissue and cell type specific, and the expression of the CYP26 enzymes is believed to regulate the formation of critical atRA concentration gradients in various tissue types. Yet, very little data exists on direct disease associations of altered CYP26 expression or activity. Nevertheless, data is emerging describing a variety of human genetic variations in the CYP26 enzymes that are associated with different pathologies. Interestingly, some of these genetic variants result in increased activity of the CYP26 enzymes potentially leading to complex gene-environment interactions due to variability in dietary intake of retinoids. This review highlights the current knowledge of structure-function of CYP26 enzymes and focuses on their role in human retinoid metabolism in different tissues.

Exome-chip association analysis of intracranial aneurysms

OBJECTIVE

To investigate to what extent low-frequency genetic variants (with minor allele frequencies <5%) affect the risk of intracranial aneurysms (IAs).

METHODS

One thousand fifty-six patients with IA and 2,097 population-based controls from the Netherlands were genotyped with the Illumina HumanExome BeadChip. After quality control (QC) of samples and single nucleotide variants (SNVs), we conducted a single variant analysis using the Fisher exact test. We also performed the variable threshold (VT) test and the sequence kernel association test (SKAT) at different minor allele count (MAC) thresholds of >5 and >0 to test the hypothesis that multiple variants within the same gene are associated with IA risk. Significant results were tested in a replication cohort of 425 patients with IA and 311 controls, and results of the 2 cohorts were combined in a meta-analysis.

RESULTS

After QC, 995 patients with IA and 2,080 controls remained for further analysis. The single variant analysis comprising 46,534 SNVs did not identify significant loci at the genome-wide level. The gene-based tests showed a statistically significant association for fibulin 2 (FBLN2) (best p = 1 × 10^-6 for the VT test, MAC >5). Associations were not statistically significant in the independent but smaller replication cohort (p > 0.57) but became slightly stronger in a meta-analysis of the 2 cohorts (best p = 4.8 × 10^-7 for the SKAT, MAC ≥1).

CONCLUSION

Gene-based tests indicated an association for FBLN2, a gene encoding an extracellular matrix protein implicated in vascular wall remodeling, but independent validation in larger cohorts is warranted. We did not identify any significant associations for single low-frequency genetic variants.