Construction and optimization of multi-platform precision pathways for precision medicine.

In the enduring challenge against disease, advancements in medical technology have empowered clinicians with novel diagnostic platforms. Whilst in some cases, a single test may provide a confident diagnosis, often additional tests are required. However, to strike a balance between diagnostic accuracy and cost-effectiveness, one must rigorously construct the clinical pathways.

Viral Integration Plays a Minor Role in the Development and Prognostication of Oral Squamous Cell Carcinoma.

Viruses are well known drivers of several human malignancies. A causative factor for oral cavity squamous cell carcinoma (OSCC) in patients with limited exposure to traditional risk factors, including tobacco use, is yet to be identified. Our study aimed to comprehensively evaluate the role of viral drivers in OSCC patients with low cumulative exposure to traditional risk factors.

Whole genome analysis reveals the genomic complexity in metastatic cutaneous squamous cell carcinoma.

Metastatic cutaneous squamous cell carcinoma (CSCC) is a highly morbid disease requiring radical surgery and adjuvant therapy, which is associated with a poor prognosis. Yet, compared to other advanced malignancies, relatively little is known of the genomic landscape of metastatic CSCC. We have previously reported the mutational signatures and mutational patterns of CCCTC-binding factor (CTCF) regions in metastatic CSCC.

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.

Whole genome duplication in oral squamous cell carcinoma in patients younger than 50 years: implications for prognosis and adverse clinicopathological factors.

Introduction: Oral squamous cell carcinoma (OSCC) in the young (<50 years), without known carcinogenic risk factors, is on the rise globally. Whole genome duplication (WGD) has been shown to occur at higher rates in cancers without an identifiable carcinogenic agent. We aimed to evaluate the prevalence of WGD in a cohort of OSCC patients under the age of 50 years.

Oral Squamous Cell Carcinoma in Young Patients Show Higher Rates of Amplification: Implications for Novel Personalized Therapy.

There is an increasing worldwide incidence of patients under 50 years of age presenting with oral squamous cell carcinoma (OSCC). The molecular mechanisms driving disease in this emerging cohort remain unclear, limiting impactful treatment options for these patients. To identify common clinically actionable targets in this cohort, we used whole genome and transcriptomic sequencing of OSCC patient samples from 26 individuals under 50 years of age.

Genome-wide analysis in Drosophila reveals diet-by-gene interactions and uncovers diet-responsive genes.

Genetic and environmental factors play a major role in metabolic health. However, they do not act in isolation, as a change in an environmental factor such as diet may exert different effects based on an individual's genotype. Here, we sought to understand how such gene-diet interactions influenced nutrient storage and utilization, a major determinant of metabolic disease.

Transcriptomic analysis of Nodal - and BMP- associated genes during development to the juvenile seastar in Parvulastra exigua (Asterinidae).

The molecular mechanisms underlying development of the pentameral body of adult echinoderms are poorly understood but are important to solve with respect to evolution of a unique body plan that contrasts with the bilateral body plan of other deuterostomes. As Nodal and BMP2/4 signalling is involved in axis formation in larvae and development of the echinoderm body plan, we used the developmental transcriptome generated for the asterinid seastar Parvulastra exigua to investigate the temporal expression patterns of Nodal and BMP2/4 genes from the embryo and across metamorphosis to the juvenile. For echinoderms, the Asteroidea represents the basal-type body architecture with a distinct (separated) ray structure.

Transcriptomic analysis of sea star development through metamorphosis to the highly derived pentameral body plan with a focus on neural transcription factors.

The Echinodermata is characterized by a secondarily evolved pentameral body plan. While the evolutionary origin of this body plan has been the subject of debate, the molecular mechanisms underlying its development are poorly understood. We assembled a de novo developmental transcriptome from the embryo through metamorphosis in the sea star Parvulastra exigua.

Transcriptional downregulation of MHC class I and melanoma de- differentiation in resistance to PD-1 inhibition.

Transcriptomic signatures designed to predict melanoma patient responses to PD-1 blockade have been reported but rarely validated. We now show that intra-patient heterogeneity of tumor responses to PD-1 inhibition limit the predictive performance of these signatures. We reasoned that resistance mechanisms will reflect the tumor microenvironment, and thus we examined PD-1 inhibitor resistance relative to T-cell activity in 94 melanoma tumors collected at baseline and at time of PD-1 inhibitor progression.

Analysis of the Whole-Exome Sequencing of Tumor and Circulating Tumor DNA in Metastatic Melanoma.

The use of circulating tumor DNA (ctDNA) to monitor cancer progression and response to therapy has significant potential but there is only limited data on whether this technique can detect the presence of low frequency subclones that may ultimately confer therapy resistance. In this study, we sought to evaluate whether whole-exome sequencing (WES) of ctDNA could accurately profile the mutation landscape of metastatic melanoma. We used WES to identify variants in matched, tumor-derived genomic DNA (gDNA) and plasma-derived ctDNA isolated from a cohort of 10 metastatic cutaneous melanoma patients.

Evaluating stably expressed genes in single cells.

Background: Single-cell RNA-seq (scRNA-seq) profiling has revealed remarkable variation in transcription, suggesting that expression of many genes at the single-cell level is intrinsically stochastic and noisy. Yet, on the cell population level, a subset of genes traditionally referred to as housekeeping genes (HKGs) are found to be stably expressed in different cell and tissue types. It is therefore critical to question whether stably expressed genes (SEGs) can be identified on the single-cell level, and if so, how can their expression stability be assessed? We have previously proposed a computational framework for ranking expression stability of genes in single cells for scRNA-seq data normalization and integration.

Melanoma Explorer: a web application to allow easy reanalysis of publicly available and clinically annotated melanoma omics data sets.

