Place of concordance-discordance model in evaluating NGS performance.

INTRODUCTION: Ideally, evaluating NGS performance requires a gold standard; in its absence, concordance between replicates is often used as substitute standard. However, the appropriateness of the concordance-discordance criterion has been rarely evaluated. This study analyzes the relationship between the probability of discordance and the probability of error under different conditions. METHODS: This study used a conditional probability approach under conditional dependence then conditional independence between two sequencing results and compares the probabilities of discordance and error in different theoretical conditions of sensitivity, specificity, and correlation between replicates, then on real results of sequencing genome NA12878. The study examines also covariate effects on discordance and error using generalized additive models with smooth functions. RESULTS: With 99% sensitivity and 99.9% specificity under conditional independence, the probability of error for a positive concordant pair of calls is 0.1%. With additional hypotheses of 0.1% prevalence and 0.9 correlation between replicates, the probability of error for a positive concordant pair is 47.4%. With real data, the estimated sensitivity, specificity, and correlation between tests for variants are around 98.98%, 99.996%, and 93%, respectively, and the error rate for positive concordant calls approximates 2.5%. In covariate effect analyses, the effects' functional form are close between discordance and error models, though the parts of deviance explained by the covariates differ between discordance and error models. CONCLUSION: With conditional independence of two sequencing results, the concordance-discordance criterion seems acceptable as substitute standard. However, with high correlation, the criterion becomes questionable because a high percentage of false concordant results appears among concordant results.

Familial transmission of chromoanagenesis leads to unpredictable unbalanced rearrangements through meiotic recombination.

Chromoanagenesis is a cellular mechanism that leads to complex chromosomal rearrangements (CCR) during a single catastrophic event. It may result in loss and/or gain of genetic material and may be responsible for various phenotypes. These rearrangements are usually sporadic. However, some familial cases have been reported. Here, we studied six families in whom an asymptomatic or paucisymptomatic parent transmitted a CCR to its offspring in an unbalanced manner. The rearrangements were characterized by karyotyping, fluorescent in situ hybridization, chromosomal microarray (CMA) and/or whole genome sequencing (WGS) in the carrier parents and offspring. We then hypothesized meiosis-pairing figures between normal and abnormal parental chromosomes that may have led to the formation of new unbalanced rearrangements through meiotic recombination. Our work indicates that chromoanagenesis might be associated with a normal phenotype and normal fertility, even in males, and that WGS may be the only way to identify these events when there is no imbalance. Subsequently, the CCR can be transmitted to the next generation in an unbalanced and unpredictable manner following meiotic recombination. Thereby, prenatal diagnosis using CMA should be proposed to these families to detect any pathogenic imbalances in the offspring.

Complete characterisation of two new large Xq28 duplications involving F8 using whole genome sequencing in patients without haemophilia A.

INTRODUCTION: Depending on the location of insertion of the gained region, F8 duplications can have variable clinical impacts from benign impact to severe haemophilia A phenotype. AIM: To characterize two large Xq28 duplications involving F8 incidentally detected by chromosome microarray analysis (CMA) in two patients presenting severe intellectual disability but no history of bleeding disorder. METHODS: Whole genome sequencing (WGS) was performed in order to characterize the two large Xq28 duplications at nucleotide level. RESULTS: In patient 1, a 60-73 kb gained region encompassing the exons 23-26 of F8 and SMIM9 was inserted at the int22h-2 locus following a non-homologous recombination between int22h-1 and int22h-2. We hypothesized that two independent events, micro-homology-mediated break-induced replication (MMBIR) and break-induced replication (BIR), could be involved in this rearrangement. In patient 2, the CMA found duplication from 101 to 116-kb long encompassing the exons 16-26 of F8 and SMIM9. The WGS analysis identified a more complex rearrangement with the presence of three genomic junctions. Due to the multiple micro-homologies observed at breakpoints, a replication-based mechanism such as fork stalling and template switching (FoSTeS) was greatly suspected. In both cases, these complex rearrangements preserved an intact copy of the F8. CONCLUSION: This study highlights the value of WGS to characterize the genomic junction at the nucleotide level and ultimately better describe the molecular mechanisms involved in Xq28 structural variations. It also emphasizes the importance of specifying the structure of the genomic gain in order to improve genotype-phenotype correlation and genetic counselling.

