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

Transcriptome-Wide Association Study Reveals New Molecular Interactions Associated with Melanoma Pathogenesis

A transcriptome-wide association study (TWAS) was conducted on genome-wide association study (GWAS) summary statistics of malignant melanoma of skin (UK Biobank dataset) and The Cancer Genome Atlas-Skin Cutaneous Melanoma (TCGA-SKCM) gene expression weights to identify melanoma susceptibility genes. The GWAS included 2465 cases and 449,799 controls, while the gene expression testing was conducted on 103 cases. Afterward, a gene enrichment analysis was applied to identify significant TWAS associations. The melanoma's gene-microRNA (miRNA) regulatory network was constructed from the TWAS genes and their corresponding miRNAs. At last, a disease enrichment analysis was conducted on the corresponding miRNAs. The TWAS detected 27 genes associated with melanoma with p-values less than 0.05 (the top three genes are LOC389458 (RBAK), C16orf73 (MEIOB), and EIF3CL). After the joint/conditional test, one gene (AMIGO1) was dropped, resulting in 26 significant genes. The Gene Ontology (GO) biological process associated the extended gene set (76 genes) with protein K11-linked ubiquitination and regulation of cell cycle phase transition. K11-linked ubiquitin chains regulate cell division. Interestingly, the extended gene set was related to different skin cancer subtypes. Moreover, the enriched pathways were nsp1 from SARS-CoV-2 that inhibit translation initiation in the host cell, cell cycle, translation factors, and DNA repair pathways full network. The gene-miRNA regulatory network identified 10 hotspot genes with the top three: TP53, BRCA1, and MDM2; and four hotspot miRNAs: mir-16, mir-15a, mir-125b, and mir-146a. Melanoma was among the top ten diseases associated with the corresponding (106) miRNAs. Our results shed light on melanoma pathogenesis and biologically significant molecular interactions.

Sleep, physical activity, and sedentary behaviors in relation to overall cancer and site-specific cancer risk: A prospective cohort study

Large prospective studies are required to better elucidate the associations of physical activity, sedentary behaviors (SBs), and sleep with overall cancer and site-specific cancer risk, accounting for the interactions with genetic predisposition. The study included 360,271 individuals in UK Biobank. After a median follow-up of 12.52 years, we found higher total physical activity (TPA) level and higher sleep scores were related to reduced risk of cancer while higher SB level showed a positive association with cancer. Compared with high TPA-healthy sleep group and low SB-healthy sleep group, low TPA-poor sleep group and high SB-poor sleep group had the highest risk for overall cancer, breast cancer, and lung cancer. Adherence to a more active exercise pattern was associated with a lower risk of cancer irrespective of genetic risk. Our study suggests that improving the quality of sleep and developing physical activity habits might yield benefits in mitigating the cancer risk.

Variant ranking pipeline for complex familial disorders

Identifying genetic susceptibility factors for complex disorders remains a challenging task. To analyze collections of small and large pedigrees where genetic heterogeneity is likely, but biological commonalities are plausible, we have developed a weights-based pipeline to prioritize variants and genes. The Weights-based vAriant Ranking in Pedigrees (WARP) pipeline prioritizes variants using 5 weights: disease incidence rate, number of cases in a family, genome fraction shared amongst cases in a family, allele frequency and variant deleteriousness. Weights, except for the population allele frequency weight, are normalized between 0 and 1. Weights are combined multiplicatively to produce family-specific-variant weights that are then averaged across all families in which the variant is observed to generate a multifamily weight. Sorting multifamily weights in descending order creates a ranked list of variants and genes for further investigation. WARP was validated using familial melanoma sequence data from the European Genome-phenome Archive. The pipeline identified variation in known germline melanoma genes POT1, MITF and BAP1 in 4 out of 13 families (31%). Analysis of the other 9 families identified several interesting genes, some of which might have a role in melanoma. WARP provides an approach to identify disease predisposing genes in studies with small and large pedigrees.

Pathogenic mechanisms and etiologic aspects of Mycobacterium avium subspecies paratuberculosis as an infectious cause of cutaneous melanoma

Infectious etiologies have previously been proposed as causes of both melanoma and non-melanoma skin cancer. This exploratory overview explains and presents the evidence for the hypothesis that a microorganism excreted in infected ruminant animal feces, Mycobacterium avium subspecies paratuberculosis (MAP), is the cause of some cases of cutaneous melanoma (CM). Occupational, residential, and recreational contact with MAP-contaminated feces, soil, sand, and natural bodies of water may confer a higher rate of CM. Included in our hypothesis are possible reasons for the differing rates and locations of CM in persons with white versus nonwhite skin, why CM develops underneath nails and in vulvar skin, why canine melanoma is an excellent model for human melanoma, and why the Bacille Calmette-Guérin (BCG) vaccine has demonstrated efficacy in the prevention and treatment of CM. The pathogenic mechanisms and etiologic aspects of MAP, as a transmittable agent underlying CM risk, are carefully deliberated in this paper. Imbalances in gut and skin bacteria, genetic risk factors, and vaccine prevention/therapy are also discussed, while acknowledging that the evidence for a causal association between MAP exposure and CM remains circumstantial.

The potential impact of dietary choices on melanoma risk: an anti-inflammatory diet

The role of inflammation in the aetiology of cancer is recognized. However, no study yet examined the association between an anti-inflammatory diet and cutaneous melanoma and explored whether it could be modified by genetic variations in cyclooxygenase-2 (COX-2), a key enzyme in inflammation. A case-control study was conducted in the IDI-IRCCS hospital in Rome, Italy with 273 cases of primary cutaneous melanoma and 269 controls frequency matched to cases. Information on socio-demographic and pigmentary characteristics, medical history, sun exposure and dietary habits were collected for all subjects. The - 765G > C polymorphism was identified in DNA extracted from blood samples. An anti-inflammatory diet score was created. Logistic regression models were fitted to obtain odds ratios (ORs) and 95% confidence intervals (CIs). A high anti-inflammatory diet score (≥ 8 anti-inflammatory dietary items) was associated with a decreased risk of cutaneous melanoma (OR: 0.29; 95%CI: 0.17-0.49, Ptrend < 0.0001) after adjusting for sex, age, education, number of common nevi, skin photo-type, solar lentigines and sunburns in childhood. COX-2 -765 G > C polymorphism was not an independent risk factor for cutaneous melanoma. Although interaction between - 765G > C genotypes and anti-inflammatory diet score was not statistically significant (p = 0.25), when stratified by -765 G > C genotypes the effect of the anti-inflammatory diet was slightly more pronounced for participants carrying - 765GG (OR: 0.17; 95%CI: 0.06-0.47, Ptrend < 0.001). Our study findings suggest that adherence to an anti-inflammatory diet is associated with a decreased risk of developing cutaneous melanoma. These results suggest the potential impact of dietary choices on melanoma risk.

iGWAS: Image-based genome-wide association of self-supervised deep phenotyping of retina fundus images

Existing imaging genetics studies have been mostly limited in scope by using imaging-derived phenotypes defined by human experts. Here, leveraging new breakthroughs in self-supervised deep representation learning, we propose a new approach, image-based genome-wide association study (iGWAS), for identifying genetic factors associated with phenotypes discovered from medical images using contrastive learning. Using retinal fundus photos, our model extracts a 128-dimensional vector representing features of the retina as phenotypes. After training the model on 40,000 images from the EyePACS dataset, we generated phenotypes from 130,329 images of 65,629 British White participants in the UK Biobank. We conducted GWAS on these phenotypes and identified 14 loci with genome-wide significance (p<5×10-8 and intersection of hits from left and right eyes). We also did GWAS on the retina color, the average color of the center region of the retinal fundus photos. The GWAS of retina colors identified 34 loci, 7 are overlapping with GWAS of raw image phenotype. Our results establish the feasibility of this new framework of genomic study based on self-supervised phenotyping of medical images.

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.

