SysBiolPGWAS: simplifying post-GWAS analysis through the use of computational technologies and integration of diverse omics datasets.

MOTIVATION: Post-genome-wide association studies (pGWAS) analysis is designed to decipher the functional consequences of significant single-nucleotide polymorphisms (SNPs) in the era of GWAS. This can be translated into research insights and clinical benefits such as the effectiveness of strategies for disease screening, treatment and prevention. However, the setup of pGWAS (pGWAS) tools can be quite complicated, and it mostly requires big data. The challenge however is, scientists are required to have sufficient experience with several of these technically complex and complicated tools in order to complete the pGWAS analysis. RESULTS: We present SysBiolPGWAS, a pGWAS web application that provides a comprehensive functionality for biologists and non-bioinformaticians to conduct several pGWAS analyses to overcome the above challenges. It provides unique functionalities for analysis involving multi-omics datasets and visualization using various bioinformatics tools. SysBiolPGWAS provides access to individual pGWAS tools and a novel custom pGWAS pipeline that integrates several individual pGWAS tools and data. The SysBiolPGWAS app was developed to be a one-stop shop for pGWAS analysis. It targets researchers in the area of the human genome and performs its analysis mainly in the autosomal chromosomes. AVAILABILITY AND IMPLEMENTATION: SysBiolPGWAS web app was developed using JavaScript/TypeScript web frameworks and is available at: https://spgwas.waslitbre.org/. All codes are available in this GitHub repository https://github.com/covenant-university-bioinformatics.

Machine learning approach to gene essentiality prediction: a review.

Essential genes are critical for the growth and survival of any organism. The machine learning approach complements the experimental methods to minimize the resources required for essentiality assays. Previous studies revealed the need to discover relevant features that significantly classify essential genes, improve on the generalizability of prediction models across organisms, and construct a robust gold standard as the class label for the train data to enhance prediction. Findings also show that a significant limitation of the machine learning approach is predicting conditionally essential genes. The essentiality status of a gene can change due to a specific condition of the organism. This review examines various methods applied to essential gene prediction task, their strengths, limitations and the factors responsible for effective computational prediction of essential genes. We discussed categories of features and how they contribute to the classification performance of essentiality prediction models. Five categories of features, namely, gene sequence, protein sequence, network topology, homology and gene ontology-based features, were generated for Caenorhabditis elegans to perform a comparative analysis of their essentiality prediction capacity. Gene ontology-based feature category outperformed other categories of features majorly due to its high correlation with the genes' biological functions. However, the topology feature category provided the highest discriminatory power making it more suitable for essentiality prediction. The major limiting factor of machine learning to predict essential genes conditionality is the unavailability of labeled data for interest conditions that can train a classifier. Therefore, cooperative machine learning could further exploit models that can perform well in conditional essentiality predictions. SHORT ABSTRACT: Identification of essential genes is imperative because it provides an understanding of the core structure and function, accelerating drug targets' discovery, among other functions. Recent studies have applied machine learning to complement the experimental identification of essential genes. However, several factors are limiting the performance of machine learning approaches. This review aims to present the standard procedure and resources available for predicting essential genes in organisms, and also highlight the factors responsible for the current limitation in using machine learning for conditional gene essentiality prediction. The choice of features and ML technique was identified as an important factor to predict essential genes effectively.

Plant genomics in Africa: present and prospects.

Plants are the world's most consumed goods. They are of high economic value and bring many health benefits. In most countries in Africa, the supply and quality of food will rise to meet the growing population's increasing demand. Genomics and other biotechnology tools offer the opportunity to improve subsistence crops and medicinal herbs in the continent. Significant advances have been made in plant genomics, which have enhanced our knowledge of the molecular processes underlying both plant quality and yield. The sequencing of complex genomes of African plant species, facilitated by the continuously evolving next-generation sequencing technologies and advanced bioinformatics approaches, has provided new opportunities for crop improvement. This review summarizes the achievements of genome sequencing projects of endemic African plants in the last two decades. We also present perspectives and challenges for future plant genomic studies that will accelerate important plant breeding programs for African communities. These challenges include a lack of basic facilities, a lack of sequencing and bioinformatics facilities, and a lack of skills to design genomics studies. However, it is imperative to state that African countries have become key players in the plant genome revolution and genome derived-biotechnology. Therefore, African governments should invest in public plant genomics research and applications, establish bioinformatics platforms and training programs, and stimulate university and industry partnerships to fully deploy plant genomics, particularly in the fields of agriculture and medicine.

