A novel integrative pipeline is presented for discovery of potential cancer-susceptibility regions (PCSRs) by calculating the number of altered genes at each chromosomal region, using expression microarray datasets of different human cancers (HCs). Our novel approach comprises primarily predicting PCSRs followed by identification of key genes in these regions to obtain potential regions harboring new cancer-associated variants. In addition to finding new cancer causal variants, another advantage in prediction of such risk regions is simultaneous study of different types of genomic variants in line with focusing on specific chromosomal regions. Using this pipeline we extracted numbers of regions with highly altered expression levels in cancer condition. Regulatory networks were also constructed for different types of cancers following the identification of altered mRNA and microRNAs. Interestingly, results showed that GAPDH, LIFR, ZEB2, mir-21, mir-30a, mir-141 and mir-200c, all located at PCSRs, are common altered factors in constructed networks. We found a number of clusters of altered mRNAs and miRNAs on predicted PCSRs (e.g.12p13.31) and their common regulators including KLF4 and SOX10. Large scale prediction of risk regions based on transcriptome data can open a window in comprehensive study of cancer risk factors and the other human diseases.
References
[1]
Nam EJ, Yoon H, Kim SW, Kim H, Kim YT, et al. (2008) MicroRNA expression profiles in serous ovarian carcinoma. Clinical Cancer Research 14: 2690–2695. doi: 10.1158/1078-0432.ccr-07-1731
[2]
Catto JW, Alcaraz A, Bjartell AS, De Vere White R, Evans CP, et al. (2011) MicroRNA in prostate, bladder, and kidney cancer: a systematic review. European urology 59: 671–681. doi: 10.1016/j.eururo.2011.01.044
[3]
Alanazi I, Ebrahimie E, Hoffmann P, Adelson DL (2013) Combined gene expression and proteomic analysis of EGF induced apoptosis in A431 cells suggests multiple pathways trigger apoptosis. Apoptosis: 1–15.
[4]
Huang L, Yu D, Wu C, Zhai K, Jiang G, et al. (2012) Copy number variation at 6q13 functions as a long-range regulator and is associated with pancreatic cancer risk. Carcinogenesis 33: 94–100. doi: 10.1093/carcin/bgr228
[5]
Hoggart CJ, Whittaker JC, De Iorio M, Balding DJ (2008) Simultaneous analysis of all SNPs in genome-wide and re-sequencing association studies. PLoS genetics 4: e1000130. doi: 10.1371/journal.pgen.1000130
[6]
Bakhtiarizadeh MR, Ebrahimi M, Ebrahimie E (2011) Discovery of EST-SSRs in Lung Cancer: Tagged ESTs with SSRs Lead to Differential Amino Acid and Protein Expression Patterns in Cancerous Tissues. PLoS One 6: e27118. doi: 10.1371/journal.pone.0027118
[7]
Merup M, Juliusson G, Wu X, Jansson M, Stellan B, et al. (2009) Amplification of multiple regions of chromosome 12, including 12q13–15, in chronic lymphocytic leukaemia. European journal of haematology 58: 174–180. doi: 10.1111/j.1600-0609.1997.tb00944.x
[8]
Schafmayer C, Buch S, V？lzke H, von Sch？nfels W, Egberts JH, et al. (2009) Investigation of the colorectal cancer susceptibility region on chromosome 8q24. 21 in a large German case-control sample. International Journal of Cancer 124: 75–80. doi: 10.1002/ijc.23872
[9]
Fehringer G, Liu G, Pintilie M, Sykes J, Cheng D, et al. (2012) Association of the 15q25 and 5p15 Lung Cancer Susceptibility Regions with Gene Expression in Lung Tumor Tissue. Cancer Epidemiology Biomarkers & Prevention 21: 1097–1104. doi: 10.1158/1055-9965.epi-11-1123-t
[10]
Bolton KL, Tyrer J, Song H, Ramus SJ, Notaridou M, et al. (2010) Common variants at 19p13 are associated with susceptibility to ovarian cancer. Nature genetics 42: 880–884.
[11]
Wu C, Miao X, Huang L, Che X, Jiang G, et al. (2012) Genome-wide association study identifies five loci associated with susceptibility to pancreatic cancer in Chinese populations. Nat Genet 44: 62–66. doi: 10.1038/ng.1020
[12]
Lan Q, Hsiung CA, Matsuo K, Hong Y-C, Seow A, et al. (2012) Genome-wide association analysis identifies new lung cancer susceptibility loci in never-smoking women in Asia. Nature genetics 44: 1330–1335.
[13]
Cropp CD, Simpson CL, Wahlfors T, Ha N, George A, et al. (2011) Genome-wide linkage scan for prostate cancer susceptibility in Finland: Evidence for a novel locus on 2q37. 3 and confirmation of signal on 17q21-q22. International Journal of Cancer 129: 2400–2407. doi: 10.1002/ijc.25906
[14]
Gudmundsson J, Sulem P, Gudbjartsson DF, Masson G, Agnarsson BA, et al.. (2012) A study based on whole-genome sequencing yields a rare variant at 8q24 associated with prostate cancer. Nature genetics.
[15]
Xu J, Mo Z, Ye D, Wang M, Liu F, et al.. (2012) Genome-wide association study in Chinese men identifies two new prostate cancer risk loci at 9q31. 2 and 19q13. 4. Nature genetics.
[16]
Monzon FA, Alvarez K, Peterson L, Truong L, Amato RJ, et al. (2011) Chromosome 14q loss defines a molecular subtype of clear-cell renal cell carcinoma associated with poor prognosis. Mod Pathol 24: 1470–1479. doi: 10.1038/modpathol.2011.107
[17]
Freedman ML, Monteiro ANA, Gayther SA, Coetzee GA, Risch A, et al. (2011) Principles for the post-GWAS functional characterization of cancer risk loci. Nat Genet 43: 513–518. doi: 10.1038/ng.840
[18]
Letessier A, Sircoulomb F, Ginestier C, Cervera N, Monville F, et al. (2006) Frequency, prognostic impact, and subtype association of 8p12, 8q24, 11q13, 12p13, 17q12, and 20q13 amplifications in breast cancers. BMC cancer 6: 245.
[19]
Tenesa A, Farrington SM, Prendergast JG, Porteous ME, Walker M, et al. (2008) Genome-wide association scan identifies a colorectal cancer susceptibility locus on 11q23 and replicates risk loci at 8q24 and 18q21. Nature genetics 40: 631–637. doi: 10.1038/ng.133
[20]
Goode EL, Chenevix-Trench G, Song H, Ramus SJ, Notaridou M, et al. (2010) A genome-wide association study identifies susceptibility loci for ovarian cancer at 2q31 and 8q24. Nature genetics 42: 874–879. doi: 10.1038/ng.668
[21]
Fan B, Dachrut S, Coral H, Yuen ST, Chu KM, et al. (2012) Integration of DNA copy number alterations and transcriptional expression analysis in human gastric cancer. PLoS One 7: e29824. doi: 10.1371/journal.pone.0029824
[22]
Hobbs MR, Jones BB, Otterud BE, Leppert M, Kriesel JD (2008) Identification of a herpes simplex labialis susceptibility region on human chromosome 21. Journal of Infectious Diseases 197: 340–346. doi: 10.1086/525540
[23]
Chen Z-J, Zhao H, He L, Shi Y, Qin Y, et al. (2010) Genome-wide association study identifies susceptibility loci for polycystic ovary syndrome on chromosome 2p16. 3, 2p21 and 9q33. 3. Nature genetics 43: 55–59. doi: 10.1038/ng.732
[24]
Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature 447: 661–678.
[25]
Couturier J, Vielh P, Salmon R, Dutrillaux B (2006) Trisomy and tetrasomy for long arm of chromosome 1 in near-diploid human endometrial adenocarcinomas. International Journal of Cancer 38: 17–19. doi: 10.1002/ijc.2910380104
[26]
Polascik TJ, Cairns P, Chang WY, Schoenberg MP, Sidransky D (1995) Distinct regions of allelic loss on chromosome 4 in human primary bladder carcinoma. Cancer Research 55: 5396–5399.
[27]
Liu TX, Becker MW, Jelinek J, Wu W-S, Deng M, et al. (2006) Chromosome 5q deletion and epigenetic suppression of the gene encoding α-catenin (CTNNA1) in myeloid cell transformation. Nature medicine 13: 78–83. doi: 10.1038/nm1512
[28]
Nana-Sinkam S, Croce C (2012) Clinical applications for microRNAs in cancer. Clinical Pharmacology & Therapeutics.
[29]
Wang D, Huang J, Hu Z (2012) RNA helicase DDX5 regulates MicroRNA expression and contributes to cytoskeletal reorganization in basal breast cancer cells. Molecular & Cellular Proteomics 11..
[30]
Wee EJH, Peters K, Nair SS, Hulf T, Stein S, et al.. (2012) Mapping the regulatory sequences controlling 93 breast cancer-associated miRNA genes leads to the identification of two functional promoters of the Hsa-mir-200b cluster, methylation of which is associated with metastasis or hormone receptor status in advanced breast cancer. Oncogene.
[31]
Mazurek A, Luo W, Krasnitz A, Hicks J, Powers RS, et al. (2012) DDX5 Regulates DNA Replication and Is Required for Cell Proliferation in a Subset of Breast Cancer Cells. Cancer Discovery 2: 812–825. doi: 10.1158/2159-8290.cd-12-0116
[32]
Hosseinpour B, Bakhtiarizadeh MR, Khosravi P, Ebrahimie E (2013) Predicting distinct organization of transcription factor binding sites on the promoter regions; a new genome-based approach to expand human embryonic stem cell regulatory network. Gene.
[33]
Zhu Y, McAvoy S, Kuhn R, Smith D (2006) RORA, a large common fragile site gene, is involved in cellular stress response. Oncogene 25: 2901–2908. doi: 10.1038/sj.onc.1209314
[34]
Shum C, Lau S, Tsoi L, Chan L, Yam J, et al.. (2012) Krüppel-like factor 4 (KLF4) suppresses neuroblastoma cell growth and determines non-tumorigenic lineage differentiation. Oncogene.
[35]
Bockhorn J, Dalton R, Nwachukwu C, Huang S, Prat A, et al. (2013) MicroRNA-30c inhibits human breast tumour chemotherapy resistance by regulating TWF1 and IL-11. Nat Commun 4: 1393. doi: 10.1038/ncomms2393
[36]
Weng Y, Chen Y, Chen J, Liu Y, Bao T (2013) Common genetic variants in the microRNA biogenesis pathway are associated with malignant peripheral nerve sheath tumor risk in a Chinese population. Cancer Epidemiology.
[37]
Guenthoer J, Diede SJ, Tanaka H, Chai X, Hsu L, et al. (2012) Assessment of palindromes as platforms for DNA amplification in breast cancer. Genome research 22: 232–245. doi: 10.1101/gr.117226.110
[38]
Chan W-Y, Wong N, Chan A, Chow J, Lee J (1998) Consistent copy number gain in chromosome 12 in primary diffuse large cell lymphomas of the stomach. The American journal of pathology 152: 11.
[39]
Balatti V, Bottoni A, Palamarchuk A, Alder H, Rassenti LZ, et al. (2012) NOTCH1 mutations in CLL associated with trisomy 12. Blood 119: 329–331. doi: 10.1182/blood-2011-10-386144
[40]
Chen R-W, Saunders PA, Wei H, Li Z, Seth P, et al. (1999) Involvement of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and p53 in neuronal apoptosis: evidence that GAPDH is upregulated by p53. The Journal of neuroscience 19: 9654–9662.
[41]
Yoshihara K, Tajima A, Adachi S, Quan J, Sekine M, et al. (2011) Germline copy number variations in BRCA1-associated ovarian cancer patients. Genes, Chromosomes and Cancer 50: 167–177. doi: 10.1002/gcc.20841
[42]
Izumi K, Conlin LK, Berrodin D, Fincher C, Wilkens A, et al. (2012) Duplication 12p and Pallister–Killian syndrome: A case report and review of the literature toward defining a Pallister–Killian syndrome minimal critical region. American Journal of Medical Genetics Part A 158: 3033–3045. doi: 10.1002/ajmg.a.35500
[43]
Karambataki M, Malousi A, Maglaveras N, Kouidou S (2010) Synonymous polymorphisms at splicing regulatory sites are associated with CpGs in neurodegenerative disease-related genes. Neuromolecular medicine 12: 260–269. doi: 10.1007/s12017-009-8111-0
[44]
Thayanithy V, Sarver AL, Kartha RV, Li L, Angstadt AY, et al. (2012) Perturbation of 14q32 miRNAs-cMYC gene network in osteosarcoma. Bone 50: 171–181. doi: 10.1016/j.bone.2011.10.012
[45]
Bakhtiarizadeh MR, Moradi-Shahrbabak M, Ebrahimie E (2013) Underlying functional genomics of fat deposition in adipose tissue. Gene.
[46]
Nikitin A, Egorov S, Daraselia N, Mazo I (2003) Pathway studio—the analysis and navigation of molecular networks. Bioinformatics 19: 2155–2157. doi: 10.1093/bioinformatics/btg290
[47]
Hosseinpour B, HajiHoseini V, Kashfi R, Ebrahimie E, Hemmatzadeh F (2012) Protein interaction network of Arabidopsis thaliana female gametophyte development identifies novel proteins and relations. PLoS One 7: e49931. doi: 10.1371/journal.pone.0049931
[48]
Huang DW, Sherman BT, Lempicki RA (2008) Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protocols 4: 44–57. doi: 10.1038/nprot.2008.211
[49]
Zambelli F, Pesole G, Pavesi G (2009) Pscan: finding over-represented transcription factor binding site motifs in sequences from co-regulated or co-expressed genes. Nucleic acids research 37: W247–W252. doi: 10.1093/nar/gkp464
