Park, Peter J.

Mutations in BRCA1 and BRCA2 genes, whether inherited or somatically acquired, cause homologous recombination deficiency (HRD) in tumor cells. The timing of HRD onset in the emerging tumor lineage is unknown. Here, we present HRDTimer, an algorithm to infer the onset of HRD-driven mutagenesis prior to cancer diagnosis. We estimate that HRD arises at 34% of SBS1-based molecular time---corresponding to a median of 8.3 years (IQR 7.1--10.4) prior to diagnosis in triple-negative breast cancers, and 15.0 years (IQR 12.0--20.6) in ER-positive breast cancers. Bulk sequencing reveals accelerated SBS1 accumulation following neoplastic transformation compared to normal tissue, influencing the estimated age of HRD onset. Single-cell duplex sequencing confirms SBS1 acceleration in tumors and further shows that non-tumor cells largely lack the HRD signature, indicating that HRD is rare in pre-malignant cells, even in BRCA1/2 mutation carriers. Together, our analysis pinpoints the onset of HRD before diagnosis, defining a window for detection and potential interception.

Germline mutations modulate the risk of developing schizophrenia (SCZ). Much less is known about the role of mosaic somatic mutations in the context of SCZ. Deep (239×) whole-genome sequencing (WGS) of brain neurons from 61 SCZ cases and 25 controls postmortem identified mutations occurring during prenatal neurogenesis. SCZ cases showed increased somatic variants in open chromatin, with increased mosaic CpG transversions (CpG>GpG) and T>G mutations at transcription factor binding sites (TFBSs) overlapping open chromatin, a result not seen in controls. Some of these variants alter gene expression, including SCZ risk genes and genes involved in neurodevelopment. Although these mutational processes can reflect a difference in factors indirectly involved in disease, increased somatic mutations at developmental TFBSs could also potentially contribute to SCZ.

Abstract Background Breast cancers in women with a germline BRCA1/2 mutation (gBRCAm) have homologous recombination deficiency (HRD) and are sensitive to therapies causing double-strand DNA breaks. In TBCRC 031 (INFORM), both neoadjuvant cisplatin and doxorubicin/cyclophosphamide (“AC”) in gBRCAm carriers resulted in a complete pathologic response in 18% and 26% respectively in patients with newly diagnosed HER2–negative breast cancer. Herein, we describe molecular features from tumor whole exome (WES) and transcriptome sequencing (RNAseq) associated with response. Methods TBCRC 031 (the INFORM Trial - NCT01670500) was a randomized phase II neoadjuvant trial comparing the efficacy of cisplatin versus AC in gBRCAm carriers with stage I-III HER2–negative breast cancer. Of 118 patients enrolled, 92 patients provided fresh frozen research biopsies, collected prior to chemotherapy initiation, which were subjected to WES and RNAseq. Variants were called using GATK best practices and intratumoral heterogeneity was inferred from mutations and variant allele frequencies using Mutant Allele Tumor Heterogeneity (MATH) scores. Mutational Signature 3 (Sig3) and Genomic Instability Score (GIS) were calculated with SigMA and scarHRD, respectively. RNAseq data were utilized to perform differential gene expression and functional analyses, while cellular deconvolution was performed with CIBERSORTx trained against breast cancer single cell data. Patients were grouped according to their residual cancer burden (RCB) as responders (RCB-0 or 1), or non-responders (RCB-2 or 3). P-values ≤ 0.05 were considered statistically significant and when appropriate, adjustment for multiple testing was performed using a false discovery rate ≤ 5%. Results Of the 92 patients with gBRCAm, 59 (64%) had triple-negative and 33 (36%) were hormone-receptor positive HER2-negative breast cancer, 40 (43%) were classified as responders and 52 (57%) as non-responders. WES analysis revealed that indices of HRD were high (Sig3 exposure predominant and GIS ≥ 42) across most samples irrespective of receptor status and not significantly different among responders and non-responders. In contrast, responders exhibited lower levels of intratumor heterogeneity than non-responders (median MATH 42.9 vs. 33.5, p = 0.01). 223 genes were differentially expressed between responders and non-responders following control for tumor hormone receptor status, BRCA1/2 mutation, and menopausal status. Pathways identified as significantly enriched in upregulated genes were indicative of intrinsic immune activation in responders (e.g., T-cell activation, immune response-signaling, and regulation of leukocyte activation). Cellular deconvolution of the tumor microenvironment confirmed that responders presented a higher proportion of T-cells (p = 0.025) and myeloid cells (p = 0.003) in the tumor samples, while perivascular-like cells were enriched in non-responders (p = 0.034). Conclusion In this analysis of the largest cohort of treatment-naïve gBRCAm breast cancer to date, we found that response to cytotoxic chemotherapy is associated with a transcriptional program of intrinsic immune activation and an increased population of intratumoral T-cells and myeloid cells. Lower levels of intratumor heterogeneity and higher immune activation were associated with response to chemotherapy gBRCAm carriers while HRD scores were not. Citation Format: Felipe Batalini, Doga C. Gulhan, Xanthi Lida Katopodi, Nikolas Kalavros, Antuan Tran, Dimitra Karagkouni, Emily Stern-Gatof, Stuart Schnitt, Judy E. Garber, Gerburg M. Wulf, Peter J. Park, Ioannis Vlachos, Nadine Tung. Intratumor heterogeneity and intrinsic immune activation are associated with response to chemotherapy in BRCA-related breast cancers [abstract]. In: Proceedings of the 2022 San Antonio Breast Cancer Symposium; 2022 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2023;83(5 Suppl):Abstract nr P5-12-03.

