Comprehensive Computational Analysis of GWAS Loci Identifies CCR2 as a Candidate Gene for Celiac Disease Pathogenesis.

Celiac disease (CD) is a gluten intolerance disorder with known genetic contribution. The recent fine mapping and genome-wide association studies (GWAS) have identified up to 57 non-HLA CD susceptibility SNPs, majority of which are non-coding variants lacking any functional annotation. Therefore, we adopted multidimensional computational approach for uncovering the plausible mechanisms through which these GWAS SNPs are connected to CD pathogenesis. At initial phase, we identified that 25 (43.85%) out of 57 CD-SNPs lies in evolutionarily constrained genetic element regions. In follow-up phases, through computational (CADD, GWAVA, and FATHMM algorithms) deleterious intensity measurements, we have discovered that 42 (3.94%) out of 1065 variants (57 CD-lead and 1008-linked SNPs; r2 ≥ 0.8) are differentially deleterious in nature to CD. Further functional scrutinization of these CD variants by public domain eQTL mapping, gene expression, knockout mouse model, and pathway analyses revealed that deleterious SNPs of CCR2 gene influences its expression levels and may also elicit a cascade of T-cell-mediated immunological events leading to intestinal gluten intolerance in genetically susceptible individuals. This study demonstrates the utility of integrated in silico analysis of annotations, gene expression, and pathways in prioritizing the potential complex disease variants from large-scale open source genomic data. J. Cell. Biochem. 118: 2193-2207, 2017. © 2017 Wiley Periodicals, Inc.

Exome Sequencing Identifies the Extremely Rare ITGAV and FN1 Variants in Early Onset Inflammatory Bowel Disease Patients.

Background: Molecular diagnosis of early onset inflammatory bowel disease (IBD) is very important for adopting suitable treatment strategies. Owing to the sparse data available, this study aims to identify the molecular basis of early onset IBD in Arab patients. Methods: A consanguineous Arab family with monozygotic twins presenting early onset IBD was screened by whole exome sequencing (WES). The variants functional characterization was performed by a series of computational biology methods. The IBD variants were further screened in in-house whole exome data of 100 Saudi cohorts ensure their rare prevalence in the population. Results: Genetic screening has identified the digenic autosomal recessive mode of inheritance of ITGAV (G58V) and FN1 (G313V) variants in IBD twins with early onset IBD. Findings from pathogenicity predictions, stability and molecular dynamics have confirmed the deleterious nature of both variants on structural features of the corresponding proteins. Functional biology data suggested that both genes show abundant expression in gastrointestinal tract and immune organs, involved in immune cell restriction, regulation of different immune related pathways. Data from knockout mouse models for ITGAV gene has revealed that the dysregulated expression of this gene impacts intestinal immune homeostasis. The defective ITGAV and FN1 involved in integrin pathway, are likely to induce intestinal inflammation by disturbing immune homeostasis. Conclusions: Our findings provide novel insights into the molecular etiology of pediatric onset IBD and may likely pave way in developing genomic medicine.

Complex Inheritance of Rare Missense Variants in PAK2, TAP2, and PLCL1 Genes in a Consanguineous Arab Family With Multiple Autoimmune Diseases Including Celiac Disease.

Background: Autoimmune diseases (AIDs) share a common molecular etiology and often present overlapping clinical presentations. Thus, this study aims to explore the complex molecular basis of AID by whole exome sequencing and computational biology analysis. Methods: Molecular screening of the consanguineous AID family and the computational biology characterization of the potential variants were performed. The potential variants were searched against the exome data of 100 healthy individuals and 30 celiac disease patients. Result: A complex inheritance pattern of PAK2 (V43A), TAP2 (F468Y), and PLCL1 (V473I) genetic variants was observed in the three probands of the AID family. The PAK2 variant (V43A) is a novel one, but TAP2 (F468Y) and PLCL1 (V473I) variants are extremely rare in local Arab (SGHP and GME) and global (gnomAD) databases. All these variants were localized in functional domains, except for the PAK2 variant (V43A) and were predicted to alter the structural (secondary structure elements, folding, active site confirmation, stability, and solvent accessibility) and functional (gene expression) features. Therefore, it is reasonable to postulate that the dysregulation of PAK2, TAP2, and PLCL1 genes is likely to elicit autoimmune reactions by altering antigen processing and presentation, T cell receptor signaling, and immunodeficiency pathways. Conclusion: Our findings highlight the importance of exploring the alternate inheritance patterns in families presenting complex autoimmune diseases, where classical genetic models often fail to explain their molecular basis. These findings may have potential implications for developing personalized therapies for complex disease patients.

