Metastatic recurrence in colorectal cancer arises from residual EMP1+ cells.

Around 30-40% of patients with colorectal cancer (CRC) undergoing curative resection of the primary tumour will develop metastases in the subsequent years1. Therapies to prevent disease relapse remain an unmet medical need. Here we uncover the identity and features of the residual tumour cells responsible for CRC relapse. An analysis of single-cell transcriptomes of samples from patients with CRC revealed that the majority of genes associated with a poor prognosis are expressed by a unique tumour cell population that we named high-relapse cells (HRCs). We established a human-like mouse model of microsatellite-stable CRC that undergoes metastatic relapse after surgical resection of the primary tumour. Residual HRCs occult in mouse livers after primary CRC surgery gave rise to multiple cell types over time, including LGR5+ stem-like tumour cells2-4, and caused overt metastatic disease. Using Emp1 (encoding epithelial membrane protein 1) as a marker gene for HRCs, we tracked and selectively eliminated this cell population. Genetic ablation of EMP1high cells prevented metastatic recurrence and mice remained disease-free after surgery. We also found that HRC-rich micrometastases were infiltrated with T cells, yet became progressively immune-excluded during outgrowth. Treatment with neoadjuvant immunotherapy eliminated residual metastatic cells and prevented mice from relapsing after surgery. Together, our findings reveal the cell-state dynamics of residual disease in CRC and anticipate that therapies targeting HRCs may help to avoid metastatic relapse.

TGFß drives immune evasion in genetically reconstituted colon cancer metastasis.

Most patients with colorectal cancer die as a result of the disease spreading to other organs. However, no prevalent mutations have been associated with metastatic colorectal cancers. Instead, particular features of the tumour microenvironment, such as lack of T-cell infiltration, low type 1 T-helper cell (TH1) activity and reduced immune cytotoxicity or increased TGFß levels predict adverse outcomes in patients with colorectal cancer. Here we analyse the interplay between genetic alterations and the tumour microenvironment by crossing mice bearing conditional alleles of four main colorectal cancer mutations in intestinal stem cells. Quadruple-mutant mice developed metastatic intestinal tumours that display key hallmarks of human microsatellite-stable colorectal cancers, including low mutational burden, T-cell exclusion and TGFß-activated stroma. Inhibition of the PD-1-PD-L1 immune checkpoint provoked a limited response in this model system. By contrast, inhibition of TGFß unleashed a potent and enduring cytotoxic T-cell response against tumour cells that prevented metastasis. In mice with progressive liver metastatic disease, blockade of TGFß signalling rendered tumours susceptible to anti-PD-1-PD-L1 therapy. Our data show that increased TGFß in the tumour microenvironment represents a primary mechanism of immune evasion that promotes T-cell exclusion and blocks acquisition of the TH1-effector phenotype. Immunotherapies directed against TGFß signalling may therefore have broad applications in treating patients with advanced colorectal cancer.

Targeting metastasis-initiating cells through the fatty acid receptor CD36.

The fact that the identity of the cells that initiate metastasis in most human cancers is unknown hampers the development of antimetastatic therapies. Here we describe a subpopulation of CD44bright cells in human oral carcinomas that do not overexpress mesenchymal genes, are slow-cycling, express high levels of the fatty acid receptor CD36 and lipid metabolism genes, and are unique in their ability to initiate metastasis. Palmitic acid or a high-fat diet specifically boosts the metastatic potential of CD36+ metastasis-initiating cells in a CD36-dependent manner. The use of neutralizing antibodies to block CD36 causes almost complete inhibition of metastasis in immunodeficient or immunocompetent orthotopic mouse models of human oral cancer, with no side effects. Clinically, the presence of CD36+ metastasis-initiating cells correlates with a poor prognosis for numerous types of carcinomas, and inhibition of CD36 also impairs metastasis, at least in human melanoma- and breast cancer-derived tumours. Together, our results indicate that metastasis-initiating cells particularly rely on dietary lipids to promote metastasis.

GEMC1 is a critical regulator of multiciliated cell differentiation.

