Predicting the impact of rare variants on RNA splicing in CAGI6.

Variants which disrupt splicing are a frequent cause of rare disease that have been under-ascertained clinically. Accurate and efficient methods to predict a variant's impact on splicing are needed to interpret the growing number of variants of unknown significance (VUS) identified by exome and genome sequencing. Here, we present the results of the CAGI6 Splicing VUS challenge, which invited predictions of the splicing impact of 56 variants ascertained clinically and functionally validated to determine splicing impact. The performance of 12 prediction methods, along with SpliceAI and CADD, was compared on the 56 functionally validated variants. The maximum accuracy achieved was 82% from two different approaches, one weighting SpliceAI scores by minor allele frequency, and one applying the recently published Splicing Prediction Pipeline (SPiP). SPiP performed optimally in terms of sensitivity, while an ensemble method combining multiple prediction tools and information from databases exceeded all others for specificity. Several challenge methods equalled or exceeded the performance of SpliceAI, with ultimate choice of prediction method likely to depend on experimental or clinical aims. One quarter of the variants were incorrectly predicted by at least 50% of the methods, highlighting the need for further improvements to splicing prediction methods for successful clinical application.

Novel Alu insertion in the ZEB2 gene causing Mowat-Wilson syndrome.

Alu elements are short, interspersed elements located throughout the genome, playing a role in human diversity, and occasionally causing genetic diseases. Here, we report a novel Alu insertion causing Mowat-Wilson syndrome, a rare neurodevelopmental disorder, in an 8-year-old boy displaying the typical clinical features for Mowat-Wilson syndrome. The variant was not initially detected in genome sequencing data, but through deep phenotyping, which pointed to only one plausible candidate gene, manual inspection of genome sequencing alignment data enabled us to identify a de novo heterozygous Alu insertion in exon 8 of the ZEB2 gene. Nanopore long-read sequencing confirmed the Alu insertion, leading to the formation of a premature stop codon and likely haploinsufficiency of ZEB2. This underscores the importance of deep phenotyping and mobile element insertion analysis in uncovering genetic causes of monogenic disorders as these elements might be overlooked in standard next-generation sequencing protocols.

Aortic dissection in a young male with persistent ductus arteriosus and a novel variant in MYLK.

Pathogenic variants in several genes involved in the function or regulation of smooth muscle cells (SMC) are known to predispose to congenital heart disease and thoracic aortic aneurysm and dissection (TAAD). Variants in MYLK are primarily known to predispose to TAAD, but a growing body of evidence points toward MYLK also playing an essential role in the regulation of SMC contraction outside the aorta. In this case report, we present a patient with co-occurrence of persistent ductus arteriosus (PDA) and thoracic aortic dissection. Genetic analyses revealed a novel splice acceptor variant (c.3986-1G > A) in MYLK, which segregated with disease in the family. RNA-analyses on fibroblasts showed that the variant induced skipping of exon 24, which resulted in an in-frame deletion of 101 amino acids. These findings suggest that MYLK-associated disease could include a broader phenotypic spectrum than isolated TAAD, including PDA and obstructive pulmonary disease. Genetic analyses could be considered in families with TAAD and PDA or obstructive pulmonary disease.

Germline (epi)genetics reveals high predisposition in females: a 5-year, nationwide, prospective Wilms tumour cohort.

BACKGROUND: Studies suggest that Wilms tumours (WT) are caused by underlying genetic (5%-10%) and epigenetic (2%-29%) mechanisms, yet studies covering both aspects are sparse. METHODS: We performed prospective whole-genome sequencing of germline DNA in Danish children diagnosed with WT from 2016 to 2021, and linked genotypes to deep phenotypes. RESULTS: Of 24 patients (58% female), 3 (13%, all female) harboured pathogenic germline variants in WT risk genes (FBXW7, WT1 and REST). Only one patient had a family history of WT (3 cases), segregating with the REST variant. Epigenetic testing revealed one (4%) additional patient (female) with uniparental disomy of chromosome 11 and Beckwith-Wiedemann syndrome (BWS). We observed a tendency of higher methylation of the BWS-related imprinting centre 1 in patients with WT than in healthy controls. Three patients (13%, all female) with bilateral tumours and/or features of BWS had higher birth weights (4780 g vs 3575 g; p=0.002). We observed more patients with macrosomia (>4250 g, n=5, all female) than expected (OR 9.98 (95% CI 2.56 to 34.66)). Genes involved in early kidney development were enriched in our constrained gene analysis, including both known (WT1, FBXW7) and candidate (CTNND1, FRMD4A) WT predisposition genes. WT predisposing variants, BWS and/or macrosomia (n=8, all female) were more common in female patients than male patients (p=0.01). CONCLUSION: We find that most females (57%) and 33% of all patients with WT had either a genetic or another indicator of WT predisposition. This emphasises the need for scrutiny when diagnosing patients with WT, as early detection of underlying predisposition may impact treatment, follow-up and genetic counselling.

