Explore the vast range of titles available.

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.

Intellectual disability (ID) is a common morbid condition with a wide range of etiologies. The list of monogenic forms of ID has increased rapidly in recent years thanks to the implementation of genomic sequencing techniques. In this study, we describe the phenotypic and genetic findings of 68 families (105 patients) all with novel ID-related variants. In addition to established ID genes, including ones for which we describe unusual mutational mechanism, some of these variants represent the first confirmatory disease–gene links following previous reports ( TRAK1 , GTF3C3 , SPTBN4 and NKX6 - 2 ), some of which were based on single families. Furthermore, we describe novel variants in 14 genes that we propose as novel candidates ( ANKHD1 , ASTN2 , ATP13A1 , FMO4 , MADD , MFSD11 , NCKAP1 , NFASC , PCDHGA10 , PPP1R21 , SLC12A2 , SLK , STK32C and ZFAT ). We highlight MADD and PCDHGA10 as particularly compelling candidates in which we identified biallelic likely deleterious variants in two independent ID families each. We also highlight NCKAP1 as another compelling candidate in a large family with autosomal dominant mild intellectual disability that fully segregates with a heterozygous truncating variant. The candidacy of NCKAP1 is further supported by its biological function, and our demonstration of relevant expression in human brain. Our study expands the locus and allelic heterogeneity of ID and demonstrates the power of positional mapping to reveal unusual mutational mechanisms.

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. 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. We used a data set of ~7,100 patients who underwent genomic testing for various Mendelian disorders to estimate the carrier frequency. 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. 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 Med18 12, 1244–1249.

Background Ciliopathies 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. Results We 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. Conclusions Our 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.

Purpose Establishing links between Mendelian phenotypes and genes enables the proper interpretation of variants therein. Autozygome, a rich source of homozygous variants, has been successfully utilized for the high throughput identification of novel autosomal recessive disease genes. Here, we highlight the utility of the autozygome for the high throughput confirmation of previously published tentative links to diseases. Methods Autozygome and exome analysis of patients with suspected Mendelian phenotypes. All variants were classified according to the American College of Medical Genetics and Genomics guidelines. Results We highlight 30 published candidate genes ( ACTL6B , ADAM22 , AGTPBP1 , APC , C12orf4 , C3orf17 (NEPRO) , CENPF , CNPY3 , COL27A1 , DMBX1 , FUT8 , GOLGA2 , KIAA0556 , LENG8 , MCIDAS , MTMR9 , MYH11 , QRSL1 , RUBCN , SLC25A42 , SLC9A1 , TBXT , TFG , THUMPD1 , TRAF3IP2 , UFC1 , UFM1 , WDR81 , XRCC2 , ZAK ) in which we identified homozygous likely deleterious variants in patients with compatible phenotypes. We also identified homozygous likely deleterious variants in 18 published candidate genes ( ABCA2 , ARL6IP1 , ATP8A2 , CDK9 , CNKSR1 , DGAT1 , DMXL2 , GEMIN4 , HCN2 , HCRT , MYO9A , PARS2 , PLOD3 , PREPL , SCLT1 , STX3 , TXNRD2 , WIPI2 ) although the associated phenotypes are sufficiently different from the original reports that they represent phenotypic expansion or potentially distinct allelic disorders. Conclusions Our results should facilitate the timely relabeling of these candidate disease genes in relevant databases to improve the yield of clinical genomic sequencing.

It is well documented that drug responses are related to Absorption, Distribution, Metabolism, and Excretion (ADME) characteristics of individual patients. Several studies have identified genetic variability in pharmacogenes, that are either directly responsible for or are associated with ADME, giving rise to individualized treatments. Our objective was to provide a comprehensive overview of pharmacogenetic variation in the Saudi population. We mined next generation sequencing (NGS) data from 11,889 unrelated Saudi nationals, to determine the presence and frequencies of known functional SNP variants in 8 clinically relevant pharmacogenes ( CYP2C9 , CYP2C19 , CYP3A5 , CYP4F2 , VKORC1 , DPYD , TPMT and NUDT15 ), recommended by the Clinical Pharmacogenetics Implementation Consortium (CPIC), and collectively identified 82 such star alleles. Functionally significant pharmacogenetic variants were prevalent especially in CYP genes (excluding CYP3A5 ), with 10–44.4% of variants predicted to be inactive or to have decreased activity. In CYP3A5 , inactive alleles (87.5%) were the most common. Only 1.8%, 0.7% and 0.7% of NUDT15 , TPMT and DPYD variants respectively, were predicted to affect gene activity. In contrast, VKORC1 was found functionally, to be highly polymorphic with 53.7% of Saudi individuals harboring variants predicted to result in decreased activity and 31.3% having variants leading to increased metabolic activity. Furthermore, among the 8 pharmacogenes studied, we detected six rare variants with an aggregated frequency of 1.1%, that among several other ethnicities, were uniquely found in Saudi population. Similarly, within our cohort, the 8 pharmacogenes yielded forty-six novel variants predicted to be deleterious. Based upon our findings, 99.2% of individuals from the Saudi population carry at least one actionable pharmacogenetic variant.

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.

