Point of Care Exome Sequencing Reveals Allelic and Phenotypic Heterogeneity Underlying Mendelian disease in Qatar.

The effectiveness of next generation sequencing at solving genetic disease has motivated the rapid adoption of this technology into clinical practice around the world. In this study, we use whole exome sequencing (WES) to assess 48 patients with Mendelian disease from 30 serial families as part of the "Qatar Mendelian Disease pilot program" - a coordinated multi-center effort to build capacity and clinical expertise in genetic medicine in Qatar. By enrolling whole families (parents plus available siblings), we demonstrate significantly improved discriminatory power for candidate variant identification over trios for both de novo and recessive inheritance patterns. For the same index cases, we further demonstrate that even in the absence of families, variant prioritization is improved up to 8-fold when a modest set of population-matched controls is used vs large public databases, stressing the poor representation of Middle Eastern alleles in presently available databases. Our in-house pipeline identified candidate disease variants in 27 of 30 families (90%), 23 of which (85%) harbor novel pathogenic variants in known disease genes, pointing to significant allelic heterogeneity and founder mutations underlying Mendelian disease in the Middle East. For 6 of these families, the clinical presentation was only partially explained by the candidate gene, suggesting phenotypic expansion of known syndromes. Our pilot study demonstrates the utility of WES for Middle Eastern populations, the dramatic improvement in variant prioritization conferred by enrolling population-matched controls and/or enrolling additional unaffected siblings at the point-of-care, and 25 novel disease-causing alleles, relevant to newborn and premarital screening panels in regional populations.

A Novel STK4 Mutation Impairs T Cell Immunity Through Dysregulation of Cytokine-Induced Adhesion and Chemotaxis Genes.

PURPOSE: Human serine/threonine kinase 4 (STK4) deficiency is a rare, autosomal recessive genetic disorder leading to combined immunodeficiency; however, the extent to which immune signaling and host defense are impaired is unclear. We assessed the functional consequences of a novel, homozygous nonsense STK4 mutation (NM_006282.2:c.871C > T, p.Arg291*) identified in a pediatric patient by comparing his innate and adaptive cell-mediated and humoral immune responses with those of three heterozygous relatives and unrelated controls. METHODS: The genetic etiology was verified by whole genome and Sanger sequencing. STK4 gene and protein expression was measured by quantitative RT-PCR and immunoblotting, respectively. Cellular abnormalities were assessed by high-throughput RT-RCR, RNA-Seq, ELISA, and flow cytometry. Antibody responses were assessed by ELISA and phage immunoprecipitation-sequencing. RESULTS: The patient exhibited partial loss of STK4 expression and complete loss of STK4 function combined with recurrent viral and bacterial infections, notably persistent Epstein-Barr virus viremia and pulmonary tuberculosis. Cellular and molecular analyses revealed abnormal fractions of T cell subsets, plasmacytoid dendritic cells, and NK cells. The transcriptional responses of the patient's whole blood and PBMC samples indicated dysregulated interferon signaling, impaired T cell immunity, and increased T cell apoptosis as well as impaired regulation of cytokine-induced adhesion and leukocyte chemotaxis genes. Nonetheless, the patient had detectable vaccine-specific antibodies and IgG responses to various pathogens, consistent with a normal CD19 + B cell fraction, albeit with a distinctive antibody repertoire, largely driven by herpes virus antigens. CONCLUSION: Patients with STK4 deficiency can exhibit broad impairment of immune function extending beyond lymphoid cells.

Point-of-care whole-exome sequencing of idiopathic male infertility.

