Amplicon-based Deep Sequencing Research Articles

Description of a novel splice site variant in UBA1 gene causing VEXAS syndrome.

The vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic (VEXAS) syndrome is a complex immune disorder consequence of somatic UBA1 variants. Most reported pathogenic UBA1 variants are missense or splice site mutations directly impairing the translational start site at p. Met41, with recent studies showing that these variants are frequent causes of recurrent inflammation in older individuals. Here we aimed to characterize a novel UBA1 variant found in two patients clinically presenting with VEXAS syndrome. Patients' data were collected from direct assessments and from their medical charts. Genomics analyses were performed by both Sanger and amplicon-based deep sequencing, mRNA studies were performed by both cDNA subcloning and mRNA sequencing. We report a novel, somatic variant in a canonical splice site of the UBA1 gene (c.346-2A>G), which was identified in two unrelated adult male patients with late-onset, unexplained inflammatory manifestations including recurrent fever, Sweet syndrome-like neutrophilic dermatosis, and lung inflammation responsive only to glucocorticoids. RNA analysis from patients' samples demonstrated aberrant mRNA splicing leading to multiple in-frame transcripts, including a transcript retaining the full sequence of intron 4 and a different transcript with the deletion of the first 15 nucleotides of exon 5. Here we describe the abnormal UBA1 transcription as a consequence of the novel c.346-2A>G variant identified in two patients with clinical features compatible with VEXAS syndrome. Overall, these results further demonstrate the expanding spectrum of variants in UBA1 leading to pathology and support for a complete gene evaluation in those candidate patients for VEXAS syndrome.

Spanish cohort of VEXAS syndrome: clinical manifestations, outcome of treatments and novel evidences about UBA1 mosaicism

BackgroundThe vacuoles, E1-enzyme, X linked, autoinflammatory and somatic (VEXAS) syndrome is an adult-onset autoinflammatory disease (AID) due to postzygotic UBA1 variants.ObjectivesTo investigate the presence of VEXAS syndrome among patients with...

Distal arthrogryposis in a girl arising from a novel TNNI2 variant inherited from paternal somatic mosaicism.

TNNI2 at 11p15.5 encodes troponin I2, fast skeletal type, which is a member of the troponin I gene family and a component of the troponin complex. Distal arthrogryposis (DA) is characterized by congenital limb contractures without primary neurological or muscular effects. DA is inherited in an autosomal dominant fashion and is clinically and genetically heterogeneous. Exome sequencing identified a causative variant in TNNI2 [NM_003282.4:c.532T>C p.(Phe178Leu)] in a Japanese girl with typical DA2b. Interestingly, the familial study using Sanger sequencing suggested a mosaic variant in her healthy father. Subsequent targeted amplicon-based deep sequencing detected the TNNI2 variant with variant allele frequencies of 9.4-17.7% in genomic DNA derived from peripheral blood leukocytes, saliva, hair, and nails in the father. We confirmed a disease-causing variant in TNNI2 in the proband inherited from her asymptomatic father with its somatic variant. Our case demonstrates that careful clinical and genetic evaluation is required in DA.

Detecting genomic mosaicism in “de novo” genetic epilepsy by amplicon-based deep sequencing

To investigate the occurrence of mosaicism in epilepsy probands and their parents using amplicon-based deep sequencing (ADS). Patients were recruited from the outpatient of Peking University First Hospital. Two hundred and sixty-four probands with pathogenic variants tested by next-generation sequencing (NGS) were enrolled. Mosaic variants were detected in seventeen disease-associated genes from 20 probands, 5 paternal, and 6 maternal parents. The frequency of mosaicism was 11.74% (31/264). Mosaicism in 11 genes was identified from 20 probands with the mutant allelic fractions (MAFs) of 12.95-38.00% in autosomal dominant genes. Five paternal mosaicisms were identified in genes with a MAF of 6.30-20.99%, and six maternal mosaic individuals with a MAF of 2.07-21.90%. Only four mosaic parents had milder seizure history. The affected sibling had the same phenotype consistent with that of the proband, who inherited the variant of SLC1A2 or STXBP1 from their unaffected mosaic mothers, respectively. Mosaic phenomenon is not rare in families with epilepsy. Phenotypes of mosaic parents were milder or normal. Mosaicism detection is helpful to identify the mutation origin and it provides a theoretical basis for prenatal diagnosis of family reproduction. ADS is a reliable way of mosaicism detection for clinical application.

