The application of next-generation sequencing in the autozygosity mapping of human recessive diseases

Genomics of rare diseases in the Greater Middle East

The Greater Middle East (GME) represents a concentrated region of unparalleled genetic diversity, characterized by an abundance of distinct alleles, founder mutations and extensive autozygosity driven by high consanguinity rates. These genetic hallmarks present a unique, yet vastly untapped resource for genomic research on Mendelian diseases. Despite this immense potential, the GME continues to face substantial challenges in comprehensive data collection and analysis. This Perspective highlights the region's unique position as a natural laboratory for genetic discovery and explores the challenges that have stifled progress thus far. Importantly, we propose strategic solutions, advocating for an all-inclusive research approach. With targeted investment and focused efforts, the latent genetic wealth in the GME can be transformed into a global hub for genomic research that will redefine and advance our understanding of the human genome.

Bi-allelic variants in WDR47 cause a complex neurodevelopmental syndrome

Brain development requires the coordinated growth of structures and cues that are essential for forming neural circuits and cognitive functions. The corpus callosum, the largest interhemispheric connection, is formed by the axons of callosal projection neurons through a series of tightly regulated cellular events, including neuronal specification, migration, axon extension and branching. Defects in any of those steps can lead to a range of disorders known as syndromic corpus callosum dysgenesis (CCD). We report five unrelated families carrying bi-allelic variants in WDR47 presenting with CCD together with other neuroanatomical phenotypes such as microcephaly and enlarged ventricles. Using in vitro and in vivo mouse models and complementation assays, we show that WDR47 is required for survival of callosal neurons by contributing to the maintenance of mitochondrial and microtubule homeostasis. We further propose that severity of the CCD phenotype is determined by the degree of the loss of function caused by the human variants. Taken together, we identify WDR47 as a causative gene of a new neurodevelopmental syndrome characterized by corpus callosum abnormalities and other neuroanatomical malformations.

Large-scale genomic investigation of pediatric cholestasis reveals a novel hepatorenal ciliopathy caused by PSKH1 mutations

PURPOSE

Pediatric cholestasis is the phenotypic expression of clinically and genetically heterogeneous disorders of bile acid synthesis and flow. Although a growing number of monogenic causes of pediatric cholestasis have been identified, the majority of cases remain undiagnosed molecularly.

METHODS

In a cohort of 299 pediatric participants (279 families) with intrahepatic cholestasis, we performed exome sequencing as a first-tier diagnostic test.

RESULTS

A likely causal variant was identified in 135 families (48.56%). These comprise 135 families that harbor variants spanning 37 genes with established or tentative links to cholestasis. In addition, we propose a novel candidate gene (PSKH1) (HGNC:9529) in 4 families. PSKH1 was particularly compelling because of strong linkage in 3 consanguineous families who shared a novel hepatorenal ciliopathy phenotype. Two of the 4 families shared a founder homozygous variant, whereas the third and fourth had different homozygous variants in PSKH1. PSKH1 encodes a putative protein serine kinase of unknown function. Patient fibroblasts displayed abnormal cilia that are long and show abnormal transport. A homozygous Pskh1 mutant mouse faithfully recapitulated the human phenotype and displayed abnormally long cilia. The phenotype could be rationalized by the loss of catalytic activity observed for each recombinant PSKH1 variant using in vitro kinase assays.

CONCLUSION

Our results support the use of genomics in the workup of pediatric cholestasis and reveal PSKH1-related hepatorenal ciliopathy as a novel candidate monogenic form.

NanoRanger enables rapid single-base-pair resolution of genomic disorders

BACKGROUND

Delineating base-resolution breakpoints of complex rearrangements is crucial for an accurate clinical understanding of pathogenic variants and for carrier screening within family networks or the broader population. However, despite advances in genetic testing using short-read sequencing (SRS), this task remains costly and challenging.

METHODS

This study addresses the challenges of resolving missing disease-causing breakpoints in complex genomic disorders with suspected homozygous rearrangements by employing multiple long-read sequencing (LRS) strategies, including a novel and efficient strategy named nanopore-based rapid acquisition of neighboring genomic regions (NanoRanger). NanoRanger does not require large amounts of ultrahigh-molecular-weight DNA and stands out for its ease of use and rapid acquisition of large genomic regions of interest with deep coverage.

FINDINGS

We describe a cohort of 16 familial cases, each harboring homozygous rearrangements that defied breakpoint determination by SRS and optical genome mapping (OGM). NanoRanger identified the breakpoints with single-base-pair resolution, enabling accurate determination of the carrier status of unaffected family members as well as the founder nature of these genomic lesions and their frequency in the local population. The resolved breakpoints revealed that repetitive DNA, gene regulatory elements, and transcription activity contribute to genome instability in these novel recessive rearrangements.

CONCLUSIONS

Our data suggest that NanoRanger greatly improves the success rate of resolving base-resolution breakpoints of complex genomic disorders and expands access to LRS for the benefit of patients with Mendelian disorders.

FUNDING

M.L. is supported by KAUST Baseline Award no. BAS/1/1080-01-01 and KAUST Research Translation Fund Award no. REI/1/4742-01.

Autozygome-guided exome-first study in a consanguineous cohort with early-onset retinal disease uncovers an isolated RIMS2 phenotype and a retina-enriched RIMS2 isoform

Leber congenital amaurosis (LCA) and early-onset retinal degeneration (EORD) are inherited retinal diseases (IRD) characterized by early-onset vision impairment. Herein, we studied 15 Saudi families by whole exome sequencing (WES) and run-of-homozygosity (ROH) detection via AutoMap in 12/15 consanguineous families. This revealed (likely) pathogenic variants in 11/15 families (73%). A potential founder variant was found in RPGRIP1. Homozygous pathogenic variants were identified in known IRD genes (ATF6, CRB1, CABP4, RDH12, RIMS2, RPGRIP1, SPATA7). We established genotype-driven clinical reclassifications for ATF6, CABP4, and RIMS2. Specifically, we observed isolated IRD in the individual with the novel RIMS2 variant, and we found a retina-enriched RIMS2 isoform conserved but not annotated in mouse. The latter illustrates potential different phenotypic consequences of pathogenic variants depending on the particular tissue/cell-type specific isoforms they affect. Lastly, a compound heterozygous genotype in GUCY2D in one non-consanguineous family was demonstrated, and homozygous variants in novel candidate genes ATG2B and RUFY3 were found in the two remaining consanguineous families. Reporting these genes will allow to validate them in other IRD cohorts. Finally, the missing heritability of the two unsolved IRD cases may be attributed to variants in non-coding regions or structural variants that remained undetected, warranting future WGS studies.

Founder mutations and rare disease in the Arab world

Founder mutations are disease-causing variants that occur frequently in geographically or culturally isolated groups whose shared ancestor(s) carried the pathogenic variant. While some disease alleles may vanish from the genetic pool due to natural selection, variants with weaker effects may survive for a long time, thereby enhancing the prevalence of some rare diseases. These are predominantly autosomal recessive diseases but can also be autosomal dominant traits with late-onset or mild phenotypes. Cultural practices, such as endogamy and consanguinity, in these isolated groups lead to higher prevalence of such rare diseases compared to the rest of the population and worldwide. In this Perspective, we define population isolates and the underlying genetic mechanisms for accumulating founder mutations. We also discuss the current and potential scientific, clinical and public-health implications of studying founder mutations in population isolates around the world, with a particular focus on the Arab population.

