The choroid plexus links innate immunity to CSF dysregulation in hydrocephalus.

The choroid plexus (ChP) is the blood-cerebrospinal fluid (CSF) barrier and the primary source of CSF. Acquired hydrocephalus, caused by brain infection or hemorrhage, lacks drug treatments due to obscure pathobiology. Our integrated, multi-omic investigation of post-infectious hydrocephalus (PIH) and post-hemorrhagic hydrocephalus (PHH) models revealed that lipopolysaccharide and blood breakdown products trigger highly similar TLR4-dependent immune responses at the ChP-CSF interface. The resulting CSF "cytokine storm", elicited from peripherally derived and border-associated ChP macrophages, causes increased CSF production from ChP epithelial cells via phospho-activation of the TNF-receptor-associated kinase SPAK, which serves as a regulatory scaffold of a multi-ion transporter protein complex. Genetic or pharmacological immunomodulation prevents PIH and PHH by antagonizing SPAK-dependent CSF hypersecretion. These results reveal the ChP as a dynamic, cellularly heterogeneous tissue with highly regulated immune-secretory capacity, expand our understanding of ChP immune-epithelial cell cross talk, and reframe PIH and PHH as related neuroimmune disorders vulnerable to small molecule pharmacotherapy.

The Psychiatric Cell Map Initiative: A Convergent Systems Biological Approach to Illuminating Key Molecular Pathways in Neuropsychiatric Disorders.

Although gene discovery in neuropsychiatric disorders, including autism spectrum disorder, intellectual disability, epilepsy, schizophrenia, and Tourette disorder, has accelerated, resulting in a large number of molecular clues, it has proven difficult to generate specific hypotheses without the corresponding datasets at the protein complex and functional pathway level. Here, we describe one path forward-an initiative aimed at mapping the physical and genetic interaction networks of these conditions and then using these maps to connect the genomic data to neurobiology and, ultimately, the clinic. These efforts will include a team of geneticists, structural biologists, neurobiologists, systems biologists, and clinicians, leveraging a wide array of experimental approaches and creating a collaborative infrastructure necessary for long-term investigation. This initiative will ultimately intersect with parallel studies that focus on other diseases, as there is a significant overlap with genes implicated in cancer, infectious disease, and congenital heart defects.

Pioneering a Global Cure for Chronic Hepatitis C Virus Infection.

This year's Lasker∼Debakey Clinical Medical Research Award honors Ralf Bartenschlager, Charles Rice, and Michael Sofia, pioneers in the development of curative and safe therapies for the 170 million people with hepatitis C virus infection.

Liver transplantation: from inception to clinical practice.

The 2012 Lasker-DeBakey Clinical Medical Research Award will be conferred on Thomas Starzl of the University of Pittsburgh School of Medicine in Pittsburgh, Pennsylvania, USA and Roy Calne of the University of Cambridge in Cambridge, UK. They are recognized for pioneering the development of liver transplantation, an intervention that saves 20,000 lives world-wide each year.

Pathogenic variants in autism gene KATNAL2 cause hydrocephalus and disrupt neuronal connectivity by impairing ciliary microtubule dynamics.

Enlargement of the cerebrospinal fluid (CSF)-filled brain ventricles (cerebral ventriculomegaly), the cardinal feature of congenital hydrocephalus (CH), is increasingly recognized among patients with autism spectrum disorders (ASD). KATNAL2, a member of Katanin family microtubule-severing ATPases, is a known ASD risk gene, but its roles in human brain development remain unclear. Here, we show that nonsense truncation of Katnal2 (Katnal2Δ17) in mice results in classic ciliopathy phenotypes, including impaired spermatogenesis and cerebral ventriculomegaly. In both humans and mice, KATNAL2 is highly expressed in ciliated radial glia of the fetal ventricular-subventricular zone as well as in their postnatal ependymal and neuronal progeny. The ventriculomegaly observed in Katnal2Δ17 mice is associated with disrupted primary cilia and ependymal planar cell polarity that results in impaired cilia-generated CSF flow. Further, prefrontal pyramidal neurons in ventriculomegalic Katnal2Δ17 mice exhibit decreased excitatory drive and reduced high-frequency firing. Consistent with these findings in mice, we identified rare, damaging heterozygous germline variants in KATNAL2 in five unrelated patients with neurosurgically treated CH and comorbid ASD or other neurodevelopmental disorders. Mice engineered with the orthologous ASD-associated KATNAL2 F244L missense variant recapitulated the ventriculomegaly found in human patients. Together, these data suggest KATNAL2 pathogenic variants alter intraventricular CSF homeostasis and parenchymal neuronal connectivity by disrupting microtubule dynamics in fetal radial glia and their postnatal ependymal and neuronal descendants. The results identify a molecular mechanism underlying the development of ventriculomegaly in a genetic subset of patients with ASD and may explain persistence of neurodevelopmental phenotypes in some patients with CH despite neurosurgical CSF shunting.

