Novel findings with reassessment of exome data: implications for validation testing and interpretation of genomic data

Isolation of single cells from human hepatoblastoma tissues for whole-exome sequencing

By combining single-cell processing with whole-exome sequencing, we have developed single-cell whole-exome sequencing to investigate the mechanisms of hepatoblastoma development and to provide potential targets and therapeutic approaches for clinical treatment. In the following protocol, we outline the steps involved in single-cell sorting, whole-genome amplification, amplification uniformity estimation, and whole-exome library construction. In addition to the cells we use, this protocol is also suitable for other cell lines and cell types.

Challenges of genetic diagnosis of inborn errors of metabolism in a major tertiary care center in Lebanon

Background: Inborn errors of metabolism are rare genetic disorders; however, these are prevalent in countries with high consanguinity rates, like Lebanon. Patients are suspected, based on a combination of clinical and biochemical features; however, the final confirmation relies on genetic testing. Using next generation sequencing, as a new genetic investigational tool, carries several challenges for the physician, the geneticist, and the families. Methods: In this retrospective study, we analyzed the clinical, biochemical, and genetic profile of inborn errors of metabolism suspected patients, seen at a major tertiary care center in Lebanon, between 2015 and 2018. Genetic testing was performed using next generation sequencing. Genotype-phenotype correlation and diagnostic yield of each testing modality were studied. Results: Out of 211 patients genetically tested, 126 were suspected to have an inborn error of metabolism. The diagnostic yield of next generation sequencing reached 64.3%. Single gene testing was requested in 53%, whole exome sequencing in 36% and gene panels in 10%. Aminoacid disorders were mostly diagnosed followed by storage disorders, organic acidemias and mitochondrial diseases. Targeted testing was performed in 77% of aminoacid and organic acid disorders and half of suspected storage disorders. Single gene sequencing was positive in 75%, whereas whole exome sequencing diagnostic yield for complex cases, like mitochondrial disorders, reached 49%. Good clinical and biochemical correlation allowed the interpretation of variants of unknown significance and negative mutations as well as therapeutic management of most patients. Conclusion: Tailoring the choice of test modality, by next generation sequencing, to the category of suspected inborn errors of metabolism may lead to rapid diagnosis, shortcutting the cost of repeated testing. Whole exome sequencing as a first-tier investigation may be considered mainly for suspected mitochondrial diseases, whereas targeted sequencing can be offered upon suspicion of a specific enzyme deficiency. Timing and modality of gene test remain challenging, in view of the cost incurred by families.

MYH7 variants cause complex congenital heart disease

MYH7, encoding the myosin heavy chain sarcomeric β-myosin heavy chain, is a common cause of both hypertrophic and dilated cardiomyopathy. Additionally, families with left ventricular noncompaction cardiomyopathy (LVNC) and congenital heart disease (CHD), typically septal defects or Ebstein anomaly, have been identified to have heterozygous pathogenic variants in MHY7. One previous case of single ventricle CHD with heart failure due to a MYH7 variant has been identified. Herein, we present a single center's experience of complex CHD due to MYH7 variants. Three probands with a history of CHD, LVNC, and/or arrhythmias were identified to have MYH7 variants through multigene panel testing or exome sequencing. These three patients collectively had 12 affected family members, four with a history of Ebstein anomaly and seven with a history of LVNC. These findings suggest a wider phenotypic spectrum in MYH7-related CHD than previously understood. Further investigation into the possible role of MYH7 in CHD and mechanism of disease is necessary to fully delineate the phenotypic spectrum of MYH7-related cardiac disease. MYH7 should be considered for families with multiple individuals with complex CHD in the setting of a family history of LVNC or arrhythmias.

Re-analysis of genomic data: An overview of the mechanisms and complexities of clinical adoption

Re-analyzing genomic information from a patient suspected of having an underlying genetic condition can improve the diagnostic yield of sequencing tests, potentially providing significant benefits to the patient and to the health care system. Although a significant number of studies have shown the clinical potential of re-analysis, less work has been performed to characterize the mechanisms responsible for driving the increases in diagnostic yield. Complexities surrounding re-analysis have also emerged. The terminology itself represents a challenge because "re-analysis" can refer to a range of different concepts. Other challenges include the increased workload that re-analysis demands of curators, adequate reimbursement pathways for clinical and diagnostic services, and the development of systems to handle large volumes of data. Re-analysis also raises ethical implications for patients and families, most notably when re-classification of a variant alters diagnosis, treatment, and prognosis. This review highlights the possibilities and complexities associated with the re-analysis of existing clinical genomic data. We propose a terminology that builds on the foundation presented in a recent statement from the American College of Medical Genetics and Genomics and describes each re-analysis process. We identify mechanisms for increasing diagnostic yield and provide perspectives on the range of challenges that must be addressed by health care systems and individual patients.