PCDH19-related epilepsy in mosaic males: The phenotypic implication of genotype and variant allele frequency.

To analyze the genotypes and phenotypes of mosaic male patients with PCDH19-related epilepsy (PCDH19-RE) and explore the correlation between genotype, variant allele frequency (VAF), and phenotypic severity. Clinical data and peripheral blood samples of 11 male mosaic patients were collected and analyzed in our study. The VAF of the PCDH19 gene from peripheral blood was quantified using amplicon-based deep sequencing. Additional 20 mosaic male patients with PCDH19-RE were collected from the published literature, with 10 patients whose VAFs of the PCDH19 gene were available for analytic purposes. In our cohort of 11 patients, 10 variants were identified, and four were novel. The VAF of the PCDH19 gene from peripheral blood ranged from 27 to 90%. The median seizure onset age was 6 months (range: 4-9 months). Clinical manifestations included cluster seizures (100%), fever sensitivity (73%), focal seizures (91%), developmental delay/intellectual disability (DD/ID, 82%), and autistic features (45%). Thirty-one mosaic male patients collected from our cohort and the literature developed seizures mostly (87%) within one year of age. Variant types included missense variants (42%), truncating variants (52%), splice variants (3%), and whole PCDH19 deletion (3%). Among 21 patients with a definite VAF from our cohort and the literature, nine had a low VAF ( ≤ 50%) and 12 had a high VAF (> 50%). Seventy-five percent of variants from the high VAF group were missense, whereas 89% of those from the low VAF group were truncations. The median seizure onset age was 6 months in the low VAF group and 9 months in the high VAF group (p = 0.018). Forty-four percent (4/9) of patients from the low VAF group achieved seizure-free for ≥1 year, whereas none of the 12 patients from the high VAF group did (p = 0.021). DD/ID was present in 83% (10/12) of the high VAF group and 56% (5/9) of the low VAF group (p = 0.331). The predominant variant types were truncating and missense variants. Missense variants tended to have higher VAFs. Patients with a high VAF were more likely to have a more severe epileptic phenotype. Our findings shed light on the phenotypic implications of VAF in mosaic males with PCDH19-RE.

Elucidating Hexanucleotide Repeat Number and Methylation within the X-Linked Dystonia-Parkinsonism (XDP)-Related SVA Retrotransposon in TAF1 with Nanopore Sequencing.

Background: X-linked dystonia-parkinsonism (XDP) is an adult-onset neurodegenerative disorder characterized by progressive dystonia and parkinsonism. It is caused by a SINE-VNTR-Alu (SVA) retrotransposon insertion in the TAF1 gene with a polymorphic (CCCTCT)n domain that acts as a genetic modifier of disease onset and expressivity. Methods: Herein, we used Nanopore sequencing to investigate SVA genetic variability and methylation. We used blood-derived DNA from 96 XDP patients for amplicon-based deep Nanopore sequencing and validated it with fragment analysis which was performed using fluorescence-based PCR. To detect methylation from blood- and brain-derived DNA, we used a Cas9-targeted approach. Results: High concordance was observed for hexanucleotide repeat numbers detected with Nanopore sequencing and fragment analysis. Within the SVA locus, there was no difference in genetic variability other than variations of the repeat motif between patients. We detected high CpG methylation frequency (MF) of the SVA and flanking regions (mean MF = 0.94, SD = ±0.12). Our preliminary results suggest only subtle differences between the XDP patient and the control in predicted enhancer sites directly flanking the SVA locus. Conclusions: Nanopore sequencing can reliably detect SVA hexanucleotide repeat numbers, methylation and, lastly, variation in the repeat motif.

