Quality control (QC) becomes more important in pre-processing analysis of high dimensional omics data. Several routine QC processes became a standard process in omics data analysis. The standard QC analysis includes calculating quality-related measures, checking the consistency among samples, detecting outlying observations and so forth. QC analysis tends to be more important in the era of high dimensional omics data. Although several QC analysis tools providing simple graphical display have been developed by many researchers, they usually require a subjective decision on QC. Here, we propose high-dimensional data quality control (HidQC) plot which is a simple and efficient QC tool for handling high dimensional omics data. HidQC plot primarily focuses on identifying samples of poor quality by conducting a contrast analysis for the between/within group distances and the summary distances. HidQC plot checks the quality by investigating the consistency of samples for each group. Unlike other QC plots, HidQC plot provides the p-value of each sample based on the permutation test, which can be used as a more objective criterion to determine whether to use the sample or not. We applied HidQC plot to MicroArray Quality Control (MAQC) project 1 data to demonstrate its usefulness.

Since most of the cancer markers that have been reported are obtained directly from cancer tissues, it is difficult to use them for early diagnosis of cancer without surgery. Thus, development of markers that can be detected by blood is crucial for making early diagnosis of cancer easier. One of the most feasible types of markers that can be detected by blood is a protein marker. Here, we focus on building prediction methods using the protein markers for early diagnosis of cancer. To develop a prediction model with high prediction ability, it is critical to choose appropriate markers first. Here, we consider a stepwise selection method using area under the receiver operating characteristic curve (Step-AUC) in order to construct a multi-protein prediction model. We showed that the performance of Step-AUC highly depends on the tuning parameter. We compared our proposed Step-AUC method to stepwise selection using information criteria and support vector machine recursive feature extraction (SVM-RFE). We observed that Step-AUC and stepwise selection using Bayesian information criteria (Step-BIC) perform better than other methods. The importance of each marker can be chosen using a new stepwise selection consistency (SSC) measure. The final models include the markers with high SSC measures. We applied our stepwise procedure to pancreatic cancer data and found two markers of interest.

An important issue in genetic association studies is gene-gene interaction (epistasis) underlying common complex diseases that are affected by multiple genetic factors. Although many methods have been proposed to analyze epistasis, the interpretation of the identified gene-gene interactions is not straightforward. In order to efficiently provide the statistical interpretation and biological evidences of gene-gene interactions, we developed the VizEpis, a tool for visualizing of gene-gene interactions in genetic association analysis and mapping of epistatic interaction to the biological evidence from public interaction databases. Using interaction network and circular plot, the VizEpis provides to explore the interaction network integrated with biological evidences in epigenetic regulation, splicing, transcription, translation and post-translation level. To aid statistical interaction in genotype level, the VizEpis provides pairwise checkerboard, forest, and funnel. VizEpis provides for the user-specified variants with statistical interaction to find biological evidences.

Pancreatic ductal adenocarcinoma (PDAC) is the most common type of pancreatic cancer, which has dismal prognosis because of its silent early symptoms, high metastatic potential, and resistance to conventional therapies. Although a PDAC patient who is diagnosed at an early stage would have a substantial increase in chance of survival, the survival rate is poor because there is no efficient non-invasive diagnostic test in the early stage. In this study, we developed an efficient prediction models to detect PDAC in its early stages. Our prediction models use both mRNA and miRNA expression data from 104 PDAC tissues and 17 normal pancreatic tissues using microarray technology. After quality control, we built prediction models based on support vector machine (SVM) from mRNA and miRNA expressions for detecting early PDAC. To prevent over-fitting effect, we conducted leave-one-out cross validation (LOOCV) and 5-fold cross validation (CV). For independent validation of prediction models, we performed evaluation on independent datasets from Gene Expression Omnibus (GEO). After the validation, we identified 28 single markers and 231 combinations of markers with powerful prediction performance. In addition, the marker candidates are annotated with cancer pathways using gene ontology analysis. Our prediction models for PDAC may have potential for early diagnosis of PDAC.

