In this study, we report the experience of the only reference clinical next-generation sequencing lab in Saudi Arabia with the first 1000 families who span a wide-range of suspected Mendelian phenotypes. A total of 1019 tests were performed in the period of March 2016–December 2016 comprising 972 solo (index only), 14 duo (parents or affected siblings only), and 33 trio (index and parents). Multigene panels accounted for 672 tests, while whole exome sequencing (WES) represented the remaining 347 tests. Pathogenic or likely pathogenic variants that explain the clinical indications were identified in 34% (27% in panels and 43% in exomes), spanning 279 genes and including 165 novel variants. While recessive mutations dominated the landscape of solved cases (71% of mutations, and 97% of which are homozygous), a substantial minority (27%) were solved on the basis of dominant mutations. The highly consanguineous nature of the study population also facilitated homozygosity for many private mutations (only 32.5% of the recessive mutations are founder), as well as the first instances of recessive inheritance of previously assumed strictly dominant disorders (involving ITPR1, VAMP1, MCTP2, and TBP). Surprisingly, however, dual molecular diagnosis was only observed in 1.5% of cases. Finally, we have encountered candidate variants in 75 genes (ABHD6, ACY3, ADGRB2, ADGRG7, AGTPBP1, AHNAK2, AKAP6, ASB3, ATXN1L, C17orf62, CABP1, CCDC186, CCP110, CLSTN2, CNTN3, CNTN5, CTNNA2, CWC22, DMAP1, DMKN, DMXL1, DSCAM, DVL2, ECI1, EP400, EPB41L5, FBXL22, GAP43, GEMIN7, GIT1, GRIK4, GRSF1, GTRP1, HID1, IFNL1, KCNC4, LRRC52, MAP7D3, MCTP2, MED26, MPP7, MRPS35, MTDH, MTMR9, NECAP2, NPAT, NRAP, PAX7, PCNX, PLCH2, PLEKHF1, PTPN12, QKI, RILPL2, RIMKLA, RIMS2, RNF213, ROBO1, SEC16A, SIAH1, SIRT2, SLAIN2, SLC22A20, SMDT1, SRRT, SSTR1, ST20, SYT9, TSPAN6, UBR4, VAMP4, VPS36, WDR59, WDYHV1, and WHSC1) not previously linked to human phenotypes and these are presented to accelerate post-publication matchmaking. Two of these genes were independently mutated in more than one family with similar phenotypes, which substantiates their link to human disease (AKAP6 in intellectual disability and UBR4 in early dementia). If the novel candidate disease genes in this cohort are independently confirmed, the yield of WES will have increased to 83%, which suggests that most “negative” clinical exome tests are unsolved due to interpretation rather than technical limitations.Electronic supplementary materialThe online version of this article (doi:10.1007/s00439-017-1821-8) contains supplementary material, which is available to authorized users.

BackgroundCiliopathies are clinically diverse disorders of the primary cilium. Remarkable progress has been made in understanding the molecular basis of these genetically heterogeneous conditions; however, our knowledge of their morbid genome, pleiotropy, and variable expressivity remains incomplete.ResultsWe applied genomic approaches on a large patient cohort of 371 affected individuals from 265 families, with phenotypes that span the entire ciliopathy spectrum. Likely causal mutations in previously described ciliopathy genes were identified in 85% (225/265) of the families, adding 32 novel alleles. Consistent with a fully penetrant model for these genes, we found no significant difference in their “mutation load” beyond the causal variants between our ciliopathy cohort and a control non-ciliopathy cohort. Genomic analysis of our cohort further identified mutations in a novel morbid gene TXNDC15, encoding a thiol isomerase, based on independent loss of function mutations in individuals with a consistent ciliopathy phenotype (Meckel-Gruber syndrome) and a functional effect of its deficiency on ciliary signaling. Our study also highlighted seven novel candidate genes (TRAPPC3, EXOC3L2, FAM98C, C17orf61, LRRCC1, NEK4, and CELSR2) some of which have established links to ciliogenesis. Finally, we show that the morbid genome of ciliopathies encompasses many founder mutations, the combined carrier frequency of which accounts for a high disease burden in the study population.ConclusionsOur study increases our understanding of the morbid genome of ciliopathies. We also provide the strongest evidence, to date, in support of the classical Mendelian inheritance of Bardet-Biedl syndrome and other ciliopathies.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-016-1099-5) contains supplementary material, which is available to authorized users.

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.