Benchmark Software and Data for Evaluating CRISPR-Cas9 Experimental Pipelines Through the Assessment of a Calibration Screen.

Genome-wide genetic screens using CRISPR-guide RNA libraries are widely performed in mammalian cells to functionally characterize individual genes and for the discovery of new anticancer therapeutic targets. As the effectiveness of such powerful and precise tools for cancer pharmacogenomics is emerging, tools and methods for their quality assessment are becoming increasingly necessary. Here, we provide an R package and a high-quality reference data set for the assessment of novel experimental pipelines through which a single calibration experiment has been executed: a screen of the HT-29 human colorectal cancer cell line with a commercially available genome-wide library of single-guide RNAs. This package and data allow experimental researchers to benchmark their screens and produce a quality-control report, encompassing several quality and validation metrics. The R code used for processing the reference data set, for its quality assessment, as well as to evaluate the quality of a user-provided screen, and to reproduce the figures presented in this article is available at https://github.com/DepMap-Analytics/HT29benchmark. The reference data is publicly available on FigShare.

ALFI: Cell cycle phenotype annotations of label-free time-lapse imaging data from cultured human cells.

Detecting and tracking multiple moving objects in a video is a challenging task. For living cells, the task becomes even more arduous as cells change their morphology over time, can partially overlap, and mitosis leads to new cells. Differently from fluorescence microscopy, label-free techniques can be easily applied to almost all cell lines, reducing sample preparation complexity and phototoxicity. In this study, we present ALFI, a dataset of images and annotations for label-free microscopy, made publicly available to the scientific community, that notably extends the current panorama of expertly labeled data for detection and tracking of cultured living nontransformed and cancer human cells. It consists of 29 time-lapse image sequences from HeLa, U2OS, and hTERT RPE-1 cells under different experimental conditions, acquired by differential interference contrast microscopy, for a total of 237.9 hours. It contains various annotations (pixel-wise segmentation masks, object-wise bounding boxes, tracking information). The dataset is useful for testing and comparing methods for identifying interphase and mitotic events and reconstructing their lineage, and for discriminating different cellular phenotypes.

Results and lessons learned from the sbv IMPROVER metagenomics diagnostics for inflammatory bowel disease challenge.

A growing body of evidence links gut microbiota changes with inflammatory bowel disease (IBD), raising the potential benefit of exploiting metagenomics data for non-invasive IBD diagnostics. The sbv IMPROVER metagenomics diagnosis for inflammatory bowel disease challenge investigated computational metagenomics methods for discriminating IBD and nonIBD subjects. Participants in this challenge were given independent training and test metagenomics data from IBD and nonIBD subjects, which could be wither either raw read data (sub-challenge 1, SC1) or processed Taxonomy- and Function-based profiles (sub-challenge 2, SC2). A total of 81 anonymized submissions were received between September 2019 and March 2020. Most participants' predictions performed better than random predictions in classifying IBD versus nonIBD, Ulcerative Colitis (UC) versus nonIBD, and Crohn's Disease (CD) versus nonIBD. However, discrimination between UC and CD remains challenging, with the classification quality similar to the set of random predictions. We analyzed the class prediction accuracy, the metagenomics features by the teams, and computational methods used. These results will be openly shared with the scientific community to help advance IBD research and illustrate the application of a range of computational methodologies for effective metagenomic classification.

Early anteroposterior regionalisation of human neural crest is shaped by a pro-mesodermal factor.

The neural crest (NC) is an important multipotent embryonic cell population and its impaired specification leads to various developmental defects, often in an anteroposterior (A-P) axial level-specific manner. The mechanisms underlying the correct A-P regionalisation of human NC cells remain elusive. Recent studies have indicated that trunk NC cells, the presumed precursors of childhood tumour neuroblastoma, are derived from neuromesodermal-potent progenitors of the postcranial body. Here we employ human embryonic stem cell differentiation to define how neuromesodermal progenitor (NMP)-derived NC cells acquire a posterior axial identity. We show that TBXT, a pro-mesodermal transcription factor, mediates early posterior NC/spinal cord regionalisation together with WNT signalling effectors. This occurs by TBXT-driven chromatin remodelling via its binding in key enhancers within HOX gene clusters and other posterior regulator-associated loci. This initial posteriorisation event is succeeded by a second phase of trunk HOX gene control that marks the differentiation of NMPs toward their TBXT-negative NC/spinal cord derivatives and relies predominantly on FGF signalling. Our work reveals a previously unknown role of TBXT in influencing posterior NC fate and points to the existence of temporally discrete, cell type-dependent modes of posterior axial identity control.

