Signed Distance Correlation (SiDCo): an online implementation of distance correlation and partial distance correlation for data-driven network analysis.

MOTIVATION: There is a need for easily accessible implementations that measure the strength of both linear and non-linear relationships between metabolites in biological systems as an approach for data-driven network development. While multiple tools implement linear Pearson and Spearman methods, there are no such tools that assess distance correlation. RESULTS: We present here SIgned Distance COrrelation (SiDCo). SiDCo is a GUI platform for calculation of distance correlation in omics data, measuring linear and non-linear dependencies between variables, as well as correlation between vectors of different lengths, e.g. different sample sizes. By combining the sign of the overall trend from Pearson's correlation with distance correlation values, we further provide a novel "signed distance correlation" of particular use in metabolomic and lipidomic analyses. Distance correlations can be selected as one-to-one or one-to-all correlations, showing relationships between each feature and all other features one at a time or in combination. Additionally, we implement "partial distance correlation," calculated using the Gaussian Graphical model approach adapted to distance covariance. Our platform provides an easy-to-use software implementation that can be applied to the investigation of any dataset. AVAILABILITY AND IMPLEMENTATION: The SiDCo software application is freely available at https://complimet.ca/sidco. Supplementary help pages are provided at https://complimet.ca/sidco. Supplementary Material shows an example of an application of SiDCo in metabolomics.

Differential Expression Feature Extraction (DEFE): A Case Study in Wheat FHB RNA-Seq Data Analysis.

In differential gene expression data analysis, one objective is to identify groups of co-expressed genes from a large dataset in order to detect the association between such a group of genes and an experimental condition. This is often done through a clustering approach, such as k-means or bipartition hierarchical clustering, based on particular similarity measures in the grouping process. In such a dataset, the gene differential expression itself is an innate attribute that can be used in the feature extraction process. For example, in a dataset consisting of multiple treatments versus their controls, the expression of a gene in each treatment would have three possible behaviors, upregulated, downregulated, or unchanged. We present in this chapter, a differential expression feature extraction (DEFE) method by using a string consisting of three numerical values at each character to denote such behavior, i.e., 1 = up, 2 = down, and 0 = unchanged, which results in up to 3B differential expression patterns across all B comparisons. This approach has been successfully applied in many research projects, and among these, we demonstrate the strength of DEFE in a case study on RNA-sequencing (RNA-seq) data analysis of wheat challenged with the phytopathogenic fungus, Fusarium graminearum. Combinations of multiple schemes of DEFE patterns revealed groups of genes putatively associated with resistance or susceptibility to FHB.

Dementia with Lewy bodies post-mortem brains reveal differentially methylated CpG sites with biomarker potential.

Dementia with Lewy bodies (DLB) is a common form of dementia with known genetic and environmental interactions. However, the underlying epigenetic mechanisms which reflect these gene-environment interactions are poorly studied. Herein, we measure genome-wide DNA methylation profiles of post-mortem brain tissue (Broadmann area 7) from 15 pathologically confirmed DLB brains and compare them with 16 cognitively normal controls using Illumina MethylationEPIC arrays. We identify 17 significantly differentially methylated CpGs (DMCs) and 17 differentially methylated regions (DMRs) between the groups. The DMCs are mainly located at the CpG islands, promoter and first exon regions. Genes associated with the DMCs are linked to "Parkinson's disease" and "metabolic pathway", as well as the diseases of "severe intellectual disability" and "mood disorders". Overall, our study highlights previously unreported DMCs offering insights into DLB pathogenesis with the possibility that some of these could be used as biomarkers of DLB in the future.

METAbolomics data Balancing with Over-sampling Algorithms (META-BOA): an online resource for addressing class imbalance.

MOTIVATION: Class imbalance, or unequal sample sizes between classes, is an increasing concern in machine learning for metabolomic and lipidomic data mining, which can result in overfitting for the over-represented class. Numerous methods have been developed for handling class imbalance, but they are not readily accessible to users with limited computational experience. Moreover, there is no resource that enables users to easily evaluate the effect of different over-sampling algorithms. RESULTS: METAbolomics data Balancing with Over-sampling Algorithms (META-BOA) is a web-based application that enables users to select between four different methods for class balancing, followed by data visualization and classification of the sample to observe the augmentation effects. META-BOA outputs a newly balanced dataset, generating additional samples in the minority class, according to the user's choice of Synthetic Minority Over-sampling Technique (SMOTE), Borderline-SMOTE (BSMOTE), Adaptive Synthetic (ADASYN) or Random Over-Sampling Examples (ROSE). To present the effect of over-sampling on the data META-BOA further displays both principal component analysis and t-distributed stochastic neighbor embedding visualization of data pre- and post-over-sampling. Random forest classification is utilized to compare sample classification in both the original and balanced datasets, enabling users to select the most appropriate method for their further analyses. AVAILABILITY AND IMPLEMENTATION: META-BOA is available at https://complimet.ca/meta-boa. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