De novo variants implicate chromatin modification, transcriptional regulation, and retinoic acid signaling in syndromic craniosynostosis.

Craniosynostosis (CS) is the most common congenital cranial anomaly. Several Mendelian forms of syndromic CS are well described, but a genetic etiology remains elusive in a substantial fraction of probands. Analysis of exome sequence data from 526 proband-parent trios with syndromic CS identified a marked excess (observed 98, expected 33, p = 4.83 × 10-20) of damaging de novo variants (DNVs) in genes highly intolerant to loss-of-function variation (probability of LoF intolerance > 0.9). 30 probands harbored damaging DNVs in 21 genes that were not previously implicated in CS but are involved in chromatin modification and remodeling (4.7-fold enrichment, p = 1.1 × 10-11). 17 genes had multiple damaging DNVs, and 13 genes (CDK13, NFIX, ADNP, KMT5B, SON, ARID1B, CASK, CHD7, MED13L, PSMD12, POLR2A, CHD3, and SETBP1) surpassed thresholds for genome-wide significance. A recurrent gain-of-function DNV in the retinoic acid receptor alpha (RARA; c.865G>A [p.Gly289Arg]) was identified in two probands with similar CS phenotypes. CS risk genes overlap with those identified for autism and other neurodevelopmental disorders, are highly expressed in cranial neural crest cells, and converge in networks that regulate chromatin modification, gene transcription, and osteoblast differentiation. Our results identify several CS loci and have major implications for genetic testing and counseling.

Pleiotropic role of TRAF7 in skull-base meningiomas and congenital heart disease.

While somatic variants of TRAF7 (Tumor necrosis factor receptor-associated factor 7) underlie anterior skull-base meningiomas, here we report the inherited mutations of TRAF7 that cause congenital heart defects. We show that TRAF7 mutants operate in a dominant manner, inhibiting protein function via heterodimerization with wild-type protein. Further, the shared genetics of the two disparate pathologies can be traced to the common origin of forebrain meninges and cardiac outflow tract from the TRAF7-expressing neural crest. Somatic and inherited mutations disrupt TRAF7-IFT57 interactions leading to cilia degradation. TRAF7-mutant meningioma primary cultures lack cilia, and TRAF7 knockdown causes cardiac, craniofacial, and ciliary defects in Xenopus and zebrafish, suggesting a mechanistic convergence for TRAF7-driven meningiomas and developmental heart defects.

TRIM71 mutations cause a neurodevelopmental syndrome featuring ventriculomegaly and hydrocephalus.

Congenital hydrocephalus (CH), characterized by cerebral ventriculomegaly, is one of the most common reasons for pediatric brain surgery. Recent studies have implicated lin-41 (lineage variant 41)/TRIM71 (tripartite motif 71) as a candidate CH risk gene, however, TRIM71 variants have not been systematically examined in a large patient cohort or conclusively linked with an OMIM syndrome. Through cross-sectional analysis of the largest assembled cohort of patients with cerebral ventriculomegaly, including neurosurgically-treated CH (totaling 2,697 parent-proband trios and 8,091 total exomes), we identified 13 protein-altering de novo variants (DNVs) in TRIM71 in unrelated children exhibiting variable ventriculomegaly, CH, developmental delay, dysmorphic features, and other structural brain defects including corpus callosum dysgenesis and white matter hypoplasia. Eight unrelated patients were found to harbor arginine variants, including two recurrent missense DNVs, at homologous positions in RPXGV motifs of different NHL domains. Seven additional patients with rare, damaging, unphased or transmitted variants of uncertain significance were also identified. NHL-domain variants of TRIM71 exhibited impaired binding to the canonical TRIM71 target CDKN1A; other variants failed to direct the subcellular localization of TRIM71 to processing bodies. Single-cell transcriptomic analysis of human embryos revealed expression of TRIM71 in early first-trimester neural stem cells of the brain. These data show TRIM71 is essential for human brain morphogenesis and that TRIM71 mutations cause a novel neurodevelopmental syndrome featuring ventriculomegaly and CH.

