Exome sequencing unravels genetic variants associated with chronic kidney disease in Saudi Arabian patients.

The use of genetic testing within nephrology is increasing and its diagnostic yield depends on the methods utilized, patient selection criteria, and population characteristics. We performed exome sequencing (ES) analysis on 102 chronic kidney disease (CKD) patients with likely genetic kidney disease. Patients had diverse CKD subtypes with/without consanguinity, positive family history, and possible hereditary renal syndrome with extra-renal abnormalities or progressive kidney disease of unknown etiology. The identified genetic variants associated with the observed kidney phenotypes were then confirmed and reported. End-stage kidney disease was reported in 51% of the cohort and a family history of kidney disease in 59%, while known consanguinity was reported in 54%. Pathogenic/likely pathogenic variants were identified in 43 patients with a diagnostic yield of 42%, and clinically associated variants of unknown significance (VUS) were identified in further 21 CKD patients (21%). A total of eight novel predicted pathogenic variants and eight VUS were detected. The clinical utility of ES within the nephrology clinic was demonstrated allowing patient management to be disease-specific. In this cohort, ES detected a diagnostic molecular abnormality in 42% of patients with CKD phenotypes. Positive family history and high rates of consanguinity likely contributed to this high diagnostic yield.

Prenatal exome sequencing and chromosomal microarray analysis in fetal structural anomalies in a highly consanguineous population reveals a propensity of ciliopathy genes causing multisystem phenotypes.

Fetal abnormalities are detected in 3% of all pregnancies and are responsible for approximately 20% of all perinatal deaths. Chromosomal microarray analysis (CMA) and exome sequencing (ES) are widely used in prenatal settings for molecular genetic diagnostics with variable diagnostic yields. In this study, we aimed to determine the diagnostic yield of trio-ES in detecting the cause of fetal abnormalities within a highly consanguineous population. In families with a history of congenital anomalies, a total of 119 fetuses with structural anomalies were recruited and DNA from invasive samples were used together with parental DNA samples for trio-ES and CMA. Data were analysed to determine possible underlying genetic disorders associated with observed fetal phenotypes. The cohort had a known consanguinity of 81%. Trio-ES led to diagnostic molecular genetic findings in 59 fetuses (with pathogenic/likely pathogenic variants) most with multisystem or renal abnormalities. CMA detected chromosomal abnormalities compatible with the fetal phenotype in another 7 cases. Monogenic ciliopathy disorders with an autosomal recessive inheritance were the predominant cause of multisystem fetal anomalies (24/59 cases, 40.7%) with loss of function variants representing the vast majority of molecular genetic abnormalities. Heterozygous de novo pathogenic variants were found in four fetuses. A total of 23 novel variants predicted to be associated with the phenotype were detected. Prenatal trio-ES and CMA detected likely causative molecular genetic defects in a total of 55% of families with fetal anomalies confirming the diagnostic utility of trio-ES and CMA as first-line genetic test in the prenatal diagnosis of multisystem fetal anomalies including ciliopathy syndromes.

Residual risk for additional recessive diseases in consanguineous couples.

PURPOSE: Consanguineous couples are typically counseled based on familial pathogenic variants identified in affected children. The residual risk for additional autosomal recessive (AR) variants, however, remains largely understudied. METHODS: First, we surveyed pedigrees of 1,859 consanguineous families for evidence of more than one AR disease. Second, we mined our database of 1,773 molecularly tested consanguineous families to identify those with more than one AR disease. Finally, we surveyed 88 women from consanguineous unions who have undergone targeted prenatal testing for a familial AR variant and followed the pregnancy outcome (n = 144). RESULTS: We found suggestive evidence of more than one AR disease in 1.94% of consanguineous pedigrees surveyed. Of 1,773 molecularly characterized consanguineous families, 2.93% had evidence of at least two AR diseases (3.54% for first cousin or closer and 2.72% for second cousin or more distant). Furthermore, we found that in 2.78% of pregnancies negative for the familial variant, the pregnancy outcome was a child with a different AR disease. CONCLUSION: Our results show that when counseling consanguineous couples for a familial AR variant, ~3% residual risk for additional AR variants should be discussed. This suggests that a broader testing strategy in consanguineous couples should be considered.

The genotypic and phenotypic spectrum of pycnodysostosis in Saudi Arabia: Novel variants and clinical findings.

