Unravelling undiagnosed rare disease cases by HiFi long-read genome sequencing

Solve-RD is a pan-European rare disease (RD) research program that aims to identify disease-causing genetic variants in previously undiagnosed RD families. We utilised 10-fold coverage HiFi long-read sequencing (LRS) for detecting causative structural variants (SVs), single nucleotide variants (SNVs), insertion-deletions (InDels), and short tandem repeat (STR) expansions in extensively studied RD families without clear molecular diagnoses. Our cohort includes 293 individuals from 114 genetically undiagnosed RD families selected by European Rare Disease Network (ERN) experts. Of these, 21 families were affected by so-called 'unsolvable' syndromes for which genetic causes remain unknown, and 93 families with at least one individual affected by a rare neurological, neuromuscular, or epilepsy disorder without genetic diagnosis despite extensive prior testing. Clinical interpretation and orthogonal validation of variants in known disease genes yielded thirteen novel genetic diagnoses due to de novo and rare inherited SNVs, InDels, SVs, and STR expansions. In an additional four families, we identified a candidate disease-causing SV affecting several genes including an MCF2 / FGF13 fusion and PSMA3 deletion. However, no common genetic cause was identified in any of the 'unsolvable' syndromes. Taken together, we found (likely) disease-causing genetic variants in 13.0% of previously unsolved families and additional candidate disease-causing SVs in another 4.3% of these families. In conclusion, our results demonstrate the added value of HiFi long-read genome sequencing in undiagnosed rare diseases.

Titin copy number variations associated with dominant inherited phenotypes

BACKGROUND

Titinopathies are caused by mutations in the titin gene (TTN). Titin is the largest known human protein; its gene has the longest coding phase with 364 exons. Titinopathies are very complex neuromuscular pathologies due to the variable age of onset of symptoms, the great diversity of pathological and muscular impairment patterns (cardiac, skeletal muscle or mixed) and both autosomal dominant and recessive modes of transmission. Until now, only few CNVs in TTN have been reported without clear genotype-phenotype associations.

METHODS

Our study includes eight families with dominant titinopathies. We performed next-generation sequencing or comparative genomic hybridisation array analyses and found CNVs in the TTN gene. We characterised these CNVs by RNA sequencing (RNAseq) analyses in six patients' muscles and performed genotype-phenotype inheritance association study by combining the clinical and biological data of these eight families.

RESULTS

Seven deletion-type CNVs in the TTN gene were identified among these families. Genotype and RNAseq results showed that five deletions do not alter the reading frame and one is out-of-reading frame. The main phenotype identified was distal myopathy associated with contractures. The analysis of morphological, clinical and genetic data and imaging let us draw new genotype-phenotype associations of titinopathies.

CONCLUSION

Identifying TTN CNVs will further increase diagnostic sensitivity in these complex neuromuscular pathologies. Our cohort of patients enabled us to identify new deletion-type CNVs in the TTN gene, with unexpected autosomal dominant transmission. This is valuable in establishing new genotype-phenotype associations of titinopathies, mainly distal myopathy in most of the patients.

Genomes in clinical care

In the era of precision medicine, genome sequencing (GS) has become more affordable and the importance of genomics and multi-omics in clinical care is increasingly being recognized. However, how to scale and effectively implement GS on an institutional level remains a challenge for many. Here, we present Genome First and Ge-Med, two clinical implementation studies focused on identifying the key pillars and processes that are required to make routine GS and predictive genomics a reality in the clinical setting. We describe our experience and lessons learned for a variety of topics including test logistics, patient care processes, data reporting, and infrastructure. Our model of providing clinical care and comprehensive genomic analysis from a single source may be used by other centers with a similar structure to facilitate the implementation of omics-based personalized health concepts in medicine.

Mobile element insertions in rare diseases: a comparative benchmark and reanalysis of 60,000 exome samples

