Multi-omic analysis suggests tumor suppressor genes evolved specific promoter features to optimize cancer resistance

Tumor suppressor genes (TSGs) exhibit distinct evolutionary features. We speculated that TSG promoters could have evolved specific features that facilitate their tumor-suppressing functions. We found that the promoter CpG dinucleotide frequencies of TSGs are significantly higher than that of non-cancer genes across vertebrate genomes, and positively correlated with gene expression across tissue types. The promoter CpG dinucleotide frequencies of all genes gradually increase with gene age, for which young TSGs have been subject to a stronger evolutionary pressure. Transcription-related features, namely chromatin accessibility, methylation and ZNF263-, SP1-, E2F4- and SP2-binding elements, are associated with gene expression. Moreover, higher promoter CpG dinucleotide frequencies and chromatin accessibility are positively associated with the ability of TSGs to resist downregulation during tumorigenesis. These results were successfully validated with independent datasets. In conclusion, TSGs evolved specific promoter features that optimized cancer resistance through achieving high expression in normal tissues and resistance to downregulation during tumorigenesis.

[1]  M. Esteller.  Epigenetics in cancer. , 2008, The New England journal of medicine.

[2]  Chao Cheng,et al.  Relating gene expression evolution with CpG content changes , 2014, BMC Genomics.

[3]  A. Jemal,et al.  Cancer statistics, 2016 , 2016, CA: a cancer journal for clinicians.

[4]  Charles Elkan,et al.  Fitting a Mixture Model By Expectation Maximization To Discover Motifs In Biopolymer , 1994, ISMB.

[5]  S. Yi,et al.  DNA methylation and evolution of duplicate genes , 2014, Proceedings of the National Academy of Sciences.

[6]  Anton Nekrutenko,et al.  Evolutionary dynamics of oncogenes and tumor suppressor genes: higher intensities of purifying selection than other genes. , 2003, Molecular biology and evolution.

[7]  Macarena Morillo-Huesca,et al.  Crosstalk between chromatin structure, cohesin activity and transcription , 2019, Epigenetics & Chromatin.

[8]  R Holliday,et al.  DNA methylation and mutation. , 1993, Mutation research.

[9]  D. Brutlag,et al.  A genome-wide analysis of CpG dinucleotides in the human genome distinguishes two distinct classes of promoters , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[10]  Peter A. Jones,et al.  The fundamental role of epigenetic events in cancer , 2002, Nature Reviews Genetics.

[11]  M. Buck,et al.  Chromatin accessibility: a window into the genome , 2014, Epigenetics & Chromatin.

[12]  Michael Gribskov,et al.  Combining evidence using p-values: application to sequence homology searches , 1998, Bioinform..

[13]  Lukas Burger,et al.  CG dinucleotides enhance promoter activity independent of DNA methylation , 2019, Genome research.

[14]  M. Esteller.  CpG island hypermethylation and tumor suppressor genes: a booming present, a brighter future , 2002, Oncogene.

[15]  William Stafford Noble,et al.  Quantifying similarity between motifs , 2007, Genome Biology.

[16]  K. J. Fryxell,et al.  CpG mutation rates in the human genome are highly dependent on local GC content. , 2005, Molecular biology and evolution.

[17]  Michael B. Stadler,et al.  Distribution, silencing potential and evolutionary impact of promoter DNA methylation in the human genome , 2007, Nature Genetics.

[18]  K. Kinzler,et al.  Cancer Genome Landscapes , 2013, Science.

[19]  Laurent Duret,et al.  Biased gene conversion and the evolution of mammalian genomic landscapes. , 2009, Annual review of genomics and human genetics.

[20]  Jun Yu,et al.  Oncogenes without a Neighboring Tumor-Suppressor Gene Are More Prone to Amplification , 2017, Molecular biology and evolution.

[21]  B. Miotto,et al.  Emerging Concept in DNA Methylation: Role of Transcription Factors in Shaping DNA Methylation Patterns , 2015, Journal of cellular physiology.

[22]  J. Dufour,et al.  Tumor suppressor and hepatocellular carcinoma. , 2008, World journal of gastroenterology.

[23]  Jun Kawai,et al.  Heterotachy in Mammalian Promoter Evolution , 2006, PLoS genetics.

[24]  L. Duret,et al.  GC-content evolution in mammalian genomes: the biased gene conversion hypothesis. , 2001, Genetics.

[25]  Yong Wang,et al.  An evaluation of new criteria for CpG islands in the human genome as gene markers , 2004, Bioinform..

[26]  A. Riggs,et al.  Epigenetic Silencing May Aid Evolution by Gene Duplication , 2003, Journal of Molecular Evolution.

[27]  Robert S. Illingworth,et al.  CpG islands influence chromatin structure via the CpG-binding protein Cfp1 , 2010, Nature.

[28]  J. Thomson,et al.  CpG island chromatin is shaped by recruitment of ZF-CxxC proteins. , 2013, Cold Spring Harbor perspectives in biology.

[29]  S. Yi,et al.  What Signatures Dominantly Associate with Gene Age? , 2016, Genome biology and evolution.

[30]  Laurence D. Hurst,et al.  The evolution of isochores , 2001, Nature Reviews Genetics.

[31]  Bradley E. Bernstein,et al.  GC-Rich Sequence Elements Recruit PRC2 in Mammalian ES Cells , 2010, PLoS genetics.

[32]  A. Galecki,et al.  Genetic Loci That Influence Cause of Death in a Heterogeneous Mouse Stock , 2004, The journals of gerontology. Series A, Biological sciences and medical sciences.

[33]  A. Bird,et al.  CpG islands and the regulation of transcription. , 2011, Genes & development.

[34]  Nansheng Chen,et al.  Using RepeatMasker to Identify Repetitive Elements in Genomic Sequences , 2009, Current protocols in bioinformatics.

[35]  Steven J. M. Jones,et al.  Comprehensive Characterization of Cancer Driver Genes and Mutations , 2018, Cell.

[36]  M. Bulyk,et al.  Polycomb-like proteins link the PRC2 complex to CpG islands , 2017, Nature.

[37]  Zhongming Zhao,et al.  TSGene 2.0: an updated literature-based knowledgebase for tumor suppressor genes , 2015, Nucleic Acids Res..

[38]  Doron Lancet,et al.  GeneHancer: genome-wide integration of enhancers and target genes in GeneCards , 2017, Database J. Biol. Databases Curation.

[39]  J. Böhme,et al.  Gene conversion of major histocompatibility complex genes is associated with CpG-rich regions , 1999, Immunogenetics.

[40]  F. Marshall.  Promoter Hypermethylation of Tumor Suppressor Genes in Urine From Kidney Cancer Patients , 2004 .

[41]  Simon Kasif,et al.  Genomewide Analysis of PRC1 and PRC2 Occupancy Identifies Two Classes of Bivalent Domains , 2008, PLoS genetics.

[42]  P. Park,et al.  CpG Islands Recruit a Histone H3 Lysine 36 Demethylase , 2010, Molecular cell.

[43]  H. Ellegren,et al.  Evolutionary Consequences of DNA Methylation on the GC Content in Vertebrate Genomes , 2015, G3: Genes, Genomes, Genetics.

[44]  A. Cheng,et al.  Oncogenes expand during evolution to withstand somatic amplification , 2018, Annals of oncology : official journal of the European Society for Medical Oncology.