Patients with CYLD cutaneous syndrome (CCS; syn. Brooke-Spiegler syndrome) carry germline mutations in the tumor suppressor CYLD and develop multiple skin tumors with diverse histophenotypes. Here, we comprehensively profile the genomic landscape of 42 benign and malignant tumors across 13 individuals from four multigenerational families and discover recurrent mutations in epigenetic modifiers DNMT3A and BCOR in 29% of benign tumors. Multi-level and microdissected sampling strikingly reveal that many clones with different DNMT3A mutations exist in these benign tumors, suggesting that intra-tumor heterogeneity is common. Integrated genomic, methylation and transcriptomic profiling in selected tumors suggest that isoform-specific DNMT3A2 mutations are associated with dysregulated methylation. Phylogenetic and mutational signature analyses confirm cylindroma pulmonary metastases from primary skin tumors. These findings contribute to existing paradigms of cutaneous tumorigenesis and metastasis.