Species-wide genetic incompatibility analysis identifies immune genes as hot spots of deleterious epistasis.

Intraspecific genetic incompatibilities prevent the assembly of specific alleles into single genotypes and influence genome- and species-wide patterns of sequence variation. A common incompatibility in plants is hybrid necrosis, characterized by autoimmune responses due to epistatic interactions between natural genetic variants. By systematically testing thousands of F1 hybrids of Arabidopsis thaliana strains, we identified a small number of incompatibility hot spots in the genome, often in regions densely populated by nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes. In several cases, these immune receptor loci interact with each other, suggestive of conflict within the immune system. A particularly dangerous locus is a highly variable cluster of NLR genes, DM2, which causes multiple independent incompatibilities with genes that encode a range of biochemical functions, including NLRs. Our findings suggest that deleterious interactions of immune receptors limit the combinations of favorable disease resistance alleles accessible to plant genomes.

Swarm Learning for decentralized and confidential clinical machine learning.

Fast and reliable detection of patients with severe and heterogeneous illnesses is a major goal of precision medicine1,2. Patients with leukaemia can be identified using machine learning on the basis of their blood transcriptomes3. However, there is an increasing divide between what is technically possible and what is allowed, because of privacy legislation4,5. Here, to facilitate the integration of any medical data from any data owner worldwide without violating privacy laws, we introduce Swarm Learning-a decentralized machine-learning approach that unites edge computing, blockchain-based peer-to-peer networking and coordination while maintaining confidentiality without the need for a central coordinator, thereby going beyond federated learning. To illustrate the feasibility of using Swarm Learning to develop disease classifiers using distributed data, we chose four use cases of heterogeneous diseases (COVID-19, tuberculosis, leukaemia and lung pathologies). With more than 16,400 blood transcriptomes derived from 127 clinical studies with non-uniform distributions of cases and controls and substantial study biases, as well as more than 95,000 chest X-ray images, we show that Swarm Learning classifiers outperform those developed at individual sites. In addition, Swarm Learning completely fulfils local confidentiality regulations by design. We believe that this approach will notably accelerate the introduction of precision medicine.

ZSCAN10 deficiency causes a neurodevelopmental disorder with characteristic oto-facial malformations.

Neurodevelopmental disorders are major indications for genetic referral and have been linked to more than 1500 loci including genes encoding transcriptional regulators. The dysfunction of transcription factors often results in characteristic syndromic presentations; however, at least half of these patients lack a genetic diagnosis. The implementation of machine learning approaches has the potential to aid in the identification of new disease genes and delineate associated phenotypes. Next generation sequencing was performed in seven affected individuals with neurodevelopmental delay and dysmorphic features. Clinical characterization included reanalysis of available neuroimaging datasets and 2D portrait image analysis with GestaltMatcher. The functional consequences of ZSCAN10 loss were modelled in mouse embryonic stem cells (mESCs), including a knockout and a representative ZSCAN10 protein truncating variant. These models were characterized by gene expression and western blot analyses, chromatin immunoprecipitation and quantitative PCR (ChIP-qPCR) and immunofluorescence staining. Zscan10 knockout mouse embryos were generated and phenotyped. We prioritized bi-allelic ZSCAN10 loss-of-function variants in seven affected individuals from five unrelated families as the underlying molecular cause. RNA-sequencing analyses in Zscan10-/- mESCs indicated dysregulation of genes related to stem cell pluripotency. In addition, we established in mESCs the loss-of-function mechanism for a representative human ZSCAN10 protein truncating variant by showing alteration of its expression levels and subcellular localization, interfering with its binding to DNA enhancer targets. Deep phenotyping revealed global developmental delay, facial asymmetry and malformations of the outer ear as consistent clinical features. Cerebral MRI showed dysplasia of the semicircular canals as an anatomical correlate of sensorineural hearing loss. Facial asymmetry was confirmed as a clinical feature by GestaltMatcher and was recapitulated in the Zscan10 mouse model along with inner and outer ear malformations. Our findings provide evidence of a novel syndromic neurodevelopmental disorder caused by bi-allelic loss-of-function variants in ZSCAN10.

