Large-scale next generation sequencing based analysis of SLCO1B1 pharmacogenetics variants in the Saudi population.

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.

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 variants in HOPS complex subunit VPS41 cause cerebellar ataxia and abnormal membrane trafficking.

Membrane trafficking is a complex, essential process in eukaryotic cells responsible for protein transport and processing. Deficiencies in vacuolar protein sorting (VPS) proteins, key regulators of trafficking, cause abnormal intracellular segregation of macromolecules and organelles and are linked to human disease. VPS proteins function as part of complexes such as the homotypic fusion and vacuole protein sorting (HOPS) tethering complex, composed of VPS11, VPS16, VPS18, VPS33A, VPS39 and VPS41. The HOPS-specific subunit VPS41 has been reported to promote viability of dopaminergic neurons in Parkinson's disease but to date has not been linked to human disease. Here, we describe five unrelated families with nine affected individuals, all carrying homozygous variants in VPS41 that we show impact protein function. All affected individuals presented with a progressive neurodevelopmental disorder consisting of cognitive impairment, cerebellar atrophy/hypoplasia, motor dysfunction with ataxia and dystonia, and nystagmus. Zebrafish disease modelling supports the involvement of VPS41 dysfunction in the disorder, indicating lysosomal dysregulation throughout the brain and providing support for cerebellar and microglial abnormalities when vps41 was mutated. This provides the first example of human disease linked to the HOPS-specific subunit VPS41 and suggests the importance of HOPS complex activity for cerebellar function.

Autozygome and high throughput confirmation of disease genes candidacy.

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.

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

This corrects the article DOI: 10.1038/gim.2017.22.

Recessive VARS2 mutation underlies a novel syndrome with epilepsy, mental retardation, short stature, growth hormone deficiency, and hypogonadism.

BACKGROUND: Most mitochondrial and cytoplasmic aminoacyl-tRNA synthetases (aaRSs) are encoded by nuclear genes. Syndromic disorders resulting from mutation of aaRSs genes display significant phenotypic heterogeneity. We expand aaRSs-related phenotypes through characterization of the clinical and molecular basis of a novel autosomal-recessive syndrome manifesting severe mental retardation, ataxia, speech impairment, epilepsy, short stature, microcephaly, hypogonadism, and growth hormone deficiency. RESULTS: A G>A variant in exon 29 of VARS2 (c.3650G>A) (NM_006295) was identified in the index case. This homozygous variant was confirmed by Sanger sequencing and segregated with disease in the family studied. The c.3650G>A change results in alteration of arginine to histidine at residue 1217 (R1217H) of the mature protein and is predicted to be pathogenic. CONCLUSIONS: These findings contribute to a growing list of aaRSs disorders, broadens the spectrum of phenotypes attributable to VARS2 mutations, and provides new insight into genotype-phenotype correlations among the mitochondrial synthetase genes.

Correction to: Recessive VARS2 mutation underlies a novel syndrome with epilepsy, mental retardation, short stature, growth hormone deficiency, and hypogonadism.

CORRECTION: After publication of the article [1], it has been brought to our attention that there is a nomenclature issue with this article. At the time of acceptance, the VARS2 mutation was considered equivalent to the VARS2 mutation. However, this has changed so that VARS now only refers to shorter mitochondrial sequence of valyl-tRNA synthesase containing 1093 amino acids. "Therefore, in the context of this article, every usage of "VARS2" should be replaced with "VARS" when referring to the causative variant".

Identification of a novel genetic locus underlying tremor and dystonia.

BACKGROUND: Five affected individuals with syndromic tremulous dystonia, spasticity, and white matter disease from a consanguineous extended family covering a period of over 24 years are presented. A positional cloning approach utilizing genome-wide linkage, homozygozity mapping and whole exome sequencing was used for genetic characterization. The impact of a calmodulin-binding transcription activator 2, (CAMTA2) isoform 2, hypomorphic mutation on mRNA and protein abundance was studied using fluorescent reporter expression cassettes. Human brain sub-region cDNA libraries were used to study the expression pattern of CAMTA2 transcript variants. RESULTS: Linkage analysis and homozygozity mapping localized the disease allele to a 2.1 Mb interval on chromosome 17 with a LOD score of 4.58. Whole exome sequencing identified a G>A change in the transcript variant 2 5'UTR of CAMTA2 that was only 6 bases upstream of the translation start site (c.-6G > A) (NM_001171166.1) and segregated with disease in an autosomal recessive manner. Transfection of wild type and mutant 5'UTR-linked fluorescent reporters showed no impact upon mRNA levels but a significant reduction in the protein fluorescent activity implying translation inhibition. CONCLUSIONS: Mutation of CAMTA2 resulting in post-transcriptional inhibition of its own gene activity likely underlies a novel syndromic tremulous dystonia.