Mobile element insertions (MEIs) are a known cause of genetic disease but have been underexplored due to technical limitations of genetic testing methods. Various bioinformatic tools have been developed to identify MEIs in Next Generation Sequencing data. However, most tools have been developed specifically for genome sequencing (GS) data rather than exome sequencing (ES) data, which remains more widely used for routine diagnostic testing. In this study, we benchmarked six MEI detection tools (ERVcaller, MELT, Mobster, SCRAMble, TEMP2 and xTea) on ES data and on GS data from publicly available genomic samples (HG002, NA12878). For all the tools we evaluated sensitivity and precision of different filtering strategies. Results show that there were substantial differences in tool performance between ES and GS data. MELT performed best with ES data and its combination with SCRAMble increased substantially the detection rate of MEIs. By applying both tools to 10,890 ES samples from Solve-RD and 52,624 samples from Radboudumc we were able to diagnose 10 patients who had remained undiagnosed by conventional ES analysis until now. Our study shows that MELT and SCRAMble can be used reliably to identify clinically relevant MEIs in ES data. This may lead to an additional diagnosis for 1 in 3000 to 4000 patients in routine clinical ES.

Diagnostic genome sequencing improves diagnostic yield: a prospective single-centre study in 1000 patients with inherited eye diseases

PURPOSE

Genome sequencing (GS) is expected to reduce the diagnostic gap in rare disease genetics. We aimed to evaluate a scalable framework for genome-based analyses 'beyond the exome' in regular care of patients with inherited retinal degeneration (IRD) or inherited optic neuropathy (ION).

METHODS

PCR-free short-read GS was performed on 1000 consecutive probands with IRD/ION in routine diagnostics. Complementary whole-blood RNA-sequencing (RNA-seq) was done in a subset of 74 patients. An open-source bioinformatics analysis pipeline was optimised for structural variant (SV) calling and combined RNA/DNA variation interpretation.

RESULTS

A definite genetic diagnosis was established in 57.4% of cases. For another 16.7%, variants of uncertain significance were identified in known IRD/ION genes, while the underlying genetic cause remained unresolved in 25.9%. SVs or alterations in non-coding genomic regions made up for 12.7% of the observed variants. The RNA-seq studies supported the classification of two unclear variants.

CONCLUSION

GS is feasible in clinical practice and reliably identifies causal variants in a substantial proportion of individuals. GS extends the diagnostic yield to rare non-coding variants and enables precise determination of SVs. The added diagnostic value of RNA-seq is limited by low expression levels of the major IRD disease genes in blood.

Clinical trio genome sequencing facilitates the interpretation of variants in cancer predisposition genes in paediatric tumour patients

The prevalence of pathogenic and likely pathogenic (P/LP) variants in genes associated with cancer predisposition syndromes (CPS) is estimated to be 8-18% for paediatric cancer patients. In more than half of the carriers, the family history is unsuspicious for CPS. Therefore, broad genetic testing could identify germline predisposition in additional children with cancer resulting in important implications for themselves and their families. We thus evaluated clinical trio genome sequencing (TGS) in a cohort of 72 paediatric patients with solid cancers other than retinoblastoma or CNS-tumours. The most prevalent cancer types were sarcoma (n = 26), neuroblastoma (n = 15), and nephroblastoma (n = 10). Overall, P/LP variants in CPS genes were identified in 18.1% of patients (13/72) and P/LP variants in autosomal-dominant CPS genes in 9.7% (7/72). Genetic evaluation would have been recommended for the majority of patients with P/LP variants according to the Jongmans criteria. Four patients (5.6%, 4/72) carried P/LP variants in autosomal-dominant genes known to be associated with their tumour type. With the immediate information on variant inheritance, TGS facilitated the identification of a de novo P/LP in NF1, a gonadosomatic mosaic in WT1 and two pathogenic variants in one patient (DICER1 and PALB2). TGS allows a more detailed characterization of structural variants with base-pair resolution of breakpoints which can be relevant for the interpretation of copy number variants. Altogether, TGS allows comprehensive identification of children with a CPS and supports the individualised clinical management of index patients and high-risk relatives.

