Beyond the exome: What's next in diagnostic testing for Mendelian conditions.

Despite advances in clinical genetic testing, including the introduction of exome sequencing (ES), more than 50% of individuals with a suspected Mendelian condition lack a precise molecular diagnosis. Clinical evaluation is increasingly undertaken by specialists outside of clinical genetics, often occurring in a tiered fashion and typically ending after ES. The current diagnostic rate reflects multiple factors, including technical limitations, incomplete understanding of variant pathogenicity, missing genotype-phenotype associations, complex gene-environment interactions, and reporting differences between clinical labs. Maintaining a clear understanding of the rapidly evolving landscape of diagnostic tests beyond ES, and their limitations, presents a challenge for non-genetics professionals. Newer tests, such as short-read genome or RNA sequencing, can be challenging to order, and emerging technologies, such as optical genome mapping and long-read DNA sequencing, are not available clinically. Furthermore, there is no clear guidance on the next best steps after inconclusive evaluation. Here, we review why a clinical genetic evaluation may be negative, discuss questions to be asked in this setting, and provide a framework for further investigation, including the advantages and disadvantages of new approaches that are nascent in the clinical sphere. We present a guide for the next best steps after inconclusive molecular testing based upon phenotype and prior evaluation, including when to consider referral to research consortia focused on elucidating the underlying cause of rare unsolved genetic disorders.

Beyond the exome: what's next in diagnostic testing for Mendelian conditions.

Despite advances in clinical genetic testing, including the introduction of exome sequencing (ES), more than 50% of individuals with a suspected Mendelian condition lack a precise molecular diagnosis. Clinical evaluation is increasingly undertaken by specialists outside of clinical genetics, often occurring in a tiered fashion and typically ending after ES. The current diagnostic rate reflects multiple factors, including technical limitations, incomplete understanding of variant pathogenicity, missing genotype-phenotype associations, complex gene-environment interactions, and reporting differences between clinical labs. Maintaining a clear understanding of the rapidly evolving landscape of diagnostic tests beyond ES, and their limitations, presents a challenge for non-genetics professionals. Newer tests, such as short-read genome or RNA sequencing, can be challenging to order and emerging technologies, such as optical genome mapping and long-read DNA or RNA sequencing, are not available clinically. Furthermore, there is no clear guidance on the next best steps after inconclusive evaluation. Here, we review why a clinical genetic evaluation may be negative, discuss questions to be asked in this setting, and provide a framework for further investigation, including the advantages and disadvantages of new approaches that are nascent in the clinical sphere. We present a guide for the next best steps after inconclusive molecular testing based upon phenotype and prior evaluation, including when to consider referral to a consortium such as GREGoR, which is focused on elucidating the underlying cause of rare unsolved genetic disorders.

Transcriptional and functional consequences of alterations to MEF2C and its topological organization in neuronal models.

Point mutations and structural variants that directly disrupt the coding sequence of MEF2C have been associated with a spectrum of neurodevelopmental disorders (NDDs). However, the impact of MEF2C haploinsufficiency on neurodevelopmental pathways and synaptic processes is not well understood, nor are the complex mechanisms that govern its regulation. To explore the functional changes associated with structural variants that alter MEF2C expression and/or regulation, we generated an allelic series of 204 isogenic human induced pluripotent stem cell (hiPSC)-derived neural stem cells and glutamatergic induced neurons. These neuronal models harbored CRISPR-engineered mutations that involved direct deletion of MEF2C or deletion of the boundary points for topologically associating domains (TADs) and chromatin loops encompassing MEF2C. Systematic profiling of mutation-specific alterations, contrasted to unedited controls that were exposed to the same guide RNAs for each edit, revealed that deletion of MEF2C caused differential expression of genes associated with neurodevelopmental pathways and synaptic function. We also discovered significant reduction in synaptic activity measured by multielectrode arrays (MEAs) in neuronal cells. By contrast, we observed robust buffering against MEF2C regulatory disruption following deletion of a distal 5q14.3 TAD and loop boundary, whereas homozygous loss of a proximal loop boundary resulted in down-regulation of MEF2C expression and reduced electrophysiological activity on MEA that was comparable to direct gene disruption. Collectively, these studies highlight the considerable functional impact of MEF2C deletion in neuronal cells and systematically characterize the complex interactions that challenge a priori predictions of regulatory consequences from structural variants that disrupt three-dimensional genome organization.

