An association test of the spatial distribution of rare missense variants within protein structures identifies Alzheimer's disease-related patterns.

More than 90% of genetic variants are rare in most modern sequencing studies, such as the Alzheimer's Disease Sequencing Project (ADSP) whole-exome sequencing (WES) data. Furthermore, 54% of the rare variants in ADSP WES are singletons. However, both single variant and unit-based tests are limited in their statistical power to detect an association between rare variants and phenotypes. To best use missense rare variants and investigate their biological effect, we examine their association with phenotypes in the context of protein structures. We developed a protein structure-based approach, protein optimized kernel evaluation of missense nucleotides (POKEMON), which evaluates rare missense variants based on their spatial distribution within a protein rather than their allele frequency. The hypothesis behind this test is that the three-dimensional spatial distribution of variants within a protein structure provides functional context to power an association test. POKEMON identified three candidate genes (TREM2, SORL1, and EXOC3L4) and another suggestive gene from the ADSP WES data. For TREM2 and SORL1, two known Alzheimer's disease (AD) genes, the signal from the spatial cluster is stable even if we exclude known AD risk variants, indicating the presence of additional low-frequency risk variants within these genes. EXOC3L4 is a novel AD risk gene that has a cluster of variants primarily shared by case subjects around the Sec6 domain. This cluster is also validated in an independent replication data set and a validation data set with a larger sample size.

Transmission of Carbapenem-Resistant Klebsiella pneumoniae in US Hospitals.

BACKGROUND: Carbapenem-resistant Klebsiella pneumoniae (CRKp) is the most prevalent carbapenem-resistant Enterobacterales in the United States. We evaluated CRKp clustering in patients in US hospitals. METHODS: From April 2016 to August 2017, 350 patients with clonal group 258 CRKp were enrolled in the Consortium on Resistance Against Carbapenems in Klebsiella and other Enterobacteriaceae, a prospective, multicenter, cohort study. A maximum likelihood tree was constructed using RAxML. Static clusters shared ≤21 single-nucleotide polymorphisms (SNP) and a most recent common ancestor. Dynamic clusters incorporated SNP distance, culture timing, and rates of SNP accumulation and transmission using the R program TransCluster. RESULTS: Most patients were admitted from home (n = 150, 43%) or long-term care facilities (n = 115, 33%). Urine (n = 149, 43%) was the most common isolation site. Overall, 55 static and 47 dynamics clusters were identified involving 210 of 350 (60%) and 194 of 350 (55%) patients, respectively. Approximately half of static clusters were identical to dynamic clusters. Static clusters consisted of 33 (60%) intrasystem and 22 (40%) intersystem clusters. Dynamic clusters consisted of 32 (68%) intrasystem and 15 (32%) intersystem clusters and had fewer SNP differences than static clusters (8 vs 9; P = .045; 95% confidence interval [CI]: -4 to 0). Dynamic intersystem clusters contained more patients than dynamic intrasystem clusters (median [interquartile range], 4 [2, 7] vs 2 [2, 2]; P = .007; 95% CI: -3 to 0). CONCLUSIONS: Widespread intrasystem and intersystem transmission of CRKp was identified in hospitalized US patients. Use of different methods for assessing genetic similarity resulted in only minor differences in interpretation.

Novel loci for Alzheimer's disease identified by a genome-wide association study in Ashkenazi Jews.

INTRODUCTION: Most Alzheimer's disease (AD) loci have been discovered in individuals with European ancestry (EA). METHODS: We applied principal component analysis using Gaussian mixture models and an Ashkenazi Jewish (AJ) reference genome-wide association study (GWAS) data set to identify Ashkenazi Jews ascertained in GWAS (n = 42,682), whole genome sequencing (WGS, n = 16,815), and whole exome sequencing (WES, n = 20,504) data sets. The association of AD was tested genome wide (GW) in the GWAS and WGS data sets and exome wide (EW) in all three data sets (EW). Gene-based analyses were performed using aggregated rare variants. RESULTS: In addition to apolipoprotein E (APOE), GW analyses (1355 cases and 1661 controls) revealed associations with TREM2 R47H (p = 9.66 × 10-9 ), rs541586606 near RAB3B (p = 5.01 × 10-8 ), and rs760573036 between SPOCK3 and ANXA10 (p = 6.32 × 10-8 ). In EW analyses (1504 cases and 2047 controls), study-wide significant association was observed with rs1003710 near SMAP2 (p = 1.91 × 10-7 ). A significant gene-based association was identified with GIPR (p = 7.34 × 10-7 ). DISCUSSION: Our results highlight the efficacy of founder populations for AD genetic studies.