A novel SMARCC1 BAFopathy implicates neural progenitor epigenetic dysregulation in human hydrocephalus.

Hydrocephalus, characterized by cerebral ventriculomegaly, is the most common disorder requiring brain surgery in children. Recent studies have implicated SMARCC1, a component of the BRG1-associated factor (BAF) chromatin remodelling complex, as a candidate congenital hydrocephalus gene. However, SMARCC1 variants have not been systematically examined in a large patient cohort or conclusively linked with a human syndrome. Moreover, congenital hydrocephalus-associated SMARCC1 variants have not been functionally validated or mechanistically studied in vivo. Here, we aimed to assess the prevalence of SMARCC1 variants in an expanded patient cohort, describe associated clinical and radiographic phenotypes, and assess the impact of Smarcc1 depletion in a novel Xenopus tropicalis model of congenital hydrocephalus. To do this, we performed a genetic association study using whole-exome sequencing from a cohort consisting of 2697 total ventriculomegalic trios, including patients with neurosurgically-treated congenital hydrocephalus, that total 8091 exomes collected over 7 years (2016-23). A comparison control cohort consisted of 1798 exomes from unaffected siblings of patients with autism spectrum disorder and their unaffected parents were sourced from the Simons Simplex Collection. Enrichment and impact on protein structure were assessed in identified variants. Effects on the human fetal brain transcriptome were examined with RNA-sequencing and Smarcc1 knockdowns were generated in Xenopus and studied using optical coherence tomography imaging, in situ hybridization and immunofluorescence. SMARCC1 surpassed genome-wide significance thresholds, yielding six rare, protein-altering de novo variants localized to highly conserved residues in key functional domains. Patients exhibited hydrocephalus with aqueductal stenosis; corpus callosum abnormalities, developmental delay, and cardiac defects were also common. Xenopus knockdowns recapitulated both aqueductal stenosis and cardiac defects and were rescued by wild-type but not patient-specific variant SMARCC1. Hydrocephalic SMARCC1-variant human fetal brain and Smarcc1-variant Xenopus brain exhibited a similarly altered expression of key genes linked to midgestational neurogenesis, including the transcription factors NEUROD2 and MAB21L2. These results suggest de novo variants in SMARCC1 cause a novel human BAFopathy we term 'SMARCC1-associated developmental dysgenesis syndrome', characterized by variable presence of cerebral ventriculomegaly, aqueductal stenosis, developmental delay and a variety of structural brain or cardiac defects. These data underscore the importance of SMARCC1 and the BAF chromatin remodelling complex for human brain morphogenesis and provide evidence for a 'neural stem cell' paradigm of congenital hydrocephalus pathogenesis. These results highlight utility of trio-based whole-exome sequencing for identifying pathogenic variants in sporadic congenital structural brain disorders and suggest whole-exome sequencing may be a valuable adjunct in clinical management of congenital hydrocephalus patients.

Sites of regulated phosphorylation that control K-Cl cotransporter activity.

Modulation of intracellular chloride concentration ([Cl(-)](i)) plays a fundamental role in cell volume regulation and neuronal response to GABA. Cl(-) exit via K-Cl cotransporters (KCCs) is a major determinant of [Cl(-)](I); however, mechanisms governing KCC activities are poorly understood. We identified two sites in KCC3 that are rapidly dephosphorylated in hypotonic conditions in cultured cells and human red blood cells in parallel with increased transport activity. Alanine substitutions at these sites result in constitutively active cotransport. These sites are highly phosphorylated in plasma membrane KCC3 in isotonic conditions, suggesting that dephosphorylation increases KCC3's intrinsic transport activity. Reduction of WNK1 expression via RNA interference reduces phosphorylation at these sites. Homologous sites are phosphorylated in all human KCCs. KCC2 is partially phosphorylated in neonatal mouse brain and dephosphorylated in parallel with KCC2 activation. These findings provide insight into regulation of [Cl(-)](i) and have implications for control of cell volume and neuronal function.

Molecular genetics of human developmental neurocranial anomalies: towards "precision surgery".

