Genome-wide association study identifies novel susceptibility loci for cutaneous squamous cell carcinoma

Genome wide association study and meta-analysis identified multiple new risk loci for freckles in 4813 Chinese individuals

Freckle is a prevalent pigmentary dermatosis with an obvious hereditary component. Dozens of freckles risk loci have been discovered through research on multiple traits or other diseases, rather than as an independent trait. To discover novel variants associated with freckles, we performed GWAS and meta-analysis in 4813 Chinese individuals. We conducted GWAS and meta-analysis of two cohorts: 197 patients and 1603 controls (Cohort I), and 336 patients and 2677 controls (Cohort II), both from China. Then we performed linkage disequilibrium (LD) analysis, eQTL study, and enrichment analysis with association results for functional implications. Finally, we discovered 59 new SNPs and 13 novel susceptibility genes associated with freckles (Pmeta <5 × 10-8), which has enriched the genetic research on freckles.

Mapping and annotating genomic loci to prioritize genes and implicate distinct polygenic adaptations for skin color

Evidence for adaptation of human skin color to regional ultraviolet radiation suggests shared and distinct genetic variants across populations. However, skin color evolution and genetics in East Asians are understudied. We quantified skin color in 48,433 East Asians using image analysis and identified associated genetic variants and potential causal genes for skin color as well as their polygenic interplay with sun exposure. This genome-wide association study (GWAS) identified 12 known and 11 previously unreported loci and SNP-based heritability was 23-24%. Potential causal genes were determined through the identification of nonsynonymous variants, colocalization with gene expression in skin tissues, and expression levels in melanocytes. Genomic loci associated with pigmentation in East Asians substantially diverged from European populations, and we detected signatures of polygenic adaptation. This large GWAS for objectively quantified skin color in an East Asian population improves understanding of the genetic architecture and polygenic adaptation of skin color and prioritizes potential causal genes.

Metastatic patterns and treatment options for head and neck cutaneous squamous cell carcinoma (Review)

According to current predictions, one-fifth of all Americans will develop skin cancer during their lifetime. Cutaneous squamous cell carcinoma (cSCC) most commonly occurs in the head and neck region, which is the area of the body with the highest level of sun exposure. High-risk head and neck cSCC (HNcSCC) is a broad category with numerous high-risk factors that are associated with unfavorable results. In cSCC staging systems, clinical and tumor traits that are likely to result in poor outcomes are identified. Metastasis occurs in ~2.5% of patients with cSCC, most often in the local lymph nodes, and there is some indication that lymph node metastasis has a distinct pattern based on the tumor site. Current findings on tumor molecular targets have suggested the use of systemic treatments, particularly immunotherapy (such as cemiplimab, pembrolizumab and nivolumab), over radiotherapy or chemotherapy for this type of metastasis. However, when used simultaneously with immunotherapy, radiotherapy may be beneficial in the treatment of metastatic HNcSCC by improving the efficacy of immunotherapy. The present review aims to assess the existing literature on metastatic HNcSCC pathways and treatment options, in order to define current and future directions. Notably, there is an urgent need to identify patients who may benefit from local or systemic cancer treatments. The treatment of lymph node metastasis presents a therapeutic challenge and requires comprehensive management.

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.

Identifying somatic fingerprints of cancers defined by germline and environmental risk factors

Numerous studies over the past generation have identified germline variants that increase specific cancer risks. Simultaneously, a revolution in sequencing technology has permitted high-throughput annotations of somatic genomes characterizing individual tumors. However, examining the relationship between germline variants and somatic alteration patterns is hugely challenged by the large numbers of variants in a typical tumor, the rarity of most individual variants, and the heterogeneity of tumor somatic fingerprints. In this article, we propose statistical methodology that frames the investigation of germline-somatic relationships in an interpretable manner. The method uses meta-features embodying biological contexts of individual somatic alterations to implicitly group rare mutations. Our team has used this technique previously through a multilevel regression model to diagnose with high accuracy tumor site of origin. Herein, we further leverage topic models from computational linguistics to achieve interpretable lower-dimensional embeddings of the meta-features. We demonstrate how the method can identify distinctive somatic profiles linked to specific germline variants or environmental risk factors. We illustrate the method using The Cancer Genome Atlas whole-exome sequencing data to characterize somatic tumor fingerprints in breast cancer patients with germline BRCA1/2 mutations and in head and neck cancer patients exposed to human papillomavirus.

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.

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.

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.

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.

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).