The transition from genomics to phenomics in personalized population health

Modern health care faces several serious challenges, including an ageing population and its inherent burden of chronic diseases, rising costs and marginal quality metrics. By assessing and optimizing the health trajectory of each individual using a data-driven personalized approach that reflects their genetics, behaviour and environment, we can start to address these challenges. This assessment includes longitudinal phenome measures, such as the blood proteome and metabolome, gut microbiome composition and function, and lifestyle and behaviour through wearables and questionnaires. Here, we review ongoing large-scale genomics and longitudinal phenomics efforts and the powerful insights they provide into wellness. We describe our vision for the transformation of the current health care from disease-oriented to data-driven, wellness-oriented and personalized population health.

Autosomal recessive IL12RB1 mutation: A case report of a Sudanese child and his father

Introduction

Mendelian susceptibility to mycobacterial disease (MSMD) is a rare inherited condition characterized by selective susceptibility to weakly virulent mycobacteria, such as substrains of the bacille Calmette–Guérin (BCG) vaccine and different environmental mycobacteria.

Case presentation

A 7-year-old Sudanese boy was referred to the immunology clinic with a suspected diagnosis of MSMD. This followed multiple presentations with disseminated tuberculosis and typhoid fever. Genetic testing surprisingly revealed pathogenic homozygous variants in IL12RB1 Exon 9, c.913A>T (p. Lys305*) in both the patient and his father, with a completely healthy asymptomatic carrier mother who is not blood related to the patient’s father.

Conclusion

It is challenging to diagnose MSMD, especially in developing countries where health systems are poor and have limited resources. Family history and genetic tests may help in early MSMD treatment and avoiding disease complications.

PMEL is mutated in oculocutaneous albinism

Oculocutaneous albinism (OCA) is a group of Mendelian disorders characterized by hypopigmentation of skin, hair and pigmented ocular structures. While much of the genetic heterogeneity of OCA has been resolved, many patients still lack a molecular diagnosis following exome sequencing. Here, we report a consanguineous family in which the index patient presented with OCA and Hirschsprung disease but tested negative for known genetic causes of OCA. Instead, he was found to have a homozygous presumptive loss of function variant in PMEL. PMEL encodes a scaffolding protein that is essential for the normal maturation of melanosomes and normal deposition of the melanin pigment therein. Numerous PMEL vertebrate ortholog mutants have been reported and all were characterized by conspicuous pigmentary abnormalities. We suggest that the patient we report is the first human equivalent of PMEL loss of function.

Case Report: A Novel Homozygous Missense Variant of FBN3 Supporting It Is a New Candidate Gene Causative of a Bardet–Biedl Syndrome–Like Phenotype

Fibrillin proteins are extracellular matrix glycoproteins assembling into microfibrils. FBN1, FBN2, and FBN3 encode the human fibrillins and mutations in FBN1 and FBN2 cause connective tissue disorders called fibrillinopathies, affecting cardiovascular, dermal, skeletal, and ocular tissues. Recently, mutations of the less characterized fibrillin family member, FBN3, have been associated in a single family with Bardet–Biedl syndrome (BBS). Here, we report on a patient born from two first cousins and affected by developmental delay, cognitive impairment, obesity, dental and genital anomalies, and brachydactyly/syndactyly. His phenotype was very similar to that reported in the previous FBN3-mutated family and fulfilled BBS clinical diagnostic criteria, although lacking polydactyly, the most recurrent clinical feature, as the previous siblings described. A familial SNP-array and proband’s WES were performed prioritizing candidate variants on the sole patient’s runs of homozygosity. This analysis disclosed a novel homozygous missense variant in FBN3 (NM_032447:c.5434A>G; NP_115823:p.Ile1812Val; rs115948457), inherited from the heterozygous parents. This study further supports that FBN3 is a candidate gene for a BBS-like syndrome characterized by developmental delay, cognitive impairment, obesity, dental, genital, and skeletal anomalies. Anyway, additional studies are necessary to investigate the exact role of the gene and possible interactions between FBN3 and BBS proteins.

Clinical, genetic, and functional characterization of the glycine receptor β-subunit A455P variant in a family affected by hyperekplexia syndrome

Hyperekplexia is a rare neurological disorder characterized by exaggerated startle responses affecting newborns with the hallmark characteristics of hypertonia, apnea, and noise or touch-induced nonepileptic seizures. The genetic causes of the disease can vary, and several associated genes and mutations have been reported to affect glycine receptors (GlyRs); however, the mechanistic links between GlyRs and hyperekplexia are not yet understood. Here, we describe a patient with hyperekplexia from a consanguineous family. Extensive genetic screening using exome sequencing coupled with autozygome analysis and iterative filtering supplemented by in silico prediction identified that the patient carries the homozygous missense mutation A455P in GLRB, which encodes the GlyR β-subunit. To unravel the physiological and molecular effects of A455P on GlyRs, we used electrophysiology in a heterologous system as well as immunocytochemistry, confocal microscopy, and cellular biochemistry. We found a reduction in glycine-evoked currents in N2A cells expressing the mutation compared to WT cells. Western blot analysis also revealed a reduced amount of GlyR β protein both in cell lysates and isolated membrane fractions. In line with the above observations, coimmunoprecipitation assays suggested that the GlyR α1-subunit retained coassembly with βA455P to form membrane-bound heteromeric receptors. Finally, structural modeling showed that the A455P mutation affected the interaction between the GlyR β-subunit transmembrane domain 4 and the other helices of the subunit. Taken together, our study identifies and validates a novel loss-of-function mutation in GlyRs whose pathogenicity is likely to cause hyperekplexia in the affected individual.

Mitochondrial "dysmorphology" in variant classification

Mitochondrial disorders are challenging to diagnose. Exome sequencing has greatly enhanced the diagnostic precision of these disorders although interpreting variants of uncertain significance (VUS) remains a formidable obstacle. Whether specific mitochondrial morphological changes can aid in the classification of these variants is unknown. Here, we describe two families (four patients), each with a VUS in a gene known to affect the morphology of mitochondria through a specific role in the fission-fusion balance. In the first, the missense variant in MFF, encoding a fission factor, was associated with impaired fission giving rise to a characteristically over-tubular appearance of mitochondria. In the second, the missense variant in DNAJA3, which has no listed OMIM phenotype, was associated with fragmented appearance of mitochondria consistent with its published deficiency states. In both instances, the highly specific phenotypes allowed us to upgrade the classification of the variants. Our results suggest that, in select cases, mitochondrial "dysmorphology" can be helpful in interpreting variants to reach a molecular diagnosis.

ASTL is mutated in female infertility

Female infertility is a relatively common phenotype with a growing number of single gene causes although these account for only a minority of cases. Here, we report a consanguineous family in which adult females who are homozygous for a truncating variant in ASTL display markedly reduced fertility in a pattern strikingly similar to Astl^-/- female mice. ASTL encodes ovastacin, which is known to trigger zona pellucida hardening (ZPH) as part of the cortical reaction upon fertilization. ZPH is required for normal early embryonic development and its absence can be caused by pathogenic variants in other zona pellucida proteins that result in a similar infertility phenotype in humans and mouse. This is the first report of ASTL-related infertility in humans and suggests that the inclusion of ASTL in female infertility gene panels is warranted.

Whole exome sequencing in 17 consanguineous Iranian pedigrees expands the mutational spectrum of inherited retinal dystrophies