A cross-disorder dosage sensitivity map of the human genome.

Rare copy-number variants (rCNVs) include deletions and duplications that occur infrequently in the global human population and can confer substantial risk for disease. In this study, we aimed to quantify the properties of haploinsufficiency (i.e., deletion intolerance) and triplosensitivity (i.e., duplication intolerance) throughout the human genome. We harmonized and meta-analyzed rCNVs from nearly one million individuals to construct a genome-wide catalog of dosage sensitivity across 54 disorders, which defined 163 dosage sensitive segments associated with at least one disorder. These segments were typically gene dense and often harbored dominant dosage sensitive driver genes, which we were able to prioritize using statistical fine-mapping. Finally, we designed an ensemble machine-learning model to predict probabilities of dosage sensitivity (pHaplo & pTriplo) for all autosomal genes, which identified 2,987 haploinsufficient and 1,559 triplosensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will provide broad utility for human disease research and clinical genetics.

Further delineation of the clinical spectrum of KAT6B disorders and allelic series of pathogenic variants.

PURPOSE: Genitopatellar syndrome and Say-Barber-Biesecker-Young-Simpson syndrome are caused by variants in the KAT6B gene and are part of a broad clinical spectrum called KAT6B disorders, whose variable expressivity is increasingly being recognized. METHODS: We herein present the phenotypes of 32 previously unreported individuals with a molecularly confirmed diagnosis of a KAT6B disorder, report 24 new pathogenic KAT6B variants, and review phenotypic information available on all published individuals with this condition. We also suggest a classification of clinical subtypes within the KAT6B disorder spectrum. RESULTS: We demonstrate that cerebral anomalies, optic nerve hypoplasia, neurobehavioral difficulties, and distal limb anomalies other than long thumbs and great toes, such as polydactyly, are more frequently observed than initially reported. Intestinal malrotation and its serious consequences can be present in affected individuals. Additionally, we identified four children with Pierre Robin sequence, four individuals who had increased nuchal translucency/cystic hygroma prenatally, and two fetuses with severe renal anomalies leading to renal failure. We also report an individual in which a pathogenic variant was inherited from a mildly affected parent. CONCLUSION: Our work provides a comprehensive review and expansion of the genotypic and phenotypic spectrum of KAT6B disorders that will assist clinicians in the assessment, counseling, and management of affected individuals.

Biallelic mutation of FBXL7 suggests a novel form of Hennekam syndrome.

Hennekam lymphangiectasia-lymphedema syndrome is an autosomal recessive disorder characterized by congenital lymphedema, intestinal lymphangiectasia, facial dysmorphism, and variable intellectual disability. Known disease genes include CCBE1, FAT4, and ADAMTS3. In a patient with clinically diagnosed Hennekam syndrome but without mutations or copy-number changes in the three known disease genes, we identified a homozygous single-exon deletion affecting FBXL7. Specifically, exon 3, which encodes the F-box domain and several leucine-rich repeats of FBXL7, is eliminated. Our analyses of databases representing >100,000 control individuals failed to identify biallelic loss-of-function variants in FBXL7. Published studies in Drosophila indicate Fbxl7 interacts with Fat, of which human FAT4 is an ortholog, and mutation of either gene yields similar morphological consequences. These data suggest that FBXL7 may be the fourth gene for Hennekam syndrome, acting via a shared pathway with FAT4.

The Genetics of Pneumothorax.

