Genome-wide association meta-analysis identifies pleiotropic risk loci for aerodigestive squamous cell cancers

Identifying therapeutic targets for cancer among 2074 circulating proteins and risk of nine cancers

Circulating proteins can reveal key pathways to cancer and identify therapeutic targets for cancer prevention. We investigate 2,074 circulating proteins and risk of nine common cancers (bladder, breast, endometrium, head and neck, lung, ovary, pancreas, kidney, and malignant non-melanoma) using cis protein Mendelian randomisation and colocalization. We conduct additional analyses to identify adverse side-effects of altering risk proteins and map cancer risk proteins to drug targets. Here we find 40 proteins associated with common cancers, such as PLAUR and risk of breast cancer [odds ratio per standard deviation increment: 2.27, 1.88-2.74], and with high-mortality cancers, such as CTRB1 and pancreatic cancer [0.79, 0.73-0.85]. We also identify potential adverse effects of protein-altering interventions to reduce cancer risk, such as hypertension. Additionally, we report 18 proteins associated with cancer risk that map to existing drugs and 15 that are not currently under clinical investigation. In sum, we identify protein-cancer links that improve our understanding of cancer aetiology. We also demonstrate that the wider consequence of any protein-altering intervention on well-being and morbidity is required to interpret any utility of proteins as potential future targets for therapeutic prevention.

Low expression of PINK1 and PARK2 predicts poor prognosis in patients with esophageal squamous cell carcinoma

BACKGROUND

The Parkinson's disease (PD) gene family expression is strongly linked to tumor development and progression; PINK1 and PARK2 are essential members of the PD gene family. However, the relationship between PINK1 and PARK2 and esophageal squamous cell carcinoma (ESCC) remains unknown. This research aims to clarify the prognostic value of PINK1 and PARK2 in ESCC.

METHODS

PINK1 and PARK2 protein levels in 232 ESCC specimens, and 125 matched adjacent normal tissues were detected by immunohistochemistry. The relationship between PINK1 and PARK2 protein expression and clinicopathological features were analyzed. Kaplan-Meier survival analysis was performed to estimate the prognostic value of the PINK1 and PARK2 proteins in patients. Cox univariate and multivariate analyses were used to assess the risk factors affecting the OS for patients with ESCC.

RESULTS

PINK1 and PARK2 had low expression in ESCC. Patients with low PINK1 had worse differentiation and advanced T and TNM stages. Lower PARK2 expression was linked to lymph node metastases and an advanced TNM stage. Furthermore, reduced PINK1 and PARK2 levels were associated with a poor prognosis for ESCC. Cox univariate and multivariate analyses revealed that PINK1, PARK2, and tumor size were closely associated with the prognosis of patients with ESCC, and PARK2 was an independent risk factor for patients with ESCC. Finally, the PINK1 and PARK2 proteins were closely related and shared the same signal pathway.

CONCLUSIONS

PINK1 and PARK2 could work as tumor suppressors in ESCC and are likely to become new treatment targets for ESCC.

The Long Non-Coding RNA ANRIL in Cancers

ANRIL (Antisense Noncoding RNA in the INK4 Locus), a long non-coding RNA encoded in the human chromosome 9p21 region, is a critical factor for regulating gene expression by interacting with multiple proteins and miRNAs. It has been found to play important roles in various cellular processes, including cell cycle control and proliferation. Dysregulation of ANRIL has been associated with several diseases like cancers and cardiovascular diseases, for instance. Understanding the oncogenic role of ANRIL and its potential as a diagnostic and prognostic biomarker in cancer is crucial. This review provides insights into the regulatory mechanisms and oncogenic significance of the 9p21 locus and ANRIL in cancer.

Functional Genetic Variants in TGFβ1 and TGFβR1 in miRNA-Binding Sites Predict Outcomes in Patients with HPV-positive Oropharyngeal Squamous Cell Carcinoma

PURPOSE

TGFβ1 and TGFβ receptor 1 (TGFβR1) participate in regulation of the host's immune system and inflammatory responses and may serve as prognostic biomarkers for human papillomavirus (HPV)-associated oropharyngeal squamous cell carcinoma (OPSCC).

EXPERIMENTAL DESIGN