Recent trio-based whole-exome sequencing studies of congenital hydrocephalus and nonsyndromic craniosynostosis have identified multiple novel disease genes that have illuminated the pathogenesis of these disorders and shed new insight into the genetic regulation of human brain and skull development. Continued study of these and other historically understudied developmental anomalies has the potential to replace the current antiquated, anatomically based disease classification systems with a molecular nomenclature that may increase precision for genetic counseling, prognostication, and surgical treatment stratification-including when not to operate. Data will also inform future clinical trials, catalyze the development of targeted therapies, and generate infrastructure and publicly available data sets relevant for other related nonsurgical neurodevelopmental and neuropsychiatric diseases.

Enhanced Ca2+ signaling, mild primary aldosteronism, and hypertension in a familial hyperaldosteronism mouse model (Cacna1hM1560V/+ ).

Gain-of-function mutations in the CACNA1H gene (encoding the T-type calcium channel CaV3.2) cause autosomal-dominant familial hyperaldosteronism type IV (FH-IV) and early-onset hypertension in humans. We used CRISPR/Cas9 to generate Cacna1hM1560V/+ knockin mice as a model of the most common FH-IV mutation, along with corresponding knockout mice (Cacna1h-/- ). Adrenal morphology of both Cacna1hM1560V/+ and Cacna1h-/- mice was normal. Cacna1hM1560V/+ mice had elevated aldosterone:renin ratios (a screening parameter for primary aldosteronism). Their adrenal Cyp11b2 (aldosterone synthase) expression was increased and remained elevated on a high-salt diet (relative autonomy, characteristic of primary aldosteronism), but plasma aldosterone was only elevated in male animals. The systolic blood pressure of Cacna1hM1560V/+ mice was 8 mmHg higher than in wild-type littermates and remained elevated on a high-salt diet. Cacna1h-/- mice had elevated renal Ren1 (renin-1) expression but normal adrenal Cyp11b2 levels, suggesting that in the absence of CaV3.2, stimulation of the renin-angiotensin system activates alternative calcium entry pathways to maintain normal aldosterone production. On a cellular level, Cacna1hM1560V/+ adrenal slices showed increased baseline and peak intracellular calcium concentrations in the zona glomerulosa compared to controls, but the frequency of calcium spikes did not rise. We conclude that FH-IV, on a molecular level, is caused by elevated intracellular Ca2+ concentrations as a signal for aldosterone production in adrenal glomerulosa cells. We demonstrate that a germline Cacna1h gain-of-function mutation is sufficient to cause mild primary aldosteronism, whereas loss of CaV3.2 channel function can be compensated for in a chronic setting.

Integrated mutational landscape analysis of uterine leiomyosarcomas.

Uterine leiomyosarcomas (uLMS) are aggressive tumors arising from the smooth muscle layer of the uterus. We analyzed 83 uLMS sample genetics, including 56 from Yale and 27 from The Cancer Genome Atlas (TCGA). Among them, a total of 55 Yale samples including two patient-derived xenografts (PDXs) and 27 TCGA samples have whole-exome sequencing (WES) data; 10 Yale and 27 TCGA samples have RNA-sequencing (RNA-Seq) data; and 11 Yale and 10 TCGA samples have whole-genome sequencing (WGS) data. We found recurrent somatic mutations in TP53, MED12, and PTEN genes. Top somatic mutated genes included TP53, ATRX, PTEN, and MEN1 genes. Somatic copy number variation (CNV) analysis identified 8 copy-number gains, including 5p15.33 (TERT), 8q24.21 (C-MYC), and 17p11.2 (MYOCD, MAP2K4) amplifications and 29 copy-number losses. Fusions involving tumor suppressors or oncogenes were deetected, with most fusions disrupting RB1, TP53, and ATRX/DAXX, and one fusion (ACTG2-ALK) being potentially targetable. WGS results demonstrated that 76% (16 of 21) of the samples harbored chromoplexy and/or chromothripsis. Clinically actionable mutational signatures of homologous-recombination DNA-repair deficiency (HRD) and microsatellite instability (MSI) were identified in 25% (12 of 48) and 2% (1 of 48) of fresh frozen uLMS, respectively. Finally, we found olaparib (PARPi; P = 0.002), GS-626510 (C-MYC/BETi; P < 0.000001 and P = 0.0005), and copanlisib (PIK3CAi; P = 0.0001) monotherapy to significantly inhibit uLMS-PDXs harboring derangements in C-MYC and PTEN/PIK3CA/AKT genes (LEY11) and/or HRD signatures (LEY16) compared to vehicle-treated mice. These findings define the genetic landscape of uLMS and suggest that a subset of uLMS may benefit from existing PARP-, PIK3CA-, and C-MYC/BET-targeted drugs.

