You Guo

Differential expression analysis for individual cancer samples based on robust within-sample relative gene expression orderings across multiple profiling platforms

The highly stable within-sample relative expression orderings (REOs) of gene pairs in a particular type of human normal tissue are widely reversed in the cancer condition. Based on this finding, we have recently proposed an algorithm named RankComp to detect differentially expressed genes (DEGs) for individual disease samples measured by a particular platform. In this paper, with 461 normal lung tissue samples separately measured by four commonly used platforms, we demonstrated that tens of millions of gene pairs with significantly stable REOs in normal lung tissue can be consistently detected in samples measured by different platforms. However, about 20% of stable REOs commonly detected by two different platforms (e.g., Affymetrix and Illumina platforms) showed inconsistent REO patterns due to the differences in probe design principles. Based on the significantly stable REOs (FDR<0.01) for normal lung tissue consistently detected by the four platforms, which tended to have large rank differences, RankComp detected averagely 1184, 1335 and 1116 DEGs per sample with averagely 96.51%, 95.95% and 94.78% precisions in three evaluation datasets with 25, 57 and 58 paired lung cancer and normal samples, respectively. Individualized pathway analysis revealed some common and subtype-specific functional mechanisms of lung cancer. Similar results were observed for colorectal cancer. In conclusion, based on the cross-platform significantly stable REOs for a particular normal tissue, differentially expressed genes and pathways in any disease sample measured by any of the platforms can be readily and accurately detected, which could be further exploited for dissecting the heterogeneity of cancer.

Common DNA methylation alterations of Alzheimer’s disease and aging in peripheral whole blood

Alzheimers disease (AD) is a common aging-related neurodegenerative illness. Recently, many studies have tried to identify AD- or aging-related DNA methylation (DNAm) biomarkers from peripheral whole blood (PWB). However, the origin of PWB biomarkers is still controversial. In this study, by analyzing 2565 DNAm profiles for PWB and brain tissue, we showed that aging-related DNAm CpGs (Age-CpGs) and AD-related DNAm CpGs (AD-CpGs) observable in PWB both mainly reflected DNAm alterations intrinsic in leukocyte subtypes rather than methylation differences introduced by the increased ratio of myeloid to lymphoid cells during aging or AD progression. The PWB Age-CpGs and AD-CpGs significantly overlapped 107 sites (P-value = 2.61×10−12) and 97 had significantly concordant methylation alterations in AD and aging (P-value < 2.2×10−16), which were significantly enriched in nervous system development, neuron differentiation and neurogenesis. More than 60.8% of these 97 concordant sites were found to be significantly correlated with age in normal peripheral CD4+ T cells and CD14+ monocytes as well as in four brain regions, and 44 sites were also significantly differentially methylated in different regions of AD brain tissue. Taken together, the PWB DNAm alterations related to both aging and AD could be exploited for identification of AD biomarkers.

An individualized prognostic signature and multi‑omics distinction for early stage hepatocellular carcinoma patients with surgical resection

Previously reported prognostic signatures for predicting the prognoses of postsurgical hepatocellular carcinoma (HCC) patients are commonly based on predefined risk scores, which are hardly applicable to samples measured by different laboratories. To solve this problem, using gene expression profiles of 170 stage I/II HCC samples, we identified a prognostic signature consisting of 20 gene pairs whose within-sample relative expression orderings (REOs) could robustly predict the disease-free survival and overall survival of HCC patients. This REOs-based prognostic signature was validated in two independent datasets. Functional enrichment analysis showed that the patients with high-risk of recurrence were characterized by the activations of pathways related to cell proliferation and tumor microenvironment, whereas the low-risk patients were characterized by the activations of various metabolism pathways. We further investigated the distinct epigenomic and genomic characteristics of the two prognostic groups using The Cancer Genome Atlas samples with multi-omics data. Epigenetic analysis showed that the transcriptional differences between the two prognostic groups were significantly concordant with DNA methylation alternations. The signaling network analysis identified several key genes (e.g. TP53, MYC) with epigenomic or genomic alternations driving poor prognoses of HCC patients. These results help us understand the multi-omics mechanisms determining the outcomes of HCC patients.

Circumvent the uncertainty in the applications of transcriptional signatures to tumor tissues sampled from different tumor sites

The expression measurements of thousands of genes are correlated with the proportions of tumor epithelial cell (PTEC) in clinical samples. Thus, for a tumor diagnostic or prognostic signature based on a summarization of expression levels of the signature genes, the risk score for a patient may dependent on the tumor tissues sampled from different tumor sites with diverse PTEC for the same patient. Here, we proposed that the within-samples relative expression orderings (REOs) based gene pairs signatures should be insensitive to PTEC variations. Firstly, by analysis of paired tumor epithelial cell and stromal cell microdissected samples from 27 cancer patients, we showed that above 80% of gene pairs had consistent REOs between the two cells, indicating these REOs would be independent of PTEC in cancer tissues. Then, by simulating tumor tissues with different PTEC using each of the 27 paired samples, we showed that about 90% REOs of gene pairs in tumor epithelial cells were maintained in tumor samples even when PTEC decreased to 30%. Especially, the REOs of gene pairs with larger expression differences in tumor epithelial cells tend to be more robust against PTEC variations. Finally, as a case study, we developed a gene pair signature which could robustly distinguish colorectal cancer tissues with various PTEC from normal tissues. We concluded that the REOs-based signatures were robust against PTEC variations.

Identifying clinically relevant drug resistance genes in drug-induced resistant cancer cell lines and post-chemotherapy tissues.

