- Open Access
Investigation of anti-cancer mechanisms by comparative analysis of naked mole rat and rat
BMC Systems Biology volume 7, Article number: S5 (2013)
The naked mole rats (NMRs) are small-sized underground rodents with plenty of unusual traits. Their life expectancy can be up to thirty years, more than seven times longer than laboratory rat. Furthermore, they are resistant to both congenital and experimentally induced cancer genesis. These peculiar physiological and pathological characteristics allow them to become a suitable model for cancer and aging research.
In this paper, we carried out a genome-wide comparative analysis of rat and NMR using the recently published genome sequence of NMR. First, we identified all the rat-NMR orthologous genes and specific genes within each of them. The expanded and contracted numbers of protein families in NMR were also analyzed when compared to rat. Seven cancer-related protein families appeared to be significantly expanded, whereas several receptor families were found to be contracted in NMR. We then chose those rat genes that were inexistent in NMR and adopted KEGG pathway database to investigate the metabolic processes in which their proteins may be involved. These genes were significantly enriched in two rat cancer pathways, "Pathway in cancer" and "Bladder cancer". In the rat "Pathway in cancer", 9 out of 14 paths leading to evading apoptosis appeared to be affected in NMR. In addition, a significant number of other NMR-missing genes enriched in several cancer-related pathways have been known to be related to a variety of cancers, implying that many of them may be also related to tumorigenesis in mammals. Finally, investigation of sequence variations among orthologous proteins between rat and NMR revealed that significant fragment insertions/deletions within important functional domains were present in some NMR proteins, which might lead to expressional and/or functional changes of these genes in different species.
Overall, this study provides insights into understanding the possible anti-cancer mechanisms of NMR as well as searching for new cancer-related candidate genes.
The naked mole rats (NMRs, Heterocephalus glaber) are mouse-sized subterranean rodents native to East Africa . They have an exceptional set of physiological traits that make them adapt to living in the underground of droughty desert. They are becoming one of the most extraordinary organisms known to science .
NMRs are the longest-lived rodent known till now and their maximum lifespan can be up to thirty years . By contrast, other similar-sized rodents such as mouse possess a life expectancy of only four years, which is far less than that of NMR. Previously published studies have indicated that the longevity of NMR was possibly because of the negligible decrement of age-related physiological characteristic along with their lifetime, such as declining fertility and mortality rate .
Besides delayed senescence, NMRs are remarkably resistant to both congenital and experimentally induced cancer genesis . Cancer is a group of complex polygenic diseases that commonly affect lots of vertebrates, and is constantly considered to be an inevitable accompanied by senescence. Cancer is the second dominant cause of mortality in the world, which cause 7.6 million of death estimated by World Health Organization . It has been recognized for quite a long time that cancer genesis is closely related to tumour suppressor genes and oncogenes . Identifying the function of the added genes may bring us in another way to explore the regulatory network of the cancer process. In addition, the mechanisms of cancer resistance present in NMR are not thoroughly clear. Thus, identification of NMR genes closely implicated in cancer may provide us effective clues for delineating causes of cancer proneness and studying anti-cancer properties for mammalian organisms.
NMRs possess several other special physiological characteristics as well. Although NMRs belong to the order of rodentia, they are actually poikilothermal animals whose body temperatures vary continuously following the environment . Furthermore, NMRs are insensitive to certain types of pain  and acid , and are well adapted to the underground surrounding at an extremely low oxygen concentration (10%-15%) .
Recently, using high-throughput next generation sequencing techniques, the genome of NMR has been sequenced. These excellent resources provide great opportunities for understanding the exceptional characteristics of NMR and improve biological and biomedical studies. Previously, some genes have been identified to be related to some of these unusual characteristics, e.g., the telomerase reverse transcriptase (TERT) gene and some other genes, which may be involved in extended longevity mechanisms of NMR . However, investigation of the genomic information of NMR at a systems-biology level is still lacking, which may provide additional information to uncover the molecular mechanisms for the extraordinary traits (e.g., anti-cancer) of NMR.
In this paper, a comparative genomics study was carried out to explore the genes that were either common between rat and NMR, or specific to each of them. We divided these genes into three groups: common genes, genes only present in NMR and genes only present in rat. We then used the Pfam database to identify the events of gain or loss of different protein families between these two species. In addition, the Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used to study the rat pathways in which the NMR-missing genes participate. A significant number of these genes were found to be enriched in the pathways related to the exceptional characteristics of NMR (such as cancer pathways), many of which have been previously reported to be associated with various cancers. Finally, we analyzed the sequence variations (such as domain insertion/deletion) of orthologous proteins to investigate the potentially expressional and/or functional alternations of them between rat and NMR. Overall, our data not only help unveil the cancer resistance mechanisms of NMR but provide insights into identifying new cancer-related genes.
Genome, database and resources
The complete set of annotated rat and NMR protein sequences were obtained from the UniProt database (http://www.uniprot.org/). For those genes with alternative splicing variants, proteins with the smallest PE value (which means the most possibility for the existence of the proteins) and the longest length were chosen to represent the gene-encoding protein sequences. A total of 20835 and 21553 proteins corresponding to their genes were finally obtained for rat and NMR, respectively.
The file containing the whole pathways of rat was downloaded from the KEGG database (http://www.genome.jp/kegg/pathway.html). The Online Mendelian Inheritance in Man (OMIM) database (http://www.ncbi.nlm.nih.gov/omim) was used to analyze the relationship between cancers and human orthologs of genes absent in NMR. Furthermore, the expression data of these genes in rat tissues were obtained from the Gene Expression Atlas (GXA) database (http://www.ebi.ac.uk/gxa/), which was used to identify whether or not the expression levels of these genes were related to cancer development.
Identification of orthologous genes between rat and NMR
To analyze the orthologous gene pairs between rat and NMR, we employed the complete set of annotated proteins of one organism as queries to search for orthologs in the other species via BLASTP with a cut-off of E-value ≤ 1e-6. Orthologous genes were further defined as bidirectional best hits.
On the basis of identified orthologs between rat and NMR, we dissected all the genes into three classifications:
Class I: the Shared genes, which were shared between rat and NMR;
Class II: the NMR-missing genes, which were absent in NMR but present in rat;
Class III: the NMR unique genes, which were found in NMR but missing in rat.
Considering that the conserved domains in a protein could provide information for its function and evolutionary dynamics, we used the Pfam database  (http://pfam.sanger.ac.uk/search), which collected a large collection of protein families, to search for gain or loss events of different protein families between these two species.
All the proteins of rat and NMR were searched against Pfam database with a cut-off of E-value ≤ 1e-5. For each protein, if two or more Pfam families were available, only the one with the smallest E-value was selected. The number of each protein family in rat and NMR was then calculated respectively.
KEGG pathway analysis of Class II genes
We further dissected the pathways containing Class II genes using the KEGG database resource, which is a collection of manually curated pathway maps according to current knowledge on protein-protein interactions .
First, each gene of Class II was mapped to their pathways. The p-value of the enrichment of NMR-missing genes in each pathway was then calculated by hypergeometric distribution test. Moreover, considering that KEGG pathways were composed of nodes which were actually modules including single gene or multiple functionally similar genes, we further analyzed the number and percentage of the nodes containing Class II genes in each of the enriched pathways.
Sequence variation analysis of Class I genes
Although Class I genes were considered as orthologous genes between rat and NMR, sequence variations had been previously observed for certain proteins. For example, the glutathione peroxidase 1 (GPx1), which is highly expressed in mouse liver and kidney, has an early stop codon in NMR. Such a variation results in the lack of the C-terminal part and may be related to an order of magnitude lower activity in NMR tissues . Thus, it would be useful to further study the orthologous genes between the two organisms for potential changes with regard to their function and/or regulation.
To systematically investigate such deletion or insertion events, we analyzed the BLASTP sequence alignment results for each Class I gene in NMR, and focused on gap-related parameters, such as alignment length, number of mismatches, and percentage of identical matches, to calculate the lengths of sequence insertions/deletions. Protein pairs of rat and NMR proteins were chosen at a cut-off of gap length>25 and percentage of mismatches <10%. To avoid incorrect protein annotation of NMR, the NMR proteins with significant fragment deletion were searched against the NMR genome using TBLASTN for further verification of the absence of these segments. Finally, all selected proteins were searched against the Conserved Domains Database (http://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml) to identify their functional domains.
Distribution of orthologous genes in rat and NMR
Figure 1 shows the Venn diagram between rat and NMR genes. A set of 15408 genes (Class I) were found in both rat and NMR, which occupied 73.9% of the total predicted NMR genes. On the other hand, 6145 rat genes (Class II) lacked orthologs in NMR. In contrast, 5427 genes (Class III) appeared to be unique in NMR. Nevertheless, considering that the number of NMR orthologous genes(Class I) is only about 71.5% of that of the rat total genes, it is possible that some genes may be mis-annotated in the currently released version of NMR genome.
Pfam analysis of rat and NMR genes
All the annotated proteins of rat and NMR were searched against Pfam database (containing 13672 families) for the classification into different protein families. 2442 and 2523 protein families (with 2416 overlapped families) were obtained in rat and NMR, respectively, indicating that the two species shared almost the same protein families.
We further analyzed the gain or loss events for each protein family in NMR. Table 1 shows the top 15 families sorted by the number of expanded genes (gain) in NMR. Among them, two distinct cancer-related protein families: "Melanoma-associated antigen family" and "Mortality factor 4 families" had significant gene expansion in NMR. A total of 56 and 15 genes were found for the two families respectively in NMR, which were 27 and 11 more than those in rat. Besides, several other families shown in Table 1 have also been reported to be related to cancer, such as the protein kinase C (PKC) family and several heat shock proteins (HSP70 and HSP90).
The PKC family, which possessed 53 proteins in NMR, had variable roles in tumour biology depending on the intracellular localizations and cell types. PKCs were generally abnormally regulated in the cancers of the breast, prostate, kidney and liver , and remained as a possible target for cancer prevention and therapy . Here, a total of 27 additional PKC members were identified in NMR, implying that these new PKCs may play an important role in preventing NMR from cancer.
HSP proteins are a group of functionally related proteins regulating protein folding and unfolding reactions. HSP70 proteins were reported to be overexpressed in the malignant melanoma . On the other hand, HSP90 proteins were also implicated to be involved in breast cancer progression because of its overexpression in breast cancer cell lines and association with survival of breast cancer . Thus, HSP70 and HSP90 proteins have been considered as the useful targets for cancer therapy [20, 21]. In this study, the significant expansions of members of these two families in NMR were consistent with their potential roles in cancer prevention and may provide clues for the anti-cancer trait of NMR.
Table 2 shows the top 15 protein families with contracted gene numbers (loss) in NMR. Among them, several protein families have been previously verified to be involved in cancer development. For example, the cadherin family is a class of type-1 trans-membrane proteins. Members of this family may play important roles in cell-cell adhesion in different organ systems . All kinds of abnormal expression of cadherin family proteins had been reported as a widespread phenomenon in mammary cancer and these proteins had been frequently implicated in tumour progression . The loss of so many genes in this family may also contribute to the complex cancer-resistant mechanisms of NMR. In addition, a significant loss of members of some receptor protein families, such as the vomeronasal organ receptor (VOR), olfactory receptor and class C G-protein-coupled receptor, may be important for some unique features of the most unusual rodent.
Pathway analysis of NMR-missing genes
Among 216 rat pathways in the KEGG database, 16 of them were found to be significantly (p-value < 0.01) enriched by the NMR-missing genes (group II genes) (Table 3). For example, at least 31 rat genes (10.7%) in the "Neuroactive ligand-receptor interaction" pathway were not detected in NMR. The loss of these genes might be one of the major reasons why NMR is insensitive to pain and acid . In addition, NMR-missing genes were also significantly enriched in two known cancer pathways, "Bladder cancer" and "Pathways in cancer" (with the p-value 2.37e-3 and 1.17e-3, respectively).
Among the rest of pathways shown in Table 3 five of them were thought to be cancer-related, which include "Cytokine-cytokine receptor interaction" , "p53 signalling pathway", "Apoptosis", "Natural killer cell mediated cytotoxicity" , "Wnt signalling pathway" " and "Notch signalling pathway" . It is possible that several of these NMR-missing genes are associated with cancer development in other mammals including humans and could be considered as candidate cancer-related genes.
As KEGG pathways are composed of nodes which may have single or multiple functionally similar genes, we also calculated the percentage of the nodes which contain at least one NMR-missing gene in each of these pathways, and obtained almost the same enriched pathways, such as "Pathways in cancer", "Neuroactive ligand-receptor interaction" and "Oxidative phosphorylation", implying that the pathway enrichment of NMR-missing genes was significant at both gene and node levels (Supplementary Table 1 [see additional file 1]).
Analysis of cancer-related genes that were absent in NMR
To investigate the potential mechanisms of the anti-cancer aspects of NMR, three cancer-related pathways, including "pathways in cancer", "MAPK (mitogen-activated protein kinase) signalling pathway" and "Wnt signalling pathway", were chosen as examples for further analysis.
1) Pathways in cancer
Twenty-nine genes in this pathway were not detected in NMR (Supplementary Table 2 [see additional file 1]). These genes were found to be strongly related to cancer. Half of them correspond to various phenotypes of cancer (e.g., leukemia, lung cancer and adrenal cortical carcinoma) based on the information retrieved from OMIM database, including several well-studied carcinogenesis genes (Bcl2, Casp8, Fas). Moreover, some important proto-oncogenes, such as Myc and Hras1, were also absent. In fact, it is well known that proto-oncogenes are normal genes that could become the oncogenes because of their overexpression or mutations. The loss of the proto-oncogenes in NMR cells may also contribute to cancer resistance.
Among all 29 NMR-missing genes in this pathway, 19 of them (65.5%) were previously reported to display differential gene expression levels between cancer and normal tissues. Thus, the absence of these genes might play important roles in suppressing cancer.
Figure 2 shows the main part of the map of "pathways in cancer". One of the key mechanisms for cancer development is "Evading apoptosis", which is a crucial oncogenic property of cancer cells. Among the 14 paths leading to evasion of apoptosis illustrated in this map, 9 of them appeared to be influenced by direct connection to the nodes which contain NMR-missing genes. Particularly, the nodes of Survivin, Bcl-2, Bcl-XL and Mtor have lost several functionally similar genes in NMR. It is possible that these changes in NMR might block the major routes to evade apoptosis, and then induce the early programmed death of cancer cells as an important cancer barrier.
Recently, a two-tier anti-cancer mechanism associated with contact inhibition regulated by p16Ink4a and p27Kip1 has been reported in NMR . However, rat cells were found to only have contact inhibition regulated by p27Kip1. This is consistent with our results as rats only have p27Kip1 gene. Thus, NMR appeared to have additional unique protective mechanisms for cancer resistance.
2) MAPK signalling pathway
The MAPK signalling pathway is an important pathway how proteins in the cytoplasm communicate the signals from the receptors on the cell membrane to the nucleus. It is in the central of a molecular metabolic network that mediates cell differentiation and proliferation. In mammalian cells, the MAPK pathway contains three major groups of proteins, including Erk (extracellular signal-regulated kinase), p38 kinase and Sapk (stress activated protein kinase). These proteins are abnormally regulated in various diseases, including cancer and inflammation.
In this study, 19 genes were found to be absent in NMR (Figure 3 and Supplementary Table 3 [see additional file 1]). By retrieving the OMIM database and GXA database, we found that 7 (36.8%) of them have been reported to be related to various phenotypes of cancer in other mammals and 10 (52.6%) of them had differential expression levels between cancer and normal tissues.
In this pathway, three proto-oncogenes, Hras1, Myc and Pdgfb (also present in the "Pathways in cancer" pathway), were absent in NMR. It has been previously shown that when one of these genes was mutated, the activity of their enzymes could be stuck in the "on" or "off" position, which was an essential step during the development of many cancers . Recently, it has also been reported that NMR and rat cells acted totally opposed if transfected with Hras1. NMR cell cycle came to an abrupt end as the presence of abnormal chromatin material and anaphase bridges and, while transfected rat cell grew rapidly and formed tumours eventually . Therefore, the loss of these genes might also play a significant role in cancer resistance.
Although the phenotypes of cancer could not be found for other NMR-missing genes, some of them have been demonstrated to be related to the survival of cancer cells, such as Jund and Park genes. Jund is an AP-1 family member involved in various biological processes such as cell apoptosis and tumour metastasis, and could regulate survival of tumour cells in prostate cancer . Prak is a protein kinase, which was previously shown to be implicated in the suppression of skin carcinogenesis . Further experiments are needed to investigate the relationship between these genes and tumorigenesis.
3) Wnt signalling pathway
The "Wnt signalling pathway" is a conserved protein-protein interaction network that regulates cell fate decisions and cell-cell communication. This pathway plays a significant role in maintaining stability of internal environment by regulating cell niche in vivo. Abnormal regulation of this pathway could lead to neoplastic proliferation which is involved in the progress of cancer cells.
Seventeen genes were found to be absent in the NMR genome (Figure 4 and Supplementary Table 4 [see additional file 1]). Among them, 4 (25%) genes have corresponding phenotypes of cancer whereas 12 (75%) genes were differentially expressed between cancer and normal tissues. Molecular-level studies have indicated that approximately 90% of the activating mutations of genes in this pathway could cause different cancers such as colorectal cancer .
Several well-studied cancer-related genes, such as Myc, Rhoa, Lef1 and Rac1, were absent in this pathway of NMR. Myc is a well-known proto-oncogene and has been frequently used to induce tumour formation in a lot of animal experiments of cancer research. Rhoa has been deeply studied and proved as a cancer-regulated gene, which controlled metastasis of tumour cells, acted as a regulator of male hormone activity in prostate cancer cells , and triggered a particular microvesicle signalling pathway in cancer cells . Lef1 protein could interact with a number of other proteins, such as Ctbp and Nlk. These interactions were thought to be responsible for the invasion and growth of prostate cancer . Rac1 was found to be associated with DNA transcription. Previous studies have reported that activation of Rac1 mediated Twist1-induced cancer cell migration . On the other hand, 12 of the NMR-missing genes in this pathway, such as Dkk4, Sox17 and Ccnd3, have not been reported to be related to any disease including cancer. It is possible that some of them are also involved in cancer formation and could be further experimentally verified.
Sequence variation analysis of Class I genes
Based on the BLASTP sequence alignment results, we found that 6349 (41.2%) of the orthologs had no gaps and 12142 (78.8%) orthologs only possessed a small (≤5 amino acids) gaps, suggesting that most of the orthologous proteins between rat and NMR had rare insertions/deletions during evolution. On the other hand, 439 orthologous proteins showed significant segment insertions/deletions whose length was more than 25 amino acids. Further analysis of these inserted and deleted fragments in these proteins revealed that many of them contained conserved sites, including functionally active sites (Supplementary Table 5 [see additional file 1]). For example, we observed that parts of the specific RNA/DNA binding site and the specific cytokine receptor motif were deleted in the NMR Fusip1 and Nrcam proteins, respectively.
Other domains affected by the insertion/deletion of certain segments included ATP-binding, Ca2+ binding and some other metal catalytic binding sites. For example, a 30-amino-acid-long sequence fragment was found to be inserted into the putative catalytic site in NMR Ship1 when compared with its ortholog in rat. It has been previously demonstrated that the phosphate domain of Ship1 was essential for catalytic activity in vivo  and the loss of Ship1 could promote leukemogenesis in a virus-infected mouse model . We suspect the insertion of such a long segment of Ship1 would change the function or expression of this gene. On the other hand, wrong annotation of some NMR proteins could not be excluded. Further studies are required to verify the presence of the sequence variations and their influence on the regulation or function of these proteins.
In this paper, a comparative genomics study was carried out to investigate the genes that were either common between rat and NMR, or specific to each of them. The majority of genes were shared by the two rodents, whereas each organism had a significant part of unique genes. Seven cancer-related protein families, such as melanoma-associated antigen family, protein kinase C family and HSP family were found to be significantly expanded. Further analysis of the genes absent in NMR indicated that the majority of them have been shown to be linked to many forms of cancer. Finally, some conserved functional domains were found to be possibly influenced by the insertion or deletion of certain fragments in NMR, which may change the expression or function of some of these genes. These results may provide important clues about the molecular mechanisms of cancer resistance of NMR and help identify new cancer-related genes  in mammals. As future topics, it is important to study such complex mechanisms from the viewpoints of network [43–46] and dynamics [47, 48] by further incorporating the expression data.
naked mole rat
Kyoto Encyclopaedia of Genes and Genomes
National Center for Biotechnology Information
Gene Expression Atlas
Conserved Domain Database
Online Mendelian Inheritance in Man
protein kinase C
heat shock proteins
telomerase reverse transcriptase
glutathione peroxidase 1
vomeronasal organ receptor
van der Horst G, Maree L, et al.: Sperm structure and motility in the eusocial naked mole-rat, Heterocephalus glaber: a case of degenerative orthogenesis in the absence of sperm competition. BMC Evol Biol. 2011, 11: 351-10.1186/1471-2148-11-351.
Buffenstein R: The naked mole-rat: a new long-living model for human aging research. J Gerontol A Biol Sci Med Sci. 2005, 60 (11): 1369-1377. 10.1093/gerona/60.11.1369.
Deacon RM, Dulu TD, Patel NB: Naked mole-rats: Behavioural phenotyping and comparison with C57BL/6 mice. Behav Brain Res. 2012, 20 231 (1): 193-200.
Buffenstein R: Negligible senescence in the longest living rodent, the naked mole-rat: insights from a successfully aging species. J Comp Physiol B. 2008, 178 (4): 439-45. 10.1007/s00360-007-0237-5.
Edrey YH, Hanes M, et al.: Successful aging and sustained good health in the naked mole rat: a long-lived mammalian model for biogerontology and biomedical research. ILAR J. 2011, 52 (1): 41-53. 10.1093/ilar.52.1.41.
World Health Organization: The Global Burden of Disease: 2004 Update. 2008, Geneva: World Health Organization
Helmig S, Schneider J: Oncogene and tumor-suppressor gene products as serum biomarkers in occupational-derived lung cancer. Expert Rev Mol Diagn. 2007, 7 (5): 555-568. 10.1586/14737184.108.40.2065.
Woodley R, Buffenstein R: Thermogenic changes with chronic cold exposure in the naked mole-rat (Heterocephalus glaber). Comp Biochem Physiol A Mol Integr Physiol. 2002, 133 (3): 827-834. 10.1016/S1095-6433(02)00199-X.
Park TJ, et al.: Selective inflammatory pain insensitivity in the African naked mole-rat (Heterocephalus glaber). PLoS Biol. 2008, 6 (1): e13-10.1371/journal.pbio.0060013.
Edrey YH, Park TJ, Kang H, Biney A, Buffenstein R: Endocrine function and neurobiology of the longest-living rodent, the naked mole-rat. Exp Gerontol. 2011, 46: 116-123. 10.1016/j.exger.2010.09.005.
Peterson BL, Larson J, Buffenstein R, Park TJ, Fall CP: Blunted Neuronal Calcium Response to Hypoxia in Naked Mole-Rat Hippocampus. PLoS ONE. 2012, 7 (2): e31568-10.1371/journal.pone.0031568. doi:10.1371/journal.pone.0031568
Kim EB, et al.: Genome sequencing reveals insights into physiology and longevity of the naked mole rat. Nature. 2011, 479 (7372): 223-7. 10.1038/nature10533.
Punta M, Coggill PC, et al.: The Pfam protein families database. Nucleic Acids Res. 2012, 40 (Database issue): D290-301.
Jitao Z, Stefan W: KEGGgraph: a graph approach to KEGG PATHWAY in R and bioconductor. Bioinformatics. 2009, 25 (11): 1470-1471. 10.1093/bioinformatics/btp167.
Kasaikina MV, Lobanov AV, Malinouski MY, Lee BC, Seravalli J, Fomenko DE, Turanov AA, Finney L, Vogt S, Park TJ, Miller RA, Hatfield DL, Gladyshev VN: Reduced utilization of selenium by naked mole rats due to a specific defect in GPx1 expression. J Biol Chem. 2011, 286 (19): 17005-14. 10.1074/jbc.M110.216267.
Griner EM, Kazanietz MG: Protein kinase C and other diacylglycerol effectors in cancer. Nat Rev Cancer. 2007, 7 (4): 281-294. 10.1038/nrc2110.
Ali AS, Ali S, et al.: Exploitation of protein kinase C: a useful target for cancer therapy. Cancer Treat Rev. 2009, 35 (1): 1-8. 10.1016/j.ctrv.2008.07.006.
Ricaniadis N, Kataki A, et al.: Long-term prognostic significance of HSP-70, c-myc and HLA-DR expression in patients with malignant melanoma. Eur J Surg Oncol. 2001, 27 (1): 88-93. 10.1053/ejso.1999.1018.
Pick E, Kluger Y, et al.: High HSP90 expression is associated with decreased survival in breast cancer. Cancer Res. 2007, 67 (7): 2932-2937. 10.1158/0008-5472.CAN-06-4511.
Galluzzi L, Giordanetto F, et al.: Targeting HSP70 for cancer therapy. Mol Cell. 2009, 36 (2): 176-177. 10.1016/j.molcel.2009.10.003.
Mahalingam D, Swords R, et al.: Targeting HSP90 for cancer therapy. Br J Cancer. 2009, 100 (10): 1523-1529. 10.1038/sj.bjc.6605066.
Berx G, van Roy F: Involvement of members of the cadherin superfamily in cancer. Cold Spring Harb Perspect Biol. 2009, 1 (6): a003129-10.1101/cshperspect.a003129.
Gama A, Schmitt F: Cadherin cell adhesion system in canine mammary cancer: a review. Vet Med Int. 2012, 2012: 357187-
Smith ES, Omerbasic D, et al.: The molecular basis of acid insensitivity in the African naked mole-rat. Science. 2011, 334 (6062): 1557-1560. 10.1126/science.1213760.
Van Dyke AL, Cote ML, et al.: Cytokine and cytokine receptor single-nucleotide polymorphisms predict risk for non-small cell lung cancer among women. Cancer Epidemiol Biomarkers Prev. 2009, 18 (6): 1829-1840. 10.1158/1055-9965.EPI-08-0962.
Mirzayans R, Andrais B, et al.: New Insights into p53 Signalling and Cancer Cell Response to DNA Damage: Implications for Cancer Therapy. J Biomed Biotechnol. 2012, 2012: 170325-
Hu W, Zheng RR, et al.: Effects of bortezomib in sensitizing human prostate cancer cell lines to NK-mediated cytotoxicity. Asian J Androl. 2012, 14 (5): 695-702. 10.1038/aja.2012.59.
Sastre-Perona A, Santisteban P: Role of the wnt pathway in thyroid cancer. Front Endocrinol (Lausanne). 2012, 3: 31-
Xia J, Li Y, et al.: Arsenic Trioxide Inhibits Cell Growth and Induces Apoptosis through Inactivation of Notch Signalling Pathway in Breast Cancer. Int J Mol Sci. 2012, 13 (8): 9627-9641.
Liang S, Mele J, et al.: Resistance to experimental tumorigenesis in cells of a long-lived mammal, the naked mole-rat (Heterocephalus glaber). Aging Cell. 2010, 9 (4): 626-635. 10.1111/j.1474-9726.2010.00588.x.
Martínez Q, Rodríguez G, et al.: LEOPARD Syndrome: Clinical Features and Gene Mutations. Mol Syndromol. 2012, 3: 145-157.
Seluanov A, Hine C, Azpurua J, Feigenson M, Bozzella M, Mao ZY, et al.: Hypersensitivity to contact inhibition provides a clue to cancer resistance of naked mole-rat. Proc Natl Acad Sci USA. 2009, 106: 19352-19357. 10.1073/pnas.0905252106.
Zerbini LF, de Vasconcellos JF, et al.: Jund-mediated repression of GADD45alpha and gamma regulates escape from cell death in prostate cancer. Cell Cycle. 2011, 10 (15): 2583-2591. 10.4161/cc.10.15.16057.
Yoshizuka N, Lai M, et al.: Prak suppresses oncogenic ras-induced hematopoietic cancer development by antagonizing the JNK pathway. Mol Cancer Res. 2012, 10 (6): 810-820. 10.1158/1541-7786.MCR-11-0576.
Giles RH, van Es JH, et al.: Caught up in a Wnt storm: Wnt signaling in cancer. Biochim Biophys Acta. 2003, 1653 (1): 1-24.
Schmidt LJ, Duncan K, et al.: RhoA as a mediator of clinically relevant androgen action in prostate cancer cells. Mol Endocrinol. 2012, 26 (5): 716-735. 10.1210/me.2011-1130.
Li B, Antonyak MA, et al.: RhoA triggers a specific signaling pathway that generates transforming microvesicles in cancer cells. Oncogene. 2012, doi: 10.1038/onc.2011.636
Li YR, Wang LG, et al.: LEF1 in Androgen-Independent Prostate Cancer: Regulation of Androgen Receptor Expression, Prostate Cancer Growth, and Invasion. Cancer Res. 2009, 69 (8): 3332-3338. 10.1158/0008-5472.CAN-08-3380.
Yang WH, Lan HY, et al.: RAC1 activation mediates Twist1-induced cancer cell migration. Nat Cell Biol. 2012, 14 (4): 366-374. 10.1038/ncb2455.
Conde C, Gloire G, et al.: Enzymatic and non-enzymatic activities of SHIP-1 in signal transduction and cancer. Biochem Pharmacol. 2011, 82 (10): 1320-1334. 10.1016/j.bcp.2011.05.031.
Lakhanpal GK, Vecchiarelli-Federico LM, Li YJ, Cui JW, Bailey ML, Spaner DE, et al.: The inositol phosphatase SHIP-1 is negatively regulated by Fli-1 and its loss accelerates leukemogenesis. Blood. 2010, 116: 428-36. 10.1182/blood-2009-10-250217.
Yang Z, Yan Z, Chen L: In silico identification of novel cancer-related genes by comparative genomics of naked mole rat and rat. Systems Biology (ISB), 2012 IEEE 6th International Conference. 2012, 285-290. 18-20 Aug. 2012, doi: 10.1109/ISB.2012.6314151
Song W, Wang J, Yang Y, Jing NH, Zhang XS, Chen L, Wu JR: Rewiring drug-activated p53-regulatory network from suppressing to promoting tumorigenesis. Journal of Molecular Cell Biology. 2012, 4 (4): 197-206. 10.1093/jmcb/mjs029.
He D, Liu Z, Honda M, Kaneko S, Chen L: Coexpression network analysis in chronic hepatitis B and C hepatic lesion reveals distinct patterns of disease progression to hepatocellular carcinoma. Journal of Molecular Cell Biology. 2012, 4 (3): 140-152. 10.1093/jmcb/mjs011.
Wen Z, Liu Z, Liu Z, Zhang Y, Chen L: An integrated approach to identify causal network modules of complex diseases with application to colorectal cancer. Journal of the American Medical Informatics Association. 2012, doi:10.1136/amiajnl-2012-001168
Wang J, Sun Y, Zheng S, Zhang X, Zhou H, Chen L: APG: an Active Protein-Gene Network Model to Quantify Regulatory Signals in Complex Biological Systems. Scientific Reports. 2013, DOI: 10.1038/srep01097
Chen L, Liu R, Liu Z, Li M, Aihara K: Detecting early-warning signals for sudden deterioration of complex diseases by dynamical network biomarkers. Scientific Reports. 2012, 2: 342-
Zhu H, Shyama Prasad Rao R, Liu J, Zeng T, Chen L: Reconstructing dynamic gene regulatory networks from sample-based transcriptional data. Nucleic Acids Research. 2012, doi:10.1093/nar/gks860
A preliminary version of this paper was published in the proceedings of IEEE ISB2012. We are also very grateful to KEGG database because of the pathway maps drawn by it.
The publication of this article is funded by grants from the National Natural Science Foundation of China (91029301, 61134013, 31171233, 61072149), a grant from Chinese Academy of Sciences (CAS) (2012OHTP10), the Chief Scientist Program of Shanghai Institutes for Biological Sciences of CAS (2009CSP002), and the FIRST program from JSPS initiated by CSTP.
This article has been published as part of BMC Systems Biology Volume 7 Supplement 2, 2013: Selected articles from The 6th International Conference of Computational Biology. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcsystbiol/supplements/7/S2.
The authors declare that they have no competing interests.
YZ and LC designed the study. ZY carried out computational studies, including comparative genomics, sequence alignment, Pfam and KEGG analysis and wrote the manuscript. YZ and LC analyzed the data and edited the manuscript. All authors read and approved the final manuscript.