Genome-wide somatic variant calling using localized colored de Bruijn graphs

ONCOLINER: A new solution for monitoring, improving, and harmonizing somatic variant calling across genomic oncology centers

The characterization of somatic genomic variation associated with the biology of tumors is fundamental for cancer research and personalized medicine, as it guides the reliability and impact of cancer studies and genomic-based decisions in clinical oncology. However, the quality and scope of tumor genome analysis across cancer research centers and hospitals are currently highly heterogeneous, limiting the consistency of tumor diagnoses across hospitals and the possibilities of data sharing and data integration across studies. With the aim of providing users with actionable and personalized recommendations for the overall enhancement and harmonization of somatic variant identification across research and clinical environments, we have developed ONCOLINER. Using specifically designed mosaic and tumorized genomes for the analysis of recall and precision across somatic SNVs, insertions or deletions (indels), and structural variants (SVs), we demonstrate that ONCOLINER is capable of improving and harmonizing genome analysis across three state-of-the-art variant discovery pipelines in genomic oncology.

The Open Pediatric Cancer Project

Background

In 2019, the Open Pediatric Brain Tumor Atlas (OpenPBTA) was created as a global, collaborative open-science initiative to genomically characterize 1,074 pediatric brain tumors and 22 patient-derived cell lines. Here, we extend the OpenPBTA to create the Open Pediatric Cancer (OpenPedCan) Project, a harmonized open-source multi-omic dataset from 6,112 pediatric cancer patients with 7,096 tumor events across more than 100 histologies. Combined with RNA-Seq from the Genotype-Tissue Expression (GTEx) and The Cancer Genome Atlas (TCGA), OpenPedCan contains nearly 48,000 total biospecimens (24,002 tumor and 23,893 normal specimens).

Findings

We utilized Gabriella Miller Kids First (GMKF) workflows to harmonize WGS, WXS, RNA-seq, and Targeted Sequencing datasets to include somatic SNVs, InDels, CNVs, SVs, RNA expression, fusions, and splice variants. We integrated summarized CPTAC whole cell proteomics and phospho-proteomics data, miRNA-Seq data, and have developed a methylation array harmonization workflow to include m-values, beta-vales, and copy number calls. OpenPedCan contains reproducible, dockerized workflows in GitHub, CAVATICA, and Amazon Web Services (AWS) to deliver harmonized and processed data from over 60 scalable modules which can be leveraged both locally and on AWS. The processed data are released in a versioned manner and accessible through CAVATICA or AWS S3 download (from GitHub), and queryable through PedcBioPortal and the NCI's pediatric Molecular Targets Platform. Notably, we have expanded PBTA molecular subtyping to include methylation information to align with the WHO 2021 Central Nervous System Tumor classifications, allowing us to create research- grade integrated diagnoses for these tumors.

Conclusions

OpenPedCan data and its reproducible analysis module framework are openly available and can be utilized and/or adapted by researchers to accelerate discovery, validation, and clinical translation.

Increased Synthetic Cytotoxicity of Combinatorial Chemoradiation Therapy in Homologous Recombination Deficient Tumors

PURPOSE

Homologous recombination deficient (HRD) tumors are exquisitely sensitive to platinum-based chemotherapy and when combined with radiation therapy (RT), leads to improved overall survival in multiple cancer types. Whether a subset of tumors with distinct molecular characteristics demonstrate increased benefit from cisplatin and RT (c-RT) is unclear. We hypothesized that HRD tumors, whether associated with BRCA mutations or genomic scars of HRD, exhibit exquisite sensitivity to c-RT, and that HRD may be a significant driver of c-RT benefit.

METHODS AND MATERIALS

Sensitivity to c-RT was examined using isogenic and sporadic breast cancer cell lines. HRD was assessed using 4 assays: RT-induced Rad51 foci, a DR-GFP reporter assay, a genomic scar score (large-scale state transitions [LST]), and clonogenic survival assays. Whole-genome sequencing of 4 breast tumors from a phase 2 clinical trial of neoadjuvant c-RT in triple-negative breast cancer was performed and HRD was defined using HRDetect.

RESULTS

BRCA1/2 deficient cell lines displayed functional HRD based on the Rad51 functional assay, with c-RT to RT or cisplatin interaction ratios (IR) of 1.11 and 26.84 for the BRCA1 isogenic pair at 2 μM cisplatin and 6 Gy, respectively. The highest LST lines demonstrated HRD and synthetic cytotoxicity to c-RT with IR at 2 Gy and cisplatin 20 μM of 7.50, and the lowest LST line with IR of 0.65. Of 4 evaluable patients in the phase 2 trial, one achieved a pathologic complete response with corresponding HRD based on multiple genomic scar scores including HRDetect and LST scores, compared with patients without a pathologic complete response.

CONCLUSIONS

HRD breast cancers, whether identified by BRCA1/2 mutation status, functional tests, or mutational signatures, appear to be significantly more sensitive to c-RT compared with isogenic controls or tumors without HRD mutational signatures. HRD tumors may be exquisitely sensitive to c-RT which warrants further clinical investigation to guide a precision oncology approach.

Reference-free structural variant detection in microbiomes via long-read co-assembly graphs

MOTIVATION

The study of bacterial genome dynamics is vital for understanding the mechanisms underlying microbial adaptation, growth, and their impact on host phenotype. Structural variants (SVs), genomic alterations of 50 base pairs or more, play a pivotal role in driving evolutionary processes and maintaining genomic heterogeneity within bacterial populations. While SV detection in isolate genomes is relatively straightforward, metagenomes present broader challenges due to the absence of clear reference genomes and the presence of mixed strains. In response, our proposed method rhea, forgoes reference genomes and metagenome-assembled genomes (MAGs) by encompassing all metagenomic samples in a series (time or other metric) into a single co-assembly graph. The log fold change in graph coverage between successive samples is then calculated to call SVs that are thriving or declining.

RESULTS

We show rhea to outperform existing methods for SV and horizontal gene transfer (HGT) detection in two simulated mock metagenomes, particularly as the simulated reads diverge from reference genomes and an increase in strain diversity is incorporated. We additionally demonstrate use cases for rhea on series metagenomic data of environmental and fermented food microbiomes to detect specific sequence alterations between successive time and temperature samples, suggesting host advantage. Our approach leverages previous work in assembly graph structural and coverage patterns to provide versatility in studying SVs across diverse and poorly characterized microbial communities for more comprehensive insights into microbial gene flux.

AVAILABILITY AND IMPLEMENTATION

rhea is open source and available at: https://github.com/treangenlab/rhea.

Adenoid cystic carcinoma of the Bartholin's gland is underpinned by MYB- and MYBL1- rearrangements

OBJECTIVE

Adenoid cystic carcinoma (AdCC) of the Bartholin's gland (AdCC-BG) is a very rare gynecologic vulvar malignancy. AdCC-BGs are slow-growing but locally aggressive and are associated with high recurrence rates. Here we sought to characterize the molecular underpinning of AdCC-BGs.

METHODS

AdCC-BGs (n = 6) were subjected to a combination of RNA-sequencing, targeted DNA-sequencing, reverse-transcription PCR, fluorescence in situ hybridization (FISH) and MYB immunohistochemistry (IHC). Clinicopathologic variables, somatic mutations, copy number alterations and chimeric transcripts were assessed.

RESULTS

All six AdCC-BGs were biphasic, composed of ductal and myoepithelial cells. Akin to salivary gland and breast AdCCs, three AdCC-BGs had the MYB::NFIB fusion gene with varying breakpoints, all of which were associated with MYB overexpression by IHC. Two AdCC-BGs were underpinned by MYBL1 fusion genes with different gene partners, including MYBL1::RAD51B and MYBL1::EWSR1 gene fusions, and showed MYB protein expression. Although the final AdCC-BG studied had MYB protein overexpression, no gene fusion was identified. AdCC-BGs harbored few additional somatic genetic alterations, and only few mutations in cancer-related genes were identified, including GNAQ, GNAS, KDM6A, AKT1 and BCL2, none of which were recurrent. Two AdCC-BGs, both with a MYB::NFIB fusion gene, developed metastatic disease.

CONCLUSIONS

AdCC-BGs constitute a convergent phenotype, whereby activation of MYB or MYBL1 can be driven by the MYB::NFIB fusion gene or MYBL1 rearrangements. Our observations further support the notion that AdCCs, irrespective of organ site, constitute a genotypic-phenotypic correlation. Assessment of MYB or MYBL1 rearrangements may be used as an ancillary marker for the diagnosis of AdCC-BGs.

Integrated approach to generate artificial samples with low tumor fraction for somatic variant calling benchmarking

BACKGROUND