A genetic influence on spontaneous pneumothoraces-those occurring without a traumatic or iatrogenic cause-is supported by several lines of evidence: 1) pneumothorax can cluster in families (i.e., familial spontaneous pneumothorax), 2) mutations in the FLCN gene have been found in both familial and sporadic cases, and 3) pneumothorax is a known complication of several genetic syndromes. Herein, we review known genetic contributions to both sporadic and familial pneumothorax. We summarize the pneumothorax-associated genetic syndromes, including Birt-Hogg-Dubé syndrome, Marfan syndrome, vascular (type IV) Ehlers-Danlos syndrome, alpha-1 antitrypsin deficiency, tuberous sclerosis complex/lymphangioleiomyomatosis, Loeys-Dietz syndrome, cystic fibrosis, homocystinuria, and cutis laxa, among others. At times, pneumothorax is their herald manifestation. These syndromes have serious potential extrapulmonary complications (e.g., malignant renal tumors in Birt-Hogg-Dubé syndrome), and surveillance and/or treatment is available for most disorders; thus, establishing a diagnosis is critical. To facilitate this, we provide an algorithm to guide the clinician in discerning which cases of spontaneous pneumothorax may have a genetic or familial contribution, which cases warrant genetic testing, and which cases should prompt an evaluation by a geneticist.

Genetic architecture of laterality defects revealed by whole exome sequencing.

Aberrant left-right patterning in the developing human embryo can lead to a broad spectrum of congenital malformations. The causes of most laterality defects are not known, with variants in established genes accounting for <20% of cases. We sought to characterize the genetic spectrum of these conditions by performing whole-exome sequencing of 323 unrelated laterality cases. We investigated the role of rare, predicted-damaging variation in 1726 putative laterality candidate genes derived from model organisms, pathway analyses, and human phenotypes. We also evaluated the contribution of homo/hemizygous exon deletions and gene-based burden of rare variation. A total of 28 candidate variants (26 rare predicted-damaging variants and 2 hemizygous deletions) were identified, including variants in genes known to cause heterotaxy and primary ciliary dyskinesia (ACVR2B, NODAL, ZIC3, DNAI1, DNAH5, HYDIN, MMP21), and genes without a human phenotype association, but with prior evidence for a role in embryonic laterality or cardiac development. Sanger validation of the latter variants in probands and their parents revealed no de novo variants, but apparent transmitted heterozygous (ROCK2, ISL1, SMAD2), and hemizygous (RAI2, RIPPLY1) variant patterns. Collectively, these variants account for 7.1% of our study subjects. We also observe evidence for an excess burden of rare, predicted loss-of-function variation in PXDNL and BMS1- two genes relevant to the broader laterality phenotype. These findings highlight potential new genes in the development of laterality defects, and suggest extensive locus heterogeneity and complex genetic models in this class of birth defects.

Familial pneumothorax: towards precision medicine.

One in 10 patients suffering from primary spontaneous pneumothoraces has a family history of the disorder. Such familial pneumothoraces can occur in isolation, but can also be the presentation of serious genetic disorders with life-threatening vascular or cancerous complications. As the pneumothorax frequently precedes the more dangerous complications by many years, it provides an opportunity to intervene in a focused manner, permitting the practice of precision medicine. In this review, we will discuss the clinical manifestations and underlying biology of the genetic causes of familial pneumothorax.

Homozygous and hemizygous CNV detection from exome sequencing data in a Mendelian disease cohort.

We developed an algorithm, HMZDelFinder, that uses whole exome sequencing (WES) data to identify rare and intragenic homozygous and hemizygous (HMZ) deletions that may represent complete loss-of-function of the indicated gene. HMZDelFinder was applied to 4866 samples in the Baylor-Hopkins Center for Mendelian Genomics (BHCMG) cohort and detected 773 HMZ deletion calls (567 homozygous or 206 hemizygous) with an estimated sensitivity of 86.5% (82% for single-exonic and 88% for multi-exonic calls) and precision of 78% (53% single-exonic and 96% for multi-exonic calls). Out of 773 HMZDelFinder-detected deletion calls, 82 were subjected to array comparative genomic hybridization (aCGH) and/or breakpoint PCR and 64 were confirmed. These include 18 single-exon deletions out of which 8 were exclusively detected by HMZDelFinder and not by any of seven other CNV detection tools examined. Further investigation of the 64 validated deletion calls revealed at least 15 pathogenic HMZ deletions. Of those, 7 accounted for 17-50% of pathogenic CNVs in different disease cohorts where 7.1-11% of the molecular diagnosis solved rate was attributed to CNVs. In summary, we present an algorithm to detect rare, intragenic, single-exon deletion CNVs using WES data; this tool can be useful for disease gene discovery efforts and clinical WES analyses.