Inherited retinal dystrophies (IRDs) constitute one of the most heterogeneous groups of Mendelian human disorders. Using autozygome-guided next-generation sequencing methods in 17 consanguineous pedigrees of Iranian descent with isolated or syndromic IRD, we identified 17 distinct genomic variants in 11 previously-reported disease genes. Consistent with a recessive inheritance pattern, as suggested by pedigrees, variants discovered in our study were exclusively bi-allelic and mostly in a homozygous state (in 15 families out of 17, or 88%). Out of the 17 variants identified, 5 (29%) were never reported before. Interestingly, two mutations (GUCY2D:c.564dup, p.Ala189ArgfsTer130 and TULP1:c.1199G > A, p.Arg400Gln) were also identified in four separate pedigrees (two pedigrees each). In addition to expanding the mutational spectrum of IRDs, our findings confirm that the traditional practice of endogamy in the Iranian population is a prime cause for the appearance of IRDs.

PLXNA2 as a candidate gene in patients with intellectual disability

Intellectual disability (ID) is one of the most common disabilities in humans. In an effort to contribute to the expanding genetic landscape of ID, we describe a novel autosomal recessive ID candidate gene. Combined autozygome/exome analysis was performed in two unrelated consanguineous families with ID. Each of the two families had a novel homozygous likely deleterious variant in PLXNA2 and displayed the core phenotype of ID. PLXNA2 belongs to a family of transmembrane proteins that function as semaphorin receptors. Sema5A-PlexinA2 is known to regulate brain development in mouse, and Plxna2-/- mice display defective associative learning, sociability, and sensorimotor gating. We note the existence of variability in the phenotype among the three patients, including the existence of variable degree of ID, ranging from borderline intellectual functioning to moderate-severe ID, and the presence of cardiac anomalies in only one of the patients. We propose incomplete penetrance as a possible explanation of the observed difference in phenotypes. Future cases will be needed to support the proposed link between PLXNA2 and ID in humans.

An Overview of Molecular Genetic Diagnosis Techniques

Our understanding of genetic disease(s) has increased exponentially since the completion of human genome sequencing and the development of numerous techniques to detect genetic variants. These techniques have not only allowed us to diagnose genetic disease, but in so doing, also provide increased understanding of the pathogenesis of these diseases to aid in developing appropriate therapeutic options. Additionally, the advent of next-generation or massively parallel sequencing (NGS/MPS) is increasingly being used in the clinical setting, as it can detect a number of abnormalities from point mutations to chromosomal rearrangements as well as aberrations within the transcriptome. In this article, we will discuss the use of multiple techniques that are used in genetic diagnosis. © 2020 by John Wiley & Sons, Inc.

Aarskog-Scott syndrome: phenotypic and genetic heterogeneity

Aarskog-Scott syndrome (AAS) is a rare developmental disorder which primarily affects males and has a relative prevalence of 1 in 25,000 in the general population. AAS patients usually present with developmental complications including short stature and facial, skeletal and urogenital anomalies. The spectrum of genotype-phenotype correlations in AAS is unclear and mutations of the FGD1 gene on the proximal short arm of chromosome X account for only 20% of the incidence of the disorder. Failure to identify pathogenic variants in patients referred for FGD1 screening suggests heterogeneity underlying pathophysiology of the condition. Furthermore, overlapping features of AAS with several other developmental disorders increase the complexity of diagnosis. Cytoskeletal signaling may be involved in the pathophysiology of AAS. The FGD1 protein family has a role in activation of CDC42 (Cell Division Control protein 42 homolog) which has a core function in remodeling of extracellular matrix and the transcriptional activation of many modulators of development. Therefore, mutations in components in the EGFR1 (Epidermal Growth Factor Receptor 1) signaling pathway, to which CDC42 belongs, may contribute to pathophysiology. Parallel sequencing strategies (so-called next generation sequencing or high throughput sequencing) enables simultaneous production of millions of sequencing reads that enormously facilitate cost-effective identification of cryptic mutations in heterogeneous monogenic disorders. Here we review the source of phenotypic and genetic heterogeneity in the context of AAS and discuss the applicability of next generation sequencing for identification of novel mutations underlying AAS.

Mapping the genetic basis of diabetes mellitus in the Australian Burmese cat (Felis catus)

Diabetes mellitus, a common endocrinopathy affecting domestic cats, shares many clinical and pathologic features with type 2 diabetes in humans. In Australia and Europe, diabetes mellitus is almost four times more common among Burmese cats than in other breeds. As a genetically isolated population, the diabetic Australian Burmese cat provides a spontaneous genetic model for studying diabetes mellitus in humans. Studying complex diseases in pedigreed breeds facilitates tighter control of confounding factors including population stratification, allelic frequencies and environmental heterogeneity. We used the feline SNV array and whole genome sequence data to undertake a genome wide-association study and runs of homozygosity analysis, of a case–control cohort of Australian and European Burmese cats. Our results identified diabetes-associated haplotypes across chromosomes A3, B1 and E1 and selective sweeps across the Burmese breed on chromosomes B1, B3, D1 and D4. The locus on chromosome B1, common to both analyses, revealed coding and splice region variants in candidate genes, ANK1, EPHX2 and LOX2, implicated in diabetes mellitus and lipid dysregulation. Mapping this condition in Burmese cats has revealed a polygenic spectrum, implicating loci linked to pancreatic beta cell dysfunction, lipid dysregulation and insulin resistance in the pathogenesis of diabetes mellitus in the Burmese cat.

Homozygosity mapping and next generation sequencing for the genetic diagnosis of hereditary ataxia and spastic paraplegia in consanguineous families

INTRODUCTION

Genetic inheritance plays key roles in patients with ataxia and/or spastic paraplegia in consanguineous families. This study aims to clarify the genetic spectrum of patients with autosomal recessive hereditary ataxia and spastic paraplegias (AR-HA/HSPs) in consanguineous families.

METHODS

A total of 36 AR-HA/HSPs consanguineous pedigrees from China were recruited into this study. Next generation sequencing (NGS), guided by homozygosity mapping (HM), was applied to identify the pathogenic variants in known genes or novel candidate genes.

RESULTS

We totally made molecular diagnosis in 47.2% (17/36) of AR-HA/HSPs families. Among them, 13 AR-HAs carried pathogenic variants in SETX (n = 4), SACS (n = 2), STUB1, HSD17B4, NEU1, ADCK3, TPP1, PLA2G6 and MTCL1, while four AR-HSPs carried pathogenic variants in SPG11, ZFYVE26, ATP13A2 and ABCD1. One homozygous nonsense mutation in MRPS27 was identified in an AR-HA family, which was potentially a novel candidate gene of AR-HA.

CONCLUSION

HM and NGS can serve as an efficient molecular diagnostic tool for AR-HA/HSPs in consanguineous families. Our findings provide a better understanding of genetic architecture of AR-HA/HSPs in consanguinity and broaden the clinical-genetic spectrum of the disease.

Analysis of transcript-deleterious variants in Mendelian disorders: implications for RNA-based diagnostics

Background

At least 50% of patients with suspected Mendelian disorders remain undiagnosed after whole-exome sequencing (WES), and the extent to which non-coding variants that are not captured by WES contribute to this fraction is unclear. Whole transcriptome sequencing is a promising supplement to WES, although empirical data on the contribution of RNA analysis to the diagnosis of Mendelian diseases on a large scale are scarce.