The landscape of genetic diseases in Saudi Arabia based on the first 1000 diagnostic panels and exomes.

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.

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.

A first-line diagnostic assay for limb-girdle muscular dystrophy and other myopathies.

BACKGROUND: Fifty random genetically unstudied families (limb-girdle muscular dystrophy (LGMD)/myopathy) were screened with a gene panel incorporating 759 OMIM genes associated with neurological disorders. Average coverage of the CDS and 10 bp flanking regions of genes was 99 %. All families were referred to the Neurosciences Clinic of King Faisal Specialist Hospital and Research Centre, Saudi Arabia. Patients presented with muscle weakness affecting the pelvic and shoulder girdle. Muscle biopsy in all cases showed dystrophic or myopathic changes. Our main objective was to evaluate a neurological gene panel as a first-line diagnostic test for LGMD/myopathies. RESULTS: Our panel identified the mutation in 76 % of families (38/50; 11 novel). Thirty-four families had mutations in LGMD-related genes with four others having variants not typically associated with LGMD. The majority of cases had recessive inheritance with homoallelic pathogenic variants (97.4 %, 37/38), as expected considering the high rate of consanguinity in the study population. In one case, we detected a heterozygous mutation in DNAJB responsible for LGMD-1E. Our cohort included seven different subtypes of LGMD2. Mutations of DYSF were the most commonly identified cause of disease followed by that in CAPN3 and FKRP. Non-LGMD myopathies were due to mutations in genes associated with congenital disorder of glycosylation (ALG2), rigid spine muscular dystrophy 1 (SEPN1), inclusion body myopathy2/Nonaka myopathy (GNE), and neuropathy (WNK1). Whole exome sequencing (WES) of patients who remained undiagnosed with the neurological panel did not improve our diagnostic yield. CONCLUSIONS: Our neurological panel achieved a high clinical sensitivity (76 %) and is an effective first-line laboratory test in patients with LGMD and other myopathies. This sensitive, cost-effective, and rapid assay significantly assists clinical practice especially in these phenotypically and genetically heterogeneous disorders. Moreover, the application of the American College of Medical Genetics (ACMG) and Association for Molecular Pathology (AMP) guidelines applied in the classification of variant pathogenecity provides a clear interpretation for physicians on the relevance of such findings.

A Novel Heterozygous De Novo MORC2 Missense Variant Causes an Early Onset and Severe Neurodevelopmental Disorder.

Microrchidia CW-type zinc finger protein 2 (MORC2) is an ATPase-containing nuclear protein which regulates transcription through chromatin remodelling and epigenetic silencing. MORC2 may have a role in the development of neurones, and dominant variants in this gene have recently been linked with disorders including Charcot-Marie-Tooth type 2Z disease, spinal muscular atrophy and, more recently, a neurodevelopmental syndrome consisting of developmental delay, impaired growth, dysmorphic facies, and axonal neuropathy (DIGFAN), presenting with hypotonia, microcephaly, brain atrophy, intellectual disability, hearing loss, faltering growth, and craniofacial dysmorphism. Notably, variants in MORC2 have shown clinical features overlapping with those of Cockayne and Leigh syndromes. Here, we report a case of MORC2-related DIGFAN syndrome in a female infant caused by a novel heterozygous de novo variant. The condition was early onset and severe, further expanding the range of genotypes associated with this disorder. Clinical features included unilateral hearing loss, developmental delay and regression within the first year of life, microcephaly, severe feeding difficulties, and faltering growth, resulting in death at 13 months of age.

The clinical utility of rapid exome sequencing in a consanguineous population.

