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.

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.

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.

Genomic and phenotypic delineation of congenital microcephaly.

PURPOSE: Congenital microcephaly (CM) is an important birth defect with long term neurological sequelae. We aimed to perform detailed phenotypic and genomic analysis of patients with Mendelian forms of CM. METHODS: Clinical phenotyping, targeted or exome sequencing, and autozygome analysis. RESULTS: We describe 150 patients (104 families) with 56 Mendelian forms of CM. Our data show little overlap with the genetic causes of postnatal microcephaly. We also show that a broad definition of primary microcephaly -as an autosomal recessive form of nonsyndromic CM with severe postnatal deceleration of occipitofrontal circumference-is highly sensitive but has a limited specificity. In addition, we expand the overlap between primary microcephaly and microcephalic primordial dwarfism both clinically (short stature in >52% of patients with primary microcephaly) and molecularly (e.g., we report the first instance of CEP135-related microcephalic primordial dwarfism). We expand the allelic and locus heterogeneity of CM by reporting 37 novel likely disease-causing variants in 27 disease genes, confirming the candidacy of ANKLE2, YARS, FRMD4A, and THG1L, and proposing the candidacy of BPTF, MAP1B, CCNH, and PPFIBP1. CONCLUSION: Our study refines the phenotype of CM, expands its genetics heterogeneity, and informs the workup of children born with this developmental brain defect.

Human knockout research: new horizons and opportunities.

Although numerous approaches have been pursued to understand the function of human genes, Mendelian genetics has by far provided the most compelling and medically actionable dataset. Biallelic loss-of-function (LOF) mutations are observed in the majority of autosomal recessive Mendelian disorders, representing natural human knockouts and offering a unique opportunity to study the physiological and developmental context of these genes. The restriction of such context to 'disease' states is artificial, however, and the recent ability to survey entire human genomes for biallelic LOF mutations has revealed a surprising landscape of knockout events in 'healthy' individuals, sparking interest in their role in phenotypic diversity beyond disease causation. As I discuss in this review, the potentially wide implications of human knockout research warrant increased investment and multidisciplinary collaborations to overcome existing challenges and reap its benefits.

Identification of a Recognizable Progressive Skeletal Dysplasia Caused by RSPRY1 Mutations.

Skeletal dysplasias are highly variable Mendelian phenotypes. Molecular diagnosis of skeletal dysplasias is complicated by their extreme clinical and genetic heterogeneity. We describe a clinically recognizable autosomal-recessive disorder in four affected siblings from a consanguineous Saudi family, comprising progressive spondyloepimetaphyseal dysplasia, short stature, facial dysmorphism, short fourth metatarsals, and intellectual disability. Combined autozygome/exome analysis identified a homozygous frameshift mutation in RSPRY1 with resulting nonsense-mediated decay. Using a gene-centric "matchmaking" system, we were able to identify a Peruvian simplex case subject whose phenotype is strikingly similar to the original Saudi family and whose exome sequencing had revealed a likely pathogenic homozygous missense variant in the same gene. RSPRY1 encodes a hypothetical RING and SPRY domain-containing protein of unknown physiological function. However, we detect strong RSPRY1 protein localization in murine embryonic osteoblasts and periosteal cells during primary endochondral ossification, consistent with a role in bone development. This study highlights the role of gene-centric matchmaking tools to establish causal links to genes, especially for rare or previously undescribed clinical entities.

A lethal neonatal phenotype of mitochondrial short-chain enoyl-CoA hydratase-1 deficiency.

Short-chain enoyl-CoA hydratase (SCEH) is a mitochondrial enzyme involved in the oxidation of fatty acids and the catabolic pathway of valine and, to a lesser extent, isoleucine. Deficiency of this enzyme was recently shown to cause an early childhood Leigh syndrome phenotype. The few reported patients were compound heterozygotes for two missense or missense with truncating variants in ECHS1 that encodes SCEH. We describe two siblings with severe refractory lactic acidosis and death within the first 2 days of life. Following negative clinical whole-exome and whole-genome sequencing, we resorted to autozygome/exome analysis on research basis and identified a homozygous splice site mutation (c.88+5G>A) in the two cases. Analysis of cDNA confirmed complete replacement of the normal transcript with an aberrant transcript (r.88_89ins 88+1_88+11) predicting premature truncation of the protein [p.(Ala31Glufs*23)]. Furthermore, quantitative reverse transcriptase polymerase chain reaction (RTPCR) showed marked reduction in ECHS1, most likely nonsense-mediated decay (NMD)-mediated. This is the first report of homozygosity for a truncating mutation in ECHS1, which may explain the severe phenotype. Our report highlights the need to consider SCEH deficiency in patients with lethal neonatal lactic acidosis, and the potentially limited sensitivity of untargeted genomic sequencing towards non-canonical splicing mutations, which may explain at least some of the 'negative' cases on clinical exome/genome sequencing.