Results

Here, we describe our experience with transcript-deleterious variants (TDVs) based on a cohort of 5647 families with suspected Mendelian diseases. We first interrogate all families for which the respective Mendelian phenotype could be mapped to a single locus to obtain an unbiased estimate of the contribution of TDVs at 18.9%. We examine the entire cohort and find that TDVs account for 15% of all “solved” cases. We compare the results of RT-PCR to in silico prediction. Definitive results from RT-PCR are obtained from blood-derived RNA for the overwhelming majority of variants (84.1%), and only a small minority (2.6%) fail analysis on all available RNA sources (blood-, skin fibroblast-, and urine renal epithelial cells-derived), which has important implications for the clinical application of RNA-seq. We also show that RNA analysis can establish the diagnosis in 13.5% of 155 patients who had received “negative” clinical WES reports. Finally, our data suggest a role for TDVs in modulating penetrance even in otherwise highly penetrant Mendelian disorders.

Conclusions

Our results provide much needed empirical data for the impending implementation of diagnostic RNA-seq in conjunction with genome sequencing.

Biallelic Mutations in Tetratricopeptide Repeat Domain 26 (Intraflagellar Transport 56) Cause Severe Biliary Ciliopathy in Humans

BACKGROUND AND AIMS

The clinical consequences of defective primary cilium (ciliopathies) are characterized by marked phenotypic and genetic heterogeneity. Although fibrocystic liver disease is an established ciliopathy phenotype, severe neonatal cholestasis is rarely recognized as such.

APPROACH AND RESULTS

We describe seven individuals from seven families with syndromic ciliopathy clinical features, including severe neonatal cholestasis (lethal in one and necessitating liver transplant in two). Positional mapping revealed a single critical locus on chromosome 7. Whole-exome sequencing revealed three different homozygous variants in Tetratricopeptide Repeat Domain 26 (TTC26) that fully segregated with the phenotype. TTC26 (intraflagellar transport [IFT] 56/DYF13) is an atypical component of IFT-B complex, and deficiency of its highly conserved orthologs has been consistently shown to cause defective ciliary function in several model organisms. We show that cilia in TTC26-mutated patient cells display variable length and impaired function, as indicated by dysregulated sonic hedgehog signaling, abnormal staining for IFT-B components, and transcriptomic clustering with cells derived from individuals with closely related ciliopathies. We also demonstrate a strong expression of Ttc26 in the embryonic mouse liver in a pattern consistent with its proposed role in the normal development of the intrahepatic biliary system.

CONCLUSIONS

In addition to establishing a TTC26-related ciliopathy phenotype in humans, our results highlight the importance of considering ciliopathies in the differential diagnosis of severe neonatal cholestasis even in the absence of more typical features.

Bi-allelic Variants in RALGAPA1 Cause Profound Neurodevelopmental Disability, Muscular Hypotonia, Infantile Spasms, and Feeding Abnormalities

Ral (Ras-like) GTPases play an important role in the control of cell migration and have been implicated in Ras-mediated tumorigenicity. Recently, variants in RALA were also described as a cause of intellectual disability and developmental delay, indicating the relevance of this pathway to neuropediatric diseases. Here, we report the identification of bi-allelic variants in RALGAPA1 (encoding Ral GTPase activating protein catalytic alpha subunit 1) in four unrelated individuals with profound neurodevelopmental disability, muscular hypotonia, feeding abnormalities, recurrent fever episodes, and infantile spasms . Dysplasia of corpus callosum with focal thinning of the posterior part and characteristic facial features appeared to be unifying findings. RalGAPA1 was absent in the fibroblasts derived from two affected individuals suggesting a loss-of-function effect of the RALGAPA1 variants. Consequently, RalA activity was increased in these cell lines, which is in keeping with the idea that RalGAPA1 deficiency causes a constitutive activation of RalA. Additionally, levels of RalGAPB, a scaffolding subunit of the RalGAP complex, were dramatically reduced, indicating a dysfunctional RalGAP complex. Moreover, RalGAPA1 deficiency clearly increased cell-surface levels of lipid raft components in detached fibroblasts, which might indicate that anchorage-dependence of cell growth signaling is disturbed. Our findings indicate that the dysregulation of the RalA pathway has an important impact on neuronal function and brain development. In light of the partially overlapping phenotype between RALA- and RALGAPA1-associated diseases, it appears likely that dysregulation of the RalA signaling pathway leads to a distinct group of genetic syndromes that we suggest could be named RALopathies.

Confirming the recessive inheritance of PERP-related erythrokeratoderma

Erythrokeratoderma (EK) is heterogeneous clinical entity characterized by excessive scaling with resulting erythrokeratotic plaques. Several genes have been linked to EK and they encode a number of proteins that are important for the integrity of the keratinocyte layer of the epidermis. PERP is a transcription factor that is activated by both p53 and p63. However, its deficiency in a mouse model appears to only recapitulate p63-mediated role in skin development and organization. We report an extended multiplex consanguineous family in which an EK phenotype with a striking similarity to that observed in Perp^-/- mice, is mapped to an autozygous region on chromosome 6 that spans PERP. Whole-exome sequencing revealed a novel variant in PERP that fully segregated with the phenotype. Functional analysis of patient- and control-derived keratinocytes revealed a deleterious effect of the identified variant on the intracellular localization of PERP. A previous report showed that PERP mutation causes a dominant form of keratoderma but a single patient in that report with a homozygous variant in PERP suggests that recessive inheritance is also possible. Our results, therefore, support the establishment of an autosomal recessive PERP-related EK phenotype in humans.

AUDACITY: A comprehensive approach for the detection and classification of Runs of Homozygosity in medical and population genomics

Runs of Homozygosity (RoHs) are popular among geneticists as the footprint of demographic processes, evolutionary forces and inbreeding in shaping our genome, and are known to confer risk of Mendelian and complex diseases. Notwithstanding growing interest in their study, there is unmet need for reliable and rapid methods for genomic analyses in large data sets. AUDACITY is a tool integrating novel RoH detection algorithm and autozygosity prediction score for prioritization of mutation-surrounding regions. It processes data in VCF file format, and outperforms existing methods in identifying RoHs of any size. Simulations and analysis of real exomes/genomes show its potential to foster future RoH studies in medical and population genomics.

Multiple Family Members With Delayed Cord Separtion and Combined Immunodeficiency With Novel Mutation in IKBKB

Background: Inhibitor of kappa kinase 2 (IKK2) deficiency is a recently described combined immunodeficiency. It undermines the nuclear factor-kappa B (NF-κB) activation pathway. Methods: The clinical and immunological data of four patients diagnosed with combined immunodeficiency (CID) from two related Saudi families were collected. Autozygosity mapping of all available members and whole exome sequencing of the index case were performed to define the genetic etiology. Results: The patients had early onset (2-4 months of age) severe infections caused by viruses, bacteria, mycobacteria, and fungi. They all had hypogammaglobulinemia and low absolute lymphocyte count. Their lymphocytes failed to respond to PHA mitogen stimulation. A novel homozygous non-sense mutation in the IKBKB gene, c.850C>T (p. Arg284*) was identified in the index patient and segregated with the disease in the rest of the family. He underwent hematopoietic stem cell transplantation (HSCT) from a fully matched sibling with no conditioning. The other three patients succumbed to their disease. Interestingly, all patients had delayed umbilical cord separation. Conclusion: IKK2 deficiency causes CID with high mortality. Immune reconstitution with HSCT should be considered as early as possible. Delayed umbilical cord separation in CID patients may be a clue to IKK2 deficiency.