Primary immunodeficiency diseases: Genomic approaches delineate heterogeneous Mendelian disorders.

BACKGROUND: Primary immunodeficiency diseases (PIDDs) are clinically and genetically heterogeneous disorders thus far associated with mutations in more than 300 genes. The clinical phenotypes derived from distinct genotypes can overlap. Genetic etiology can be a prognostic indicator of disease severity and can influence treatment decisions. OBJECTIVE: We sought to investigate the ability of whole-exome screening methods to detect disease-causing variants in patients with PIDDs. METHODS: Patients with PIDDs from 278 families from 22 countries were investigated by using whole-exome sequencing. Computational copy number variant (CNV) prediction pipelines and an exome-tiling chromosomal microarray were also applied to identify intragenic CNVs. Analytic approaches initially focused on 475 known or candidate PIDD genes but were nonexclusive and further tailored based on clinical data, family history, and immunophenotyping. RESULTS: A likely molecular diagnosis was achieved in 110 (40%) unrelated probands. Clinical diagnosis was revised in about half (60/110) and management was directly altered in nearly a quarter (26/110) of families based on molecular findings. Twelve PIDD-causing CNVs were detected, including 7 smaller than 30 Kb that would not have been detected with conventional diagnostic CNV arrays. CONCLUSION: This high-throughput genomic approach enabled detection of disease-related variants in unexpected genes; permitted detection of low-grade constitutional, somatic, and revertant mosaicism; and provided evidence of a mutational burden in mixed PIDD immunophenotypes.

Increased bone turnover, osteoporosis, progressive tibial bowing, fractures, and scoliosis in a patient with a final-exon SATB2 frameshift mutation.

Haploinsufficiency of SATB2 causes cleft palate, intellectual disability with deficient speech, facial and dental abnormalities, and other variable features known collectively as SATB2-associated syndrome. This phenotype was accompanied by osteoporosis, fractures, and tibial bowing in two previously reported adult patients; each possessed SATB2 mutations either predicted or demonstrated to escape nonsense-mediated decay, suggesting that the additional bone defects result from a dominant negative effect and/or age-dependent penetrance. These hypotheses remain to be confirmed, as do the specific downstream defects causing bone abnormalities. We report a SATB2 mutation (c.2018dupA; p.(H673fs)) in a 15-year-old patient whose SATB2-associated syndrome phenotype is accompanied by osteoporosis, fractures, progressive tibial bowing, and scoliosis. As this homeodomain-disrupting and predicted truncating mutation resides within the final exon of SATB2, escape from nonsense-mediated decay is likely. Thus, we provide further evidence of bone phenotypes beyond those typically associated with SATB2-associated syndrome in individuals with potential dominant-negative SATB2 alleles, as well as evidence for age-dependence of bone features. Elevations in alkaline phosphatase, urinary N-telopeptide/creatinine ratio, and osteocalcin in the patient indicate increased bone turnover. We propose surveillance and treatment with osteoclast inhibitors to prevent fractures and to slow progressive bone deformities. © 2016 Wiley Periodicals, Inc.

A visual and curatorial approach to clinical variant prioritization and disease gene discovery in genome-wide diagnostics.

