Osteolytic lesions (OL) characterize symptomatic multiple myeloma. The mechanisms of how malignant plasma cells (PC) cause OL in one region while others show no signs of bone destruction despite subtotal infiltration remain unknown. We report on a single-cell RNA sequencing (scRNA-seq) study of PC obtained prospectively from random bone marrow aspirates (BM) and paired imaging-guided biopsies of OL. We analyze 148,630 PC from 24 different locations in 10 patients and observe vast inter- and intra-patient heterogeneity based on scRNA-seq analyses. Beyond the limited evidence for spatial heterogeneity from whole-exome sequencing, we find an additional layer of complexity by integrated analysis of anchored scRNA-seq datasets from the BM and OL. PC from OL are characterized by differentially expressed genes compared to PC from BM, including upregulation of genes associated with myeloma bone disease like DKK1, HGF and TIMP-1 as well as recurrent downregulation of JUN/FOS, DUSP1 and HBB. Assessment of PC from longitudinally collected samples reveals transcriptional changes after induction therapy. Our study contributes to the understanding of destructive myeloma bone disease.

In ultraviolet (UV) radiation–exposed skin, mutations fuel clonal cell growth. The relationship between UV exposure and the accumulation of clonal mutations (CMs) and the correlation between CMs and skin cancer risk are largely unexplored. We characterized 450 individual-matched sun-exposed (SE) and non-SE (NE) normal human skin samples. The number and relative contribution of CMs were significantly different between SE and NE areas. Furthermore, we identified hotspots in TP53, NOTCH1, and GRM3 where mutations were significantly associated with UV exposure. In the normal skin from patients with cutaneous squamous cell carcinoma, we found that the cancer burden was associated with the UV-induced mutations, with the difference mostly conferred by the low-frequency CMs. These findings provide previously unknown information on UV’s carcinogenic effect and pave the road for future development of quantitative assessment of subclinical UV damage and skin cancer risk.

The link between carcinogen exposure and cancer immunogenicity is unclear. Single exposure to 12-dimethylbenz[a]anthracene (DMBA) at puberty accelerated spontaneous breast carcinogenesis in mouse mammary tumor virus-polyoma middle tumor-antigen transgenic (MMTV-PyMT^tg or PyMT) and MMTV-Her2/neu^tg (Her2) mice. Paradoxically, DMBA-treated PyMT and Her2 animals were protected from metastasis. CD8⁺ T cells significantly infiltrated DMBA-exposed breast cancers. CD8⁺ T cell depletion resulted in severe lung and liver metastasis in DMBA-treated PyMT mice. Besides increasing tumor mutational burden, DMBA exposure up-regulated Chemokine (C-C motif) ligand 21 (CCL21) in cancer cells and heightened antigen presentation. CCL21 injection suppressed breast cancer growth, and CCL21 receptor deletion attenuated T cell immunity against cancer metastasis in DMBA-treated PyMT animals. CCL21 expression correlated with increased mutational burden and cytolytic activity across human cancers. Higher CCL21 levels correlated with increased CD8⁺ T cell infiltrates in human breast cancer and predicted lower breast cancer distant recurrence rate. Collectively, carcinogen exposure induces immune-activating factors within cancer cells that promote CD8⁺ T cell immunity against metastasis.

Next-Generation Sequencing (NGS) technologies have rapidly advanced our understanding of human variation in cancer. To accurately translate the raw sequencing data into practical knowledge, annotation tools, algorithms and pipelines must be developed that keep pace with the rapidly evolving technology. Currently, a challenge exists in accurately annotating multi-nucleotide variants (MNVs). These tandem substitutions, when affecting multiple nucleotides within a single protein codon of a gene, result in a translated amino acid involving all nucleotides in that codon. Most existing variant callers report a MNV as individual single-nucleotide variants (SNVs), often resulting in multiple triplet codon sequences and incorrect amino acid predictions. To correct potentially misannotated MNVs among reported SNVs, a primary challenge resides in haplotype phasing which is to determine whether the neighboring SNVs are co-located on the same chromosome. Here we describe MAC (Multi-Nucleotide Variant Annotation Corrector), an integrative pipeline developed to correct potentially mis-annotated MNVs. MAC was designed as an application that only requires a SNV file and the matching BAM file as data inputs. Using an example data set containing 3024 SNVs and the corresponding whole-genome sequencing BAM files, we show that MAC identified eight potentially mis-annotated SNVs, and accurately updated the amino acid predictions for seven of the variant calls. MAC can identify and correct amino acid predictions that result from MNVs affecting multiple nucleotides within a single protein codon, which cannot be handled by most existing SNV-based variant pipelines. The MAC software is freely available and represents a useful tool for the accurate translation of genomic sequence to protein function.