The CDKN2A/p16(INK) (4a) 5'UTR sequence and translational regulation: impact of novel variants predisposing to melanoma

TranSNPs: A class of functional SNPs affecting mRNA translation potential revealed by fraction-based allelic imbalance

Few studies have explored the association between SNPs and alterations in mRNA translation potential. We developed an approach to identify SNPs that can mark allele-specific protein expression levels and could represent sources of inter-individual variation in disease risk. Using MCF7 cells under different treatments, we performed polysomal profiling followed by RNA sequencing of total or polysome-associated mRNA fractions and designed a computational approach to identify SNPs showing a significant change in the allelic balance between total and polysomal mRNA fractions. We identified 147 SNPs, 39 of which located in UTRs. Allele-specific differences at the translation level were confirmed in transfected MCF7 cells by reporter assays. Exploiting breast cancer data from TCGA we identified UTR SNPs demonstrating distinct prognosis features and altering binding sites of RNA-binding proteins. Our approach produced a catalog of tranSNPs, a class of functional SNPs associated with allele-specific translation and potentially endowed with prognostic value for disease risk.

Current Practice in Bicistronic IRES Reporter Use: A Systematic Review

Bicistronic reporter assays have been instrumental for transgene expression, understanding of internal ribosomal entry site (IRES) translation, and identification of novel cap-independent translational elements (CITE). We observed a large methodological variability in the use of bicistronic reporter assays and data presentation or normalization procedures. Therefore, we systematically searched the literature for bicistronic IRES reporter studies and analyzed methodological details, data visualization, and normalization procedures. Two hundred fifty-seven publications were identified using our search strategy (published 1994-2020). Experimental studies on eukaryotic adherent cell systems and the cell-free translation assay were included for further analysis. We evaluated the following methodological details for 176 full text articles: the bicistronic reporter design, the cell line or type, transfection methods, and time point of analyses post-transfection. For the cell-free translation assay, we focused on methods of in vitro transcription, type of translation lysate, and incubation times and assay temperature. Data can be presented in multiple ways: raw data from individual cistrons, a ratio of the two, or fold changes thereof. In addition, many different control experiments have been suggested when studying IRES-mediated translation. In addition, many different normalization and control experiments have been suggested when studying IRES-mediated translation. Therefore, we also categorized and summarized their use. Our unbiased analyses provide a representative overview of bicistronic IRES reporter use. We identified parameters that were reported inconsistently or incompletely, which could hamper data reproduction and interpretation. On the basis of our analyses, we encourage adhering to a number of practices that should improve transparency of bicistronic reporter data presentation and improve methodological descriptions to facilitate data replication.

Variants in the 5'UTR reduce SHOX expression and contribute to SHOX haploinsufficiency

SHOX haploinsufficiency causes 70-90% of Léri-Weill dyschondrosteosis (LWD) and 2-10% of idiopathic short stature (ISS). Deletions removing the entire gene or enhancers and point mutations in the coding region represent a well-established cause of haploinsufficiency. During diagnostic genetic testing on ISS/LWD patients, in addition to classic SHOX defects, five 5'UTR variants (c.-58G > T, c.-55C > T, c.-51G > A, c.-19G > A, and c.-9del), were detected whose pathogenetic role was unclear and were thus classified as VUS (Variants of Uncertain Significance). The purpose of the present study was to investigate the role of these noncoding variations in SHOX haploinsufficiency. The variants were tested for their ability to interfere with correct gene expression of a regulated reporter gene (luciferase assay). The negative effect on the mRNA splicing predicted in silico for c.-19G > A was assayed in vitro through a minigene splicing assay. The luciferase assay showed that c.-51G > A, c.-19G > A, and c.-9del significantly reduce luciferase activity by 60, 35, and 40% at the homozygous state. Quantification of the luciferase mRNA showed that c.-51G > A and c.-9del might interfere with the correct SHOX expression mainly at the post-transcriptional level. The exon trapping assay demonstrated that c.-19G > A determines the creation of a new branch site causing an aberrant mRNA splicing. In conclusion, this study allowed us to reclassify two of the 5'UTR variants identified during SHOX diagnostic screening as likely pathogenic, one remains as a VUS, and two as likely benign variants. This analysis for the first time expands the spectrum of the genetic causes of SHOX haploinsufficiency to noncoding variations in the 5'UTR.

Insights into Genetic Susceptibility to Melanoma by Gene Panel Testing: Potential Pathogenic Variants in ACD, ATM, BAP1, and POT1

The contribution of recently established or candidate susceptibility genes to melanoma missing heritability has yet to be determined. Multigene panel testing could increase diagnostic yield and better define the role of candidate genes. We characterized 273 CDKN2A/ARF and CDK4-negative probands through a custom-designed targeted gene panel that included CDKN2A/ARF, CDK4, ACD, BAP1, MITF, POT1, TERF2IP, ATM, and PALB2. Co-segregation, loss of heterozygosity (LOH)/protein expression analysis, and splicing characterization were performed to improve variant classification. We identified 16 (5.9%) pathogenic and likely pathogenic variants in established high/medium penetrance cutaneous melanoma susceptibility genes (BAP1, POT1, ACD, MITF, and TERF2IP), including two novel variants in BAP1 and 4 in POT1. We also found four deleterious and five likely deleterious variants in ATM (3.3%). Thus, including potentially deleterious variants in ATM increased the diagnostic yield to about 9%. Inclusion of rare variants of uncertain significance would increase the overall detection yield to 14%. At least 10% of melanoma missing heritability may be explained through panel testing in our population. To our knowledge, this is the highest frequency of putative ATM deleterious variants reported in melanoma families, suggesting a possible role in melanoma susceptibility, which needs further investigation.

