There is a lack of robust generalizable predictive biomarkers of response to immune checkpoint blockade in multiple types of cancer. We develop hDirect-MAP, an algorithm that maps T cells into a shared high-dimensional (HD) expression space of diverse T cell functional signatures in which cells group by the common T cell phenotypes rather than dimensional reduced features or a distorted view of these features. Using projection-free single-cell modeling, hDirect-MAP first removed a large group of cells that did not contribute to response and then clearly distinguished T cells into response-specific subpopulations that were defined by critical T cell functional markers of strong differential expression patterns. We found that these grouped cells cannot be distinguished by dimensional-reduction algorithms but are blended by diluted expression patterns. Moreover, these identified response-specific T cell subpopulations enabled a generalizable prediction by their HD metrics. Tested using five single-cell RNA-seq or mass cytometry datasets from basal cell carcinoma, squamous cell carcinoma and melanoma, hDirect-MAP demonstrated common response-specific T cell phenotypes that defined a generalizable and accurate predictive biomarker.

SUMMARY: Mutation is the key for a variant of concern (VOC) to overcome selective pressures, but this process is still unclear. Understanding the association of the mutational process with VOCs is an unmet need. Motivation: Here, we developed VOC-alarm, a method to predict VOCs and their caused COVID surges, using mutations of about 5.7 million SARS-CoV-2 complete sequences. We found that VOCs rely on lineage-level entropy value of mutation numbers to compete with other variants, suggestive of the importance of population-level mutations in the virus evolution. Thus, we hypothesized that VOCs are a result of a mutational process across the globe. Results: Analyzing the mutations from January 2020 to December 2021, we simulated the mutational process by estimating the pace of evolution, and thus divided the time period, January 2020-March 2022, into eight stages. We predicted Alpha, Delta, Delta Plus (AY.4.2) and Omicron (B.1.1.529) by their mutational entropy values in the Stages I, III, V and VII with accelerated paces, respectively. In late November 2021, VOC-alarm alerted that Omicron strongly competed with Delta and Delta plus to become a highly transmissible variant. Using simulated data, VOC-alarm also predicted that Omicron could lead to another COVID surge from January 2022 to March 2022. AVAILABILITY AND IMPLEMENTATION: Our software implementation is available at https://github.com/guangxujin/VOC-alarm. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

There is an urgent public health need to better understand Severe Acute Respiratory Syndrome (SARS)-CoV-2/COVID-19, particularly how sequences of the viruses could lead to diverse incidence and mortality of COVID-19 in different countries. However, because of its unknown ancestors and hosts, elucidating the genetic variations of the novel coronavirus, SARS-CoV-2, has been difficult. Without needing to know ancestors, we identified an uneven distribution of local genome similarities among the viruses categorized by geographic regions, and it was strongly correlated with incidence and mortality. To ensure unbiased and origin-independent analyses, we used a pairwise comparison of local genome sequences of virus genomes by Basic Local Alignment Search Tool (BLAST). We found a strong statistical correlation between dominance of the SARS-CoV-2 in distributions of uneven similarities and the incidence and mortality of illness. Genomic annotation of the BLAST hits also showed that viruses from geographic regions with severe infections tended to have more dynamic genomic regions in the SARS-CoV-2 receptor-binding domain (RBD) and receptor-binding motif (RBM) of the spike protein (S protein). Dynamic domains in the S protein were also confirmed by a canyon region of mismatches coincident with RBM and RBD, without hits of alignments of 100% matching. Thus, our origin-independent analysis suggests that the dynamic and unstable SARS-CoV-2-RBD could be the main reason for diverse incidence and mortality of COVID-19 infection.

MOTIVATION: Accurate detection of somatic mutations is a crucial step toward understanding cancer. Various tools have been developed to detect somatic mutations from cancer genome sequencing data by mapping reads to a universal reference genome and inferring likelihoods from complex statistical models. However, read mapping is frequently obstructed by mismatches between germline and somatic mutations on a read and the reference genome. Previous attempts to develop personalized genome tools are not compatible with downstream statistical models for somatic mutation detection. RESULTS: We present PRESM, a tool that builds personalized reference genomes by integrating germline mutations into the reference genome. The aforementioned obstacle is circumvented by using a two-step germline substitution procedure, maintaining positional fidelity using an innovative workaround. Reads derived from tumor tissue can be positioned more accurately along a personalized reference than a universal reference due to the reduced genetic distance between the subject (tumor genome) and the target (the personalized genome). Application of PRESM's personalized genome reduced false-positive (FP) somatic mutation calls by as much as 55.5%, and facilitated the discovery of a novel somatic point mutation on a germline insertion in PDE1A, a phosphodiesterase associated with melanoma. Moreover, all improvements in calling accuracy were achieved without parameter optimization, as PRESM itself is parameter-free. Hence, similar increases in read mapping and decreases in the FP rate will persist when PRESM-built genomes are applied to any user-provided dataset. AVAILABILITY AND IMPLEMENTATION: The software is available at https://github.com/precisionomics/PRESM. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SUMMARY: Systematic studies of drug repositioning require the integration of multi-level drug data, including basic chemical information (such as SMILES), drug targets, target-related signaling pathways, clinical trial information and Food and Drug Administration (FDA)-approval information, to predict new potential indications of existing drugs. Currently available databases, however, lack query support for multi-level drug information and thus are not designed to support drug repositioning studies. DrugMap Central (DMC), an online tool, is developed to help fill the gap. DMC enables the users to integrate, query, visualize, interrogate, and download multi-level data of known drugs or compounds quickly for drug repositioning studies all within one system. AVAILABILITY: DMC is accessible at http://r2d2drug.org/DMC.aspx. CONTACT: [email protected].