Transcriptional control of gene expression by microRNAs.

MicroRNAs (miRNAs) control gene expression in animals and plants. Like another class of small RNAs, siRNAs, they affect gene expression posttranscriptionally. While siRNAs in addition act in transcriptional gene silencing, a role of miRNAs in transcriptional regulation has been less clear. We show here that in moss Physcomitrella patens mutants without a DICER-LIKE1b gene, maturation of miRNAs is normal but cleavage of target RNAs is abolished and levels of these transcripts are drastically reduced. These mutants accumulate miRNA:target-RNA duplexes and show hypermethylation of the genes encoding target RNAs, leading to gene silencing. This pathway occurs also in the wild-type upon hormone treatment. We propose that initiation of epigenetic silencing by DNA methylation depends on the ratio of the miRNA and its target RNA.

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.

Bi-allelic loss-of-function variants in BCAS3 cause a syndromic neurodevelopmental disorder.

BCAS3 microtubule-associated cell migration factor (BCAS3) is a large, highly conserved cytoskeletal protein previously proposed to be critical in angiogenesis and implicated in human embryogenesis and tumorigenesis. Here, we established BCAS3 loss-of-function variants as causative for a neurodevelopmental disorder. We report 15 individuals from eight unrelated families with germline bi-allelic loss-of-function variants in BCAS3. All probands share a global developmental delay accompanied by pyramidal tract involvement, microcephaly, short stature, strabismus, dysmorphic facial features, and seizures. The human phenotype is less severe compared with the Bcas3 knockout mouse model and cannot be explained by angiogenic defects alone. Consistent with being loss-of-function alleles, we observed absence of BCAS3 in probands' primary fibroblasts. By comparing the transcriptomic and proteomic data based on probands' fibroblasts with those of the knockout mouse model, we identified similar dysregulated pathways resulting from over-representation analysis, while the dysregulation of some proposed key interactors could not be confirmed. Together with the results from a tissue-specific Drosophila loss-of-function model, we demonstrate a vital role for BCAS3 in neural tissue development.

Bi-allelic loss-of-function variants in KIF21A cause severe fetal akinesia with arthrogryposis multiplex.

BACKGROUND: Fetal akinesia (FA) results in variable clinical presentations and has been associated with more than 166 different disease loci. However, the underlying molecular cause remains unclear in many individuals. We aimed to further define the set of genes involved. METHODS: We performed in-depth clinical characterisation and exome sequencing on a cohort of 23 FA index cases sharing arthrogryposis as a common feature. RESULTS: We identified likely pathogenic or pathogenic variants in 12 different established disease genes explaining the disease phenotype in 13 index cases and report 12 novel variants. In the unsolved families, a search for recessive-type variants affecting the same gene was performed; and in five affected fetuses of two unrelated families, a homozygous loss-of-function variant in the kinesin family member 21A gene (KIF21A) was found. CONCLUSION: Our study underlines the broad locus heterogeneity of FA with well-established and atypical genotype-phenotype associations. We describe KIF21A as a new factor implicated in the pathogenesis of severe neurogenic FA sequence with arthrogryposis of multiple joints, pulmonary hypoplasia and facial dysmorphisms. This hypothesis is further corroborated by a recent report on overlapping phenotypes observed in Kif21a null piglets.

Detection of circulating cell-free HPV DNA of 13 HPV types for patients with cervical cancer as potential biomarker to monitor therapy response and to detect relapse.