Association of mitochondrial variants and haplogroups identified by whole exome sequencing with Alzheimer's disease.

INTRODUCTION: Findings regarding the association between mitochondrial DNA (mtDNA) variants and Alzheimer's disease (AD) are inconsistent. METHODS: We developed a pipeline for accurate assembly and variant calling in mitochondrial genomes embedded within whole exome sequences (WES) from 10,831 participants from the Alzheimer's Disease Sequencing Project (ADSP). Association of AD risk was evaluated with each mtDNA variant and variants located in 1158 nuclear genes related to mitochondrial function using the SCORE test. Gene-based tests were performed using SKAT-O. RESULTS: Analysis of 4220 mtDNA variants revealed study-wide significant association of AD with a rare MT-ND4L variant (rs28709356 C>T; minor allele frequency = 0.002; P = 7.3 × 10-5 ) as well as with MT-ND4L in a gene-based test (P = 6.71 × 10-5 ). Significant association was also observed with a MT-related nuclear gene, TAMM41, in a gene-based test (P = 2.7 × 10-5 ). The expression of TAMM41 was lower in AD cases than controls (P = .00046) or mild cognitive impairment cases (P = .03). DISCUSSION: Significant findings in MT-ND4L and TAMM41 provide evidence for a role of mitochondria in AD.

Whole exome sequencing study identifies novel rare and common Alzheimer's-Associated variants involved in immune response and transcriptional regulation.

The Alzheimer's Disease Sequencing Project (ADSP) undertook whole exome sequencing in 5,740 late-onset Alzheimer disease (AD) cases and 5,096 cognitively normal controls primarily of European ancestry (EA), among whom 218 cases and 177 controls were Caribbean Hispanic (CH). An age-, sex- and APOE based risk score and family history were used to select cases most likely to harbor novel AD risk variants and controls least likely to develop AD by age 85 years. We tested ~1.5 million single nucleotide variants (SNVs) and 50,000 insertion-deletion polymorphisms (indels) for association to AD, using multiple models considering individual variants as well as gene-based tests aggregating rare, predicted functional, and loss of function variants. Sixteen single variants and 19 genes that met criteria for significant or suggestive associations after multiple-testing correction were evaluated for replication in four independent samples; three with whole exome sequencing (2,778 cases, 7,262 controls) and one with genome-wide genotyping imputed to the Haplotype Reference Consortium panel (9,343 cases, 11,527 controls). The top findings in the discovery sample were also followed-up in the ADSP whole-genome sequenced family-based dataset (197 members of 42 EA families and 501 members of 157 CH families). We identified novel and predicted functional genetic variants in genes previously associated with AD. We also detected associations in three novel genes: IGHG3 (p = 9.8 × 10-7), an immunoglobulin gene whose antibodies interact with ß-amyloid, a long non-coding RNA AC099552.4 (p = 1.2 × 10-7), and a zinc-finger protein ZNF655 (gene-based p = 5.0 × 10-6). The latter two suggest an important role for transcriptional regulation in AD pathogenesis.

Correction: Whole exome sequencing study identifies novel rare and common Alzheimer's-Associated variants involved in immune response and transcriptional regulation.

A correction to this paper has been published and can be accessed via a link at the top of the paper.

PopCluster: an algorithm to identify genetic variants with ethnicity-dependent effects.

MOTIVATION: Over the last decade, more diverse populations have been included in genome-wide association studies. If a genetic variant has a varying effect on a phenotype in different populations, genome-wide association studies applied to a dataset as a whole may not pinpoint such differences. It is especially important to be able to identify population-specific effects of genetic variants in studies that would eventually lead to development of diagnostic tests or drug discovery. RESULTS: In this paper, we propose PopCluster: an algorithm to automatically discover subsets of individuals in which the genetic effects of a variant are statistically different. PopCluster provides a simple framework to directly analyze genotype data without prior knowledge of subjects' ethnicities. PopCluster combines logistic regression modeling, principal component analysis, hierarchical clustering and a recursive bottom-up tree parsing procedure. The evaluation of PopCluster suggests that the algorithm has a stable low false positive rate (∼4%) and high true positive rate (>80%) in simulations with large differences in allele frequencies between cases and controls. Application of PopCluster to data from genetic studies of longevity discovers ethnicity-dependent heterogeneity in the association of rs3764814 (USP42) with the phenotype. AVAILABILITY AND IMPLEMENTATION: PopCluster was implemented using the R programming language, PLINK and Eigensoft software, and can be found at the following GitHub repository: https://github.com/gurinovich/PopCluster with instructions on its installation and usage. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Pseudomonas poae-Associated Fatal Septic Transfusion Reaction, Peoria, Illinois, USA, 2017.