Netpro2vec: A Graph Embedding Framework for Biomedical Applications.

The ever-increasing importance of structured data in different applications, especially in the biomedical field, has driven the need for reducing its complexity through projections into a more manageable space. The latest methods for learning features on graphs focus mainly on the neighborhood of nodes and edges. Methods capable of providing a representation that looks beyond the single node neighborhood are kernel graphs. However, they produce handcrafted features unaccustomed with a generalized model. To reduce this gap, in this work we propose a neural embedding framework, based on probability distribution representations of graphs, named Netpro2vec. The goal is to look at basic node descriptions other than the degree, such as those induced by the Transition Matrix and Node Distance Distribution. Netpro2vec provides embeddings completely independent from the task and nature of the data. The framework is evaluated on synthetic and various real biomedical network datasets through a comprehensive experimental classification phase and is compared to well-known competitors.

SRSF1-dependent inhibition of C9ORF72-repeat RNA nuclear export: genome-wide mechanisms for neuroprotection in amyotrophic lateral sclerosis.

BACKGROUND: Loss of motor neurons in amyotrophic lateral sclerosis (ALS) leads to progressive paralysis and death. Dysregulation of thousands of RNA molecules with roles in multiple cellular pathways hinders the identification of ALS-causing alterations over downstream changes secondary to the neurodegenerative process. How many and which of these pathological gene expression changes require therapeutic normalisation remains a fundamental question. METHODS: Here, we investigated genome-wide RNA changes in C9ORF72-ALS patient-derived neurons and Drosophila, as well as upon neuroprotection taking advantage of our gene therapy approach which specifically inhibits the SRSF1-dependent nuclear export of pathological C9ORF72-repeat transcripts. This is a critical study to evaluate (i) the overall safety and efficacy of the partial depletion of SRSF1, a member of a protein family involved itself in gene expression, and (ii) a unique opportunity to identify neuroprotective RNA changes. RESULTS: Our study shows that manipulation of 362 transcripts out of 2257 pathological changes, in addition to inhibiting the nuclear export of repeat transcripts, is sufficient to confer neuroprotection in C9ORF72-ALS patient-derived neurons. In particular, expression of 90 disease-altered transcripts is fully reverted upon neuroprotection leading to the characterisation of a human C9ORF72-ALS disease-modifying gene expression signature. These findings were further investigated in vivo in diseased and neuroprotected Drosophila transcriptomes, highlighting a list of 21 neuroprotective changes conserved with 16 human orthologues in patient-derived neurons. We also functionally validated the high neuroprotective potential of one of these disease-modifying transcripts, demonstrating that inhibition of ALS-upregulated human KCNN1-3 (Drosophila SK) voltage-gated potassium channel orthologs mitigates degeneration of human motor neurons and Drosophila motor deficits. CONCLUSIONS: Strikingly, the partial depletion of SRSF1 leads to expression changes in only a small proportion of disease-altered transcripts, indicating that not all RNA alterations need normalization and that the gene therapeutic approach is safe in the above preclinical models as it does not disrupt globally gene expression. The efficacy of this intervention is also validated at genome-wide level with transcripts modulated in the vast majority of biological processes affected in C9ORF72-ALS. Finally, the identification of a characteristic signature with key RNA changes modified in both the disease state and upon neuroprotection also provides potential new therapeutic targets and biomarkers.

Defining the signalling determinants of a posterior ventral spinal cord identity in human neuromesodermal progenitor derivatives.

The anteroposterior axial identity of motor neurons (MNs) determines their functionality and vulnerability to neurodegeneration. Thus, it is a crucial parameter in the design of strategies aiming to produce MNs from human pluripotent stem cells (hPSCs) for regenerative medicine/disease modelling applications. However, the in vitro generation of posterior MNs corresponding to the thoracic/lumbosacral spinal cord has been challenging. Although the induction of cells resembling neuromesodermal progenitors (NMPs), the bona fide precursors of the spinal cord, offers a promising solution, the progressive specification of posterior MNs from these cells is not well defined. Here, we determine the signals guiding the transition of human NMP-like cells toward thoracic ventral spinal cord neurectoderm. We show that combined WNT-FGF activities drive a posterior dorsal pre-/early neural state, whereas suppression of TGFß-BMP signalling pathways promotes a ventral identity and neural commitment. Based on these results, we define an optimised protocol for the generation of thoracic MNs that can efficiently integrate within the neural tube of chick embryos. We expect that our findings will facilitate the comparison of hPSC-derived spinal cord cells of distinct axial identities.

