Whole-exome sequencing of Nigerian benign prostatic hyperplasia reveals increased alterations in apoptotic pathways

BACKGROUND

Through whole-exome sequencing of 60 formalin-fixed paraffin-embedded Nigerian (NGRn) benign prostatic hyperplasia (BPH) samples, we identified germline and somatic alterations in apoptotic pathways impacting BPH development and progression. Prostate enlargement is a common occurrence in male aging; however, this enlargement can lead to lower urinary tract symptoms that negatively impact quality of life. This impact is disproportionately present in men of African ancestry. BPH pathophysiology is poorly understood and studies examining non-European populations are lacking.

METHODS

In this study, NGRn BPH, normal prostate, and prostate cancer (PCa) tumor samples were sequenced and compared to characterize genetic alterations in NGRn BPH.

RESULTS

Two hundred and two nonbenign, ClinVar-annotated germline variants were present in NGRn BPH samples. Six genes [BRCA1 (92%), HSD3B1 (85%), TP53 (37%), PMS2 (23%), BARD1 (20%), and BRCA2 (17%)] were altered in at least 10% of samples; however, compared to NGRn normal and tumor, the frequency of alterations in BPH samples showed no significant differences at the gene or variant level. BRCA2_rs11571831 and TP53_rs1042522 germline alterations had a statistically significant co-occurrence interaction in BPH samples. In at least two BPH samples, 173 genes harbored somatic variants known to be clinically actionable. Three genes (COL18A1, KIF16B, and LRP1) showed a statistically significant (p < 0.05) higher frequency in BPH. NGRn BPH also had five gene pairs (PKD1/KIAA0100, PKHD1/PKD1, DNAH9/LRP1B, NWD1/DCHS2, and TCERG1/LMTK2) with statistically significant co-occurring interactions. Two hundred and seventy-nine genes contained novel somatic variants in NGRn BPH. Three genes (CABP1, FKBP1C, and RP11-595B24.2) had a statistically significant (p < 0.05) higher alteration frequency in NGRn BPH and three were significantly higher in NGRn tumor (CACNA1A, DMKN, and CACNA2D2). Pairwise Fisher's exact tests showed 14 gene pairs with statistically significant (p < 0.05) interactions and four interactions approaching significance (p < 0.10). Mutational patterns in NGRn BPH were similar to COSMIC (Catalog of Somatic Mutations in Cancer) signatures associated with aging and dysfunctional DNA damage repair.

CONCLUSIONS

NGRn BPH contained significant germline alteration interactions (BRCA2_rs11571831 and TP53_rs1042522) and increased somatic alteration frequencies (LMTK2, LRP1, COL18A1, CABP1, and FKBP1C) that impact apoptosis. Normal prostate development is maintained by balancing apoptotic and proliferative activity. Dysfunction in either mechanism can lead to abnormal prostate growth. This work is the first to examine genomic sequencing in NGRn BPH and provides data that fill known gaps in the understanding BPH and how it impacts men of African ancestry.

Somatic variants as a cause of drug-resistant epilepsy including mesial temporal lobe epilepsy with hippocampal sclerosis

OBJECTIVE

The contribution of somatic variants to epilepsy has recently been demonstrated, particularly in the etiology of malformations of cortical development. The aim of this study was to determine the diagnostic yield of somatic variants in genes that have been previously associated with a somatic or germline epilepsy model, ascertained from resected brain tissue from patients with multidrug-resistant focal epilepsy.

METHODS

Forty-two patients were recruited across three categories: (1) malformations of cortical development, (2) mesial temporal lobe epilepsy with hippocampal sclerosis, and (3) nonlesional focal epilepsy. Participants were subdivided based on histopathology of the resected brain. Paired blood- and brain-derived DNA samples were sequenced using high-coverage targeted next generation sequencing to high depth (585× and 1360×, respectively). Variants were identified using Genome Analysis ToolKit (GATK4) MuTect-2 and confirmed using high-coverage Amplicon-EZ sequencing.

RESULTS

Sequence data on 41 patients passed quality control. Four somatic variants were validated following amplicon sequencing: within CBL, ALG13, MTOR, and FLNA. The diagnostic yield across 41 patients was 10%, 9% in mesial temporal lobe epilepsy with hippocampal sclerosis and 20% in malformations of cortical development.

SIGNIFICANCE

This study provides novel insights into the etiology of mesial temporal lobe epilepsy with hippocampal sclerosis, highlighting a potential pathogenic role of somatic variants in CBL and ALG13. We also report candidate diagnostic somatic variants in FLNA in focal cortical dysplasia, while providing further insight into the importance of MTOR and related genes in focal cortical dysplasia. This work demonstrates the potential molecular diagnostic value of variants in both germline and somatic epilepsy genes.

Enrichment of Immune Dysregulation Disorders in Adult Patients with Human Inborn Errors of Immunity

Human inborn errors of immunity (IEI) comprise a group of diseases resulting from molecular variants that compromise innate and adaptive immunity. Clinical features of IEI patients are dominated by susceptibility to a spectrum of infectious diseases, as well as autoimmune, autoinflammatory, allergic, and malignant phenotypes that usually appear in childhood, which is when the diagnosis is typically made. However, some IEI patients are identified in adulthood due to symptomatic delay of the disease or other reasons that prevent the request for a molecular study. The application of next-generation sequencing (NGS) as a diagnostic technique has given rise to an ever-increasing identification of IEI-monogenic causes, thus improving the diagnostic yield and facilitating the possibility of personalized treatment. This work was a retrospective study of 173 adults with IEI suspicion that were sequenced between 2005 and 2023. Sanger, targeted gene-panel, and whole exome sequencing were used for molecular diagnosis. Disease-causing variants were identified in 44 of 173 (25.43%) patients. The clinical phenotype of these 44 patients was mostly related to infection susceptibility (63.64%). An enrichment of immune dysregulation diseases was found when cohorts with molecular diagnosis were compared to those without. Immune dysregulation disorders, group 4 from the International Union of Immunological Societies Expert Committee (IUIS), were the most prevalent among these adult patients. Immune dysregulation as a new item in the Jeffrey Model Foundation warning signs for adults significantly increases the sensitivity for the identification of patients with an IEI-producing molecular defect.

Somatic EPAS1 Variants in Pheochromocytoma and Paraganglioma in Patients With Sickle Cell Disease

CONTEXT