In the United States, fatal transfusion-transmitted infections from red blood cell units are rare. Although this pattern mostly reflects how inhospitable refrigerated red blood cell units are to contaminant growth, fatalities caused by microorganisms that can grow at storage temperature (4°C), but not in standard clinical blood cultures at 37°C, are probably underestimated. We analyzed a fatal red blood cell transfusion in Peoria, Illinois, USA, that occurred in 2017. Samples from the patient's whole blood and the red blood cell unit remained culture-negative during the investigation, despite direct visualization of gram-negative bacilli within the unit immediately after transfusion. We identified the bacteria as Pseudomonas poae, a nonpathogenic pseudomonad carrying multiple cold-shock domain protein genes, and confirmed its cold tolerance and inability to grow at 37°C. Our work indicates transfusion reaction workups need to include testing for psychrophilic organisms, which could explain the cause of other apparently culture-negative transfusion reactions.

A rare missense variant of CASP7 is associated with familial late-onset Alzheimer's disease.

INTRODUCTION: The genetic architecture of Alzheimer's disease (AD) is only partially understood. METHODS: We conducted an association study for AD using whole sequence data from 507 genetically enriched AD cases (i.e., cases having close relatives affected by AD) and 4917 cognitively healthy controls of European ancestry (EA) and 172 enriched cases and 179 controls of Caribbean Hispanic ancestry. Confirmation of top findings from stage 1 was sought in two family-based genome-wide association study data sets and in a whole genome-sequencing data set comprising members from 42 EA and 115 Caribbean Hispanic families. RESULTS: We identified associations in EAs with variants in 12 novel loci. The most robust finding is a rare CASP7 missense variant (rs116437863; P = 2.44 × 10-10) which improved when combined with results from stage 2 data sets (P = 1.92 × 10-10). DISCUSSION: Our study demonstrated that an enriched case design can strengthen genetic signals, thus allowing detection of associations that would otherwise be missed in a traditional case-control study.

A long noncoding RNA from the HBS1L-MYB intergenic region on chr6q23 regulates human fetal hemoglobin expression.

The HBS1L-MYB intergenic region (chr6q23) regulates erythroid cell proliferation, maturation, and fetal hemoglobin (HbF) expression. An enhancer element within this locus, highlighted by a 3-bp deletion polymorphism (rs66650371), is known to interact with the promoter of the neighboring gene, MYB, to increase its expression, thereby regulating HbF production. RNA polymerase II binding and a 50-bp transcript from this enhancer region reported in ENCODE datasets suggested the presence of a long noncoding RNA (lncRNA). We characterized a novel 1283bp transcript (HMI-LNCRNA; chr6:135,096,362-135,097,644; hg38) that was transcribed from the enhancer region of MYB. Within erythroid cells, HMI-LNCRNA was almost exclusively present in nucleus, and was much less abundant than the mRNA for MYB. HMI-LNCRNA expression was significantly higher in erythroblasts derived from cultured adult peripheral blood CD34+ cells which expressed more HBB, compared to erythroblasts from cultured cord blood CD34+ cells which expressed much more HBG. Down-regulation of HMI-LNCRNA in HUDEP-2 cells, which expressed mostly HBB, significantly upregulated HBG expression both at the mRNA (200-fold) and protein levels, and promoted erythroid maturation. No change was found in the expression of BCL11A and other key transcription factors known to modulate HBG expression. HMI-LNCRNA plays an important role in regulating HBG expression, and its downregulation can result in a significant increase in HbF. HMI-LNCRNA might be a potential therapeutic target for HbF induction treatment in sickle cell disease and ß-thalassemia.

A phased SNP-based classification of sickle cell anemia HBB haplotypes.

