Phenotypic screen of sixty-eight colorectal cancer cell lines identifies CEACAM6 and CEACAM5 as markers of acid resistance.

Elevated cancer metabolism releases lactic acid and CO2 into the under-perfused tumor microenvironment, resulting in extracellular acidosis. The surviving cancer cells must adapt to this selection pressure; thus, targeting tumor acidosis is a rational therapeutic strategy to manage tumor growth. However, none of the major approved treatments are based explicitly on disrupting acid handling, signaling, or adaptations, possibly because the distinction between acid-sensitive and acid-resistant phenotypes is not clear. Here, we report pH-related phenotypes of sixty-eight colorectal cancer (CRC) cell lines by measuring i) extracellular acidification as a readout of acid production by fermentative metabolism and ii) growth of cell biomass over a range of extracellular pH (pHe) levels as a measure of the acid sensitivity of proliferation. Based on these measurements, CRC cell lines were grouped along two dimensions as "acid-sensitive"/"acid-resistant" versus "low metabolic acid production"/"high metabolic acid production." Strikingly, acid resistance was associated with the expression of CEACAM6 and CEACAM5 genes coding for two related cell-adhesion molecules, and among pH-regulating genes, of CA12. CEACAM5/6 protein levels were strongly induced by acidity, with a further induction under hypoxia in a subset of CRC lines. Lack of CEACAM6 (but not of CEACAM5) reduced cell growth and their ability to differentiate. Finally, CEACAM6 levels were strongly increased in human colorectal cancers from stage II and III patients, compared to matched samples from adjacent normal tissues. Thus, CEACAM6 is a marker of acid-resistant clones in colorectal cancer and a potential motif for targeting therapies to acidic regions within the tumors.

Goblet cell differentiation subgroups in colorectal cancer.

The poor prognosis of relatively undifferentiated cancers has long been recognized, suggesting that selection against differentiation and in favor of uncontrolled growth is one of the most powerful drivers of cancer progression. Goblet cells provide the mucous surface of the gut, and when present in colorectal cancers (CRC), the cancers are called mucinous. We have used the presence of MUC2, the main mucous product of goblet cells, and an associated gene product, TFF3, to classify a large panel of nearly 80 CRC-derived cell lines into five categories based on their levels of MUC2 and TFF3 expression. We have then shown that these five patterns of expression can be easily identified in the direct analysis of tumor specimens allowing a much finer characterization of CRCs with respect to the presence of goblet cell differentiation. In particular, about 30% of all CRCs fall into the category of expressing TFF3 but not MUC2, which has not previously been acknowledged. Using the cell line data, we suggest that there are up to 12 genes (MUC2, TFF3, ATOH1, SPDEF, CDX1, CDX2, GATA6, HES1, ETS2, OLFM4, TOX3, and LGR5) that may be involved in selection against goblet cell differentiation in CRC by changes in methylation rather than mutations. Of these, LGR5, which is particularly associated with lack of goblet cell features, may function in the control of differentiation rather than direct control of cell growth, as has so far mostly been assumed. These results emphasize the importance of methylation changes in driving cancer progression.

Two contrasting papers stimulate a commentary on the origins of tumour immunology, current cancer immunotherapies, and the future potential for cancer immunotherapy.

Two early papers expressing conflicting views on the occurrence of effective immune attack against cancers stimulate an analysis of the gradual development of an understanding of tumour biology. This understanding has led to the development of the strikingly effective check point blocking and CART anti-cancer immunotherapies, and the promise of more widely applicable therapies based on T cell attracting genetically engineered monoclonal antibodies.

Polygenic inheritance, GWAS, polygenic risk scores, and the search for functional variants.

The reconciliation between Mendelian inheritance of discrete traits and the genetically based correlation between relatives for quantitative traits was Fisher's infinitesimal model of a large number of genetic variants, each with very small effects, whose causal effects could not be individually identified. The development of genome-wide genetic association studies (GWAS) raised the hope that it would be possible to identify single polymorphic variants with identifiable functional effects on complex traits. It soon became clear that, with larger and larger GWAS on more and more complex traits, most of the significant associations had such small effects, that identifying their individual functional effects was essentially hopeless. Polygenic risk scores that provide an overall estimate of the genetic propensity to a trait at the individual level have been developed using GWAS data. These provide useful identification of groups of individuals with substantially increased risks, which can lead to recommendations of medical treatments or behavioral modifications to reduce risks. However, each such claim will require extensive investigation to justify its practical application. The challenge now is to use limited genetic association studies to find individually identifiable variants of significant functional effect that can help to understand the molecular basis of complex diseases and traits, and so lead to improved disease prevention and treatment. This can best be achieved by 1) the study of rare variants, often chosen by careful candidate assessment, and 2) the careful choice of phenotypes, often extremes of a quantitative variable, or traits with relatively high heritability.