Exploring the Genetic Landscape of Retinal Diseases in North-Western Pakistan Reveals a High Degree of Autozygosity and a Prevalent Founder Mutation in ABCA4

Variants in more than 271 different genes have been linked to hereditary retinal diseases, making comprehensive genomic approaches mandatory for accurate diagnosis. We explored the genetic landscape of retinal disorders in consanguineous families from North-Western Pakistan, harboring a population of approximately 35 million inhabitants that remains relatively isolated and highly inbred (~50% consanguinity). We leveraged on the high degree of consanguinity by applying genome-wide high-density single-nucleotide polymorphism (SNP) genotyping followed by targeted Sanger sequencing of candidate gene(s) lying inside autozygous intervals. In addition, we performed whole-exome sequencing (WES) on at least one proband per family. We identified 7 known and 4 novel variants in a total of 10 genes (ABCA4, BBS2, CNGA1, CNGA3, CNGB3, MKKS, NMNAT1, PDE6B, RPE65, and TULP1) previously known to cause inherited retinal diseases. In spite of all families being consanguineous, compound heterozygosity was detected in one family. All homozygous pathogenic variants resided in autozygous intervals ≥2.0 Mb in size. Putative founder variants were observed in the ABCA4 (NM_000350.2:c.214G>A; p.Gly72Arg; ten families) and NMNAT1 genes (NM_022787.3:c.25G>A; p.Val9Met; two families). We conclude that geographic isolation and sociocultural tradition of intrafamilial mating in North-Western Pakistan favor both the clinical manifestation of rare “generic” variants and the prevalence of founder mutations.

MDH1 deficiency is a metabolic disorder of the malate-aspartate shuttle associated with early onset severe encephalopathy

The reversible oxidation of L-malate to oxaloacetate is catalyzed by NAD(H)-dependent malate dehydrogenase (MDH). MDH plays essential roles in the malate-aspartate shuttle and the tricarboxylic acid cycle. These metabolic processes are important in mitochondrial NADH supply for oxidative phosphorylation. Recently, bi-allelic mutations in mitochondrial MDH2 were identified in patients with global developmental delay, epilepsy and lactic acidosis. We now report two patients from an extended consanguineous family with a deleterious variant in the cytosolic isoenzyme of MDH (MDH1). The homozygous missense variant in the NAD⁺-binding domain of MDH1 led to severely diminished MDH protein expression. The patients presented with global developmental delay, epilepsy and progressive microcephaly. Both patients had normal concentrations of plasma amino acids, acylcarnitines, lactate, and urine organic acids. To identify the metabolic consequences of MDH1 deficiency, untargeted metabolomics was performed on dried blood spots (DBS) from the patients and in MDH1 knockout HEK293 cells that were generated by Crispr/Cas9. Increased levels of glutamate and glycerol-3-phosphate were found in DBS of both patients. In MDH1 KO HEK293 cells, increased levels of glycerol-3-phosphate were also observed, as well as increased levels of aspartate and decreased levels of fumarate. The consistent finding of increased concentrations of glycerol-3-phosphate may represent a compensatory mechanism to enhance cytosolic oxidation of NADH by the glycerol-P-shuttle. In conclusion, MDH1 deficiency is a new metabolic defect in the malate-aspartate shuttle characterized by a severe neurodevelopmental phenotype with elevated concentrations of glycerol-3-phosphate as a potential biomarker.

Lessons Learned from Large-Scale, First-Tier Clinical Exome Sequencing in a Highly Consanguineous Population

We report the results of clinical exome sequencing (CES) on >2,200 previously unpublished Saudi families as a first-tier test. The predominance of autosomal-recessive causes allowed us to make several key observations. We highlight 155 genes that we propose to be recessive, disease-related candidates. We report additional mutational events in 64 previously reported candidates (40 recessive), and these events support their candidacy. We report recessive forms of genes that were previously associated only with dominant disorders and that have phenotypes ranging from consistent with to conspicuously distinct from the known dominant phenotypes. We also report homozygous loss-of-function events that can inform the genetics of complex diseases. We were also able to deduce the likely causal variant in most couples who presented after the loss of one or more children, but we lack samples from those children. Although a similar pattern of mostly recessive causes was observed in the prenatal setting, the higher proportion of loss-of-function events in these cases was notable. The allelic series presented by the wealth of recessive variants greatly expanded the phenotypic expression of the respective genes. We also make important observations about dominant disorders; these observations include the pattern of de novo variants, the identification of 74 candidate dominant, disease-related genes, and the potential confirmation of 21 previously reported candidates. Finally, we describe the influence of a predominantly autosomal-recessive landscape on the clinical utility of rapid sequencing (Flash Exome). Our cohort's genotypic and phenotypic data represent a unique resource that can contribute to improved variant interpretation through data sharing.

Bi-allelic Mutations in FAM149B1 Cause Abnormal Primary Cilium and a Range of Ciliopathy Phenotypes in Humans

Ciliopathies are clinical disorders of the primary cilium with widely recognized phenotypic and genetic heterogeneity. In two Arab consanguineous families, we mapped a ciliopathy phenotype that most closely matches Joubert syndrome (hypotonia, developmental delay, typical facies, oculomotor apraxia, polydactyly, and subtle posterior fossa abnormalities) to a single locus in which a founder homozygous truncating variant in FAM149B1 was identified by exome sequencing. We subsequently identified a third Arab consanguineous multiplex family in which the phenotype of Joubert syndrome/oral-facial-digital syndrome (OFD VI) was found to co-segregate with the same founder variant in FAM149B1. Independently, autozygosity mapping and exome sequencing in a consanguineous Turkish family with Joubert syndrome highlighted a different homozygous truncating variant in the same gene. FAM149B1 encodes a protein of unknown function. Mutant fibroblasts were found to have normal ciliogenesis potential. However, distinct cilia-related abnormalities were observed in these cells: abnormal accumulation IFT complex at the distal tips of the cilia, which assumed bulbous appearance, increased length of the primary cilium, and dysregulated SHH signaling. We conclude that FAM149B1 is required for normal ciliary biology and that its deficiency results in a range of ciliopathy phenotypes in humans along the spectrum of Joubert syndrome.

PUS7 mutations impair pseudouridylation in humans and cause intellectual disability and microcephaly

Pseudouridylation is the most common post-transcriptional modification, wherein uridine is isomerized into 5-ribosyluracil (pseudouridine, Ψ). The resulting increase in base stacking and creation of additional hydrogen bonds are thought to enhance RNA stability. Pseudouridine synthases are encoded in humans by 13 genes, two of which are linked to Mendelian diseases: PUS1 and PUS3. Very recently, PUS7 mutations were reported to cause intellectual disability with growth retardation. We describe two families in which two different homozygous PUS7 mutations (missense and frameshift deletion) segregate with a phenotype comprising intellectual disability and progressive microcephaly. Short stature and hearing loss were variable in these patients. Functional characterization of the two mutations confirmed that both result in decreased levels of Ψ₁₃ in tRNAs. Furthermore, the missense variant of the S. cerevisiae ortholog failed to complement the growth defect of S. cerevisiae pus7Δ trm8Δ mutants. Our results confirm that PUS7 is a bona fide Mendelian disease gene and expand the list of human diseases caused by impaired pseudouridylation.