Background SLCO1B1 plays an important role in mediating hepatic clearance of many different drugs including statins, angiotensin-converting enzyme inhibitors, chemotherapeutic agents and antibiotics. Several variants in SLCO1B1 have been shown to have a clinically significant impact, in relation to efficacy of these medications. This study provides a comprehensive overview of SLCO1B1 variation in Saudi individuals, one of the largest Arab populations in the Middle East. Methods The dataset of 11,889 (9,961 exomes and 1,928 pharmacogenetic gene panel) Saudi nationals, was used to determine the presence and frequencies of SLCO1B1 variants, as described by the Clinical Pharmacogenetic Implementation Consortium (CPIC). Results We identified 141 previously described SNPs, of which rs2306283 (50%) and rs4149056 (28%), were the most common. In addition, we observed six alleles [ *15 (24.7%) followed by *20 (8.04%), *14 (5.86%), *5 (3.84%), *31 (0.21%) and *9 (0.03%)] predicted to be clinically actionable. Allele diplotype to phenotype conversion revealed 41 OATP1B1 diplotypes. We estimated the burden of rare, and novel predicted deleterious variants, resulting from 17 such alterations. Conclusions The data we present, from one of the largest Arab cohorts studied to date, provides the most comprehensive overview of SLCO1B1 variants, and the subsequent OATP1B1 activity of this ethnic group, which thus far remains relatively underrepresented in available international genomic databases. We believe that the presented data provides a basis for further clinical investigations and the application of personalized statin drug therapy guidance in Arabs.

Severe combined immunodeficiency disease (SCID) is the most severe form of primary immunodeficiency disorders (PID). T-cell receptor excision circle (TREC) copy number analysis is an efficient tool for population-based newborn screening (NBS) for SCID and other T cell lymphopenias. We sought to assess the incidence of SCID among Saudi newborn population and examine the feasibility of using targeted next generation sequencing PID gene panel (T-NGS PID) on DNA isolated from dried blood spots (DBSs) in routine NBS programs as a mutation screening tool for samples with low TREC count. Punches from 8,718 DBS collected on Guthrie cards were processed anonymously for the TREC assay. DNA was extracted from samples with confirmed low TREC count, then screened for 22q11.2 deletion syndrome by real-time polymerase chain reaction and for mutations in PID-related genes by T-NGS PID panel. Detected mutations were confirmed by Sanger sequencing. Sixteen out of the 8,718 samples were confirmed to have low TREC copy number. Autosomal recessive mutations in were confirmed in three samples. Two additional samples were positive for the 22q11.2 deletion syndrome. In this study, we provide evidence for high incidence of SCID among Saudi population (1/2,906 live births) and demonstrate the feasibility of using T-NGS PID panel on DNA extracted from DBSs as a new reliable, rapid, and cost-effective mutation screening method for newborns with low TREC assay, which can be implemented as part of NBS programs for SCID.

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder with genetic and clinical heterogeneity. The interplay of de novo and inherited rare variants has been suspected in the development of ASD. Here, we applied whole exome sequencing (WES) on 19 trios from singleton Saudi families with ASD. We developed an analysis pipeline that allows capturing both de novo and inherited rare variants predicted to be deleterious. A total of 47 unique rare variants were detected in 17 trios including 38 which are newly discovered. The majority were either autosomal recessive or X-linked. Our pipeline uncovered variants in 15 ASD-candidate genes, including 5 ( GLT8D1 , HTATSF1 , OR6C65 , ITIH6 and DDX26B ) that have not been reported in any human condition. The remaining variants occurred in genes formerly associated with ASD or other neurological disorders. Examples include SUMF1 , KDM5B and MXRA5 (Known-ASD genes), PRODH2 and KCTD21 (implicated in schizophrenia), as well as USP9X and SMS (implicated in intellectual disability). Consistent with expectation and previous studies, most of the genes implicated herein are enriched for biological processes pertaining to neuronal function. Our findings underscore the private and heterogeneous nature of the genetic architecture of ASD even in a population with high consanguinity rates.