Somatic EPAS1 variants account for 5% to 8% of all pheochromocytoma and paragangliomas (PPGL) but are detected in over 90% of PPGL in patients with congenital cyanotic heart disease, where hypoxemia may select for EPAS1 gain-of-function variants. Sickle cell disease (SCD) is an inherited hemoglobinopathy associated with chronic hypoxia and there are isolated reports of PPGL in patients with SCD, but a genetic link between the conditions has yet to be established.

OBJECTIVE

To determine the phenotype and EPAS1 variant status of patients with PPGL and SCD.

METHODS

Records of 128 patients with PPGL under follow-up at our center from January 2017 to December 2022 were screened for SCD diagnosis. For identified patients, clinical data and biological specimens were obtained, including tumor, adjacent non-tumor tissue and peripheral blood. Sanger sequencing of exons 9 and 12 of EPAS1, followed by amplicon next-generation sequencing of identified variants was performed on all samples.

RESULTS

Four patients with both PPGL and SCD were identified. Median age at PPGL diagnosis was 28 years. Three tumors were abdominal paragangliomas and 1 was a pheochromocytoma. No germline pathogenic variants in PPGL-susceptibility genes were identified in the cohort. Genetic testing of tumor tissue detected unique EPAS1 variants in all 4 patients. Variants were not detected in the germline, and 1 variant was detected in lymph node tissue of a patient with metastatic disease.

CONCLUSION

We propose that somatic EPAS1 variants may be acquired through exposure to chronic hypoxia in SCD and drive PPGL development. Future work is needed to further characterize this association.

Genomic Alterations, Gene Expression Profiles and Functional Enrichment of Normal-Karyotype Acute Myeloid Leukaemia Based on Targeted Next-Generation Sequencing

Characterising genomic variants is paramount in understanding the pathogenesis and heterogeneity of normal-karyotype acute myeloid leukaemia (AML-NK). In this study, clinically significant genomic biomarkers were ascertained using targeted DNA sequencing and RNA sequencing on eight AML-NK patients' samples collected at disease presentation and after complete remission. In silico and Sanger sequencing validations were performed to validate variants of interest, and they were followed by the performance of functional and pathway enrichment analyses for overrepresentation analysis of genes with somatic variants. Somatic variants involving 26 genes were identified and classified as follows: 18/42 (42.9%) as pathogenic, 4/42 (9.5%) as likely pathogenic, 4/42 (9.5%) as variants of unknown significance, 7/42 (16.7%) as likely benign and 9/42 (21.4%) as benign. Nine novel somatic variants were discovered, of which three were likely pathogenic, in the CEBPA gene with significant association with its upregulation. Transcription misregulation in cancer tops the affected pathways involving upstream genes (CEBPA and RUNX1) that were deregulated in most patients during disease presentation and were closely related to the most enriched molecular function gene ontology category, DNA-binding transcription activator activity RNA polymerase II-specific (GO:0001228). In summary, this study elucidated putative variants and their gene expression profiles along with functional and pathway enrichment in AML-NK patients.

Interpreting variants in genes affected by clonal hematopoiesis in population data

Reference population databases like the Genome Aggregation Database (gnomAD) have improved our ability to interpret the human genome. Variant frequencies and frequency-derived tools (such as depletion scores) have become fundamental to variant interpretation and the assessment of variant-gene-disease relationships. Clonal hematopoiesis (CH) obstructs variant interpretation as somatic variants that provide proliferative advantage will affect variant frequencies, depletion scores, and downstream filtering. Further, default filtering of variants or genes associated with CH risks filtering bona fide germline variants as variants associated with CH can also cause Mendelian conditions. Here, we provide our insights on interpreting population variant data in genes affected by clonal hematopoiesis, as well as recommendations for careful review of 36 established CH genes associated with neurodevelopmental conditions.

Frequent somatic mosaicism in T lymphocyte subsets in individuals with and without multiple sclerosis

Background

Somatic variants are variations in an individual's genome acquired after the zygotic stadium and result from mitotic errors or not (fully) repaired DNA damage.

Objectives

To investigate whether somatic mosaicism in T lymphocyte subsets is enriched early in multiple sclerosis (MS).

Methods

We identified somatic variants with variant allele fractions ≥1% across the whole exome in CD4⁺ and CD8⁺ T lymphocytes of 21 treatment-naive MS patients with <5 years of disease duration and 16 partially age-matched healthy controls. We investigated the known somatic STAT3 variant p.Y640F in peripheral blood in a larger cohort of 446 MS patients and 259 controls.

Results

All subjects carried 1-142 variants in CD4⁺ or CD8⁺ T lymphocytes. Variants were more common, more abundant, and increased with age in CD8⁺ T lymphocytes. Somatic variants were common in the genes DNMT3A and especially STAT3. Overall, the presence or abundance of somatic variants, including the STAT3 p.Y640F variant, did not differ between MS patients and controls.

Conclusions

Somatic variation in T lymphocyte subsets is widespread in both control individuals and MS patients. Somatic mosaicism in T lymphocyte subsets is not enriched in early MS and thus unlikely to contribute to MS risk, but future research needs to address whether a subset of variants influences disease susceptibility.

Mosaic RASopathies: A review of disorders caused by somatic pathogenic variants in the genes of the RAS/MAPK pathway

Mosaic RASopathies are a heterogeneous group of diseases characterized by the presence at birth or early onset of congenital anomalies, cutaneous and vascular anomalies, segmental overgrowth, and increased cancer risk. They are caused by somatic pathogenic variants of the genes belonging the RAt Sarcoma Mitogen-activated protein kinase (RAS/MAPK) pathway causing its hyperactivation. Here, we review the clinical and molecular characteristics of this heterogeneous group of diseases, including the possibilities of molecular diagnosis and new therapeutic perspectives.

Somatic variants in diverse genes leads to a spectrum of focal cortical malformations

