eDGAR: a database of Disease-Gene Associations with annotated Relationships among genes

Identification of ferroptosis-related genes and potential drugs in osteoarthritis

BACKGROUND

Osteoarthritis (OA) is a common chronic degenerative joint disease in orthopedics, and ferroptosis is a newly identified mode of cell death present in OA. Inhibition of inflammatory cytokine expression and modulation of chondrocyte ferroptosis related pathways may be novel strategies for the treatment of OA. The purpose of this work was to uncover prospective biomarkers and molecular processes of ferroptosis in OA, as well as to better understand the molecular mechanisms of ferroptosis in OA treated with resveratrol.

MATERIAL AND METHODS

We obtained OA gene expression profiles from the Gene Expression Omnibus (GEO) database. OA-expressed ferroptosis-related genes were identified using Genecards data, differential gene analysis, and weighted gene co-expression network analysis. Enrichment analysis was utilized to identify signaling pathways and molecular mechanisms linked with ferroptosis in OA, while immune infiltration analysis indicated immune cell infiltration in OA. The action targets of resveratrol were taken from the TCM database to determine the therapeutic targets of resveratrol for the treatment of OA. To validate the molecular process, molecular docking was performed using the therapeutic targets' enrichment analysis. Finally, in vitro investigations confirmed the molecular mechanism of ferroptosis in resveratrol-treated OA.

RESULTS

Bioinformatic analysis identified 462 OA ferroptosis gene sets, with GPX4, TFRC, SLC7A11, EGFR, and IL1B serving as significant hub genes. Enrichment analysis revealed that ferroptosis was also linked to animal mitophagy, the FoxO signaling pathway, the Toll-like receptor signaling pathway, the PI3K-Akt signaling pathway, inflammation, immune response activation, and cellular autophagy. The immune infiltration data revealed that T_cells_CD4_memory_resting, T_cells_CD4_memory_activated, NK_cells_activated, and Mast_cells_activated were considerably infiltrated in OA. Resveratrol ameliorated OA via modulating autophagy and ferroptosis via GPX4, TFRC, SLC7A11, EGFR, and IL1B, according to a mechanistic study.

CONCLUSION

We discovered the mechanism of GPX4, TFRC, SLC7A11, and EGFR, IL1B ferroptosis-related genes in OA, and preliminary evidence suggests that resveratrol improves OA by regulating ferroptosis and immunological processes, which may give a new route for OA treatment.

Loop Catalog: a comprehensive HiChIP database of human and mouse samples

