Characterization and mitigation of artifacts derived from NGS library preparation due to structure-specific sequences in the human genome.

BACKGROUND: Hybridization capture-based targeted next generation sequencing (NGS) is gaining importance in routine cancer clinical practice. DNA library preparation is a fundamental step to produce high-quality sequencing data. Numerous unexpected, low variant allele frequency calls were observed in libraries using sonication fragmentation and enzymatic fragmentation. In this study, we investigated the characteristics of the artifact reads induced by sonication and enzymatic fragmentation. We also developed a bioinformatic algorithm to filter these sequencing errors. RESULTS: We used pairwise comparisons of somatic single nucleotide variants (SNVs) and insertions and deletions (indels) of the same tumor DNA samples prepared using both ultrasonic and enzymatic fragmentation protocols. Our analysis revealed that the number of artifact variants was significantly greater in the samples generated using enzymatic fragmentation than using sonication. Most of the artifacts derived from the sonication-treated libraries were chimeric artifact reads containing both cis- and trans-inverted repeat sequences of the genomic DNA. In contrast, chimeric artifact reads of endonuclease-treated libraries contained palindromic sequences with mismatched bases. Based on these distinctive features, we proposed a mechanistic hypothesis model, PDSM (pairing of partial single strands derived from a similar molecule), by which these sequencing errors derive from ultrasonication and enzymatic fragmentation library preparation. We developed a bioinformatic algorithm to generate a custom mutation "blacklist" in the BED region to reduce errors in downstream analyses. CONCLUSIONS: We first proposed a mechanistic hypothesis model (PDSM) of sequencing errors caused by specific structures of inverted repeat sequences and palindromic sequences in the natural genome. This new hypothesis predicts the existence of chimeric reads that could not be explained by previous models, and provides a new direction for further improving NGS analysis accuracy. A bioinformatic algorithm, ArtifactsFinder, was developed and used to reduce the sequencing errors in libraries produced using sonication and enzymatic fragmentation.

Genotype-phenotype of autosomal dominant polycystic kidney disease in Malta.

BACKGROUND: Autosomal dominant polycystic kidney disease (ADPKD) is characterized by the development of multiple renal cysts causing kidney enlargement and end-stage renal disease (ESRD) in half the patients by 60 years of age. The aim of the study was to determine the genetic aetiology in Maltese patients clinically diagnosed with ADPKD and correlate the clinical features. METHODS: A total of 60 patients over 18 years of age clinically diagnosed with ADPKD were studied using a customized panel of genes that had sufficient evidence of disease diagnosis using next generation sequencing (NGS). The genes studied were PKD1, PKD2, GANAB, DNAJB11, PKHD1 and DZIP1L. Selected variants were confirmed by bidirectional Sanger sequencing with specifically designed primers. Cases where no clinically significant variant was identified by the customized gene panel were then studied by Whole Exome Sequencing (WES). Microsatellite analysis was performed to determine the origin of an identified recurrent variant in the PKD2 gene. Clinical features were studied for statistical correlation with genetic results. RESULTS: Genetic diagnosis was reached in 49 (82%) of cases studied. Pathogenic/likely pathogenic variants PKD1 and PKD2 gene were found in 25 and in 23 cases respectively. The relative proportion of genetically diagnosed PKD1:PKD2 cases was 42:38. A pathogenic variant in the GANAB gene was identified in 1 (2%) case. A potentially significant heterozygous likely pathogenic variant was identified in PKHD1 in 1 (2%) case. Potentially significant variants of uncertain significance were seen in 4 (7%) cases of the study cohort. No variants in DNAJB11 and DZIP1L were observed. Whole exome sequencing (WES) added the diagnostic yield by 10% over the gene panel analysis. Overall no clinically significant variant was detected in 6 (10%) cases of the study population by a customized gene panel and WES. One recurrent variant the PKD2 c.709+1G > A was observed in 19 (32%) cases. Microsatellite analysis showed that all variant cases shared the same haplotype indicating that their families may have originated from a common ancestor and confirmed it to be a founder variant in the Maltese population. The rate of decline in eGFR was steeper and progression to ESRD was earlier in cases with PKD1 variants when compared to cases with PKD2 variants. Cases segregating truncating variants in PKD1 showed a significantly earlier onset of ESRD and this was significantly worse in cases with frameshift variants. Overall extrarenal manifestations were commoner in cases segregating truncating variants in PKD1. CONCLUSIONS: This study helps to show that a customized gene panel is the first-line method of choice for studying patients with ADPKD followed by WES which increased the detection of variants present in the PKD1 pseudogene region. A founder variant in the PKD2 gene was identified in our Maltese cohort with ADPKD. Phenotype of patients with ADPKD is significantly related to the genotype confirming the important role of molecular investigations in the diagnosis and prognosis of polycystic kidney disease. Moreover, the findings also highlight the variability in the clinical phenotype and indicate that other factors including epigenetic and environmental maybe be important determinants in Autosomal Dominant Polycystic Kidney Disease.

