Search

Neurofibromatosis type 1 (NF1) is an autosomal dominant genetic disorder characterized by the predisposition to develop tumors such as malignant peripheral nerve sheath tumors (MPNSTs) which represents the primary cause of death for NF1-affected patients. Regardless of the high incidence and mortality, the molecular mechanisms underneath MPNST growth and metastatic progression remain poorly understood. In this proof-of-concept study, we performed somatic whole-exome sequencing (WES) to profile the genomic alterations in four samples from a patient with NF1-associated MPNST, consisting of a benign plexiform neurofibroma, a primary MPNST, and metastases from lung and skin tissues. By comparing genomic patterns, we identified a high level of variability across samples with distinctive genetic changes which allow for the definition of profiles of the early phase with respect to the late metastatic stages. Pathogenic and likely pathogenic variants were abundant in the primary tumor, whereas the metastatic samples exhibited a high level of copy-number variations (CNVs), highlighting a possible genomic instability in the late phases. The most known MPNST-related genes, such as TP53 and SUZ12, were identified in CNVs observed within the primary tumor. Pathway analysis of altered early genes in MPNST pointed to a potential role in cell motility, division and metabolism. Moreover, we employed survival analysis with the TCGA sarcoma genomic dataset on 262 affected patients, in order to corroborate the predictive significance of the identified early and metastatic MPNST driver genes. Specifically, the expression changes related to the mutated genes, such as in RBMX, PNPLA6 and AGAP2, were associated with reduced patient survival, distinguishing them as potential prognostic biomarkers. This study underlines the relevance of integrating genomic results with clinical information for early diagnosis and prognostic understanding of tumor aggressiveness. [ABSTRACT FROM AUTHOR]

PURA syndrome is a congenital developmental disorder caused by de novo mutations in the PURA gene, which encodes a DNA/RNA-binding protein essential for transcriptional and translational regulation. We present the case of an 11-year-old patient with a de novo frameshift variant in the PURA gene, identified through whole exome sequencing (WES). In addition to the classical PURA deficiency phenotype, our patient exhibited pronounced sialorrhea and seizures, which were effectively treated with the ketogenic diet (KD). Our integrative approach, combining a literature review and bioinformatics data, has led to the first documented clinical case showing improvement in both sialorrhea and seizures with KD treatment, a phenomenon not previously reported. Although a direct relationship between the de novo PURA mutation and the KD was not established, we identified a novel frameshift deletion associated with a new clinical phenotype. [ABSTRACT FROM AUTHOR]

The purpose of this study was to assess the added diagnostic value of whole genome sequencing (WGS) for patients with inherited retinal diseases (IRDs) who remained undiagnosed after whole exome sequencing (WES). WGS was performed for index patients in 66 families. The datasets were analyzed according to GATK's guidelines. Additionally, DeepVariant was complemented by GATK's workflow, and a novel structural variant pipeline was developed. Overall, a molecular diagnosis was established in 19/66 (28.8%) index patients. Pathogenic deletions and one deep-intronic variant contributed to the diagnostic yield in 4/19 and 1/19 index patients, respectively. The remaining diagnoses (14/19) were attributed to exonic variants that were missed during WES analysis due to bioinformatic limitations, newly described loci, or unclear pathogenicity. The added diagnostic value of WGS equals 5/66 (9.6%) for our cohort, which is comparable to previous studies. This figure would decrease further to 1/66 (1.5%) with a standardized and reliable copy number variant workflow during WES analysis. Given the higher costs and limited added value, the implementation of WGS as a first-tier assay for inherited eye disorders in a diagnostic laboratory remains untimely. Instead, progress in bioinformatic tools and communication between diagnostic and clinical teams have the potential to ameliorate diagnostic yields. [ABSTRACT FROM AUTHOR]

The evolving landscape of clinical genetics is becoming increasingly relevant in the field of nephrology. HNF1B-associated renal disease presents with a diverse array of renal and extrarenal manifestations, prominently featuring cystic kidney disease and diabetes mellitus. For the genetic analyses, whole exome sequencing (WES) and multiplex ligation-dependent probe amplification (MLPA) were performed. Bioinformatics analysis was performed with Ingenuity Clinical Insights software (Qiagen). The patient's electronic record was utilized after receiving informed consent. In this report, we present seven cases of HNF1B-associated kidney disease, each featuring distinct genetic abnormalities and displaying diverse extrarenal manifestations. Over 12 years, the mean decline in eGFR averaged −2.22 ± 0.7 mL/min/1.73 m2. Diabetes mellitus was present in five patients, kidney dysplastic lesions in six patients, pancreatic dysplasia, hypomagnesemia and abnormal liver function tests in three patients each. This case series emphasizes the phenotypic variability and the fast decline in kidney function associated with HNF-1B-related disease. Additionally, it underscores that complex clinical presentations may have a retrospectively straightforward explanation through the use of diverse genetic analytical tools. [ABSTRACT FROM AUTHOR]