Molecular Diagnostic Outcomes from 700 Cases: What Can We Learn from a Retrospective Analysis of Clinical Exome Sequencing?

Clinical exome sequencing (CES) aids in the diagnosis of rare genetic disorders. Herein, we report the molecular diagnostic yield and spectrum of genetic alterations contributing to disease in 700 pediatric cases analyzed at the Children's Hospital of Philadelphia. The overall diagnostic yield was 23%, with three cases having more than one molecular diagnosis and 2.6% having secondary/additional findings. A candidate gene finding was reported in another 8.4% of cases. The clinical indications with the highest diagnostic yield were neurodevelopmental disorders (including seizures), whereas immune- and oncology-related indications were negatively associated with molecular diagnosis. The rapid expansion of knowledge regarding the genome's role in human disease necessitates reanalysis of CES samples. To capture these new discoveries, a subset of cases (n = 240) underwent reanalysis, with an increase in diagnostic yield. We describe our experience reporting CES results in a pediatric setting, including reporting of secondary findings, reporting newly discovered genetic conditions, and revisiting negative test results. Finally, we highlight the challenges associated with implementing critical updates to the CES workflow. Although these updates are necessary, they demand an investment of time and resources from the laboratory. In summary, these data demonstrate the clinical utility of exome sequencing and reanalysis, while highlighting the critical considerations for continuous improvement of a CES test in a clinical laboratory.

Palmoplantar keratoderma with deafness phenotypic variability in a patient with an inherited GJB2 frameshift variant and novel missense variant

Background

Variants in the GJB2 gene encoding the gap junction protein connexin‐26 (Cx26) can cause autosomal recessive nonsyndromic hearing loss or a variety of phenotypically variable autosomal dominant disorders that effect skin and hearing, such as palmoplantar keratoderma (PPK) with deafness and keratitis–ichthyosis–deafness (KID) syndrome. Here, we report a patient with chronic mucocutaneous candidiasis, hyperkeratosis with resorption of the finger tips, profound bilateral sensorineural hearing loss, and normal hair and ocular examination. Exome analysis identified a novel missense variant in GJB2 (NM_004004.5:c.101T>A, p.Met34Lys) that was inherited from a mosaic unaffected parent in the setting of a well‐reported GJB2 loss of function variant (NM_004004.5:c.35delG, p.Gly12Valfs*2) on the other allele.

Method

Rat epidermal keratinocytes were transfected with cDNA encoding wildtype Cx26 and/or the Met34Lys mutant of Cx26. Fixed cells were immunolabeled in order to assess the subcellular location of the Cx26 mutant and cell images were captured.

Results

Expression in rat epidermal keratinocytes revealed that the Met34Lys mutant was retained in the endoplasmic reticulum, unlike wildtype Cx26, and failed to reach the plasma membrane to form gap junctions. Additionally, the Met34Lys mutant acted dominantly to wildtype Cx26, restricting its delivery to the cell surface.

Conclusion

Overall, we show the p.Met34Lys variant is a novel dominant acting variant causing PPK with deafness. The presence of a loss a function variant on the other allele creates a more severe clinical phenotype, with some features reminiscent of KID syndrome.

A flexible computational pipeline for research analyses of unsolved clinical exome cases

Exome sequencing has enabled molecular diagnoses for rare disease patients but often with initial diagnostic rates of ~25-30%. Here we develop a robust computational pipeline to rank variants for reassessment of unsolved rare disease patients. A comprehensive web-based patient report is generated in which all deleterious variants can be filtered by gene, variant characteristics, OMIM disease and Phenolyzer scores, and all are annotated with an ACMG classification and links to ClinVar. The pipeline ranked 21/34 previously diagnosed variants as top, with 26 in total ranked ≤7th, 3 ranked ≥13th; 5 failed the pipeline filters. Pathogenic/likely pathogenic variants by ACMG criteria were identified for 22/145 unsolved cases, and a previously undefined candidate disease variant for 27/145. This open access pipeline supports the partnership between clinical and research laboratories to improve the diagnosis of unsolved exomes. It provides a flexible framework for iterative developments to further improve diagnosis.