PURPOSE: Nonobstructive azoospermia (NOA) affects 1% of the male population; however, despite state-of-the-art clinical assessment, for most patients the cause is unknown. We capitalized on an analysis of multiplex families in the Middle East to identify highly penetrant genetic causes. METHODS: We used whole-exome sequencing (WES) in 8 consanguineous families and combined newly discovered genes with previously reported ones to create a NOA gene panel, which was used to identify additional variants in 75 unrelated idiopathic NOA subjects and 74 fertile controls. RESULTS: In five of eight families, we identified rare deleterious recessive variants in CCDC155, NANOS2, SPO11, TEX14, and WNK3 segregating with disease. These genes, which are novel to human NOA, have remarkable testis-specific expression, and murine functional evidence supports roles for them in spermatogenesis. Among 75 unrelated NOA subjects, we identified 4 (~5.3%) with additional recessive variants in these newly discovered genes and 6 with deleterious variants in previously reported NOA genes, yielding an overall genetic etiology for 13.3% subjects versus 0 fertile controls (p = 0.001). CONCLUSION: NOA affects millions of men, many of whom remain idiopathic despite extensive laboratory evaluation. The genetic etiology for a substantial fraction of these patients (>50% familial and >10% sporadic) may be discovered by WES at the point of care.

Two hits in one: whole genome sequencing unveils LIG4 syndrome and urofacial syndrome in a case report of a child with complex phenotype.

BACKGROUND: Ligase IV syndrome, a hereditary disease associated with compromised DNA damage response mechanisms, and Urofacial syndrome, caused by an impairment of neural cell signaling, are both rare genetic disorders, whose reports in literature are limited. We describe the first case combining both disorders in a specific phenotype. CASE PRESENTATION: We report a case of a 7-year old girl presenting with a complex phenotype characterized by multiple congenital abnormalities and dysmorphic features, microcephaly, short stature, combined immunodeficiency and severe vesicoureteral reflux. Whole Genome Sequencing was performed and a novel ligase IV homozygous missense c.T1312C/p.Y438H mutation was detected, and is believed to be responsible for most of the clinical features of the child, except vesicoureteral reflux which has not been previously described for ligase IV deficiency. However, we observed a second rare damaging (nonsense) homozygous mutation (c.C2125T/p.R709X) in the leucine-rich repeats and immunoglobulin-like domains 2 gene that encodes a protein implicated in neural cell signaling and oncogenesis. Interestingly, this mutation has recently been reported as pathogenic and causing urofacial syndrome, typically displaying vesicoureteral reflux. Thus, this second mutation completes the missing genetic explanation for this intriguing clinical puzzle. We verified that both mutations fit an autosomal recessive inheritance model due to extensive consanguinity. CONCLUSIONS: We successfully identified a novel ligase IV mutation, causing ligase IV syndrome, and an additional rare leucine-rich repeats and immunoglobulin-like domains 2 gene nonsense mutation, in the context of multiple autosomal recessive conditions due to extensive consanguinity. This work demonstrates the utility of Whole Genome Sequencing data in clinical diagnosis in such cases where the combination of multiple rare phenotypes results in very intricate clinical pictures. It also reports a novel causative mutation and a clinical phenotype, which will help in better defining the essential features of both ligase IV and leucine-rich repeats and immunoglobulin-like domains 2 deficiency syndromes.

Predicted deleterious variants in the human genome relevant to gene therapy with adeno-associated virus vectors.

Based on the observation that humans have variable responses of gene expression with the same dose of an adeno-associated vector, we hypothesized that there are deleterious variants in genes coding for processes required for adeno-associated virus (AAV)-mediated gene transfer/expression that may hamper or enhance the effectiveness of AAV-mediated gene therapy. To assess this hypothesis, we evaluated 69,442 whole genome sequences from three populations (European, African/African American, and Qatari) for predicted deleterious variants in 62 genes known to play a role in AAV-mediated gene transfer/expression. The analysis identified 5,564 potentially deleterious mutations of which 27 were classified as common based on an allele frequency ≥1% in at least one population studied. Many of these deleterious variants are predicated to prevent while others enhance effective AAV gene transfer/expression, and several are linked to known hereditary disorders. The data support the hypothesis that, like other drugs, human genetic variability contributes to the person-to-person effectiveness of AAV gene therapy and the screening for genetic variability should be considered as part of future clinical trials.