Genome-wide association studies significantly increased the number of hypertension risk variants; however, most of them focused on European societies. There is lack of such studies in developing countries, including Pakistan. The lack of research studies and the high prevalence of hypertension in the Pakistani community prompted us to design this study. Aldosterone synthase (CYP11B2) was thoroughly studied in different ethnic groups; however, no such study has been conducted in the Pashtun population of Khyber Pakhtunkhwa, Pakistan. In essential hypertension, the aldosterone synthase gene (CYP11B2) plays a significant role. Aldosterone synthesis is affected by both hereditary and environmental factors. Aldosterone synthase (encoded by the CYP11B2 gene) controls the conversion of deoxycorticosterone to aldosterone and, thus, has genetic influences. Polymorphisms in the CYP11B2 gene are linked to an increased risk of hypertension. Previous research on the polymorphism of the aldosterone synthase (CYP11B2) gene and its relationship to hypertension produced inconclusive results. The present study investigates the relationship between CYP11B2 gene polymorphism and hypertension in Pakistan's Pashtun population. We used the nascent exome sequencing method to identify variants associated with hypertension. The research was divided into two phases. In phase one, DNA samples from 200 adult hypertension patients (of age ≥ 30 years) and 200 controls were pooled (n = 200/pool) and subjected to Exome Sequencing. In the second phase, the WES reported SNPs were genotyped using the Mass ARRAY technique to verify and confirm the association between WES-identified SNPs and hypertension. WES identified a total of eight genetic variants in the CYP11B2 gene. The chi-square test and logistic regression analysis were used to estimate the minor allele frequencies (MAFs) and chosen SNPs relationships with hypertension. The frequency of minor allele T was found to be higher in cases compared to the control (42% vs. 30%: p = 0.001) for rs1799998 of CYP11B2 gene, while no significant results (p > 0.05) were observed for the remaining SNPs; rs4536, rs4537, rs4545, rs4543, rs4539, rs4546 and rs6418 showed no positive association with HTN in the studied population (all p > 0.05). Our study findings suggest that rs1799998 increases susceptibly to HTN in the Pashtun population of KP, Pakistan. [ABSTRACT FROM AUTHOR]

Incontinentia pigmenti (IP) is a rare skin disease combined with anomalies of the teeth, eyes, and central nervous system (CNS). Mutations of the IKBKG gene are responsible for IP. Among the most frequent CNS abnormalities found in IP using magnetic resonance imaging (MRI) are corpus callosum (CC) abnormalities. The aim of the study was to determine the presence of CC abnormalities, their relationship with the IKBKG mutations, and the possible presence of mutations of other genes. A group of seven IP patients was examined. Analyses of the IKBKG gene and the X-chromosome inactivation pattern were performed, as well as MRI and whole exome sequencing (WES) with the focus on the genes relevant for neurodegeneration. WES analysis showed IKBKG mutation in all examined patients. A patient who had a mutation of a gene other than IKBKG was excluded from further study. Four of the seven patients had clinically diagnosed CNS anomalies; two out of four had MRI-diagnosed CC anomalies. The simultaneous presence of IKBKG mutation and CC abnormalities and the absence of other mutations indicate that IKBKG may be the cause of CC abnormalities and should be included in the list of genes responsible for CC abnormalities. [ABSTRACT FROM AUTHOR]

