Evaluation of tissue- and plasma-derived tumor mutational burden (TMB) and genomic alterations of interest in CheckMate 848, a study of nivolumab combined with ipilimumab and nivolumab alone in patients with advanced or metastatic solid tumors with high TMB.

BACKGROUND: An accumulation of somatic mutations in tumors leads to increased neoantigen levels and antitumor immune response. Tumor mutational burden (TMB) reflects the rate of somatic mutations in the tumor genome, as determined from tumor tissue (tTMB) or blood (bTMB). While high tTMB is a biomarker of immune checkpoint inhibitor (ICI) treatment efficacy, few studies have explored the clinical utility of bTMB, a less invasive alternative for TMB assessment. Establishing the correlation between tTMB and bTMB would provide insight into whether bTMB is a potential substitute for tTMB. We explored the tumor genomes of patients enrolled in CheckMate 848 with measurable TMB. The correlation between tTMB and bTMB, and the factors affecting it, were evaluated. METHODS: In the phase 2 CheckMate 848 (NCT03668119) study, immuno-oncology-naïve patients with advanced, metastatic, or unresectable solid tumors and tTMB-high or bTMB-high (≥10 mut/Mb) were prospectively randomized 2:1 to receive nivolumab plus ipilimumab or nivolumab monotherapy. Tissue and plasma DNA sequencing was performed using the Foundation Medicine FoundationOne CDx and bTMB Clinical Trial Assays, respectively. tTMB was quantified from coding variants, insertions, and deletions, and bTMB from somatic base substitutions. Correlations between tTMB and bTMB were determined across samples and with respect to maximum somatic allele frequency (MSAF). Assay agreement and variant composition were also evaluated. RESULTS: A total of 1,438 and 1,720 unique tissue and blood samples, respectively, were obtained from 1,954 patients and included >100 screened disease ontologies, with 1,017 unique pairs of tTMB and bTMB measurements available for assessment. Median tTMB and bTMB were 3.8 and 3.5 mut/Mb, respectively. A significant correlation between tTMB and bTMB (r=0.48, p<0.0001) was observed across all sample pairs, which increased to r=0.54 (p<0.0001) for samples with MSAF≥1%. Assay concordance was highest for samples with MSAF≥10% across multiple disease ontologies and observed for both responders and non-responders to ICI therapy. The variants contributing to tTMB and bTMB were similar. CONCLUSIONS: We observed that tTMB and bTMB had a statistically significant correlation, particularly for samples with high MSAF, and that this correlation applied across disease ontologies. Further investigation into the clinical utility of bTMB is warranted.

Misannotated Multi-Nucleotide Variants in Public Cancer Genomics Datasets Lead to Inaccurate Mutation Calls with Significant Implications.

Although next-generation sequencing is widely used in cancer to profile tumors and detect variants, most somatic variant callers used in these pipelines identify variants at the lowest possible granularity, single-nucleotide variants (SNV). As a result, multiple adjacent SNVs are called individually instead of as a multi-nucleotide variants (MNV). With this approach, the amino acid change from the individual SNV within a codon could be different from the amino acid change based on the MNV that results from combining SNV, leading to incorrect conclusions about the downstream effects of the variants. Here, we analyzed 10,383 variant call files (VCF) from the Cancer Genome Atlas (TCGA) and found 12,141 incorrectly annotated MNVs. Analysis of seven commonly mutated genes from 178 studies in cBioPortal revealed that MNVs were consistently missed in 20 of these studies, whereas they were correctly annotated in 15 more recent studies. At the BRAF V600 locus, the most common example of MNV, several public datasets reported separate BRAF V600E and BRAF V600M variants instead of a single merged V600K variant. VCFs from the TCGA Mutect2 caller were used to develop a solution to merge SNV to MNV. Our custom script used the phasing information from the SNV VCF and determined whether SNVs were at the same codon and needed to be merged into MNV before variant annotation. This study shows that institutions performing NGS sequencing for cancer genomics should incorporate the step of merging MNV as a best practice in their pipelines. SIGNIFICANCE: Identification of incorrect mutation calls in TCGA, including clinically relevant BRAF V600 and KRAS G12, will influence research and potentially clinical decisions.

Modeling performance of sample collection sites using whole exome sequencing metrics.

Although next-generation sequencing assays are routinely carried out using samples from cancer trials, the sequencing data are not always of the required quality. There is a need to evaluate the performance of tissue collection sites and provide feedback about the quality of next-generation sequencing data. This study used a modeling approach based on whole exome sequencing quality control (QC) metrics to evaluate the relative performance of sites participating in the Bristol Myers Squibb Immuno-Oncology clinical trials sample collection. We identified several events for the sample swap. Overall, most sites performed well and few showed poor performance. These findings can increase awareness of sample failure and improve the quality of samples.

