Genome-wide association study identifies 14 novel risk alleles associated with basal cell carcinoma

PTGES2 and RNASET2 identified as novel potential biomarkers and therapeutic targets for basal cell carcinoma: insights from proteome-wide mendelian randomization, colocalization, and MR-PheWAS analyses

Introduction

Basal cell carcinoma (BCC) is the most common skin cancer, lacking reliable biomarkers or therapeutic targets for effective treatment. Genome-wide association studies (GWAS) can aid in identifying drug targets, repurposing existing drugs, predicting clinical trial side effects, and reclassifying patients in clinical utility. Hence, the present study investigates the association between plasma proteins and skin cancer to identify effective biomarkers and therapeutic targets for BCC.

Methods

Proteome-wide mendelian randomization was performed using inverse-variance-weight and Wald Ratio methods, leveraging 1 Mb cis protein quantitative trait loci (cis-pQTLs) in the UK Biobank Pharma Proteomics Project (UKB-PPP) and the deCODE Health Study, to determine the causal relationship between plasma proteins and skin cancer and its subtypes in the FinnGen R10 study and the SAIGE database of Lee lab. Significant association with skin cancer and its subtypes was defined as a false discovery rate (FDR) < 0.05. pQTL to GWAS colocalization analysis was executed using a Bayesian model to evaluate five exclusive hypotheses. Strong colocalization evidence was defined as a posterior probability for shared causal variants (PP.H4) of ≥0.85. Mendelian randomization-Phenome-wide association studies (MR-PheWAS) were used to evaluate potential biomarkers and therapeutic targets for skin cancer and its subtypes within a phenome-wide human disease category.

Results

PTGES2, RNASET2, SF3B4, STX8, ENO2, and HS3ST3B1 (besides RNASET2, five other plasma proteins were previously unknown in expression quantitative trait loci (eQTL) and methylation quantitative trait loci (mQTL)) were significantly associated with BCC after FDR correction in the UKB-PPP and deCODE studies. Reverse MR showed no association between BCC and these proteins. PTGES2 and RNASET2 exhibited strong evidence of colocalization with BCC based on a posterior probability PP.H4 >0.92. Furthermore, MR-PheWAS analysis showed that BCC was the most significant phenotype associated with PTGES2 and RNASET2 among 2,408 phenotypes in the FinnGen R10 study. Therefore, PTGES2 and RNASET2 are highlighted as effective biomarkers and therapeutic targets for BCC within the phenome-wide human disease category.

Conclusion

The study identifies PTGES2 and RNASET2 plasma proteins as novel, reliable biomarkers and therapeutic targets for BCC, suggesting more effective clinical application strategies for patients.

Retrospective clinical study analysis of skin adverse reactions related to epidermal growth factor receptor inhibitors

Background

Epidermal growth factor receptor inhibitors (EGFRIs) represent a cornerstone in the targeted therapy of malignant tumors. While effective, dermatological adverse events (dAEs) associated with EGFRIs pose a significant challenge, often necessitating treatment discontinuation due to their severity and potential to impede the continuity of cancer therapy. Despite extensive research, the specific mechanisms and predictors of these adverse events remain poorly understood, particularly in diverse populations. This gap in knowledge underscores the need for targeted studies to better predict and manage these events, enhancing patient outcomes and adherence to life-saving therapies.

Methods

This observational study was conducted at The First Affiliated Hospital of Guangxi Medical University, covering cancer patients treated with EGFRIs from 2020 to 2022. We analyzed clinical data including patient demographics, treatment specifics, and the development and timing of dAEs. The study employed SPSS 26.0 software for data analysis, focusing on the incidence of dAEs and factors influencing their occurrence. We used Kaplan-Meier and Cox regression methods to establish a predictive model for dAEs, tracking their onset and impact on treatment continuity.

Results

In our study of 120 patients treated with EGFR inhibitors at The First Affiliated Hospital of Guangxi Medical University, we found a high prevalence of dAEs, with 84.2% of patients experiencing such effects. The most common manifestations were papulopustular rashes, observed as pustules in 52.5% and papules in 57.4% of cases, followed by nail lesions in 62.4% of patients, oral or other mucosal ulcers in 34.7%, and hair changes in 26.7%. The median incubation time (MIT) for dAEs was 5 weeks. We identified drug type, ethnicity, and occupation as statistically significant risk factors (P<0.05 for all) that influenced the MIT, which the Cox regression model further identified as protective factors. Nomograms were developed to assess the risk of dAEs, although it is important to note that these models have only been internally validated, lacking external validation data at this stage.

Conclusions

The study highlights the high incidence of EGFRIs-associated dAEs, with specific dermatological manifestations posing significant challenges in cancer therapy. The identification of drug type, ethnicity, and occupation as influential factors on the MIT for dAEs informs clinical decisions. Our prediction model serves as a practical tool for evaluating the risk of developing dAEs over time, aiming to optimize patient management and mitigate treatment interruptions.

Comparisons between wrinkles and photo-ageing detected and self-reported by the participant or identified by trained assessors reveal insights from Chinese individuals in the Singapore/Malaysia Cross-sectional Genetics Epidemiology Study (SMCGES) cohort

BACKGROUND

Changes develop on the facial skin as a person ages. Other than chronological time, it has been discovered that gender, ethnicity, air pollution, smoking, nutrition, and sun exposure are notable risk factors that influence the development of skin ageing phenotypes such as wrinkles and photo-ageing. These risk factors can be quantified through epidemiological collection methods. We previously studied wrinkles and photo-ageing in detail using photo-numeric scales. The analysis was performed on the ethnic Chinese skin by three trained assessors. Recent studies have shown that it is possible to use self-reported data to identify skin-related changes including skin colour and skin cancer. In order to investigate the association between risk factors and skin ageing phenotypic outcomes in large-scale epidemiological studies, it would be useful to evaluate whether it is also possible for participants to self-report signs of ageing on their skin.

AIM

We have previously identified several validated photo-numeric scales for wrinkling and photo-ageing to use on ethnic Chinese skin. Using these scales, our trained assessors grade wrinkling and photo-ageing with moderately high inter-assessor concordance and agreement. The main objective of this study involves letting participants grade self-reported wrinkling and photo-ageing using these same scales. We aim to compare the concordance and agreement between signs of skin ageing by the participant and signs of ageing identified by our assessors.

METHOD

Three trained assessors studied facial photo-ageing on 1081 ethnic Chinese young adults from the Singapore/Malaysia Cross-sectional Genetics Epidemiology Study (SMCGES) cohort. Self-reported facial photo-ageing data by the same 1081 participants were also collated and the two sets of data are compared.

RESULTS

Here, we found that self-reported signs of photo-ageing are concordant with photo-ageing detected by our assessors. This finding is consistent whether photo-ageing is evaluated through studying wrinkle variations (Spearman's rank correlation (ρ) value: 0.246-0.329) or through studying dyspigmentation patterns (Spearman's rank correlation (ρ) value 0.203-0.278). When studying individual wrinkles, both participants and assessors often detect the presence of the same wrinkle (Spearman's rank correlation (ρ) value 0.249-0.366). A weak-to-fair level of agreement between both participants and assessors (Cohen's kappa (κ) values: 0.041-0.233) persists and is statistically significant after accounting for agreements due to chance. Both the participant and the assessor are largely consistent in evaluating the extent of photo-ageing (area under curve (AUC) values 0.689-0.769) and in discerning between the presence or absence of a given facial wrinkle (area under curve (AUC) values 0.601-0.856).

CONCLUSION

When we analyse the overall appearance of the face, our results show that signs of photo-ageing identified by the participant are concordant with signs of photo-ageing identified by our assessors. When we focused our analysis on specific areas of the face, we found that participants were more likely to identify and self-report the same wrinkles that our assessors have also detected. Here, we found that self-reported signs of skin ageing provide a satisfactory approximation to the signs of skin ageing identified by our assessors. The ability to use self-reported signs of skin ageing should also be evaluated on scales beyond the ones discussed in this study. Currently, there are not as many photo-numeric scales for quantifying dyspigmentation patterns as there are for quantifying wrinkle variations. As Chinese skin is known to become dyspigmented more easily with age, more photo-numeric scales need to be developed and properly validated.