MultiEditR: The first tool for the detection and quantification of RNA editing from Sanger sequencing demonstrates comparable fidelity to RNA-seq

We present MultiEditR (Multiple Edit Deconvolution by Inference of Traces in R), the first algorithm specifically designed to detect and quantify RNA editing from Sanger sequencing (z.umn.edu/multieditr). Although RNA editing is routinely evaluated by measuring the heights of peaks from Sanger sequencing traces, the accuracy and precision of this approach has yet to be evaluated against gold standard next-generation sequencing methods. Through a comprehensive comparison to RNA sequencing (RNA-seq) and amplicon-based deep sequencing, we show that MultiEditR is accurate, precise, and reliable for detecting endogenous and programmable RNA editing.

Assessment of the gene mosaicism burden in blood and its implications for immune disorders

There are increasing evidences showing the contribution of somatic genetic variants to non-cancer diseases. However, their detection using massive parallel sequencing methods still has important limitations. In addition, the relative importance and dynamics of somatic variation in healthy tissues are not fully understood. We performed high-depth whole-exome sequencing in 16 samples from patients with a previously determined pathogenic somatic variant for a primary immunodeficiency and tested different variant callers detection ability. Subsequently, we explored the load of somatic variants in the whole blood of these individuals and validated it by amplicon-based deep sequencing. Variant callers allowing low frequency read thresholds were able to detect most of the variants, even at very low frequencies in the tissue. The genetic load of somatic coding variants detectable in whole blood is low, ranging from 1 to 2 variants in our dataset, except for one case with 17 variants compatible with clonal haematopoiesis under genetic drift. Because of the ability we demonstrated to detect this type of genetic variation, and its relevant role in disorders such as primary immunodeficiencies, we suggest considering this model of gene mosaicism in future genetic studies and considering revisiting previous massive parallel sequencing data in patients with negative results.

A novel CRISPR-based malaria diagnostic capable of Plasmodium detection, species differentiation, and drug-resistance genotyping

BackgroundCRISPR-based diagnostics are a new class of highly sensitive and specific assays with multiple applications in infectious disease diagnosis. SHERLOCK, or Specific High-Sensitivity Enzymatic Reporter UnLOCKing, is one such CRISPR-based diagnostic that combines recombinase polymerase pre-amplification, CRISPR-RNA base-pairing, and LwCas13a activity for nucleic acid detection.MethodsWe developed SHERLOCK assays capable of detecting all Plasmodium species known to cause human malaria and species-specific detection of P. vivax and P. falciparum, the species responsible for the majority of malaria cases worldwide. We further tested these assays using a diverse panel of clinical samples from the Democratic Republic of the Congo, Uganda, and Thailand and pools of Anopheles mosquitoes from Thailand. In addition, we developed a prototype SHERLOCK assay capable of detecting the dihydropteroate synthetase (dhps) single nucleotide variant A581G associated with P. falciparum sulfadoxine resistance.FindingsThe suite of Plasmodium assays achieved analytical sensitivities ranging from 2•5-18•8 parasites per reaction when tested against laboratory strain genomic DNA. When compared to real-time PCR, the P. falciparum assay achieved 94% sensitivity and 94% specificity during testing of 123 clinical samples. Compared to amplicon-based deep sequencing, the dhps SHERLOCK assay achieved 73% sensitivity and 100% specificity when applied to a panel of 43 clinical samples, with false-negative calls only at lower parasite densities.InterpretationThese novel SHERLOCK assays demonstrate the versatility of CRISPR-based diagnostics and their potential as a new generation of molecular tools for malaria diagnosis and surveillance.FundingNational Institutes of Health (T32GM007092, R21AI148579, K24AI134990, R01AI121558, UL1TR002489, P30CA016086)

