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

Protein structural insights into a rare PCSK9 gain-of-function variant (R496W) causing familial hypercholesterolemia in a Saudi family: whole exome sequencing and computational analysis

Familial hypercholesterolemia (FH) is a globally underdiagnosed genetic condition associated with premature cardiovascular death. The genetic etiology data on Arab FH patients is scarce. Therefore, this study aimed to identify the genetic basis of FH in a Saudi family using whole exome sequencing (WES) and multidimensional bioinformatic analysis. Our WES findings revealed a rare heterozygous gain-of-function variant (R496W) in the exon 9 of the PCSK9 gene as a causal factor for FH in this family. This variant was absent in healthy relatives of the proband and 200 healthy normolipidemic controls from Saudi Arabia. Furthermore, this variant has not been previously reported in various regional and global population genomic variant databases. Interestingly, this variant is classified as "likely pathogenic" (PP5) based on the variant interpretation guidelines of the American College of Medical Genetics (ACMG). Computational functional characterization suggested that this variant could destabilize the native PCSK9 protein and alter its secondary and tertiary structural features. In addition, this variant was predicted to negatively influence its ligand-binding ability with LDLR and Alirocumab antibody molecules. This rare PCSK9 (R496W) variant is likely to expand our understanding of the genetic basis of FH in Saudi Arabia. This study also provides computational structural insights into the genotype-protein phenotype relationship of PCSK9 pathogenic variants and contributes to the development of personalized medicine for FH patients in the future.

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.

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.

Identification of new susceptibility loci associated with rheumatoid arthritis

OBJECTIVES

The present study aimed to discover novel susceptibility loci associated with risk of rheumatoid arthritis (RA).

METHODS

We performed a new genome-wide association study (GWAS) in Chinese subjects (1027 RA cases and 2879 controls) and further conducted an expanded meta-analysis with previous GWAS summary data and replication studies. The functional roles of the associated loci were interrogated using publicly available databases. Dual-luciferase reporter and cytokine assay were also used for exploring variant function.

RESULTS

We identified five new susceptibility loci (IL12RB2, BOLL-PLCL1, CCR2, TCF7 and IQGAP1; p_meta <5.00E-08) with same effect direction in each study cohort. The sensitivity analyses showed that the genetic association of at least three loci was reliable and robust. All these lead variants are expression quantitative trait loci and overlapped with epigenetic marks in immune cells. Furthermore, genes within the five loci are genetically associated with risk of other autoimmune diseases, and genes within four loci are known functional players in autoimmunity, which supports the validity of our findings. The reporter assay showed that the risk allele of rs8030390 in IQGAP1 have significantly increased reporter activity in HEK293T cells. In addition, the cytokine assay found that the risk allele of rs244672 in TCF7 was most significantly associated with increased plasma IL-17A levels in healthy controls. Finally, identified likely causal genes in these loci significantly interacted with RA drug targets.

CONCLUSION

This study identified novel RA risk loci and highlighted that comprehensive genetic study can provide important information for RA pathogenesis and drug therapy.

Molecular modelling and dynamic simulations of sequestosome 1 (SQSTM1) missense mutations linked to Paget disease of bone

The Paget disease (PDB; OMIM is 167250) is a chronic bone disease caused by pathogenic mutations in Sequestome1/p62 (SQSTM1) gene. This study has aimed to interpret the relationship of PDB linked SQSTM1 mutations with protein structure and its molecular dynamic features. The disease causative missense mutations were initially collected, and then analyzed for their, exonic and domain distribution, impact on secondary and tertiary structures, and their ability on protein-ligand interactions, using a combination of systems biology approaches. Our results show that most PDB linked SQSTM1 missense mutations affect amino acid residues clustered within or near the UBA domain (aa 389-434), which participates in the ubiquitination of substrates. We also report that the majority mutations occurred in α-helices over β-strands but their effects on the secondary structure were mostly neutral. Global tertiary structure deviations were minimal; however, at amino acid residue level minor structural changes were evident. The molecular dynamics simulation analysis showed that both PB1 and UBA domains were under constant structural fluctuations resulting in closed form conformation of SQSMT1 protein structure, when it is bound to PRKCI ligand. We also found salt bridge conformation changes in the UBA domain of SQSTM1 mutants when they bound to the PRKCI interactor protein. This finding suggests the possibility that mutations in SQSTM1 could impair its ability to ubiquitinate the substrates, eventually affecting autophagy and apoptosis, especially in mature osteoclasts. This study presents the additional insight into structure and function relationship between SQSTM1 mutations and PDB pathogenesis. Communicated by Ramaswamy H. Sarma.

pCADD: SNV prioritisation in Sus scrofa

BACKGROUND

