Yin and Yang of disease genes and death genes between reciprocally scale-free biological networks

Biological networks often show a scale-free topology with node degree following a power-law distribution. Lethal genes tend to form functional hubs, whereas non-lethal disease genes are located at the periphery. Uni-dimensional analyses, however, are flawed. We created and investigated two distinct scale-free networks; a protein–protein interaction (PPI) and a perturbation sensitivity network (PSN). The hubs of both networks exhibit a low molecular evolutionary rate (P < 8 × 10−12, P < 2 × 10−4) and a high codon adaptation index (P < 2 × 10−16, P < 2 × 10−8), indicating that both hubs have been shaped under high evolutionary selective pressure. Moreover, the topologies of PPI and PSN are inversely proportional: hubs of PPI tend to be located at the periphery of PSN and vice versa. PPI hubs are highly enriched with lethal genes but not with disease genes, whereas PSN hubs are highly enriched with disease genes and drug targets but not with lethal genes. PPI hub genes are enriched with essential cellular processes, but PSN hub genes are enriched with environmental interaction processes, having more TATA boxes and transcription factor binding sites. It is concluded that biological systems may balance internal growth signaling and external stress signaling by unifying the two opposite scale-free networks that are seemingly opposite to each other but work in concert between death and disease.

[1]  Jung Kyoon Choi,et al.  Epigenetic regulation and the variability of gene expression , 2008, Nature Genetics.

[2]  A. Barabasi,et al.  The human disease network , 2007, Proceedings of the National Academy of Sciences.

[3]  A. Barabasi,et al.  Lethality and centrality in protein networks , 2001, Nature.

[4]  Ju Han Kim,et al.  Genomic characterization of perturbation sensitivity , 2007, ISMB/ECCB.

[5]  Yves A. Lussier,et al.  Network models of genome-wide association studies uncover the topological centrality of protein interactions in complex diseases , 2013, J. Am. Medical Informatics Assoc..

[6]  L. Hurst The Ka/Ks ratio: diagnosing the form of sequence evolution. , 2002, Trends in genetics : TIG.

[7]  B. Pugh,et al.  Identification and Distinct Regulation of Yeast TATA Box-Containing Genes , 2004, Cell.

[8]  A. E. Hirsh,et al.  Protein dispensability and rate of evolution , 2001, Nature.

[9]  M. Gouy,et al.  Codon catalog usage and the genome hypothesis. , 1980, Nucleic acids research.

[10]  J. Bähler,et al.  Meta-analysis of genome regulation and expression variability across hundreds of environmental and genetic perturbations in fission yeast. , 2010, Molecular bioSystems.

[11]  Dmitrij Frishman,et al.  Complex principal component and correlation structure of 16 yeast genomic variables. , 2011, Molecular biology and evolution.

[12]  A. E. Hirsh,et al.  Functional genomic analysis of the rates of protein evolution. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[13]  D. Kültz,et al.  Molecular and evolutionary basis of the cellular stress response. , 2005, Annual review of physiology.

[14]  Stefan Wuchty,et al.  Peeling the yeast protein network , 2005, Proteomics.

[15]  A. E. Hirsh,et al.  Evolutionary Rate in the Protein Interaction Network , 2002, Science.

[16]  Jaak Vilo,et al.  Building and analysing genome-wide gene disruption networks , 2002, ECCB.

[17]  Yudong D. He,et al.  Functional Discovery via a Compendium of Expression Profiles , 2000, Cell.

[18]  Scott A. Rifkin,et al.  Genetic Properties Influencing the Evolvability of Gene Expression , 2007, Science.

[19]  A. Barabasi,et al.  Network biology: understanding the cell's functional organization , 2004, Nature Reviews Genetics.

[20]  J. Bähler,et al.  Tuning gene expression to changing environments: from rapid responses to evolutionary adaptation , 2008, Nature Reviews Genetics.

[21]  M. DePamphilis,et al.  HUMAN DISEASE , 1957, The Ulster Medical Journal.

[22]  Alan F. Scott,et al.  McKusick's Online Mendelian Inheritance in Man (OMIM®) , 2008, Nucleic Acids Res..

[23]  G Humphreys,et al.  Codon usage can affect efficiency of translation of genes in Escherichia coli. , 1984, Nucleic acids research.

[24]  M. Kimura Evolutionary Rate at the Molecular Level , 1968, Nature.

[25]  Yongjin Li,et al.  Discovering disease-genes by topological features in human protein-protein interaction network , 2006, Bioinform..

[26]  G. Hong,et al.  Nucleic Acids Research , 2015, Nucleic Acids Research.

[27]  Massimo Marchiori,et al.  Error and attacktolerance of complex network s , 2004 .

[28]  Albert-László Barabási,et al.  Error and attack tolerance of complex networks , 2000, Nature.

[29]  P. Sharp,et al.  The codon Adaptation Index--a measure of directional synonymous codon usage bias, and its potential applications. , 1987, Nucleic acids research.