Communicating computational workflows in a regulatory environment.

The volume of nucleic acid sequence data has exploded recently, amplifying the challenge of transforming data into meaningful information. Processing data can require an increasingly complex ecosystem of customized tools, which increases difficulty in communicating analyses in an understandable way yet is of sufficient detail to enable informed decisions or repeats. This can be of particular interest to institutions and companies communicating computations in a regulatory environment. BioCompute Objects (BCOs; an instance of pipeline documentation that conforms to the IEEE 2791-2020 standard) were developed as a standardized mechanism for analysis reporting. A suite of BCOs is presented, representing interconnected elements of a computation modeled after those that might be found in a regulatory submission but are shared publicly - in this case a pipeline designed to identify viral contaminants in biological manufacturing, such as for vaccines.

Consumption of sucralose- and acesulfame-potassium-containing diet soda alters the relative abundance of microbial taxa at the species level: findings of two pilot studies.

Sucralose and acesulfame-potassium consumption alters gut microbiota in rodents, with unclear effects in humans. We examined effects of three-times daily sucralose- and acesulfame-potassium-containing diet soda consumption for 1 (n = 17) or 8 (n = 8) weeks on gut microbiota composition in young adults. After 8 weeks of diet soda consumption, the relative abundance of Proteobacteria, specifically Enterobacteriaceae, increased; and, increased abundance of two Proteobacteria taxa was also observed after 1 week of diet soda consumption compared with sparkling water. In addition, three taxa in the Bacteroides genus increased following 1 week of diet soda consumption compared with sparkling water. The clinical relevance of these findings and effects of sucralose and acesulfame-potassium consumption on human gut microbiota warrant further investigation in larger studies. Clinical trial registration: NCT02877186 and NCT03125356.

The Glycan Structure Dictionary-a dictionary describing commonly used glycan structure terms.

Recent technological advances in glycobiology have resulted in a large influx of data and the publication of many papers describing discoveries in glycoscience. However, the terms used in describing glycan structural features are not standardized, making it difficult to harmonize data across biomolecular databases, hampering the harvesting of information across studies and hindering text mining and curation efforts. To address this shortcoming, the Glycan Structure Dictionary has been developed as a reference dictionary to provide a standardized list of widely used glycan terms that can help in the curation and mapping of glycan structures described in publications. Currently, the dictionary has 190 glycan structure terms with 297 synonyms linked to 3,332 publications. For a term to be included in the dictionary, it must be present in at least 2 peer-reviewed publications. Synonyms, annotations, and cross-references to GlyTouCan, GlycoMotif, and other relevant databases and resources are also provided when available. The purpose of this effort is to facilitate biocuration, assist in the development of text mining tools, improve the harmonization of search, and browse capabilities in glycoinformatics resources and help to map glycan structures to function and disease. It is also expected that authors will use these terms to describe glycan structures in their manuscripts over time. A mechanism is also provided for researchers to submit terms for potential incorporation. The dictionary is available at https://wiki.glygen.org/Glycan_structure_dictionary.

Worldwide Glycoscience Informatics Infrastructure: The GlySpace Alliance.

The GlySpace Alliance was formed in 2018 among the principal investigators of three major glycoscience portals: Glyco@Expasy, GlyCosmos, and GlyGen, representing Europe, Asia, and the United States, respectively. While each of these portals has its unique user interface, the aim is to provide the same basic data set of glycan-related omics data. These portals will be introduced with the aim to enable users to find their target information in the most efficient manner, in particular, in terms of the chemical structures of glycans and their functions.

Multi-omics approaches to studying gastrointestinal microbiome in the context of precision medicine and machine learning.