Update of penetrance estimates in Birt-Hogg-Dubé syndrome.

BACKGROUND: Birt-Hogg-Dubé (BHD) syndrome is a rare genetic syndrome caused by pathogenic or likely pathogenic germline variants in the FLCN gene. Patients with BHD syndrome have an increased risk of fibrofolliculomas, pulmonary cysts, pneumothorax and renal cell carcinoma. There is debate regarding whether colonic polyps should be added to the criteria. Previous risk estimates have mostly been based on small clinical case series. METHODS: A comprehensive review was conducted to identify studies that had recruited families carrying pathogenic or likely pathogenic variants in FLCN. Pedigree data were requested from these studies and pooled. Segregation analysis was used to estimate the cumulative risk of each manifestation for carriers of FLCN pathogenic variants. RESULTS: Our final dataset contained 204 families that were informative for at least one manifestation of BHD (67 families informative for skin manifestations, 63 for lung, 88 for renal carcinoma and 29 for polyps). By age 70 years, male carriers of the FLCN variant have an estimated 19% (95% CI 12% to 31%) risk of renal tumours, 87% (95% CI 80% to 92%) of lung involvement and 87% (95% CI 78% to 93%) of skin lesions, while female carriers had an estimated 21% (95% CI 13% to 32%) risk of renal tumours, 82% (95% CI 73% to 88%) of lung involvement and 78% (95% CI 67% to 85%) of skin lesions. The cumulative risk of colonic polyps by age 70 years old was 21% (95% CI 8% to 45%) for male carriers and 32% (95% CI 16% to 53%) for female carriers. CONCLUSIONS: These updated penetrance estimates, based on a large number of families, are important for the genetic counselling and clinical management of BHD syndrome.

Clinical implications of genetic testing in familial intermediate and late-onset colorectal cancer.

The genetic background of familial, late-onset colorectal cancer (CRC) (i.e., onset > age 50 years) has not been studied as thoroughly as other subgroups of familial CRC, and the proportion of families with a germline genetic predisposition to CRC remains to be defined. To define the contribution of known or suggested CRC predisposition genes to familial late-onset CRC, we analyzed 32 well-established or candidate CRC predisposition genes in 75 families with late-onset CRC. We identified pathogenic or likely pathogenic variants in five patients in MSH6 (n = 1), MUTYH (monoallelic; n = 2) and NTHL1 (monoallelic; n = 2). In addition, we identified a number of variants of unknown significance in particular in the lower penetrant Lynch syndrome-associated mismatch repair (MMR) gene MSH6 (n = 6). In conclusion, screening using a comprehensive cancer gene panel in families with accumulation of late-onset CRC appears not to have a significant clinical value due to the low level of high-risk pathogenic variants detected. Our data suggest that only patients with abnormal MMR immunohistochemistry (IHC) or microsatellite instability (MSI) analyses, suggestive of Lynch syndrome, or a family history indicating another cancer predisposition syndrome should be prioritized for such genetic evaluations. Variants in MSH6 and MUTYH have previously been proposed to be involved in digenic or oligogenic hereditary predisposition to CRC. Accumulation of variants in MSH6 and monoallelic, pathogenic variants in MUTYH in our study indicates that digenic or oligogenic inheritance might be involved in late-onset CRC and warrants further studies of complex types of inheritance.

Breast cancer survival in Nordic BRCA2 mutation carriers-unconventional association with oestrogen receptor status.