Large-scale genome-wide association study of 398,238 women unveils seven novel loci associated with high-grade serous epithelial ovarian cancer risk

Background

Nineteen genomic regions have been associated with high-grade serous ovarian cancer (HGSOC). We used data from the Ovarian Cancer Association Consortium (OCAC), Consortium of Investigators of Modifiers of BRCA1/BRCA2 (CIMBA), UK Biobank (UKBB), and FinnGen to identify novel HGSOC susceptibility loci and develop polygenic scores (PGS).

Methods

We analyzed >22 million variants for 398,238 women. Associations were assessed separately by consortium and meta-analysed. OCAC and CIMBA data were used to develop PGS which were trained on FinnGen data and validated in UKBB and BioBank Japan.

Results

Eight novel variants were associated with HGSOC risk. An interesting discovery biologically was finding that TP53 3'-UTR SNP rs78378222 was associated with HGSOC (per T allele relative risk (RR)=1.44, 95%CI:1.28-1.62, P=1.76×10-9). The optimal PGS included 64,518 variants and was associated with an odds ratio of 1.46 (95%CI:1.37-1.54) per standard deviation in the UKBB validation (AUROC curve=0.61, 95%CI:0.59-0.62).

Conclusions

This study represents the largest GWAS for HGSOC to date. The results highlight that improvements in imputation reference panels and increased sample sizes can identify HGSOC associated variants that previously went undetected, resulting in improved PGS. The use of updated PGS in cancer risk prediction algorithms will then improve personalized risk prediction for HGSOC.

Polygenic risk scores, radiation treatment exposures and subsequent cancer risk in childhood cancer survivors

Survivors of childhood cancer are at increased risk for subsequent cancers attributable to the late effects of radiotherapy and other treatment exposures; thus, further understanding of the impact of genetic predisposition on risk is needed. Combining genotype data for 11,220 5-year survivors from the Childhood Cancer Survivor Study and the St Jude Lifetime Cohort, we found that cancer-specific polygenic risk scores (PRSs) derived from general population, genome-wide association study, cancer loci identified survivors of European ancestry at increased risk of subsequent basal cell carcinoma (odds ratio per s.d. of the PRS: OR = 1.37, 95% confidence interval (CI) = 1.29-1.46), female breast cancer (OR = 1.42, 95% CI = 1.27-1.58), thyroid cancer (OR = 1.48, 95% CI = 1.31-1.67), squamous cell carcinoma (OR = 1.20, 95% CI = 1.00-1.44) and melanoma (OR = 1.60, 95% CI = 1.31-1.96); however, the association for colorectal cancer was not significant (OR = 1.19, 95% CI = 0.94-1.52). An investigation of joint associations between PRSs and radiotherapy found more than additive increased risks of basal cell carcinoma, and breast and thyroid cancers. For survivors with radiotherapy exposure, the cumulative incidence of subsequent cancer by age 50 years was increased for those with high versus low PRS. These findings suggest a degree of shared genetic etiology for these malignancy types in the general population and survivors, which remains evident in the context of strong radiotherapy-related risk.

Risk Factors and Innovations in Risk Assessment for Melanoma, Basal Cell Carcinoma, and Squamous Cell Carcinoma

Skin cancer is the most frequently diagnosed cancer globally and is preventable. Various risk factors contribute to different types of skin cancer, including melanoma, basal cell carcinoma, and squamous cell carcinoma. These risk factors encompass both extrinsic, such as UV exposure and behavioral components, and intrinsic factors, especially involving genetic predisposition. However, the specific risk factors vary among the skin cancer types, highlighting the importance of precise knowledge to facilitate appropriate early diagnosis and treatment for at-risk individuals. Better understanding of the individual risk factors has led to the development of risk scores, allowing the identification of individuals at particularly high risk. These advances contribute to improved prevention strategies, emphasizing the commitment to mitigating the impact of skin cancer.

Comparisons between eyebags, droopy eyelids, and eyebrow positioning identified by photo-numeric scales or identified by written descriptive scales: Insights from the Singapore/Malaysia cross-sectional genetics epidemiology study (SMCGES) cohort

BACKGROUND

We evaluate skin sagging phenotypes (eyebags, droopy eyelids, low eyebrow positioning) using written descriptive scales and photo-numeric scales. We also study how anti-ageing interventions and digital screen time influence skin sagging.

AIM

We compare the two phenotype assessment methods with each other.

METHOD

Skin sagging and personal lifestyle data obtained from 2885 ethnic Chinese young adults from the Singapore/Malaysia cross-sectional genetics epidemiology study (SMCGES) cohort were collated and compared.

RESULTS

Significant correlations (p-value < 0.001) between written descriptive scales and photo-numeric scales were observed for eyebags (0.25) and eyebrow positioning (0.08). Significant correlations (p-value < 0.001) were observed after combining both scales for eyebags (0.38), droopy eyelids (0.30), and eyebrow positioning (0.30). Anti-ageing interventions are associated with delayed progression of eyebags from 18-45 years old, droopy eyelids from 31-45 years old, and eyebrow positioning from 35-40 years old. Significantly lower (p-value < 0.02) eyebrow positioning is associated with both <1 and 1-3 h of screen time stratified by age.

CONCLUSION

Written descriptive scales provide comparable results to photo-numeric scales. However, validating and adapting photo-numeric scales for different populations identifies phenotypes better. Anti-ageing interventions are beneficial at different age ranges. Screen time is associated with skin sagging in young (18-30 years old) participants.

Multi-ancestry genome-wide meta-analysis identifies novel basal cell carcinoma loci and shared genetic effects with squamous cell carcinoma

Basal cell carcinoma (BCC) is one of the most common malignancies worldwide, yet its genetic determinants are incompletely defined. We perform a European ancestry genome-wide association (GWA) meta-analysis and a Hispanic/Latino ancestry GWA meta-analysis and meta-analyze both in a multi-ancestry GWAS meta-analysis of BCC, totaling 50,531 BCC cases and 762,234 controls from four cohorts (GERA, Mass-General Brigham Biobank, UK Biobank, and 23andMe research cohort). Here we identify 122 BCC-associated loci, of which 36 were novel, and subsequently fine-mapped these associations. We also identify an association of the well-known pigment gene SLC45A2 as well as associations at RCC2 and CLPTM1L with BCC in Hispanic/Latinos. We examine these BCC loci for association with cutaneous squamous cell carcinoma (cSCC) in 16,407 SCC cases and 762,486 controls of European ancestry, and 33 SNPs show evidence of association. Our study findings provide important insights into the genetic basis of BCC and cSCC susceptibility.

Leveraging Single-Cell Populations to Uncover the Genetic Basis of Complex Traits

The ease and throughput of single-cell genomics have steadily improved, and its current trajectory suggests that surveying single-cell populations will become routine. We discuss the merger of quantitative genetics with single-cell genomics and emphasize how this synergizes with advantages intrinsic to plants. Single-cell population genomics provides increased detection resolution when mapping variants that control molecular traits, including gene expression or chromatin accessibility. Additionally, single-cell population genomics reveals the cell types in which variants act and, when combined with organism-level phenotype measurements, unveils which cellular contexts impact higher-order traits. Emerging technologies, notably multiomics, can facilitate the measurement of both genetic changes and genomic traits in single cells, enabling single-cell genetic experiments. The implementation of single-cell genetics will advance the investigation of the genetic architecture of complex molecular traits and provide new experimental paradigms to study eukaryotic genetics.

Systematic review on gene-sun exposure interactions in skin cancer

BACKGROUND

The risk of skin cancer is determined by environmental factors like ultraviolet radiation (UVR), personal habits like time spent outdoors and genetic factors. This review aimed to survey existing studies in gene-environment (GxE) interaction on skin cancer risk, and report on GxE effect estimates.

METHODS

We searched Embase, Medline (Ovid) and Web of Science (Core Collection) and included only primary research that reported on GxE on the risk of the three most common types of skin cancer: basal cell carcinoma (BCC), squamous cell carcinoma (SCC) and melanoma. Quality assessment followed the Newcastle-Ottawa Scale. Meta-analysis was not possible because no two studies examined the same interaction. This review was registered on PROSPERO (CRD42021238064).