Variants in NAA15 cause pediatric hypertrophic cardiomyopathy

The NatA N-acetyltransferase complex is important for cotranslational protein modification and regulation of multiple cellular processes. The NatA complex includes the core components of NAA10, the catalytic subunit, and NAA15, the auxiliary component. Both NAA10 and NAA15 have been associated with neurodevelopmental disorders with overlapping clinical features, including variable intellectual disability, dysmorphic facial features, and, less commonly, congenital anomalies such as cleft lip or palate. Cardiac arrhythmias, including long QT syndrome, ventricular tachycardia, and ventricular fibrillation were among the first reported cardiac manifestations in patients with NAA10-related syndrome. Recently, three individuals with NAA10-related syndrome have been reported to also have hypertrophic cardiomyopathy (HCM). The general and cardiac phenotypes of NAA15-related syndrome are not as well described as NAA10-related syndrome. Congenital heart disease, including ventricular septal defects, and arrhythmias, such as ectopic atrial tachycardia, have been reported in a small proportion of patients with NAA15-related syndrome. Given the relationship between NAA10 and NAA15, we propose that HCM is also likely to occur in NAA15-related disorder. We present two patients with pediatric HCM found to have NAA15-related disorder via exome sequencing, providing the first evidence that variants in NAA15 can cause HCM.

JARID2 haploinsufficiency is associated with a clinically distinct neurodevelopmental syndrome

PURPOSE

JARID2, located on chromosome 6p22.3, is a regulator of histone methyltransferase complexes that is expressed in human neurons. So far, 13 individuals sharing clinical features including intellectual disability (ID) were reported with de novo heterozygous deletions in 6p22-p24 encompassing the full length JARID2 gene (OMIM 601594). However, all published individuals to date have a deletion of at least one other adjoining gene, making it difficult to determine if JARID2 is the critical gene responsible for the shared features. We aim to confirm JARID2 as a human disease gene and further elucidate the associated clinical phenotype.

METHODS

Chromosome microarray analysis, exome sequencing, and an online matching platform (GeneMatcher) were used to identify individuals with single-nucleotide variants or deletions involving JARID2.

RESULTS

We report 16 individuals in 15 families with a deletion or single-nucleotide variant in JARID2. Several of these variants are likely to result in haploinsufficiency due to nonsense-mediated messenger RNA (mRNA) decay. All individuals have developmental delay and/or ID and share some overlapping clinical characteristics such as facial features with those who have larger deletions involving JARID2.

CONCLUSION

We report that JARID2 haploinsufficiency leads to a clinically distinct neurodevelopmental syndrome, thus establishing gene-disease validity for the purpose of diagnostic reporting.

Interplay between probe design and test performance: overlap between genomic regions of interest, capture regions and high quality reference calls influence performance of WES-based assays

Whole exome sequencing (WES)-based assays undergo rigorous validation before being implemented in diagnostic laboratories. This validation process generates experimental evidence that allows laboratories to predict the performance of the intended assay. The NA12878 Genome in a Bottle (GIAB) HapMap reference sample is commonly used for validation in diagnostic laboratories. We investigated what data points should be taken into consideration when validating WES-based assays using the GIAB reference in a diagnostic setting. We delineate specific factors that require special consideration and identify OMIM genes associated with diseases that may 'bypass' validation. Four replicates of the NA12878 sample were sequenced at the CHEO Genetics Diagnostic Laboratory on a NextSeq 500; the data were analyzed using the bcbio_nexgen v1.1.2 pipeline. The hap.py validation engine, Real Time Genomics vcfeval tool, and high confidence (HC) variant calls in HC regions available for the GIAB sample were used to validate the obtained variant calls. The same validation process was then used to evaluate variant calls obtained for the same sample by two other clinical diagnostic laboratories. We showed that variant calls in NA12878 can be confidently measured only in the regions that intersect between the GIAB HC regions and the target regions of exome capture. Of the 4139 (as of October 2019) OMIM genes associated with a phenotype and having a known molecular basis of disease, 84 were fully outside of the GIAB HC regions and many of the remaining OMIM genes were only partially covered by the HC regions. A significant proportion of variants identified in the NA12878 sample outside of the HC regions have unknown (UNK) status due to the absence of HC reference alleles. Verification of such calls is possible either by an alternative truth set or by orthogonal testing. Similarly, many variants outside of exome capture regions, if not accounted for, will be deemed false negatives due to insufficient probe coverage. Our results demonstrate the importance of the intersection between genomic regions of interest, capture regions, and the high confidence regions. If not considered, false and ambiguous variant calls could have a negative impact on diagnostic accuracy of the intended WES-based diagnostic assay and increase the need for confirmatory testing. To enable laboratories to identify 'problematic' regions and optimize validation efforts, we have made our VCF and BED files available in UCSC Genome Browser: NA12878 WES Benchmark. Relevant genes and genome annotations are evolving, we implemented a general purpose algorithm to cross-reference OMIM genes with the genomic regions of interest that can be applied to capture genes/regions outside HC regions (see repository of data material section).