Deciphering the roles of aryl hydrocarbon receptor (AHR) in regulating carcinogenesis

Aryl hydrocarbon receptor (AHR) is a ligand-dependent receptor that belongs to the superfamily of basic helix-loop-helix (bHLH) transcription factors. The activation of the canonical AHR signaling pathway is known to induce the expression of cytochrome P450 enzymes, facilitating the detoxification metabolism in the human body. Additionally, AHR could interact with various signaling pathways such as epidermal growth factor receptor (EGFR), signal transducer and activator of transcription 3 (STAT3), hypoxia-inducible factor-1α (HIF-1α), nuclear factor ekappa B (NF-κβ), estrogen receptor (ER), and androgen receptor (AR) signaling pathways. Over the past 30 years, several studies have reported that various chemical, physical, or biological agents, such as tobacco, hydrocarbon compounds, industrial and agricultural chemical wastes, drugs, UV, viruses, and other toxins, could affect AHR expression or activity, promoting cancer development. Thus, it is valuable to overview how these factors regulate AHR-mediated carcinogenesis. Current findings have reported that many compounds could act as AHR ligands to drive the expressions of AHR-target genes, such as CYP1A1, CYP1B1, MMPs, and AXL, and other targets that exert a pro-proliferation or anti-apoptotic effect, like XIAP. Furthermore, some other physical and chemical agents, such as UV and 3-methylcholanthrene, could promote AHR signaling activities, increasing the signaling activities of a few oncogenic pathways, such as the phosphatidylinositol 3-kinase/protein kinase B (PI3K/AKT) and mitogen-activated protein kinase/extracellular signal-regulated kinase (MAPK/ERK) pathways. Understanding how various factors regulate AHR-mediated carcinogenesis processes helps clinicians and scientists plan personalized therapeutic strategies to improve anti-cancer treatment efficacy. As many studies that have reported the roles of AHR in regulating carcinogenesis are preclinical or observational clinical studies that did not explore the detailed mechanisms of how different chemical, physical, or biological agents promote AHR-mediated carcinogenesis processes, future studies should focus on conducting large-scale and functional studies to unravel the underlying mechanism of how AHR interacts with different factors in regulating carcinogenesis processes.

Patterns of item nonresponse behaviour to survey questionnaires are systematic and associated with genetic loci

Response to survey questionnaires is vital for social and behavioural research, and most analyses assume full and accurate response by participants. However, nonresponse is common and impedes proper interpretation and generalizability of results. We examined item nonresponse behaviour across 109 questionnaire items in the UK Biobank (N = 360,628). Phenotypic factor scores for two participant-selected nonresponse answers, 'Prefer not to answer' (PNA) and 'I don't know' (IDK), each predicted participant nonresponse in follow-up surveys (incremental pseudo-R2 = 0.056), even when controlling for education and self-reported health (incremental pseudo-R2 = 0.046). After performing genome-wide association studies of our factors, PNA and IDK were highly genetically correlated with one another (rg = 0.73 (s.e. = 0.03)) and with education (rg,PNA = -0.51 (s.e. = 0.03); rg,IDK = -0.38 (s.e. = 0.02)), health (rg,PNA = 0.51 (s.e. = 0.03); rg,IDK = 0.49 (s.e. = 0.02)) and income (rg,PNA = -0.57 (s.e. = 0.04); rg,IDK = -0.46 (s.e. = 0.02)), with additional unique genetic associations observed for both PNA and IDK (P < 5 × 10-8). We discuss how these associations may bias studies of traits correlated with item nonresponse and demonstrate how this bias may substantially affect genome-wide association studies. While the UK Biobank data are deidentified, we further protected participant privacy by avoiding exploring non-response behaviour to single questions, assuring that no information can be used to associate results with any particular respondents.