BACKGROUND: The natural history of breast cancer among BRCA2 carriers has not been clearly established. In a previous study from Iceland, positive ER status was a negative prognostic factor. We sought to identify factors that predicted survival after invasive breast cancer in an expanded cohort of BRCA2 carriers. METHODS: We studied 608 women with invasive breast cancer and a pathogenic BRCA2 mutation (variant) from four Nordic countries. Information on prognostic factors and treatment was retrieved from health records and by analysis of archived tissue specimens. Hazard ratios (HR) were estimated for breast cancer-specific survival using Cox regression. RESULTS: About 77% of cancers were ER-positive, with the highest proportion (83%) in patients under 40 years. ER-positive breast cancers were more likely to be node-positive (59%) than ER-negative cancers (34%) (P < 0.001). The survival analysis included 584 patients. Positive ER status was protective in the first 5 years from diagnosis (multivariate HR = 0.49; 95% CI 0.26-0.93, P = 0.03); thereafter, the effect was adverse (HR = 1.91; 95% CI 1.07-3.39, P = 0.03). The adverse effect of positive ER status was limited to women who did not undergo endocrine treatment (HR = 2.36; 95% CI 1.26-4.44, P = 0.01) and patients with intact ovaries (HR = 1.99; 95% CI 1.11-3.59, P = 0.02). CONCLUSIONS: The adverse effect of a positive ER status in BRCA2 carriers with breast cancer may be contingent on exposure to ovarian hormones.

The BRCA1 c. 5096G>A p.Arg1699Gln (R1699Q) intermediate risk variant: breast and ovarian cancer risk estimation and recommendations for clinical management from the ENIGMA consortium.

BACKGROUND: We previously showed that the BRCA1 variant c.5096G>A p.Arg1699Gln (R1699Q) was associated with an intermediate risk of breast cancer (BC) and ovarian cancer (OC). This study aimed to assess these cancer risks for R1699Q carriers in a larger cohort, including follow-up of previously studied families, to further define cancer risks and to propose adjusted clinical management of female BRCA1*R1699Q carriers. METHODS: Data were collected from 129 BRCA1*R1699Q families ascertained internationally by ENIGMA (Evidence-based Network for the Interpretation of Germline Mutant Alleles) consortium members. A modified segregation analysis was used to calculate BC and OC risks. Relative risks were calculated under both monogenic model and major gene plus polygenic model assumptions. RESULTS: In this cohort the cumulative risk of BC and OC by age 70 years was 20% and 6%, respectively. The relative risk for developing cancer was higher when using a model that included the effects of both the R1699Q variant and a residual polygenic component compared with monogenic model (for BC 3.67 vs 2.83, and for OC 6.41 vs 5.83). CONCLUSION: Our results confirm that BRCA1*R1699Q confers an intermediate risk for BC and OC. Breast surveillance for female carriers based on mammogram annually from age 40 is advised. Bilateral salpingo-oophorectomy should be considered based on family history.

Human papillomavirus-related carcinoma with adenoid cystic-like features of the sinonasal tract: clinical and morphological characterization of six new cases.

AIMS: Human papillomavirus (HPV) is known as causative for squamous cell carcinoma (SCC) of the oropharynx, but is also found not infrequently in carcinomas of the sinonasal tract. Recently, a subset of these carcinomas was recognized to harbour HPV33 and have a significant morphological overlap with adenoid cystic carcinoma (ACC), a rare and aggressive carcinoma originating in the minor salivary glands. Termed 'HPV-related carcinoma with ACC-like features', only nine cases have been reported. To clarify the occurrence of these tumours we screened a large material for the presence of HPV-related ACC-like carcinoma. The identified tumours were characterized immunohistochemically and with fluorescence in-situ hybridization, and clinicopathological information for all cases is presented. METHODS AND RESULTS: Forty-seven candidate cases were screened for presence of HPV. Six cases were identified and genotyped as HPV types 33, 35, and 56. All six cases had areas of dysplastic mucosal lining and showed remarkable heterogeneous morphologies. MYB, MYBL1, and NFIB genes were intact and, interestingly, staining for MYB protein was largely negative in contrast to what was found in ACC. One patient experienced a local recurrence 11 years after initial treatment and the remaining five patients were alive without evidence of disease. CONCLUSION: We report six new cases of HPV-related ACC-like carcinoma and found that, although in a small material, the prognosis for these patients seems more favourable than for ACC. For the distinction between ACC and HPV-related ACC-like carcinoma, p16, MYB immunohistochemistry or investigation of MYB, MYBL1 and NFIB gene status are valuable.