Loss-of-function of Endothelin receptor type A results in Oro-Oto-Cardiac syndrome

Craniofacial morphogenesis is regulated in part by signaling from the Endothelin receptor type A (EDNRA). Pathogenic variants in EDNRA signaling pathway components EDNRA, GNAI3, PCLB4, and EDN1 cause Mandibulofacial Dysostosis with Alopecia (MFDA), Auriculocondylar syndrome (ARCND) 1, 2, and 3, respectively. However, cardiovascular development is normal in MFDA and ARCND individuals, unlike Ednra knockout mice. One explanation may be that partial EDNRA signaling remains in MFDA and ARCND, as mice with reduced, but not absent, EDNRA signaling also lack a cardiovascular phenotype. Here we report an individual with craniofacial and cardiovascular malformations mimicking the Ednra -/- mouse phenotype, including a distinctive micrognathia with microstomia and a hypoplastic aortic arch. Exome sequencing found a novel homozygous missense variant in EDNRA (c.1142A>C; p.Q381P). Bioluminescence resonance energy transfer assays revealed that this amino acid substitution in helix 8 of EDNRA prevents recruitment of G proteins to the receptor, abrogating subsequent receptor activation by its ligand, Endothelin-1. This homozygous variant is thus the first reported loss-of-function EDNRA allele, resulting in a syndrome we have named Oro-Oto-Cardiac Syndrome. Further, our results illustrate that EDNRA signaling is required for both normal human craniofacial and cardiovascular development, and that limited EDNRA signaling is likely retained in ARCND and MFDA individuals. This work illustrates a straightforward approach to identifying the functional consequence of novel genetic variants in signaling molecules associated with malformation syndromes.

NKX2-6 related congenital heart disease: Biallelic homeodomain-disrupting variants and truncus arteriosus

Congenital heart defects (CHD) are the most common birth defect and are both clinically and genetically heterogeneous. Truncus arteriosus (TA), characterized by a single arterial vessel arising from both ventricles giving rise to the coronary, pulmonary and systemic arteries, is rare and only responsible for 1% of all CHD. Two consanguineous families with TA were previously identified to have homozygous nonsense variants within the gene NKX2-6. NKX2-6 is a known downstream target of TBX1, an important transcriptional regulator implicated in the cardiac phenotype of 22q11.2 microdeletion syndrome. Herein, we report two siblings with TA presumably caused by compound heterozygous NKX2-6 variants without a history of consanguinity. Two in-house cohorts with conotruncal defects (CTD) were sequenced for variants in NKX2-6 and no additional cases of biallelic NKX2-6 variants were identified. The similar phenotype of these cases, and the clustering of variants that likely result in a truncated protein that disrupts the homeobox domain, suggest that biallelic loss of function for NKX2-6 is a rare genetic etiology for TA in particular, and possibly other types of CHD.

Limitations of exome sequencing in detecting rare and undiagnosed diseases