RESULTS

In total 260 records were identified after exclusion of duplicates. Fifteen studies were included in the final synthesis-12 used candidate gene approach. We found some evidence of GxE interactions with sun exposure, notably, with MC1R, CAT and NOS1 genes in melanoma, HAL and IL23A in BCC and HAL and XRCC1 in SCC.

CONCLUSION

Sun exposure seems to interact with genes involved in pigmentation, oxidative stress and immunosuppression, indicating that excessive UV exposure might exhaust oxidative defence and repair systems differentially, dependent on genetic make-up. Further research is warranted to better understand skin cancer epidemiology and develop sun exposure recommendations. A genome-wide approach is recommended as it might uncover unknown disease pathways dependent on UV radiation.

Design and quality control of large-scale two-sample Mendelian randomization studies

Background

Mendelian randomization (MR) studies are susceptible to metadata errors (e.g. incorrect specification of the effect allele column) and other analytical issues that can introduce substantial bias into analyses. We developed a quality control (QC) pipeline for the Fatty Acids in Cancer Mendelian Randomization Collaboration (FAMRC) that can be used to identify and correct for such errors.

Methods

We collated summary association statistics from fatty acid and cancer genome-wide association studies (GWAS) and subjected the collated data to a comprehensive QC pipeline. We identified metadata errors through comparison of study-specific statistics to external reference data sets (the National Human Genome Research Institute-European Bioinformatics Institute GWAS catalogue and 1000 genome super populations) and other analytical issues through comparison of reported to expected genetic effect sizes. Comparisons were based on three sets of genetic variants: (i) GWAS hits for fatty acids, (ii) GWAS hits for cancer and (iii) a 1000 genomes reference set.

Results

We collated summary data from 6 fatty acid and 54 cancer GWAS. Metadata errors and analytical issues with the potential to introduce substantial bias were identified in seven studies (11.6%). After resolving metadata errors and analytical issues, we created a data set of 219 842 genetic associations with 90 cancer types, generated in analyses of 566 665 cancer cases and 1 622 374 controls.

Conclusions

