Next Article in Journal
Evolution and Structure of Proteins and Proteomes
Next Article in Special Issue
Medicine in the Light of Evolution
Previous Article in Journal
The Emerging Role of Cohesin in the DNA Damage Response
Previous Article in Special Issue
Evolutionary Emergence of Drug Resistance in Candida Opportunistic Pathogens
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 9
Issue 12
10.3390/genes9120582
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Selective Pressures on Human Cancer Genes along the Evolution of Mammals

by
Alberto Vicens
1,2 and
David Posada
1,2,3,*
1
Department of Biochemistry, Genetics and Immunology, University of Vigo, 36310 Vigo, Spain
2
Biomedical Research Center (CINBIO), University of Vigo, 36310 Vigo, Spain
3
Galicia Sur Health Research Institute, 36310 Vigo, Spain
*
Author to whom correspondence should be addressed.
Genes 2018, 9(12), 582; https://doi.org/10.3390/genes9120582
Submission received: 30 October 2018 / Revised: 21 November 2018 / Accepted: 21 November 2018 / Published: 28 November 2018
(This article belongs to the Special Issue Evolutionary Medicine)
Browse Figures
Versions Notes

Abstract

:
Cancer is a disease driven by both somatic mutations that increase survival and proliferation of cell lineages and the evolution of genes associated with cancer risk in populations. Several genes associated with cancer in humans, hereafter cancer genes, show evidence of germline positive selection among species. Taking advantage of a large collection of mammalian genomes, we systematically looked for signatures of germline positive selection in 430 cancer genes available in COSMIC. We identified 40 cancer genes with a robust signal of positive selection in mammals. We found evidence for fewer selective constraints—higher number of non-synonymous substitutions per non-synonymous site to the number of synonymous substitutions per synonymous site (dN/dS)—and higher incidence of positive selection—more positively selected sites—in cancer genes bearing germline and recessive mutations that predispose to cancer. This finding suggests a potential association between relaxed selection, positive selection, and risk of hereditary cancer. On the other hand, we did not find significant differences in terms of tissue or gene type. Human cancer genes under germline positive selection in mammals are significantly enriched in the processes of DNA repair, with high presence of Fanconi anaemia/Breast Cancer A (FA/BRCA) pathway components and T cell proliferation genes. We also show that the inferred positively selected sites in the two genes with the strongest signal of positive selection, i.e., BRCA2 and PTPRC, are in regions of functional relevance, which could be relevant to cancer susceptibility.

1. Introduction

Cancer is a genomic disease caused by mutations in genes that control normal cell functions, in particular growth and division. A fundamental goal of cancer genomics is identifying mutations that confer a selective advantage to the cell and increase survival and proliferation, the so-called driver mutations, as well as the genes carrying the driver mutations in each tumor, known as driver genes or cancer genes. Although traditionally the focus has been put on somatic driver mutations, those that appear during an individual lifetime as cells divide and grow, there are also germline mutations in the human population that predispose to cancer [1]. Today, more than 500 human cancer genes have been identified, of which approximately 90% contain somatic mutations and 20% bear germline mutations [2,3,4]
Several studies have identified a number of human cancer genes undergoing germline positive selection among species. Clark et al. revealed a strong evidence of positive selection on oncogenes and tumor suppressor genes in the chimpanzee lineage [5]. Nielsen et al. identified an elevated number of tumor suppressor and apoptosis genes under strong positive selection in humans and chimpanzees [6]. Subsequent genome-wide screenings in mammals unveiled positive selection in genes with roles in immunity and reproduction, but also related with apoptosis and cancer [7,8]. Given the recurrent observation of positive selection acting on human cancer genes across species, several authors have proposed that evolutionary pressures affecting organismal fitness, such as sexual selection, pathogen–host interactions, parent–offspring conflict, or maternal–fetal conflict, could lead to increased cancer risk in humans as a pleiotropic effect [6,7,8,9,10]. Another intriguing matter is whether the germline evolution of human cancer genes is influenced by their different characteristics, like genetic type, main tissue, or inheritance mode. For example, in hominoids it was observed that tumor suppressor genes (TSGs) tend to accumulate more non-synonymous substitutions than oncogenes, suggesting that the former are subjected to more relaxed purifying selection due to their recessive effects [11].
To address these questions in more detail, we carried out a comprehensive analysis of the evolution of 430 human cancer genes in mammals, after applying stringent criteria for evolutionary analysis. Using the ratio of number of non-synonymous substitutions per non-synonymous site to the number of synonymous substitutions per synonymous site (dN/dS), we identified 40 human cancer genes under putative positive selection in mammals. These genes are functionally enriched in DNA repair and immunity and are associated with hereditary cancer and recessive effects.

2. Materials and Methods

2.1. Cancer Genes