Background: The reduction in next-generation sequencing (NGS) costs allows for using this method for newborn screening for monogenic diseases (MDs). In this report, we describe a clinical case of a newborn participating in the EXAMEN project (ClinicalTrials.gov Identifier: NCT05325749). Methods: The child presented with convulsive syndrome on the third day of life. Generalized convulsive seizures were accompanied by electroencephalographic patterns corresponding to epileptiform activity. Proband WES expanded to trio sequencing was performed. Results: A differential diagnosis was made between symptomatic (dysmetabolic, structural, infectious) neonatal seizures and benign neonatal seizures. There were no data in favor of the dysmetabolic, structural, or infectious nature of seizures. Molecular karyotyping and whole exome sequencing were not informative. Trio WES revealed a de novo variant in the KCNJ9 gene (1:160087612T > C, p.Phe326Ser, NM_004983), for which, according to the OMIM database, no association with the disease has been described to date. Three-dimensional modeling was used to predict the structure of the KCNJ9 protein using the known structure of its homologs. According to the predictions, Phe326Ser change possibly disrupts the hydrophobic contacts with the valine side chain. Destabilization of the neighboring structures may undermine the formation of GIRK2/GIRK3 tetramers necessary for their proper functioning. Conclusions: We believe that the identified variant may be the cause of the disease in this patient but further studies, including the search for other patients with the KCNJ9 variants, are needed. [ABSTRACT FROM AUTHOR]

Primary congenital hypothyroidism (CH) is a common neonatal endocrine disorder characterized by elevated concentrations of thyroid stimulating hormone (TSH) and low concentrations of free thyroxine (FT4). PAX8 and NKX2-1 are important transcription factors involved in thyroid development. In this study, we detected three novel variants in PAX8 (c.149A > C and c.329G > A) and NKX2-1 (c.706A > G) by whole exome sequencing (WES) in three unrelated CH patients with variable phenotypes. The results of Western blot and immunofluorescence analysis showed that the three variants had no effect on protein expression and subcellular localization. However, the results of the electrophoretic mobility shift assay (EMSA) and dual-luciferase reporter assay suggested that the three variants in PAX8 and NKX2-1 both affected their DNA-binding ability and reduced their transactivation capacity. Moreover, a dominant-negative effect in K236E−NKX2-1 was identified by dual-luciferase reporter assay. To sum up, our findings extend our knowledge of the current mutation spectrum of PAX8 and NKX2-1 and provide important information for diagnosing, treating, and preventing CH in these families. [ABSTRACT FROM AUTHOR]

Our aim was to analyze the phenotypic-genetic correlations in a patient diagnosed with early onset corticobasal syndrome with progressive non-fluent aphasia (CBS-PNFA), characterized by predominant apraxia of speech, accompanied by prominent right-sided upper-limb limb-kinetic apraxia, alien limb phenomenon, synkinesis, myoclonus, mild cortical sensory loss, and right-sided hemispatial neglect. Whole-exome sequencing (WES) identified rare single heterozygous variants in ATP7B (c.3207C>A), SORL1 (c.352G>A), SETX (c.2385_2387delAAA), and FOXP1 (c.1762G>A) genes. The functional analysis revealed that the deletion in the SETX gene changed the splicing pattern, which was accompanied by lower SETX mRNA levels in the patient's fibroblasts, suggesting loss-of-function as the underlying mechanism. In addition, the patient's fibroblasts demonstrated altered mitochondrial architecture with decreased connectivity, compared to the control individuals. This is the first association of the CBS-PNFA phenotype with the most common ATP7B pathogenic variant p.H1069Q, previously linked to Wilson's disease, and early onset Parkinson's disease. This study expands the complex clinical spectrum related to variants in well-known disease genes, such as ATP7B, SORL1, SETX, and FOXP1, corroborating the hypothesis of oligogenic inheritance. To date, the FOXP1 gene has been linked exclusively to neurodevelopmental speech disorders, while our study highlights its possible relevance for adult-onset progressive apraxia of speech, which guarantees further study. [ABSTRACT FROM AUTHOR]

Neurocytomas are rare low-grade brain tumors predominantly affecting young adults, but their cellular origin and molecular pathogenesis is largely unknown. We previously reported a sellar neurocytoma that secreted excess arginine vasopressin causing syndrome of inappropriate anti-diuretic hormone (SIADH). Whole exome sequencing in 21 neurocytoma tumor tissues identified somatic mutations in the plant homeodomain finger protein 14 (PHF14) in 3/21 (14%) tumors. Of these mutations, two were missense mutations and 4 caused splicing site losses, resulting in PHF14 dysfunction. Employing shRNA-mediated knockdown and CRISPR/Cas9-based knockout approaches, we demonstrated that loss of PHF14 increased proliferation and colony formation in five different human, mouse and rat mesenchymal and differentiated cell lines. Additionally, we demonstrated that PHF14 depletion resulted in upregulation of platelet derived growth factor receptor-alpha (PDGFRα) mRNA and protein in neuroblastoma SHSY-5Y cells and led to increased sensitivity to treatment with the PDGFR inhibitor Sunitinib. Furthermore, in a neurocytoma primary culture harboring splicing loss PHF14 mutations, overexpression of wild-type PHF14 and sunitinib treatment inhibited cell proliferation. Nude mice, inoculated with PHF14 knockout SHSY-5Y cells developed earlier and larger tumors than control cell-inoculated mice and Sunitinib administration caused greater tumor suppression in mice harboring PHF-14 knockout than control SHSY-5Y cells. Altogether our studies identified mutations of PHF14 in 14% of neurocytomas, demonstrate it can serve as an alternative pathway for certain cancerous behavior, and suggest a potential role for Sunitinib treatment in some patients with residual/recurrent neurocytoma. [ABSTRACT FROM AUTHOR]