The human gastrointestinal (gut) microbiome plays a critical role in maintaining host health and has been increasingly recognized as an important factor in precision medicine. High-throughput sequencing technologies have revolutionized -omics data generation, facilitating the characterization of the human gut microbiome with exceptional resolution. The analysis of various -omics data, including metatranscriptomics, metagenomics, glycomics, and metabolomics, holds potential for personalized therapies by revealing information about functional genes, microbial composition, glycans, and metabolites. This multi-omics approach has not only provided insights into the role of the gut microbiome in various diseases but has also facilitated the identification of microbial biomarkers for diagnosis, prognosis, and treatment. Machine learning algorithms have emerged as powerful tools for extracting meaningful insights from complex datasets, and more recently have been applied to metagenomics data via efficiently identifying microbial signatures, predicting disease states, and determining potential therapeutic targets. Despite these rapid advancements, several challenges remain, such as key knowledge gaps, algorithm selection, and bioinformatics software parametrization. In this mini-review, our primary focus is metagenomics, while recognizing that other -omics can enhance our understanding of the functional diversity of organisms and how they interact with the host. We aim to explore the current intersection of multi-omics, precision medicine, and machine learning in advancing our understanding of the gut microbiome. A multidisciplinary approach holds promise for improving patient outcomes in the era of precision medicine, as we unravel the intricate interactions between the microbiome and human health.

Modeling and integration of N-glycan biomarkers in a comprehensive biomarker data model.