Comparison of mRNA splicing assay protocols across multiple laboratories: recommendations for best practice in standardized clinical testing.

BACKGROUND: Accurate evaluation of unclassified sequence variants in cancer predisposition genes is essential for clinical management and depends on a multifactorial analysis of clinical, genetic, pathologic, and bioinformatic variables and assays of transcript length and abundance. The integrity of assay data in turn relies on appropriate assay design, interpretation, and reporting. METHODS: We conducted a multicenter investigation to compare mRNA splicing assay protocols used by members of the ENIGMA (Evidence-Based Network for the Interpretation of Germline Mutant Alleles) consortium. We compared similarities and differences in results derived from analysis of a panel of breast cancer 1, early onset (BRCA1) and breast cancer 2, early onset (BRCA2) gene variants known to alter splicing (BRCA1: c.135-1G>T, c.591C>T, c.594-2A>C, c.671-2A>G, and c.5467+5G>C and BRCA2: c.426-12_8delGTTTT, c.7988A>T, c.8632+1G>A, and c.9501+3A>T). Differences in protocols were then assessed to determine which elements were critical in reliable assay design. RESULTS: PCR primer design strategies, PCR conditions, and product detection methods, combined with a prior knowledge of expected alternative transcripts, were the key factors for accurate splicing assay results. For example, because of the position of primers and PCR extension times, several isoforms associated with BRCA1, c.594-2A>C and c.671-2A>G, were not detected by many sites. Variation was most evident for the detection of low-abundance transcripts (e.g., BRCA2 c.8632+1G>A Δ19,20 and BRCA1 c.135-1G>T Δ5q and Δ3). Detection of low-abundance transcripts was sometimes addressed by using more analytically sensitive detection methods (e.g., BRCA2 c.426-12_8delGTTTT ins18bp). CONCLUSIONS: We provide recommendations for best practice and raise key issues to consider when designing mRNA assays for evaluation of unclassified sequence variants.

The evolutionary impact of childhood cancer on the human gene pool.

Germline pathogenic variants associated with increased childhood mortality must be subject to natural selection. Here, we analyze publicly available germline genetic metadata from 4,574 children with cancer [11 studies; 1,083 whole exome sequences (WES), 1,950 whole genome sequences (WGS), and 1,541 gene panel] and 141,456 adults [125,748 WES and 15,708 WGS]. We find that pediatric cancer predisposition syndrome (pCPS) genes [n = 85] are highly constrained, harboring only a quarter of the loss-of-function variants that would be expected. This strong indication of selective pressure on pCPS genes is found across multiple lines of germline genomics data from both pediatric and adult cohorts. For six genes [ELP1, GPR161, VHL and SDHA/B/C], a clear lack of mutational constraint calls the pediatric penetrance and/or severity of associated cancers into question. Conversely, out of 23 known pCPS genes associated with biallelic risk, two [9%, DIS3L2 and MSH2] show significant constraint, indicating that they may monoallelically increase childhood cancer risk. In summary, we show that population genetic data provide empirical evidence that heritable childhood cancer leads to natural selection powerful enough to have significantly impacted the present-day gene pool.

Actionability and familial uptake following opportunistic genomic screening in a pediatric cancer cohort.

The care for patients with serious conditions is increasingly guided by genomic medicine, and genomic medicine may equally transform care for healthy individual if genomic population screening is implemented. This study examines the medical impact of opportunistic genomic screening (OGS) in a cohort of patients undergoing comprehensive genomic germline DNA testing for childhood cancer, including the impact on their relatives. Medical actionability and uptake after cascade testing in the period following disclosure of OGS results was quantified. A secondary finding was reported to 19/595 (3.2%) probands primarily in genes related to cardiovascular and lipid disorders. After a mean follow up time of 1.6 years (Interquartile range (IQR): 0.57-1.92 yrs.) only 12 (63%) of these variants were found to be medically actionable. Clinical follow up or treatment was planned in 16 relatives, and as in the probands, the prescribed treatment was primarily betablockers or cholesterol lowering therapy. No invasive procedures or implantation of medical devices were performed in probands or relatives, and no reproductive counseling was requested. After an average of 1.6 years of follow-up 2.25 relatives per family with an actionable finding had been tested. This real-world experience of OGS grants new insight into the practical implementation effects and derived health care demands of genotype-first screening. The resulting health care effect and impact on demand for genetic counseling and workup in relatives extends beyond the effect in the probands.

