Quality of life in children with Down syndrome and its association with parent and child demographic characteristics: Parent-reported measures.

BACKGROUND: This study aims to explore the association between the quality of life (QoL) in children with Down syndrome (DS) and its relationship with demographic characteristics of both parents and children. The investigation encompasses five domains: physical and psychological well-being, autonomy and parental relationship, social well-being, and peers, as well as school and the learning environment. METHOD: An online questionnaire, the KIDSCREEN-27, was used to measure the QoL of 112 families with DS in Saudi Arabia, referred to as "Parent-Reported Measures." Descriptive statistics were analyzed using the Statistical Package for Social Sciences. RESULTS: The study found that the QoL of children with DS showed high scores in the psychological well-being, autonomy, parental relations, school, and learning environment domains. However, the physical and social well-being and peer domains had lower scores, although still considered "good scores." Family income had a positively significant influence on all QoL domains. Specifically, higher family income was associated with better QoL outcomes, except for social well-being. Parental age was found to influence psychological well-being, while parental education and the relationship between the parent and child influenced social well-being. Lastly, the child's gender was found to have an impact on the school and learning environment domain. CONCLUSION: The study highlights the importance of understanding the impact of the demographic variability of children with DS and their parents on the QoL of their children. It emphasizes the need to address the needs of families with lower incomes and the importance of parental education and relationships with their children in improving social well-being. The findings could aid policymakers and healthcare providers in improving the QoL for families with children who have DS.

Identification of a de novo LRP1 mutation in a Saudi family with Tetralogy of Fallot.

BACKGROUND: Tetralogy of Fallot (TOF) is a rare, complex congenital heart defect caused by genetic and environmental interactions that results in abnormal heart development during the early stages of pregnancy. Genetic basis of TOF in Saudi populations is not yet studied. Therefore, the objective of this study is to screen for the molecular defects causing TOF in Saudi patients. METHODS: A family with non-syndromic TOF was recruited from the Western region of Saudi Arabia. Whole exome sequencing (WES) was performed on the proband and her parents. The identified candidate variant was verified by sanger sequencing. Also, different computational biology tools were used to figure out how candidate variants affect the structure and function of candidate protein involved in TOF. RESULTS: A novel heterozygous de novo mutation in LRP1 (p. G3311D) gene was identified in the index case. Also, this variant was absent in the in-house exome sequencing data of 80 healthy Saudi individuals. This variant was predicted to be likely pathogenic, as it negatively affects the biophysical chemical properties and stability of the protein. Furthermore, functional biology data from knock out mouse models confirms that molecular defects in LRP1 gene leads to cardiac defects and lethality. This variant was not previously reported in both Arab and global population genetic databases. CONCLUSION: The findings in this study postulate that the LRP1 variant has a role in TOF pathogenesis and facilitate accurate diagnosis as well as the understanding of underlying molecular mechanisms and pathophysiology of the disease.

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.

Identification of miRNA-mRNA-TFs regulatory network and crucial pathways involved in asthma through advanced systems biology approaches.

Asthma is a life-threatening and chronic inflammatory lung disease that is posing a true global health challenge. The genetic basis of the disease is fairly well examined. However, the molecular crosstalk between microRNAs (miRNAs), target genes, and transcription factors (TFs) networks and their contribution to disease pathogenesis and progression is not well explored. Therefore, this study was aimed at dissecting the molecular network between mRNAs, miRNAs, and TFs using robust computational biology approaches. The transcriptomic data of bronchial epithelial cells of severe asthma patients and healthy controls was studied by different systems biology approaches like differentially expressed gene detection, functional enrichment, miRNA-target gene pairing, and mRNA-miRNA-TF molecular networking. We detected the differential expression of 1703 (673 up-and 1030 down-regulated) genes and 71 (41 up-and 30 down-regulated) miRNAs in the bronchial epithelial cells of asthma patients. The DEGs were found to be enriched in key pathways like IL-17 signaling (KEGG: 04657), Th1 and Th2 cell differentiation (KEGG: 04658), and the Th17 cell differentiation (KEGG: 04659) (p-values = 0.001). The results from miRNAs-target gene pairs-transcription factors (TFs) have detected the key roles of 3 miRs (miR-181a-2-3p; miR-203a-3p; miR-335-5p), 6 TFs (TFAM, FOXO1, GFI1, IRF2, SOX9, and HLF) and 32 miRNA target genes in eliciting autoimmune reactions in bronchial epithelial cells of the respiratory tract. Through systemic implementation of comprehensive system biology tools, this study has identified key miRNAs, TFs, and miRNA target gene pairs as potential tissue-based asthma biomarkers.

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.

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.

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.

