Background: Congenital heart disease (CHD) affects ~ 1% of live births and is the most common birth defect. Although the genetic contribution to the CHD has been long suspected, it has only been well established recently. De novo variants are estimated to contribute to approximately 8% of sporadic CHD. Methods: CHD is genetically heterogeneous, making pathway enrichment analysis an effective approach to explore and statistically validate CHD-associated genes. In this study, we performed novel gene and pathway enrichment analyses of high-impact de novo variants in the recently published whole-exome sequencing (WES) data generated from a cohort of CHD 2645 parent-offspring trios to identify new CHD-causing candidate genes and mutations. We performed rigorous variant- and gene-level filtrations to identify potentially damaging variants, followed by enrichment analyses and gene prioritization. Results: Our analyses revealed 23 novel genes that are likely to cause CHD, including HSP90AA1, ROCK2, IQGAP1, and CHD4, and sharing biological functions, pathways, molecular interactions, and properties with known CHD-causing genes. Conclusions: Ultimately, these findings suggest novel genes that are likely to be contributing to CHD pathogenesis.

Computational analyses of human patient exomes aim to filter out as many nonpathogenic genetic variants (NPVs) as possible, without removing the true disease-causing mutations. This involves comparing the patient’s exome with public databases to remove reported variants inconsistent with disease prevalence, mode of inheritance, or clinical penetrance. However, variants frequent in a given exome cohort, but absent or rare in public databases, have also been reported and treated as NPVs, without rigorous exploration. We report the generation of a blacklist of variants frequent within an in-house cohort of 3,104 exomes. This blacklist did not remove known pathogenic mutations from the exomes of 129 patients and decreased the number of NPVs remaining in the 3,104 individual exomes by a median of 62%. We validated this approach by testing three other independent cohorts of 400, 902, and 3,869 exomes. The blacklist generated from any given cohort removed a substantial proportion of NPVs (11–65%). We analyzed the blacklisted variants computationally and experimentally. Most of the blacklisted variants corresponded to false signals generated by incomplete reference genome assembly, location in low-complexity regions, bioinformatic misprocessing, or limitations inherent to cohort-specific private alleles (e.g., due to sequencing kits, and genetic ancestries). Finally, we provide our precalculated blacklists, together with ReFiNE, a program for generating customized blacklists from any medium-sized or large in-house cohort of exome (or other next-generation sequencing) data via a user-friendly public web server. This work demonstrates the power of extracting variant blacklists from private databases as a specific in-house but broadly applicable tool for optimizing exome analysis.

Gain-of-function (GOF) variants give rise to increased/novel protein functions whereas loss-of-function (LOF) variants lead to diminished protein function. Experimental approaches for identifying GOF and LOF are generally slow and costly, whilst available computational methods have not been optimized to discriminate between GOF and LOF variants. We have developed LoGoFunc, a machine learning method for predicting pathogenic GOF, pathogenic LOF, and neutral genetic variants, trained on a broad range of gene-, protein-, and variant-level features describing diverse biological characteristics. LoGoFunc outperforms other tools trained solely to predict pathogenicity for identifying pathogenic GOF and LOF variants and is available athttps://itanlab.shinyapps.io/goflof/.

Attention deficit/hyperactivity disorder (ADHD) is a common and heritable phenotype frequently accompanied by insomnia, anxiety, and depression. Here, using a reverse phenotyping approach, we report heterozygous coding variations in the core circadian clock gene cryptochrome 1 in 15 unrelated multigenerational families with combined ADHD and insomnia. The variants led to functional alterations in the circadian molecular rhythms, providing a mechanistic link to the behavioral symptoms. One variant, CRY1Δ11 c.1657+3A>C, is present in approximately 1% of Europeans, therefore standing out as a diagnostic and therapeutic marker. We showed by exome sequencing in an independent cohort of patients with combined ADHD and insomnia that 8 of 62 patients and 0 of 369 controls carried CRY1Δ11. Also, we identified a variant, CRY1Δ6 c.825+1G>A, that shows reduced affinity for BMAL1/CLOCK and causes an arrhythmic phenotype. Genotype-phenotype correlation analysis revealed that this variant segregated with ADHD and delayed sleep phase disorder (DSPD) in the affected family. Finally, we found in a phenome-wide association study involving 9438 unrelated adult Europeans that CRY1Δ11 was associated with major depressive disorder, insomnia, and anxiety. These results defined a distinctive group of circadian psychiatric phenotypes that we propose to designate as “circiatric” disorders.