The generation of multiciliated cells (MCCs) is required for the proper function of many tissues, including the respiratory tract, brain, and germline. Defects in MCC development have been demonstrated to cause a subclass of mucociliary clearance disorders termed reduced generation of multiple motile cilia (RGMC). To date, only two genes, Multicilin (MCIDAS) and cyclin O (CCNO) have been identified in this disorder in humans. Here, we describe mice lacking GEMC1 (GMNC), a protein with a similar domain organization as Multicilin that has been implicated in DNA replication control. We have found that GEMC1-deficient mice are growth impaired, develop hydrocephaly with a high penetrance, and are infertile, due to defects in the formation of MCCs in the brain, respiratory tract, and germline. Our data demonstrate that GEMC1 is a critical regulator of MCC differentiation and a candidate gene for human RGMC or related disorders.

Distribution of ADAT-Dependent Codons in the Human Transcriptome.

Nucleotide modifications in the anticodons of transfer RNAs (tRNA) play a central role in translation efficiency, fidelity, and regulation of translation, but, for most of these modifications, the details of their function remain unknown. The heterodimeric adenosine deaminases acting on tRNAs (ADAT2-ADAT3, or ADAT) are enzymes present in eukaryotes that convert adenine (A) to inosine (I) in the first anticodon base (position 34) by hydrolytic deamination. To explore the influence of ADAT activity on mammalian translation, we have characterized the human transcriptome and proteome in terms of frequency and distribution of ADAT-related codons. Eight different tRNAs can be modified by ADAT and, once modified, these tRNAs will recognize NNC, NNU and NNA codons, but not NNG codons. We find that transcripts coding for proteins highly enriched in these eight amino acids (ADAT-aa) are specifically enriched in NNC, NNU and NNA codons. We also show that the proteins most enriched in ADAT-aa are composed preferentially of threonine, alanine, proline, and serine (TAPS). We propose that the enrichment in ADAT-codons in these proteins is due to the similarities in the codons that correspond to TAPS.

htSeqTools: high-throughput sequencing quality control, processing and visualization in R.

UNLABELLED: We provide a Bioconductor package with quality assessment, processing and visualization tools for high-throughput sequencing data, with emphasis in ChIP-seq and RNA-seq studies. It includes detection of outliers and biases, inefficient immuno-precipitation and overamplification artifacts, de novo identification of read-rich genomic regions and visualization of the location and coverage of genomic region lists. AVAILABILITY: www.bioconductor.org.

A mathematical framework to determine the temporal sequence of somatic genetic events in cancer.

Human cancer is caused by the accumulation of genetic alterations in cells. Of special importance are changes that occur early during malignant transformation because they may result in oncogene addiction and represent promising targets for therapeutic intervention. Here we describe a computational approach, called Retracing the Evolutionary Steps in Cancer (RESIC), to deduce the temporal sequence of genetic events during tumorigenesis from cross-sectional genomic data of tumors at their fully transformed stage. When applied to a dataset of 70 advanced colorectal cancers, our algorithm accurately predicts the sequence of APC, KRAS, and TP53 mutations previously defined by analyzing tumors at different stages of colon cancer formation. We further validate the method with glioblastoma and leukemia sample data and then apply it to complex integrated genomics databases, finding that high-level EGFR amplification appears to be a late event in primary glioblastomas. RESIC represents the first evolutionary mathematical approach to identify the temporal sequence of mutations driving tumorigenesis and may be useful to guide the validation of candidate genes emerging from cancer genome surveys.

Cellular Senescence Is Immunogenic and Promotes Antitumor Immunity.

Cellular senescence is a stress response that activates innate immune cells, but little is known about its interplay with the adaptive immune system. Here, we show that senescent cells combine several features that render them highly efficient in activating dendritic cells (DC) and antigen-specific CD8 T cells. This includes the release of alarmins, activation of IFN signaling, enhanced MHC class I machinery, and presentation of senescence-associated self-peptides that can activate CD8 T cells. In the context of cancer, immunization with senescent cancer cells elicits strong antitumor protection mediated by DCs and CD8 T cells. Interestingly, this protection is superior to immunization with cancer cells undergoing immunogenic cell death. Finally, the induction of senescence in human primary cancer cells also augments their ability to activate autologous antigen-specific tumor-infiltrating CD8 lymphocytes. Our study indicates that senescent cancer cells can be exploited to develop efficient and protective CD8-dependent antitumor immune responses. SIGNIFICANCE: Our study shows that senescent cells are endowed with a high immunogenic potential-superior to the gold standard of immunogenic cell death. We harness these properties of senescent cells to trigger efficient and protective CD8-dependent antitumor immune responses. See related article by Chen et al., p. 432. This article is highlighted in the In This Issue feature, p. 247.