While exome sequencing (ES) is commonly the final diagnostic step in clinical genetics, it may miss diagnoses. To clarify the limitations of ES, we investigated the diagnostic yield of genetic tests beyond ES in our Undiagnosed Diseases Network (UDN) participants. We reviewed the yield of additional genetic testing including genome sequencing (GS), copy number variant (CNV), noncoding variant (NCV), repeat expansion (RE), or methylation testing in UDN cases with nondiagnostic ES results. Overall, 36/54 (67%) of total diagnoses were based on clinical findings and coding variants found by ES and 3/54 (6%) were based on clinical findings only. The remaining 15/54 (28%) required testing beyond ES. Of these, 7/15 (47%) had NCV, 6/15 (40%) CNV, and 2/15 (13%) had a RE or a DNA methylation disorder. Thus 18/54 (33%) of diagnoses were not solved exclusively by ES. Several methods were needed to detect and/or confirm the functional effects of the variants missed by ES, and in some cases by GS. These results indicate that tests to detect elusive variants should be considered after nondiagnostic preliminary steps. Further studies are needed to determine the cost-effectiveness of tests beyond ES that provide diagnoses and insights to possible treatment.

Imprinted genes in clinical exome sequencing: Review of 538 cases and exploration of mouse-human conservation in the identification of novel human disease loci

Human imprinting disorders cause a range of dysmorphic and neurocognitive phenotypes, and they may elude traditional molecular diagnosis such exome sequencing. The discovery of novel disorders related to imprinted genes has lagged behind traditional Mendelian disorders because current diagnostic technology, especially unbiased testing, has limited utility in their discovery. To identify novel imprinting disorders, we reviewed data for every human gene hypothesized to be imprinted, identified each mouse ortholog, determined its imprinting status in the mouse, and analyzed its function in humans and mice. We identified 17 human genes that are imprinted in both humans and mice, and have functional data in mice or humans to suggest that dysregulated expression would lead to an abnormal phenotype in humans. These 17 genes, along with known imprinted genes, were preferentially flagged 538 clinical exome sequencing tests. The identified genes were: DIRAS3 [1p31.3], TP73 [1p36.32], SLC22A3 [6q25.3], GRB10 [7p12.1], DDC [7p12.2], MAGI2 [7q21.11], PEG10 [7q21.3], PPP1R9A [7q21.3], CALCR [7q21.3], DLGAP2 [8p23.3], GLIS3 [9p24.2], INPP5F [10q26.11], ANO1 [11q13.3], SLC38A4 [12q13.11], GATM [15q21.1], PEG3 [19q13.43], and NLRP2 [19q13.42]. In the 538 clinical cases, eight cases (1.7%) reported variants in a causative known imprinted gene. There were 367/758 variants (48.4%) in imprinted genes that were not known to cause disease, but none of those variants met the criteria for clinical reporting. Imprinted disorders play a significant role in human disease, and additional human imprinted disorders remain to be discovered. Therefore, evolutionary conservation is a potential tool to identify novel genes involved in human imprinting disorders and to identify them in clinical testing.

A highly sensitive and specific workflow for detecting rare copy-number variants from exome sequencing data

Background

Exome sequencing (ES) is a first-tier diagnostic test for many suspected Mendelian disorders. While it is routine to detect small sequence variants, it is not a standard practice in clinical settings to detect germline copy-number variants (CNVs) from ES data due to several reasons relating to performance. In this work, we comprehensively characterized one of the most sensitive ES-based CNV tools, ExomeDepth, against SNP array, a standard of care test in clinical settings to detect genome-wide CNVs.

Methods

We propose a modified ExomeDepth workflow by excluding exons with low mappability prior to variant calling to drastically reduce the false positives originating from the repetitive regions of the genome, and an iterative variant calling framework to assess the reproducibility. We used a cohort of 307 individuals with clinical ES data and clinical SNP array to estimate the sensitivity and false discovery rate of the CNV detection using exome sequencing. Further, we performed targeted testing of the STRC gene in 1972 individuals. To reduce the number of variants for downstream analysis, we performed a large-scale iterative variant calling process with random control cohorts to assess the reproducibility of the CNVs.

Results

The modified workflow presented in this paper reduced the number of total variants identified by one third while retaining a higher sensitivity of 97% and resulted in an improved false discovery rate of 11.4% compared to the default ExomeDepth pipeline. The exclusion of exons with low mappability removes 4.5% of the exons, including a subset of exons (0.6%) in disease-associated genes which are intractable by short-read next-generation sequencing (NGS). Results from the reproducibility analysis showed that the clinically reported variants were reproducible 100% of the time and that the modified workflow can be used to rank variants from high to low confidence. Targeted testing of 30 CNVs identified in STRC, a challenging gene to ascertain by NGS, showed a 100% validation rate.