Simple Summary: The micropapillary component in lung adenocarcinoma has a valid predictive role for patient prognosis. To investigate targeted clinical strategies, we collected 31 stage I lung adenocarcinoma samples and performed microdissection to separate micropapillary and non-micropapillary components, followed by whole exome sequencing. We focused on the genomic features, evolutionary advantages, and associated clinical implications of the micropapillary component, and we proposed and validated the possible presence of essential mutations, TP53 and ZNF469, in the formation of the micropapillary component. These features may be associated with the formation of high-grade invasive patterns in this component. Background: Micropapillary components are observed in a considerable proportion of ground-glass opacities (GGOs) and contribute to the poor prognosis of patients with invasive lung adenocarcinoma (LUAD). However, the underlying mutational processes related to the presence of micropapillary components remain obscure, limiting the development of clinical interventions. Methods: We collected 31 GGOs, which were separated into paired micropapillary and non-micropapillary components using microdissection. Whole-exome sequencing (WES) was performed on the GGO components, and bioinformatics analysis was conducted to reveal the genomic features of the micropapillary component in invasive LUAD. Results: The micropapillary component had more genomic variations, including tumor mutation burden, intratumoral heterogeneity, and copy number variation. We also observed the enrichment of AID/APOBEC mutation signatures and an increased activation of the RTK/Ras, Notch, and Wnt oncogenic pathways within the micropapillary component. A phylogenetic analysis further suggested that ERBB2/3/4, NCOR1/2, TP53, and ZNF469 contributed to the micropapillary component's progression during the early invasion of LUAD, a finding that was validated in the TCGA cohort. Conclusions: Our results revealed specific mutational characteristics of the micropapillary component of invasive LUAD in an Asian population. These characteristics were associated with the formation of high-grade invasive patterns. These preliminary findings demonstrated the potential of targeting the micropapillary component in patients with early-stage LUAD. [ABSTRACT FROM AUTHOR]

Cardiovascular disease (CVD) is one of the biggest threats to public health worldwide. Identifying key genetic contributors to CVD enables clinicians to assess the most effective treatment course and prognosis, as well as potentially inform family members. This often involves either whole exome sequencing (WES) or targeted panel analysis of known pathogenic genes. In the future, tailored or personalized therapeutic strategies may be implemented, such as gene therapy. With the recent revolution in deep sequencing technologies, we know that up to 90% of the human genome is transcribed, despite only 2% of the 6 billion DNA bases coding for proteins. The long non-coding RNA (lncRNA) "genes" make up an important and significant fraction of this "dark matter" of the genome. We highlight how, despite lncRNA genes exceeding that of classical protein-coding genes by number, the "non-coding" human genome is neglected when looking for genetic components of disease. WES platforms and pathogenic gene panels still do not cover even characterized lncRNA genes that are functionally involved in the pathophysiology of CVD. We suggest that the importance of lncRNAs in disease causation and progression be taken as seriously as that of pathogenic protein variants and mutations, and that this is maybe a new area of attention for clinical geneticists. [ABSTRACT FROM AUTHOR]

