Fumarate hydratase variant prevalence and manifestations among individuals receiving germline testing.

BACKGROUND: Germline variants in fumarate hydratase (FH) are associated with autosomal dominant (AD) hereditary leiomyomatosis and renal cell cancer (HLRCC) and autosomal recessive (AR) fumarase deficiency (FMRD). The prevalence and cancer penetrance across different FH variants remain unclear. METHODS: A database containing 120,061 records from individuals undergoing cancer germline testing was obtained. FH variants were classified into 3 categories: AD HLRCC variants, AR FMRD variants, and variants of unknown significance (VUSs). Individuals with variants from these categories were compared with those with negative genetic testing. RESULTS: FH variants were detected in 1.3% of individuals (AD HLRCC, 0.3%; AR FMRD, 0.4%; VUS, 0.6%). The rate of AD HLRCC variants discovered among reportedly asymptomatic individuals without a clear indication for HLRCC testing was 1 in 2668 (0.04%). In comparison with those with negative genetic testing, the renal cell carcinoma (RCC) prevalence was elevated with AD HLRCC variants (17.0% vs 4.5%; P < .01) and VUSs (6.4% vs 4.5%; P = .02) but not with AR FMRD variants. CONCLUSIONS: The prevalence of HLRCC discovered incidentally on germline testing is similar to recent population carrier estimates, and this suggests that this is a relatively common cancer syndrome. Compared with those with negative genetic testing, those with VUSs had an elevated risk of RCC, whereas those with AR FMRD variants did not.

Whole-genome sequencing of phenotypically distinct inflammatory breast cancers reveals similar genomic alterations to non-inflammatory breast cancers.

BACKGROUND: Inflammatory breast cancer (IBC) has a highly invasive and metastatic phenotype. However, little is known about its genetic drivers. To address this, we report the largest cohort of whole-genome sequencing (WGS) of IBC cases. METHODS: We performed WGS of 20 IBC samples and paired normal blood DNA to identify genomic alterations. For comparison, we used 23 matched non-IBC samples from the Cancer Genome Atlas Program (TCGA). We also validated our findings using WGS data from the International Cancer Genome Consortium (ICGC) and the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium. We examined a wide selection of genomic features to search for differences between IBC and conventional breast cancer. These include (i) somatic and germline single-nucleotide variants (SNVs), in both coding and non-coding regions; (ii) the mutational signature and the clonal architecture derived from these SNVs; (iii) copy number and structural variants (CNVs and SVs); and (iv) non-human sequence in the tumors (i.e., exogenous sequences of bacterial origin). RESULTS: Overall, IBC has similar genomic characteristics to non-IBC, including specific alterations, overall mutational load and signature, and tumor heterogeneity. In particular, we observed similar mutation frequencies between IBC and non-IBC, for each gene and most cancer-related pathways. Moreover, we found no exogenous sequences of infectious agents specific to IBC samples. Even though we could not find any strongly statistically distinguishing genomic features between the two groups, we did find some suggestive differences in IBC: (i) The MAST2 gene was more frequently mutated (20% IBC vs. 0% non-IBC). (ii) The TGF ß pathway was more frequently disrupted by germline SNVs (50% vs. 13%). (iii) Different copy number profiles were observed in several genomic regions harboring cancer genes. (iv) Complex SVs were more frequent. (v) The clonal architecture was simpler, suggesting more homogenous tumor-evolutionary lineages. CONCLUSIONS: Whole-genome sequencing of IBC manifests a similar genomic architecture to non-IBC. We found no unique genomic alterations shared in just IBCs; however, subtle genomic differences were observed including germline alterations in TGFß pathway genes and somatic mutations in the MAST2 kinase that could represent potential therapeutic targets.

Molecular medicine tumor board: whole-genome sequencing to inform on personalized medicine for a man with advanced prostate cancer.