Pooled testing with replication as a mass testing strategy for the COVID-19 pandemics.

During the COVID-19 pandemic it is essential to test as many people as possible, in order to detect early outbreaks of the infection. Present testing solutions are based on the extraction of RNA from patients using oropharyngeal and nasopharyngeal swabs, and then testing with real-time PCR for the presence of specific RNA filaments identifying the virus. This approach is limited by the availability of reactants, trained technicians and laboratories. One of the ways to speed up the testing procedures is a group testing, where the swabs of multiple patients are grouped together and tested. In this paper we propose to use the group testing technique in conjunction with an advanced replication scheme in which each patient is allocated in two or more groups to reduce the total numbers of tests and to allow testing of even larger numbers of people. Under mild assumptions, a 13 ×  average reduction of tests can be achieved compared to individual testing without delay in time.

Learning from Metabolic Networks: Current Trends and Future Directions for Precision Medicine.

BACKGROUND: Systems biology and network modeling represent, nowadays, the hallmark approaches for the development of predictive and targeted-treatment based precision medicine. The study of health and disease as properties of the human body system allows the understanding of the genotype-phenotype relationship through the definition of molecular interactions and dependencies. In this scenario, metabolism plays a central role as its interactions are well characterized and it is considered an important indicator of the genotype- phenotype associations. In metabolic systems biology, the genome-scale metabolic models are the primary scaffolds to integrate multi-omics data as well as cell-, tissue-, condition- specific information. Modeling the metabolism has both investigative and predictive values. Several methods have been proposed to model systems, which involve steady-state or kinetic approaches, and to extract knowledge through machine and deep learning. METHODS: This review collects, analyzes, and compares the suitable data and computational approaches for the exploration of metabolic networks as tools for the development of precision medicine. To this extent, we organized it into three main sections: "Data and Databases", "Methods and Tools", and "Metabolic Networks for medicine". In the first one, we have collected the most used data and relative databases to build and annotate metabolic models. In the second section, we have reported the state-of-the-art methods and relative tools to reconstruct, simulate, and interpret metabolic systems. Finally, we have reported the most recent and innovative studies that exploited metabolic networks to study several pathological conditions, not only those directly related to metabolism. CONCLUSION: We think that this review can be a guide to researchers of different disciplines, from computer science to biology and medicine, in exploring the power, challenges and future promises of the metabolism as predictor and target of the so-called P4 medicine (predictive, preventive, personalized and participatory).

Correction to: Clustering analysis of tumor metabolic networks.

An amendment to this paper has been published and can be accessed via the original article.

Clustering analysis of tumor metabolic networks.

BACKGROUND: Biological networks are representative of the diverse molecular interactions that occur within cells. Some of the commonly studied biological networks are modeled through protein-protein interactions, gene regulatory, and metabolic pathways. Among these, metabolic networks are probably the most studied, as they directly influence all physiological processes. Exploration of biochemical pathways using multigraph representation is important in understanding complex regulatory mechanisms. Feature extraction and clustering of these networks enable grouping of samples obtained from different biological specimens. Clustering techniques separate networks depending on their mutual similarity. RESULTS: We present a clustering analysis on tissue-specific metabolic networks for single samples from three primary tumor sites: breast, lung, and kidney cancer. The metabolic networks were obtained by integrating genome scale metabolic models with gene expression data. We performed network simplification to reduce the computational time needed for the computation of network distances. We empirically proved that networks clustering can characterize groups of patients in multiple conditions. CONCLUSIONS: We provide a computational methodology to explore and characterize the metabolic landscape of tumors, thus providing a general methodology to integrate analytic metabolic models with gene expression data. This method represents a first attempt in clustering large scale metabolic networks. Moreover, this approach gives the possibility to get valuable information on what are the effects of different conditions on the overall metabolism.

