Variation in histone configurations correlates with gene expression across nine inbred strains of mice.

The Diversity Outbred (DO) mice and their inbred founders are widely used models of human disease. However, although the genetic diversity of these mice has been well documented, their epigenetic diversity has not. Epigenetic modifications, such as histone modifications and DNA methylation, are important regulators of gene expression and, as such, are a critical mechanistic link between genotype and phenotype. Therefore, creating a map of epigenetic modifications in the DO mice and their founders is an important step toward understanding mechanisms of gene regulation and the link to disease in this widely used resource. To this end, we performed a strain survey of epigenetic modifications in hepatocytes of the DO founders. We surveyed four histone modifications (H3K4me1, H3K4me3, H3K27me3, and H3K27ac), as well as DNA methylation. We used ChromHMM to identify 14 chromatin states, each of which represents a distinct combination of the four histone modifications. We found that the epigenetic landscape is highly variable across the DO founders and is associated with variation in gene expression across strains. We found that epigenetic state imputed into a population of DO mice recapitulated the association with gene expression seen in the founders, suggesting that both histone modifications and DNA methylation are highly heritable mechanisms of gene expression regulation. We illustrate how DO gene expression can be aligned with inbred epigenetic states to identify putative cis-regulatory regions. Finally, we provide a data resource that documents strain-specific variation in the chromatin state and DNA methylation in hepatocytes across nine widely used strains of laboratory mice.

A longitudinal circulating tumor DNA-based model associated with survival in metastatic non-small-cell lung cancer.

One of the great challenges in therapeutic oncology is determining who might achieve survival benefits from a particular therapy. Studies on longitudinal circulating tumor DNA (ctDNA) dynamics for the prediction of survival have generally been small or nonrandomized. We assessed ctDNA across 5 time points in 466 non-small-cell lung cancer (NSCLC) patients from the randomized phase 3 IMpower150 study comparing chemotherapy-immune checkpoint inhibitor (chemo-ICI) combinations and used machine learning to jointly model multiple ctDNA metrics to predict overall survival (OS). ctDNA assessments through cycle 3 day 1 of treatment enabled risk stratification of patients with stable disease (hazard ratio (HR) = 3.2 (2.0-5.3), P < 0.001; median 7.1 versus 22.3 months for high- versus low-intermediate risk) and with partial response (HR = 3.3 (1.7-6.4), P < 0.001; median 8.8 versus 28.6 months). The model also identified high-risk patients in an external validation cohort from the randomized phase 3 OAK study of ICI versus chemo in NSCLC (OS HR = 3.73 (1.83-7.60), P = 0.00012). Simulations of clinical trial scenarios employing our ctDNA model suggested that early ctDNA testing outperforms early radiographic imaging for predicting trial outcomes. Overall, measuring ctDNA dynamics during treatment can improve patient risk stratification and may allow early differentiation between competing therapies during clinical trials.

Comprehensive Genomic Profiling (CGP)-Informed Personalized Molecular Residual Disease (MRD) Detection: An Exploratory Analysis from the PREDATOR Study of Metastatic Colorectal Cancer (mCRC) Patients Undergoing Surgical Resection.

A majority of patients with metastatic colorectal cancer (mCRC) experience recurrence post curative-intent surgery. The addition of adjuvant chemotherapy has shown to provide limited survival benefits when applied to all patients. Therefore, a biomarker to assess molecular residual disease (MRD) accurately and guide treatment selection is highly desirable for high-risk patients. This feasibility study evaluated the prognostic value of a tissue comprehensive genomic profiling (CGP)-informed, personalized circulating tumor DNA (ctDNA) assay (FoundationOne®Tracker) (Foundation Medicine, Inc., Cambridge, MA, USA) by correlating MRD status with clinical outcomes. ctDNA analysis was performed retrospectively on plasma samples from 69 patients with resected mCRC obtained at the MRD and the follow-up time point. Tissue CGP identified potentially actionable alterations in 54% (37/69) of patients. MRD-positivity was significantly associated with lower disease-free survival (DFS) (HR: 4.97, 95% CI: 2.67−9.24, p < 0.0001) and overall survival (OS) (HR: 27.05, 95% CI: 3.60−203.46, p < 0.0001). Similarly, ctDNA positive status at the follow-up time point correlated with a marked reduction in DFS (HR: 8.78, 95% CI: 3.59−21.49, p < 0.0001) and OS (HR: 20.06, 95% CI: 2.51−160.25, p < 0.0001). The overall sensitivity and specificity at the follow-up time point were 69% and 100%, respectively. Our results indicate that MRD detection using the tissue CGP-informed ctDNA assay is prognostic of survival outcomes in patients with resected mCRC. The concurrent MRD detection and identification of actionable alterations has the potential to guide perioperative clinical decision-making.

Changes in Circulating Tumor DNA Reflect Clinical Benefit Across Multiple Studies of Patients With Non-Small-Cell Lung Cancer Treated With Immune Checkpoint Inhibitors.