Assessing computational genomics skills: Our experience in the H3ABioNet African bioinformatics network.

The H3ABioNet pan-African bioinformatics network, which is funded to support the Human Heredity and Health in Africa (H3Africa) program, has developed node-assessment exercises to gauge the ability of its participating research and service groups to analyze typical genome-wide datasets being generated by H3Africa research groups. We describe a framework for the assessment of computational genomics analysis skills, which includes standard operating procedures, training and test datasets, and a process for administering the exercise. We present the experiences of 3 research groups that have taken the exercise and the impact on their ability to manage complex projects. Finally, we discuss the reasons why many H3ABioNet nodes have declined so far to participate and potential strategies to encourage them to do so.

Correction: Heuristic-enabled active machine learning: A case study of predicting essential developmental stage and immune response genes in Drosophila melanogaster.

[This corrects the article DOI: 10.1371/journal.pone.0288023.].

An in-silico analysis of OGT gene association with diabetes mellitus.

O-GlcNAcylation is a nutrient-sensing post-translational modification process. This cycling process involves two primary proteins: the O-linked N-acetylglucosamine transferase (OGT) catalysing the addition, and the glycoside hydrolase OGA (O-GlcNAcase) catalysing the removal of the O-GlCNAc moiety on nucleocytoplasmic proteins. This process is necessary for various critical cellular functions. The O-linked N-acetylglucosamine transferase (OGT) gene produces the OGT protein. Several studies have shown the overexpression of this protein to have biological implications in metabolic diseases like cancer and diabetes mellitus (DM). This study retrieved 159 SNPs with clinical significance from the SNPs database. We probed the functional effects, stability profile, and evolutionary conservation of these to determine their fit for this research. We then identified 7 SNPs (G103R, N196K, Y228H, R250C, G341V, L367F, and C845S) with predicted deleterious effects across the four tools used (PhD-SNPs, SNPs&Go, PROVEAN, and PolyPhen2). Proceeding with this, we used ROBETTA, a homology modelling tool, to model the proteins with these point mutations and carried out a structural bioinformatics method- molecular docking- using the Glide model of the Schrodinger Maestro suite. We used a previously reported inhibitor of OGT, OSMI-1, as the ligand for these mutated protein models. As a result, very good binding affinities and interactions were observed between this ligand and the active site residues within 4Å of OGT. We conclude that these mutation points may be used for further downstream analysis as drug targets for treating diabetes mellitus.

Mathematical modeling of malaria transmission dynamics in humans with mobility and control states.

Malaria importation is one of the hypothetical drivers of malaria transmission dynamics across the globe. Several studies on malaria importation focused on the effect of the use of conventional malaria control strategies as approved by the World Health Organization (WHO) on malaria transmission dynamics but did not capture the effect of the use of traditional malaria control strategies by vigilant humans. In order to handle the aforementioned situation, a novel system of Ordinary Differential Equations (ODEs) was developed comprising the human and the malaria vector compartments. Analysis of the system was carried out to assess its quantitative properties. The novel computational algorithm used to solve the developed system of ODEs was implemented and benchmarked with the existing Runge-Kutta numerical solution method. Furthermore, simulations of different vigilant conditions useful to control malaria were carried out. The novel system of malaria models was well-posed and epidemiologically meaningful based on its quantitative properties. The novel algorithm performed relatively better in terms of model simulation accuracy than Runge-Kutta. At the best model-fit condition of 98% vigilance to the use of conventional and traditional malaria control strategies, this study revealed that malaria importation has a persistent impact on malaria transmission dynamics. In lieu of this, this study opined that total vigilance to the use of the WHO-approved and traditional malaria management tools would be the most effective control strategy against malaria importation.

In silico Structure Prediction, Molecular Docking, and Dynamic Simulation of Plasmodium falciparum AP2-I Transcription Factor.

