Integrating expression-related SNPs into genome-wide gene- and pathway-based analyses identified novel lung cancer susceptibility genes

Identification of genetically predicted DNA methylation markers associated with non-small cell lung cancer risk among 34,964 cases and 448,579 controls

BACKGROUND

Although the associations between genetic variations and lung cancer risk have been explored, the epigenetic consequences of DNA methylation in lung cancer development are largely unknown. Here, the genetically predicted DNA methylation markers associated with non-small cell lung cancer (NSCLC) risk by a two-stage case-control design were investigated.

METHODS

The genetic prediction models for methylation levels based on genetic and methylation data of 1595 subjects from the Framingham Heart Study were established. The prediction models were applied to a fixed-effect meta-analysis of screening data sets with 27,120 NSCLC cases and 27,355 controls to identify the methylation markers, which were then replicated in independent data sets with 7844 lung cancer cases and 421,224 controls. Also performed was a multi-omics functional annotation for the identified CpGs by integrating genomics, epigenomics, and transcriptomics and investigation of the potential regulation pathways.

RESULTS

Of the 29,894 CpG sites passing the quality control, 39 CpGs associated with NSCLC risk (Bonferroni-corrected p ≤ 1.67 × 10-6 ) were originally identified. Of these, 16 CpGs remained significant in the validation stage (Bonferroni-corrected p ≤ 1.28 × 10-3 ), including four novel CpGs. Multi-omics functional annotation showed nine of 16 CpGs were potentially functional biomarkers for NSCLC risk. Thirty-five genes within a 1-Mb window of 12 CpGs that might be involved in regulatory pathways of NSCLC risk were identified.

CONCLUSIONS

Sixteen promising DNA methylation markers associated with NSCLC were identified. Changes of the methylation level at these CpGs might influence the development of NSCLC by regulating the expression of genes nearby.

PLAIN LANGUAGE SUMMARY

The epigenetic consequences of DNA methylation in lung cancer development are still largely unknown. This study used summary data of large-scale genome-wide association studies to investigate the associations between genetically predicted levels of methylation biomarkers and non-small cell lung cancer risk at the first time. This study looked at how well larotrectinib worked in adult patients with sarcomas caused by TRK fusion proteins. These findings will provide a unique insight into the epigenetic susceptibility mechanisms of lung cancer.

Regulatory network and targeted interventions for CCDC family in tumor pathogenesis

CCDC (coiled-coil domain-containing) is a coiled helix domain that exists in natural proteins. There are about 180 CCDC family genes, encoding proteins that are involved in intercellular transmembrane signal transduction and genetic signal transcription, among other functions. Alterations in expression, mutation, and DNA promoter methylation of CCDC family genes have been shown to be associated with the pathogenesis of many diseases, including primary ciliary dyskinesia, infertility, and tumors. In recent studies, CCDC family genes have been found to be involved in regulation of growth, invasion, metastasis, chemosensitivity, and other biological behaviors of malignant tumor cells in various cancer types, including nasopharyngeal carcinoma, lung cancer, colorectal cancer, and thyroid cancer. In this review, we summarize the involvement of CCDC family genes in tumor pathogenesis and the relevant upstream and downstream molecular mechanisms. In addition, we summarize the potential of CCDC family genes as tumor therapy targets. The findings discussed here help us to further understand the role and the therapeutic applications of CCDC family genes in tumors.

Potential Involvement of NSD1, KRT24 and ACACA in the Genetic Predisposition to Colorectal Cancer

Simple Summary

Methods used for the identification of hereditary cancer genes have evolved in parallel to technological progress; however, much of the genetic predisposition to cancer remains unexplained. A new in silico method based on Knudson’s two-hit hypothesis recently identified ~50 putative cancer predisposing genes, but their actual association with cancer has not yet been validated. In our study, we aimed to assess the involvement of these genes in familial/early-onset colorectal cancer (CRC) using different lines of evidence. Our results indicated that most of those genes were not associated with a genetic predisposition to CRC, but suggested a possible association for NSD1, KRT24 and ACACA.

Abstract

The ALFRED (Allelic Loss Featuring Rare Damaging) in silico method was developed to identify cancer predisposition genes through the identification of somatic second hits. By applying ALFRED to ~10,000 tumor exomes, 49 candidate genes were identified. We aimed to assess the causal association of the identified genes with colorectal cancer (CRC) predisposition. Of the 49 genes, NSD1, HDAC10, KRT24, ACACA and TP63 were selected based on specific criteria relevant for hereditary CRC genes. Gene sequencing was performed in 736 patients with familial/early onset CRC or polyposis without germline pathogenic variants in known genes. Twelve (predicted) damaging variants in 18 patients were identified. A gene-based burden test in 1596 familial/early-onset CRC patients, 271 polyposis patients, 543 TCGA CRC patients and >134,000 controls (gnomAD, non-cancer), revealed no clear association with CRC for any of the studied genes. Nevertheless, (non-significant) over-representation of disruptive variants in NSD1, KRT24 and ACACA in CRC patients compared to controls was observed. A somatic second hit was identified in one of 20 tumors tested, corresponding to an NSD1 carrier. In conclusion, most genes identified through the ALFRED in silico method were not relevant for CRC predisposition, although a possible association was detected for NSD1, KRT24 and ACACA.