Tuberculosis-targeted next-generation sequencing and machine learning: An ultrasensitive diagnostic strategy for paucibacillary pulmonary tuberculosis and tuberculous meningitis.

BACKGROUND: Existing diagnostic approaches for paucibacillary tuberculosis (TB) are limited by the low sensitivity of testing methods and difficulty in obtaining suitable samples. METHODS: An ultrasensitive TB diagnostic strategy was established, integrating efficient and specific TB targeted next-generation sequencing and machine learning models, and validated in clinical cohorts to test plasma cfDNA, cerebrospinal fluid (CSF) DNA collected from tuberculous meningitis (TBM) and pediatric pulmonary TB (PPTB) patients. RESULTS: In the detection of 227 samples, application of the specific thresholds of CSF DNA (AUC = 0.974) and plasma cfDNA (AUC = 0.908) yielded sensitivity of 97.01 % and the specificity of 95.65 % in CSF samples and sensitivity of 82.61 % and specificity of 86.36 % in plasma samples, respectively. In the analysis of 44 paired samples from TBM patients, our strategy had a high concordance of 90.91 % (40/44) in plasma cfDNA and CSF DNA with both sensitivity of 95.45 % (42/44). In the PPTB patient, the sensitivity of the TB diagnostic strategy yielded higher sensitivity on plasma specimen than Xpert assay on gastric lavage (28.57 % VS. 15.38 %). CONCLUSIONS: Our TB diagnostic strategy provides greater detection sensitivity for paucibacillary TB, while plasma cfDNA as an easily collected specimen, could be an appropriate sample type for PTB and TBM diagnosis.

Novel mutations of the ABCA12, KRT1 and ST14 genes in three unrelated newborns showing congenital ichthyosis.

BACKGROUND: Congenital ichthyosis (CI) is a heterogeneous group of genetic disorders characterized by generalized dry skin, scaling and hyperkeratosis, often associated to erythroderma. They are rare diseases, with overall incidence of 6.7 in 100,000. Clinical manifestations are due to mutations in genes mostly involved in skin barrier formation. Based on clinical presentation, CI is distinguished in non-syndromic and syndromic forms. To date, mutations of more than 50 genes have been associated to different types of CI. CASES PRESENTATION: We report on three Italian unrelated newborns showing clinical signs compatible with different forms of CI of variable severity, namely Harlequin ichtyosis (HI), epidermolytic ichtyosis (EI) and autosomal recessive ichtyosis with hypotrichosis (ARIH). Target next generation sequencing (NGS) analysis identified three novel mutations of the ABCA12, KRT1 and ST14 genes, respectively associated to such congenital ichtyoses, not reported in literature. Genomic investigation allowed to provide the more appropriate management to each patient, based on an individualized approach. CONCLUSIONS: Our report highlights the wide genetic heterogeneity and phenotypic variability of CI. It expands the current knowledge on such diseases, widening their genomic database, and providing a better clinical characterization. Furthermore, it underlines the clinical relevance of NGS, which is essential to address the management of patients. Indeed, it may guide towards the most adequate approach, preventing clinical obstinacy for subjects with more severe forms and unfavorable outcomes (together with the support, in such situations, of bioethicists included within the multidisciplinary care team), as well as reassuring families in those with milder course and favorable evolution.