BRCA1 R1699Q variant displaying ambiguous functional abrogation confers intermediate breast and ovarian cancer risk.

BACKGROUND: Clinical classification of rare sequence changes identified in the breast cancer susceptibility genes BRCA1 and BRCA2 is essential for appropriate genetic counselling of individuals carrying these variants. We previously showed that variant BRCA1 c.5096G>A p.Arg1699Gln in the BRCA1 transcriptional transactivation domain demonstrated equivocal results from a series of functional assays, and proposed that this variant may confer low to moderate risk of cancer. METHODS: Measures of genetic risk (report of family history, segregation) were assessed for 68 BRCA1 c.5096G>A p.Arg1699Gln (R1699Q) families recruited through family cancer clinics, comparing results with 34 families carrying the previously classified pathogenic BRCA1 c.5095C>T p.Arg1699Trp (R1699W) mutation at the same residue, and to 243 breast cancer families with no BRCA1 pathogenic mutation (BRCA-X). RESULTS: Comparison of BRCA1 carrier prediction scores of probands using the BOADICEA risk prediction tool revealed that BRCA1 c.5096G>A p.Arg1699Gln variant carriers had family histories that were less 'BRCA1-like' than BRCA1 c.5095C>T p.Arg1699Trp mutation carriers (p<0.00001), but more 'BRCA1-like' than BRCA-X families (p=0.0004). Further, modified segregation analysis of the subset of 30 families with additional genotyping showed that BRCA1 c.5096G >A p.Arg1699Gln had reduced penetrance compared with the average truncating BRCA1 mutation penetrance (p=0.0002), with estimated cumulative risks to age 70 of breast or ovarian cancer of 24%. CONCLUSIONS: Our results provide substantial evidence that the BRCA1 c.5096G>A p.Arg1699Gln (R1699Q) variant, demonstrating ambiguous functional deficiency across multiple assays, is associated with intermediate risk of breast and ovarian cancer, highlighting challenges for risk modelling and clinical management of patients of this and other potential moderate-risk variants.

Homozygous splice variant (c.1741-6G>A) of the COL6A1 gene in three patients with Ullrich congenital muscular dystrophy.

The three major collagen VI genes: COL6A1, COL6A2, and COL6A3 encode microfibrillar components of extracellular matrices in multiple tissues including muscles and tendons. Pathogenic variants in the collagen VI genes cause collagen VI-related dystrophies representing a continuum of conditions from Bethlem myopathy at the milder end to Ullrich congenital muscular dystrophy at the more severe end. Here we describe a pathogenic variant in the COL6A1 gene (NM_001848.3; c.1741-6G>A) found in homozygosity in three patients with Ullrich congenital muscular dystrophy. The patients suffered from severe muscle impairment characterised by proximal weakness, distal hyperlaxity, joint contractures, wheelchair-dependency, and use of nocturnal non-invasive ventilation. The pathogenicity was verified by RNA analyses showing that the variant induced aberrant splicing leading to a frameshift and loss of function. The analyses were in line with immunocytochemistry studies of patient-derived skin fibroblasts and muscle tissue demonstrating impaired secretion of collagen VI into the extracellular matrix. Thereby, we add the variant c.1741-6G>A to the list of pathogenic, recessive, splice variants in COL6A1 causing Ullrich congenital muscular dystrophy. The variant is listed in ClinVar as of "uncertain significance" and "likely benign" and may presumably have been overlooked in other patients.

The effect of a single SMARCA4 exon deletion on RNA splicing: Implications for variant classification.

BACKGROUND: Exon deletions are generally considered pathogenic, particularly when they are located out of frame. Here, we describe a pediatric, female patient presenting with hypercalcemia and a small cell carcinoma of the ovary, hypercalcemic type, and carrying a germline de novo SMARCA4 exon 14 deletion. METHODS: The SMARCA4 deletion was identified by whole genome sequencing, and the effect on the RNA level was examined by gel- and capillary electrophoresis and nanopore sequencing. RESULTS: The deletion was in silico predicted to be truncating, but RNA analysis revealed two major transcripts with deletion of exon 14 alone or exon 14 through 15, where the latter was located in-frame. Because the patient's phenotype matched that of other patients with pathogenic germline variants in SMARCA4, the deletion was classified as likely pathogenic. CONCLUSION: We propose to include RNA analysis in classification of single-exon deletions, especially if located outside of known functional domains, as this can identify any disparate effects on the RNA and DNA level, which may have implications for variant classification using the American College of Medical Genetics and Genomics guidelines.

