A Cosine Similarity-Based Method to Infer Variability of Chromatin Accessibility at the Single-Cell Level

Computational methods to explore chromatin state dynamics

The human genome is marked by several singular and combinatorial histone modifications that shape the different states of chromatin and its three-dimensional organization. Genome-wide mapping of these marks as well as histone variants and open chromatin regions is commonly carried out via profiling DNA-protein binding or via chromatin accessibility methods. After the generation of epigenomic datasets in a cell type, statistical models can be used to annotate the noncoding regions of DNA and infer the combinatorial histone marks or chromatin states (CS). These methods involve partitioning the genome and labeling individual segments based on their CS patterns. Chromatin labels enable the systematic discovery of genomic function and activity and can label the gene body, promoters or enhancers without using other genomic maps. CSs are dynamic and change under different cell conditions, such as in normal, preneoplastic or tumor cells. This review aims to explore the available computational tools that have been developed to capture CS alterations under two or more cellular conditions.

A fast and efficient algorithm for DNA sequence similarity identification

DNA sequence similarity analysis is necessary for enormous purposes including genome analysis, extracting biological information, finding the evolutionary relationship of species. There are two types of sequence analysis which are alignment-based (AB) and alignment-free (AF). AB is effective for small homologous sequences but becomes NP-hard problem for long sequences. However, AF algorithms can solve the major limitations of AB. But most of the existing AF methods show high time complexity and memory consumption, less precision, and less performance on benchmark datasets. To minimize these limitations, we develop an AF algorithm using a 2D \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k-mer$$\end{document}k-mer count matrix inspired by the CGR approach. Then we shrink the matrix by analyzing the neighbors and then measure similarities using the best combinations of pairwise distance (PD) and phylogenetic tree methods. We also dynamically choose the value of k for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k-mer$$\end{document}k-mer. We develop an efficient system for finding the positions of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k-mer$$\end{document}k-mer in the count matrix. We apply our system in six different datasets. We achieve the top rank for two benchmark datasets from AFproject, 100% accuracy for two datasets (16 S Ribosomal, 18 Eutherian), and achieve a milestone for time complexity and memory consumption in comparison to the existing study datasets (HEV, HIV-1). Therefore, the comparative results of the benchmark datasets and existing studies demonstrate that our method is highly effective, efficient, and accurate. Thus, our method can be used with the top level of authenticity for DNA sequence similarity measurement.

Targets preliminary screening for the fresh natural drug molecule based on Cosine-correlation and similarity-comparison of local network

Background

Chinese herbal medicine is made up of hundreds of natural drug molecules and has played a major role in traditional Chinese medicine (TCM) for several thousand years. Therefore, it is of great significance to study the target of natural drug molecules for exploring the mechanism of treating diseases with TCM. However, it is very difficult to determine the targets of a fresh natural drug molecule due to the complexity of the interaction between drug molecules and targets. Compared with traditional biological experiments, the computational method has the advantages of less time and low cost for targets screening, but it remains many great challenges, especially for the molecules without social ties.

Methods

This study proposed a novel method based on the Cosine-correlation and Similarity-comparison of Local Network (CSLN) to perform the preliminary screening of targets for the fresh natural drug molecules and assign weights to them through a trained parameter.

Results

The performance of CSLN is superior to the popular drug-target-interaction (DTI) prediction model GRGMF on the gold standard data in the condition that is drug molecules are the objects for training and testing. Moreover, CSLN showed excellent ability in checking the targets screening performance for a fresh-natural-drug-molecule (scenario simulation) on the TCMSP (13 positive samples in top20), meanwhile, Western-Blot also further verified the accuracy of CSLN.

Conclusions

In summary, the results suggest that CSLN can be used as an alternative strategy for screening targets of fresh natural drug molecules.

Detecting molecular subtypes from multi-omics datasets using SUMO

We present a data integration framework that uses non-negative matrix factorization of patient-similarity networks to integrate continuous multi-omics datasets for molecular subtyping. It is demonstrated to have the capability to handle missing data without using imputation and to be consistently among the best in detecting subtypes with differential prognosis and enrichment of clinical associations in a large number of cancers. When applying the approach to data from individuals with lower-grade gliomas, we identify a subtype with a significantly worse prognosis. Tumors assigned to this subtype are hypomethylated genome wide with a gain of AP-1 occupancy in demethylated distal enhancers. The tumors are also enriched for somatic chromosome 7 (chr7) gain, chr10 loss, and other molecular events that have been suggested as diagnostic markers for "IDH wild type, with molecular features of glioblastoma" by the cIMPACT-NOW consortium but have yet to be included in the World Health Organization (WHO) guidelines.