UV-radiation and MC1R germline mutations are risk factors for the development of conventional and spitzoid melanomas in children and adolescents

BACKGROUND

Genomic characterisation has led to an improved understanding of adult melanoma. However, the aetiology of melanoma in children is still unclear and identifying the correct diagnosis and therapeutic strategies remains challenging.

METHODS

Exome sequencing of matched tumour-normal pairs from 26 paediatric patients was performed to study the mutational spectrum of melanomas. The cohort was grouped into different categories: spitzoid melanoma (SM), conventional melanoma (CM), and other melanomas (OT).

FINDINGS

In all patients with CM (n = 10) germline variants associated with melanoma were found in low to moderate melanoma risk genes: in 8 patients MC1R variants, in 2 patients variants in MITF, PTEN and BRCA2. Somatic BRAF mutations were detected in 60% of CMs, homozygous deletions of CDKN2A in 20%, TERTp mutations in 30%. In the SM group (n = 12), 5 patients carried at least one MC1R variant; somatic BRAF mutations were detected in 8.3%, fusions in 25% of the cases. No SM showed a homozygous CDKN2A deletion nor a TERTp mutation. In 81.8% of the CM/SM cases the UV damage signatures SBS7 and/or DBS1 were detected. The patient with melanoma arising in giant congenital nevus (CNM) demonstrated the characteristic NRAS Q61K mutation.

INTERPRETATION

UV-radiation and MC1R germline variants are risk factors in the development of conventional and spitzoid paediatric melanomas. Paediatric CMs share genomic similarities with adult CMs while the SMs differ genetically from the CM group. Consistent genetic characterization of all paediatric melanomas will potentially lead to better subtype differentiation, treatment, and prevention in the future.

FUNDING

Found in Acknowledgement.

Combined loss of CDH1 and downstream regulatory sequences drive early-onset diffuse gastric cancer and increase penetrance of hereditary diffuse gastric cancer

BACKGROUND

Germline CDH1 pathogenic or likely pathogenic variants cause hereditary diffuse gastric cancer (HDGC). Once a genetic cause is identified, stomachs' and breasts' surveillance and/or prophylactic surgery is offered to asymptomatic CDH1 carriers, which is life-saving. Herein, we characterized an inherited mechanism responsible for extremely early-onset gastric cancer and atypical HDGC high penetrance.

METHODS

Whole-exome sequencing (WES) re-analysis was performed in an unsolved HDGC family. Accessible chromatin and CDH1 promoter interactors were evaluated in normal stomach by ATAC-seq and 4C-seq, and functional analysis was performed using CRISPR-Cas9, RNA-seq and pathway analysis.

RESULTS

We identified a germline heterozygous 23 Kb CDH1-TANGO6 deletion in a family with eight diffuse gastric cancers, six before age 30. Atypical HDGC high penetrance and young cancer-onset argued towards a role for the deleted region downstream of CDH1, which we proved to present accessible chromatin, and CDH1 promoter interactors in normal stomach. CRISPR-Cas9 edited cells mimicking the CDH1-TANGO6 deletion display the strongest CDH1 mRNA downregulation, more impacted adhesion-associated, type-I interferon immune-associated and oncogenic signalling pathways, compared to wild-type or CDH1-deleted cells. This finding solved an 18-year family odyssey and engaged carrier family members in a cancer prevention pathway of care.

CONCLUSION

In this work, we demonstrated that regulatory elements lying down-stream of CDH1 are part of a chromatin network that control CDH1 expression and influence cell transcriptome and associated signalling pathways, likely explaining high disease penetrance and very young cancer-onset. This study highlights the importance of incorporating scientific-technological updates and clinical guidelines in routine diagnosis, given their impact in timely genetic diagnosis and disease prevention.

Detection and Characterization of a De Novo Alu Retrotransposition Event Causing NKX2-1-Related Disorder

BACKGROUND

Heterozygous NKX2-1 loss-of-function variants cause combinations of hyperkinetic movement disorders (MDs, particularly childhood-onset chorea), pulmonary dysfunction, and hypothyroidism. Mobile element insertions (MEIs) are potential disease-causing structural variants whose detection in routine diagnostics remains challenging.