PURPOSE: Molecular profiling of cancer is increasingly common as part of routine care in oncology, and germline and somatic profiling may provide insights and actionable targets for men with metastatic prostate cancer. However, all reported cases are of deidentified individuals without full medical and genomic data available in the public domain. PATIENT AND METHODS: We present a case of whole-genome tumor and germline sequencing in a patient with advanced prostate cancer, who has agreed to make his genomic and clinical data publicly available. RESULTS: We describe an 84-year-old Caucasian male with a Gleason 10 oligometastastic hormone-sensitive prostate cancer. Whole-genome sequencing provided insights into his tumor's underlying mutational processes and the development of an SPOP mutation. It also revealed an androgen-receptor dependency of his cancer which was reflected in his durable response to radiation and hormonal therapy. Potentially actionable genomic lesions in the tumor were identified through a personalized medicine approach for potential future therapy, but at the moment, he remains in remission, illustrating the hormonal sensitivity of his SPOP-driven prostate cancer. We also placed this patient in the context of a large prostate-cancer cohort from the PCAWG (Pan-cancer Analysis of Whole Genomes) group. In this comparison, the patient's cancer appears typical in terms of the number and type of somatic mutations, but it has a somewhat larger contribution from the mutational process associated with aging. CONCLUSION: We combined the expertise of medical oncology and genomics approaches to develop a molecular tumor board to integrate the care and study of this patient, who continues to have an outstanding response to his combined modality treatment. This identifiable case potentially helps overcome barriers to clinical and genomic data sharing.

Author Correction: Retrospective evaluation of whole exome and genome mutation calls in 746 cancer samples.

Correction to this paper has been published: https://doi.org/10.1038/s41467-020-20128-w.

Retrospective evaluation of whole exome and genome mutation calls in 746 cancer samples.

The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) curated consensus somatic mutation calls using whole exome sequencing (WES) and whole genome sequencing (WGS), respectively. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, which aggregated whole genome sequencing data from 2,658 cancers across 38 tumour types, we compare WES and WGS side-by-side from 746 TCGA samples, finding that ~80% of mutations overlap in covered exonic regions. We estimate that low variant allele fraction (VAF < 15%) and clonal heterogeneity contribute up to 68% of private WGS mutations and 71% of private WES mutations. We observe that ~30% of private WGS mutations trace to mutations identified by a single variant caller in WES consensus efforts. WGS captures both ~50% more variation in exonic regions and un-observed mutations in loci with variable GC-content. Together, our analysis highlights technological divergences between two reproducible somatic variant detection efforts.

An integrative ENCODE resource for cancer genomics.

ENCODE comprises thousands of functional genomics datasets, and the encyclopedia covers hundreds of cell types, providing a universal annotation for genome interpretation. However, for particular applications, it may be advantageous to use a customized annotation. Here, we develop such a custom annotation by leveraging advanced assays, such as eCLIP, Hi-C, and whole-genome STARR-seq on a number of data-rich ENCODE cell types. A key aspect of this annotation is comprehensive and experimentally derived networks of both transcription factors and RNA-binding proteins (TFs and RBPs). Cancer, a disease of system-wide dysregulation, is an ideal application for such a network-based annotation. Specifically, for cancer-associated cell types, we put regulators into hierarchies and measure their network change (rewiring) during oncogenesis. We also extensively survey TF-RBP crosstalk, highlighting how SUB1, a previously uncharacterized RBP, drives aberrant tumor expression and amplifies the effect of MYC, a well-known oncogenic TF. Furthermore, we show how our annotation allows us to place oncogenic transformations in the context of a broad cell space; here, many normal-to-tumor transitions move towards a stem-like state, while oncogene knockdowns show an opposing trend. Finally, we organize the resource into a coherent workflow to prioritize key elements and variants, in addition to regulators. We showcase the application of this prioritization to somatic burdening, cancer differential expression and GWAS. Targeted validations of the prioritized regulators, elements and variants using siRNA knockdowns, CRISPR-based editing, and luciferase assays demonstrate the value of the ENCODE resource.

Using sigLASSO to optimize cancer mutation signatures jointly with sampling likelihood.

Multiple mutational processes drive carcinogenesis, leaving characteristic signatures in tumor genomes. Determining the active signatures from a full repertoire of potential ones helps elucidate mechanisms of cancer development. This involves optimally decomposing the counts of cancer mutations, tabulated according to their trinucleotide context, into a linear combination of known signatures. Here, we develop sigLASSO (a software tool at github.com/gersteinlab/siglasso) to carry out this optimization efficiently. sigLASSO has four key aspects: (1) It jointly optimizes the likelihood of sampling and signature fitting, by explicitly factoring multinomial sampling into the objective function. This is particularly important when mutation counts are low and sampling variance is high (e.g., in exome sequencing). (2) sigLASSO uses L1 regularization to parsimoniously assign signatures, leading to sparse and interpretable solutions. (3) It fine-tunes model complexity, informed by data scale and biological priors. (4) Consequently, sigLASSO can assess model uncertainty and abstain from making assignments in low-confidence contexts.