Human axial progenitors generate trunk neural crest cells in vitro.

The neural crest (NC) is a multipotent embryonic cell population that generates distinct cell types in an axial position-dependent manner. The production of NC cells from human pluripotent stem cells (hPSCs) is a valuable approach to study human NC biology. However, the origin of human trunk NC remains undefined and current in vitro differentiation strategies induce only a modest yield of trunk NC cells. Here we show that hPSC-derived axial progenitors, the posteriorly-located drivers of embryonic axis elongation, give rise to trunk NC cells and their derivatives. Moreover, we define the molecular signatures associated with the emergence of human NC cells of distinct axial identities in vitro. Collectively, our findings indicate that there are two routes toward a human post-cranial NC state: the birth of cardiac and vagal NC is facilitated by retinoic acid-induced posteriorisation of an anterior precursor whereas trunk NC arises within a pool of posterior axial progenitors.

Glycosphingolipid metabolic reprogramming drives neural differentiation.

Neural development is accomplished by differentiation events leading to metabolic reprogramming. Glycosphingolipid metabolism is reprogrammed during neural development with a switch from globo- to ganglio-series glycosphingolipid production. Failure to execute this glycosphingolipid switch leads to neurodevelopmental disorders in humans, indicating that glycosphingolipids are key players in this process. Nevertheless, both the molecular mechanisms that control the glycosphingolipid switch and its function in neurodevelopment are poorly understood. Here, we describe a self-contained circuit that controls glycosphingolipid reprogramming and neural differentiation. We find that globo-series glycosphingolipids repress the epigenetic regulator of neuronal gene expression AUTS2. AUTS2 in turn binds and activates the promoter of the first and rate-limiting ganglioside-producing enzyme GM3 synthase, thus fostering the synthesis of gangliosides. By this mechanism, the globo-AUTS2 axis controls glycosphingolipid reprogramming and neural gene expression during neural differentiation, which involves this circuit in neurodevelopment and its defects in neuropathology.

A web-oriented software for the optimization of pooled experiments in NGS for detection of rare mutations.

BACKGROUND: The cost per patient of next generation sequencing for detection of rare mutations may be significantly reduced using pooled experiments. Recently, some techniques have been proposed for the planning of pooled experiments and for the optimal allocation of patients into pools. However, the lack of a user friendly resource for planning the design of pooled experiments forces the scientists to do frequent, complex and long computations. RESULTS: OPENDoRM is a powerful collection of novel mathematical algorithms usable via an intuitive graphical user interface. It enables researchers to speed up the planning of their routine experiments, as well as, to support scientists without specific bioinformatics expertises. Users can automatically carry out analysis in terms of costs associated with the optimal allocation of patients in pools. They are also able to choose between three distinct pooling mathematical methods, each of which also suggests the optimal configuration for the submitted experiment. Importantly, in order to keep track of the performed experiments, users can save and export the results of their experiments in standard tabular and charts contents. CONCLUSION: OPENDoRM is a freely available web-oriented application for the planning of pooled NGS experiments, available at: http://www-labgtp.na.icar.cnr.it/OPENDoRM. Its easy and intuitive graphical user interface enables researchers to plan theirs experiments using novel algorithms, and to interactively visualize the results.

Cold adaptation shapes the robustness of metabolic networks in Drosophila melanogaster.

When ectotherms are exposed to low temperatures, they enter a cold-induced coma (chill coma) that prevents resource acquisition, mating, oviposition, and escape from predation. There is substantial variation in time taken to recover from chill coma both within and among species, and this variation is correlated with habitat temperatures such that insects from cold environments recover more quickly. This suggests an adaptive response, but the mechanisms underlying variation in recovery times are unknown, making it difficult to decisively test adaptive hypotheses. We use replicated lines of Drosophila melanogaster selected in the laboratory for fast (hardy) or slow (susceptible) chill-coma recovery times to investigate modifications to metabolic profiles associated with cold adaptation. We measured metabolite concentrations of flies before, during, and after cold exposure using nuclear magnetic resonance (NMR) spectroscopy to test the hypotheses that hardy flies maintain metabolic homeostasis better during cold exposure and recovery, and that their metabolic networks are more robust to cold-induced perturbations. The metabolites of cold-hardy flies were less cold responsive and their metabolic networks during cold exposure were more robust, supporting our hypotheses. Metabolites involved in membrane lipid synthesis, tryptophan metabolism, oxidative stress, energy balance, and proline metabolism were altered by selection on cold tolerance. We discuss the potential significance of these alterations.