BACKGROUND: HPV-related cervical cancer (CC) is the fourth most frequent cancer in women worldwide. Cell-free tumour DNA is a potent biomarker to detect treatment response, residual disease, and relapse. We investigated the potential use of cell-free circulating HPV-DNA (cfHPV-DNA) in plasma of patients with CC. METHODS: cfHPV-DNA levels were measured using a highly sensitive next-generation sequencing-based approach targeting a panel of 13 high-risk HPV types. RESULTS: Sequencing was performed in 69 blood samples collected from 35 patients, of which 26 were treatment-naive when the first liquid biopsy sample was retrieved. cfHPV-DNA was successfully detected in 22/26 (85%) cases. A significant correlation between tumour burden and cfHPV-DNA levels was observed: cfHPV-DNA was detectable in all treatment-naive patients with advanced-stage disease (17/17, FIGO IB3-IVB) and in 5/9 patients with early-stage disease (FIGO IA-IB2). Sequential samples revealed a decrease of cfHPV-DNA levels in 7 patients corresponding treatment response and an increase in a patient with relapse. CONCLUSIONS: In this proof-of-concept study we demonstrated the potential of cfHPV-DNA as a biomarker for therapy monitoring in patients with primary and recurrent CC. Our findings facilitate the development of a sensitive and precise, non-invasive, inexpensive, and easily accessible tool in CC diagnosis, therapy monitoring and follow-up.

Genetic and functional analysis of chymotrypsin-like protease (CTRL) in chronic pancreatitis.

BACKGROUND: Genetic predisposition is crucial in the pathogenesis of early-onset chronic pancreatitis (CP). So far, several genetic alterations have been identified as risk factors, predominantly in genes encoding digestive enzymes. However, many early-onset CP cases have no identified underlying cause. Chymotrypsins are a family of serine proteases that can cleave trypsinogen and lead to its degradation. Because genetic alterations in the chymotrypsins CTRC, CTRB1, and CTRB2 are associated with CP, we genetically and functionally investigated chymotrypsin-like protease (CTRL) as a potential risk factor. METHODS: We screened 1005 non-alcoholic CP patients and 1594 controls for CTRL variants by exome sequencing. We performed Western blots and activity assays to analyse secretion and proteolytic activity. We measured BiP mRNA expression to investigate the potential impact of identified alterations on endoplasmic reticulum (ER) stress. RESULTS: We identified 13 heterozygous non-synonymous CTRL variants: five exclusively in patients and three only in controls. Functionality was unchanged in 6/13 variants. Four alterations showed normal secretion but reduced (p.G20S, p.G56S, p.G61S) or abolished (p.S208F) activity. Another three variants (p.C201Y, p.G215R and p.C220G) were not secreted and already showed reduced or no activity intracellularly. However, intracellular retention did not lead to ER stress. CONCLUSION: We identified several CTRL variants, some showing potent effects on protease function and secretion. We observed these effects in variants found in patients and controls, and CTRL loss-of-function variants were not significantly more common in patients than controls. Therefore, CTRL is unlikely to play a relevant role in the development of CP.

First case of adult onset neuronal intranuclear inclusion disease with both typical radiological signs and NOTCH2NLC repeat expansions in a Caucasian individual.

BACKGROUND AND PURPOSE: Adult onset neuronal intranuclear inclusion disease (NIID) is a rare neurodegenerative disorder with a heterogeneous clinical presentation that can mimic stroke and various forms of dementia. To date, it has been described almost exclusively in Asian individuals. METHODS: This case presentation includes magnetic resonance imaging (MRI) of the neurocranium, histology by skin biopsy, and long-read genome sequencing. RESULTS: A 75-year-old Caucasian female presented with paroxysmal encephalopathy twice within a 14-month period. Brain MRI revealed high-intensity signals at the cerebral corticomedullary junction (diffusion-weighted imaging) and the paravermal area (fluid-attenuated inversion recovery), a typical distribution observed in adult onset NIID. The diagnosis was corroborated by skin biopsy, which demonstrated eosinophilic intranuclear inclusion bodies, and confirmed by long-read genome sequencing, showing an expansion of the GGC repeat in exon 1 of NOTCH2NLC. CONCLUSIONS: Our case proves adult onset NOTCH2NLC-GGC-positive NIID with typical findings on MRI and histology in a Caucasian patient and underscores the need to consider this diagnosis in non-Asian individuals.