This study included 1,013 patients with incident OPSCC, of whom 489 had tumor HPV16 status determined. All patients were genotyped for two functional polymorphisms: TGFβ1 rs1800470 and TGFβR1 rs334348. Univariate and multivariate Cox regression models were performed to evaluate associations between the polymorphisms and overall survival (OS), disease-specific survival (DSS), and disease-free survival (DFS).

RESULTS

Patients with TGFβ1 rs1800470 CT or CC genotype had 70%-80% reduced risks of OS, DSS, and DFS compared with patients with TT genotype, and patients with TGFβR1 rs334348 GA or GG genotype had 30%-40% reduced risk of OS, DSS, and DFS compared with patients with AA genotype. Furthermore, among patients with HPV-positive (HPV+) OPSCC, the same patterns were observed but the risk reductions were greater: up to 80%-90% for TGFβ1 rs1800470 CT or CC genotype and 70%-85% for TGFβR1 rs334348 GA or GG genotype. The risk reductions were still greater (up to 17 to 25 times reduced) for patients with both TGFβ1 rs1800470 CT or CC genotype and TGFβR1 rs334348 GA or GG genotype compared with patients with both TGFβ1 rs1800470 TT genotype and TGFβR1 rs334348 AA genotype among patients with HPV+ OPSCC.

CONCLUSIONS

Our findings indicate that TGFβ1 rs1800470 and TGFβR1 rs334348 may individually or jointly modify risks of death and recurrence in patients with OPSCC, particularly those with HPV+ OPSCC undergoing definitive radiotherapy, and may serve as prognostic biomarkers, which could lead to better personalized treatment and improved prognosis.

Pan-cancer and cross-population genome-wide association studies dissect shared genetic backgrounds underlying carcinogenesis

Integrating genomic data of multiple cancers allows de novo cancer grouping and elucidating the shared genetic basis across cancers. Here, we conduct the pan-cancer and cross-population genome-wide association study (GWAS) meta-analysis and replication studies on 13 cancers including 250,015 East Asians (Biobank Japan) and 377,441 Europeans (UK Biobank). We identify ten cancer risk variants including five pleiotropic associations (e.g., rs2076295 at DSP on 6p24 associated with lung cancer and rs2525548 at TRIM4 on 7q22 nominally associated with six cancers). Quantifying shared heritability among the cancers detects positive genetic correlations between breast and prostate cancer across populations. Common genetic components increase the statistical power, and the large-scale meta-analysis of 277,896 breast/prostate cancer cases and 901,858 controls identifies 91 newly genome-wide significant loci. Enrichment analysis of pathways and cell types reveals shared genetic backgrounds across said cancers. Focusing on genetically correlated cancers can contribute to enhancing our insights into carcinogenesis.

Computational screening and analysis of deleterious nsSNPs in human p14ARF (CDKN2A gene) protein using molecular dynamic simulation approach

Cyclin-dependent kinase inhibitor 2 A (CDKN2A) gene belongs to the cyclin-dependent kinase family that code for two transcripts (p16INK4A and p14ARF), both work as tumor suppressors proteins. The mutation that occurs in the p14ARF protein can lead to different types of cancers. Single nucleotide polymorphisms (SNPs) are an important type of genetic alteration that can lead to different types of diseases. In this study, we applied the computational strategy on human p14ARF protein to identify the potential deleterious nsSNPs and check their impact on the structure, function, and protein stability. We applied more than ten prediction tools to screen the retrieved 288 nsSNPs, consequently extracting four deleterious nsSNPs i.e., rs139725688 (R10G), rs139725688 (R21W), rs374360796 (F23L) and rs747717236 (L124R). Homology modeling, conservation and conformational analysis of mutant models were performed to examine the divergence of these variants from the native p14ARF structure. All-atom molecular dynamics simulation revealed a significant impact of these mutations on protein stability, compactness, globularity, solvent accessibility and secondary structure elements. Protein-protein interactions indicated that p14ARF operates as a hub linking clusters of different proteins and that changes in p14ARF may result in the disassociation of numerous signal cascades. Our current study is the first survey of computational analysis on p14ARF protein that determines the association of these nsSNPs with the altered function of p14ARF protein and leads to the development of various types of cancers. This research proposes the described functional SNPs as possible targets for proteomic investigations, diagnostic procedures, and treatments.Communicated by Ramaswamy H. Sarma.