GEMC1 and MCIDAS interactions with SWI/SNF complexes regulate the multiciliated cell-specific transcriptional program.

Multiciliated cells (MCCs) project dozens to hundreds of motile cilia from their apical surface to promote the movement of fluids or gametes in the mammalian brain, airway or reproductive organs. Differentiation of MCCs requires the sequential action of the Geminin family transcriptional activators, GEMC1 and MCIDAS, that both interact with E2F4/5-DP1. How these factors activate transcription and the extent to which they play redundant functions remains poorly understood. Here, we demonstrate that the transcriptional targets and proximal proteomes of GEMC1 and MCIDAS are highly similar. However, we identified distinct interactions with SWI/SNF subcomplexes; GEMC1 interacts primarily with the ARID1A containing BAF complex while MCIDAS interacts primarily with BRD9 containing ncBAF complexes. Treatment with a BRD9 inhibitor impaired MCIDAS-mediated activation of several target genes and compromised the MCC differentiation program in multiple cell based models. Our data suggest that the differential engagement of distinct SWI/SNF subcomplexes by GEMC1 and MCIDAS is required for MCC-specific transcriptional regulation and mediated by their distinct C-terminal domains.

Metabolic rewiring induced by ranolazine improves melanoma responses to targeted therapy and immunotherapy.

Resistance of melanoma to targeted therapy and immunotherapy is linked to metabolic rewiring. Here, we show that increased fatty acid oxidation (FAO) during prolonged BRAF inhibitor (BRAFi) treatment contributes to acquired therapy resistance in mice. Targeting FAO using the US Food and Drug Administration-approved and European Medicines Agency-approved anti-anginal drug ranolazine (RANO) delays tumour recurrence with acquired BRAFi resistance. Single-cell RNA-sequencing analysis reveals that RANO diminishes the abundance of the therapy-resistant NGFRhi neural crest stem cell subpopulation. Moreover, by rewiring the methionine salvage pathway, RANO enhances melanoma immunogenicity through increased antigen presentation and interferon signalling. Combination of RANO with anti-PD-L1 antibodies strongly improves survival by increasing antitumour immune responses. Altogether, we show that RANO increases the efficacy of targeted melanoma therapy through its effects on FAO and the methionine salvage pathway. Importantly, our study suggests that RANO could sensitize BRAFi-resistant tumours to immunotherapy. Since RANO has very mild side-effects, it might constitute a therapeutic option to improve the two main strategies currently used to treat metastatic melanoma.

A differentiation-based microRNA signature identifies leiomyosarcoma as a mesenchymal stem cell-related malignancy.

Smooth muscle (SM) is a spontaneously contractile tissue that provides physical support and function to organs such as the uterus. Uterine smooth muscle-related neoplasia comprise common well-differentiated benign lesions called leiomyomas (ULM), and rare, highly aggressive and pleomorphic tumors named leiomyosarcomas (ULMS). MicroRNAs (miRNAs) are small non-coding RNAs that play essential roles in normal cellular development and tissue homeostasis that can be used to accurately subclassify different tumor types. Here, we demonstrate that miRNAs are required for full smooth muscle cell (SMC) differentiation of bone marrow-derived human mesenchymal stem cells (hMSCs). We also report a miRNA signature associated with this process. Moreover, we show that this signature, along with miRNA profiles for ULMS and ULM, are able to subclassify tumors of smooth muscle origin along SM differentiation. Using multiple computational analyses, we determined that ULMS are more similar to hMSCs as opposed to ULM, which are linked with more mature SMCs and myometrium. Furthermore, a comparison of the SM differentiation and ULMS miRNA signatures identified miRNAs strictly associated with SM maturation or transformation, as well as those modulated in both processes indicating a possible dual role. These results support separate origins and/or divergent transformation pathways for ULM and ULMS, resulting in drastically different states of differentiation. In summary, this work expands on our knowledge of the regulation of SM differentiation and sarcoma pathogenesis.

Identification and drug-induced reversion of molecular signatures of Alzheimer's disease onset and progression in AppNL-G-F, AppNL-F, and 3xTg-AD mouse models.

BACKGROUND: In spite of many years of research, our understanding of the molecular bases of Alzheimer's disease (AD) is still incomplete, and the medical treatments available mainly target the disease symptoms and are hardly effective. Indeed, the modulation of a single target (e.g., ß-secretase) has proven to be insufficient to significantly alter the physiopathology of the disease, and we should therefore move from gene-centric to systemic therapeutic strategies, where AD-related changes are modulated globally. METHODS: Here we present the complete characterization of three murine models of AD at different stages of the disease (i.e., onset, progression and advanced). We combined the cognitive assessment of these mice with histological analyses and full transcriptional and protein quantification profiling of the hippocampus. Additionally, we derived specific Aß-related molecular AD signatures and looked for drugs able to globally revert them. RESULTS: We found that AD models show accelerated aging and that factors specifically associated with Aß pathology are involved. We discovered a few proteins whose abundance increases with AD progression, while the corresponding transcript levels remain stable, and showed that at least two of them (i.e., lfit3 and Syt11) co-localize with Aß plaques in the brain. Finally, we found two NSAIDs (dexketoprofen and etodolac) and two anti-hypertensives (penbutolol and bendroflumethiazide) that overturn the cognitive impairment in AD mice while reducing Aß plaques in the hippocampus and partially restoring the physiological levels of AD signature genes to wild-type levels. CONCLUSIONS: The characterization of three AD mouse models at different disease stages provides an unprecedented view of AD pathology and how this differs from physiological aging. Moreover, our computational strategy to chemically revert AD signatures has shown that NSAID and anti-hypertensive drugs may still have an opportunity as anti-AD agents, challenging previous reports.

Genomic organization of eukaryotic tRNAs.

BACKGROUND: Surprisingly little is known about the organization and distribution of tRNA genes and tRNA-related sequences on a genome-wide scale. While tRNA gene complements are usually reported in passing as part of genome annotation efforts, and peculiar features such as the tandem arrangements of tRNA gene in Entamoeba histolytica have been described in some detail, systematic comparative studies are rare and mostly restricted to bacteria. We therefore set out to survey the genomic arrangement of tRNA genes and pseudogenes in a wide range of eukaryotes to identify common patterns and taxon-specific peculiarities. RESULTS: In line with previous reports, we find that tRNA complements evolve rapidly and tRNA gene and pseudogene locations are subject to rapid turnover. At phylum level, the distributions of the number of tRNA genes and pseudogenes numbers are very broad, with standard deviations on the order of the mean. Even among closely related species we observe dramatic changes in local organization. For instance, 65% and 87% of the tRNA genes and pseudogenes are located in genomic clusters in zebrafish and stickleback, resp., while such arrangements are relatively rare in the other three sequenced teleost fish genomes. Among basal metazoa, Trichoplax adherens has hardly any duplicated tRNA gene, while the sea anemone Nematostella vectensis boasts more than 17000 tRNA genes and pseudogenes. Dramatic variations are observed even within the eutherian mammals. Higher primates, for instance, have 616 +/- 120 tRNA genes and pseudogenes of which 17% to 36% are arranged in clusters, while the genome of the bushbaby Otolemur garnetti has 45225 tRNA genes and pseudogenes of which only 5.6% appear in clusters. In contrast, the distribution is surprisingly uniform across plant genomes. Consistent with this variability, syntenic conservation of tRNA genes and pseudogenes is also poor in general, with turn-over rates comparable to those of unconstrained sequence elements. Despite this large variation in abundance in Eukarya we observe a significant correlation between the number of tRNA genes, tRNA pseudogenes, and genome size. CONCLUSIONS: The genomic organization of tRNA genes and pseudogenes shows complex lineage-specific patterns characterized by an extensive variability that is in striking contrast to the extreme levels of sequence-conservation of the tRNAs themselves. The comprehensive analysis of the genomic organization of tRNA genes and pseudogenes in Eukarya provides a basis for further studies into the interplay of tRNA gene arrangements and genome organization in general.

Differential Distribution of the wlaN and cgtB Genes, Associated with Guillain-Barré Syndrome, in Campylobacter jejuni Isolates from Humans, Broiler Chickens, and Wild Birds.