OBJECTIVE

To establish the molecular diagnosis of two first-degree relatives with clinically suspected NKX2-1-related disorder who had negative NKX2-1 Sanger (SS), whole-exome (WES), and whole-genome (WGS) sequencing.

METHODS

The proband's WES was analyzed for MEIs. A candidate MEI in NKX2-1 underwent optimized SS after plasmid cloning. Functional studies exploring NKX2-1 haploinsufficiency at RNA and protein levels were performed.

RESULTS

A 347-bp AluYa5 insertion with a 65-bp poly-A tail followed by a 16-bp duplication of the pre-insertion wild-type sequence in exon 3 of NKX2-1 (ENST00000354822.7:c.556_557insAlu541_556dup) segregated with the disease phenotype.

CONCLUSIONS

We identified a de novo exonic AluYa5 insertion causing NKX2-1-related disorder in SS/WES/WGS-negative cases, suggesting that MEI analysis of short-read sequencing data or targeted long-read sequencing could unmask the molecular diagnosis of unsolved MD cases. © 2022 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Heterozygous UCHL1 loss-of-function variants cause a neurodegenerative disorder with spasticity, ataxia, neuropathy, and optic atrophy

PURPOSE

Biallelic variants in UCHL1 have been associated with a progressive early-onset neurodegenerative disorder, autosomal recessive spastic paraplegia type 79. In this study, we investigated heterozygous UCHL1 variants on the basis of results from cohort-based burden analyses.

METHODS

Gene-burden analyses were performed on exome and genome data of independent cohorts of patients with hereditary ataxia and spastic paraplegia from Germany and the United Kingdom in a total of 3169 patients and 33,141 controls. Clinical data of affected individuals and additional independent families were collected and evaluated. Patients' fibroblasts were used to perform mass spectrometry-based proteomics.

RESULTS

UCHL1 was prioritized in both independent cohorts as a candidate gene for an autosomal dominant disorder. We identified a total of 34 cases from 18 unrelated families, carrying 13 heterozygous loss-of-function variants (15 families) and an inframe insertion (3 families). Affected individuals mainly presented with spasticity (24/31), ataxia (28/31), neuropathy (11/21), and optic atrophy (9/17). The mass spectrometry-based proteomics showed approximately 50% reduction of UCHL1 expression in patients' fibroblasts.

CONCLUSION

Our bioinformatic analysis, in-depth clinical and genetic workup, and functional studies established haploinsufficiency of UCHL1 as a novel disease mechanism in spastic ataxia.

Exome-wide analysis of copy number variation shows association of the human leukocyte antigen region with asthma in UK Biobank

BACKGROUND

The role of copy number variants (CNVs) in susceptibility to asthma is not well understood. This is, in part, due to the difficulty of accurately measuring CNVs in large enough sample sizes to detect associations. The recent availability of whole-exome sequencing (WES) in large biobank studies provides an unprecedented opportunity to study the role of CNVs in asthma.

METHODS

We called common CNVs in 49,953 individuals in the first release of UK Biobank WES using ClinCNV software. CNVs were tested for association with asthma in a stage 1 analysis comprising 7098 asthma cases and 36,578 controls from the first release of sequencing data. Nominally-associated CNVs were then meta-analysed in stage 2 with an additional 17,280 asthma cases and 115,562 controls from the second release of UK Biobank exome sequencing, followed by validation and fine-mapping.

RESULTS