Plasmodium falciparum Apicomplexan Apetala 2 Invasion (PfAP2-I) transcription factor (TF) is a protein that regulates the expression of a subset of gene families involved in P. falciparum red blood cell (RBC) invasion. Inhibiting PfAP2-I TF with small molecules represents a potential new antimalarial therapeutic target to combat drug resistance, which this study aims to achieve. The 3D model structure of PfAP2-I was predicted ab initio using ROBETTA prediction tool and was validated using Save server 6.0 and MolProbity. Computed Atlas of Surface Topography of proteins (CASTp) 3.0 was used to predict the active sites of the PfAP2-I modeled structure. Pharmacophore modeling of the control ligand and PfAP2-I modeled structure was carried out using the Pharmit server to obtain several compounds used for molecular docking analysis. Molecular docking and postdocking studies were conducted using AutoDock vina and Discovery studio. The designed ligands' toxicity predictions and in silico drug-likeness were performed using the SwissADME predictor and OSIRIS Property Explorer. The modeled protein structure from the ROBETTA showed a validation result of 96.827 for ERRAT, 90.2% of the amino acid residues in the most favored region for the Ramachandran plot, and MolProbity score of 1.30 in the 98th percentile. Five (5) best hit compounds from molecular docking analysis were selected based on their binding affinity (between -8.9 and -11.7 Kcal/mol) to the active site of PfAP2-I and were considered for postdocking studies. For the absorption, distribution, metabolism, elimination, and toxicity (ADMET) properties, compound MCULE-7146940834 had the highest drug score (0.63) and drug-likeness (6.76). MCULE-7146940834 maintained a stable conformation within the flexible protein's active site during simulation. The good, estimated binding energies, drug-likeness, drug score, and molecular dynamics simulation interaction observed for MCULE-7146940834 against PfAP2-I show that MCULE-7146940834 can be considered a lead candidate for PfAP2-I inhibition. Experimental validations should be carried out to ascertain the efficacy of these predicted best hit compounds.

Polygenic risk scores for disease risk prediction in Africa: current challenges and future directions.

Early identification of genetic risk factors for complex diseases can enable timely interventions and prevent serious outcomes, including mortality. While the genetics underlying many Mendelian diseases have been elucidated, it is harder to predict risk for complex diseases arising from the combined effects of many genetic variants with smaller individual effects on disease aetiology. Polygenic risk scores (PRS), which combine multiple contributing variants to predict disease risk, have the potential to influence the implementation for precision medicine. However, the majority of existing PRS were developed from European data with limited transferability to African populations. Notably, African populations have diverse genetic backgrounds, and a genomic architecture with smaller haplotype blocks compared to European genomes. Subsequently, growing evidence shows that using large-scale African ancestry cohorts as discovery for PRS development may generate more generalizable findings. Here, we (1) discuss the factors contributing to the poor transferability of PRS in African populations, (2) showcase the novel Africa genomic datasets for PRS development, (3) explore the potential clinical utility of PRS in African populations, and (4) provide insight into the future of PRS in Africa.

Strengthening Bioinformatics and Genomics Analyses Skills in Africa for Attainment of the Sustainable Development Goals: Report of the 2nd Conference of the Nigerian Bioinformatics and Genomics Network.

The second conference of the Nigerian Bioinformatics and Genomics Network (NBGN21) was held from October 11 to October 13, 2021. The event was organized by the Nigerian Bioinformatics and Genomics Network. A 1-day genomic analysis workshop on genome-wide association study and polygenic risk score analysis was organized as part of the conference. It was organized primarily as a research capacity building initiative to empower Nigerian researchers to take a leading role in this cutting-edge field of genomic data science. The theme of the conference was "Leveraging Bioinformatics and Genomics for the attainments of the Sustainable Development Goals." The conference used a hybrid approach-virtual and in-person. It served as a platform to bring together 235 registered participants mainly from Nigeria and virtually, from all over the world. NBGN21 had four keynote speakers and four leading Nigerian scientists received awards for their contributions to genomics and bioinformatics development in Nigeria. A total of 100 travel fellowships were awarded to delegates within Nigeria. A major topic of discussion was the application of bioinformatics and genomics in the achievement of the Sustainable Development Goals (SDG3-Good Health and Well-Being, SDG4-Quality Education, and SDG 15-Life on Land [Biodiversity]). In closing, most of the NBGN21 conference participants were interviewed and interestingly they agreed that bioinformatics and genomic analysis of African genomes are vital in identifying population-specific genetic variants that confer susceptibility to different diseases that are endemic in Africa. The knowledge of this can empower African healthcare systems and governments for timely intervention, thereby enhancing good health and well-being.