In this large MR collaboration, 11.6% of included studies were affected by a substantial metadata error or analytical issue. By increasing the integrity of collated summary data prior to their analysis, our protocol can be used to increase the reliability of downstream MR analyses. Our pipeline is available to other researchers via the CheckSumStats package (https://github.com/MRCIEU/CheckSumStats).

Advancing Cancer Research: Current Knowledge on Cutaneous Neoplasia

Skin cancers require a multidisciplinary approach. The updated guidelines introduce new insights into the management of these diseases. Melanoma (MM), the third most common skin cancer, a malignant melanocytic tumor, which is classified into four major histological subtypes, continues to have the potential to be a lethal disease. The mortality-incidence ratio is higher in Eastern European countries compared to Western European countries, which shows the need for better prevention and early detection in Eastern European countries. Basal cell carcinoma (BCC) and squamous cell carcinoma (cSCC) remain the top two skin cancers, and their incidence continues to grow. The gold standard in establishing the diagnosis and establishing the histopathological subtype in BCC and SCC is a skin biopsy. Sebaceous carcinoma (SeC) is an uncommon and potentially aggressive cutaneous malignancy showing sebaceous differentiation. It accounts for 0.7% of skin cancers and 3-6.7% of cancer-related deaths. Due to the rapid extension to the regional lymph nodes, SeC requires early treatment. The main treatment for sebaceous carcinoma is surgical treatment, including Mohs micrographic surgery, which has the advantage of complete margin evaluation and low recurrence rates. Primary cutaneous lymphomas (PCLs) are a heterogeneous group of lymphoproliferative diseases, with no evidence of extracutaneous determination at the moment of the diagnosis. PCLs have usually a very different evolution, prognosis, and treatment compared to the lymphomas that may secondarily involve the skin. The aim of our review is to summarize the important changes in the approach to treating melanoma, non-melanoma skin, cutaneous T and B cell lymphomas, and other types of skin cancers. For all skin cancers, optimal patient management requires a multidisciplinary approach including dermatology, medical oncology, and radiation oncology.

A multi-phenotype analysis reveals 19 susceptibility loci for basal cell carcinoma and 15 for squamous cell carcinoma

Basal cell carcinoma and squamous cell carcinoma are the most common skin cancers, and have genetic overlap with melanoma, pigmentation traits, autoimmune diseases, and blood biochemistry biomarkers. In this multi-trait genetic analysis of over 300,000 participants from Europe, Australia and the United States, we reveal 78 risk loci for basal cell carcinoma (19 previously unknown and replicated) and 69 for squamous cell carcinoma (15 previously unknown and replicated). The previously unknown risk loci are implicated in cancer development and progression (e.g. CDKL1), pigmentation (e.g. TPCN2), cardiometabolic (e.g. FADS2), and immune-regulatory pathways for innate immunity (e.g. IFIH1), and HIV-1 viral load modulation (e.g. CCR5). We also report an optimised polygenic risk score for effective risk stratification for keratinocyte cancer in the Canadian Longitudinal Study of Aging (794 cases and 18139 controls), which could facilitate skin cancer surveillance e.g. in high risk subpopulations such as transplantees.

Implementation of individualised polygenic risk score analysis: a test case of a family of four

Background

Polygenic risk scores (PRS) have been widely applied in research studies, showing how population groups can be stratified into risk categories for many common conditions. As healthcare systems consider applying PRS to keep their populations healthy, little work has been carried out demonstrating their implementation at an individual level.

Case presentation

We performed a systematic curation of PRS sources from established data repositories, selecting 15 phenotypes, comprising an excess of 37 million SNPs related to cancer, cardiovascular, metabolic and autoimmune diseases. We tested selected phenotypes using whole genome sequencing data for a family of four related individuals. Individual risk scores were given percentile values based upon reference distributions among 1000 Genomes Iberians, Europeans, or all samples. Over 96 billion allele effects were calculated in order to obtain the PRS for each of the individuals analysed here.

Conclusions

Our results highlight the need for further standardisation in the way PRS are developed and shared, the importance of individual risk assessment rather than the assumption of inherited averages, and the challenges currently posed when translating PRS into risk metrics.

A systematic review of skin ageing genes: gene pleiotropy and genes on the chromosomal band 16q24.3 may drive skin ageing

Skin ageing is the result of intrinsic genetic and extrinsic lifestyle factors. However, there is no consensus on skin ageing phenotypes and ways to quantify them. In this systematic review, we first carefully identified 56 skin ageing phenotypes from multiple literature sources and sought the best photo-numeric grading scales to evaluate them. Next, we conducted a systematic review on all 44 Genome-wide Association Studies (GWAS) on skin ageing published to date and identified genetic risk factors (2349 SNPs and 366 genes) associated with skin ageing. We identified 19 promising SNPs found to be significantly (p-Value < 1E-05) associated with skin ageing phenotypes in two or more independent studies. Here we show, using enrichment analyses strategies and gene expression data, that (1) pleiotropy is a recurring theme among skin ageing genes, (2) SNPs associated with skin ageing phenotypes are mostly located in a small handful of 44 pleiotropic and hub genes (mostly on the chromosome band 16q24.3) and 32 skin colour genes. Since numerous genes on the chromosome band 16q24.3 and skin colour genes show pleiotropy, we propose that (1) genes traditionally identified to contribute to skin colour have more than just skin pigmentation roles, and (2) further progress towards understand the development of skin pigmentation requires understanding the contributions of genes on the chromosomal band 16q24.3. We anticipate our systematic review to serve as a hub to locate primary literature sources pertaining to the genetics of skin ageing and to be a starting point for more sophisticated work examining pleiotropic genes, hub genes, and skin ageing phenotypes.

Polygenic risk score as a determinant of risk of keratinocyte cancer in an Australian population-based cohort

BACKGROUND

Keratinocyte cancer (KC) risk is determined by genetic and environmental factors. Genetic risk can be quantified by polygenic risk scores (PRS), which sum the combined effects of single nucleotide polymorphisms (SNPs).

OBJECTIVES

Our objective here was to evaluate the contribution of the summed genetic score to predict the KC risk in the phenotypically well-characterised Nambour population.

METHODS

We used PLINK v1.90 to calculate PRS for 432 cases, 566 controls, using 78 genome-wide independent SNPs that are associated with KC risk. We assessed the association between PRS and KC using logistic regression, stratifying the cohort into 3 risk groups (high 20%, intermediate 60%, low 20%).

RESULTS

The fully adjusted model including traditional risk factors (phenotypic and sun exposure-related), showed a significant 50% increase in odds of KC per standard deviation of PRS (odds ratio (OR) =1.51; 95% confidence interval (CI) =1.30-1.76, P=5.75 × 10^-8 ). Those in the top 20% PRS had over three times the risk of KC of those in the lowest 20% (OR=3.45; 95% CI=2.18-5.50, P=1.5×10^-7 ) and higher absolute risk of KC per 100 person-years of 2.96 compared with 1.34. Area under the ROC curve increased from 0.72 to 0.74 on adding PRS to the fully adjusted model.

CONCLUSIONS

These results show that PRS can enhance the prediction of KC above traditional risk factors.

Single-cell analysis of human basal cell carcinoma reveals novel regulators of tumor growth and the tumor microenvironment

How basal cell carcinoma (BCC) interacts with its tumor microenvironment to promote growth is unclear. We use singe-cell RNA sequencing to define the human BCC ecosystem and discriminate between normal and malignant epithelial cells. We identify spatial biomarkers of tumors and their surrounding stroma that reinforce the heterogeneity of each tissue type. Combining pseudotime, RNA velocity-PAGA, cellular entropy, and regulon analysis in stromal cells reveals a cancer-specific rewiring of fibroblasts, where STAT1, TGF-β, and inflammatory signals induce a noncanonical WNT5A program that maintains the stromal inflammatory state. Cell-cell communication modeling suggests that tumors respond to the sudden burst of fibroblast-specific inflammatory signaling pathways by producing heat shock proteins, whose expression we validated in situ. Last, dose-dependent treatment with an HSP70 inhibitor suppresses in vitro vismodegib-resistant BCC cell growth, Hedgehog signaling, and in vivo tumor growth in a BCC mouse model, validating HSP70's essential role in tumor growth and reinforcing the critical nature of tumor microenvironment cross-talk in BCC progression.

Genetic polymorphism in Methylenetetrahydrofolate Reductase chloride transport protein 6 (MTHFR CLCN6) gene is associated with keratinocyte skin cancer in a cohort of renal transplant recipients

Background

Renal transplant recipients (RTRs) are at increased risk of keratinocyte cancer (KC), especially cutaneous squamous cell carcinoma (cSCC). Previous studies identified a genetic variant of the Methylenetetrahydrofolate Reductase (MTHFR) gene, C677T, which conferred a risk for diagnosis of cSCC in Irish RTRs.

Objective

We sought to find further genetic variation in MTHFR and overlap genes that may be associated with a diagnosis of KC in RTRs.

Methods

Genotyping of a combined RTR population (n = 821) from two centres, Ireland (n = 546) and the USA (n = 275), was performed. This included 290 RTRs with KC and 444 without. Eleven single nucleotide polymorphisms (SNPs) in the MTHFR gene and seven in the overlap gene MTHFR Chloride transport protein 6 (CLCN6) were evaluated and association explored by time to event analysis (from transplant to first KC) using Cox proportional hazards model.

Results

Polymorphism at MTHFR CLCN6 (rs9651118) was significantly associated with KC in RTRs (HR 1.50, 95% CI 1.17–1.91, p < 0.00061) and cSCC (HR 1.63, 95% CI 1.14–2.34, p = 0.007). A separate SNP, MTHFR C677T, was also significantly associated with KC in the Irish population (HR 1.31, 95% CI 1.05–1.63, p = 0.016), but not American RTRs.

Conclusions

We report the association of a SNP in the MTHFR overlap gene, CLCN6 and KC in a combined RTR population. While the exact function of CLCN6 is not known, it is proposed to be involved in folate availability. Future applications could include incorporation in a polygenic risk score for KC in RTRs to help identify those at increased risk beyond traditional risk factor assessment.

Genome-wide association study of actinic keratosis identifies new susceptibility loci implicated in pigmentation and immune regulation pathways

Actinic keratosis (AK) is a common precancerous cutaneous neoplasm that arises on chronically sun-exposed skin. AK susceptibility has a moderate genetic component, and although a few susceptibility loci have been identified, including IRF4, TYR, and MC1R, additional loci have yet to be discovered. We conducted a genome-wide association study of AK in non-Hispanic white participants of the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort (n = 63,110, discovery cohort), with validation in the Mass-General Brigham (MGB) Biobank cohort (n = 29,130). We identified eleven loci (P < 5 × 10-8), including seven novel loci, of which four novel loci were validated. In a meta-analysis (GERA + MGB), one additional novel locus, TRPS1, was identified. Genes within the identified loci are implicated in pigmentation (SLC45A2, IRF4, BNC2, TYR, DEF8, RALY, HERC2, and TRPS1), immune regulation (FOXP1 and HLA-DQA1), and cell signaling and tissue remodeling (MMP24) pathways. Our findings provide novel insight into the genetics and pathogenesis of AK susceptibility.

The Immunogenetics of Non-melanoma Skin Cancer

Non-melanoma skin cancer (NMSC) is the most common malignancy seen in Caucasians and includes basal cell carcinoma (BCC) and squamous cell carcinoma (SCC). The incidence of NMSC is showing an increasing trend which is attributed to the increased use of sunbeds, recreational sun exposure, aging population, and partly to improved screening and reporting. Ultraviolet (UV) radiation plays the most crucial role in the pathogenesis of both BCC and SCC by inducing DNA damage and mutagenic photoproducts. Other risk factors are fair skin, old age, genetic predisposition, immunosuppression, ionizing radiation, organic chemicals, and HPV infection. The role of genomic instability, genetic mutations/aberrations, and host immunity has been fairly illustrated in several studies. This chapter aims to discuss these aspects of NMSC in detail.

Polygenic Risk Scores Stratify Keratinocyte Cancer Risk among Solid Organ Transplant Recipients with Chronic Immunosuppression in a High Ultraviolet Radiation Environment

Solid organ transplant recipients (SOTRs) have elevated risks for basal cell carcinoma (BCC) and squamous cell carcinoma (SCC), especially in high UVR environments. We assessed whether polygenic risk scores can improve the prediction of BCC and SCC risks and multiplicity over and above the traditional risk factors in SOTRs in a high UV setting. We built polygenic risk scores for BCC (n = 594,881) and SCC (n = 581,431) using UK Biobank and 23andMe datasets, validated them in the Australian QSkin Sun and Health Study cohort (n > 6,300), and applied them in SOTRs in the skin tumor in allograft recipients cohort from Queensland, Australia, a high UV environment. About half of the SOTRs with a high genetic risk developed BCC (absolute risk = 45.45%, 95% confidence interval = 33.14-58.19%) and SCC (absolute risk = 44.12%, 95% confidence interval = 32.08-56.68%). For both cancers, SOTRs in the top quintile were at >3-fold increased risk relative to those in the bottom quintile. The respective polygenic risk scores improved risk predictions by 2% for BCC (area under the curve = 0.77 vs. 0.75, P = 0.0691) and SCC (area under the curve = 0.84 vs. 0.82, P = 0.0260), over and above the established risk factors, and 19.03% (for BCC) and 18.10% (for SCC) of the SOTRs were reclassified in a high/medium/low risk scenario. The polygenic risk scores also added predictive accuracy for tumor multiplicity (BCC R2 = 0.21 vs. 0.19, P = 3.2 × 10-3; SCC R2 = 0.30 vs. 0.27, P = 4.6 × 10-4).

Impact of FOXP1 rs2687201 genetic variant on the susceptibility to HCV-related hepatocellular carcinoma in Egyptians

Hepatocellular carcinoma (HCC) constitutes a challenging health problem in Egypt due to the high incidence of hepatitis C virus (HCV) infection. Improved understanding of genetic mechanisms underlying the individual predisposition to HCC will lead to enhancements in the early diagnosis, treatment, and prevention of this disease. Transcription factor forkhead box P1 (FOXP1) is involved in the cellular processes of proliferation, differentiation, metabolism, and longevity. In addition, it has been implicated in hepatic tumorigenesis. The present study explored the association of C/A single-nucleotide polymorphism in the FOXP1 gene (rs2687201) with HCC susceptibility in HCV Egyptian patients. The study included 108 patients with HCV-dependant HCC, 86 HCV patients, and 80- age and gender-matched healthy controls. rs2687201 genotyping was performed by allelic discrimination method using TaqMan real-time PCR assays while FOXP1 gene expression and protein level were determined using qRT-PCR and enzyme-linked immunoassay, respectively. Our results revealed a significant association between FOXP1 rs2687201 and HCC risk where (A) allele was significantly more frequent in patients with HCC compared to controls (odds ratio [OR]: 1.88, 95% confidence interval [CI]: 1.17-3.04, p = 0.01) and to HCV patients (OR: 1.85, 95% CI: 1.62-2.94, p = 0.012). Furthermore, FOXP1 gene and protein expression levels were remarkably higher in (CA + AA) than in CC genotype carriers in a dominant model. The (CA + AA) genotype displayed a significantly shorter overall survival than the CC genotype in HCC patients. In conclusion, FOXP1 gene polymorphism rs2687201 is significantly associated with HCC, but not with HCV infection, in Egyptian patients.

Polyunsaturated Fatty Acid Levels and the Risk of Keratinocyte Cancer: A Mendelian Randomization Analysis

BACKGROUND

Keratinocyte cancer is the commonest cancer, imposing a high economic burden on the health care system. Observational studies have shown mixed associations between polyunsaturated fatty acids (PUFA) and keratinocyte cancer, basal cell carcinoma (BCC), and squamous cell carcinoma (SCC). We explored whether genetically predicted PUFA levels are associated with BCC and SCC risks.

METHODS

We conducted a two-sample Mendelian randomization study using PUFA level genome-wide association studies (GWAS) from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (n > 8,000), and the meta-analysis GWASs from UKB, 23andMe, and Qskin for BCC (n = 651,138) and SCC (n = 635,331) risk.

RESULTS

One SD increase in genetically predicted levels of linoleic acid [OR = 0.94, 95% confidence interval (CI) = 0.91-0.97, P = 1.4 × 10^-4] and alpha-linolenic acid (OR = 0.91, 95% CI = 0.86-0.96, P = 5.1 × 10^-4) was associated with a reduced BCC risk, while arachidonic acid (OR = 1.04, 95% CI = 1.02-1.06, P = 3.2 × 10^-4) and eicosapentaenoic acid (OR = 1.10, 95% CI = 1.04-1.16, P = 1.5 × 10^-3) were associated with an increased BCC risk.

CONCLUSIONS

Higher genetically predicted levels of linoleic acid and alpha-linolenic acid were associated with a reduced BCC risk, but arachidonic acid and eicosapentaenoic acid were associated with a higher BCC risk.

IMPACT

PUFA-related diet and supplementation could influence BCC etiology.

Review of the Molecular Genetics of Basal Cell Carcinoma; Inherited Susceptibility, Somatic Mutations, and Targeted Therapeutics

Simple Summary

Basal cell carcinoma is the most common human cancer worldwide. The molecular basis of BCC involves an interplay of inherited genetic susceptibility and somatic mutations, commonly induced by exposure to UV radiation. In this review, we outline the currently known germline and somatic mutations implicated in the pathogenesis of BCC with particular attention paid toward affected molecular pathways. We also discuss polymorphisms and associated phenotypic traits in addition to active areas of BCC research. We finally provide a brief overview of existing non-surgical treatments and emerging targeted therapeutics for BCC such as Hedgehog pathway inhibitors, immune modulators, and histone deacetylase inhibitors.

Abstract

Basal cell carcinoma (BCC) is a significant public health concern, with more than 3 million cases occurring each year in the United States, and with an increasing incidence. The molecular basis of BCC is complex, involving an interplay of inherited genetic susceptibility, including single nucleotide polymorphisms and genetic syndromes, and sporadic somatic mutations, often induced by carcinogenic exposure to UV radiation. This review outlines the currently known germline and somatic mutations implicated in the pathogenesis of BCC, including the key molecular pathways affected by these mutations, which drive oncogenesis. With advances in next generation sequencing and our understanding of the molecular genetics of BCC, established and emerging targeted therapeutics are offering new avenues for the non-surgical treatment of BCC. These agents, including Hedgehog pathway inhibitors, immune modulators, and histone deacetylase inhibitors, will also be discussed.

Citrus Consumption and Risk of Non-Melanoma Skin Cancer in the UK Biobank

Background: Non-melanoma skin cancer (NMSC) incidence has been dramatically increasing worldwide. Psoralen, a known photocarcinogen, is naturally abundant in citrus products, leading to the hypothesis that high citrus consumption may increase NMSC risk.Methods: We fitted age- and multivariable-adjusted logistic regression models to evaluate the association between citrus consumption and NMSC risk among 197,372 UKBB participants. A total of 9,613 NMSC cases were identified using International Classification of Disease 10 codes. Citrus consumption data were collected via five rounds of 24-hour recall questionnaires.Results: We found no association between high total citrus consumption and NMSC risk, although a slightly elevated NMSC risk was observed among participants who consumed >0 to half a serving of total citrus per day (OR [95% CI] = 1.08 [1.01-1.16]). There was no association between individual citrus products and NMSC risk.Conclusion: High citrus consumption was not associated with an increased risk of NMSC in our UKBB sample. Further studies are needed to clarify these associations.Supplemental data for this article is available online at https://doi.org/10.1080/01635581.2021.1952439 .

Characterization of intermediate-sized insertions using whole-genome sequencing data and analysis of their functional impact on gene expression

Intermediate-sized insertions are one of the structural variants contributing to genome diversity. However, due to technical difficulties in identifying them, their importance in disease pathogenicity and gene expression regulation remains unclear. We used whole-genome sequencing data of 174 Japanese samples to characterize intermediate-sized insertions using a highly-accurate insertion calling method (IMSindel software and joint-call recovery) and obtained a catalogue of 4254 insertions. We constructed an imputation panel comprising of insertions and SNVs from all samples, and conducted imputation of intermediate-sized insertions for 82 publicly-available Japanese samples. Positive Predictive Value of imputation, evaluated using Nanopore long-read sequencing data, was 97%. Subsequent eQTL analysis predicted 128 (~ 3.0%) insertions as causative for gene expression level changes. Enrichment analysis of causal insertions for genome regulatory elements showed significant associations with CTCF-binding sites, super-enhancers, and promoters. Among 17 causal insertions found in the same causal set with GWAS hits, there were insertions associated with changes in expression of cancer-related genes such as BRCA1, ZNF222, and ABCB10. Analysis of insertions sequences revealed that 461 insertions were short tandem duplications frequently found in early-replicating regions of genome. Furthermore, comparison of functional importance of intermediate-sized insertions with that of intermediate-sized deletions detected in the same sample set in our previous study showed that insertions were more frequent in genic regions, and proportion of functional candidates was smaller in insertions. Here, we characterize a high-confidence set of intermediate-sized insertions and indicate their importance in gene expression regulation. Our results emphasize the importance of considering intermediate-sized insertions in trait association studies.

Epigenetic Alterations in Keratinocyte Carcinoma

Basal cell carcinoma (BCC) and squamous cell carcinoma (SCC) are both derived from epidermal keratinocytes but are phenotypically diverse. To improve the understanding of keratinocyte carcinogenesis, it is critical to understand epigenetic alterations, especially those that govern gene expression. We examined changes to the enhancer-associated histone acetylation mark H3K27ac by mapping matched tumor-normal pairs from 11 patients (five with BCC and six with SCC) undergoing Mohs surgery. Our analysis uncovered cancer-specific enhancers on the basis of differential H3K27ac peaks between matched tumor-normal pairs. We also uncovered biological pathways potentially altered in keratinocyte carcinoma, including enriched epidermal development and Wnt signaling pathways enriched in BCCs and enriched immune response and cell activation pathways in SCCs. We also observed enrichment of transcription factors that implicated SMAD and JDP2 in BCC pathogenesis and FOXP1 in SCC pathogenesis. On the basis of these findings, we prioritized three loci with putative regulation events (FGFR2 enhancer in BCC, intragenic regulation of FOXP1 in SCC, and WNT5A promoter in both subtypes) and validated our findings with published gene expression data. Our findings highlight unique and shared epigenetic alterations in histone modifications and potential regulators for BCCs and SCCs that likely impact the divergent oncogenic pathways, paving the way for targeted drug discoveries.

Common genetic polymorphisms contribute to the association between chronic lymphocytic leukaemia and non-melanoma skin cancer

BACKGROUND

Epidemiological studies have demonstrated a positive association between chronic lymphocytic leukaemia (CLL) and non-melanoma skin cancer (NMSC). We hypothesized that shared genetic risk factors between CLL and NMSC could contribute to the association observed between these diseases.

METHODS

We examined the association between (i) established NMSC susceptibility loci and CLL risk in a meta-analysis including 3100 CLL cases and 7667 controls and (ii) established CLL loci and NMSC risk in a study of 4242 basal cell carcinoma (BCC) cases, 825 squamous cell carcinoma (SCC) cases and 12802 controls. Polygenic risk scores (PRS) for CLL, BCC and SCC were constructed using established loci. Logistic regression was used to estimate odds ratios (ORs) and 95% confidence intervals (CIs).

RESULTS

Higher CLL-PRS was associated with increased BCC risk (OR4th-quartile-vs-1st-quartile = 1.13, 95% CI: 1.02-1.24, Ptrend = 0.009), even after removing the shared 6p25.3 locus. No association was observed with BCC-PRS and CLL risk (Ptrend = 0.68). These findings support a contributory role for CLL in BCC risk, but not for BCC in CLL risk. Increased CLL risk was observed with higher SCC-PRS (OR4th-quartile-vs-1st-quartile = 1.22, 95% CI: 1.08-1.38, Ptrend = 1.36 × 10-5), which was driven by shared genetic susceptibility at the 6p25.3 locus.

CONCLUSION

These findings highlight the role of pleiotropy regarding the pathogenesis of CLL and NMSC and shows that a single pleiotropic locus, 6p25.3, drives the observed association between genetic susceptibility to SCC and increased CLL risk. The study also provides evidence that genetic susceptibility for CLL increases BCC risk.

Genome-wide association study in almost 195,000 individuals identifies 50 previously unidentified genetic loci for eye color

Human eye color is highly heritable, but its genetic architecture is not yet fully understood. We report the results of the largest genome-wide association study for eye color to date, involving up to 192,986 European participants from 10 populations. We identify 124 independent associations arising from 61 discrete genomic regions, including 50 previously unidentified. We find evidence for genes involved in melanin pigmentation, but we also find associations with genes involved in iris morphology and structure. Further analyses in 1636 Asian participants from two populations suggest that iris pigmentation variation in Asians is genetically similar to Europeans, albeit with smaller effect sizes. Our findings collectively explain 53.2% (95% confidence interval, 45.4 to 61.0%) of eye color variation using common single-nucleotide polymorphisms. Overall, our study outcomes demonstrate that the genetic complexity of human eye color considerably exceeds previous knowledge and expectations, highlighting eye color as a genetically highly complex human trait.

Loss-of-Function Variants in the Tumor-Suppressor Gene PTPN14 Confer Increased Cancer Risk

The success of genome-wide association studies (GWAS) in identifying common, low-penetrance variant-cancer associations for the past decade is undisputed. However, discovering additional high-penetrance cancer mutations in unknown cancer predisposing genes requires detection of variant-cancer association of ultra-rare coding variants. Consequently, large-scale next-generation sequence data with associated phenotype information are needed. Here, we used genotype data on 166,281 Icelanders, of which, 49,708 were whole-genome sequenced and 408,595 individuals from the UK Biobank, of which, 41,147 were whole-exome sequenced, to test for association between loss-of-function burden in autosomal genes and basal cell carcinoma (BCC), the most common cancer in Caucasians. A total of 25,205 BCC cases and 683,058 controls were tested. Rare germline loss-of-function variants in PTPN14 conferred substantial risks of BCC (OR, 8.0; P = 1.9 × 10^-12), with a quarter of carriers getting BCC before age 70 and over half in their lifetime. Furthermore, common variants at the PTPN14 locus were associated with BCC, suggesting PTPN14 as a new, high-impact BCC predisposition gene. A follow-up investigation of 24 cancers and three benign tumor types showed that PTPN14 loss-of-function variants are associated with high risk of cervical cancer (OR, 12.7, P = 1.6 × 10^-4) and low age at diagnosis. Our findings, using power-increasing methods with high-quality rare variant genotypes, highlight future prospects for new discoveries on carcinogenesis. SIGNIFICANCE: This study identifies the tumor-suppressor gene PTPN14 as a high-impact BCC predisposition gene and indicates that inactivation of PTPN14 by germline sequence variants may also lead to increased risk of cervical cancer.

Genetic and functional interaction network analysis reveals global enrichment of regulatory T cell genes influencing basal cell carcinoma susceptibility

BACKGROUND

Basal cell carcinoma (BCC) of the skin is the most common form of human cancer, with more than 90% of tumours presenting with clear genetic activation of the Hedgehog pathway. However, polygenic risk factors affecting mechanisms such as DNA repair and cell cycle checkpoints or which modulate the tumour microenvironment or host immune system play significant roles in determining whether genetic mutations culminate in BCC development. We set out to define background genetic factors that play a role in influencing BCC susceptibility via promoting or suppressing the effects of oncogenic drivers of BCC.

METHODS

We performed genome-wide association studies (GWAS) on 17,416 cases and 375,455 controls. We subsequently performed statistical analysis by integrating data from population-based genetic studies of multi-omics data, including blood- and skin-specific expression quantitative trait loci and methylation quantitative trait loci, thereby defining a list of functionally relevant candidate BCC susceptibility genes from our GWAS loci. We also constructed a local GWAS functional interaction network (consisting of GWAS nearest genes) and another functional interaction network, consisting specifically of candidate BCC susceptibility genes.

RESULTS

A total of 71 GWAS loci and 46 functional candidate BCC susceptibility genes were identified. Increased risk of BCC was associated with the decreased expression of 26 susceptibility genes and increased expression of 20 susceptibility genes. Pathway analysis of the functional candidate gene regulatory network revealed strong enrichment for cell cycle, cell death, and immune regulation processes, with a global enrichment of genes and proteins linked to TReg cell biology.

CONCLUSIONS

Our genome-wide association analyses and functional interaction network analysis reveal an enrichment of risk variants that function in an immunosuppressive regulatory network, likely hindering cancer immune surveillance and effective antitumour immunity.

Characterization of Hepatitis B Virus Integrations Identified in Hepatocellular Carcinoma Genomes

Hepatocellular carcinoma (HCC) is a leading cause of cancer-related mortality. Almost half of HCC cases are associated with hepatitis B virus (HBV) infections, which often lead to HBV sequence integrations in the human genome. Accurate identification of HBV integration sites at a single nucleotide resolution is critical for developing a better understanding of the cancer genome landscape and of the disease itself. Here, we performed further analyses and characterization of HBV integrations identified by our recently reported VIcaller platform in recurrent or known HCC genes (such as TERT, MLL4, and CCNE1) as well as non-recurrent cancer-related genes (such as CSMD2, NKD2, and RHOU). Our pathway enrichment analysis revealed multiple pathways involving the alcohol dehydrogenase 4 gene, such as the metabolism pathways of retinol, tyrosine, and fatty acid. Further analysis of the HBV integration sites revealed distinct patterns involving the integration upper breakpoints, integrated genome lengths, and integration allele fractions between tumor and normal tissues. Our analysis also implies that the VIcaller method has diagnostic potential through discovering novel clonal integrations in cancer-related genes. In conclusion, although VIcaller is a hypothesis free virome-wide approach, it can still be applied to accurately identify genome-wide integration events of a specific candidate virus and their integration allele fractions.

Disease risk scores for skin cancers

We trained and validated risk prediction models for the three major types of skin cancer- basal cell carcinoma (BCC), squamous cell carcinoma (SCC), and melanoma-on a cross-sectional and longitudinal dataset of 210,000 consented research participants who responded to an online survey covering personal and family history of skin cancer, skin susceptibility, and UV exposure. We developed a primary disease risk score (DRS) that combined all 32 identified genetic and non-genetic risk factors. Top percentile DRS was associated with an up to 13-fold increase (odds ratio per standard deviation increase >2.5) in the risk of developing skin cancer relative to the middle DRS percentile. To derive lifetime risk trajectories for the three skin cancers, we developed a second and age independent disease score, called DRSA. Using incident cases, we demonstrated that DRSA could be used in early detection programs for identifying high risk asymptotic individuals, and predicting when they are likely to develop skin cancer. High DRSA scores were not only associated with earlier disease diagnosis (by up to 14 years), but also with more severe and recurrent forms of skin cancer.

Programmed cell death 1 (PDCD1) gene haplotypes and susceptibility of patients to basal cell carcinoma

Programmed death-1 (PD-1), as an immunoinhibitory receptor encoded by programmed cell death-1 (PDCD1) gene, has a pivotal role in tolerance to self-antigens. Mutations of PDCD1 may participate in susceptibility to basal cell carcinoma (BCC) as the most common of skin cancer. We studied the impacts of two single nucleotide polymorphisms (SNPs) within PDCD1 and their haplotypes in BCC susceptibility in an Iranian population. The blood samples were collected from 210 BCC and 220 healthy individuals. After the extraction of genomic DNA, the genotypes and alleles of PD1.1 G/A (rs36084323) and PD1.6 G/A (rs10204525) SNPs were determined by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). Four haplotypes were estimated by these SNPs. Our data revealed that genotype and allele frequencies of PD1.1 and PD1.6 polymorphisms in BCC patients were similar to those in healthy individuals. The results of estimated haplotypes for PDCD1 indicated that GG and AA haplotypes of PDCD1 had protective effects on BCC susceptibility (OR = 0.7, 95% CI = 0.51-0.96, p = 0.03 and OR = 0.57, 95% CI = 0.35-0.91, p = 0.02, respectively), while GA and AG haplotypes served as the risk factors for developing BCC (OR = 1.76, 95% CI = 1.09-2.84, p = 0.02 and OR = 3.87, 95% CI = 1.95-7.69, p = <0.001, respectively). Based on these findings, frequency distributions of PDCD1 haplotypes have important roles in the determination of BCC development in the Iranian population. However, larger multicenter studies are required to confirm this conclusion.