Five of 189 CNVs were associated with asthma in stage 2, including a deletion overlapping the HLA-DQA1 and HLA-DQB1 genes, a duplication of CHROMR/PRKRA, deletions within MUC22 and TAP2, and a duplication in FBRSL1. The HLA-DQA1, HLA-DQB1, MUC22 and TAP2 genes all reside within the human leukocyte antigen (HLA) region on chromosome 6. In silico analyses demonstrated that the deletion overlapping HLA-DQA1 and HLA-DQB1 is likely to be an artefact arising from under-mapping of reads from non-reference HLA haplotypes, and that the CHROMR/PRKRA and FBRSL1 duplications represent presence/absence of pseudogenes within the HLA region. Bayesian fine-mapping of the HLA region suggested that there are two independent asthma association signals. The variants with the largest posterior inclusion probability in the two credible sets were an amino acid change in HLA-DQB1 (glutamine to histidine at residue 253) and a multi-allelic amino acid change in HLA-DRB1 (presence/absence of serine, glycine or leucine at residue 11).

CONCLUSIONS

At least two independent loci characterised by amino acid changes in the HLA-DQA1, HLA-DQB1 and HLA-DRB1 genes are likely to account for association of SNPs and CNVs in this region with asthma. The high divergence of haplotypes in the HLA can give rise to spurious CNVs, providing an important, cautionary tale for future large-scale analyses of sequencing data.

Transcript-Specific Loss-of-Function Variants in VPS16 Are Enriched in Patients With Dystonia

Background and Objectives

Our objective was to improve rare variant interpretation using statistical measures as well as publicly accessible annotation of expression levels and tissue specificity of different splice isoforms. We describe rare VPS16 variants observed in patients with dystonia and patients without dystonia, elaborate on our interpretation of VPS16 variants affecting different transcripts, and provide detailed clinical description of the movement disorder caused by VPS16 variants.

Methods

In-house exome and genome data sets (n = 11,539) were screened for rare heterozygous missense and putative loss-of-function (pLoF) variants in VPS16. Using pext (proportion expressed across transcripts) values from the Genome Aggregation Database (gnomAD), we differentiated variants affecting weakly and highly expressed exons/transcripts and applied statistical measures to systematically identify disease-associated genetic variation among patients with dystonia (n = 280).

Results

Six different heterozygous pLoFs in VPS16 transcripts were identified in 13 individuals. Three of these pLoFs occurred in 9 individuals with different phenotypes, and 3 pLoFs were identified in 4 unrelated individuals with early-onset dystonia. Although pLoFs were enriched in the dystonia cohort (n = 280; p = 2.04 × 10⁻⁴; 4/280 cases vs 9/11,259 controls; Fisher exact test), it was not exome-wide significant. According to the pext values in gnomAD, all 3 pLoFs observed in the patients with dystonia were located in the highly expressed canonical transcript ENST00000380445.3, whereas 2 of 3 pLoFs detected in 8 individuals without dystonia were located in the first exon of the noncanonical transcript ENST00000380443.3 that is weakly expressed across all tissues. Taking these biological implications into account, pLoFs involving the canonical transcript were exome-wide significantly enriched in patients with dystonia (p = 1.67 × 10⁻⁶; 4/280 cases vs 1/11,259 controls; Fisher exact test). All VPS16 patients showed mild progressive dystonia with writer's cramp as the presenting symptom between age 7 and 34 years (mean 20 years) that often progressed to generalized dystonia and was even accompanied by hyperkinetic movements and myoclonus in 1 patient.

Discussion

Our data provide strong evidence for VPS16 pLoFs to be implicated in dystonia and knowledge on exon resolution expression levels as well as statistical measures proved to be useful for variant interpretation.

Enhanced epigenetic profiling of classical human monocytes reveals a specific signature of healthy aging in the DNA methylome

The impact of healthy aging on molecular programming of immune cells is poorly understood. Here, we report comprehensive characterization of healthy aging in human classical monocytes, with a focus on epigenomic, transcriptomic, and proteomic alterations, as well as the corresponding proteomic and metabolomic data for plasma, using healthy cohorts of 20 young and 20 older males (~27 and ~64 years old on average). For each individual, we performed eRRBS-based DNA methylation profiling, which allowed us to identify a set of age-associated differentially methylated regions (DMRs) – a novel, cell-type specific signature of aging in DNA methylome. Hypermethylation events were associated with H3K27me3 in the CpG islands near promoters of lowly-expressed genes, while hypomethylated DMRs were enriched in H3K4me1 marked regions and associated with age-related increase of expression of the corresponding genes, providing a link between DNA methylation and age-associated transcriptional changes in primary human cells.

