Mathematical Modeling: Tuning the MYC response

  1. Ying Zheng
  2. David Levens  Is a corresponding author
  1. National Cancer Institute, United States

Cancer develops when cells in the body gain mutations that enable them to divide rapidly and form a tumor. Some of these mutations can increase the expression of particular genes in cells. For example, a gene that encodes a transcription factor protein called MYC is upregulated in most types of tumor, where it supports several aspects of tumor development.

Transcription factors can bind to promoter regions at the start of genes to regulate the first step in gene expression, a process called transcription. MYC can bind to virtually all promoters in the genome and trigger widespread increases in gene activity. However, it is not clear whether MYC binding directly increases the transcription of all the genes (Levens, 2013; Lin et al., 2012; Nie et al., 2012; Wolf et al., 2015), or whether MYC only regulates specific sets of genes that then increase the expression of other genes (Kress et al., 2015).

Previous studies have shown that producing more MYC in cell lines does not increase the expression of all genes equally (Lin et al., 2012; Sabò et al., 2014; Walz et al., 2014). Now, in eLife, Elmar Wolf at the University of Würzburg and colleagues – including Francesca Lorenzin as first author – show that changing the amount of MYC in cells does not change the amount of MYC bound at each gene in the same way (Lorenzin et al., 2016). This prompted the researchers to seek an explanation for the relationship between MYC levels and gene expression across all genes.

Lorenzin et al. – who are based at the University of Würzburg and the Max-Delbrück-Center for Molecular Medicine – manipulated the levels of MYC in a human cell line called U20S, which normally has less of this protein compared to many cancer cell lines. Upregulating MYC in U2OS cells to match the levels seen in cancer resulted in MYC molecules being bound to virtually all promoters.

Further experiments indicated that MYC binding to any promoter can alter the activity of the corresponding gene, but that MYC is better at binding to some promoters than others. At normal levels of MYC, promoters where MYC binds weakly (i.e. those with a low binding affinity) will have few MYC molecules bound to them. However, under the same conditions, high-affinity promoters may already be saturated with MYC molecules. If cells make more MYC, the number of MYC molecules binding to the lower-affinity promoters would increase, but the high-affinity promoters would be largely unaffected. Therefore, the increases in MYC production observed in cancer cells can alter the transcription of some genes more than others. Because every cancer starts with a different pattern of gene expression and a different level of MYC, this may result in each tumor having a distinct transcription profile.

What determines how strongly MYC binds to a promoter? Lorenzin et al. modeled MYC binding across the genome, taking into account that MYC is known to bind strongly to a DNA motif called the E-box (Guo et al., 2014), and less tightly to other DNA sequences (Figure 1). But the model based on E-boxes alone did not explain how MYC is distributed across the genome. Lorenzin et al. then considered the possibility that interactions between MYC and other proteins may help to draw MYC to the promoters of active genes. For example, a protein called WDR5 was speculated to help MYC target genes (Thomas et al., 2015). Indeed, including interactions between MYC and WDR5 markedly improved the model for many (but not all) genes. Moreover, impairing this interaction had a greater impact on the expression of genes with high affinity for MYC than those with low affinity.

The ability of MYC to bind to individual promoters is influenced by interactions with DNA and other proteins.

The MYC protein forms a dimer with another protein called MAX, which allows it to bind to virtually any promoter in the genome and increase gene expression by activating transcription (green arrow). However, the ability of MYC to bind to individual promoters (binding affinity) varies. MYC prefers to bind to DNA motifs called E-boxes, but will also bind to other DNA sequences with lower affinity. A protein called WDR5, which is associated with the DNA of active genes, can help to recruit MYC to particular promoters (beige arrow). MYC also interacts with other DNA-associated proteins – including TRAPP, Paf1, BRD4 and PTEF-b – during transcription, which may also help MYC to bind to certain promoters with higher affinity.

Besides WDR5, MYC interacts with many other proteins and protein complexes during transcription (Figure 1). Do these interactions also help to recruit and retain MYC at promoters, and if so, how might they be incorporated into the model? Perhaps MYC interacts with different partners in different stages of transcription to ensure that the process happens in an efficient and orderly manner. It should be noted that some of these partners might trigger the degradation or modification of MYC leading to its ejection from certain promoters, which would change its apparent binding affinity, and in turn change the outputs of these promoters.

Finally, normal cells contain less MYC than the cells used by Lorenzin et al. and so the majority of promoters in a normal cell are not saturated with MYC molecules. In cancer cells, the genes with the highest affinity for MYC, such as the ribosomal protein genes, are saturated with MYC and are also the most sensitive to its inhibition or reduction. Understanding exactly how MYC binding influences gene expression promises to aid efforts to develop new therapeutic strategies that inhibit its cancer-promoting actions, while preserving its normal responsibilities in cells.

References

Article and author information

Author details

  1. Ying Zheng

    Laboratory of Pathology, Center for Cancer Research, National Cancer Institute, Bethesda, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0674-2438
  2. David Levens

    Laboratory of Pathology, Center for Cancer Research, National Cancer Institute, Bethesda, United States
    For correspondence
    levensd@mail.nih.gov
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-7616-922X

Publication history

  1. Version of Record published:

Copyright

© 2016, Zheng et al.

This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.

Metrics

  • 1,527
    views
  • 330
    downloads
  • 2
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Ying Zheng
  2. David Levens
(2016)
Mathematical Modeling: Tuning the MYC response
eLife 5:e18871.
https://doi.org/10.7554/eLife.18871

Further reading

    1. Cancer Biology
    2. Computational and Systems Biology
    Rosalyn W Sayaman, Masaru Miyano ... Mark LaBarge
    Research Article

    Effects from aging in single cells are heterogenous, whereas at the organ- and tissue-levels aging phenotypes tend to appear as stereotypical changes. The mammary epithelium is a bilayer of two major phenotypically and functionally distinct cell lineages: luminal epithelial and myoepithelial cells. Mammary luminal epithelia exhibit substantial stereotypical changes with age that merit attention because these cells are the putative cells-of-origin for breast cancers. We hypothesize that effects from aging that impinge upon maintenance of lineage fidelity increase susceptibility to cancer initiation. We generated and analyzed transcriptomes from primary luminal epithelial and myoepithelial cells from younger <30 (y)ears old and older >55y women. In addition to age-dependent directional changes in gene expression, we observed increased transcriptional variance with age that contributed to genome-wide loss of lineage fidelity. Age-dependent variant responses were common to both lineages, whereas directional changes were almost exclusively detected in luminal epithelia and involved altered regulation of chromatin and genome organizers such as SATB1. Epithelial expression of gap junction protein GJB6 increased with age, and modulation of GJB6 expression in heterochronous co-cultures revealed that it provided a communication conduit from myoepithelial cells that drove directional change in luminal cells. Age-dependent luminal transcriptomes comprised a prominent signal that could be detected in bulk tissue during aging and transition into cancers. A machine learning classifier based on luminal-specific aging distinguished normal from cancer tissue and was highly predictive of breast cancer subtype. We speculate that luminal epithelia are the ultimate site of integration of the variant responses to aging in their surrounding tissue, and that their emergent phenotype both endows cells with the ability to become cancer-cells-of-origin and represents a biosensor that presages cancer susceptibility.

    1. Cancer Biology
    Jae Hun Shin, Jooyoung Park ... Alfred LM Bothwell
    Research Article

    Metastasis is the leading cause of cancer-related mortality. Paneth cells provide stem cell niche factors in homeostatic conditions, but the underlying mechanisms of cancer stem cell niche development are unclear. Here, we report that Dickkopf-2 (DKK2) is essential for the generation of cancer cells with Paneth cell properties during colon cancer metastasis. Splenic injection of Dkk2 knockout (KO) cancer organoids into C57BL/6 mice resulted in a significant reduction of liver metastases. Transcriptome analysis showed reduction of Paneth cell markers such as lysozymes in KO organoids. Single-cell RNA sequencing analyses of murine metastasized colon cancer cells and patient samples identified the presence of lysozyme positive cells with Paneth cell properties including enhanced glycolysis. Further analyses of transcriptome and chromatin accessibility suggested hepatocyte nuclear factor 4 alpha (HNF4A) as a downstream target of DKK2. Chromatin immunoprecipitation followed by sequencing analysis revealed that HNF4A binds to the promoter region of Sox9, a well-known transcription factor for Paneth cell differentiation. In the liver metastatic foci, DKK2 knockout rescued HNF4A protein levels followed by reduction of lysozyme positive cancer cells. Taken together, DKK2-mediated reduction of HNF4A protein promotes the generation of lysozyme positive cancer cells with Paneth cell properties in the metastasized colon cancers.