Pigmentary traits, sun exposure, and risk of non-Hodgkin's lymphoma/chronic lymphocytic leukemia: A study within the French E3N prospective cohort

To investigate whether risk factors for keratinocyte carcinomas (KCs), namely pigmentary traits and sun exposure, are associated with risk of non‐Hodgkin's lymphomas (NHL) and chronic lymphocytic leukemia (CLL). E3N is a prospective cohort of French women aged 40–65 years at inclusion in 1990. Cancer data were collected at baseline and updated every 2–3 years. Hazard Ratios (HRs) and 95% confidence intervals (CIs) for associations between pigmentary traits and sun exposure, and risk of CLL/NHL were estimated using Cox models. With a median follow‐up of 24 years, 622 incident cases of CLL/NHL were ascertained among the 92,097 included women. The presence of nevi was associated with CLL/NHL risk: HR for “many or very many nevi” relative to “no nevi”: 1.56 [1.15; 2.11]. Such association with number of nevi appears to be mostly limited to risk of CLL: HR for “many or very many nevi”: 3.00 [1.38; 6.52]; versus 1.32 [0.94; 1.84] for NHL. Women whose skin was highly sensitive to sunburn also had a higher risk of CLL: HR = 1.96 [1.21; 3.18], while no increase in risk of NHL was observed. Skin or hair color, number of freckles, and average daily ultraviolet (UV) dose during spring and summer in location of residence at birth or at inclusion (kJ/m²) were not associated with CLL/NHL risk. Some pigmentary traits (presence of nevi and skin sensitivity), but not sun exposure, were associated with CLL/NHL. These observations suggest that CLL may share some constitutional risk factors with keratinocyte cancers.