BACKGROUND: Sickle cell anemia causes severe complications and premature death. Five common ß-globin gene cluster haplotypes are each associated with characteristic fetal hemoglobin (HbF) levels. As HbF is the major modulator of disease severity, classifying patients according to haplotype is useful. The first method of haplotype classification used restriction fragment length polymorphisms (RFLPs) to detect single nucleotide polymorphisms (SNPs) in the ß-globin gene cluster. This is labor intensive, and error prone. METHODS: We used genome-wide SNP data imputed to the 1000 Genomes reference panel to obtain phased data distinguishing parental alleles. RESULTS: We successfully haplotyped 813 sickle cell anemia patients previously classified by RFLPs with a concordance >98%. Four SNPs (rs3834466, rs28440105, rs10128556, and rs968857) marking four different restriction enzyme sites unequivocally defined most haplotypes. We were able to assign a haplotype to 86% of samples that were either partially or misclassified using RFLPs. CONCLUSION: Phased data using only four SNPs allowed unequivocal assignment of a haplotype that was not always possible using a larger number of RFLPs. Given the availability of genome-wide SNP data, our method is rapid and does not require high computational resources.

The genetic basis of asymptomatic codon 8 frame-shift (HBB:c25_26delAA) ß(0) -thalassaemia homozygotes.

Two 21-year old dizygotic twin men of Iraqi descent were homozygous for HBB codon 8, deletion of two nucleotides (-AA) frame-shift ß(0) -thalassaemia mutation (FSC8; HBB:c25_26delAA). Both were clinically well, had splenomegaly, and were never transfused. They had mild microcytic anaemia (Hb 120-130 g/l) and 98% of their haemoglobin was fetal haemoglobin (HbF). Both were carriers of Hph α-thalassaemia mutation. On the three major HbF quantitative trait loci (QTL), the twins were homozygous for G>A HBG2 Xmn1 site at single nucleotide polymorphism (SNP) rs7482144, homozygous for 3-bp deletion HBS1L-MYB intergenic polymorphism (HMIP) at rs66650371, and heterozygous for the A>C BCL11A intron 2 polymorphism at rs766432. These findings were compared with those found in 22 other FSC8 homozygote patients: four presented with thalassaemia intermedia phenotype, and 18 were transfusion dependent. The inheritance of homozygosity for HMIP 3-bp deletion at rs66650371 and heterozygosity for Hph α-thalassaemia mutation was found in the twins and not found in any of the other 22 patients. Further studies are needed to uncover likely additional genetic variants that could contribute to the exceptionally high HbF levels and mild phenotype in these twins.

A candidate transacting modulator of fetal hemoglobin gene expression in the Arab-Indian haplotype of sickle cell anemia.

Fetal hemoglobin (HbF) levels are higher in the Arab-Indian (AI) ß-globin gene haplotype of sickle cell anemia compared with African-origin haplotypes. To study genetic elements that effect HbF expression in the AI haplotype we completed whole genome sequencing in 14 Saudi AI haplotype sickle hemoglobin homozygotes-seven selected for low HbF (8.2% ± 1.3%) and seven selected for high HbF (23.5% ± 2.6%). An intronic single nucleotide polymorphism (SNP) in ANTXR1, an anthrax toxin receptor (chromosome 2p13), was associated with HbF. These results were replicated in two independent Saudi AI haplotype cohorts of 120 and 139 patients, but not in 76 Saudi Benin haplotype, 894 African origin haplotype and 44 AI haplotype patients of Indian origin, suggesting that this association is effective only in the Saudi AI haplotype background. ANTXR1 variants explained 10% of the HbF variability compared with 8% for BCL11A. These two genes had independent, additive effects on HbF and together explained about 15% of HbF variability in Saudi AI sickle cell anemia patients. ANTXR1 was expressed at mRNA and protein levels in erythroid progenitors derived from induced pluripotent stem cells (iPSCs) and CD34+ cells. As CD34+ cells matured and their HbF decreased ANTXR1 expression increased; as iPSCs differentiated and their HbF increased, ANTXR1 expression decreased. Along with elements in cis to the HbF genes, ANTXR1 contributes to the variation in HbF in Saudi AI haplotype sickle cell anemia and is the first gene in trans to HBB that is associated with HbF only in carriers of the Saudi AI haplotype. Am. J. Hematol. 91:1118-1122, 2016. © 2016 Wiley Periodicals, Inc.