Mutation in PLK4, encoding a master regulator of centriole formation, defines a novel locus for primordial dwarfism.

BACKGROUND: Primordial dwarfism (PD) is a heterogeneous clinical entity characterised by severe prenatal and postnatal growth deficiency. Despite the recent wave of disease gene discovery, the causal mutations in many PD patients remain unknown. OBJECTIVE: To describe a PD family that maps to a novel locus. METHODS: Clinical, imaging and laboratory phenotyping of a new family with PD followed by autozygosity mapping, linkage analysis and candidate gene sequencing. RESULTS: We describe a multiplex consanguineous Saudi family in which two full siblings and one half-sibling presented with classical features of Seckel syndrome in addition to optic nerve hypoplasia. We were able to map the phenotype to a single novel locus on 4q25-q28.2, in which we identified a five base-pair deletion in PLK4, which encodes a master regulator of centriole duplication. CONCLUSIONS: Our discovery further confirms the role of genes involved in centriole biology in the pathogenesis of PD.

Mutation in ADAT3, encoding adenosine deaminase acting on transfer RNA, causes intellectual disability and strabismus.

BACKGROUND: Intellectual disability (ID) is one of the most common forms of disability worldwide, displaying a wide range of aetiologies and affecting nearly 2% of the global population. OBJECTIVE: To describe a novel autosomal recessive form of ID with strabismus and its underlying aetiology. MATERIALS AND METHODS: Autozygosity mapping, linkage analysis and exome sequencing were performed in a large multiplex consanguineous family that segregates ID and strabismus. Exome sequencing was independently performed in three other consanguineous families segregating the same disease. Direct sequencing of the resulting candidate gene was performed in four additional families with the same phenotype. RESULTS: A single missense mutation was identified in ADAT3 in all studied families on an ancient ancestral haplotype. This gene encodes one of two eukaryotic proteins that are necessary for the deamination of adenosine at position 34 to inosine in t-RNA. Our results show the first human mutation in the t-RNA editing machinery and expand the landscape of pathways involved in the pathogenesis of ID.

Beyond the exome: utility of long-read whole genome sequencing in exome-negative autosomal recessive diseases.

BACKGROUND: Long-read whole genome sequencing (lrWGS) has the potential to address the technical limitations of exome sequencing in ways not possible by short-read WGS. However, its utility in autosomal recessive Mendelian diseases is largely unknown. METHODS: In a cohort of 34 families in which the suspected autosomal recessive diseases remained undiagnosed by exome sequencing, lrWGS was performed on the Pacific Bioscience Sequel IIe platform. RESULTS: Likely causal variants were identified in 13 (38%) of the cohort. These include (1) a homozygous splicing SV in TYMS as a novel candidate gene for lethal neonatal lactic acidosis, (2) a homozygous non-coding SV that we propose impacts STK25 expression and causes a novel neurodevelopmental disorder, (3) a compound heterozygous SV in RP1L1 with complex inheritance pattern in a family with inherited retinal disease, (4) homozygous deep intronic variants in LEMD2 and SNAP91 as novel candidate genes for neurodevelopmental disorders in two families, and (5) a promoter SNV in SLC4A4 causing non-syndromic band keratopathy. Surprisingly, we also encountered causal variants that could have been identified by short-read exome sequencing in 7 families. The latter highlight scenarios that are especially challenging at the interpretation level. CONCLUSIONS: Our data highlight the continued need to address the interpretation challenges in parallel with efforts to improve the sequencing technology itself. We propose a path forward for the implementation of lrWGS sequencing in the setting of autosomal recessive diseases in a way that maximizes its utility.