Comprehensive Gene Panel Testing for Hearing Loss in Children: Understanding Factors Influencing Diagnostic Yield.

OBJECTIVE: To evaluate factors influencing the diagnostic yield of comprehensive gene panel testing (CGPT) for hearing loss (HL) in children and to understand the characteristics of undiagnosed probands. STUDY DESIGN: This was a retrospective cohort study of 474 probands with childhood-onset HL who underwent CGPT between 2016 and 2020 at a single center. Main outcomes and measures included the association between clinical variables and diagnostic yield and the genetic and clinical characteristics of undiagnosed probands. RESULTS: The overall diagnostic yield was 44% (209/474) with causative variants involving 41 genes. While the diagnostic yield was high in the probands with congenital, bilateral, and severe HL, it was low in those with unilateral, noncongenital, or mild HL; cochlear nerve deficiency; preterm birth; neonatal intensive care unit admittance; certain ancestry; and developmental delay. Follow-up studies on 49 probands with initially inconclusive CGPT results changed the diagnostic status to likely positive or negative outcomes in 39 of them (80%). Reflex to exome sequencing on 128 undiagnosed probands by CGPT revealed diagnostic findings in 8 individuals, 5 of whom had developmental delays. The remaining 255 probands were undiagnosed, with 173 (173/255) having only a single variant in the gene(s) associated with autosomal recessive HL and 28% (48/173) having a matched phenotype. CONCLUSION: CGPT efficiently identifies the genetic etiologies of HL in children. CGPT-undiagnosed probands may benefit from follow-up studies or expanded testing.

Developmental genomics of limb malformations: Allelic series in association with gene dosage effects contribute to the clinical variability.

Genetic heterogeneity, reduced penetrance, and variable expressivity, the latter including asymmetric body axis plane presentations, have all been described in families with congenital limb malformations (CLMs). Interfamilial and intrafamilial heterogeneity highlight the complexity of the underlying genetic pathogenesis of these developmental anomalies. Family-based genomics by exome sequencing (ES) and rare variant analyses combined with whole-genome array-based comparative genomic hybridization were implemented to investigate 18 families with limb birth defects. Eleven of 18 (61%) families revealed explanatory variants, including 7 single-nucleotide variant alleles and 3 copy number variants (CNVs), at previously reported "disease trait associated loci": BHLHA9, GLI3, HOXD cluster, HOXD13, NPR2, and WNT10B. Breakpoint junction analyses for all three CNV alleles revealed mutational signatures consistent with microhomology-mediated break-induced replication, a mechanism facilitated by Alu/Alu-mediated rearrangement. Homozygous duplication of BHLHA9 was observed in one Turkish kindred and represents a novel contributory genetic mechanism to Gollop-Wolfgang Complex (MIM: 228250), where triplication of the locus has been reported in one family from Japan (i.e., 4n = 2n + 2n versus 4n = 3n + 1n allelic configurations). Genes acting on limb patterning are sensitive to a gene dosage effect and are often associated with an allelic series. We extend an allele-specific gene dosage model to potentially assist, in an adjuvant way, interpretations of interconnections among an allelic series, clinical severity, and reduced penetrance of the BHLHA9-related CLM spectrum.

Variant-specific changes in RAC3 function disrupt corticogenesis in neurodevelopmental phenotypes.