Whole MYBPC3 NGS sequencing as a molecular strategy to improve the efficiency of molecular diagnosis of patients with hypertrophic cardiomyopathy.

Hypertrophic cardiomyopathy (HCM) is the most common heritable cardiomyopathy, historically believed to affect 1 of 500 people. MYBPC3 pathogenic variations are the most frequent cause of familial HCM and more than 90% of them introduce a premature termination codon. The current study aims to determine the prevalence of deep intronic MYBPC3 pathogenic variations that could lead to splice mutations. To improve molecular diagnosis, a next-generation sequencing (NGS) workflow based on whole MYBPC3 sequencing of a cohort of 93 HCM patients, for whom no putatively causative point mutations were identified after NGS sequencing of a panel of 48 cardiomyopathy-causing genes, was performed. Our approach led us to reconsider the molecular diagnosis of six patients of the cohort (6.5%). These HCM probands were carriers of either a new large MYBPC3 rearrangement or splice intronic variations (five cases). Four pathogenic intronic variations, including three novel ones, were detected. Among them, the prevalence of one of them (NM_000256.3:c.1927+ 600 C>T) was estimated at about 0.35% by the screening of 1,040 unrelated HCM individuals. This study suggests that deep MYBPC3 splice mutations account for a significant proportion of HCM cases (6.5% of this cohort). Consequently, NGS sequencing of MYBPC3 intronic sequences have to be performed systematically.

Development of a new expanded next-generation sequencing panel for genetic diseases involved in dyslipidemia.

The aim of this study was to provide an efficient tool: reliable, able to increase the molecular diagnosis performance, to facilitate the detection of copy number variants (CNV), to assess genetic risk scores (wGRS) and to offer the opportunity to explore candidate genes. Custom SeqCap EZ libraries, NextSeq500 sequencing and a homemade pipeline enable the analysis of 311 dyslipidemia-related genes. In the training group (48 DNA from patients with a well-established molecular diagnosis), this next-generation sequencing (NGS) workflow showed an analytical sensitivity >99% (n = 532 variants) without any false negative including a partial deletion of one exon. In the prospective group, from 25 DNA from patients without prior molecular analyses, 18 rare variants were identified in the first intention panel genes, allowing the diagnosis of monogenic dyslipidemia in 11 patients. In six other patients, the analysis of minor genes and wGRS determination provided a hypothesis to explain the dyslipidemia. Remaining data from the whole NGS workflow identified four patients with potentially deleterious variants. This NGS process gives a major opportunity to accede to an enhanced understanding of the genetic of dyslipidemia by simultaneous assessment of multiple genetic determinants.

Whole genome paired-end sequencing elucidates functional and phenotypic consequences of balanced chromosomal rearrangement in patients with developmental disorders.

BACKGROUND: Balanced chromosomal rearrangements associated with abnormal phenotype are rare events, but may be challenging for genetic counselling, since molecular characterisation of breakpoints is not performed routinely. We used next-generation sequencing to characterise breakpoints of balanced chromosomal rearrangements at the molecular level in patients with intellectual disability and/or congenital anomalies. METHODS: Breakpoints were characterised by a paired-end low depth whole genome sequencing (WGS) strategy and validated by Sanger sequencing. Expression study of disrupted and neighbouring genes was performed by RT-qPCR from blood or lymphoblastoid cell line RNA. RESULTS: Among the 55 patients included (41 reciprocal translocations, 4 inversions, 2 insertions and 8 complex chromosomal rearrangements), we were able to detect 89% of chromosomal rearrangements (49/55). Molecular signatures at the breakpoints suggested that DNA breaks arose randomly and that there was no major influence of repeated elements. Non-homologous end-joining appeared as the main mechanism of repair (55% of rearrangements). A diagnosis could be established in 22/49 patients (44.8%), 15 by gene disruption (KANSL1, FOXP1, SPRED1, TLK2, MBD5, DMD, AUTS2, MEIS2, MEF2C, NRXN1, NFIX, SYNGAP1, GHR, ZMIZ1) and 7 by position effect (DLX5, MEF2C, BCL11B, SATB2, ZMIZ1). In addition, 16 new candidate genes were identified. Systematic gene expression studies further supported these results. We also showed the contribution of topologically associated domain maps to WGS data interpretation. CONCLUSION: Paired-end WGS is a valid strategy and may be used for structural variation characterisation in a clinical setting.