Transcription factor chromatin profiling genome-wide using uliCUT&RUN in single cells and individual blastocysts

Determining chromatin-associated protein localization across the genome has provided insight into the functions of DNA-binding proteins and their connections to disease. However, established protocols requiring large quantities of cell or tissue samples currently limit applications for clinical and biomedical research in this field. Furthermore, most technologies have been optimized to assess abundant histone protein localization, prohibiting the investigation of nonhistone protein localization in low cell numbers. We recently described a protocol to profile chromatin-associated protein localization in as low as one cell: ultra-low-input cleavage under targets and release using nuclease (uliCUT&RUN). Optimized from chromatin immunocleavage and CUT&RUN, uliCUT&RUN is a tethered enzyme-based protocol that utilizes a combination of recombinant protein, antibody recognition and stringent purification to selectively target proteins of interest and isolate the associated DNA. Performed in native conditions, uliCUT&RUN profiles protein localization to chromatin with low input and high precision. Compared with other profiling technologies, uliCUT&RUN can determine nonhistone protein chromatin occupancies in low cell numbers, permitting the investigation into the molecular functions of a range of DNA-binding proteins within rare samples. From sample preparation to sequencing library submission, the uliCUT&RUN protocol takes <2 d to perform, with the accompanying data analysis timeline dependent on experience level.

Single Cell Omics for Musculoskeletal Research

PURPOSE OF REVIEW

The ability to analyze the molecular events occurring within individual cells as opposed to populations of cells is revolutionizing our understanding of musculoskeletal tissue development and disease. Single cell studies have the great potential of identifying cellular subpopulations that work in a synchronized fashion to regenerate and repair damaged tissues during normal homeostasis. In addition, such studies can elucidate how these processes break down in disease as well as identify cellular subpopulations that drive the disease. This review highlights three emerging technologies: single cell RNA sequencing (scRNA-seq), Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq), and Cytometry by Time-Of-Flight (CyTOF) mass cytometry.

RECENT FINDINGS

Technological and bioinformatic tools to analyze the transcriptome, epigenome, and proteome at the individual cell level have advanced rapidly making data collection relatively easy; however, understanding how to access and interpret the data remains a challenge for many scientists. It is, therefore, of paramount significance to educate the musculoskeletal community on how single cell technologies can be used to answer research questions and advance translation. This article summarizes talks given during a workshop on "Single Cell Omics" at the 2020 annual meeting of the Orthopedic Research Society. Studies that applied scRNA-seq, ATAC-seq, and CyTOF mass cytometry to cartilage development and osteoarthritis are reviewed. This body of work shows how these cutting-edge tools can advance our understanding of the cellular heterogeneity and trajectories of lineage specification during development and disease.

Lung organoid simulations for modelling and predicting the effect of mutations on SARS-CoV-2 infectivity

The global pandemic caused by the SARS-CoV-2 virus continues to spread. Infection with SARS- CoV-2 causes COVID-19, a disease of variable severity. Mutation has already altered the SARS-CoV-2 genome from its original reported sequence and continued mutation is highly probable. These mutations can: (i) have no significant impact (they are silent), (ii) result in a complete loss or reduction of infectivity, or (iii) induce increase in infectivity. Physical generation, for research purposes, of viral mutations that could enhance infectivity are controversial and highly regulated. The primary purpose of this project was to evaluate the ability of the DeepNEU machine learning stem-cell simulation platform to enable rapid and efficient assessment of the potential impact of viral loss-of-function (LOF) and gain-of-function (GOF) mutations on SARS-CoV-2 infectivity. Our data suggest that SARS-CoV-2 infection can be simulated in human alveolar type lung cells. Simulation of infection in these lung cells can be used to model and assess the impact of LOF and GOF mutations in the SARS-CoV2 genome. We have also created a four- factor infectivity measure: the DeepNEU Case Fatality Rate (dnCFR). dnCFR can be used to assess infectivity based on the presence or absence of the key viral proteins (NSP3, Spike-RDB, N protein, and M protein). dnCFR was used in this study, not to only assess the impact of different mutations on SARS-CoV2 infectivity, but also to categorize the effects of mutations as loss of infectivity or gain of infectivity events.

Profiling Chromatin Accessibility at Single-cell Resolution