We analyzed 131 human brains (44 neurotypical, 19 with Tourette syndrome, 9 with schizophrenia, and 59 with autism) for somatic mutations after whole genome sequencing to a depth of more than 200×. Typically, brains had 20 to 60 detectable single-nucleotide mutations, but ~6% of brains harbored hundreds of somatic mutations. Hypermutability was associated with age and damaging mutations in genes implicated in cancers and, in some brains, reflected in vivo clonal expansions. Somatic duplications, likely arising during development, were found in ~5% of normal and diseased brains, reflecting background mutagenesis. Brains with autism were associated with mutations creating putative transcription factor binding motifs in enhancer-like regions in the developing brain. The top-ranked affected motifs corresponded to MEIS (myeloid ectopic viral integration site) transcription factors, suggesting a potential link between their involvement in gene regulation and autism.

Mutations provide an enduring record Somatic mutations pepper our cells with change, but because they are not in the germline, they do not propagate to the next generation. Bizzotto et al. leveraged data on the distribution of somatic mutations in adults to take a backward look at the earliest moments of human development. Calculation of cellular lineages on the basis of shared somatic mutations shows the number of cells from which the body will develop when the human embryo gastrulates. The lineage for forebrain cells is identifiable, as are the asymmetrical fates spun out of many of the gastrula cells. Science , this issue p. 1249

Abstract Primary central nervous system (CNS) lymphoma (PCNSL) is a diffuse large B-cell lymphoma (DLBCL) confined to the CNS. A genome-wide gene expression comparison between PCNSL and non-CNS DLBCL was performed, the latter consisting of both nodal and extranodal DLBCL (nDLBCL and enDLBCL), to identify a “CNS signature.” Pathway analysis with the program SigPathway revealed that PCNSL is characterized notably by significant differential expression of multiple extracellular matrix (ECM) and adhesion-related pathways. The most significantly up-regulated gene is the ECM-related osteopontin (SPP1). Expression at the protein level of ECM-related SPP1 and CHI3L1 in PCNSL cells was demonstrated by immunohistochemistry. The alterations in gene expression can be interpreted within several biologic contexts with implications for PCNSL, including CNS tropism (ECM and adhesion-related pathways, SPP1, DDR1), B-cell migration (CXCL13, SPP1), activated B-cell subtype (MUM1), lymphoproliferation (SPP1, TCL1A, CHI3L1), aggressive clinical behavior (SPP1, CHI3L1, MUM1), and aggressive metastatic cancer phenotype (SPP1, CHI3L1). The gene expression signature discovered in our study may represent a true “CNS signature” because we contrasted PCNSL with wide-spectrum non-CNS DLBCL on a genomic scale and performed an in-depth bioinformatic analysis.