In animal breeding, identification of causative genetic variants is of major importance and high economical value. Usually, the number of candidate variants exceeds the number of variants that can be validated. One way of prioritizing probable candidates is by evaluating their potential to have a deleterious effect, e.g. by predicting their consequence. Due to experimental difficulties to evaluate variants that do not cause an amino-acid substitution, other prioritization methods are needed. For human genomes, the prediction of deleterious genomic variants has taken a step forward with the introduction of the combined annotation dependent depletion (CADD) method. In theory, this approach can be applied to any species. Here, we present pCADD (p for pig), a model to score single nucleotide variants (SNVs) in pig genomes.

RESULTS

To evaluate whether pCADD captures sites with biological meaning, we used transcripts from miRNAs and introns, sequences from genes that are specific for a particular tissue, and the different sites of codons, to test how well pCADD scores differentiate between functional and non-functional elements. Furthermore, we conducted an assessment of examples of non-coding and coding SNVs, which are causal for changes in phenotypes. Our results show that pCADD scores discriminate between functional and non-functional sequences and prioritize functional SNVs, and that pCADD is able to score the different positions in a codon relative to their redundancy. Taken together, these results indicate that based on pCADD scores, regions with biological relevance can be identified and distinguished according to their rate of adaptation.

CONCLUSIONS

We present the ability of pCADD to prioritize SNVs in the pig genome with respect to their putative deleteriousness, in accordance to the biological significance of the region in which they are located. We created scores for all possible SNVs, coding and non-coding, for all autosomes and the X chromosome of the pig reference sequence Sscrofa11.1, proposing a toolbox to prioritize variants and evaluate sequences to highlight new sites of interest to explain biological functions that are relevant to animal breeding.

Identification of key regulatory genes connected to NF-κB family of proteins in visceral adipose tissues using gene expression and weighted protein interaction network

Obesity is connected to the activation of chronic inflammatory pathways in both adipocytes and macrophages located in adipose tissues. The nuclear factor (NF)-κB is a central molecule involved in inflammatory pathways linked to the pathology of different complex metabolic disorders. Investigating the gene expression data in the adipose tissue would potentially unravel disease relevant gene interactions. The present study is aimed at creating a signature molecular network and at prioritizing the potential biomarkers interacting with NF-κB family of proteins in obesity using system biology approaches. The dataset GSE88837 associated with obesity was downloaded from Gene Expression Omnibus (GEO) database. Statistical analysis represented the differential expression of a total of 2650 genes in adipose tissues (p = <0.05). Using concepts like correlation, semantic similarity, and theoretical graph parameters we narrowed down genes to a network of 23 genes strongly connected with NF-κB family with higher significance. Functional enrichment analysis revealed 21 of 23 target genes of NF-κB were found to have a critical role in the pathophysiology of obesity. Interestingly, GEM and PPP1R13L were predicted as novel genes which may act as potential target or biomarkers of obesity as they occur with other 21 target genes with known obesity relationship. Our study concludes that NF-κB and prioritized target genes regulate the inflammation in adipose tissues through several molecular signaling pathways like NF-κB, PI3K-Akt, glucocorticoid receptor regulatory network, angiogenesis and cytokine pathways. This integrated system biology approaches can be applied for elucidating functional protein interaction networks of NF-κB protein family in different complex diseases. Our integrative and network-based approach for finding therapeutic targets in genomic data could accelerate the identification of novel drug targets for obesity.

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.

Protein phenotype diagnosis of autosomal dominant calmodulin mutations causing irregular heart rhythms

The life-threatening group of irregular cardiac rhythmic disorders also known as Cardiac Arrhythmias (CA) are caused by mutations in highly conserved Calmodulin (CALM/CaM) genes. Herein, we present a multidimensional approach to diagnose changes in phenotypic, stability, and Ca²⁺ ion binding properties of CA-causing mutations. Mutation pathogenicity was determined by diverse computational machine learning approaches. We further modeled the mutations in 3D protein structure and analyzed residue level phenotype plasticity. We have also examined the influence of torsion angles, number of H-bonds, and free energy dynamics on the stability, near-native simulation dynamic potential of residue fluctuations in protein structures, Ca²⁺ ion binding potentials, of CaM mutants. Our study recomends to use M-CAP method for measuring the pathogenicity of CA causing CaM variants. Interestingly, most CA-causing variants we analyzed, exists in either third (V/H-96, S/I-98, V-103) or fourth (G/V-130, V/E/H-132, H-134, P-136, G-141, and L-142) EF-hands located in carboxyl domains of the CaM molecule. We observed that the minor structural fluctuations caused by these variants are likely tolerable owing to the highly flexible nature of calmodulin's globular domains. However, our molecular docking results supports that these variants disturb the affinity of CaM toward Ca²⁺ ions and corroborate previous findings from functional studies. Taken together, these computational findings can explain the molecular reasons for subtle changes in structure, flexibility, and stability aspects of mutant CaM molecule. Our comprehensive molecular scanning approach demonstrates the utility of computational methods in quick preliminary screening of CA- CaM mutations before undertaking time consuming and complicated functional laboratory assays.