Campylobacter jejuni causes campylobacteriosis, a bacterial gastroenteritis with high incidence worldwide. Moreover, C. jejuni infection can trigger the polyneuropathic disorder denominated Guillain-Barré syndrome (GBS). The C. jejuni strains that can elicit GBS carry either wlaN or cgtB, coding both genes for a ß-1,3-galactosyltransferase enzyme that is required for the production of sialylated lipooligosaccharide (LOSSIAL). We described a differential prevalence of the genes wlaN and cgtB in C. jejuni isolates from three different ecological niches: humans, broiler chickens, and wild birds. The distribution of both genes, which is similar between broiler chicken and human isolates and distinct when compared to the wild bird isolates, suggests a host-dependent distribution. Moreover, the prevalence of the wlaN and cgtB genes seems to be restricted to some clonal complexes. Gene sequencing identified the presence of new variants of the G- homopolymeric tract within the wlaN gene. Furthermore, we detected two variants of a G rich region within the cgtB gene, suggesting that, similarly to wlaN, the G-tract in the cgtB gene mediates the phase variation control of cgtB expression. Caco-2 cell invasion assays indicate that there is no evident correlation between the production of LOSSIAL and the ability to invade eukaryotic cells.

Tousled-Like Kinases Suppress Innate Immune Signaling Triggered by Alternative Lengthening of Telomeres.

The Tousled-like kinases 1 and 2 (TLK1/2) control histone deposition through the ASF1 histone chaperone and influence cell cycle progression and genome maintenance, yet the mechanisms underlying TLK-mediated genome stability remain uncertain. Here, we show that TLK loss results in severe chromatin decompaction and altered genome accessibility, particularly affecting heterochromatic regions. Failure to maintain heterochromatin increases spurious transcription of repetitive elements and induces features of alternative lengthening of telomeres (ALT). TLK depletion culminates in a cGAS-STING-TBK1-mediated innate immune response that is independent of replication-stress signaling and attenuated by the depletion of factors required to produce extra-telomeric DNA. Analysis of human cancers reveals that chromosomal instability correlates with high TLK2 and low STING levels in many cohorts. Based on these findings, we propose that high TLK levels contribute to immune evasion in chromosomally unstable and ALT+ cancers.

Eleven grand challenges in single-cell data science.

The recent boom in microfluidics and combinatorial indexing strategies, combined with low sequencing costs, has empowered single-cell sequencing technology. Thousands-or even millions-of cells analyzed in a single experiment amount to a data revolution in single-cell biology and pose unique data science problems. Here, we outline eleven challenges that will be central to bringing this emerging field of single-cell data science forward. For each challenge, we highlight motivating research questions, review prior work, and formulate open problems. This compendium is for established researchers, newcomers, and students alike, highlighting interesting and rewarding problems for the coming years.

Homology-based annotation of non-coding RNAs in the genomes of Schistosoma mansoni and Schistosoma japonicum.

BACKGROUND: Schistosomes are trematode parasites of the phylum Platyhelminthes. They are considered the most important of the human helminth parasites in terms of morbidity and mortality. Draft genome sequences are now available for Schistosoma mansoni and Schistosoma japonicum. Non-coding RNA (ncRNA) plays a crucial role in gene expression regulation, cellular function and defense, homeostasis, and pathogenesis. The genome-wide annotation of ncRNAs is a non-trivial task unless well-annotated genomes of closely related species are already available. RESULTS: A homology search for structured ncRNA in the genome of S. mansoni resulted in 23 types of ncRNAs with conserved primary and secondary structure. Among these, we identified rRNA, snRNA, SL RNA, SRP, tRNAs and RNase P, and also possibly MRP and 7SK RNAs. In addition, we confirmed five miRNAs that have recently been reported in S. japonicum and found two additional homologs of known miRNAs. The tRNA complement of S. mansoni is comparable to that of the free-living planarian Schmidtea mediterranea, although for some amino acids differences of more than a factor of two are observed: Leu, Ser, and His are overrepresented, while Cys, Meth, and Ile are underrepresented in S. mansoni. On the other hand, the number of tRNAs in the genome of S. japonicum is reduced by more than a factor of four. Both schistosomes have a complete set of minor spliceosomal snRNAs. Several ncRNAs that are expected to exist in the S. mansoni genome were not found, among them the telomerase RNA, vault RNAs, and Y RNAs. CONCLUSION: The ncRNA sequences and structures presented here represent the most complete dataset of ncRNA from any lophotrochozoan reported so far. This data set provides an important reference for further analysis of the genomes of schistosomes and indeed eukaryotic genomes at large.

