Identification of Susceptibility Loci for Cutaneous Squamous Cell Carcinoma

Association between lifetime smoking and cutaneous squamous cell carcinoma: A 2-sample Mendelian randomization study

Background/Purpose

Cutaneous squamous cell carcinoma (cSCC) is one of the most common malignancies worldwide. While several environmental risk factors for cSCC are well established, there is conflicting evidence on cigarette smoking (and its potential causal effect) and cSCC risk. Furthermore, it is unclear if these potential associations represent causal, modifiable risk factors for cSCC development. This study aims to assess the nature of the associations between cigarette smoking traits (smoking initiation, amount smoked, and lifetime smoking exposure) and cSCC risk using two-sample Mendelian randomization analyses.

Methods

Genetic instruments, based on common genetic variants associated with cigarette smoking traits (P < 5 × 10-8), were derived from published genome-wide association studies (GWASs). For cSCC, we used GWAS summary statistics from the Kaiser Permanente GERA cohort (7701 cSCC cases and 60,167 controls; all non-Hispanic Whites).

Results

We found modest evidence that genetically determined lifetime smoking was associated with cSCC (inverse-variance weighted method: OR[95% CI] = 1.47[1.09-1.98]; P = .012), suggesting it may be a causal risk factor for cSCC. We did not detect any evidence of association between genetically determined smoking initiation or amount smoked and cSCC risk.

Conclusion

Study findings highlight the importance of smoking prevention and may support risk-stratified cSCC screening strategies based on carcinogen exposure and other genetic and clinical information.

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.

Basonuclin-2 regulates extracellular matrix production and degradation

Epithelial-mesenchymal transition is essential for tissue patterning and organization. It involves both regulation of cell motility and alterations in the composition and organization of the ECM-a complex environment of proteoglycans and fibrous proteins essential for tissue homeostasis, signaling in response to chemical and biomechanical stimuli, and is often dysregulated under conditions such as cancer, fibrosis, and chronic wounds. Here, we demonstrate that basonuclin-2 (BNC2), a mesenchymal-expressed gene, that is, strongly associated with cancer and developmental defects across genome-wide association studies, is a novel regulator of ECM composition and degradation. We find that at endogenous levels, BNC2 controls the expression of specific collagens, matrix metalloproteases, and other matrisomal components in breast cancer cells, and in fibroblasts that are primarily responsible for the production and processing of the ECM within the tumour microenvironment. In so doing, BNC2 modulates the motile and invasive properties of cancers, which likely explains the association of high BNC2 expression with increasing cancer grade and poor patient prognosis.

Identification of Genetic Risk Factors for Keratinocyte Cancer in Immunosuppressed Solid Organ Transplant Recipients: A Case-Control Study

Because of long-term immunosuppression, solid organ transplant recipients are at increased risk for keratinocyte cancer. We matched solid organ transplant patients (n = 150), cases with keratinocyte cancers and tumor-free controls, considering the most important risk factors for keratinocyte cancer in solid organ transplant recipients. Using whole exome data of germline DNA from this patient cohort, we identified several genetic loci associated with the occurrence of multiple keratinocyte cancers. We found one genome-wide significant association of a common single nucleotide polymorphism located in EXOC3 (rs72698504). In addition, we found several variants with a p-value of less than 10^-5 associated with the number of keratinocyte cancers. These variants were located in the genes CYB561, WASHC1, PITRM1-AS1, MUC8, ABI3BP, and THBS2-AS1. Using whole exome sequencing data, we performed groupwise tests for rare missense variants in our dataset and found robust associations (p < 10^-6, Burden Zeggini test) between MC1R, EPHA8, EPO, MYCT1, ADGRG3, and MGME1 and keratinocyte cancer. Thus, overall, we detected genes involved in pigmentation/UV protection, tumor suppression, immunomodulation, intracellular traffic, and response to UV as genetic risk factors for multiple keratinocyte cancers in solid organ transplant recipients. We also grouped selected genes to pathways and found a selection of genes involved in the "cellular response to UV" to be significantly associated with multiple keratinocyte cancers.

Molecular Alterations in Cutaneous Squamous Cell Carcinoma in Immunocompetent and Immunosuppressed Hosts-A Systematic Review

The characterization of cutaneous squamous cell carcinoma (cSCC) at the molecular level is lacking in the current literature due to the high mutational burden of this disease. Immunosuppressed patients afflicted with cSCC experience considerable morbidity and mortality. In this article, we review the molecular profile of cSCC among the immunosuppressed and immunocompetent populations at the genetic, epigenetic, transcriptomic, and proteometabolomic levels, as well as describing key differences in the tumor immune microenvironment between these two populations. We feature novel biomarkers from the recent literature which may serve as potential targets for therapy.

Microbiota, Oxidative Stress, and Skin Cancer: An Unexpected Triangle

Mounting evidence indicates that the microbiota, the unique combination of micro-organisms residing in a specific environment, plays an essential role in the development of a wide range of human diseases, including skin cancer. Moreover, a persistent imbalance of microbial community, named dysbiosis, can also be associated with oxidative stress, a well-known emerging force involved in the pathogenesis of several human diseases, including cutaneous malignancies. Although their interplay has been somewhat suggested, the connection between microbiota, oxidative stress, and skin cancer is a largely unexplored field. In the present review, we discuss the current knowledge on these topics, suggesting potential therapeutic strategies.

Epidemiology and survival outcomes of patients with orbital region non-cutaneous squamous cell carcinoma: a population-based analysis

Background

Non-cutaneous squamous cell carcinoma (ncSCC) of the orbital region is very rare. Thus, its epidemiological characteristics and prognosis are poorly understood. The aim of the study was to assess the epidemiological characteristics and survival outcomes of ncSCC of the orbital region.

Methods

Incidence and demographic data on ncSCC of the orbital region were extracted from the Surveillance, Epidemiology, and End Results (SEER) database and analyzed. The chi-square test was used to calculate the differences between groups. Univariate and multivariate Cox regression analyses were performed to determine the independent prognostic factors for disease-specific survival (DSS) and overall survival (OS).

Results