Validating newly discovered biomarkers in large, publicly available data sets is often difficult and requires specialized computer programming skills. Melanoma Explorer is a web application that enables easy interrogation of melanoma omics data sets that are freely available in online data repositories with a point-and-click interface. Two use cases are demonstrated.

DCARS: differential correlation across ranked samples.

Motivation: Genes act as a system and not in isolation. Thus, it is important to consider coordinated changes of gene expression rather than single genes when investigating biological phenomena such as the aetiology of cancer. We have developed an approach for quantifying how changes in the association between pairs of genes may inform the outcome of interest called Differential Correlation across Ranked Samples (DCARS).

Quantitative Performance Evaluator for Proteomics (QPEP): Web-based Application for Reproducible Evaluation of Proteomics Preprocessing Methods.

Tandem mass spectrometry is one of the most popular techniques for quantitation of proteomes. There exists a large variety of options in each stage of data preprocessing that impact the bias and variance of the summarized protein-level values. Using a newly released data set satisfying a replicated Latin squares design, a diverse set of performance metrics has been developed and implemented in a web-based application, Quantitative Performance Evaluator for Proteomics (QPEP).

Nodal and BMP expression during the transition to pentamery in the sea urchin Heliocidaris erythrogramma: insights into patterning the enigmatic echinoderm body plan.

Background: The molecular mechanisms underlying the development of the unusual echinoderm pentameral body plan and their likeness to mechanisms underlying the development of the bilateral plans of other deuterostomes are of interest in tracing body plan evolution. In this first study of the spatial expression of genes associated with Nodal and BMP2/4 signalling during the transition to pentamery in sea urchins, we investigate Heliocidaris erythrogramma, a species that provides access to the developing adult rudiment within days of fertilization.

Results: BMP2/4, and the putative downstream genes, Six1/2, Eya, Tbx2/3 and Msx were expressed in the earliest morphological manifestation of pentamery during development, the five hydrocoele lobes.

Differential distribution improves gene selection stability and has competitive classification performance for patient survival.

A consistent difference in average expression level, often referred to as differential expression (DE), has long been used to identify genes useful for classification. However, recent cancer studies have shown that when transcription factors or epigenetic signals become deregulated, a change in expression variability (DV) of target genes is frequently observed. This suggests that assessing the importance of genes by either differential expression or variability alone potentially misses sets of important biomarkers that could lead to improved predictions and treatments.

PD-L1 Negative Status is Associated with Lower Mutation Burden, Differential Expression of Immune-Related Genes, and Worse Survival in Stage III Melanoma.

Purpose: Understanding why some melanomas test negative for PD-L1 by IHC may have implications for the application of anti-PD-1 therapies in melanoma management. This study sought to determine somatic mutation and gene expression patterns associated with tumor cell PD-L1 expression, or lack thereof, in stage III metastatic melanoma to better define therapeutically relevant patient subgroups.

Experimental Design: IHC for PD-L1 was assessed in 52 American Joint Committee on Cancer stage III melanoma lymph node specimens and compared with specimen-matched comprehensive clinicopathologic, genomic, and transcriptomic data.

Targeting activating mutations of EZH2 leads to potent cell growth inhibition in human melanoma by derepression of tumor suppressor genes.

The epigenetic modifier EZH2 is part of the polycomb repressive complex that suppresses gene expression via histone methylation. Activating mutations in EZH2 are found in a subset of melanoma that contributes to disease progression by inactivating tumor suppressor genes. In this study we have targeted EZH2 with a specific inhibitor (GSK126) or depleted EZH2 protein by stable shRNA knockdown.

Combining BET and HDAC inhibitors synergistically induces apoptosis of melanoma and suppresses AKT and YAP signaling.

Histone acetylation marks have an important role in controlling gene expression and are removed by histone deacetylases (HDACs). These marks are read by bromodomain and extra-terminal (BET) proteins and novel inhibitiors of these proteins are currently in clinical development. Inhibitors of HDAC and BET proteins have individually been shown to cause apoptosis and reduce growth of melanoma cells.

Transcriptomic analysis of Nodal- and BMP-associated genes during juvenile development of the sea urchin Heliocidaris erythrogramma.

Understanding the unusual radial body plan of echinoderms and its relationship to the bilateral plan of other deuterostomes remains a challenge. The molecular processes of embryonic and early larval development in sea urchins are well characterised, but those giving rise to the adult and its radial body remain poorly studied. We used the developmental transcriptome generated for Heliocidaris erythrogramma, a species that forms the juvenile soon after gastrulation, to investigate changes in gene expression underlying radial body development.

ClassifyR: an R package for performance assessment of classification with applications to transcriptomics.

Unlabelled: Although a large collection of classification software packages exist in R, a new generic framework for linking custom classification functions with classification performance measures is needed. A generic classification framework has been designed and implemented as an R package in an object oriented style. Its design places emphasis on parallel processing, reproducibility and extensibility.

Wnt inhibitory factor 1 (WIF1) is a marker of osteoblastic differentiation stage and is not silenced by DNA methylation in osteosarcoma.

Wnt pathway targeting is of high clinical interest for treating bone loss disorders such as osteoporosis. These therapies inhibit the action of negative regulators of osteoblastic Wnt signaling. The report that Wnt inhibitory factor 1 (WIF1) was epigenetically silenced via promoter DNA methylation in osteosarcoma (OS) raised potential concerns for such treatment approaches.

Detection and classification of peaks in 5' cap RNA sequencing data.

Background: The large-scale sequencing of 5' cap enriched cDNA promises to reveal the diversity of transcription initiation across entire genomes. The process of transcription is noisy, and there is often no single, exact start site. This creates the need for a fast and simple method of identifying transcription start peaks based on this type of data.