Genomic DNA was isolated from peripheral blood mononuclear cells obtained from consenting patients (n = 5) using QIAGEN FlexiGene chemistry on an AutoGen FLEX-STAR (Hollister, MA, USA). Briefly, red blood cells were lysed; white cells were pelleted by centrifugation, lysed, and treated with protease. DNA was precipitated with isopropanol, washed with 70% ethanol, and suspended in Tris-EDTA (TE) buffer. Paired-end libraries were prepared using the Bravo liquid handler (Agilent, Santa Clara, CA, USA) and whole exome capture was performed using an Agilent V4+UTR or Agilent V5+UTR kit. Concentration and size distribution of the libraries were determined on Qubit (Invitrogen, Carlsbad, CA, USA) and Agilent Bioanalyzer DNA 1000 chips. Exome libraries were loaded onto TruSeq Rapid run paired end flow cells and sequenced on an Illumina HiSeq 2500 sequencer using TruSeq Rapid SBS kit version 1. Base calling was performed using Illumina's RTA version 1.17.21.3. Filtering for variants predicted to be deleterious (missense, indel, stop codon change or changes at splice junctions) was accomplished using Ingenuity Variant Analysis version 5.2 settings for dominant variants. Variants with a call quality of at least 100 for cases and 50 for controls, outside the top 5% of exonically variable 100 base windows in healthy public genomes, and with an allele frequency of 0.5% or less in NHLBI exomes, ExAC database, or gnomAD were kept.