HLA-A*3101 and carbamazepine-induced hypersensitivity reactions in Europeans.

BACKGROUND: Carbamazepine causes various forms of hypersensitivity reactions, ranging from maculopapular exanthema to severe blistering reactions. The HLA-B*1502 allele has been shown to be strongly correlated with carbamazepine-induced Stevens-Johnson syndrome and toxic epidermal necrolysis (SJS-TEN) in the Han Chinese and other Asian populations but not in European populations. METHODS: We performed a genomewide association study of samples obtained from 22 subjects with carbamazepine-induced hypersensitivity syndrome, 43 subjects with carbamazepine-induced maculopapular exanthema, and 3987 control subjects, all of European descent. We tested for an association between disease and HLA alleles through proxy single-nucleotide polymorphisms and imputation, confirming associations by high-resolution sequence-based HLA typing. We replicated the associations in samples from 145 subjects with carbamazepine-induced hypersensitivity reactions. RESULTS: The HLA-A*3101 allele, which has a prevalence of 2 to 5% in Northern European populations, was significantly associated with the hypersensitivity syndrome (P=3.5×10(-8)). An independent genomewide association study of samples from subjects with maculopapular exanthema also showed an association with the HLA-A*3101 allele (P=1.1×10(-6)). Follow-up genotyping confirmed the variant as a risk factor for the hypersensitivity syndrome (odds ratio, 12.41; 95% confidence interval [CI], 1.27 to 121.03), maculopapular exanthema (odds ratio, 8.33; 95% CI, 3.59 to 19.36), and SJS-TEN (odds ratio, 25.93; 95% CI, 4.93 to 116.18). CONCLUSIONS: The presence of the HLA-A*3101 allele was associated with carbamazepine-induced hypersensitivity reactions among subjects of Northern European ancestry. The presence of the allele increased the risk from 5.0% to 26.0%, whereas its absence reduced the risk from 5.0% to 3.8%. (Funded by the U.K. Department of Health and others.).

Direct detection of indirect transmission of Streptobacillus moniliformis rat bite fever infection.

We describe the evaluation of culture-negative synovial fluid from a 3-year-old boy by PCR and electrospray ionization followed by mass spectrometry (PCR/ESI-MS). Our patient developed a diffuse rash and fever with systemic signs and symptoms of sepsis, but four sets of blood cultures obtained prior to initiation of antibiotics were negative. After 1 week of illness, he developed right-knee swelling. Analysis of synovial fluid was consistent with infection, but cultures of specimens obtained following initiation of antimicrobial treatment were negative for growth. PCR/ESI-MS detected Streptobacillus moniliformis in the synovial fluid sample. Our patient completed an appropriate course of antibiotic treatment and remained completely asymptomatic in follow-up evaluation. This unique case suggests that PCR/ESI-MS may be a useful diagnostic tool for direct detection of unusual or unexpected pathogens directly from clinical specimens, particularly when samples have been obtained from patients following initiation of antibiotic therapy.

Ureaplasma parvum prosthetic joint infection detected by PCR.

We describe the first reported case of Ureaplasma parvum prosthetic joint infection (PJI) detected by PCR. Ureaplasma species do not possess a cell wall and are usually associated with colonization and infection of mucosal surfaces (not prosthetic material). U. parvum is a relatively new species name for certain serovars of Ureaplasma urealyticum, and PCR is useful for species determination. Our patient presented with late infection of his right total knee arthroplasty. Intraoperative fluid and tissue cultures and pre- and postoperative synovial fluid cultures were all negative. To discern the pathogen, we employed PCR coupled with electrospray ionization mass spectrometry (PCR/ESI-MS). Our patient's failure to respond to empirical antimicrobial treatment and our previous experience with PCR/ESI-MS in culture-negative cases of infection prompted us to use this approach over other diagnostic modalities. PCR/ESI-MS detected U. parvum in all samples. U. parvum-specific PCR testing was performed on all synovial fluid samples to confirm the U. parvum detection.

Identification of occult Fusobacterium nucleatum central nervous system infection by use of PCR-electrospray ionization mass spectrometry.

Anaerobic bacteria are often difficult to detect, especially after the initiation of antibiotics. We describe the application of PCR-electrospray ionization mass spectrometry (PCR/ESI-MS) using a sample of cerebrospinal fluid to identify an anaerobic Gram-negative bacillus, Fusobacterium nucleatum, in a patient with "culture-negative" meningitis and cerebral abscesses.