GMMchi: gene expression clustering using Gaussian mixture modeling.

BACKGROUND: Cancer evolution consists of a stepwise acquisition of genetic and epigenetic changes, which alter the gene expression profiles of cells in a particular tissue and result in phenotypic alterations acted upon by natural selection. The recurrent appearance of specific genetic lesions across individual cancers and cancer types suggests the existence of certain "driver mutations," which likely make up the major contribution to tumors' selective advantages over surrounding normal tissue and as such are responsible for the most consequential aspects of the cancer cells' gene expression patterns and phenotypes. We hypothesize that such mutations are likely to cluster with specific dichotomous shifts in the expression of the genes they most closely control, and propose GMMchi, a Python package that leverages Gaussian Mixture Modeling to detect and characterize bimodal gene expression patterns across cancer samples, as a tool to analyze such correlations using 2 × 2 contingency table statistics. RESULTS: Using well-defined simulated data, we were able to confirm the robust performance of GMMchi, reaching 85% accuracy with a sample size of n = 90. We were also able to demonstrate a few examples of the application of GMMchi with respect to its capacity to characterize background florescent signals in microarray data, filter out uninformative background probe sets, as well as uncover novel genetic interrelationships and tumor characteristics. Our approach to analysing gene expression analysis in cancers provides an additional lens to supplement traditional continuous-valued statistical analysis by maximizing the information that can be gathered from bulk gene expression data. CONCLUSIONS: We confirm that GMMchi robustly and reliably extracts bimodal patterns from both colorectal cancer (CRC) cell line-derived microarray and tumor-derived RNA-Seq data and verify previously reported gene expression correlates of some well-characterized CRC phenotypes. AVAILABILITY: The Python package GMMchi and our cell line microarray data used in this paper is available for downloading on GitHub at https://github.com/jeffliu6068/GMMchi .

The fine-scale genetic structure of the British population.

Fine-scale genetic variation between human populations is interesting as a signature of historical demographic events and because of its potential for confounding disease studies. We use haplotype-based statistical methods to analyse genome-wide single nucleotide polymorphism (SNP) data from a carefully chosen geographically diverse sample of 2,039 individuals from the United Kingdom. This reveals a rich and detailed pattern of genetic differentiation with remarkable concordance between genetic clusters and geography. The regional genetic differentiation and differing patterns of shared ancestry with 6,209 individuals from across Europe carry clear signals of historical demographic events. We estimate the genetic contribution to southeastern England from Anglo-Saxon migrations to be under half, and identify the regions not carrying genetic material from these migrations. We suggest significant pre-Roman but post-Mesolithic movement into southeastern England from continental Europe, and show that in non-Saxon parts of the United Kingdom, there exist genetically differentiated subgroups rather than a general 'Celtic' population.

Single-molecule DNA-mapping and whole-genome sequencing of individual cells.

To elucidate cellular diversity and clonal evolution in tissues and tumors, one must resolve genomic heterogeneity in single cells. To this end, we have developed low-cost, mass-producible micro-/nanofluidic chips for DNA extraction from individual cells. These chips have modules that collect genomic DNA for sequencing or map genomic structure directly, on-chip, with denaturation-renaturation (D-R) optical mapping [Marie R, et al. (2013) Proc Natl Acad Sci USA 110:4893-4898]. Processing of single cells from the LS174T colorectal cancer cell line showed that D-R mapping of single molecules can reveal structural variation (SV) in the genome of single cells. In one experiment, we processed 17 fragments covering 19.8 Mb of the cell's genome. One megabase-large fragment aligned well to chromosome 19 with half its length, while the other half showed variable alignment. Paired-end single-cell sequencing supported this finding, revealing a region of complexity and a 50-kb deletion. Sequencing struggled, however, to detect a 20-kb gap that D-R mapping showed clearly in a megabase fragment that otherwise mapped well to the reference at the pericentromeric region of chromosome 4. Pericentromeric regions are complex and show substantial sequence homology between different chromosomes, making mapping of sequence reads ambiguous. Thus, D-R mapping directly, from a single molecule, revealed characteristics of the single-cell genome that were challenging for short-read sequencing.