NUP214 deficiency causes severe encephalopathy and microcephaly in humans

Nuclear pore complex (NPC) is a fundamental component of the nuclear envelope and is key to the nucleocytoplasmic transport. Mutations in several NUP genes that encode individual components of NPC known as nucleoporins have been identified in recent years among patients with static encephalopathies characterized by developmental delay and microcephaly. We describe a multiplex consanguineous family in which four affected members presented with severe neonatal hypotonia, profound global developmental delay, progressive microcephaly and early death. Autozygome and linkage analysis revealed that this phenotype is linked to a founder disease haplotype (chr9:127,113,732-135,288,807) in which whole exome sequencing revealed the presence of a novel homozygous missense variant in NUP214. Functional analysis of patient-derived fibroblasts recapitulated the dysmorphic phenotype of nuclei that was previously described in NUP214 knockdown cells. In addition, the typical rim staining of NUP214 is largely displaced, further supporting the deleterious effect of the variant. Our data expand the list of NUP genes that are mutated in encephalopathy disorders in humans.

Identification and Analysis of Genes Associated with Inherited Retinal Diseases

Inherited retinal diseases (IRDs) display a very high degree of clinical and genetic heterogeneity, which poses challenges in finding the underlying defects in known IRD-associated genes and in identifying novel IRD-associated genes. Knowledge on the molecular and clinical aspects of IRDs has increased tremendously in the last decade. Here, we outline the state-of-the-art techniques to find the causative genetic variants, with special attention for next-generation sequencing which can combine molecular diagnostics and retinal disease gene identification. An important aspect is the functional assessment of rare variants with RNA and protein effects which can only be predicted in silico. We therefore describe the in vitro assessment of putative splice defects in human embryonic kidney cells. In addition, we outline the use of stem cell technology to generate photoreceptor precursor cells from patients' somatic cells which can subsequently be used for RNA and protein studies. Finally, we outline the in silico methods to interpret the causality of variants associated with inherited retinal disease and the registry of these variants.

The many faces of peroxisomal disorders: Lessons from a large Arab cohort

Defects in the peroxisomes biogenesis and/or function result in peroxisomal disorders. In this study, we describe the largest Arab cohort to date (72 families) of clinically, biochemically and molecularly characterized patients with peroxisomal disorders. At the molecular level, we identified 43 disease-causing variants, half of which are novel. The founder nature of many of the variants allowed us to calculate the minimum disease burden for these disorders in our population ~1:30 000, which is much higher than previous estimates in other populations. Clinically, we found an interesting trend toward genotype/phenotype correlation in terms of long-term survival. Nearly half (40/75) of our peroxisomal disorders patients had documented survival beyond 1 year of age. Most unusual among the long-term survivors was a multiplex family in which the affected members presented as adults with non-specific intellectual disability and epilepsy. Other unusual presentations included the very recently described peroxisomal fatty acyl-CoA reductase 1 disorder as well as CRD, spastic paraparesis, white matter (CRSPW) syndrome. We conclude that peroxisomal disorders are highly heterogeneous in their clinical presentation. Our data also confirm the demonstration that milder forms of Zellweger spectrum disorders cannot be ruled out by the "gold standard" very long chain fatty acids assay, which highlights the value of a genomics-first approach in these cases.

Identification of novel loci for pediatric cholestatic liver disease defined by KIF12, PPM1F, USP53, LSR, and WDR83OS pathogenic variants

PURPOSE

Genetic testing in pediatric cholestasis can be very informative but genetic causes have not been fully characterized.

METHODS

Exome sequencing and positional mapping in seven families with cholestatic liver disease and negative clinical testing for known disease genes.

RESULTS

KIF12, which encodes a microtubule motor protein with a tentative role in cell polarity, was found to harbor three homozygous likely deleterious variants in three families with sclerosing cholangitis. KIF12 expression is dependent on HNF-1β, deficiency which is known to cause bile duct dysmorphogenesis associated with loss of KIF12 expression. In another extended family, we mapped an apparently novel syndrome of sclerosing cholangitis, short stature, hypothyroidism, and abnormal tongue pigmentation in two cousins to a homozygous variant in PPM1F (POPX2), a regulator of kinesin-mediated ciliary transport. In the fifth family, a syndrome of normal gamma glutamyltransferase (GGT) cholestasis and hearing loss was found to segregate with a homozygous truncating variant in USP53, which encodes an interactor with TJP2. In the sixth family, we mapped a novel syndrome of transient neonatal cholestasis, intellectual disability, and short stature to a homozygous variant in LSR, an important regulator of liver development. In the last family of three affected siblings, a novel syndrome of intractable itching, hypercholanemia, short stature, and intellectual disability was mapped to a single locus that contains a homozygous truncating variant in WDR83OS (C19orf56), known to interact with ATP13A2 and BSEP.

CONCLUSION

Our results expand the genetic heterogeneity of pediatric cholestatic liver disease and highlight the vulnerability of bile homeostasis to a wide range of molecular perturbations.

Biallelic UFM1 and UFC1 mutations expand the essential role of ufmylation in brain development

The post-translational modification of proteins through the addition of UFM1, also known as ufmylation, plays a critical developmental role as revealed by studies in animal models. The recent finding that biallelic mutations in UBA5 (the E1-like enzyme for ufmylation) cause severe early-onset encephalopathy with progressive microcephaly implicates ufmylation in human brain development. More recently, a homozygous UFM1 variant was proposed as a candidate aetiology of severe early-onset encephalopathy with progressive microcephaly. Here, we establish a locus for severe early-onset encephalopathy with progressive microcephaly based on two families, and map the phenotype to a novel homozygous UFM1 mutation. This mutation has a significantly diminished capacity to form thioester intermediates with UBA5 and with UFC1 (the E2-like enzyme for ufmylation), with resulting impaired ufmylation of cellular proteins. Remarkably, in four additional families where eight children have severe early-onset encephalopathy with progressive microcephaly, we identified two biallelic UFC1 mutations, which impair UFM1-UFC1 intermediate formation with resulting widespread reduction of cellular ufmylation, a pattern similar to that observed with UFM1 mutation. The striking resemblance between UFM1- and UFC1-related clinical phenotype and biochemical derangements strongly argues for an essential role for ufmylation in human brain development. The hypomorphic nature of UFM1 and UFC1 mutations and the conspicuous depletion of biallelic null mutations in the components of this pathway in human genome databases suggest that it is necessary for embryonic survival, which is consistent with the embryonic lethal nature of knockout models for the orthologous genes.

Expanding the phenome and variome of skeletal dysplasia

PURPOSE

To describe our experience with a large cohort (411 patients from 288 families) of various forms of skeletal dysplasia who were molecularly characterized.

METHODS

Detailed phenotyping and next-generation sequencing (panel and exome).

RESULTS

Our analysis revealed 224 pathogenic/likely pathogenic variants (54 (24%) of which are novel) in 123 genes with established or tentative links to skeletal dysplasia. In addition, we propose 5 genes as candidate disease genes with suggestive biological links (WNT3A, SUCO, RIN1, DIP2C, and PAN2). Phenotypically, we note that our cohort spans 36 established phenotypic categories by the International Skeletal Dysplasia Nosology, as well as 18 novel skeletal dysplasia phenotypes that could not be classified under these categories, e.g., the novel C3orf17-related skeletal dysplasia. We also describe novel phenotypic aspects of well-known disease genes, e.g., PGAP3-related Toriello-Carey syndrome-like phenotype. We note a strong founder effect for many genes in our cohort, which allowed us to calculate a minimum disease burden for the autosomal recessive forms of skeletal dysplasia in our population (7.16E-04), which is much higher than the global average.