Immunogenetic Determinants of Susceptibility to Head and Neck Cancer in the Million Veteran Program Cohort

Increasing rates of human papillomavirus (HPV)-driven oropharyngeal cancer (OPC) have largely offset declines in tobacco-associated head and neck squamous cell carcinoma (HNSCC) at non-OPC sites. Host immunity is an important modulator of HPV infection, persistence, and clearance, and also of immune evasion in both virally- and nonvirally-driven cancers. However, the association between collective known cancer-related immune gene variants and HNSCC susceptibility has not been fully characterized. Here, we conducted a genetic association study in the multiethnic Veterans Affairs Million Veteran Program cohort, evaluating 16,050 variants in 1,576 immune genes in 4,012 HNSCC cases (OPC = 1,823; non-OPC = 2,189) and 16,048 matched controls. Significant polymorphisms were further examined in a non-Hispanic white (NHW) validation cohort (OPC = 1,206; non-OPC = 955; controls = 4,507). For overall HNSCC susceptibility in NHWs, we discovered and validated a novel 9q31.1 SMC2 association and replicated the known 6p21.32 HLA-DQ-DR association. Six loci/genes for overall HNSCC susceptibility were selectively enriched in African-Americans (6p21.32 HLA-G, 9q21.33 GAS1, 11q12.2 CD6, 11q23.2 NCAM1/CD56, 17p13.1 CD68, 18q22.2 SOCS6); all 6 genes function in antigen-presenting regulation and T-cell activation. Two additional loci (10q26 DMBT1, 15q22.2 TPM1) were uncovered for non-OPC susceptibility, and three loci (11q24 CRTAM, 16q21 CDH5, 18q12.1 CDH2) were identified for HPV-positive OPC susceptibility. This study underscores the role of immune gene variants in modulating susceptibility for both HPV-driven and non-HPV-driven HNSCC. Additional large studies, particularly in racially diverse populations, are needed to further validate the associations and to help elucidate other potential immune factors and mechanisms that may underlie HNSCC risk.

SIGNIFICANCE

Several inherited variations in immune system genes are significantly associated with susceptibility to head and neck cancer, which could help improve personalized cancer risk estimates.

Integrative Multiomics to Dissect the Lung Transcriptional Landscape of Pulmonary Arterial Hypertension

Pulmonary arterial hypertension (PAH) remains an incurable and often fatal disease despite currently available therapies. Multiomics systems biology analysis can shed new light on PAH pathobiology and inform translational research efforts. Using RNA sequencing on the largest PAH lung biobank to date (96 disease and 52 control), we aim to identify gene co-expression network modules associated with PAH and potential therapeutic targets. Co-expression network analysis was performed to identify modules of co-expressed genes which were then assessed for and prioritized by importance in PAH, regulatory role, and therapeutic potential via integration with clinicopathologic data, human genome-wide association studies (GWAS) of PAH, lung Bayesian regulatory networks, single-cell RNA-sequencing data, and pharmacotranscriptomic profiles. We identified a co-expression module of 266 genes, called the pink module, which may be a response to the underlying disease process to counteract disease progression in PAH. This module was associated not only with PAH severity such as increased PVR and intimal thickness, but also with compensated PAH such as lower number of hospitalizations, WHO functional class and NT-proBNP. GWAS integration demonstrated the pink module is enriched for PAH-associated genetic variation in multiple cohorts. Regulatory network analysis revealed that BMPR2 regulates the main target of FDA-approved riociguat, GUCY1A2, in the pink module. Analysis of pathway enrichment and pink hub genes (i.e. ANTXR1 and SFRP4) suggests the pink module inhibits Wnt signaling and epithelial-mesenchymal transition. Cell type deconvolution showed the pink module correlates with higher vascular cell fractions (i.e. myofibroblasts). A pharmacotranscriptomic screen discovered ubiquitin-specific peptidases (USPs) as potential therapeutic targets to mimic the pink module signature. Our multiomics integrative study uncovered a novel gene subnetwork associated with clinicopathologic severity, genetic risk, specific vascular cell types, and new therapeutic targets in PAH. Future studies are warranted to investigate the role and therapeutic potential of the pink module and targeting USPs in PAH.