BATL: Bayesian annotations for targeted lipidomics.

MOTIVATION: Bioinformatic tools capable of annotating, rapidly and reproducibly, large, targeted lipidomic datasets are limited. Specifically, few programs enable high-throughput peak assessment of liquid chromatography-electrospray ionization tandem mass spectrometry data acquired in either selected or multiple reaction monitoring modes. RESULTS: We present here Bayesian Annotations for Targeted Lipidomics, a Gaussian naïve Bayes classifier for targeted lipidomics that annotates peak identities according to eight features related to retention time, intensity, and peak shape. Lipid identification is achieved by modeling distributions of these eight input features across biological conditions and maximizing the joint posterior probabilities of all peak identities at a given transition. When applied to sphingolipid and glycerophosphocholine selected reaction monitoring datasets, we demonstrate over 95% of all peaks are rapidly and correctly identified. AVAILABILITY AND IMPLEMENTATION: BATL software is freely accessible online at https://complimet.ca/batl/ and is compatible with Safari, Firefox, Chrome and Edge. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Metabolomics and computational analysis of the role of monoamine oxidase activity in delirium and SARS-COV-2 infection.

Delirium is an acute change in attention and cognition occurring in ~ 65% of severe SARS-CoV-2 cases. It is also common following surgery and an indicator of brain vulnerability and risk for the development of dementia. In this work we analyzed the underlying role of metabolism in delirium-susceptibility in the postoperative setting using metabolomic profiling of cerebrospinal fluid and blood taken from the same patients prior to planned orthopaedic surgery. Distance correlation analysis and Random Forest (RF) feature selection were used to determine changes in metabolic networks. We found significant concentration differences in several amino acids, acylcarnitines and polyamines linking delirium-prone patients to known factors in Alzheimer's disease such as monoamine oxidase B (MAOB) protein. Subsequent computational structural comparison between MAOB and angiotensin converting enzyme 2 as well as protein-protein docking analysis showed that there potentially is strong binding of SARS-CoV-2 spike protein to MAOB. The possibility that SARS-CoV-2 influences MAOB activity leading to the observed neurological and platelet-based complications of SARS-CoV-2 infection requires further investigation.

Identification of transcriptional subtypes in lung adenocarcinoma and squamous cell carcinoma through integrative analysis of microarray and RNA sequencing data.

Classification of tumors into subtypes can inform personalized approaches to treatment including the choice of targeted therapies. The two most common lung cancer histological subtypes, lung adenocarcinoma and lung squamous cell carcinoma, have been previously divided into transcriptional subtypes using microarray data, and corresponding signatures were subsequently used to classify RNA-seq data. Cross-platform unsupervised classification facilitates the identification of robust transcriptional subtypes by combining vast amounts of publicly available microarray and RNA-seq data. However, cross-platform classification is challenging because of intrinsic differences in data generated using the two gene expression profiling technologies. In this report, we show that robust gene expression subtypes can be identified in integrated data representing over 3500 normal and tumor lung samples profiled using two widely used platforms, Affymetrix HG-U133 Plus 2.0 Array and Illumina HiSeq RNA sequencing. We tested and analyzed consensus clustering for 384 combinations of data processing methods. The agreement between subtypes identified in single-platform and cross-platform normalized data was then evaluated using a variety of statistics. Results show that unsupervised learning can be achieved with combined microarray and RNA-seq data using selected preprocessing, cross-platform normalization, and unsupervised feature selection methods. Our analysis confirmed three lung adenocarcinoma transcriptional subtypes, but only two consistent subtypes in squamous cell carcinoma, as opposed to four subtypes previously identified. Further analysis showed that tumor subtypes were associated with distinct patterns of genomic alterations in genes coding for therapeutic targets. Importantly, by integrating quantitative proteomics data, we were able to identify tumor subtype biomarkers that effectively classify samples on the basis of both gene and protein expression. This study provides the basis for further integrative data analysis across gene and protein expression profiling platforms.