Post-zygotically acquired genetic variants, or somatic variants, that arise during cortical development have emerged as important causes of focal epilepsies, particularly those due to malformations of cortical development. Pathogenic somatic variants have been identified in many genes within the PI3K-AKT-mTOR-signalling pathway in individuals with hemimegalencephaly and focal cortical dysplasia (type II), and more recently in SLC35A2 in individuals with focal cortical dysplasia (type I) or non-dysplastic epileptic cortex. Given the expanding role of somatic variants across different brain malformations, we sought to delineate the landscape of somatic variants in a large cohort of patients who underwent epilepsy surgery with hemimegalencephaly or focal cortical dysplasia. We evaluated samples from 123 children with hemimegalencephaly (n = 16), focal cortical dysplasia type I and related phenotypes (n = 48), focal cortical dysplasia type II (n = 44), or focal cortical dysplasia type III (n = 15). We performed high-depth exome sequencing in brain tissue-derived DNA from each case and identified somatic single nucleotide, indel and large copy number variants. In 75% of individuals with hemimegalencephaly and 29% with focal cortical dysplasia type II, we identified pathogenic variants in PI3K-AKT-mTOR pathway genes. Four of 48 cases with focal cortical dysplasia type I (8%) had a likely pathogenic variant in SLC35A2. While no other gene had multiple disease-causing somatic variants across the focal cortical dysplasia type I cohort, four individuals in this group had a single pathogenic or likely pathogenic somatic variant in CASK, KRAS, NF1 and NIPBL, genes previously associated with neurodevelopmental disorders. No rare pathogenic or likely pathogenic somatic variants in any neurological disease genes like those identified in the focal cortical dysplasia type I cohort were found in 63 neurologically normal controls (P = 0.017), suggesting a role for these novel variants. We also identified a somatic loss-of-function variant in the known epilepsy gene, PCDH19, present in a small number of alleles in the dysplastic tissue from a female patient with focal cortical dysplasia IIIa with hippocampal sclerosis. In contrast to focal cortical dysplasia type II, neither focal cortical dysplasia type I nor III had somatic variants in genes that converge on a unifying biological pathway, suggesting greater genetic heterogeneity compared to type II. Importantly, we demonstrate that focal cortical dysplasia types I, II and III are associated with somatic gene variants across a broad range of genes, many associated with epilepsy in clinical syndromes caused by germline variants, as well as including some not previously associated with radiographically evident cortical brain malformations.

Targeted Sequencing of Plasma-Derived vs. Urinary cfDNA from Patients with Triple-Negative Breast Cancer

In breast cancer, the genetic profiling of circulating cell-free DNA (cfDNA) from blood plasma was shown to have good potential for clinical use. In contrast, only a few studies were performed investigating urinary cfDNA. In this pilot study, we analyzed plasma-derived and matching urinary cfDNA samples obtained from 15 presurgical triple-negative breast cancer patients. We used a targeted next-generation sequencing approach to identify and compare genetic alterations in both body fluids. The cfDNA concentration was higher in urine compared to plasma, but there was no significant correlation between matched samples. Bioinformatical analysis revealed a total of 3339 somatic breast-cancer-related variants (VAF ≥ 3%), whereof 1222 vs. 2117 variants were found in plasma-derived vs. urinary cfDNA, respectively. Further, 431 shared variants were found in both body fluids. Throughout the cohort, the recovery rate of plasma-derived mutations in matching urinary cfDNA was 47% and even 63% for pathogenic variants only. The most frequently occurring pathogenic and likely pathogenic mutated genes were NF1, CHEK2, KMT2C and PTEN in both body fluids. Notably, a pathogenic CHEK2 (T519M) variant was found in all 30 samples. Taken together, our results indicated that body fluids appear to be valuable sources bearing complementary information regarding the genetic tumor profile.

Molecular Profiling of Tumor Tissue in Mexican Patients with Colorectal Cancer

Colorectal cancer is a heterogeneous disease with multiple genomic changes that influence the clinical management of patients; thus, the search for new molecular targets remains necessary. The aim of this study was to identify genetic variants in tumor tissues from Mexican patients with colorectal cancer, using massive parallel sequencing. A total of 4813 genes were analyzed in tumoral DNA from colorectal cancer patients, using the TruSight One Sequencing panel. From these, 192 variants with clinical associations were found distributed in 168 different genes, of which 46 variants had not been previous reported in the literature or databases, although genes harboring those variants had already been described in colorectal cancer. Enrichment analysis of the affected genes was performed using Reactome software; pathway over-representation showed significance for disease, signal transduction, and immune system subsets in all patients, while exclusive subsets such as DNA repair, autophagy, and RNA metabolism were also found. Those characteristics, whether individual or shared, could give tumors specific capabilities for survival, aggressiveness, or response to treatment. Our results can be useful for future investigations targeting specific characteristics of tumors in colorectal cancer patients. The identification of exclusive or common pathways in colorectal cancer patients could be important for better diagnosis and personalized cancer treatment.

Unique molecular characteristics of NAFLD-associated liver cancer accentuate β-catenin/TNFRSF19-mediated immune evasion

BACKGROUND & AIMS

The hepatic manifestation of the metabolic syndrome, non-alcoholic fatty liver disease (NAFLD), can lead to the development of hepatocellular carcinoma (HCC). Despite a strong causative link, NAFLD-HCC is often underrepresented in systematic genome explorations.

METHODS

Herein, tumor-normal pairs from 100 patients diagnosed with NAFLD-HCC were subject to next-generation sequencing. Bioinformatic analyses were performed to identify key genomic, epigenomic and transcriptomic events associated with the pathogenesis of NAFLD-HCC. Establishment of primary patient-derived NAFLD-HCC culture was used as a representative human model for downstream in vitro investigations of the underlying CTNNB1 S45P driver mutation. A syngeneic immunocompetent mouse model was used to further test the involvement of CTNNB1^mutand TNFRSF19 in reshaping the tumor microenvironment.

RESULTS

Mutational processes operative in the livers of patients with NAFLD inferred susceptibility to tumor formation through defective DNA repair pathways. Dense promoter mutations and dysregulated transcription factors accentuated activated transcriptional regulation in NAFLD-HCC, in particular the enrichment of MAZ-MYC activities. Somatic events common in HCCs arising from NAFLD and viral hepatitis B infection underscore similar driver pathways, although an incidence shift highlights CTNNB1^mut dominance in NAFLD-HCC (33%). Immune exclusion correlated evidently with CTNNB1^mut. Chromatin immunoprecipitation-sequencing integrated with transcriptome and immune profiling revealed a unique transcriptional axis, wherein CTNNB1^mut leads to an upregulation of TNFRSF19 which subsequently represses senescence-associated secretory phenotype-like cytokines (including IL6 and CXCL8). This phenomenon could be reverted by the Wnt-modulator ICG001.

CONCLUSIONS

The unique mutational processes in the livers of patients with NAFLD and NAFLD-HCC allude to a "field effect" involving a gain-of-function role of CTNNB1 mutations in immune exclusion.

LAY SUMMARY

The increasing prevalence of metabolic syndrome in adult populations means that NAFLD is poised to be the major cause of liver cancer in the 21st century. We showed a strong "field effect" in the livers of patients with NAFLD, wherein activated β-catenin was involved in reshaping the tumor-immune microenvironment.