Genetics of the human face: Identification of large-effect single gene variants.

To discover specific variants with relatively large effects on the human face, we have devised an approach to identifying facial features with high heritability. This is based on using twin data to estimate the additive genetic value of each point on a face, as provided by a 3D camera system. In addition, we have used the ethnic difference between East Asian and European faces as a further source of face genetic variation. We use principal components (PCs) analysis to provide a fine definition of the surface features of human faces around the eyes and of the profile, and chose upper and lower 10% extremes of the most heritable PCs for looking for genetic associations. Using this strategy for the analysis of 3D images of 1,832 unique volunteers from the well-characterized People of the British Isles study and 1,567 unique twin images from the TwinsUK cohort, together with genetic data for 500,000 SNPs, we have identified three specific genetic variants with notable effects on facial profiles and eyes.

Somatic selection of poorly differentiating variant stem cell clones could be a key to human ageing.

Any replicating system in which heritable variants with differing replicative potentials can arise is subject to a Darwinian evolutionary process. The continually replicating adult tissue stem cells that control the integrity of many tissues of long-lived, multicellular, complex vertebrate organisms, including humans, constitute such a replicating system. Our suggestion is that somatic selection for mutations (or stable epigenetic changes) that cause an increased rate of adult tissue stem cell proliferation, and their long-term persistence, at the expense of normal differentiation, is a major key to the ageing process. Once an organism has passed the reproductive age, there is no longer any significant counterselection at the organismal level to this inevitable cellular level Darwinian process.

Cellular polarity modulates drug resistance in primary colorectal cancers via orientation of the multidrug resistance protein ABCB1.

Colonic epithelial cells are highly polarised with a lumen-facing apical membrane, termed the brush border, and a basal membrane in contact with the underlying extracellular matrix (ECM). This polarity is often maintained in cancer tissue in the form of neoplastic glands and has prognostic value. We compared the cellular polarity of several ex vivo spheroid colonic cancer cultures with their parental tumours and found that those grown as non-attached colonies exhibited apical brush border proteins on their outer cellular membranes. Transfer of these cultures to an ECM, such as collagen, re-established the centralised apical polarity observed in vivo. The multidrug resistance protein ABCB1 also became aberrantly polarised to outer colony membranes in suspension cultures, unlike cultures grown in collagen, where it was polarised to central lumens. This polarity switch was dependent on the presence of serum or selected serum components, including epidermal growth factor (EGF), transforming growth factor-ß1 (TGF-ß1) and insulin-like growth factor-1 (IGF-1). The apical/basal orientation of primary cancer colon cultures cultured in collagen/serum was modulated by α2ß1 integrin signalling. The polarisation of ABCB1 in colonies significantly altered drug uptake and sensitivity, as the outward polarisation of ABCB1 in suspension colonies effluxed substrates more effectively than ECM-grown colonies with ABCB1 polarised to central lumens. Thus, serum-free suspension colonies were more resistant to a variety of anti-cancer drugs than ECM-grown colonies. In conclusion, the local stroma, or absence thereof, can have profound effects on the sensitivity of colorectal cultures to drugs that are ABCB1 substrates. © 2018 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Stromal uptake and transmission of acid is a pathway for venting cancer cell-generated acid.