Novel MYO1D Missense Variant Identified Through Whole Exome Sequencing and Computational Biology Analysis Expands the Spectrum of Causal Genes of Laterality Defects.

Laterality defects (LDs) or asymmetrically positioned organs are a group of rare developmental disorders caused by environmental and/or genetic factors. However, the exact molecular pathophysiology of LD is not yet fully characterised. In this context, studying Arab population presents an ideal opportunity to discover the novel molecular basis of diseases owing to the high rate of consanguinity and genetic disorders. Therefore, in the present study, we studied the molecular basis of LD in Arab patients, using next-generation sequencing method. We discovered an extremely rare novel missense variant in MYO1D gene (Pro765Ser) presenting with visceral heterotaxy and left isomerism with polysplenia syndrome. The proband in this index family has inherited this homozygous variant from her heterozygous parents following the autosomal recessive pattern. This is the first report to show MYO1D genetic variant causing left-right axis defects in humans, besides previous known evidence from zebrafish, frog and Drosophila models. Moreover, our multilevel bioinformatics-based structural (protein variant structural modelling, divergence, and stability) analysis has suggested that Ser765 causes minor structural drifts and stability changes, potentially affecting the biophysical and functional properties of MYO1D protein like calmodulin binding and microfilament motor activities. Functional bioinformatics analysis has shown that MYO1D is ubiquitously expressed across several human tissues and is reported to induce severe phenotypes in knockout mouse models. In conclusion, our findings show the expanded genetic spectrum of LD, which could potentially pave way for the novel drug target identification and development of personalised medicine for high-risk families.

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 differential analysis of uterine leiomyomas and leiomyosarcomas through weighted gene network and pathway tracing approaches.

Uterine smooth muscular neoplastic growths like benign leiomyomas (UL) and metastatic leiomyosarcomas (ULMS) share similar clinical symptoms, radiological and histological appearances making their clinical distinction a difficult task. Therefore, the objective of this study is to identify key genes and pathways involved in transformation of UL to ULMS through molecular differential analysis. Global gene expression profiles of 25 ULMS, 25 UL, and 29 myometrium (Myo) tissues generated on Affymetrix U133A 2.0 human genome microarrays were analyzed by deploying robust statistical, molecular interaction network, and pathway enrichment methods. The comparison of expression signals across Myo vs UL, Myo vs ULMS, and UL vs ULMS groups identified 249, 1037, and 716 significantly expressed genes, respectively (p ≤ 0.05). The analysis of 249 DEGs from Myo vs UL confirms multistage dysregulation of various key pathways in extracellular matrix, collagen, cell contact inhibition, and cytokine receptors transform normal myometrial cells to benign leiomyomas (p value ≤ 0.01). The 716 DEGs between UL vs ULMS were found to affect cell cycle, cell division related Rho GTPases and PI3K signaling pathways triggering uncontrolled growth and metastasis of tumor cells (p value ≤ 0.01). Integration of gene networking data, with additional parameters like estimation of mutation burden of tumors and cancer driver gene identification, has led to the finding of 4 hubs (JUN, VCAN, TOP2A, and COL1A1) and 8 bottleneck genes (PIK3R1, MYH11, KDR, ESR1, WT1, CCND1, EZH2, and CDKN2A), which showed a clear distinction in their distribution pattern among leiomyomas and leiomyosarcomas. This study provides vital clues for molecular distinction of UL and ULMS which could further assist in identification of specific diagnostic markers and therapeutic targets.Abbreviations UL: Uterine Leiomyomas; ULMS: Uterine Leiomyosarcoma; Myo: Myometrium; DEGs: Differential Expressed Genes; RMA: Robust Multiarray Average; DC: Degree of Centrality; BC: Betweenness of Centrality; CGC: Cancer Gene Census; FDR: False Discovery Rate; TCGA: Cancer Genome Atlas; BP: Biological Process; CC: Cellular Components; MF: Molecular Function; PPI: Protein-Protein Interaction.

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.

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.

