A comparative mRNA- and miRNA transcriptomics reveals novel molecular signatures associated with metastatic prostate cancers.

Background: Prostate cancer (PC) is a fatally aggressive urogenital cancer killing millions of men, globally. Thus, this study aims to identify key miRNAs, target genes, and drug targets associated with prostate cancer metastasis. Methods: The miRNA and mRNA expression datasets of 148 prostate tissue biopsies (39 tumours and 109 normal tissues), were analysed by differential gene expression analysis, protein interactome mapping, biological pathway analysis, miRNA-mRNA networking, drug target analysis, and survival curve analysis. Results: The dysregulated expression of 53 miRNAs and their 250 target genes involved in Hedgehog, ErbB, and cAMP signalling pathways connected to cell growth, migration, and proliferation of prostate cancer cells was detected. The subsequent miRNA-mRNA network and expression status analysis have helped us in narrowing down their number to 3 hub miRNAs (hsa-miR-455-3p, hsa-miR-548c-3p, and hsa-miR-582-5p) and 9 hub genes (NFIB, DICER1, GSK3B, DCAF7, FGFR1OP, ABHD2, NACC2, NR3C1, and FGF2). Further investigations with different systems biology methods have prioritized NR3C1, ABHD2, and GSK3B as potential genes involved in prostate cancer metastasis owing to their high mutation load and expression status. Interestingly, down regulation of NR3C1 seems to improve the prostate cancer patient survival rate beyond 150 months. The NR3C1, ABHD2, and GSK3B genes are predicted to be targeted by hsa-miR-582-5p, besides some antibodies, PROTACs and inhibitory molecules. Conclusion: This study identified key miRNAs (miR-548c-3p and miR-582-5p) and target genes (NR3C1, ABHD2, and GSK3B) as potential biomarkers for metastatic prostate cancers from large-scale gene expression data using systems biology approaches.

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.

Bioinformatics insights into the genes and pathways on severe COVID-19 pathology in patients with comorbidities.

Background: Coronavirus disease (COVID-19) infection is known for its severe clinical pathogenesis among individuals with pre-existing comorbidities. However, the molecular basis of this observation remains elusive. Thus, this study aimed to map key genes and pathway alterations in patients with COVID-19 and comorbidities using robust systems biology approaches. Methods: The publicly available genome-wide transcriptomic datasets from 120 COVID-19 patients, 281 patients suffering from different comorbidities (like cardiovascular diseases, atherosclerosis, diabetes, and obesity), and 252 patients with different infectious diseases of the lung (respiratory syncytial virus, influenza, and MERS) were studied using a range of systems biology approaches like differential gene expression, gene ontology (GO), pathway enrichment, functional similarity, mouse phenotypic analysis and drug target identification. Results: By cross-mapping the differentially expressed genes (DEGs) across different datasets, we mapped 274 shared genes to severe symptoms of COVID-19 patients or with comorbidities alone. GO terms and functional pathway analysis highlighted genes in dysregulated pathways of immune response, interleukin signaling, FCGR activation, regulation of cytokines, chemokines secretion, and leukocyte migration. Using network topology parameters, phenotype associations, and functional similarity analysis with ACE2 and TMPRSS2-two key receptors for this virus-we identified 17 genes with high connectivity (CXCL10, IDO1, LEPR, MME, PTAFR, PTGS2, MAOB, PDE4B, PLA2G2A, COL5A1, ICAM1, SERPINE1, ABCB1, IL1R1, ITGAL, NCAM1 and PRKD1) potentially contributing to the clinical severity of COVID-19 infection in patients with comorbidities. These genes are predicted to be tractable and/or with many existing approved inhibitors, modulators, and enzymes as drugs. Conclusion: By systemic implementation of computational methods, this study identified potential candidate genes and pathways likely to confer disease severity in COVID-19 patients with pre-existing comorbidities. Our findings pave the way to develop targeted repurposed therapies in COVID-19 patients.

Genotype-protein phenotype characterization of NOD2 and IL23R missense variants associated with inflammatory bowel disease: A paradigm from molecular modelling, dynamics, and docking simulations.