Detection of mobile elements insertions for routine clinical diagnostics in targeted sequencing data

Background

Targeted sequencing approaches such as gene panel or exome sequencing have become standard of care for the diagnosis of rare and common genetic disease. The detection and interpretation of point mutations, small insertions and deletions, and even exon‐level copy number variants are well established in clinical genetic testing. Other types of genetic variation such as mobile elements insertions (MEIs) are technically difficult to detect. In addition, their downstream clinical interpretation is more complex compared to point mutations due to a larger genomic footprint that can not only predict a clear loss of protein function but might disturb gene regulation and splicing even when located within the non‐coding regions. As a consequence, the contribution of MEIs to disease and tumor development remains largely unexplored in routine diagnostics.

Methods

In this study, we investigated the occurrence of MEIs in 7,693 exome datasets from individuals with rare diseases and healthy relatives as well as 788 cancer patients analyzed by panel sequencing.

Results

We present several exemplary cases highlighting the diagnostic value of MEIs and propose a strategy for the detection, prioritization, and clinical interpretation of MEIs in routine clinical diagnostics.

Conclusion

In this paper, we state that detection and interpretation of MEIs in clinical practice in targeted NGS data can be performed relatively easy despite the fact that MEIs very rarely occur in coding parts of the human genome. Large scale reanalysis of MEIs in existing cohorts may solve otherwise unsolvable cases.

New euchromatic variant dup(11)(p15.3p15.1) transmitted through two generations defined by low coverage whole genome sequencing

We report on a 14-year old boy, his father, and his paternal uncle, all three carriers of a duplication of chromosomal region 11p15.3-p15.1. The aberration was transmitted by the grandmother, who is carrier of a balanced insertion 46,XX,ins(14;11)(q32.1;p15.3p15.1). In order to determine the precise molecular basis of this structural variant, we performed low-coverage whole genome sequencing on the boy's father. This approach allowed precise determination of the genomic breakpoints and revealed a duplication of 6.9 Mb, centromeric to the Beckwith-Wiedemann/Silver-Russell syndrome critical region in 11p15.5, that inserted in inverse orientation into 14q32.12 (according to HGVS nomenclature: NC_000014.8:g.92871000_92871001ins[NC_000011.9:g.12250642_19165928inv;T]). To our knowledge, this is the first report of a duplication of 11p15.3-p15.1 involving more than 40 genes and transmitted through two generations without apparent clinical effects.

Efficient and flexible Integration of variant characteristics in rare variant association studies using integrated nested Laplace approximation

Rare variants are thought to play an important role in the etiology of complex diseases and may explain a significant fraction of the missing heritability in genetic disease studies. Next-generation sequencing facilitates the association of rare variants in coding or regulatory regions with complex diseases in large cohorts at genome-wide scale. However, rare variant association studies (RVAS) still lack power when cohorts are small to medium-sized and if genetic variation explains a small fraction of phenotypic variance. Here we present a novel Bayesian rare variant Association Test using Integrated Nested Laplace Approximation (BATI). Unlike existing RVAS tests, BATI allows integration of individual or variant-specific features as covariates, while efficiently performing inference based on full model estimation. We demonstrate that BATI outperforms established RVAS methods on realistic, semi-synthetic whole-exome sequencing cohorts, especially when using meaningful biological context, such as functional annotation. We show that BATI achieves power above 70% in scenarios in which competing tests fail to identify risk genes, e.g. when risk variants in sum explain less than 0.5% of phenotypic variance. We have integrated BATI, together with five existing RVAS tests in the 'Rare Variant Genome Wide Association Study' (rvGWAS) framework for data analyzed by whole-exome or whole genome sequencing. rvGWAS supports rare variant association for genes or any other biological unit such as promoters, while allowing the analysis of essential functionalities like quality control or filtering. Applying rvGWAS to a Chronic Lymphocytic Leukemia study we identified eight candidate predisposition genes, including EHMT2 and COPS7A.