Inferring Potential Cancer Driving Synonymous Variants

Synonymous single nucleotide variants (sSNVs) are often considered functionally silent, but a few cases of cancer-causing sSNVs have been reported. From available databases, we collected four categories of sSNVs: germline, somatic in normal tissues, somatic in cancerous tissues, and putative cancer drivers. We found that screening sSNVs for recurrence among patients, conservation of the affected genomic position, and synVep prediction (synVep is a machine learning-based sSNV effect predictor) recovers cancer driver variants (termed proposed drivers) and previously unknown putative cancer genes. Of the 2.9 million somatic sSNVs found in the COSMIC database, we identified 2111 proposed cancer driver sSNVs. Of these, 326 sSNVs could be further tagged for possible RNA splicing effects, RNA structural changes, and affected RBP motifs. This list of proposed cancer driver sSNVs provides computational guidance in prioritizing the experimental evaluation of synonymous mutations found in cancers. Furthermore, our list of novel potential cancer genes, galvanized by synonymous mutations, may highlight yet unexplored cancer mechanisms.

Management of BRCA Tumour Testing in an Integrated Molecular Tumour Board Multidisciplinary Model

Tumour testing of the BRCA1/2 genes is routinely performed in patients with different cancer histological subtypes. To accurately identify patients with tumour-detected germline pathogenic variants (PVs) is a relevant issue currently under investigation. This study aims at evaluating the performance of the tumour-to-germline diagnostic flowchart model defined at our Institutional Molecular Tumour Board (MTB). Results from tumour BRCA sequencing of 641 consecutive unselected cancer patients were discussed during weekly MTB meetings with the early involvement of clinical geneticists for appropriate referral to genetic counselling. The overall tumour detection rate of BRCA1/2 PVs was 8.7% (56/641), ranging from 24.4% (31/127) in high-grade ovarian cancer to 3.9% (12/304) in tumours not associated with germline BRCA1/2 PVs. Thirty-seven patients with PVs (66%) were evaluated by a clinical geneticist, and in 24 of them (64.9%), germline testing confirmed the presence of the PV in blood. Nine of these patients (37.5%) were not eligible for germline testing according to the criteria in use at our institution. Cascade testing was subsequently performed on 18 relatives. The tumour-to-germline diagnostic pipeline, developed in the framework of our institutional MTB, compared with guideline-based germline testing following genetic counselling, proved to be effective in identifying a higher number of germline BRCA PVs carriers.

Validation of an NGS Panel Designed for Detection of Actionable Mutations in Tumors Common in Latin America

Next-generation sequencing (NGS) is progressively being used in clinical practice. However, several barriers preclude using this technology for precision oncology in most Latin American countries. To overcome some of these barriers, we have designed a 25-gene panel that contains predictive biomarkers for most current and near-future available therapies in Chile and Latin America. Library preparation was optimized to account for low DNA integrity observed in formalin-fixed paraffin-embedded tissue. The workflow includes an automated bioinformatic pipeline that accounts for the underrepresentation of Latin Americans in genome databases. The panel detected small insertions, deletions, and single nucleotide variants down to allelic frequencies of 0.05 with high sensitivity, specificity, and reproducibility. The workflow was validated in 272 clinical samples from several solid tumor types, including gallbladder (GBC). More than 50 biomarkers were detected in these samples, mainly in BRCA1/2, KRAS, and PIK3CA genes. In GBC, biomarkers for PARP, EGFR, PIK3CA, mTOR, and Hedgehog signaling inhibitors were found. Thus, this small NGS panel is an accurate and sensitive method that may constitute a more cost-efficient alternative to multiple non-NGS assays and costly, large NGS panels. This kind of streamlined assay with automated bioinformatics analysis may facilitate the implementation of precision medicine in Latin America.

Somatic copy number variants in neuropsychiatric disorders

Copy number variants (CNVs) have been implicated in neuropsychiatric disorders, with rare-inherited and de novo CNVs (dnCNVs) having large effects on disease liability. Recent studies started exploring a class of dnCNVs that occur post-zygotically, and are therefore present in some but not all cells of the body. Analogous to conditional mutations in animal models, the presence of risk mutations in a fraction of cells has the potential to enlighten how damaging mutations affect cell-type/cell-circuit specific pathologies leading to neuropsychiatric manifestations. Although mosaic CNVs appear to contribute to a modest fraction of risk (0.3-0.5%), expanding our insights about them with more sensitive experimental and statistical methods, has the potential to help clarify mechanisms of neuropsychiatric disease.

Adult-onset autoinflammation caused by somatic mutations in UBA1: A Dutch case series of patients with VEXAS

BACKGROUND

A novel autoinflammatory syndrome was recently described in male patients who harbored somatic mutations in the X-chromosomal UBA1 gene. These patients were characterized by adult-onset, treatment-refractory inflammation with fever, cytopenia, dysplastic bone marrow, vacuoles in myeloid and erythroid progenitor cells, cutaneous and pulmonary inflammation, chondritis, and vasculitis, which is abbreviated as VEXAS.

OBJECTIVE

This study aimed to (retrospectively) diagnose VEXAS in patients who had previously been registered as having unclassified autoinflammation. We furthermore aimed to describe clinical experiences with this multifaceted, complex disease.

METHODS

A systematic reanalysis of whole-exome sequencing data from a cohort of undiagnosed patients with autoinflammation from academic hospitals in The Netherlands was performed. When no sequencing data were available, targeted Sanger sequencing was applied in cases with high clinical suspicion of VEXAS.

RESULTS

A total of 12 male patients who carried mutations in UBA1 were identified. These patients presented with adult-onset (mean age 67 years, range 47-79 years) autoinflammation with systemic symptoms, elevated inflammatory parameters, and multiorgan involvement, most typically involving the skin and bone marrow. Novel features of VEXAS included interstitial nephritis, cardiac involvement, stroke, and intestinal perforation related to treatment with tocilizumab. Although many types of treatment were initiated, most patients became treatment-refractory, with a high mortality rate of 50%.

CONCLUSION

VEXAS should be considered in the differential diagnosis of males with adult-onset autoinflammation characterized by systemic symptoms and multiorgan involvement. Early diagnosis can prevent unnecessary diagnostic procedures and provide better prognostic information and more suitable treatment options, including stem cell transplantation.

Quantitative and Qualitative QC of Next-Generation Sequencing for Detecting Somatic Variants: An Example of Detecting Clonal Hematopoiesis of Indeterminate Potential