Gene signature to predict prognostic survival of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) has a high incidence and poor prognosis and is the second most fatal cancer, and certain HCC patients also show high heterogeneity. This study developed a prognostic model for predicting clinical outcomes of HCC. RNA and microRNA (miRNA) sequencing data of HCC were obtained from the cancer genome atlas. RNA dysregulation between HCC tumors and adjacent normal liver tissues was examined by DESeq algorithms. Survival analysis was conducted to determine the basic prognostic indicators. We identified competing endogenous RNA (ceRNA) containing 15,364 pairs of mRNA–long noncoding RNA (lncRNA). An imbalanced ceRNA network comprising 8 miRNAs, 434 mRNAs, and 81 lncRNAs was developed using hypergeometric test. Functional analysis showed that these RNAs were closely associated with biosynthesis. Notably, 53 mRNAs showed a significant prognostic correlation. The least absolute shrinkage and selection operator’s feature selection detected four characteristic genes (SAPCD2, DKC1, CHRNA5, and UROD), based on which a four-gene independent prognostic signature for HCC was constructed using Cox regression analysis. The four-gene signature could stratify samples in the training, test, and external validation sets (p <0.01). Five-year survival area under ROC curve (AUC) in the training and validation sets was greater than 0.74. The current prognostic gene model exhibited a high stability and accuracy in predicting the overall survival (OS) of HCC patients.

Mitochondrially localized MPZL3 emerges as a signaling hub of mammalian physiology

MPZL3 is a nuclear-encoded, mitochondrially localized, immunoglobulin-like V-type protein that functions as a key regulator of epithelial cell differentiation, lipid metabolism, ROS production, glycemic control, and energy expenditure. Recently, MPZL3 has surfaced as an important modulator of sebaceous gland function and of hair follicle cycling, an organ transformation process that is also governed by peripheral clock gene activity and PPARγ. Given the phenotype similarities and differences between Mpzl3 and Pparγ knockout mice, we propose that MPZL3 serves as a signaling hub that is regulated by core clock gene products and/or PPARγ to translate signals from these nuclear transcription factors to the mitochondria to modulate circadian and metabolic regulation. Conservation between murine and human MPZL3 suggests that human MPZL3 may have similarly complex functions in health and disease. We summarize current knowledge and discuss future directions to elucidate the full spectrum of MPZL3 functions in mammalian physiology.

A cross-tissue transcriptome-wide association study identifies novel susceptibility genes for lung cancer in Chinese populations

Although dozens of susceptibility loci have been identified for lung cancer in genome-wide association studies (GWASs), the susceptibility genes and underlying mechanisms remain unclear. In this study, we conducted a cross-tissue transcriptome-wide association study (TWAS) with UTMOST based on summary statistics from 13 327 lung cancer cases and 13 328 controls and the genetic-expression matrix over 44 human tissues in the Genotype-Tissue Expression (GTEx) project. After further evaluating the associations in each tissue, we revealed 6 susceptibility genes in known loci and identified 12 novel ones. Among those, five novel genes, including DCAF16 (Pcross-tissue = 2.57 × 10-5, PLung = 2.89 × 10-5), CBL (Pcross-tissue = 5.08 × 10-7, PLung = 1.82 × 10-4), ATR (Pcross-tissue = 1.45 × 10-5, PLung = 9.68 × 10-5), GYPE (Pcross-tissue = 1.45 × 10-5, PLung = 2.17 × 10-3) and PARD3 (Pcross-tissue = 5.79 × 10-6, PLung = 4.05 × 10-3), were significantly associated with the risk of lung cancer in both cross-tissue and lung tissue models. Further colocalization analysis indicated that rs7667864 (C > A) and rs2298650 (G > T) drove the GWAS association signals at 4p15.31-32 (OR = 1.09, 95%CI: 1.04-1.12, PGWAS = 5.54 × 10-5) and 11q23.3 (OR = 1.08, 95%CI: 1.04-1.13, PGWAS = 5.55 × 10-5), as well as the expression of DCAF16 (βGTEx = 0.24, PGTEx = 9.81 × 10-15; βNJLCC = 0.29, PNJLCC = 3.84 × 10-8) and CBL (βGTEx = -0.17, PGTEx = 2.82 × 10-8; βNJLCC = -0.32, PNJLCC = 2.61 × 10-7) in lung tissue. Functional annotations and phenotype assays supported the carcinogenic effect of these novel susceptibility genes in lung carcinogenesis.