Identification of mobile retrocopies during genetic testing: Consequences for routine diagnosis.

Human retrocopies, that is messenger RNA transcripts benefitting from the long interspersed element 1 machinery for retrotransposition, may have specific consequences for genomic testing. Next genetration sequencing (NGS) techniques allow the detection of such mobile elements but they may be misinterpreted as genomic duplications or be totally overlooked. We report eight observations of retrocopies detected during diagnostic NGS analyses of targeted gene panels, exome, or genome sequencing. For seven cases, while an exons-only copy number gain was called, read alignment inspection revealed a depth of coverage shift at every exon-intron junction where indels were also systematically called. Moreover, aberrant chimeric read pairs spanned entire introns or were paired with another locus for terminal exons. The 8th retrocopy was present in the reference genome and thus showed a normal NGS profile. We emphasize the existence of retrocopies and strategies to accurately detect them at a glance during genetic testing and discuss pitfalls for genetic testing.

The epilepsy phenotypic spectrum associated with a recurrent CUX2 variant.

OBJECTIVE: Cut homeodomain transcription factor CUX2 plays an important role in dendrite branching, spine development, and synapse formation in layer II to III neurons of the cerebral cortex. We identify a recurrent de novo CUX2 p.Glu590Lys as a novel genetic cause for developmental and epileptic encephalopathy (DEE). METHODS: The de novo p.Glu590Lys variant was identified by whole-exome sequencing (n = 5) or targeted gene panel (n = 4). We performed electroclinical and imaging phenotyping on all patients. RESULTS: The cohort comprised 7 males and 2 females. Mean age at study was 13 years (0.5-21.0). Median age at seizure onset was 6 months (2 months to 9 years). Seizure types at onset were myoclonic, atypical absence with myoclonic components, and focal seizures. Epileptiform activity on electroencephalogram was seen in 8 cases: generalized polyspike-wave (6) or multifocal discharges (2). Seizures were drug resistant in 7 or controlled with valproate (2). Six patients had a DEE: myoclonic DEE (3), Lennox-Gastaut syndrome (2), and West syndrome (1). Two had a static encephalopathy and genetic generalized epilepsy, including absence epilepsy in 1. One infant had multifocal epilepsy. Eight had severe cognitive impairment, with autistic features in 6. The p.Glu590Lys variant affects a highly conserved glutamine residue in the CUT domain predicted to interfere with CUX2 binding to DNA targets during neuronal development. INTERPRETATION: Patients with CUX2 p.Glu590Lys display a distinctive phenotypic spectrum, which is predominantly generalized epilepsy, with infantile-onset myoclonic DEE at the severe end and generalized epilepsy with severe static developmental encephalopathy at the milder end of the spectrum. Ann Neurol 2018;83:926-934.

Structure-activity relationship study: Mechanism of cyto-genotoxicity of Nitropyrazole-derived high energy density materials family.