Gluconeogenesis and PEPCK are critical components of healthy aging and dietary restriction life extension.

High glucose diets are unhealthy, although the mechanisms by which elevated glucose is harmful to whole animal physiology are not well understood. In Caenorhabditis elegans, high glucose shortens lifespan, while chemically inflicted glucose restriction promotes longevity. We investigated the impact of glucose metabolism on aging quality (maintained locomotory capacity and median lifespan) and found that, in addition to shortening lifespan, excess glucose negatively impacts locomotory healthspan. Conversely, disrupting glucose utilization by knockdown of glycolysis-specific genes results in large mid-age physical improvements via a mechanism that requires the FOXO transcription factor DAF-16. Adult locomotory capacity is extended by glycolysis disruption, but maximum lifespan is not, indicating that limiting glycolysis can increase the proportion of life spent in mobility health. We also considered the largely ignored role of glucose biosynthesis (gluconeogenesis) in adult health. Directed perturbations of gluconeogenic genes that specify single direction enzymatic reactions for glucose synthesis decrease locomotory healthspan, suggesting that gluconeogenesis is needed for healthy aging. Consistent with this idea, overexpression of the central gluconeogenic gene pck-2 (encoding PEPCK) increases health measures via a mechanism that requires DAF-16 to promote pck-2 expression in specific intestinal cells. Dietary restriction also features DAF-16-dependent pck-2 expression in the intestine, and the healthspan benefits conferred by dietary restriction require pck-2. Together, our results describe a new paradigm in which nutritional signals engage gluconeogenesis to influence aging quality via DAF-16. These data underscore the idea that promotion of gluconeogenesis might be an unappreciated goal for healthy aging and could constitute a novel target for pharmacological interventions that counter high glucose consequences, including diabetes.

C. elegans Heterochromatin Factor SET-32 Plays an Essential Role in Transgenerational Establishment of Nuclear RNAi-Mediated Epigenetic Silencing.

The dynamic process by which nuclear RNAi engages a transcriptionally active target, before the repressive state is stably established, remains largely a mystery. Here, we found that the onset of exogenous dsRNA-induced nuclear RNAi in C. elegans is a transgenerational process, and it requires a putative histone methyltransferase (HMT), SET-32. By developing a CRISPR-based genetic approach, we found that silencing establishment at the endogenous targets of germline nuclear RNAi also requires SET-32. Although SET-32 and two H3K9 HMTs, MET-2 and SET-25, are dispensable for the maintenance of silencing, they do contribute to transcriptional repression in mutants that lack the germline nuclear Argonaute protein HRDE-1, suggesting a conditional role of heterochromatin in the maintenance phase. Our study indicates that (1) establishment and maintenance of siRNA-guided transcriptional repression are two distinct processes with different genetic requirements and (2) the rate-limiting step of the establishment phase is a transgenerational, chromatin-based process.

The spatial and temporal dynamics of nuclear RNAi-targeted retrotransposon transcripts in Caenorhabditis elegans.

Nuclear RNA interference provides a unique approach to the study of RNA-mediated transgenerational epigenetic inheritance. A paradox in the field is that expression of target loci is necessary for the initiation and maintenance of their silencing. How expression and repression are coordinated during animal development is poorly understood. To resolve this gap, we took imaging, deep-sequencing and genetic approaches towards delineating the developmental regulation and subcellular localization of RNA transcripts of two representative endogenous targets, the LTR retrotransposons Cer3 and Cer8. By examining wild-type worms and a collection of mutant strains, we found that the expression and silencing cycle of Cer3 and Cer8 is coupled with embryonic and germline development. Strikingly, endogenous targets exhibit a hallmark of nuclear enrichment of their RNA transcripts. In addition, germline and somatic repressions of Cer3 have different genetic requirements for three heterochromatin enzymes, MET-2, SET-25 and SET-32, in conjunction with the nuclear Argonaute protein HRDE-1. These results provide the first comprehensive cellular and developmental characterization of nuclear RNAi activities throughout the animal reproductive cycle.

MORC-1 Integrates Nuclear RNAi and Transgenerational Chromatin Architecture to Promote Germline Immortality.