Variants in RAC3, encoding a small GTPase RAC3 which is critical for the regulation of actin cytoskeleton and intracellular signal transduction, are associated with a rare neurodevelopmental disorder with structural brain anomalies and facial dysmorphism. We investigated a cohort of 10 unrelated participants presenting with global psychomotor delay, hypotonia, behavioural disturbances, stereotyped movements, dysmorphic features, seizures and musculoskeletal abnormalities. MRI of brain revealed a complex pattern of variable brain malformations, including callosal abnormalities, white matter thinning, grey matter heterotopia, polymicrogyria/dysgyria, brainstem anomalies and cerebellar dysplasia. These patients harboured eight distinct de novo RAC3 variants, including six novel variants (NM_005052.3): c.34G > C p.G12R, c.179G > A p.G60D, c.186_188delGGA p.E62del, c.187G > A p.D63N, c.191A > G p.Y64C and c.348G > C p.K116N. We then examined the pathophysiological significance of these novel and previously reported pathogenic variants p.P29L, p.P34R, p.A59G, p.Q61L and p.E62K. In vitro analyses revealed that all tested RAC3 variants were biochemically and biologically active to variable extent, and exhibited a spectrum of different affinities to downstream effectors including p21-activated kinase 1. We then focused on the four variants p.Q61L, p.E62del, p.D63N and p.Y64C in the Switch II region, which is essential for the biochemical activity of small GTPases and also a variation hot spot common to other Rho family genes, RAC1 and CDC42. Acute expression of the four variants in embryonic mouse brain using in utero electroporation caused defects in cortical neuron morphology and migration ending up with cluster formation during corticogenesis. Notably, defective migration by p.E62del, p.D63N and p.Y64C were rescued by a dominant negative version of p21-activated kinase 1. Our results indicate that RAC3 variants result in morphological and functional defects in cortical neurons during brain development through variant-specific mechanisms, eventually leading to heterogeneous neurodevelopmental phenotypes.

Phenotypic and mutational spectrum of ROR2-related Robinow syndrome.

Robinow syndrome is characterized by a triad of craniofacial dysmorphisms, disproportionate-limb short stature, and genital hypoplasia. A significant degree of phenotypic variability seems to correlate with different genes/loci. Disturbances of the noncanonical WNT-pathway have been identified as the main cause of the syndrome. Biallelic variants in ROR2 cause an autosomal recessive form of the syndrome with distinctive skeletal findings. Twenty-two patients with a clinical diagnosis of autosomal recessive Robinow syndrome were screened for variants in ROR2 using multiple molecular approaches. We identified 25 putatively pathogenic ROR2 variants, 16 novel, including single nucleotide variants and exonic deletions. Detailed phenotypic analyses revealed that all subjects presented with a prominent forehead, hypertelorism, short nose, abnormality of the nasal tip, brachydactyly, mesomelic limb shortening, short stature, and genital hypoplasia in male patients. A total of 19 clinical features were present in more than 75% of the subjects, thus pointing to an overall uniformity of the phenotype. Disease-causing variants in ROR2, contribute to a clinically recognizable autosomal recessive trait phenotype with multiple skeletal defects. A comprehensive quantitative clinical evaluation of this cohort delineated the phenotypic spectrum of ROR2-related Robinow syndrome. The identification of exonic deletion variant alleles further supports the contention of a loss-of-function mechanism in the etiology of the syndrome.

High prevalence of multilocus pathogenic variation in neurodevelopmental disorders in the Turkish population.

Neurodevelopmental disorders (NDDs) are clinically and genetically heterogenous; many such disorders are secondary to perturbation in brain development and/or function. The prevalence of NDDs is > 3%, resulting in significant sociocultural and economic challenges to society. With recent advances in family-based genomics, rare-variant analyses, and further exploration of the Clan Genomics hypothesis, there has been a logarithmic explosion in neurogenetic "disease-associated genes" molecular etiology and biology of NDDs; however, the majority of NDDs remain molecularly undiagnosed. We applied genome-wide screening technologies, including exome sequencing (ES) and whole-genome sequencing (WGS), to identify the molecular etiology of 234 newly enrolled subjects and 20 previously unsolved Turkish NDD families. In 176 of the 234 studied families (75.2%), a plausible and genetically parsimonious molecular etiology was identified. Out of 176 solved families, deleterious variants were identified in 218 distinct genes, further documenting the enormous genetic heterogeneity and diverse perturbations in human biology underlying NDDs. We propose 86 candidate disease-trait-associated genes for an NDD phenotype. Importantly, on the basis of objective and internally established variant prioritization criteria, we identified 51 families (51/176 = 28.9%) with multilocus pathogenic variation (MPV), mostly driven by runs of homozygosity (ROHs) - reflecting genomic segments/haplotypes that are identical-by-descent. Furthermore, with the use of additional bioinformatic tools and expansion of ES to additional family members, we established a molecular diagnosis in 5 out of 20 families (25%) who remained undiagnosed in our previously studied NDD cohort emanating from Turkey.

