A Perspective on Protein Structure Prediction Using Quantum Computers

Despite the recent advancements by deep learning methods such as AlphaFold2, in silico protein structure prediction remains a challenging problem in biomedical research. With the rapid evolution of quantum computing, it is natural to ask whether quantum computers can offer some meaningful benefits for approaching this problem. Yet, identifying specific problem instances amenable to quantum advantage and estimating the quantum resources required are equally challenging tasks. Here, we share our perspective on how to create a framework for systematically selecting protein structure prediction problems that are amenable for quantum advantage, and estimate quantum resources for such problems on a utility-scale quantum computer. As a proof-of-concept, we validate our problem selection framework by accurately predicting the structure of a catalytic loop of the Zika Virus NS3 Helicase, on quantum hardware.

Sphingolipid biosynthesis is essential for metabolic rewiring during TH17 cell differentiation

T helper 17 (TH17) cells are implicated in autoimmune diseases, and several metabolic processes are shown to be important for their development and function. In this study, we report an essential role for sphingolipids synthesized through the de novo pathway in TH17 cell development. Deficiency of SPTLC1, a major subunit of serine palmitoyl transferase enzyme complex that catalyzes the first and rate-limiting step of de novo sphingolipid synthesis, impaired glycolysis in differentiating TH17 cells by increasing intracellular reactive oxygen species (ROS) through enhancement of nicotinamide adenine dinucleotide phosphate oxidase 2 activity. Increased ROS leads to impaired activation of mammalian target of rapamycin C1 and reduced expression of hypoxia-inducible factor 1-alpha and c-Myc-induced glycolytic genes. SPTLCI deficiency protected mice from developing experimental autoimmune encephalomyelitis and experimental T cell transfer colitis. Our results thus show a critical role for de novo sphingolipid biosynthetic pathway in shaping adaptive immune responses with implications in autoimmune diseases.

Contemporary Antiretroviral Therapy Dysregulates Iron Transport and Augments Mitochondrial Dysfunction in HIV-Infected Human Microglia and Neural-Lineage Cells

HIV-associated cognitive dysfunction during combination antiretroviral therapy (cART) involves mitochondrial dysfunction, but the impact of contemporary cART on chronic metabolic changes in the brain and in latent HIV infection is unclear. We interrogated mitochondrial function in a human microglia (hμglia) cell line harboring inducible HIV provirus and in SH-SY5Y cells after exposure to individual antiretroviral drugs or cART, using the MitoStress assay. cART-induced changes in protein expression, reactive oxygen species (ROS) production, mitochondrial DNA copy number, and cellular iron were also explored. Finally, we evaluated the ability of ROS scavengers or plasmid-mediated overexpression of the antioxidant iron-binding protein, Fth1, to reverse mitochondrial defects. Contemporary antiretroviral drugs, particularly bictegravir, depressed multiple facets of mitochondrial function by 20-30%, with the most pronounced effects in latently infected HIV+ hμglia and SH-SY5Y cells. Latently HIV-infected hμglia exhibited upregulated glycolysis. Increases in total and/or mitochondrial ROS, mitochondrial DNA copy number, and cellular iron accompanied mitochondrial defects in hμglia and SH-SY5Y cells. In SH-SY5Y cells, cART reduced mitochondrial iron-sulfur-cluster-containing supercomplex and subunit expression and increased Nox2 expression. Fth1 overexpression or pre-treatment with N-acetylcysteine prevented cART-induced mitochondrial dysfunction. Contemporary cART impairs mitochondrial bioenergetics in hμglia and SH-SY5Y cells, partly through cellular iron accumulation; some effects differ by HIV latency.

Design and Synthesis of Thiazole Scaffold-Based Small Molecules as Anticancer Agents Targeting the Human Lactate Dehydrogenase A Enzyme

A new series of thiazole central scaffold-based small molecules of hLDHA inhibitors were designed using an in silico approach. Molecular docking analysis of designed molecules with hLDHA (PDB ID: 1I10) demonstrates that Ala 29, Val 30, Arg 98, Gln 99, Gly 96, and Thr 94 possessed strong interaction with the compounds. Compounds 8a, 8b, and 8d showed good binding affinity (-8.1 to -8.8 kcal/mol), whereas an additional interaction of NO₂ at the ortho position in compounds 8c with Gln 99 through hydrogen bonding enhanced the affinity to -9.8 kcal/mol. Selected high-scored compounds were synthesized and screened for hLDHA inhibitory activities and in vitro anticancer activity in six cancer cell lines. Biochemical enzyme inhibition assays showed the highest hLDHA inhibitory activity observed with compounds 8b, 8c, and 8l. Compounds 8b, 8c, 8j, 8l, and 8m depicted significant anticancer activities, exhibiting IC₅₀ values in the range of 1.65-8.60 μM in HeLa and SiHa cervical cancer cell lines. Compounds 8j and 8m exhibited notable anticancer activity with IC₅₀ values of 7.90 and 5.15 μM, respectively, in liver cancer cells (HepG2). Interestingly, compounds 8j and 8m did not induce noticeable toxicity in the human embryonic kidney cells (HEK293). Insilico absorption, distribution, metabolism, and excretion profiling demonstrates that the compounds possess drug-likeness, and results may pave the way for the development of novel thiazole-based biologically active small molecules for therapeutics.

CNV-ClinViewer: enhancing the clinical interpretation of large copy-number variants online

MOTIVATION