High-throughput sequencing (HTS) has become the gold standard approach for variant analysis in cancer research. However, somatic variants may occur at low fractions due to contamination from normal cells or tumor heterogeneity; this poses a significant challenge for standard HTS analysis pipelines. The problem is exacerbated in scenarios with minimal tumor DNA, such as circulating tumor DNA in plasma. Assessing sensitivity and detection of HTS approaches in such cases is paramount, but time-consuming and expensive: specialized experimental protocols and a sufficient quantity of samples are required for processing and analysis. To overcome these limitations, we propose a new computational approach specifically designed for the generation of artificial datasets suitable for this task, simulating ultra-deep targeted sequencing data with low-fraction variants and demonstrating their effectiveness in benchmarking low-fraction variant calling.

RESULTS

Our approach enables the generation of artificial raw reads that mimic real data without relying on pre-existing data by using NEAT, a fine-grained read simulator that generates artificial datasets using models learned from multiple different datasets. Then, it incorporates low-fraction variants to simulate somatic mutations in samples with minimal tumor DNA content. To prove the suitability of the created artificial datasets for low-fraction variant calling benchmarking, we used them as ground truth to evaluate the performance of widely-used variant calling algorithms: they allowed us to define tuned parameter values of major variant callers, considerably improving their detection of very low-fraction variants.

CONCLUSIONS

Our findings highlight both the pivotal role of our approach in creating adequate artificial datasets with low tumor fraction, facilitating rapid prototyping and benchmarking of algorithms for such dataset type, as well as the important need of advancing low-fraction variant calling techniques.

APOBEC3 mutagenesis drives therapy resistance in breast cancer

Acquired genetic alterations commonly drive resistance to endocrine and targeted therapies in metastatic breast cancer 1-7 , however the underlying processes engendering these diverse alterations are largely uncharacterized. To identify the mutational processes operant in breast cancer and their impact on clinical outcomes, we utilized a well-annotated cohort of 3,880 patient samples with paired tumor-normal sequencing data. The mutational signatures associated with apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3) enzymes were highly prevalent and enriched in post-treatment compared to treatment-naïve hormone receptor-positive (HR+) cancers. APOBEC3 mutational signatures were independently associated with shorter progression-free survival on antiestrogen plus CDK4/6 inhibitor combination therapy in patients with HR+ metastatic breast cancer. Whole genome sequencing (WGS) of breast cancer models and selected paired primary-metastatic samples demonstrated that active APOBEC3 mutagenesis promoted resistance to both endocrine and targeted therapies through characteristic alterations such as RB1 loss-of-function mutations. Evidence of APOBEC3 activity in pre-treatment samples illustrated a pervasive role for this mutational process in breast cancer evolution. The study reveals APOBEC3 mutagenesis to be a frequent mediator of therapy resistance in breast cancer and highlights its potential as a biomarker and target for overcoming resistance.

Converging and evolving immuno-genomic routes toward immune escape in breast cancer

The interactions between tumor and immune cells along the course of breast cancer progression remain largely unknown. Here, we extensively characterize multiple sequential and parallel multiregion tumor and blood specimens of an index patient and a cohort of metastatic triple-negative breast cancers. We demonstrate that a continuous increase in tumor genomic heterogeneity and distinct molecular clocks correlated with resistance to treatment, eventually allowing tumors to escape from immune control. TCR repertoire loses diversity over time, leading to convergent evolution as breast cancer progresses. Although mixed populations of effector memory and cytotoxic single T cells coexist in the peripheral blood, defects in the antigen presentation machinery coupled with subdued T cell recruitment into metastases are observed, indicating a potent immune avoidance microenvironment not compatible with an effective antitumor response in lethal metastatic disease. Our results demonstrate that the immune responses against cancer are not static, but rather follow dynamic processes that match cancer genomic progression, illustrating the complex nature of tumor and immune cell interactions.

Validation and benchmarking of targeted panel sequencing for cancer genomic profiling

OBJECTIVES

To validate a large next-generation sequencing (NGS) panel for comprehensive genomic profiling and improve patient access to more effective precision oncology treatment strategies.

METHODS

OncoPanScan was designed by targeting 825 cancer-related genes to detect a broad range of genomic alterations. A practical validation strategy was used to evaluate the assay's analytical performance, involving 97 tumor specimens with 25 paired blood specimens, 10 engineered cell lines, and 121 artificial reference DNA samples.

RESULTS

Overall, 1107 libraries were prepared and the sequencing failure rate was 0.18%. Across alteration classes, sensitivity ranged from 0.938 to more than 0.999, specificity ranged from 0.889 to more than 0.999, positive predictive value ranged from 0.867 to more than 0.999, repeatability ranged from 0.908 to more than 0.999, and reproducibility ranged from 0.832 to more than 0.999. The limit of detection for variants was established based on variant frequency, while for tumor mutation burden and microsatellite instability, it was based on tumor content, resulting in a minimum requirement of 20% tumor content. Benchmarking variant calls against validated NGS assays revealed that variations in the dry-bench processes were the primary cause of discordances.

CONCLUSIONS

This study presents a detailed validation framework and empirical recommendations for large panel validation and elucidates the sources of discordant alteration calls by comparing with "gold standard measures."

Evolution of structural rearrangements in prostate cancer intracranial metastases

Intracranial metastases in prostate cancer are uncommon but clinically aggressive. A detailed molecular characterization of prostate cancer intracranial metastases would improve our understanding of their pathogenesis and the search for new treatment strategies. We evaluated the clinical and molecular characteristics of 36 patients with metastatic prostate cancer to either the dura or brain parenchyma. We performed whole genome sequencing (WGS) of 10 intracranial prostate cancer metastases, as well as WGS of primary prostate tumors from men who later developed metastatic disease (n = 6) and nonbrain prostate cancer metastases (n = 36). This first whole genome sequencing study of prostate intracranial metastases led to several new insights. First, there was a higher diversity of complex structural alterations in prostate cancer intracranial metastases compared to primary tumor tissues. Chromothripsis and chromoplexy events seemed to dominate, yet there were few enrichments of specific categories of structural variants compared with non-brain metastases. Second, aberrations involving the AR gene, including AR enhancer gain were observed in 7/10 (70%) of intracranial metastases, as well as recurrent loss of function aberrations involving TP53 in 8/10 (80%), RB1 in 2/10 (20%), BRCA2 in 2/10 (20%), and activation of the PI3K/AKT/PTEN pathway in 8/10 (80%). These alterations were frequently present in tumor tissues from other sites of disease obtained concurrently or sequentially from the same individuals. Third, clonality analysis points to genomic factors and evolutionary bottlenecks that contribute to metastatic spread in patients with prostate cancer. These results describe the aggressive molecular features underlying intracranial metastasis that may inform future diagnostic and treatment approaches.

Identifying novel mechanisms of biallelic TP53 loss refines poor outcome for patients with multiple myeloma

Biallelic TP53 inactivation is the most important high-risk factor associated with poor survival in multiple myeloma. Classical biallelic TP53 inactivation has been defined as simultaneous mutation and copy number loss in most studies; however, numerous studies have demonstrated that other factors could lead to the inactivation of TP53. Here, we hypothesized that novel biallelic TP53 inactivated samples existed in the multiple myeloma population. A random forest regression model that exploited an expression signature of 16 differentially expressed genes between classical biallelic TP53 and TP53 wild-type samples was subsequently established and used to identify novel biallelic TP53 samples from monoallelic TP53 groups. The model reflected high accuracy and robust performance in newly diagnosed relapsed and refractory populations. Patient survival of classical and novel biallelic TP53 samples was consistently much worse than those with mono-allelic or wild-type TP53 status. We also demonstrated that some predicted biallelic TP53 samples simultaneously had copy number loss and aberrant splicing, resulting in overexpression of high-risk transcript variants, leading to biallelic inactivation. We discovered that splice site mutation and overexpression of the splicing factor MED18 were reasons for aberrant splicing. Taken together, our study unveiled the complex transcriptome of TP53, some of which might benefit future studies targeting abnormal TP53.

In Silico Evaluation of Variant Calling Methods for Bacterial Whole-Genome Sequencing Assays

Identification and analysis of clinically relevant strains of bacteria increasingly relies on whole-genome sequencing. The downstream bioinformatics steps necessary for calling variants from short-read sequences are well-established but seldom validated against haploid genomes. We devised an in silico workflow to introduce single nucleotide polymorphisms (SNP) and indels into bacterial reference genomes, and computationally generate sequencing reads based on the mutated genomes. We then applied the method to Mycobacterium tuberculosis H37Rv, Staphylococcus aureus NCTC 8325, and Klebsiella pneumoniae HS11286, and used the synthetic reads as truth sets for evaluating several popular variant callers. Insertions proved especially challenging for most variant callers to correctly identify, relative to deletions and single nucleotide polymorphisms. With adequate read depth, however, variant callers that use high quality soft-clipped reads and base mismatches to perform local realignment consistently had the highest precision and recall in identifying insertions and deletions ranging from1 to 50 bp. The remaining variant callers had lower recall values associated with identification of insertions greater than 20 bp.

Clinicopathologic and genomic features of lobular like invasive mammary carcinoma: is it a distinct entity?

This study describes "lobular-like invasive mammary carcinomas" (LLIMCas), a group of low- to intermediate-grade invasive mammary carcinomas with discohesive, diffusely infiltrative cells showing retained circumferential membranous immunoreactivity for both E-cadherin and p120. We analyzed the clinical-pathologic features of 166 LLIMCas compared to 104 classical invasive lobular carcinomas (ILCs) and 100 grade 1 and 2 invasive ductal carcinomas (IDCs). Tumor size and pT stage of LLIMCas were intermediate between IDCs and ILCs, and yet often underestimated on imaging and showed frequent positive margins on the first resection. Despite histomorphologic similarities to classical ILC, the discohesion in LLIMCa was independent of E-cadherin/p120 immunophenotypic alteration. An exploratory, hypothesis-generating analysis of the genomic features of 14 randomly selected LLIMCas and classical ILCs (7 from each category) was performed utilizing an FDA-authorized targeted capture sequencing assay (MSK-IMPACT). None of the seven LLIMCas harbored CDH1 loss-of-function mutations, and none of the CDH1 alterations detected in two of the LLIMCas was pathogenic. In contrast, all seven ILCs harbored CDH1 loss-of-function mutations coupled with the loss of heterozygosity of the CDH1 wild-type allele. Four of the six evaluable LLIMCas were positive for CDH1 promoter methylation, which may partially explain the single-cell infiltrative morphology seen in LLIMCa. Further studies are warranted to better define the molecular basis of the discohesive cellular morphology in LLIMCa. Until more data becomes available, identifying LLIMCas and distinguishing them from typical IDCs and ILCs would be justified. In patients with LLIMCas, preoperative MRI should be entertained to guide surgical management.

OpenPBTA: The Open Pediatric Brain Tumor Atlas

Pediatric brain and spinal cancers are collectively the leading disease-related cause of death in children; thus, we urgently need curative therapeutic strategies for these tumors. To accelerate such discoveries, the Children's Brain Tumor Network (CBTN) and Pacific Pediatric Neuro-Oncology Consortium (PNOC) created a systematic process for tumor biobanking, model generation, and sequencing with immediate access to harmonized data. We leverage these data to establish OpenPBTA, an open collaborative project with over 40 scalable analysis modules that genomically characterize 1,074 pediatric brain tumors. Transcriptomic classification reveals universal TP53 dysregulation in mismatch repair-deficient hypermutant high-grade gliomas and TP53 loss as a significant marker for poor overall survival in ependymomas and H3 K28-mutant diffuse midline gliomas. Already being actively applied to other pediatric cancers and PNOC molecular tumor board decision-making, OpenPBTA is an invaluable resource to the pediatric oncology community.

NR4A3 Expression Is Consistently Absent in Acinic Cell Carcinomas of the Breast: A Potential Nosologic Shift

Acinic cell carcinoma (AciCC) is a tumor that is recognized in both the breast and salivary glands. Recently, the recurrent genomic rearrangement, t(4;9)(q13;q31) was identified in salivary AciCC that results in constitutive upregulation of the nuclear transcription factor NR4A3, which can be detected by immunohistochemistry. In this study, we sought to evaluate NR4A3 expression in breast AciCC using immunohistochemistry. Strong and diffuse nuclear staining was considered a positive result. Sixteen AciCCs were studied, including 8 pure AciCCs and 8 AciCCs admixed with other types (invasive carcinoma of no special type in 5 cases and metaplastic carcinoma in 3 cases). All 16 AciCCs showed negative results for NR4A3 expression. Four cases with available material were evaluated for rearrangements of the NR4A3 gene by fluorescence in situ hybridization and no rearrangements were observed. Whole-genome sequencing of 1 AciCC revealed a TP53 splice-site mutation, high levels of genomic instability, and genomic features of homologous recombination DNA repair defects; a structural variant analysis of this case did not reveal the presence of a t(4;9) rearrangement. We conclude that breast AciCCs consistently lack NR4A3 rearrangement or overexpression, unlike most of the salivary AciCCs, and that consistent with previous results, breast AciCCs are associated with genomic alterations more similar to those seen in triple-negative breast carcinomas than salivary gland AciCCs. These results suggest that unlike other salivary gland-like tumors that occur in the breast, the molecular underpinnings of the salivary gland and breast AciCCs are different and that the salivary gland and breast AciCCs likely represent distinct entities.

satmut_utils: a simulation and variant calling package for multiplexed assays of variant effect

The impact of millions of individual genetic variants on molecular phenotypes in coding sequences remains unknown. Multiplexed assays of variant effect (MAVEs) are scalable methods to annotate relevant variants, but existing software lacks standardization, requires cumbersome configuration, and does not scale to large targets. We present satmut_utils as a flexible solution for simulation and variant quantification. We then benchmark MAVE software using simulated and real MAVE data. We finally determine mRNA abundance for thousands of cystathionine beta-synthase variants using two experimental methods. The satmut_utils package enables high-performance analysis of MAVEs and reveals the capability of variants to alter mRNA abundance.

The motif composition of variable number tandem repeats impacts gene expression

Understanding the impact of DNA variation on human traits is a fundamental question in human genetics. Variable number tandem repeats (VNTRs) make up ∼3% of the human genome but are often excluded from association analysis owing to poor read mappability or divergent repeat content. Although methods exist to estimate VNTR length from short-read data, it is known that VNTRs vary in both length and repeat (motif) composition. Here, we use a repeat-pangenome graph (RPGG) constructed on 35 haplotype-resolved assemblies to detect variation in both VNTR length and repeat composition. We align population-scale data from the Genotype-Tissue Expression (GTEx) Consortium to examine how variations in sequence composition may be linked to expression, including cases independent of overall VNTR length. We find that 9422 out of 39,125 VNTRs are associated with nearby gene expression through motif variations, of which only 23.4% are accessible from length. Fine-mapping identifies 174 genes to be likely driven by variation in certain VNTR motifs and not overall length. We highlight two genes, CACNA1C and RNF213, that have expression associated with motif variation, showing the utility of RPGG analysis as a new approach for trait association in multiallelic and highly variable loci.

Challenges in neoantigen-directed therapeutics

A fundamental prerequisite for the efficacy of cancer immunotherapy is the presence of functional, antigen-specific T cells within the tumor. Neoantigen-directed therapy is a promising strategy that aims at targeting the host's immune response against tumor-specific antigens, thereby eradicating cancer cells. Initial forays have been made in clinical environments utilizing vaccines and adoptive cell therapy; however, many challenges lie ahead. We provide an in-depth overview of the current state of the field with an emphasis on in silico neoantigen discovery and the clinical aspects that need to be addressed to unlock the full potential of this therapy.

Mutational Landscape of Bladder Cancer in Mexican Patients: KMT2D Mutations and chr11q15.5 Amplifications Are Associated with Muscle Invasion

Bladder cancer (BC) is the most common neoplasm of the urinary tract, which originates in the epithelium that covers the inner surface of the bladder. The molecular BC profile has led to the development of different classifications of non-muscle invasive bladder cancer (NMIBC) and muscle-invasive bladder cancer (MIBC). However, the genomic BC landscape profile of the Mexican population, including NMIBC and MIBC, is unknown. In this study, we aimed to identify somatic single nucleotide variants (SNVs) and copy number variations (CNVs) in Mexican patients with BC and their associations with clinical and pathological characteristics. We retrospectively evaluated 37 patients treated between 2012 and 2021 at the National Cancer Institute-Mexico (INCan). DNA samples were obtained from paraffin-embedded tumor tissues and exome sequenced. Strelka2 and Lancet packages were used to identify SNVs and insertions or deletions. FACETS was used to determine CNVs. We found a high frequency of mutations in TP53 and KMT2D, gains in 11q15.5 and 19p13.11-q12, and losses in 7q11.23. STAG2 mutations and 1q11.23 deletions were also associated with NMIBC and low histologic grade.

Comprehensive Genomic Profiling of Cell-Free Circulating Tumor DNA Detects Response to Ribociclib Plus Letrozole in a Patient with Metastatic Breast Cancer

Analysis of cell-free circulating tumor DNA obtained by liquid biopsy is a non-invasive approach that may provide clinically actionable information when conventional tissue biopsy is inaccessible or infeasible. Here, we followed a patient with hormone receptor-positive and human epidermal growth factor receptor (HER) 2-negative breast cancer who developed bone metastases seven years after mastectomy. We analyzed circulating cell-free DNA (cfDNA) extracted from plasma using high-depth massively parallel sequencing targeting 468 cancer-associated genes, and we identified a clonal hotspot missense mutation in the PIK3CA gene (3:178952085, A > G, H1047R) and amplification of the CCND1 gene. Whole-exome sequencing revealed that both alterations were present in the primary tumor. After treatment with ribociclib plus letrozole, the genetic abnormalities were no longer detected in cfDNA. These results underscore the clinical utility of combining liquid biopsy and comprehensive genomic profiling to monitor treatment response in patients with metastasized breast cancer.

Triple MAPK inhibition salvaged a relapsed post-BCMA CAR-T cell therapy multiple myeloma patient with a BRAF V600E subclonal mutation

Background

Multiple Myeloma (MM) is a progressive plasma cell neoplasm characterized by heterogeneous clonal expansion. Despite promising response rates achieved with anti-BCMA CAR-T cell therapy, patients may still relapse and there are currently no clear therapeutic options in post-CAR-T settings. In this report, we present a case of a post-BCMA CAR-T relapsed/refractory (RR) MM patient with skin extramedullary disease (EMD) in which a novel MAPK inhibition combinatorial strategy was implemented based on next-generation sequencing and in vitro experiments.

Case presentation

A 61-year-old male with penta-refractory MM penta- (IgA lambda), ISS stage 3 with hyperdiploidy, gain of 1q21 and del13 was treated with anti-BCMA CAR-T cell therapy, achieving a best response of VGPR. He progressed after 6 months and was salvaged for a short period with autologous stem cell transplantation. Eventually, he progressed with extramedullary disease manifested as subcutaneous nodules. Based on whole-exome sequencing, we identified a BRAF (V600E) dominant subclone in both bone marrow and cutaneous plasmacytoma. Following in vitro experiments, and according to our previous studies, we implemented a triple MAPK inhibition strategy under which the patient achieved a very good partial response for 110 days, which allowed to bridge him to subsequent clinical trials and eventually achieve a stringent complete response (sCR).

Conclusion

Here, we show the applicability, effectiveness, and tolerability the triple MAPK inhibition strategy in the context of post-BCMA CAR-T failure in specific subset of patients. The triple therapy could bridge our hospice bound RRMM patient with BRAF (V600E) to further therapeutic options where sCR was achieved. We will further evaluate triple MAPK inhibition in patients with BRAF V600E in a precision medicine clinical trial launching soon.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13045-022-01330-3.

SETD2 Haploinsufficiency Enhances Germinal Center-Associated AICDA Somatic Hypermutation to Drive B-cell Lymphomagenesis

SETD2 is the sole histone methyltransferase responsible for H3K36me3, with roles in splicing, transcription initiation, and DNA damage response. Homozygous disruption of SETD2 yields a tumor suppressor effect in various cancers. However, SETD2 mutation is typically heterozygous in diffuse large B-cell lymphomas. Here we show that heterozygous Setd2 deficiency results in germinal center (GC) hyperplasia and increased competitive fitness, with reduced DNA damage checkpoint activity and apoptosis, resulting in accelerated lymphomagenesis. Impaired DNA damage sensing in Setd2-haploinsufficient germinal center B (GCB) and lymphoma cells associated with increased AICDA-induced somatic hypermutation, complex structural variants, and increased translocations including those activating MYC. DNA damage was selectively increased on the nontemplate strand, and H3K36me3 loss was associated with greater RNAPII processivity and mutational burden, suggesting that SETD2-mediated H3K36me3 is required for proper sensing of cytosine deamination. Hence, Setd2 haploinsufficiency delineates a novel GCB context-specific oncogenic pathway involving defective epigenetic surveillance of AICDA-mediated effects on transcribed genes.

SIGNIFICANCE

Our findings define a B cell-specific oncogenic effect of SETD2 heterozygous mutation, which unleashes AICDA mutagenesis of nontemplate strand DNA in the GC reaction, resulting in lymphomas with heavy mutational burden. GC-derived lymphomas did not tolerate SETD2 homozygous deletion, pointing to a novel context-specific therapeutic vulnerability. This article is highlighted in the In This Issue feature, p. 1599.

ESR1 hotspot mutations in endometrial stromal sarcoma with high-grade transformation and endocrine treatment

High-grade endometrial stromal sarcomas (HGESSs) are more aggressive and have higher rates of resistance to endocrine therapy than low-grade endometrial stromal sarcomas (LGESSs). The pathogenesis of hormonal resistance in these lesions has yet to be defined. Here we sought to histologically and genetically characterize 3 LGESSs and their recurrences that underwent histologic high-grade transformation following endocrine therapy. For this, DNA from primary tumors and select subsequent recurrences were subject to massively parallel sequencing targeting 468 cancer-related genes. Somatic mutation analyses were performed using validated bioinformatics methods. In addition, RNA from each case was evaluated for the presence of gene fusions using targeted RNA-sequencing. All patients initially presented with LGESS, developed HGESS recurrences, and received at least 2 lines of hormonal suppressive therapy. Gene fusions classically described as associated with LGESS were identified in all 3 cases, including JAZF1-PHF1, EPC1-PHF1 and JAZF1-SUZ12 fusions for Cases 1, 2 and 3, respectively. Targeted sequencing analysis revealed that none of the primary LGESS, however the HGESS recurrences of Cases 1 and 3, and the LGESS and HGESS recurrences of Case 2 post endocrine treatment harbored ESR1 p.Y537S hotspot mutations. These ESR1 ligand-binding domain mutations have been found as a mechanism of acquired endocrine resistance in breast cancer. Also, a reduction in estrogen receptor (ER) expression was observed in recurrences. Our findings suggest that the ESR1 p.Y537S hotspot mutation in LGESS with histologic high-grade transformation may be associated with endocrine resistance in these lesions. Furthermore, our data suggest that genetic analyses may be performed in recurrent LGESS following hormonal therapy, development of high-grade morphology, and/or altered/diminished ER expression.

Loss of MMR and TGFBR2 Increases the Susceptibility to Microbiota-Dependent Inflammation-Associated Colon Cancer

BACKGROUND AND AIMS

Mutations in DNA mismatch repair (MMR) genes are causative in Lynch syndrome and a significant proportion of sporadic colorectal cancers (CRCs). MMR-deficient (dMMR) CRCs display increased mutation rates, with mutations frequently accumulating at short repetitive DNA sequences throughout the genome (microsatellite instability). The TGFBR2 gene is one of the most frequently mutated genes in dMMR CRCs. Therefore, we generated an animal model to study how the loss of both TGFBR2 signaling impacts dMMR-driven intestinal tumorigenesis in vivo and explore the impact of the gut microbiota.

METHODS

We generated VCMsh2/Tgfbr2 mice in which Msh2loxP and Tgfbr2loxP alleles are inactivated by Villin-Cre recombinase in the intestinal epithelium. VCMsh2/Tgfbr2 mice were analyzed for their rate of intestinal cancer development and for the mutational spectra and gene expression profiles of tumors. In addition, we assessed the impact of chemically induced chronic inflammation and gut microbiota composition on colorectal tumorigenesis.

RESULTS

VCMsh2/Tgfbr2 mice developed small intestinal adenocarcinomas and CRCs with histopathological features highly similar to CRCs in Lynch syndrome patients. The CRCs in VCMsh2/Tgfbr2 mice were associated with the presence of colitis and displayed genetic and histological features that resembled inflammation-associated CRCs in human patients. The development of CRCs in VCMsh2/Tgfbr2 mice was strongly modulated by the gut microbiota composition, which in turn was impacted by the TGFBR2 status of the tumors.

CONCLUSIONS

Our results demonstrate a synergistic interaction between MMR and TGFBR2 inactivation in inflammation-associated colon tumorigenesis and highlight the crucial impact of the gut microbiota on modulating the incidence of inflammation-associated CRCs.

RAD51B Harbors Germline Mutations Associated With Pancreatic Ductal Adenocarcinoma

PURPOSE

Genetic alterations in many components of the homologous recombination, DNA damage response, and repair (HR-DDR) pathway are involved in the hereditary cancer syndromes, including familial pancreatic cancer. HR-DDR genes beyond BRCA1, BRCA2, ATM, and PALB2 may also mutate and confer the HR-DDR deficiency in pancreatic ductal adenocarcinoma (PDAC).

METHODS

We conducted a study to examine the genetic alterations using a companion diagnostic 15-gene HR-DDR panel in PDACs. HR-DDR gene mutations were identified and characterized by whole-exome sequencing and whole-genome sequencing. Different HR-DDR gene mutations are associated with variable homologous recombination deficiency (HRD) scores.

RESULTS

Eight of 50 PDACs with at least one HR-DDR gene mutation were identified. One tumor with BRCA2 mutations is associated with a high HRD score. However, another tumor with a CHEK2 mutation is associated with a zero HRD score. Notably, four of eight PDACs in this study harbor a RAD51B gene mutation. All four RAD51B gene mutations were germline mutations. However, currently, RAD51B is not the gene panel for germline tests.

CONCLUSION

The finding in this study thus supports including RAD51B in the germline test of HR-DDR pathway genes.

Shedding light on the hidden human proteome expands immunopeptidome in cancer

Unrestrained cellular growth and immune escape of a tumor are associated with the incidental errors of the genome and transcriptome. Advances in next-generation sequencing have identified thousands of genomic and transcriptomic aberrations that generate variant peptides that assemble the hidden proteome, further expanding the immunopeptidome. Emerging next-generation sequencing technologies and a number of computational methods estimated the abundance of immune infiltration from bulk transcriptome have advanced our understanding of tumor microenvironments. Here, we will characterize several major types of tumor-specific antigens arising from single-nucleotide variants, insertions and deletions, gene fusion, alternative splicing, RNA editing and non-coding RNAs. Finally, we summarize the current state-of-the-art computational and experimental approaches or resources and provide an integrative pipeline for the identification of candidate tumor antigens. Together, the systematic investigation of the hidden proteome in cancer will help facilitate the development of effective and durable immunotherapy targets for cancer.

Achieving robust somatic mutation detection with deep learning models derived from reference data sets of a cancer sample

BACKGROUND

Accurate detection of somatic mutations is challenging but critical in understanding cancer formation, progression, and treatment. We recently proposed NeuSomatic, the first deep convolutional neural network-based somatic mutation detection approach, and demonstrated performance advantages on in silico data.

RESULTS

In this study, we use the first comprehensive and well-characterized somatic reference data sets from the SEQC2 consortium to investigate best practices for using a deep learning framework in cancer mutation detection. Using the high-confidence somatic mutations established for a cancer cell line by the consortium, we identify the best strategy for building robust models on multiple data sets derived from samples representing real scenarios, for example, a model trained on a combination of real and spike-in mutations had the highest average performance.

CONCLUSIONS

The strategy identified in our study achieved high robustness across multiple sequencing technologies for fresh and FFPE DNA input, varying tumor/normal purities, and different coverages, with significant superiority over conventional detection approaches in general, as well as in challenging situations such as low coverage, low variant allele frequency, DNA damage, and difficult genomic regions.

Detection of minor variants in Mycobacterium tuberculosis whole genome sequencing data

The study of genetic minority variants is fundamental to the understanding of complex processes such as evolution, fitness, transmission, virulence, heteroresistance and drug tolerance in Mycobacterium tuberculosis (Mtb). We evaluated the performance of the variant calling tool LoFreq to detect de novo as well as drug resistance conferring minor variants in both in silico and clinical Mtb next generation sequencing (NGS) data. The in silico simulations demonstrated that LoFreq is a conservative variant caller with very high precision (≥96.7%) over the entire range of depth of coverage tested (30x to1000x), independent of the type and frequency of the minor variant. Sensitivity increased with increasing depth of coverage and increasing frequency of the variant, and was higher for calling insertion and deletion (indel) variants than for single nucleotide polymorphisms (SNP). The variant frequency limit of detection was 0.5% and 3% for indel and SNP minor variants, respectively. For serial isolates from a patient with DR-TB; LoFreq successfully identified all minor Mtb variants in the Rv0678 gene (allele frequency as low as 3.22% according to targeted deep sequencing) in whole genome sequencing data (median coverage of 62X). In conclusion, LoFreq can successfully detect minor variant populations in Mtb NGS data, thus limiting the need for filtering of possible false positive variants due to sequencing error. The observed performance statistics can be used to determine the limit of detection in existing whole genome sequencing Mtb data and guide the required depth of future studies that aim to investigate the presence of minor variants.

The histologic phenotype of lung cancers is associated with transcriptomic features rather than genomic characteristics

Histology plays an essential role in therapeutic decision-making for lung cancer patients. However, the molecular determinants of lung cancer histology are largely unknown. We conduct whole-exome sequencing and microarray profiling on 19 micro-dissected tumor regions of different histologic subtypes from 9 patients with lung cancers of mixed histology. A median of 68.9% of point mutations and 83% of copy number aberrations are shared between different histologic components within the same tumors. Furthermore, different histologic components within the tumors demonstrate similar subclonal architecture. On the other hand, transcriptomic profiling reveals shared pathways between the same histologic subtypes from different patients, which is supported by the analyses of the transcriptomic data from 141 cell lines and 343 lung cancers of different histologic subtypes. These data derived from mixed histologic subtypes in the setting of identical genetic background and exposure history support that the histologic fate of lung cancer cells is associated with transcriptomic features rather than the genomic profiles in most tumors.

The molecular determinants of lung cancer histologic subtypes are not well understood. Here the authors analyze lung cancers of mixed histology and find that histologic subtypes are associated with transcriptomic features rather than genomic profiles in most tumors.

Calling small variants using universality with Bayes-factor-adjusted odds ratios

The application of next-generation sequencing in research and particularly in clinical routine requires highly accurate variant calling. Here we describe UVC, a method for calling small variants of germline or somatic origin. By unifying opposite assumptions with sublation, we discovered the following two empirical laws to improve variant calling: allele fraction at high sequencing depth is inversely proportional to the cubic root of variant-calling error rate, and odds ratios adjusted with Bayes factors can model various sequencing biases. UVC outperformed other variant callers on the GIAB germline truth sets, 192 scenarios of in silico mixtures simulating 192 combinations of tumor/normal sequencing depths and tumor/normal purities, the GIAB somatic truth sets derived from physical mixture, and the SEQC2 somatic reference sets derived from the breast-cancer cell-line HCC1395. UVC achieved 100% concordance with the manual review conducted by multiple independent researchers on a Qiagen 71-gene-panel dataset derived from 16 patients with colon adenoma. UVC outperformed other unique molecular identifier (UMI)-aware variant callers on the datasets used for publishing these variant callers. Performance was measured with sensitivity-specificity trade off for called variants. The improved variant calls generated by UVC from previously published UMI-based sequencing data provided additional insight about DNA damage repair. UVC is open-sourced under the BSD 3-Clause license at https://github.com/genetronhealth/uvc and quay.io/genetronhealth/gcc-6-3-0-uvc-0-6-0-441a694.

Whole genome and exome sequencing reference datasets from a multi-center and cross-platform benchmark study

With the rapid advancement of sequencing technologies, next generation sequencing (NGS) analysis has been widely applied in cancer genomics research. More recently, NGS has been adopted in clinical oncology to advance personalized medicine. Clinical applications of precision oncology require accurate tests that can distinguish tumor-specific mutations from artifacts introduced during NGS processes or data analysis. Therefore, there is an urgent need to develop best practices in cancer mutation detection using NGS and the need for standard reference data sets for systematically measuring accuracy and reproducibility across platforms and methods. Within the SEQC2 consortium context, we established paired tumor-normal reference samples and generated whole-genome (WGS) and whole-exome sequencing (WES) data using sixteen library protocols, seven sequencing platforms at six different centers. We systematically interrogated somatic mutations in the reference samples to identify factors affecting detection reproducibility and accuracy in cancer genomes. These large cross-platform/site WGS and WES datasets using well-characterized reference samples will represent a powerful resource for benchmarking NGS technologies, bioinformatics pipelines, and for the cancer genomics studies.

Benchmarking pipelines for subclonal deconvolution of bulk tumour sequencing data

Intratumour heterogeneity provides tumours with the ability to adapt and acquire treatment resistance. The development of more effective and personalised treatments for cancers, therefore, requires accurate characterisation of the clonal architecture of tumours, enabling evolutionary dynamics to be tracked. Many methods exist for achieving this from bulk tumour sequencing data, involving identifying mutations and performing subclonal deconvolution, but there is a lack of systematic benchmarking to inform researchers on which are most accurate, and how dataset characteristics impact performance. To address this, we use the most comprehensive tumour genome simulation tool available for such purposes to create 80 bulk tumour whole exome sequencing datasets of differing depths, tumour complexities, and purities, and use these to benchmark subclonal deconvolution pipelines. We conclude that i) tumour complexity does not impact accuracy, ii) increasing either purity or purity-corrected sequencing depth improves accuracy, and iii) the optimal pipeline consists of Mutect2, FACETS and PyClone-VI. We have made our benchmarking datasets publicly available for future use.

Morphological and genomic characteristics of breast cancers occurring in individuals with Lynch Syndrome

PURPOSE

Lynch syndrome (LS) is defined by germline pathogenic mutations involving DNA Mismatch Repair (MMR) genes and linked with the development of MMR-deficient (MMRd) colon and endometrial cancers. Whether breast cancers (BC) developing in context of LS are causally related to MMR deficiency (MMRd), remains controversial. Thus, we explored the morphological and genomic characteristics of BCs occurring in LS individuals.

EXPERIMENTAL DESIGN

A retrospective analysis of 20,110 cancer patients who underwent multigene panel genetic testing was performed to identify individuals with a likely pathogenic/pathogenic germline variant in MLH1, MSH2, MSH6 or PMS2 who developed BCs. The histological characteristics and immunohistochemical (IHC) assessment of BCs for MMR proteins and programmed death-ligand 1 (PD-L1) expression were assessed on cases with available materials. DNA samples from paired tumors and blood were sequenced with MSK-IMPACT ({greater than or equal to}468 key cancer genes). MSI status was assessed utilizing MSISensor. Mutational signatures were defined using SigMA.

RESULTS

272 LS individuals were identified, 13 (5%) of whom had primary BCs. The majority of BCs (92%) were hormone receptor positive tumors. Five (42%) of 12 BCs displayed loss of MMR proteins by IHC. Four (36%) of 11 BCs subjected to tumor-normal sequencing showed dominant microsatellite instability mutational signatures, high tumor mutational burden and indeterminate (27%) or high MSISensor scores (9%). One patient with metastatic MMRd BC received anti-PD1 therapy and achieved a robust and durable response.

CONCLUSIONS

A subset of BCs developing in LS individuals are etiologically linked to MMRd and may benefit from anti-PD1/PD-L1 immunotherapy.

Recurrence biomarkers of triple negative breast cancer treated with neoadjuvant chemotherapy and anti-EGFR antibodies

To find metastatic recurrence biomarkers of triple-negative breast cancer (TNBC) treated by neoadjuvant chemotherapy and anti-EGFR antibodies (NAT), we evaluated tumor genomic, transcriptomic, and immune features, using MSK-IMPACT assay, gene arrays, Nanostring technology, and TIL assessment on H&E. Six patients experienced a rapid fatal recurrence (RR) and other 6 had later non-fatal recurrences (LR). Before NAT, RR had low expression of 6 MHC class I and 13 MHC class II genes but were enriched in upregulated genes involved in the cell cycle-related pathways. Their TIL number before NAT in RR was very low (<5%) and did not increase after treatment. In post-NAT residual tumors, RR cases showed high expression of SOX2 and CXCR4. Our results indicate that high expression of cell cycle genes, combined with cold immunological phenotype, may predict strong TNBC resistance to NAT and rapid progression after it. This biomarker combination is worth validation in larger studies.

Toward best practice in cancer mutation detection with whole-genome and whole-exome sequencing

Clinical applications of precision oncology require accurate tests that can distinguish true cancer-specific mutations from errors introduced at each step of next-generation sequencing (NGS). To date, no bulk sequencing study has addressed the effects of cross-site reproducibility, nor the biological, technical and computational factors that influence variant identification. Here we report a systematic interrogation of somatic mutations in paired tumor-normal cell lines to identify factors affecting detection reproducibility and accuracy at six different centers. Using whole-genome sequencing (WGS) and whole-exome sequencing (WES), we evaluated the reproducibility of different sample types with varying input amount and tumor purity, and multiple library construction protocols, followed by processing with nine bioinformatics pipelines. We found that read coverage and callers affected both WGS and WES reproducibility, but WES performance was influenced by insert fragment size, genomic copy content and the global imbalance score (GIV; G > T/C > A). Finally, taking into account library preparation protocol, tumor content, read coverage and bioinformatics processes concomitantly, we recommend actionable practices to improve the reproducibility and accuracy of NGS experiments for cancer mutation detection.

Mesonephric and mesonephric-like carcinomas of the female genital tract: molecular characterization including cases with mixed histology and matched metastases

Mesonephric carcinoma of the cervix is a rare tumor derived from Wolffian remnants. Mesonephric-like carcinomas of the ovary and endometrium, while morphologically similar, do not have obvious Wolffian derivation. Here, we sought to characterize the repertoire of genetic alterations in primary mesonephric and mesonephric-like carcinomas, in the distinct histologic components of mixed cases, as well as in matched primary tumors and metastases. DNA from microdissected tumor and normal tissue from mesonephric carcinomas (cervix, n = 8) and mesonephric-like carcinomas (ovarian n = 15, endometrial n = 13) were subjected to sequencing targeting 468 cancer-related genes. The histologically distinct components of four cases with mixed histology and four primary tumors and their matched metastases were microdissected and analyzed separately. Mesonephric-like carcinomas were underpinned by somatic KRAS mutations (25/28, 89%) akin to mesonephric carcinomas (8/8, 100%), but also harbored genetic alterations more frequently reported in Müllerian tumors. Mesonephric-like carcinomas that lacked KRAS mutations harbored NRAS (n = 2, ovary) or BRAF (n = 1, endometrium) hotspot mutations. PIK3CA mutations were identified in both mesonephric-like (8/28, 28%) and mesonephric carcinomas (2/8, 25%). Only mesonephric-like tumors harbored CTNNB1 hotspot (4/28, 14%) and PTEN (3/13, 23%) mutations. Copy number analysis revealed frequent gains of chromosomes 1q and 10 in both mesonephric (87% 1q; 50% chromosome 10) and mesonephric-like tumors (89% 1q; 43% chromosome 10). Chromosome 12 gains were more frequent in ovarian mesonephric-like carcinomas, and losses of chromosome 9 were more frequent in mesonephric than in mesonephric-like carcinomas (both p = 0.01, Fisher's exact test). The histologically distinct components of four mixed cases were molecularly related and shared similar patterns of genetic alterations. The progression from primary to metastatic lesions involved the acquisition of additional mutations, and/or shifts from subclonal to clonal mutations. Our findings suggest that mesonephric-like carcinomas are derived from a Müllerian substrate with differentiation along Wolffian/mesonephric lines.

Somatic variant analysis of linked-reads sequencing data with Lancet

SUMMARY

We present a new version of the popular somatic variant caller, Lancet, that supports the analysis of linked-reads sequencing data. By seamlessly integrating barcodes and haplotype read assignments within the colored De Bruijn graph local-assembly framework, Lancet computes a barcode-aware coverage and identifies variants that disagree with the local haplotype structure.

AVAILABILITY AND IMPLEMENTATION

Lancet is implemented in C++ and available for academic and non-commercial research purposes as an open-source package at https://github.com/nygenome/lancet.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Targeting Neoepitopes to Treat Solid Malignancies: Immunosurgery

Successful outcome of immune checkpoint blockade in patients with solid cancers is in part associated with a high tumor mutational burden (TMB) and the recognition of private neoantigens by T-cells. The quality and quantity of target recognition is determined by the repertoire of ‘neoepitope’-specific T-cell receptors (TCRs) in tumor-infiltrating lymphocytes (TIL), or peripheral T-cells. Interferon gamma (IFN-γ), produced by T-cells and other immune cells, is essential for controlling proliferation of transformed cells, induction of apoptosis and enhancing human leukocyte antigen (HLA) expression, thereby increasing immunogenicity of cancer cells. TCR αβ-dependent therapies should account for tumor heterogeneity and availability of the TCR repertoire capable of reacting to neoepitopes and functional HLA pathways. Immunogenic epitopes in the tumor-stroma may also be targeted to achieve tumor-containment by changing the immune-contexture in the tumor microenvironment (TME). Non protein-coding regions of the tumor-cell genome may also contain many aberrantly expressed, non-mutated tumor-associated antigens (TAAs) capable of eliciting productive anti-tumor immune responses. Whole-exome sequencing (WES) and/or RNA sequencing (RNA-Seq) of cancer tissue, combined with several layers of bioinformatic analysis is commonly used to predict possible neoepitopes present in clinical samples. At the ImmunoSurgery Unit of the Champalimaud Centre for the Unknown (CCU), a pipeline combining several tools is used for predicting private mutations from WES and RNA-Seq data followed by the construction of synthetic peptides tailored for immunological response assessment reflecting the patient’s tumor mutations, guided by MHC typing. Subsequent immunoassays allow the detection of differential IFN-γ production patterns associated with (intra-tumoral) spatiotemporal differences in TIL or peripheral T-cells versus TIL. These bioinformatics tools, in addition to histopathological assessment, immunological readouts from functional bioassays and deep T-cell ‘adaptome’ analyses, are expected to advance discovery and development of next-generation personalized precision medicine strategies to improve clinical outcomes in cancer in the context of i) anti-tumor vaccination strategies, ii) gauging mutation-reactive T-cell responses in biological therapies and iii) expansion of tumor-reactive T-cells for the cellular treatment of patients with cancer.