Analysis and Simulation of Glioblastoma Cell Lines-Derived Extracellular Vesicles Metabolome.

Glioblastoma (GBM) is one of the most aggressive cancers of the central nervous system. Despite current advances in non-invasive imaging and the advent of novel therapeutic modalities, patient survival remains very low. There is a critical need for the development of effective biomarkers for GBM diagnosis and therapeutic monitoring. Extracellular vesicles (EVs) produced by GBM tumors have been shown to play an important role in cellular communication and modulation of the tumor microenvironment. As GBM-derived EVs contain specific "molecular signatures" of their parental cells and are able to transmigrate across the blood-brain barrier into biofluids such as the blood and cerebrospinal fluid (CSF), they are considered as a valuable source of potential diagnostic biomarkers. Given the relatively harsh extracellular environment of blood and CSF, EVs have to endure and adapt to different conditions. The ability of EVs to adjust and function depends on their lipid bilayer, metabolic content and enzymes and transport proteins. The knowledge of EVs metabolic characteristics and adaptability is essential for their utilization as diagnostic and therapeutic tools. The main aim of this study was to determine the metabolome of small EVs or exosomes derived from different GBM cells and compare to the metabolic profile of their parental cells using NMR spectroscopy. In addition, a possible flux of metabolic processes in GBM-derived EVs was simulated using constraint-based modeling from published proteomics information. Our results showed a clear difference between the metabolic profiles of GBM cells, EVs and media. Machine learning analysis of EV metabolomics, as well as flux simulation, supports the notion of active metabolism within EVs, including enzymatic reactions and the transfer of metabolites through the EV membrane. These results are discussed in the context of novel GBM diagnostics and therapeutic monitoring.

Mycobiome Sequencing and Analysis Applied to Fungal Community Profiling of the Lower Respiratory Tract During Fungal Pathogenesis.

Invasive fungal infections are an increasingly important cause of human morbidity and mortality. We generated a next-generation sequencing (NGS)-based method designed to detect a wide range of fungi and applied it to analysis of the fungal microbiome (mycobiome) of the lung during fungal infection. Internal transcribed spacer 1 (ITS1) amplicon sequencing and a custom analysis pipeline detected 96% of species from three mock communities comprised of potential fungal lung pathogens with good recapitulation of the expected species distributions (Pearson correlation coefficients r = 0.63, p = 0.004; r = 0.71, p < 0.001; r = 0.62, p = 0.002). We used this pipeline to analyze mycobiomes of bronchoalveolar lavage (BAL) specimens classified as culture-negative (n = 50) or culture-positive (n = 39) for Blastomyces dermatitidis/gilchristii, the causative agent of North America blastomycosis. Detected in 91.4% of the culture-positive samples, Blastomyces dominated (>50% relative abundance) the mycobiome in 68.6% of these culture-positive samples but was absent in culture-negative samples. To overcome any bias in relative abundance due to between-sample variation in fungal biomass, an abundance-weighting calculation was used to normalize the data by accounting for sample-specific PCR cycle number and PCR product concentration data utilized during sample preparation. After normalization, there was a statistically significant greater overall abundance of ITS1 amplicon in the Blastomyces-culture-positive samples versus culture-negative samples. Moreover, the normalization revealed a greater biomass of yeast and environmental fungi in several Blastomyces-culture-positive samples than in the culture-negative samples. Successful detection of Coccidioides, Scedosporium, Phaeoacremonium, and Aspergillus in 6 additional culture-positive BALs by ITS1 amplicon sequencing demonstrates the ability of this method to detect a broad range of fungi from clinical specimens, suggesting that it may be a potentially useful adjunct to traditional fungal microbiological testing for the diagnosis of respiratory mycoses.

Database of resistance related metabolites in Wheat Fusarium head blight Disease (MWFD).

Database URL: https://bioinfo.nrc.ca/mwfd/index.php.

NKT Cell-Deficient Mice Harbor an Altered Microbiota That Fuels Intestinal Inflammation during Chemically Induced Colitis.