Saudi Familial Hypercholesterolemia Patients With Rare LDLR Stop Gain Variant Showed Variable Clinical Phenotype and Resistance to Multiple Drug Regimen.

Familial hypercholesterolemia (FH), a well-known lipid disease caused by inherited genetic defects in cholesterol uptake and metabolism is underdiagnosed in many countries including Saudi Arabia. The present study aims to identify the molecular basis of severe clinical manifestations of FH patients from unrelated Saudi consanguineous families. Two Saudi families with multiple FH patients fulfilling the combined FH diagnostic criteria of Simon Broome Register, and the Dutch Lipid Clinic Network (DLCN) were recruited. LipidSeq, a targeted resequencing panel for monogenic dyslipidemias, was used to identify causative pathogenic mutation in these two families and in 92 unrelated FH cases. Twelve FH patients from two unrelated families were sharing a very rare, pathogenic and founder LDLR stop gain mutation i.e., c.2027delG (p.Gly676Alafs*33) in both the homozygous or heterozygous states, but not in unrelated patients. Based on the variant zygosity, a marked phenotypic heterogeneity in terms of LDL-C levels, clinical presentations and resistance to anti-lipid treatment regimen (ACE inhibitors, ß-blockers, ezetimibe, statins) of the FH patients was observed. This loss-of-function mutation is predicted to alter the free energy dynamics of the transcribed RNA, leading to its instability. Protein structural mapping has predicted that this non-sense mutation eliminates key functional domains in LDLR, which are essential for the receptor recycling and LDL particle binding. In conclusion, by combining genetics and structural bioinformatics approaches, this study identified and characterized a very rare FH causative LDLR pathogenic variant determining both clinical presentation and resistance to anti-lipid drug treatment.

Whole exome sequencing identifies rare biallelic ALMS1 missense and stop gain mutations in familial Alström syndrome patients.

Alström syndrome (AS, OMIM ID 203800) is a rare childhood multiorgan disorder, which is widely studied in non-Arab ethnic patients. The clinical and molecular basis of AS and the mode of disease inheritance in consanguineous Arab populations is not well investigated. Therefore, to identify the molecular basis of AS in familial forms, the present study performed whole exome sequencing of 5 AS patients belonging to 2 different Bedouin families from Saudi Arabia. The present study identified the AS causative rare biallelic mutations in ALMS gene:T376S in exon 5 and S909* in exon 8 for family A and an R2721* in exon 10 (R2721*) for family B. ALMS1 targeted genetic sequencing of healthy population controls and family members has confirmed its extremely rare frequency and autosomal recessive mode of inheritance. The truncating mutations S909* and R2721* could cause the loss of CC domains and ALMS motif on C-terminal end of the protein and creates unstable protein, which eventually undergoes intracellular degradation. The premature protein truncating mutations described in our study may eventually provide further insight into the functional domains of the ALMS1 protein and contribute to the understanding of the phenotypic spectrum of AS. Whole exome sequencing based molecular diagnosis is expected to rule out ambiguity surrounding clinical diagnosis of suspected AS cases.

Whole exome sequencing of a Saudi family and systems biology analysis identifies CPED1 as a putative causative gene to Celiac Disease.

Celiac disease (CD) is a gastrointestinal disorder whose genetic basis is not fully understood. Therefore, we studied a Saudi family with two CD affected siblings to discover the causal genetic defect. Through whole exome sequencing (WES), we identified that both siblings have inherited an extremely rare and deleterious CPED1 genetic variant (c.241 A > G; p.Thr81Ala) segregating as autosomal recessive mutation, suggesting its putative causal role in the CD. Saudi population specific minor allele frequency (MAF) analysis has confirmed its extremely rare prevalence in homozygous condition (MAF is 0.0004). The Sanger sequencing analysis confirmed the absence of this homozygous variant in 100 sporadic Saudi CD cases. Genotype-Tissue Expression (GTEx) data has revealed that CPED1 is abundantly expressed in gastrointestinal mucosa. By using a combination of systems biology approaches like protein 3D modeling, stability analysis and nucleotide sequence conservation analysis, we have further established that this variant is deleterious to the structural and functional aspects of CPED1 protein. To the best of our knowledge, this variant has not been previously reported in CD or any other gastrointestinal disease. The cell culture and animal model studies could provide further insight into the exact role of CPED1 p.Thr81Ala variant in the pathophysiology of CD. In conclusion, by using WES and systems biology analysis, present study for the first-time reports CPED1 as a potential causative gene for CD in a Saudi family with potential implications to both disease diagnosis and genetic counseling.