Advances in cutaneous squamous cell carcinoma

Human malignancies arise predominantly in tissues of epithelial origin, where the stepwise transformation from healthy epithelium to premalignant dysplasia to invasive neoplasia involves sequential dysregulation of biological networks that govern essential functions of epithelial homeostasis. Cutaneous squamous cell carcinoma (cSCC) is a prototype epithelial malignancy, often with a high tumour mutational burden. A plethora of risk genes, dominated by UV-induced sun damage, drive disease progression in conjunction with stromal interactions and local immunomodulation, enabling continuous tumour growth. Recent studies have identified subpopulations of SCC cells that specifically interact with the tumour microenvironment. These advances, along with increased knowledge of the impact of germline genetics and somatic mutations on cSCC development, have led to a greater appreciation of the complexity of skin cancer pathogenesis and have enabled progress in neoadjuvant immunotherapy, which has improved pathological complete response rates. Although measures for the prevention and therapeutic management of cSCC are associated with clinical benefit, the prognosis remains poor for advanced disease. Elucidating how the genetic mechanisms that drive cSCC interact with the tumour microenvironment is a current focus in efforts to understand, prevent and treat cSCC.

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.

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.

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.

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.

Interpretation of the role of germline and somatic non-coding mutations in cancer: expression and chromatin conformation informed analysis

BACKGROUND

There has been extensive scrutiny of cancer driving mutations within the exome (especially amino acid altering mutations) as these are more likely to have a clear impact on protein functions, and thus on cell biology. However, this has come at the neglect of systematic identification of regulatory (non-coding) variants, which have recently been identified as putative somatic drivers and key germline risk factors for cancer development. Comprehensive understanding of non-coding mutations requires understanding their role in the disruption of regulatory elements, which then disrupt key biological functions such as gene expression.

MAIN BODY

We describe how advancements in sequencing technologies have led to the identification of a large number of non-coding mutations with uncharacterized biological significance. We summarize the strategies that have been developed to interpret and prioritize the biological mechanisms impacted by non-coding mutations, focusing on recent annotation of cancer non-coding variants utilizing chromatin states, eQTLs, and chromatin conformation data.

CONCLUSION

We believe that a better understanding of how to apply different regulatory data types into the study of non-coding mutations will enhance the discovery of novel mechanisms driving cancer.

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.

Natural Killer cells demonstrate distinct eQTL and transcriptome-wide disease associations, highlighting their role in autoimmunity

Natural Killer cells are innate lymphocytes with central roles in immunosurveillance and are implicated in autoimmune pathogenesis. The degree to which regulatory variants affect Natural Killer cell gene expression is poorly understood. Here we perform expression quantitative trait locus mapping of negatively selected Natural Killer cells from a population of healthy Europeans (n = 245). We find a significant subset of genes demonstrate expression quantitative trait loci specific to Natural Killer cells and these are highly informative of human disease, in particular autoimmunity. A Natural Killer cell transcriptome-wide association study across five common autoimmune diseases identifies further novel associations at 27 genes. In addition to these cis observations, we find novel master-regulatory regions impacting expression of trans gene networks at regions including 19q13.4, the Killer cell Immunoglobulin-like Receptor region, GNLY, MC1R and UVSSA. Our findings provide new insights into the unique biology of Natural Killer cells, demonstrating markedly different expression quantitative trait loci from other immune cells, with implications for disease mechanisms.

Natural Killer cells are key mediators of anti-tumour immunosurveillance and anti-viral immunity. Here, the authors map regulatory genetic variation in primary Natural Killer cells, providing new insights into their role in human health and disease.

Benzo[a]pyrene—Environmental Occurrence, Human Exposure, and Mechanisms of Toxicity