Pycnodysostosis is characterized by short stature, osteosclerosis, acro-osteolysis, increased tendency of fractures, and distinctive dysmorphic features. It is a rare autosomal recessive disease caused by biallelic CTSK mutations. The clinical details of 18 patients from Saudi Arabia were reviewed. Short stature, osteopetrosis, acro-osteolysis, and distinctive facial dysmorphism were documented in all cases. Our results highlight the significant complications associated with this disease. The large anterior fontanelle is one of the cardinal signs of this disease; however, half of our patients had small fontanelles and a quarter had craniosynostosis, which caused optic nerve compression. Sleep apnea was of the major complications in three patients. Bone fracture can be a presenting symptom, and in our patients it mainly occurred after the age of 3 years. Bone marrow suppression was seen in a single patient of our cohort who was misdiagnosed initially with malignant osteopetrosis. In this study, we also describe two novel (c.5G > A [p.Trp2Ter], c.538G > A [p.Gly180Ser]) and two reported (c.244-29 A > G, c.830C > T [p.Ala277Val]) CTSK mutations. Our results indicate that the recurrent intronic variant, c.244-29 A > G is likely to be a founder mutation, as it was found in 78% (14/18 patients) of our cohort belonging to the same tribe.

Genotype-phenotype correlation of 33 patients with maple syrup urine disease.

Maple syrup urine disease (MSUD) is a rare autosomal recessive inherited disorder due to defects in the branched-chain α-ketoacid dehydrogenase complex (BCKDC). MSUD varies in severity and its clinical spectrum is quite broad, ranging from mild to severe phenotypes. Thirty-three MSUD patients were recruited into this study for molecular genetic variant profiling and genotype-phenotype correlation. Except for one patient, all other patients presented with the classic neonatal form of the disease. Seventeen different variants were detected where nine were novel. The detected variants spanned across the entire BCKDHA, BCKDHB and DBT genes. All variants were in homozygous forms. The commonest alterations were nonsense and frameshift variants, followed by missense variants. For the prediction of variant's pathogenicity, we used molecular modeling and several in silico tools including SIFT, Polyphen2, Condel, and Provean. In addition, six other tools were used for the prediction of the conservation of the variants' sites including Eigen-PC, GERP++, SiPhy, PhastCons vertebrates and primates, and PhyloP100 rank scores. Herein, we presented a comprehensive characterization of a large cohort of patients with MSUD. The clinical severity of the variants' phenotypes was well correlated with the genotypes. The study underscores the importance of the use of in silico analysis of MSUD genotypes for the prediction of the clinical outcomes in patients with MSUD.

Identification of TMC1 as a relatively common cause for nonsyndromic hearing loss in the Saudi population.

Hearing loss (HL) is the most common sensory disorder worldwide and genetic factors contribute to approximately half of congenital HL cases. HL is subject to extensive genetic heterogeneity, rendering molecular diagnosis difficult. Mutations of the transmembrane channel-like 1 (TMC1) gene cause hearing defects in humans and mice. The precise function of TMC1 protein in the inner ear is unknown, although it is predicted to be involved in functional maturation of cochlear hair cells. TMC1 mutations result in autosomal recessive (DFNB7/11) and sometimes dominant (DFNA36) nonsyndromic HL. Mutations in TMC1 are responsible for a significant portion of HL, particularly in consanguineous populations. To evaluate the importance of TMC1 mutations in the Saudi population, we used a combination of autozygome-guided candidate gene mutation analysis and targeted next generation sequencing in 366 families with HL previously shown to lack mutations in GJB2. We identified 12 families that carried five causative TMC1 mutations; including three novel (c.362+3A > G; c.758C > T [p.Ser253Phe]; c.1396_1398delACC [p.Asn466del]) and two reported mutations (c.100C > T [p.Arg34Ter]; c.1714G > A [p.Asp572Asn]). Each of the identified recessive mutation was classified as severe, by both age of onset and severity of HL. Similarly, consistent with the previously reported dominant variant p.Asp572Asn, the HL phenotype was progressive. Eight families in our cohort were found to share the pathogenic p.Arg34Ter mutation and linkage disequilibrium was observed between p.Arg34Ter and SNPs investigated. Our results indicate that TMC1 mutations account for about 3.3% (12/366) of Saudi HL cases and that the recurrent TMC1 mutation p.Arg34Ter is likely to be a founder mutation.

Bialleleic PKD1 mutations underlie early-onset autosomal dominant polycystic kidney disease in Saudi Arabian families.