The Transcriptomic Analysis of Circulating Immune Cells in a Celiac Family Unveils Further Insights Into Disease Pathogenesis

Celiac disease (CD), the most common chronic enteropathy worldwide, is triggered and sustained by a dysregulated immune response to dietary gluten in genetically susceptible individuals. Up to date either the role of environmental factors and the pathways leading to mucosal damage have been only partially unraveled. Therefore, we seized the unique opportunity to study a naturally-occurring experimental model of a family composed of both parents suffering from CD (one on a gluten-free diet) and two non-celiac daughters. The control group consisted in four unrelated cases, two celiac and two non-celiac subjects, all matching with family members for both disease status and genetic susceptibility. In this privileged setting, we sought to investigate gene expression in peripheral blood mononuclear cells (PBMCs), a population known to mirror the immune response state within the gut. To this purpose, PBMCs were obtained from the four biopsied-proven CD patients and the four non-celiac cases. Each group included two family members and two unrelated control subjects. After RNA purification and cDNA synthesis, each sample underwent a microarray screen on a whole-transcriptome scale, and the hybridization results were visualized by hierarchical clustering. Differentially expressed genes (DEG) were partitioned into clusters displaying comparable regulations among samples. These clusters were subjected to both functional and pathway analysis by using the Kyoto Encyclopedia of Genes and Genomes. Interestingly, on a global gene expression level, the family members clustered together, regardless of their disease status. A relevant fraction of DEG belonged to a limited number of pathways, and could be differentiated based on disease status: active CD vs. treated CD and CD vs. controls. These pathways were mainly involved in immune function regulation, cell-cell junctions, protein targeting and degradation, exosome trafficking, and signal transduction. Worth of noting, a small group of genes mapping on the male-specific region of the Y chromosome, and previously linked to cardiovascular risk, was found to be strongly upregulated in the active CD case belonging to the family, who suddenly died of a heart attack. Our results provide novel information on CD pathogenesis and may be useful in identifying new therapeutic tools and risk factors associated with this condition.

Computational Protein Phenotype Characterization of IL10RA Mutations Causative to Early Onset Inflammatory Bowel Disease (IBD)

The deleterious amino acid substitution mutations in IL-10 receptor alpha gene are most frequently reported in several autoimmune diseases including early onset-inflammatory bowel disease (IBD). Despite the important role of IL-10 RA in maintaining immune homeostasis, the specific structural and functional implications of these mutations on protein phenotype, stability, ligand binding and post translational characteristics is not well explored. Therefore, this study performed the multidimensional computational analysis of IL10RA missense variations causative to pediatric or early onset inflammatory bowel disease (<5 years of age). Our computational algorithmic screening identified the deleterious nature of p. W45G, p. Y57C, p. W69G, p.T84I, p.Y91C, p.R101W, p.R117C, and p.R117H, IBD causative IL10-RA mutations. The sensitivity and specificity analysis of different computational methods showed that CADD outperform SIFT, PolyPhen 2.0, FATHMM, LRT, MetaLR, MetaSVM, PROVEAN and Condel in predicting the pathogenicity of IL10RA mutations. Our three-dimensional protein modeling assays showed that the point mutations cause major drifts in the structural plasticity of IL10 RA molecule and negatively influence its stability. Findings from molecular docking analysis have shown that these point mutations decrease the binding affinity of IL10RA toward IL10 and may likely to disturb the IL10 signaling pathway. This study provides an easy frame work for phenotypic characterization of mutant IL10RA molecule in terms of structure, flexibility and stability aspects. Our approach may also add a new dimension to conventional functional biology assays in quickly studying IL10 RA mutations and also for designing and developing inhibitors for mutant IL10RA molecule.

Understanding Celiac Disease From Genetics to the Future Diagnostic Strategies

Celiac disease (CD) is an autoimmune disorder characterized by the permanent inflammation of the small bowel, triggered by the ingestion of gluten. It is associated with a number of symptoms, the most common being gastrointestinal. The prevalence of this illness worldwide is 1%. One of the main problems of CD is its difficulty to be diagnosed due to the various presentations of the disease. Besides, in many cases, CD is asymptomatic. Celiac disease is a multifactorial disease, HLA-DQ2 and HLA-DQ8 haplotypes are predisposition factors. Nowadays, molecular markers are being studied as diagnostic tools. In this review, we explore CD from its basic concept, manifestations, types, current and future methods of diagnosis, and associated disorders. Before addressing the therapeutic approaches, we also provide a brief overview of CD genetics and treatment.