Genetic Variants in ARHGEF6 Cause Congenital Anomalies of the Kidneys and Urinary Tract in Humans, Mice, and Frogs.

BACKGROUND: About 40 disease genes have been described to date for isolated CAKUT, the most common cause of childhood CKD. However, these genes account for only 20% of cases. ARHGEF6, a guanine nucleotide exchange factor that is implicated in biologic processes such as cell migration and focal adhesion, acts downstream of integrin-linked kinase (ILK) and parvin proteins. A genetic variant of ILK that causes murine renal agenesis abrogates the interaction of ILK with a murine focal adhesion protein encoded by Parva , leading to CAKUT in mice with this variant. METHODS: To identify novel genes that, when mutated, result in CAKUT, we performed exome sequencing in an international cohort of 1265 families with CAKUT. We also assessed the effects in vitro of wild-type and mutant ARHGEF6 proteins, and the effects of Arhgef6 deficiency in mouse and frog models. RESULTS: We detected six different hemizygous variants in the gene ARHGEF6 (which is located on the X chromosome in humans) in eight individuals from six families with CAKUT. In kidney cells, overexpression of wild-type ARHGEF6 -but not proband-derived mutant ARHGEF6 -increased active levels of CDC42/RAC1, induced lamellipodia formation, and stimulated PARVA-dependent cell spreading. ARHGEF6-mutant proteins showed loss of interaction with PARVA. Three-dimensional Madin-Darby canine kidney cell cultures expressing ARHGEF6-mutant proteins exhibited reduced lumen formation and polarity defects. Arhgef6 deficiency in mouse and frog models recapitulated features of human CAKUT. CONCLUSIONS: Deleterious variants in ARHGEF6 may cause dysregulation of integrin-parvin-RAC1/CDC42 signaling, thereby leading to X-linked CAKUT.

De novo mutations in the BMP signaling pathway in lambdoid craniosynostosis.

Lambdoid craniosynostosis (CS) is a congenital anomaly resulting from premature fusion of the cranial suture between the parietal and occipital bones. Predominantly sporadic, it is the rarest form of CS and its genetic etiology is largely unexplored. Exome sequencing of 25 kindreds, including 18 parent-offspring trios with sporadic lambdoid CS, revealed a marked excess of damaging (predominantly missense) de novo mutations that account for ~ 40% of sporadic cases. These mutations clustered in the BMP signaling cascade (P = 1.6 × 10-7), including mutations in genes encoding BMP receptors (ACVRL1 and ACVR2A), transcription factors (SOX11, FOXO1) and a transcriptional co-repressor (IFRD1), none of which have been implicated in other forms of CS. These missense mutations are at residues critical for substrate or target sequence recognition and many are inferred to cause genetic gain-of-function. Additionally, mutations in transcription factor NFIX were implicated in syndromic craniosynostosis affecting diverse sutures. Single cell RNA sequencing analysis of the mouse lambdoid suture identified enrichment of mutations in osteoblast precursors (P = 1.6 × 10-6), implicating perturbations in the balance between proliferation and differentiation of osteoprogenitor cells in lambdoid CS. The results contribute to the growing knowledge of the genetics of CS, have implications for genetic counseling, and further elucidate the molecular etiology of premature suture fusion.

Mutations in ASPRV1 Cause Dominantly Inherited Ichthyosis.

The discovery of genetic causes of inherited skin disorders has been pivotal to the understanding of epidermal differentiation, function, and renewal. Here we show via exome sequencing that mutations in ASPRV1 (aspartic peptidase retroviral-like 1) cause a dominant Mendelian disorder featuring palmoplantar keratoderma and lamellar ichthyosis, a phenotype that has otherwise been exclusively recessive. ASPRV1 encodes a mammalian-specific and stratified epithelia-specific protease important in processing of filaggrin, a critical component of the uppermost epidermal layer. Three different heterozygous ASPRV1 missense mutations in four unrelated ichthyosis kindreds segregate with disease and disrupt protein residues within close proximity to each other and autocatalytic cleavage sites. Expression of mutant ASPRV1 proteins demonstrates that all three mutations alter ASPRV1 auto-cleavage and filaggrin processing, a function vital to epidermal barrier integrity.