Chromatin-wide profiling of DYRK1A reveals a role as a gene-specific RNA polymerase II CTD kinase.

DYRK1A is a dosage-sensitive protein kinase that fulfills key roles during development and in tissue homeostasis, and its dysregulation results in human pathologies. DYRK1A is present in both the nucleus and cytoplasm of mammalian cells, although its nuclear function remains unclear. Genome-wide analysis of DYRK1A-associated loci reveals that the kinase is recruited preferentially to promoters of genes actively transcribed by RNA polymerase II (RNAPII), which are functionally associated with translation, RNA processing, and cell cycle. DYRK1A-bound promoter sequences are highly enriched in a conserved palindromic motif, which is necessary to drive DYRK1A-dependent transcriptional activation. DYRK1A phosphorylates the C-terminal domain (CTD) of RNAPII at Ser2 and Ser5. Depletion of DYRK1A results in reduced association of RNAPII at the target promoters as well as hypophosphorylation of the RNAPII CTD along the target gene bodies. These results are consistent with DYRK1A being a transcriptional regulator by acting as a CTD kinase.

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.

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.

A single center experience of prenatal parent-fetus trio exome sequencing for pregnancies with congenital anomalies.

OBJECTIVES: To examine the diagnostic yield of trio exome sequencing in fetuses with multiple structural defects with no pathogenic findings in cytogenetic and microarray analyses. METHODS: We recruited 51 fetuses with two or more defects, non-immune fetal hydrops or fetal akinesia deformation syndrome|or fetal akinesia deformation sequence (FADS). Trio exome sequencing was performed on DNA from chorionic villi samples and parental blood. Detection of genomic variation and prioritization of clinically relevant variants was performed according to in-house standard operating procedures. RESULTS: Median maternal and gestational age was 32.0 years and 21.0 weeks, respectively. Forty-three (84.3%) fetuses had two or more affected organ systems. The remaining fetuses had isolated fetal hydrops or FADS. In total, the exome analysis established the genetic cause for the clinical abnormalities in 22 (43.1%, 95% CI 29.4%-57.8%) pregnancies. CONCLUSIONS: In fetuses with multiple defects, hydrops or FADS and normal standard genetic results, trio exome sequencing has the potential to identify genetic anomalies in more than 40% of 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.

Mutations in DCHS1 cause mitral valve prolapse.

Mitral valve prolapse (MVP) is a common cardiac valve disease that affects nearly 1 in 40 individuals. It can manifest as mitral regurgitation and is the leading indication for mitral valve surgery. Despite a clear heritable component, the genetic aetiology leading to non-syndromic MVP has remained elusive. Four affected individuals from a large multigenerational family segregating non-syndromic MVP underwent capture sequencing of the linked interval on chromosome 11. We report a missense mutation in the DCHS1 gene, the human homologue of the Drosophila cell polarity gene dachsous (ds), that segregates with MVP in the family. Morpholino knockdown of the zebrafish homologue dachsous1b resulted in a cardiac atrioventricular canal defect that could be rescued by wild-type human DCHS1, but not by DCHS1 messenger RNA with the familial mutation. Further genetic studies identified two additional families in which a second deleterious DCHS1 mutation segregates with MVP. Both DCHS1 mutations reduce protein stability as demonstrated in zebrafish, cultured cells and, notably, in mitral valve interstitial cells (MVICs) obtained during mitral valve repair surgery of a proband. Dchs1(+/-) mice had prolapse of thickened mitral leaflets, which could be traced back to developmental errors in valve morphogenesis. DCHS1 deficiency in MVP patient MVICs, as well as in Dchs1(+/-) mouse MVICs, result in altered migration and cellular patterning, supporting these processes as aetiological underpinnings for the disease. Understanding the role of DCHS1 in mitral valve development and MVP pathogenesis holds potential for therapeutic insights for this very common disease.