Inflammatory bowel disease (IBD) is a gastrointestinal disease with an underlying contribution of genetic, microbial, environment, immunity factors. The coding region risk markers identified by IBD genome wide association studies have not been well characterized at protein phenotype level. Therefore, this study is conducted to characterize the role of NOD2 (Arg675Trp and Gly908Arg) and IL23R (Gly149Arg and Arg381Gln) missense variants on the structural and functional features of corresponding proteins. Thus, we used different variant pathogenicity assays, molecular modelling, secondary structure, stability, molecular dynamics, and molecular docking analysis methods. Our findings suggest that SIFT, Polyphen, GREP++, PhyloP, SiPhy and REVEL methods are very sensitive in determining pathogenicity of NOD2 and IL23R missense variants. We have also noticed that all the tested missense variants could potentially alter secondary (α-helices, ß-strands, and coils) and tertiary (residue level deviations) structural features. Moreover, our molecular dynamics (MD) simulation findings have simulated that NOD2 (Arg675Trp and Gly908Arg) and IL23R (Gly149Arg and Arg381Gln) variants creates rigid local structures comprising the protein flexibility and conformations. These predictions are corroborated by molecular docking results, where we noticed that NOD2 and IL23R missense variants induce molecular interaction deformities with RIPK2 and JAK2 ligand molecules, respectively. These functional alterations could potentially alter the signal transduction pathway cascade involved in inflammation and autoimmunity. Drug library searches and findings from docking studies have identified the inhibitory effects of Tacrolimus and Celecoxib drugs on NOD2 and IL23R variant forms, underlining their potential to contribute to personalized medicine for IBD. The present study supports the utilization of computational methods as primary filters (pre-in vitro and in vivo) in studying the disease potential mutations in the context of genptype-protein phenotype characteristics.

Identification of a Rare Exon 19 Skipping Mutation in ALMS1 Gene in Alström Syndrome Patients From Two Unrelated Saudi Families.

Background: Alström syndrome (AS) is a very rare childhood disorder characterized by cardiomyopathy, progressive hearing loss and blindness. Inherited genetic variants of ALMS1 gene are the known molecular cause of this disease. The objective of this study was to characterize the genetic basis and understand the genotype-phenotype relationship in Saudi AS patients. Methods: Clinical phenotyping and whole-exome sequencing (WES) analysis were performed on six AS patients belonging to two unrelated consanguineous Saudi families. Sanger sequencing was performed to determine the mode of inheritance of ALMS1 variant in first-degree family relatives and also to ensure its rare prevalence in 100 healthy population controls. Results: We identified that Alström patients from both the families were sharing a very rare ALMS1, 3'-splice site acceptor (c.11873-2 A>T) variant, which skips entire exon-19 and shortens the protein by 80 amino acids. This disease variant was inherited by AS patients in autosomal recessive mode and is not yet reported in any population-specific genetic databases. AS patients carrying this mutation showed heterogeneity in clinical presentations. Computational analysis of the mutant centroid structure of ALMS1 mRNA revealed that exon-19 skipping enlarges the hairpin loop and decreases the free energy, eventually affecting its folding pattern, stability, and function. Hence, we propose c.11873-2A as an AS causative potential founder mutation in Saudi Arabia because it is found in two families lacking a common lineage. Conclusions: We conclude that WES analysis potentially helps in clinical phenotyping, early diagnosis, and better clinical management of Alström patients showing variable clinical expressivity.

Molecular profiling of lamellar ichthyosis pathogenic missense mutations on the structural and stability aspects of TGM1 protein.

Lamellar ichthyosis (LI) is a rare inherited disease where affected infants present a extensive skin scaling characterized by hyperkeratosis. Inherited mutations in the Transglutaminase 1 (TGM1) protein is one of the known causative genetic factor for the LI. The main objective of this study is to explore the impact of LI causative missense mutations on the structural and stability aspects of TGM1 protein using structural modeling, molecular docking and molecular dynamics approaches. By testing all LI causative TMG1 mutations against multiple stability prediction methods, we found that L362R and L388P mutations positioned in the Transglut_core domain were most destabilizing to the stability of TGM1 protein. These 2 mutations were 3D protein modeled and further analyzed by molecular docking and dynamic simulation methods. Molecular docking of these TGM1 mutant structures with chitosan, a natural polyphenolic compound and known inducer for transglutaminase enzyme, has shown stable molecular interactions between the native TGM1-chitosan and TGM1(L388P)-chitosan complex, when compared to the TGM1(L362R)-chitosan complex. Interestingly, molecular dynamics analysis have also yielded similar findings, where L388P-chitosan complex is shown to develop B-sheets and attain better stability, whereas TGM1-L362R complex possessed coils over the simulation period, pointing its highly destabilizing behavior on the protein structure. This study concludes that missense mutations in Transglut_core domain of the TGM1 are deleterious to the stability and structural changes of TGM1 protein and also suggest that chitosan molecule could act as a natural activator against few pathogenic TGM1 mutations. Communicated by Ramaswamy H. Sarma.

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.