The overall incidence of ncSCC in the orbital region from 1975 to 2019 was 0.68/1,000,000, and the incidence showed an increasing trend during this period. A total of 1,265 patients with ncSCC of the orbital region (mean age, 65.3 years) were identified in the SEER database. Of these, 65.1% were aged ≥60 years, 87.4% were White, and 73.5% were male. The conjunctiva (74.5%) was the most common primary site, followed by the orbit (12.1%), lacrimal apparatus (10.8%), and overlapping lesion of the eye and adnexa (2.7%). Multivariate Cox regression analysis revealed that age, primary site, SEER summary stage, and surgery were independent prognostic factors for DSS, whereas age, sex, marital status, primary site, SEER summary stage, and surgery were independent prognostic factors for OS.

Conclusions

The incidence of ncSCC in the orbital region has increased over the past 40 years. It usually affects White people, men, and people aged ≥60 years, and its most common site is the conjunctiva. Orbital SCC has worse survival outcomes than SCC of other sites in the orbital region. Surgery is the independent protective treatment for ncSCC of the orbital region.

Interventions After First Post-Transplant Cutaneous Squamous Cell Carcinoma: A Proposed Decision Framework

Cutaneous squamous cell carcinoma (CSCC) is a major cause of morbidity and mortality after organ transplant. Many patients subsequently develop multiple CSCC following a first CSCC, and the risk of metastasis and death is significantly increased compared to the general population. Post-transplant CSCC represents a disease at the interface of dermatology and transplant medicine. Both systemic chemoprevention and modulation of immunosuppression are frequently employed in patients with multiple CSCC, yet there is little consensus on their use after first CSCC to reduce risk of subsequent tumors. While relatively few controlled trials have been undertaken, extrapolation of observational data suggests the most effective interventions may be at the time of first CSCC. We review the need for intervention after a first post-transplant CSCC and evidence for use of various approaches as secondary prevention, before discussing barriers preventing engagement with this approach and finally highlight areas for future research. Close collaboration between specialties to ensure prompt deployment of these interventions after a first CSCC may improve patient outcomes.

Inflammatory Bowel Disease Risk Variants Are Associated with an Increased Risk of Skin Cancer

BACKGROUND

Inflammatory bowel disease is associated with an increased risk of skin cancer. The aims of this study were to determine whether IBD susceptibility variants are also associated with skin cancer susceptibility and if such risk is augmented by use of immune-suppressive therapy.

METHODS

The discovery cohort included participants in the UK Biobank. The validation cohort included participants in the Michigan Genomics Initiative. The primary outcome of interest was skin cancer, subgrouped into nonmelanoma skin cancers (NMSC) and melanoma skin cancers (MSC). Multivariable logistic regression with matched controls (3 controls:1 case) was performed to identify genomic predictors of skin malignancy in the discovery cohort. Variants with P < .05 were tested for replication in the validation cohort. Validated Single nucleotide polymorphisms were then evaluated for effect modification by immune-suppressive medications.

RESULTS

The discovery cohort included 10,247 cases of NMSC and 1883 cases of MSC. The validation cohort included 7334 cases of NMSC and 3304 cases of MSC. Twenty-nine variants were associated with risk of NMSC in the discovery cohort, of which 5 replicated in the validation cohort (increased risk, rs7773324-A [DUSP22; IRF4], rs2476601-G [PTPN22], rs1847472-C [BACH2], rs72810983-A [CPEB4]; decreased risk, rs6088765-G [PROCR; MMP24]). Twelve variants were associated with risk of MSC in the discovery cohort, of which 4 were replicated in the validation cohort (increased risk, rs61839660-T [IL2RA]; decreased risk, rs17391694-C [GIPC2; MGC27382], rs6088765-G [PROCR; MMP24], and rs1728785-C [ZFP90]). No effect modification was observed.

CONCLUSIONS

The results of this study highlight shared genetic susceptibility across IBD and skin cancer, with increased risk of NMSC in those who carry risk variants in IRF4, PTPN22, CPEB4, and BACH2 and increased risk of MSC in those who carry a risk variant in IL2RA.

Evolutionary Origins of Metabolic Reprogramming in Cancer

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

Novel genetic associations with five aesthetic facial traits: A genome-wide association study in the Chinese population

Background: The aesthetic facial traits are closely related to life quality and strongly influenced by genetic factors, but the genetic predispositions in the Chinese population remain poorly understood.

Methods: A genome-wide association studies (GWAS) and subsequent validations were performed in 26,806 Chinese on five facial traits: widow’s peak, unibrow, double eyelid, earlobe attachment, and freckles. Functional annotation was performed based on the expression quantitative trait loci (eQTL) variants, genome-wide polygenic scores (GPSs) were developed to represent the combined polygenic effects, and single nucleotide polymorphism (SNP) heritability was presented to evaluate the contributions of the variants.