BACKGROUND: Polycystic kidney disease (PKD) is one of the most common genetic renal diseases and may be inherited in an autosomal dominant or autosomal recessive pattern. Pathogenic variants in two major genes, PKD1 and PKD2, and two rarer genes, GANAB and DNAJB11, cause autosomal dominant PKD (ADPKD). Early onset and severe PKD can occur with PKD1 and PKD2 pathogenic variants and such phenotypes may be modified by second alleles inherited in trans. Homozygous or compound heterozygous hypomorphic PKD1 variants may also cause a moderate to severe disease PKD phenotype. METHODS: Targeted renal gene panel followed by Sanger sequencing of PKD1 gene were employed to investigate molecular causes in early onset PKD patients. RESULTS: In this study, we report four consanguineous Saudi Arabian families with early onset PKD which were associated with biallelic variants in PKD1 gene. CONCLUSIONS: Our findings confirm that PKD1 alleles may combine to produce severe paediatric onset PKD mimicking the more severe autosomal recessive ciliopathy syndromes associated with PKD. Screening of parents of such children may also reveal subclinical PKD phenotypes.

ROR1 is essential for proper innervation of auditory hair cells and hearing in humans and mice.

Hair cells of the inner ear, the mechanosensory receptors, convert sound waves into neural signals that are passed to the brain via the auditory nerve. Little is known about the molecular mechanisms that govern the development of hair cell-neuronal connections. We ascertained a family with autosomal recessive deafness associated with a common cavity inner ear malformation and auditory neuropathy. Via whole-exome sequencing, we identified a variant (c.2207G>C, p.R736T) in ROR1 (receptor tyrosine kinase-like orphan receptor 1), cosegregating with deafness in the family and absent in ethnicity-matched controls. ROR1 is a tyrosine kinase-like receptor localized at the plasma membrane. At the cellular level, the mutation prevents the protein from reaching the cellular membrane. In the presence of WNT5A, a known ROR1 ligand, the mutated ROR1 fails to activate NF-κB. Ror1 is expressed in the inner ear during development at embryonic and postnatal stages. We demonstrate that Ror1 mutant mice are severely deaf, with preserved otoacoustic emissions. Anatomically, mutant mice display malformed cochleae. Axons of spiral ganglion neurons show fasciculation defects. Type I neurons show impaired synapses with inner hair cells, and type II neurons display aberrant projections through the cochlear sensory epithelium. We conclude that Ror1 is crucial for spiral ganglion neurons to innervate auditory hair cells. Impairment of ROR1 function largely affects development of the inner ear and hearing in humans and mice.

Inclusion of joint laxity, recurrent patellar dislocation, and short distal ulnae as a feature of Van Den Ende-Gupta syndrome: a case report.

BACKGROUND: Van Den Ende-Gupta Syndrome (VDEGS) is an extremely rare autosomal recessive syndrome with less than 20 reported families (approximately 40 patients) in the worldwide literature. CASE PRESENTATION: We have assessed one consanguineous Saudi family with typical features of VDEGS. Two siblings were affected with almost identical features; including blepharophimosis, arachnodactyly, flexion contractures of the elbows, camptodactyly, slender ribs, hooked lateral clavicular ends, and bilateral radial head dislocations. Both patients had several unusual features; including joint laxity, flat feet, recurrent patellar dislocations, and bilateral short distal ulnae. Full sequencing of SCARF2 revealed a homozygous mutation c.773G > A (p. Cys258Tyr) in both affected children. The parents (both with no abnormalities) were heterozygous for the same mutation. CONCLUSION: Joint laxity, recurrent patellar dislocations, and short distal ulnae should be included as part of the clinical spectrum of VDEGS.

Delineating the phenotypic spectrum of hyperphosphatasia with mental retardation syndrome 4 in 14 patients of Middle-Eastern origin.

Hyperphosphatasia with mental retardation syndrome 4 (HPMRS4) is a rare autosomal recessive condition caused by an impairment of glycosylphophatidylinositol biosynthesis. The cardinal features of HPMRS4 include; characteristic facial features, severe intellectual disability and various neurologic abnormalities. We report here detailed clinical, biochemical, and molecular findings of 14 patients clinically suspected to have HPMRS4, from three Middle-Eastern Countries; Saudi Arabia, Qatar, and Oman. All patients in our series presented with the cardinal features pointing to HPMRS4 and with an elevated alkaline phosphatase level. Five patients had megalocornea, which have been reported recently in an Arab patient. Additionally, fracture, bilateral coxa valga, camptodactyly, truncal obesity, and hyperpigmented macules of the upper thigh, each was seen once and was not described before with HPMRS4. Additional clinical and radiological findings are described, supporting the novel clinical and radiological findings recently described in Egyptian patients. The utilization of homozygosity mapping coupled with PGAP3 sequencing and whole exome sequencing facilitated the mutation detection in these patients. These missense mutations include c.320C > T (p.S107 L), c.850C > T (p.H284Y), and c.851A > G (p.H284R) in the PGAP3 gene. We believe that the recurrent mutations identified in our cohort may represent founder mutations in big tribes from a certain geographical region of Saudi Arabia, Qatar, and Oman. Therefore, in case of a clinical suspicion of HPMRS4 in these populations, targeted genetic testing for the identified mutations should be performed first to expedite the genetic diagnosis.

