scMerge leverages factor analysis, stable expression, and pseudoreplication to merge multiple single-cell RNA-seq datasets.

Concerted examination of multiple collections of single-cell RNA sequencing (RNA-seq) data promises further biological insights that cannot be uncovered with individual datasets. Here we present scMerge, an algorithm that integrates multiple single-cell RNA-seq datasets using factor analysis of stably expressed genes and pseudoreplicates across datasets. Using a large collection of public datasets, we benchmark scMerge against published methods and demonstrate that it consistently provides improved cell type separation by removing unwanted factors; scMerge can also enhance biological discovery through robust data integration, which we show through the inference of development trajectory in a liver dataset collection.

bcGST-an interactive bias-correction method to identify over-represented gene-sets in boutique arrays.

MOTIVATION: Gene annotation and pathway databases such as Gene Ontology and Kyoto Encyclopaedia of Genes and Genomes are important tools in Gene-Set Test (GST) that describe gene biological functions and associated pathways. GST aims to establish an association relationship between a gene-set of interest and an annotation. Importantly, GST tests for over-representation of genes in an annotation term. One implicit assumption of GST is that the gene expression platform captures the complete or a very large proportion of the genome. However, this assumption is neither satisfied for the increasingly popular boutique array nor the custom designed gene expression profiling platform. Specifically, conventional GST is no longer appropriate due to the gene-set selection bias induced during the construction of these platforms. RESULTS: We propose bcGST, a bias-corrected GST by introducing bias-correction terms in the contingency table needed for calculating the Fisher's Exact Test. The adjustment method works by estimating the proportion of genes captured on the array with respect to the genome in order to assist filtration of annotation terms that would otherwise be falsely included or excluded. We illustrate the practicality of bcGST and its stability through multiple differential gene expression analyses in melanoma and the Cancer Genome Atlas cancer studies. AVAILABILITY AND IMPLEMENTATION: The bcGST method is made available as a Shiny web application at http://shiny.maths.usyd.edu.au/bcGST/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Cross-Platform Omics Prediction procedure: a statistical machine learning framework for wider implementation of precision medicine.

In this modern era of precision medicine, molecular signatures identified from advanced omics technologies hold great promise to better guide clinical decisions. However, current approaches are often location-specific due to the inherent differences between platforms and across multiple centres, thus limiting the transferability of molecular signatures. We present Cross-Platform Omics Prediction (CPOP), a penalised regression model that can use omics data to predict patient outcomes in a platform-independent manner and across time and experiments. CPOP improves on the traditional prediction framework of using gene-based features by selecting ratio-based features with similar estimated effect sizes. These components gave CPOP the ability to have a stable performance across datasets of similar biology, minimising the effect of technical noise often generated by omics platforms. We present a comprehensive evaluation using melanoma transcriptomics data to demonstrate its potential to be used as a critical part of a clinical screening framework for precision medicine. Additional assessment of generalisation was demonstrated with ovarian cancer and inflammatory bowel disease studies.

Distinct Molecular Profiles and Immunotherapy Treatment Outcomes of V600E and V600K BRAF-Mutant Melanoma.

PURPOSE: BRAF V600E and V600K melanomas have distinct clinicopathologic features, and V600K appear to be less responsive to BRAFi±MEKi. We investigated mechanisms for this and explored whether genotype affects response to immunotherapy. EXPERIMENTAL DESIGN: Pretreatment formalin-fixed paraffin-embedded tumors from patients treated with BRAFi±MEKi underwent gene expression profiling and DNA sequencing. Molecular results were validated using The Cancer Genome Atlas (TCGA) data. An independent cohort of V600E/K patients treated with anti-PD-1 immunotherapy was examined. RESULTS: Baseline tissue and clinical outcome with BRAFi±MEKi were studied in 93 patients (78 V600E, 15 V600K). V600K patients had numerically less tumor regression (median, -31% vs. -52%, P = 0.154) and shorter progression-free survival (PFS; median, 5.7 vs. 7.1 months, P = 0.15) compared with V600E. V600K melanomas had lower expression of the ERK pathway feedback regulator dual-specificity phosphatase 6, confirmed with TCGA data (116 V600E, 17 V600K). Pathway analysis showed V600K had lower expression of ERK and higher expression of PI3K-AKT genes than V600E. Higher mutational load was observed in V600K, with a higher proportion of mutations in PIK3R1 and tumor-suppressor genes. In patients treated with anti-PD-1, V600K (n = 19) had superior outcomes than V600E (n = 84), including response rate (53% vs. 29%, P = 0.059), PFS (median, 19 vs. 2.7 months, P = 0.049), and overall survival (20.4 vs. 11.7 months, P = 0.081). CONCLUSIONS: BRAF V600K melanomas appear to benefit less from BRAFi±MEKi than V600E, potentially due to less reliance on ERK pathway activation and greater use of alternative pathways. In contrast, these melanomas have higher mutational load and respond better to immunotherapy.