CONCLUSION

By expanding the phenotypic, allelic, and locus heterogeneity of skeletal dysplasia in humans, we hope our study will improve the diagnostic rate of patients with these conditions.

Genotype-phenotype correlations, dystonia and disease progression in spinocerebellar ataxia type 14

Background: Spinocerebellar ataxia type 14 is a rare form of autosomal dominant cerebellar ataxia caused by mutations in protein kinase Cγ gene. Clinically, it presents with a slowly progressive, mainly pure cerebellar ataxia.

Methods: Using next generation sequencing, we screened 194 families with autosomal dominant cerebellar ataxia and normal polyglutamine repeats. In‐depth phenotyping was performed using validated clinical rating scales neuroimaging and electrophysiological investigations.

Results: We identified 25 individuals from 13 families carrying pathogenic mutations in protein kinase Cγ gene. A total of 10 unique protein kinase Cγ gene mutations have been confirmed of which 5 are novel and 5 were previously described. Our data suggest that the age at onset is highly variable; disease course is slowly progressive and rarely associated with severe disability. However, one third of patients presented with a complex ataxia comprising severe focal and/or task‐induced dystonia, peripheral neuropathy, parkinsonism, myoclonus, and pyramidal syndrome. The most complex phenotype is related to a missense mutation in the catalytic domain in exon 11.

Conclusion: We present one of the largest genetically confirmed spinocerebellar ataxia type 14 cohorts contributing novel variants and clinical characterisation. We show that although protein kinase Cγ gene mutations present mainly as slowly progressive pure ataxia, more than a third of cases had a complex phenotype. Overall, our case series extends the phenotype and suggests that protein kinase Cγ gene mutations should be considered in patients with slowly progressive autosomal dominant cerebellar ataxia, particularly when myoclonus, dystonia, or mild cognitive impairment are present in the absence of polyglutamine expansion. © 2018 The Authors. Movement Disorders published by Wiley Periodicals, Inc. on behalf of International Parkinson and Movement Disorder Society.

Genetic investigation of 93 families with microphthalmia or posterior microphthalmos

Microphthalmia is a developmental eye defect that is highly variable in severity and in its potential for systemic association. Despite the discovery of many disease genes in microphthalmia, at least 50% of patients remain undiagnosed genetically. Here, we describe a cohort of 147 patients (93 families) from our highly consanguineous population with various forms of microphthalmia (including the distinct entity of posterior microphthalmos) that were investigated using a next-generation sequencing multi-gene panel (i-panel) as well as whole exome sequencing and molecular karyotyping. A potentially causal mutation was identified in the majority of the cohort with microphthalmia (61%) and posterior microphthalmos (82%). The identified mutations (55 point mutations, 15 of which are novel) spanned 24 known disease genes, some of which have not or only very rarely been linked to microphthalmia (PAX6, SLC18A2, DSC3 and CNKSR1). Our study has also identified interesting candidate variants in 2 genes that have not been linked to human diseases (MYO10 and ZNF219), which we present here as novel candidates for microphthalmia. In addition to revealing novel phenotypic aspects of microphthalmia, this study expands its allelic and locus heterogeneity and highlights the need for expanded testing of patients with this condition.

Genetics of Obesity in Consanguineous Populations: Toward Precision Medicine and the Discovery of Novel Obesity Genes

OBJECTIVE

Consanguinity has been instrumental in the elucidation of many Mendelian genetic diseases. Here, the unique advantage of consanguineous populations was considered in the quest for genes causing obesity.

METHODS

PubMed was searched for articles relevant to consanguinity and obesity published between 1995 and 2016. Some earlier articles of interest were also consulted.

RESULTS

Although obesity is the most heritable disorder, even in outbred populations, only 2% to 5% of severe obesity cases have so far been proven to be caused by single gene mutations. In some highly consanguineous populations, a remarkably higher proportion of obesity cases because of known and novel monogenic variants has been identified (up to 30%).

CONCLUSIONS

Combining the power conferred by consanguinity with current large-capacity sequencing techniques should bring new genetic factors and molecular mechanisms to the fore, unveiling a large part of the yet-elusive neurohumoral circuitry involved in the regulation of energy homeostasis and appetite. Importantly, the undertaking of such initiatives is destined to unfold novel targets for the development of precision medicine relevant to different forms of obesity.

A null mutation in MICU2 causes abnormal mitochondrial calcium homeostasis and a severe neurodevelopmental disorder

Mitochondrial calcium homeostasis is a tightly controlled process that is required for a variety of cellular functions. The mitochondrial calcium uniporter complex plays a critical role in this process. MICU2 is a major component of the mitochondrial calcium uniporter complex and its deficiency has been shown to impair mitochondrial calcium [Ca2+]m homeostasis although the exact mechanism remains unclear. We used exome sequencing, positional mapping, and functional characterization of MICU2 deficiency to investigate the role of MICU2 in calcium homeostasis. Using combined autozygome/exome analysis, a homozygous truncating mutation in MICU2 was found to fully segregate with a neurodevelopmental disorder in the form of severe cognitive impairment, spasticity, and white matter involvement in a multiplex consanguineous family. Patient-derived MICU2-deficient cells displayed impaired [Ca2+]m homeostasis, with associated increase in mitochondrial sensitivity to oxidative stress, and abnormal regulation of inner mitochondrial membrane potential. This is the first demonstration of MICU2 deficiency in humans, which we suggest causes a distinct neurodevelopmental phenotype secondary to impaired mitochondrial calcium uniporter-mediated regulation of intracellular calcium homeostasis.

Misleading presentation of haemoglobin electrophoresis data

Haemoglobinopathies are common in sub-Saharan Africa. As such haemoglobin electrophoresis are required to inform clinical decision making. However, haemoglobin electrophoresis is an assay that detects protein at either alkaline or acidic pH. Such assays do not interrogate gene sequences but rather the product of a gene. As many post-transcriptional and post-translational modifications impact the final output of the gene (i.e. protein), presentation of such protein-based assay must accurately reflect the technique employed. It is grossly misleading and scientifically inaccurate to report cellulose acetate and/or citrate agar haemoglobin electrophoresis results as 'haemoglobin genotype'. We propose a new paradigm in which haemoglobin electrophoresis data would be presented as 'haemoglobin phenotype' at a specified pH.

FUNDING

None declared.

Whole exome sequencing as a diagnostic tool for patients with ciliopathy-like phenotypes