SLC5A1 Mutations in Saudi Arabian Patients With Congenital Glucose-Galactose Malabsorption.

Congenital glucose-galactose malabsorption (cGGM) is a rare autosomal recessive disorder, caused by mutations in the SLC5A1 gene, encoding the sodium/glucose cotransporter 1, which may result in severe life-threatening osmotic diarrhea due to the accumulation of unabsorbed sugars in the intestinal lumen. If treated early with elimination of glucose and galactose from the diet, patients usually recover and develop normally. We present clinical and molecular data from 16 unrelated cGGM diagnosed Saudi patients from consanguineous families with majority of them having previous positive family history of cGGM. Sanger sequencing for the full coding regions of SLC5A1 for all patients resulted in the identification of 4 allelic variants in a homozygous state. Two mutations are novel; c.265G>A (p.G89R) and c.1304 G>A (p.G435D), and 2 have been previously reported to cause cGGM, c.765 C>G (p.C255W) and c.1136 G>A (p.R379Q). This is the first report delineating the clinical and molecular basis of cGGM in patients from this region.

Two novel LHX3 mutations in patients with combined pituitary hormone deficiency including cervical rigidity and sensorineural hearing loss.

BACKGROUND: Congenital combined pituitary hormone deficiency (CPHD) is a rare heterogeneous group of conditions. CPHD-type 3 (CPHD3; MIM# 221750) is caused by recessive mutations in LHX3, a LIM-homeodomain transcription factor gene. The isoforms of LHX3 are critical for pituitary gland formation and specification of the anterior pituitary hormone-secreting cell types. They also play distinct roles in the development of neuroendocrine and auditory systems. CASE PRESENTATION: Here, we summarize the clinical, endocrinological, radiological and molecular features of three patients from two unrelated families. Clinical evaluation revealed severe CPHD coupled with cervical vertebral malformations (rigid neck, scoliosis), mild developmental delay and moderate sensorineural hearing loss (SNHL). The patients were diagnosed with CPHD3 based on the array of hormone deficiencies and other associated syndromic symptoms, suggestive of targeted LHX3 gene sequencing. A novel missense mutation c.437G > T (p. Cys146Phe) and a novel nonsense mutation c.466C > T (p. Arg156Ter), both in homozygous forms, were found. The altered Cys146 resides in the LIM2 domain of the encoded protein and is a phylogenetically conserved residue, which mediates LHX3 transcription factor binding with a zinc cation. The p. Arg156Ter is predicted to result in a severely truncated protein, lacking the DNA binding homeodomain. CONCLUSIONS: Considering genotype/phenotype correlation, we suggest that the presence of SNHL and limited neck rotation should be considered in the differential diagnosis of CPHD3 to facilitate molecular diagnosis. This report describes the first LHX3 mutations from Saudi patients and highlights the importance of combining molecular diagnosis with the clinical findings. In addition, it also expands the knowledge of LHX3-related CPHD3 phenotype and the allelic spectrum for this gene.

Clinical genomics can facilitate countrywide estimation of autosomal recessive disease burden.

BACKGROUND: Most autosomal recessive diseases are rare, but they collectively account for a substantial proportion of disease burden, especially in consanguineous populations. Estimation of this disease burden, however, is hampered by many factors, including lack of countrywide registries. Establishing carrier frequency can be a practical surrogate to estimate disease burden, although the requirement of a large representative cohort may be challenging. PURPOSE: We propose that the application of clinical genomics in the diagnostic setting offers a unique opportunity to estimate carrier frequency in the population as a secondary benefit. METHODS: We used a data set of ~7,100 patients who underwent genomic testing for various Mendelian disorders to estimate the carrier frequency. RESULTS: We were able to calculate the frequency of 259 confirmed founder recessive mutations. We found the corresponding disease burden to be, at minimum, ~7 per 1,000 children born to first-cousin parents, with disorders related to intellectual disability and vision impairment being the most common. CONCLUSION: Our approach can be utilized to inform the design of new policies for the prevention of genetic disorders and highlights an important secondary benefit of clinical genomics.Genet Med 18 12, 1244-1249.

Variation in DNAH1 may contribute to primary ciliary dyskinesia.

BACKGROUND: Primary Ciliary Dyskinesia (PCD) is a genetically heterogeneous ciliopathy caused by ultrastructural defects in ciliary or flagellar structure and is characterized by a number of clinical symptoms including recurrent respiratory infections progressing to permanent lung damage and infertility. CASE PRESENTATION: Here we describe our search to delineate the molecular basis in two affected sisters with clinically diagnosed PCD from a consanguineous Saudi Arabian family, in which all known genes have been excluded. A homozygosity mapping-based approach was utilized that ultimately identified one single affected-shared region of homozygosity using 10 additional unaffected family members. A plausible candidate gene was directly sequenced and analyzed for mutations. A novel homozygous missense aberration (p.Lys1154Gln) was identified in both sisters in the DNAH1 gene that segregated completely with the disease phenotype. Further confirmation of this interesting variant was provided by exome-wide analysis in the proband. CONCLUSION: Molecular variation in DNAH1 may play a role in PCD and its potential contribution should be considered in patients where all known genes are excluded.

Tracheobronchial anomalies in a patient with Schimke immuno-osseous dysplasia (SIOD).

Schimke Immuno-osseous Dysplasia (SIOD) is a rare genetic disorder with multiple systemic manifestations. Pulmonary manifestations have been described but not well characterized. They are believed to be secondary to decreased elasticity, and include emphysema, pulmonary hypertension and bronchiectasis. We describe a 24-year-old female patient with SIOD with tracheobronchial anomalies not reported before, including tracheomalacia and a tracheal bronchus with Epstein-Barr virus (EBV) related leiomyoma causing endobronchial obstruction. Such anomalies, in addition to the difficult upper airway associated with SIOD present specific challenges during management. This case reports documents tracheobronchial abnormalities that have not been described before in SIOD.

Novel PHKG2 mutation causing GSD IX with prominent liver disease: report of three cases and review of literature.

Glycogen storage disease type IX (GSD IX) is a common form of glycogenosis due to mutations in PHKA1, PHKA2, or PHKB and PHKG2 genes resulting in the deficiency of phosphorylase kinase. The first two genes are X-linked while the latter two follow an autosomal recessive inheritance. The majority of cases of GSD IX are attributed to defects in PHKA2 which usually cause a mild disease. We report three patients with PHKG2-related GSD IX presenting with significant hepatic involvement, fibrosis, and cirrhosis. Interestingly, the homozygosity mapping resolved a dilemma about an erroneously normal phosphorylase kinase activity in patient 1. The novel mutation found in all the three patients (p.G220E) affects the catalytic subunit of the phosphorylase kinase. Increasing evidence shows that patients with PHKG2 mutations have a severe hepatic phenotype within the heterogeneous GSD IX disorder. Therefore, defect in PHKG2 should be considered in patients with suspected glycogenosis associated with significant liver fibrosis and cirrhosis.

Clinical characteristics and long-term management for patients with vitamin D-dependent rickets type II: a retrospective study at a single center in Saudi Arabia.

Introduction: Hereditary Vitamin D-dependent rickets type II (HVDDR-type II) is a rare autosomal recessive disorder caused by molecular variation in the gene encoding the vitamin D receptor (VDR). This study aims to evaluate phenotype and genotype characteristics and long-term follow-up of the largest group of patients with (HVDDR-type II) in Saudi Arabia. Methodology: We conducted a retrospective chart review to collect the clinical, biochemical, and genetic data for all HVDDR-type II patients currently receiving treatment at King Faisal Specialist Hospital & Research Centre, Riyadh, Saudi Arabia. Results: A total of 42 patients, 57.1% female, and 42.9% male were included in the study. Seven patients were treated with high doses of oral calcium, while 35 patients were treated with IV calcium infusion. The median age at presentation was 15.5 months. Alopecia was found in 97.6%, 21.4% presented with bowing legs, 14.3% with delayed walking, 9.5% with seizure, and 2.4% presented with respiratory failure, while a family history of the disease was positive in 71.4% of total patients. Molecular genetic testing of the VDR gene in our cohort identified six different gene variants c.885 C>A (p.Tyr295Ter), c.88 C>T (p.Arg30Ter), c.1036G>A (p.Val346Met), c.820C>T (p.Arg274Cys), c.803 T>C (p.Ile268Thr), and c.2T>G (p.Met1?). Conclusion: We are describing the largest cohort of patients with HVDDR-type II, their clinical biochemical findings, and the most prevalent genetic variants in our population.

Identification of Variants Underlying Phenylalanine Hydroxylase Deficiency in Saudi Arabia.

Background: Deleterious mutations in the human gene phenylalanine hydroxylase (PAH) encoding the phenylalanine hydroxylase enzyme give rise to classic phenylketonuria and hyperphenylalaninemia. Our study was designed to characterize the spectrum of variants in the PAH gene in Saudi patients. Materials and Methods: We screened a cohort of 72 Saudi patients with clinical and biochemical diagnoses of hyperphenylalaninemia at the largest tertiary care center in Saudi Arabia; the King Faisal Specialist Hospital and Research Center (KFSH&RC), Riyadh. All patient's charts were reviewed under an approved study by Institutional Review Board. Results: Twenty-one different PAH variants were identified among the 144 PAH alleles assessed by targeted gene sequencing. Within the studied cohort, 60 of 72 patients had homozygous mutations with the the remaining 12 being compound heterozygotes. The most prevalent of the disease alleles identified in this study was the p.(Arg252Trp) mutation, which accounted for 38 of 144 alleles (26.4%). With the high incidence of genetic disorders in the population, religiously permissible preventive reproductive measures are a priority in our practice. Prenatal diagnoses carried out on four fetuses revealed two that were homozygous for PAH pathogenic variants. In addition, pre-implantation genetic diagnoses were initiated for 19 families. Eight of these families completed more than one full cycle of treatment, from which one healthy newborn was delivered. Conclusions: This study describes the spectrum of PAH variants in the Saudi population and highlights the molecular heterogeneity underlying phenylketonuria and hyperphenylalaninemia. These results add to the existing knowledge about PAH variants in Middle Eastern Countries. These results can be further translated to provide: informed counseling; cascade carrier testing in extended family members; and pre-marital screening.

Wild-type S100A3 and S100A13 restore calcium homeostasis and mitigate mitochondrial dysregulation in pulmonary fibrosis patient-derived cells.

Patients with digenic S100A3 and S100A13 mutations exhibited an atypical and progressive interstitial pulmonary fibrosis, with impaired intracellular calcium homeostasis and mitochondrial dysfunction. Here we provide direct evidence of a causative effect of the mutation on receptor mediated calcium signaling and calcium store responses in control cells transfected with mutant S100A3 and mutant S100A13. We demonstrate that the mutations lead to increased mitochondrial mass and hyperpolarization, both of which were reversed by transfecting patient-derived cells with the wild type S100A3 and S100A13, or extracellular treatment with the recombinant proteins. In addition, we demonstrate increased secretion of inflammatory mediators in patient-derived cells and in control cells transfected with the mutant-encoding constructs. These findings indicate that treatment of patients' cells with recombinant S100A3 and S100A13 proteins is sufficient to normalize most of cellular responses, and may therefore suggest the use of these recombinant proteins in the treatment of this devastating disease.

Genetic association analysis of 77,539 genomes reveals rare disease etiologies.

The genetic etiologies of more than half of rare diseases remain unknown. Standardized genome sequencing and phenotyping of large patient cohorts provide an opportunity for discovering the unknown etiologies, but this depends on efficient and powerful analytical methods. We built a compact database, the 'Rareservoir', containing the rare variant genotypes and phenotypes of 77,539 participants sequenced by the 100,000 Genomes Project. We then used the Bayesian genetic association method BeviMed to infer associations between genes and each of 269 rare disease classes assigned by clinicians to the participants. We identified 241 known and 19 previously unidentified associations. We validated associations with ERG, PMEPA1 and GPR156 by searching for pedigrees in other cohorts and using bioinformatic and experimental approaches. We provide evidence that (1) loss-of-function variants in the Erythroblast Transformation Specific (ETS)-family transcription factor encoding gene ERG lead to primary lymphoedema, (2) truncating variants in the last exon of transforming growth factor-ß regulator PMEPA1 result in Loeys-Dietz syndrome and (3) loss-of-function variants in GPR156 give rise to recessive congenital hearing impairment. The Rareservoir provides a lightweight, flexible and portable system for synthesizing the genetic and phenotypic data required to study rare disease cohorts with tens of thousands of participants.