MOTIVATION: For the past few decades, many statistical methods in genome-wide association studies (GWAS) have been developed to identify SNP-SNP interactions for case-control studies. However, there has been less work for prospective cohort studies, involving the survival time. Recently, Gui et al. (2011) proposed a novel method, called Surv-MDR, for detecting gene-gene interactions associated with survival time. Surv-MDR is an extension of the multifactor dimensionality reduction (MDR) method to the survival phenotype by using the log-rank test for defining a binary attribute. However, the Surv-MDR method has some drawbacks in the sense that it needs more intensive computations and does not allow for a covariate adjustment. In this article, we propose a new approach, called Cox-MDR, which is an extension of the generalized multifactor dimensionality reduction (GMDR) to the survival phenotype by using a martingale residual as a score to classify multi-level genotypes as high- and low-risk groups. The advantages of Cox-MDR over Surv-MDR are to allow for the effects of discrete and quantitative covariates in the frame of Cox regression model and to require less computation than Surv-MDR. RESULTS: Through simulation studies, we compared the power of Cox-MDR with those of Surv-MDR and Cox regression model for various heritability and minor allele frequency combinations without and with adjusting for covariate. We found that Cox-MDR and Cox regression model perform better than Surv-MDR for low minor allele frequency of 0.2, but Surv-MDR has high power for minor allele frequency of 0.4. However, when the effect of covariate is adjusted for, Cox-MDR and Cox regression model perform much better than Surv-MDR. We also compared the performance of Cox-MDR and Surv-MDR for a real data of leukemia patients to detect the gene-gene interactions with the survival time. CONTACT: [email protected]; [email protected].

BACKGROUND: Because common complex diseases are affected by multiple genes and environmental factors, it is essential to investigate gene-gene and/or gene-environment interactions to understand genetic architecture of complex diseases. After the great success of large scale genome-wide association (GWA) studies using the high density single nucleotide polymorphism (SNP) chips, the study of gene-gene interaction becomes a next challenge. Multifactor dimensionality reduction (MDR) analysis has been widely used for the gene-gene interaction analysis. In practice, however, it is not easy to perform high order gene-gene interaction analyses via MDR in genome-wide level because it requires exploring a huge search space and suffers from a computational burden due to high dimensionality. RESULTS: We propose dimensional reduction analysis, Gene-MDR analysis for the fast and efficient high order gene-gene interaction analysis. The proposed Gene-MDR method is composed of two-step applications of MDR: within- and between-gene MDR analyses. First, within-gene MDR analysis summarizes each gene effect via MDR analysis by combining multiple SNPs from the same gene. Second, between-gene MDR analysis then performs interaction analysis using the summarized gene effects from within-gene MDR analysis. We apply the Gene-MDR method to bipolar disorder (BD) GWA data from Wellcome Trust Case Control Consortium (WTCCC). The results demonstrate that Gene-MDR is capable of detecting high order gene-gene interactions associated with BD. CONCLUSION: By reducing the dimension of genome-wide data from SNP level to gene level, Gene-MDR efficiently identifies high order gene-gene interactions. Therefore, Gene-MDR can provide the key to understand complex disease etiology.

Most common complex traits are affected by multiple genes and/or environmental factors. To understand genetic architecture of complex traits, the investigation of gene-gene and gene-environment interactions can be essential. However, conducting gene-gene interaction using genome-wide data requires exploring a huge search space and suffers from a computation burden due to high dimensionality of genetic data. To identify gene-gene interaction more efficiently, we propose a gene-based reduction method which first summarizes the gene effect by combining multiple single nucleotide polymorphism (SNP) and then performs the gene-gene interaction via the summarized gene effect. By reducing the search space from SNPs to gene, our gene-based method becomes efficient and fast for identifying gene-gene interaction in genome wide association studies. The gene-based reduction method is illustrated by hypertension data from a Korean population.