Stringent toxicological tests have to be performed prior to the industrial development of alternative chemicals particularly high energy dense materials (HEDMs) such as explosives. The properties (e.g., power, stability) of these compounds are constantly being improved, the current axis of research being the nitration of nitrogen heterocycles leading to HEDMs such as nitropyrazole-derived molecules. However, except for 3,4,5-trinitropyrazole (3,4,5-TNP), which was shown to be highly toxic in mice, the toxicological impact of these HEDMs has so far not been investigated. Furthermore, as industrials are strongly advised to develop alternative safety testing assays to in vivo experiments, we herein focused on determining the cytotoxic and genotoxic effects of seven Nitropyrazole-derived HEDMs on three rodent cell lines (mouse embryonic BALB/3T3 clone A31 cells, Chinese hamster ovary cells CHO-K1 and mouse lymphoma L5178Y TK +/- clone (3.7.2C) cells), two human fibroblast lines (CRC05, PFS04062) and on the human hepatic HepaRG model (both in proliferative and differentiated cells). A stronger cytotoxic effect was observed for 1,3-dinitropyrazole (1, 3-DNP) and 3,4,5-TNP in all cell lines, though differentiated HepaRG cells clearly displayed fewer likely due to the metabolism and elimination of these molecules by their functional biotransformation pathways. At the mechanistic level, the sub-chronic cytotoxic and genotoxic effects were linked to ROS/RNS production (experimental assays), HA2.X and to transcriptomic data highlighting the increase in DNA repair mechanisms.

Transcriptional regulation of CRMP5 controls neurite outgrowth through Sox5.

Transcriptional regulation of proteins involved in neuronal polarity is a key process that underlies the ability of neurons to transfer information in the central nervous system. The Collapsin Response Mediator Protein (CRMP) family is best known for its role in neurite outgrowth regulation conducting to neuronal polarity and axonal guidance, including CRMP5 that drives dendrite differentiation. Although CRMP5 is able to control dendritic development, the regulation of its expression remains poorly understood. Here we identify a Sox5 consensus binding sequence in the putative promoter sequence upstream of the CRMP5 gene. By luciferase assays we show that Sox5 increases CRMP5 promoter activity, but not if the putative Sox5 binding site is mutated. We demonstrate that Sox5 can physically bind to the CRMP5 promoter DNA in gel mobility shift and chromatin immunoprecipitation assays. Using a combination of real-time RT-PCR and quantitative immunocytochemistry, we provide further evidence for a Sox5-dependent upregulation of CRMP5 transcription and protein expression in N1E115 cells: a commonly used cell line model for neuronal differentiation. Furthermore, we report that increasing Sox5 levels in this neuronal cell line inhibits neurite outgrowth. This inhibition requires CRMP5 because CRMP5 knockdown prevents the Sox5-dependent effect. We confirm the physiological relevance of the Sox5-CRMP5 pathway in the regulation of neurite outgrowth using mouse primary hippocampal neurons. These findings identify Sox5 as a critical modulator of neurite outgrowth through the selective activation of CRMP5 expression.

Centralization errors in comparative genomic hybridization array analysis of pituitary tumor samples.

Reliable interpretation of comparative genomic hybridization array (aCGH) results requires centralization and normalization of the data. We evaluated the reliability of aCGH centralization by comparing aCGH results (with classical centralization-normalization steps) to fluorescence in situ hybridization (FISH) results. In addition, we propose a method to correct centralization bias. Sixty-six pituitary tumors were analyzed (Agilent aCGH + SNP 4 × 180K microarray). For each tumor, the FISH-based log2 (ratios) of a subset of chromosomes were compared with the corresponding aCGH raw log2 (ratios). With our new normalization-centralization process, this difference was added to all log2 (ratios), before performing loess regression on non-altered probes only. Finally, the mean log2 (ratio) and the percentage of normal probes were compared between CGHnormaliter and our new FISH-based method. For 11 tumors, FISH results and raw CGH log2 (ratios) differed significantly. In addition, nine tumors showed discrepancies between results generated by CGHnormaliter and our new-method. Such discrepancies seemed to occur with tumours with many abnormalities (0%-40% normal probes), rather than in those tumours with fewer abnormalities (31%-100% normal probes). Five tumors had too few normal probes to allow normalization. In these tumors, which can exhibit many changes in DNA copy number, we found that centralization bias was frequent and uncorrected by current normalization methods. Therefore, an external control for centralization, such as FISH analysis, is required to insure reliable interpretation of aCGH data.

Toxicokinetics and tolerance of a high energy material 3,4,5-trinitropyrazole (TNP) in mice.