BACKGROUND

Because next-generation sequencing (NGS) for detecting somatic mutations has been adopted in clinical fields, both qualitative and quantitative QC of the somatic variants through whole coding regions detected by NGS is crucial. However, specific applications or guidelines, especially for quantitative QC, are currently insufficient. Our goal was to devise a practical approach for both quantitative and qualitative QC using an example of detecting clonal hematopoiesis of indeterminate potential (CHIP).

METHODS

We applied the QC scheme using commercial reference materials and in-house QC materials (IQCM) composed of haplotype map and cancer cell lines for monitoring CHIP.

RESULTS

This approach efficiently validated a customized CHIP NGS assay. Accuracy, analytical sensitivity, analytical specificity, qualitative precision (concordance), and limit of detection achieved were 99.87%, 98.53%, 100.00%, 100.00%, and 1.00%, respectively. The quantitative precision analysis also had a higher CV percentage at a lower alternative read depth (R2 = 0.749∼0.858). Use of IQCM ensured more than 100-fold reduction in the cost per run compared with that achieved using commercial reference materials.

CONCLUSION

Our approach determined the general analytical performance of NGS for detecting CHIP and recognized limitations such as lower precision at a lower level of variant burden. This approach could also be theoretically expanded to a general NGS assay for detecting somatic variants. Considering the reliable NGS results and cost-effectiveness, we propose the use of IQCM for QC of NGS assays at clinical laboratories.

Epithelial Cells of Deep Infiltrating Endometriosis Harbor Mutations in Cancer Driver Genes

Endometriosis is an inflammatory condition manifested by the presence of endometrial-like tissue outside of the uterine cavity. The most common clinical presentations of endometriosis are dysmenorrhea, infertility, and severe pelvic pain. Few hypotheses attempt to explain the pathogenesis of endometriosis; however, none of the theories have been fully confirmed or considered universal. We examined somatic mutations in eutopic endometrium samples, deep endometriotic nodules and peripheral blood from 13 women with deep endometriosis of the rectovaginal space. Somatic variants were identified in laser microdissected samples using next-generation sequencing. A custom panel of 1296 cancer-related genes was employed, and selected genes representing cancer drivers and non-drivers for endometrial and ovarian cancer were thoroughly investigated. All 59 detected somatic variants were of low mutated allele frequency (<10%). In deep ectopic lesions, detected variants were significantly more often located in cancer driver genes, whereas in eutopic endometrium, there was no such distribution. Our results converge with other reports, where cancer-related mutations were found in endometriosis without cancer, particularly recurrent KRAS mutations. Genetic alterations located in ectopic endometriotic nodules could contribute to their formation; nevertheless, to better understand the pathogenesis of this disease, more research in this area must be performed.

Somatic and de novo Germline Variants of MEDs in Human Neural Tube Defects

Background

Neural tube defects (NTDs) are among the most common and severe congenital defects in humans. Their genetic etiology is complex and remains poorly understood. The Mediator complex (MED) plays a vital role in neural tube development in animal models. However, no studies have yet examined the role of its human homolog in the etiology of NTDs.

Methods

In this study, 48 pairs of neural lesion site and umbilical cord tissues from NTD and 21 case-parent trios were involved in screening for NTD-related somatic and germline de novo variants. A series of functional cell assays were performed. We generated a Med12 p.Arg1784Cys knock-in mouse using CRISPR/Cas9 technology to validate the human findings.

Results

One somatic variant, MED12 p.Arg1782Cys, was identified in the lesion site tissue from an NTD fetus. This variant was absent in any other normal tissue from different germ layers of the same case. In 21 case-parent trios, one de novo stop-gain variant, MED13L p.Arg1760^∗, was identified. Cellular functional studies showed that MED12 p.Arg1782Cys decreased MED12 protein level and affected the regulation of MED12 on the canonical-WNT signaling pathway. The Med12 p.Arg1784Cys knock-in mouse exhibited exencephaly and spina bifida.

Conclusion

These findings provide strong evidence that functional variants of MED genes are associated with the etiology of some NTDs. We demonstrated a potentially important role for somatic variants in the occurrence of NTDs. Our study is the first study in which an NTD-related variant identified in humans was validated in mice using CRISPR/Cas9 technology.

Somatic TARDBP variants as a cause of semantic dementia

The aetiology of late-onset neurodegenerative diseases is largely unknown. Here we investigated whether de novo somatic variants for semantic dementia can be detected, thereby arguing for a more general role of somatic variants in neurodegenerative disease. Semantic dementia is characterized by a non-familial occurrence, early onset (<65 years), focal temporal atrophy and TDP-43 pathology. To test whether somatic variants in neural progenitor cells during brain development might lead to semantic dementia, we compared deep exome sequencing data of DNA derived from brain and blood of 16 semantic dementia cases. Somatic variants observed in brain tissue and absent in blood were validated using amplicon sequencing and digital PCR. We identified two variants in exon one of the TARDBP gene (L41F and R42H) at low level (1-3%) in cortical regions and in dentate gyrus in two semantic dementia brains, respectively. The pathogenicity of both variants is supported by demonstrating impaired splicing regulation of TDP-43 and by altered subcellular localization of the mutant TDP-43 protein. These findings indicate that somatic variants may cause semantic dementia as a non-hereditary neurodegenerative disease, which might be exemplary for other late-onset neurodegenerative disorders.

The Potential of Five Immune-Related Prognostic Genes to Predict Survival and Response to Immune Checkpoint Inhibitors for Soft Tissue Sarcomas Based on Multi-Omic Study

Low response rates to immunotherapy have been reported in soft tissue sarcoma (STS). There are few predictive biomarkers of response, and the tumor immune microenvironment associated with progression and prognosis remains unclear in STS. Gene expression data from the Cancer Genome Atlas were used to identify the immune-related prognostic genes (IRPGs) and construct the immune gene-related prognostic model (IGRPM). The tumor immune microenvironment was characterized to reveal differences between patients with different prognoses. Furthermore, somatic mutation data and DNA methylation data were analyzed to understand the underlying mechanism leading to different prognoses. The IGRPM was constructed using five IRPGs (IFIH1, CTSG, STC2, SECTM1, and BIRC5). Two groups (high- and low-risk patients) were identified based on the risk score. Low-risk patients with higher overall survival time had higher immune scores, more immune cell infiltration (e.g., CD8 T cell and activated natural killer cells), higher expression of immune-stimulating molecules, higher stimulating cytokines and corresponding receptors, higher innate immunity molecules, and stronger antigen-presenting capacity. However, inhibition of immunity was observed in low-risk patients owing to the higher expression of immune checkpoint molecules and inhibiting cytokines. High-risk patients had high tumor mutation burden, which did not significantly influence survival. Gene set enrichment analysis further revealed that pathways of cell cycle and cancers were activated in high-risk patients. DNA methylation analysis indicated that relative high methylation was associated with better overall survival. Finally, the age, mitotic counts, and risk scores were independent prognostic factors for STS. Five IRPGs performed well in risk stratification of patients and are candidate biomarkers for predicting response to immunotherapy. Differences observed through the multi-omic study of patients with different prognoses may reveal the underlying mechanism of the development and progression of STS, and thereby improve treatment.