Photodynamic Therapy for Basal Cell Carcinoma: The Clinical Context for Future Research Priorities

Photodynamic therapy (PDT) is an established treatment option for low-risk basal cell carcinoma (BCC). BCC is the most common human cancer and also a convenient cancer in which to study PDT treatment. This review clarifies challenges to researchers evident from the clinical use of PDT in BCC treatment. It outlines the context of PDT and how PDT treatments for BCC have been developed hitherto. The sections examine the development of systemic and subsequently topical photosensitizers, light delivery regimens, and the use of PDT in different patient populations and subtypes of BCC. The outcomes of topical PDT are discussed in comparison with alternative treatments, and topical PDT applications in combination and adjuvant therapy are considered. The intention is to summarize the clinical relevance and expose areas of research need in the BCC context, ultimately to facilitate improvements in PDT treatment.

Cancer PRSweb: An Online Repository with Polygenic Risk Scores for Major Cancer Traits and Their Evaluation in Two Independent Biobanks

To facilitate scientific collaboration on polygenic risk scores (PRSs) research, we created an extensive PRS online repository for 35 common cancer traits integrating freely available genome-wide association studies (GWASs) summary statistics from three sources: published GWASs, the NHGRI-EBI GWAS Catalog, and UK Biobank-based GWASs. Our framework condenses these summary statistics into PRSs using various approaches such as linkage disequilibrium pruning/p value thresholding (fixed or data-adaptively optimized thresholds) and penalized, genome-wide effect size weighting. We evaluated the PRSs in two biobanks: the Michigan Genomics Initiative (MGI), a longitudinal biorepository effort at Michigan Medicine, and the population-based UK Biobank (UKB). For each PRS construct, we provide measures on predictive performance and discrimination. Besides PRS evaluation, the Cancer-PRSweb platform features construct downloads and phenome-wide PRS association study results (PRS-PheWAS) for predictive PRSs. We expect this integrated platform to accelerate PRS-related cancer research.