Selecting the appropriate hurdles and endpoints for pentilludin, a novel antiaddiction pharmacotherapeutic targeting the receptor type protein tyrosine phosphatase D

Substance use disorders provide challenges for development of effective medications. Use of abused substances is likely initiated, sustained and "quit" by complex brain and pharmacological mechanisms that have both genetic and environmental determinants. Medical utilities of prescribed stimulants and opioids provide complex challenges for prevention: how can we minimize their contribution to substance use disorders while retaining medical benefits for pain, restless leg syndrome, attention deficit hyperactivity disorder, narcolepsy and other indications. Data required to support assessments of reduced abuse liability and resulting regulatory scheduling differs from information required to support licensing of novel prophylactic or therapeutic anti-addiction medications, adding further complexity and challenges. I describe some of these challenges in the context of our current efforts to develop pentilludin as a novel anti-addiction therapeutic for a target that is strongly supported by human and mouse genetic and pharmacologic studies, the receptor type protein tyrosine phosphatase D (PTPRD).

MDM4: What do we know about the association between its polymorphisms and cancer?

MDM4 is an important p53-negative regulator, consequently, it is involved in cell proliferation, DNA repair, and apoptosis regulation. MDM4 overexpression and amplification are described to lead to cancer formation, metastasis, and poor disease prognosis. Several MDM4 SNPs are in non-coding regions, and some affect the MDM4 regulation by disrupting the micro RNA binding site in 3'UTR (untranslated region). Here, we gathered several association studies with different MDM4 SNPs and populations to understand the relationship between its SNPs and solid tumor risk. Many studies failed to replicate their results regarding different populations, cancer types, and risk genotypes, leading to conflicting conclusions. We suggested that distinct haplotype patterns in different populations might affect the association between MDM4 SNPs and cancer risk. Thus, we propose to investigate some linkage SNPs in specific haplotypes to provide informative MDM4 markers for association studies with cancer.

Exploring the genetic overlap between twelve psychiatric disorders

The widespread comorbidity among psychiatric disorders demonstrated in epidemiological studies1-5 is mirrored by non-zero, positive genetic correlations from large-scale genetic studies6-10. To identify shared biological processes underpinning this observed phenotypic and genetic covariance and enhance molecular characterization of general psychiatric disorder liability11-13, we used several strategies aimed at uncovering pleiotropic, that is, cross-trait-associated, single-nucleotide polymorphisms (SNPs), genes and biological pathways. We conducted cross-trait meta-analysis on 12 psychiatric disorders to identify pleiotropic SNPs. The meta-analytic signal was driven by schizophrenia, hampering interpretation and joint biological characterization of the cross-trait meta-analytic signal. Subsequent pairwise comparisons of psychiatric disorders identified substantial pleiotropic overlap, but mainly among pairs of psychiatric disorders, and mainly at less stringent P-value thresholds. Only annotations related to evolutionarily conserved genomic regions were significant for multiple (9 out of 12) psychiatric disorders. Overall, identification of shared biological mechanisms remains challenging due to variation in power and genetic architecture between psychiatric disorders.

Approach to Localized Squamous Cell Cancer of the Esophagus

OPINION STATEMENT

Esophageal cancer is a leading cause of cancer deaths worldwide, with an increasing incidence in recent decades. The majority of esophageal cancers are squamous cell carcinoma. The 5-year survival rate of esophageal squamous cell carcinoma (ESCC) is poor, and there remains globally a pressing need for novel treatments that improve patient outcomes and quality of life. In this review, we discuss management of localized ESCC with an update on relevant newly published literature, including targeted therapy and novel biomarkers. The standard treatment approach for locally advanced, resectable ESCC is currently chemoradiation with or without surgery. Here we discuss different approaches to endoscopic resection, surgery, and radiation therapy. Although the typical chemotherapy regimen is a combination of a platinum with a fluoropyrimidine or paclitaxel, different regimens are being evaluated. With the landscape of immunotherapy rapidly evolving, at the forefront of new treatments for ESCC is immunotherapy and other targeted agents. Ultimately, the treatment approach should be individualized to each patient.

Genetic Analysis of Lung Cancer and the Germline Impact on Somatic Mutation Burden

BACKGROUND

Germline genetic variation contributes to lung cancer (LC) susceptibility. Previous genome-wide association studies (GWAS) have implicated susceptibility loci involved in smoking behaviors and DNA repair genes, but further work is required to identify susceptibility variants.

