Identification of shared and unique susceptibility pathways among cancers of the lung, breast, and prostate from genome-wide association studies and tissue-specific protein interactions

The relationship between 11 different polygenic longevity scores, parental lifespan, and disease diagnosis in the UK Biobank

Large-scale genome-wide association studies (GWAS) strongly suggest that most traits and diseases have a polygenic component. This observation has motivated the development of disease-specific "polygenic scores (PGS)" that are weighted sums of the effects of disease-associated variants identified from GWAS that correlate with an individual's likelihood of expressing a specific phenotype. Although most GWAS have been pursued on disease traits, leading to the creation of refined "Polygenic Risk Scores" (PRS) that quantify risk to diseases, many GWAS have also been pursued on extreme human longevity, general fitness, health span, and other health-positive traits. These GWAS have discovered many genetic variants seemingly protective from disease and are often different from disease-associated variants (i.e., they are not just alternative alleles at disease-associated loci) and suggest that many health-positive traits also have a polygenic basis. This observation has led to an interest in "polygenic longevity scores (PLS)" that quantify the "risk" or genetic predisposition of an individual towards health. We derived 11 different PLS from 4 different available GWAS on lifespan and then investigated the properties of these PLS using data from the UK Biobank (UKB). Tests of association between the PLS and population structure, parental lifespan, and several cancerous and non-cancerous diseases, including death from COVID-19, were performed. Based on the results of our analyses, we argue that PLS are made up of variants not only robustly associated with parental lifespan, but that also contribute to the genetic architecture of disease susceptibility, morbidity, and mortality.

Investigating the tissue specificity and prognostic impact of cis-regulatory cancer risk variants

The tissue-specific incidence of cancers and their genetic basis are poorly understood. Although prior studies have shown global correlation across tissues for cancer risk single-nucleotide polymorphisms (SNPs) identified through genome-wide association studies (GWAS), any shared functional regulation of gene expression on a per SNP basis has not been well characterized. We set to quantify cis-mediated gene regulation and tissue sharing for SNPs associated with eight common cancers. We identify significant tissue sharing for individual SNPs and global enrichment for breast, colorectal, and Hodgkin lymphoma cancer risk SNPs in multiple tissues. In addition, we observe increasing tissue sharing for cancer risk SNPs overlapping with super-enhancers for breast cancer and Hodgkin lymphoma providing further evidence of tissue specificity. Finally, for genes under cis-regulation by breast cancer SNPs, we identify a phenotype characterized by low expression of tumor suppressors and negative regulators of the WNT pathway associated with worse freedom from progression and overall survival in patients who eventually develop breast cancer. Our results introduce a paradigm for functionally annotating individual cancer risk SNPs and will inform the design of future translational studies aimed to personalize assessment of inherited cancer risk across tissues.

Rare variants in the endocytic pathway are associated with Alzheimer's disease, its related phenotypes, and functional consequences

Late-onset Alzheimer’s disease (LOAD) is the most common type of dementia causing irreversible brain damage to the elderly and presents a major public health challenge. Clinical research and genome-wide association studies have suggested a potential contribution of the endocytic pathway to AD, with an emphasis on common loci. However, the contribution of rare variants in this pathway to AD has not been thoroughly investigated. In this study, we focused on the effect of rare variants on AD by first applying a rare-variant gene-set burden analysis using genes in the endocytic pathway on over 3,000 individuals with European ancestry from three large whole-genome sequencing (WGS) studies. We identified significant associations of rare-variant burden within the endocytic pathway with AD, which were successfully replicated in independent datasets. We further demonstrated that this endocytic rare-variant enrichment is associated with neurofibrillary tangles (NFTs) and age-related phenotypes, increasing the risk of obtaining severer brain damage, earlier age-at-onset, and earlier age-of-death. Next, by aggregating rare variants within each gene, we sought to identify single endocytic genes associated with AD and NFTs. Careful examination using NFTs revealed one significantly associated gene, ANKRD13D. To identify functional associations, we integrated bulk RNA-Seq data from over 600 brain tissues and found two endocytic expression genes (eGenes), HLA-A and SLC26A7, that displayed significant influences on their gene expressions. Differential expressions between AD patients and controls of these three identified genes were further examined by incorporating scRNA-Seq data from 48 post-mortem brain samples and demonstrated distinct expression patterns across cell types. Taken together, our results demonstrated strong rare-variant effect in the endocytic pathway on AD risk and progression and functional effect of gene expression alteration in both bulk and single-cell resolution, which may bring more insight and serve as valuable resources for future AD genetic studies, clinical research, and therapeutic targeting.

Late-onset Alzheimer’s disease (LOAD) is the most common type of dementia and a leading cause of death in the world. Clinical and genetic studies have suggested the potential contribution of the cellular transportation pathway to AD with an emphasis on common variants. In this study, we investigated the effect of rare variants within the cellular transportation pathway and examined three large datasets with over 3,000 individuals with European ancestry. We reported enrichment of rare deleterious variants in the cellular transportation pathway in AD patients from all three datasets. We also observed an elevation of rare deleterious variants in this pathway was associated with individuals with severer brain damages (AD progression), earlier age-at-onset, and earlier age-of-death. By aggregating rare variants in each gene from the cellular transportation pathway, we revealed one gene in which rare variants were significantly associated with the progression of AD. By integrating gene expression data from brain tissues, we identified two additional genes whose rare-variant effect displayed significant influences on gene expression. Taken together, our results demonstrated that rare-variant effect in the cellular transportation pathway is strongly associated with the risk and the progression of AD, which may serve as future clinical and therapeutic targets.

