CIViCdb 2022: evolution of an open-access cancer variant interpretation knowledgebase.

CIViC (Clinical Interpretation of Variants in Cancer; civicdb.org) is a crowd-sourced, public domain knowledgebase composed of literature-derived evidence characterizing the clinical utility of cancer variants. As clinical sequencing becomes more prevalent in cancer management, the need for cancer variant interpretation has grown beyond the capability of any single institution. CIViC contains peer-reviewed, published literature curated and expertly-moderated into structured data units (Evidence Items) that can be accessed globally and in real time, reducing barriers to clinical variant knowledge sharing. We have extended CIViC's functionality to support emergent variant interpretation guidelines, increase interoperability with other variant resources, and promote widespread dissemination of structured curated data. To support the full breadth of variant interpretation from basic to translational, including integration of somatic and germline variant knowledge and inference of drug response, we have enabled curation of three new Evidence Types (Predisposing, Oncogenic and Functional). The growing CIViC knowledgebase has over 300 contributors and distributes clinically-relevant cancer variant data currently representing >3200 variants in >470 genes from >3100 publications.

Phase 1/dose expansion trial of brentuximab vedotin and lenalidomide in relapsed or refractory diffuse large B-cell lymphoma.

New therapies are needed for patients with relapsed/refractory (rel/ref) diffuse large B-cell lymphoma (DLBCL) who do not benefit from or are ineligible for stem cell transplant and chimeric antigen receptor therapy. The CD30-targeted, antibody-drug conjugate brentuximab vedotin (BV) and the immunomodulator lenalidomide (Len) have demonstrated promising activity as single agents in this population. We report the results of a phase 1/dose expansion trial evaluating the combination of BV/Len in rel/ref DLBCL. Thirty-seven patients received BV every 21 days, with Len administered continuously for a maximum of 16 cycles. The maximum tolerated dose of the combination was 1.2 mg/kg BV with 20 mg/d Len. BV/Len was well tolerated with a toxicity profile consistent with their use as single agents. Most patients required granulocyte colony-stimulating factor support because of neutropenia. The overall response rate was 57% (95% CI, 39.6-72.5), complete response rate, 35% (95% CI, 20.7-52.6); median duration of response, 13.1 months; median progression-free survival, 10.2 months (95% CI, 5.5-13.7); and median overall survival, 14.3 months (95% CI, 10.2-35.6). Response rates were highest in patients with CD30+ DLBCL (73%), but they did not differ according to cell of origin (P = .96). NK cell expansion and phenotypic changes in CD8+ T-cell subsets in nonresponders were identified by mass cytometry. BV/Len represents a potential treatment option for patients with rel/ref DLBCL. This combination is being further explored in a phase 3 study (registered on https://clinicaltrials.org as NCT04404283). This trial was registered on https://clinicaltrials.gov as NCT02086604.

A comparative analysis of RAS variants in patients with disorders of somatic mosaicism.

PURPOSE: RAS genes (HRAS, KRAS, and NRAS) are commonly found to be mutated in cancers, and activating RAS variants are also found in disorders of somatic mosaicism (DoSM). A survey of the mutational spectrum of RAS variants in DoSM has not been performed. METHODS: A total of 938 individuals with suspected DoSM underwent high-sensitivity clinical next-generation sequencing-based testing. We investigated the mutational spectrum and genotype-phenotype associations of mosaic RAS variants. RESULTS: In this article, we present a series of individuals with DoSM with RAS variants. Classic hotspots, including Gly12, Gly13, and Gln61 constituted the majority of RAS variants observed in DoSM. Furthermore, we present 12 individuals with HRAS and KRAS in-frame duplication/insertion (dup/ins) variants in the switch II domain. Among the 18.3% individuals with RAS in-frame dup/ins variants, clinical findings were mainly associated with vascular malformations. Hotspots were associated with a broad phenotypic spectrum, including vascular tumors, vascular malformations, nevoid proliferations, segmental overgrowth, digital anomalies, and combinations of these. The median age at testing was higher and the variant allelic fraction was lower in individuals with in-frame dup/ins variants than those in individuals with mosaic RAS hotspots. CONCLUSION: Our work provides insight into the allelic and clinical heterogeneity of mosaic RAS variants in nonmalignant conditions.

Large scale genotype- and phenotype-driven machine learning in Von Hippel-Lindau disease.

Von Hippel-Lindau (VHL) disease is a hereditary cancer syndrome where individuals are predisposed to tumor development in the brain, adrenal gland, kidney, and other organs. It is caused by pathogenic variants in the VHL tumor suppressor gene. Standardized disease information has been difficult to collect due to the rarity and diversity of VHL patients. Over 4100 unique articles published until October 2019 were screened for germline genotype-phenotype data. Patient data were translated into standardized descriptions using Human Genome Variation Society gene variant nomenclature and Human Phenotype Ontology terms and has been manually curated into an open-access knowledgebase called Clinical Interpretation of Variants in Cancer. In total, 634 unique VHL variants, 2882 patients, and 1991 families from 427 papers were captured. We identified relationship trends between phenotype and genotype data using classic statistical methods and spectral clustering unsupervised learning. Our analyses reveal earlier onset of pheochromocytoma/paraganglioma and retinal angiomas, phenotype co-occurrences and genotype-phenotype correlations including hotspots. It confirms existing VHL associations and can be used to identify new patterns and associations in VHL disease. Our database serves as an aggregate knowledge translation tool to facilitate sharing information about the pathogenicity of VHL variants.

Standards for the classification of pathogenicity of somatic variants in cancer (oncogenicity): Joint recommendations of Clinical Genome Resource (ClinGen), Cancer Genomics Consortium (CGC), and Variant Interpretation for Cancer Consortium (VICC).

PURPOSE: Several professional societies have published guidelines for the clinical interpretation of somatic variants, which specifically address diagnostic, prognostic, and therapeutic implications. Although these guidelines for the clinical interpretation of variants include data types that may be used to determine the oncogenicity of a variant (eg, population frequency, functional, and in silico data or somatic frequency), they do not provide a direct, systematic, and comprehensive set of standards and rules to classify the oncogenicity of a somatic variant. This insufficient guidance leads to inconsistent classification of rare somatic variants in cancer, generates variability in their clinical interpretation, and, importantly, affects patient care. Therefore, it is essential to address this unmet need. METHODS: Clinical Genome Resource (ClinGen) Somatic Cancer Clinical Domain Working Group and ClinGen Germline/Somatic Variant Subcommittee, the Cancer Genomics Consortium, and the Variant Interpretation for Cancer Consortium used a consensus approach to develop a standard operating procedure (SOP) for the classification of oncogenicity of somatic variants. RESULTS: This comprehensive SOP has been developed to improve consistency in somatic variant classification and has been validated on 94 somatic variants in 10 common cancer-related genes. CONCLUSION: The comprehensive SOP is now available for classification of oncogenicity of somatic variants.

Heterozygous HMGB1 loss-of-function variants are associated with developmental delay and microcephaly.

13q12.3 microdeletion syndrome is a rare cause of syndromic intellectual disability. Identification and genetic characterization of patients with 13q12.3 microdeletion syndrome continues to expand the phenotypic spectrum associated with it. Previous studies identified four genes within the approximately 300 Kb minimal critical region including two candidate protein coding genes: KATNAL1 and HMGB1. To date, no patients carrying a sequence-level variant or a single gene deletion in HMGB1 or KATNAL1 have been described. Here we report six patients with loss-of-function variants involving HMGB1 and who had phenotypic features similar to the previously described 13q12.3 microdeletion syndrome cases. Common features included developmental delay, language delay, microcephaly, obesity and dysmorphic features. In silico analyses suggest that HMGB1 is likely to be intolerant to loss-of-function, and previous in vitro data are in line with the role of HMGB1 in neurodevelopment. These results strongly suggest that haploinsufficiency of the HMGB1 gene may play a critical role in the pathogenesis of the 13q12.3 microdeletion syndrome.

Single-agent ibrutinib in relapsed or refractory follicular lymphoma: a phase 2 consortium trial.

Most patients with follicular lymphoma (FL) experience multiple relapses necessitating subsequent lines of therapy. Ibrutinib, a Bruton tyrosine kinase (BTK) inhibitor approved for the treatment of several B-cell malignancies, showed promising activity in FL in a phase 1 study. We report the results of a phase 2 trial evaluating ibrutinib in recurrent FL. Forty patients with recurrent FL were treated with ibrutinib 560 mg/d until progression or intolerance. The primary end point was overall response rate (ORR). Exploratory analyses included correlations of outcome with recurrent mutations identified in a cancer gene panel that used next-generation sequencing in pretreatment biopsies from 31 patients and results of early interim positron emission tomography/computed tomography scans in 20 patients. ORR was 37.5% with a complete response rate of 12.5%, median progression-free survival (PFS) of 14 months, and 2-year PFS of 20.4%. Response rates were significantly higher among patients whose disease was sensitive to rituximab (52.6%) compared with those who were rituximab refractory (16.7%) (P = .04). CARD11 mutations were present in 16% of patients (5 of 31) and predicted resistance to ibrutinib with only wild-type patients responding (P = .002). Maximum standardized uptake value at cycle 1 day 8 correlated with response and PFS. Ibrutinib was well-tolerated with a toxicity profile similar to labeled indications. Ibrutinib is a well-tolerated treatment with modest activity in relapsed FL. Evaluation of BTK inhibitors in earlier lines of therapy may be warranted on the basis of improved response rates in rituximab-sensitive disease. Somatic mutations such as CARD11 may have an impact on response to ibrutinib, may inform clinical decisions, and should be evaluated in larger data sets. This trial was registered at www.clinicaltrials.gov as #NCT01849263.

Standard operating procedure for somatic variant refinement of sequencing data with paired tumor and normal samples.

PURPOSE: Following automated variant calling, manual review of aligned read sequences is required to identify a high-quality list of somatic variants. Despite widespread use in analyzing sequence data, methods to standardize manual review have not been described, resulting in high inter- and intralab variability. METHODS: This manual review standard operating procedure (SOP) consists of methods to annotate variants with four different calls and 19 tags. The calls indicate a reviewer's confidence in each variant and the tags indicate commonly observed sequencing patterns and artifacts that inform the manual review call. Four individuals were asked to classify variants prior to, and after, reading the SOP and accuracy was assessed by comparing reviewer calls with orthogonal validation sequencing. RESULTS: After reading the SOP, average accuracy in somatic variant identification increased by 16.7% (p value = 0.0298) and average interreviewer agreement increased by 12.7% (p value < 0.001). Manual review conducted after reading the SOP did not significantly increase reviewer time. CONCLUSION: This SOP supports and enhances manual somatic variant detection by improving reviewer accuracy while reducing the interreviewer variability for variant calling and annotation.

Recurrent somatic mutations affecting B-cell receptor signaling pathway genes in follicular lymphoma.

Follicular lymphoma (FL) is the most common form of indolent non-Hodgkin lymphoma, yet it remains only partially characterized at the genomic level. To improve our understanding of the genetic underpinnings of this incurable and clinically heterogeneous disease, whole-exome sequencing was performed on tumor/normal pairs from a discovery cohort of 24 patients with FL. Using these data and mutations identified in other B-cell malignancies, 1716 genes were sequenced in 113 FL tumor samples from 105 primarily treatment-naive individuals. We identified 39 genes that were mutated significantly above background mutation rates. CREBBP mutations were associated with inferior PFS. In contrast, mutations in previously unreported HVCN1, a voltage-gated proton channel-encoding gene and B-cell receptor signaling modulator, were associated with improved PFS. In total, 47 (44.8%) patients harbor mutations in the interconnected B-cell receptor (BCR) and CXCR4 signaling pathways. Histone gene mutations were more frequent than previously reported (identified in 43.8% of patients) and often co-occurred (17.1% of patients). A novel, recurrent hotspot was identified at a posttranslationally modified residue in the histone H2B family. This study expands the number of mutated genes described in several known signaling pathways and complexes involved in lymphoma pathogenesis (BCR, Notch, SWitch/sucrose nonfermentable (SWI/SNF), vacuolar ATPases) and identified novel recurrent mutations (EGR1/2, POU2AF1, BTK, ZNF608, HVCN1) that require further investigation in the context of FL biology, prognosis, and treatment.

Adapting crowdsourced clinical cancer curation in CIViC to the ClinGen minimum variant level data community-driven standards.

Harmonization of cancer variant representation, efficient communication, and free distribution of clinical variant-associated knowledge are central problems that arise with increased usage of clinical next-generation sequencing. The Clinical Genome Resource (ClinGen) Somatic Working Group (WG) developed a minimal variant level data (MVLD) representation of cancer variants, and has an ongoing collaboration with Clinical Interpretations of Variants in Cancer (CIViC), an open-source platform supporting crowdsourced and expert-moderated cancer variant curation. Harmonization between MVLD and CIViC variant formats was assessed by formal field-by-field analysis. Adjustments to the CIViC format were made to harmonize with MVLD and support ClinGen Somatic WG curation activities, including four new features in CIViC: (1) introduction of an assertions feature for clinical variant assessment following the Association of Molecular Pathologists (AMP) guidelines, (2) group-level curation tracking for organizations, enabling member transparency, and curation effort summaries, (3) introduction of ClinGen Allele Registry IDs to CIViC, and (4) mapping of CIViC assertions into ClinVar submission with automated submissions. A generalizable workflow utilizing MVLD and new CIViC features is outlined for use by ClinGen Somatic WG task teams for curation and submission to ClinVar, and provides a model for promoting harmonization of cancer variant representation and efficient distribution of this information.

DGIdb 2.0: mining clinically relevant drug-gene interactions.

The Drug-Gene Interaction Database (DGIdb, www.dgidb.org) is a web resource that consolidates disparate data sources describing drug-gene interactions and gene druggability. It provides an intuitive graphical user interface and a documented application programming interface (API) for querying these data. DGIdb was assembled through an extensive manual curation effort, reflecting the combined information of twenty-seven sources. For DGIdb 2.0, substantial updates have been made to increase content and improve its usefulness as a resource for mining clinically actionable drug targets. Specifically, nine new sources of drug-gene interactions have been added, including seven resources specifically focused on interactions linked to clinical trials. These additions have more than doubled the overall count of drug-gene interactions. The total number of druggable gene claims has also increased by 30%. Importantly, a majority of the unrestricted, publicly-accessible sources used in DGIdb are now automatically updated on a weekly basis, providing the most current information for these sources. Finally, a new web view and API have been developed to allow searching for interactions by drug identifiers to complement existing gene-based search functionality. With these updates, DGIdb represents a comprehensive and user friendly tool for mining the druggable genome for precision medicine hypothesis generation.

Neoantigen Landscape Supports Feasibility of Personalized Cancer Vaccine for Follicular Lymphoma.

Personalized cancer vaccines designed to target neoantigens represent a promising new treatment paradigm in oncology. In contrast to classical idiotype vaccines, we hypothesized that 'polyvalent' vaccines could be engineered for the personalized treatment of follicular lymphoma (FL) using neoantigen discovery by combined whole exome sequencing (WES) and RNA sequencing (RNA-Seq). Fifty-eight tumor samples from 57 patients with FL underwent WES and RNA-Seq. Somatic and B-cell clonotype neoantigens were predicted and filtered to identify high-quality neoantigens. B-cell clonality was determined by alignment of B-cell receptor (BCR) CDR3 regions from RNA-Seq data, grouping at the protein level, and comparison to the BCR repertoire from healthy individuals using RNA-Seq data. An average of 52 somatic mutations per patient (range: 2-172) were identified, and two or more (median: 15) high-quality neoantigens were predicted for 56 of 58 FL samples. The predicted neoantigen peptides were composed of missense mutations (77%), indels (9%), gene fusions (3%), and BCR sequences (11%). Building off of these preclinical analyses, we initiated a pilot clinical trial using personalized neoantigen vaccination combined with PD-1 blockade in patients with relapsed or refractory FL (#NCT03121677). Synthetic long peptide (SLP) vaccines targeting predicted high-quality neoantigens were successfully synthesized for and administered to all four patients enrolled. Initial results demonstrate feasibility, safety, and potential immunologic and clinical responses. Our study suggests that a genomics-driven personalized cancer vaccine strategy is feasible for patients with FL, and this may overcome prior challenges in the field.

Knockdown of Hspa9, a del(5q31.2) gene, results in a decrease in hematopoietic progenitors in mice.

Heterozygous deletions spanning chromosome 5q31.2 occur frequently in the myelodysplastic syndromes (MDS) and are highly associated with progression to acute myeloid leukemia (AML) when p53 is mutated. Mutagenesis screens in zebrafish and mice identified Hspa9 as a del(5q31.2) candidate gene that may contribute to MDS and AML pathogenesis, respectively. To test whether HSPA9 haploinsufficiency recapitulates the features of ineffective hematopoiesis observed in MDS, we knocked down the expression of HSPA9 in primary human hematopoietic cells and in a murine bone marrow-transplantation model using lentivirally mediated gene silencing. Knockdown of HSPA9 in human cells significantly delayed the maturation of erythroid precursors, but not myeloid or megakaryocytic precursors, and suppressed cell growth by 6-fold secondary to an increase in apoptosis and a decrease in the cycling of cells compared with control cells. Erythroid precursors, B lymphocytes, and the bone marrow progenitors c-kit(+)/lineage(-)/Sca-1(+) (KLS) and megakaryocyte/erythrocyte progenitor (MEP) were significantly reduced in a murine Hspa9-knockdown model. These abnormalities suggest that cooperating gene mutations are necessary for del(5q31.2) MDS cells to gain clonal dominance in the bone marrow. Our results demonstrate that Hspa9 haploinsufficiency alters the hematopoietic progenitor pool in mice and contributes to abnormal hematopoiesis.

Accessing clinical-grade genomic classification data through the ClinGen Data Platform.

The Clinical Genome Resource (ClinGen) serves as an authoritative resource on the clinical relevance of genes and variants. In order to support our curation activities and to disseminate our findings to the community, we have developed a Data Platform of informatics resources backed by standardized data models. In this workshop we demonstrate our publicly available resources including curation interfaces, (Variant Curation Interface, CIViC), supporting infrastructure (Allele Registry, Genegraph), and data models (SEPIO, GA4GH VRS, VA).

Mutations associated with progression in follicular lymphoma predict inferior outcomes at diagnosis: Alliance A151303.

Follicular lymphoma (FL) is clinically heterogeneous, with select patients tolerating extended watch-and-wait, whereas others require prompt treatment, suffer progression of disease within 24 months of treatment (POD24), and/or experience aggressive histologic transformation (t-FL). Because our understanding of the relationship between genetic alterations in FL and patient outcomes remains limited, we conducted a clinicogenomic analysis of 370 patients with FL or t-FL (from Cancer and Leukemia Group B/Alliance trials 50402/50701/50803, or real-world cohorts from Washington University School of Medicine, Cleveland Clinic, or University of Miami). FL subsets by grade, stage, watch-and-wait, or POD24 status did not differ by mutation burden, whereas mutation burden was significantly higher in relapsed/refractory (rel/ref) FL and t-FL than in newly diagnosed (dx) FL. Nonetheless, mutation burden in dx FL was not associated with frontline progression-free survival (PFS). CREBBP was the only gene more commonly mutated in FL than in t-FL yet mutated CREBBP was associated with shorter frontline PFS in FL. Mutations in 20 genes were more common in rel/ref FL or t-FL than in dx FL, including 6 significantly mutated genes (SMGs): STAT6, TP53, IGLL5, B2M, SOCS1, and MYD88. We defined a mutations associated with progression (MAP) signature as ≥2 mutations in these 7 genes (6 rel/ref FL or t-FL SMGs plus CREBBP). Patients with dx FL possessing a MAP signature had shorter frontline PFS, revealing a 7-gene set offering insight into FL progression risk potentially more generalizable than the m7-Follicular Lymphoma International Prognostic Index (m7-FLIPI), which had modest prognostic value in our cohort. Future studies are warranted to validate the poor prognosis associated with a MAP signature in dx FL, potentially facilitating novel trials specifically in this high-risk subset of patients.

Ultra-Deep Sequencing Reveals the Mutational Landscape of Classical Hodgkin Lymphoma.

The malignant Hodgkin and Reed Sternberg (HRS) cells of classical Hodgkin lymphoma (cHL) are scarce in affected lymph nodes, creating a challenge to detect driver somatic mutations. As an alternative to cell purification techniques, we hypothesized that ultra-deep exome sequencing would allow genomic study of HRS cells, thereby streamlining analysis and avoiding technical pitfalls. To test this, 31 cHL tumor/normal pairs were exome sequenced to approximately 1,000× median depth of coverage. An orthogonal error-corrected sequencing approach verified >95% of the discovered mutations. We identified mutations in genes novel to cHL including: CDH5 and PCDH7, novel stop gain mutations in IL4R, and a novel pattern of recurrent mutations in pathways regulating Hippo signaling. As a further application of our exome sequencing, we attempted to identify expressed somatic single-nucleotide variants (SNV) in single-nuclei RNA sequencing (snRNA-seq) data generated from a patient in our cohort. Our snRNA analysis identified a clear cluster of cells containing a somatic SNV identified in our deep exome data. This cluster has differentially expressed genes that are consistent with genes known to be dysregulated in HRS cells (e.g., PIM1 and PIM3). The cluster also contains cells with an expanded B-cell clonotype further supporting a malignant phenotype. This study provides proof-of-principle that ultra-deep exome sequencing can be utilized to identify recurrent mutations in HRS cells and demonstrates the feasibility of snRNA-seq in the context of cHL. These studies provide the foundation for the further analysis of genomic variants in large cohorts of patients with cHL. SIGNIFICANCE: Our data demonstrate the utility of ultra-deep exome sequencing in uncovering somatic variants in Hodgkin lymphoma, creating new opportunities to define the genes that are recurrently mutated in this disease. We also show for the first time the successful application of snRNA-seq in Hodgkin lymphoma and describe the expression profile of a putative cluster of HRS cells in a single patient.

Impact of a 40-Gene Targeted Panel Test on Physician Decision Making for Patients With Acute Myeloid Leukemia.

PURPOSE: Physicians treating hematologic malignancies increasingly order targeted sequencing panels to interrogate recurrently mutated genes. The precise impact of these panels on clinical decision making is not well understood. METHODS: Here, we report our institutional experience with a targeted 40-gene panel (MyeloSeq) that is used to generate a report for both genetic variants and variant allele frequencies for the treating physician (the limit of mutation detection is approximately one AML cell in 50). RESULTS: In total, 346 sequencing reports were generated for 325 patients with suspected hematologic malignancies over an 8-month period (August 2018 to April 2019). To determine the influence of genomic data on clinical care for patients with acute myeloid leukemia (AML), we analyzed 122 consecutive reports from 109 patients diagnosed with AML and surveyed the treating physicians with a standardized questionnaire. The panel was ordered most commonly at diagnosis (61.5%), but was also used to assess response to therapy (22.9%) and to detect suspected relapse (15.6%). The panel was ordered at multiple timepoints during the disease course for 11% of patients. Physicians self-reported that 50 of 114 sequencing reports (44%) influenced clinical care decisions in 44 individual patients. Influences were often nuanced and extended beyond identifying actionable genetic variants with US Food and Drug Administration-approved drugs. CONCLUSION: This study provides insights into how physicians are currently using multigene panels capable of detecting relatively rare AML cells. The most influential way to integrate these tools into clinical practice will be to perform prospective clinical trials that assess patient outcomes in response to genomically driven interventions.

CIViCpy: A Python Software Development and Analysis Toolkit for the CIViC Knowledgebase.

PURPOSE: Precision oncology depends on the matching of tumor variants to relevant knowledge describing the clinical significance of those variants. We recently developed the Clinical Interpretations for Variants in Cancer (CIViC; civicdb.org) crowd-sourced, expert-moderated, and open-access knowledgebase. CIViC provides a structured framework for evaluating genomic variants of various types (eg, fusions, single-nucleotide variants) for their therapeutic, prognostic, predisposing, diagnostic, or functional utility. CIViC has a documented application programming interface for accessing CIViC records: assertions, evidence, variants, and genes. Third-party tools that analyze or access the contents of this knowledgebase programmatically must leverage this application programming interface, often reimplementing redundant functionality in the pursuit of common analysis tasks that are beyond the scope of the CIViC Web application. METHODS: To address this limitation, we developed CIViCpy (civicpy.org), a software development kit for extracting and analyzing the contents of the CIViC knowledgebase. CIViCpy enables users to query CIViC content as dynamic objects in Python. We assess the viability of CIViCpy as a tool for advancing individualized patient care by using it to systematically match CIViC evidence to observed variants in patient cancer samples. RESULTS: We used CIViCpy to evaluate variants from 59,437 sequenced tumors of the American Association for Cancer Research Project GENIE data set. We demonstrate that CIViCpy enables annotation of > 1,200 variants per second, resulting in precise variant matches to CIViC level A (professional guideline) or B (clinical trial) evidence for 38.6% of tumors. CONCLUSION: The clinical interpretation of genomic variants in cancers requires high-throughput tools for interoperability and analysis of variant interpretation knowledge. These needs are met by CIViCpy, a software development kit for downstream applications and rapid analysis. CIViCpy is fully documented, open-source, and available free online.

A harmonized meta-knowledgebase of clinical interpretations of somatic genomic variants in cancer.

Precision oncology relies on accurate discovery and interpretation of genomic variants, enabling individualized diagnosis, prognosis and therapy selection. We found that six prominent somatic cancer variant knowledgebases were highly disparate in content, structure and supporting primary literature, impeding consensus when evaluating variants and their relevance in a clinical setting. We developed a framework for harmonizing variant interpretations to produce a meta-knowledgebase of 12,856 aggregate interpretations. We demonstrated large gains in overlap between resources across variants, diseases and drugs as a result of this harmonization. We subsequently demonstrated improved matching between a patient cohort and harmonized interpretations of potential clinical significance, observing an increase from an average of 33% per individual knowledgebase to 57% in aggregate. Our analyses illuminate the need for open, interoperable sharing of variant interpretation data. We also provide a freely available web interface (search.cancervariants.org) for exploring the harmonized interpretations from these six knowledgebases.