Proliferation and invasion of cancer cells require favorable pH, yet potentially toxic quantities of acid are produced metabolically. Membrane-bound transporters extrude acid from cancer cells, but little is known about the mechanisms that handle acid once it is released into the poorly perfused extracellular space. Here, we studied acid handling by myofibroblasts (colon cancer-derived Hs675.T, intestinal InMyoFib, embryonic colon-derived CCD-112-CoN), skin fibroblasts (NHDF-Ad), and colorectal cancer (CRC) cells (HCT116, HT29) grown in monoculture or coculture. Expression of the acid-loading transporter anion exchanger 2 (AE2) (SLC4A2 product) was detected in myofibroblasts and fibroblasts, but not in CRC cells. Compared with CRC cells, Hs675.T and InMyoFib myofibroblasts had very high capacity to absorb extracellular acid. Acid uptake into CCD-112-CoN and NHDF-Ad cells was slower and comparable to levels in CRC cells, but increased alongside SLC4A2 expression under stimulation with transforming growth factor ß1 (TGFß1), a cytokine involved in cancer-stroma interplay. Myofibroblasts and fibroblasts are connected by gap junctions formed by proteins such as connexin-43, which allows the absorbed acid load to be transmitted across the stromal syncytium. To match the stimulatory effect on acid uptake, cell-to-cell coupling in NHDF-Ad and CCD-112-CoN cells was strengthened with TGFß1. In contrast, acid transmission was absent between CRC cells, even after treatment with TGFß1. Thus, stromal cells have the necessary molecular apparatus for assembling an acid-venting route that can improve the flow of metabolic acid through tumors. Importantly, the activities of stromal AE2 and connexin-43 do not place an energetic burden on cancer cells, allowing resources to be diverted for other activities.

Myofibroblasts are distinguished from activated skin fibroblasts by the expression of AOC3 and other associated markers.

Pericryptal myofibroblasts in the colon and rectum play an important role in regulating the normal colorectal stem cell niche and facilitating tumor progression. Myofibroblasts previously have been distinguished from normal fibroblasts mostly by the expression of α smooth muscle actin (αSMA). We now have identified AOC3 (amine oxidase, copper containing 3), a surface monoamine oxidase, as a new marker of myofibroblasts by showing that it is the target protein of the myofibroblast-reacting mAb PR2D3. The normal and tumor tissue distribution and the cell line reactivity of AOC3 match that expected for myofibroblasts. We have shown that the surface expression of AOC3 is sensitive to digestion by trypsin and collagenase and that anti-AOC3 antibodies can be used for FACS sorting of myofibroblasts obtained by nonenzymatic procedures. Whole-genome microarray mRNA-expression profiles of myofibroblasts and skin fibroblasts revealed four additional genes that are significantly differentially expressed in these two cell types: NKX2-3 and LRRC17 in myofibroblasts and SHOX2 and TBX5 in skin fibroblasts. TGFß substantially down-regulated AOC3 expression in myofibroblasts but in skin fibroblasts it dramatically increased the expression of αSMA. A knockdown of NKX2-3 in myofibroblasts caused a decrease of myofibroblast-related gene expression and increased expression of the fibroblast-associated gene SHOX2, suggesting that NKX2-3 is a key mediator for maintaining myofibroblast characteristics. Our results show that colorectal myofibroblasts, as defined by the expression of AOC3, NKX2-3, and other markers, are a distinctly different cell type from TGFß-activated fibroblasts.

A Mathematician's Odyssey.

In this overview of my research, I have aimed to give the background as to how I came to be involved in my various areas of interest, with an emphasis on the early phases of my career, which largely determined my future directions. I had the enormous good fortune to have worked under two of the most outstanding scientists of the twentieth century, R.A. Fisher and Joshua Lederberg. From mathematics and statistics, I went to population genetics and the early use of computers for modeling and simulation. Molecular biology took me into the laboratory and eventually to somatic cell genetics and human gene mapping. One chance encounter led me into the HLA field and another led me into research on cancer, especially colorectal cancer. On the way, I became a champion of the Human Genome Project and of the need for scientists to help promote the public understanding of science.

The CDX1-microRNA-215 axis regulates colorectal cancer stem cell differentiation.