Immune checkpoint inhibitors, which unleash a patient's own T cells to kill tumors, are revolutionizing cancer treatment. Several independent studies suggest that higher non-synonymous mutational burden assessed by whole exome sequencing (WES) in tumors is associated with improved objective response, durable clinical benefit, and progression-free survival in immune checkpoint inhibitors treatment. Next-generation sequencing (NGS) is a promising technology being used in the clinic to direct patient treatment. Cancer genome WES poses a unique challenge due to tumor heterogeneity and sequencing artifacts introduced by formalin-fixed, paraffin-embedded (FFPE) tissue. In order to evaluate the data interoperability of WES data from different sources to survey tumor mutational landscape, we compared WES data of several tumor/normal matched samples from five commercial vendors. A large data discrepancy was observed from vendors' self-reported data. Independent data analysis from vendors' raw NGS data shows that whole exome sequencing data from qualified vendors can be combined and analyzed uniformly to derive comparable quantitative estimates of tumor mutational burden. [ABSTRACT FROM AUTHOR]

Simple Summary: Copy-number variations (CNVs) have important clinical implications for several diseases and cancers. We reviewed 50 popular CNV calling tools and included 11 tools for benchmarking in a reference cohort encompassing 39 whole genome sequencing (WGS) samples paired current clinical standard—SNP-array based CNV calling. For nine samples we also performed whole exome sequencing (WES), to address the effect of sequencing protocol on CNV calling. Furthermore, we included Gold Standard reference sample NA12878, and tested 12 samples with CNVs confirmed by multiplex ligation-dependent probe amplification (MLPA). Tool performance varied greatly in the number of called CNVs and bias for CNV lengths. Some tools had near-perfect recall of CNVs from arrays for some samples, but poor precision. We suggest combining the best tools also based on different methodologies: GATK gCNV, Lumpy, DELLY, and cn.MOPS. Copy-number variations (CNVs) have important clinical implications for several diseases and cancers. Relevant CNVs are hard to detect because common structural variations define large parts of the human genome. CNV calling from short-read sequencing would allow single protocol full genomic profiling. We reviewed 50 popular CNV calling tools and included 11 tools for benchmarking in a reference cohort encompassing 39 whole genome sequencing (WGS) samples paired current clinical standard—SNP-array based CNV calling. Additionally, for nine samples we also performed whole exome sequencing (WES), to address the effect of sequencing protocol on CNV calling. Furthermore, we included Gold Standard reference sample NA12878, and tested 12 samples with CNVs confirmed by multiplex ligation-dependent probe amplification (MLPA). Tool performance varied greatly in the number of called CNVs and bias for CNV lengths. Some tools had near-perfect recall of CNVs from arrays for some samples, but poor precision. Several tools had better performance for NA12878, which could be a result of overfitting. We suggest combining the best tools also based on different methodologies: GATK gCNV, Lumpy, DELLY, and cn.MOPS. Reducing the total number of called variants could potentially be assisted by the use of background panels for filtering of frequently called variants. [ABSTRACT FROM AUTHOR]

Mutations in the HSD17B3 gene cause HSD17B3 deficiency and result in 46, XY Disorders of Sex Development (46, XY DSD). The diagnosis of 46, XY DSD is very challenging and not rarely is confirmed only at older ages, when an affected XY female presents with primary amenorrhea or develops progressive virilization. The patient described in this paper represents a case of discrepancies between non-invasive prenatal testing (NIPT) and ultrasound based fetal sex determination detected during prenatal screening. Exome sequencing was performed on the cell free fetal DNA (cffDNA), amniotic fluid, and the parents. Libraries were generated according to the manufacturer's protocols using TruSight One Kits (Illumina Inc., San Diego, CA, USA). Sequencing was carried out on NEXT Seq 500 (Illumina) to mean sequencing depth of at least 100×. A panel of sexual disease genes was used in order to search for a causative variant. The finding of a mutation (c.645 A>T, p.Glu215Asp) in HSD17B3 gene in amniotic fluid as well as in cffDNA and both parents supported the hypothesis of the HSD17B3 deficiency. In conclusion, we used clinical exome sequencing and non-invasive prenatal detection, providing a solution for NIPT of a single-gene disorder. Early genetic diagnoses are useful for patients and clinicians, contribute to clinical knowledge of DSD, and are invaluable for genetic counseling of couples contemplating future pregnancies. [ABSTRACT FROM AUTHOR]

Synchronous primary malignancies occur in a small proportion of head and neck squamous cell carcinoma (HNSCC) patients. Here, we analysed three synchronous primaries and a recurrence from one patient by comparing the genomic and transcriptomic profiles among the tumour samples and determining the recurrence origin. We found remarkable levels of heterogeneity among the primary tumours, and through the patterns of shared mutations, we traced the origin of the recurrence. Interestingly, the patient carried germline variants that might have predisposed him to carcinogenesis, together with a history of alcohol and tobacco consumption. The mutational signature analysis confirmed the impact of alcohol exposure, with Signature 16 present in all tumour samples. Characterisation of immune cell infiltration highlighted an immunosuppressive environment in all samples, which exceeded the potential activity of T cells. Studies such as the one described here have important clinical value and contribute to personalised treatment decisions for patients with synchronous primaries and matched recurrences. [ABSTRACT FROM AUTHOR]