PET/MRI and genetic intrapatient heterogeneity in head and neck cancers.

PURPOSE: The relation between functional imaging and intrapatient genetic heterogeneity remains poorly understood. The aim of our study was to investigate spatial sampling and functional imaging by FDG-PET/MRI to describe intrapatient tumour heterogeneity. METHODS: Six patients with oropharyngeal cancer were included in this pilot study. Two tumour samples per patient were taken and sequenced by next-generation sequencing covering 327 genes relevant in head and neck cancer. Corresponding regions were delineated on pretherapeutic FDG-PET/MRI images to extract apparent diffusion coefficients and standardized uptake values. RESULTS: Samples were collected within the primary tumour (n = 3), within the primary tumour and the involved lymph node (n = 2) as well as within two independent primary tumours (n = 1). Genetic heterogeneity of the primary tumours was limited and most driver gene mutations were found ubiquitously. Slightly increasing heterogeneity was found between primary tumours and lymph node metastases. One private predicted driver mutation within a primary tumour and one in a lymph node were found. However, the two independent primary tumours did not show any shared mutations in spite of a clinically suspected field cancerosis. No conclusive correlation between genetic heterogeneity and heterogeneity of PET/MRI-derived parameters was observed. CONCLUSION: Our limited data suggest that single sampling might be sufficient in some patients with oropharyngeal cancer. However, few driver mutations might be missed and, if feasible, spatial sampling should be considered. In two independent primary tumours, both lesions should be sequenced. Our data with a limited number of patients do not support the concept that multiparametric PET/MRI features are useful to guide biopsies for genetic tumour characterization.

Description of Citrobacter cronae sp. nov., isolated from human rectal swabs and stool samples.

Nine independent Gram-negative bacterial strains were isolated from rectal swabs or stool samples of immunocompromised patients from two different wards of a university hospital. All isolates were phylogenetically analysed based on their 16S rRNA gene sequence, housekeeping gene recN, multilocus sequence analysis of concatenated partial fusA, leuS, pyrG and rpoB sequences, and by whole genome sequencing data. The analysed strains of the new species cluster together and form a separate branch with Citrobacter werkmanii NBRC105721T as the most closely related species. An average nucleotide identity value of 95.9-96% and computation of digital DNA-DNA hybridization values separate the new species from all other type strains of the genus Citrobacter. Biochemical characteristics further delimit the isolates from closely related Citrobacter type strains. As a result of the described data, a new Citrobacter species is introduced, for which the name Citrobacter cronae sp. nov. is proposed. The type strain is Tue2-1T with a G+C DNA content of 52.2 mol%.

Tissue-specific DNA methylation loss during ageing and carcinogenesis is linked to chromosome structure, replication timing and cell division rates.

DNA methylation is an epigenetic mechanism known to affect gene expression and aberrant DNA methylation patterns have been described in cancer. However, only a small fraction of differential methylation events target genes with a defined role in cancer, raising the question of how aberrant DNA methylation contributes to carcinogenesis. As recently a link has been suggested between methylation patterns arising in ageing and those arising in cancer, we asked which aberrations are unique to cancer and which are the product of normal ageing processes. We therefore compared the methylation patterns between ageing and cancer in multiple tissues. We observed that hypermethylation preferentially occurs in regulatory elements, while hypomethylation is associated with structural features of the chromatin. Specifically, we observed consistent hypomethylation of late-replicating, lamina-associated domains. The extent of hypomethylation was stronger in cancer, but in both ageing and cancer it was proportional to the replication timing of the region and the cell division rate of the tissue. Moreover, cancer patients who displayed more hypomethylation in late-replicating, lamina-associated domains had higher expression of cell division genes. These findings suggest that different cell division rates contribute to tissue- and cancer type-specific DNA methylation profiles.

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.