Replication of GWAS Coding SNPs Implicates MMEL1 as a Potential Susceptibility Locus among Saudi Arabian Celiac Disease Patients

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.

Genome-Wide Association Study-Guided Exome Rare Variant Burden Analysis Identifies IL1R1 and CD3E as Potential Autoimmunity Risk Genes for Celiac Disease

Celiac disease (CeD) is a multifactorial autoimmune enteropathy characterized by the overactivation of the immune system in response to dietary gluten. The molecular etiology of CeD is still not well-understood. Therefore, this study aims to identify potential candidate genes involved in CeD pathogenesis by applying multilayered system biology approaches. Initially, we identified rare coding variants shared between the affected siblings in two rare Arab CeD families by whole-exome sequencing (WES). Then we used the STRING database to construct a protein network of rare variants and genome-wide association study (GWAS) loci to explore their molecular interactions in CeD. Furthermore, the hub genes identified based on network topology parameters were subjected to a series of computational validation analyses like pathway enrichment, gene expression, knockout mouse model, and variant pathogenicity predictions. Our findings have shown the absence of rare variants showing classical Mendelian inheritance in both families. However, interactome analysis of rare WES variants and GWAS loci has identified a total of 11 hub genes. The multidimensional computational analysis of hub genes has prioritized IL1R1 for family A and CD3E for family B as potential genes. These genes were connected to CeD pathogenesis pathways of T-cell selection, cytokine signaling, and adaptive immune response. Future multi-omics studies may uncover the roles of IL1R1 and CD3E in gluten sensitivity. The present investigation lays forth a novel approach integrating next-generation sequencing (NGS) of familial cases, GWAS, and computational analysis for solving the complex genetic architecture of CeD.

Investigating the association of polymorphisms of ANKRD55 and MMEL1 with susceptibility to multiple sclerosis in Iranian population

OBJECTIVE

Multiple sclerosis (MS) is an autoimmune neurological disability in which immune cells attack the myelin sheaths that protect nerve fibers. The pathogenesis of the disease involves both complex genetic effects as well as multifaceted gene-environment interactions. In the present study, we examined the association of two Single nucleotide polymorphisms (SNPs) in ANKRD55 (rs6859219) and MMEL (rs3748816) with MS in the Iranian population. ANKRD55 is specifically expressed in human peripheral blood mononuclear cells and CD4 + T cells, while MMEL1is involved in the degradation of both neuropeptides and β-amyloid.

METHODS

In this case-control study, 110 patients with MS and 110 matched healthy controls were enrolled. The Participants were genotyped for ANKRD55 and MMEL1 SNPs using PCR-RFLP and Real-Time TaqMan SNP Genotyping respectively. The results were finally analyzed using SPSS software version 22.

RESULTS

Our results did not show significant differences in allelic frequencies of two SNPs among cases and controls (P-Value >0.05). However, for ANKRD55 (rs6859219), CA genotype was shown to have a protective effect (p = 0.035 and OR = 0.55), while CC genotype was a susceptive genotype to MS (p = 0.036 and OR = 1.8). There was no significant difference in genotypic frequencies of SNP rs3748816 in MMEL1.

CONCLUSION

We could successfully replicate the association of ANKRD55 (rs6859219) with susceptibility to MS in the Iranian population. Our result can provide an insight into better understanding the pathogenesis of MS and also improve the genetic counseling for patients affected with multiple sclerosis in Iran.

Celiac disease susceptibility: The genome and beyond

Celiac Disease (CeD) is an immune-mediated complex disease that is triggered by the ingestion of gluten and develops in genetically susceptible individuals. It has been known for a long time that the Human Leucocyte Antigen (HLA) molecules DQ2 and DQ8 are necessary, although not sufficient, for the disease development, and therefore other susceptibility genes and (epi)genetic events must participate in CeD pathogenesis. The advances in Genomics during the last 15 years have made CeD one of the immune-related disorders with the best-characterized genetic component. In the present work, we will first review the main Genome-Wide Association Studies (GWAS) carried out in the disorder, and emphasize post-GWAS discoveries, including diverse integrative strategies, SNP prioritization approaches, and insights into the Microbiome through the host Genomics. Second, we will explore CeD-related Epigenetics and Epigenomics, mostly focusing on the emerging knowledge of the celiac methylome, and the vast but yet under-explored non-coding RNA (ncRNA) landscape. We conclude that much has been done in the field although there are still completely unvisited areas in the post-Genomics of CeD. Chromatin conformation and accessibility, and Epitranscriptomics are promising domains that need to be unveiled to complete the big picture of the celiac Genome.

Exploring celiac disease candidate pathways by global gene expression profiling and gene network cluster analysis