Until recently, few molecular signatures of drug resistance identified in drug-induced resistant cancer cell models can be translated into clinical practice. Here, we defined differentially expressed genes (DEGs) between pre-chemotherapy colorectal cancer (CRC) tissue samples of non-responders and responders for 5-fluorouracil and oxaliplatin-based therapy as clinically relevant drug resistance genes (CRG5-FU/L-OHP). Taking CRG5-FU/L-OHP as reference, we evaluated the clinical relevance of several types of genes derived from HCT116 CRC cells with resistance to 5-fluorouracil and oxaliplatin, respectively. The results revealed that DEGs between parental and resistant cells, when both were treated with the corresponding drug for a certain time, were significantly consistent with the CRG5-FU/L-OHP as well as the DEGs between the post-chemotherapy CRC specimens of responders and non-responders. This study suggests a novel strategy to extract clinically relevant drug resistance genes from both drug-induced resistant cell models and post-chemotherapy cancer tissue specimens.

Discriminating cancer-related and cancer-unrelated chemoradiation-response genes for locally advanced rectal cancers

For patients with locally advanced rectal cancer (LARC) treated with preoperation chemoradiation (pCRT), identifying differentially expressed (DE) genes between non-responders and responders is a common approach for investigating mechanisms of chemoradiation resistance. However, some of such DE genes might be irrelevant to cancer itself but simply reflect the pharmacokinetic differences of the normal tissues. In this study, we adopted the RankComp algorithm to identify DE genes for each of LARC sample compared with its own normal state. Then, we identified genes with significantly different deregulation frequencies between the non-responders and responders, defined as cancer-related pCRT-response genes. Pathway enrichment and protein-protein interaction analyses showed that these genes specifically and intensively interacted with currently known effective genes of pCRT, involving in DNA replication, cell cycle and DNA repair. In contrast, after excluding the cancer-related pCRT-response genes, the other DE genes between non-responders and responders were enriched in many pathways of drug and protein metabolisms and transports, and interacted with both the known effective genes and pharmacokinetic genes. Hence, these two types of DE genes should be distinguished for investigating mechanisms of pCRT response in LARCs.

A qualitative signature for predicting pathological response to neoadjuvant chemoradiation in locally advanced rectal cancers

BACKGROUND AND PURPOSE The standard therapy for locally advanced rectal cancers (LARCs) is neoadjuvant chemoradiation (nCRT) followed by surgical resection. Pathological response to nCRT varies among patients, and it remains a challenge to predict pathological response to nCRT in LARCs. MATERIAL AND METHODS Using 42 samples as the training cohort, we searched a signature by screening the gene pairs whose within-sample relative expression orderings are significantly correlated with the pathological response. The signature was validated in both a public cohort of 46 samples and a cohort of 33 samples measured at our laboratory. RESULTS A signature consisting of 27 gene pairs was identified in the training cohort with an accuracy of 92.86% and an area under the receiver operating characteristic curve (AUC) of 0.95. The accuracy was 89.13% for the public test cohort and 90.91% for the private test cohort, with AUC being 0.95 and 0.91, respectively. Furthermore, the signature was used to predict disease-free survival benefits from 5Fu-based chemotherapy in 285 locally advanced colorectal cancers. CONCLUSIONS The signature consisting of 27 gene pairs can robustly predict clinical response of LARCs to nCRT.

Quantitative or qualitative transcriptional diagnostic signatures? A case study for colorectal cancer

BackgroundDue to experimental batch effects, the application of a quantitative transcriptional signature for disease diagnoses commonly requires inter-sample data normalization, which would be hardly applicable under common clinical settings. Many cancers might have qualitative differences with the non-cancer states in the gene expression pattern. Therefore, it is reasonable to explore the power of qualitative diagnostic signatures which are robust against experimental batch effects and other random factors.ResultsFirstly, using data of technical replicate samples from the MicroArray Quality Control (MAQC) project, we demonstrated that the low-throughput PCR-based technologies also exist large measurement variations for gene expression even when the samples were measured in the same test site. Then, we demonstrated the critical limitation of low stability for classifiers based on quantitative transcriptional signatures in applications to individual samples through a case study using a support vector machine and a naïve Bayesian classifier to discriminate colorectal cancer tissues from normal tissues. To address this problem, we identified a signature consisting of three gene pairs for discriminating colorectal cancer tissues from non-cancer (normal and inflammatory bowel disease) tissues based on within-sample relative expression orderings (REOs) of these gene pairs. The signature was well verified using 22 independent datasets measured by different microarray and RNA_seq platforms, obviating the need of inter-sample data normalization.ConclusionsSubtle quantitative information of gene expression measurements tends to be unstable under current technical conditions, which will introduce uncertainty to clinical applications of the quantitative transcriptional diagnostic signatures. For diagnosis of disease states with qualitative transcriptional characteristics, the qualitative REO-based signatures could be robustly applied to individual samples measured by different platforms.

Evaluating hepatocellular carcinoma cell lines for tumor samples using within-sample relative expression orderings of genes

Concerns are raised about the representativeness of cell lines for tumours due to the culture environment and misidentification. Liver is a major metastatic destination of many cancers, which might further confuse the origin of hepatocellular carcinoma cell lines. Therefore, it is of crucial importance to understand how well they can represent hepatocellular carcinoma.

A qualitative signature for early diagnosis of hepatocellular carcinoma based on relative expression orderings

Currently, using biopsy specimens to confirm suspicious liver lesions of early hepatocellular carcinoma are not entirely reliable because of insufficient sampling amount and inaccurate sampling location. It is necessary to develop a signature to aid early hepatocellular carcinoma diagnosis using biopsy specimens even when the sampling location is inaccurate.