METHODS

To identify LC susceptibility loci, a family history-based genome-wide association by proxy (GWAx) of LC (48 843 European proxy LC patients, 195 387 controls) was combined with a previous LC GWAS (29 266 patients, 56 450 controls) by meta-analysis. Colocalization was used to explore candidate genes and overlap with existing traits at discovered susceptibility loci. Polygenic risk scores (PRS) were tested within an independent validation cohort (1 666 LC patients vs 6 664 controls) using variants selected from the LC susceptibility loci and a novel selection approach using published GWAS summary statistics. Finally, the effects of the LC PRS on somatic mutational burden were explored in patients whose tumor resections have been profiled by exome (n = 685) and genome sequencing (n = 61). Statistical tests were 2-sided.

RESULTS

The GWAx-GWAS meta-analysis identified 8 novel LC loci. Colocalization implicated DNA repair genes (CHEK1), metabolic genes (CYP1A1), and smoking propensity genes (CHRNA4 and CHRNB2). PRS analysis demonstrated that these variants, as well as subgenome-wide significant variants related to expression quantitative trait loci and/or smoking propensity, assisted in LC genetic risk prediction (odds ratio = 1.37, 95% confidence interval = 1.29 to 1.45; P < .001). Patients with higher genetic PRS loads of smoking-related variants tended to have higher mutation burdens in their lung tumors.

CONCLUSIONS

This study has expanded the number of LC susceptibility loci and provided insights into the molecular mechanisms by which these susceptibility variants contribute to LC development.

Functional studies of lung cancer GWAS beyond association

Fourteen years after the first genome-wide association study (GWAS) of lung cancer was published, approximately forty-five genomic loci have now been significantly associated with lung cancer risk. While functional characterization was performed for several of these loci, a comprehensive summary of current molecular understanding of lung cancer risk has been lacking. Further, many novel computational and experimental tools now became available to accelerate the functional assessment of disease-associated variants, moving beyond locus-by-locus approaches. In this review, we first highlight the heterogeneity of lung cancer GWAS findings across histological subtypes, ancestries, and smoking status, which poses unique challenges to follow-up studies. We then summarize the published lung cancer post-GWAS studies for each risk-associated locus to assess the current understanding of biological mechanisms beyond the initial statistical association. We further summarize strategies for GWAS functional follow-up studies considering cutting-edge functional genomics tools and providing a catalog of available resources relevant to lung cancer. Overall, we aim to highlight the importance of integrating computational and experimental approaches to draw biological insights from the lung cancer GWAS results beyond association.

Common variants in breast cancer risk loci predispose to distinct tumor subtypes

BACKGROUND

Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER) status, but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear.

METHODS

Among 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes.

RESULTS

Eighty-five of 173 variants were associated with at least one tumor feature (false discovery rate < 5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at p < 0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions.

CONCLUSION

This report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction.

Genetic Predisposition to Persistent Human Papillomavirus-Infection and Virus-Induced Cancers

Human papillomaviruses (HPVs) are the most common sexually transmitted pathogens worldwide and among the more than 200 identified HPV types, approximately 15 high risk (HR-HPV) types are oncogenic, being strongly associated with the development of cervical cancer, anogenital cancers and an increasing fraction of head and neck squamous cell carcinomas (HNSCC). HPV-associated cervix cancer accounts for 83% of HPV-attributable cancers, and more than two-thirds of those cases occur in developing countries. Despite the high frequency of HPV infections, in most cases, the virus is cleared by the host immune response and only a small proportion of infected individuals develop persistent infections that can result in malignant transformation, indicating that other elements, including biological, genetic and environmental factors may influence the individual susceptibility to HPV-associated cancers. Previous studies have quantified that heritability, in the form of genetic variants, common in the general population, is implicated in nearly 30% of cervical cancers and a large number of studies conducted across various populations have identified genetic variants that appear to be associated with genes that predispose or protect the host to HPV infections thereby affecting individual susceptibility to HPV-associated cancers. In this article, we provide an overview of gene association studies on HPV-associated cancers with emphasis on genome-wide association study (GWAS) that have identified novel genetic factors linked to HPV infection or HPV-associated cancers.