Data mining in psychiatric research.

Mathematical sciences and computational methods have found new applications in fields like medicine over the last few decades. Modern data acquisition and data analysis protocols have been of great assistance to medical researchers and clinical scientists. Especially in psychiatry, technology and science have made new computational methods available to assist the development of predictive modeling and to identify diseases more accurately. Data mining (or knowledge discovery) aims to extract information from large datasets and solve challenging tasks, like patient assessment, early mental disease diagnosis, and drug efficacy assessment. Accurate and fast data analysis methods are very important, especially when dealing with severe psychiatric diseases like schizophrenia. In this paper, we focus on computational methods related to data analysis and more specifically to data mining. Then, we discuss some related research in the field of psychiatry.

Classification of cancer cell death with spectral dimensionality reduction and generalized eigenvalues.

OBJECTIVE: Accurate cell death discrimination is a time consuming and expensive process that can only be performed in biological laboratories. Nevertheless, it is very useful and arises in many biological and medical applications. METHODS AND MATERIAL: Raman spectra are collected for 84 samples of A549 cell line (human lung cancer epithelia cells) that has been exposed to toxins to simulate the necrotic and apoptotic death. The proposed data mining approach for the multiclass cell death discrimination problem uses a multiclass regularized generalized eigenvalue algorithm for classification (multiReGEC), together with a dimensionality reduction algorithm based on spectral clustering. RESULTS: The proposed algorithmic scheme can classify A549 lung cancer cells from three different classes (apoptotic death, necrotic death and control cells) with 97.78%± 0.047 accuracy versus 92.22 ± 0.095 without the proposed feature selection preprocessing. The spectrum areas depicted by the algorithm corresponds to the 〉C O bond from the lipids and the lipid bilayer. This chemical structure undergoes different change of state based on cell death type. Further evidence of the validity of the technique is obtained through the successful classification of 7 cell spectra that undergo hyperthermic treatment. CONCLUSIONS: In this study we propose a fast and automated way of processing Raman spectra for cell death discrimination, using a feature selection algorithm that not only enhances the classification accuracy, but also gives more insight in the undergoing cell death process.

Prediction of the responsiveness to pharmacological chaperones: lysosomal human alpha-galactosidase, a case of study.

BACKGROUND: The pharmacological chaperones therapy is a promising approach to cure genetic diseases. It relies on substrate competitors used at sub-inhibitory concentration which can be administered orally, reach difficult tissues and have low cost. Clinical trials are currently carried out for Fabry disease, a lysosomal storage disorder caused by inherited genetic mutations of alpha-galactosidase. Regrettably, not all genotypes respond to these drugs. RESULTS: We collected the experimental data available in literature on the enzymatic activity of ninety-six missense mutants of lysosomal alpha-galactosidase measured in the presence of pharmacological chaperones. We associated with each mutation seven features derived from the analysis of 3D-structure of the enzyme, two features associated with their thermo-dynamic stability and four features derived from sequence alone. Structural and thermodynamic analysis explains why some mutants of human lysosomal alpha-galactosidase cannot be rescued by pharmacological chaperones: approximately forty per cent of the non responsive cases examined can be correctly associated with a negative prognostic feature. They include mutations occurring in the active site pocket, mutations preventing disulphide bridge formation and severely destabilising mutations. Despite this finding, prediction of mutations responsive to pharmacological chaperones cannot be achieved with high accuracy relying on combinations of structure- and thermodynamic-derived features even with the aid of classical and state of the art statistical learning methods.We developed a procedure to predict responsive mutations with an accuracy as high as 87%: the method scores the mutations by using a suitable position-specific substitution matrix. Our approach is of general applicability since it does not require the knowledge of 3D-structure but relies only on the sequence. CONCLUSIONS: Responsiveness to pharmacological chaperones depends on the structural/functional features of the disease-associated protein, whose complex interplay is best reflected on sequence conservation by evolutionary pressure. We propose a predictive method which can be applied to screen novel mutations of alpha galactosidase. The same approach can be extended on a genomic scale to find candidates for therapy with pharmacological chaperones among proteins with unknown tertiary structures.