Somatic variants in new candidate genes identified in focal cortical dysplasia type II

OBJECTIVE

Focal cortical dysplasia type II (FCDII) is a malformation of cortex development commonly found in children with drug-resistant epilepsy. FCDII has been associated with somatic mutations in mammalian target of rapamycin (mTOR)-related pathway genes and an upregulation of mTOR. Somatic mutations were found in 10%-63% of FCDII samples; the frequency of the mutant allele was 0.93%-33.5%. This study aimed to find new candidate genes involved in FCDII.

METHODS

We collected resected FCD lesions, perilesional brain tissues, and peripheral blood from 17 children with pathologically confirmed FCDII. We performed whole exome sequencing and followed a set of screening and analysis strategies to identify potentially deleterious somatic variants (PDSVs) in brain-expressed genes. We performed site-specific amplicon sequencing to validate the results. We also performed an in vitro functional study on an IRS1 variant.

RESULTS

In six of 17 samples, we identified seven PDSVs in seven genes, including two frameshift variants and five missense variants. The frequencies of the variant allele were 1.29%-5.50%. The genes were MTOR, TSC2, IRS1, RAB6B, RALA, HTR6, and ZNF337. PDSVs in IRS1, RAB6B, ZNF337, RALA, and HTR6 had not been previously associated with FCD. In one lesion, two PDSVs were found in two genes. In a transfected cell line, we demonstrated that the c.1791dupG (identified in FCDII from Patient 1) led to a truncated IRS1 and significant mTOR hyperactivation compared to cells that carried wild-type IRS1. mTOR was also activated in FCDII tissue from Patient 1.

SIGNIFICANCE

Seven PDSVs were identified in FCDII lesions in six of 17 children. Five variant genes had not been previously associated with cortical malformations. We demonstrated that the IRS1 variant led to mTOR hyperactivation in vitro. Although functional experiments are needed, the results provide evidence for novel candidate genes in the pathogenesis of FCDII.

A practical framework and online tool for mutational signature analyses show inter-tissue variation and driver dependencies