A protocol for annotation of total body photography for machine learning to analyze skin phenotype and lesion classification

Introduction

Artificial Intelligence (AI) has proven effective in classifying skin cancers using dermoscopy images. In experimental settings, algorithms have outperformed expert dermatologists in classifying melanoma and keratinocyte cancers. However, clinical application is limited when algorithms are presented with 'untrained' or out-of-distribution lesion categories, often misclassifying benign lesions as malignant, or misclassifying malignant lesions as benign. Another limitation often raised is the lack of clinical context (e.g., medical history) used as input for the AI decision process. The increasing use of Total Body Photography (TBP) in clinical examinations presents new opportunities for AI to perform holistic analysis of the whole patient, rather than a single lesion. Currently there is a lack of existing literature or standards for image annotation of TBP, or on preserving patient privacy during the machine learning process.

Methods

This protocol describes the methods for the acquisition of patient data, including TBP, medical history, and genetic risk factors, to create a comprehensive dataset for machine learning. 500 patients of various risk profiles will be recruited from two clinical sites (Australia and Spain), to undergo temporal total body imaging, complete surveys on sun behaviors and medical history, and provide a DNA sample. This patient-level metadata is applied to image datasets using DICOM labels. Anonymization and masking methods are applied to preserve patient privacy. A two-step annotation process is followed to label skin images for lesion detection and classification using deep learning models. Skin phenotype characteristics are extracted from images, including innate and facultative skin color, nevi distribution, and UV damage. Several algorithms will be developed relating to skin lesion detection, segmentation and classification, 3D mapping, change detection, and risk profiling. Simultaneously, explainable AI (XAI) methods will be incorporated to foster clinician and patient trust. Additionally, a publicly released dataset of anonymized annotated TBP images will be released for an international challenge to advance the development of new algorithms using this type of data.

Conclusion

The anticipated results from this protocol are validated AI-based tools to provide holistic risk assessment for individual lesions, and risk stratification of patients to assist clinicians in monitoring for skin cancer.

Genetic predisposition to childhood obesity does not influence the risk of developing skin cancer in adulthood

The relationship between body mass index (BMI) and melanoma and other skin cancers remains unclear. The objective of this study was to employ the Mendelian randomization (MR) approach to evaluate the effects of genetically predicted childhood adiposity on the risk of developing skin cancer later in life. Two-sample MR analyses were conducted using summary data from genome-wide association study (GWAS) meta-analyses of childhood BMI, melanoma, cutaneous squamous cell carcinoma (cSCC), and basal cell carcinoma (BCC). We used the inverse-variance-weighted (IVW) methods to obtain a pooled estimate across all genetic variants for childhood BMI. We performed multiple sensitivity analyses to evaluate the potential influence of various assumptions on our findings. We found no evidence that genetically predicted childhood BMI was associated with risks of developing melanoma, cSCC, or BCC in adulthood (OR, 95% CI: melanoma: 1.02 (0.93-1.13), cSCC 0.94 (0.79-1.11), BCC 0.97 (0.84-1.12)). Our findings do not support the conclusions from observational studies that childhood BMI is associated with increased risks of melanoma, cSCC, or BCC in adulthood. Intervening on childhood adiposity will not reduce the risk of common skin cancers later in life.

Genetic Analysis of Perceived Youthfulness Reveals Differences in How Men's and Women's Age Is Assessed

Skin aging is a natural process that occurs over time but can be accelerated by sun exposure. Measuring skin age in a large population can provide insight into the extent of skin damage from sun exposure and skin cancer risk. Understanding the genetics of skin aging, within and across sexes (males and females), could improve our understanding of the genetic drivers of both skin aging and skin cancer. We used UK Biobank data to examine the genetic overlap between perceived youthfulness and traits relevant to actinic photoaging. Our GWAS identified 22 genome-wide significant loci for women and 43 for men. The genetic correlation (rg) between perceived youthfulness in men and women was significantly less than unity (rg = 0.75, 95% confidence interval = 0.69-0.80), suggesting a gene-by-sex interaction. In women, perceived youthfulness was modestly correlated with keratinocyte cancer (rg = -0.19) and skin tanning (rg = 0.18). In men, perceived youthfulness was correlated with male-pattern baldness (rg = -0.23). This suggests that the genetic architecture of perceived youthfulness may differ between sexes, with genes influencing skin tanning and skin cancer susceptibility driving the difference in women, whereas genes influencing male-pattern baldness and other puberty-related traits drive the difference in men. We recommend that future genetic analysis of skin aging include a sex-stratified component.

The Skin Molecular Ecosystem Holds the Key to Nevogenesis and Melanomagenesis

Early detection of melanoma is critical to good patient outcomes, but we still know little about the mechanisms of early melanoma development. Normal epidermis has many of the sequence variants and genetic architecture disruptions found in both benign nevi, melanomas, and other skin cancers, yet continues to behave more or less normally. One hypothesis is that many melanocytes in this context are "tumor competent" but are regulated by the microenvironment provided by the surrounding keratinocytes to inhibit progress to nevi or melanoma. There is evidence of accumulating disorder in several measures of the genomic and epigenomic landscape from normal skin through nevi to melanoma that may be key to promoting nevogenesis and melanomagenesis.

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.

Novel ancestry-specific primary open-angle glaucoma loci and shared biology with vascular mechanisms and cell proliferation

Primary open-angle glaucoma (POAG), a leading cause of irreversible blindness globally, shows disparity in prevalence and manifestations across ancestries. We perform meta-analysis across 15 biobanks (of the Global Biobank Meta-analysis Initiative) (n = 1,487,441: cases = 26,848) and merge with previous multi-ancestry studies, with the combined dataset representing the largest and most diverse POAG study to date (n = 1,478,037: cases = 46,325) and identify 17 novel significant loci, 5 of which were ancestry specific. Gene-enrichment and transcriptome-wide association analyses implicate vascular and cancer genes, a fifth of which are primary ciliary related. We perform an extensive statistical analysis of SIX6 and CDKN2B-AS1 loci in human GTEx data and across large electronic health records showing interaction between SIX6 gene and causal variants in the chr9p21.3 locus, with expression effect on CDKN2A/B. Our results suggest that some POAG risk variants may be ancestry specific, sex specific, or both, and support the contribution of genes involved in programmed cell death in POAG pathogenesis.

Active natural compounds perturb the melanoma risk-gene network

Cutaneous melanoma is an aggressive type of skin cancer with a complex genetic landscape caused by the malignant transformation of melanocytes. This study aimed at providing an in silico network model based on the systematic profiling of the melanoma-associated genes considering germline mutations, somatic mutations, and genome-wide association study signals accounting for a total of 232 unique melanoma risk genes. A protein-protein interaction network was constructed using the melanoma risk genes as seeds and evaluated to describe the functional landscape in which the melanoma genes operate within the cellular milieu. Not only were the majority of the melanoma risk genes able to interact with each other at the protein level within the core of the network, but this showed significant enrichment for genes whose expression is altered in human melanoma specimens. Functional annotation showed the melanoma risk network to be significantly associated with processes related to DNA metabolism and telomeres, DNA damage and repair, cellular ageing, and response to radiation. We further explored whether the melanoma risk network could be used as an in silico tool to predict the efficacy of anti-melanoma phytochemicals, that are considered active molecules with potentially less systemic toxicity than classical cytotoxic drugs. A significant portion of the melanoma risk network showed differential expression when SK-MEL-28 human melanoma cells were exposed to the phytochemicals harmine and berberine chloride. This reinforced our hypothesis that the network modeling approach not only provides an alternative way to identify molecular pathways relevant to disease but it may also represent an alternative screening approach to prioritize potentially active compounds.