How distinct transcriptional programs are enacted to generate cellular heterogeneity and plasticity, and enable complex fate decisions are important open questions. One key regulator is the cell’s epigenome state that drives distinct transcriptional programs by regulating chromatin accessibility. Genome-wide chromatin accessibility measurements can impart insights into regulatory sequences (in)accessible to DNA-binding proteins at a single-cell resolution. This review outlines molecular methods and bioinformatic tools for capturing cell-to-cell chromatin variation using single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) in a scalable fashion. It also covers joint profiling of chromatin with transcriptome/proteome measurements, computational strategies to integrate multi-omic measurements, and predictive bioinformatic tools to infer chromatin accessibility from single-cell transcriptomic datasets. Methodological refinements that increase power for cell discovery through robust chromatin coverage and integrate measurements from multiple modalities will further expand our understanding of gene regulation during homeostasis and disease.

Single-cell ATAC-seq signal extraction and enhancement with SCATE

Single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) is the state-of-the-art technology for analyzing genome-wide regulatory landscapes in single cells. Single-cell ATAC-seq data are sparse and noisy, and analyzing such data is challenging. Existing computational methods cannot accurately reconstruct activities of individual cis-regulatory elements (CREs) in individual cells or rare cell subpopulations. We present a new statistical framework, SCATE, that adaptively integrates information from co-activated CREs, similar cells, and publicly available regulome data to substantially increase the accuracy for estimating activities of individual CREs. We demonstrate that SCATE can be used to better reconstruct the regulatory landscape of a heterogeneous sample.

Application of information theoretical approaches to assess diversity and similarity in single-cell transcriptomics

Single-cell transcriptomics offers a powerful way to reveal the heterogeneity of individual cells. To date, many information theoretical approaches have been proposed to assess diversity and similarity, and characterize the latent heterogeneity in transcriptome data. Diversity implies gene expression variations and can facilitate the identification of signature genes; while, similarity unravels co-expression patterns for cell type clustering. In this review, we summarized 16 measures of information theory used for evaluating diversity and similarity in single-cell transcriptomic data, provide references and shed light on selected theoretical properties when there is a need to select proper measurements in general cases. We further provide an R package assembling discussed approaches to improve the researchers own single-cell transcriptome study. At last, we prospected further applications of diversity and similarity measures in support of depicting heterogeneity in single-cell multi-omics data.

Joint profiling of chromatin accessibility and CAR-T integration site analysis at population and single-cell levels

Chimeric antigen receptor (CAR)-T immunotherapy has yielded impressive results in several B cell malignancies, establishing itself as a powerful means to redirect the natural properties of T lymphocytes. In this strategy, the T cell genome is modified by the integration of lentiviral vectors encoding CAR that direct tumor cell killing. However, this therapeutic approach is often limited by the extent of CAR-T cell expansion in vivo. A major outstanding question is whether or not CAR-T integration itself enhances the proliferative competence of individual T cells by rewiring their regulatory landscape. To address this question, it is critical to define the identity of an individual CAR-T cell and simultaneously chart where the CAR-T vector integrates into the genome. Here, we report the development of a method called EpiVIA (https://github.com/VahediLab/epiVIA) for the joint profiling of the chromatin accessibility and lentiviral integration site analysis at the population and single-cell levels. We validate our technique in clonal cells with previously defined integration sites and further demonstrate the ability to measure lentiviral integration sites and chromatin accessibility of host and viral genomes at the single-cell resolution in CAR-T cells. We anticipate that EpiVIA will enable the single-cell deconstruction of gene regulation during CAR-T therapy, leading to the discovery of cellular factors associated with durable treatment.

A multi-level similarity measure for the retrieval of the common CT imaging signs of lung diseases

The common CT imaging signs of lung diseases (CISLs) which frequently appear in lung CT images are widely used in the diagnosis of lung diseases. Computer-aided diagnosis (CAD) based on the CISLs can improve radiologists' performance in the diagnosis of lung diseases. Since similarity measure is important for CAD, we propose a multi-level method to measure the similarity between the CISLs. The CISLs are characterized in the low-level visual scale, mid-level attribute scale, and high-level semantic scale, for a rich representation. The similarity at multiple levels is calculated and combined in a weighted sum form as the final similarity. The proposed multi-level similarity method is capable of computing the level-specific similarity and optimal cross-level complementary similarity. The effectiveness of the proposed similarity measure method is evaluated on a dataset of 511 lung CT images from clinical patients for CISLs retrieval. It can achieve about 80% precision and take only 3.6 ms for the retrieval process. The extensive comparative evaluations on the same datasets are conducted to validate the advantages on retrieval performance of our multi-level similarity measure over the single-level measure and the two-level similarity methods. The proposed method can have wide applications in radiology and decision support. Graphical abstract.