Radiogenomics in head and neck cancer: correlation of radiomic heterogeneity and somatic mutations in TP53, FAT1 and KMT2D

PURPOSE

Genetic tumour profiles and radiomic features can be used to complement clinical information in head and neck squamous cell carcinoma (HNSCC) patients. Radiogenomics imply the potential to investigate complementarity or interrelations of radiomic and genomic features, and prognostic factors might be determined. The aim of our study was to explore radiogenomics in HNSCC.

METHODS

For 20 HNSCC patients treated with primary radiochemotherapy, next-generation sequencing (NGS) of tumour and corresponding normal tissue was performed. In total, 327 genes were investigated by panel sequencing. Radiomic features were extracted from computed tomography data. A hypothesis-driven approach was used for radiogenomic correlations of selected image-based heterogeneity features and well-known driver gene mutations in HNSCC.

RESULTS

The most frequently mutated driver genes in our cohort were TP53 (involved in cell cycle control), FAT1 (Wnt signalling, cell-cell contacts, migration) and KMT2D (chromatin modification). Radiomic features of heterogeneity did not correlate significantly with somatic mutations in TP53 or KMT2D. However, somatic mutations in FAT1 and smaller primary tumour volumes were associated with reduced radiomic intra-tumour heterogeneity.

CONCLUSION

The landscape of somatic variants in our cohort is well in line with previous reports. An association of somatic mutations in FAT1 with reduced radiomic tumour heterogeneity could potentially elucidate the previously described favourable outcomes of these patients. Larger studies are needed to validate this exploratory data in the future.

Allele balance bias identifies systematic genotyping errors and false disease associations

In recent years, next‐generation sequencing (NGS) has become a cornerstone of clinical genetics and diagnostics. Many clinical applications require high precision, especially if rare events such as somatic mutations in cancer or genetic variants causing rare diseases need to be identified. Although random sequencing errors can be modeled statistically and deep sequencing minimizes their impact, systematic errors remain a problem even at high depth of coverage. Understanding their source is crucial to increase precision of clinical NGS applications. In this work, we studied the relation between recurrent biases in allele balance (AB), systematic errors, and false positive variant calls across a large cohort of human samples analyzed by whole exome sequencing (WES). We have modeled the AB distribution for biallelic genotypes in 987 WES samples in order to identify positions recurrently deviating significantly from the expectation, a phenomenon we termed allele balance bias (ABB). Furthermore, we have developed a genotype callability score based on ABB for all positions of the human exome, which detects false positive variant calls that passed state‐of‐the‐art filters. Finally, we demonstrate the use of ABB for detection of false associations proposed by rare variant association studies. Availability: https://github.com/Francesc-Muyas/ABB.

A statistical approach to detection of copy number variations in PCR-enriched targeted sequencing data

BACKGROUND

Multiplex polymerase chain reaction (PCR) is a common enrichment technique for targeted massive parallel sequencing (MPS) protocols. MPS is widely used in biomedical research and clinical diagnostics as the fast and accurate tool for the detection of short genetic variations. However, identification of larger variations such as structure variants and copy number variations (CNV) is still being a challenge for targeted MPS. Some approaches and tools for structural variants detection were proposed, but they have limitations and often require datasets of certain type, size and expected number of amplicons affected by CNVs. In the paper, we describe novel algorithm for high-resolution germinal CNV detection in the PCR-enriched targeted sequencing data and present accompanying tool.

RESULTS

We have developed a machine learning algorithm for the detection of large duplications and deletions in the targeted sequencing data generated with PCR-based enrichment step. We have performed verification studies and established the algorithm's sensitivity and specificity. We have compared developed tool with other available methods applicable for the described data and revealed its higher performance.

CONCLUSION

We showed that our method has high specificity and sensitivity for high-resolution copy number detection in targeted sequencing data using large cohort of samples.