The transcription factor caudal-type homeobox 1 (CDX1) is a key regulator of differentiation in the normal colon and in colorectal cancer (CRC). CDX1 activates the expression of enterocyte genes, but it is not clear how the concomitant silencing of stem cell genes is achieved. MicroRNAs (miRNAs) are important mediators of gene repression and have been implicated in tumor suppression and carcinogenesis, but the roles of miRNAs in differentiation, particularly in CRC, remain poorly understood. Here, we identified microRNA-215 (miR-215) as a direct transcriptional target of CDX1 by using high-throughput small RNA sequencing to profile miRNA expression in two pairs of CRC cell lines: CDX1-low HCT116 and HCT116 with stable CDX1 overexpression, and CDX1-high LS174T and LS174T with stable CDX1 knockdown. Validation of candidate miRNAs identified by RNA-seq in a larger cell-line panel revealed miR-215 to be most significantly correlated with CDX1 expression. Quantitative ChIP-PCR and promoter luciferase assays confirmed that CDX1 directly activates miR-215 transcription. miR-215 expression is depleted in FACS-enriched cancer stem cells compared with unsorted samples. Overexpression of miR-215 in poorly differentiated cell lines causes a decrease in clonogenicity, whereas miR-215 knockdown increases clonogenicity and impairs differentiation in CDX1-high cell lines. We identified the genome-wide targets of miR-215 and found that miR-215 mediates the repression of cell cycle and stemness genes downstream of CDX1. In particular, the miR-215 target gene BMI1 has been shown to promote stemness and self-renewal and to vary inversely with CDX1. Our work situates miR-215 as a link between CDX1 expression and BMI1 repression that governs differentiation in CRC.

Recommendations of the 2006 Human Variome Project meeting.

Lists of variations in genomic DNA and their effects have been kept for some time and have been used in diagnostics and research. Although these lists have been carefully gathered and curated, there has been little standardization and coordination, complicating their use. Given the myriad possible variations in the estimated 24,000 genes in the human genome, it would be useful to have standard criteria for databases of variation. Incomplete collection and ascertainment of variants demonstrates a need for a universally accessible system. These and other problems led to the World Heath Organization-cosponsored meeting on June 20-23, 2006 in Melbourne, Australia, which launched the Human Variome Project. This meeting addressed all areas of human genetics relevant to collection of information on variation and its effects. Members of each of eight sessions (the clinic and phenotype, the diagnostic laboratory, the research laboratory, curation and collection, informatics, relevance to the emerging world, integration and federation and funding and sustainability) developed a number of recommendations that were then organized into a total of 96 recommendations to act as a foundation for future work worldwide. Here we summarize the background of the project, the meeting and its recommendations.

Rapidly derived colorectal cancer cultures recapitulate parental cancer characteristics and enable personalized therapeutic assays.

We have developed a simple procedure for deriving pure cultures of growing cancer cells from colorectal cancers, including material refrigerated overnight, for pathological characterization and cytotoxicity assays. Forty-six cancers were processed and cultures set up under varying culture conditions. Use of a Rho kinase (ROCK1) inhibitor markedly increased culture survival, resulting in 80% of samples growing in culture for at least 1 month and beyond. Overnight refrigeration of samples before culture initiation had little effect on success rates, paving the way for cultures to be established for samples collected over wide geographical areas, such as those for clinical trials. Primary cultures demonstrated good correlation for differentiation markers compared to parent cancers, and were highly dynamic in 3D culture. In Matrigel, many colonies formed central lumens, indicating the presence of stem-like cells. Viable colonies in these cultures recapitulated the in vivo generation of carcinoembryonic antigen (CEA)-positive necrotic/apoptotic debris, much of which was derived from abnormal vacuolated dynamic 'bubble cells' that have not previously been described. Although bubble cells morphologically resembled signet ring cells, a rare cancer subtype, immunostaining suggested that they were most likely derived from terminally differentiated enterocytes. Micro-assays showed that drug toxicity could be measured in these cultures within hours and with sensitivity down to a few hundred cells. Primary cultures derived by our method provide valid in vitro avatars for studying the pathology of cancers in vitro and are amenable to pre-clinical drug testing, paving the way for personalized cancer treatment.

Direct and immune mediated antibody targeting of ERBB receptors in a colorectal cancer cell-line panel.