DAAM2 Variants Cause Nephrotic Syndrome via Actin Dysregulation.

The discovery of >60 monogenic causes of nephrotic syndrome (NS) has revealed a central role for the actin regulators RhoA/Rac1/Cdc42 and their effectors, including the formin INF2. By whole-exome sequencing (WES), we here discovered bi-allelic variants in the formin DAAM2 in four unrelated families with steroid-resistant NS. We show that DAAM2 localizes to the cytoplasm in podocytes and in kidney sections. Further, the variants impair DAAM2-dependent actin remodeling processes: wild-type DAAM2 cDNA, but not cDNA representing missense variants found in individuals with NS, rescued reduced podocyte migration rate (PMR) and restored reduced filopodia formation in shRNA-induced DAAM2-knockdown podocytes. Filopodia restoration was also induced by the formin-activating molecule IMM-01. DAAM2 also co-localizes and co-immunoprecipitates with INF2, which is intriguing since variants in both formins cause NS. Using in vitro bulk and TIRF microscopy assays, we find that DAAM2 variants alter actin assembly activities of the formin. In a Xenopus daam2-CRISPR knockout model, we demonstrate actin dysregulation in vivo and glomerular maldevelopment that is rescued by WT-DAAM2 mRNA. We conclude that DAAM2 variants are a likely cause of monogenic human SRNS due to actin dysregulation in podocytes. Further, we provide evidence that DAAM2-associated SRNS may be amenable to treatment using actin regulating compounds.

Mechanisms of Congenital Heart Disease Caused by NAA15 Haploinsufficiency.

[Figure: see text].

Unexpected role of SIX1 variants in craniosynostosis: expanding the phenotype of SIX1-related disorders.

BACKGROUND: Pathogenic heterozygous SIX1 variants (predominantly missense) occur in branchio-otic syndrome (BOS), but an association with craniosynostosis has not been reported. METHODS: We investigated probands with craniosynostosis of unknown cause using whole exome/genome (n=628) or RNA (n=386) sequencing, and performed targeted resequencing of SIX1 in 615 additional patients. Expression of SIX1 protein in embryonic cranial sutures was examined in the Six1nLacZ/+ reporter mouse. RESULTS: From 1629 unrelated cases with craniosynostosis we identified seven different SIX1 variants (three missense, including two de novo mutations, and four nonsense, one of which was also present in an affected twin). Compared with population data, enrichment of SIX1 loss-of-function variants was highly significant (p=0.00003). All individuals with craniosynostosis had sagittal suture fusion; additionally four had bilambdoid synostosis. Associated BOS features were often attenuated; some carrier relatives appeared non-penetrant. SIX1 is expressed in a layer basal to the calvaria, likely corresponding to the dura mater, and in the mid-sagittal mesenchyme. CONCLUSION: Craniosynostosis is associated with heterozygous SIX1 variants, with possible enrichment of loss-of-function variants compared with classical BOS. We recommend screening of SIX1 in craniosynostosis, particularly when sagittal±lambdoid synostosis and/or any BOS phenotypes are present. These findings highlight the role of SIX1 in cranial suture homeostasis.

A genome-wide case-only test for the detection of digenic inheritance in human exomes.

Whole-exome sequencing (WES) has facilitated the discovery of genetic lesions underlying monogenic disorders. Incomplete penetrance and variable expressivity suggest a contribution of additional genetic lesions to clinical manifestations and outcome. Some monogenic disorders may therefore actually be digenic. However, only a few digenic disorders have been reported, all discovered by candidate gene approaches applied to at least one locus. We propose here a two-locus genome-wide test for detecting digenic inheritance in WES data. This approach uses the gene as the unit of analysis and tests all pairs of genes to detect pairwise gene × gene interactions underlying disease. It is a case-only method, which has several advantages over classic case-control tests, in particular by avoiding recruitment of controls. Our simulation studies based on real WES data identified two major sources of type I error inflation in this case-only test: linkage disequilibrium and population stratification. Both were corrected by specific procedures. Moreover, our case-only approach is more powerful than the corresponding case-control test for detecting digenic interactions in various population stratification scenarios. Finally, we confirmed the potential of our unbiased, genome-wide approach by successfully identifying a previously reported digenic lesion in patients with craniosynostosis. Our case-only test is a powerful and timely tool for detecting digenic inheritance in WES data from patients.