A systematic study and literature review of parental somatic mosaicism of FBN1 pathogenic variants in Marfan syndrome

BackgroundA proportion of de novo variants in patients affected by genetic disorders, particularly those with autosomal dominant (AD) inheritance, could be the consequence of somatic mosaicism in one of the...

Abstract PD10-02: The influence of pathogenic variants in breast cancer predisposition genes on secondary breast cancer events in a prospectively collected cohort of breast cancer patients

Abstract Background: About 10-20% of breast cancers (BC) are hereditary. Patients (pts) with BC and pathogenic variants (PVs) in the BRCA1 or BRCA2 gene are at increased risk of contralateral breast cancer (CBC). However, risks of CBC and recurrent BC associated with PVs in other established BC predisposition genes (e.g. PALB2, CHEK2, ATM) are less well established. We performed research genetic testing of over 4000 pts with BC treated at Mayo Clinic. Because most pts were unaware of their mutation status prior to selection of therapy, this cohort provides a unique opportunity to evaluate the natural history of BC recurrence after treatment of the index BC in pts with a deleterious mutation in predisposition genes. Methods: Women diagnosed and treated for BC at Mayo Clinic between 2001-2016 were enrolled in an IRB-approved research cohort. Pts were not selected for age of diagnosis or family history. Coding regions and consensus splice sites of established BC predisposition genes were subjected to amplicon-based deep sequencing using the QIAseq method in germline DNA extracted from peripheral blood or buccal samples. PVs in the genes were identified using GATK haplotype caller and Vardict. Patient, tumor, oncologic management and follow up data were reviewed. Herein we report on the differences in the incidence of local recurrence and CBC, among carriers and non- carriers and by predisposition gene using the cumulative incidence estimator and Fine and Gray regression to account for competing risks. Results: Of 4272 women with BC (median age 57, range 21-94), a PV in a BC predisposition gene was identified in 287 (6.7%). The most common PVs were seen in CHEK2 (67, 1.6%), BRCA1 (62, 1.5%), BRCA2 (60, 1.4%), ATM (47, 1.1%); 8 (0.2%) patients had more than one PV. Pts with any PVs presented at younger age (median 51 vs 58 years, p<0.001), were more likely to have stage IV disease (4.5% vs 2.3%, p=0.02), estrogen receptor (ER) negative BC (19.5% vs 13.5%, p=0.005), and bilateral cancer (8.8% vs 4.3%, p<0.001). ER- BC proportion varied by gene - 38.9% in pts with PVs associated with ER- BC (BRCA1, BARD1, BRIP1, RAD51C, RAD51D) and 11.9% in those with PVs associated with ER+ BC (BRCA2, ATM, CDH1, CHEK2), p<0.001. NCCN criteria for genetic testing were assessed for 3763, with 69.9% of those with PVs qualifying for testing compared to 46.3% of those without PVs (p<0.001). Among 3973 pts (249 with PVs, 3724 without) with unilateral stage 0-III cancer, 168 pts developed local recurrence, 131 pts developed CBC and 248 pts developed distant disease with median 5 years of follow-up. Patients with any PV had higher rates of local recurrence after breast conserving surgery (HR 2.3, 95% CI: 1.2-4.2, p=0.01) but overall did not have significantly higher risk of CBC if contralateral prophylactic mastectomy was not performed (HR 1.7, 95% CI: 0.9-3.2, p=0.11). Specifically, pts with PVs in genes other than BRCA1 or BRCA2 did not have a higher risk of CBC compared to pts with no PV (HR 1.0, 95% CI: 0.4-2.7, p=0.96); the comparison of CBC risk in BRCA1 or BRCA2 versus those with no PV was HR 2.9 (95% CI: 1.3-6.5, p=0.01), while BRCA1 or BRCA2 compared to those with other PVs was HR 2.7 (95% CI: 0.8-9.3, p=0.13). In pts with CHEK2 PVs, the probability of CBC at 10-years follow-up was 5.4%, which was similar to pts without a mutation (5.3%), p=0.47. Conclusion: While PVs in known BC predisposition genes increase risk of BC development, risk of CBC in carriers of PVs in non-BRCA predisposition genes with BC are not as high as seen with BRCA1 or 2 PV carriers. This information regarding future BC will be helpful in personalizing decisions on the extent of surgery in carriers of PVs in BC predisposition genes. Citation Format: Judy C. Boughey, Tanya L. Hoskin, Siddhartha Yadav, Fergus J. Couch. The influence of pathogenic variants in breast cancer predisposition genes on secondary breast cancer events in a prospectively collected cohort of breast cancer patients [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PD10-02.

Evaluation of circulating tumor DNA as a biomarker in pancreatic cancer with liver metastasis.

Pancreatic cancer is an aggressive, solid tumor, with a grave prognosis. Despite surgical treatment in patients with pancreatic cancer, the rate of recurrence is high. In addition, although tumor biomarkers are frequently used to confirm advanced pancreatic cancer, this is not accurate and the biomarkers currently used cannot indicate prognosis. This study sought to evaluate circulating tumor DNA as a tumor biomarker to prognosticate pancreatic cancer. Patients with advanced pancreatic cancer and liver metastasis (N = 104) were included, and blood samples were collected from all patients. The mutant allele frequency was measured using amplicon-based deep sequencing on a cell-free DNA panel covering 14 genes with > 240 hot spots. In patients with advanced pancreatic cancer, 50% (N = 52) had detectable ctDNA levels, with TP53 (45%, N = 47) and KRAS (42.3%, N = 44) mutations the most common. Patients with detectable circulating tumor DNA levels also had significantly worse overall survival and progression free survival than ctDNA negative patients (8.4 vs 16 months, P<0.0001 for overall survival; 3.2 vs 7.9 months, P<0.0001 for progression-free survival). In a multivariate analysis, ctDNA status was independently associated with overall survival and progression-free survival (HR = 3.1, 95%CI = 1.9–5.0, P<0.0001; HR 2.6, 95%CI = 1.7–4.0, P<0.0001, respectively). Moreover, circulating tumor DNA significantly correlated with a higher number of liver metastases, the presence of lung and/or peritoneal metastases, tumor burden, and higher carbohydrate antigen 19–9 levels. This study supports the use of circulating tumor DNA as an independent prognostic marker for advanced pancreatic cancer.

A somatic activating KRAS variant identified in an affected lesion of a patient with Gorham-Stout disease.

Gorham-Stout disease (GSD), a rare disorder of unknown etiology, is characterized by massive osteolysis that is associated with proliferation and dilation of lymphatic vessels. Variants in cancer-associated genes have been described in complex lymphatic anomalies. To explore the pathogenesis of GSD, we performed the amplicon-based deep sequencing on 50 cancer-related genes to assay affected tissues from the six patients with GSD. In one patient, a somatic activating KRAS c.182A > G variant (p.Q61R) was detected in 1% of the tissue sample. Conversely, the mutant allele was not detected in uninvolved normal skin and blood samples. Histopathology of the patient's tissue sample showed proliferation of abnormal lymphatic and blood vascular endothelial cells, osteoclasts, and activated macrophages. The activating KRAS variant is a known 'hotspot' variant, frequently identified in several types of human cancer. This is the first report of identifying a pathogenic variant in a patient with GSD. This finding may set the stage for elucidation of pathophysiology and the development of novel therapies for GSD.

Serial circulating tumour DNA analysis for locally advanced rectal cancer treated with preoperative therapy: prediction of pathological response and postoperative recurrence

BackgroundThe “watch-and-wait” approach is a common treatment option amongst patients with locally advanced rectal cancer (LARC). However, the diagnostic sensitivity of clinical modalities, such as colonoscopy and magnetic resonance imaging to determine pathological response, is not high. We analysed the clinical utility of circulating tumour DNA (ctDNA) of patients with LARC to predict response to preoperative therapy and postoperative recurrence.MethodsA serial ctDNA analysis of 222 plasma samples from 85 patients with LARC was performed using amplicon-based deep sequencing on a cell-free DNA panel covering 14 genes with over 240 hotspots.ResultsctDNA was detected in 57.6% and 22.3% of samples at baseline and after preoperative treatment, respectively, which was significantly different (P = 0.0003). Change in ctDNA was an independent predictor of complete response to preoperative therapy (P = 0.0276). In addition, postoperative ctDNA and carcinoembryonic antigen (CEA) were independent prognostic markers for risk of recurrence after surgery (ctDNA, P = 0.0127 and CEA, P = 0.0105), with a combined analysis having cumulative effects on recurrence-free survival (P = 1.0 × 10–16).ConclusionsSerial ctDNA analysis may offer clinically useful predictive and prognostic markers for response to preoperative therapy and postoperative recurrence in patients with LARC.

Abstract A37: Early change in circulating tumor DNA as a potential predictor of response to chemotherapy in patients with metastatic colorectal cancer

Abstract The impact of the change in ctDNA levels after chemotherapy on the clinical outcomes of patients with metastatic colorectal cancer (mCRC) remains unclear. The present study evaluated the clinical implications of the early change in ctDNA levels as a predictor of objective response and the clinical outcome in mCRC patients who received chemotherapy. We investigated the effects of after/before ratio of ctDNA levels 2 and 8 weeks after initiation of second-line chemotherapy on objective response rate (ORR), progression-free survival (PFS), and overall survival (OS). ctDNA was detected using amplicon-based deep sequencing with a molecular barcode that includes &amp;gt;240 hotspot mutations in 14 colon cancer-related genes. Patients with lower ctDNA level (≤ 50%) 8 weeks after initiation of chemotherapy compared to baseline showed significantly longer PFS and OS than the patients with higher (&amp;gt; 50%) ctDNA level in multivariate analysis. In patients achieving partial response or stable disease, the after/before ratio of ctDNA level 8 weeks after initiation of chemotherapy was significantly lower than those in patients with progressive disease. The present study suggests that an early change in the ctDNA level might serve as a biomarker to predict the chemotherapeutic efficacy and clinical outcomes in patients with mCRC. Citation Format: Hitoshi Zembutsu, Hiroki Osumi, Eiji Shinozaki, Kensei Yamaguchi. Early change in circulating tumor DNA as a potential predictor of response to chemotherapy in patients with metastatic colorectal cancer [abstract]. In: Proceedings of the AACR Special Conference on Advances in Liquid Biopsies; Jan 13-16, 2020; Miami, FL. Philadelphia (PA): AACR; Clin Cancer Res 2020;26(11_Suppl):Abstract nr A37.

I-Silence, Please! An Alternative for Gene Disruption via Adenine Base Editors.

i-Silence, Please! An Alternative for Gene Disruption via Adenine Base Editors.

Early change in circulating tumor DNA as a potential predictor of response to chemotherapy in patients with metastatic colorectal cancer

The impact of ctDNA changes after chemotherapy on the clinical outcomes of patients with metastatic colorectal cancer (mCRC) remains unclear. The present study evaluated the clinical implications of the early change in ctDNA levels as a predictor of objective response and clinical outcome in mCRC patients who received chemotherapy. We investigated the effects of after/before ratio of ctDNA levels 2 and 8 weeks after initiation of second-line chemotherapy, on objective response rate (ORR), progression-free survival (PFS), and overall survival (OS). ctDNA was detected using amplicon-based deep sequencing with a molecular barcode encompassing >240 hotspot mutations in 14 colon cancer-related genes. In multivariate analysis, as compared to baseline, patients with lower ctDNA level (≤50%) 8 weeks after initiation of chemotherapy showed significantly longer PFS and OS than the patients with higher (>50%) ctDNA level. In patients achieving a partial response or stable disease, the after/before ratio of ctDNA level 8 weeks after initiation of chemotherapy was significantly lower than those in patients with progressive disease. The present study suggests that an early change in the ctDNA level might serve as a biomarker to predict the chemotherapeutic efficacy and clinical outcomes in patients with mCRC.

Sequence characterization of RET in 117 Chinese Hirschsprung disease families identifies a large burden of de novo and parental mosaic mutations

BackgroundHirschsprung disease (HSCR) is an inherited congenital disorder characterized by the absence of enteric ganglia in the distal part of the gut. RET is the major causative gene and contains > 80% of all known disease-causing mutations.ResultsTo determine the incidence of RET pathogenic variants, be they Mendelian inherited, mosaic in parents or true de novo variants (DNVs) in 117 Chinese families, we used high-coverage NGS and droplet digital polymerase chain reaction (ddPCR) to identify 15 (12.8%) unique RET coding variants (7 are novel); one was inherited from a heterozygous unaffected mother, 11 were DNVs (73.3%), and 3 full heterozygotes were inherited from parental mosaicism (2 paternal, 1 maternal): two clinically unaffected parents were identified by NGS and confirmed by ddPCR, with mutant allele frequency (13–27%) that was the highest in hair, lowest in urine and similar in blood and saliva. An extremely low-level paternal mosaicism (0.03%) was detected by ddPCR in blood. Six positive-controls were examined to compare the mosaicism detection limit and sensitivity of NGS, amplicon-based deep sequencing and ddPCR.ConclusionOur findings expand the clinical and molecular spectrum of RET variants in HSCR and reveal a high frequency of RET DNVs in the Chinese population.

Utility of Targeted, Amplicon-Based Deep Sequencing To Detect Resistance to First-Line Tuberculosis Drugs in Botswana.

Multidrug-resistant tuberculosis (TB) is an alarming threat, and targeted deep sequencing (DS) may be an effective method for rapid identification of drug-resistant profiles, including detection of heteroresistance. We evaluated the sensitivity and specificity of targeted DS versus phenotypic drug susceptibility testing (pDST) among patients starting first-line anti-TB therapy in Botswana. Overall, we found high concordance between DS and pDST. Lower sensitivity of DS, which targets established high-confidence resistance variants, was observed for detecting isoniazid resistance among HIV-infected patients.

Identification of single nucleotide variants using position-specific error estimation in deep sequencing data

BackgroundTargeted deep sequencing is a highly effective technology to identify known and novel single nucleotide variants (SNVs) with many applications in translational medicine, disease monitoring and cancer profiling. However, identification of SNVs using deep sequencing data is a challenging computational problem as different sequencing artifacts limit the analytical sensitivity of SNV detection, especially at low variant allele frequencies (VAFs).MethodsTo address the problem of relatively high noise levels in amplicon-based deep sequencing data (e.g. with the Ion AmpliSeq technology) in the context of SNV calling, we have developed a new bioinformatics tool called AmpliSolve. AmpliSolve uses a set of normal samples to model position-specific, strand-specific and nucleotide-specific background artifacts (noise), and deploys a Poisson model-based statistical framework for SNV detection.ResultsOur tests on both synthetic and real data indicate that AmpliSolve achieves a good trade-off between precision and sensitivity, even at VAF below 5% and as low as 1%. We further validate AmpliSolve by applying it to the detection of SNVs in 96 circulating tumor DNA samples at three clinically relevant genomic positions and compare the results to digital droplet PCR experiments.ConclusionsAmpliSolve is a new tool for in-silico estimation of background noise and for detection of low frequency SNVs in targeted deep sequencing data. Although AmpliSolve has been specifically designed for and tested on amplicon-based libraries sequenced with the Ion Torrent platform it can, in principle, be applied to other sequencing platforms as well. AmpliSolve is freely available at https://github.com/dkleftogi/AmpliSolve.