Conclusions

In summary, we introduced a modification to the default ExomeDepth workflow to reduce the false positives originating from the repetitive regions of the genome, created a large-scale iterative variant calling framework for reproducibility, and provided recommendations for implementation in clinical settings.

Increased diagnostic yield by reanalysis of data from a hearing loss gene panel

Background

Congenital hearing loss affects approximately 1–2 infants out of every 1000, with 50% of the cases resulting from genetic factors. Targeted gene panels have been widely used for genetic diagnosis of hearing loss. This study aims to reveal new diagnoses via reanalyzing historical data of a multigene panel, and exam the reasons for new diagnoses.

Methods

A total of 210 samples were enlisted, including clinical reports and sequencing data of patients with congenital/prelingual hearing loss who were referred to clinical genetic testing from October 2014 to June 2017. All variants listed on the original clinical reports were reinterpreted according to the standards and guidelines recommended by the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP). Expanded analysis of raw data were performed in undiagnosed cases.

Results

Re-analysis resulted in nine new diagnoses, improving the overall diagnostic rate from 39 to 43%. New diagnoses were attributed to newly published clinical evidence in the literature, adoption of new interpretation guidelines and expanded analysis range.

Conclusion

This work demonstrates benefits of reanalysis of targeted gene panel data, indicating that periodical reanalysis should be performed in clinical practice.

Electronic supplementary material

The online version of this article (10.1186/s12920-019-0531-6) contains supplementary material, which is available to authorized users.

Rapid and accurate interpretation of clinical exomes using Phenoxome: a computational phenotype-driven approach

Clinical exome sequencing (CES) has become the preferred diagnostic platform for complex pediatric disorders with suspected monogenic etiologies. Despite rapid advancements, the major challenge still resides in identifying the casual variants among the thousands of variants detected during CES testing, and thus establishing a molecular diagnosis. To improve the clinical exome diagnostic efficiency, we developed Phenoxome, a robust phenotype-driven model that adopts a network-based approach to facilitate automated variant prioritization. Phenoxome dissects the phenotypic manifestation of a patient in concert with their genomic profile to filter and then prioritize variants that are likely to affect the function of the gene (potentially pathogenic variants). To validate our method, we have compiled a clinical cohort of 105 positive patient samples that represent a wide range of genetic heterogeneity. Phenoxome identifies the causative variants within the top 5, 10, or 25 candidates in more than 50%, 71%, or 88% of these exomes, respectively. Furthermore, we show that our method is optimized for clinical testing by outperforming the current state-of-art method. We have demonstrated the performance of Phenoxome using a clinical cohort and showed that it enables rapid and accurate interpretation of clinical exomes. Phenoxome is available at https://phenoxome.chop.edu/ .

A comprehensive iterative approach is highly effective in diagnosing individuals who are exome negative

PURPOSE

Sixty to seventy-five percent of individuals with rare and undiagnosed phenotypes remain undiagnosed after exome sequencing (ES). With standard ES reanalysis resolving 10-15% of the ES negatives, further approaches are necessary to maximize diagnoses in these individuals.

METHODS

In 38 ES negative patients an individualized genomic-phenotypic approach was employed utilizing (1) phenotyping; (2) reanalyses of FASTQ files, with innovative bioinformatics; (3) targeted molecular testing; (4) genome sequencing (GS); and (5) conferring of clinical diagnoses when pathognomonic clinical findings occurred.

RESULTS

Certain and highly likely diagnoses were made in 18/38 (47%) individuals, including identifying two new developmental disorders. The majority of diagnoses (>70%) were due to our bioinformatics, phenotyping, and targeted testing identifying variants that were undetected or not prioritized on prior ES. GS diagnosed 3/18 individuals with structural variants not amenable to ES. Additionally, tentative diagnoses were made in 3 (8%), and in 5 individuals (13%) candidate genes were identified. Overall, diagnoses/potential leads were identified in 26/38 (68%).

CONCLUSIONS

Our comprehensive approach to ES negatives maximizes the ES and clinical data for both diagnoses and candidate gene identification, without GS in the majority. This iterative approach is cost-effective and is pertinent to the current conundrum of ES negatives.