Germline-expressed endogenous small interfering RNAs (endo-siRNAs) transmit multigenerational epigenetic information to ensure fertility in subsequent generations. In Caenorhabditis elegans, nuclear RNAi ensures robust inheritance of endo-siRNAs and deposition of repressive H3K9me3 marks at target loci. How target silencing is maintained in subsequent generations is poorly understood. We discovered that morc-1 is essential for transgenerational fertility and acts as an effector of endo-siRNAs. Unexpectedly, morc-1 is dispensable for siRNA inheritance but is required for target silencing and maintenance of siRNA-dependent chromatin organization. A forward genetic screen identified mutations in met-1, which encodes an H3K36 methyltransferase, as potent suppressors of morc-1(-) and nuclear RNAi mutant phenotypes. Further analysis of nuclear RNAi and morc-1(-) mutants revealed a progressive, met-1-dependent enrichment of H3K36me3, suggesting that robust fertility requires repression of MET-1 activity at nuclear RNAi targets. Without MORC-1 and nuclear RNAi, MET-1-mediated encroachment of euchromatin leads to detrimental decondensation of germline chromatin and germline mortality.

Decoupling the downstream effects of germline nuclear RNAi reveals that H3K9me3 is dispensable for heritable RNAi and the maintenance of endogenous siRNA-mediated transcriptional silencing in Caenorhabditis elegans.

BACKGROUND: Germline nuclear RNAi in C. elegans is a transgenerational gene-silencing pathway that leads to H3K9 trimethylation (H3K9me3) and transcriptional silencing at the target genes. H3K9me3 induced by either exogenous double-stranded RNA (dsRNA) or endogenous siRNA (endo-siRNA) is highly specific to the target loci and transgenerationally heritable. Despite these features, the role of H3K9me3 in siRNA-mediated transcriptional silencing and inheritance of the silencing state at native target genes is unclear. In this study, we took combined genetic and whole-genome approaches to address this question. RESULTS: Here we demonstrate that siRNA-mediated H3K9me3 requires combined activities of three H3K9 histone methyltransferases: MET-2, SET-25, and SET-32. set-32 single, met-2 set-25 double, and met-2 set-25;set-32 triple mutant adult animals all exhibit prominent reductions in H3K9me3 throughout the genome, with met-2 set-25;set-32 mutant worms losing all detectable H3K9me3 signals. Surprisingly, loss of high-magnitude H3K9me3 at the native nuclear RNAi targets has no effect on the transcriptional silencing state. In addition, the exogenous dsRNA-induced transcriptional silencing and heritable RNAi at oma-1, a well-established nuclear RNAi reporter gene, are completely resistant to the loss of H3K9me3. CONCLUSIONS: Nuclear RNAi-mediated H3K9me3 in C. elegans requires multiple histone methyltransferases, including MET-2, SET-25, and SET-32. H3K9me3 is not essential for dsRNA-induced heritable RNAi or the maintenance of endo-siRNA-mediated transcriptional silencing in C. elegans. We propose that siRNA-mediated transcriptional silencing in C. elegans can be maintained by an H3K9me3-independent mechanism.

A transgenerational role of the germline nuclear RNAi pathway in repressing heat stress-induced transcriptional activation in C. elegans.

BACKGROUND: Environmental stress-induced transgenerational epigenetic effects have been observed in various model organisms and human. The capacity and mechanism of such phenomena are poorly understood. In C. elegans, siRNA mediates transgenerational gene silencing through the germline nuclear RNAi pathway. This pathway is also required to maintain the germline immortality when C. elegans is under heat stress. However, the underlying molecular mechanism is unknown. In this study, we investigated the impact of heat stress on chromatin, transcription, and siRNAs at the whole-genome level, and whether any of the heat-induced effects is transgenerationally heritable in either the wild-type or the germline nuclear RNAi mutant animals. RESULTS: We performed 12-generation temperature-shift experiments using the wild-type C. elegans and a mutant strain that lacks the germline-specific nuclear Argonaute protein HRDE-1/WAGO-9. By examining the mRNA, small RNA, RNA polymerase II, and H3K9 trimethylation profiles at the whole-genome level, we revealed an epigenetic role of HRDE-1 in repressing heat stress-induced transcriptional activation of over 280 genes. Many of these genes are in or near LTR (long-terminal repeat) retrotransposons. Strikingly, for some of these genes, the heat stress-induced transcriptional activation in the hrde-1 mutant intensifies in the late generations under the heat stress and is heritable for at least two generations after the mutant animals are shifted back to lower temperature. hrde-1 mutation also leads to siRNA expression changes of many genes. This effect on siRNA is dependent on both the temperature and generation. CONCLUSIONS: Our study demonstrated that a large number of the endogenous targets of the germline nuclear RNAi pathway in C. elegans are sensitive to heat-induced transcriptional activation. This effect at certain genomic loci including LTR retrotransposons is transgenerational. Germline nuclear RNAi antagonizes this temperature effect at the transcriptional level and therefore may play a key role in heat stress response in C. elegans.