Causal relationships between risk of venous thromboembolism and 18 cancers: a bidirectional Mendelian randomization analysis

BACKGROUND

People with cancer experience high rates of venous thromboembolism (VTE). Risk of subsequent cancer is also increased in people experiencing their first VTE. The causal mechanisms underlying this association are not completely understood, and it is unknown whether VTE is itself a risk factor for cancer.

METHODS

We used data from large genome-wide association study meta-analyses to perform bidirectional Mendelian randomization analyses to estimate causal associations between genetic liability to VTE and risk of 18 different cancers.

RESULTS

We found no conclusive evidence that genetic liability to VTE was causally associated with an increased incidence of cancer, or vice versa. We observed an association between liability to VTE and pancreatic cancer risk [odds ratio for pancreatic cancer: 1.23 (95% confidence interval: 1.08-1.40) per log-odds increase in VTE risk, P = 0.002]. However, sensitivity analyses revealed this association was predominantly driven by a variant proxying non-O blood group, with inadequate evidence to suggest a causal relationship.

CONCLUSIONS

These findings do not support the hypothesis that genetic liability to VTE is a cause of cancer. Existing observational epidemiological associations between VTE and cancer are therefore more likely to be driven by pathophysiological changes which occur in the setting of active cancer and anti-cancer treatments. Further work is required to explore and synthesize evidence for these mechanisms.

Genetic testing for familial melanoma

While the average lifetime risk of melanoma worldwide is approximately 3%, those with inherited high-penetrance mutations face an increased lifetime risk of 52-84%. In countries of low melanoma incidence, such as in Southern Europe, familial melanoma genetic testing may be warranted when there are two first degree relatives with a melanoma diagnosis. Testing criteria for high incidence countries such as USA, or with very-high incidence, such as Australia and New Zealand, would require a threshold of 3 to 4 affected family members. A mutation in the most common gene associated with familial melanoma, CDKN2A, is identified in approximately 10-40% of those meeting testing criteria. However, the use of multi-gene panels covering additional less common risk genes can significantly increase the diagnostic yield. Currently, genetic testing for familial melanoma is typically conducted by qualified genetic counsellors, however with increasing demand on testing services and high incidence rate in certain countries, a mainstream model should be considered. With appropriate training, dermatologists are well placed to identify high risk individuals and offer melanoma genetic test in dermatology clinics. Genetic testing should be given in conjunction with pre- and post-test consultation. Informed patient consent should cover possible results, the limitations and implications of testing including inconclusive results, and potential for genetic discrimination. Previous studies reporting on participant outcomes of genetic testing for familial melanoma have found significant improvements in both sun protective behavior and screening frequency in mutation carriers.

Natural hair color and skin cancers: A two-sample Mendelian randomization study

Previous observational studies have indicated an association between hair color and the risk of melanoma and keratinocyte skin cancer (KSC); however, different hair colors show inconsistent effects on skin cancers. Here, we conducted a two-sample Mendelian randomization (MR) study to evaluate the causal relationship between natural hair color and skin cancers by using 211 single nucleotide polymorphisms as genetic instruments from a genome-wide meta-analysis of 360,270 individuals of European ancestry. Light hair colors (red, blonde, and light brown) were associated with high levels of cutaneous melanoma (CM) and KSC (CM-inverse variance weighted [IVW] odds ratio [OR]-red: 1.034, 95% confidence interval [CI]: 1.025-1.044, P < 0.001; OR-blonde: 1.008, 95% CI: 1.003-1.014, P = 0.003; OR-light brown: 1.006, 95% CI: 1.002-1.011, P = 0.009; KSC-IVW OR-red: 1.078, 95% CI: 1.053-1.103, P < 0.001; OR-blonde: 1.024, 95% CI: 1.009-1.040, P = 0.002; OR-light brown: 1.018, 95% CI: 1.004-1.033, P = 0.01). However, dark brown hair showed an inverse causal relationship with skin cancers (CM IVW OR: 0.987, 95% CI: 0.984-0.990, P < 0.001; KSC IVW OR: 0.979, 95% CI: 0.970-0.988, P < 0.001). Black hair was associated with a decreased risk of KSC (IVW OR: 0.954, 95% CI: 0.913-0.997, P = 0.036) but showed no causal relationship with CM. The present study provides strong MR evidence of a causal association between hair color and skin cancer. Secondary MR analyses enhances result robustness by replicating findings, exploring gender-specific effects, and providing a more comprehensive understanding of the complex relationship between hair color and skin cancers. More large-scale MR studies or randomized controlled trials are required to further investigate the mechanisms of the association between hair color and skin cancers.

The ATM Ser49Cys Variant Effects ATM Function as a Regulator of Oncogene-Induced Senescence

An apical component of the cell cycle checkpoint and DNA damage repair response is the ataxia-telangiectasia mutated (ATM) Ser/Thr protein kinase. A variant of ATM, Ser49Cys (rs1800054; minor allele frequency = 0.011), has been associated with an elevated risk of melanoma development; however, the functional consequence of this variant is not defined. ATM-dependent signalling in response to DNA damage has been assessed in a panel of patient-derived lymphoblastoid lines and primary human melanocytic cell strains heterozygous for the ATM Ser49Cys variant allele. The ATM Ser49Cys allele appears functional for acute p53-dependent signalling in response to DNA damage. Expression of the variant allele did reduce the efficacy of oncogene expression in inducing senescence. These findings demonstrate that the ATM 146C>G Ser49Cys allele has little discernible effect on the acute response to DNA damage but has reduced function observed in the chronic response to oncogene over-expression. Analysis of melanoma, naevus and skin colour genomics and GWAS analyses have demonstrated no association of this variant with any of these outcomes. The modest loss of function detected suggest that the variant may act as a modifier of other variants of ATM/p53-dependent signalling.

Assessing the genetic risk of nodular melanoma using a candidate gene approach

BACKGROUND

Nodular melanoma (NM) is a challenge to diagnose early due to its rapid growth and more atypical clinical presentation, making it the largest contributor to melanoma mortality.

OBJECTIVES

Our study aim was to perform a rare-variant allele (RVA) analysis of whole-exome sequencing of patients with NM and non-NM (minor allele frequency ≤ 1% non-Finnish European) for a set of 500 candidate genes potentially implicated in melanoma.

METHODS

This study recruited 131 participants with NM and 194 with non-NM from South-east Queensland and patients with NM from Victoria to perform a comparative analysis of possible genetic differences or similarities between the two melanoma cohorts.

RESULTS