Exome sequencing and metabolomic analysis of a chronic kidney disease and hearing loss patient family revealed RMND1 mutation induced sphingolipid metabolism defects.

Mitochondrial disorders (MIDs) shows overlapping clinical presentations owing to the genetic and metabolic defects of mitochondria. However, specific relationship between inherited mutations in nuclear encoded mitochondrial proteins and their functional impacts in terms of metabolic defects in patients is not yet well explored. Therefore, using high throughput whole exome sequencing (WES), we screened a chronic kidney disease (CKD) and sensorineural hearing loss (SNHL) patient, and her family members to ascertain the mode of inheritance of the mutation, and healthy population controls to establish its rare frequency. The impact of mutation on biophysical characteristics of the protein was further studied by mapping it in 3D structure. Furthermore, LC-MS tandem mass spectrophotometry based untargeted metabolomic profiling was done to study the fluctuations in plasma metabolites relevant to disease causative mutations and kidney damage. We identified a very rare homozygous c.631G > A (p.Val211Met) pathogenic mutation in RMND1 gene in the proband, which is inherited in an autosomal recessive fashion. This gene is involved in the mitochondrial translational pathways and contribute in mitochondrial energy metabolism. The p.Val211Met mutation is found to disturb the structural orientation (RMSD is -2.95 Å) and stability (ΔΔG is -0.552 Kcal/mol) of the RMND1 protein. Plasma metabolomics analysis revealed the aberrant accumulation of metabolites connected to lipid and amino acid metabolism pathways. Of these metabolites, pathway networking has discovered ceramide, a metabolite of sphingolipids, which plays a role in different signaling cascades including mitochondrial membrane biosynthesis, is highly elevated in this patient. This study suggests that genetic defects in RMND1 gene alters the mitochondrial energy metabolism leading to the accumulation of ceramide, and subsequently promote dysregulated apoptosis and tissue necrosis in kidneys.

In Silico Approach to Investigate the Structural and Functional Attributes of Familial Hypercholesterolemia Variants Reported in the Saudi Population.

Familial hypercholesterolemia (FH) is a metabolic disorder that leads primarily to premature cardiovascular diseases, the main cause of mortality in Saudi Arabia (SA). FH is underreported and underdiagnosed in SA with statistical evidence of high expected prevalence in such a consanguineous community. Lacking knowledge of which and how these alterations are actually impacting lipid metabolism is one of the main reasons why FH is insufficiently diagnosed in the region. The aim of this study was to develop a fast prediction approach using an integrated bioinformatics method for future screening of the potential causative variants from national registries. A total of 21 variants were detected with majority rate in LDLR (81%). Variants were classified based on the type of mutation. Missense variants resulting in amino acid changes, c.1429G>A (p.D477N), c.1474G>A (p.D492N), c.1731G>T (p.W577C), and c.1783C>T (p.R595W) in LDLR gene, in addition to c.9835A>G (p.S3279G) in APOB, were shown to be deleterious by concordant analysis. Furthermore, functional interaction deformities showed a significant loss and gain of energies in the mutated proteins. These findings will help in distinguishing the most harmful mutations needed to be screened for clinically diagnosed FH patients in SA. Such computational research is necessary to avoid time consumption and the usage of expensive biological experiments. This can be a fast track to facilitate the future filtering and screening of causative mutations from national registries.

Interpreting the Mechanism of APOE (p.Leu167del) Mutation in the Incidence of Familial Hypercholesterolemia; An In-silico Approach.