Benchmark of site- and structure-specific quantitative tissue N-glycoproteomics for discovery of potential N-glycoprotein markers: a case study of pancreatic cancer

Pancreatic cancer is a highly malignant tumor of the digestive tract that is difficult to diagnose and treat. It is more common in developed countries and has become one of the main causes of death in some countries and regions. Currently, pancreatic cancer generally has a poor prognosis, partly due to the lack of symptoms in the early stages of pancreatic cancer. Therefore, most cases are diagnosed at advanced stage. With the continuous in-depth research of glycoproteomics in precision medical diagnosis, there have been some reports on quantitative analysis of cancer-related cells, plasma or tissues to find specific biomarkers for targeted therapy. This research is based on the developed complete N-linked glycopeptide database search engine GPSeeker, combined with liquid-mass spectrometry and stable diethyl isotope labeling, providing a benchmark of site- and structure-specific quantitative tissue N-glycoproteomics for discovery of potential N-glycoprotein markers. With spectrum-level FDR ≤1%, 20,038 intact N-Glycopeptides corresponding to 4518 peptide backbones, 228 N-glycan monosaccharide compositions 1026 N-glycan putative structures, 4460 N-glycosites and 3437 intact N-glycoproteins were identified. With the criteria of ≥1.5-fold change and p value<0.05, 52 differentially expressed intact N-glycopeptides (DEGPs) were found in pancreatic cancer tussues relative to control, where 38 up-regulated and 14 down-regulated, respectively.

Integration of GWAS and eQTL Analysis to Identify Risk Loci and Susceptibility Genes for Gastric Cancer

Genome-wide association studies (GWAS) have identified several susceptibility loci for gastric cancer (GC), but the majority of identified single-nucleotide polymorphisms (SNPs) fall within the non-coding region and are likely to exert their biological function by modulating gene expression. To systematically estimate expression-associated SNPs (eSNPs) that confer genetic predisposition to GC, we evaluated the associations of 314,203 stomach tissue-specific eSNPs with GC risk in three GWAS datasets (2,631 cases and 4,373 controls). Subsequently, we conducted a gene-based analysis to calculate the cumulative effect of eSNPs through sequence kernel association combined test and Sherlock integrative analysis. At the SNP-level, we identified two novel variants (rs836545 at 7p22.1 and rs1892252 at 6p22.2) associated with GC risk. The risk allele carriers of rs836545-T and rs1892252-G exhibited higher expression levels of DAGLB (P = 3.70 × 10^–18) and BTN3A2 (P = 3.20 × 10^–5), respectively. Gene-based analyses identified DAGLB and FBXO43 as novel susceptibility genes for GC. DAGLB and FBXO43 were significantly overexpressed in GC tissues than in their adjacent tissues (P = 5.59 × 10^–7 and P = 3.90 × 10^–6, respectively), and high expression level of these two genes was associated with an unfavorable prognosis of GC patients (P = 1.30 × 10^–7 and P = 7.60 × 10^–3, respectively). Co-expression genes with these two novel genes in normal stomach tissues were significantly enriched in several cancer-related pathways, including P53, MAPK and TGF-beta pathways. In summary, our findings confirm the importance of eSNPs in dissecting the genetic basis of GC, and the identified eSNPs and relevant genes will provide new insight into the genetic and biological basis for the mechanism of GC development.

Association Analysis of Driver Gene-Related Genetic Variants Identified Novel Lung Cancer Susceptibility Loci with 20,871 Lung Cancer Cases and 15,971 Controls

BACKGROUND

A substantial proportion of cancer driver genes (CDG) are also cancer predisposition genes. However, the associations between genetic variants in lung CDGs and the susceptibility to lung cancer have rarely been investigated.

METHODS

We selected expression-related single-nucleotide polymorphisms (eSNP) and nonsynonymous variants of lung CDGs, and tested their associations with lung cancer risk in two large-scale genome-wide association studies (20,871 cases and 15,971 controls of European descent). Conditional and joint association analysis was performed to identify independent risk variants. The associations of independent risk variants with somatic alterations in lung CDGs or recurrently altered pathways were investigated using data from The Cancer Genome Atlas (TCGA) project.

RESULTS

We identified seven independent SNPs in five lung CDGs that were consistently associated with lung cancer risk in discovery (P < 0.001) and validation (P < 0.05) stages. Among these loci, rs78062588 in TPM3 (1q21.3) was a new lung cancer susceptibility locus (OR = 0.86, P = 1.65 × 10-6). Subgroup analysis by histologic types further identified nine lung CDGs. Analysis of somatic alterations found that in lung adenocarcinomas, rs78062588[C] allele (TPM3 in 1q21.3) was associated with elevated somatic copy number of TPM3 (OR = 1.16, P = 0.02). In lung adenocarcinomas, rs1611182 (HLA-A in 6p22.1) was associated with truncation mutations of the transcriptional misregulation in cancer pathway (OR = 0.66, P = 1.76 × 10-3).