Benzo[a]pyrene (B[a]P) is the main representative of polycyclic aromatic hydrocarbons (PAHs), and has been repeatedly found in the air, surface water, soil, and sediments. It is present in cigarette smoke as well as in food products, especially when smoked and grilled. Human exposure to B[a]P is therefore common. Research shows growing evidence concerning toxic effects induced by this substance. This xenobiotic is metabolized by cytochrome P450 (CYP P450) to carcinogenic metabolite: 7β,8α-dihydroxy-9α,10α-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene (BPDE), which creates DNA adducts, causing mutations and malignant transformations. Moreover, B[a]P is epigenotoxic, neurotoxic, and teratogenic, and exhibits pro-oxidative potential and causes impairment of animals’ fertility. CYP P450 is strongly involved in B[a]P metabolism, and it is simultaneously expressed as a result of the association of B[a]P with aromatic hydrocarbon receptor (AhR), playing an essential role in the cancerogenic potential of various xenobiotics. In turn, polymorphism of CYP P450 genes determines the sensitivity of the organism to B[a]P. It was also observed that B[a]P facilitates the multiplication of viruses, which may be an additional problem with the widespread COVID-19 pandemic. Based on publications mainly from 2017 to 2022, this paper presents the occurrence of B[a]P in various environmental compartments and human surroundings, shows the exposure of humans to this substance, and describes the mechanisms of its toxicity.

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.

Rosmarinus officinalis L. Leaf Extracts and Their Metabolites Inhibit the Aryl Hydrocarbon Receptor (AhR) Activation In Vitro and in Human Keratinocytes: Potential Impact on Inflammatory Skin Diseases and Skin Cancer

Aryl hydrocarbon receptor (AhR) activation by environmental agents and microbial metabolites is potentially implicated in a series of skin diseases. Hence, it would be very important to identify natural compounds that could inhibit the AhR activation by ligands of microbial origin as 6-formylindolo[3,2-b]carbazole (FICZ), indirubin (IND) and pityriazepin (PZ) or the prototype ligand 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Five different dry Rosmarinus officinalis L. extracts (ROEs) were assayed for their activities as antagonists of AhR ligand binding with guinea pig cytosol in the presence of [³H]TCDD. The methanolic ROE was further assayed towards CYP1A1 mRNA induction using RT-PCR in human keratinocytes against TCDD, FICZ, PZ, and IND. The isolated metabolites, carnosic acid, carnosol, 7-O-methyl-epi-rosmanol, 4′,7-O-dimethylapigenin, and betulinic acid, were assayed for their agonist and antagonist activity in the presence and absence of TCDD using the gel retardation assay (GRA). All assayed ROE extracts showed similar dose-dependent activities with almost complete inhibition of AhR activation by TCDD at 100 ppm. The methanol ROE at 10 ppm showed 99%, 50%, 90%, and 85% inhibition against TCDD, FICZ, IND, and PZ, respectively, in human keratinocytes. Most assayed metabolites exhibited dose-dependent antagonist activity. ROEs inhibit AhR activation by TCDD and by the Malassezia metabolites FICZ, PZ, and IND. Hence, ROE could be useful for the prevention or treatment of skin diseases mediated by activation of AhR.

The Aryl Hydrocarbon Receptor in the Pathogenesis of Environmentally-Induced Squamous Cell Carcinomas of the Skin

Cutaneous squamous cell carcinoma (SCC) is one of the most frequent malignancies in humans and academia as well as public authorities expect a further increase of its incidence in the next years. The major risk factor for the development of SCC of the general population is the repeated and unprotected exposure to ultraviolet (UV) radiation. Another important risk factor, in particular with regards to occupational settings, is the chronic exposure to polycyclic aromatic hydrocarbons (PAH) which are formed during incomplete combustion of organic material and thus can be found in coal tar, creosote, bitumen and related working materials. Importantly, both exposomal factors unleash their carcinogenic potential, at least to some extent, by activating the aryl hydrocarbon receptor (AHR). The AHR is a ligand-dependent transcription factor and key regulator in xenobiotic metabolism and immunity. The AHR is expressed in all cutaneous cell-types investigated so far and maintains skin integrity. We and others have reported that in response to a chronic exposure to environmental stressors, in particular UV radiation and PAHs, an activation of AHR and downstream signaling pathways critically contributes to the development of SCC. Here, we summarize the current knowledge about AHR's role in skin carcinogenesis and focus on its impact on defense mechanisms, such as DNA repair, apoptosis and anti-tumor immune responses. In addition, we discuss the possible consequences of a simultaneous exposure to different AHR-stimulating environmental factors for the development of cutaneous SCC.