Metabolic and Metabo-Clinical Signatures of Type 2 Diabetes, Obesity, Retinopathy, and Dyslipidemia.

Macro- and microvascular complications of type 2 diabetes (T2D), obesity, and dyslipidemia share common metabolic pathways. In this study, using a total of 1,300 metabolites from 996 Qatari adults (57% with T2D) and 1,159 metabolites from an independent cohort of 2,618 individuals from the Qatar BioBank (11% with T2D), we identified 373 metabolites associated with T2D, obesity, retinopathy, dyslipidemia, and lipoprotein levels, 161 of which were novel. Novel metabolites included phospholipids, sphingolipids, lysolipids, fatty acids, dipeptides, and metabolites of the urea cycle and xanthine, steroid, and glutathione metabolism. The identified metabolites enrich pathways of oxidative stress, lipotoxicity, glucotoxicity, and proteolysis. Second, we identified 15 patterns we defined as "metabo-clinical signatures." These are clusters of patients with T2D who group together based on metabolite levels and reveal the same clustering in two or more clinical variables (obesity, LDL, HDL, triglycerides, and retinopathy). These signatures revealed metabolic pathways associated with different clinical patterns and identified patients with extreme (very high/low) clinical variables associated with extreme metabolite levels in specific pathways. Among our novel findings are the role of N-acetylmethionine in retinopathy in conjunction with dyslipidemia and the possible roles of N-acetylvaline and pyroglutamine in association with high cholesterol levels and kidney function.

Novel MYO5B mutation in microvillous inclusion disease of Syrian ancestry.

Microvillus inclusion disease (MVID) is a rare autosomal recessive condition characterized by a lack of microvilli on the surface of enterocytes, resulting in severe, life-threatening diarrhea that could lead to mortality within the first year of life. We identify two unrelated families, each with one child presenting with severe MVID from birth. Using trio whole-exome sequencing, we observed that the two families share a novel nonsense variant (Glu1589*) in the MYO5B gene, a type Vb myosin motor protein in which rare damaging mutations were previously described to cause MVID. This founder mutation was very rare in public databases and is likely specific to patients of Syrian ancestry. We present a detailed account of both patients' clinical histories to fully characterize the effect of this variant and expand the genotype-phenotype databases for MVID patients from the Middle East.

Whole-exome sequencing identifies common and rare variant metabolic QTLs in a Middle Eastern population.

Metabolomics-genome-wide association studies (mGWAS) have uncovered many metabolic quantitative trait loci (mQTLs) influencing human metabolic individuality, though predominantly in European cohorts. By combining whole-exome sequencing with a high-resolution metabolomics profiling for a highly consanguineous Middle Eastern population, we discover 21 common variant and 12 functional rare variant mQTLs, of which 45% are novel altogether. We fine-map 10 common variant mQTLs to new metabolite ratio associations, and 11 common variant mQTLs to putative protein-altering variants. This is the first work to report common and rare variant mQTLs linked to diseases and/or pharmacological targets in a consanguineous Arab cohort, with wide implications for precision medicine in the Middle East.

Exome sequencing-based identification of novel type 2 diabetes risk allele loci in the Qatari population.