Functional biology of the Steel syndrome founder allele and evidence for clan genomics derivation of COL27A1 pathogenic alleles worldwide.

Previously we reported the identification of a homozygous COL27A1 (c.2089G>C; p.Gly697Arg) missense variant and proposed it as a founder allele in Puerto Rico segregating with Steel syndrome (STLS, MIM #615155); a rare osteochondrodysplasia characterized by short stature, congenital bilateral hip dysplasia, carpal coalitions, and scoliosis. We now report segregation of this variant in five probands from the initial clinical report defining the syndrome and an additional family of Puerto Rican descent with multiple affected adult individuals. We modeled the orthologous variant in murine Col27a1 and found it recapitulates some of the major Steel syndrome associated skeletal features including reduced body length, scoliosis, and a more rounded skull shape. Characterization of the in vivo murine model shows abnormal collagen deposition in the extracellular matrix and disorganization of the proliferative zone of the growth plate. We report additional COL27A1 pathogenic variant alleles identified in unrelated consanguineous Turkish kindreds suggesting Clan Genomics and identity-by-descent homozygosity contributing to disease in this population. The hypothesis that carrier states for this autosomal recessive osteochondrodysplasia may contribute to common complex traits is further explored in a large clinical population cohort. Our findings augment our understanding of COL27A1 biology and its role in skeletal development; and expand the functional allelic architecture in this gene underlying both rare and common disease phenotypes.

Biallelic GRM7 variants cause epilepsy, microcephaly, and cerebral atrophy.

OBJECTIVE: Defects in ion channels and neurotransmitter receptors are implicated in developmental and epileptic encephalopathy (DEE). Metabotropic glutamate receptor 7 (mGluR7), encoded by GRM7, is a presynaptic G-protein-coupled glutamate receptor critical for synaptic transmission. We previously proposed GRM7 as a candidate disease gene in two families with neurodevelopmental disorders (NDDs). One additional family has been published since. Here, we describe three additional families with GRM7 biallelic variants and deeply characterize the associated clinical neurological and electrophysiological phenotype and molecular data in 11 affected individuals from six unrelated families. METHODS: Exome sequencing and family-based rare variant analyses on a cohort of 220 consanguineous families with NDDs revealed three families with GRM7 biallelic variants; three additional families were identified through literature search and collaboration with a clinical molecular laboratory. RESULTS: We compared the observed clinical features and variants of 11 affected individuals from the six unrelated families. Identified novel deleterious variants included two homozygous missense variants (c.2671G>A:p.Glu891Lys and c.1973G>A:p.Arg685Gln) and one homozygous stop-gain variant (c.1975C>T:p.Arg659Ter). Developmental delay, neonatal- or infantile-onset epilepsy, and microcephaly were universal. Three individuals had hypothalamic-pituitary-axis dysfunction without pituitary structural abnormality. Neuroimaging showed cerebral atrophy and hypomyelination in a majority of cases. Two siblings demonstrated progressive loss of myelination by 2 years in both and an acquired microcephaly pattern in one. Five individuals died in early or late childhood. CONCLUSION: Detailed clinical characterization of 11 individuals from six unrelated families demonstrates that rare biallelic GRM7 pathogenic variants can cause DEEs, microcephaly, hypomyelination, and cerebral atrophy.

Biallelic and De Novo Variants in DONSON Reveal a Clinical Spectrum of Cell Cycle-opathies with Microcephaly, Dwarfism and Skeletal Abnormalities.

Co-occurrence of primordial dwarfism and microcephaly together with particular skeletal findings are seen in a wide range of Mendelian syndromes including microcephaly micromelia syndrome (MMS, OMIM 251230), microcephaly, short stature, and limb abnormalities (MISSLA, OMIM 617604), and microcephalic primordial dwarfisms (MPDs). Genes associated with these syndromes encode proteins that have crucial roles in DNA replication or in other critical steps of the cell cycle that link DNA replication to cell division. We identified four unrelated families with five affected individuals having biallelic or de novo variants in DONSON presenting with a core phenotype of severe short stature (z score < -3 SD), additional skeletal abnormalities, and microcephaly. Two apparently unrelated families with identical homozygous c.631C > T p.(Arg211Cys) variant had clinical features typical of Meier-Gorlin syndrome (MGS), while two siblings with compound heterozygous c.346delG p.(Asp116Ile*62) and c.1349A > G p.(Lys450Arg) variants presented with Seckel-like phenotype. We also identified a de novo c.683G > T p.(Trp228Leu) variant in DONSON in a patient with prominent micrognathia, short stature and hypoplastic femur and tibia, clinically diagnosed with Femoral-Facial syndrome (FFS, OMIM 134780). Biallelic variants in DONSON have been recently described in individuals with microcephalic dwarfism. These studies also demonstrated that DONSON has an essential conserved role in the cell cycle. Here we describe novel biallelic and de novo variants that are associated with MGS, Seckel-like phenotype and FFS, the last of which has not been associated with any disease gene to date.

Paralog Studies Augment Gene Discovery: DDX and DHX Genes.

Members of a paralogous gene family in which variation in one gene is known to cause disease are eight times more likely to also be associated with human disease. Recent studies have elucidated DHX30 and DDX3X as genes for which pathogenic variant alleles are involved in neurodevelopmental disorders. We hypothesized that variants in paralogous genes encoding members of the DExD/H-box RNA helicase superfamily might also underlie developmental delay and/or intellectual disability (DD and/or ID) disease phenotypes. Here we describe 15 unrelated individuals who have DD and/or ID, central nervous system (CNS) dysfunction, vertebral anomalies, and dysmorphic features and were found to have probably damaging variants in DExD/H-box RNA helicase genes. In addition, these individuals exhibit a variety of other tissue and organ system involvement including ocular, outer ear, hearing, cardiac, and kidney tissues. Five individuals with homozygous (one), compound-heterozygous (two), or de novo (two) missense variants in DHX37 were identified by exome sequencing. We identified ten total individuals with missense variants in three other DDX/DHX paralogs: DHX16 (four individuals), DDX54 (three individuals), and DHX34 (three individuals). Most identified variants are rare, predicted to be damaging, and occur at conserved amino acid residues. Taken together, these 15 individuals implicate the DExD/H-box helicases in both dominantly and recessively inherited neurodevelopmental phenotypes and highlight the potential for more than one disease mechanism underlying these disorders.

The Genomics of Arthrogryposis, a Complex Trait: Candidate Genes and Further Evidence for Oligogenic Inheritance.

Arthrogryposis is a clinical finding that is present either as a feature of a neuromuscular condition or as part of a systemic disease in over 400 Mendelian conditions. The underlying molecular etiology remains largely unknown because of genetic and phenotypic heterogeneity. We applied exome sequencing (ES) in a cohort of 89 families with the clinical sign of arthrogryposis. Additional molecular techniques including array comparative genomic hybridization (aCGH) and Droplet Digital PCR (ddPCR) were performed on individuals who were found to have pathogenic copy number variants (CNVs) and mosaicism, respectively. A molecular diagnosis was established in 65.2% (58/89) of families. Eleven out of 58 families (19.0%) showed evidence for potential involvement of pathogenic variation at more than one locus, probably driven by absence of heterozygosity (AOH) burden due to identity-by-descent (IBD). RYR3, MYOM2, ERGIC1, SPTBN4, and ABCA7 represent genes, identified in two or more families, for which mutations are probably causative for arthrogryposis. We also provide evidence for the involvement of CNVs in the etiology of arthrogryposis and for the idea that both mono-allelic and bi-allelic variants in the same gene cause either similar or distinct syndromes. We were able to identify the molecular etiology in nine out of 20 families who underwent reanalysis. In summary, our data from family-based ES further delineate the molecular etiology of arthrogryposis, yielded several candidate disease-associated genes, and provide evidence for mutational burden in a biological pathway or network. Our study also highlights the importance of reanalysis of individuals with unsolved diagnoses in conjunction with sequencing extended family members.

Exome Sequencing of a Primary Ovarian Insufficiency Cohort Reveals Common Molecular Etiologies for a Spectrum of Disease.

CONTEXT: Primary ovarian insufficiency (POI) encompasses a spectrum of premature menopause, including both primary and secondary amenorrhea. For 75% to 90% of individuals with hypergonadotropic hypogonadism presenting as POI, the molecular etiology is unknown. Common etiologies include chromosomal abnormalities, environmental factors, and congenital disorders affecting ovarian development and function, as well as syndromic and nonsyndromic single gene disorders suggesting POI represents a complex trait. OBJECTIVE: To characterize the contribution of known disease genes to POI and identify molecular etiologies and biological underpinnings of POI. DESIGN, SETTING, AND PARTICIPANTS: We applied exome sequencing (ES) and family-based genomics to 42 affected female individuals from 36 unrelated Turkish families, including 31 with reported parental consanguinity. RESULTS: This analysis identified likely damaging, potentially contributing variants and molecular diagnoses in 16 families (44%), including 11 families with likely damaging variants in known genes and five families with predicted deleterious variants in disease genes (IGSF10, MND1, MRPS22, and SOHLH1) not previously associated with POI. Of the 16 families, 2 (13%) had evidence for potentially pathogenic variants at more than one locus. Absence of heterozygosity consistent with identity-by-descent mediated recessive disease burden contributes to molecular diagnosis in 15 of 16 (94%) families. GeneMatcher allowed identification of additional families from diverse genetic backgrounds. CONCLUSIONS: ES analysis of a POI cohort further characterized locus heterogeneity, reaffirmed the association of genes integral to meiotic recombination, demonstrated the likely contribution of genes involved in hypothalamic development, and documented multilocus pathogenic variation suggesting the potential for oligogenic inheritance contributing to the development of POI.

Identifying Genes Whose Mutant Transcripts Cause Dominant Disease Traits by Potential Gain-of-Function Alleles.

Premature termination codon (PTC)-bearing transcripts are often degraded by nonsense-mediated decay (NMD) resulting in loss-of-function (LoF) alleles. However, not all PTCs result in LoF mutations, i.e., some such transcripts escape NMD and are translated to truncated peptide products that result in disease due to gain-of-function (GoF) effects. Since the location of the PTC is a major factor determining transcript fate, we hypothesized that depletion of protein-truncating variants (PTVs) within the gene region predicted to escape NMD in control databases could provide a rank for genic susceptibility for disease through GoF versus LoF. We developed an NMD escape intolerance score to rank genes based on the depletion of PTVs that would render them able to escape NMD using the Atherosclerosis Risk in Communities Study (ARIC) and the Exome Aggregation Consortium (ExAC) control databases, which was further used to screen the Baylor-Center for Mendelian Genomics disease database. This analysis revealed 1,996 genes significantly depleted for PTVs that are predicted to escape from NMD, i.e., PTVesc; further studies provided evidence that revealed a subset as candidate genes underlying Mendelian phenotypes. Importantly, these genes have characteristically low pLI scores, which can cause them to be overlooked as candidates for dominant diseases. Collectively, we demonstrate that this NMD escape intolerance score is an effective and efficient tool for gene discovery in Mendelian diseases due to production of truncated or altered proteins. More importantly, we provide a complementary analytical tool to aid identification of genes associated with dominant traits through a mechanism distinct from LoF.

Phenotypic expansion illuminates multilocus pathogenic variation.

PURPOSE: Multilocus variation-pathogenic variants in two or more disease genes-can potentially explain the underlying genetic basis for apparent phenotypic expansion in cases for which the observed clinical features extend beyond those reported in association with a "known" disease gene. METHODS: Analyses focused on 106 patients, 19 for whom apparent phenotypic expansion was previously attributed to variation at known disease genes. We performed a retrospective computational reanalysis of whole-exome sequencing data using stringent Variant Call File filtering criteria to determine whether molecular diagnoses involving additional disease loci might explain the observed expanded phenotypes. RESULTS: Multilocus variation was identified in 31.6% (6/19) of families with phenotypic expansion and 2.3% (2/87) without phenotypic expansion. Intrafamilial clinical variability within two families was explained by multilocus variation identified in the more severely affected sibling. CONCLUSION: Our findings underscore the role of multiple rare variants at different loci in the etiology of genetically and clinically heterogeneous cohorts. Intrafamilial phenotypic and genotypic variability allowed a dissection of genotype-phenotype relationships in two families. Our data emphasize the critical role of the clinician in diagnostic genomic analyses and demonstrate that apparent phenotypic expansion may represent blended phenotypes resulting from pathogenic variation at more than one locus.

Mutations in the mitochondrial ribosomal protein MRPS22 lead to primary ovarian insufficiency.

Primary ovarian insufficiency (POI) is characterized by amenorrhea and loss or dysfunction of ovarian follicles prior to the age of 40. POI has been associated with autosomal recessive mutations in genes involving hormonal signaling and folliculogenesis, however, the genetic etiology of POI most often remains unknown. Here we report MRPS22 homozygous missense variants c.404G>A (p.R135Q) and c.605G>A (p.R202H) identified in four females from two independent consanguineous families as a novel genetic cause of POI in adolescents. Both missense mutations identified in MRPS22 are rare, occurred in highly evolutionarily conserved residues, and are predicted to be deleterious to protein function. In contrast to prior reports of mutations in MRPS22 associated with severe mitochondrial disease, the POI phenotype is far less severe. Consistent with this genotype-phenotype correlation, mitochondrial defects in oxidative phosphorylation or rRNA levels were not detected in fibroblasts derived from the POI patients, suggesting a non-bioenergetic or tissue-specific mitochondrial defect. Furthermore, we demonstrate in a Drosophila model that mRpS22 deficiency specifically in somatic cells of the ovary had no effect on fertility, whereas flies with mRpS22 deficiency specifically in germ cells were infertile and agametic, demonstrating a cell autonomous requirement for mRpS22 in germ cell development. These findings collectively identify that MRPS22, a component of the small mitochondrial ribosome subunit, is critical for ovarian development and may therefore provide insight into the pathophysiology and treatment of ovarian dysfunction.

Expanding the Phenotype of Homozygous KCNMA1 Mutations; Dyskinesia, Epilepsy, Intellectual Disability, Cerebellar and Corticospinal Tract Atrophy

BACKGROUND: The KCNMA1 gene encodes the α-subunit of the large conductance, voltage, and calcium-sensitive potassium channel (BK channels) that plays a critical role in neuronal excitability. Heterozygous mutations in KCNMA1 were first illustrated in a large family with generalized epilepsy and paroxysmal nonkinesigenic dyskinesia. Recent research has established homozygous KCNMA1 mutations accountable for the phenotype of cerebellar atrophy, developmental delay, and seizures. CASE REPORT: Here, we report the case of a patient with a novel homozygous truncating mutation in KCNMA1 (p.Arg458Ter) presenting with both the loss- and gain-of-function phenotype with paroxysmal dyskinesia, epilepsy, intellectual delay, and corticospinal­cerebellar tract atrophy. CONCLUSION: This report extends the KNCMA1 mutation phenotype with a patient who carries a novel frameshift variant, presenting with both the gain- and loss-of-function phenotypes along with spinal tract involvement as a novel characteristic.

A biallelic ANTXR1 variant expands the anthrax toxin receptor associated phenotype to tooth agenesis.

Tooth development is regulated by multiple genetic pathways, which ultimately drive the complex interactions between the oral epithelium and mesenchyme. Disruptions at any time point during this process may lead to failure of tooth development, also known as tooth agenesis (TA). TA is a common craniofacial abnormality in humans and represents the failure to develop one or more permanent teeth. Many genes and potentially subtle variants in these genes contribute to the TA phenotype. We report the clinical and genetic impact of a rare homozygous ANTXR1 variant (c.1312C>T), identified by whole exome sequencing (WES), in a consanguineous Turkish family with TA. Mutations in ANTXR1 have been associated with GAPO (growth retardation, alopecia, pseudoanodontia, and optic atrophy) syndrome and infantile hemangioma, however no clinical characteristics associated with these conditions were observed in our study family. We detected the expression of Antxr1 in oral and dental tissues of developing mouse embryos, further supporting a role for this gene in tooth development. Our findings implicate ANTXR1 as a candidate gene for isolated TA, suggest the involvement of specific hypomorphic alleles, and expand the previously known ANTXR1-associated phenotypes.

WNT Signaling Perturbations Underlie the Genetic Heterogeneity of Robinow Syndrome.

Locus heterogeneity characterizes a variety of skeletal dysplasias often due to interacting or overlapping signaling pathways. Robinow syndrome is a skeletal disorder historically refractory to molecular diagnosis, potentially stemming from substantial genetic heterogeneity. All current known pathogenic variants reside in genes within the noncanonical Wnt signaling pathway including ROR2, WNT5A, and more recently, DVL1 and DVL3. However, ∼70% of autosomal-dominant Robinow syndrome cases remain molecularly unsolved. To investigate this missing heritability, we recruited 21 families with at least one family member clinically diagnosed with Robinow or Robinow-like phenotypes and performed genetic and genomic studies. In total, four families with variants in FZD2 were identified as well as three individuals from two families with biallelic variants in NXN that co-segregate with the phenotype. Importantly, both FZD2 and NXN are relevant protein partners in the WNT5A interactome, supporting their role in skeletal development. In addition to confirming that clustered -1 frameshifting variants in DVL1 and DVL3 are the main contributors to dominant Robinow syndrome, we also found likely pathogenic variants in candidate genes GPC4 and RAC3, both linked to the Wnt signaling pathway. These data support an initial hypothesis that Robinow syndrome results from perturbation of the Wnt/PCP pathway, suggest specific relevant domains of the proteins involved, and reveal key contributors in this signaling cascade during human embryonic development. Contrary to the view that non-allelic genetic heterogeneity hampers gene discovery, this study demonstrates the utility of rare disease genomic studies to parse gene function in human developmental pathways.

REST Final-Exon-Truncating Mutations Cause Hereditary Gingival Fibromatosis.

Hereditary gingival fibromatosis (HGF) is the most common genetic form of gingival fibromatosis that develops as a slowly progressive, benign, localized or generalized enlargement of keratinized gingiva. HGF is a genetically heterogeneous disorder and can be transmitted either as an autosomal-dominant or autosomal-recessive trait or appear sporadically. To date, four loci (2p22.1, 2p23.3-p22.3, 5q13-q22, and 11p15) have been mapped to autosomes and one gene (SOS1) has been associated with the HGF trait observed to segregate in a dominant inheritance pattern. Here we report 11 individuals with HGF from three unrelated families. Whole-exome sequencing (WES) revealed three different truncating mutations including two frameshifts and one nonsense variant in RE1-silencing transcription factor (REST) in the probands from all families and further genetic and genomic analyses confirmed the WES-identified findings. REST is a transcriptional repressor that is expressed throughout the body; it has different roles in different cellular contexts, such as oncogenic and tumor-suppressor functions and hematopoietic and cardiac differentiation. Here we show the consequences of germline final-exon-truncating mutations in REST for organismal development and the association with the HGF phenotype.

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.

De novo missense variants in PPP1CB are associated with intellectual disability and congenital heart disease.

Intellectual disabilities are genetically heterogeneous and can be associated with congenital anomalies. Using whole-exome sequencing (WES), we identified five different de novo missense variants in the protein phosphatase-1 catalytic subunit beta (PPP1CB) gene in eight unrelated individuals who share an overlapping phenotype of dysmorphic features, macrocephaly, developmental delay or intellectual disability (ID), congenital heart disease, short stature, and skeletal and connective tissue abnormalities. Protein phosphatase-1 (PP1) is a serine/threonine-specific protein phosphatase involved in the dephosphorylation of a variety of proteins. The PPP1CB gene encodes a PP1 subunit that regulates the level of protein phosphorylation. All five altered amino acids we observed are highly conserved among the PP1 subunit family, and all are predicted to disrupt PP1 subunit binding and impair dephosphorylation. Our data suggest that our heterozygous de novo PPP1CB pathogenic variants are associated with syndromic intellectual disability.