CONCLUSIONS

Genetic variants can regulate functions of lung CDGs and influence lung cancer susceptibility.

IMPACT

Our findings might help unravel biological mechanisms underlying lung cancer susceptibility.

Silencing DVL3 defeats MTX resistance and attenuates stemness via Notch Signaling Pathway in colorectal cancer

Chemoresistance and recurrence of colorectal cancer are mainly caused by the existence of cancer stem-like cells. Dishevelled-3 (DVL3) is a common member of both Wnt/β-catenin pathway and the Notch signaling pathway, which were involved in cancer progression, chemoresistance and even maintenance of stem cell-like properties. However, the underlying biological function of DVL3 still remains unclear. In this study, we proposed DVL3 was simultaneously involved in Methotrexate (MTX) resistance and Colorectal cancer (CRC) stemness by bioinformatic analysis. We also demonstrated DVL3 knockdown sensitized CRC cells to MTX and inhibited their stem cell-like properties by functional experiments. As for the mechanism, DVL3 silencing attenuated the activated Notch signaling by down-regulating Notch intracellular domain (NICD) as well as its downstream targets. Additionally, we demonstrated that CRC cancer tissues had greater DVL3 expression than adjacent non-cancer tissues and patients' overall survival was closely associated with DVL3 according to the data in our clinical center. Accordingly, our data suggest that DVL3 is a key regulator in CRC stemness and chemoresistance and targeting DVL3 could be a potential strategy for CRC therapy.

Systematic analyses of regulatory variants in DNase I hypersensitive sites identified two novel lung cancer susceptibility loci

DNase I hypersensitive sites (DHS) are abundant in regulatory elements, such as promoter, enhancer and transcription factor binding sites. Many studies have revealed that disease-associated variants were concentrated in DHS-related regions. However, limited studies are available on the roles of DHS-related variants in lung cancer. In this study, we performed a large-scale case-control study with 20 871 lung cancer cases and 15 971 controls to evaluate the associations between regulatory genetic variants in DHS and lung cancer susceptibility. The expression quantitative trait loci (eQTL) analysis and pathway-enrichment analysis were performed to identify the possible target genes and pathways. In addition, we performed motif-based analysis to explore the lung-cancer-related motifs using sequence kernel association test. Two novel variants, rs186332 in 20q13.3 (C>T, odds ratio [OR] = 1.17, 95% confidence interval [95% CI]: 1.10-1.24, P = 8.45 × 10-7) and rs4839323 in 1p13.2 (T>C, OR = 0.92, 95% CI: 0.89-0.95, P = 1.02 × 10-6) showed significant association with lung cancer risk. The eQTL analysis suggested that these two SNPs might regulate the expression of MRGBP and SLC16A1, respectively. What's more, the expression of both MRGBP and SLC16A1 was aberrantly elevated in lung tumor tissues. The motif-based analysis identified 10 motifs related to the risk of lung cancer (P < 1.71 × 10-4). Our findings suggested that variants in DHS might modify lung cancer susceptibility through regulating the expression of surrounding genes. This study provided us a deeper insight into the roles of DHS-related genetic variants for lung cancer.

Identification of an early diagnostic biomarker of lung adenocarcinoma based on co-expression similarity and construction of a diagnostic model

Background

The purpose of this study was to achieve early and accurate diagnosis of lung cancer and long-term monitoring of the therapeutic response.

Methods

We downloaded GSE20189 from GEO database as analysis data. We also downloaded human lung adenocarcinoma RNA-seq transcriptome expression data from the TCGA database as validation data. Finally, the expression of all of the genes underwent z test normalization. We used ANOVA to identify differentially expressed genes specific to each stage, as well as the intersection between them. Two methods, correlation analysis and co-expression network analysis, were used to compare the expression patterns and topological properties of each stage. Using the functional quantification algorithm, we evaluated the functional level of each significantly enriched biological function under different stages. A machine-learning algorithm was used to screen out significant functions as features and to establish an early diagnosis model. Finally, survival analysis was used to verify the correlation between the outcome and the biomarkers that we found.

Results

We screened 12 significant biomarkers that could distinguish lung cancer patients with diverse risks. Patients carrying variations in these 12 genes also presented a poor outcome in terms of survival status compared with patients without variations.

Conclusions

We propose a new molecular-based noninvasive detection method. According to the expression of the stage-specific gene set in the peripheral blood of patients with lung cancer, the difference in the functional level is quantified to realize the early diagnosis and prediction of lung cancer.