Detecting Rewiring Events in Protein-Protein Interaction Networks Based on Transcriptomic Data

Proteins rarely carry out their cellular functions in isolation. Instead, eukaryotic proteins engage in about six interactions with other proteins on average. The aggregated protein interactome of an organism forms a “hairy ball”-type protein-protein interaction (PPI) network. Yet, in a typical human cell, only about half of all proteins are expressed at a particular time. Hence, it has become common practice to prune the full PPI network to the subset of expressed proteins. If RNAseq data is available, one can further resolve the specific protein isoforms present in a cell or tissue. Here, we review various approaches, software tools and webservices that enable users to construct context-specific or tissue-specific PPI networks and how these are rewired between two cellular conditions. We illustrate their different functionalities on the example of the interactions involving the human TNR6 protein. In an outlook, we describe how PPI networks may be integrated with epigenetic data or with data on the activity of splicing factors.

Cross-cancer evaluation of polygenic risk scores for 16 cancer types in two large cohorts

Even distinct cancer types share biological hallmarks. Here, we investigate polygenic risk score (PRS)-specific pleiotropy across 16 cancers in European ancestry individuals from the Genetic Epidemiology Research on Adult Health and Aging cohort (16,012 cases, 50,552 controls) and UK Biobank (48,969 cases, 359,802 controls). Within cohorts, each PRS is evaluated in multivariable logistic regression models against all other cancer types. Results are then meta-analyzed across cohorts. Ten positive and one inverse cross-cancer associations are found after multiple testing correction. Two pairs show bidirectional associations; the melanoma PRS is positively associated with oral cavity/pharyngeal cancer and vice versa, whereas the lung cancer PRS is positively associated with oral cavity/pharyngeal cancer, and the oral cavity/pharyngeal cancer PRS is inversely associated with lung cancer. Overall, we validate known, and uncover previously unreported, patterns of pleiotropy that have the potential to inform investigations of risk prediction, shared etiology, and precision cancer prevention strategies.

Pan-cancer study detects genetic risk variants and shared genetic basis in two large cohorts

Deciphering the shared genetic basis of distinct cancers has the potential to elucidate carcinogenic mechanisms and inform broadly applicable risk assessment efforts. Here, we undertake genome-wide association studies (GWAS) and comprehensive evaluations of heritability and pleiotropy across 18 cancer types in two large, population-based cohorts: the UK Biobank (408,786 European ancestry individuals; 48,961 cancer cases) and the Kaiser Permanente Genetic Epidemiology Research on Adult Health and Aging cohorts (66,526 European ancestry individuals; 16,001 cancer cases). The GWAS detect 21 genome-wide significant associations independent of previously reported results. Investigations of pleiotropy identify 12 cancer pairs exhibiting either positive or negative genetic correlations; 25 pleiotropic loci; and 100 independent pleiotropic variants, many of which are regulatory elements and/or influence cross-tissue gene expression. Our findings demonstrate widespread pleiotropy and offer further insight into the complex genetic architecture of cross-cancer susceptibility.

Association between integrin-linked kinase and hyperthermia in oral squamous cell carcinoma

The present study aimed to observe the effect of the biological functions of integrin-linked kinase (ILK) silencing combined with hyperthermia on Tca8113 cells. Lentivirus-mediated short hairpin RNA (shRNA)-targeting ILK was transfected into oral squamous cell carcinoma (OSCC) Tca8113 cells and, combined with hyperthermia, several experimental methods were used to detect their biological behavior in vitro. On the basis of in vitro experiments, Tca8113 cells were transplanted into nude mice models, and ILK-shRNA-lentivirus was injected into the nude mice transplanted tumor and combined with hyperthermia. Tumor morphology and the associated protein expression changes were determined. Subsequent to ILK silencing combined with hyperthermia, the growth, migration and proliferation of Tca8113 cells were significantly inhibited. Flow cytometry revealed that the cells were blocked in the S phase, and western blot analysis demonstrated that ILK, phosphorylated (p)-RAC-alpha serine/threonine-protein kinase (Akt), p-glycogen synthase kinase-3β and p-heat shock factor 1 protein expression levels were significantly decreased, while apoptosis-associated protein B-cell lymphoma-2-associated X protein expression and the efficacy of hypothermia were significantly increased. By ILK silencing combined with hyperthermia, a significant therapeutic effect on transplanted tumors was observed in nude mice. Immunohistochemistry revealed the same results as the in vitro experiments. ILK silencing combined with hyperthermia can inhibit the growth, proliferation and migration of Tca8113 cells, promote Tca8113 cell apoptosis, inhibit the phosphatidylinositol-3-kinase/Akt signaling pathway and increase hyperthermia sensitivity; the combination therapy exhibits a synergistic sensitizing effect. Therefore, ILK silencing combined with hypothermia may serve as a novel combination therapy strategy against OSCC.