NKT cells are unconventional T cells that respond to self and microbe-derived lipid and glycolipid Ags presented by the CD1d molecule. Invariant NKT (iNKT) cells influence immune responses in numerous diseases. Although only a few studies have examined their role during intestinal inflammation, it appears that iNKT cells protect from Th1-mediated inflammation but exacerbate Th2-mediated inflammation. Studies using iNKT cell-deficient mice and chemically induced dextran sodium sulfate (DSS) colitis have led to inconsistent results. In this study, we show that CD1d-deficient mice, which lack all NKT cells, harbor an altered intestinal microbiota that is associated with exacerbated intestinal inflammation at steady-state and following DSS treatment. This altered microbiota, characterized by increased abundance of the bacterial phyla Proteobacteria, Deferribacteres, and TM7, among which the mucin-eating Mucispirillum, as well as members of the genus Prevotella and segmented filamentous bacteria, was transmissible upon fecal transplant, along with the procolitogenic phenotype. Our results also demonstrate that this proinflammatory microbiota influences iNKT cell function upon activation during DSS colitis. Collectively, alterations of the microbiota have a major influence on colitis outcome and therefore have to be accounted for in such experimental settings and in studies focusing on iNKT cells.

Global Analysis of the Fungal Microbiome in Cystic Fibrosis Patients Reveals Loss of Function of the Transcriptional Repressor Nrg1 as a Mechanism of Pathogen Adaptation.

The microbiome shapes diverse facets of human biology and disease, with the importance of fungi only beginning to be appreciated. Microbial communities infiltrate diverse anatomical sites as with the respiratory tract of healthy humans and those with diseases such as cystic fibrosis, where chronic colonization and infection lead to clinical decline. Although fungi are frequently recovered from cystic fibrosis patient sputum samples and have been associated with deterioration of lung function, understanding of species and population dynamics remains in its infancy. Here, we coupled high-throughput sequencing of the ribosomal RNA internal transcribed spacer 1 (ITS1) with phenotypic and genotypic analyses of fungi from 89 sputum samples from 28 cystic fibrosis patients. Fungal communities defined by sequencing were concordant with those defined by culture-based analyses of 1,603 isolates from the same samples. Different patients harbored distinct fungal communities. There were detectable trends, however, including colonization with Candida and Aspergillus species, which was not perturbed by clinical exacerbation or treatment. We identified considerable inter- and intra-species phenotypic variation in traits important for host adaptation, including antifungal drug resistance and morphogenesis. While variation in drug resistance was largely between species, striking variation in morphogenesis emerged within Candida species. Filamentation was uncoupled from inducing cues in 28 Candida isolates recovered from six patients. The filamentous isolates were resistant to the filamentation-repressive effects of Pseudomonas aeruginosa, implicating inter-kingdom interactions as the selective force. Genome sequencing revealed that all but one of the filamentous isolates harbored mutations in the transcriptional repressor NRG1; such mutations were necessary and sufficient for the filamentous phenotype. Six independent nrg1 mutations arose in Candida isolates from different patients, providing a poignant example of parallel evolution. Together, this combined clinical-genomic approach provides a high-resolution portrait of the fungal microbiome of cystic fibrosis patient lungs and identifies a genetic basis of pathogen adaptation.

Regulation of obesity-related insulin resistance with gut anti-inflammatory agents.

Obesity has reached epidemic proportions, but little is known about its influence on the intestinal immune system. Here we show that the gut immune system is altered during high-fat diet (HFD) feeding and is a functional regulator of obesity-related insulin resistance (IR) that can be exploited therapeutically. Obesity induces a chronic phenotypic pro-inflammatory shift in bowel lamina propria immune cell populations. Reduction of the gut immune system, using beta7 integrin-deficient mice (Beta7(null)), decreases HFD-induced IR. Treatment of wild-type HFD C57BL/6 mice with the local gut anti-inflammatory, 5-aminosalicyclic acid (5-ASA), reverses bowel inflammation and improves metabolic parameters. These beneficial effects are dependent on adaptive and gut immunity and are associated with reduced gut permeability and endotoxemia, decreased visceral adipose tissue inflammation, and improved antigen-specific tolerance to luminal antigens. Thus, the mucosal immune system affects multiple pathways associated with systemic IR and represents a novel therapeutic target in this disease.

Mapping the cellular response to small molecules using chemogenomic fitness signatures.

Genome-wide characterization of the in vivo cellular response to perturbation is fundamental to understanding how cells survive stress. Identifying the proteins and pathways perturbed by small molecules affects biology and medicine by revealing the mechanisms of drug action. We used a yeast chemogenomics platform that quantifies the requirement for each gene for resistance to a compound in vivo to profile 3250 small molecules in a systematic and unbiased manner. We identified 317 compounds that specifically perturb the function of 121 genes and characterized the mechanism of specific compounds. Global analysis revealed that the cellular response to small molecules is limited and described by a network of 45 major chemogenomic signatures. Our results provide a resource for the discovery of functional interactions among genes, chemicals, and biological processes.