We retrieved 574 single-copy genes associated with human cancer from the Cancer Gene Census (CGC) project [4] of the COSMIC repository (https://cancer.sanger.ac.uk/cosmic, accession date: 5 March 2018). We only collected genes classified into Tier 1, which refers to genes with a documented activity relevant to cancer. The catalogue of the retrieved cancer genes, along with information about their function and associated mutations, is shown in Table S1.

2.2. Sequence Data Collection

We downloaded a representative human DNA sequence for each cancer gene from the Ensembl Genes database (Release 91, human genome version GRCh38.p12) using BioMart (accession date: 8 March 2018). For each human gene, we chose a single isoform on the basis of the following criteria, in the indicated order: GENCODE validation, APRIS annotation as principal 1, best transcript support level (TSL), and longer transcript. We discarded 22 genes whose best TSL was less than 1 (Table S2). Using the selected human isoform as a reference, we downloaded the corresponding orthologues from 32 mammalian genomes (Table S3, Figure S1) using the Bioconductor package BiomaRt [12]. When more than one ortholog was obtained for a given species, we chose the isoform with the best orthology confidence score. We discarded 17 genes for which less than 15 mammal orthologues were found (Table S2), so the final number of retrieved ortholog groups was 535 (Table S4).

2.3. Multiple Sequence Alignment

We aligned the coding sequences for each ortholog group using MACSE [13], a program that accounts for frameshifts and stop codons. The resulting multiple sequence alignments were further refined with TrimAl [14], removing taxa and sites with more than 60% gaps across rows and columns, respectively. After trimming, we discarded 71 genes that contained less than 10 orthologues (Table S2), in order to maximize the statistical power for the selection analyses, ending up with a list of 464 genes.

2.4. Estimation of Phylogenetic Trees

We inferred maximum likelihood (ML) gene trees for the 464 genes using Randomized Acelerated Maximum Likelihood-Next generation (RAxML-NG). All reconstructions were performed using the general time reversible substitution model [15] with gamma-distributed rate variation among sites [16]. For each gene, we obtained 10 starting trees using randomized stepwise addition parsimony. We assessed nodal support using 100 bootstrap replicates [17]. To minimize the impact of estimation errors and incomplete lineage sorting in subsequent analyses, we discarded 27 genes whose estimated tree topologies were quite distinct (normalized Robinson–Foulds (RF) distance ≥ 0.6) from a species tree assembled for 19 mammals with a well-known phylogenetic position (Figure S1). We calculated the RF distances with ETE3-compare [18].

2.5. Codon-Based Selection Models

We estimated nonsynonymous (dN) and synonymous (dS) substitution rates using the program codeml of the PAML package v4.9c [19] for 437 mammal gene trees. Because dS saturation decreases the power of detecting positive selection in codon-based models [20], we further discarded seven genes with an estimated dS > 15 (Table S2). To estimate the global dN/dS ratios for each of the remaining 430 genes, we used the one-ratio (M0) model, which assumes the same dN/dS for all branches in the gene tree and across sites. To identify genes under putative positive selection, we carried out three different tests. First, we compared different site-models using two likelihood ratio tests (LRTs): M1a (neutral) versus M2a (selection), and M8 (beta selection) versus M8a (beta neutral) [21,22]. The resulting p-values were adjusted for multiple testing using the Benjamini–Hochberg procedure [23] with a family-wise significance level of 0.05. In addition, we also tested for evidence of episodic positive selection using BUSTED [24], as implemented in Hyphy [25]. In order to be very stringent, only genes inferred to be under positive selection by the three tests were finally considered to be positively selected genes (PSGs). For those genes in which the LRT was significant, we considered as positively selected sites (PSS) those with a Bayes Empirical Bayes (BEB) posterior probability > 0.95 of having a dN/dS > 1 under both M2a and M8 [26].

2.6. Gene Ontology Enrichment Analysis

To identify enriched Gene Ontology (GO) terms in the PSGs, we used GOrilla [27]. We compared the list of 40 PSGs with a background list of the 574 cancer genes from CGC. We searched for significant GO terms (p-value < 0.01) in the three available ontologies: biological process, cellular component, and molecular function.

2.7. Pathogenic Germline Mutations

We retrieved the list of pathogenic germline variants for Breast cancer type 2 susceptibility (BRCA2) protein from the study published by the TCGA PanCanAtlas Germline Working Group [28].

2.8. Comparison of dN/dS Ratios across COSMIC Categories

We compared the M0 dN/dS ratios obtained across four different CGC–COSMIC classifications: mutation type, inheritance, tissue type, and cancer role. To test for significant dN/dS and proportion of PSGs differences between and among categories, we performed ANOVAs (for multiple comparisons) and t-tests (for pairwise comparisons) using the ggpubr package for R [29]. We adjusted the p-value for multiple pairwise comparisons using the Benjamin–Hochberg procedure. To compare the proportion of genes under putative positive selection across groups, we applied the chi-squared test (p < 0.05) function (chisq.test) implemented in R.

3. Results

After multiple processing steps and stringent criteria (see Materials and Methods), we finally assessed the selective pressures along the mammal phylogeny on 430 human cancer genes. Multiple sequence alignments included 11–32 taxa and were 108–4984 nt long (Table S5).

3.1. Long-Term Selective Pressures on Human Cancer Genes

The mean dN/dS for the 430 cancer genes examined was 0.122. The LRTs among site-specific dN/dS models were significant for 56 (M2a) and 61 (M8) genes, while BUSTED was significant for 357 genes, in all three tests after correcting for multiple testing (p-adj < 0.05) (Table S5). Within these, 40 genes were identified with M2, M8, and BUSTED, and therefore considered to be PSGs (Table 1; Figure 1). All these genes showed at least one PSS under M2a or M8 (BEB > 0.95) (Table S5).

3.2. Comparison of Selection Estimates across Functional Categories

We compared global dN/dS values across COSMIC categories. Genes bearing only germline mutations (i.e., associated with hereditary cancer) showed significantly higher dN/dS estimates than genes with only somatic mutations (i.e., associated with sporadic cancer) or with both somatic and germline mutations (Figure 2A), mainly due to a significantly increase in dN (Figure S2). We also observed higher dN/dS values for cancer genes associated with recessive mutations than for cancer genes with dominant mutations (Figure 2B), again due to a significantly increase in dN (Figure S2). We noticed that these two mutational categories are not independent, as 33 out of 34 (97%) genes with germline mutations are associated with recessive inheritance. On the other hand, the global dN/dS estimates were not significantly different among tissue types (epithelial, lymphoid, mesenchymal, and others) (Figure 2C) or cancer role (fusion genes, oncogenes, and TSGs) (Figure 2D).
We also compared the proportion of PSGs across COSMIC categories, observing a significant increase in the germline (Figure 3A) and recessive categories (Figure 3B). We did not detect significant differences in the proportion of PSGs among tissue types (Figure 3C) or cancer role (Figure 3D).
Because a longer gene is more likely to have more spurious PSS than a shorter gene, we assessed whether the significant patterns just described were influenced by differences in protein length among categories. For this, we performed three statistical analyses. First, we compared protein length with global dN/dS, without detecting a significant correlation (Figure S3A). Second, we compared the protein length of PSGs and non-PSGs, again without observing significant differences (Figure S3B). Third, we compared the protein length among mutational categories (Figure S3C). Here, genes carrying germline mutations were not different from genes carrying only somatic mutations, but genes bearing both mutations were significantly larger; anyway, this did not interfere with our selection analyses. In addition, we contrasted the number of PSSs, normalized by sequence length, across COSMIC categories for both M2a and M8 models. We did not find significant differences in the proportion of PSSs for any functional category, regardless of the model (Figure S4).

3.3. Functional Enrichment of Positively Selected Cancer Genes

The 40 PSGs were enriched in biological processes associated with DNA repair (p-value = 2.63 × 10−4) and the regulation of T cell proliferation (p-value = 4.1 × 10−4). We did not identify significant enrichment of GO terms for molecular function or cellular component. Among DNA repair genes, we detected a high presence of members of the Fanconi Anemia Complementation Group (FANCC, FANCD2, FANCG, FANCA, and FANCE), which participate in homologous recombination DNA repair [30], and genes involved in double-strand break repair (BRCA2, BRIP1, BARD, BLM, CHEK2, and PALB2). In addition to the enrichment of genes involved in regulation of T cell proliferation, we observed a high proportion of PSGs related to immunity (Table 1).

3.4. Functional Relevance of Positively Selected Sites in Cancer Genes

To understand the functional relevance of PSSs, we mapped their location in the two genes with the highest number of PSS, namely, BRCA2 and PTPRC. BRCA2 is a gene involved in double-strand break repair whose deficiency leads to hereditary breast and ovarian cancer [30]. We identified 16 and 25 PSSs under M2a and M8, respectively (Table S5). These two sets were nested, so we mapped the 25 M8 residues on human BRCA2 protein scheme (Figure 4A). We found that three selected residues (positions 711, 748, and 800) are located in the binding region of nucleophosmin (NPM) [31]. Six PSSs (1158, 1646, 1708, 1913, 2035, and 2037) are distributed along the BRC repeats that bind to RAD51 [32]. Among these, residues 1646 and 1708 sit in the interaction region with the polymerase Eta (POLH) [33]. We noticed that, in position 1913, 11 out of 27 species, including humans, have a Cys residue, whereas eight species have a His. Because Cys and His residues are often involved in specific functions within protein structures [34], replacements in this position could be functionally relevant. Four PSSs (2530, 2572, 2574, and 2884) locate within the interaction region with SEM1, a gene involved in DNA damage repair and cell cycle progression [35]. Among these, residues 2530, 2572, and 2574 cluster in the helical subdomain that interacts with FANCD2 [36], a partner of the FA/BRCA complex (also a PSG). Residue 2884 sits in the first Oligosaccharide binding (OB domain. Several PSSs were found accumulated in a disordered segment of the C-terminal region (3363–408) with no documented activity. We observed some PSSs mapped in close proximity to natural variants predisposing to human cancer [28] (Figure 4A). The stop-gained variant Y792*, associated with pancreatic adenocarcinoma (PAAD), is close to three PSSs. The PSS 1158 is in close proximity to the stop-gained variant Q1037*, also associated with PAAD. The selected residue 1646 is flanked by the stop-gained variants S1630* and Y1655*, associated with ovarian cancer and head-neck squamous carcinoma, respectively. The PSS in 1708 maps close to the pathogenic mutation Y1762* associated with ovarian cancer. In the intervening regions between the BRC6-BRC7-BRC8 repeats, we detected a clustering of three PSSs (1913, 2035, and 2037) with two pathogenic variants predisposing to breast cancer (E1953* and S1955*) and the mutation K2013* related to ovarian cancer. Three PSSs (2530, 2572 and 2574) and three pathogenic variants (R2494* associated with bladder urothelial carcinoma, R2520* with ovarian cancer, and W2626 with rectum adenocarcinoma) were found co-localized in the helical domain.
PTPRC encodes a tyrosine phosphatase also known as CD45 that regulates T- and B-cell antigen receptor signaling and is associated with oncogenic transformation through changes in expression [37]. We detected 36 and 52 PSSs under M2a and M8, respectively (Table 1). Remarkably, the 52 M8 PSSs concentrate on the extracellular region involved in T cell receptor activation [38], whereas the cytoplasmic segment, which contains the phosphatase domains, seems to be under strong purifying selection (no PSSs; Figure 4B). Within the extracellular region (positions 26–577 of human PTPRC), the PSSs cluster in the cysteine-rich (CR) domain and across the three Fibronectin type 3 (FN3) domains (Figure 4C). We did not find information about germline variants in PTPRC associated with cancer.

4. Discussion

4.1. Cancer Genes Show Relatively Low dN/dS Values

Because human cancer genes are generally involved in essential cellular functions such as DNA repair, regulation of cell cycle, and apoptosis, strong purifying selection removing deleterious germline mutations is expected, resulting in a dN/dS << 1. On average, the estimated dN/dS ratio across cancer genes (0.12) was somehow lower than previous estimates obtained from mammalian genomes, which yielded values between 0.15 and 0.22 [6,39,40,41]. On the other hand, Thomas et al. obtained a lower dN/dS for cancer-related genes (dN/dS = 0.079) than for other disease-related genes (dN/dS = 0.101) or for non-disease-related genes (dN/dS = 0.100), when comparing humans and rodents [42]. Our results also concord with Blekhman et al., who found significantly lower dN/dS values (dN/dS = 0.061) in human genes associated with cancer compared to genes involved in Mendelian (dN/dS = 0.133) and complex diseases (dN/dS = 0.203) [43].

4.2. Positive Selection on Human Cancer Genes is Associated with Hereditary Cancer and Recessive Mutations

Human cancer genes bearing only germline mutations in COSMIC yielded higher dN/dS ratios and higher proportion of PSGs than human cancer genes with only somatic mutations. This result suggests that genes associated with hereditary cancer have less selective constraints than those genes related to sporadic cancer. Indeed, it is expected that slightly deleterious (and advantageous) mutations under weak purifying selection reach higher frequencies in the populations than under strong purifying selection [44]. Therefore, it is expected that variants under relaxed purifying selection are more likely associated with hereditary cancer than with sporadic cancer, as mutations in the latter are expected to be highly deleterious and are therefore removed by purifying selection. In addition, since cancer is a complex and often a late-onset disease, where each allele contributes to a small fraction of cancer risk, it is plausible that variants directly associated with cancer susceptibility are under relaxed selective constraints [43,45,46]. COSMIC does not provide information about the exact number of germline and somatic mutations in each gene category, so we could not evaluate the impact of mutation burden on selective constraints.
We detected that genes associated with recessive mutations show significantly higher dN/dS values and were more often positively selected than genes associated with dominant mutations. It is expected that genes with dominant effects are subjected to stronger purifying selection than those with recessive effects, as deleterious mutations in the former will have immediate disadvantages in heterozygosity. It has been reported that genes with recessive disease mutations have higher dN/dS than genes with dominant disease mutations [43,47]. This result is not independent from the one just discussed, as almost all cancer genes with germline mutations in our dataset were associated with recessive inheritance. Because genes associated with germline mutations and dominant inheritance are unrepresented in our dataset, we were unable to test the potential interaction between mutation type and inheritance mode.

4.3. Lack of Variation in Selection across Tissues or Cancer Gene Role

Our analysis did not reveal significant differences in the dN/dS ratios or the proportion of PSGs among tissue types (epithelial, lymphoid, mesenchymal, and others) or cancer role (fusion genes, oncogenes, and TSGs). This might suggest that selective pressures on human cancer genes along the mammal lineage are not directly related to cancer. Since blood and bone marrow cancers are promoted by alterations in the immune system, and genes involved in immunity often undergo fast adaptation [6,7,8], we would have expected a stronger signal of positive selection on genes associated with lymphoid cancers. A previous study on the evolution of cancer genes in hominoids identified a higher dN/dS in TSGs with respect to oncogenes [11]. Although the difference was not statistically significant, we observed a higher average dN/dS for TSGs than for oncogenes. These patterns might be attributed to TSGs being more enriched in genes with recessive mutations (109 out of 206, 52.9%, TSGs carry recessive mutations) than oncogenes (15 out of 208, 7.2%). Nonetheless, the proportion of PSGs was very similar between fusion genes and TSGs (Figure 3C), while the former included very few genes with recessive mutations (8 out of 225, 3.5%). Therefore, although the inheritance factor might influence the strength of purifying selection, positive selection is likely driven by other properties.

4.4. Signalling Pathways and Biological Functions of Cancer Genes under Positive Selection

Within the putative list of PSGs, we found an enrichment of genes involved in DNA repair, with a high presence of genes involved in the Fanconi Anemia (FA)/BRCA pathway. Evidence of adaptive evolution in some components of the FA/BRCA pathway (BRCA2, CHEK2, FANCC, FANCB, FANCD2, and FANCE) has been previously identified in mammals [48]. In addition to these genes, we also identified signatures of positive selection in FANCA and FANCG. The systematic positive selection observed on the FA/BRCA complex might suggest a mechanism of coevolution to maintain the interactions among partners of this network [48,49]. Positive selection on the FA/BRCA complex could be driven by different selective pressures. On one hand, positive selection on this DNA repair pathway could favor a molecular mechanism of tumor resistance to counteract the increased cancer risk associated with longevity [50]. This hypothesis is supported by the signature of positive selection of some FA/BRCA components (BRCA2, FANCA, FANCE, and FANCL) identified in long-lived and cancer-resistant species [51,52]. On the other hand, germline variants in the FA/BRCA repair pathway have been associated with hereditary breast–ovarian cancer and Fanconi Anemia in humans [30,32]. Therefore, it is possible that, at least a portion of selected alleles in the FA/BRCA pathway could be pleiotropic and have deleterious late-onset effects like higher cancer risk. Future biochemical characterization of mutants and genetic association studies will help to better understand the consequences of positive selection on the components of the FA/BRCA pathway.
As mentioned in the introduction, antagonistic pleiotropy between positive selection associated with organismal fitness and increased cancer risk could be a general mechanism behind the long-term molecular adaptation of human cancer genes [6,7,8,10]. Immune response, placentation, and spermatogenesis, expected to be shaped by pathogen–host coevolution, maternal–fetal interactions, and sexual selection, respectively, are biological processes also often associated with positive selection on cancer genes [9,10]. Interestingly, we observed a high proportion of immunity-related genes under positive selection, with enrichment in the process of T cell proliferation, which would support the hypothesis of pathogen–host interactions driving adaptive changes of cancer-related genes. At the same time, the identification of several positively selected genes expressed in testis (such as BRIP1, BUB1B, KTN1, and RANBP2) is also concordant with the hypothesis that the genetic pathways of spermatogenesis, which evolve in response to sexual selection and intrasexual conflict, often coincide with those used by cancer cells to increase their survival and replication [6,9,10].

4.5. Functional Relevance of Residues under Positive Selection in Cancer Genes

We identified 25 PSSs in BRCA2, where some positions mapped close to natural variants associated with cancer. A similar signature of positive selection in BRCA2 was previously identified by O’Connell [48]. Although no PSS matched with pathogenic variants, which is not surprising since cancer is not a driver of species adaptation, they might be involved in cancer susceptibility or tumor resistance, especially when located in regions of functional importance. Future biochemical studies will help to assess the evolutionary significance of variants in PSSs. It is worth to mention that, although we focused on the potential effect of mutations in PSSs, it is expected that most of the mutations associated with cancer risk fall on sites under strong purifying selection during species evolution. Therefore, in the future, it would be interesting to consider the most conserved domains of human cancer genes as potential candidate disease regions.

Data Availability

The data and code required to reproduce the results of this study are available in https://github.com/avicens/cancer_genes_selection.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/9/12/582/s1, Figure S1: Mammal phylogenetic tree of mammalian species assembled for this study, Figure S2: Global dN and dS values according to mutation type (A,B) and inheritance (C,D), Figure S3: (A) Correlation between protein sequence length and global dN/dS estimates; (B) Comparison of protein length between positively selected and not selected genes; (C) Comparison of protein length among mutation-type categories. Statistical p-values are shown for multiple and pairwise comparisons. Mean of sequence length is indicated below each category in (B,C), Figure S4: Proportion of positively selected sites across COSMIC categories: (A) mutation type; (B) genetic dominance; (C) tissue type; and (D) cancer role. The number of genes in each category is indicated within each square. Significance levels for chi-square tests are indicated below each plot: non-significant (ns) and significant p-value < 0.05 (*). Abbreviated categories: Som + Ger: Genes bearing both somatic and germline mutations; TSG: tumor suppressor genes. Table S1; List of Cancer Gene Census of COSMIC genes used in the study. Table S2: Genes discarded from the analysis in the different pre-processing steps, Table S3: Genes discarded from the analysis in the different pre-processing steps, Table S4: Protein Ensembl accessions of all orthologues retrieved in this study, Table S5: Evolutionary parameters estimated for each gene.

Author Contributions

Conceptualization, A.V. and D.P.; Methodology, A.V. and D.P.; Formal Analysis, A.V..; Investigation, A.V. and D.P.; Writing-Original Draft Preparation, A.V.; Writing-Review & Editing, D.P.; Supervision, D.P.

Funding

This research was funded by the Spanish Government (Juan de la Cierva postdoctoral fellowship IJCI-2016-29550 to A.V and research grant BFU2015-63774-P to D.P) and the European Research Council (ERC-617457-PHYLOCANCER to D.P.).

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Bodmer, W.; Tomilson, I. Rare genetic variants and the risk of cancer. Curr. Opin. Genet. Dev. 2010, 20, 262–267. [Google Scholar] [CrossRef] [PubMed]
  2. Bailey, M.H.; Tokheim, C.; Porta-Pardo, E.; Sengupta, S.; Bertrand, D.; Weerasinghe, A.; Colaprico, A.; Wendl, M.C.; Kim, J.; Reardon, B. Comprehensive characterization of cancer driver genes and mutations. Cell 2018, 173, 371–385.e18. [Google Scholar] [CrossRef]
  3. Martincorena, I.; Raine, K.M.; Gerstung, M.; Dawson, K.J.; Haase, K.; Van Loo, P.; Davies, H.; Stratton, M.R.; Campbell, P.J. Universal patterns of selection in cancer and somatic tissues. bioRxiv 2017. [Google Scholar] [CrossRef] [PubMed]
  4. Sondka, Z.; Bamford, S.; Cole, C.G.; Ward, S.A.; Dunham, I.; Forbes, S.A. The COSMIC Cancer Gene Census: describing genetic dysfunction across all human cancers. Nat. Rev. Cancer 2018. [Google Scholar] [CrossRef] [PubMed]
  5. Clark, A.G.; Glanowski, S.; Nielsen, R.; Thomas, P.D.; Kejariwal, A.; Todd, M.A.; Tanenbaum, D.M.; Civello, D.; Lu, F.; Murphy, B.; et al. Inferring nonneutral evolution from human-chimp-mouse orthologous gene trios. Science 2003, 302, 1960–1963. [Google Scholar] [CrossRef] [PubMed]
  6. Nielsen, R.; Bustamante, C.; Clark, A.G.; Glanowski, S.; Sackton, T.B.; Hubisz, M.J.; Fledel-Alon, A.; Tanenbaum, D.M.; Civello, D.; White, T.J.; et al. A scan for positively selected genes in the genomes of humans and chimpanzees. PLoS Biol. 2005, 3, e170. [Google Scholar] [CrossRef]
  7. Kosiol, C.; Vinař, T.; da Fonseca, R.R.; Hubisz, M.J.; Bustamante, C.D.; Nielsen, R.; Siepel, A. Patterns of positive selection in six mammalian genomes. PLoS Genet. 2008, 4, e1000144. [Google Scholar] [CrossRef]
  8. Da Fonseca, R.R.; Kosiol, C.; Vinař, T.; Siepel, A.; Nielsen, R. Positive selection on apoptosis related genes. FEBS Lett. 2010, 584, 469–476. [Google Scholar] [CrossRef]
  9. Kleene, K.C. Sexual selection, genetic conflict, selfish genes, and the atypical patterns of gene expression in spermatogenic cells. Dev. Biol. 2005, 277, 16–26. [Google Scholar] [CrossRef]
  10. Crespi, B.J.; Summers, K. Positive selection in the evolution of cancer. Biol. Rev. Camb. Philos. Soc. 2006, 81, 407–424. [Google Scholar] [CrossRef]
  11. Kang, L.; Michalak, P. The evolution of cancer-related genes in hominoids. J. Mol. Evol. 2015, 80, 37–41. [Google Scholar] [CrossRef] [PubMed]
  12. Durinck, S.; Moreau, Y.; Kasprzyk, A.; Davis, S.; De Moor, B.; Brazma, A.; Huber, W. BioMart and Bioconductor: a powerful link between biological databases and microarray data analysis. Bioinformatics 2005, 21, 3439–3440. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Ranwez, V.; Harispe, S.; Delsuc, F.; Douzery, E.J.P. MACSE: Multiple Alignment of Coding SEquences accounting for frameshifts and stop codons. PLoS ONE 2011, 6, e22594. [Google Scholar] [CrossRef] [PubMed]
  14. Capella-Gutiérrez, S.; Silla-Martínez, J.M.; Gabaldón, T. trimAl: A tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 2009, 25, 1972–1973. [Google Scholar] [CrossRef] [PubMed]
  15. Rodríguez, F.; Oliver, J.L.; Marín, A.; Medina, J.R. The general stochastic model of nucleotide substitution. J. Theor. Biol. 1990, 142, 485–501. [Google Scholar] [CrossRef]
  16. Yang, Z. Maximum-likelihood estimation of phylogeny from DNA sequences when substitution rates differ over sites. Mol. Biol. Evol. 1993, 10, 1396–1401. [Google Scholar] [CrossRef]
  17. Felsenstein, J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 1985, 39, 783–791. [Google Scholar] [CrossRef] [PubMed]
  18. Huerta-Cepas, J.; Serra, F.; Bork, P. ETE 3: reconstruction, analysis, and visualization of phylogenomic data. Mol. Biol. Evol. 2016, 33, 1635–1638. [Google Scholar] [CrossRef]
  19. Yang, Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24, 1586–1591. [Google Scholar] [CrossRef]
  20. Gharib, W.H.; Robinson-Rechavi, M. The branch-site test of positive selection is surprisingly robust but lacks power under synonymous substitution saturation and variation in GC. Mol. Biol. Evol. 2013, 30, 1675–1686. [Google Scholar] [CrossRef]
  21. Yang, Z.; Nielsen, R.; Goldman, N.; Pedersen, A. Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics 2000, 155, 431–449. [Google Scholar] [PubMed]
  22. Wong, W.S.W.; Yang, Z.; Goldman, N.; Nielsen, R. Accuracy and power of statistical methods for detecting adaptive evolution in protein coding sequences and for identifying positively selected sites. Genetics 2004, 168, 1041–1051. [Google Scholar] [CrossRef] [PubMed]
  23. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  24. Murrell, B.; Weaver, S.; Smith, M.D.; Wertheim, J.O.; Murrell, S.; Aylward, A.; Eren, K.; Pollner, T.; Martin, D.P.; Smith, D.M.; et al. Gene-wide identification of episodic selection. Mol. Biol. Evol. 2015, 32, 1365–1371. [Google Scholar] [CrossRef]
  25. Pond, S.L.K.; Frost, S.D.W.; Muse, S. V HyPhy: hypothesis testing using phylogenies. Bioinformatics 2005, 21, 676–679. [Google Scholar] [CrossRef] [PubMed]
  26. Yang, Z.; Wong, S.W.; Nielsen, R. Bayes empirical bayes inference of amino acid sites under positive selection. Mol. Biol. Evol. 2005, 22, 1107–1118. [Google Scholar] [CrossRef] [PubMed]
  27. Eden, E.; Navon, R.; Steinfeld, I.; Lipson, D.; Yakhini, Z. GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics 2009, 10, 48. [Google Scholar] [CrossRef]
  28. Huang, K.; Mashl, R.J.; Wu, Y.; Ritter, D.I.; Wang, J.; Oh, C.; Paczkowska, M.; Reynolds, S.; Wyczalkowski, M.A.; Oak, N.; et al. Pathogenic germline variants in 10,389 adult cancers. Cell 2018, 173, 355–370. [Google Scholar] [CrossRef]
  29. Team, R. core R: A language and Environment for statistical computing. Available online: http://www.r-project.org (accessed on 5 May 2018).
  30. Bogliolo, M.; Surrallés, J. The Fanconi Anemia/BRCA pathway: FANCD2 at the crossroad between repair and checkpoint responses to DNA damage. In Madame Curie Bioscience Database; Landes Bioscience: Austin, TX, USA, 2013. [Google Scholar]
  31. Wang, H.-F.; Takenaka, K.; Nakanishi, A.; Miki, Y. BRCA2 and nucleophosmin coregulate centrosome amplification and form a complex with the ρ effector kinase ROCK2. Cancer Res. 2011, 71, 68–77. [Google Scholar] [CrossRef]
  32. Roy, R.; Chun, J.; Powell, S.N. BRCA1 and BRCA2: Different roles in a common pathway of genome protection. Nat. Rev. Cancer 2012, 12, 68–78. [Google Scholar] [CrossRef]
  33. Buisson, R.; Niraj, J.; Pauty, J.; Maity, R.; Zhao, W.; Coulombe, Y.; Sung, P.; Masson, J.-Y. Breast cancer proteins PALB2 and BRCA2 stimulate polymerase η in recombination-associated DNA synthesis at blocked replication forks. Cell Rep. 2014, 6, 553–564. [Google Scholar] [CrossRef] [PubMed]
  34. Dayhoff, M.; Schwartz, R.; Orcutt, B. Atlas of Protein Sequence and Structure; Dayhoff, M., Ed.; National b.: Washington, DC, USA, 1978. [Google Scholar]
  35. Marston, N.J.; Richards, W.J.; Hughes, D.; Bertwistle, D.; Marshall, C.J.; Ashworth, A. Interaction between the product of the breast cancer susceptibility gene BRCA2 and DSS1, a protein functionally conserved from yeast to mammals. Mol. Cell. Biol. 1999, 19, 4633–4642. [Google Scholar] [CrossRef] [PubMed]
  36. Hussain, S.; Wilson, J.B.; Medhurst, A.L.; Hejna, J.; Witt, E.; Ananth, S.; Davies, A.; Masson, J.-Y.; Moses, R.; West, S.C.; et al. Direct interaction of FANCD2 with BRCA2 in DNA damage response pathways. Hum. Mol. Genet. 2004, 13, 1241–1248. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Yuan, L.; Zeng, G.; Chen, L.; Wang, G.; Wang, X.; Cao, X.; Lu, M.; Liu, X.; Qian, G.; Xiao, Y.; et al. Identification of key genes and pathways in human clear cell renal cell carcinoma (ccRCC) by co-expression analysis. Int. J. Biol. Sci. 2018, 14, 266–279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Chang, V.T.; Fernandes, R.A.; Ganzinger, K.A.; Lee, S.F.; Siebold, C.; McColl, J.; Jönsson, P.; Palayret, M.; Harlos, K.; Coles, C.H.; et al. Initiation of T cell signaling by CD45 segregation at “close contacts”. Nat. Immunol. 2016, 17, 574–582. [Google Scholar] [CrossRef] [PubMed]
  39. Romiguier, J.; Ranwez, V.; Douzery, E.J.P.; Galtier, N. Genomic evidence for large, long-lived ancestors to placental mammals. Mol. Biol. Evol. 2013, 30, 5–13. [Google Scholar] [CrossRef] [PubMed]
  40. Figuet, E.; Nabholz, B.; Bonneau, M.; Mas Carrio, E.; Nadachowska-Brzyska, K.; Ellegren, H.; Galtier, N. Life history traits, protein evolution, and the nearly neutral theory in amniotes. Mol. Biol. Evol. 2016, 33, 1517–1527. [Google Scholar] [CrossRef] [PubMed]
  41. Bustamante, C.D.; Fledel-Alon, A.; Williamson, S.; Nielsen, R.; Todd Hubisz, M.; Glanowski, S.; Tanenbaum, D.M.; White, T.J.; Sninsky, J.J.; Hernandez, R.D.; et al. Natural selection on protein-coding genes in the human genome. Nature 2005, 437, 1153–1157. [Google Scholar] [CrossRef]
  42. Thomas, M.A.; Weston, B.; Joseph, M.; Wu, W.; Nekrutenko, A.; Tonellato, P.J. Evolutionary dynamics of oncogenes and tumor suppressor genes: Higher intensities of purifying selection than other genes. Mol. Biol. Evol. 2003, 20, 964–968. [Google Scholar] [CrossRef]
  43. Blekhman, R.; Man, O.; Herrmann, L.; Boyko, A.R.; Indap, A.; Kosiol, C.; Bustamante, C.D.; Teshima, K.M.; Przeworski, M. Natural selection on genes that underlie human disease susceptibility. Curr. Biol. 2008, 18, 883–889. [Google Scholar] [CrossRef]
  44. Eyre-Walker, A.; Keightley, P.D.; Smith, N.G.C.; Gaffney, D. Quantifying the slightly deleterious mutation model of molecular evolution. Mol. Biol. Evol. 2002, 19, 2142–2149. [Google Scholar] [CrossRef] [PubMed]
  45. Wright, A.; Charlesworth, B.; Rudan, I.; Carothers, A.; Campbell, H. A polygenic basis for late-onset disease. Trends Genet. 2003, 19, 97–106. [Google Scholar] [CrossRef]
  46. Spataro, N.; Rodríguez, J.A.; Navarro, A.; Bosch, E. Properties of human disease genes and the role of genes linked to Mendelian disorders in complex disease aetiology. Hum. Mol. Genet. 2017, 26, ddw405. [Google Scholar] [CrossRef] [PubMed]
  47. Furney, S.J.; Albà, M.M.; López-Bigas, N. Differences in the evolutionary history of disease genes affected by dominant or recessive mutations. BMC Genomics 2006, 7, 165. [Google Scholar] [CrossRef]
  48. O’Connell, M.J. Selection and the cell cycle: Positive darwinian selection in a well-known DNA damage response pathway. J. Mol. Evol. 2010, 71, 444–457. [Google Scholar] [CrossRef] [PubMed]
  49. Qian, W.; Zhou, H.; Tang, K. Recent coselection in human populations revealed by protein–protein interaction network. Genome Biol. Evol. 2015, 7, 136–153. [Google Scholar] [CrossRef] [PubMed]
  50. Tollis, M.; Schiffman, J.D.; Boddy, A.M. Evolution of cancer suppression as revealed by mammalian comparative genomics. Curr. Opin. Genet. Dev. 2017, 42, 40–47. [Google Scholar] [CrossRef]
  51. Morgan, C.C.; Mc Cartney, A.M.; Donoghue, M.T.A.; Loughran, N.B.; Spillane, C.; Teeling, E.C.; O’Connell, M.J. Molecular adaptation of telomere associated genes in mammals. BMC Evol. Biol. 2013, 13. [Google Scholar] [CrossRef]
  52. Zhang, G.; Cowled, C.; Shi, Z.; Huang, Z.; Bishop-Lilly, K.A.; Fang, X.; Wynne, J.W.; Xiong, Z.; Baker, M.L.; Zhao, W.; et al. Comparative analysis of bat genomes provides insight into the evolution of flight and immunity. Science 2013, 339, 456–460. [Google Scholar] [CrossRef]
Genes 09 00582 g001 550
Figure 1. Distribution of global number of non-synonymous substitutions per non-synonymous site to the number of synonymous substitutions per synonymous site (dN/dS) estimates. Dashed, red, and blue vertical lines indicate mean dN/dS for all (0.122), positively selected (0.290), and not positively selected genes (0.105), respectively.
Figure 1. Distribution of global number of non-synonymous substitutions per non-synonymous site to the number of synonymous substitutions per synonymous site (dN/dS) estimates. Dashed, red, and blue vertical lines indicate mean dN/dS for all (0.122), positively selected (0.290), and not positively selected genes (0.105), respectively.
Genes 09 00582 g001
Genes 09 00582 g002 550
Figure 2. Global dN/dS across COSMIC categories. (A) Mutation type; (B) inheritance; (C) tissue type, and (D) cancer role; p-values are shown for each comparison: (ns): p-value > 0.01; (**): p-value < 0.01; (***): p-value < 0.001.
Figure 2. Global dN/dS across COSMIC categories. (A) Mutation type; (B) inheritance; (C) tissue type, and (D) cancer role; p-values are shown for each comparison: (ns): p-value > 0.01; (**): p-value < 0.01; (***): p-value < 0.001.
Genes 09 00582 g002
Genes 09 00582 g003 550
Figure 3. Number of positively selected genes across COSMIC categories. (A) Mutation type; (B) inheritance; (C) tissue type, and (D) cancer role; p-values for chi-squared tests are shown underneath (significance p < 0.05).
Figure 3. Number of positively selected genes across COSMIC categories. (A) Mutation type; (B) inheritance; (C) tissue type, and (D) cancer role; p-values for chi-squared tests are shown underneath (significance p < 0.05).
Genes 09 00582 g003
Genes 09 00582 g004 550
Figure 4. Positive selected sites in BRCA2 and PTPRC. (A) Mapping of positively selected sites (PSSs; red) and cancer predisposition variants (missensense and stop-gained; blue) in human BRCA2. Proteins interacting with different domains of BRCA2 are indicated underneath. Domain abbreviations: H: Helical domain; OB: Oligonucleotide Binding; NLS: Nuclear Localization Sequence; (B) Site-specific posterior probabilities of dN/dS > 1 in PTPRC, under the selection model M2a, are represented on the y-axis. The domain scheme is represented in the x-axis. The dashed line indicates the significance level for positive selected sites (Bayes Empirical Bayes (BEB) prob = 0.95). Domain abbreviations: CR: cysteine-rich; FN3 type 3: fibronectin 3; TM: transmembrane; (C) Crystal structure of the extracellular region of PTPRC (PDB ID: 5fmv, Chang et al., 2016), showing CR and FN3 domains with different colors. PSSs are shown as red spheres.
Figure 4. Positive selected sites in BRCA2 and PTPRC. (A) Mapping of positively selected sites (PSSs; red) and cancer predisposition variants (missensense and stop-gained; blue) in human BRCA2. Proteins interacting with different domains of BRCA2 are indicated underneath. Domain abbreviations: H: Helical domain; OB: Oligonucleotide Binding; NLS: Nuclear Localization Sequence; (B) Site-specific posterior probabilities of dN/dS > 1 in PTPRC, under the selection model M2a, are represented on the y-axis. The domain scheme is represented in the x-axis. The dashed line indicates the significance level for positive selected sites (Bayes Empirical Bayes (BEB) prob = 0.95). Domain abbreviations: CR: cysteine-rich; FN3 type 3: fibronectin 3; TM: transmembrane; (C) Crystal structure of the extracellular region of PTPRC (PDB ID: 5fmv, Chang et al., 2016), showing CR and FN3 domains with different colors. PSSs are shown as red spheres.
Genes 09 00582 g004
Table
Table 1. List of cancer genes showing evidence of positive selection.
Table 1. List of cancer genes showing evidence of positive selection.
GeneFunctiondN/dS
IL2T cell proliferation and regulation of the immune response.0.748
FASApoptosis0.601
FCRL4B cell receptor signaling0.601
NUTM2AUnknown0.574
PALB2Tumor necrosis factor, apoptosis0.507
PDCD1LG2T cell proliferation; immune response0.506
FANCGFanconi Anemia (FA) group; DNA repair0.500
BRCA2Double-strand break repair and/or homologous recombination0.468
CD274T cell effector regulation; attenuation of anti-tumor immunity0.465
FANCCF.A. group; DNA repair0.445
CASP8Protease inhibitor; apoptosis0.363
PTPRCProtein phosphatase; receptor; immune response0.361
FANCD2F.A. group; DNA repair0.348
BARD1Control of the cell cycle in response to DNA damage0.333
ERCC5DNA repair0.321
NCOA4Androgen receptor signaling0.312
NINCentrosome localization0.307
BRIP1Double-strand break repair and/or homologous recombination0.286
COL1A1Collagen component0.276
CD79BB cell differentiation and activation0.275
BLMBasic helix-loop transcription factor0.264
CD79AB cell differentiation and activation0.253
PMS2DNA binding protein0.225
KTN1Kinesin-driven vesicle motility; cadherin binding0.211
PRF1Apoptosis; immune response0.197
SETChaperone; phosphatase inhibitor0.180
ARHGEF12Regulation of RhoA GTPase0.178
CHEK2Checkpoint-mediated cell cycle arrest, activation of DNA repair and apoptosis0.175
PTPRBProtein phosphatase; receptor; angiogenesis0.156
SS18Chromatin-binding protein; transcription regulation0.153
FLT3Regulation of apoptotic process0.144
COL2A1Collagen component 0.132
MLLT6Nucleic acid binding; zinc finger transcription factor0.130
KDM6ATranscription factor; chromatin remodeling0.120
POU2AF1Transcriptional coactivator; immune response0.106
MED12Nucleic acid binding; transcription cofactor0.097
RBM15RNA binding protein0.088
RABEP1Membrane fusion; apoptosis0.086
BRAFTransduction of mitogenic signals; apoptosis0.079
PICALMVesicle coat protein0.068

Share and Cite

MDPI and ACS Style

Vicens, A.; Posada, D. Selective Pressures on Human Cancer Genes along the Evolution of Mammals. Genes 2018, 9, 582. https://doi.org/10.3390/genes9120582

AMA Style

Vicens A, Posada D. Selective Pressures on Human Cancer Genes along the Evolution of Mammals. Genes. 2018; 9(12):582. https://doi.org/10.3390/genes9120582

Chicago/Turabian Style

Vicens, Alberto, and David Posada. 2018. "Selective Pressures on Human Cancer Genes along the Evolution of Mammals" Genes 9, no. 12: 582. https://doi.org/10.3390/genes9120582

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop