Biostatistical Aspects of Whole Genome Sequencing Studies: Preprocessing and Quality Control.

Rapid advances in high-throughput DNA sequencing technologies have enabled large-scale whole genome sequencing (WGS) studies. Before performing association analysis between phenotypes and genotypes, preprocessing and quality control (QC) of the raw sequence data need to be performed. Because many biostatisticians have not been working with WGS data so far, we first sketch Illumina's short-read sequencing technology. Second, we explain the general preprocessing pipeline for WGS studies. Third, we provide an overview of important QC metrics, which are applied to WGS data: on the raw data, after mapping and alignment, after variant calling, and after multisample variant calling. Fourth, we illustrate the QC with the data from the GENEtic SequencIng Study Hamburg-Davos (GENESIS-HD), a study involving more than 9000 human whole genomes. All samples were sequenced on an Illumina NovaSeq 6000 with an average coverage of 35× using a PCR-free protocol. For QC, one genome in a bottle (GIAB) trio was sequenced in four replicates, and one GIAB sample was successfully sequenced 70 times in different runs. Fifth, we provide empirical data on the compression of raw data using the DRAGEN original read archive (ORA). The most important quality metrics in the application were genetic similarity, sample cross-contamination, deviations from the expected Het/Hom ratio, relatedness, and coverage. The compression ratio of the raw files using DRAGEN ORA was 5.6:1, and compression time was linear by genome coverage. In summary, the preprocessing, joint calling, and QC of large WGS studies are feasible within a reasonable time, and efficient QC procedures are readily available.

ROS1 genomic rearrangements are rare actionable drivers in microsatellite stable colorectal cancer.

c-Ros oncogene 1, receptor tyrosine kinase (ROS1) genomic rearrangements have been reported previously in rare cases of colorectal cancer (CRC), yet little is known about the frequency, molecular characteristics, and therapeutic vulnerabilities of ROS1-driven CRC. We analyzed a clinical dataset of 40 589 patients with CRC for ROS1 genomic rearrangements and their associated genomic characteristics (Foundation Medicine, Inc [FMI]). We moreover report the disease course and treatment response of an index patient with ROS1-rearranged metastatic CRC. ROS1 genomic rearrangements were identified in 34 (0.08%) CRC samples. GOPC-ROS1 was the most common ROS1 fusion identified (11 samples), followed by TTC28-ROS1 (3 samples). Four novel 5' gene partners of ROS1 were identified (MCM9, SRPK1, EPHA6, P4HA1). Contrary to previous reports on fusion-positive CRC, ROS1-rearrangements were found exclusively in microsatellite stable (MSS) CRCs. KRAS mutations were significantly less abundant in ROS1-rearranged vs ROS1 wild type cases. The index patient presented with chemotherapy-refractory metastatic right-sided colon cancer harboring GOPC-ROS1. Molecularly targeted treatment with crizotinib induced a rapid and sustained partial response. After 15 months on crizotinib disseminated tumor progression occurred and KRAS Q61H emerged in tissue and liquid biopsies. ROS1 rearrangements define a small, yet therapeutically actionable molecular subgroup of MSS CRC. In summary, the high prevalence of GOPC-ROS1 and noncanonical ROS1 fusions pose diagnostic challenges. We advocate NGS-based comprehensive molecular profiling of MSS CRCs that are wild type for RAS and BRAF and patient enrollment in precision trials.

Genomic landscape of pleural and peritoneal mesothelioma tumours.

BACKGROUND: Malignant pleural and peritoneal mesotheliomas are rare malignancies with unacceptable poor prognoses and limited treatment options. The genomic landscape is mainly characterised by the loss of tumour suppressor genes and mutations in DNA repair genes. Currently, data from next-generation sequencing (NGS) of mesothelioma tumours is restricted to a limited number of cases; moreover, data comparing molecular features of mesothelioma from the pleural and peritoneal origin with NGS are lacking. METHODS: We analysed 1113 pleural mesothelioma and 355 peritoneal mesothelioma samples. All tumours were sequenced with the FoundationOne® or FoundationOne®CDx assay for detection of substitutions, insertion-deletions, copy-number alterations and selected rearrangements in at least 324 cancer genes. RESULTS: This analysis revealed alterations in 19 genes with an overall prevalence of at least 2%. Alterations in BAP1, CDKN2A, CDKN2B, NF2, MTAP, TP53 and SETD2 occurred with a prevalence of at least 10%. Peritoneal, compared to pleural mesothelioma, was characterised by a lower prevalence of alterations in CDKN2A, CDKN2B and MTAP. Moreover, we could define four distinct subgroups according to alterations in BAP1 and CDKN2A/B. Alterations in Hedgehog pathway-related genes (PTCH1/2 and SUFU) and Hippo pathway-related gene (NF2) as well as KRAS, EGFR, PDGFRA/B, ERBB2 and FGFR3 were detected in both cohorts. CONCLUSION: Here, we report the molecular aberrations from the largest cohort of patients with mesothelioma. This analysis identified a proportion of patients with targetable alterations and suggests that molecular profiling can identify new treatment options for patients with mesothelioma.

Molecular and immunophenotypic characterization of SMARCB1 (INI1) - deficient intrathoracic Neoplasms.

The switch/sucrose-non-fermenting (SWI/SNF) complex is an ATP-dependent chromatin remodeling complex that plays important roles in DNA repair, transcription and cell differentiation. This complex consists of multiple subunits and is of particular interest in thoracic malignancies due to frequent subunit alteration of SMARCA4 (BRG1). Much less is known about SMARCB1 (INI1) deficient intrathoracic neoplasms, which are rare, often misclassified and understudied. In a retrospective analysis of 1479 intrathoracic malignant neoplasms using immunohistochemistry for INI1 (SMARCB1) on tissue micro arrays (TMA) and a search through our hospital sarcoma database, we identified in total nine intrathoracic, INI1 deficient cases (n = 9). We characterized these cases further by additional immunohistochemistry, broad targeted genomic analysis, methylation profiling and correlated them with clinical and radiological data. This showed that genomic SMARCB1 together with tumor suppressor alterations drive tumorigenesis in some of these cases, rather than epigenetic changes such as DNA methylation. A proper diagnostic classification, however, remains challenging. Intrathoracic tumors with loss or alteration of SMARCB1 (INI1) are highly aggressive and remain often underdiagnosed due to their rarity, which leads to false diagnostic interpretations. A better understanding of these tumors and proper diagnosis is important for better patient care as clinical trials and more targeted therapeutic options are emerging.

ERBB2-amplified lobular breast carcinoma exhibits concomitant CDK12 co-amplification associated with poor prognostic features.

Most invasive lobular breast carcinomas (ILBCs) are luminal-type carcinomas with an HER2-negative phenotype (ERBB2 or HER2 un-amplified) and CDH1 mutations. Rare variants include ERBB2-amplified subtypes associated with an unfavorable prognosis and less response to anti-HER2 targeted therapies. We analyzed the clinicopathological and molecular features of ERBB2-amplified ILBC and compared these characteristics with ERBB2-unamplified ILBC. A total of 253 patients with ILBC were analyzed. Paraffin-embedded formalin-fixed tumor samples from 250 of these patients were added to a tissue microarray. Protein expression of prognostic, stem cell and breast-specific markers was tested by immunohistochemistry (IHC). Hybrid capture-based comprehensive genomic profiling (CGP) was performed for 10 ILBCs that were either fluorescent in situ hybridization (FISH) or IHC positive for HER2 amplification/overexpression and 10 ILBCs that were either FISH or IHC negative. Results were compared with a CGP database of 44,293 invasive breast carcinomas. The CGP definition of ERBB2 amplification was five copies or greater. A total of 17 of 255 ILBC (5%) were ERBB2 amplified. ERBB2-amplified ILBC had higher tumor stage (p < 0.0001), more frequent positive nodal status (p = 0.00022), more distant metastases (p = 0.012), and higher histological grade (p < 0.0001), and were more often hormone receptor negative (p < 0.001) and more often SOX10 positive (p = 0.005). ERBB2 short variant sequence mutations were more often detected in ERBB2-unamplified tumors (6/10, p = 0.027), whereas CDH1 mutations/copy loss were frequently present in both subgroups (9/10 and 7/10, respectively). Amplification of pathogenic genes were more common in HER2-positive ILBC (p = 0.0009). CDK12 gene amplification (≥6 copies) was detected in 7 of 10 ERBB2-amplified ILBC (p = 0.018). There were no CDK12 gene amplifications reported in 44,293 invasive breast carcinomas in the FMI Insights CGP database. ERBB2-amplified ILBC is a distinct molecular subgroup with frequent coamplification of CDK12, whereas ERBB2 sequence mutations occur only in ERBB2-unamplified ILBC. CDK12/ERBB2 co-amplification may explain the poor prognosis and therapy resistance of ERBB2-amplified ILBC.

Molecular and functional landscape of malignant serous effusions for precision oncology.

Personalized treatment for patients with advanced solid tumors critically depends on the deep characterization of tumor cells from patient biopsies. Here, we comprehensively characterize a pan-cancer cohort of 150 malignant serous effusion (MSE) samples at the cellular, molecular, and functional level. We find that MSE-derived cancer cells retain the genomic and transcriptomic profiles of their corresponding primary tumors, validating their use as a patient-relevant model system for solid tumor biology. Integrative analyses reveal that baseline gene expression patterns relate to global ex vivo drug sensitivity, while high-throughput drug-induced transcriptional changes in MSE samples are indicative of drug mode of action and acquired treatment resistance. A case study exemplifies the added value of multi-modal MSE profiling for patients who lack genetically stratified treatment options. In summary, our study provides a functional multi-omics view on a pan-cancer solid tumor cohort and underlines the feasibility and utility of MSE-based precision oncology.

Reporting of somatic variants in clinical cancer care: recommendations of the Swiss Society of Molecular Pathology.

Somatic variant testing through next-generation sequencing (NGS) is well integrated into Swiss molecular pathology laboratories and has become a standard diagnostic method for numerous indications in cancer patient care. Currently, there is a wide variation in reporting practices within our country, and as patients move between different hospitals, it is increasingly necessary to standardize NGS reports to ease their reinterpretation. Additionally, as many different stakeholders-oncologists, hematologists, geneticists, pathologists, and patients-have access to the NGS report, it needs to contain comprehensive and detailed information in order to answer the questions of experts and avoid misinterpretation by non-experts. In 2017, the Swiss Institute of Bioinformatics conducted a survey to assess the differences in NGS reporting practices across ten pathology institutes in Switzerland. The survey examined 68 reporting items and identified 48 discrepancies. Based on these findings, the Swiss Society of Molecular Pathology initiated a Delphi method to reach a consensus on a set of recommendations for NGS reporting. Reports should include clinical information about the patient and the diagnosis, technical details about the sample and the test performed, and a list of all clinically relevant variants and variants of uncertain significance. In the absence of a consensus on an actionability scheme, the five-class pathogenicity scheme proposed by the ACMG/AMP guideline must be included in the reports. The Swiss Society of Molecular Pathology recognizes the importance of including clinical actionability in the report and calls on the European community of molecular pathologists and oncologists to reach a consensus on this issue.

Altered cytoplasmic and nuclear ADP-ribosylation levels analyzed with an improved ADP-ribose binder are a prognostic factor in renal cell carcinoma.

ADP-ribosylation (ADPR) of proteins is catalyzed by ADP-ribosyltransferases, which are targeted by inhibitors (i.e. poly(ADP-ribose) polymerase inhibitors [PARPi]). Although renal cell carcinoma (RCC) cells are sensitive in vitro to PARPi, studies on the association between ADPR levels and somatic loss of function mutations in DNA damage repair genes are currently missing. Here we observed, in two clear cell RCC (ccRCC) patient cohorts (n = 257 and n = 241) stained with an engineered ADP-ribose binding macrodomain (eAf1521), that decreased cytoplasmic ADPR (cyADPR) levels significantly correlated with late tumor stage, high-ISUP (the International Society of Urological Pathology) grade, presence of necrosis, dense lymphocyte infiltration, and worse patient survival (p < 0.01 each). cyADPR proved to be an independent prognostic factor (p = 0.001). Comparably, absence of nuclear ADPR staining in ccRCC correlated with absence of PARP1 staining (p < 0.01) and worse patient outcome (p < 0.05). In papillary RCC the absence of cyADPR was also significantly associated with tumor progression and worse patient outcome (p < 0.05 each). To interrogate whether the ADPR status could be associated with genetic alterations in DNA repair, chromatin remodeling, and histone modulation, we performed DNA sequence analysis and identified a significant association of increased ARID1A mutations in ccRCCcyADPR+++/PARP1+ compared with ccRCCcyADPR-/PARP1- (31% versus 4%; p < 0.05). Collectively, our data suggest the prognostic value of nuclear and cytoplasmic ADPR levels in RCC that might be further influenced by genetic alterations.

Acquired resistance to anti-PD1 therapy in patients with NSCLC associates with immunosuppressive T cell phenotype.

Immune checkpoint inhibitor treatment has the potential to prolong survival in non-small cell lung cancer (NSCLC), however, some of the patients develop resistance following initial response. Here, we analyze the immune phenotype of matching tumor samples from a cohort of NSCLC patients showing good initial response to immune checkpoint inhibitors, followed by acquired resistance at later time points. By using imaging mass cytometry and whole exome and RNA sequencing, we detect two patterns of resistance¨: One group of patients is characterized by reduced numbers of tumor-infiltrating CD8+ T cells and reduced expression of PD-L1 after development of resistance, whereas the other group shows high CD8+ T cell infiltration and high expression of PD-L1 in addition to markedly elevated expression of other immune-inhibitory molecules. In two cases, we detect downregulation of type I and II IFN pathways following progression to resistance, which could lead to an impaired anti-tumor immune response. This study thus captures the development of immune checkpoint inhibitor resistance as it progresses and deepens our mechanistic understanding of immunotherapy response in NSCLC.

Unravelling homologous recombination repair deficiency and therapeutic opportunities in soft tissue and bone sarcoma.

Defects in homologous recombination repair (HRR) in tumors correlate with poor prognosis and metastases development. Determining HRR deficiency (HRD) is of major clinical relevance as it is associated with therapeutic vulnerabilities and remains poorly investigated in sarcoma. Here, we show that specific sarcoma entities exhibit high levels of genomic instability signatures and molecular alterations in HRR genes, while harboring a complex pattern of chromosomal instability. Furthermore, sarcomas carrying HRDness traits exhibit a distinct SARC-HRD transcriptional signature that predicts PARP inhibitor sensitivity in patient-derived sarcoma cells. Concomitantly, HRDhigh sarcoma cells lack RAD51 nuclear foci formation upon DNA damage, further evidencing defects in HRR. We further identify the WEE1 kinase as a therapeutic vulnerability for sarcomas with HRDness and demonstrate the clinical benefit of combining DNA damaging agents and inhibitors of DNA repair pathways ex vivo and in the clinic. In summary, we provide a personalized oncological approach to treat sarcoma patients successfully.

MTPpilot: An Interactive Software for Visualization of Next-Generation Sequencing Results in Molecular Tumor Boards.

PURPOSE: Comprehensive targeted next-generation sequencing (NGS) panels are routinely used in modern molecular cancer diagnostics. In molecular tumor boards, the detected genomic alterations are often discussed to decide the next treatment options for patients with cancer. With the increasing size and complexity of NGS panels, the discussion of these results becomes increasingly complex, especially if they are reported in a text-based form, as it is the standard in current molecular pathology. METHODS: We have developed the Molecular Tumor Profiling pilot (MTPpilot) webservice using HTML, PHP, JavaScript, and MySQL to support the clinical discussion of NGS results at molecular tumor boards. RESULTS: MTPpilot integrates various public genome, network, and cancer mutation databases with interactive visualization tools to assess the functional impact of mutations and support clinical decision making at tumor boards. CONCLUSION: MTPpilot is tailored for discussion of NGS gene panel results at molecular tumor boards. It is freely available as a webservice at MTPpilot.

Comparison of calling pipelines for whole genome sequencing: an empirical study demonstrating the importance of mapping and alignment.

Rapid advances in high-throughput DNA sequencing technologies have enabled the conduct of whole genome sequencing (WGS) studies, and several bioinformatics pipelines have become available. The aim of this study was the comparison of 6 WGS data pre-processing pipelines, involving two mapping and alignment approaches (GATK utilizing BWA-MEM2 2.2.1, and DRAGEN 3.8.4) and three variant calling pipelines (GATK 4.2.4.1, DRAGEN 3.8.4 and DeepVariant 1.1.0). We sequenced one genome in a bottle (GIAB) sample 70 times in different runs, and one GIAB trio in triplicate. The truth set of the GIABs was used for comparison, and performance was assessed by computation time, F1 score, precision, and recall. In the mapping and alignment step, the DRAGEN pipeline was faster than the GATK with BWA-MEM2 pipeline. DRAGEN showed systematically higher F1 score, precision, and recall values than GATK for single nucleotide variations (SNVs) and Indels in simple-to-map, complex-to-map, coding and non-coding regions. In the variant calling step, DRAGEN was fastest. In terms of accuracy, DRAGEN and DeepVariant performed similarly and both superior to GATK, with slight advantages for DRAGEN for Indels and for DeepVariant for SNVs. The DRAGEN pipeline showed the lowest Mendelian inheritance error fraction for the GIAB trios. Mapping and alignment played a key role in variant calling of WGS, with the DRAGEN outperforming GATK.

Prognostic impact of tumour mutational burden in resected stage I and II lung adenocarcinomas from a European Thoracic Oncology Platform Lungscape cohort.

BACKGROUND: The primary objective of this study is to evaluate tumor mutational burden (TMB), its associations with selected clinicopathological and molecular characteristics as well as its clinical significance, in a retrospective cohort of surgically resected stage I-II lung adenocarcinomas, subset of the ETOP Lungscape cohort. METHODS: TMB was evaluated on tumor DNA extracted from resected primary lung adenocarcinomas, based on FoundationOne®CDx (F1CDx) genomic profiling, centrally performed at the University Hospital Zurich. The F1CDx test sequences the complete exons of 324 cancer-related genes and detects substitutions, insertions and deletions (indels), copy number alterations and gene rearrangements. In addition, the genomic biomarkers TMB and microsatellite instability (MSI) are analyzed. RESULTS: In the Lungscape cohort, TMB was assessed in 78 surgically resected lung adenocarcinomas from two Swiss centers (62 % males, 55 %/45 % stage I/II). Median TMB was 7.6 Muts/Mb, with TMB high (≥10 Muts/Mb) in 40 % of cases (95 %CI:29 %-52 %). The most frequently mutated genes were TP53/KRAS/EGFR/MLL2 detected in 58 %/38 %/33 %/30 % of samples, respectively. TMB was significantly higher among males (TMB high: 50 % vs 23 % in females, p = 0.032), as well as among current/former smokers (TMB high: 44 % vs 8 % in never smokers, p = 0.023). Furthermore, TMB was significantly higher in TP53 mutated than in non-mutated patients (TMB high: 60 % vs 12 %, p < 0.001), while it was higher in EGFR non-mutated patients compared to EGFR mutated (TMB high: 48 % vs 23 %, p = 0.049). At a median follow-up time of 56.1 months (IQR:38.8-72.0), none of the three outcome variables (OS, RFS, TTR) differed significantly by TMB status (all p-values > 5 %). This was also true when adjusting for clinicopathological characteristics. CONCLUSIONS: While presence of TP53 mutations and absence of EGFR mutations are associated with high TMB, increased TMB had no significant prognostic impact in patients with resected stage I/II lung adenocarcinoma beyond T and N classification, in both unadjusted and adjusted analyses.

Clinicopathological and Genomic Profiles of Atypical Fibroxanthoma and Pleomorphic Dermal Sarcoma Identify Overlapping Signatures with a High Mutational Burden.

Atypical fibroxanthoma (AFX) and pleomorphic dermal sarcoma (PDS) are rare tumors developing in chronically sun-exposed skin. Clinicopathological features are similar, but they differ in prognosis, while PDS has a more aggressive course with a higher risk for local recurrence and metastases. In current clinical practice, they are diagnosed by exclusion using immunohistochemistry. Thus, stringent diagnostic criteria and correct differentiation are critical in management and treatment for optimal outcomes. This retrospective single-center study collected clinicopathological data and tumor samples of 10 AFX and 18 PDS. Extracted genomic DNA from tumor specimens was analyzed by a next-generation sequencing (NGS) platform (FoundationOne-CDx™). Among 65 identified mutations, TP53 inactivating mutations were observed in all tumor specimens. In both AFX and PDS, the known pathogenic gene alterations in CDKN2A, TERT promoter, and NOTCH1 were frequently present, along with high mutational burden and stable Micro-Satellite Instability status. The mutational profiles differed only in ASXL1, which was only present in AFX. Further differences were identified in likely pathogenic and unknown gene alterations. Similarities in their genomic signatures could help to distinguish them from other malignancies, but they are not distinguishable between each other using the FoundationOne-CDx™ NGS panel. Therefore, histological criteria to determine diagnosis remain valid. For further insight, performing deep tumor profiling may be necessary.

Harmonization and Standardization of Panel-Based Tumor Mutational Burden Measurement: Real-World Results and Recommendations of the Quality in Pathology Study.

INTRODUCTION: Tumor mutational burden (TMB) is a quantitative assessment of the number of somatic mutations within a tumor genome. Immunotherapy benefit has been associated with TMB assessed by whole-exome sequencing (wesTMB) and gene panel sequencing (psTMB). The initiatives of Quality in Pathology (QuIP) and Friends of Cancer Research have jointly addressed the need for harmonization among TMB testing options in tissues. This QuIP study identifies critical sources of variation in psTMB assessment. METHODS: A total of 20 samples from three tumor types (lung adenocarcinoma, head and neck squamous cell carcinoma, and colon adenocarcinoma) with available WES data were analyzed for psTMB using six panels across 15 testing centers. Interlaboratory and interplatform variation, including agreement on variant calling and TMB classification, were investigated. Bridging factors to transform psTMB to wesTMB values were empirically derived. The impact of germline filtering was evaluated. RESULTS: Sixteen samples had low interlaboratory and interpanel psTMB variation, with 87.7% of pairwise comparisons revealing a Spearman's ρ greater than 0.6. A wesTMB cut point of 199 missense mutations projected to psTMB cut points between 7.8 and 12.6 mutations per megabase pair; the corresponding psTMB and wesTMB classifications agreed in 74.9% of cases. For three-tier classification with cut points of 100 and 300 mutations, agreement was observed in 76.7%, weak misclassification in 21.8%, and strong misclassification in 1.5% of cases. Confounders of psTMB estimation included fixation artifacts, DNA input, sequencing depth, genome coverage, and variant allele frequency cut points. CONCLUSIONS: This study provides real-world evidence that all evaluated panels can be used to estimate TMB in a routine diagnostic setting and identifies important parameters for reliable tissue TMB assessment that require careful control. As complex or composite biomarkers beyond TMB are likely playing an increasing role in therapy prediction, the efforts by QuIP and Friends of Cancer Research also delineate a general framework and blueprint for the evaluation of such assays.

Trithiocarbonates: exploration of a new head group for HDAC inhibitors.

Inhibition of histone deacetylases class I/II enzymes is a new, promising approach for cancer therapy. In the present study, we disclose a new structural class of HDAC inhibitors with the trithiocarbonate motif. A clear structure-activity-relationship was obtained for the cap-linker motif and the putative Zn(2+) complexing head group. Selected analogs display potent inhibition of HDAC enzymatic activity and a cellular potency comparable to that of suberoylanilide hydroxamic acid (SAHA), recently approved for treatment of patients with advanced cutaneous T-cell lymphoma.