Computational and drug target analysis of functional single nucleotide polymorphisms associated with Haemoglobin Subunit Beta (HBB) gene.

There is overwhelming evidence implicating Haemoglobin Subunit Beta (HBB) protein in the onset of beta thalassaemia. In this study for the first time, we used a combined SNP informatics and computer algorithms such as Neural network, Bayesian network, and Support Vector Machine to identify deleterious non-synonymous Single Nucleotide Polymorphisms (nsSNPs) present in the HBB gene. Our findings highlight three major mutation points (R31G, W38S, and Q128P) within the HBB gene sequence that have significant statistical and computational associations with the onset of beta thalassaemia. The dynamic simulation study revealed that R31G, W38S, and Q128P elicited high structural perturbation and instability, however, the wild type protein was considerably stable. Ten compounds with therapeutic potential against HBB were also predicted by structure-based virtual screening. Interestingly, the instability caused by the mutations was reversed upon binding to a ligand. This study has been able to predict potential deleterious mutants that can be further explored in the understanding of the pathological basis of beta thalassaemia and the design of tailored inhibitors.

The Nigerian Bioinformatics and Genomics Network (NBGN): a collaborative platform to advance bioinformatics and genomics in Nigeria.

Africa plays a central importance role in the human origins, and disease susceptibility, agriculture and biodiversity conservation. Nigeria as the most populous and most diverse country in Africa, owing to its 250 ethnic groups and over 500 different native languages is imperative to any global genomic initiative. The newly inaugurated Nigerian Bioinformatics and Genomics Network (NBGN) becomes necessary to facilitate research collaborative activities and foster opportunities for skills' development amongst Nigerian bioinformatics and genomics investigators. NBGN aims to advance and sustain the fields of genomics and bioinformatics in Nigeria by serving as a vehicle to foster collaboration, provision of new opportunities for interactions between various interdisciplinary subfields of genomics, computational biology and bioinformatics as this will provide opportunities for early career researchers. To provide the foundation for sustainable collaborations, the network organises conferences, workshops, trainings and create opportunities for collaborative research studies and internships, recognise excellence, openly share information and create opportunities for more Nigerians to develop the necessary skills to exceed in genomics and bioinformatics. NBGN currently has attracted more than 650 members around the world. Research collaborations between Nigeria, Africa and the West will grow and all stakeholders, including funding partners, African scientists, researchers across the globe, physicians and patients will be the eventual winners. The exponential membership growth and diversity of research interests of NBGN just within weeks of its establishment and the unanticipated attendance of its activities suggest the significant importance of the network to bioinformatics and genomics research in Nigeria.

Potential Anti-Cancer Flavonoids Isolated From Caesalpinia bonduc Young Twigs and Leaves: Molecular Docking and In Silico Studies.

Tyrosine kinase (TK), vascular endothelial growth factor (VEGF), and matrix metalloproteinases (MMP) are important cancer therapeutic target proteins. Based on reported anti-cancer and cytotoxic activities of Caesalpinia bonduc, this study isolated phytochemicals from young twigs and leaves of C bonduc and identified the interaction between them and cancer target proteins (TK, VEGF, and MMP) in silico. AutoDock Vina, iGEMDOCK, and analysis of pharmacokinetic and pharmacodynamic properties of the isolated bioactives as therapeutic molecules were performed. Seven phytochemicals (7-hydroxy-4'-methoxy-3,11-dehydrohomoisoflavanone, 4,4'-dihydroxy-2'-methoxy-chalcone, 7,4'-dihydroxy-3,11-dehydrohomoisoflavanone, luteolin, quercetin-3-methyl, kaempferol-3-O-ß-d-xylopyranoside and kaempferol-3-O-α-l-rhamnopyranosyl-(1 → 2)-ß-D-xylopyranoside) were isolated. Molecular docking analysis showed that the phytochemicals displayed strong interactions with the proteins compared with their respective drug inhibitors. Pharmacokinetic and pharmacodynamic properties of the compounds were promising suggesting that they can be developed as putative lead compounds for developing new anti-cancer drugs.