Pathogenic copy-number variants (CNVs) can cause a heterogeneous spectrum of rare and severe disorders. However, most CNVs are benign and are part of natural variation in human genomes. CNV pathogenicity classification, genotype-phenotype analyses, and therapeutic target identification are challenging and time-consuming tasks that require the integration and analysis of information from multiple scattered sources by experts.

RESULTS

Here, we introduce the CNV-ClinViewer, an open-source web application for clinical evaluation and visual exploration of CNVs. The application enables real-time interactive exploration of large CNV datasets in a user-friendly designed interface and facilitates semi-automated clinical CNV interpretation following the ACMG guidelines by integrating the ClassifCNV tool. In combination with clinical judgment, the application enables clinicians and researchers to formulate novel hypotheses and guide their decision-making process. Subsequently, the CNV-ClinViewer enhances for clinical investigators' patient care and for basic scientists' translational genomic research.

AVAILABILITY AND IMPLEMENTATION

The web application is freely available at https://cnv-ClinViewer.broadinstitute.org and the open-source code can be found at https://github.com/LalResearchGroup/CNV-clinviewer.

Phosphorylation-mediated regulation of the Bacillus anthracis phosphoglycerate mutase by the Ser/Thr protein kinase PrkC

Bacillus anthracis Ser/Thr protein kinase PrkC is necessary for phenotypic memory and spore germination, and the loss of PrkC-dependent phosphorylation events affect the spore development. During sporulation, Bacillus sp. can store 3-Phosphoglycerate (3-PGA) that will be required at the onset of germination when ATP will be necessary. The Phosphoglycerate mutase (Pgm) catalyzes the isomerization of 2-PGA and 3-PGA and is important for spore germination as a key metabolic enzyme that maintains 3-PGA pool at later events. Therefore, regulation of Pgm is important for an efficient spore germination process and metabolic switching. While the increased expression of Pgm in B. anthracis decreases spore germination efficiency, it remains unexplored if PrkC could directly influence Pgm activity. Here, we report the phosphorylation and regulation of Pgm by PrkC and its impact on Pgm stability and catalytic activity. Mass spectrometry revealed Pgm phosphorylation on seven threonine residues. In silico mutational analysis highlighted the role of Thr⁴⁵⁹ residue towards metal and substrate binding. Altogether, we demonstrated that PrkC-mediated Pgm phosphorylation negatively regulates its activity that is essential to maintain Pgm in its apo-like isoform before germination. This study advances the role of Pgm regulation that represents an important switch for B. anthracis resumption of metabolism and spore germination.

Genomics to Notebook (g2nb): extending the electronic notebook to address the challenges of bioinformatics analysis

We present Genomics to Notebook (g2nb), an environment that combines the JupyterLab notebook system with widely-used bioinformatics platforms. Galaxy, GenePattern, and the JavaScript versions of IGV and Cytoscape are currently available within g2nb. The analyses and visualizations within those platforms are presented as cells in a notebook, making thousands of genomics methods available within the notebook metaphor and allowing notebooks to contain workflows utilizing multiple software packages on remote servers, all without the need for programming. The g2nb environment is, to our knowledge, the only notebook-based system that incorporates multiple bioinformatics analysis platforms into a notebook interface.

Abstract 2073: Genomics to Notebook (g2nb): Extending the electronic notebook to address the needs of cancer bioinformatics

As the availability of genomic data and analysis tools from large-scale cancer initiatives continues to increase, with single-cell studies adding new dimensions to the potential scientific insights, the need has become more urgent for a software environment that supports the rapid pace of cancer data science. The The Jupyter Notebook environment has become the de facto medium for this purpose due to its ease in combining scientific exposition with executable code to form a single reproducible “research narrative” document. However, analyses are often compute-intensive, requiring more resources than are frequently available within a notebook environment running on a desktop or laptop computer. Additionally, thousands of tools, modules, and plugins are readily available on integrative software platforms such as Galaxy, GenePattern, and Cytoscape, outside the notebook paradigm. Finally, many biomedical investigators lack the programming expertise required to fully realize the benefits and utility of the notebook metaphor. To address these issues, we have released Genomics to Notebook (g2nb), which builds on JupyterLab to add access to bioinformatics platforms and other functionality for the non-programmer through the components described below, while retaining all programmatic features of JupyterLab.The g2nb environment incorporates bioinformatics software platforms within the notebook interface, allowing a single notebook to contain a workflow spanning multiple tools and servers. When run, the entire analysis appears to execute seamlessly within the notebook. To achieve this, we developed a new analysis cell type that provides an interface within the notebook to tools that are hosted on a remote Galaxy or GenePattern server. Analysis cells present a web form-like interface, similar to that of the original platforms, requiring an investigator to provide only the input parameters and data. The popular visualization tools Cytoscape and Integrative Genomics Viewer (IGV) are also supported in their web-based formats as notebook cells, with additional platforms and visualizers added regularly. The g2nb environment is freely available at the g2nb workspace, http://g2nb.org, where scientists can use all of the g2nb functionality with only a web browser. Those who wish to use g2nb locally can use the provided Docker container or install the packages via the conda or pip package managers. The online workspace also includes a library of featured genomic analysis notebooks, including templates for common analysis tasks as well as cancer-specific research scenarios and compute-intensive methods. Scientists can easily copy these notebooks, use them as is, or adapt them for their research purposes.