Celiac disease (CeD) is a gastrointestinal autoimmune disorder, whose specific molecular basis is not yet fully interpreted. Therefore, in this study, we compared the global gene expression profile of duodenum tissues from CeD patients, both at the time of disease diagnosis and after two years of the gluten-free diet. A series of advanced systems biology approaches like differential gene expression, protein–protein interactions, gene network-cluster analysis were deployed to annotate the candidate pathways relevant to CeD pathogenesis. The duodenum tissues from CeD patients revealed the differential expression of 106 up- and 193 down-regulated genes. The pathway enrichment of differentially expressed genes (DEGs) highlights the involvement of biological pathways related to loss of cell division regulation (cell cycle, p53 signalling pathway), immune system processes (NOD-like receptor signalling pathway, Th1, and Th2 cell differentiation, IL-17 signalling pathway) and impaired metabolism and absorption (mineral and vitamin absorptions and drug metabolism) in celiac disease. The molecular dysfunctions of these 3 biological events tend to increase the number of intraepithelial lymphocytes (IELs) and villous atrophy of the duodenal mucosa promoting the development of CeD. For the first time, this study highlights the involvement of aberrant cell division, immune system, absorption, and metabolism pathways in CeD pathophysiology and presents potential novel therapeutic opportunities.

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.

The genetic association study of TP53 polymorphisms in Saudi obese patients

Obesity is a multifactorial metabolic disorder characterized by low grade chronic inflammation. Rare and novel mutations in genes which are vital in several key pathways have been reported to alter the energy expenditure which regulates body weight. The TP53 or p53 gene plays a prominent role in regulating various metabolic activities such as glycolysis, lipolysis, and glycogen synthesis. Recent genome-wide association studies reported that tumor suppressor gene p53 variants play a critical role in the predisposition of type 2 diabetes and obesity. Till date, no reports are available from the Arabian population; hence the present study was intended to assess the association between p53 variants with risk of obesity development in the Saudi population. We have selected three p53 polymorphisms, rs1642785 (C > G), and rs9894946 (A > G), and rs1042522 (Pro72Arg; C > G) and assessed their association with obesity risk in the Saudi population. Phenotypic and biochemical parameters were also evaluated to check their association with p53 genotypes and obesity. Genotyping was carried out on 136 obese and 122 normal samples. We observed that there is significantly increased prevalence p52 Pro72Arg (rs1042522) polymorphism in obese persons when compared to controls at GG genotype in overall comparison (OR: 2.169, 95% CI: 1.086-4.334, p = 0.02716). Male obese subjects showed three-fold higher risk at GG genotype (OR: 3.275, 95% CI: 1.230-8.716, p = 0.01560) and two-fold risk at G allele (OR: 1.827, 95% CI: 1.128-2.958, p = 0.01388) of p53 variant Pro72Arg respectively. This variant has also shown significant influence on cholesterol, LDL level, and random insulin levels in obese subjects (p ≤ 0.05). In conclusion, p53 Pro72Arg variant is highly prevalent among obese individuals and may act as a genetic modifier for obesity development among Saudis.

Health related quality of life among Saudi children and adolescents with celiac disease

BACKGROUND/AIMS

Celiac disease (CD) is an immune-mediated enteropathy triggered by gluten ingestion in genetically susceptible individuals. This study reports on the quality of life (QOL) of Saudi children and adolescents with CD.

PATIENTS AND METHODS

This is a case control study that included Saudi patients with CD, aged 9-18 years, who attended CD Clinic at KAU between February 2017 and July 2018. The study was conducted using the Short-Form (SF-36) questionnaire for all candidates, CD-specific QOL questionnaire for the CD patients group, and CD screen questionnaire for the control group.

RESULTS

Overall, 354 subjects were studied (111 CD patients and 243 control). Female subjects constituted 56.8% of both patient and control groups. In the generic SF-36 questionnaire, QOL was comparable between patients and controls in all domains except for the general health domain, which showed difference in favor of the controls (55.01 ± 26.41 and 62.96 ± 18.16, P = 0.005). We also found that males have lower QOL scores in the domains comprising health change (P = 0.02), physical functioning (P = 0.04, role functioning/emotional and emotional well-being (P = 0.049). The CD-specific QOL showed excellent and good scores for 79.3% of patients overall in the tested domains. Patients with poor adherence showed poor QOL in all generic (SF-36) domains but not in CD-specific domains.

CONCLUSION

The HRQOL for Saudi CD children on GFD is generally comparable to the healthy control with exception of the general health domain. Adherence to the GFD improves the generic (SF-36) QOL domains.

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 400^th and 556^th 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.

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; r²  ≥ 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.