A supervised approach for identifying discriminating genotype patterns and its application to breast cancer data

MOTIVATION Large-scale association studies, investigating the genetic determinants of a phenotype of interest, are producing increasing amounts of genomic variation data on human cohorts. A fundamental challenge in these studies is the detection of genotypic patterns that discriminate individuals exhibiting the phenotype under study from individuals that do not possess it. The difficulty stems from the large number of single nucleotide polymorphism (SNP) combinations that have to be tested. The discrimination problem becomes even more involved when additional high-throughput data, such as gene expression data, are available for the same cohort. RESULTS We have developed a graph theoretic approach for identifying discriminating patterns (DPs) for a given phenotype in a genotyped population. The method is based on representing the SNP data as a bipartite graph of individuals and their SNP states, and identifying fully connected subgraphs of this graph that relate individuals enriched for a given phenotypic group. The method can handle additional data types such as expression profiles of the genotyped population. It is reminiscent of biclustering approaches with the crucial difference that its search process is guided by the phenotype under consideration in a supervised manner. We tested our approach in simulations and on real data. In simulations, our method was able to retrieve planted patterns with high success rate. We then applied our approach to a dataset of 72 breast cancer patients with available gene expression profiles, genotyped over 695 SNPs. We detected several DPs that were highly significant with respect to various clinical phenotypes, and investigated the groups of patients and the groups of genes they defined. We found the patient groups to be highly enriched for other phenotypes and to display expression coherency among their profiles. The gene groups displayed functional coherency and involved genes with known role in cancer, providing additional support to their involvement. AVAILABILITY The program is available upon request.

[1]  Dorit S. Hochbaum,et al.  Approximating Clique and Biclique Problems , 1998, J. Algorithms.

[2]  Y. Benjamini,et al.  Resampling-based false discovery rate controlling multiple test procedures for correlated test statistics , 1999 .

[3]  C. Sing,et al.  Complex adaptive systems and human health: the influence of common genotypes of the apolipoprotein E (ApoE) gene polymorphism and age on the relational order within a field of lipid metabolism traits , 2000, Human Genetics.

[4]  R. Tibshirani,et al.  Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[5]  J. Haines,et al.  Association of single-nucleotide polymorphisms of the tau gene with late-onset Parkinson disease. , 2001, JAMA.

[6]  R. Waterland,et al.  Mannose 6‐phosphate/insulin‐like growth factor 2 receptor (M6P/IGF2R) variants in American and Japanese populations , 2001, Human mutation.

[7]  R. Clough,et al.  M6P/IGF2R loss of heterozygosity in head and neck cancer associated with poor patient prognosis , 2003, BMC Cancer.

[8]  Dror Rawitz,et al.  Approximating Element-Weighted Vertex Deletion Problems for the Complete k-Partite Property , 2002, J. Algorithms.

[9]  Roded Sharan,et al.  Discovering statistically significant biclusters in gene expression data , 2002, ISMB.

[10]  R. Nicholson,et al.  Structure-function analysis of LIV-1, the breast cancer-associated protein that belongs to a new subfamily of zinc transporters. , 2003, The Biochemical journal.

[11]  Joseph T. Chang,et al.  Spectral biclustering of microarray data: coclustering genes and conditions. , 2003, Genome research.

[12]  R. Tibshirani,et al.  Repeated observation of breast tumor subtypes in independent gene expression data sets , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[13]  Y. Yatabe,et al.  Reduced Expression of the let-7 MicroRNAs in Human Lung Cancers in Association with Shortened Postoperative Survival , 2004, Cancer Research.

[14]  Jason H. Moore,et al.  STUDENTJAMA. The challenges of whole-genome approaches to common diseases. , 2004, JAMA.

[15]  D. Thomas,et al.  Statistical Methods in GeneticEpidemiology , 2014 .

[16]  Y. Ritov,et al.  Acetylcholinesterase/paraoxonase genotype and expression predict anxiety scores in Health, Risk Factors, Exercise Training, and Genetics study. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[17]  M. Ellis,et al.  The mannose 6-phosphate/insulin-like growth factor 2 receptor (M6P/IGF2R), a putative breast tumor suppressor gene , 1998, Breast Cancer Research and Treatment.

[18]  G. Church,et al.  Modular epistasis in yeast metabolism , 2005, Nature Genetics.

[19]  P. Donnelly,et al.  Genome-wide strategies for detecting multiple loci that influence complex diseases , 2005, Nature Genetics.

[20]  T. Deguchi,et al.  Association of Lys173Arg polymorphism with CYP11B2 expression in normal adrenal glands and aldosterone-producing adenomas. , 2005, The Journal of clinical endocrinology and metabolism.

[21]  L. Altucci,et al.  Tumor-selective action of HDAC inhibitors involves TRAIL induction in acute myeloid leukemia cells , 2005, Nature Medicine.

[22]  A. Tsalenko,et al.  Genetic variation in putative regulatory loci controlling gene expression in breast cancer , 2006, Proceedings of the National Academy of Sciences of the United States of America.