Computational Identification of Metabolic Pathways of Plasmodium falciparum using the k-Shortest Path Algorithm.

Plasmodium falciparum, a malaria pathogen, has shown substantial resistance to treatment coupled with poor response to some vaccines thereby requiring urgent, holistic, and broad approach to prevent this endemic disease. Understanding the biology of the malaria parasite has been identified as a vital approach to overcome the threat of malaria. This study is aimed at identifying essential proteins unique to malaria parasites using a reconstructed iPfa genome-scale metabolic model (GEM) of the 3D7 strain of Plasmodium falciparum by filling gaps in the model with nineteen (19) metabolites and twenty-three (23) reactions obtained from the MetaCyc database. Twenty (20) currency metabolites were removed from the network because they have been identified to produce shortcuts that are biologically infeasible. The resulting modified iPfa GEM was a model using the k-shortest path algorithm to identify possible alternative metabolic pathways in glycolysis and pentose phosphate pathways of Plasmodium falciparum. Heuristic function was introduced for the optimal performance of the algorithm. To validate the prediction, the essentiality of the reactions in the reconstructed network was evaluated using betweenness centrality measure, which was applied to every reaction within the pathways considered in this study. Thirty-two (32) essential reactions were predicted among which our method validated fourteen (14) enzymes already predicted in the literature. The enzymatic proteins that catalyze these essential reactions were checked for homology with the host genome, and two (2) showed insignificant similarity, making them possible drug targets. In conclusion, the application of the intelligent search technique to the metabolic network of P. falciparum predicts potential biologically relevant alternative pathways using graph theory-based approach.

In Silico Knockout Screening of Plasmodium falciparum Reactions and Prediction of Novel Essential Reactions by Analysing the Metabolic Network.

Malaria is an infectious disease that affects close to half a million individuals every year and Plasmodium falciparum is a major cause of malaria. The treatment of this disease could be done effectively if the essential enzymes of this parasite are specifically targeted. Nevertheless, the development of the parasite in resisting existing drugs now makes discovering new drugs a core responsibility. In this study, a novel computational model that makes the prediction of new and validated antimalarial drug target cheaper, easier, and faster has been developed. We have identified new essential reactions as potential targets for drugs in the metabolic network of the parasite. Among the top seven (7) predicted essential reactions, four (4) have been previously identified in earlier studies with biological evidence and one (1) has been with computational evidence. The results from our study were compared with an extensive list of seventy-seven (77) essential reactions with biological evidence from a previous study. We present a list of thirty-one (31) potential candidates for drug targets in Plasmodium falciparum which includes twenty-four (24) new potential candidates for drug targets.

Genomic Interventions in Medicine.

Lately, the term "genomics" has become ubiquitous in many scientific articles. It is a rapidly growing aspect of the biomedical sciences that studies the genome. The human genome contains a torrent of information that gives clues about human origin, evolution, biological function, and diseases. In a bid to demystify the workings of the genome, the Human Genome Project (HGP) was initiated in 1990, with the chief goal of sequencing the approximately 3 billion nucleotide base pairs of the human DNA. Since its completion in 2003, the HGP has opened new avenues for the application of genomics in clinical practice. This review attempts to overview some milestone discoveries that paved way for the initiation of the HGP, remarkable revelations from the HGP, and how genomics is influencing a paradigm shift in routine clinical practice. It further highlights the challenges facing the implementation of genomic medicine, particularly in Africa. Possible solutions are also discussed.

Development of Bioinformatics Infrastructure for Genomics Research.