Severe epidermolysis bullosa simplex phenotype caused by codominant mutations p.Ile377Thr in keratin 14 and p.Gly138Glu in keratin 5

Epidermolysis bullosa simplex (EBS) is a rare skin disease usually inherited in an autosomal dominant pattern. EBS is resulting from mutations in keratin 5 (KRT5) and keratin 14 (KRT14) genes encoding the keratins 5 and 14 proteins expressed in the keratinocytes of the basal layer of the epidermis. To date, seven pathogenic mutations have been reported to be responsible for EBS in the Canadian population from the province of Quebec: p.Pro25Leu, p.Leu150Pro, p.Met327Thr and p.Arg559X in KRT5; p.Arg125Ser, p.Ile377Thr and p.Ile412Phe in KRT14. Here, we present a novel French-Canadian patient diagnosed with EBS confined to the soles but presenting a severe complication form including blisters, hyperkeratosis, skin erosions and toenail abnormalities. Mutation screening was performed by direct sequencing of the entire coding regions of KRT5 and KRT14 genes and revealed the previously reported missense heterozygous mutation c. 1130T > C in KRT14 (p.Ile377Thr). Furthermore, this patient is carrying a second mutation in KRT5, c.413G > A (p.Gly138Glu), which has been linked to an increased risk of basal cell carcinoma in the literature. We suspect an impact of the p.Gly138Glu variant on the EBS phenotype severity of the studied patient. The pathogenicity and consequences of both genetic variations were simulated by in silico tools.

Delineation of clinical and biological factors associated with cutaneous squamous cell carcinoma among patients with chronic lymphocytic leukemia

BACKGROUND

The incidence of cutaneous squamous cell carcinoma (SCC) in patients with chronic lymphocytic leukemia (CLL) is significantly higher compared with age- and sex-matched controls.

OBJECTIVE

To evaluate the association of factors associated with SCC risk.

METHODS

Clinical CLL and SCC data were obtained from Mayo Clinic CLL Resource and self-reported questionnaires among patients with newly diagnosed CLL. We computed the CLL International Prognostic Index (CLL-IPI) from CLL prognostic factors, and a polygenic risk score from SCC susceptibility variants. We used Cox regression to estimate hazard ratios (HRs) and 95% confidence intervals (CIs).