A unified haplotype-based method for accurate and comprehensive variant calling

Almost all haplotype-based variant callers were designed specifically for detecting common germline variation in diploid populations, and give suboptimal results in other scenarios. Here we present Octopus, a variant caller that uses a polymorphic Bayesian genotyping model capable of modeling sequencing data from a range of experimental designs within a unified haplotype-aware framework. Octopus combines sequencing reads and prior information to phase-called genotypes of arbitrary ploidy, including those with somatic mutations. We show that Octopus accurately calls germline variants in individuals, including single nucleotide variants, indels and small complex replacements such as microinversions. Using a synthetic tumor data set derived from clean sequencing data from a sample with known germline haplotypes and observed mutations in a large cohort of tumor samples, we show that Octopus is more sensitive to low-frequency somatic variation, yet calls considerably fewer false positives than other methods. Octopus also outputs realigned evidence BAM files to aid validation and interpretation.

2-kupl: mapping-free variant detection from DNA-seq data of matched samples

BACKGROUND

The detection of genome variants, including point mutations, indels and structural variants, is a fundamental and challenging computational problem. We address here the problem of variant detection between two deep-sequencing (DNA-seq) samples, such as two human samples from an individual patient, or two samples from distinct bacterial strains. The preferred strategy in such a case is to align each sample to a common reference genome, collect all variants and compare these variants between samples. Such mapping-based protocols have several limitations. DNA sequences with large indels, aggregated mutations and structural variants are hard to map to the reference. Furthermore, DNA sequences cannot be mapped reliably to genomic low complexity regions and repeats.

RESULTS

We introduce 2-kupl, a k-mer based, mapping-free protocol to detect variants between two DNA-seq samples. On simulated and actual data, 2-kupl achieves higher accuracy than other mapping-free protocols. Applying 2-kupl to prostate cancer whole exome sequencing data, we identify a number of candidate variants in hard-to-map regions and propose potential novel recurrent variants in this disease.

CONCLUSIONS

We developed a mapping-free protocol for variant calling between matched DNA-seq samples. Our protocol is suitable for variant detection in unmappable genome regions or in the absence of a reference genome.

Massively parallel sequencing analysis of 68 gastric-type cervical adenocarcinomas reveals mutations in cell cycle-related genes and potentially targetable mutations

Gastric-type cervical adenocarcinoma (GCA) is an aggressive type of endocervical adenocarcinoma characterized by mucinous morphology, gastric-type mucin, lack of association with human papillomavirus (HPV) and resistance to chemo/radiotherapy. We characterized the landscape of genetic alterations in a large cohort of GCAs, and compared it with that of usual-type HPV-associated endocervical adenocarcinomas (UEAs), pancreatic adenocarcinomas (PAs) and intestinal-type gastric adenocarcinomas (IGAs). GCAs (n = 68) were subjected to massively parallel sequencing targeting 410-468 cancer-related genes. Somatic mutations and copy number alterations (CNAs) were determined using validated bioinformatics methods. Mutational data for UEAs (n = 21), PAs (n = 178), and IGAs (n = 148) from The Cancer Genome Atlas (TCGA) were obtained from cBioPortal. GCAs most frequently harbored somatic mutations in TP53 (41%), CDKN2A (18%), KRAS (18%), and STK11 (10%). Potentially targetable mutations were identified in ERBB3 (10%), ERBB2 (8%), and BRAF (4%). GCAs displayed low levels of CNAs with no recurrent amplifications or homozygous deletions. In contrast to UEAs, GCAs harbored more frequent mutations affecting cell cycle-related genes including TP53 (41% vs 5%, p < 0.01) and CDKN2A (18% vs 0%, p = 0.01), and fewer PIK3CA mutations (7% vs 33%, p = 0.01). TP53 mutations were less prevalent in GCAs compared to PAs (41% vs 56%, p < 0.05) and IGAs (41% vs 57%, p < 0.05). GCAs showed a higher frequency of STK11 mutations than PAs (10% vs 2%, p < 0.05) and IGAs (10% vs 1%, p < 0.05). GCAs harbored more frequent mutations in ERBB2 and ERBB3 (9% vs 1%, and 10% vs 0.5%, both p < 0.01) compared to PAs, and in CDKN2A (18% vs 1%, p < 0.05) and KRAS (18% vs 6%, p < 0.05) compared to IGAs. GCAs harbor recurrent somatic mutations in cell cycle-related genes and in potentially targetable genes, including ERBB2/3. Mutations in genes such as STK11 may be used as supportive evidence to help distinguish GCAs from other adenocarcinomas with similar morphology in metastatic sites.

Molecular landscape and subtype-specific therapeutic response of nasopharyngeal carcinoma revealed by integrative pharmacogenomics

Nasopharyngeal carcinoma (NPC) is a malignant head and neck cancer type with high morbidity in Southeast Asia, however the pathogenic mechanism of this disease is poorly understood. Using integrative pharmacogenomics, we find that NPC subtypes maintain distinct molecular features, drug responsiveness, and graded radiation sensitivity. The epithelial carcinoma (EC) subtype is characterized by activations of microtubule polymerization and defective mitotic spindle checkpoint related genes, whereas sarcomatoid carcinoma (SC) and mixed sarcomatoid-epithelial carcinoma (MSEC) subtypes exhibit enriched epithelial-mesenchymal transition (EMT) and invasion promoting genes, which are well correlated with their morphological features. Furthermore, patient-derived organoid (PDO)-based drug test identifies potential subtype-specific treatment regimens, in that SC and MSEC subtypes are sensitive to microtubule inhibitors, whereas EC subtype is more responsive to EGFR inhibitors, which is synergistically enhanced by combining with radiotherapy. Through combinational chemoradiotherapy (CRT) screening, effective CRT regimens are also suggested for patients showing less sensitivity to radiation. Altogether, our study provides an example of applying integrative pharmacogenomics to establish a personalized precision oncology for NPC subtype-guided therapies.

Clonal relationship and directionality of progression of synchronous endometrial and ovarian carcinomas in patients with DNA mismatch repair-deficiency associated syndromes

Sporadic synchronous endometrial (ECs) and ovarian cancers (OCs), although clinically considered to be independent primaries, have been shown to be clonally related and likely constitute metastases from each other. We sought to define whether synchronous ECs/OCs in patients with DNA mismatch repair (MMR)-deficiency syndromes would be clonally related. We subjected synchronous ECs/OCs from four patients (LS3-LS6) with clinically confirmed Lynch syndrome (LS) and one patient with constitutional mismatch repair-deficiency syndrome (CMMRD) to massively parallel sequencing targeting 468 cancer-related genes. Somatic mutations, copy number alterations (CNAs), clonal relatedness and clonal decomposition analyses were performed using previously described bioinformatics methods. All synchronous ECs/OCs analyzed were considered independent primaries based on clinicopathologic criteria. Sequencing analysis revealed that the ECs/OCs of three cases (LS2-CMMRD, L3, L4) harbored similar repertoires of somatic mutations and CNAs and were clonally related. In these three cases, a subset of subclonal mutations in the EC became clonal in the OC, suggesting that the EC was likely the substrate from which the OC developed. LS5's EC/OC had distinct mutational profiles but shared specific CNAs. In contrast, LS6's EC/OC harbored distinct somatic mutations and lacked CNAs, consistent with each tumor constituting an independent primary lesion. In LS5 and LS6, PTEN mutations and PTEN loss of protein expression were found to be restricted to the EC. Finally, re-analysis of sequencing data of sporadic synchronous ECs/OCs supported the observations made in the current study that the directionality of progression is likely from the endometrium to the ovary. In conclusion, contrary to sporadic synchronous ECs/OCs, which are almost invariably clonally related, ECs/OCs simultaneously involving the uterus and ovary in LS patients may represent distinct primary tumors. A subset of MMR-deficiency syndrome-related synchronous ECs/OCs, however, may originate from a single primary tumor at variance with their clinical diagnosis, with the endometrium being the likeliest site of origin.

Genetic characterisation of adult primary pleomorphic uterine rhabdomyosarcoma and comparison with uterine carcinosarcoma

AIMS

To characterise the genetic alterations in adult primary uterine rhabdomyosarcomas (uRMSs) and to investigate whether these tumours are genetically distinct from uterine carcinosarcomas (UCSs).

METHODS AND RESULTS