Exposome and Skin. Part 2. The Influential Role of the Exposome, Beyond UVR, in Actinic Keratosis, Bowen's Disease and Squamous Cell Carcinoma: A Proposal

Actinic keratosis (AK) is the main risk factor for the development of cutaneous invasive squamous cell carcinoma (SCC). It represents the first sign of severe chronic ultraviolet radiation exposure, which has a clear significant effect. Nevertheless, the skin is exposed to many other exposome factors which should be thoroughly considered. Our aim was to assess the impact of exposome factors other than ultraviolet radiation (UVR) on the etiopathology of AK and Bowen's disease (BD) and progression of AK to SCC and to design tailored prevention strategies. We performed an exhaustive literature search in September 2021 through PubMed on the impact of exposome factors other than UVR on AK, BD and SCC. We conducted several parallel searches combining terms of the following topics: AK, BD, SCC and microbiome, hormones, nutrition, alcohol, tobacco, viral infections, chemical contaminants and air pollution. Notably, skin microbiome studies have shown how Staphylococcus aureus infections are associated with AK and AK-to-SCC progression by the production of chronic inflammation. Nutritional studies have demonstrated how a caloric restriction in fat intake, oral nicotinamide and moderate consumption of wine significantly reduce the number of premalignant keratoses and SCC. Regarding lifestyle factors, both alcohol and smoking are associated with the development of SCC in a dose-dependent manner. Relevant environmental factors are viral infections and chemical contaminants. Human papillomavirus infections induce deregulation of cellular proliferation and are associated with AK, BD and SCC. In addition to outdoor jobs, occupations such as industrial processing and farming also increase the risk of developing keratoses and SCC. The exposome of AK will undoubtedly help the understanding of its etiopathology and possible progression to SCC and will serve as a basis to design tailored prevention strategies.

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.

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.

A UVB-responsive common variant at chromosome band 7p21.1 confers tanning response and melanoma risk via regulation of the aryl hydrocarbon receptor, AHR

Genome-wide association studies (GWASs) have identified a melanoma-associated locus on chromosome band 7p21.1 with rs117132860 as the lead SNP and a secondary independent signal marked by rs73069846. rs117132860 is also associated with tanning ability and cutaneous squamous cell carcinoma (cSCC). Because ultraviolet radiation (UVR) is a key environmental exposure for all three traits, we investigated the mechanisms by which this locus contributes to melanoma risk, focusing on cellular response to UVR. Fine-mapping of melanoma GWASs identified four independent sets of candidate causal variants. A GWAS region-focused Capture-C study of primary melanocytes identified physical interactions between two causal sets and the promoter of the aryl hydrocarbon receptor (AHR). Subsequent chromatin state annotation, eQTL, and luciferase assays identified rs117132860 as a functional variant and reinforced AHR as a likely causal gene. Because AHR plays critical roles in cellular response to dioxin and UVR, we explored links between this SNP and AHR expression after both 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and ultraviolet B (UVB) exposure. Allele-specific AHR binding to rs117132860-G was enhanced following both, consistent with predicted weakened AHR binding to the risk/poor-tanning rs117132860-A allele, and allele-preferential AHR expression driven from the protective rs117132860-G allele was observed following UVB exposure. Small deletions surrounding rs117132860 introduced via CRISPR abrogates AHR binding, reduces melanocyte cell growth, and prolongs growth arrest following UVB exposure. These data suggest AHR is a melanoma susceptibility gene at the 7p21.1 risk locus and rs117132860 is a functional variant within a UVB-responsive element, leading to allelic AHR expression and altering melanocyte growth phenotypes upon exposure.