BACKGROUND: Apolipoprotein E (APOE) gene is a ligand protein in humans which mediates the metabolism of cholesterol by binding to the low-density lipoprotein receptor (LDLR). P.Leu167del mutation in APOE gene was recently connected with Familial Hypercholesterolemia, a condition associated with premature cardiovascular disease. The consequences of this mutation on the protein structure and its receptor binding capacity remain largely unknown. OBJECTIVE: The current study aims to further decipher the underlying mechanism of this mutation using advanced software-based algorithms. The consequences of disrupting the leucine zipper by this mutation was studied at the structural and functional level of the APOE protein. METHODS: 3D protein modeling for both APOE and LDLR (wild types), along with APOE (p.Leu167del) mutant type were generated using homology modeling template-based alignment. Structural deviation analysis was performed to evaluate the spatial orientation and the stability of the mutant APOE structure. Molecular docking analysis simulating APOE-LDLR protein interaction was carried out, in order to evaluate the impact of the mutation on the binding affinity. RESULT: Structural deviation analysis for APOE mutated model showed low degree of deviance scoring root-mean-square deviation, (RMSD) = 0.322 Å. Whereas Docking simulation revealed an enhanced molecular interaction towards the LDLR with an estimation of +171.03 kJ/mol difference in binding free energy. CONCLUSION: This in-silico study suggests that p.Leu167del is causing the protein APOE to associate strongly with its receptor, LDLR. This gain-of-function is likely hindering the ability of LDLR to be effectively recycled back to the surface of the hepatocytes to clear cholesterol from the circulation therefore leading to FH.

The Spectrum of Familial Hypercholesterolemia (FH) in Saudi Arabia: Prime Time for Patient FH Registry.

BACKGROUND: Familial hypercholesterolemia (FH) is a life-threatening inherited condition. Untreated patients have the risk to develop raised plasma levels of cholesterol, atherosclerosis and cardiovascular disease (CVD). If diagnosed and treated early in life, the pathological consequences due to atherosclerosis could be avoided and patients with FH can have an anticipated normal life. Mounting evidence suggests that FH is underdiagnosed and undertreated in all populations. The underlying molecular basis of FH is the presence of mutations in one or more genes in the low-density lipoprotein receptor (LDLR), apolipoprotein B (APOB) or proprotein convertase subtilisin/kexin 9 (PCSK9). However, their prevalence is largely unknown in Saudi Arabia but given the high rates of consanguinity, the prevalence appears to be higher. Furthermore, the high prevalence of obesity and diabetes mellitus in Saudi Arabia increases the vascular disease burden in FH cases by adding additional CVD risk factors. OBJECTIVE: This article explores the spectrum of FH-causing mutations in the highly consanguineous Saudi community, the need for establishing the Saudi FH registry, the challenges in creating gene databases, and cascade screening. CONCLUSION: The establishment of FH registry and genetic testing should raise awareness not only among healthcare professionals, but the general population as well. It also helps to provide the best treatment regimen in a cost effective manner to this under-recognised population of FH patients.

Whole exome sequencing of a consanguineous family identifies the possible modifying effect of a globally rare AK5 allelic variant in celiac disease development among Saudi patients.

Celiac disease (CD), a multi-factorial auto-inflammatory disease of the small intestine, is known to occur in both sporadic and familial forms. Together HLA and Non-HLA genes can explain up to 50% of CD's heritability. In order to discover the missing heritability due to rare variants, we have exome sequenced a consanguineous Saudi family presenting CD in an autosomal recessive (AR) pattern. We have identified a rare homozygous insertion c.1683_1684insATT, in the conserved coding region of AK5 gene that showed classical AR model segregation in this family. Sequence validation of 200 chromosomes each of sporadic CD cases and controls, revealed that this extremely rare (EXac MAF 0.000008) mutation is highly penetrant among general Saudi populations (MAF is 0.62). Genotype and allelic distribution analysis have indicated that this AK5 (c.1683_1684insATT) mutation is negatively selected among patient groups and positively selected in the control group, in whom it may modify the risk against CD development [p<0.002]. Our observation gains additional support from computational analysis which predicted that Iso561 insertion shifts the existing H-bonds between 400th and 556th amino acid residues lying near the functional domain of adenylate kinase. This shuffling of amino acids and their H-bond interactions is likely to disturb the secondary structure orientation of the polypeptide and induces the gain-of-function in nucleoside phosphate kinase activity of AK5, which may eventually down-regulates the reactivity potential of CD4+ T-cells against gluten antigens. Our study underlines the need to have population-specific genome databases to avoid false leads and to identify true candidate causal genes for the familial form of celiac disease.