Culture negative empyema in a critically ill child: an opportunity for rapid molecular diagnostics.

BACKGROUND: Nucleic acid amplification technologies (NAAT) are advancing our ability to make rapid molecular diagnoses in patients with serious culture negative infections. This is the first report of PCR coupled to electrospray ionization mass spectrometry use in the evaluation of complicated community acquired pneumonia in a pediatric patient. CASE PRESENTATION: We present a case of culture negative empyema in a critically ill, Caucasian, 2-year-old girl who was treated with broad-spectrum empiric antibiotics, in which the length of stay was prolonged by adverse effects of the empiric antibiotic treatment. PCR coupled to electrospray ionization mass spectrometry was applied to culture negative fluid and tissue samples from the patient in order to determine the etiology of the empyema. CONCLUSIONS: Using this method, Streptococcus mitis/viridans was identified as the pathogen. A retrospective review of cases of empyema in children at our institution found that 87.5% of cases were negative for identification of a pathogen and antibiotics were administered to 100% of cases prior to collecting pleural fluid for culture. Understanding the role of Streptococcus mitis/viridans group in the etiology of empyema using an advanced NAAT coupled with mass spectrometry can enlighten clinicians as to the impact of this pathogen in community acquired pneumonia and help assist with antibiotic stewardship.

MitoH3: Mitochondrial Haplogroup and Homoplasmic/Heteroplasmic Variant Calling Pipeline for Alzheimer's Disease Sequencing Project.

Background: Mitochondrial DNA (mtDNA) is a double-stranded circular DNA and has multiple copies in each cell. Excess heteroplasmy, the coexistence of distinct variants in copies of mtDNA within a cell, may lead to mitochondrial impairments. Accurate determination of heteroplasmy in whole-genome sequencing (WGS) data has posed a significant challenge because mitochondria carrying heteroplasmic variants cannot be distinguished during library preparation. Moreover, sequencing errors, contamination, and nuclear mtDNA segments can reduce the accuracy of heteroplasmic variant calling. Objective: To efficiently and accurately call mtDNA homoplasmic and heteroplasmic variants from the large-scale WGS data generated from the Alzheimer's Disease Sequencing Project (ADSP), and test their association with Alzheimer's disease (AD). Methods: In this study, we present MitoH3-a comprehensive computational pipeline for calling mtDNA homoplasmic and heteroplasmic variants and inferring haplogroups in the ADSP WGS data. We first applied MitoH3 to 45 technical replicates from 6 subjects to define a threshold for detecting heteroplasmic variants. Then using the threshold of 5% ≤variant allele fraction≤95%, we further applied MitoH3 to call heteroplasmic variants from a total of 16,113 DNA samples with 6,742 samples from cognitively normal controls and 6,183 from AD cases. Results: This pipeline is available through the Singularity container engine. For 4,311 heteroplasmic variants identified from 16,113 samples, no significant variant count difference was observed between AD cases and controls. Conclusions: Our streamlined pipeline, MitoH3, enables computationally efficient and accurate analysis of a large number of samples.

A comparative study of structural variant calling in WGS from Alzheimer's disease families.

Detecting structural variants (SVs) in whole-genome sequencing poses significant challenges. We present a protocol for variant calling, merging, genotyping, sensitivity analysis, and laboratory validation for generating a high-quality SV call set in whole-genome sequencing from the Alzheimer's Disease Sequencing Project comprising 578 individuals from 111 families. Employing two complementary pipelines, Scalpel and Parliament, for SV/indel calling, we assessed sensitivity through sample replicates (N = 9) with in silico variant spike-ins. We developed a novel metric, D-score, to evaluate caller specificity for deletions. The accuracy of deletions was evaluated by Sanger sequencing. We generated a high-quality call set of 152,301 deletions of diverse sizes. Sanger sequencing validated 114 of 146 detected deletions (78.1%). Scalpel excelled in accuracy for deletions ≤100 bp, whereas Parliament was optimal for deletions >900 bp. Overall, 83.0% and 72.5% of calls by Scalpel and Parliament were validated, respectively, including all 11 deletions called by both Parliament and Scalpel between 101 and 900 bp. Our flexible protocol successfully generated a high-quality deletion call set and a truth set of Sanger sequencing-validated deletions with precise breakpoints spanning 1-17,000 bp.