The Meaning of "Cause" in Genetics

Causation has multiple distinct meanings in genetics. One reason for this is meaning slippage between two concepts of the gene: Mendelian and molecular. Another reason is that a variety of genetic methods address different kinds of causal relationships. Some genetic studies address causes of traits in individuals, which can only be assessed when single genes follow predictable inheritance patterns that reliably cause a trait. A second sense concerns the causes of trait differences within a population. Whereas some single genes can be said to cause population-level differences, most often these claims concern the effects of many genes. Polygenic traits can be understood using heritability estimates, which estimate the relative influences of genetic and environmental differences to trait differences within a population. Attempts to understand the molecular mechanisms underlying polygenic traits have been developed, although causal inference based on these results remains controversial. Genetic variation has also recently been leveraged as a randomizing factor to identify environmental causes of trait differences. This technique-Mendelian randomization-offers some solutions to traditional epidemiological challenges, although it is limited to the study of environments with known genetic influences.

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.

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 .

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.

Extra-anogenital giant cutaneous squamous cell carcinomas require multidisciplinary management

Cutaneous squamous cell carcinoma (cSCC) is a highly prevalent cancer and the majority of cSCC have a good prognosis. However, a subset of cSCC can progress to advanced disease. We present the first reported case of a giant cSCC located on the breast. In addition, a systematic literature search on extra-anogenital giant (EAG) cSCC was performed using Pubmed and MEDLINE databases. Thirty-one articles could be retained which were relevant for this review. A total of 42 well-described cases were identified. Median age at presentation was 70 years (range 9-100 years). Twenty-four (57%) patients were male, eighteen were female (43%). The estimated median delay of treatment was 12 months (range 1 to >240 months). In 27 patients (64%) the giant cSCC was localized on the neck, face or scalp, 6 on the thoracic wall or back (14%), 4 on the lower limb (10%), 2 on the hip or buttock (5%), 2 on the upper limb (5%), one (2%) on the breast. Tumor stage at presentation was T2, T3 and T4 in respectively 26 (62%), 11 (26%) and 5 (12%) cases. Lymph node metastases were identified in 1 patient (2%) and distant metastases in another patient (2%). In 34 out of 42 cases (81%) primary radical surgical excision was performed, 3 received radiotherapy, 2 chemotherapy and 3 palliative care. In the cases with reported follow-up, four patients (4/30: 13%) died of disease. The treatment of EAG cSCC poses many problems, making a multidisciplinary approach of paramount importance.

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.

Neuropeptides and neurohormones in immune, inflammatory and cellular responses to ultraviolet radiation

Humans are exposed to varying amounts of ultraviolet radiation (UVR) through sunlight. UVR penetrates into human skin leading to release of neuropeptides, neurotransmitters and neuroendocrine hormones. These messengers released from local sensory nerves, keratinocytes, Langerhans cells (LCs), mast cells, melanocytes and endothelial cells (ECs) modulate local and systemic immune responses, mediate inflammation and promote differing cell biologic effects. In this review, we will focus on both animal and human studies that elucidate the roles of calcitonin gene-related peptide (CGRP), substance P (SP), nerve growth factor (NGF), nitric oxide and proopiomelanocortin (POMC) derivatives in mediating immune and inflammatory effects of exposure to UVR as well as other cell biologic effects of UVR exposure.

LINC00641 Inhibits the Development of Cutaneous Squamous Cell Carcinoma By Downregulating miR-424 in A431 Cells