Drosophila Larval Brain Neoplasms Present Tumour-Type Dependent Genome Instability.

Single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) are found at different rates in human cancer. To determine if these genetic lesions appear in Drosophila tumors we have sequenced the genomes of 17 malignant neoplasms caused by mutations in l(3)mbt, brat, aurA, or lgl We have found CNVs and SNPs in all the tumors. Tumor-linked CNVs range between 11 and 80 per sample, affecting between 92 and 1546 coding sequences. CNVs are in average less frequent in l(3)mbt than in brat lines. Nearly half of the CNVs fall within the 10 to 100Kb range, all tumor samples contain CNVs larger that 100 Kb and some have CNVs larger than 1Mb. The rates of tumor-linked SNPs change more than 20-fold depending on the tumor type: at late time points brat, l(3)mbt, and aurA and lgl lines present median values of SNPs/Mb of exome of 0.16, 0.48, and 3.6, respectively. Higher SNP rates are mostly accounted for by C > A transversions, which likely reflect enhanced oxidative stress conditions in the affected tumors. Both CNVs and SNPs turn over rapidly. We found no evidence for selection of a gene signature affected by CNVs or SNPs in the cohort. Altogether, our results show that the rates of CNVs and SNPs, as well as the distribution of CNV sizes in this cohort of Drosophila tumors are well within the range of those reported for human cancer. Genome instability is therefore inherent to Drosophila malignant neoplastic growth at a variable extent that is tumor type dependent.

Tousled-like kinases stabilize replication forks and show synthetic lethality with checkpoint and PARP inhibitors.

DNA sequence and epigenetic information embedded in chromatin must be faithfully duplicated and transmitted to daughter cells during cell division. However, how chromatin assembly and DNA replication are integrated remains unclear. We examined the contribution of the Tousled-like kinases 1 and 2 (TLK1/TLK2) to chromatin assembly and maintenance of replication fork integrity. We show that TLK activity is required for DNA replication and replication-coupled nucleosome assembly and that lack of TLK activity leads to replication fork stalling and the accumulation of single-stranded DNA, a phenotype distinct from ASF1 depletion. Consistent with these results, sustained TLK depletion gives rise to replication-dependent DNA damage and p53-dependent cell cycle arrest in G1. We find that deficient replication-coupled de novo nucleosome assembly renders replication forks unstable and highly dependent on the ATR and CHK1 checkpoint kinases, as well as poly(adenosine 5'-diphosphate-ribose) polymerase (PARP) activity, to avoid collapse. Human cancer data revealed frequent up-regulation of TLK genes and an association with poor patient outcome in multiple types of cancer, and depletion of TLK activity leads to increased replication stress and DNA damage in a panel of cancer cells. Our results reveal a critical role for TLKs in chromatin replication and suppression of replication stress and identify a synergistic lethal relationship with checkpoint signaling and PARP that could be exploited in treatment of a broad range of cancers.

Mex3a Marks a Slowly Dividing Subpopulation of Lgr5+ Intestinal Stem Cells.

Highly proliferative Lgr5+ stem cells maintain the intestinal epithelium and are thought to be largely homogeneous. Although quiescent intestinal stem cell (ISC) populations have been described, the identity and features of such a population remain controversial. Here we report unanticipated heterogeneity within the Lgr5+ ISC pool. We found that expression of the RNA-binding protein Mex3a labels a slowly cycling subpopulation of Lgr5+ ISCs that contribute to all intestinal lineages with distinct kinetics. Single-cell transcriptome profiling revealed that Lgr5+ cells adopt two discrete states, one of which is defined by a Mex3a expression program and relatively low levels of proliferation genes. During homeostasis, Mex3a+ cells continually shift into the rapidly dividing, self-renewing ISC pool. Chemotherapy and radiation preferentially target rapidly dividing Lgr5+ cells but spare the Mex3a-high/Lgr5+ population, helping to promote regeneration of the intestinal epithelium following toxic insults. Thus, Mex3a defines a reserve-like ISC population within the Lgr5+ compartment.