Three tumours originally diagnosed as primary adult pleomorphic uRMS were subjected to massively parallel sequencing targeting 468 cancer-related genes and RNA-sequencing. Mutational profiles were compared with those of UCSs (n = 57) obtained from The Cancer Genome Atlas. Sequencing data analyses were performed using validated bioinformatic approaches. Pathogenic TP53 mutations and high levels of genomic instability were detected in the three cases. uRMS1 harboured a likely pathogenic YTHDF2-FOXR1 fusion. uRMS2 harboured a PPP2R1A hotspot mutation and amplification of multiple genes, including WHSC1L1, FGFR1, MDM2, and CCNE1, whereas uRMS3 harboured an FBXW7 hotspot mutation and an ANKRD11 homozygous deletion. Hierarchical clustering of somatic mutations and copy number alterations revealed that these tumours initially diagnosed as pleomorphic uRMSs and UCSs were similar. Subsequent comprehensive pathological re-review of the three uRMSs revealed previously unidentified minute pan-cytokeratin-positive atypical glands in one case (uRMS3), favouring its reclassification as UCS with extensive rhabdomyosarcomatous overgrowth.

CONCLUSIONS

Adult pleomorphic uRMSs harbour TP53 mutations and high levels of copy number alterations. Our findings underscore the challenge in discriminating between uRMS and UCS with rhabdomyosarcomatous differentiation.

TERT promoter hotspot mutations and gene amplification in metaplastic breast cancer

Metaplastic breast cancers (MBCs) are characterized by complex genomes, which seem to vary according to their histologic subtype. TERT promoter hotspot mutations and gene amplification are rare in common forms of breast cancer, but present in a subset of phyllodes tumors. Here, we sought to determine the frequency of genetic alterations affecting TERT in a cohort of 60 MBCs with distinct predominant metaplastic components (squamous, 23%; spindle, 27%; osseous, 8%; chondroid, 42%), and to compare the repertoire of genetic alterations of MBCs according to the presence of TERT promoter hotspot mutations or gene amplification. Forty-four MBCs were subjected to: whole-exome sequencing (WES; n = 27) or targeted sequencing of 341-468 cancer-related genes (n = 17); 16 MBCs were subjected to Sanger sequencing of the TERT promoter, TP53 and selected exons of PIK3CA, HRAS, and BRAF. TERT promoter hotspot mutations (n = 9) and TERT gene amplification (n = 1) were found in 10 of the 60 MBCs analyzed, respectively. These TERT alterations were less frequently found in MBCs with predominant chondroid differentiation than in other MBC subtypes (p = 0.01, Fisher's exact test) and were mutually exclusive with TP53 mutations (p < 0.001, CoMEt). In addition, a comparative analysis of the MBCs subjected to WES or targeted cancer gene sequencing (n = 44) revealed that MBCs harboring TERT promoter hotspot mutations or gene amplification (n = 6) more frequently harbored PIK3CA than TERT wild-type MBCs (n = 38; p = 0.001; Fisher's exact test). In conclusion, TERT somatic genetic alterations are found in a subset of TP53 wild-type MBCs with squamous/spindle differentiation, highlighting the genetic diversity of these cancers.

Moss enables high sensitivity single-nucleotide variant calling from multiple bulk DNA tumor samples

Intra-tumor heterogeneity renders the identification of somatic single-nucleotide variants (SNVs) a challenging problem. In particular, low-frequency SNVs are hard to distinguish from sequencing artifacts. While the increasing availability of multi-sample tumor DNA sequencing data holds the potential for more accurate variant calling, there is a lack of high-sensitivity multi-sample SNV callers that utilize these data. Here we report Moss, a method to identify low-frequency SNVs that recur in multiple sequencing samples from the same tumor. Moss provides any existing single-sample SNV caller the ability to support multiple samples with little additional time overhead. We demonstrate that Moss improves recall while maintaining high precision in a simulated dataset. On multi-sample hepatocellular carcinoma, acute myeloid leukemia and colorectal cancer datasets, Moss identifies new low-frequency variants that meet manual review criteria and are consistent with the tumor's mutational signature profile. In addition, Moss detects the presence of variants in more samples of the same tumor than reported by the single-sample caller. Moss' improved sensitivity in SNV calling will enable more detailed downstream analyses in cancer genomics.

Genomic characterization of small cell carcinomas of the uterine cervix

Small cell carcinoma (SCC) of the uterine cervix is a rare and aggressive form of neuroendocrine carcinoma, which resembles small cell lung cancer (SCLC) in its histology and poor survival rate. Here, we sought to define the genetic underpinning of SCCs of the uterine cervix and compare their mutational profiles with those of human papillomavirus (HPV)‐positive head and neck squamous cell carcinomas, HPV‐positive cervical carcinomas, and SCLCs using publicly available data. Using a combination of whole‐exome and targeted massively parallel sequencing, we found that the nine uterine cervix SCCs, which were HPV18‐positive (n = 8) or HPV16‐positive (n = 1), harbored a low mutation burden, few copy number alterations, and other than TP53 in two cases no recurrently mutated genes. The majority of mutations were likely passenger missense mutations, and only few affected previously described cancer‐related genes. Using RNA‐sequencing, we identified putative viral integration sites on 18q12.3 and on 8p22 in two SCCs of the uterine cervix. The overall nonsilent mutation rate of uterine cervix SCCs was significantly lower than that of SCLCs, HPV‐driven cervical adeno‐ and squamous cell carcinomas, or HPV‐positive head and neck squamous cell carcinomas. Unlike SCLCs, which are reported to harbor almost universal TP53 and RB1 mutations and a dominant tobacco smoke‐related signature 4, uterine cervix SCCs rarely harbored mutations affecting these genes (2/9, 22% TP53; 0% RB1) and displayed a dominant aging (67%) or APOBEC mutational signature (17%), akin to HPV‐driven cancers, including cervical adeno‐ and squamous cell carcinomas and head and neck squamous cell carcinomas. Taken together, in contrast to SCLCs, which are characterized by highly recurrent TP53 and RB1 alterations, uterine cervix SCCs were positive for HPV leading to inactivation of the suppressors p53 and RB, suggesting that these SCCs are convergent phenotypes.

Molecular characterization of high-grade serous ovarian cancers occurring in younger and older women

OBJECTIVES

To determine if the mutational landscapes and genomic features of homologous recombination DNA repair defects (HRD) vary between younger and older patients with high-grade serous ovarian cancer (HGSOC).

METHODS

Younger and older women were defined as bottom and top age quartiles, respectively. HGSOCs from 15 younger (median 49 years, range 35-53) and 15 older women (median 72 years, range 70-87) were subjected to whole-exome sequencing (WES). For validation, HGSOC WES data were obtained from The Cancer Genome Atlas (TCGA), including 38 younger (median 45 years, range 34-50) and 30 older women (median 74 years, range 68-84). Mutational profiles, BRCA1/2 status, genomic HRD features, and for TCGA cases RNA-sequencing-based HRD transcriptomic signatures were assessed.

RESULTS

In the institutional cohort, pathogenic germline BRCA1/2 mutations were more frequent in younger (5/15) than older women (0/15, p = 0.042). No somatic BRCA1/2 mutations were identified. HGSOCs from older patients preferentially displayed aging-related mutational signatures and, in contrast to younger patients, harbored CCNE1 amplifications (3/15, 20%). In the TCGA cohort, pathogenic germline BRCA1 (younger 8/38, older 0/30, p = 0.007) but not BRCA2 mutations (young 3/38, older 4/30, p = 0.691) were more frequent in younger patients. Again, no somatic BRCA1/2 mutations were identified. HGSOCs from younger women more frequently displayed genomic features of HRD (all, p < 0.05), a significant HRD gene-signature enrichment, but less frequently CCNE1 amplification (p = 0.05). Immunoreactive CLOVAR subtypes were more common in HGSOCs from younger women, and proliferative subtypes in HGSOCs from older women (p = 0.041).

CONCLUSIONS

HGSOC patients diagnosed at an older age less frequently harbor pathogenic BRCA1 germline mutations and genomic features of HRD than younger women. Individualized treatment options, particularly pertaining to use of PARP inhibitors, in older women may be warranted.

Structural variant detection in cancer genomes: computational challenges and perspectives for precision oncology

Cancer is generally characterized by acquired genomic aberrations in a broad spectrum of types and sizes, ranging from single nucleotide variants to structural variants (SVs). At least 30% of cancers have a known pathogenic SV used in diagnosis or treatment stratification. However, research into the role of SVs in cancer has been limited due to difficulties in detection. Biological and computational challenges confound SV detection in cancer samples, including intratumor heterogeneity, polyploidy, and distinguishing tumor-specific SVs from germline and somatic variants present in healthy cells. Classification of tumor-specific SVs is challenging due to inconsistencies in detected breakpoints, derived variant types and biological complexity of some rearrangements. Full-spectrum SV detection with high recall and precision requires integration of multiple algorithms and sequencing technologies to rescue variants that are difficult to resolve through individual methods. Here, we explore current strategies for integrating SV callsets and to enable the use of tumor-specific SVs in precision oncology.