A significant proportion of colorectal cancer (CRC) patients are resistant to anti-ERBB1 [avian erythroblastic leukemia viral (v-erb-b) oncogene homolog, receptor for EGF] monoclonal antibodies (Mabs). We evaluated both immune and nonimmune effects of cetuximab (anti-ERBB1 Mab), trastuzumab (anti-ERBB2 Mab), pertuzumab (anti-ERBB2 Mab), and lapatinib (dual ERBB1 and ERBB2 tyrosine kinase inhibitor) in a large well-characterized panel of 64 CRC cell lines to find response predictive tumor characteristics. There was a significant correlation between the direct effects of cetuximab and lapatinib. Both agents were associated (P = 0.0004) with "triple' wild-type status in KRAS, BRAF, and PIK3CA exon 20. Most cell lines were resistant to the direct effects of anti-ERBB2 Mabs, suggesting that the effects of lapatinib might mainly be through ERBB1. Microarray mRNA expression profiles of sensitive and resistant cell lines showed that although ERBB1 receptor or ligand levels did not associate with cetuximab sensitivity, high levels of ERBB2 (P = 0.036) and amphiregulin (P = 0.026) predicted sensitivity to lapatinib. However, higher ERBB1 expression predicted susceptibility to cetuximab-induced antibody-dependent cellular cytotoxicity and occurred independently of KRAS/BRAF/PIK3CA mutations (P = 0.69). Lapatinib may be an effective alternative therapy to cetuximab in triple wild-type tumors. Microarray analysis provides suggestive biomarkers for resistance. ERBB1 levels, independent of mutation status, predict immune killing. Therefore, anti-ERBB1 antibodies may be considered in CRC tumors with higher ERBB1 expression and favorable FcγR polymorphisms.

Hypoxia and lineage specification of cell line-derived colorectal cancer stem cells.

Hypoxia is an important regulator of normal and cancer stem cell (CSC) differentiation. Colorectal CSCs from SW1222, LS180, and CCK81 colorectal cancer-derived cell lines are able to differentiate into complex 3D lumen-containing structures in normoxia, whereas in hypoxia, they form undifferentiated dense colonies that have reduced expression of the enterocyte differentiation marker CDX1, lack goblet cell formation, and have increased expression of BMI1 and activated Notch1. Hypoxia increases the clonogenicity of CSCs, which is cumulative as each round of hypoxia enriches for more CSCs. The hypoxic phenotype is reversible, because cells from hypoxic-dense colonies are able to reform differentiated structures when regrown in normoxia. We show that CDX1 is able to stimulate the generation of lumens even in hypoxia and has a negative feedback on BMI1 expression. Knockdown of CDX1 reduces lumen formation but does not affect goblet cell formation, suggesting that enterocytes and goblet cells form from different progenitor cells. Notch inhibition by dibenzazepine (DBZ) allowed CSCs to form goblet cells in both normoxia and hypoxia. Finally, we show that Hif1α, but not CA9, is an important mediator of the effects of hypoxia on the clonogenicity and differentiation of CSCs. In summary, hypoxia maintains the stem-like phenotype of colorectal cell line-derived CSCs and prevents differentiation of enterocytes and goblet cells by regulating CDX1 and Notch1, suggesting that this regulation is an important component of how hypoxia controls the switch between stemness and differentiation in CSCs.

Replication error deficient and proficient colorectal cancer gene expression differences caused by 3'UTR polyT sequence deletions.

Replication error deficient (RER+) colorectal cancers are a distinct subset of colorectal cancers, characterized by inactivation of the DNA mismatch repair system. These cancers are typically pseudodiploid, accumulate mutations in repetitive sequences as a result of their mismatch repair deficiency, and have distinct pathologies. Regulatory sequences controlling all aspects of mRNA processing, especially including message stability, are found in the 3'UTR sequence of most genes. The relevant sequences are typically A/U-rich elements or U repeats. Microarray analysis of 14 RER+ (deficient) and 16 RER- (proficient) colorectal cancer cell lines confirms a striking difference in expression profiles. Analysis of the incidence of mononucleotide repeat sequences in the 3'UTRs, 5'UTRs, and coding sequences of those genes most differentially expressed in RER+ versus RER- cell lines has shown that much of this differential expression can be explained by the occurrence of a massive enrichment of genes with 3'UTR T repeats longer than 11 base pairs in the most differentially expressed genes. This enrichment was confirmed by analysis of two published consensus sets of RER differentially expressed probesets for a large number of primary colorectal cancers. Sequence analysis of the 3'UTRs of a selection of the most differentially expressed genes shows that they all contain deletions in these repeats in all RER+ cell lines studied. These data strongly imply that deregulation of mRNA stability through accumulation of mutations in repetitive regulatory 3'UTR sequences underlies the striking difference in expression profiles between RER+ and RER- colorectal cancers.