Molecular biomarkers measure discrete components of biological processes that can contribute to disorders when impaired. Great interest exists in discovering early cancer biomarkers to improve outcomes. Biomarkers represented in a standardized data model, integrated with multi-omics data, may improve the understanding and use of novel biomarkers such as glycans and glycoconjugates. Among altered components in tumorigenesis, N-glycans exhibit substantial biomarker potential, when analyzed with their protein carriers. However, such data are distributed across publications and databases of diverse formats, which hamper their use in research and clinical application. Mass spectrometry measures of 50 N-glycans on 7 serum proteins in liver disease were integrated (as a panel) into a cancer biomarker data model, providing a unique identifier, standard nomenclature, links to glycan resources, and accession and ontology annotations to standard protein, gene, disease, and biomarker information. Data provenance was documented with a standardized United States Food and Drug Administration-supported BioCompute Object. Using the biomarker data model allows the capture of granular information, such as glycans with different levels of abundance in cirrhosis, hepatocellular carcinoma, and transplant groups. Such representation in a standardized data model harmonizes glycomics data in a unified framework, making glycan-protein biomarker data exploration more available to investigators and to other data resources. The biomarker data model we describe can be used by researchers to describe their novel glycan and glycoconjugate biomarkers; it can integrate N-glycan biomarker data with multi-source biomedical data and can foster discovery and insight within a unified data framework for glycan biomarker representation, thereby making the data FAIR (Findable, Accessible, Interoperable, Reusable) (https://www.go-fair.org/fair-principles/).

Cardiomyocyte-specific regression of nitrosative stress-mediated S-Nitrosylation of IKKγ alleviates pathological cardiac hypertrophy.

IKKγ prototypically promotes NFκBp65 activity by regulating the assembly of the IKK holocomplex. In hypertrophied cardiomyocytes, the p65-p300 complex-induced regenerative efforts are neutralized by the p53-p300 complex-mediated apoptotic load resulting in compromised cardiac function. The present study reports that nitrosative stress leads to S-Nitrosylation of IKKγ in hypertrophied cardiomyocytes in a pre-clinical model. Using a cardiomyocyte-targeted nanoconjugate, IKKγ S-Nitrosylation-resistant mutant plasmids were delivered to the pathologically hypertrophied heart that resulted in improved cardiac function by amelioration of cardiomyocyte apoptosis and simultaneous induction of their cell cycle re-entry machinery. Mechanistically, in IKKγ S-Nitrosyl mutant-transfected hypertrophied cells, increased IKKγ-p300 binding downregulated the binding of p53 and p65 with p300. This shifted the binding preference of p65 from p300 to HDAC1 resulting in upregulated expression of cyclin D1 and CDK2 via the p27/pRb pathway. This approach has therapeutic advantage over mainstream anti-hypertrophic remedies which concomitantly reduce the regenerative prowess of resident cardiomyocytes during hypertrophy upon downregulation of myocyte apoptosis. Therefore, cardiomyocyte-targeted delivery of IKKγ S-Nitrosyl mutants during hypertrophy can be exploited as a novel strategy to re-muscularize the diseased heart.

Higher levels of Bifidobacteria and tumor necrosis factor in children with drug-resistant epilepsy are associated with anti-seizure response to the ketogenic diet.

BACKGROUND: Recently, studies have suggested a role for the gut microbiota in epilepsy. Gut microbial changes during ketogenic diet (KD) treatment of drug-resistant epilepsy have been described. Inflammation is associated with certain types of epilepsy and specific inflammation markers decrease during KD. The gut microbiota plays an important role in the regulation of the immune system and inflammation. METHODS: 28 children with drug-resistant epilepsy treated with the ketogenic diet were followed in this observational study. Fecal and serum samples were collected at baseline and three months after dietary intervention. FINDINGS: We identified both gut microbial and inflammatory changes during treatment. KD had a general anti-inflammatory effect. Novel bioinformatics and machine learning approaches identified signatures of specific Bifidobacteria and TNF (tumor necrosis factor) associated with responders before starting KD. During KD, taxonomic and inflammatory profiles between responders and non-responders were more similar than at baseline. INTERPRETATION: Our results suggest that children with drug-resistant epilepsy are more likely to benefit from KD treatment when specific Bifidobacteria and TNF are elevated. We here present a novel signature of interaction of the gut microbiota and the immune system associated with anti-epileptic response to KD treatment. This signature could be used as a prognostic biomarker to identify potential responders to KD before starting treatment. Our findings may also contribute to the development of new anti-seizure therapies by targeting specific components of the gut microbiota. FUNDING: This study was supported by the Swedish Brain Foundation, Margarethahemmet Society, Stiftelsen Sunnerdahls Handikappfond, Linnea & Josef Carlssons Foundation, and The McCormick Genomic & Proteomic Center.

Communicating regulatory high-throughput sequencing data using BioCompute Objects.

This project demonstrates the use of the IEEE 2791-2020 Standard (BioCompute Objects [BCO]) to enable the complete and concise communication of results from next generation sequencing (NGS) analysis. One arm of a clinical trial was replicated using synthetically generated data made to resemble real biological data and then two independent analyses were performed. The first simulated a pharmaceutical regulatory submission to the US Food and Drug Administration (FDA) including analysis of results and a BCO. The second simulated an FDA review that included an independent analysis of the submitted data. Of the 118 simulated patient samples generated, 117 (99.15%) were in agreement in the two analyses. This process exemplifies how a template BCO (tBCO), including a verification kit, facilitates transparency and reproducibility, thereby reinforcing confidence in the regulatory submission process.

Whole Genome Variant Dataset for Enriching Studies across 18 Different Cancers.

Whole genome sequencing (WGS) has helped to revolutionize biology, but the computational challenge remains for extracting valuable inferences from this information. Here, we present the cancer-associated variants from the Cancer Genome Atlas (TCGA) WGS dataset. This set of data will allow cancer researchers to further expand their analysis beyond the exomic regions of the genome to the entire genome. A total of 1342 WGS alignments available from the consortium were processed with VarScan2 and deposited to the NCI Cancer Cloud. The sample set covers 18 different cancers and reveals 157,313,519 pooled (non-unique) cancer-associated single-nucleotide variations (SNVs) across all samples. There was an average of 117,223 SNVs per sample, with a range from 1111 to 775,470 and a standard deviation of 163,273. The dataset was incorporated into BigQuery, which allows for fast access and cross-mapping, which will allow researchers to enrich their current studies with a plethora of newly available genomic data.

Enhancing the interoperability of glycan data flow between ChEBI, PubChem and GlyGen.

Glycans play a vital role in health, disease, bioenergy, biomaterials and bio-therapeutics. As a result, there is keen interest to identify and increase glycan data in bioinformatics databases like ChEBI and PubChem, and connecting them to resources at the EMBL-EBI and NCBI to facilitate access to important annotations at a global level. GlyTouCan is a comprehensive archival database that contains glycans obtained primarily through batch upload from glycan repositories, glycoprotein databases and individual laboratories. In many instances, the glycan structures deposited in GlyTouCan may not be fully defined or have supporting experimental evidence and citations. Databases like ChEBI and PubChem were designed to accommodate complete atomistic structures with well-defined chemical linkages. As a result, they cannot easily accommodate the structural ambiguity inherent in glycan databases. Consequently, there is a need to improve the organization of glycan data coherently to enhance connectivity across the major NCBI, EMBL-EBI and glycoscience databases. This paper outlines a workflow developed in collaboration between GlyGen, ChEBI and PubChem to improve the visibility and connectivity of glycan data across these resources. GlyGen hosts a subset of glycans (~29,000) from the GlyTouCan database and has submitted valuable glycan annotations to the PubChem database and integrated over 10,500 (including ambiguously defined) glycans into the ChEBI database. The integrated glycans were prioritized based on links to PubChem and connectivity to glycoprotein data. The pipeline provides a blueprint for how glycan data can be harmonized between different resources. The current PubChem, ChEBI and GlyTouCan mappings can be downloaded from GlyGen (https://data.glygen.org).

Mice with dysfunctional TGF-ß signaling develop altered intestinal microbiome and colorectal cancer resistant to 5FU.

Emerging data show a rise in colorectal cancer (CRC) incidence in young men and women that is often chemoresistant. One potential risk factor is an alteration in the microbiome. Here, we investigated the role of TGF-ß signaling on the intestinal microbiome and the efficacy of chemotherapy for CRC induced by azoxymethane and dextran sodium sulfate in mice. We used two genotypes of TGF-ß-signaling-deficient mice (Smad4+/- and Smad4+/-Sptbn1+/-), which developed CRC with similar phenotypes and had similar alterations in the intestinal microbiome. Using these mice, we evaluated the intestinal microbiome and determined the effect of dysfunctional TGF-ß signaling on the response to the chemotherapeutic agent 5-Fluoro-uracil (5FU) after induction of CRC. Using shotgun metagenomic sequencing, we determined gut microbiota composition in mice with CRC and found reduced amounts of beneficial species of Bacteroides and Parabacteroides in the mutants compared to the wild-type (WT) mice. Furthermore, the mutant mice with CRC were resistant to 5FU. Whereas the abundances of E. boltae, B.dorei, Lachnoclostridium sp., and Mordavella sp. were significantly reduced in mice with CRC, these species only recovered to basal amounts after 5FU treatment in WT mice, suggesting that the alterations in the intestinal microbiome resulting from compromised TGF-ß signaling impaired the response to 5FU. These findings could have implications for inhibiting the TGF-ß pathway in the treatment of CRC or other cancers.

COVID-19 biomarkers and their overlap with comorbidities in a disease biomarker data model.

In response to the COVID-19 outbreak, scientists and medical researchers are capturing a wide range of host responses, symptoms and lingering postrecovery problems within the human population. These variable clinical manifestations suggest differences in influential factors, such as innate and adaptive host immunity, existing or underlying health conditions, comorbidities, genetics and other factors-compounding the complexity of COVID-19 pathobiology and potential biomarkers associated with the disease, as they become available. The heterogeneous data pose challenges for efficient extrapolation of information into clinical applications. We have curated 145 COVID-19 biomarkers by developing a novel cross-cutting disease biomarker data model that allows integration and evaluation of biomarkers in patients with comorbidities. Most biomarkers are related to the immune (SAA, TNF-∝ and IP-10) or coagulation (D-dimer, antithrombin and VWF) cascades, suggesting complex vascular pathobiology of the disease. Furthermore, we observe commonality with established cancer biomarkers (ACE2, IL-6, IL-4 and IL-2) as well as biomarkers for metabolic syndrome and diabetes (CRP, NLR and LDL). We explore these trends as we put forth a COVID-19 biomarker resource (https://data.oncomx.org/covid19) that will help researchers and diagnosticians alike.

Bioinformatics tools developed to support BioCompute Objects.

Developments in high-throughput sequencing (HTS) result in an exponential increase in the amount of data generated by sequencing experiments, an increase in the complexity of bioinformatics analysis reporting and an increase in the types of data generated. These increases in volume, diversity and complexity of the data generated and their analysis expose the necessity of a structured and standardized reporting template. BioCompute Objects (BCOs) provide the requisite support for communication of HTS data analysis that includes support for workflow, as well as data, curation, accessibility and reproducibility of communication. BCOs standardize how researchers report provenance and the established verification and validation protocols used in workflows while also being robust enough to convey content integration or curation in knowledge bases. BCOs that encapsulate tools, platforms, datasets and workflows are FAIR (findable, accessible, interoperable and reusable) compliant. Providing operational workflow and data information facilitates interoperability between platforms and incorporation of future dataset within an HTS analysis for use within industrial, academic and regulatory settings. Cloud-based platforms, including High-performance Integrated Virtual Environment (HIVE), Cancer Genomics Cloud (CGC) and Galaxy, support BCO generation for users. Given the 100K+ userbase between these platforms, BioCompute can be leveraged for workflow documentation. In this paper, we report the availability of platform-dependent and platform-independent BCO tools: HIVE BCO App, CGC BCO App, Galaxy BCO API Extension and BCO Portal. Community engagement was utilized to evaluate tool efficacy. We demonstrate that these tools further advance BCO creation from text editing approaches used in earlier releases of the standard. Moreover, we demonstrate that integrating BCO generation within existing analysis platforms greatly streamlines BCO creation while capturing granular workflow details. We also demonstrate that the BCO tools described in the paper provide an approach to solve the long-standing challenge of standardizing workflow descriptions that are both human and machine readable while accommodating manual and automated curation with evidence tagging. Database URL:  https://www.biocomputeobject.org/resources.

Streamlined Subpopulation, Subtype, and Recombination Analysis of HIV-1 Half-Genome Sequences Generated by High-Throughput Sequencing.

High-throughput sequencing (HTS) has been widely used to characterize HIV-1 genome sequences. There are no algorithms currently that can directly determine genotype and quasispecies population using short HTS reads generated from long genome sequences without additional software. To establish a robust subpopulation, subtype, and recombination analysis workflow, we amplified the HIV-1 3'-half genome from plasma samples of 65 HIV-1-infected individuals and sequenced the entire amplicon (∼4,500 bp) by HTS. With direct analysis of raw reads using HIVE-hexahedron, we showed that 48% of samples harbored 2 to 13 subpopulations. We identified various subtypes (17 A1s, 4 Bs, 27 Cs, 6 CRF02_AGs, and 11 unique recombinant forms) and defined recombinant breakpoints of 10 recombinants. These results were validated with viral genome sequences generated by single genome sequencing (SGS) or the analysis of consensus sequence of the HTS reads. The HIVE-hexahedron workflow is more sensitive and accurate than just evaluating the consensus sequence and also more cost-effective than SGS.IMPORTANCE The highly recombinogenic nature of human immunodeficiency virus type 1 (HIV-1) leads to recombination and emergence of quasispecies. It is important to reliably identify subpopulations to understand the complexity of a viral population for drug resistance surveillance and vaccine development. High-throughput sequencing (HTS) provides improved resolution over Sanger sequencing for the analysis of heterogeneous viral subpopulations. However, current methods of analysis of HTS reads are unable to fully address accurate population reconstruction. Hence, there is a dire need for a more sensitive, accurate, user-friendly, and cost-effective method to analyze viral quasispecies. For this purpose, we have improved the HIVE-hexahedron algorithm that we previously developed with in silico short sequences to analyze raw HTS short reads. The significance of this study is that our standalone algorithm enables a streamlined analysis of quasispecies, subtype, and recombination patterns from long HIV-1 genome regions without the need of additional sequence analysis tools. Distinct viral populations and recombination patterns identified by HIVE-hexahedron are further validated by comparison with sequences obtained by single genome sequencing (SGS).

GlyGen data model and processing workflow.

SUMMARY: Glycoinformatics plays a major role in glycobiology research, and the development of a comprehensive glycoinformatics knowledgebase is critical. This application note describes the GlyGen data model, processing workflow and the data access interfaces featuring programmatic use case example queries based on specific biological questions. The GlyGen project is a data integration, harmonization and dissemination project for carbohydrate and glycoconjugate-related data retrieved from multiple international data sources including UniProtKB, GlyTouCan, UniCarbKB and other key resources. AVAILABILITY AND IMPLEMENTATION: GlyGen web portal is freely available to access at https://glygen.org. The data portal, web services, SPARQL endpoint and GitHub repository are also freely available at https://data.glygen.org, https://api.glygen.org, https://sparql.glygen.org and https://github.com/glygener, respectively. All code is released under license GNU General Public License version 3 (GNU GPLv3) and is available on GitHub https://github.com/glygener. The datasets are made available under Creative Commons Attribution 4.0 International (CC BY 4.0) license. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

OncoMX: A Knowledgebase for Exploring Cancer Biomarkers in the Context of Related Cancer and Healthy Data.

PURPOSE: The purpose of OncoMX1 knowledgebase development was to integrate cancer biomarker and relevant data types into a meta-portal, enabling the research of cancer biomarkers side by side with other pertinent multidimensional data types. METHODS: Cancer mutation, cancer differential expression, cancer expression specificity, healthy gene expression from human and mouse, literature mining for cancer mutation and cancer expression, and biomarker data were integrated, unified by relevant biomedical ontologies, and subjected to rule-based automated quality control before ingestion into the database. RESULTS: OncoMX provides integrated data encompassing more than 1,000 unique biomarker entries (939 from the Early Detection Research Network [EDRN] and 96 from the US Food and Drug Administration) mapped to 20,576 genes that have either mutation or differential expression in cancer. Sentences reporting mutation or differential expression in cancer were extracted from more than 40,000 publications, and healthy gene expression data with samples mapped to organs are available for both human genes and their mouse orthologs. CONCLUSION: OncoMX has prioritized user feedback as a means of guiding development priorities. By mapping to and integrating data from several cancer genomics resources, it is hoped that OncoMX will foster a dynamic engagement between bioinformaticians and cancer biomarker researchers. This engagement should culminate in a community resource that substantially improves the ability and efficiency of exploring cancer biomarker data and related multidimensional data.

Bioinformatics and machine learning in gastrointestinal microbiome research and clinical application.

The scientific community currently defines the human microbiome as all the bacteria, viruses, fungi, archaea, and eukaryotes that occupy the human body. When considering the variable locations, composition, diversity, and abundance of our microbial symbionts, the sheer volume of microorganisms reaches hundreds of trillions. With the onset of next generation sequencing (NGS), also known as high-throughput sequencing (HTS) technologies, the barriers to studying the human microbiome lowered significantly, making in-depth microbiome research accessible. Certain locations on the human body, such as the gastrointestinal, oral, nasal, and skin microbiomes have been heavily studied through community-focused projects like the Human Microbiome Project (HMP). In particular, the gastrointestinal microbiome (GM) has received significant attention due to links to neurological, immunological, and metabolic diseases, as well as cancer. Though HTS technologies allow deeper exploration of the GM, data informing the functional characteristics of microbiota and resulting effects on human function or disease are still sparse. This void is compounded by microbiome variability observed among humans through factors like genetics, environment, diet, metabolic activity, and even exercise; making GM research inherently difficult to study. This chapter describes an interdisciplinary approach to GM research with the goal of mitigating the hindrances of translating findings into a clinical setting. By applying tools and knowledge from microbiology, metagenomics, bioinformatics, machine learning, predictive modeling, and clinical study data from children with treatment-resistant epilepsy, we describe a proof-of-concept approach to clinical translation and precision application of GM research.