Miniature short hairpin RNA screens to characterize antiproliferative drugs.

The application of new proteomics and genomics technologies support a view in which few drugs act solely by inhibiting a single cellular target. Indeed, drug activity is modulated by complex, often incompletely understood cellular mechanisms. Therefore, efforts to decipher mode of action through genetic perturbation such as RNAi typically yields "hits" that fall into several categories. Of particular interest to the present study, we aimed to characterize secondary activities of drugs on cells. Inhibiting a known target can result in clinically relevant synthetic phenotypes. In one scenario, drug perturbation could, for example, improperly activate a protein that normally inhibits a particular kinase. In other cases, additional, lower affinity targets can be inhibited as in the example of inhibition of c-Kit observed in Bcr-Abl-positive cells treated with Gleevec. Drug transport and metabolism also play an important role in the way any chemicals act within the cells. Finally, RNAi per se can also affect cell fitness by more general off-target effects, e.g., via the modulation of apoptosis or DNA damage repair. Regardless of the root cause of these unwanted effects, understanding the scope of a drug's activity and polypharmacology is essential for better understanding its mechanism(s) of action, and such information can guide development of improved therapies. We describe a rapid, cost-effective approach to characterize primary and secondary effects of small-molecules by using small-scale libraries of virally integrated short hairpin RNAs. We demonstrate this principle using a "minipool" composed of shRNAs that target the genes encoding the reported protein targets of approved drugs. Among the 28 known reported drug-target pairs, we successfully identify 40% of the targets described in the literature and uncover several unanticipated drug-target interactions based on drug-induced synthetic lethality. We provide a detailed protocol for performing such screens and for analyzing the data. This cost-effective approach to mammalian knockdown screens, combined with the increasing maturation of RNAi technology will expand the accessibility of similar approaches in academic settings.

A phenotypic screening platform to identify small molecule modulators of Chlamydomonas reinhardtii growth, motility and photosynthesis.

Chemical biology, the interfacial discipline of using small molecules as probes to investigate biology, is a powerful approach of developing specific, rapidly acting tools that can be applied across organisms. The single-celled alga Chlamydomonas reinhardtii is an excellent model system because of its photosynthetic ability, cilia-related motility and simple genetics. We report the results of an automated fitness screen of 5,445 small molecules and subsequent assays on motility/phototaxis and photosynthesis. Cheminformatic analysis revealed active core structures and was used to construct a naïve Bayes model that successfully predicts algal bioactive compounds.

A comprehensive platform for highly multiplexed mammalian functional genetic screens.

BACKGROUND: Genome-wide screening in human and mouse cells using RNA interference and open reading frame over-expression libraries is rapidly becoming a viable experimental approach for many research labs. There are a variety of gene expression modulation libraries commercially available, however, detailed and validated protocols as well as the reagents necessary for deconvolving genome-scale gene screens using these libraries are lacking. As a solution, we designed a comprehensive platform for highly multiplexed functional genetic screens in human, mouse and yeast cells using popular, commercially available gene modulation libraries. The Gene Modulation Array Platform (GMAP) is a single microarray-based detection solution for deconvolution of loss and gain-of-function pooled screens. RESULTS: Experiments with specially constructed lentiviral-based plasmid pools containing ~78,000 shRNAs demonstrated that the GMAP is capable of deconvolving genome-wide shRNA "dropout" screens. Further experiments with a larger, ~90,000 shRNA pool demonstrate that equivalent results are obtained from plasmid pools and from genomic DNA derived from lentivirus infected cells. Parallel testing of large shRNA pools using GMAP and next-generation sequencing methods revealed that the two methods provide valid and complementary approaches to deconvolution of genome-wide shRNA screens. Additional experiments demonstrated that GMAP is equivalent to similar microarray-based products when used for deconvolution of open reading frame over-expression screens. CONCLUSION: Herein, we demonstrate four major applications for the GMAP resource, including deconvolution of pooled RNAi screens in cells with at least 90,000 distinct shRNAs. We also provide detailed methodologies for pooled shRNA screen readout using GMAP and compare next-generation sequencing to GMAP (i.e. microarray) based deconvolution methods.