Citation Format: Michael M. Reich, Thorin Tabor, John Liefeld, Jayadev Joshi, Forrest Kim, Helga Thorvaldsdottir, Daniel Blankenberg, Jill P. Mesirov. Genomics to Notebook (g2nb): Extending the electronic notebook to address the needs of cancer bioinformatics [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 2073.

Evaluating homophily in networks via HONTO (HOmophily Network TOol): a case study of chromosomal interactions in human PPI networks

SUMMARY

It has been observed in different kinds of networks, such as social or biological ones, a typical behavior inspired by the general principle 'similarity breeds connections'. These networks are defined as homophilic as nodes belonging to the same class preferentially interact with each other. In this work, we present HONTO (HOmophily Network TOol), a user-friendly open-source Python3 package designed to evaluate and analyze homophily in complex networks. The tool takes in input from the network along with a partition of its nodes into classes and yields a matrix whose entries are the homophily/heterophily z-score values. To complement the analysis, the tool also provides z-score values of nodes that do not interact with any other node of the same class. Homophily/heterophily z-scores values are presented as a heatmap allowing a visual at-a-glance interpretation of results.

AVAILABILITY AND IMPLEMENTATION

Tool's source code is available at https://github.com/cumbof/honto under the MIT license, installable as a package from PyPI (pip install honto) and conda-forge (conda install -c conda-forge honto), and has a wrapper for the Galaxy platform available on the official Galaxy ToolShed (Blankenberg et al., 2014) at https://toolshed.g2.bx.psu.edu/view/fabio/honto.

Analysis of nested alternate open reading frames and their encoded proteins

Transcriptional and post-transcriptional mechanisms diversify the proteome beyond gene number, while maintaining a sequence relationship between original and altered proteins. A new mechanism breaks this paradigm, generating novel proteins by translating alternative open reading frames (Alt-ORFs) within canonical host mRNAs. Uniquely, ‘alt-proteins’ lack sequence homology with host ORF-derived proteins. We show global amino acid frequencies, and consequent biochemical characteristics of Alt-ORFs nested within host ORFs (nAlt-ORFs), are genetically-driven, and predicted by summation of frequencies of hundreds of encompassing host codon-pairs. Analysis of 101 human nAlt-ORFs of length ≥150 codons confirms the theoretical predictions, revealing an extraordinarily high median isoelectric point (pI) of 11.68, due to anomalous charged amino acid levels. Also, nAlt-ORF proteins exhibit a >2-fold preference for reading frame 2 versus 3, predicted mitochondrial and nuclear localization, and elevated codon adaptation index indicative of natural selection. Our results provide a theoretical and conceptual framework for exploration of these largely unannotated, but potentially significant, alternative ORFs and their encoded proteins.

Delivery of ceramide phosphoethanolamine lipids to the cleavage furrow through the endocytic pathway is essential for male meiotic cytokinesis

Cell division, wherein 1 cell divides into 2 daughter cells, is fundamental to all living organisms. Cytokinesis, the final step in cell division, begins with the formation of an actomyosin contractile ring, positioned midway between the segregated chromosomes. Constriction of the ring with concomitant membrane deposition in a specified spatiotemporal manner generates a cleavage furrow that physically separates the cytoplasm. Unique lipids with specific biophysical properties have been shown to localize to intercellular bridges (also called midbody) connecting the 2 dividing cells; however, their biological roles and delivery mechanisms remain largely unknown. In this study, we show that ceramide phosphoethanolamine (CPE), the structural analog of sphingomyelin, has unique acyl chain anchors in Drosophila spermatocytes and is essential for meiotic cytokinesis. The head group of CPE is also important for spermatogenesis. We find that aberrant central spindle and contractile ring behavior but not mislocalization of phosphatidylinositol phosphates (PIPs) at the plasma membrane is responsible for the male meiotic cytokinesis defect in CPE-deficient animals. Further, we demonstrate the enrichment of CPE in multivesicular bodies marked by Rab7, which in turn localize to cleavage furrow. Volume electron microscopy analysis using correlative light and focused ion beam scanning electron microscopy shows that CPE-enriched Rab7 positive endosomes are juxtaposed on contractile ring material. Correlative light and transmission electron microscopy reveal Rab7 positive endosomes as a multivesicular body-like organelle that releases its intraluminal vesicles in the vicinity of ingressing furrows. Genetic ablation of Rab7 or Rab35 or expression of dominant negative Rab11 results in significant meiotic cytokinesis defects. Further, we show that Rab11 function is required for localization of CPE positive endosomes to the cleavage furrow. Our results imply that endosomal delivery of CPE to ingressing membranes is crucial for meiotic cytokinesis.

During cytokinesis, it is known that unique lipids with specific biophysical properties are localized to intercellular bridges connecting the two dividing cells, but how does this occur? This study shows that multivesicular bodies deliver sphingolipids to the ingressing membranes during male meiotic cytokinesis.

Evaluation of vertical and horizontal distribution of particulate matter near an urban roadway using an unmanned aerial vehicle

Measurement of traffic emissions has gained a lot of interest in recent times due to its contribution to urban pollution. This paper reports the outcome from an unmanned aerial vehicle (UAV) based measurement of PM concentration near an urban roadway at Kolkata, India. A total of 54 flights were carried out for simultaneous measurements of PM₁, PM_2.5 and PM₁₀ mass concentration and meteorological parameters in vertical as well as in horizontal direction. Results for the vertical flight up to 100 m showed that the PM₁, PM_2.5 and PM₁₀ concentrations at higher altitudes are less (mean; 24.6, 39.9 and 103.8 μg m^-3) compared to the respective ground level concentrations (mean; 26.3, 50.4 and 201.9 μg m^-3). For all the three particle sizes, the majority of the cases of higher PM concentration at higher altitudes happened during the evening flight. Low mixing height and low wind speed are suggested to be the reasons for the poor dispersion of pollutants in the evening. While there was a 7-10% fall of fine particles (PM₁ and PM_2.5) mass concentrations up to 90 m away from the road, no trend could be seen for PM₁₀. The random forest model to predict the UAV/Ground concentration ratio showed high accuracy (R² = 0.82-0.95) for all three particle sizes. This is an important finding from this study, which shows how UAV measurement data can be used to generate models that can predict the higher altitude concentrations from the ground based measurements.

The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2022 update

Galaxy is a mature, browser accessible workbench for scientific computing. It enables scientists to share, analyze and visualize their own data, with minimal technical impediments. A thriving global community continues to use, maintain and contribute to the project, with support from multiple national infrastructure providers that enable freely accessible analysis and training services. The Galaxy Training Network supports free, self-directed, virtual training with >230 integrated tutorials. Project engagement metrics have continued to grow over the last 2 years, including source code contributions, publications, software packages wrapped as tools, registered users and their daily analysis jobs, and new independent specialized servers. Key Galaxy technical developments include an improved user interface for launching large-scale analyses with many files, interactive tools for exploratory data analysis, and a complete suite of machine learning tools. Important scientific developments enabled by Galaxy include Vertebrate Genome Project (VGP) assembly workflows and global SARS-CoV-2 collaborations.

GenomicSuperSignature facilitates interpretation of RNA-seq experiments through robust, efficient comparison to public databases

Millions of transcriptomic profiles have been deposited in public archives, yet remain underused for the interpretation of new experiments. We present a method for interpreting new transcriptomic datasets through instant comparison to public datasets without high-performance computing requirements. We apply Principal Component Analysis on 536 studies comprising 44,890 human RNA sequencing profiles and aggregate sufficiently similar loading vectors to form Replicable Axes of Variation (RAV). RAVs are annotated with metadata of originating studies and by gene set enrichment analysis. Functionality to associate new datasets with RAVs, extract interpretable annotations, and provide intuitive visualization are implemented as the GenomicSuperSignature R/Bioconductor package. We demonstrate the efficient and coherent database search, robustness to batch effects and heterogeneous training data, and transfer learning capacity of our method using TCGA and rare diseases datasets. GenomicSuperSignature aids in analyzing new gene expression data in the context of existing databases using minimal computing resources.

Re-Expression of ERα and AR in Receptor Negative Endocrine Cancers via GSK3 Inhibition

DNA methylation, catalyzed by DNA methyltransferase (DNMT), is a well-characterized epigenetic modification in cancer cells. In particular, promoter hypermethylation of AR and ESR1 results in loss of expression on Androgen Receptor (AR) and Estrogen Receptor (ER), respectively, and is associated with a hormone refractory state. We now report that Glycogen Synthase Kinase 3 (GSK3) phosphorylates DNMT1 at S714, which is localized to a 62 amino acid region referred to as auto-inhibitory linker, which functions to occlude the DNA from the active site of DNMT1 to prevent the methylation of unmethylated DNA. Molecular Dynamics simulation indicates that phosphorylation at S714 resulted in conformational rearrangement of the autoinhibitory domain that inactivated its ability to block the methylation of unmethylated DNA and resulted in enhanced DNA binding. Treatment with a novel and more selective inhibitor of GSK3 resulted in decreased methylation of the promoter region of genes encoding the Androgen Receptor (AR) and Estrogen Receptor alpha (ERa) and re-expression of the AR and ERa in AR negative prostate cancer and ER negative breast cancer cells, respectively. As a result, concurrent treatment with the GSK3 inhibitor resulted in responsiveness of AR negative prostate cancer and ER negative breast cancer cells to inhibitors of the AR or ER, respectively, in in vitro and in vivo experimental models.

Expanding the Galaxy's reference data

Summary

Properly and effectively managing reference datasets is an important task for many bioinformatics analyses. Refgenie is a reference asset management system that allows users to easily organize, retrieve and share such datasets. Here, we describe the integration of refgenie into the Galaxy platform. Server administrators are able to configure Galaxy to make use of reference datasets made available on a refgenie instance. In addition, a Galaxy Data Manager tool has been developed to provide a graphical interface to refgenie's remote reference retrieval functionality. A large collection of reference datasets has also been made available using the CVMFS (CernVM File System) repository from GalaxyProject.org, with mirrors across the USA, Canada, Europe and Australia, enabling easy use outside of Galaxy.

Availability and implementation

The ability of Galaxy to use refgenie assets was added to the core Galaxy framework in version 22.01, which is available from https://github.com/galaxyproject/galaxy under the Academic Free License version 3.0. The refgenie Data Manager tool can be installed via the Galaxy ToolShed, with source code managed at https://github.com/BlankenbergLab/galaxy-tools-blankenberg/tree/main/data_managers/data_manager_refgenie_pull and released using an MIT license. Access to existing data is also available through CVMFS, with instructions at https://galaxyproject.org/admin/reference-data-repo/. No new data were generated or analyzed in support of this research.

Using Galaxy to Perform Large-Scale Interactive Data Analyses-An Update

Modern biology continues to become increasingly computational. Datasets are becoming progressively larger, more complex, and more abundant. The computational savviness necessary to analyze these data creates an ongoing obstacle for experimental biologists. Galaxy (galaxyproject.org) provides access to computational biology tools in a web-based interface. It also provides access to major public biological data repositories, allowing private data to be combined with public datasets. Galaxy is hosted on high-capacity servers worldwide and is accessible for free, with an option to be installed locally. This article demonstrates how to employ Galaxy to perform biologically relevant analyses on publicly available datasets. These protocols use both standard and custom tools, serving as a tutorial and jumping-off point for more intensive and/or more specific analyses using Galaxy. © 2021 Wiley Periodicals LLC. Basic Protocol 1: Finding human coding exons with highest SNP density Basic Protocol 2: Calling peaks for ChIP-seq data Basic Protocol 3: Compare datasets using genomic coordinates Basic Protocol 4: Working with multiple alignments Basic Protocol 5: Single cell RNA-seq.

SimText: A text mining framework for interactive analysis and visualization of similarities among biomedical entities

Summary

Literature exploration in PubMed on a large number of biomedical entities (e.g. genes, diseases or experiments) can be time-consuming and challenging, especially when assessing associations between entities. Here, we describe SimText, a user-friendly toolset that provides customizable and systematic workflows for the analysis of similarities among a set of entities based on text. SimText can be used for (i) text collection from PubMed and extraction of words with different text mining approaches, and (ii) interactive analysis and visualization of data using unsupervised learning techniques in an interactive app.

Availability and implementation

We developed SimText as an open-source R software and integrated it into Galaxy (https://usegalaxy.eu), an online data analysis platform with supporting self-learning training material available at https://training.galaxyproject.org. A command-line version of the toolset is available for download from GitHub (https://github.com/dlal-group/simtext) or as Docker image (https://hub.docker.com/r/dlalgroup/simtext/tags.).

Supplementary information

Supplementary data are available at Bioinformatics online.

A single-cell RNA-sequencing training and analysis suite using the Galaxy framework

Background

The vast ecosystem of single-cell RNA-sequencing tools has until recently been plagued by an excess of diverging analysis strategies, inconsistent file formats, and compatibility issues between different software suites. The uptake of 10x Genomics datasets has begun to calm this diversity, and the bioinformatics community leans once more towards the large computing requirements and the statistically driven methods needed to process and understand these ever-growing datasets.

Results

Here we outline several Galaxy workflows and learning resources for single-cell RNA-sequencing, with the aim of providing a comprehensive analysis environment paired with a thorough user learning experience that bridges the knowledge gap between the computational methods and the underlying cell biology. The Galaxy reproducible bioinformatics framework provides tools, workflows, and trainings that not only enable users to perform 1-click 10x preprocessing but also empower them to demultiplex raw sequencing from custom tagged and full-length sequencing protocols. The downstream analysis supports a range of high-quality interoperable suites separated into common stages of analysis: inspection, filtering, normalization, confounder removal, and clustering. The teaching resources cover concepts from computer science to cell biology. Access to all resources is provided at the singlecell.usegalaxy.eu portal.

Conclusions

The reproducible and training-oriented Galaxy framework provides a sustainable high-performance computing environment for users to run flexible analyses on both 10x and alternative platforms. The tutorials from the Galaxy Training Network along with the frequent training workshops hosted by the Galaxy community provide a means for users to learn, publish, and teach single-cell RNA-sequencing analysis.

No more business as usual: Agile and effective responses to emerging pathogen threats require open data and open analytics

The current state of much of the Wuhan pneumonia virus (severe acute respiratory syndrome coronavirus 2 [SARS-CoV-2]) research shows a regrettable lack of data sharing and considerable analytical obfuscation. This impedes global research cooperation, which is essential for tackling public health emergencies and requires unimpeded access to data, analysis tools, and computational infrastructure. Here, we show that community efforts in developing open analytical software tools over the past 10 years, combined with national investments into scientific computational infrastructure, can overcome these deficiencies and provide an accessible platform for tackling global health emergencies in an open and transparent manner. Specifically, we use all SARS-CoV-2 genomic data available in the public domain so far to (1) underscore the importance of access to raw data and (2) demonstrate that existing community efforts in curation and deployment of biomedical software can reliably support rapid, reproducible research during global health crises. All our analyses are fully documented at https://github.com/galaxyproject/SARS-CoV-2.

The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2020 update

Galaxy (https://galaxyproject.org) is a web-based computational workbench used by tens of thousands of scientists across the world to analyze large biomedical datasets. Since 2005, the Galaxy project has fostered a global community focused on achieving accessible, reproducible, and collaborative research. Together, this community develops the Galaxy software framework, integrates analysis tools and visualizations into the framework, runs public servers that make Galaxy available via a web browser, performs and publishes analyses using Galaxy, leads bioinformatics workshops that introduce and use Galaxy, and develops interactive training materials for Galaxy. Over the last two years, all aspects of the Galaxy project have grown: code contributions, tools integrated, users, and training materials. Key advances in Galaxy's user interface include enhancements for analyzing large dataset collections as well as interactive tools for exploratory data analysis. Extensions to Galaxy's framework include support for federated identity and access management and increased ability to distribute analysis jobs to remote resources. New community resources include large public servers in Europe and Australia, an increasing number of regional and local Galaxy communities, and substantial growth in the Galaxy Training Network.

A crowdsourced set of curated structural variants for the human genome

A high quality benchmark for small variants encompassing 88 to 90% of the reference genome has been developed for seven Genome in a Bottle (GIAB) reference samples. However a reliable benchmark for large indels and structural variants (SVs) is more challenging. In this study, we manually curated 1235 SVs, which can ultimately be used to evaluate SV callers or train machine learning models. We developed a crowdsourcing app-SVCurator-to help GIAB curators manually review large indels and SVs within the human genome, and report their genotype and size accuracy. SVCurator displays images from short, long, and linked read sequencing data from the GIAB Ashkenazi Jewish Trio son [NIST RM 8391/HG002]. We asked curators to assign labels describing SV type (deletion or insertion), size accuracy, and genotype for 1235 putative insertions and deletions sampled from different size bins between 20 and 892,149 bp. 'Expert' curators were 93% concordant with each other, and 37 of the 61 curators had at least 78% concordance with a set of 'expert' curators. The curators were least concordant for complex SVs and SVs that had inaccurate breakpoints or size predictions. After filtering events with low concordance among curators, we produced high confidence labels for 935 events. The SVCurator crowdsourced labels were 94.5% concordant with the heuristic-based draft benchmark SV callset from GIAB. We found that curators can successfully evaluate putative SVs when given evidence from multiple sequencing technologies.

Software engineering for scientific big data analysis

The increasing complexity of data and analysis methods has created an environment where scientists, who may not have formal training, are finding themselves playing the impromptu role of software engineer. While several resources are available for introducing scientists to the basics of programming, researchers have been left with little guidance on approaches needed to advance to the next level for the development of robust, large-scale data analysis tools that are amenable to integration into workflow management systems, tools, and frameworks. The integration into such workflow systems necessitates additional requirements on computational tools, such as adherence to standard conventions for robustness, data input, output, logging, and flow control. Here we provide a set of 10 guidelines to steer the creation of command-line computational tools that are usable, reliable, extensible, and in line with standards of modern coding practices.

Child Weight Gain Trajectories Linked To Oral Microbiota Composition

Gut and oral microbiota perturbations have been observed in obese adults and adolescents; less is known about their influence on weight gain in young children. Here we analyzed the gut and oral microbiota of 226 two-year-olds with 16S rRNA gene sequencing. Weight and length were measured at seven time points and used to identify children with rapid infant weight gain (a strong risk factor for childhood obesity), and to derive growth curves with innovative Functional Data Analysis (FDA) techniques. We showed that growth curves were associated negatively with diversity, and positively with the Firmicutes-to-Bacteroidetes ratio, of the oral microbiota. We also demonstrated an association between the gut microbiota and child growth, even after controlling for the effect of diet on the microbiota. Lastly, we identified several bacterial genera that were associated with child growth patterns. These results suggest that by the age of two, the oral microbiota of children with rapid infant weight gain may have already begun to establish patterns often seen in obese adults. They also suggest that the gut microbiota at age two, while strongly influenced by diet, does not harbor obesity signatures many researchers identified in later life stages.

The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update

Galaxy (homepage: https://galaxyproject.org, main public server: https://usegalaxy.org) is a web-based scientific analysis platform used by tens of thousands of scientists across the world to analyze large biomedical datasets such as those found in genomics, proteomics, metabolomics and imaging. Started in 2005, Galaxy continues to focus on three key challenges of data-driven biomedical science: making analyses accessible to all researchers, ensuring analyses are completely reproducible, and making it simple to communicate analyses so that they can be reused and extended. During the last two years, the Galaxy team and the open-source community around Galaxy have made substantial improvements to Galaxy's core framework, user interface, tools, and training materials. Framework and user interface improvements now enable Galaxy to be used for analyzing tens of thousands of datasets, and >5500 tools are now available from the Galaxy ToolShed. The Galaxy community has led an effort to create numerous high-quality tutorials focused on common types of genomic analyses. The Galaxy developer and user communities continue to grow and be integral to Galaxy's development. The number of Galaxy public servers, developers contributing to the Galaxy framework and its tools, and users of the main Galaxy server have all increased substantially.

Recommendations for the packaging and containerizing of bioinformatics software

Software Containers are changing the way scientists and researchers develop, deploy and exchange scientific software. They allow labs of all sizes to easily install bioinformatics software, maintain multiple versions of the same software and combine tools into powerful analysis pipelines. However, containers and software packages should be produced under certain rules and standards in order to be reusable, compatible and easy to integrate into pipelines and analysis workflows. Here, we presented a set of recommendations developed by the BioContainers Community to produce standardized bioinformatics packages and containers. These recommendations provide practical guidelines to make bioinformatics software more discoverable, reusable and transparent.  They are aimed to guide developers, organisations, journals and funders to increase the quality and sustainability of research software.

Recommendations for the packaging and containerizing of bioinformatics software

Software Containers are changing the way scientists and researchers develop, deploy and exchange scientific software. They allow labs of all sizes to easily install bioinformatics software, maintain multiple versions of the same software and combine tools into powerful analysis pipelines. However, containers and software packages should be produced under certain rules and standards in order to be reusable, compatible and easy to integrate into pipelines and analysis workflows. Here, we presented a set of recommendations developed by the BioContainers Community to produce standardized bioinformatics packages and containers. These recommendations provide practical guidelines to make bioinformatics software more discoverable, reusable and transparent.  They are aimed to guide developers, organisations, journals and funders to increase the quality and sustainability of research software.

Biology Needs Evolutionary Software Tools: Let's Build Them Right

Research in population genetics and evolutionary biology has always provided a computational backbone for life sciences as a whole. Today evolutionary and population biology reasoning are essential for interpretation of large complex datasets that are characteristic of all domains of today's life sciences ranging from cancer biology to microbial ecology. This situation makes algorithms and software tools developed by our community more important than ever before. This means that we, developers of software tool for molecular evolutionary analyses, now have a shared responsibility to make these tools accessible using modern technological developments as well as provide adequate documentation and training.

The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2016 update

High-throughput data production technologies, particularly ‘next-generation’ DNA sequencing, have ushered in widespread and disruptive changes to biomedical research. Making sense of the large datasets produced by these technologies requires sophisticated statistical and computational methods, as well as substantial computational power. This has led to an acute crisis in life sciences, as researchers without informatics training attempt to perform computation-dependent analyses. Since 2005, the Galaxy project has worked to address this problem by providing a framework that makes advanced computational tools usable by non experts. Galaxy seeks to make data-intensive research more accessible, transparent and reproducible by providing a Web-based environment in which users can perform computational analyses and have all of the details automatically tracked for later inspection, publication, or reuse. In this report we highlight recently added features enabling biomedical analyses on a large scale.

Online resources for genomic analysis using high-throughput sequencing

The availability of high-throughput sequencing has created enormous possibilities for scientific discovery. However, the massive amount of data being generated has resulted in a severe informatics bottleneck. A large number of tools exist for analyzing next-generation sequencing (NGS) data, yet often there remains a disconnect between these research tools and the ability of many researchers to use them. As a consequence, several online resources and communities have been developed to assist researchers with both the management and the analysis of sequencing data sets. Here we describe the use and applications of common file formats for coding and storing genomic data, consider several web-accessible open-source resources for the visualization and analysis of NGS data, and provide examples of typical analyses with links to further detailed exercises.

Controlling for contamination in re-sequencing studies with a reproducible web-based phylogenetic approach

Polymorphism discovery is a routine application of next-generation sequencing technology where multiple samples are sent to a service provider for library preparation, subsequent sequencing, and bioinformatic analyses. The decreasing cost and advances in multiplexing approaches have made it possible to analyze hundreds of samples at a reasonable cost. However, because of the manual steps involved in the initial processing of samples and handling of sequencing equipment, cross-contamination remains a significant challenge. It is especially problematic in cases where polymorphism frequencies do not adhere to diploid expectation, for example, heterogeneous tumor samples, organellar genomes, as well as during bacterial and viral sequencing. In these instances, low levels of contamination may be readily mistaken for polymorphisms, leading to false results. Here we describe practical steps designed to reliably detect contamination and uncover its origin, and also provide new, Galaxy-based, readily accessible computational tools and workflows for quality control. All results described in this report can be reproduced interactively on the web as described at http://usegalaxy.org/contamination.

Wrangling Galaxy's reference data

Summary: The Galaxy platform has developed into a fully featured collaborative workbench, with goals of inherently capturing provenance to enable reproducible data analysis, and of making it straightforward to run one’s own server. However, many Galaxy platform tools rely on the presence of reference data, such as alignment indexes, to function efficiently. Until now, the building of this cache of data for Galaxy has been an error-prone manual process lacking reproducibility and provenance.

The Galaxy Data Manager framework is an enhancement that changes the management of Galaxy’s built-in data cache from a manual procedure to an automated graphical user interface (GUI) driven process, which contains the same openness, reproducibility and provenance that is afforded to Galaxy’s analysis tools. Data Manager tools allow the Galaxy administrator to download, create and install additional datasets for any type of reference data in real time.

Availability and implementation: The Galaxy Data Manager framework is implemented in Python and has been integrated as part of the core Galaxy platform. Individual Data Manager tools can be defined locally or installed from a ToolShed, allowing the Galaxy community to define additional Data Manager tools as needed, with full versioning and dependency support.

Contact:dan@bx.psu.edu. or anton@bx.psu.edu

Supplementary information:Supplementary data is available at Bioinformatics online.

Dissemination of scientific software with Galaxy ToolShed

The proliferation of web-based integrative analysis frameworks has enabled users to perform complex analyses directly through the web. Unfortunately, it also revoked the freedom to easily select the most appropriate tools. To address this, we have developed Galaxy ToolShed.

Analysis of next-generation sequencing data using Galaxy

The extraordinary throughput of next-generation sequencing (NGS) technology is outpacing our ability to analyze and interpret the data. This chapter will focus on practical informatics methods, strategies, and software tools for transforming NGS data into usable information through the use of a web-based platform, Galaxy. The Galaxy interface is explored through several different types of example analyses. Instructions for running one's own Galaxy server on local hardware or on cloud computing resources are provided. Installing new tools into a personal Galaxy instance is also demonstrated.

Using Galaxy to perform large-scale interactive data analyses

Innovations in biomedical research technologies continue to provide experimental biologists with novel and increasingly large genomic and high-throughput data resources to be analyzed. As creating and obtaining data has become easier, the key decision faced by many researchers is a practical one: where and how should an analysis be performed? Datasets are large and analysis tool set-up and use is riddled with complexities outside of the scope of core research activities. The authors believe that Galaxy provides a powerful solution that simplifies data acquisition and analysis in an intuitive Web application, granting all researchers access to key informatics tools previously only available to computational specialists working in Unix-based environments. We will demonstrate through a series of biomedically relevant protocols how Galaxy specifically brings together (1) data retrieval from public and private sources, for example, UCSC's Eukaryote and Microbial Genome Browsers, (2) custom tools (wrapped Unix functions, format standardization/conversions, interval operations), and 3rd-party analysis tools.

An encyclopedia of mouse DNA elements (Mouse ENCODE)

To complement the human Encyclopedia of DNA Elements (ENCODE) project and to enable a broad range of mouse genomics efforts, the Mouse ENCODE Consortium is applying the same experimental pipelines developed for human ENCODE to annotate the mouse genome.

Making whole genome multiple alignments usable for biologists

SUMMARY

Here we describe a set of tools implemented within the Galaxy platform designed to make analysis of multiple genome alignments truly accessible for biologists. These tools are available through both a web-based graphical user interface and a command-line interface.

AVAILABILITY AND IMPLEMENTATION

This open-source toolset was implemented in Python and has been integrated into the online data analysis platform Galaxy (public web access: http://usegalaxy.org; download: http://getgalaxy.org). Additional help is available as a live supplement from http://usegalaxy.org/u/dan/p/maf.

CONTACT

james.taylor@emory.edu; anton@bx.psu.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Integrating diverse databases into an unified analysis framework: a Galaxy approach

Recent technological advances have lead to the ability to generate large amounts of data for model and non-model organisms. Whereas, in the past, there have been a relatively small number of central repositories that serve genomic data, an increasing number of distinct specialized data repositories and resources have been established. Here, we describe a generic approach that provides for the integration of a diverse spectrum of data resources into a unified analysis framework, Galaxy (http://usegalaxy.org). This approach allows the simplified coupling of external data resources with the data analysis tools available to Galaxy users, while leveraging the native data mining facilities of the external data resources. DATABASE URL: http://usegalaxy.org.

Manipulation of FASTQ data with Galaxy

Summary: Here, we describe a tool suite that functions on all of the commonly known FASTQ format variants and provides a pipeline for manipulating next generation sequencing data taken from a sequencing machine all the way through the quality filtering steps.

Availability and Implementation: This open-source toolset was implemented in Python and has been integrated into the online data analysis platform Galaxy (public web access: http://usegalaxy.org; download: http://getgalaxy.org). Two short movies that highlight the functionality of tools described in this manuscript as well as results from testing components of this tool suite against a set of previously published files are available at http://usegalaxy.org/u/dan/p/fastq

Contact:james.taylor@emory.edu; anton@bx.psu.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

Galaxy: a web-based genome analysis tool for experimentalists

High-throughput data production has revolutionized molecular biology. However, massive increases in data generation capacity require analysis approaches that are more sophisticated, and often very computationally intensive. Thus, making sense of high-throughput data requires informatics support. Galaxy (http://galaxyproject.org) is a software system that provides this support through a framework that gives experimentalists simple interfaces to powerful tools, while automatically managing the computational details. Galaxy is distributed both as a publicly available Web service, which provides tools for the analysis of genomic, comparative genomic, and functional genomic data, or a downloadable package that can be deployed in individual laboratories. Either way, it allows experimentalists without informatics or programming expertise to perform complex large-scale analysis with just a Web browser.

28-way vertebrate alignment and conservation track in the UCSC Genome Browser

This article describes a set of alignments of 28 vertebrate genome sequences that is provided by the UCSC Genome Browser. The alignments can be viewed on the Human Genome Browser (March 2006 assembly) at http://genome.ucsc.edu, downloaded in bulk by anonymous FTP from http://hgdownload.cse.ucsc.edu/goldenPath/hg18/multiz28way, or analyzed with the Galaxy server at http://g2.bx.psu.edu. This article illustrates the power of this resource for exploring vertebrate and mammalian evolution, using three examples. First, we present several vignettes involving insertions and deletions within protein-coding regions, including a look at some human-specific indels. Then we study the extent to which start codons and stop codons in the human sequence are conserved in other species, showing that start codons are in general more poorly conserved than stop codons. Finally, an investigation of the phylogenetic depth of conservation for several classes of functional elements in the human genome reveals striking differences in the rates and modes of decay in alignability. Each functional class has a distinctive period of stringent constraint, followed by decays that allow (for the case of regulatory regions) or reject (for coding regions and ultraconserved elements) insertions and deletions.

A framework for collaborative analysis of ENCODE data: making large-scale analyses biologist-friendly

The standardization and sharing of data and tools are the biggest challenges of large collaborative projects such as the Encyclopedia of DNA Elements (ENCODE). Here we describe a compact Web application, Galaxy2(ENCODE), that effectively addresses these issues. It provides an intuitive interface for the deposition and access of data, and features a vast number of analysis tools including operations on genomic intervals, utilities for manipulation of multiple sequence alignments, and molecular evolution algorithms. By providing a direct link between data and analysis tools, Galaxy2(ENCODE) allows addressing biological questions that are beyond the reach of existing software. We use Galaxy2(ENCODE) to show that the ENCODE regions contain >2000 unannotated transcripts under strong purifying selection that are likely functional. We also show that the ENCODE regions are representative of the entire genome by estimating the rate of nucleotide substitution and comparing it to published data. Although each of these analyses is complex, none takes more than 15 min from beginning to end. Finally, we demonstrate how new tools can be added to Galaxy2(ENCODE) with almost no effort. Every section of the manuscript is supplemented with QuickTime screencasts. Galaxy2(ENCODE) and the screencasts can be accessed at http://g2.bx.psu.edu.

Galaxy: a platform for interactive large-scale genome analysis

Accessing and analyzing the exponentially expanding genomic sequence and functional data pose a challenge for biomedical researchers. Here we describe an interactive system, Galaxy, that combines the power of existing genome annotation databases with a simple Web portal to enable users to search remote resources, combine data from independent queries, and visualize the results. The heart of Galaxy is a flexible history system that stores the queries from each user; performs operations such as intersections, unions, and subtractions; and links to other computational tools. Galaxy can be accessed at http://g2.bx.psu.edu.