HiChIP enables cost-effective and high-resolution profiling of chromatin loops. To leverage the increasing number of HiChIP datasets, we developed Loop Catalog (https://loopcatalog.lji.org), a web-based database featuring loop calls from 1000+ distinct human and mouse HiChIP samples from 152 studies plus 44 high-resolution Hi-C samples. We demonstrate its utility for interpreting GWAS and eQTL variants through SNP-to-gene linking, identifying enriched sequence motifs and motif pairs, and generating regulatory networks and 2D representations of chromatin structure. Our catalog spans over 4.19M unique loops, and with embedded analysis modules, constitutes an important resource for the field.

The landscape of structural variation in aye-ayes (Daubentonia madagascariensis)

Aye-ayes (Daubentonia madagascariensis) are one of the 25 most critically endangered primate species in the world. Endemic to Madagascar, their small and highly fragmented populations make them particularly vulnerable to both genetic disease and anthropogenic environmental changes. Over the past decade, conservation genomic efforts have largely focused on inferring and monitoring population structure based on single nucleotide variants to identify and protect critical areas of genetic diversity. However, the recent release of a highly contiguous genome assembly allows, for the first time, for the study of structural genomic variation (deletions, duplications, insertions, and inversions) which are likely to impact a substantial proportion of the species' genome. Based on whole-genome, short-read sequencing data from 14 individuals, >1,000 high-confidence autosomal structural variants were detected, affecting ~240 kb of the aye-aye genome. The majority of these variants (>85%) were deletions shorter than 200 bp, consistent with the notion that longer structural mutations are often associated with strongly deleterious fitness effects. For example, two deletions longer than 850 bp located within disease-linked genes were predicted to impose substantial fitness deficits owing to a resulting frameshift and gene fusion, respectively; whereas several other major effect variants outside of coding regions are likely to impact gene regulatory landscapes. Taken together, this first glimpse into the landscape of structural variation in aye-ayes will enable future opportunities to advance our understanding of the traits impacting the fitness of this endangered species, as well as allow for enhanced evolutionary comparisons across the full primate clade.

Landscape of genomic structural variations in Indian population-based cohorts: Deeper insights into their prevalence and clinical relevance

Structural variations (SV) are large (>50 base pairs) genomic rearrangements comprising deletions, duplications, insertions, inversions, and translocations. Studying SVs is important because they play active and critical roles in regulating gene expression, determining disease predispositions, and identifying population-specific differences among individuals of diverse ancestries. However, SV discoveries in the Indian population using whole-genome sequencing (WGS) have been limited. In this study, using short-read WGS having an average 42X depth of coverage, we identify and characterize 36,210 SVs from 529 individuals enrolled in population-based cohorts in India. These SVs include 24,574 deletions, 2,913 duplications, 8,710 insertions, and 13 inversions; 1.26% (456 out of 36,210) of the identified SVs can potentially impact the coding regions of genes. Furthermore, 56 of these SVs are highly intolerant to loss-of-function changes to the mapped genes, and five SVs impacting ADAMTS17, CCDC40, and RHCE are common in our study individuals. Seven rare SVs significantly impact dosage sensitivity of genes known to be associated with various clinical phenotypes. Most of the SVs in our study are rare and heterozygous. This fine-scale SV discovery in the underrepresented Indian population provides valuable insights that extend beyond Eurocentric human genetic studies.

Open Genes-a new comprehensive database of human genes associated with aging and longevity

The Open Genes database was created to enhance and simplify the search for potential aging therapy targets. We collected data on 2402 genes associated with aging and developed convenient tools for searching and comparing gene features. A comprehensive description of genes has been provided, including lifespan-extending interventions, age-related changes, longevity associations, gene evolution, associations with diseases and hallmarks of aging, and functions of gene products. For each experiment, we presented the necessary structured data for evaluating the experiment's quality and interpreting the study's findings. Our goal was to stay objective and precise while connecting a particular gene to human aging. We distinguished six types of studies and 12 criteria for adding genes to our database. Genes were classified according to the confidence level of the link between the gene and aging. All the data collected in a database are provided both by an API and a user interface. The database is publicly available on a website at https://open-genes.org/.

ioSearch: An approach for identifying interacting multiomics biomarkers using a novel algorithm with application on breast cancer data sets

Identification of biomarkers by integrating multiple omics together is important because complex diseases occur due to an intricate interplay of various genetic materials. Traditional single-omics association tests neither explore this crucial interomics dependence nor identify moderately weak signals due to the multiple-testing burden. Conversely, multiomics data integration imparts complementary information but suffers from an increased multiple-testing burden, data diversity inherent with different omics features, high-dimensionality, and so forth. Most of the available methods address subtype classification using dimension-reduction techniques to circumvent the sample size issue but interacting multiomics biomarker identification methods are unavailable. We propose a two-step model that first investigates phenotype-omics association using logistic regression. Then, selects disease-associated omics using sparse principal components which explores the interrelationship of multiple variables from two omics in a multivariate multiple regression framework. On the basis of this model, we developed a multiomics biomarker identification algorithm, interacting omics search (ioSearch), that jointly tests the effect of multiple omics with disease and between-omics associations by using pathway information that subsequently reduces the multiple-testing burden. Further, inference in terms of p values potentially makes it an easily interpretable biomarker identification tool. Extensive simulation demonstrates ioSearch as statistically powerful with a controlled Type-I error rate. Its application to publicly available breast cancer data sets identified relevant omics features in important pathways.

Unleashing the potential of CRISPR multiplexing: Harnessing Cas12 and Cas13 for precise gene modulation in eye diseases

Gene therapy is a flourishing field with the potential to revolutionize the treatment of genetic diseases. The emergence of CRISPR-Cas9 has significantly advanced targeted and efficient genome editing. Although CRISPR-Cas9 has demonstrated promising potential applications in various genetic disorders, it faces limitations in simultaneously targeting multiple genes. Novel CRISPR systems, such as Cas12 and Cas13, have been developed to overcome these challenges, enabling multiplexing and providing unique advantages. Cas13, in particular, targets mRNA instead of genomic DNA, permitting precise gene expression control and mitigating off-target effects. This review investigates the potential of Cas12 and Cas13 in ocular gene therapy applications, such as suppression of inflammation and cell death. In addition, the capabilities of Cas12 and Cas13 are explored in addressing potential targets related with disease mechanisms such as aberrant isoforms, mitochondrial genes, cis-regulatory sequences, modifier genes, and long non-coding RNAs. Anatomical accessibility and relative immune privilege of the eye provide an ideal organ system for evaluating these novel techniques' efficacy and safety. By targeting multiple genes concurrently, CRISPR-Cas12 and Cas13 systems hold promise for treating a range of ocular disorders, including glaucoma, retinal dystrophies, and age-related macular degeneration. Nonetheless, additional refinement is required to ascertain the safety and efficacy of these approaches in ocular disease treatments. Thus, the development of Cas12 and Cas13 systems marks a significant advancement in gene therapy, offering the potential to devise effective treatments for ocular disorders.

A large-scale targeted proteomics of serum and tissue shows the utility of classifying high grade and low grade meningioma tumors

BACKGROUND

Meningiomas are the most prevalent primary brain tumors. Due to their increasing burden on healthcare, meningiomas have become a pivot of translational research globally. Despite many studies in the field of discovery proteomics, the identification of grade-specific markers for meningioma is still a paradox and requires thorough investigation. The potential of the reported markers in different studies needs further verification in large and independent sample cohorts to identify the best set of markers with a better clinical perspective.

METHODS

A total of 53 fresh frozen tumor tissue and 51 serum samples were acquired from meningioma patients respectively along with healthy controls, to validate the prospect of reported differentially expressed proteins and claimed markers of Meningioma mined from numerous manuscripts and knowledgebases. A small subset of Glioma/Glioblastoma samples were also included to investigate inter-tumor segregation. Furthermore, a simple Machine Learning (ML) based analysis was performed to evaluate the classification accuracy of the list of proteins.

RESULTS

A list of 15 proteins from tissue and 12 proteins from serum were found to be the best segregator using a feature selection-based machine learning strategy with an accuracy of around 80% in predicting low grade (WHO grade I) and high grade (WHO grade II and WHO grade III) meningiomas. In addition, the discriminant analysis could also unveil the complexity of meningioma grading from a segregation pattern, which leads to the understanding of transition phases between the grades.

CONCLUSIONS

The identified list of validated markers could play an instrumental role in the classification of meningioma as well as provide novel clinical perspectives in regard to prognosis and therapeutic targets.

XGDAG: explainable gene-disease associations via graph neural networks

MOTIVATION

Disease gene prioritization consists in identifying genes that are likely to be involved in the mechanisms of a given disease, providing a ranking of such genes. Recently, the research community has used computational methods to uncover unknown gene-disease associations; these methods range from combinatorial to machine learning-based approaches. In particular, during the last years, approaches based on deep learning have provided superior results compared to more traditional ones. Yet, the problem with these is their inherent black-box structure, which prevents interpretability.

RESULTS

We propose a new methodology for disease gene discovery, which leverages graph-structured data using graph neural networks (GNNs) along with an explainability phase for determining the ranking of candidate genes and understanding the model's output. Our approach is based on a positive-unlabeled learning strategy, which outperforms existing gene discovery methods by exploiting GNNs in a non-black-box fashion. Our methodology is effective even in scenarios where a large number of associated genes need to be retrieved, in which gene prioritization methods often tend to lose their reliability.

AVAILABILITY AND IMPLEMENTATION

The source code of XGDAG is available on GitHub at: https://github.com/GiDeCarlo/XGDAG. The data underlying this article are available at: https://www.disgenet.org/, https://thebiogrid.org/, https://doi.org/10.1371/journal.pcbi.1004120.s003, and https://doi.org/10.1371/journal.pcbi.1004120.s004.

A novel gluten knowledge base of potential biomedical and health-related interactions extracted from the literature: using machine learning and graph analysis methodologies to reconstruct the bibliome

BACKGROUND

In return for their nutritional properties and broad availability, cereal crops have been associated with different alimentary disorders and symptoms, with the majority of the responsibility being attributed to gluten. Therefore, the research of gluten-related literature data continues to be produced at ever-growing rates, driven in part by the recent exploratory studies that link gluten to non-traditional diseases and the popularity of gluten-free diets, making it increasingly difficult to access and analyse practical and structured information. In this sense, the accelerated discovery of novel advances in diagnosis and treatment, as well as exploratory studies, produce a favourable scenario for disinformation and misinformation.

OBJECTIVES

Aligned with, the European Union strategy "Delivering on EU Food Safety and Nutrition in 2050" which emphasizes the inextricable links between imbalanced diets, the increased exposure to unreliable sources of information and misleading information, and the increased dependency on reliable sources of information; this paper presents GlutKNOIS, a public and interactive literature-based database that reconstructs and represents the experimental biomedical knowledge extracted from the gluten-related literature. The developed platform includes different external database knowledge, bibliometrics statistics and social media discussion to propose a novel and enhanced way to search, visualise and analyse potential biomedical and health-related interactions in relation to the gluten domain.

METHODS

For this purpose, the presented study applies a semi-supervised curation workflow that combines natural language processing techniques, machine learning algorithms, ontology-based normalization and integration approaches, named entity recognition methods, and graph knowledge reconstruction methodologies to process, classify, represent and analyse the experimental findings contained in the literature, which is also complemented by data from the social discussion.

RESULTS

and Conclusions: In this sense, 5,814 documents were manually annotated and 7,424 were fully automatically processed to reconstruct the first online gluten-related knowledge database of evidenced health-related interactions that produce health or metabolic changes based on the literature. In addition, the automatic processing of the literature combined with the knowledge representation methodologies proposed has the potential to assist in the revision and analysis of years of gluten research. The reconstructed knowledge base is public and accessible at https://sing-group.org/glutknois/.

Topology Driven Analysis of Protein - Protein Interactome for Prioritizing Key Comorbid Genes via Sub Graph Based Average Path Length Centrality

In gene-based therapies, local perturbations associated with one disease can lead to comorbidity as it influences the pathways involved with the other diseases. The key genes orchestrating the common biological mechanisms are need to be prioritized for addressing the challenges introduced by the cross talks between disease modules. Here, a local centrality measure named Sub graph based Average Path length Double Specific Betweenness centrality (SAPDSB) for prioritizing the comorbid genes via Protein-Protein Interaction Network (PPIN) analysis is presented. This approach can be used to identify putative biomarkers which can be repurposed for the management of comorbidity. Proposed network based topological measure is designed specifically to prioritize the comorbid genes that are most likely to be present in the overlap of disease modules. In order to attain this, the estimated average path length of the seed network which holds Protein-Protein Interactions (PPIs) of the disease genes is exploited. Prioritized comorbid genes are further pruned using centrality-based cut-off values and specificity scores. The biological significance of the resultant genes is corroborated with connectivity analysis using leave-one-out method, pathway enrichment analysis and a comparative analysis using single disease-based gene prioritization tools. For performance analysis, proposed approach is tested using case studies involving common diseases and rare neurodegenerative diseases. For case study1, diseases such as Diabetes, Carcinoma and Alzheimer's are considered in a pairwise manner while for case study2, Amyotrophic Lateral Sclerosis (ALS) and Spinal Muscular Atrophy (SMA) are considered. As outcome, prioritized candidate genes and biological pathways associated with respective disease pairs have been found. The associations from top 10 candidate genes in different disease pair combinations of Diabetes-Carcinoma-Alzheimer's revealed common genes like CREBBP, TP53, HSP90AA1 and the common pathway namely p53 pathway feedback loops 2. Out of the pathways retrieved from the top 10 genes associated with ALS-SMA disease pair, 60% of unique pathways are found to be leading to both diseases and its comorbidities. Comparative analysis of the proposed method with recent similar approach also reported a clear degree of benefits in performance.

miRNA-mRNA-protein dysregulated network in COPD in women

Rationale: Chronic obstructive pulmonary disease (COPD) is a complex disease caused by a multitude of underlying mechanisms, and molecular mechanistic modeling of COPD, especially at a multi-molecular level, is needed to facilitate the development of molecular diagnostic and prognostic tools and efficacious treatments. Objectives: To investigate the miRNA-mRNA-protein dysregulated network to facilitate prediction of biomarkers and disease subnetwork in COPD in women. Measurements and Results: Three omics data blocks (mRNA, miRNA, and protein) collected from BAL cells from female current-smoker COPD patients, smokers with normal lung function, and healthy never-smokers were integrated with miRNA-mRNA-protein regulatory networks to construct a COPD-specific dysregulated network. Furthermore, downstream network topology, literature annotation, and functional enrichment analysis identified both known and novel disease-related biomarkers and pathways. Both abnormal regulations in miRNA-induced mRNA transcription and protein translation repression play roles in COPD. Finally, the let-7-AIFM1-FKBP1A pathway is highlighted in COPD pathology. Conclusion: For the first time, a comprehensive miRNA-mRNA-protein dysregulated network of primary immune cells from the lung related to COPD in females was constructed to elucidate specific biomarkers and disease pathways. The multi-omics network provides a new molecular insight from a multi-molecular aspect and highlights dysregulated interactions. The highlighted let-7-AIFM1-FKBP1A pathway also indicates new hypotheses of COPD pathology.

Validating methods for testing natural molecules on molecular pathways of interest in silico and in vitro

Differentially expressed genes can serve as drug targets and are used to predict drug response and disease progression. In silico drug analysis based on the expression of these genetic biomarkers allows the detection of putative therapeutic agents, which could be used to reverse a pathological gene expression signature. Indeed, a set of bioinformatics tools can increase the accuracy of drug discovery, helping in biomarker identification. Once a drug target is identified, in vitro cell line models of disease are used to evaluate and validate the therapeutic potential of putative drugs and novel natural molecules. This study describes the development of efficacious PCR primers that can be used to identify gene expression of specific genetic pathways, which can lead to the identification of natural molecules as therapeutic agents in specific molecular pathways. For this study, genes involved in health conditions and processes were considered. In particular, the expression of genes involved in obesity, xenobiotics metabolism, endocannabinoid pathway, leukotriene B4 metabolism and signaling, inflammation, endocytosis, hypoxia, lifespan, and neurotrophins were evaluated. Exploiting the expression of specific genes in different cell lines can be useful in in vitro to evaluate the therapeutic effects of small natural molecules.

Artificial Intelligence, Healthcare, Clinical Genomics, and Pharmacogenomics Approaches in Precision Medicine

Precision medicine has greatly aided in improving health outcomes using earlier diagnosis and better prognosis for chronic diseases. It makes use of clinical data associated with the patient as well as their multi-omics/genomic data to reach a conclusion regarding how a physician should proceed with a specific treatment. Compared to the symptom-driven approach in medicine, precision medicine considers the critical fact that all patients do not react to the same treatment or medication in the same way. When considering the intersection of traditionally distinct arenas of medicine, that is, artificial intelligence, healthcare, clinical genomics, and pharmacogenomics—what ties them together is their impact on the development of precision medicine as a field and how they each contribute to patient-specific, rather than symptom-specific patient outcomes. This study discusses the impact and integration of these different fields in the scope of precision medicine and how they can be used in preventing and predicting acute or chronic diseases. Additionally, this study also discusses the advantages as well as the current challenges associated with artificial intelligence, healthcare, clinical genomics, and pharmacogenomics.

Integrated Analysis of Tissue-Specific Gene Expression in Diabetes by Tensor Decomposition Can Identify Possible Associated Diseases

In the field of gene expression analysis, methods of integrating multiple gene expression profiles are still being developed and the existing methods have scope for improvement. The previously proposed tensor decomposition-based unsupervised feature extraction method was improved by introducing standard deviation optimization. The improved method was applied to perform an integrated analysis of three tissue-specific gene expression profiles (namely, adipose, muscle, and liver) for diabetes mellitus, and the results showed that it can detect diseases that are associated with diabetes (e.g., neurodegenerative diseases) but that cannot be predicted by individual tissue expression analyses using state-of-the-art methods. Although the selected genes differed from those identified by the individual tissue analyses, the selected genes are known to be expressed in all three tissues. Thus, compared with individual tissue analyses, an integrated analysis can provide more in-depth data and identify additional factors, namely, the association with other diseases.

Distinct Roles of NANOS1 and NANOS3 in the Cell Cycle and NANOS3-PUM1-FOXM1 Axis to Control G2/M Phase in a Human Primordial Germ Cell Model

Nanos RNA-binding proteins are critical factors of germline development throughout the animal kingdom and their dysfunction causes infertility. During evolution, mammalian Nanos paralogues adopted divergent roles in germ cell biology. However, the molecular basis behind this divergence, such as their target mRNAs, remains poorly understood. Our RNA-sequencing analysis in a human primordial germ cell model-TCam-2 cell line revealed distinct pools of genes involved in the cell cycle process downregulated upon NANOS1 and NANOS3 overexpression. We show that NANOS1 and NANOS3 proteins influence different stages of the cell cycle. Namely, NANOS1 is involved in the G1/S and NANOS3 in the G2/M phase transition. Many of their cell cycle targets are known infertility and cancer-germ cell genes. Moreover, NANOS3 in complex with RNA-binding protein PUM1 causes 3′UTR-mediated repression of FOXM1 mRNA encoding a transcription factor crucial for G2/M phase transition. Interestingly, while NANOS3 and PUM1 act as post-transcriptional repressors of FOXM1, FOXM1 potentially acts as a transcriptional activator of NANOS3, PUM1, and itself. Finally, by utilizing publicly available RNA-sequencing datasets, we show that the balance between FOXM1-NANOS3 and FOXM1-PUM1 expression levels is disrupted in testis cancer, suggesting a potential role in this disease.

Exploring the molecular mechanism of action of Polygonum capitatum Buch-Ham. ex D. Don for the treatment of bacterial prostatitis based on network pharmacology and experimental verification

ETHNOPHARMACOLOGY RELEVANCE

Polygonum capitatum Buch-Ham. ex D. Don (CNPC2009), a traditional Miao-national herbal medicine, has been widely used with considerable therapeutic efficacy in the treatment of various urologic disorders including prostatitis. However, the molecular mechanism of action (MOA) remains unclear.

AIM OF THE STUDY

In this study, UPLC-Q-Exactive-MS and Network pharmacological methods were used to explore the underlying molecular MOA of Polygonum capitatum Buch-Ham. Ex D. Don (P.capitatum) for the treatment bacterial prostatitis (BP).

MATERIALS AND METHODS

The UPLC-Q-Exactive-MS technique was used to identify the chemical components of P. capitatum. Databases such as SwissTargetPrediction, Gene Cards, and OMIM were used to predict the targets of P. capitatum for the treatment of BP. The Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) was used to analyze the protein-protein interaction (PPI) and construct a PPI network, and the Metascape was used for Gene Ontology (GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. In addition, experimental treatment of Escherichia coli (E.coli)-induced BP was verified.

RESULTS

A total of 31 molecular components were identified by UPLC-Q-Exactive-MS. Network pharmacology revealed that P. capitatum may act on the AKT1, PI3K, MTO, EGFR and other targets through active components such as Gallic acid, Quercetin, Luteolin, Protocatechuic Acid, Kaempferol and thereby regulate PI3K-AKT, ErbB, AMPK, HIF-1, and other signaling pathways to intervene in the pathological mechanism of BP. Verification through experimental results showed that compared with the model group, treatment with P. capitatum could significantly inhibit bacterial growth in prostate tissues, lowered the prostate index, down-regulated the levels of inflammatory mediators(IL-1β, IL-6, and TNF-α) in prostate tissues, and down-regulate the protein expression and mRNA expression levels of AKT and PI3K.

CONCLUSION

This study preliminarily revealed the MOA of P. capitatum for treating BP with multiple components, multiple targets, and multiple pathways, especially affecting the PI3K-AKT signaling pathways.

Neuronal and cardiac toxicity of pharmacological compounds identified through transcriptomic analysis of human pluripotent stem cell-derived embryoid bodies

Concurrent with the '3R' principle, the embryonic stem cell test (EST) using mouse embryonic stem cells, developed in 2000, remains the solely accepted in vitro method for embryotoxicity testing. However, the scope and implementation of EST for embryotoxicity screening, compliant with regulatory requirements, are limited. This is due to its technical complexity, long testing period, labor-intensive methodology, and limited endpoint data, leading to misclassification of embryotoxic potential. In this study, we used human induced pluripotent stem cell (hiPSC)-derived embryoid bodies (EB) as an in vitro model to investigate the embryotoxic effects of a carefully selected set of pharmacological compounds. Morphology, viability, and differentiation potential were investigated after exposing EBs to folic acid, all-trans-retinoic acid, dexamethasone, and valproic acid for 15 days. The results showed that the compounds differentially repressed cell growth, compromised morphology, and triggered apoptosis in the EBs. Further, transcriptomics was employed to compare subtle temporal changes between treated and untreated cultures. Gene ontology and pathway analysis revealed that dysregulation of a large number of genes strongly correlated with impaired neuroectoderm and cardiac mesoderm formation. This aberrant gene expression pattern was associated with several disorders of the brain like mental retardation, multiple sclerosis, stroke and of the heart like dilated cardiomyopathy, ventricular tachycardia, and ventricular arrhythmia. Lastly, these in vitro findings were validated using in ovo chick embryo model. Taken together, pharmacological compound or drug-induced defective EB development from hiPSCs could potentially be used as a suitable in vitro platform for embryotoxicity screening.

MOSES: A New Approach to Integrate Interactome Topology and Functional Features for Disease Gene Prediction

Disease gene prediction is to date one of the main computational challenges of precision medicine. It is still uncertain if disease genes have unique functional properties that distinguish them from other non-disease genes or, from a network perspective, if they are located randomly in the interactome or show specific patterns in the network topology. In this study, we propose a new method for disease gene prediction based on the use of biological knowledge-bases (gene-disease associations, genes functional annotations, etc.) and interactome network topology. The proposed algorithm called MOSES is based on the definition of two somewhat opposing sets of genes both disease-specific from different perspectives: warm seeds (i.e., disease genes obtained from databases) and cold seeds (genes far from the disease genes on the interactome and not involved in their biological functions). The application of MOSES to a set of 40 diseases showed that the suggested putative disease genes are significantly enriched in their reference disease. Reassuringly, known and predicted disease genes together, tend to form a connected network module on the human interactome, mitigating the scattered distribution of disease genes which is probably due to both the paucity of disease-gene associations and the incompleteness of the interactome.

A computational approach for identification of core modules from a co-expression network and GWAS data

This protocol describes the application of the "omnigenic" model of the genetic architecture of complex traits to identify novel "core" genes influencing a disease-associated phenotype. Core genes are hypothesized to directly regulate disease and may serve as therapeutic targets. This protocol leverages GWAS data, a co-expression network, and publicly available data, including the GTEx database and the International Mouse Phenotyping Consortium Database, to identify modules enriched for genes with "core-like" characteristics. For complete details on the use and execution of this protocol, please refer to Sabik et al. (2020).

An overview of germline variations in genes of primary immunodeficiences through integrative analysis of ClinVar, HGMD® and dbSNP databases

Primary immunodeficiencies (PID) are a diverse group of genetic disorders caused by inadequate development and function of immune system. Identifying genetic etiology is important for genetic counselling and treatment decisions. Clinical relevance of genetic variants is a complex problem depending on gene-specific and variant specific genotype-phenotype interactions. To address this challenge, we aimed to characterize the pathogenic landscape of PID genes by combining the analysis of germline variations reported in ClinVar and HGMD^® and identification of damaging variations available in dbSNP. We generated a joint ClinVar/HGMD database, which included 111,940 variants, among them 32,452 were classified as pathogenic/likely pathogenic. From a total of 5,415,794 bi- or multiallelic variants in PID genes recorded in dbSNP, we retrieved 38,291 high impact (HI) biallelic variants with presumably disruptive impact in the protein, of them 25,500 variants were not present in ClinVar/HGMD. Using a functional prediction algorithm, we additionally identified 28,507 deleterious and 56,016 neutral missense variants among dbSNP variants and created a collection of damaging and neutral variations in PID genes, not currently present in ClinVar/HGMD, with their allele frequencies and mappings to protein domains. The distribution of pathogenic variants from ClinVar/HGMD, HI variants and deleterious missense variants from dbSNP was analyzed in the context of hereditary pattern and gene specific metrics, such as pLI and haploinsufficiency. Our report summarized data on complex gene-specific variability in PID genes and might be useful for the identification of the most promising variants and gene regions for further study.

Advancing clinical genomics and precision medicine with GVViZ: FAIR bioinformatics platform for variable gene-disease annotation, visualization, and expression analysis

Background

Genetic disposition is considered critical for identifying subjects at high risk for disease development. Investigating disease-causing and high and low expressed genes can support finding the root causes of uncertainties in patient care. However, independent and timely high-throughput next-generation sequencing data analysis is still a challenge for non-computational biologists and geneticists.

Results

In this manuscript, we present a findable, accessible, interactive, and reusable (FAIR) bioinformatics platform, i.e., GVViZ (visualizing genes with disease-causing variants). GVViZ is a user-friendly, cross-platform, and database application for RNA-seq-driven variable and complex gene-disease data annotation and expression analysis with a dynamic heat map visualization. GVViZ has the potential to find patterns across millions of features and extract actionable information, which can support the early detection of complex disorders and the development of new therapies for personalized patient care. The execution of GVViZ is based on a set of simple instructions that users without a computational background can follow to design and perform customized data analysis. It can assimilate patients’ transcriptomics data with the public, proprietary, and our in-house developed gene-disease databases to query, easily explore, and access information on gene annotation and classified disease phenotypes with greater visibility and customization. To test its performance and understand the clinical and scientific impact of GVViZ, we present GVViZ analysis for different chronic diseases and conditions, including Alzheimer’s disease, arthritis, asthma, diabetes mellitus, heart failure, hypertension, obesity, osteoporosis, and multiple cancer disorders. The results are visualized using GVViZ and can be exported as image (PNF/TIFF) and text (CSV) files that include gene names, Ensembl (ENSG) IDs, quantified abundances, expressed transcript lengths, and annotated oncology and non-oncology diseases.

Conclusions

We emphasize that automated and interactive visualization should be an indispensable component of modern RNA-seq analysis, which is currently not the case. However, experts in clinics and researchers in life sciences can use GVViZ to visualize and interpret the transcriptomics data, making it a powerful tool to study the dynamics of gene expression and regulation. Furthermore, with successful deployment in clinical settings, GVViZ has the potential to enable high-throughput correlations between patient diagnoses based on clinical and transcriptomics data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-021-00336-1.

Mapping OMIM Disease-Related Variations on Protein Domains Reveals an Association Among Variation Type, Pfam Models, and Disease Classes

Human genome resequencing projects provide an unprecedented amount of data about single-nucleotide variations occurring in protein-coding regions and often leading to observable changes in the covalent structure of gene products. For many of these variations, links to Online Mendelian Inheritance in Man (OMIM) genetic diseases are available and are reported in many databases that are collecting human variation data such as Humsavar. However, the current knowledge on the molecular mechanisms that are leading to diseases is, in many cases, still limited. For understanding the complex mechanisms behind disease insurgence, the identification of putative models, when considering the protein structure and chemico-physical features of the variations, can be useful in many contexts, including early diagnosis and prognosis. In this study, we investigate the occurrence and distribution of human disease–related variations in the context of Pfam domains. The aim of this study is the identification and characterization of Pfam domains that are statistically more likely to be associated with disease-related variations. The study takes into consideration 2,513 human protein sequences with 22,763 disease-related variations. We describe patterns of disease-related variation types in biunivocal relation with Pfam domains, which are likely to be possible markers for linking Pfam domains to OMIM diseases. Furthermore, we take advantage of the specific association between disease-related variation types and Pfam domains for clustering diseases according to the Human Disease Ontology, and we establish a relation among variation types, Pfam domains, and disease classes. We find that Pfam models are specific markers of patterns of variation types and that they can serve to bridge genes, diseases, and disease classes. Data are available as Supplementary Material for 1,670 Pfam models, including 22,763 disease-related variations associated to 3,257 OMIM diseases.

Network bioinformatics analysis provides insight into drug repurposing for COVID-19

The COVID-19 disease caused by the SARS-CoV-2 virus is a health crisis worldwide. While developing novel drugs and vaccines is long, repurposing existing drugs against COVID-19 can yield treatments with known preclinical, pharmacokinetic, pharmacodynamic, and toxicity profiles, which can rapidly enter clinical trials. In this study, we present a novel network-based drug repurposing platform to identify candidates for the treatment of COVID-19. At the time of the initial outbreak, knowledge about SARS-CoV-2 was lacking, but based on its similarity with other viruses, we sought to identify repurposing candidates to be tested rapidly at the clinical or preclinical levels. We first analyzed the genome sequence of SARS-CoV-2 and confirmed SARS as the closest virus by genome similarity, followed by MERS and other human coronaviruses. Using text mining and database searches, we obtained 34 COVID-19-related genes to seed the construction of a molecular network where our module detection and drug prioritization algorithms identified 24 disease-related human pathways, five modules, and 78 drugs to repurpose. Based on clinical knowledge, we re-prioritized 30 potentially repurposable drugs against COVID-19 (including pseudoephedrine, andrographolide, chloroquine, abacavir, and thalidomide). Our work shows how in silico repurposing analyses can yield testable candidates to accelerate the response to novel disease outbreaks.

Approaches to Identify and Characterise the Post-Transcriptional Roles of lncRNAs in Cancer

It is becoming increasingly evident that the non-coding genome and transcriptome exert great influence over their coding counterparts through complex molecular interactions. Among non-coding RNAs (ncRNA), long non-coding RNAs (lncRNAs) in particular present increased potential to participate in dysregulation of post-transcriptional processes through both RNA and protein interactions. Since such processes can play key roles in contributing to cancer progression, it is desirable to continue expanding the search for lncRNAs impacting cancer through post-transcriptional mechanisms. The sheer diversity of mechanisms requires diverse resources and methods that have been developed and refined over the past decade. We provide an overview of computational resources as well as proven low-to-high throughput techniques to enable identification and characterisation of lncRNAs in their complex interactive contexts. As more cancer research strategies evolve to explore the non-coding genome and transcriptome, we anticipate this will provide a valuable primer and perspective of how these technologies have matured and will continue to evolve to assist researchers in elucidating post-transcriptional roles of lncRNAs in cancer.

An emerging role of chromatin-interacting RNA-binding proteins in transcription regulation

Transcription factors (TFs) are well-established key factors orchestrating gene transcription, and RNA-binding proteins (RBPs) are mainly thought to participate in post-transcriptional control of gene. In fact, these two steps are functionally coupled, offering a possibility for reciprocal communications between transcription and regulatory RNAs and RBPs. Recently, a series of exploratory studies, utilizing functional genomic strategies, have revealed that RBPs are prevalently involved in transcription control genome-wide through their interactions with chromatin. Here, we present a refined census of RBPs to grope for such an emerging role and discuss the global view of RBP-chromatin interactions and their functional diversities in transcription regulation.

High-depth African genomes inform human migration and health

The African continent is regarded as the cradle of modern humans and African genomes contain more genetic variation than those from any other continent, yet only a fraction of the genetic diversity among African individuals has been surveyed¹. Here we performed whole-genome sequencing analyses of 426 individuals-comprising 50 ethnolinguistic groups, including previously unsampled populations-to explore the breadth of genomic diversity across Africa. We uncovered more than 3 million previously undescribed variants, most of which were found among individuals from newly sampled ethnolinguistic groups, as well as 62 previously unreported loci that are under strong selection, which were predominantly found in genes that are involved in viral immunity, DNA repair and metabolism. We observed complex patterns of ancestral admixture and putative-damaging and novel variation, both within and between populations, alongside evidence that Zambia was a likely intermediate site along the routes of expansion of Bantu-speaking populations. Pathogenic variants in genes that are currently characterized as medically relevant were uncommon-but in other genes, variants denoted as 'likely pathogenic' in the ClinVar database were commonly observed. Collectively, these findings refine our current understanding of continental migration, identify gene flow and the response to human disease as strong drivers of genome-level population variation, and underscore the scientific imperative for a broader characterization of the genomic diversity of African individuals to understand human ancestry and improve health.

Highlighting Human Enzymes Active in Different Metabolic Pathways and Diseases: The Case Study of EC 1.2.3.1 and EC 2.3.1.9

Enzymes are key proteins performing the basic functional activities in cells. In humans, enzymes can be also responsible for diseases, and the molecular mechanisms underlying the genotype to phenotype relationship are under investigation for diagnosis and medical care. Here, we focus on highlighting enzymes that are active in different metabolic pathways and become relevant hubs in protein interaction networks. We perform a statistics to derive our present knowledge on human metabolic pathways (the Kyoto Encyclopaedia of Genes and Genomes (KEGG)), and we found that activity aldehyde dehydrogenase (NAD(+)), described by Enzyme Commission number EC 1.2.1.3, and activity acetyl-CoA C-acetyltransferase (EC 2.3.1.9) are the ones most frequently involved. By associating functional activities (EC numbers) to enzyme proteins, we found the proteins most frequently involved in metabolic pathways. With our analysis, we found that these proteins are endowed with the highest numbers of interaction partners when compared to all the enzymes in the pathways and with the highest numbers of predicted interaction sites. As specific enzyme protein test cases, we focus on Alpha-Aminoadipic Semialdehyde Dehydrogenase (ALDH7A1, EC 2.3.1.9) and Acetyl-CoA acetyltransferase, cytosolic and mitochondrial (gene products of ACAT2 and ACAT1, respectively; EC 2.3.1.9). With computational approaches we show that it is possible, by starting from the enzyme structure, to highlight clues of their multiple roles in different pathways and of putative mechanisms promoting the association of genes to disease.

Precise breakpoint detection in a patient with 9p- syndrome

We present a case of 9p- syndrome with a complex chromosomal event originally characterized by the classical karyotype approach as 46,XX,der(9)t(9;13)(p23;q13). We used advanced technologies (Bionano Genomics genome imaging and 10× Genomics sequencing) to characterize the location of the translocation and accompanying deletion on Chromosome 9 and duplication on Chromosome 13 with single-nucleotide breakpoint resolution. The translocation breakpoint was at Chr 9:190938 and Chr 13:50850492, the deletion at Chr 9:1-190938, and the duplication at Chr 13:50850492-114364328. We identified genes in the deletion and duplication regions that are known to be associated with this patient's phenotype (e.g., ZIC2 in dysmorphic facial features, FOXD4 in developmental delay, RNASEH2B in developmental delay, and PCDH9 in autism). Our results indicate that clinical genomic assessment of individuals with complex karyotypes can be refined to a single-base-pair resolution when utilizing Bionano and 10× Genomics sequencing. With the 10× Genomics data, we were also able to characterize other variation (e.g., loss of function) throughout the remainder of the patient's genome. Overall, the Bionano and 10× technologies complemented each other and provided important insight into our patient with 9p- syndrome. Altogether, these results indicate that newer technologies can identify precise genomic variants associated with unique patient phenotypes that permit discovery of novel genotype-phenotype correlations and therapeutic strategies.

Biomedical data and computational models for drug repositioning: a comprehensive review

Drug repositioning can drastically decrease the cost and duration taken by traditional drug research and development while avoiding the occurrence of unforeseen adverse events. With the rapid advancement of high-throughput technologies and the explosion of various biological data and medical data, computational drug repositioning methods have been appealing and powerful techniques to systematically identify potential drug-target interactions and drug-disease interactions. In this review, we first summarize the available biomedical data and public databases related to drugs, diseases and targets. Then, we discuss existing drug repositioning approaches and group them based on their underlying computational models consisting of classical machine learning, network propagation, matrix factorization and completion, and deep learning based models. We also comprehensively analyze common standard data sets and evaluation metrics used in drug repositioning, and give a brief comparison of various prediction methods on the gold standard data sets. Finally, we conclude our review with a brief discussion on challenges in computational drug repositioning, which includes the problem of reducing the noise and incompleteness of biomedical data, the ensemble of various computation drug repositioning methods, the importance of designing reliable negative samples selection methods, new techniques dealing with the data sparseness problem, the construction of large-scale and comprehensive benchmark data sets and the analysis and explanation of the underlying mechanisms of predicted interactions.

PanglaoDB: a web server for exploration of mouse and human single-cell RNA sequencing data

Single-cell RNA sequencing is an increasingly used method to measure gene expression at the single cell level and build cell-type atlases of tissues. Hundreds of single-cell sequencing datasets have already been published. However, studies are frequently deposited as raw data, a format difficult to access for biological researchers due to the need for data processing using complex computational pipelines. We have implemented an online database, PanglaoDB, accessible through a user-friendly interface that can be used to explore published mouse and human single cell RNA sequencing studies. PanglaoDB contains pre-processed and pre-computed analyses from more than 1054 single-cell experiments covering most major single cell platforms and protocols, based on more than 4 million cells from a wide range of tissues and organs. The online interface allows users to query and explore cell types, genetic pathways and regulatory networks. In addition, we have established a community-curated cell-type marker compendium, containing more than 6000 gene-cell-type associations, as a resource for automatic annotation of cell types.

NPInter v4.0: an integrated database of ncRNA interactions

Noncoding RNAs (ncRNAs) play crucial regulatory roles in a variety of biological circuits. To document regulatory interactions between ncRNAs and biomolecules, we previously created the NPInter database (http://bigdata.ibp.ac.cn/npinter). Since the last version of NPInter was issued, a rapidly growing number of studies have reported novel interactions and accumulated numerous high-throughput interactome data. We have therefore updated NPInter to its fourth edition in which are integrated 600 000 new experimentally identified ncRNA interactions. ncRNA-DNA interactions derived from ChIRP-seq data and circular RNA interactions have been included in the database. Additionally, disease associations were annotated to the interacting molecules. The database website has also been redesigned with a more user-friendly interface and several additional functional modules. Overall, NPInter v4.0 now provides more comprehensive data and services for researchers working on ncRNAs and their interactions with other biomolecules.

VarSite: Disease variants and protein structure

VarSite is a web server mapping known disease-associated variants from UniProt and ClinVar, together with natural variants from gnomAD, onto protein 3D structures in the Protein Data Bank. The analyses are primarily image-based and provide both an overview for each human protein, as well as a report for any specific variant of interest. The information can be useful in assessing whether a given variant might be pathogenic or benign. The structural annotations for each position in the protein include protein secondary structure, interactions with ligand, metal, DNA/RNA, or other protein, and various measures of a given variant's possible impact on the protein's function. The 3D locations of the disease-associated variants can be viewed interactively via the 3dmol.js JavaScript viewer, as well as in RasMol and PyMOL. Users can search for specific variants, or sets of variants, by providing the DNA coordinates of the base change(s) of interest. Additionally, various agglomerative analyses are given, such as the mapping of disease and natural variants onto specific Pfam or CATH domains. The server is freely accessible to all at: https://www.ebi.ac.uk/thornton-srv/databases/VarSite.

A pediatric perspective on genomics and prevention in the twenty-first century

We present evidence from diverse disciplines and populations to identify the current and emerging role of genomics in prevention from both medical and public health perspectives as well as key challenges and potential untoward consequences of increasing the role of genomics in these endeavors. We begin by comparing screening in healthy populations (newborn screening), with testing in symptomatic populations, which may incidentally identify secondary findings and at-risk relatives. Emerging evidence suggests that variants in genes subject to the reporting of secondary findings are more common than expected in patients who otherwise would not meet the criteria for testing and population testing for variants in these genes may more precisely identify discrete populations to target for various prevention strategies starting in childhood. Conversely, despite its theoretical promise, recent studies attempting to demonstrate benefits of next-generation sequencing for newborn screening have instead demonstrated numerous barriers and pitfalls to this approach. We also examine the special cases of pharmacogenomics and polygenic risk scores as examples of ways genomics can contribute to prevention amongst a broader population than that affected by rare Mendelian disease. We conclude with unresolved questions which will benefit from future investigations of the role of genomics in disease prevention.

Human gene and disease associations for clinical-genomics and precision medicine research

We are entering the era of personalized medicine in which an individual's genetic makeup will eventually determine how a doctor can tailor his or her therapy. Therefore, it is becoming critical to understand the genetic basis of common diseases, for example, which genes predispose and rare genetic variants contribute to diseases, and so on. Our study focuses on helping researchers, medical practitioners, and pharmacists in having a broad view of genetic variants that may be implicated in the likelihood of developing certain diseases. Our focus here is to create a comprehensive database with mobile access to all available, authentic and actionable genes, SNPs, and classified diseases and drugs collected from different clinical and genomics databases worldwide, including Ensembl, GenCode, ClinVar, GeneCards, DISEASES, HGMD, OMIM, GTR, CNVD, Novoseek, Swiss-Prot, LncRNADisease, Orphanet, GWAS Catalog, SwissVar, COSMIC, WHO, and FDA. We present a new cutting-edge gene-SNP-disease-drug mobile database with a smart phone application, integrating information about classified diseases and related genes, germline and somatic mutations, and drugs. Its database includes over 59 000 protein-coding and noncoding genes; over 67 000 germline SNPs and over a million somatic mutations reported for over 19 000 protein-coding genes located in over 1000 regions, published with over 3000 articles in over 415 journals available at the PUBMED; over 80 000 ICDs; over 123 000 NDCs; and over 100 000 classified gene-SNP-disease associations. We present an application that can provide new insights into the information about genetic basis of human complex diseases and contribute to assimilating genomic with phenotypic data for the availability of gene-based designer drugs, precise targeting of molecular fingerprints for tumor, appropriate drug therapy, predicting individual susceptibility to disease, diagnosis, and treatment of rare illnesses are all a few of the many transformations expected in the decade to come.

TS-GOEA: a web tool for tissue-specific gene set enrichment analysis based on gene ontology

BACKGROUND

The Gene Ontology (GO) knowledgebase is the world's largest source of information on the functions of genes. Since the beginning of GO project, various tools have been developed to perform GO enrichment analysis experiments. GO enrichment analysis has become a commonly used method of gene function analysis. Existing GO enrichment analysis tools do not consider tissue-specific information, although this information is very important to current research.

RESULTS

In this paper, we built an easy-to-use web tool called TS-GOEA that allows users to easily perform experiments based on tissue-specific GO enrichment analysis. TS-GOEA uses strict threshold statistical method for GO enrichment analysis, and provides statistical tests to improve the reliability of the analysis results. Meanwhile, TS-GOEA provides tools to compare different experimental results, which is convenient for users to compare the experimental results. To evaluate its performance, we tested the genes associated with platelet disease with TS-GOEA.

CONCLUSIONS

TS-GOEA is an effective GO analysis tool with unique features. The experimental results show that our method has better performance and provides a useful supplement for the existing GO enrichment analysis tools. TS-GOEA is available at http://120.77.47.2:5678.

Debutant iOS app and gene-disease complexities in clinical genomics and precision medicine

Background

The last decade has seen a dramatic increase in the availability of scientific data, where human-related biological databases have grown not only in count but also in volume, posing unprecedented challenges in data storage, processing, analysis, exchange, and curation. Next generation sequencing (NGS) advancements have facilitated and accelerated the process of identifying genetic variations. Adopting NGS with Whole-Genome and RNA sequencing in a diagnostic context has the potential to improve disease-risk detection in support of precision medicine and drug discovery. Several bioinformatics pipelines have been developed to strengthen variant interpretation by efficiently processing and analyzing sequence data, whereas many published results show how genomics data can be proactively incorporated into medical practices and improve utilization of clinical information. To utilize the wealth of genomics and health, there is a crucial need to generate appropriate gene-disease annotation repositories accessed through modern technology.

Results

Our focus here is to create a comprehensive database with mobile access to actionable genes and classified diseases, considered the foundation for clinical genomics and precision medicine. We present a publicly available iOS app, PAS-Gen, which invites global users to freely download it on iPhone and iPad devices, quickly adopt its easy to use interface, and search for genes and related diseases. PAS-Gen was developed using Swift, XCODE, and PHP scripting that uses Web and MySQL database servers, which includes over 59,000 protein-coding and non-coding genes, and over 90,000 classified gene-disease associations. PAS-Gen is founded on the clinical and scientific premise that easier healthcare and genomics data sharing will accelerate future medical discoveries.

Conclusions

We present a cutting-edge gene-disease database with a smart phone application, integrating information on classified diseases and related genes. The PAS-Gen app will assist researchers, medical practitioners, and pharmacists by providing a broad and view of genes that may be implicated in the likelihood of developing certain diseases. This tool with accelerate users’ abilities to understand the genetic basis of human complex diseases and by assimilating genomic and phenotypic data will support future work to identify gene-specific designer drugs, target precise molecular fingerprints for tumors, suggest appropriate drug therapies, predict individual susceptibility to disease, and diagnose and treat rare illnesses.

CAGI SickKids challenges: Assessment of phenotype and variant predictions derived from clinical and genomic data of children with undiagnosed diseases

Whole-genome sequencing (WGS) holds great potential as a diagnostic test. However, the majority of patients currently undergoing WGS lack a molecular diagnosis, largely due to the vast number of undiscovered disease genes and our inability to assess the pathogenicity of most genomic variants. The CAGI SickKids challenges attempted to address this knowledge gap by assessing state-of-the-art methods for clinical phenotype prediction from genomes. CAGI4 and CAGI5 participants were provided with WGS data and clinical descriptions of 25 and 24 undiagnosed patients from the SickKids Genome Clinic Project, respectively. Predictors were asked to identify primary and secondary causal variants. In addition, for CAGI5, groups had to match each genome to one of three disorder categories (neurologic, ophthalmologic, and connective), and separately to each patient. The performance of matching genomes to categories was no better than random but two groups performed significantly better than chance in matching genomes to patients. Two of the ten variants proposed by two groups in CAGI4 were deemed to be diagnostic, and several proposed pathogenic variants in CAGI5 are good candidates for phenotype expansion. We discuss implications for improving in silico assessment of genomic variants and identifying new disease genes.

PhenPath: a tool for characterizing biological functions underlying different phenotypes

Background

Many diseases are associated with complex patterns of symptoms and phenotypic manifestations. Parsimonious explanations aim at reconciling the multiplicity of phenotypic traits with the perturbation of one or few biological functions. For this, it is necessary to characterize human phenotypes at the molecular and functional levels, by exploiting gene annotations and known relations among genes, diseases and phenotypes. This characterization makes it possible to implement tools for retrieving functions shared among phenotypes, co-occurring in the same patient and facilitating the formulation of hypotheses about the molecular causes of the disease.

Results

We introduce PhenPath, a new resource consisting of two parts: PhenPathDB and PhenPathTOOL.

The former is a database collecting the human genes associated with the phenotypes described in Human Phenotype Ontology (HPO) and OMIM Clinical Synopses. Phenotypes are then associated with biological functions and pathways by means of NET-GE, a network-based method for functional enrichment of sets of genes. The present version considers only phenotypes related to diseases. PhenPathDB collects information for 18 OMIM Clinical synopses and 7137 HPO phenotypes, related to 4292 diseases and 3446 genes. Enrichment of Gene Ontology annotations endows some 87.7, 86.9 and 73.6% of HPO phenotypes with Biological Process, Molecular Function and Cellular Component terms, respectively. Furthermore, 58.8 and 77.8% of HPO phenotypes are also enriched for KEGG and Reactome pathways, respectively. Based on PhenPathDB, PhenPathTOOL analyzes user-defined sets of phenotypes retrieving diseases, genes and functional terms which they share. This information can provide clues for interpreting the co-occurrence of phenotypes in a patient.

Conclusions

The resource allows finding molecular features useful to investigate diseases characterized by multiple phenotypes, and by this, it can help researchers and physicians in identifying molecular mechanisms and biological functions underlying the concomitant manifestation of phenotypes. The resource is freely available at http://phenpath.biocomp.unibo.it.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5868-x) contains supplementary material, which is available to authorized users.

Active repurposing of drug candidates for melanoma based on GWAS, PheWAS and a wide range of omics data

Background

Drug repurposing is a swift, safe, and cheap drug discovery method. Melanoma disorders present low survival and high mortality rates and are challenging to diagnose and treat. Moreover, there is a high volume of worldwide investigations that are attempting to find melanoma-related genes of influence, which can be identified as responsive targets for reliable treatment.

Method

In this study, we used a wide range of data analyses to analyze over 1100 genes and proteins of influence with respect to cutaneous malignant melanoma. Our analysis included various investigational results from genome- and phenome-wide association studies (GWAS and PheWAS, respectively), biomedical, transcriptomic, and metabolomic datasets. We then researched the DrugBank for potential melanoma targets from the selected list. We excluded known melanoma targets to obtain a list of druggable proteins. We performed a precise analysis of the drugs’ pathogenesis and checked the expression profiles of the selected drugs having high associations with known anti-melanoma drugs.

Result

We found 35 drugs that interacted with 20 unique targets. These drugs appear to have high melanoma treatment potentials. We confirmed our results with previous studies and found supporting references for 30 of these drugs.

In conclusion, this investigation can be applied to various diseases for the efficient and economical repurposing of various drug compounds. For further validation, the results may be applicable for in vivo tests and clinical trials.

Protein-Protein Interaction Network Analysis Reveals Several Diseases Highly Associated with Polycystic Ovarian Syndrome

Based on clinical observations, women with polycystic ovarian syndrome (PCOS) are prone to developing several other diseases, such as metabolic and cardiovascular diseases. However, the molecular association between PCOS and these diseases remains poorly understood. Recent studies showed that the information from protein-protein interaction (PPI) network analysis are useful in understanding the disease association in detail. This study utilized this approach to deepen the knowledge on the association between PCOS and other diseases. A PPI network for PCOS was constructed using PCOS-related proteins (PCOSrp) obtained from PCOSBase. MCODE was used to identify highly connected regions in the PCOS network, known as subnetworks. These subnetworks represent protein families, where their molecular information is used to explain the association between PCOS and other diseases. Fisher's exact test and comorbidity data were used to identify PCOS-disease subnetworks. Pathway enrichment analysis was performed on the PCOS-disease subnetworks to identify significant pathways that are highly involved in the PCOS-disease associations. Migraine, schizophrenia, depressive disorder, obesity, and hypertension, along with twelve other diseases, were identified to be highly associated with PCOS. The identification of significant pathways, such as ribosome biogenesis, antigen processing and presentation, and mitophagy, suggest their involvement in the association between PCOS and migraine, schizophrenia, and hypertension.

100 Years of evolving gene-disease complexities and scientific debutants

It's been over 100 years since the word `gene' is around and progressively evolving in several scientific directions. Time-to-time technological advancements have heavily revolutionized the field of genomics, especially when it's about, e.g. triple code development, gene number proposition, genetic mapping, data banks, gene-disease maps, catalogs of human genes and genetic disorders, CRISPR/Cas9, big data and next generation sequencing, etc. In this manuscript, we present the progress of genomics from pea plant genetics to the human genome project and highlight the molecular, technical and computational developments. Studying genome and epigenome led to the fundamentals of development and progression of human diseases, which includes chromosomal, monogenic, multifactorial and mitochondrial diseases. World Health Organization has classified, standardized and maintained all human diseases, when many academic and commercial online systems are sharing information about genes and linking to associated diseases. To efficiently fathom the wealth of this biological data, there is a crucial need to generate appropriate gene annotation repositories and resources. Our focus has been how many gene-disease databases are available worldwide and which sources are authentic, timely updated and recommended for research and clinical purposes. In this manuscript, we have discussed and compared 43 such databases and bioinformatics applications, which enable users to connect, explore and, if possible, download gene-disease data.

Functional and Structural Features of Disease-Related Protein Variants

Modern sequencing technologies provide an unprecedented amount of data of single-nucleotide variations occurring in coding regions and leading to changes in the expressed protein sequences. A significant fraction of these single-residue variations is linked to disease onset and collected in public databases. In recent years, many scientific studies have been focusing on the dissection of salient features of disease-related variations from different perspectives. In this work, we complement previous analyses by updating a dataset of disease-related variations occurring in proteins with 3D structure. Within this dataset, we describe functional and structural features that can be of interest for characterizing disease-related variations, including major chemico-physical properties, the strength of association to disease of variation types, their effect on protein stability, their location on the protein structure, and their distribution in Pfam structural/functional protein models. Our results support previous findings obtained in different data sets and introduce Pfam models as possible fingerprints of patterns of disease related single-nucleotide variations.

GM604 regulates developmental neurogenesis pathways and the expression of genes associated with amyotrophic lateral sclerosis

Background

Amyotrophic lateral sclerosis (ALS) is currently an incurable disease without highly effective pharmacological treatments. The peptide drug GM604 (GM6 or Alirinetide) was developed as a candidate ALS therapy, which has demonstrated safety and good drug-like properties with a favorable pharmacokinetic profile. GM6 is hypothesized to bolster neuron survival through the multi-target regulation of developmental pathways, but mechanisms of action are not fully understood.

Methods

This study used RNA-seq to evaluate transcriptome responses in SH-SY5Y neuroblastoma cells following GM6 treatment (6, 24 and 48 h).

Results

We identified 2867 protein-coding genes with expression significantly altered by GM6 (FDR < 0.10). Early (6 h) responses included up-regulation of Notch and hedgehog signaling components, with increased expression of developmental genes mediating neurogenesis and axon growth. Prolonged GM6 treatment (24 and 48 h) altered the expression of genes contributing to cell adhesion and the extracellular matrix. GM6 further down-regulated the expression of genes associated with mitochondria, inflammatory responses, mRNA processing and chromatin organization. GM6-increased genes were located near GC-rich motifs interacting with C2H2 zinc finger transcription factors, whereas GM6-decreased genes were located near AT-rich motifs associated with helix-turn-helix homeodomain factors. Such motifs interacted with a diverse network of transcription factors encoded by GM6-regulated genes (STAT3, HOXD11, HES7, GLI1). We identified 77 ALS-associated genes with expression significantly altered by GM6 treatment (FDR < 0.10), which were known to function in neurogenesis, axon guidance and the intrinsic apoptosis pathway.

Conclusions

Our findings support the hypothesis that GM6 acts through developmental-stage pathways to influence neuron survival. Gene expression responses were consistent with neurotrophic effects, ECM modulation, and activation of the Notch and hedgehog neurodevelopmental pathways. This multifaceted mechanism of action is unique among existing ALS drug candidates and may be applicable to multiple neurodegenerative diseases.

Electronic supplementary material

The online version of this article (10.1186/s40035-018-0135-7) contains supplementary material, which is available to authorized users.