Background: Cutaneous squamous cell carcinoma (CSCC) is the most deadly disease among nonmelanoma skin cancers. LINC00641 plays a role in various cancers, but its role in CSCC has not been reported so far. Methods and Materials: The expression of LINC00641 and miR-424 in cells was detected by RT-qPCR. CCK-8 and colony formation assay were used to detect the proliferation of cells. Western blot was used to detect the expression levels of proliferation-, invasion-, and migration-related proteins. Wound Healing and Transwell experiments detected the ability of cell invasion and migration. In animal experiments, a tumor-bearing model was established in nude mice, and tumor volume was measured and photographed. The expression levels of proliferation-, invasion-, and migration-related proteins were detected by western blot. Results: The expression of LINC00641 was significantly decreased in CSCC cell lines. The overexpression of LINC00641 at the cellular level inhibited the proliferation and migration of CSCC cell line A431 by downregulating the expression of miR-424. The overexpression of LINC00641 in animals inhibited the tumor volume of nude mice by downregulating the expression of miR-424 to inhibit the expression of proliferation- and migration-related proteins. Conclusion: LINC00641 inhibits the development of CSCC by downregulating miR-424.

Differential Molecular Expression Patterns Associated with Metastasis in Cutaneous Squamous Cell Carcinoma: A Systematic Review and Meta-Analysis

The majority of cutaneous squamous cell carcinomas are treated by surgical removal; however, approximately 4% of tumors will metastasize. Molecular expression testing may improve accuracy in estimating the prognosis and defining the mechanisms important in the disease progression, which may impact response to therapy. Using PubMed (MEDLINE) and EMBASE, a systematic review was performed to evaluate studies published from January 2005 to August 2019 reporting tumor protein or RNA expression along with either outcomes (metastasis or death) or a comparison of primary with metastatic tumor samples. Inclusion criteria were met by 45 studies containing 81 comparisons of 44 distinct proteins and 25 microRNAs. On meta-analysis of studies analyzing primary tumor samples in terms of later outcomes, high primary tumor expression of PD-L1 (OR = 2.34, 95% confidence interval = 1.09-5.02, P = 0.030), EGFR (OR = 2.57, 95% confidence interval = 1.24-5.33, P = 0.011), and podoplanin (OR = 2.33, 95% confidence interval = 1.00-5.41, P = 0.049) conferred increased odds for metastasis. In comparison, metastatic tissue was more likely to have a high expression of PD-L1 than primary tissue (OR = 3.13, 95% confidence interval = 1.00-9.75, P = 0.049). Further studies are needed to confirm whether testing for PD-L1, EGFR, and podoplanin expression aids in cutaneous squamous cell carcinomas prognostic estimation of metastasis or death or predicts response to therapy.

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.

No Evidence for a Causal Relationship Between Cancers and Parkinson's Disease

BACKGROUND

Epidemiological data suggest that cancer patients have a reduced risk of subsequent Parkinson's disease (PD) development, but the prevalence of PD in melanoma patients is often reported to be increased. Causal relationships between cancers and PD have not been fully explored.

OBJECTIVE

To study causal relationship between different cancers and PD.

METHODS

We used GWAS summary statistics of 15 different types of cancers and two-sample Mendelian randomization to study the causal relationship with PD.

RESULTS

There was no evidence to support a causal relationship between the studied cancers and PD. We also performed reverse analyses between PD and cancers with available full summary statistics (melanoma, breast, prostate, endometrial and keratinocyte cancers) and did not find evidence of causal relationship.

CONCLUSION

We found no evidence to support a causal relationship between cancers and PD and the previously reported associations could be a result of genetic pleiotropy, shared biology or biases.

A nomogram predicting subclinical extension of cutaneous squamous cell carcinoma in Chinese individuals: A retrospective study