The high-energy compound 3,4,5-trinitropyrazole (TNP) was developed as an alternative to other less energetic and more sensitive explosive materials, in particular 1-methyl-2,4,6-trinitrobenzene (TNT). However, the level of toxicity of TNP remains understudied. Here using an in vivo CD1 mouse model, we mimicked an acute exposure (24 h) to TNP, given either orally or intravenously, and determined the maximum administrable doses (190 mg/kg and 11 mg/kg, respectively), as well as the lethal dose for 50% (LD50) of female or male mice (390 mg/kg for both) treated intravenously with TNP alone. Several metabolites including nitroso-dinitro-pyrazole, hydroxylamino-dinitro-pyrazole, hydroxyl-dinitro-pyrazole and amino-dinitro-pyrazole were identified in urine. TNP is quickly metabolized and eliminated via urine as two main amino-dinitro-pyrazole metabolites. A comparison of the transcriptomic effects of TNP and TNT after 10 days exposure enabled us to demonstrate no major induction of transcripts involved both in cell death mechanisms (apoptosis, necrosis, autophagy) and physiological pathways (glycolysis, ATP production). Finally, subchronic exposure to TNP was replicated in female mice, fed 16.8-52.8 mg/kg/day of TNP for one month, to study the impact on cellular functions. Although blood TNP levels remained high, a lower rate of TNP accumulation in the liver and lungs were observed than during an acute exposure. Conversely, cellular stress functions explored using the RT2 Profiler™ PCR Array Mouse Molecular Toxicology PathwayFinder remained unaltered after this chronic exposure. These findings demonstrate that TNP can be rapidly eliminated in vivo without accumulating in vital organs.

What Does This Mutation Mean? The Tools and Pitfalls of Variant Interpretation in Lymphoid Malignancies.

High throughput sequencing (HTS) is increasingly important in determining cancer diagnoses, with subsequent prognostic and therapeutic implications. The biology of cancer is becoming increasingly deciphered and it is clear that therapy needs to be individually tailored. Whilst translational research plays an important role in lymphoid malignancies, few guidelines exist to guide biologists and routine laboratories through this constantly evolving field. In this article, we review the challenges of interpreting HTS in lymphoid malignancies and provide a toolkit to interpret single nucleotide variants obtained from HTS. We define the pre-analytical issues such as sequencing DNA obtained from formalin-fixed and paraffin-embedded tissue (FFPE), the acquisition of germline DNA, or the bioinformatic pitfalls, the analytical issues encountered and how to manage them. We describe the main constitutional and cancer databases, their characteristics and limitations, with an emphasis on variant interpretation in lymphoid malignancies. Finally, we discuss the challenges of predictions that one can make using in silico or in vitro modelling, pharmacogenomic screening, and the limits of those prediction tools. This description of the current status in genomic interpretation highlights the need for new large databases and international collaboration in the lymphoma field.

An original phylogenetic approach identified mitochondrial haplogroup T1a1 as inversely associated with breast cancer risk in BRCA2 mutation carriers.

INTRODUCTION: Individuals carrying pathogenic mutations in the BRCA1 and BRCA2 genes have a high lifetime risk of breast cancer. BRCA1 and BRCA2 are involved in DNA double-strand break repair, DNA alterations that can be caused by exposure to reactive oxygen species, a main source of which are mitochondria. Mitochondrial genome variations affect electron transport chain efficiency and reactive oxygen species production. Individuals with different mitochondrial haplogroups differ in their metabolism and sensitivity to oxidative stress. Variability in mitochondrial genetic background can alter reactive oxygen species production, leading to cancer risk. In the present study, we tested the hypothesis that mitochondrial haplogroups modify breast cancer risk in BRCA1/2 mutation carriers. METHODS: We genotyped 22,214 (11,421 affected, 10,793 unaffected) mutation carriers belonging to the Consortium of Investigators of Modifiers of BRCA1/2 for 129 mitochondrial polymorphisms using the iCOGS array. Haplogroup inference and association detection were performed using a phylogenetic approach. ALTree was applied to explore the reference mitochondrial evolutionary tree and detect subclades enriched in affected or unaffected individuals. RESULTS: We discovered that subclade T1a1 was depleted in affected BRCA2 mutation carriers compared with the rest of clade T (hazard ratio (HR) = 0.55; 95% confidence interval (CI), 0.34 to 0.88; P = 0.01). Compared with the most frequent haplogroup in the general population (that is, H and T clades), the T1a1 haplogroup has a HR of 0.62 (95% CI, 0.40 to 0.95; P = 0.03). We also identified three potential susceptibility loci, including G13708A/rs28359178, which has demonstrated an inverse association with familial breast cancer risk. CONCLUSIONS: This study illustrates how original approaches such as the phylogeny-based method we used can empower classical molecular epidemiological studies aimed at identifying association or risk modification effects.