Molecular insights into the coding region mutations of low-density lipoprotein receptor adaptor protein 1 (LDLRAP1) linked to familial hypercholesterolemia.

BACKGROUND: Familial hypercholesterolemia (FH) is a lipid disorder caused by pathogenic mutations in LDLRAP1 gene. The present study has aimed to deepen our understanding about the pathogenicity predictions of FH causative genetic mutations, as well as their relationship to phenotype changes in LDLRAP1 protein, by utilizing multidirectional computational analysis. METHODS: FH linked LDLRAP1 mutations were mined from databases, and the prediction ability of several pathogenicity classifiers against these clinical variants, was assessed through different statistical measures. Furthermore, these mutations were 3D modelled in protein structures to assess their impact on protein phenotype changes. RESULTS: Our findings suggest that Polyphen-2, when compared with SIFT, M-CAP and CADD tools, can make better pathogenicity predictions for FH causative LDLRAP1 mutations. Through, 3D simulation and superimposition analysis of LDLRAP1 protein structures, it was found that missense mutations do not create any gross changes in the protein structure, although they could induce subtle structural changes at the level of amino acid residues. Near native molecular dynamic analysis revealed that missense mutations could induce variable degrees of fluctuation differences guiding the protein flexibility. Stability analysis showed that most missense mutations shifts the free energy equilibrium and hence they destabilize the protein. Molecular docking analysis demonstrates the molecular shifts in hydrogen and ionic bonds and Van der waals bonding properties, which further cause differences in the binding energy of LDLR-LDLRAP1 proteins. CONCLUSIONS: The diverse computational approaches used in the present study may provide a new dimension for exploring the structure-function relationship of the novel and deleterious LDLRAP1 mutations linked to FH.

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.

Dissecting the Role of NF-κb Protein Family and Its Regulators in Rheumatoid Arthritis Using Weighted Gene Co-Expression Network.

Rheumatoid arthritis (RA) is a chronic synovial autoinflammatory disease that destructs the cartilage and bone, leading to disability. The functional regulation of major immunity-related pathways like nuclear factor kappa B (NF-κB), which is involved in the chronic inflammatory reactions underlying the development of RA, remains to be explored. Therefore, this study has adopted statistical and knowledge-based systemic investigations (like gene correlation, semantic similarity, and topological parameters based on graph theory) to study the gene expression status of NF-κB protein family (NKPF) and its regulators in synovial tissues to trace the molecular pathways through which these regulators contribute to RA. A complex protein-protein interaction map (PPIM) of 2,742 genes and 37,032 interactions was constructed from differentially expressed genes (p ≤ 0.05). PPIM was further decomposed into a Regulator Allied Protein Interaction Network (RAPIN) based on the interaction between genes (5 NKPF, 31 seeds, 131 hubs, and 652 bottlenecks). Pathway network analysis has shown the RA-specific disturbances in the functional connectivity between seed genes (RIPK1, ATG7, TLR4, TNFRSF1A, KPNA1, CFLAR, SNW1, FOSB, PARVA, CX3CL1, and TRPC6) and NKPF members (RELA, RELB, NFKB2, and REL). Interestingly, these genes are known for their involvement in inflammation and immune system (signaling by interleukins, cytokine signaling in immune system, NOD-like receptor signaling, MAPK signaling, Toll-like receptor signaling, and TNF signaling) pathways connected to RA. This study, for the first time, reports that SNW1, along with other NK regulatory genes, plays an important role in RA pathogenesis and might act as potential biomarker for RA. Additionally, these genes might play important roles in RA pathogenesis, as well as facilitate the development of effective targeted therapies. Our integrative data analysis and network-based methods could accelerate the identification of novel drug targets for RA from high-throughput genomic data.

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.