ENIGMA CHEK2gether Project: A Comprehensive Study Identifies Functionally Impaired CHEK2 Germline Missense Variants Associated with Increased Breast Cancer Risk.

PURPOSE: Germline pathogenic variants in CHEK2 confer moderately elevated breast cancer risk (odds ratio, OR ∼ 2.5), qualifying carriers for enhanced breast cancer screening. Besides pathogenic variants, dozens of missense CHEK2 variants of uncertain significance (VUS) have been identified, hampering the clinical utility of germline genetic testing (GGT). EXPERIMENTAL DESIGN: We collected 460 CHEK2 missense VUS identified by the ENIGMA consortium in 15 countries. Their functional characterization was performed using CHEK2-complementation assays quantifying KAP1 phosphorylation and CHK2 autophosphorylation in human RPE1-CHEK2-knockout cells. Concordant results in both functional assays were used to categorize CHEK2 VUS from 12 ENIGMA case-control datasets, including 73,048 female patients with breast cancer and 88,658 ethnicity-matched controls. RESULTS: A total of 430/460 VUS were successfully analyzed, of which 340 (79.1%) were concordant in both functional assays and categorized as functionally impaired (N = 102), functionally intermediate (N = 12), or functionally wild-type (WT)-like (N = 226). We then examined their association with breast cancer risk in the case-control analysis. The OR and 95% CI (confidence intervals) for carriers of functionally impaired, intermediate, and WT-like variants were 2.83 (95% CI, 2.35-3.41), 1.57 (95% CI, 1.41-1.75), and 1.19 (95% CI, 1.08-1.31), respectively. The meta-analysis of population-specific datasets showed similar results. CONCLUSIONS: We determined the functional consequences for the majority of CHEK2 missense VUS found in patients with breast cancer (3,660/4,436; 82.5%). Carriers of functionally impaired missense variants accounted for 0.5% of patients with breast cancer and were associated with a moderate risk similar to that of truncating CHEK2 variants. In contrast, 2.2% of all patients with breast cancer carried functionally wild-type/intermediate missense variants with no clinically relevant breast cancer risk in heterozygous carriers.

Functional characterization of MLH1 missense variants identified in Lynch syndrome patients.

Germline mutations in the human DNA mismatch repair (MMR) genes MSH2 and MLH1 are associated with the inherited cancer disorder Lynch syndrome (LS), also known as hereditary nonpolyposis colorectal cancer or HNPCC. A proportion of MSH2 and MLH1 mutations found in suspected LS patients give rise to single amino acid substitutions. The functional consequences in regard to pathogenicity of many of these variants are unclear. We have examined the functionality of a panel of MLH1 missense mutations found in LS families, by testing the variant proteins in functional assays, addressing subcellular localization, and protein-protein interaction with the dimer partner PMS2 and the MMR-associated exonuclease 1. We show that a significant proportion of examined variant proteins have functional defects in either subcellular localization or protein-protein interactions, which is suspected to lead to the cancer phenotype observed in patients. Moreover, the obtained results correlate well with reported MMR activity and with in silico analysis for a majority of the variants.

9p21.3 Microdeletion involving CDKN2A/2B in a young patient with multiple primary cancers and review of the literature.

Germline pathogenic variants in CDKN2A predispose to various cancers, including melanoma, pancreatic cancer, and neural system tumors, whereas CDKN2B variants are associated with renal cell carcinoma. A few case reports have described heterozygous germline deletions spanning both CDKN2A and CDKN2B associated with a cancer predisposition syndrome (CPS) that constitutes a risk of cancer beyond those associated with haploinsufficiency of each gene individually, indicating an additive effect or a contiguous gene deletion syndrome. We report a young woman with a de novo germline 9p21 microdeletion involving the CDKN2A/CDKN2B genes, who developed six primary cancers since childhood, including a very rare extraskeletal osteosarcoma (eOS) at the age of 8. To our knowledge this is the first report of eOS in a patient with CDKN2A/CDKN2B deletion.