BACKGROUND: Genome-wide data are increasingly important in the clinical evaluation of human disease. However, the large number of variants observed in individual patients challenges the efficiency and accuracy of diagnostic review. Recent work has shown that systematic integration of clinical phenotype data with genotype information can improve diagnostic workflows and prioritization of filtered rare variants. We have developed visually interactive, analytically transparent analysis software that leverages existing disease catalogs, such as the Online Mendelian Inheritance in Man database (OMIM) and the Human Phenotype Ontology (HPO), to integrate patient phenotype and variant data into ranked diagnostic alternatives. METHODS: Our tool, "OMIM Explorer" ( http://www.omimexplorer.com ), extends the biomedical application of semantic similarity methods beyond those reported in previous studies. The tool also provides a simple interface for translating free-text clinical notes into HPO terms, enabling clinical providers and geneticists to contribute phenotypes to the diagnostic process. The visual approach uses semantic similarity with multidimensional scaling to collapse high-dimensional phenotype and genotype data from an individual into a graphical format that contextualizes the patient within a low-dimensional disease map. The map proposes a differential diagnosis and algorithmically suggests potential alternatives for phenotype queries--in essence, generating a computationally assisted differential diagnosis informed by the individual's personal genome. Visual interactivity allows the user to filter and update variant rankings by interacting with intermediate results. The tool also implements an adaptive approach for disease gene discovery based on patient phenotypes. RESULTS: We retrospectively analyzed pilot cohort data from the Baylor Miraca Genetics Laboratory, demonstrating performance of the tool and workflow in the re-analysis of clinical exomes. Our tool assigned to clinically reported variants a median rank of 2, placing causal variants in the top 1 % of filtered candidates across the 47 cohort cases with reported molecular diagnoses of exome variants in OMIM Morbidmap genes. Our tool outperformed Phen-Gen, eXtasy, PhenIX, PHIVE, and hiPHIVE in the prioritization of these clinically reported variants. CONCLUSIONS: Our integrative paradigm can improve efficiency and, potentially, the quality of genomic medicine by more effectively utilizing available phenotype information, catalog data, and genomic knowledge.

Hutterite-type cataract maps to chromosome 6p21.32-p21.31, cosegregates with a homozygous mutation in LEMD2, and is associated with sudden cardiac death.

BACKGROUND: Juvenile-onset cataracts are known among the Hutterites of North America. Despite being identified over 30 years ago, this autosomal recessive condition has not been mapped, and the disease gene is unknown. METHODS: We performed whole exome sequencing of three Hutterite-type cataract trios and follow-up genotyping and mapping in four extended kindreds. RESULTS: Trio exomes enabled genome-wide autozygosity mapping, which localized the disease gene to a 9.5-Mb region on chromosome 6p. This region contained two candidate variants, LEMD2 c.T38G and MUC21 c.665delC. Extended pedigrees recruited for variant genotyping revealed multiple additional relatives with juvenile-onset cataract, as well as six deceased relatives with both cataracts and sudden cardiac death. The candidate variants were genotyped in 84 family members, including 17 with cataracts; only the variant in LEMD2 cosegregated with cataracts (LOD = 9.62). SNP-based fine mapping within the 9.5 Mb linked region supported this finding by refining the cataract locus to a 0.5- to 2.9-Mb subregion (6p21.32-p21.31) containing LEMD2 but not MUC21. LEMD2 is expressed in mouse and human lenses and encodes a LEM domain-containing protein; the c.T38G missense mutation is predicted to mutate a highly conserved residue within this domain (p.Leu13Arg). CONCLUSION: We performed a genetic and genomic study of Hutterite-type cataract and found evidence for an association of this phenotype with sudden cardiac death. Using combined genetic and genomic approaches, we mapped cataracts to a small portion of chromosome 6 and propose that they result from a homozygous missense mutation in LEMD2.

Analysis of the ABCA4 genomic locus in Stargardt disease.

Autosomal recessive Stargardt disease (STGD1, MIM 248200) is caused by mutations in the ABCA4 gene. Complete sequencing of ABCA4 in STGD patients identifies compound heterozygous or homozygous disease-associated alleles in 65-70% of patients and only one mutation in 15-20% of patients. This study was designed to find the missing disease-causing ABCA4 variation by a combination of next-generation sequencing (NGS), array-Comparative Genome Hybridization (aCGH) screening, familial segregation and in silico analyses. The entire 140 kb ABCA4 genomic locus was sequenced in 114 STGD patients with one known ABCA4 exonic mutation revealing, on average, 200 intronic variants per sample. Filtering of these data resulted in 141 candidates for new mutations. Two variants were detected in four samples, two in three samples, and 20 variants in two samples, the remaining 117 new variants were detected only once. Multimodal analysis suggested 12 new likely pathogenic intronic ABCA4 variants, some of which were specific to (isolated) ethnic groups. No copy number variation (large deletions and insertions) was detected in any patient suggesting that it is a very rare event in the ABCA4 locus. Many variants were excluded since they were not conserved in non-human primates, were frequent in African populations and, therefore, represented ancestral, and not disease-associated, variants. The sequence variability in the ABCA4 locus is extensive and the non-coding sequences do not harbor frequent mutations in STGD patients of European-American descent. Defining disease-associated alleles in the ABCA4 locus requires exceptionally well characterized large cohorts and extensive analyses by a combination of various approaches.

The Alu-rich genomic architecture of SPAST predisposes to diverse and functionally distinct disease-associated CNV alleles.

Intragenic copy-number variants (CNVs) contribute to the allelic spectrum of both Mendelian and complex disorders. Although pathogenic deletions and duplications in SPAST (mutations in which cause autosomal-dominant spastic paraplegia 4 [SPG4]) have been described, their origins and molecular consequences remain obscure. We mapped breakpoint junctions of 54 SPAST CNVs at nucleotide resolution. Diverse combinations of exons are deleted or duplicated, highlighting the importance of particular exons for spastin function. Of the 54 CNVs, 38 (70%) appear to be mediated by an Alu-based mechanism, suggesting that the Alu-rich genomic architecture of SPAST renders this locus susceptible to various genome rearrangements. Analysis of breakpoint Alus further informs a model of Alu-mediated CNV formation characterized by small CNV size and potential involvement of mechanisms other than homologous recombination. Twelve deletions (22%) overlap part of SPAST and a portion of a nearby, directly oriented gene, predicting novel chimeric genes in these subjects' genomes. cDNA from a subject with a SPAST final exon deletion contained multiple SPAST:SLC30A6 fusion transcripts, indicating that SPAST CNVs can have transcriptional effects beyond the gene itself. SLC30A6 has been implicated in Alzheimer disease, so these fusion gene data could explain a report of spastic paraplegia and dementia cosegregating in a family with deletion of the final exon of SPAST. Our findings provide evidence that the Alu genomic architecture of SPAST predisposes to diverse CNV alleles with distinct transcriptional--and possibly phenotypic--consequences. Moreover, we provide further mechanistic insights into Alu-mediated copy-number change that are extendable to other loci.

Exonic duplication CNV of NDRG1 associated with autosomal-recessive HMSN-Lom/CMT4D.

PURPOSE: Copy-number variations as a mutational mechanism contribute significantly to human disease. Approximately one-half of the patients with Charcot-Marie-Tooth (CMT) disease have a 1.4 Mb duplication copy-number variation as the cause of their neuropathy. However, non-CMT1A neuropathy patients rarely have causative copy-number variations, and to date, autosomal-recessive disease has not been associated with copy-number variation as a mutational mechanism. METHODS: We performed Agilent 8 × 60 K array comparative genomic hybridization on DNA from 12 recessive Turkish families with CMT disease. Additional molecular studies were conducted to detect breakpoint junctions and to evaluate gene expression levels in a family in which we detected an intragenic duplication copy-number variation. RESULTS: We detected an ~6.25 kb homozygous intragenic duplication in NDRG1, a gene known to be causative for recessive HMSNL/CMT4D, in three individuals from a Turkish family with CMT neuropathy. Further studies showed that this intragenic copy-number variation resulted in a homozygous duplication of exons 6-8 that caused decreased mRNA expression of NDRG1. CONCLUSION: Exon-focused high-resolution array comparative genomic hybridization enables the detection of copy-number variation carrier states in recessive genes, particularly small copy-number variations encompassing or disrupting single genes. In families for whom a molecular diagnosis has not been elucidated by conventional clinical assays, an assessment for copy-number variations in known CMT genes might be considered.

Platform comparison of detecting copy number variants with microarrays and whole-exome sequencing.

Copy number variation (CNV) is a common source of genetic variation that has been implicated in many genomic disorders, Mendelian diseases, and common/complex traits. Genomic microarrays are often employed for CNV detection. More recently, whole-exome sequencing (WES) has enabled detection of clinically relevant point mutations and small insertion-deletion exome wide. We evaluated (de Ligt et al. 2013) [1] the utility of short-read WES (SOLiD 5500xl) to detect clinically relevant CNVs in DNA from 10 patients with intellectual disability and compared these results to data from three independent high-resolution microarray platforms. Calls made by the different platforms and detection software are available at dbVar under nstd84.