Ciliopathies are a group of rare disorders characterized by a high genetic and phenotypic variability, which complicates their molecular diagnosis. Hence the need to use the latest powerful approaches to faster identify the genetic defect in these patients. We applied whole exome sequencing to six consanguineous families clinically diagnosed with ciliopathy-like disease, and for which mutations in predominant Bardet-Biedl syndrome (BBS) genes had previously been excluded. Our strategy, based on first applying several filters to ciliary variants and using many of the bioinformatics tools available, allowed us to identify causal mutations in BBS2, ALMS1 and CRB1 genes in four families, thus confirming the molecular diagnosis of ciliopathy. In the remaining two families, after first rejecting the presence of pathogenic variants in common cilia-related genes, we adopted a new filtering strategy combined with prioritisation tools to rank the final candidate genes for each case. Thus, we propose CORO2B, LMO7 and ZNF17 as novel candidate ciliary genes, but further functional studies will be needed to confirm their role. Our data show the usefulness of this strategy to diagnose patients with unclear phenotypes, and therefore the success of applying such technologies to achieve a rapid and reliable molecular diagnosis, improving genetic counselling for these patients. In addition, the described pipeline also highlights the common pitfalls associated to the large volume of data we have to face and the difficulty of assigning a functional role to these changes, hence the importance of designing the most appropriate strategy according to each case.

GZF1 Mutations Expand the Genetic Heterogeneity of Larsen Syndrome

Larsen syndrome is characterized by the dislocation of large joints and other less consistent clinical findings. Heterozygous FLNB mutations account for the majority of Larsen syndrome cases, but biallelic mutations in CHST3 and B4GALT7 have been more recently described, thus confirming the existence of recessive forms of the disease. In a multiplex consanguineous Saudi family affected by severe and recurrent large joint dislocation and severe myopia, we identified a homozygous truncating variant in GZF1 through a combined autozygome and exome approach. Independently, the same approach identified a second homozygous truncating GZF1 variant in another multiplex consanguineous family affected by severe myopia, retinal detachment, and milder skeletal involvement. GZF1 encodes GDNF-inducible zinc finger protein 1, a transcription factor of unknown developmental function, which we found to be expressed in the eyes and limbs of developing mice. Global transcriptional profiling of cells from affected individuals revealed a shared pattern of gene dysregulation and significant enrichment of genes encoding matrix proteins, including P3H2, which hints at a potential disease mechanism. Our results suggest that GZF1 mutations cause a phenotype of severe myopia and significant articular involvement not previously described in Larsen syndrome.

Autozygosity reveals recessive mutations and novel mechanisms in dominant genes: implications in variant interpretation

PURPOSE

The purpose of this study is to describe recessive alleles in strictly dominant genes. Identifying recessive mutations in genes for which only dominant disease or risk alleles have been reported can expand our understanding of the medical relevance of these genes both phenotypically and mechanistically. The Saudi population is enriched for autozygosity, which enhances the homozygous occurrence of alleles, including pathogenic alleles in genes that have been associated only with a dominant inheritance pattern.

METHODS

Exome sequencing of patients from consanguineous families with likely recessive phenotypes was performed. In one family, the genotype of the deceased children was inferred from their parents due to lack of available samples.

RESULTS

We describe the identification of 11 recessive variants (5 of which are reported here for the first time) in 11 genes for which only dominant disease or risk alleles have been reported. The observed phenotypes for these recessive variants were novel (e.g., FBN2-related myopathy and CSF1R-related brain malformation and osteopetrosis), typical (e.g., ACTG2-related visceral myopathy), or an apparently healthy state (e.g., PDE11A), consistent with the corresponding mouse knockout phenotypes.

CONCLUSION

Our results show that, in the era of genomic sequencing and "reverse phenotyping," recessive variants in dominant genes should not be dismissed based on perceived "incompatibility" with the patient's phenotype before careful consideration.Genet Med advance online publication 06 April 2017.

Clinical genomics expands the morbid genome of intellectual disability and offers a high diagnostic yield

Intellectual disability (ID) is a measurable phenotypic consequence of genetic and environmental factors. In this study, we prospectively assessed the diagnostic yield of genomic tools (molecular karyotyping, multi-gene panel and exome sequencing) in a cohort of 337 ID subjects as a first-tier test and compared it with a standard clinical evaluation performed in parallel. Standard clinical evaluation suggested a diagnosis in 16% of cases (54/337) but only 70% of these (38/54) were subsequently confirmed. On the other hand, the genomic approach revealed a likely diagnosis in 58% (n=196). These included copy number variants in 14% (n=54, 15% are novel), and point mutations revealed by multi-gene panel and exome sequencing in the remaining 43% (1% were found to have Fragile-X). The identified point mutations were mostly recessive (n=117, 81%), consistent with the high consanguinity of the study cohort, but also X-linked (n=8, 6%) and de novo dominant (n=19, 13%). When applied directly on all cases with negative molecular karyotyping, the diagnostic yield of exome sequencing was 60% (77/129). Exome sequencing also identified likely pathogenic variants in three novel candidate genes (DENND5A, NEMF and DNHD1) each of which harbored independent homozygous mutations in patients with overlapping phenotypes. In addition, exome sequencing revealed de novo and recessive variants in 32 genes (MAMDC2, TUBAL3, CPNE6, KLHL24, USP2, PIP5K1A, UBE4A, TP53TG5, ATOH1, C16ORF90, SLC39A14, TRERF1, RGL1, CDH11, SYDE2, HIRA, FEZF2, PROCA1, PIANP, PLK2, QRFPR, AP3B2, NUDT2, UFC1, BTN3A2, TADA1, ARFGEF3, FAM160B1, ZMYM5, SLC45A1, ARHGAP33 and CAPS2), which we highlight as potential candidates on the basis of several lines of evidence, and one of these genes (SLC39A14) was biallelically inactivated in a potentially treatable form of hypermanganesemia and neurodegeneration. Finally, likely causal variants in previously published candidate genes were identified (ASTN1, HELZ, THOC6, WDR45B, ADRA2B and CLIP1), thus supporting their involvement in ID pathogenesis. Our results expand the morbid genome of ID and support the adoption of genomics as a first-tier test for individuals with ID.

Prenatal homozygosity mapping detects a novel mutation in CHST3 in a fetus with skeletal dysplasia and joint dislocations

In selected cases, homozygosity mapping followed by direct sequencing of one or a few carefully selected candidate genes in a prenatal setting can be beneficial to obtain diagnosis in consanguineous families.

Novel phenotypes and loci identified through clinical genomics approaches to pediatric cataract

Pediatric cataract is highly heterogeneous clinically and etiologically. While mostly isolated, cataract can be part of many multisystem disorders, further complicating the diagnostic process. In this study, we applied genomic tools in the form of a multi-gene panel as well as whole-exome sequencing on unselected cohort of pediatric cataract (166 patients from 74 families). Mutations in previously reported cataract genes were identified in 58% for a total of 43 mutations, including 15 that are novel. GEMIN4 was independently mutated in families with a syndrome of cataract, global developmental delay with or without renal involvement. We also highlight a recognizable syndrome that resembles galactosemia (a fulminant infantile liver disease with cataract) caused by biallelic mutations in CYP51A1. A founder mutation in RIC1 (KIAA1432) was identified in patients with cataract, brain atrophy, microcephaly with or without cleft lip and palate. For non-syndromic pediatric cataract, we map a novel locus in a multiplex consanguineous family on 4p15.32 where exome sequencing revealed a homozygous truncating mutation in TAPT1. We report two further candidates that are biallelically inactivated each in a single cataract family: TAF1A (cataract with global developmental delay) and WDR87 (non-syndromic cataract). In addition to positional mapping data, we use iSyTE developmental lens expression and gene-network analysis to corroborate the proposed link between the novel candidate genes and cataract. Our study expands the phenotypic, allelic and locus heterogeneity of pediatric cataract. The high diagnostic yield of clinical genomics supports the adoption of this approach in this patient group.