BloodChIP Xtra: an expanded database of comparative genome-wide transcription factor binding and gene-expression profiles in healthy human stem/progenitor subsets and leukemic cells.

The BloodChIP Xtra database (http://bloodchipXtra.vafaeelab.com/) facilitates genome-wide exploration and visualization of transcription factor (TF) occupancy and chromatin configuration in rare primary human hematopoietic stem (HSC-MPP) and progenitor (CMP, GMP, MEP) cells and acute myeloid leukemia (AML) cell lines (KG-1, ME-1, Kasumi1, TSU-1621-MT), along with chromatin accessibility and gene expression data from these and primary patient AMLs. BloodChIP Xtra features significantly more datasets than our earlier database BloodChIP (two primary cell types and two cell lines). Improved methodologies for determining TF occupancy and chromatin accessibility have led to increased availability of data for rare primary cell types across the spectrum of healthy and AML hematopoiesis. However, there is a continuing need for these data to be integrated in an easily accessible manner for gene-based queries and use in downstream applications. Here, we provide a user-friendly database based around genome-wide binding profiles of key hematopoietic TFs and histone marks in healthy stem/progenitor cell types. These are compared with binding profiles and chromatin accessibility derived from primary and cell line AML and integrated with expression data from corresponding cell types. All queries can be exported to construct TF-gene and protein-protein networks and evaluate the association of genes with specific cellular processes.

Genome-wide transcription factor-binding maps reveal cell-specific changes in the regulatory architecture of human HSPCs.

Hematopoietic stem and progenitor cells (HSPCs) rely on a complex interplay among transcription factors (TFs) to regulate differentiation into mature blood cells. A heptad of TFs (FLI1, ERG, GATA2, RUNX1, TAL1, LYL1, LMO2) bind regulatory elements in bulk CD34+ HSPCs. However, whether specific heptad-TF combinations have distinct roles in regulating hematopoietic differentiation remains unknown. We mapped genome-wide chromatin contacts (HiC, H3K27ac, HiChIP), chromatin modifications (H3K4me3, H3K27ac, H3K27me3) and 10 TF binding profiles (heptad, PU.1, CTCF, STAG2) in HSPC subsets (stem/multipotent progenitors plus common myeloid, granulocyte macrophage, and megakaryocyte erythrocyte progenitors) and found TF occupancy and enhancer-promoter interactions varied significantly across cell types and were associated with cell-type-specific gene expression. Distinct regulatory elements were enriched with specific heptad-TF combinations, including stem-cell-specific elements with ERG, and myeloid- and erythroid-specific elements with combinations of FLI1, RUNX1, GATA2, TAL1, LYL1, and LMO2. Furthermore, heptad-occupied regions in HSPCs were subsequently bound by lineage-defining TFs, including PU.1 and GATA1, suggesting that heptad factors may prime regulatory elements for use in mature cell types. We also found that enhancers with cell-type-specific heptad occupancy shared a common grammar with respect to TF binding motifs, suggesting that combinatorial binding of TF complexes was at least partially regulated by features encoded in DNA sequence motifs. Taken together, this study comprehensively characterizes the gene regulatory landscape in rare subpopulations of human HSPCs. The accompanying data sets should serve as a valuable resource for understanding adult hematopoiesis and a framework for analyzing aberrant regulatory networks in leukemic cells.

Blood-based transcriptomic signature panel identification for cancer diagnosis: benchmarking of feature extraction methods.

Liquid biopsy has shown promise for cancer diagnosis due to its minimally invasive nature and the potential for novel biomarker discovery. However, the low concentration of relevant blood-based biosources and the heterogeneity of samples (i.e. the variability of relative abundance of molecules identified), pose major challenges to biomarker discovery. Moreover, the number of molecular measurements or features (e.g. transcript read counts) per sample could be in the order of several thousand, whereas the number of samples is often substantially lower, leading to the curse of dimensionality. These challenges, among others, elucidate the importance of a robust biomarker panel identification or feature extraction step wherein relevant molecular measurements are identified prior to classification for cancer detection. In this work, we performed a benchmarking study on 12 feature extraction methods using transcriptomic profiles derived from different blood-based biosources. The methods were assessed both in terms of their predictive performance and the robustness of the biomarker panels in diagnosing cancer or stratifying cancer subtypes. While performing the comparison, the feature extraction methods are categorized into feature subset selection methods and transformation methods. A transformation feature extraction method, namely partial least square discriminant analysis, was found to perform consistently superior in terms of classification performance. As part of the benchmarking study, a generic pipeline has been created and made available as an R package to ensure reproducibility of the results and allow for easy extension of this study to other datasets (https://github.com/VafaeeLab/bloodbased-pancancer-diagnosis).

Disentangling single-cell omics representation with a power spectral density-based feature extraction.

Emerging single-cell technologies provide high-resolution measurements of distinct cellular modalities opening new avenues for generating detailed cellular atlases of many and diverse tissues. The high dimensionality, sparsity, and inaccuracy of single cell sequencing measurements, however, can obscure discriminatory information, mask cellular subtype variations and complicate downstream analyses which can limit our understanding of cell function and tissue heterogeneity. Here, we present a novel pre-processing method (scPSD) inspired by power spectral density analysis that enhances the accuracy for cell subtype separation from large-scale single-cell omics data. We comprehensively benchmarked our method on a wide range of single-cell RNA-sequencing datasets and showed that scPSD pre-processing, while being fast and scalable, significantly reduces data complexity, enhances cell-type separation, and enables rare cell identification. Additionally, we applied scPSD to transcriptomics and chromatin accessibility cell atlases and demonstrated its capacity to discriminate over 100 cell types across the whole organism and across different modalities of single-cell omics data.

Supervised application of internal validation measures to benchmark dimensionality reduction methods in scRNA-seq data.

A typical single-cell RNA sequencing (scRNA-seq) experiment will measure on the order of 20 000 transcripts and thousands, if not millions, of cells. The high dimensionality of such data presents serious complications for traditional data analysis methods and, as such, methods to reduce dimensionality play an integral role in many analysis pipelines. However, few studies have benchmarked the performance of these methods on scRNA-seq data, with existing comparisons assessing performance via downstream analysis accuracy measures, which may confound the interpretation of their results. Here, we present the most comprehensive benchmark of dimensionality reduction methods in scRNA-seq data to date, utilizing over 300 000 compute hours to assess the performance of over 25 000 low-dimension embeddings across 33 dimensionality reduction methods and 55 scRNA-seq datasets. We employ a simple, yet novel, approach, which does not rely on the results of downstream analyses. Internal validation measures (IVMs), traditionally used as an unsupervised method to assess clustering performance, are repurposed to measure how well-formed biological clusters are after dimensionality reduction. Performance was further evaluated over nearly 200 000 000 iterations of DBSCAN, a density-based clustering algorithm, showing that hyperparameter optimization using IVMs as the objective function leads to near-optimal clustering. Methods were also assessed on the extent to which they preserve the global structure of the data, and on their computational memory and time requirements across a large range of sample sizes. Our comprehensive benchmarking analysis provides a valuable resource for researchers and aims to guide best practice for dimensionality reduction in scRNA-seq analyses, and we highlight Latent Dirichlet Allocation and Potential of Heat-diffusion for Affinity-based Transition Embedding as high-performing algorithms.

A comprehensive integrated drug similarity resource for in-silico drug repositioning and beyond.

Drug similarity studies are driven by the hypothesis that similar drugs should display similar therapeutic actions and thus can potentially treat a similar constellation of diseases. Drug-drug similarity has been derived by variety of direct and indirect sources of evidence and frequently shown high predictive power in discovering validated repositioning candidates as well as other in-silico drug development applications. Yet, existing resources either have limited coverage or rely on an individual source of evidence, overlooking the wealth and diversity of drug-related data sources. Hence, there has been an unmet need for a comprehensive resource integrating diverse drug-related information to derive multi-evidenced drug-drug similarities. We addressed this resource gap by compiling heterogenous information for an exhaustive set of small-molecule drugs (total of 10 367 in the current version) and systematically integrated multiple sources of evidence to derive a multi-modal drug-drug similarity network. The resulting database, 'DrugSimDB' currently includes 238 635 drug pairs with significant aggregated similarity, complemented with an interactive user-friendly web interface (http://vafaeelab.com/drugSimDB.html), which not only enables database ease of access, search, filtration and export, but also provides a variety of complementary information on queried drugs and interactions. The integration approach can flexibly incorporate further drug information into the similarity network, providing an easily extendable platform. The database compilation and construction source-code has been well-documented and semi-automated for any-time upgrade to account for new drugs and up-to-date drug information.

Pan-Cancer Analysis Reveals Functional Similarity of Three lncRNAs across Multiple Tumors.

Long non-coding RNAs (lncRNAs) are emerging as key regulators in many biological processes. The dysregulation of lncRNA expression has been associated with many diseases, including cancer. Mounting evidence suggests lncRNAs to be involved in cancer initiation, progression, and metastasis. Thus, understanding the functional implications of lncRNAs in tumorigenesis can aid in developing novel biomarkers and therapeutic targets. Rich cancer datasets, documenting genomic and transcriptomic alterations together with advancement in bioinformatics tools, have presented an opportunity to perform pan-cancer analyses across different cancer types. This study is aimed at conducting a pan-cancer analysis of lncRNAs by performing differential expression and functional analyses between tumor and non-neoplastic adjacent samples across eight cancer types. Among dysregulated lncRNAs, seven were shared across all cancer types. We focused on three lncRNAs, found to be consistently dysregulated among tumors. It has been observed that these three lncRNAs of interest are interacting with a wide range of genes across different tissues, yet enriching substantially similar biological processes, found to be implicated in cancer progression and proliferation.

Investigation of Plasma-Derived Lipidome Profiles in Experimental Cerebral Malaria in a Mouse Model Study.

Cerebral malaria (CM), a fatal complication of Plasmodium infection that affects children, especially under the age of five, in sub-Saharan Africa and adults in South-East Asia, results from incompletely understood pathogenetic mechanisms. Increased release of circulating miRNA, proteins, lipids and extracellular vesicles has been found in CM patients and experimental mouse models. We compared lipid profiles derived from the plasma of CBA mice infected with Plasmodium berghei ANKA (PbA), which causes CM, to those from Plasmodium yoelii (Py), which does not. We previously showed that platelet-free plasma (18k fractions enriched from plasma) contains a high number of extracellular vesicles (EVs). Here, we found that this fraction produced at the time of CM differed dramatically from those of non-CM mice, despite identical levels of parasitaemia. Using high-resolution liquid chromatography-mass spectrometry (LCMS), we identified over 300 lipid species within 12 lipid classes. We identified 45 and 75 lipid species, mostly including glycerolipids and phospholipids, with significantly altered concentrations in PbA-infected mice compared to Py-infected and uninfected mice, respectively. Total lysophosphatidylethanolamine (LPE) levels were significantly lower in PbA infection compared to Py infection and controls. These results suggest that experimental CM could be characterised by specific changes in the lipid composition of the 18k fraction containing circulating EVs and can be considered an appropriate model to study the role of lipids in the pathophysiology of CM.

TDAview: an online visualization tool for topological data analysis.

SUMMARY: TDAview is an online tool for topological data analysis (TDA) and visualization. It implements the Mapper algorithm for TDA and provides extensive graph visualization options. TDAview is a user-friendly tool that allows biologists and clinicians without programming knowledge to harness the power of TDA. TDAview supports an analysis and visualization mode in which a Mapper graph is constructed based on user-specified parameters, followed by graph visualization. It can also be used in a visualization only mode in which TDAview is used for visualizing the data properties of a Mapper graph generated using other open-source software. The graph visualization options allow data exploration by graphical display of metadata variable values for nodes and edges, as well as the generation of publishable figures. TDAview can handle large datasets, with tens of thousands of data points, and thus has a wide range of applications for high-dimensional data, including the construction of topology-based gene co-expression networks. AVAILABILITY AND IMPLEMENTATION: TDAview is a free online tool available at https://voineagulab.github.io/TDAview/. The source code, usage documentation and example data are available at TDAview GitHub repository: https://github.com/Voineagulab/TDAview.

Human Group IIA Phospholipase A2-Three Decades on from Its Discovery.

Phospholipase A2 (PLA2) enzymes were first recognized as an enzyme activity class in 1961. The secreted (sPLA2) enzymes were the first of the five major classes of human PLA2s to be identified and now number nine catalytically-active structurally homologous proteins. The best-studied of these, group IIA sPLA2, has a clear role in the physiological response to infection and minor injury and acts as an amplifier of pathological inflammation. The enzyme has been a target for anti-inflammatory drug development in multiple disorders where chronic inflammation is a driver of pathology since its cloning in 1989. Despite intensive effort, no clinically approved medicines targeting the enzyme activity have yet been developed. This review catalogues the major discoveries in the human group IIA sPLA2 field, focusing on features of enzyme function that may explain this lack of success and discusses future research that may assist in realizing the potential benefit of targeting this enzyme. Functionally-selective inhibitors together with isoform-selective inhibitors are necessary to limit the apparent toxicity of previous drugs. There is also a need to define the relevance of the catalytic function of hGIIA to human inflammatory pathology relative to its recently-discovered catalysis-independent function.

Expression of long noncoding RNAs in cancer-associated fibroblasts linked to patient survival in ovarian cancer.

Cancer-associated fibroblasts (CAFs) are the most abundant cell type in the tumor microenvironment and are responsible for producing the desmoplastic reaction that is a poor prognostic factor in ovarian cancer. Long non-coding RNAs (lncRNAs) have been shown to play important roles in cancer. However, very little is known about the role of lncRNAs in the tumor microenvironment. We aimed to identify lncRNAs expressed in ovarian CAFs that were associated with patient survival and used computational approaches to predict their function. Increased expression of 9 lncRNAs and decreased expression of 1 lncRNA in ovarian CAFs were found to be associated with poorer overall survival. A "guilt-by-association" approach was used to predict the function of these lncRNAs. In particular, MIR155HG was predicted to play a role in immune response. Further investigation revealed high MIR155HG expression to be associated with higher infiltrates of immune cell subsets. In conclusion, these data indicate expression on several lncRNAs in CAFs are associated with patient survival and are likely to play an important role in regulating CAF function.

In silico prediction of physical protein interactions and characterization of interactome orphans.

Protein-protein interactions (PPIs) are useful for understanding signaling cascades, predicting protein function, associating proteins with disease and fathoming drug mechanism of action. Currently, only ∼ 10% of human PPIs may be known, and about one-third of human proteins have no known interactions. We introduce FpClass, a data mining-based method for proteome-wide PPI prediction. At an estimated false discovery rate of 60%, we predicted 250,498 PPIs among 10,531 human proteins; 10,647 PPIs involved 1,089 proteins without known interactions. We experimentally tested 233 high- and medium-confidence predictions and validated 137 interactions, including seven novel putative interactors of the tumor suppressor p53. Compared to previous PPI prediction methods, FpClass achieved better agreement with experimentally detected PPIs. We provide an online database of annotated PPI predictions (http://ophid.utoronto.ca/fpclass/) and the prediction software (http://www.cs.utoronto.ca/~juris/data/fpclass/).

Unraveling the mechanistic complexity of Alzheimer's disease through systems biology.

Alzheimer's disease (AD) is a complex, multifactorial disease that has reached global epidemic proportions. The challenge remains to fully identify its underlying molecular mechanisms that will enable development of accurate diagnostic tools and therapeutics. Conventional experimental approaches that target individual or small sets of genes or proteins may overlook important parts of the regulatory network, which limits the opportunity of identifying multitarget interventions. Our perspective is that a more complete insight into potential treatment options for AD will only be made possible through studying the disease as a system. We propose an integrative systems biology approach that we argue has been largely untapped in AD research. We present key publications to demonstrate the value of this approach and discuss the potential to intensify research efforts in AD through transdisciplinary collaboration. We highlight challenges and opportunities for significant breakthroughs that could be made if a systems biology approach is fully exploited.

3D-IncNet: Head and Neck (H&N) Primary Tumors Segmentation and Survival Prediction.

Cancer begins when healthy cells change and grow out of control, forming a mass called a tumor. Head and neck (H&N) cancers usually develop in or around the head and neck, including the mouth (oral cavity), nose and sinuses, throat (pharynx), and voice box (larynx). 4% of all cancers are H&N cancers with a very low survival rate (a five-year survival rate of 64.7%). FDG-PET/CT imaging is often used for early diagnosis and staging of H&N tumors, thus improving these patients' survival rates. This work presents a novel 3D-Inception-Residual aided with 3D depth-wise convolution and squeeze and excitation block. We introduce a 3D depth-wise convolution-inception encoder consisting of an additional 3D squeeze and excitation block and a 3D depth-wise convolution-based residual learning decoder (3D-IncNet), which not only helps to recalibrate the channel-wise features but adaptively through explicit inter-dependencies modeling but also integrate the coarse and fine features resulting in accurate tumor segmentation. We further demonstrate the effectiveness of inception-residual encoder-decoder architecture in achieving better dice scores and the impact of depth-wise convolution in lowering the computational cost. We applied random forest for survival prediction on deep, clinical, and radiomics features. Experiments are conducted on the benchmark HECKTOR21 challenge, which showed significantly better performance by surpassing the state-of-the-artwork and achieved 0.836 and 0.811 concordance index and dice scores, respectively. We made the model and code publicly available.

AICRF: ancestry inference of admixed population with deep conditional random field.

Ancestry inference of admixed populations is an important issue in anthropology and studies of gene discovery, and characterization. Usually, local ancestor inference (LAI) methods use fixed-length windows to divide chromosomes into smaller blocks. The accuracy of LAI algorithms will decrease if a window with an inappropriate length is used to infer the ancestry of admixed individuals. In this study, we first present a heuristic function to determine a proper window length for LAI methods. This heuristic is based on the distance between the ancestral populations of admixed individuals. Then we introduce a method for ancestry inference of admixed population with deep conditional random field (AICRF). AICRF uses a conditional random field (CRF) parameterized by probable extreme learning machines (PELMs) trained on reference panels where PELM is a novel probabilistic ELM classifier. This method does not require many statistical or biological parameters. We evaluate the performance of AICRF in comparison with RFMix. Experimental results show that AICRF is more accurate than RFMix with increasing admixture times.

An introduction to representation learning for single-cell data analysis.

Single-cell-resolved systems biology methods, including omics- and imaging-based measurement modalities, generate a wealth of high-dimensional data characterizing the heterogeneity of cell populations. Representation learning methods are routinely used to analyze these complex, high-dimensional data by projecting them into lower-dimensional embeddings. This facilitates the interpretation and interrogation of the structures, dynamics, and regulation of cell heterogeneity. Reflecting their central role in analyzing diverse single-cell data types, a myriad of representation learning methods exist, with new approaches continually emerging. Here, we contrast general features of representation learning methods spanning statistical, manifold learning, and neural network approaches. We consider key steps involved in representation learning with single-cell data, including data pre-processing, hyperparameter optimization, downstream analysis, and biological validation. Interdependencies and contingencies linking these steps are also highlighted. This overview is intended to guide researchers in the selection, application, and optimization of representation learning strategies for current and future single-cell research applications.

Biomarker Reproducibility Challenge: A Review of Non-Nucleotide Biomarker Discovery Protocols from Body Fluids in Breast Cancer Diagnosis.

Breast cancer has now become the most commonly diagnosed cancer, accounting for one in eight cancer diagnoses worldwide. Non-invasive diagnostic biomarkers and associated tests are superlative candidates to complement or improve current approaches for screening, early diagnosis, or prognosis of breast cancer. Biomarkers detected from body fluids such as blood (serum/plasma), urine, saliva, nipple aspiration fluid, and tears can detect breast cancer at its early stages in a minimally invasive way. The advancements in high-throughput molecular profiling (omics) technologies have opened an unprecedented opportunity for unbiased biomarker detection. However, the irreproducibility of biomarkers and discrepancies of reported markers have remained a major roadblock to clinical implementation, demanding the investigation of contributing factors and the development of standardised biomarker discovery pipelines. A typical biomarker discovery workflow includes pre-analytical, analytical, and post-analytical phases, from sample collection to model development. Variations introduced during these steps impact the data quality and the reproducibility of the findings. Here, we present a comprehensive review of methodological variations in biomarker discovery studies in breast cancer, with a focus on non-nucleotide biomarkers (i.e., proteins, lipids, and metabolites), highlighting the pre-analytical to post-analytical variables, which may affect the accurate identification of biomarkers from body fluids.

Small and Large Extracellular Vesicles Derived from Pleural Mesothelioma Cell Lines Offer Biomarker Potential.

Pleural mesothelioma, previously known as malignant pleural mesothelioma, is an aggressive and fatal cancer of the pleura, with one of the poorest survival rates. Pleural mesothelioma is in urgent clinical need for biomarkers to aid early diagnosis, improve prognostication, and stratify patients for treatment. Extracellular vesicles (EVs) have great potential as biomarkers; however, there are limited studies to date on their role in pleural mesothelioma. We conducted a comprehensive proteomic analysis on different EV populations derived from five pleural mesothelioma cell lines and an immortalized control cell line. We characterized three subtypes of EVs (10 K, 18 K, and 100 K), and identified a total of 4054 unique proteins. Major differences were found in the cargo between the three EV subtypes. We show that 10 K EVs were enriched in mitochondrial components and metabolic processes, while 18 K and 100 K EVs were enriched in endoplasmic reticulum stress. We found 46 new cancer-associated proteins for pleural mesothelioma, and the presence of mesothelin and PD-L1/PD-L2 enriched in 100 K and 10 K EV, respectively. We demonstrate that different EV populations derived from pleural mesothelioma cells have unique cancer-specific proteomes and carry oncogenic cargo, which could offer a novel means to extract biomarkers of interest for pleural mesothelioma from liquid biopsies.

Computer-Aided Diagnosis of Melanoma Subtypes Using Reflectance Confocal Images.

Lentigo maligna (LM) is an early form of pre-invasive melanoma that predominantly affects sun-exposed areas such as the face. LM is highly treatable when identified early but has an ill-defined clinical border and a high rate of recurrence. Atypical intraepidermal melanocytic proliferation (AIMP), also known as atypical melanocytic hyperplasia (AMH), is a histological description that indicates melanocytic proliferation with uncertain malignant potential. Clinically and histologically, AIMP can be difficult to distinguish from LM, and indeed AIMP may, in some cases, progress to LM. The early diagnosis and distinction of LM from AIMP are important since LM requires a definitive treatment. Reflectance confocal microscopy (RCM) is an imaging technique often used to investigate these lesions non-invasively, without biopsy. However, RCM equipment is often not readily available, nor is the associated expertise for RCM image interpretation easy to find. Here, we implemented a machine learning classifier using popular convolutional neural network (CNN) architectures and demonstrated that it could correctly classify lesions between LM and AIMP on biopsy-confirmed RCM image stacks. We identified local z-projection (LZP) as a recent fast approach for projecting a 3D image into 2D while preserving information and achieved high-accuracy machine classification with minimal computational requirements.

Extracellular vesicles from biological fluids as potential markers in castration resistant prostate cancer.

PURPOSE: Extracellular vesicles (EV) secreted from cancer cells are present in various biological fluids, carrying distinctly different cellular components compared to normal cells, and have great potential to be used as markers for disease initiation, progression, and response to treatment. This under-utilised tool provides insights into a better understanding of prostate cancer. METHODS: EV from serum and urine of healthy men and castration-resistant prostate cancer (CRPC) patients were isolated and characterised by transmission electron microscopy, particle size analysis, and western blot. Proteomic and cholesterol liquid chromatography-mass spectrometry (LC-MS) analyses were conducted. RESULTS: There was a successful enrichment of small EV/exosomes isolated from serum and urine. EV derived from biological fluids of CRPC patients had significant differences in composition when compared with those from healthy controls. Analysis of matched serum and urine samples from six prostate cancer patients revealed specific EV proteins common in both types of biological fluid for each patient. CONCLUSION: Some of the EV proteins identified from our analyses have potential to be used as CRPC markers. These markers may depict a pattern in cancer progression through non-invasive sample collection.