BACKGROUND: Type 2 diabetes (T2D) susceptibility is influenced by genetic and lifestyle factors. To date, the majority of genetic studies of T2D have been in populations of European and Asian descent. The focus of this study is on genetic variations underlying T2D in Qataris, a population with one of the highest incidences of T2D worldwide. RESULTS: Illumina HiSeq exome sequencing was performed on 864 Qatari subjects (574 T2D cases, 290 controls). Sequence kernel association test (SKAT) gene-based analysis identified an association for low frequency potentially deleterious variants in 6 genes. However, these findings were not replicated by SKAT analysis in an independent cohort of 12,699 exomes, primarly due to the absence of low frequency potentially deleterious variants in 5 of the 6 genes. Interestingly one of the genes identified, catenin beta 1 (CTNNB1, ß-catenin), is the key effector of the Wnt pathway and interacts with the nuclear receptor transcription factor 7-like 2 (TCF7L2), variants which are the most strongly associated with risk of developing T2D worldwide. Single variant analysis did not identify any associated variants, suggesting the SKAT association signal was not driven by individual variants. None of the 6 associated genes were among 634 previously described T2D genes. CONCLUSIONS: The observation that genes not previously linked to T2D in prior studies of European and Asian populations are associated with T2D in Qatar provides new insights into the complexity of T2D pathogenesis and emphasizes the importance of understudied populations when assessing genetic variation in the pathogenesis of common disorders.

The Qatar genome: a population-specific tool for precision medicine in the Middle East.

Reaching the full potential of precision medicine depends on the quality of personalized genome interpretation. In order to facilitate precision medicine in regions of the Middle East and North Africa (MENA), a population-specific genome for the indigenous Arab population of Qatar (QTRG) was constructed by incorporating allele frequency data from sequencing of 1,161 Qataris, representing 0.4% of the population. A total of 20.9 million single nucleotide polymorphisms (SNPs) and 3.1 million indels were observed in Qatar, including an average of 1.79% novel variants per individual genome. Replacement of the GRCh37 standard reference with QTRG in a best practices genome analysis workflow resulted in an average of 7* deeper coverage depth (an improvement of 23%) and 756,671 fewer variants on average, a reduction of 16% that is attributed to common Qatari alleles being present in QTRG. The benefit for using QTRG varies across ancestries, a factor that should be taken into consideration when selecting an appropriate reference for analysis.

Type 2 Diabetes Risk Allele Loci in the Qatari Population.

BACKGROUND: The prevalence of type 2 diabetes (T2D) is increasing in the Middle East. However, the genetic risk factors for T2D in the Middle Eastern populations are not known, as the majority of studies of genetic risk for T2D are in Europeans and Asians. METHODS: All subjects were ≥3 generation Qataris. Cases with T2D (n = 1,124) and controls (n = 590) were randomly recruited and assigned to the 3 known Qatari genetic subpopulations [Bedouin (Q1), Persian/South Asian (Q2) and African (Q3)]. Subjects underwent genotyping for 37 single nucleotide polymorphisms (SNPs) in 29 genes known to be associated with T2D in Europeans and/or Asian populations, and an additional 27 tag SNPs related to these susceptibility loci. Pre-study power analysis suggested that with the known incidence of T2D in adult Qataris (22%), the study population size would be sufficient to detect significant differences if the SNPs were risk factors among Qataris, assuming that the odds ratio (OR) for T2D SNPs in Qatari's is greater than or equal to the SNP with highest known OR in other populations. RESULTS: Haplotype analysis demonstrated that Qatari haplotypes in the region of known T2D risk alleles in Q1 and Q2 genetic subpopulations were similar to European haplotypes. After Benjamini-Hochberg adjustment for multiple testing, only two SNPs (rs7903146 and rs4506565), both associated with transcription factor 7-like 2 (TCF7L2), achieved statistical significance in the whole study population. When T2D subjects and control subjects were assigned to the known 3 Qatari subpopulations, and analyzed individually and with the Q1 and Q2 genetic subpopulations combined, one of these SNPs (rs4506565) was also significant in the admixed group. No other SNPs associated with T2D in all Qataris or individual genetic subpopulations. CONCLUSIONS: With the caveats of the power analysis, the European/Asian T2D SNPs do not contribute significantly to the high prevalence of T2D in the Qatari population, suggesting that the genetic risks for T2D are likely different in Qataris compared to Europeans and Asians.

Correction: Type 2 Diabetes Risk Allele Loci in the Qatari Population.

[This corrects the article DOI: 10.1371/journal.pone.0156834.].