Estimation of the carrier frequency of fumarate hydratase alterations and implications for kidney cancer risk in hereditary leiomyomatosis and renal cancer.

BACKGROUND: Hereditary leiomyomatosis and renal cancer (HLRCC) is a cancer syndrome associated with a germline mutation in fumarate hydratase (FH). The syndrome is associated with cutaneous and uterine leiomyomas, and some patients develop a lethal form of kidney cancer. This study provides estimates for the FH carrier frequency and kidney cancer penetrance. METHODS: Data sets containing sequencing data for the FH gene were used: the 1000 Genomes Project (1000GP) and the Exome Aggregation Consortium (ExAC). Alterations in the FH gene were characterized on the basis of different variant risk tiers: 1) ClinVar annotated variants, 2) loss-of-function alterations, and 3) highly impactful missense alterations. The cumulative incidence of FH alterations overall and by different world populations was evaluated in 1000GP and ExAC. A lifetime penetrance of HLRCC kidney cancer risk was generated with 3 estimates of the annual incidence. RESULTS: The overall allele frequencies of tier 1 to 3 FH alterations in the ExAC and 1000GP data sets were 2.54 × 10-3 (1 in 393) and 1.20 × 10-3 (1 in 835), respectively. There were differences in the allele frequencies of FH alterations between world populations. Based on various estimates of the percentage of kidney cancers with FH alterations, the lifetime kidney cancer penetrance for carrier estimate 3 in ExAC was 1.7% to 5.8%. CONCLUSIONS: FH alterations are common and are carried by approximately 1 in 1000 individuals according to the more conservative estimates. The lifetime kidney cancer penetrance appears lower than previously estimated. Although databases are not population cohorts, they provide a useful quantitative estimate of rare variants with low penetrance.

Passenger Mutations in More Than 2,500 Cancer Genomes: Overall Molecular Functional Impact and Consequences.

The dichotomous model of "drivers" and "passengers" in cancer posits that only a few mutations in a tumor strongly affect its progression, with the remaining ones being inconsequential. Here, we leveraged the comprehensive variant dataset from the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) project to demonstrate that-in addition to the dichotomy of high- and low-impact variants-there is a third group of medium-impact putative passengers. Moreover, we also found that molecular impact correlates with subclonal architecture (i.e., early versus late mutations), and different signatures encode for mutations with divergent impact. Furthermore, we adapted an additive-effects model from complex-trait studies to show that the aggregated effect of putative passengers, including undetected weak drivers, provides significant additional power (∼12% additive variance) for predicting cancerous phenotypes, beyond PCAWG-identified driver mutations. Finally, this framework allowed us to estimate the frequency of potential weak-driver mutations in PCAWG samples lacking any well-characterized driver alterations.

Whole-genome analysis of papillary kidney cancer finds significant noncoding alterations.

To date, studies on papillary renal-cell carcinoma (pRCC) have largely focused on coding alterations in traditional drivers, particularly the tyrosine-kinase, Met. However, for a significant fraction of tumors, researchers have been unable to determine a clear molecular etiology. To address this, we perform the first whole-genome analysis of pRCC. Elaborating on previous results on MET, we find a germline SNP (rs11762213) in this gene predicting prognosis. Surprisingly, we detect no enrichment for small structural variants disrupting MET. Next, we scrutinize noncoding mutations, discovering potentially impactful ones associated with MET. Many of these are in an intron connected to a known, oncogenic alternative-splicing event; moreover, we find methylation dysregulation nearby, leading to a cryptic promoter activation. We also notice an elevation of mutations in the long noncoding RNA NEAT1, and these mutations are associated with increased expression and unfavorable outcome. Finally, to address the origin of pRCC heterogeneity, we carry out whole-genome analyses of mutational processes. First, we investigate genome-wide mutational patterns, finding they are governed mostly by methylation-associated C-to-T transitions. We also observe significantly more mutations in open chromatin and early-replicating regions in tumors with chromatin-modifier alterations. Finally, we reconstruct cancer-evolutionary trees, which have markedly different topologies and suggested evolutionary trajectories for the different subtypes of pRCC.