Results: In total, 21 genetic associations were identified, of which ten were novel: GMDS-AS1 (rs4959669, p = 1.29 × 10⁻⁴⁹) and SPRED2 (rs13423753, p = 2.99 × 10⁻¹⁴) for widow’s peak, a previously unreported trait; FARSB (rs36015125, p = 1.96 × 10⁻²¹) for unibrow; KIF26B (rs7549180, p = 2.41 × 10⁻¹⁵), CASC2 (rs79852633, p = 4.78 × 10⁻¹¹⁾, RPGRIP1L (rs6499632, p = 9.15 × 10⁻¹¹), and PAX1 (rs147581439, p = 3.07 × 10⁻⁸) for double eyelid; ZFHX3 (rs74030209, p = 9.77 × 10⁻¹⁴) and LINC01107 (rs10211400, p = 6.25 × 10⁻¹⁰) for earlobe attachment; and SPATA33 (rs35415928, p = 1.08 × 10⁻⁸) for freckles. Functionally, seven identified SNPs tag the missense variants and six may function as eQTLs. The combined polygenic effect of the associations was represented by GPSs and contributions of the variants were evaluated using SNP heritability.

Conclusion: These identifications may facilitate a better understanding of the genetic basis of features in the Chinese population and hopefully inspire further genetic research on facial development.

Cross-ancestry genome-wide meta-analysis of 61,047 cases and 947,237 controls identifies new susceptibility loci contributing to lung cancer

To identify new susceptibility loci to lung cancer among diverse populations, we performed cross-ancestry genome-wide association studies in European, East Asian and African populations and discovered five loci that have not been previously reported. We replicated 26 signals and identified 10 new lead associations from previously reported loci. Rare-variant associations tended to be specific to populations, but even common-variant associations influencing smoking behavior, such as those with CHRNA5 and CYP2A6, showed population specificity. Fine-mapping and expression quantitative trait locus colocalization nominated several candidate variants and susceptibility genes such as IRF4 and FUBP1. DNA damage assays of prioritized genes in lung fibroblasts indicated that a subset of these genes, including the pleiotropic gene IRF4, potentially exert effects by promoting endogenous DNA damage.

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.

Genome-Wide Analysis Identifies Candidate Genes Encoding Beak Color of Duck

Beak color diversity is a broadly occurring phenomenon in birds. Here, we used ducks to identify candidate genes for yellow, black, and spotted beaks. For this, an F₂ population consisting of 275 ducks was genotyped using whole genome resequencing containing 12.6 M single-nucleotide polymorphisms (SNPs) and three beak colors. Genome-wide association studies (GWAS) was used to identify the candidate and potential SNPs for three beak colors in ducks (yellow, spotted, and black). The results showed that 2753 significant SNPs were associated with black beaks, 7462 with yellow, and 17 potential SNPs with spotted beaks. Based on SNP annotation, MITF, EDNRB2, members of the POU family, and the SLC superfamily were the candidate genes regulating pigmentation. Meanwhile, isoforms MITF-M and EDNRB2 were significantly different between black and yellow beaks. MITF and EDNRB2 likely play a synergistic role in the regulation of melanin synthesis, and their mutations contribute to phenotypic differences in beak melanin deposition among individuals. This study provides new insights into genetic factors that may influence the diversity of beak color.

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.

Citrus-Gene interaction and melanoma risk in the UK Biobank

High citrus consumption may increase melanoma risk; however, little is known about the biological mechanisms of this association, or whether it is modified by genetic variants. We conducted a genome-wide analysis of gene-citrus consumption interactions on melanoma risk among 1563 melanoma cases and 193 296 controls from the UK Biobank. Both the 2-degrees-of-freedom (df) joint test of genetic main effect and gene-environment (G-E) interaction and the standard 1-df G-E interaction test were performed. Three index SNPs (lowest P-value SNP among highly correlated variants [r²  > .6]) were identified from among the 365 genome-wide significant 2-df test results (rs183783391 on chromosome 3 [MITF], rs869329 on chromosome 9 [MTAP] and rs11446223 on chromosome 16 [DEF8]). Although all three were statistically significant for the 2-df test (4.25e-08, 1.98e-10 and 4.93e-13, respectively), none showed evidence of interaction according to the 1-df test (P = .73, .24 and .12, respectively). Eight nonindex, 2-df test significant SNPs on chromosome 16 were significant (P < .05) according to the 1-df test, providing evidence of citrus-gene interaction. Seven of these SNPs were mapped to AFG3L1P (rs199600347, rs111822773, rs113178244, rs3803683, rs73283867, rs78800020, rs73283871), and one SNP was mapped to GAS8 (rs74583214). We identified several genetic loci that may elucidate the association between citrus consumption and melanoma risk. Further studies are needed to confirm these findings.

Cutaneous squamous cell carcinoma arising in immunosuppressed patients: a systematic review of tumor profiling studies

As solid organ transplantation becomes more prevalent, more individuals are living as members of the immunosuppressed population with an elevated risk for cutaneous squamous cell carcinoma (cSCC). Although great progress has been made in understanding the pathogenesis of cSCC in general, little is known about the drivers of tumorigenesis in immunosuppressed patients and organ-transplant recipients, specifically. This systematic review sought to synthesize information regarding the genetic and epigenetic alterations as well as changes in protein and mRNA expression that place this growing population at risk for cSCC, influence treatment response, and promote tumor aggressiveness. This review will provide investigators with a framework to identify future areas of investigation and clinicians with additional insight into how to best manage these patients.

Cutaneous Squamous Cell Carcinoma: The Frontier of Cancer Immunoprevention

Cutaneous squamous cell carcinoma (cSCC) is the second most common cancer, with its incidence rising steeply. Immunosuppression is a well-established risk factor for cSCC, and this risk factor highlights the critical role of the immune system in regulating cSCC development and progression. Further highlighting the nature of cSCC as an immunological disorder, substantial evidence demonstrates a tight association between cSCC risk and age-related immunosenescence. Besides the proven efficacy of immune checkpoint blockade therapy for advanced cSCC, novel immunotherapy that targets cSCC precursor lesions has shown efficacy for cSCC prevention. Furthermore, the appreciation of the interplay between keratinocytes, commensal papillomaviruses, and the immune system has revealed the possibility for the development of a preventive cSCC vaccine. cSCC shares fundamental aspects of its origin and pathogenesis with mucosal SCCs. Therefore, advances in the field of cSCC immunoprevention will inform our approach to the management of mucosal SCCs and potentially other epithelial cancers.

Seven key hub genes identified by gene co-expression network in cutaneous squamous cell carcinoma

Background

Cutaneous squamous cell carcinoma (cSCC) often follows actinic keratosis (AK) and is the second most common skin cancer worldwide. To reduce metastasis risk, it is important to diagnose and treat cSCC early. This study aimed to identify hub genes associated with cSCC and AK.

Methods

This study used three datasets GSE45216, GSE98774, and GSE108008. We combined samples from the GSE45216 and GSE98774 datasets into the new dataset GSE45216–98774. We applied a weighted gene co-expression network analysis (WGCNA) to investigate key modules and hub genes associated with cSCC and AK. We considered the hub genes found in both the GSE45216–98774 and GSE108008 datasets as validated hub genes. We tested whether the expression of hub genes could predict patient survival outcomes in other cancers using TCGA pan-cancer data.

Results

We identified modules most relevant to cSCC and AK. Additionally, we identified and validated seven hub genes of cSCC: GATM, ARHGEF26, PTHLH, MMP1, POU2F3, MMP10 and GATA3. We did not find validated hub genes for AK. Each hub gene was significantly associated with the survival of various cancer types. Only GATA3 was significantly associated with melanoma survival.

Conclusions

We applied WGCNA to find seven hub genes that play important roles in cSCC tumorigenesis. These results provide new insights that help explain the pathogenesis of cSCC. These hub genes may become biomarkers or therapeutic targets for accurate diagnosis and treatment of cSCC in the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08604-y.

A custom capture sequence approach for oculocutaneous albinism identifies structural variant alleles at the OCA2 locus

Oculocutaneous albinism (OCA) is a heritable disorder of pigment production that manifests as hypopigmentation and altered eye development. Exon sequencing of known OCA genes is unsuccessful in producing a complete molecular diagnosis for a significant number of affected individuals. We sequenced the DNA of individuals with OCA using short-read custom capture sequencing that targeted coding, intronic, and noncoding regulatory regions of known OCA genes, and genome-wide association study-associated pigmentation loci. We identified an OCA2 complex structural variant (CxSV), defined by a 143 kb inverted segment reintroduced in intron 1, upstream of the native location. The corresponding CxSV junctions were observed in 11/390 probands screened. The 143 kb CxSV presents in one family as a copy number variant duplication for the 143 kb region. In the remaining 10/11 families, the 143 kb CxSV acquired an additional 184 kb deletion across the same region, restoring exons 3-19 of OCA2 to a copy-number neutral state. Allele-associated haplotype analysis found rare SNVs rs374519281 and rs139696407 are linked with the 143 kb CxSV in both OCA2 alleles. For individuals in which customary molecular evaluation does not reveal a biallelic OCA diagnosis, we recommend preliminary screening for these haplotype-associated rare variants, followed by junction-specific validation for the OCA2 143 kb CxSV.

Identifying Susceptibility Loci for Cutaneous Squamous Cell Carcinoma Using a Fast Sequence Kernel Association Test

Cutaneous squamous cell carcinoma (cSCC) accounts for about 20% of all skin cancers, the most common type of malignancy in the United States. Genome-wide association studies (GWAS) have successfully identified multiple genetic variants associated with the risk of cSCC. Most of these studies were single-locus-based, testing genetic variants one-at-a-time. In this article, we performed gene-based association tests to evaluate the joint effect of multiple variants, especially rare variants, on the risk of cSCC by using a fast sequence kernel association test (fastSKAT). The study included 1,710 cSCC cases and 24,304 cancer-free controls from the Nurses' Health Study, the Nurses' Health Study II and the Health Professionals Follow-up Study. We used UCSC Genome Browser to define gene units as candidate loci, and further evaluated the association between all variants within each gene unit and disease outcome. Four genes HP1BP3, DAG1, SEPT7P2, and SLFN12 were identified using Bonferroni adjusted significance level. Our study is complementary to the existing GWASs, and our findings may provide additional insights into the etiology of cSCC. Further studies are needed to validate these findings.

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.

Biophysical Compatibility of a Heterotrimeric Tyrosinase-TYRP1-TYRP2 Metalloenzyme Complex

Tyrosinase (TYR) is a copper-containing monooxygenase central to the function of melanocytes. Alterations in its expression or activity contribute to variations in skin, hair and eye color, and underlie a variety of pathogenic pigmentary phenotypes, including several forms of oculocutaneous albinism (OCA). Many of these phenotypes are linked to individual missense mutations causing single nucleotide variants and polymorphisms (SNVs) in TYR. We previously showed that two TYR homologues, TYRP1 and TYRP2, modulate TYR activity and stabilize the TYR protein. Accordingly, to investigate whether TYR, TYRP1, and TYRP2 are biophysically compatible with various heterocomplexes, we computationally docked a high-quality 3D model of TYR to the crystal structure of TYRP1 and to a high-quality 3D model of TYRP2. Remarkably, the resulting TYR-TYRP1 heterodimer was complementary in structure and energy with the TYR-TYRP2 heterodimer, with TYRP1 and TYRP2 docking to different adjacent surfaces on TYR that apposed a third realistic protein interface between TYRP1-TYRP2. Hence, the 3D models are compatible with a heterotrimeric TYR-TYRP1-TYRP2 complex. In addition, this heterotrimeric TYR-TYRP1-TYRP2 positioned the C-terminus of each folded enzymatic domain in an ideal position to allow their C-terminal transmembrane helices to form a putative membrane embedded three-helix bundle. Finally, pathogenic TYR mutations causing OCA1A, which also destabilize TYR biochemically, cluster on an unoccupied protein interface at the periphery of the heterotrimeric complex, suggesting that this may be a docking site for OCA2, an anion channel. Pathogenic OCA2 mutations result in similar phenotypes to those produced by OCA1A TYR mutations. While this complex may be difficult to detect in vitro, due to the complex environment of the vertebrate cellular membranous system, our results support the existence of a heterotrimeric complex in melanogenesis.

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.

Loss of retinoic acid receptor-related receptor alpha (Rorα) promotes the progression of UV-induced cSCC

Cutaneous squamous cell carcinoma (cSCC) is prevalent in the world, accounting for a huge part of non-melanoma skin cancer. Most cSCCs are associated with a distinct pre-cancerous lesion, the actinic keratosis (AK). However, the progression trajectory from normal skin to AK and cSCC has not been fully demonstrated yet. To identify genes involved in this progression trajectory and possible therapeutic targets for cSCC, here we constructed a UV-induced cSCC mouse model covering the progression from normal skin to AK to cSCC, which mimicked the solar UV radiation perfectly using the solar-like ratio of UVA and UVB, firstly. Then, transcriptome analysis and a series of bioinformatics analyses and cell experiments proved that Rorα is a key transcript factor during cSCC progression. Rorα could downregulate the expressions of S100a9 and Sprr2f in cSCC cells, which can inhibit the proliferation and migration in cSCC cells, but not the normal keratinocyte. Finally, further animal experiments confirmed the inhibitory effect of cSCC growth by Rorα in vivo. Our findings showed that Rorα would serve as a potential novel target for cSCC, which will facilitate the treatment of cSCC in the future.

Genetic ancestry, skin pigmentation, and the risk of cutaneous squamous cell carcinoma in Hispanic/Latino and non-Hispanic white populations

Although cutaneous squamous cell carcinoma (cSCC) is one of the most common malignancies in individuals of European ancestry, the incidence of cSCC in Hispanic/Latinos is also increasing. cSCC has both a genetic and environmental etiology. Here, we examine the role of genetic ancestry, skin pigmentation, and sun exposure in Hispanic/Latinos and non-Hispanic whites on cSCC risk. We observe an increased cSCC risk with greater European ancestry (P = 1.27 × 10⁻⁴²) within Hispanic/Latinos and with greater northern (P = 2.38 × 10⁻⁶⁵) and western (P = 2.28 × 10⁻⁴⁹) European ancestry within non-Hispanic whites. These associations are significantly, but not completely, attenuated after considering skin pigmentation-associated loci, history of actinic keratosis, and sun-protected versus sun-exposed anatomical sites. We also report an association of the well-known pigment variant Ala111Thr (rs1426654) at SLC24A5 with cSCC in Hispanic/Latinos. These findings demonstrate a strong correlation of northwestern European genetic ancestry with cSCC risk in both Hispanic/Latinos and non-Hispanic whites, largely but not entirely mediated through its impact on skin pigmentation.

Eric Jorgenson and Hélène Choquet et al. find that northwestern European genetic ancestry is associated with increased risk of cutaneous squamous cell carcinoma (cSCC) in non-Hispanic whites, and more so in Hispanic/Latinos of the US. The ancestry effect is largely, but not entirely explained by genetic determinants of skin pigmentation in both populations.

Functional annotation of melanoma risk loci identifies novel susceptibility genes

Genome-wide association study (GWAS)-identified single-nucleotide polymorphisms (SNPs) are tag SNPs located in both transcribed and non-coding regulatory DNA regions, rather than representing causal or functional variants for disease. To identify functional variants or genes for melanoma susceptibility, we used functional mapping and annotation (FUMA) to perform functional annotation of the summary statistics of 2541 significant melanoma risk SNPs (P < 5 × 10-8) identified by GWAS. The original GWAS melanoma study included 15 990 cases and 26 409 controls, representing the largest international meta-analysis of melanoma susceptibility. We prioritized 330 unique genes, including those in immune cytokine signaling pathways, from 19 loci through positional, expression quantitative trait locus, and chromatin interaction mapping. In comparison, only 38 melanoma-related genes were identified in the original meta-analysis. In addition to the well-known melanoma susceptibility genes confirmed in the meta-analysis (MC1R, CDKN2A, TERT, OCA2 and ARNT/SETDB1), we also identified additional novel genes using FUMA to map SNPs to genes. Through chromatin interaction mapping, we prioritized IFNA7, IFNA10, IFNA16, IFNA17, IFNA14, IFNA6, IFNA21, IFNA4, IFNE and IFNA5; these 10 most significant genes are all involved in immune system and cytokine signaling pathways. In the gene analysis, we identified 72 genes with a P < 2.5 × 10-6. The genes associated with melanoma risk were DEF8 (P = 1.09 × 10-57), DBNDD1 (P = 2.19 × 10-42), SPATA33 (P = 3.54 × 10-38) and MC1R (P = 1.04 × 10-36). In summary, this study identifies novel putative melanoma susceptibility genes and provides a guide for further experimental validation of functional variants and disease-related genes.

Predicting target genes of non-coding regulatory variants with IRT

SUMMARY

Interpreting genetic variants of unknown significance (VUS) is essential in clinical applications of genome sequencing for diagnosis and personalized care. Non-coding variants remain particularly difficult to interpret, despite making up a large majority of trait associations identified in genome-wide association studies (GWAS) analyses. Predicting the regulatory effects of non-coding variants on candidate genes is a key step in evaluating their clinical significance. Here, we develop a machine-learning algorithm, Inference of Connected expression quantitative trait loci (eQTLs) (IRT), to predict the regulatory targets of non-coding variants identified in studies of eQTLs. We assemble datasets using eQTL results from the Genotype-Tissue Expression (GTEx) project and learn to separate positive and negative pairs based on annotations characterizing the variant, gene and the intermediate sequence. IRT achieves an area under the receiver operating characteristic curve (ROC-AUC) of 0.799 using random cross-validation, and 0.700 for a more stringent position-based cross-validation. Further evaluation on rare variants and experimentally validated regulatory variants shows a significant enrichment in IRT identifying the true target genes versus negative controls. In gene-ranking experiments, IRT achieves a top-1 accuracy of 50% and top-3 accuracy of 90%. Salient features, including GC-content, histone modifications and Hi-C interactions are further analyzed and visualized to illustrate their influences on predictions. IRT can be applied to any VUS of interest and each candidate nearby gene to output a score reflecting the likelihood of regulatory effect on the expression level. These scores can be used to prioritize variants and genes to assist in patient diagnosis and GWAS follow-up studies.

AVAILABILITY AND IMPLEMENTATION

Codes and data used in this work are available at https://github.com/miaecle/eQTL_Trees.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Landscape of copy number aberrations in esophageal squamous cell carcinoma from a high endemic region of South Africa

Background

Esophageal squamous cell carcinoma (ESCC) is an aggressive cancer with one of the highest world incidences in the Eastern Cape region of South Africa. Several genome wide studies have been performed on ESCC cohorts from Asian countries, North America, Malawi and other parts of the world but none have been conducted on ESCC tumors from South Africa to date, where the molecular pathology and etiology of this disease remains unclear. We report here tumor associated copy number changes observed in 51 ESCC patients’ samples from the Eastern Cape province of South Africa.

Methods

We extracted tumor DNA from 51 archived ESCC specimens and interrogated tumor associated DNA copy number changes using Affymetrix® 500 K SNP array technology. The Genomic Identification of Significant Targets in Cancer (GISTIC 2.0) algorithm was applied to identify significant focal regions of gains and losses. Gains of the top recurrent cancer genes were validated by fluorescence in situ hybridization and their protein expression assessed by immunohistochemistry.

Results

Twenty-three significant focal gains were identified across samples. Gains involving the CCND1, MYC, EGFR and JAG1 loci recapitulated those described in studies on Asian and Malawian cohorts. The two most significant gains involved the chromosomal sub-bands 3q28, encompassing the TPRG1 gene and 11q13.3 including the CTTN, PPFIA1and SHANK2 genes. There was no significant homozygous loss and the most recurrent hemizygous deletion involved the B3GAT1 gene on chromosome 11q25. Focal gains on 11q13.3 in 37% of cases (19/51), consistently involved CTTN and SHANK2 genes. Twelve of these cases (23,5%), had a broader region of gain that also included the CCND1, FGF19, FGF4 and FGF3 genes. SHANK2 and CTTN are co-amplified in several cancers, these proteins interact functionally together and are involved in cell motility. Immunohistochemistry confirmed both Shank2 (79%) and cortactin (69%) protein overexpression in samples with gains of these genes. In contrast, cyclin D1 (65%) was moderately expressed in samples with CCND1 DNA gain.

Conclusions

This study reports copy number changes in a South African ESCC cohort and highlights similarities and differences with cohorts from Asia and Malawi. Our results strongly suggest a role for CTTN and SHANK2 in the pathogenesis of ESCC in South Africa.

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 meta-analysis identifies eight new susceptibility loci for cutaneous squamous cell carcinoma

Cutaneous squamous cell carcinoma (SCC) is one of the most common cancers in the United States. Previous genome-wide association studies (GWAS) have identified 14 single nucleotide polymorphisms (SNPs) associated with cutaneous SCC. Here, we report the largest cutaneous SCC meta-analysis to date, representing six international cohorts and totaling 19,149 SCC cases and 680,049 controls. We discover eight novel loci associated with SCC, confirm all previously associated loci, and perform fine mapping of causal variants. The novel SNPs occur within skin-specific regulatory elements and implicate loci involved in cancer development, immune regulation, and keratinocyte differentiation in SCC susceptibility.

The authors perform a meta-analysis of cutaneous squamous cell carcinoma, identifying causal variants within skin-specific regulatory elements.

Transcriptomic Profiles of MV4-11 and Kasumi 1 Acute Myeloid Leukemia Cell Lines Modulated by Epigenetic Modifiers Trichostatin A and 5-Azacytidine

Background: Acute myeloid leukemia (AML) is the most common form of acute leukemias in adults which is clinically and molecularly heterogeneous. Several risk and genetic factors have been widely investigated to characterize AML. However, the concomitant epigenetic factors in controlling the gene expression lead to AML transformation was not fully understood. This study was aimed to identify epigenetically regulated genes in AML cell lines induced by epigenetic modulating agents, Trichostatin A (TSA) and 5-Azacytidine (5-Aza). Materials and Methods: MV4-11 and Kasumi 1 were treated with TSA and/or 5-Aza at IC₅₀ concentration. Gene expression profiling by microarray was utilized using SurePrint G3 Human Gene Expression v3. Gene ontology and KEGG pathway annotations were analyzed by DAVID bioinformatics software using EASE enrichment score. mRNA expression of the differentially expressed genes were verified by quantitative real time PCR. Results: Gene expression analysis revealed a significant changes in the expression of 24,822, 15,720, 15,654 genes in MV4-11 and 12,598, 8828, 18,026 genes in Kasumi 1, in response to TSA, 5-Aza and combination treatments, respectively, compared to non-treated (p<0.05). 7 genes (SOCS3, TUBA1C, CCNA1, MAP3K6, PTPRC, STAT6 and RUNX1) and 4 genes (ANGPTL4, TUBB2A, ADAM12 and PTPN6) shown to be predominantly expressed in MV4-11 and Kasumi 1, respectively (EASE<0.1). The analysis also revealed phagosome pathway commonly activated in both cell lines. Conclusion: Our data showed a distinct optimal biological characteristic and pathway in different types of leukemic cell lines. These finding may help in the identification of cell-specific epigenetic biomarker in the pathogenesis of AML.

Association between Human Leukocyte Antigen Type and Keratinocyte Carcinoma Risk in Renal Transplant Recipients

Keratinocyte carcinoma (KC), defined as squamous cell carcinoma and basal cell carcinoma, is the most common malignancy among white, non-Hispanic renal transplant recipients. Although recent genome-wide association studies reported that class II HLA is associated with KC risk, epidemiologic data on HLA type and KC risk in renal transplant recipients is limited. Using an institutional cohort of white, non-Hispanic renal transplant recipients transplanted between 1993 and 2017, we examined the association between pretransplant molecular HLA types and KC risk. Posttransplant KCs were captured using the International Classification of Diseases codes and validated using pathology reports. Cox proportional hazards regression models were used to estimate hazard ratios of incident KC, squamous cell carcinoma, and basal cell carcinoma, adjusting for age, male sex, history of KC, Charlson comorbidity index, HLA mismatch, transplant type, year of transplant, and the type of immunosuppression. Among 617 subjects (mean age 53 years, 67% male), 10% developed posttransplant KC. Multivariable Cox regression analyses showed HLA-DRB1∗13 was associated with KC risk (hazard ratio, 1.84; 95% confidence interval, 1.00-3.38) and squamous cell carcinoma risk (hazard ratio, 2.24; 95% confidence interval, 1.12-4.49), whereas HLA-DRB1∗14 (hazard ratio, 2.81; 95% confidence interval, 1.14-6.91) was associated with basal cell carcinoma risk. Our findings suggest that a subset of renal transplant recipients with specific HLA polymorphisms may be at increased KC risk.

Sex-Stratified Polygenic Risk Score Identifies Individuals at Increased Risk of Basal Cell Carcinoma

The incidence of basal cell carcinoma (BCC) is higher among men than women. Susceptibility loci for BCC have been identified through genome-wide association studies, and two previous studies have found polygenic risk scores (PRS) to be significantly associated with the risk of BCC. However, to our knowledge, sex-stratified PRS analyses examining the genetic contribution to BCC risk among men and women have not been previously reported. To quantify the contribution of genetic variability on the BCC risk by sex, we derived a polygenic risk score and estimated the genetic relative risk distribution for men and women. Using 29 published single nucleotide polymorphisms, we found that the estimated relative risk of BCC increases with higher percentiles of the polygenic risk score. For men, the estimated risk of BCC is twice the average population risk at the 88th percentile, while for women, this occurs at the 99th percentile. Our findings indicate that there is a significant impact of genetic variation on the risk of developing BCC and that this impact may be greater for men than for women. Polygenic risk scores may be clinically useful tools for risk stratification, particularly in combination with other known risk factors for BCC development.

Cutaneous keratinocyte cancers of the head and neck: Epidemiology, risk factors and clinical, dermoscopic and reflectance confocal microscopic features

Keratinocyte cancers are the most common malignancy among people with European ancestry, and are very common on sun-exposed areas of the head and neck. Incidence is directly correlated with latitude and annual ultraviolet radiation incidence, although there are a number of other environmental, occupational and genetic risk factors, and keratinocyte cancers become more common at middle age. Basal cell carcinomas (BCC) are the most common, comprising 80% of keratinocyte cancers, but have a very low rate of metastases and low mortality. Squamous cell carcinomas (SCC) make up 20% of keratinocyte cancers, and have relatively infrequent metastases, at 5-16%. While there are no precursor lesions for BCC, SCC represents the final stage in a spectrum of cellular atypia and dysplasia, from actinic keratoses to in situ SCC to invasive SCC. Dermoscopy is a well-established diagnostic tool for keratinocyte cancers, and reflectance confocal microscopy is emerging as another useful diagnostic tool, particularly on functionally and cosmetically sensitive sites like the face.

Genome-wide association studies and polygenic risk scores for skin cancer: clinically useful yet?

BACKGROUND

Genome-wide association studies (GWAS) have identified thousands of susceptibility variants, although most have been associated with small individual risk estimates that offer little predictive value. However, combining multiple variants into polygenic risk scores (PRS) may be more informative. Multiple studies have developed PRS composed of GWAS-identified variants for cutaneous cancers. This review highlights data from these studies.

OBJECTIVES

To review published GWAS and PRS studies for melanoma, cutaneous squamous cell carcinoma (cSCC) and basal cell carcinoma (BCC), and discuss their potential clinical utility.

METHODS

We searched PubMed and the National Human Genome Research Institute-European Bioinformatics Institute GWAS catalogue to identify relevant studies.

RESULTS

Results from 21 GWAS (11 melanoma, 3 cSCC, 7 BCC) and 11 PRS studies are summarized. Six loci in pigmentation genes overlap between these three cancers (ASIP/RALY, IRF4, MC1R, OCA2, SLC45A2 and TYR). Additional loci overlap for cSCC/BCC and BCC/melanoma, but no other loci are shared between cSCC and melanoma. PRS for melanoma show roughly two-to-threefold increases in risk and modest improvements in risk prediction (2-7% increases). PRS are associated with twofold and threefold increases in risk of cSCC and BCC, respectively, with small improvements (2% increase) in predictive ability.

CONCLUSIONS

Existing data indicate that PRS may offer small, but potentially meaningful, improvements to risk prediction. Additional research is needed to clarify the potential utility of PRS in cutaneous carcinomas. Clinical translation will require well-powered validation studies incorporating known risk factors to evaluate PRS as tools for screening. What's already known about this topic? Over 50 susceptibility loci for melanoma, basal cell carcinoma (BCC) and cutaneous squamous cell carcinoma (cSCC) have been identified in genome-wide association studies (GWAS). Polygenic risk scores (PRS) using variants identified from GWAS have also been developed for melanoma, BCC and cSCC, and investigated with respect to clinical risk prediction. What does this study add? This review provides an overview of GWAS findings and the potential clinical utility of PRS for melanoma, BCC and cSCC.

Exploring various polygenic risk scores for skin cancer in the phenomes of the Michigan genomics initiative and the UK Biobank with a visual catalog: PRSWeb

Polygenic risk scores (PRS) are designed to serve as single summary measures that are easy to construct, condensing information from a large number of genetic variants associated with a disease. They have been used for stratification and prediction of disease risk. The primary focus of this paper is to demonstrate how we can combine PRS and electronic health records data to better understand the shared and unique genetic architecture and etiology of disease subtypes that may be both related and heterogeneous. PRS construction strategies often depend on the purpose of the study, the available data/summary estimates, and the underlying genetic architecture of a disease. We consider several choices for constructing a PRS using data obtained from various publicly-available sources including the UK Biobank and evaluate their abilities to predict not just the primary phenotype but also secondary phenotypes derived from electronic health records (EHR). This study was conducted using data from 30,702 unrelated, genotyped patients of recent European descent from the Michigan Genomics Initiative (MGI), a longitudinal biorepository effort within Michigan Medicine. We examine the three most common skin cancer subtypes in the USA: basal cell carcinoma, cutaneous squamous cell carcinoma, and melanoma. Using these PRS for various skin cancer subtypes, we conduct a phenome-wide association study (PheWAS) within the MGI data to evaluate PRS associations with secondary traits. PheWAS results are then replicated using population-based UK Biobank data and compared across various PRS construction methods. We develop an accompanying visual catalog called PRSweb that provides detailed PheWAS results and allows users to directly compare different PRS construction methods.

Potential role of TP63 in apical periodontitis development

AIM

To investigate the expression of TP63 in apical periodontitis (AP) tissues and the association of single nucleotide polymorphisms (SNPs) in the TP63 gene with AP using a case-control dataset.

METHODOLOGY

Expression of TP63 in human AP lesions (apical abscess, radicular cyst, periapical granuloma) was evaluated using immunohistochemistry. A case-control association study was performed to assess the association of TP63 polymorphisms in individuals having AP with or without associated pain. Cases were defined as subjects with deep caries and AP (n = 151) and subjects with symptomatic apical periodontitis or acute apical abscess (n = 124). Subjects without AP (n = 169) and asymptomatic (n = 196) were used as controls, respectively. Saliva samples were collected as source of genomic DNA. Twelve SNPs in the TP63 gene were selected for genotyping using Taqman chemistry in real-time PCR. Data analysis was performed using PLINK software. The Bonferroni method was applied to correct for multiple testing; α ≤ 0.004 indicates significant differences between groups.

RESULTS

TP63 expression was evident in apical abscesses and radicular cysts, while weaker expression was observed in periapical granulomas. Positive expression was observed in mononuclear cells in the granulation tissues of all AP lesions. Regarding the presence of AP, a trend for allelic association was observed for rs16864812 and rs9810322 (P = 0.04) and rs9810322 genotypes were also nominally associated with AP under a dominant model (P = 0.04). When considering the presence of periapical pain, a trend for allelic and genotypic association was observed for rs10155037 (P = 0.03). Haplotypes were also associated with AP and periapical pain (P ≤ 0.05).

CONCLUSIONS

Apical periodontitis is a complex multifactorial condition and it is likely that multiple genes and environmental effects may influence its susceptibility, progression or both. TP63 variants may play a role in AP pathogenesis and susceptibility, individually or interactively with other genes. Additional studies in other populations and functional studies are needed to improve understanding of the role of TP63 in AP.

Keratinocyte Carcinomas: Current Concepts and Future Research Priorities

Cutaneous squamous cell carcinoma (cSCC) and basal cell carcinoma (BCC) are keratinocyte carcinomas, the most frequently diagnosed cancers in fair-skinned populations. Ultraviolet radiation (UVR) is the main driving carcinogen for these tumors, but immunosuppression, pigmentary factors, and aging are also risk factors. Scientific discoveries have improved the understanding of the role of human papillomaviruses (HPV) in cSCC as well as the skin microbiome and a compromised immune system in the development of both cSCC and BCC. Genomic analyses have uncovered genetic risk variants, high-risk susceptibility genes, and somatic events that underlie common pathways important in keratinocyte carcinoma tumorigenesis and tumor characteristics that have enabled development of prediction models for early identification of high-risk individuals. Advances in chemoprevention in high-risk individuals and progress in targeted and immune-based treatment approaches have the potential to decrease the morbidity and mortality associated with these tumors. As the incidence and prevalence of keratinocyte carcinoma continue to increase, strategies for prevention, including effective sun-protective behavior, educational interventions, and reduction of tanning bed access and usage, are essential. Gaps in our knowledge requiring additional research to reduce the high morbidity and costs associated with keratinocyte carcinoma include better understanding of factors leading to more aggressive tumors, the roles of microbiome and HPV infection, prediction of response to therapies including immune checkpoint blockade, and how to tailor both prevention and treatment to individual risk factors and needs.

Complement C3a promotes proliferation, migration and stemness in cutaneous squamous cell carcinoma

Background

Complement C3 has been shown to be highly expressed in cutaneous squamous cell carcinoma (cSCC) tumour tissues and is correlated with tumour cell growth. This study aimed to investigate the mechanism of C3 in cSCC malignant transformation.

Methods

C3 expression was analysed in cSCC cell lines A431, Tca8113, SCC13, HSC‐5 and HSC‐1 and in immortalized HaCaT keratinocytes. Proliferation and migration of cSCC were determined after C3a exposure. Expression of cyclin D1, cyclin E, vascular endothelial growth factor (VEGF), pro‐matrix metalloproteinase 1 (pro‐MMP1), pro‐matrix metalloproteinase 2 (pro‐MMP2), stemness factors, GSK‐3β, and β‐catenin were analyzed. Tumour growth was examined in a murine xenograft model.

Results

C3 expression was much more highly expressed in all cSCC cell lines than in HaCaT cells. C3a treatment significantly promoted cSCC cell proliferation and migration and upregulated cyclin D1, cyclin E, VEGF, pro‐MMP1 and pro‐MMP2 expression, which were impeded by the C3aR antagonist. Moreover, the expression of stemness factors Sox‐2, Nanog, Oct‐4, c‐Myc and CD‐44 was stimulated by C3a and slowed by C3aR disruption. Knockdown of Sox‐2 by siRNA transfection suppressed cell proliferation and migration, constrained VEGF secretion and inhibited pro‐MMP1 and pro‐MMP2 expression. C3a also activated the Wnt and β‐catenin pathway in cSCC cells. Disruption of C3aR expression dampened tumour growth and the expression of Wnt‐1, β‐catenin and Sox‐2 in the xenograft model.

Conclusions

C3a enhanced cell proliferation, migration and stemness in cSCC, and this activity was correlated with activation of the Wnt and β‐catenin pathway.

Human health in relation to exposure to solar ultraviolet radiation under changing stratospheric ozone and climate

The Montreal Protocol has limited increases in the UV-B (280-315 nm) radiation reaching the Earth's surface as a result of depletion of stratospheric ozone. Nevertheless, the incidence of skin cancers continues to increase in most light-skinned populations, probably due mainly to risky sun exposure behaviour. In locations with strong sun protection programs of long duration, incidence is now reducing in younger age groups. Changes in the epidemiology of UV-induced eye diseases are less clear, due to a lack of data. Exposure to UV radiation plays a role in the development of cataracts, pterygium and possibly age-related macular degeneration; these are major causes of visual impairment world-wide. Photodermatoses and phototoxic reactions to drugs are not uncommon; management of the latter includes recognition of the risks by the prescribing physician. Exposure to UV radiation has benefits for health through the production of vitamin D in the skin and modulation of immune function. The latter has benefits for skin diseases such as psoriasis and possibly for systemic autoimmune diseases such as multiple sclerosis. The health risks of sun exposure can be mitigated through appropriate sun protection, such as clothing with both good UV-blocking characteristics and adequate skin coverage, sunglasses, shade, and sunscreen. New sunscreen preparations provide protection against a broader spectrum of solar radiation, but it is not clear that this has benefits for health. Gaps in knowledge make it difficult to derive evidence-based sun protection advice that balances the risks and benefits of sun exposure.