Mutational signatures are patterns of mutations that arise during tumorigenesis. We present an enhanced, practical framework for mutational signature analyses. Applying these methods on 3,107 whole genome sequenced (WGS) primary cancers of 21 organs reveals known signatures and nine previously undescribed rearrangement signatures. We highlight inter-organ variability of signatures and present a way of visualizing that diversity, reinforcing our findings in an independent analysis of 3,096 WGS metastatic cancers. Signatures with a high level of genomic instability are dependent on TP53 dysregulation. We illustrate how uncertainty in mutational signature identification and assignment to samples affects tumor classification, reinforcing that using multiple orthogonal mutational signature data is not only beneficial, it is essential for accurate tumor stratification. Finally, we present a reference web-based tool for cancer and experimentally-generated mutational signatures, called Signal (https://signal.mutationalsignatures.com), that also supports performing mutational signature analyses.

Somatic Variants: New Kids on the Block in Human Immunogenetics

Somatic variants are not inherited but acquired during an individual's lifetime, and individuals are increasingly considered as complex mosaics of genetically distinct cells. Whereas this concept is long-recognized in cancer, this review focuses on the growing role of somatic variants in immune cells in nonmalignant immune-related disorders, such as primary immunodeficiency and autoimmune diseases. Older case reports described somatic variants early in development, leading to large numbers of affected cells and severe phenotypes. Thanks to technological evolution, it is now feasible to detect somatic variants occurring later in life and affecting fewer cells. Hence, only recently is the scale at which somatic variants contribute to monogenic diseases being uncovered and is their contribution to complex diseases being explored systematically.

Massive computational identification of somatic variants in exonic splicing enhancers using The Cancer Genome Atlas

Owing to the development of next-generation sequencing (NGS) technologies, a large number of somatic variants have been identified in various types of cancer. However, the functional significance of most somatic variants remains unknown. Somatic variants that occur in exonic splicing enhancer (ESE) regions are thought to prevent serine and arginine-rich (SR) proteins from binding to ESE sequence motifs, which leads to exon skipping. We computationally identified somatic variants in ESEs by compiling numerous open-access datasets from The Cancer Genome Atlas (TCGA). Using somatic variants and RNA-seq data from 9635 patients across 32 TCGA projects, we identified 646 ESE-disrupting variants. The false positive rate of our method, estimated using a permutation test, was approximately 1%. Of these ESE-disrupting variants, approximately 71% were located in the binding motifs of four classical SR proteins. ESE-disrupting variants occurred in proportion to the number of somatic variants, but not necessarily in the specific genes associated with the biological processes of cancer. Existing bioinformatics tools could not predict the pathogenicity of ESE-disrupting variants identified in this study, although these variants could cause exon skipping. We demonstrated that ESE-disrupting nonsense variants tended to escape nonsense-mediated decay surveillance. Using integrated analyses of open access data, we could specifically identify ESE-disrupting variants. We have generated a powerful tool, which can handle datasets without normal samples or raw data, and thus contribute to reducing variants of uncertain significance because our statistical approach only uses the exon-junction read counts from the tumor samples.

A Canadian guideline on the use of next-generation sequencing in oncology

Rapid advancements in next-generation sequencing (ngs) technology have created an unprecedented opportunity to decipher the molecular profile of tumours to more effectively prevent, diagnose, and treat cancer. Oncologists now have the option to order molecular tests that can guide treatment decisions. However, to date, most oncologists have received limited training in genomics, and they are now faced with the challenge of understanding how such tests and their interpretation align with patient management. Guidance on how to effectively use ngs technology is therefore needed to aid oncologists in applying the results of genomic tests. The Canadian guideline presented here describes best practices and unmet needs related to ngs-based testing for somatic variants in oncology, including clinical application, assay and sample selection, bioinformatics and interpretation of reports performed by laboratories, patient communication, and clinical trials.

Risk stratification by somatic mutation burden in Ewing sarcoma

BACKGROUND

Up to one-third of patients with localized Ewing sarcoma (ES) develop recurrent disease, but current biomarkers do not accurately identify this high-risk group. Therefore, the objective of this study was to determine the utility of mutational burden in predicting outcomes in patients with localized ES.

METHODS

Clinical and genomic data from 99 patients with ES, of whom 63 had localized disease at diagnosis, were obtained from the cBioPortal for Cancer Genomics. Genomic data included the type and number of somatic mutations using cBioPortal mutation calling. Primary endpoints were overall survival (OS) and the time to progression (TTP).

RESULTS

Patients had a median number of 11 somatic mutations. Patients were stratified according to whether they had a lower or higher mutational burden if they had ≤11 or >11 mutations, respectively. Higher mutational burden was significantly associated with inferior OS and TTP, a finding that was confirmed by univariate and multivariable analyses. In patients who had localized disease at diagnosis, higher mutational burden was the only variable significantly associated with inferior OS and TTP. The presence of a mutation in either stromal antigen 2 (STAG2) or tumor protein 53 (TP53), both of which were correlated previously with shorter OS in patients with ES, were significantly associated with higher mutational burden. Upon stratifying patients who had localized disease based on a standard panel of cancer genes, higher risk stratification was correlated significantly with inferior TTP and trended toward significance with inferior OS.

CONCLUSIONS

Patients who have localized ES and a higher mutational burden have inferior OS and TTP compared with those who have lower mutation burden. The current findings suggest that the somatic mutation burden can be used to better risk stratify these patients and to guide clinical decision making.

Clinical whole-genome sequencing from routine formalin-fixed, paraffin-embedded specimens: pilot study for the 100,000 Genomes Project

PURPOSE

Fresh-frozen (FF) tissue is the optimal source of DNA for whole-genome sequencing (WGS) of cancer patients. However, it is not always available, limiting the widespread application of WGS in clinical practice. We explored the viability of using formalin-fixed, paraffin-embedded (FFPE) tissues, available routinely for cancer patients, as a source of DNA for clinical WGS.

METHODS

We conducted a prospective study using DNAs from matched FF, FFPE, and peripheral blood germ-line specimens collected from 52 cancer patients (156 samples) following routine diagnostic protocols. We compared somatic variants detected in FFPE and matching FF samples.

RESULTS

We found the single-nucleotide variant agreement reached 71% across the genome and somatic copy-number alterations (CNAs) detection from FFPE samples was suboptimal (0.44 median correlation with FF) due to nonuniform coverage. CNA detection was improved significantly with lower reverse crosslinking temperature in FFPE DNA extraction (80 °C or 65 °C depending on the methods). Our final data showed somatic variant detection from FFPE for clinical decision making is possible. We detected 98% of clinically actionable variants (including 30/31 CNAs).

CONCLUSION

We present the first prospective WGS study of cancer patients using FFPE specimens collected in a routine clinical environment proving WGS can be applied in the clinic.

Functional Interaction Between BRCA1 and DNA Repair in Yeast May Uncover a Role of RAD50, RAD51, MRE11A, and MSH6 Somatic Variants in Cancer Development

In this study, we determined if BRCA1 partners involved in DNA double-strand break (DSB) and mismatch repair (MMR) may contribute to breast and ovarian cancer development. Taking advantage the functional conservation of DNA repair pathways between yeast and human, we expressed several BRCA1 missense variants in DNA repair yeast mutants to identify functional interaction between BRCA1 and DNA repair in BRCA1-induced genome instability. The pathogenic p.C61G, pA1708E, p.M775R, and p.I1766S, and the neutral pS1512I BRCA1 variants increased intra-chromosomal recombination in the DNA-repair proficient strain RSY6. In the mre11, rad50, rad51, and msh6 deletion strains, the BRCA1 variants p.C61G, pA1708E, p.M775R, p.I1766S, and pS1215I did not increase intra-chromosomal recombination suggesting that a functional DNA repair pathway is necessary for BRCA1 variants to determine genome instability. The pathogenic p.C61G and p.I1766S and the neutral p.N132K, p.Y179C, and p.N550H variants induced a significant increase of reversion in the msh2Δ strain; the neutral p.Y179C and the pathogenic p.I1766S variant induced gene reversion also, in the msh6Δ strain. These results imply a functional interaction between MMR and BRCA1 in modulating genome instability. We also performed a somatic mutational screening of MSH6, RAD50, MRE11A, and RAD51 genes in tumor samples from 34 patients and identified eight pathogenic or predicted pathogenic rare missense variants: four in MSH6, one in RAD50, one in MRE11A, and two in RAD51. Although we found no correlation between BRCA1 status and these somatic DNA repair variants, this study suggests that somatic missense variants in DNA repair genes may contribute to breast and ovarian tumor development.

Identification of Novel Functional Variants of SIN3A and SRSF1 among Somatic Variants in Acute Myeloid Leukemia Patients

The advent of massively parallel sequencing, also called next-generation sequencing (NGS), has dramatically influenced cancer genomics by accelerating the identification of novel molecular alterations. Using a whole genome sequencing (WGS) approach, we identified somatic coding and noncoding variants that may contribute to leukemogenesis in 11 adult Korean acute myeloid leukemia (AML) patients, with serial tumor samples (primary and relapse) available for 5 of them; somatic variants were identified in 187 AML-related genes, including both novel (SIN3A, C10orf53, PTPRR, and RERGL) and well-known (NPM1, RUNX1, and CEPBA) AML-related genes. Notably, SIN3A expression shows prognostic value in AML. A newly designed method, referred to as "hot-zone" analysis, detected two putative functional noncoding variants that can alter transcription factor binding affinity near PPP1R10 and SRSF1. Moreover, the functional importance of the SRSF1 noncoding variant was further investigated by luciferase assays, which showed that the variant is critical for the regulation of gene expression leading to leukemogenesis. We expect that further functional investigation of these coding and noncoding variants will contribute to a more in-depth understanding of the underlying molecular mechanisms of AML and the development of targeted anti-cancer drugs.

Comprehensive genomic transcriptomic tumor-normal gene panel analysis for enhanced precision in patients with lung cancer

A CMS approved test for lung cancer is based on tumor-only analysis of a targeted 35 gene panel, specifically excluding the use of the patient's normal germline tissue. However, this tumor-only approach increases the risk of mistakenly identifying germline single nucleotide polymorphisms (SNPs) as somatically-derived cancer driver mutations (false positives). 621 patients with 30 different cancer types, including lung cancer, were studied to compare the precision of tumor somatic variant calling in 35 genes using tumor-only DNA sequencing versus tumor-normal DNA plus RNA sequencing. When sequencing of lung cancer was performed using tumor genomes alone without normal germline controls, 94% of variants identified were SNPs and thus false positives. Filtering for common SNPs still resulted in as high as 48% false positive variant calling. With tumor-only sequencing, 29% of lung cancer patients had a false positive variant call in at least one of twelve genes with directly targetable drugs. RNA analysis showed 18% of true somatic variants were not expressed. Thus, sequencing and analysis of both normal germline and tumor genomes is necessary for accurate identification of molecular targets. Treatment decisions based on tumor-only analysis may result in the administration of ineffective therapies while also increasing the risk of negative drug-related side effects.

ClinGen Cancer Somatic Working Group - standardizing and democratizing access to cancer molecular diagnostic data to drive translational research

A growing number of academic and community clinics are conducting genomic testing to inform treatment decisions for cancer patients (1). In the last 3-5 years, there has been a rapid increase in clinical use of next generation sequencing (NGS) based cancer molecular diagnostic (MolDx) testing (2). The increasing availability and decreasing cost of tumor genomic profiling means that physicians can now make treatment decisions armed with patient-specific genetic information. Accumulating research in the cancer biology field indicates that there is significant potential to improve cancer patient outcomes by effectively leveraging this rich source of genomic data in treatment planning (3). To achieve truly personalized medicine in oncology, it is critical to catalog cancer sequence variants from MolDx testing for their clinical relevance along with treatment information and patient outcomes, and to do so in a way that supports large-scale data aggregation and new hypothesis generation. One critical challenge to encoding variant data is adopting a standard of annotation of those variants that are clinically actionable. Through the NIH-funded Clinical Genome Resource (ClinGen) (4), in collaboration with NLM's ClinVar database and >50 academic and industry based cancer research organizations, we developed the Minimal Variant Level Data (MVLD) framework to standardize reporting and interpretation of drug associated alterations (5). We are currently involved in collaborative efforts to align the MVLD framework with parallel, complementary sequence variants interpretation clinical guidelines from the Association of Molecular Pathologists (AMP) for clinical labs (6). In order to truly democratize access to MolDx data for care and research needs, these standards must be harmonized to support sharing of clinical cancer variants. Here we describe the processes and methods developed within the ClinGen's Somatic WG in collaboration with over 60 cancer care and research organizations as well as CLIA-certified, CAP-accredited clinical testing labs to develop standards for cancer variant interpretation and sharing.

Detection of somatic variants in peripheral blood lymphocytes using a next generation sequencing multigene pan cancer panel

Next Generation Sequencing (NGS) multigene panels, which are routinely used to assess hereditary cancer risk, can detect both inherited germline variants and somatic variants in cancer-risk genes. We evaluated the frequency and distribution of likely somatic Pathogenic and Likely Pathogenic variants (PVs) detected in >220,000 individuals who underwent clinical testing with a 25-gene panel between September 2013 and March 2016. Likely somatic PVs are defined as variants with NGS read frequencies from 10% to 30%. Overall, 137 (0.06%) individuals were identified as carrying likely somatic PVs, most commonly in TP53 (73), CHEK2 (27), and ATM (20). Among this group, a second PV with a NGS read frequency consistent with a germline variant within the same gene or a different gene on the panel was detected in 21 individuals (15.3%), which is similar to the detection rate in our general testing population. Likely somatic PVs accounted for 38.8% of all PVs in TP53. In comparison, likely somatic PVs accounted for <1% of PVs in most other genes. Likely somatic PVs were more frequently identified in older individuals (p < 0.001). Additional studies are ongoing to further investigate the incidence and clinical implications of somatic variants, enabling the appropriate medical management for these patients.

Integrated analysis of oral tongue squamous cell carcinoma identifies key variants and pathways linked to risk habits, HPV, clinical parameters and tumor recurrence

Oral tongue squamous cell carcinomas (OTSCC) are a homogeneous group of tumors characterized by aggressive behavior, early spread to lymph nodes and a higher rate of regional failure. Additionally, the incidence of OTSCC among younger population (<50yrs) is on the rise; many of whom lack the typical associated risk factors of alcohol and/or tobacco exposure. We present data on single nucleotide variations (SNVs), indels, regions with loss of heterozygosity (LOH), and copy number variations (CNVs) from fifty-paired oral tongue primary tumors and link the significant somatic variants with clinical parameters, epidemiological factors including human papilloma virus (HPV) infection and tumor recurrence. Apart from the frequent somatic variants harbored in TP53, CASP8, RASA1, NOTCH and CDKN2A genes, significant amplifications and/or deletions were detected in chromosomes 6-9, and 11 in the tumors. Variants in CASP8 and CDKN2A were mutually exclusive. CDKN2A, PIK3CA, RASA1 and DMD variants were exclusively linked to smoking, chewing, HPV infection and tumor stage. We also performed a whole-genome gene expression study that identified matrix metalloproteases to be highly expressed in tumors and linked pathways involving arachidonic acid and NF-k-B to habits and distant metastasis, respectively. Functional knockdown studies in cell lines demonstrated the role of CASP8 in a HPV-negative OTSCC cell line. Finally, we identified a 38-gene minimal signature that predicts tumor recurrence using an ensemble machine-learning method. Taken together, this study links molecular signatures to various clinical and epidemiological factors in a homogeneous tumor population with a relatively high HPV prevalence.

Available resources and challenges for the clinical annotation of somatic variations

Next-generation sequencing (NGS) has become an important tool for identifying clinically relevant variants in both inherited disorders and oncology. Variants annotation that enables the creation of meaningful clinical reports often requires mining multiple publicly available databases. There are a number of such resources that have been designed to catalog and mine a plethora of germline variants or mutations. However, when analyzing tumor specimens in clinical settings, one may need to use different or ancillary resources that are specific for somatic variants or actionable mutations that may have clinical or treatment implications. The purpose of this review is to recapitulate the state of the art of somatic variation databases, which can aid in the clinical interpretation of NGS-based assays in oncology. In addition, the current need for collating various annotation sources into one-stop solutions to facilitate faster query execution and better integration into existing laboratory information systems are discussed.