Cutaneous squamous cell carcinoma usually extends beyond the visible margin. Little is known about the predictors for cutaneous squamous cell carcinoma with subclinical extension in Chinese individuals. This study aimed to construct a nomogram for predicting the probability of subclinical extension of cutaneous squamous cell carcinoma in Chinese patients.A retrospective analysis was conducted using data from Mohs micrographic surgery-treated cutaneous squamous cell carcinoma patients at a single institution between December 1, 2009 and October 31, 2019. Subclinical extension was defined as a lesion requiring ≥ 2 Mohs stages or with final safe margins of ≥ 5 mm. A nomogram predicting the probability of subclinical extension was constructed using the predictors identified in multivariable analysis.Of 274 patients included, 119 (43.4%) had subclinical extension. In multivariable analysis, male sex (odds ratio [OR], 2.45; 95% confidence interval [CI], 1.40-4.29; P = .002), lesions on mucocutaneous areas (OR, 3.71; 95% CI, 1.34-10.32; P = .012) and extremities (OR, 2.40; 95% CI, 1.20-4.78; P = .013), maximum diameter of 10 to 19 mm (OR, 14.15; 95% CI, 4.24-47.28; P < .001), and 20 to 29 mm (OR, 9.21; 95% CI, 2.80-30.29; P < .001) were associated with subclinical extension. A nomogram incorporating these 3 variables demonstrated promising predictive ability (C statistics = 0.78; 95% CI, 0.67-0.89).The nomogram incorporating sex, tumor location, and maximum diameter can provide individualized prediction for subclinical extension in Chinese patients with cutaneous squamous cell carcinoma. This information may help surgeons determine appropriate margins at the first Mohs stage.

A New Insight into the Potential Role of Tryptophan-Derived AhR Ligands in Skin Physiological and Pathological Processes

The aryl hydrocarbon receptor (AhR) plays a crucial role in environmental responses and xenobiotic metabolism, as it controls the transcription profiles of several genes in a ligand-specific and cell-type-specific manner. Various barrier tissues, including skin, display the expression of AhR. Recent studies revealed multiple roles of AhR in skin physiology and disease, including melanogenesis, inflammation and cancer. Tryptophan metabolites are distinguished among the groups of natural and synthetic AhR ligands, and these include kynurenine, kynurenic acid and 6-formylindolo[3,2-b]carbazole (FICZ). Tryptophan derivatives can affect and regulate a variety of signaling pathways. Thus, the interest in how these substances influence physiological and pathological processes in the skin is expanding rapidly. The widespread presence of these substances and potential continuous exposure of the skin to their biological effects indicate the important role of AhR and its ligands in the prevention, pathogenesis and progression of skin diseases. In this review, we summarize the current knowledge of AhR in skin physiology. Moreover, we discuss the role of AhR in skin pathological processes, including inflammatory skin diseases, pigmentation disorders and cancer. Finally, the impact of FICZ, kynurenic acid, and kynurenine on physiological and pathological processes in the skin is considered. However, the mechanisms of how AhR regulates skin function require further investigation.

The Role of Cell Adhesion Molecule 1 (CADM1) in Cutaneous Malignancies

Cell adhesion ability is one of the components to establish cell organization and shows a great contribution to human body construction consisting of various types of cells mixture to orchestrate tissue specific function. The cell adhesion molecule 1 (CADM1) is a molecule of cell adhesion with multiple functions and has been identified as a tumor suppressor gene. CADM1 has multifunctions on the pathogenesis of malignancies, and other normal cells such as immune cells. However, little is known about the function of CADM1 on cutaneous cells and cutaneous malignancies. CADM1 plays an important role in connecting cells with each other, contacting cells to deliver their signal, and acting as a scaffolding molecule for other immune cells to develop their immune responses. A limited number of studies reveal the contribution of CADM1 on the development of cutaneous malignancies. Solid cutaneous malignancies, such as cutaneous squamous cell carcinoma and malignant melanoma, reduce their CADM1 expression to promote the invasion and metastasis of the tumor. On the contrary to these cutaneous solid tumors except for Merkel cell carcinoma, cutaneous lymphomas, such as adult-T cell leukemia/lymphoma, mycosis fungoides, and Sézary syndrome, increase their CADM1 expression for the development of tumor environment. Based on the role of CADM1 in the etiology of tumor development, the theory of CADM1 contribution will desirably be applied to skin tumors according to other organ malignancies, however, the characteristics of skin as a multicomponent peripheral organ should be kept in mind to conclude their prognoses.

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.