The dChip survival analysis module for microarray data.  
Jump to Full Text  
MedLine Citation:

PMID: 21388547 Owner: NLM Status: MEDLINE 
Abstract/OtherAbstract:

BACKGROUND: Genomewide expression signatures are emerging as potential marker for overall survival and disease recurrence risk as evidenced by recent commercialization of gene expression based biomarkers in breast cancer. Similar predictions have recently been carried out using genomewide copy number alterations and microRNAs. Existing software packages for microarray data analysis provide functions to define expressionbased survival gene signatures. However, there is no software that can perform survival analysis using SNP array data or draw survival curves interactively for expressionbased sample clusters. RESULTS: We have developed the survival analysis module in the dChip software that performs survival analysis across the genome for gene expression and copy number microarray data. Built on the current dChip software's microarray analysis functions such as chromosome display and clustering, the new survival functions include interactive exploring of KaplanMeier (KM) plots using expression or copy number data, computing survival pvalues from the logrank test and Cox models, and using permutation to identify significant chromosome regions associated with survival. CONCLUSIONS: The dChip survival module provides userfriendly way to perform survival analysis and visualize the results in the context of genes and cytobands. It requires no coding expertise and only minimal learning curve for thousands of existing dChip users. The implementation in Visual C++ also enables fast computation. The software and demonstration data are freely available at http://dchipsurv.chenglilab.org. 
Authors:

Samir B Amin; Parantu K Shah; Aimin Yan; Sophia Adamia; Stéphane Minvielle; Hervé AvetLoiseau; Nikhil C Munshi; Cheng Li 
Related Documents
:

11004177  The toryz (yeck bisz) operon encodes a third respiratory trimethylamine noxide reducta... 23006437  Genomewide expressions in autologous eutopic and ectopic endometrium of fertile women ... 12101127  Spot 42 rna mediates discoordinate expression of the e. coli galactose operon. 6434517  Anaerobic and leucinedependent expression of a peptide transport gene in salmonella ty... 23719067  Current progress on gene therapy for primary immunodeficiencies. 12969137  Early aldosterone upregulated genes: new pathways for renal disease? 
Publication Detail:

Type: Journal Article; Research Support, N.I.H., Extramural; Research Support, NonU.S. Gov't Date: 20110309 
Journal Detail:

Title: BMC bioinformatics Volume: 12 ISSN: 14712105 ISO Abbreviation: BMC Bioinformatics Publication Date: 2011 
Date Detail:

Created Date: 20110401 Completed Date: 20110525 Revised Date: 20130630 
Medline Journal Info:

Nlm Unique ID: 100965194 Medline TA: BMC Bioinformatics Country: England 
Other Details:

Languages: eng Pagination: 72 Citation Subset: IM 
Affiliation:

Department of Biostatistics and Computational Biology, DanaFarber Cancer Institute and Harvard School of Public Health, 450 Brookline Ave, Boston, MA 02215, USA. 
Export Citation:

APA/MLA Format Download EndNote Download BibTex 
MeSH Terms  
Descriptor/Qualifier:

Breast Neoplasms
/
genetics Cluster Analysis Gene Expression Profiling / methods* Humans KaplanMeier Estimate Oligonucleotide Array Sequence Analysis / methods* Polymorphism, Single Nucleotide Proportional Hazards Models Software* 
Grant Support  
ID/Acronym/Agency:

1R01GM077122/GM/NIGMS NIH HHS 
Comments/Corrections 
Full Text  
Journal Information Journal ID (nlmta): BMC Bioinformatics ISSN: 14712105 Publisher: BioMed Central 
Article Information Download PDF Copyright ©2011 Amin et al; licensee BioMed Central Ltd. openaccess: Received Day: 27 Month: 9 Year: 2010 Accepted Day: 9 Month: 3 Year: 2011 collection publication date: Year: 2011 Electronic publication date: Day: 9 Month: 3 Year: 2011 Volume: 12First Page: 72 Last Page: 72 ID: 3068974 Publisher Id: 147121051272 PubMed Id: 21388547 DOI: 10.1186/147121051272 
The dChip survival analysis module for microarray data  
Samir B Amin12  Email: samir_amin@dfci.harvard.edu 
Parantu K Shah1  Email: parantu.shah@gmail.com 
Aimin Yan1  Email: jlfmssm@gmail.com 
Sophia Adamia2  Email: Sophia_Adamia@dfci.harvard.edu 
Stéphane Minvielle34  Email: stephane.minvielle@chunantes.fr 
Hervé AvetLoiseau34  Email: herve.avetloiseau@chunantes.fr 
Nikhil C Munshi25  Email: Nikhil_Munshi@dfci.harvard.edu 
Cheng Li1  Email: lch3000@gmail.com 
1Department of Biostatistics and Computational Biology, DanaFarber Cancer Institute and Harvard School of Public Health, 450 Brookline Ave, Boston, MA, 02215, USA 