PURPOSE: As immune checkpoint inhibitors (ICI) become increasingly used in frontline settings, identifying early indicators of response is needed. Recent studies suggest a role for circulating tumor DNA (ctDNA) in monitoring response to ICI, but uncertainty exists in the generalizability of these studies. Here, the role of ctDNA for monitoring response to ICI is assessed through a standardized approach by assessing clinical trial data from five independent studies. PATIENTS AND METHODS: Patient-level clinical and ctDNA data were pooled and harmonized from 200 patients across five independent clinical trials investigating the treatment of patients with non-small-cell lung cancer with programmed cell death-1 (PD-1)/programmed death ligand-1 (PD-L1)-directed monotherapy or in combination with chemotherapy. CtDNA levels were measured using different ctDNA assays across the studies. Maximum variant allele frequencies were calculated using all somatic tumor-derived variants in each unique patient sample to correlate ctDNA changes with overall survival (OS) and progression-free survival (PFS). RESULTS: We observed strong associations between reductions in ctDNA levels from on-treatment liquid biopsies with improved OS (OS; hazard ratio, 2.28; 95% CI, 1.62 to 3.20; P < .001) and PFS (PFS; hazard ratio 1.76; 95% CI, 1.31 to 2.36; P < .001). Changes in the maximum variant allele frequencies ctDNA values showed strong association across different outcomes. CONCLUSION: In this pooled analysis of five independent clinical trials, consistent and robust associations between reductions in ctDNA and outcomes were found across multiple end points assessed in patients with non-small-cell lung cancer treated with an ICI. Additional tumor types, stages, and drug classes should be included in future analyses to further validate this. CtDNA may serve as an important tool in clinical development and an early indicator of treatment benefit.

Pan-cancer Analysis of Homologous Recombination Repair-associated Gene Alterations and Genome-wide Loss-of-Heterozygosity Score.

PURPOSE: To study associations across tumor types between genome-wide loss of heterozygosity (gLOH) and alterations in homologous recombination repair (HRR)-associated genes beyond BRCA1 and BRCA2. EXPERIMENTAL DESIGN: Genomic profiling using a targeted next-generation sequencing assay examining 324-465 genes (FoundationOne, FoundationOne Heme, and FoundationOne CDx; Foundation Medicine, Inc.) was performed in a cohort of 160,790 samples across different tumor types. Zygosity predictions and gLOH status were calculated and linked with alterations in 18 HRR-associated genes (BRCA1, BRCA2, PALB2, BARD1, ATR, ATRX, ATM, BAP1, RAD51B, RAD51C, RAD51D, BRIP1, NBN, CHEK1, CHEK2, FANCA, FANCC, MRE11) and other genomic features, using Fisher's exact test and Mann-Whitney U tests. RESULTS: We identified a strong correlation between elevated gLOH and biallelic alterations in a core set of HRR-associated genes beyond BRCA1 and BRCA2, such as BARD1, PALB2, FANCC, RAD51C, and RAD51D (particularly in breast, ovarian, pancreatic, and prostate cancer). Monoallelic/heterozygous alterations in HRR-associated genes were not associated with elevated gLOH. gLOH was also independently associated with TP53 loss. Co-occurrence of TP53 loss and alterations in HRR-associated genes, and combined loss of TP53-PTEN or TP53-RB1, was associated with a higher gLOH than each of the events separately. CONCLUSIONS: Biallelic alterations in core HRR-associated genes are frequent, strongly associated with elevated gLOH, and enriched in breast, ovarian, pancreatic, and prostate cancer. This analysis could inform the design of the next generation of clinical trials examining DNA repair-targeting agents, including PARP inhibitors.

Prdm9 and Meiotic Cohesin Proteins Cooperatively Promote DNA Double-Strand Break Formation in Mammalian Spermatocytes.

Meiotic recombination is required for correct segregation of chromosomes to gametes and to generate genetic diversity. In mice and humans, DNA double-strand breaks (DSBs) are initiated by SPO11 at recombination hotspots activated by PRDM9-catalyzed histone modifications on open chromatin. However, the DSB-initiating and repair proteins are associated with a linear proteinaceous scaffold called the chromosome axis, the core of which is composed of cohesin proteins. STAG3 is a stromalin subunit common to all meiosis-specific cohesin complexes. Mutations of meiotic cohesin proteins, especially STAG3, perturb both axis formation and recombination in the mouse, prompting determination of how the processes are mechanistically related. Protein interaction and genetic analyses revealed that PRDM9 interacts with STAG3 and REC8 in cooperative relationships that promote normal levels of meiotic DSBs at recombination hotspots in spermatocytes. The efficacy of the Prdm9-Stag3 genetic interaction in promoting DSB formation depends on PRDM9-mediated histone methyltransferase activity. Moreover, STAG3 deficiency has a major effect on DSB number even in the absence of PRDM9, showing that its role is not restricted to canonical PRDM9-activated hotspots. STAG3 and REC8 promote axis localization of the DSB-promoting proteins HORMAD1, IHO1, and MEI4, as well as SPO11 activity. These results establish that PRDM9 and axis-associated cohesin complexes together coordinate and facilitate meiotic recombination by recruiting key proteins for initiation of DSBs, thereby associating activated hotspots with DSB-initiating complexes on the axis.

Uncoupling of transcriptomic and cytological differentiation in mouse spermatocytes with impaired meiosis.

Cell differentiation is driven by changes in gene expression that manifest as changes in cellular phenotype or function. Altered cellular phenotypes, stemming from genetic mutations or other perturbations, are widely assumed to directly correspond to changes in the transcriptome and vice versa. Here, we exploited the cytologically well-defined Prdm9 mutant mouse as a model of developmental arrest to test whether parallel programs of cellular differentiation and gene expression are tightly coordinated, or can be disassociated. By comparing cytological phenotype markers and transcriptomes in wild-type and mutant spermatocytes, we identified multiple instances of cellular and molecular uncoupling in Prdm9-/- mutants. Most notably, although Prdm9-/- germ cells undergo cytological arrest in a late-leptotene/zygotene stage, they nevertheless develop gene expression signatures characteristic of later developmental substages. These findings suggest that transcriptomic changes may not reliably map to cellular phenotypes in developmentally perturbed systems.