Phenotypic analysis revealed that a majority of patients diagnosed with NM were older males with a higher frequency of fair skin and red hair than is seen in the general population. The distribution of common melanoma polygenic risk scores was similar in patients with NM and non-NM, with over 28% in the highest quantile of scores. There was also a similar frequency of carriage of familial/high-penetrant melanoma gene and loss-of-function variants. We identified 39 genes by filtering 500 candidate genes based on the greatest frequency in NM compared with non-NM cases. The genes with RVAs of greatest frequency in NM included PTCH1, ARID2 and GHR. Rare variants in the SMO gene, which interacts with PTCH1 as ligand and receptor, were also identified, providing evidence that the Hedgehog pathway may contribute to NM risk. There was a cumulative effect in carrying multiple rare variants in the NM-associated genes. A 14.8-fold increased ratio for NM compared with non-NM was seen when two RVAs of the 39 genes were carried by a patient.

CONCLUSIONS

This study highlights the importance of considering frequency of RVA to identify those at risk of NM in addition to known high penetrance genes.

Genetic variants and haplotype structures of miRNA host genes in cancer and obesity

Cancer and obesity are two important public health problems. This study aimed to investigate the role of genetic variants and haplotypes of miRNA host genes in cancer and obesity. Data from the catalog of genome-wide association studies (GWAS) were used to find significant variants (index). Then, 1000-genome phase 3 data were used to find haplotypic variants (proxy) associated with these diseases. The candidate variants and haplotypes were identified from proxy and index variants. Finally, SNP function analysis was performed. All GWAS-significant cancer-associated miRNA host gene variants, including MIR4713HG, MIR663AHG, MIR99AHG and MIR4435-2HG, were also significantly associated with obesity. The rs703764 variant was common between cutaneous melanoma and obesity traits in the European population (P ≤ 5E-8). The rs2414098 variant was associated with endometrial cancer (P ≤ 5E-13), and the rs7173595 variant was associated with waist-hip ratio (P ≤ 5E-13) and new CGGCATCA haplotypic located at MIR4713HG was identified in the European population. In addition, the ATCTTGTT haplotype for rs17041868 in MIR4435-2HG was identified to be associated with obesity traits (waist-hip ratio and BMI) in the European population (P ≤ 5E-8). This study found that rs703764 is a common genetic marker between cancer and obesity. The CGGCATCA haplotype is common between endometrial cancer and waist-hip ratio. Also, ATCTTGTT haplotype is associated with obesity traits. These results indicate that the variants and haplotypes of miRNAs host genes play an important role between cancer and obesity in the European population. It is suggested to investigate the effect of these structures in other populations.Communicated by Ramaswamy H. Sarma.

FORGEdb: a tool for identifying candidate functional variants and uncovering target genes and mechanisms for complex diseases