In 1959, 63 years after the death of John Langdon Down, Jérôme Lejeune discovered trisomy 21 as the genetic reason for Down syndrome. Screening for Down syndrome has been applied since the 1960s by using maternal age as the risk parameter. Since then, several advances have been made. First trimester screening, combining maternal age, maternal serum parameters and ultrasound findings, emerged in the 1990s with a detection rate (DR) of around 90–95% and a false positive rate (FPR) of around 5%, also looking for trisomy 13 and 18. With the development of high-resolution ultrasound, around 50% of fetal anomalies are now detected in the first trimester. Non-invasive prenatal testing (NIPT) for trisomy 21, 13 and 18 is a highly efficient screening method and has been applied as a first-line or a contingent screening approach all over the world since 2012, in some countries without a systematic screening program. Concomitant with the rise in technology, the possibility of screening for other genetic conditions by analysis of cfDNA, such as sex chromosome anomalies (SCAs), rare autosomal anomalies (RATs) and microdeletions and duplications, is offered by different providers to an often not preselected population of pregnant women. Most of the research in the field is done by commercial providers, and some of the tests are on the market without validated data on test performance. This raises difficulties in the counseling process and makes it nearly impossible to obtain informed consent. In parallel with the advent of new screening technologies, an expansion of diagnostic methods has begun to be applied after invasive procedures. The karyotype has been the gold standard for decades. Chromosomal microarrays (CMAs) able to detect deletions and duplications on a submicroscopic level have replaced the conventional karyotyping in many countries. Sequencing methods such as whole exome sequencing (WES) and whole genome sequencing (WGS) tremendously amplify the diagnostic yield in fetuses with ultrasound anomalies. [ABSTRACT FROM AUTHOR]

Corneal dystrophies (CDs) are a diverse group of inherited disorders with a heterogeneous genetic background. Here, we report the identification of a novel ZNF143 heterozygous missense mutation in three individuals of the same family with clinical and pathological features that are consistent with endothelial CD. Ophthalmologic examination revealed diffuse corneal clouding and edema with decreased endothelial cell density. Pathological findings showed increased corneal thickness due to edema of basal epithelial cells and stroma, and abnormal metaplastic endothelium with stratified epithelium-like changes. Patients' metaplastic corneal endothelial cells expressed predominantly cytokerain 7, cytokeratin 19, and E-cadherin. Although Sanger sequencing did not detect any mutation associated with endothelial CDs, whole exome sequencing identified the ZNF143 c.937G>C p.(Asp313His) mutation as a candidate gene for our patients' endothelial CD. In-vitro functional studies demonstrated that mutant ZNF143 promoted the mesenchymal-to-epithelial transition; it upregulated the expression of genes associated with epithelialization in human corneal endothelial cells. Additionally, proinflammatory cytokine responsive genes were significantly enriched after mutant ZNF143 transfection, which may contribute to the severe phenotype of the three patients. These findings link a mutation in ZNF143 with endothelial CD for the first time. [ABSTRACT FROM AUTHOR]

Immune checkpoint inhibitors, which unleash a patient's own T cells to kill tumors, are revolutionizing cancer treatment. Several independent studies suggest that higher non-synonymous mutational burden assessed by whole exome sequencing (WES) in tumors is associated with improved objective response, durable clinical benefit, and progression-free survival in immune checkpoint inhibitors treatment. Next-generation sequencing (NGS) is a promising technology being used in the clinic to direct patient treatment. Cancer genome WES poses a unique challenge due to tumor heterogeneity and sequencing artifacts introduced by formalin-fixed, paraffin-embedded (FFPE) tissue. In order to evaluate the data interoperability of WES data from different sources to survey tumor mutational landscape, we compared WES data of several tumor/normal matched samples from five commercial vendors. A large data discrepancy was observed from vendors' self-reported data. Independent data analysis from vendors' raw NGS data shows that whole exome sequencing data from qualified vendors can be combined and analyzed uniformly to derive comparable quantitative estimates of tumor mutational burden.