2Department of Medical Oncology, DanaFarber Cancer Institute, 450 Brookline Ave, Boston, MA, 02215, USA 

3Hematology Department, University Hospital, Nantes, France 

4Centre de Recherche en Cancérologie, INSERM U892, Nantes, France 

5Veterans Administration Boston Healthcare System and Harvard Medical School. 1400 VFW Pkwy, West Roxbury, MA, 02132, USA 
In cancer clinical practice, predicting patient survival based on traditional tumor staging systems using clinical, histopathological and molecular markers remains an integral component in the treatment decision for patients. For example, patients with advanced disease and poor survival prognosis are subjected to more aggressive treatments. However, this conventional approach is nonspecific and has limited success in the cancer treatment. Many patients have recurrence despite having aggressive therapy based on survival risk score [^{1},^{2}].
With highthroughput cancer genomics data, we and others have reported using genomewide expression signatures to predict survival risk, and these signatures are now increasingly being used in treatment decision for several cancer types [^{3}^{}^{6}]. Survival predictions have also been carried out using genomewide copy number alterations [^{7},^{8}] and microRNAs [^{9},^{10}]. Encouraged by these results, researchers routinely analyze large sets of microarray data in relation to survival information. Common analysis tasks and endpoints include gene signatures that predict survival risk, survival difference between sample groups defined by unsupervised clustering, and survival analysis using the copynumber data of local genomic regions. Such survival analysis on a highdimensional data requires statistical programming and commandline skills, or the use of the existing software packages such as BRBArrayTools, Survival Online tool and Prediction Analysis for Microarrays (PAM) [^{11}^{}^{13}].
However, there is no specific utility that can perform survival analysis using SNP array data or draw survival curves interactively for expressionbased sample clusters. We have developed the widelyused dChip software that can efficiently process and derive gene expression and copy number data from microarray datasets (http://www.dchip.org) [^{14},^{15}], and have pioneered using SNP arrays to find chromosomal alterations such as amplification, deletion, and loss of heterozygosity (LOH) [^{16}]. Thus, the addition of survival functions will be helpful for researchers to query and correlate chromosomal regions of interest with associated survival data.
Here, we describe the survival analysis module in the dChip software that performs survival analysis across the genome for gene expression and copy number microarray data. The new survival functions include interactive exploring of KaplanMeier (KM) plots using both expression and copy number data, computing survival pvalues from the logrank test and Cox models, and using permutation to assess the survival significance of copy numbers genomewide. Researchers can also compare survival curves between sample clustering groups derived from expression data. The dChip survival module enables userfriendly, interactive survival analysis and visualization of microarray data in the context of genes and cytobands. It requires no need for coding and minimal learning curve for existing dChip users. The implementation in Visual C++ also enables fast computation for processing large data sets from studies such as the Cancer Genome Atlas (TCGA).
The survival analysis functions are implemented in dChip using Visual C++ and optimized for fast computation. The computed logrank test and Cox model statistics and pvalues are confirmed using R code. Figure 1 summarizes the preliminary raw data analysis and new workflow functions in two categories. a) those for SNP copy number data, and b) those for expressionbased sample clustering groups.
Here we will use two example data sets to demonstrate the functions: 1) performing survival analysis using SNP data, and 2) drawing KM plots using expressionbased sample clustering groups. For the first dataset [^{7}], we will discuss the following dChip analysis steps: SNP data input and normalization, plotting copy number data in the chromosome view, carrying out survival analysis using the logrank and Cox model, and the permutation function to adjust for multiple testing and assess the genomewide significance of the survival scores. For the second dataset, we will use a gene expression dataset consisting of 170 uniformly treated patients with multiple myeloma with clinical followup of more than five years (Munshi et al., manuscript in preparation). We will first perform unsupervised hierarchical clustering and define gene signatures that classify the samples into subgroups, and then compare KM curves by the logrank test among these subgroups.
We will use a 192sample microarray dataset processed on Affymetrix SNP 500 K microarray platform [^{7}] to illustrate the usage of the new analysis and visualization functions. This data set represents 192 uniformly treated patients with multiple myeloma (MM). MM is a common type of hematological cancer and characterized by malignant clonal transformation of plasma cells in bone marrow with excess production of a monoclonal immunoglobulin. Chromosomal aberrations are a hallmark of MM with specific changes (del13 and t(4;14)) giving poor prognosis and other changes (hyperdiploidy and t(11;14)) conferring better survival [^{17}].
We use the dChip software to normalize all the arrays of 192 myeloma samples and additional 10 normal blood samples to compute modelbased signal values. Normalization and modelbased signals are calculated for each of the two subarrays and combined. The median genotype call rates are 96.77% and 97.35% respectively for the 250 K Nsp and Sty subarrays. A tabdelimited sample information file is prepared with columns specifying survival outcome including survival time and event indicator (0 = alive, 1 = death). Numerical columns are marked in the column header such as "Survival(numeric)", and will be standardized and displayed above the samples in the clustering or chromosome data views. After loading the normalized data into dChip using the menu function "Analysis > Open group", we use the "Tools > Array List File" menu function to create an array list to order samples by the values of a particular sample variable. Array list files are also useful when doing survival analysis on a subset of samples, leaving other samples out of the downstream analysis.
The next step displays the data along chromosomes using the menu function "Analysis > Chromosome" to specify analysis parameters for copy number and LOH analysis (Figure 2). The copy number analysis functions are explained in detail in the dChip manual [^{14}]. When the ploidy of samples are unknown (i.e. tumor samples), we can check the option "Scale copy number mode to 2 copy" to adjust for ploidy effect samplewise, so that in hyperdiploid samples the chromosomes with normal copy numbers are estimated to have two copies rather than deletion events.
Once the data is displayed by chromosomes, we can toggle between showing individual and all chromosomes using the "Chromosome > Show All" menu. Figure 3 shows the gain and loss of copy numbers in the chromosome view, with SNPs on the rows and samples on the columns. It also uses an array list file to sort all the samples by survival time irrespective of the event indicator. We can browse the genome and quickly observe whether copy gain and loss events are associated with survival. Figure 3 shows that hyperdiploid samples with three copies in chromosome 5 tend to locate on the right side (thus longer survival time) rather than on the left side, suggesting hyperdiploidy is associated with better survival outcome of myeloma patients.
In the first survival analysis function, we use the menu function "Chromosome > Compute Score" to perform genomewide, singleSNP survival analysis by selecting the "Survival LogRank pvalue" option under "Scoring method" (Figure 4). It will perform SNPwise logrank test for association between copy numbers and survival outcome. The data variables of either event free survival or overall survival can be selected by clicking "Response variable" to open the "Select factors" dialog (Figure 4). For a SNP, all the samples are divided into three groups based on this SNP's inferred copy number: Deletion (≤ 1.5 copy), Gain (≥ 2.5 copy), and "No change" (between 1.5 and 2.5 copy). This copy thresholds are adjustable at the "Options > Score" dialog. The logrank test is then applied to the three sample groups, testing the null hypothesis that there is no survival difference between the groups.
Across the genome, we compute the logrank test score for one of every 10 consecutive SNPs to speed up the computation. This is because the inferred copy numbers are correlated for nearby SNPs when computed from 10SNP local median smoothing, leading to correlated logrank test scores for nearby SNPs. The window size of 10 is useradjustable at the dialog option "Tools > Options > Score > Use one of every 10 markers". Once the logrank test is applied to all the SNPs, dChip will display the survival scores genomewide (Figure 3). The blue curve on the right displays the negative log10 transformed pvalues from the logrank test, with larger values indicating stronger association between copy numbers and survival outcome. The vertical red line in the gray box indicates the score threshold for the genomewide significance at the 0.05 level by the MaxT permutation method (see the section below on permutation). One or more regions across the genome could have scores exceeding the threshold. Inspecting the genes within these regions could identify candidate survivalassociated genes (Figure 5).
To draw a KaplanMeier plot (KM plot) for an individual SNP, we first click on either a SNP nearby a gene of interest (the menu function "Chromosome > Find Gene" can search for genes) or a SNP from the regions with significant survival scores. Then, we use the "View > KM plot" function to draw a KM plot using the copy number data of the SNP to form Deletion, Gain, and "No change" groups as above (Figure 6). The pvalue to test the different survival rates between the groups is computed using the logrank test when there are two or more groups each with 5 or more expected number of events. The SNP (rs16915104) in this KM plot is within the transcribed region of the PLEKHA5 (PEPP2), an Xlinked human homeobox gene at chromosome 12p12.3, which encodes transcription factor with known oncogenic role in cancer and drug resistance via phosphoinositidemediated signal pathways [^{18},^{19}]. Using this approach, we can first perform genomewide screen of SNPs to locate chromosome regions of significant logrank scores, and then zoom into the peak regions and use KM plots to check the survival correlation of the SNPs and genes in the region.
The second analysis function performs univariate Cox regression for survival outcome using the inferred copy numbers of SNPs. After calculating and displaying the inferred copy numbers in the chromosome view, we use the menu function "Chromosome > Compute Score" and select the "Survival Cox Regression" option under "Scoring method" to perform univariate Cox regression for individual SNPs across the genome (Figure 4). The computing time will vary depending on the number of samples and SNPs; it takes less than 10 minutes to process the 500 K SNP, 192 sample data on a 2.4 GHz CPU with 3 GB RAM.
This function displays the absolute zscores computed from the Cox model in the gray box on the right of the copy number view (Figure 3 and 5). Negative zscores are displayed in blue, indicating less hazard or longer survival as copy number increases: copy number gains are associated with longer survival, or copy number deletions are associated with shorter survival. Positive zscores are displayed in red, indicating higher hazard or shorter survivals as copy number increases. As before, ordering samples based on survival time visualizes and confirms the association between copy number alteration events and survival variables. Other continuous or binary sample variables can also be selected as "Covariates" to perform multivariate Cox regression to assess the explanatory power of a SNP's copy number to survival in the context of these variables (Figure 4). Permutation can be applied to determine the genomewide threshold to call chromosome regions significantly associated with survival outcome.
We also extend the existing permutation function in dChip to assess the genomewide significance of the survival scores derived from the logrank test or Cox regression. The number of permutation runs is specified at the "Chromosome > Compute Score" dialog (Figure 4). The permutation is computationally intensive but also benefits from computing for only a subset of SNPs across the genome.
The permutation tests the null hypothesis that there is no chromosome region in the cancer genome whose copy number is associated with survival, and therefore any observed association is due to random chance. To simulate data sets under the null hypothesis, we can either permute survival times with censoring indicators across samples or permute chromosome region blocks within every sample. The survival scores from the simulated data sets are then compared to those from the original data set. Specifically for permuting chromosome region blocks, for each sample, whose SNPs are ordered first by chromosomes and then by positions within chromosome, we randomly partition the whole genome into K (≥ 2) blocks, and randomly switch the order of these blocks while preserving the order of SNPs within each block. In this way, the SNPs and their copy numbers in a sample are randomly relocated in blocks to new positions in the genome, while only minimally perturbing the dependence of the copy number data of neighboring SNPs. The same permutation applies to all samples using a different random partition for each sample. The survival score at each SNP locus can then be computed for the permuted data set, and the MaxT method can be applied to assess the significance of the original scores [^{20}]. The maximal survival scores from every permuted data set form the score distribution, whose 95^{th }largest value is the genomewide threshold at the 0.05 significance level to determine the chromosome regions significantly associated with survival in the original dataset. Similarly, clinical variables such as survival times and censoring indicators can be permuted together among all the samples and this achieves similar genomewide significance threshold (Additional file 1: Figure S1 and Additional file 2: Figure S2).
Unsupervised hierarchical clustering is frequently used to discover novel sample subgroups from microarray data or inspect the expressionbased samples clusters in relation to clinical variables. In particular, it is interesting to know whether sample clusters correspond to differential survival outcomes.
We develop a third analysis function in dChip to facilitate drawing KM plots based on sample clusters. We first use the "Analysis > ANOVA & Correlation" menu to specify the survival variable and event indicator, similar to Figure 4. Next, we use "Analysis > Clustering & Enrichment" menu to cluster samples using a variationfiltered gene list. Interested in whether sample clusters correspond to different survival outcomes, we can click to select a main sample cluster branch (in blue color) and use Controlclick to select and color additional sample clusters (Figure 7A). We then use the menu "View/KM Plot" to display the KaplanMeier plot and logrank pvalue based on the specified sample clusters (Figure 7B). The dChip Analysis View will also show the details of logrank test and the number of samples omitted due to missing survival data.
Gene expression changes and copynumber alterations are known to influence cancer progression and patient survival. Currently, the survival analysis of genomewide copy number data using SNP arrays and the assessment of significance requires specialized statistical software and skills. We have developed the survival analysis module within the dChip software to streamline the survival analysis and interactive visualization of SNP copy number data and expressionbased sample clusters. We use data analysis examples to show how dChip can interactively explore KM plots and find survival associated genomic regions of interest. The easy userinterface and fast computing make these dChip functions accessible to biomedical researchers.
Survival analysis using copy number data provides options to use either the logrank test or Cox regression model to compute survival association. The genomewide view after both analysis methods displays a sidecurve showing survival scores for consecutive SNPs. The results from the dChip survival functions agree with previously published results on the same dataset [^{7}]. The amplifications in chromosome 1q and deletions in 1p and 16q are associated with poor survival, whereas the copy number gain of chromosomes 5 and 11 predicts a good prognosis.
We can change the size of SNP window for copynumber smoothing from the default 10 SNPs to other numbers. Increasing the window size will speed up the computation when sample size is large. But we need to be cautious that large window sizes may miss potential survival association of SNPs, especially those SNPs near or within the gene coding regions.
Unsupervised hierarchical clustering is a standard approach to analyzing expression profiling. Researchers are often interested in whether samples clusters correspond to different eventfree survival, overall survival, or treatment response outcomes. We have also developed functions in dChip to draw sample clusterbased KM plots. A user can select two or multiple cluster nodes to compute, display, and export highquality survival figures. To our knowledge, few software packages provide such graphical interface to ease the survival analysis without coding. We find Survival Online (SO) tool by Corradi et al. [^{12}] a useful online portal for Cox regression and survival analysis using gene expression data. At present, dChipSurv provides similar analysis for SNP array data. Together, these two applications will provide complementary set of features to users in need of survival analysis using expression and copynumber microarray data.
In summary, the dChip survival module addresses the frequent need of many researchers to integrate survival data analysis under a single microarray analysis package with minimal learning curve, fast computing, and no requirement of programming skills. We will add more survival functions as well as extend them to microRNA and RNAseq data in future dChip versions.
Project name: dChip survival analysis module
Project home page: http://dchipsurv.chenglilab.org/
Operating system(s): Windows 2000 or after
Programming language: Visual C++ 2005
Source code: Available on request.
License: The software is freely available.
CL designed and implemented the software module, S Amin and CL wrote the manuscript, S Amin, PKS and AY contributed to the data analysis and manuscript preparation, S Adamia, SM, HA and NCM designed and performed the microarray experiments.
All authors read and approved the final draft.
Figure S1
Click here for additional data file (147121051272S1.DOCX)
Additional file 2
Figure S2
Click here for additional data file (147121051272S2.DOC)
We thank Cheng Li Lab members for constructive discussion. This work has been supported by NIH grant 1R01GM077122. P.K.S. has been supported by the Claudia Adams Barr Program in Innovative Basic Cancer Research.
References
Vigano A,Dorgan M,Buckingham J,Bruera E,SuarezAlmazor ME,Survival prediction in terminal cancer patients: a systematic review of the medical literaturePalliat MedYear: 200014536337410.1191/02692160070153619211064783  
Williams C,Brunskill S,Altman D,Briggs A,Campbell H,Clarke M,Glanville J,Gray A,Harris A,Johnston K,et al. Costeffectiveness of using prognostic information to select women with breast cancer for adjuvant systemic therapyHealth Technol AssessYear: 20061034iiiiv ixxi, 1204. 16959170  
Kim C,Paik S,Geneexpressionbased prognostic assays for breast cancerNat Rev Clin OncolYear: 20107634034710.1038/nrclinonc.2010.6120440284  
van de Vijver M,He Y,van't Veer L,Dai H,Hart A,Voskuil D,Schreiber G,Peterse J,Roberts C,Marton M,et al. A geneexpression signature as a predictor of survival in breast cancerN Engl J MedYear: 2002347251999200910.1056/NEJMoa02196712490681  
Decaux O,Lodé L,Magrangeas F,Charbonnel C,Gouraud W,Jézéquel P,Attal M,Harousseau J,Moreau P,Bataille R,et al. Prediction of survival in multiple myeloma based on gene expression profiles reveals cell cycle and chromosomal instability signatures in highrisk patients and hyperdiploid signatures in lowrisk patients: a study of the Intergroupe Francophone du MyélomeJ Clin OncolYear: 200826294798480510.1200/JCO.2007.13.854518591550  
Carter S,Eklund A,Kohane I,Harris L,Szallasi Z,A signature of chromosomal instability inferred from gene expression profiles predicts clinical outcome in multiple human cancersNat GenetYear: 20063891043104810.1038/ng186116921376  
AvetLoiseau H,Li C,Magrangeas F,Gouraud W,Charbonnel C,Harousseau J,Attal M,Marit G,Mathiot C,Facon T,et al. Prognostic significance of copynumber alterations in multiple myelomaJ Clin OncolYear: 200927274585459010.1200/JCO.2008.20.613619687334  
Zhang Y,Martens J,Yu J,Jiang J,Sieuwerts A,Smid M,Klijn J,Wang Y,Foekens J,Copy number alterations that predict metastatic capability of human breast cancerCancer ResYear: 20096993795380110.1158/00085472.CAN08459619336569  
Schetter A,Leung S,Sohn J,Zanetti K,Bowman E,Yanaihara N,Yuen S,Chan T,Kwong D,Au G,et al. MicroRNA expression profiles associated with prognosis and therapeutic outcome in colon adenocarcinomaJAMAYear: 2008299442543610.1001/jama.299.4.42518230780  
Ji J,Shi J,Budhu A,Yu Z,Forgues M,Roessler S,Ambs S,Chen Y,Meltzer P,Croce C,et al. MicroRNA expression, survival, and response to interferon in liver cancerN Engl J MedYear: 2009361151437144710.1056/NEJMoa090128219812400  
Simon R,Lam A,Li MC,Ngan M,Menenzes S,Zhao Y,Analysis of gene expression data using BRBArrayToolsCancer InformYear: 20073111719455231  
Corradi L,Mirisola V,Porro I,Torterolo L,Fato M,Romano P,Pfeffer U,Survival Online: a webbased service for the analysis of correlations between gene expression and clinical and followup dataBMC BioinformaticsYear: 200910Suppl 12S1010.1186/1471210510S12S1019828070  
Tibshirani R,Hastie T,Narasimhan B,Chu G,Diagnosis of multiple cancer types by shrunken centroids of gene expressionProc Natl Acad Sci USAYear: 200299106567657210.1073/pnas.08209929912011421  
Li C,Amin S,Fan J, Lin X, Liu JSAnalysis of cancer genome alterations using single nucleotide polymorphism (SNP) microarraysNew developments in biostatistics and bioinformaticsYear: 2009Singapore: Higher Education Press & World Scientific Publishing209238 Frontiers of Statistics, vol 1. full_text.  
Li C,Wong W,Modelbased analysis of oligonucleotide arrays: expression index computation and outlier detectionProc Natl Acad Sci USAYear: 2001981313610.1073/pnas.01140409811134512  
Zhao X,Li C,Paez JG,Chin K,Janne PA,Chen TH,Girard L,Minna J,Christiani D,Leo C,et al. An integrated view of copy number and allelic alterations in the cancer genome using single nucleotide polymorphism arraysCancer ResYear: 20046493060307110.1158/00085472.CAN03330815126342  
Fonseca R,Barlogie B,Bataille R,Bastard C,Bergsagel PL,Chesi M,Davies FE,Drach J,Greipp PR,Kirsch IR,et al. Genetics and cytogenetics of multiple myeloma: a workshop reportCancer ResYear: 20046441546155810.1158/00085472.CAN03287614989251  
Hofmann O,Caballero OL,Stevenson BJ,Chen YT,Cohen T,Chua R,Maher CA,Panji S,Schaefer U,Kruger A,et al. Genomewide analysis of cancer/testis gene expressionProc Natl Acad Sci USAYear: 200810551204222042710.1073/pnas.081077710519088187  
Wittig R,Nessling M,Will RD,Mollenhauer J,Salowsky R,Munstermann E,Schick M,Helmbach H,Gschwendt B,Korn B,et al. Candidate genes for crossresistance against DNAdamaging drugsCancer ResYear: 200262226698670512438269  
Westfall PH,Young SS,Resamplingbased multiple testing: examples and methods for Pvalue adjustmentYear: 1993New York: Wiley 
Figures
Article Categories:

Previous Document: Psychological approach to successful ageing predicts future quality of life in older adults.
Next Document: Fish oil rich diet in comparison to saturated fat rich diet offered protection against lipopolysacch...