RESULTS

Among 1269 patients with CLL, the median follow-up was 7 years, and SCC subsequently developed in 124 patients. Significant associations with SCC risk were history of skin cancer (HR=4.80; 95% CI: 3.37-6.83), CLL-IPI (HR=1.42; 95% CI: 1.13-1.80), and polygenic risk score (HR=2.58; 95% CI: 1.50-4.43). In a multivariable model, these factors were independent predictors (C statistic = 0.69; 95% CI: 0.62-0.76). T-cell immunosuppressive treatments were also associated with SCC risk (HR=2.29; 95% CI: 1.47-3.55; adjusted for age, sex, and prior SCC).

LIMITATIONS

The sample size decreases when combining all risk factors in a single model.

CONCLUSION

SCC risk includes history of skin cancer, an aggressive disease at time of CLL diagnosis, receiving T-cell immunosuppressive treatments, and high polygenic risk score. Future studies should develop prediction models that include these factors to improved screening guidelines.

Genetic and environmental factors underlying keratinocyte carcinoma risk

Recent large-scale GWAS and large epidemiologic studies have accelerated the discovery of genes and environmental factors that contribute to the risk of keratinocyte carcinoma (KC), which includes basal cell carcinoma (BCC) and squamous cell carcinoma (SCC). This Review summarizes the genomic regions associated with SCC and BCC risk, examines the genetic overlap between SCC and BCC, and discusses biological pathways involved in SCC and BCC development. Next, we review environmental factors that are associated with KC risk, including those that are shared between SCC and BCC as well as others that associated with only one type of KC. We conclude with a critical appraisal of current research and potential directions for future research.

Genome-wide association meta-analyses combining multiple risk phenotypes provide insights into the genetic architecture of cutaneous melanoma susceptibility

Most genetic susceptibility to cutaneous melanoma remains to be discovered. Meta-analysis genome-wide association study (GWAS) of 36,760 cases of melanoma (67% newly genotyped) and 375,188 controls identified 54 significant (P < 5 × 10-8) loci with 68 independent single nucleotide polymorphisms. Analysis of risk estimates across geographical regions and host factors suggests the acral melanoma subtype is uniquely unrelated to pigmentation. Combining this meta-analysis with GWAS of nevus count and hair color, and transcriptome association approaches, uncovered 31 potential secondary loci for a total of 85 cutaneous melanoma susceptibility loci. These findings provide insights into cutaneous melanoma genetic architecture, reinforcing the importance of nevogenesis, pigmentation and telomere maintenance, together with identifying potential new pathways for cutaneous melanoma pathogenesis.

Combined analysis of keratinocyte cancers identifies novel genome-wide loci

The keratinocyte cancers (KC), basal cell carcinoma (BCC) and squamous cell carcinoma (SCC) are the most common cancers in fair-skinned people. KC treatment represents the second highest cancer healthcare expenditure in Australia. Increasing our understanding of the genetic architecture of KC may provide new avenues for prevention and treatment. We first conducted a series of genome-wide association studies (GWAS) of KC across three European ancestry datasets from Australia, Europe and USA, and used linkage disequilibrium (LD) Score regression (LDSC) to estimate their pairwise genetic correlations. We employed a multiple-trait approach to map genes across the combined set of KC GWAS (total N = 47 742 cases, 634 413 controls). We also performed meta-analyses of BCC and SCC separately to identify trait specific loci. We found substantial genetic correlations (generally 0.5–1) between BCC and SCC suggesting overlapping genetic risk variants. The multiple trait combined KC GWAS identified 63 independent genome-wide significant loci, 29 of which were novel. Individual separate meta-analyses of BCC and SCC identified an additional 13 novel loci not found in the combined KC analysis. Three new loci were implicated using gene-based tests. New loci included common variants in BRCA2 (distinct to known rare high penetrance cancer risk variants), and in CTLA4, a target of immunotherapy in melanoma. We found shared and trait specific genetic contributions to BCC and SCC. Considering both, we identified a total of 79 independent risk loci, 45 of which are novel.

Genome-wide identification of loci associated with growth in rainbow trout

Background

Growth is a major economic production trait in aquaculture. Improvements in growth performance will reduce time and cost for fish to reach market size. However, genes underlying growth have not been fully explored in rainbow trout.

Results

A previously developed 50 K gene-transcribed SNP chip, containing ~ 21 K SNPs showing allelic imbalances potentially associated with important aquaculture production traits including body weight, muscle yield, was used for genotyping a total of 789 fish with available phenotypic data for bodyweight gain. Genotyped fish were obtained from two consecutive generations produced in the NCCCWA growth-selection breeding program. Weighted single-step GBLUP (WssGBLUP) was used to perform a genome-wide association (GWA) analysis to identify quantitative trait loci (QTL) associated with bodyweight gain. Using genomic sliding windows of 50 adjacent SNPs, 247 SNPs associated with bodyweight gain were identified. SNP-harboring genes were involved in cell growth, cell proliferation, cell cycle, lipid metabolism, proteolytic activities, chromatin modification, and developmental processes. Chromosome 14 harbored the highest number of SNPs (n = 50). An SNP window explaining the highest additive genetic variance for bodyweight gain (~ 6.4%) included a nonsynonymous SNP in a gene encoding inositol polyphosphate 5-phosphatase OCRL-1. Additionally, based on a single-marker GWA analysis, 33 SNPs were identified in association with bodyweight gain. The highest SNP explaining variation in bodyweight gain was identified in a gene coding for thrombospondin-1 (THBS1) (R² = 0.09).

Conclusion

The majority of SNP-harboring genes, including OCRL-1 and THBS1, were involved in developmental processes. Our results suggest that development-related genes are important determinants for growth and could be prioritized and used for genomic selection in breeding programs.

Predicting keratinocyte carcinoma in patients with actinic keratosis: development and internal validation of a multivariable risk-prediction model

BACKGROUND

Patients with actinic keratosis (AK) are at increased risk for developing keratinocyte carcinoma (KC) but predictive factors and their risk rates are unknown.

OBJECTIVES

To develop and internally validate a prediction model to calculate the absolute risk of a first KC in patients with AK.

METHODS

The risk-prediction model was based on the prospective population-based Rotterdam Study cohort. We hereto analysed the data of participants with at least one AK lesion at cohort baseline using a multivariable Cox proportional hazards model and included 13 a priori defined candidate predictor variables considering phenotypic, genetic and lifestyle risk factors. KCs were identified by linkage of the data with the Dutch Pathology Registry.

RESULTS

Of the 1169 AK participants at baseline, 176 (15·1%) developed a KC after a median follow-up of 1·8 years. The final model with significant predictors was obtained after backward stepwise selection and comprised the presence of four to nine AKs [hazard ratio (HR) 1·68, 95% confidence interval (CI) 1·17-2·42], 10 or more AKs (HR 2·44, 95% CI 1·65-3·61), AK localization on the upper extremities (HR 0·75, 95% CI 0·52-1·08) or elsewhere except the head (HR 1·40, 95% CI 0·98-2·01) and coffee consumption (HR 0·92, 95% CI 0·84-1·01). Evaluation of the discriminative ability of the model showed a bootstrap validated concordance index (c-index) of 0·60.

CONCLUSIONS

We showed that the risk of KC in patients with AK can be calculated with the use of four easily assessable predictor variables. Given the c-index, extension of the model with additional, currently unknown predictor variables is desirable. Linked Comment: Kim et al. Br J Dermatol 2020; 183:415-416.

Single-cell RNA-sequencing of differentiating iPS cells reveals dynamic genetic effects on gene expression

Recent developments in stem cell biology have enabled the study of cell fate decisions in early human development that are impossible to study in vivo. However, understanding how development varies across individuals and, in particular, the influence of common genetic variants during this process has not been characterised. Here, we exploit human iPS cell lines from 125 donors, a pooled experimental design, and single-cell RNA-sequencing to study population variation of endoderm differentiation. We identify molecular markers that are predictive of differentiation efficiency of individual lines, and utilise heterogeneity in the genetic background across individuals to map hundreds of expression quantitative trait loci that influence expression dynamically during differentiation and across cellular contexts.