BACKGROUND: Quantification of protein expression by means of mass spectrometry (MS) has been introduced in various proteomics studies. In particular, two label-free quantification methods, such as spectral counting and spectra feature analysis have been extensively investigated in a wide variety of proteomic studies. The cornerstone of both methods is peptide identification based on a proteomic database search and subsequent estimation of peptide retention time. However, they often suffer from restrictive database search and inaccurate estimation of the liquid chromatography (LC) retention time. Furthermore, conventional peptide identification methods based on the spectral library search algorithms such as SEQUEST or SpectraST have been found to provide neither the best match nor high-scored matches. Lastly, these methods are limited in the sense that target peptides cannot be identified unless they have been previously generated and stored into the database or spectral libraries.To overcome these limitations, we propose a novel method, namely Quantification method based on Finding the Identical Spectral set for a Homogenous peptide (Q-FISH) to estimate the peptide's abundance from its tandem mass spectrometry (MS/MS) spectra through the direct comparison of experimental spectra. Intuitively, our Q-FISH method compares all possible pairs of experimental spectra in order to identify both known and novel proteins, significantly enhancing identification accuracy by grouping replicated spectra from the same peptide targets. RESULTS: We applied Q-FISH to Nano-LC-MS/MS data obtained from human hepatocellular carcinoma (HCC) and normal liver tissue samples to identify differentially expressed peptides between the normal and disease samples. For a total of 44,318 spectra obtained through MS/MS analysis, Q-FISH yielded 14,747 clusters. Among these, 5,777 clusters were identified only in the HCC sample, 6,648 clusters only in the normal tissue sample, and 2,323 clusters both in the HCC and normal tissue samples. While it will be interesting to investigate peptide clusters only found from one sample, further examined spectral clusters identified both in the HCC and normal samples since our goal is to identify and assess differentially expressed peptides quantitatively. The next step was to perform a beta-binomial test to isolate differentially expressed peptides between the HCC and normal tissue samples. This test resulted in 84 peptides with significantly differential spectral counts between the HCC and normal tissue samples. We independently identified 50 and 95 peptides by SEQUEST, of which 24 and 56 peptides, respectively, were found to be known biomarkers for the human liver cancer. Comparing Q-FISH and SEQUEST results, we found 22 of the differentially expressed 84 peptides by Q-FISH were also identified by SEQUEST. Remarkably, of these 22 peptides discovered both by Q-FISH and SEQUEST, 13 peptides are known for human liver cancer and the remaining 9 peptides are known to be associated with other cancers. CONCLUSIONS: We proposed a novel statistical method, Q-FISH, for accurately identifying protein species and simultaneously quantifying the expression levels of identified peptides from mass spectrometry data. Q-FISH analysis on human HCC and liver tissue samples identified many protein biomarkers that are highly relevant to HCC. Q-FISH can be a useful tool both for peptide identification and quantification on mass spectrometry data analysis. It may also prove to be more effective in discovering novel protein biomarkers than SEQUEST and other standard methods.

MOTIVATION: Gene-gene interactions are important contributors to complex biological traits. Multifactor dimensionality reduction (MDR) is a method to analyze gene-gene interactions and has been applied to many genetics studies of complex diseases. In order to identify the best interaction model associated with disease susceptibility, MDR classifiers corresponding to interaction models has been constructed and evaluated as a predictor of disease status via a certain measure such as balanced accuracy (BA). It has been shown that the performance of MDR tends to depend on the choice of the evaluation measures. RESULTS: In this article, we introduce two types of new evaluation measures. First, we develop weighted BA (wBA) that utilizes the quantitative information on the effect size of each multi-locus genotype on a trait. Second, we employ ordinal association measures to assess the performance of MDR classifiers. Simulation studies were conducted to compare the proposed measures with BA, a current measure. Our results showed that the wBA and tau(b) improved the power of MDR in detecting gene-gene interactions. Noticeably, the power increment was higher when data contains the greater number of genetic markers. Finally, we applied the proposed evaluation measures to real data.

UNLABELLED: Pharmaceutical industry has been striving to reduce the costs of drug development and increase productivity. Among the many different attempts, drug repositioning (retargeting existing drugs) comes into the spotlight because of its financial efficiency. We introduce IDMap which predicts novel relationships between targets and chemicals and thus is capable of repositioning the marketed drugs by using text mining and chemical structure information. Also capable of mapping commercial chemicals to possible drug targets and vice versa, IDMap creates convenient environments for identifying the potential lead and its targets, especially in the field of drug repositioning. AVAILABILITY: IDMap executable and its user manual including color images are freely available to non-commercial users at http://www.equispharm.com/idmap