The majority of disease-associated variants identified through genome-wide association studies are located outside of protein-coding regions. Prioritizing candidate regulatory variants and gene targets to identify potential biological mechanisms for further functional experiments can be challenging. To address this challenge, we developed FORGEdb ( https://forgedb.cancer.gov/ ; https://forge2.altiusinstitute.org/files/forgedb.html ; and https://doi.org/10.5281/zenodo.10067458 ), a standalone and web-based tool that integrates multiple datasets, delivering information on associated regulatory elements, transcription factor binding sites, and target genes for over 37 million variants. FORGEdb scores provide researchers with a quantitative assessment of the relative importance of each variant for targeted functional experiments.

ATM Variant as a Cause of Hereditary Cutaneous Melanoma in a Spanish Family: Case Report

Introduction

Ataxia-Telangiectasia Mutated (ATM) is a cancer predisposition gene; carriers of germline pathogenic variants have an increased risk of developing malignancies, including breast, prostate, pancreatic, and ovarian cancer. Most ATM variants are of uncertain significance. Findings from genome-wide association studies (GWAS) suggest that ATM may be a low-risk melanoma susceptibility locus.

Case Report

We report the case of a Hispanic family whose members who have presented cutaneous melanoma have been found to be carriers for the ATM pathogenic variant c.3747-1G>C (rs730881364), one of whom was diagnosed at 24 years old.

Discussion

We describe for the first time the possible clinical association between ATM (c.3747-1G>C) and familial melanoma. In silico splice site analysis predicts that this alteration will weaken the native splice acceptor site and will result in the creation or strengthening of a novel splice acceptor site, assuming a variant that entails loss of functionality that is probably pathogenic and related to oncogenesis. However, we cannot exclude that cutaneous melanoma in both members and at an early age is the result of chance, environmental interaction, other uncontrolled external factors, or the interaction of other genetic alterations other than the ATM variant described in this study.

Analysing the causal relationship between potentially protective and risk factors and cutaneous melanoma: A Mendelian randomization study

BACKGROUND

Previous observational studies reported altered melanoma risks in relation to many potential factors, such as coffee intake, smoking habits and photodamage-related conditions. Considering the susceptibility of epidemiological studies to residual confounders, there remains uncertainty about the actual causal roles of these reported factors in melanoma aetiology.

OBJECTIVES

This study aims to investigate the causal association between cutaneous melanoma (CM) and previously reported factors: coffee intake, alcohol consumption, lifetime smoking, socioeconomic status (SES), ease of skin tanning, childhood sunburn and facial ageing, providing insight into its underlying aetiology and preventative strategies.

METHODS

We utilized a two-sample MR analysis on data from the largest meta-analysis summary statistics of confirmed cutaneous melanoma including 30,134 patients. Genetic instrumental variables were constructed by identifying single nucleotide polymorphisms (SNPs) that associate with corresponding factors. Inverse variance weighted (IVW) was the primary MR method. For sensitivity and heterogeneity, MR Egger, weighted median, simple mode, weighted mode and MR Egger intercept tests were examined.

RESULTS

Cutaneous melanoma risks were found to be elevated in association with a predisposition towards ease of skin tanning (IVW: OR = 2.842, 95% CI 2.468-3.274, p < 0.001) and with childhood sunburn history (IVW: OR = 6.317, 95% CI 4.479-8.909, p < 0.001). Repeated MR after removing potential confounders and outliers demonstrated resolved horizontal pleiotropy and statistically significant results that closely mirrored the initial findings. Other potential factors, such as coffee intake, alcohol consumption, smoking and socioeconomic status (SES), indicated insignificant effects on melanoma risk in the analysis, and therefore, our Mendelian randomization study does not support their roles in modifying melanoma risks.

CONCLUSIONS

Our extensive MR analysis provides strong evidence of the causative role of ease of skin tanning and childhood sunburn history in elevating melanoma risk. Curtailing ultraviolet radiation (UVR) exposure may be the single best preventative strategy to reduce melanoma risk.

Acceptability and appropriateness of a risk-tailored organised melanoma screening program: Qualitative interviews with key informants

INTRODUCTION

In Australia, opportunistic screening (occurring as skin checks) for the early detection of melanoma is common, and overdiagnosis is a recognised concern. Risk-tailored cancer screening is an approach to cancer control that aims to provide personalised screening tailored to individual risk. This study aimed to explore the views of key informants in Australia on the acceptability and appropriateness of risk-tailored organised screening for melanoma, and to identify barriers, facilitators and strategies to inform potential future implementation. Acceptability and appropriateness are crucial, as successful implementation will require a change of practice for clinicians and consumers.

METHODS

This was a qualitative study using semi-structured interviews. Key informants were purposively selected to ensure expertise in melanoma early detection and screening, prioritising senior or executive perspectives. Consumers were expert representatives. Data were analysed deductively using the Tailored Implementation for Chronic Diseases (TICD) checklist.

RESULTS

Thirty-six participants were interviewed (10 policy makers; 9 consumers; 10 health professionals; 7 researchers). Key informants perceived risk-tailored screening for melanoma to be acceptable and appropriate in principle. Barriers to implementation included lack of trial data, reluctance for low-risk groups to not screen, variable skill level in general practice, differing views on who to conduct screening tests, confusing public health messaging, and competing health costs. Key facilitators included the perceived opportunity to improve health equity and the potential cost-effectiveness of a risk-tailored screening approach. A range of implementation strategies were identified including strengthening the evidence for cost-effectiveness, engaging stakeholders, developing pathways for people at low risk, evaluating different risk assessment criteria and screening delivery models and targeted public messaging.

CONCLUSION

Key informants were supportive in principle of risk-tailored melanoma screening, highlighting important next steps. Considerations around risk assessment, policy and modelling the costs of current verses future approaches will help inform possible future implementation of risk-tailored population screening for melanoma.

Identification of genomic-wide genetic links between cutaneous melanoma and obesity-related physical traits via cFDR

BACKGROUND

Both epidemiological and clinical studies have suggested the comorbidity between cutaneous melanoma (CM) and obesity-related physical traits. However, it remains unclear about their shared genetic architecture.

OBJECTIVE

To determine the shared genetic architecture between CM and obesity-related physical traits through conditional false discovery rate (cFDR) analysis.

METHOD

Quantile-quantile plots were firstly built to assess the pleiotropic enrichment of shared single nucleotide polymorphisms between CM and each trait. Then, cFDR and conjunctional cFDR (ccFDR) were used to identify the shared risk loci between CM and each trait. Moreover, the functional evaluation of shared risk genes was carried out through analyses of expression quantitative trait loci (eQTL), Kyoto Encyclopedia of Genes and Genomes and gene ontology, respectively. Finally, single-cell sequence analysis was performed to locate the expression of eQTL-mapped genes in tissues.

RESULTS

Successive pleiotropic enrichment was found between CM and 5 obesity-related traits or height. 24 shared risk loci were identified between CM and 13 traits except appendicular lean mass using ccFDR analysis, with 17 novel and 4 validated loci. The functions of ccFDR-identified and eQTL-mapped genes were revealed to be mainly involved in cellular senescence, proliferation, meiotic nuclear division, cell cycle, and the metabolism of lipid, cholesterol and glucose. Single-cell sequence analysis showed that keratinocytes contribute to the occurrence and aggressiveness of CM through secreting paracrine cytokines.

CONCLUSION

Our findings demonstrate the significant genetic correlation between CM and obesity-related physical traits, which may provide a novel genetical basis for the pathogenesis and treatment of CM.

Telomere length is associated with increased risk of cutaneous melanoma: a Mendelian randomization study

RESULTS

The MR analysis using two TL GWAS datasets revealed strong and consistent evidence that long TL is causally associated with an increased risk of CM. The analysis of the Codd et al. dataset found that long TL significantly predicted an elevated risk of CM (IVW OR = 2.411, 95% CI 2.092-2.780, P = 8.05E-34). Similarly, the analysis of the Li et al. dataset yielded consistent positive results across all MR methods, providing further robustness to the causal relationship (IVW OR = 2.324, 95% CI 1.516-3.565, P = 1.11E-04). The study provides evidence for a causal association between TL and CM susceptibility, indicating that longer TL increases the risk of developing CM and providing insight into the unique telomere biology in melanoma pathogenesis. Telomere maintenance pathways may be a potential target for preventing and treating CM.

Mendelian Randomization Analysis reveals Inverse Genetic Risks between Skin Cancers and Vitiligo

Several observational studies have demonstrated a consistent pattern of decreased melanoma risk among patients with vitiligo. More recently, this finding has been supported by a suggested genetic relationship between the two entities, with certain variants significantly associated with an increased risk of melanoma, basal cell carcinoma, and squamous cell carcinoma but a decreased risk of vitiligo. We compared 48 associated variants from a recently published GWAS and identified three variants-located in the TYR, MC1R-DEF8, and RALY-EIF2S2-ASIP-AHCY-ITCH loci- that correlated with an increased risk for melanoma, basal cell carcinoma, and squamous cell carcinoma and a decreased risk for vitiligo. We then used results of skin cancers and vitiligo GWAS to compare the shared genetic properties between these two traits through an unbiased Mendelian randomization analysis. Our results suggest that the inverse genetic relationship between common skin cancers and vitiligo is broader than previously reported owing to the influence of shared genome-wide significant associations.

Colitis-associated carcinogenesis: crosstalk between tumors, immune cells and gut microbiota

Colorectal cancer (CRC) is the third most common cancer worldwide. One of the main causes of colorectal cancer is inflammatory bowel disease (IBD), which includes ulcerative colitis (UC) and Crohn's disease (CD). Intestinal epithelial cells (IECs), intestinal mesenchymal cells (IMCs), immune cells, and gut microbiota construct the main body of the colon and maintain colon homeostasis. In the development of colitis and colitis-associated carcinogenesis, the damage, disorder or excessive recruitment of different cells such as IECs, IMCs, immune cells and intestinal microbiota play different roles during these processes. This review aims to discuss the various roles of different cells and the crosstalk of these cells in transforming intestinal inflammation to cancer, which provides new therapeutic methods for chemotherapy, targeted therapy, immunotherapy and microbial therapy.

TINF2 is a major susceptibility gene in Danish patients with multiple primary melanoma

TINF2 encodes the TINF2 protein, which is a subunit in the shelterin complex critical for telomere regulation. Three recent studies have associated six truncating germline variants in TINF2 that have previously been associated with a cancer predisposition syndrome (CPS) caused by elongation of the telomeres. This has added TINF2 to the long telomere syndrome genes, together with other telomere maintenance genes such as ACD, POT1, TERF2IP, and TERT. We report a clinical study of 102 Danish patients with multiple primary melanoma (MPM) in which a germline truncating variant in TINF2 (p.(Arg265Ter)) was identified in four unrelated participants. The telomere lengths of three variant carriers were >90% percentile. In a routine diagnostic setting, the variant was identified in two more families, including an additional MPM patient and monozygotic twins with thyroid cancer and other cancer types. A total of 10 individuals from six independent families were confirmed carriers, all with cancer history, predominantly melanoma. Our findings suggest a major role of TINF2 in Danish patients with MPM. In addition to melanoma, other cancers in the six families include thyroid, renal, breast, and sarcoma, supporting a CPS in which melanoma, thyroid cancer, and sarcoma predominate. Further studies are needed to establish the full spectrum of associated cancer types and characterize lifetime cancer risk in carriers.

Unveiling the Genetic Complexity of Teratozoospermia: Integrated Genomic Analysis Reveals Novel Insights into lncRNAs' Role in Male Infertility

Male infertility is a global health issue, affecting over 20 million men worldwide. Genetic factors are crucial in various male infertility forms, including teratozoospermia. Nonetheless, the genetic causes of male infertility remain largely unexplored. In this study, we employed whole-genome sequencing and RNA expression analysis to detect differentially expressed (DE) long-noncoding RNAs (lncRNAs) in teratozoospermia, along with mutations that are exclusive to teratozoospermic individuals within these DE lncRNAs regions. Bioinformatic tools were used to assess variants' impact on lncRNA structure, function, and lncRNA-miRNA interactions. Our analysis identified 1166 unique mutations in teratozoospermic men within DE lncRNAs, distinguishing them from normozoospermic men. Among these, 64 variants in 23 lncRNAs showed potential regulatory roles, 7 variants affected 4 lncRNA structures, while 37 variants in 17 lncRNAs caused miRNA target loss or gain. Pathway Enrichment and Gene Ontology analyses of the genes targeted by the affected miRNAs revealed dysregulated pathways in teratozoospermia and a link between male infertility and cancer. This study lists novel variants and lncRNAs associated for the first time with teratozoospermia. These findings pave the way for future studies aiming to enhance diagnosis and therapy in the field of male infertility.

Genetic variants for smoking behaviour and risk of skin cancer

Observational studies have suggested that smoking may increase the risk of cutaneous squamous cell carcinoma (cSCC) while decreasing the risks of basal cell carcinoma (BCC), and melanoma. However, it remains possible that confounding by other factors may explain these associations. The aim of this investigation was to use Mendelian randomization (MR) to test whether smoking is associated with skin cancer, independently of other factors. Two-sample MR analyses were conducted to determine the causal effect of smoking measures on skin cancer risk using genome-wide association study (GWAS) summary statistics. We used the inverse-variance-weighted estimator to derive separate risk estimates across genetic instruments for all smoking measures. A genetic predisposition to smoking initiation was associated with lower risks of all skin cancer types, although none of the effect estimates reached statistical significance (OR 95% CI BCC 0.91, 0.82-1.01; cSCC 0.82, 0.66-1.01; melanoma 0.91, 0.82-1.01). Results for other measures were similar to smoking initiation with the exception of smoking intensity which was associated with a significantly reduced risk of melanoma (OR 0.67, 95% CI 0.51-0.89). Our findings support the findings of observational studies linking smoking to lower risks of melanoma and BCC. However, we found no evidence that smoking is associated with an elevated risk of cSCC; indeed, our results are most consistent with a decreased risk, similar to BCC and melanoma.

Uncovering the complex relationship between balding, testosterone and skin cancers in men

Male-pattern baldness (MPB) is related to dysregulation of androgens such as testosterone. A previously observed relationship between MPB and skin cancer may be due to greater exposure to ultraviolet radiation or indicate a role for androgenic pathways in the pathogenesis of skin cancers. We dissected this relationship via Mendelian randomization (MR) analyses, using genetic data from recent male-only meta-analyses of cutaneous melanoma (12,232 cases; 20,566 controls) and keratinocyte cancers (KCs) (up to 17,512 cases; >100,000 controls), followed by stratified MR analysis by body-sites. We found strong associations between MPB and the risk of KC, but not with androgens, and multivariable models revealed that this relationship was heavily confounded by MPB single nucleotide polymorphisms involved in pigmentation pathways. Site-stratified MR analyses revealed strong associations between MPB with head and neck squamous cell carcinoma and melanoma, suggesting that sun exposure on the scalp, rather than androgens, is the main driver. Men with less hair covering likely explains, at least in part, the higher incidence of melanoma in men residing in countries with high ambient UV.

Elsinochrome A induces cell apoptosis and autophagy in photodynamic therapy

Elsinochrome A (EA) is a perylene quinone natural photosensitizer, photosensitizer under light excitation generates reactive oxygen species (ROS) to induce apoptosis, so can be used for treating tumors, that is so-called photodynamic therapy (PDT). However, the molecular mechanism, especially related to apoptosis and autophagy, is still unclear. In this study, we aimed to explore the mechanism of EA-PDT-induced B16 cells apoptosis and autophagy. The action of EA-PDT on mitochondrial permeability transition pore (MPTP), mitochondrial membrane potential (MMP) and the mitochondrial function were researched by fluorescence technique and Extracellular Flux Analyzer. Illumina sequencing, tandem mass tags Quantitative Proteomics and Western Blot studied the mechanism at the gene and protein levels. The results indicated that EA-PDT had excellent phototoxicity in vitro. EA could bind to the mitochondria. EA-PDT for 5 min caused MPTP opening, MMP decreasing and abnormal mitochondrial function with a concentration-dependent characteristic. EA-PDT resulted in an increase intracellular ROS and the number of autophagosomes. Caspase2, caspase9 and tnf were upregulated, and bcl2, prkn, atg2, atg9 and atg10 were downregulated. Our results indicated that EA-PDT induced cell apoptosis and autophagy through the mediation of ROS/Atg/Parkin. This study can provide enlightenment for exploring potential targets of drug development for the PDT of melanoma.

Genomic characterisation of the overlap of endometriosis with 76 comorbidities identifies pleiotropic and causal mechanisms underlying disease risk

Comorbid conditions can be driven by underlying pleiotropic and causal mechanisms that can provide insights into shared molecular and biological processes contributing to disease risk. Endometriosis is a chronic condition affecting one in nine women of reproductive age and poses many challenges including lengthy diagnostic delays and limited treatment efficacy owing to poor understanding of disease aetiology. To shed light on the underlying biological mechanisms and to identify potential risk factors, we examine the epidemiological and genomic relationship between endometriosis and its comorbidities. In the UK Biobank 292 ICD10 codes were epidemiologically correlated with endometriosis diagnosis, including gynaecological, immune, infection, pain, psychiatric, cancer, gastrointestinal, urinary, bone and cardiovascular traits. A subset of the identified comorbidities (n = 76) underwent follow-up genetic analysis. Whilst Mendelian randomisation suggested causality was not responsible for most comorbid relationships, 22 traits were genetically correlated with endometriosis, including pain, gynaecological and gastrointestinal traits, suggestive of a shared genetic background. Pleiotropic genetic variants and genes were identified using gene-based and colocalisation analysis. Shared genetic risk factors and potential target genes suggest a diverse collection of biological systems are involved in these comorbid relationships including coagulation factors, development of the female reproductive tract and cell proliferation. These findings highlight the diversity of traits with epidemiological and genomic overlap with endometriosis and implicate a key role for pleiotropy in the comorbid relationships.

Preferential killing of melanoma cells by a p16-related peptide

We report the identification of a synthetic, cell-penetrating peptide able to kill human melanoma cells efficiently and selectively, while being less toxic to normal human melanocytes and nontoxic to human fibroblasts. The peptide is based on the target-binding site of the melanoma suppressor and senescence effector p16 (also known as INK4A or CDKN2A), coupled to a cell-penetrating moiety. The killing is by apoptosis and appears to act by a route other than the canonical downstream target of p16 and CDK4, the retinoblastoma (RB) protein family, as it is also effective in HeLa cells and a melanocyte line expressing HPV E7 oncogenes, which both lack any active RB. There was varying toxicity to other types of cancer cell lines, such as glioblastoma. Melanoma cell killing by a p16-derived peptide was reported once before but only at a higher concentration, while selectivity and generality were not previously tested.

Melanoma risk prediction based on a polygenic risk score and clinical risk factors

Melanoma is one of the most commonly diagnosed cancers in the Western world: third in Australia, fifth in the USA and sixth in the European Union. Predicting an individual's personal risk of developing melanoma may aid them in undertaking effective risk reduction measures. The objective of this study was to use the UK Biobank to predict the 10-year risk of melanoma using a newly developed polygenic risk score (PRS) and an existing clinical risk model. We developed the PRS using a matched case-control training dataset ( N  = 16 434) in which age and sex were controlled by design. The combined risk score was developed using a cohort development dataset ( N  = 54 799) and its performance was tested using a cohort testing dataset ( N  = 54 798). Our PRS comprises 68 single-nucleotide polymorphisms and had an area under the receiver operating characteristic curve of 0.639 [95% confidence interval (CI) = 0.618-0.661]. In the cohort testing data, the hazard ratio per SD of the combined risk score was 1.332 (95% CI = 1.263-1.406). Harrell's C-index was 0.685 (95% CI = 0.654-0.715). Overall, the standardized incidence ratio was 1.193 (95% CI = 1.067-1.335). By combining a PRS and a clinical risk score, we have developed a risk prediction model that performs well in terms of discrimination and calibration. At an individual level, information on the 10-year risk of melanoma can motivate people to take risk-reduction action. At the population level, risk stratification can allow more effective population-level screening strategies to be implemented.

Autoimmune alleles at the major histocompatibility locus modify melanoma susceptibility

Autoimmunity and cancer represent two different aspects of immune dysfunction. Autoimmunity is characterized by breakdowns in immune self-tolerance, while impaired immune surveillance can allow for tumorigenesis. The class I major histocompatibility complex (MHC-I), which displays derivatives of the cellular peptidome for immune surveillance by CD8+ T cells, serves as a common genetic link between these conditions. As melanoma-specific CD8+ T cells have been shown to target melanocyte-specific peptide antigens more often than melanoma-specific antigens, we investigated whether vitiligo- and psoriasis-predisposing MHC-I alleles conferred a melanoma-protective effect. In individuals with cutaneous melanoma from both The Cancer Genome Atlas (n = 451) and an independent validation set (n = 586), MHC-I autoimmune-allele carrier status was significantly associated with a later age of melanoma diagnosis. Furthermore, MHC-I autoimmune-allele carriers were significantly associated with decreased risk of developing melanoma in the Million Veteran Program (OR = 0.962, p = 0.024). Existing melanoma polygenic risk scores (PRSs) did not predict autoimmune-allele carrier status, suggesting these alleles provide orthogonal risk-relevant information. Mechanisms of autoimmune protection were neither associated with improved melanoma-driver mutation association nor improved gene-level conserved antigen presentation relative to common alleles. However, autoimmune alleles showed higher affinity relative to common alleles for particular windows of melanocyte-conserved antigens and loss of heterozygosity of autoimmune alleles caused the greatest reduction in presentation for several conserved antigens across individuals with loss of HLA alleles. Overall, this study presents evidence that MHC-I autoimmune-risk alleles modulate melanoma risk unaccounted for by current PRSs.

Population-based analysis of POT1 variants in a cutaneous melanoma case-control cohort

Pathogenic germline variants in the protection of telomeres 1 gene (POT1) have been associated with predisposition to a range of tumour types, including melanoma, glioma, leukaemia and cardiac angiosarcoma. We sequenced all coding exons of the POT1 gene in 2928 European-descent melanoma cases and 3298 controls, identifying 43 protein-changing genetic variants. We performed POT1-telomere binding assays for all missense and stop-gained variants, finding nine variants that impair or disrupt protein-telomere complex formation, and we further define the role of variants in the regulation of telomere length and complex formation through molecular dynamics simulations. We determine that POT1 coding variants are a minor contributor to melanoma burden in the general population, with only about 0.5% of melanoma cases carrying germline pathogenic variants in this gene, but should be screened in individuals with a strong family history of melanoma and/or multiple malignancies.

Explainable multi-task learning improves the parallel estimation of polygenic risk scores for many diseases through shared genetic basis

Many complex diseases share common genetic determinants and are comorbid in a population. We hypothesized that the co-occurrences of diseases and their overlapping genetic etiology can be exploited to simultaneously improve multiple diseases' polygenic risk scores (PRS). This hypothesis was tested using a multi-task learning (MTL) approach based on an explainable neural network architecture. We found that parallel estimations of the PRS for 17 prevalent cancers in a pan-cancer MTL model were generally more accurate than independent estimations for individual cancers in comparable single-task learning (STL) models. Such performance improvement conferred by positive transfer learning was also observed consistently for 60 prevalent non-cancer diseases in a pan-disease MTL model. Interpretation of the MTL models revealed significant genetic correlations between the important sets of single nucleotide polymorphisms used by the neural network for PRS estimation. This suggested a well-connected network of diseases with shared genetic basis.

Increase in power by obtaining 10 or more controls per case when type-1 error is small in large-scale association studies

BACKGROUND

The rule of thumb that there is little gain in statistical power by obtaining more than 4 controls per case, is based on type-1 error α = 0.05. However, association studies that evaluate thousands or millions of associations use smaller α and may have access to plentiful controls. We investigate power gains, and reductions in p-values, when increasing well beyond 4 controls per case, for small α.

METHODS

We calculate the power, the median expected p-value, and the minimum detectable odds-ratio (OR), as a function of the number of controls/case, as α decreases.

RESULTS

As α decreases, at each ratio of controls per case, the increase in power is larger than for α = 0.05. For α between 10-6 and 10-9 (typical for thousands or millions of associations), increasing from 4 controls per case to 10-50 controls per case increases power. For example, a study with power = 0.2 (α = 5 × 10-8) with 1 control/case has power = 0.65 with 4 controls/case, but with 10 controls/case has power = 0.78, and with 50 controls/case has power = 0.84. For situations where obtaining more than 4 controls per case provides small increases in power beyond 0.9 (at small α), the expected p-value can decrease by orders-of-magnitude below α. Increasing from 1 to 4 controls/case reduces the minimum detectable OR toward the null by 20.9%, and from 4 to 50 controls/case reduces by an additional 9.7%, a result which applies regardless of α and hence also applies to "regular" α = 0.05 epidemiology.

CONCLUSIONS

At small α, versus 4 controls/case, recruiting 10 or more controls/cases can increase power, reduce the expected p-value by 1-2 orders of magnitude, and meaningfully reduce the minimum detectable OR. These benefits of increasing the controls/case ratio increase as the number of cases increases, although the amount of benefit depends on exposure frequencies and true OR. Provided that controls are comparable to cases, our findings suggest greater sharing of comparable controls in large-scale association studies.

Cross-cancer pleiotropic analysis identifies three novel genetic risk loci for colorectal cancer

BACKGROUND

To understand the shared genetic basis between colorectal cancer (CRC) and other cancers and identify potential pleiotropic loci for compensating the missing genetic heritability of CRC.

METHODS

We conducted a systematic genome-wide pleiotropy scan to appraise associations between cancer-related genetic variants and CRC risk among European populations. Single nucleotide polymorphism (SNP)-set analysis was performed using data from the UK Biobank and the Study of Colorectal Cancer in Scotland (10 039 CRC cases and 30 277 controls) to evaluate the overlapped genetic regions for susceptibility of CRC and other cancers. The variant-level pleiotropic associations between CRC and other cancers were examined by CRC genome-wide association study meta-analysis and the pleiotropic analysis under composite null hypothesis (PLACO) pleiotropy test. Gene-based, co-expression and pathway enrichment analyses were performed to explore potential shared biological pathways. The interaction between novel genetic variants and common environmental factors was further examined for their effects on CRC.

RESULTS

Genome-wide pleiotropic analysis identified three novel SNPs (rs2230469, rs9277378 and rs143190905) and three mapped genes (PIP4K2A, HLA-DPB1 and RTEL1) to be associated with CRC. These genetic variants were significant expressions quantitative trait loci in colon tissue, influencing the expression of their mapped genes. Significant interactions of PIP4K2A and HLA-DPB1 with environmental factors, including smoking and alcohol drinking, were observed. All mapped genes and their co-expressed genes were significantly enriched in pathways involved in carcinogenesis.

CONCLUSION

Our findings provide an important insight into the shared genetic basis between CRC and other cancers. We revealed several novel CRC susceptibility loci to help understand the genetic architecture of CRC.

The MC1R r allele does not increase melanoma risk in MITF E318K carriers

BACKGROUND

Population-wide screening for melanoma is not cost-effective, but genetic characterization could facilitate risk stratification and targeted screening. Common Melanocortin-1 receptor (MC1R) red hair colour (RHC) variants and Microphthalmia-associated transcription factor (MITF) E318K separately confer moderate melanoma susceptibility, but their interactive effects are relatively unexplored.

OBJECTIVES

To evaluate whether MC1R genotypes differentially affect melanoma risk in MITF E318K+ vs. E318K- individuals.

MATERIALS AND METHODS

Melanoma status (affected or unaffected) and genotype data (MC1R and MITF E318K) were collated from research cohorts (five Australian and two European). In addition, RHC genotypes from E318K+ individuals with and without melanoma were extracted from databases (The Cancer Genome Atlas and Medical Genome Research Bank, respectively). χ2 and logistic regression were used to evaluate RHC allele and genotype frequencies within E318K+/- cohorts depending on melanoma status. Replication analysis was conducted on 200 000 general-population exomes (UK Biobank).

RESULTS

The cohort comprised 1165 MITF E318K- and 322 E318K+ individuals. In E318K- cases MC1R R and r alleles increased melanoma risk relative to wild type (wt), P < 0.001 for both. Similarly, each MC1R RHC genotype (R/R, R/r, R/wt, r/r and r/wt) increased melanoma risk relative to wt/wt (P < 0.001 for all). In E318K+ cases, R alleles increased melanoma risk relative to the wt allele [odds ratio (OR) 2.04 (95% confidence interval 1.67-2.49); P = 0.01], while the r allele risk was comparable with the wt allele [OR 0.78 (0.54-1.14) vs. 1.00, respectively]. E318K+ cases with the r/r genotype had a lower but not significant melanoma risk relative to wt/wt [OR 0.52 (0.20-1.38)]. Within the E318K+ cohort, R genotypes (R/R, R/r and R/wt) conferred a significantly higher risk compared with non-R genotypes (r/r, r/wt and wt/wt) (P < 0.001). UK Biobank data supported our findings that r did not increase melanoma risk in E318K+ individuals.

CONCLUSIONS

RHC alleles/genotypes modify melanoma risk differently in MITF E318K- and E318K+ individuals. Specifically, although all RHC alleles increase risk relative to wt in E318K- individuals, only MC1R R increases melanoma risk in E318K+ individuals. Importantly, in the E318K+ cohort the MC1R r allele risk is comparable with wt. These findings could inform counselling and management for MITF E318K+ individuals.

Causal relationships between risk of venous thromboembolism and 18 cancers: a bidirectional Mendelian randomisation analysis

Background

People with cancer experience high rates of venous thromboembolism (VTE). Additionally, risk of subsequent cancer is increased in people experiencing their first VTE. The causal mechanisms underlying this association are not completely understood, and it is unknown whether VTE is itself a risk factor for cancer.

Methods

We used data from large genome-wide association study meta-analyses to perform bi-directional Mendelian randomisation analyses to estimate causal associations between genetically-proxied lifetime risk of VTE and risk of 18 different cancers.

Results

We found no conclusive evidence that genetically-proxied lifetime risk of VTE was causally associated with an increased incidence of cancer, or vice-versa. We observed an association between VTE and pancreatic cancer risk (odds ratio for pancreatic cancer 1.23 (95% confidence interval 1.08 - 1.40) per log-odds increase in risk of VTE, P = 0.002). However, sensitivity analyses revealed this association was predominantly driven by a variant proxying non-O blood group, with inadequate evidence from Mendelian randomisation to suggest a causal relationship.

Conclusions

These findings do not support the hypothesis that genetically-proxied lifetime risk of VTE is a cause of cancer. Existing observational epidemiological associations between VTE and cancer are therefore more likely to be driven by pathophysiological changes which occur in the setting of active cancer and anti-cancer treatments. Further work is required to explore and synthesise evidence for these mechanisms.

The journey from melanocytes to melanoma

Over the past decade, melanoma has led the field in new cancer treatments, with impressive gains in on-treatment survival but more modest improvements in overall survival. Melanoma presents heterogeneity and transcriptional plasticity that recapitulates distinct melanocyte developmental states and phenotypes, allowing it to adapt to and eventually escape even the most advanced treatments. Despite remarkable advances in our understanding of melanoma biology and genetics, the melanoma cell of origin is still fiercely debated because both melanocyte stem cells and mature melanocytes can be transformed. Animal models and high-throughput single-cell sequencing approaches have opened new opportunities to address this question. Here, we discuss the melanocytic journey from the neural crest, where they emerge as melanoblasts, to the fully mature pigmented melanocytes resident in several tissues. We describe a new understanding of melanocyte biology and the different melanocyte subpopulations and microenvironments they inhabit, and how this provides unique insights into melanoma initiation and progression. We highlight recent findings on melanoma heterogeneity and transcriptional plasticity and their implications for exciting new research areas and treatment opportunities. The lessons from melanocyte biology reveal how cells that are present to protect us from the damaging effects of ultraviolet radiation reach back to their origins to become a potentially deadly cancer.

Targeted long-read sequencing of the Ewing sarcoma 6p25.1 susceptibility locus identifies germline-somatic interactions with EWSR1-FLI1 binding

Ewing sarcoma (EwS) is a rare bone and soft tissue malignancy driven by chromosomal translocations encoding chimeric transcription factors, such as EWSR1-FLI1, that bind GGAA motifs forming novel enhancers that alter nearby expression. We propose that germline microsatellite variation at the 6p25.1 EwS susceptibility locus could impact downstream gene expression and EwS biology. We performed targeted long-read sequencing of EwS blood DNA to characterize variation and genomic features important for EWSR1-FLI1 binding. We identified 50 microsatellite alleles at 6p25.1 and observed that EwS-affected individuals had longer alleles (>135 bp) with more GGAA repeats. The 6p25.1 GGAA microsatellite showed chromatin features of an EWSR1-FLI1 enhancer and regulated expression of RREB1, a transcription factor associated with RAS/MAPK signaling. RREB1 knockdown reduced proliferation and clonogenic potential and reduced expression of cell cycle and DNA replication genes. Our integrative analysis at 6p25.1 details increased binding of longer GGAA microsatellite alleles with acquired EWSR-FLI1 to promote Ewing sarcomagenesis by RREB1-mediated proliferation.

Joint analysis of multiple phenotypes for extremely unbalanced case-control association studies

In genome-wide association studies (GWAS) for thousands of phenotypes in biobanks, most binary phenotypes have substantially fewer cases than controls. Many widely used approaches for joint analysis of multiple phenotypes produce inflated type I error rates for such extremely unbalanced case-control phenotypes. In this research, we develop a method to jointly analyze multiple unbalanced case-control phenotypes to circumvent this issue. We first group multiple phenotypes into different clusters based on a hierarchical clustering method, then we merge phenotypes in each cluster into a single phenotype. In each cluster, we use the saddlepoint approximation to estimate the p value of an association test between the merged phenotype and a single nucleotide polymorphism (SNP) which eliminates the issue of inflated type I error rate of the test for extremely unbalanced case-control phenotypes. Finally, we use the Cauchy combination method to obtain an integrated p value for all clusters to test the association between multiple phenotypes and a SNP. We use extensive simulation studies to evaluate the performance of the proposed approach. The results show that the proposed approach can control type I error rate very well and is more powerful than other available methods. We also apply the proposed approach to phenotypes in category IX (diseases of the circulatory system) in the UK Biobank. We find that the proposed approach can identify more significant SNPs than the other viable methods we compared with.