Recent advances in automated high-resolution fluorescence microscopy and robotic handling have made the systematic and cost effective study of diverse morphological changes within a large population of cells possible under a variety of perturbations, e.g., drugs, compounds, metal catalysts, RNA interference (RNAi). Cell population-based studies deviate from conventional microscopy studies on a few cells, and could provide stronger statistical power for drawing experimental observations and conclusions. However, it is challenging to manually extract and quantify phenotypic changes from the large amounts of complex image data generated. Thus, bioimage informatics approaches are needed to rapidly and objectively quantify and analyze the image data. This paper provides an overview of the bioimage informatics challenges and approaches in image-based studies for drug and target discovery. The concepts and capabilities of image-based screening are first illustrated by a few practical examples investigating different kinds of phenotypic changes caEditorsused by drugs, compounds, or RNAi. The bioimage analysis approaches, including object detection, segmentation, and tracking, are then described. Subsequently, the quantitative features, phenotype identification, and multidimensional profile analysis for profiling the effects of drugs and targets are summarized. Moreover, a number of publicly available software packages for bioimage informatics are listed for further reference. It is expected that this review will help readers, including those without bioimage informatics expertise, understand the capabilities, approaches, and tools of bioimage informatics and apply them to advance their own studies.

MOTIVATION: Prediction of synergistic effects of drug combinations has traditionally been relied on phenotypic response data. However, such methods cannot be used to identify molecular signaling mechanisms of synergistic drug combinations. In this article, we propose an enhanced Petri-Net (EPN) model to recognize the synergistic effects of drug combinations from the molecular response profiles, i.e. drug-treated microarray data. METHODS: We addressed the downstream signaling network of the targets for the two individual drugs used in the pairwise combinations and applied EPN to the identified targeted signaling network. In EPN, drugs and signaling molecules are assigned to different types of places, while drug doses and molecular expressions are denoted by color tokens. The changes of molecular expressions caused by treatments of drugs are simulated by two actions of EPN: firing and blasting. Firing is to transit the drug and molecule tokens from one node or place to another, and blasting is to reduce the number of molecule tokens by drug tokens in a molecule node. The goal of EPN is to mediate the state characterized by control condition without any treatment to that of treatment and to depict the drug effects on molecules by the drug tokens. RESULTS: We applied EPN to our generated pairwise drug combination microarray data. The synergistic predictions using EPN are consistent with those predicted using phenotypic response data. The molecules responsible for the synergistic effects with their associated feedback loops display the mechanisms of synergism. AVAILABILITY: The software implemented in Python 2.7 programming language is available from request. CONTACT: [email protected].

BACKGROUND: Transcriptional regulation is a fundamental process in biological systems, where transcription factors (TFs) have been revealed to play crucial roles. In recent years, in addition to TFs, an increasing number of non-coding RNAs (ncRNAs) have been shown to mediate post-transcriptional processes and regulate many critical pathways in both prokaryotes and eukaryotes. On the other hand, with more and more high-throughput biological data becoming available, it is possible and imperative to quantitatively study gene regulation in a systematic and detailed manner. RESULTS: Most existing studies for inferring transcriptional regulatory interactions and the activity of TFs ignore the possible post-transcriptional effects of ncRNAs. In this work, we propose a novel framework to infer the activity of regulators including both TFs and ncRNAs by exploring the expression profiles of target genes and (post)transcriptional regulatory relationships. We model the integrated regulatory system by a set of biochemical reactions which lead to a log-bilinear problem. The inference process is achieved by an iterative algorithm, in which two linear programming models are efficiently solved. In contrast to available related studies, the effects of ncRNAs on transcription process are considered in this work, and thus more reasonable and accurate reconstruction can be expected. In addition, the approach is suitable for large-scale problems from the viewpoint of computation. Experiments on two synthesized data sets and a model system of Escherichia coli (E. coli) carbon source transition from glucose to acetate illustrate the effectiveness of our model and algorithm. CONCLUSION: Our results show that incorporating the post-transcriptional regulation of ncRNAs into system model can mine the hidden effects from the regulation activity of TFs in transcription processes and thus can uncover the biological mechanisms in gene regulation in a more accurate manner. The software for the algorithm in this paper is available upon request.

Transcriptional regulation is a fundamental process in biological systems, where transcription factors (TFs) play crucial roles. Except for TFs, an increasing number of small non-coding RNAs (ncRNAs) have been shown to mediate post-transcriptional processes in both prokaryotes and eukaryotes. In this work, we propose a novel approach to infer the activities of regulators including TFs and ncRNAs by exploring target gene expression profiles and (post) transcriptional regulatory relationships. The inference process is efficiently achieved by an iteration algorithm, in which two linear programming models are iteratively solved. In contrast to the existing works, for the first time, the effects of ncRNAs on transcription process are considered and thus more reasonable inference can be expected. Experiments on a model system of E. coli carbon source transition from glucose to acetate illustrate the effectiveness of our method.