BACKGROUND: Although pockets of bioinformatics excellence have developed in Africa, generally, large-scale genomic data analysis has been limited by the availability of expertise and infrastructure. H3ABioNet, a pan-African bioinformatics network, was established to build capacity specifically to enable H3Africa (Human Heredity and Health in Africa) researchers to analyze their data in Africa. Since the inception of the H3Africa initiative, H3ABioNet's role has evolved in response to changing needs from the consortium and the African bioinformatics community. OBJECTIVES: H3ABioNet set out to develop core bioinformatics infrastructure and capacity for genomics research in various aspects of data collection, transfer, storage, and analysis. METHODS AND RESULTS: Various resources have been developed to address genomic data management and analysis needs of H3Africa researchers and other scientific communities on the continent. NetMap was developed and used to build an accurate picture of network performance within Africa and between Africa and the rest of the world, and Globus Online has been rolled out to facilitate data transfer. A participant recruitment database was developed to monitor participant enrollment, and data is being harmonized through the use of ontologies and controlled vocabularies. The standardized metadata will be integrated to provide a search facility for H3Africa data and biospecimens. Because H3Africa projects are generating large-scale genomic data, facilities for analysis and interpretation are critical. H3ABioNet is implementing several data analysis platforms that provide a large range of bioinformatics tools or workflows, such as Galaxy, the Job Management System, and eBiokits. A set of reproducible, portable, and cloud-scalable pipelines to support the multiple H3Africa data types are also being developed and dockerized to enable execution on multiple computing infrastructures. In addition, new tools have been developed for analysis of the uniquely divergent African data and for downstream interpretation of prioritized variants. To provide support for these and other bioinformatics queries, an online bioinformatics helpdesk backed by broad consortium expertise has been established. Further support is provided by means of various modes of bioinformatics training. CONCLUSIONS: For the past 4 years, the development of infrastructure support and human capacity through H3ABioNet, have significantly contributed to the establishment of African scientific networks, data analysis facilities, and training programs. Here, we describe the infrastructure and how it has affected genomics and bioinformatics research in Africa.

Cluster analysis of Plasmodium RNA-seq time-course data identifies stage-specific co-regulated biological processes and regulatory elements.

In this study, we interpreted RNA-seq time-course data of three developmental stages of Plasmodium species by clustering genes based on similarities in their expression profile without prior knowledge of the gene function. Functional enrichment of clusters of upregulated genes at specific time-points reveals potential targetable biological processes with information on their timings. We identified common consensus sequences that these clusters shared as potential points of coordinated transcriptional control. Five cluster groups showed upregulated profile patterns of biological interest. This included two clusters from the Intraerythrocytic Developmental Cycle (cluster 4 = 16 genes, and cluster 9 = 32 genes), one from the sexual development stage (cluster 2 = 851 genes), and two from the gamete-fertilization stage in the mosquito host (cluster 4 = 153 genes, and cluster 9 = 258 genes). The IDC expressed the least numbers of genes with only 1448 genes showing any significant activity of the 5020 genes (~29%) in the experiment. Gene ontology (GO) enrichment analysis of these clusters revealed a total of 671 uncharacterized genes implicated in 14 biological processes and components associated with these stages, some of which are currently being investigated as drug targets in on-going research. Five putative transcription regulatory binding motifs shared by members of each cluster were also identified, one of which was also identified in a previous study by separate researchers. Our study shows stage-specific genes and biological processes that may be important in antimalarial drug research efforts. In addition, timed-coordinated control of separate processes may explain the paucity of factors in parasites.

Modeling of the Glycolysis Pathway in Plasmodium falciparum using Petri Nets.

Malaria is one of the deadly diseases, which affects a large number of the world's population. The Plasmodium falciparum parasite during erythrocyte stages produces its energy mainly through anaerobic glycolysis, with pyruvate being converted into lactate. The glycolysis metabolism in P. falci-parum is one of the important metabolic pathways of the parasite because the parasite is entirely dependent on it for energy. Also, several glycolytic enzymes have been proposed as drug targets. Petri nets (PNs) have been recognized as one of the important models for representing biological pathways. In this work, we built a qualitative PN model for the glycolysis pathway in P. falciparum and analyzed the model for its structural and quantitative properties using PN theory. From PlasmoCyc files, a total of 11 reactions were extracted; 6 of these were reversible and 5 were irreversible. These reactions were catalyzed by a total number of 13 enzymes. We extracted some of the essential reactions in the pathway using PN model, which are the possible drug targets without which the pathway cannot function. This model also helps to improve the understanding of the biological processes within this pathway.