Prospective cardiovascular events in patients with advanced thoracic cancer treated with immune checkpoint inhibitor.

INTRODUCTION: Myocarditis is the most lethal cardiovascular immune related adverse events with a low incidence, depending on the studies. We prospectively studied the potential interest of a systematic screening to early detect immune related myocarditis and confirm the incidence of immune-induced myocarditis in advanced lung cancer and the impact of troponin systematic screening in early detection of other major cardiovascular events (MACE). MATERIAL AND METHODS: This prospective bicentric study includes adults who received at least one dose of immune checkpoint inhibitor (ICI) for advanced lung cancer. Cardiac biomarkers dosage, ECG and transthoracic echography (TTE) were done at baseline. Diagnosis of myocarditis was based on European Society of Cardiology recommendations. MACEs were reported during the observation period. RESULTS: Among 298 patients, 5 (1.68 %) immune-induced myocarditis occurred, all being asymptomatic with at first troponin elevation, treated by corticosteroids and ICI's discontinuation. No attributable death occurred, and no specific clinical characteristics were identified with myocarditis onset. Three patients were rechallenged with ICI after troponin normalization in the absence of other therapeutic options. Recurrence occurred in 2 patients, with a re-increase of troponin and a de novo modification of the ECG. Systematic cardiovascular screening also led to 14 cardiovascular diseases detection and 11 MACEs during ICI. CONCLUSION: Systematic cardiovascular screening has uncovered slightly more immuno-induced myocarditis cases than reported previously, but without altering treatment strategies due to their subclinical nature. Additionally, it helps detecting other cardiovascular diseases in this comorbid population.

Performance comparisons between clustering models for reconstructing NGS results from technical replicates.

To improve the performance of individual DNA sequencing results, researchers often use replicates from the same individual and various statistical clustering models to reconstruct a high-performance callset. Here, three technical replicates of genome NA12878 were considered and five model types were compared (consensus, latent class, Gaussian mixture, Kamila-adapted k-means, and random forest) regarding four performance indicators: sensitivity, precision, accuracy, and F1-score. In comparison with no use of a combination model, i) the consensus model improved precision by 0.1%; ii) the latent class model brought 1% precision improvement (97%-98%) without compromising sensitivity (= 98.9%); iii) the Gaussian mixture model and random forest provided callsets with higher precisions (both >99%) but lower sensitivities; iv) Kamila increased precision (>99%) and kept a high sensitivity (98.8%); it showed the best overall performance. According to precision and F1-score indicators, the compared non-supervised clustering models that combine multiple callsets are able to improve sequencing performance vs. previously used supervised models. Among the models compared, the Gaussian mixture model and Kamila offered non-negligible precision and F1-score improvements. These models may be thus recommended for callset reconstruction (from either biological or technical replicates) for diagnostic or precision medicine purposes.

Direct Comparative Analysis of a Pharmacogenomics Panel with PacBio Hifi® Long-Read and Illumina Short-Read Sequencing.