The statistical properties of gene-set analysis

The rapid increase in loci discovered in genome-wide association studies has created a need to understand the biological implications of these results. Gene-set analysis provides a means of gaining such understanding, but the statistical properties of gene-set analysis are not well understood, which compromises our ability to interpret its results. In this Analysis article, we provide an extensive statistical evaluation of the core structure that is inherent to all gene- set analyses and we examine current implementations in available tools. We show which factors affect valid and successful detection of gene sets and which provide a solid foundation for performing and interpreting gene-set analysis.

Enrichment of putative PAX8 target genes at serous epithelial ovarian cancer susceptibility loci

BACKGROUND

Genome-wide association studies (GWAS) have identified 18 loci associated with serous ovarian cancer (SOC) susceptibility but the biological mechanisms driving these findings remain poorly characterised. Germline cancer risk loci may be enriched for target genes of transcription factors (TFs) critical to somatic tumorigenesis.

METHODS

All 615 TF-target sets from the Molecular Signatures Database were evaluated using gene set enrichment analysis (GSEA) and three GWAS for SOC risk: discovery (2196 cases/4396 controls), replication (7035 cases/21 693 controls; independent from discovery), and combined (9627 cases/30 845 controls; including additional individuals).

RESULTS

The PAX8-target gene set was ranked 1/615 in the discovery (PGSEA<0.001; FDR=0.21), 7/615 in the replication (PGSEA=0.004; FDR=0.37), and 1/615 in the combined (PGSEA<0.001; FDR=0.21) studies. Adding other genes reported to interact with PAX8 in the literature to the PAX8-target set and applying an alternative to GSEA, interval enrichment, further confirmed this association (P=0.006). Fifteen of the 157 genes from this expanded PAX8 pathway were near eight loci associated with SOC risk at P<10-5 (including six with P<5 × 10-8). The pathway was also associated with differential gene expression after shRNA-mediated silencing of PAX8 in HeyA8 (PGSEA=0.025) and IGROV1 (PGSEA=0.004) SOC cells and several PAX8 targets near SOC risk loci demonstrated in vitro transcriptomic perturbation.

CONCLUSIONS

Putative PAX8 target genes are enriched for common SOC risk variants. This finding from our agnostic evaluation is of particular interest given that PAX8 is well-established as a specific marker for the cell of origin of SOC.

Familial Lung Cancer: A Brief History from the Earliest Work to the Most Recent Studies

Lung cancer is the deadliest cancer in the United States, killing roughly one of four cancer patients in 2016. While it is well-established that lung cancer is caused primarily by environmental effects (particularly tobacco smoking), there is evidence for genetic susceptibility. Lung cancer has been shown to aggregate in families, and segregation analyses have hypothesized a major susceptibility locus for the disease. Genetic association studies have provided strong evidence for common risk variants of small-to-moderate effect. Rare and highly penetrant alleles have been identified by linkage studies, including on 6q23-25. Though not common, some germline mutations have also been identified via sequencing studies. Ongoing genomics studies aim to identify additional high penetrance germline susceptibility alleles for this deadly disease.

PI3K/Akt/mTOR Signaling and Plasma Membrane Proteins Are Implicated in Responsiveness to Adjuvant Dendritic Cell Vaccination for Metastatic Colorectal Cancer

PURPOSE

We have previously demonstrated that patients with metastatic colorectal cancer who exhibit immune responses to a dendritic cell (DC) vaccine have superior recurrence-free survival following surgery, compared with patients in whom responses do not occur. We sought to characterize the patterns of T-lymphocyte infiltration and somatic mutations in metastases that are associated with and predictive of response to the DC vaccine.

EXPERIMENTAL DESIGN

Cytotoxic, memory, and regulatory T cells in resected metastases and surrounding normal liver tissue from 22 patients (11 responders and 11 nonresponders) were enumerated by immunohistochemistry prior to vaccine administration. In conjunction with tumor sequencing, the combined multivariate and collapsing method was used to identify gene mutations that are associated with vaccine response. We also derived a response prediction score for each patient using his/her tumor genotype data and variant association effect sizes computed from the other 21 patients; greater weighting was placed on gene products with cell membrane-related functions.

RESULTS

There was no correlation between vaccine response and intratumor, peritumor, or hepatic densities of T-cell subpopulations. Associated genes were found to be enriched in the PI3K/Akt/mTOR signaling axis (P < 0.001). Applying a consistent prediction score cutoff over 22 rounds of leave-one-out cross-validation correctly inferred vaccine response in 21 of 22 patients (95%).

CONCLUSIONS

Adjuvant DC vaccination has shown promise as a form of immunotherapy for patients with metastatic colorectal cancer. Its efficacy may be influenced by somatic mutations that affect pathways involving PI3K, Akt, and mTOR, as well as tumor surface proteins. Clin Cancer Res; 23(2); 399-406. ©2016 AACR.