Combining molecular and immunohistochemical analyses of key drivers in primary melanomas: interplay between germline and somatic variations

Due to the high mutational somatic burden of Cutaneous Malignant Melanoma (CMM) a thorough profiling of the driver mutations and their interplay is necessary to explain the timing of tumorigenesis or for the identification of actionable genetic events. The aim of this study was to establish the mutation rate of some of the key drivers in melanoma tumorigenesis combining molecular analyses and/or immunohistochemistry in 93 primary CMMs from an Italian cohort also characterized for germline status, and to investigate an interplay between germline and somatic variants. BRAF mutations were present in 68% of cases, while CDKN2A germline mutations were found in 16 % and p16 loss in tissue was found in 63%. TERT promoter somatic mutations were detected in 38% of cases while the TERT –245T>C polymorphism was found in 51% of cases. NRAS mutations were found in 39% of BRAF negative or undetermined cases. NF1 was expressed in all cases analysed. MC1R variations were both considered as a dichotomous variable or scored. While a positive, although not significant association between CDKN2A germline mutations, but not MC1R variants, and BRAF somatic mutation was found, we did not observe other associations between germline and somatic events. A yet undescribed inverse correlation between TERT –245T>C polymorphism and the presence of BRAF mutation was found. It is possible to hypothesize that –245T>C polymorphism could be included in those genotypes which may influence the occurrence of BRAF mutations. Further studies are needed to investigate the role of –245T>C polymorphism as a germline predictor of BRAF somatic mutation status.

Functional analysis of a CDKN2A 5'UTR germline variant associated with pancreatic cancer development

CDKN2A coding region germline variants are associated with pancreatic adenocarcinoma (PC) susceptibility. Recently, we described functional germline 5’UTR CDKN2A variants from melanoma patients affecting the post-transcriptional regulation of p16^INK4a mRNA that is dependent, at least in part, on an Internal Ribosome Entry Site (IRES) in the 5’UTR region. Here we describe a 5’UTR c.-201_-198delinsCTTT CDKN2A variant (frequency 0.0028 based on 350 PC patients), which seems to be private to PC, since it has never been found in public databases nor in thousands of melanoma patients tested. Functional analyses confirmed IRES activity of the 5’UTR in BX-PC3 PC cells and revealed a functional impact of the identified variant. Using gene reporter assays we observed reduced translation potential in cells treated with the mTOR inhibitor Torin1, a condition that favors the assessment of IRES activity. At the endogenous gene level we quantified allelic imbalance among polysome-associated mRNAs using a patient-derived cell line heterozygous for the c.-201_-198delinsCTTT. Overall, we conclude that this very rare private variant can be considered a potential mutation, specifically associated with PC. Our data indicate that sequencing of the entire 5'UTR of CDKN2A should be included in routine screening of PC cases with suspected inherited susceptibility.

Performance of in silico tools for the evaluation of p16INK4a (CDKN2A) variants in CAGI

Correct phenotypic interpretation of variants of unknown significance for cancer-associated genes is a diagnostic challenge as genetic screenings gain in popularity in the next-generation sequencing era. The Critical Assessment of Genome Interpretation (CAGI) experiment aims to test and define the state of the art of genotype-phenotype interpretation. Here, we present the assessment of the CAGI p16INK4a challenge. Participants were asked to predict the effect on cellular proliferation of 10 variants for the p16INK4a tumor suppressor, a cyclin-dependent kinase inhibitor encoded by the CDKN2A gene. Twenty-two pathogenicity predictors were assessed with a variety of accuracy measures for reliability in a medical context. Different assessment measures were combined in an overall ranking to provide more robust results. The R scripts used for assessment are publicly available from a GitHub repository for future use in similar assessment exercises. Despite a limited test-set size, our findings show a variety of results, with some methods performing significantly better. Methods combining different strategies frequently outperform simpler approaches. The best predictor, Yang&Zhou lab, uses a machine learning method combining an empirical energy function measuring protein stability with an evolutionary conservation term. The p16INK4a challenge highlights how subtle structural effects can neutralize otherwise deleterious variants.

Cap-Independent Translational Control of Carcinogenesis

Translational regulation has been shown to play an important role in cancer and tumor progression. Despite this fact, the role of translational control in cancer is an understudied and under appreciated field, most likely due to the technological hurdles and paucity of methods available to establish that changes in protein levels are due to translational regulation. Tumors are subjected to many adverse stress conditions such as hypoxia or starvation. Under stress conditions, translation is globally downregulated through several different pathways in order to conserve energy and nutrients. Many of the proteins that are synthesized during stress in order to cope with the stress use a non-canonical or cap-independent mechanism of initiation. Tumor cells have utilized these alternative mechanisms of translation initiation to promote survival during tumor progression. This review will specifically discuss the role of cap-independent translation initiation, which relies on an internal ribosome entry site (IRES) to recruit the ribosomal subunits internally to the messenger RNA. We will provide an overview of the role of IRES-mediated translation in cancer by discussing the types of genes that use IRESs and the conditions under which these mechanisms of initiation are used. We will specifically focus on three well-studied examples: Apaf-1, p53, and c-Jun, where IRES-mediated translation has been demonstrated to play an important role in tumorigenesis or tumor progression.