BACKGROUND: Pharmacogenetics (PGx) aims to determine genetic signatures that can be used in clinical settings to individualize treatment for each patient, including anti-cancer drugs, anti-psychotics, and painkillers. Taken together, a better understanding of the impacts of genetic variants on the corresponding protein function or expression permits the prediction of the pharmacological response: responders, non-responders, and those with adverse drug reactions (ADRs). OBJECTIVE: This work provides a comparison between innovative long-read sequencing (LRS) and short-read sequencing (SRS) techniques. METHODS AND MATERIALS: The gene panel captured using PacBio HiFi® sequencing was tested on thirteen clinical samples on GENTYANE's platform. SRS, using a comprehensive pharmacogenetics panel, was performed in routine settings at the Civil Hospitals of Lyon. We focused on complex regions analysis, including copy number variations (CNVs), structural variants, repeated regions, and phasing-haplotyping for three key pharmacogenes: CYP2D6, UGT1A1, and NAT2. RESULTS: Variants and the corresponding expected star (*) alleles were reported. Although only 38.4% concordance was found for haplotype determination and 61.5% for diplotype, this did not affect the metabolism scoring. A better accuracy of LRS was obtained for the detection of the CYP2D6*5 haplotype in the presence of the duplicated wild-type CYP2D6*2 form. A total concordance was performed for UGT1A1 TA repeat detection. Direct phasing using the LRS approach allowed us to correct certain NAT2 profiles. CONCLUSIONS: Combining an optimized variant-calling pipeline and with direct phasing analysis, LRS is a robust technique for PGx analysis that can minimize the risk of mis-haplotyping.

Paired Comparison of Routine Molecular Screening of Patient Samples with Advanced Non-Small Cell Lung Cancer in Circulating Cell-Free DNA Using Three Targeted Assays.

INTRODUCTION: Progressive advanced non-small cell lung cancer (NSCLC) accounts for about 80-85% of all lung cancers. Approximately 10-50% of patients with NSCLC harbor targetable activating mutations, such as in-frame deletions in Exon 19 (Ex19del) of EGFR. Currently, for patients with advanced NSCLC, testing for sensitizing mutations in EGFR is mandatory prior to the administration of tyrosine kinase inhibitors. PATIENTS AND METHODS: Plasma was collected from patients with NSCLC. We carried out targeted NGS using the Plasma-SeqSensei™ SOLID CANCER IVD kit on cfDNA (circulating free DNA). Clinical concordance for plasma detection of known oncogenic drivers was reported. In a subset of cases, validation was carried out using an orthogonal OncoBEAMTM EGFR V2 assay, as well as with our custom validated NGS assay. Somatic alterations were filtered, removing somatic mutations attributable to clonal hematopoiesis for our custom validated NGS assay. RESULTS: In the plasma samples, driver targetable mutations were studied, with a mutant allele frequency (MAF) ranging from 0.00% (negative detection) to 82.25%, using the targeted next-generation sequencing Plasma-SeqSensei™ SOLID CANCER IVD Kit. In comparison with the OncoBEAMTM EGFR V2 kit, the EGFR concordance is 89.16% (based on the common genomic regions). The sensitivity and specificity rates based on the genomic regions (EGFR exons 18, 19, 20, and 21) were 84.62% and 94.67%. Furthermore, the observed clinical genomic discordances were present in 25% of the samples: 5% in those linked to the lower of coverage of the OncoBEAMTM EGFR V2 kit, 7% in those induced by the sensitivity limit on the EGFR with the Plasma-SeqSensei™ SOLID CANCER IVD Kit, and 13% in the samples linked to the larger KRAS, PIK3CA, BRAF coverage of the Plasma-SeqSensei™ SOLID CANCER IVD kit. Most of these somatic alterations were cross validated in our orthogonal custom validated NGS assay, used in the routine management of patients. The concordance is 82.19% in the common genomic regions (EGFR exons 18, 19, 20, 21; KRAS exons 2, 3, 4; BRAF exons 11, 15; and PIK3CA exons 10, 21). The sensitivity and specificity rates were 89.38% and 76.12%, respectively. The 32% of genomic discordances were composed of 5% caused by the limit of coverage of the Plasma-SeqSensei™ SOLID CANCER IVD kit, 11% induced by the sensitivity limit of our custom validated NGS assay, and 16% linked to the additional oncodriver analysis, which is only covered by our custom validated NGS assay. CONCLUSIONS: The Plasma-SeqSensei™ SOLID CANCER IVD kit resulted in de novo detection of targetable oncogenic drivers and resistance alterations, with a high sensitivity and accuracy for low and high cfDNA inputs. Thus, this assay is a sensitive, robust, and accurate test.