BACKGROUND: The clinical utility of exome sequencing is now well documented. Rapid exome sequencing (RES) is more resource-intensive than regular exome sequencing and is typically employed in specialized clinical settings wherein urgent molecular diagnosis is thought to influence acute management. Studies on the clinical utility of RES have been largely limited to outbred populations. METHODS: Here, we describe our experience with rapid exome sequencing (RES) in a highly consanguineous population. Clinical settings included intensive care units, prenatal cases approaching the legal cutoff for termination, and urgent transplant decisions. RESULTS: A positive molecular finding (a pathogenic or likely pathogenic variant that explains the phenotype) was observed in 80 of 189 cases (42%), while 15 (8%) and 94 (50%) received ambiguous (variant of uncertain significance (VUS)) and negative results, respectively. The consanguineous nature of the study population gave us an opportunity to observe highly unusual and severe phenotypic expressions of previously reported genes. Clinical utility was observed in nearly all (79/80) cases with positive molecular findings and included management decisions, prognostication, and reproductive counseling. Reproductive counseling is a particularly important utility in this population where the overwhelming majority (86%) of identified variants are autosomal recessive, which are more actionable in this regard than the de novo variants typically reported by RES elsewhere. Indeed, our cost-effectiveness analysis shows compelling cost savings in the study population. CONCLUSIONS: This work expands the diversity of environments in which RES has a demonstrable clinical utility.

Identification of pharmacogenetic variants from large scale next generation sequencing data in the Saudi population.

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.

Genetic Alterations in Pediatric Thyroid Cancer Using a Comprehensive Childhood Cancer Gene Panel.

CONTEXT: Pediatric differentiated thyroid cancer (DTC) differs from adult DTC in its underlying genetics and clinicopathological features. In this report, we studied these aspects in 48 cases of pediatric DTC. PATIENTS AND METHODS: We used the comprehensive Oncomine Childhood Cancer Gene panel on Ion Torrent next-generation sequencing platform. We included 48 patients (37 girls and 11 boys) with pediatric DTC (median age 17 years; range, 5-18 years) and studied the association between these genetic alterations and the clinicopathological features and outcome. RESULTS: Of 48 tumors, 33 (69%) had somatic genetic alterations that were mutually exclusive except in one tumor. BRAFV600E and RET-PTC1 were the most common, occurring in 9 different tumors (19%) each. RET-PTC3 and ETV6-NTRK3 were the next most common, with each occurring in 4 different tumors (8%). Other genetic alterations including EML4-NTRK1, EML4-ALK, NRAS, KRAS, PTEN, and CREBBP occurred once each. There were no differences between those who had mutations and those without mutations with respect to age, sex, tumor multifocality, extrathyroidal extension, vascular invasion, lymph node or distant metastasis, and American Thyroid Association response to therapy status at the last follow-up visits. Similarly, none of these factors was different between those with fusion genes vs single-point mutations vs no mutations. CONCLUSIONS: In pediatric DTC, fusion genes are more common than single-point mutations. The most common genetic alterations are RET-PTC1, BRAFV600E, RET-PTC3, and ETV6-NTRK3. Other alterations occur rarely. Genetic alterations do not correlate with the clinicopathological features or the outcome.

Integrated Analysis of Whole Exome Sequencing and Copy Number Evaluation in Parkinson's Disease.

Genetic studies of the familial forms of Parkinson's disease (PD) have identified a number of causative genes with an established role in its pathogenesis. These genes only explain a fraction of the diagnosed cases. The emergence of Next Generation Sequencing (NGS) expanded the scope of rare variants identification in novel PD related genes. In this study we describe whole exome sequencing (WES) genetic findings of 60 PD patients with 125 variants validated in 51 of these cases. We used strict criteria for variant categorization that generated a list of variants in 20 genes. These variants included loss of function and missense changes in 18 genes that were never previously linked to PD (NOTCH4, BCOR, ITM2B, HRH4, CELSR1, SNAP91, FAM174A, BSN, SPG7, MAGI2, HEPHL1, EPRS, PUM1, CLSTN1, PLCB3, CLSTN3, DNAJB9 and NEFH) and 2 genes that were previously associated with PD (EIF4G1 and ATP13A2). These genes either play a critical role in neuronal function and/or have mouse models with disease related phenotypes. We highlight NOTCH4 as an interesting candidate in which we identified a deleterious truncating and a splice variant in 2 patients. Our combined molecular approach provides a comprehensive strategy applicable for complex genetic disorders.

Whole exome sequencing reveals inherited and de novo variants in autism spectrum disorder: a trio study from Saudi families.

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.