Double-negative B cells and DNASE1L3 colocalise with microbiota in gut-associated lymphoid tissue

Intestinal homeostasis is maintained by the response of gut-associated lymphoid tissue to bacteria transported across the follicle associated epithelium into the subepithelial dome. The initial response to antigens and how bacteria are handled is incompletely understood. By iterative application of spatial transcriptomics and multiplexed single-cell technologies, we identify that the double negative 2 subset of B cells, previously associated with autoimmune diseases, is present in the subepithelial dome in health. We show that in this location double negative 2 B cells interact with dendritic cells co-expressing the lupus autoantigens DNASE1L3 and C1q and microbicides. We observe that in humans, but not in mice, dendritic cells expressing DNASE1L3 are associated with sampled bacteria but not DNA derived from apoptotic cells. We propose that fundamental features of autoimmune diseases are microbiota-associated, interacting components of normal intestinal immunity.

Apoptosis in mesenchymal stromal cells activates an immunosuppressive secretome predicting clinical response in Crohn's disease

In vivo apoptosis of human mesenchymal stromal cells (MSCs) plays a critical role in delivering immunomodulation. Yet, caspase activity not only mediates the dying process but also death-independent functions that may shape the immunogenicity of apoptotic cells. Therefore, a better characterization of the immunological profile of apoptotic MSCs (ApoMSCs) could shed light on their mechanistic action and therapeutic applications. We analyzed the transcriptomes of MSCs undergoing apoptosis and identified several immunomodulatory factors and chemokines dependent on caspase activation following Fas stimulation. The ApoMSC secretome inhibited human T cell proliferation and activation, and chemoattracted monocytes in vitro. Both immunomodulatory activities were dependent on the cyclooxygenase2 (COX2)/prostaglandin E2 (PGE2) axis. To assess the clinical relevance of ApoMSC signature, we used the peripheral blood mononuclear cells (PBMCs) from a cohort of fistulizing Crohn's disease (CD) patients who had undergone MSC treatment (ADMIRE-CD). Compared with healthy donors, MSCs exposed to patients' PBMCs underwent apoptosis and released PGE2 in a caspase-dependent manner. Both PGE2 and apoptosis were significantly associated with clinical responses to MSCs. Our findings identify a new mechanism whereby caspase activation delivers ApoMSC immunosuppression. Remarkably, such molecular signatures could implicate translational tools for predicting patients' clinical responses to MSC therapy in CD.

The transcriptional landscape of endogenous retroelements delineates esophageal adenocarcinoma subtypes

Most cancer types exhibit aberrant transcriptional activity, including derepression of retrotransposable elements (RTEs). However, the degree, specificity and potential consequences of RTE transcriptional activation may differ substantially among cancer types and subtypes. Representing one extreme of the spectrum, we characterize the transcriptional activity of RTEs in cohorts of esophageal adenocarcinoma (EAC) and its precursor Barrett's esophagus (BE) from the OCCAMS (Oesophageal Cancer Clinical and Molecular Stratification) consortium, and from TCGA (The Cancer Genome Atlas). We found exceptionally high RTE inclusion in the EAC transcriptome, driven primarily by transcription of genes incorporating intronic or adjacent RTEs, rather than by autonomous RTE transcription. Nevertheless, numerous chimeric transcripts straddling RTEs and genes, and transcripts from stand-alone RTEs, particularly KLF5- and SOX9-controlled HERVH proviruses, were overexpressed specifically in EAC. Notably, incomplete mRNA splicing and EAC-characteristic intronic RTE inclusion was mirrored by relative loss of the respective fully-spliced, functional mRNA isoforms, consistent with compromised cellular fitness. Defective RNA splicing was linked with strong transcriptional activation of a HERVH provirus on Chr Xp22.32 and defined EAC subtypes with distinct molecular features and prognosis. Our study defines distinguishable RTE transcriptional profiles of EAC, reflecting distinct underlying processes and prognosis, thus providing a framework for targeted studies.

Mechanistic insights into the interactions between cancer drivers and the tumour immune microenvironment

BACKGROUND

The crosstalk between cancer and the tumour immune microenvironment (TIME) has attracted significant interest in the latest years because of its impact on cancer evolution and response to treatment. Despite this, cancer-specific tumour-TIME interactions and their mechanistic insights are still poorly understood.

METHODS

Here, we compute the significant interactions occurring between cancer-specific genetic drivers and five anti- and pro-tumour TIME features in 32 cancer types using Lasso regularised ordinal regression. Focusing on head and neck squamous cancer (HNSC), we rebuild the functional networks linking specific TIME driver alterations to the TIME state they associate with.

RESULTS

The 477 TIME drivers that we identify are multifunctional genes whose alterations are selected early in cancer evolution and recur across and within cancer types. Tumour suppressors and oncogenes have an opposite effect on the TIME and the overall anti-tumour TIME driver burden is predictive of response to immunotherapy. TIME driver alterations predict the immune profiles of HNSC molecular subtypes, and perturbations in keratinization, apoptosis and interferon signalling underpin specific driver-TIME interactions.

CONCLUSIONS

Overall, our study delivers a comprehensive resource of TIME drivers, gives mechanistic insights into their immune-regulatory role, and provides an additional framework for patient prioritisation to immunotherapy. The full list of TIME drivers and associated properties are available at http://www.network-cancer-genes.org .

Reflected stemness as a potential driver of the tumour microenvironment

A fundamental requirement for cancer initiation is the activation of developmental programmes by mutant cells. Oncogenic signals often confer an undifferentiated, stem cell-like phenotype that supports the long-term proliferative potential of cancer cells. Although cancer is a genetically driven disease, mutations in cancer-driver genes alone are insufficient for tumour formation, and the proliferation of cells harbouring oncogenic mutations depends on their microenvironment. In this Opinion article we discuss how the reprogrammed status of cancer cells not only represents the essence of their tumorigenicity but triggers 'reflected stemness' in their surrounding normal counterparts. We propose that this reciprocal interaction underpins the establishment of the tumour microenvironment (TME).

Exome sequencing of hepatocellular carcinoma in lemurs identifies potential cancer drivers: A pilot study

Background and objectives

Hepatocellular carcinoma occurs frequently in prosimians, but the cause of these liver cancers in this group is unknown. Characterizing the genetic changes associated with hepatocellular carcinoma in prosimians may point to possible causes, treatments and methods of prevention, aiding conservation efforts that are particularly crucial to the survival of endangered lemurs. Although genomic studies of cancer in non-human primates have been hampered by a lack of tools, recent studies have demonstrated the efficacy of using human exome capture reagents across primates.

Methodology

In this proof-of-principle study, we applied human exome capture reagents to tumor-normal pairs from five lemurs with hepatocellular carcinoma to characterize the mutational landscape of this disease in lemurs.

Results

Several genes implicated in human hepatocellular carcinoma, including ARID1A, TP53 and CTNNB1, were mutated in multiple lemurs, and analysis of cancer driver genes mutated in these samples identified enrichment of genes involved with TP53 degradation and regulation. In addition to these similarities with human hepatocellular carcinoma, we also noted unique features, including six genes that contain mutations in all five lemurs. Interestingly, these genes are infrequently mutated in human hepatocellular carcinoma, suggesting potential differences in the etiology and/or progression of this cancer in lemurs and humans.

Conclusions and implications

Collectively, this pilot study suggests that human exome capture reagents are a promising tool for genomic studies of cancer in lemurs and other non-human primates.

Lay Summary

Hepatocellular carcinoma occurs frequently in prosimians, but the cause of these liver cancers is unknown. In this proof-of-principle study, we applied human DNA sequencing tools to tumor-normal pairs from five lemurs with hepatocellular carcinoma and compared the lemur mutation profiles to those of human hepatocellular carcinomas.

Comparative assessment of genes driving cancer and somatic evolution in non-cancer tissues: an update of the Network of Cancer Genes (NCG) resource

Background

Genetic alterations of somatic cells can drive non-malignant clone formation and promote cancer initiation. However, the link between these processes remains unclear and hampers our understanding of tissue homeostasis and cancer development.

Results

Here, we collect a literature-based repertoire of 3355 well-known or predicted drivers of cancer and non-cancer somatic evolution in 122 cancer types and 12 non-cancer tissues. Mapping the alterations of these genes in 7953 pan-cancer samples reveals that, despite the large size, the known compendium of drivers is still incomplete and biased towards frequently occurring coding mutations. High overlap exists between drivers of cancer and non-cancer somatic evolution, although significant differences emerge in their recurrence. We confirm and expand the unique properties of drivers and identify a core of evolutionarily conserved and essential genes whose germline variation is strongly counter-selected. Somatic alteration in even one of these genes is sufficient to drive clonal expansion but not malignant transformation.

Conclusions

Our study offers a comprehensive overview of our current understanding of the genetic events initiating clone expansion and cancer revealing significant gaps and biases that still need to be addressed. The compendium of cancer and non-cancer somatic drivers, their literature support, and properties are accessible in the Network of Cancer Genes and Healthy Drivers resource at http://www.network-cancer-genes.org/.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02607-z.

Clone competition as a mechanism to reduce tumor formation

With age, clones carrying somatic mutations in well-known cancer driver genes progressively populate adult tissues, yet cancer transformation is rare. In a recent issue of Nature, Colom et al. showed that competition between mutated clones with different fitness could act as a tumor-protective mechanism.

Immunogenomics of Colorectal Cancer Response to Checkpoint Blockade: Analysis of the KEYNOTE 177 Trial and Validation Cohorts

BACKGROUND & AIMS

Colorectal cancer (CRC) shows variable response to immune checkpoint blockade, which can only partially be explained by high tumor mutational burden (TMB). We conducted an integrated study of the cancer tissue and associated tumor microenvironment (TME) from patients treated with pembrolizumab (KEYNOTE 177 clinical trial) or nivolumab to dissect the cellular and molecular determinants of response to anti- programmed cell death 1 (PD1) immunotherapy.

METHODS

We selected multiple regions per tumor showing variable T-cell infiltration for a total of 738 regions from 29 patients, divided into discovery and validation cohorts. We performed multiregional whole-exome and RNA sequencing of the tumor cells and integrated these with T-cell receptor sequencing, high-dimensional imaging mass cytometry, detection of programmed death-ligand 1 (PDL1) interaction in situ, multiplexed immunofluorescence, and computational spatial analysis of the TME.

RESULTS

In hypermutated CRCs, response to anti-PD1 immunotherapy was not associated with TMB but with high clonality of immunogenic mutations, clonally expanded T cells, low activation of Wnt signaling, deregulation of the interferon gamma pathway, and active immune escape mechanisms. Responsive hypermutated CRCs were also rich in cytotoxic and proliferating PD1+CD8 T cells interacting with PDL1+ antigen-presenting macrophages.

CONCLUSIONS

Our study clarified the limits of TMB as a predictor of response of CRC to anti-PD1 immunotherapy. It identified a population of antigen-presenting macrophages interacting with CD8 T cells that consistently segregate with response. We therefore concluded that anti-PD1 agents release the PD1-PDL1 interaction between CD8 T cells and macrophages to promote cytotoxic antitumor activity.

O-linked mucin-type glycosylation regulates the transcriptional programme downstream of EGFR

Aberrant mucin-type O-linked glycosylation is a common occurrence in cancer where the upregulation of sialyltransferases is often seen leading to the early termination of O-glycan chains. Mucin-type O-linked glycosylation is not limited to mucins and occurs on many cell surface glycoproteins including EGFR, where the number of sites can be limited. Upon EGF ligation, EGFR induces a signaling cascade and may also translocate to the nucleus where it directly regulates gene transcription, a process modulated by Galectin-3 and MUC1 in some cancers. Here, we show that upon EGF binding, breast cancer cells carrying different O-glycans respond by transcribing different gene expression signatures. MMP10, the principal gene upregulated when cells carrying sialylated core 1 glycans were stimulated with EGF, is also upregulated in ER-positive breast carcinoma reported to express high levels of ST3Gal1 and hence mainly core 1 sialylated O-glycans. In contrast, isogenic cells engineered to carry core 2 glycans upregulate CX3CL1 and FGFBP1 and these genes are upregulated in ER-negative breast carcinomas, also known to express longer core 2 O-glycans. Changes in O-glycosylation did not significantly alter signal transduction downstream of EGFR in core 1 or core 2 O-glycan expressing cells. However, striking changes were observed in the formation of an EGFR/galectin-3/MUC1/β-catenin complex at the cell surface that is present in cells carrying short core 1-based O-glycans but absent in core 2 carrying cells.

Normal Somatic Mutations in Cancer Transformation

Gene alterations play a prominent role in driving cancer initiation and progression. However, the genetic events that occur in normal cells prior to tumorigenesis are still unknown. Recent studies have started to map somatic mutations in normal human tissues, and here we discuss their implications for our understanding of tumorigenesis.

Pan-cancer detection of driver genes at the single-patient resolution

BACKGROUND

Identifying the complete repertoire of genes that drive cancer in individual patients is crucial for precision oncology. Most established methods identify driver genes that are recurrently altered across patient cohorts. However, mapping these genes back to patients leaves a sizeable fraction with few or no drivers, hindering our understanding of cancer mechanisms and limiting the choice of therapeutic interventions.

RESULTS

We present sysSVM2, a machine learning software that integrates cancer genetic alterations with gene systems-level properties to predict drivers in individual patients. Using simulated pan-cancer data, we optimise sysSVM2 for application to any cancer type. We benchmark its performance on real cancer data and validate its applicability to a rare cancer type with few known driver genes. We show that drivers predicted by sysSVM2 have a low false-positive rate, are stable and disrupt well-known cancer-related pathways.

CONCLUSIONS

sysSVM2 can be used to identify driver alterations in patients lacking sufficient canonical drivers or belonging to rare cancer types for which assembling a large enough cohort is challenging, furthering the goals of precision oncology. As resources for the community, we provide the code to implement sysSVM2 and the pre-trained models in all TCGA cancer types ( https://github.com/ciccalab/sysSVM2 ).

Approaching Cancer Evolution from Different Angles

Dr Francesca Ciccarelli (The Francis Crick Institute, UK) and Dr James De Gregori (University of Colorado, USA) interview 3 top scientists in clinical (Dr Charles Swanton, The Francis Crick Institute, UK), molecular (Dr Kornelia Polyak, Dana-Farber Cancer Institute, USA), and evolutionary cancer research (Dr Carlo Maley, Arizona State University, USA) to discuss the current status of knowledge, the challenges, and the opportunities to move the field forward.

Identification of non-cancer cells from cancer transcriptomic data

Interactions between cancer cells and non-cancer cells composing the tumour microenvironment play a primary role in determining cancer progression and shaping the response to therapy. The qualitative and quantitative characterisation of the different cell populations in the tumour microenvironment is therefore crucial to understand its role in cancer. In recent years, many experimental and computational approaches have been developed to identify the cell populations composing heterogeneous tissue samples, such as cancer. In this review, we describe the state-of-the-art approaches for the quantification of non-cancer cells from bulk and single-cell cancer transcriptomic data, with a focus on immune cells. We illustrate the main features of these approaches and highlight their applications for the analysis of the tumour microenvironment in solid cancers. We also discuss techniques that are complementary and alternative to RNA sequencing, particularly focusing on approaches that can provide spatial information on the distribution of the cells within the tumour in addition to their qualitative and quantitative measurements. This article is part of a Special Issue entitled: Transcriptional Profiles and Regulatory Gene Networks edited by Dr. Federico Manuel Giorgi and Dr. Shaun Mahony.

Mapping the Human Kinome in Response to DNA Damage

We provide a catalog for the effects of the human kinome on cell survival in response to DNA-damaging agents, covering all major DNA repair pathways. By treating 313 kinase-deficient cell lines with ten diverse DNA-damaging agents, including seven commonly used chemotherapeutics, we identified examples of vulnerability and resistance that are kinase specific. To investigate synthetic lethal interactions, we tested the response to carmustine for 25 cell lines by establishing a phenotypic fluorescence-activated cell sorting (FACS) assay designed to validate gene-drug interactions. We show apoptosis, cell cycle changes, and DNA damage and proliferation after alkylation- or crosslink-induced damage. In addition, we reconstitute the cellular sensitivity of DYRK4, EPHB6, MARK3, and PNCK as a proof of principle for our study. Furthermore, using global phosphoproteomics on cells lacking MARK3, we provide evidence for its role in the DNA damage response. Our data suggest that cancers with inactivating mutations in kinases, including MARK3, are particularly vulnerable to alkylating chemotherapeutic agents.

Molecular Biology and Evolution of Cancer: From Discovery to Action

Cancer progression is an evolutionary process. During this process, evolving cancer cell populations encounter restrictive ecological niches within the body, such as the primary tumor, circulatory system, and diverse metastatic sites. Efforts to prevent or delay cancer evolution-and progression-require a deep understanding of the underlying molecular evolutionary processes. Herein we discuss a suite of concepts and tools from evolutionary and ecological theory that can inform cancer biology in new and meaningful ways. We also highlight current challenges to applying these concepts, and propose ways in which incorporating these concepts could identify new therapeutic modes and vulnerabilities in cancer.

IL-36 Promotes Systemic IFN-I Responses in Severe Forms of Psoriasis

Psoriasis is an immune-mediated skin disorder associated with severe systemic comorbidities. Whereas IL-36 is a key disease driver, the pathogenic role of this cytokine has mainly been investigated in skin. Thus, its effects on systemic immunity and extracutaneous disease manifestations remain poorly understood. To address this issue, we investigated the consequences of excessive IL-36 activity in circulating immune cells. We initially focused our attention on generalized pustular psoriasis (GPP), a clinical variant associated with pervasive upregulation of IL-36 signaling. By undertaking blood and neutrophil RNA sequencing, we demonstrated that affected individuals display a prominent IFN-I signature, which correlates with abnormal IL-36 activity. We then validated the association between IL-36 deregulation and IFN-I over-expression in patients with severe psoriasis vulgaris (PV). We also found that the activation of IFN-I genes was associated with extracutaneous morbidity, in both GPP and PV. Finally, we undertook mechanistic experiments, demonstrating that IL-36 acts directly on plasmacytoid dendritic cells, where it potentiates toll-like receptor (TLR)-9 activation and IFN-α production. This effect was mediated by the upregulation of PLSCR1, a phospholipid scramblase mediating endosomal TLR-9 translocation. These findings identify an IL-36/ IFN-I axis contributing to extracutaneous inflammation in psoriasis.

Patient-specific cancer genes contribute to recurrently perturbed pathways and establish therapeutic vulnerabilities in esophageal adenocarcinoma

The identification of cancer-promoting genetic alterations is challenging particularly in highly unstable and heterogeneous cancers, such as esophageal adenocarcinoma (EAC). Here we describe a machine learning algorithm to identify cancer genes in individual patients considering all types of damaging alterations simultaneously. Analysing 261 EACs from the OCCAMS Consortium, we discover helper genes that, alongside well-known drivers, promote cancer. We confirm the robustness of our approach in 107 additional EACs. Unlike recurrent alterations of known drivers, these cancer helper genes are rare or patient-specific. However, they converge towards perturbations of well-known cancer processes. Recurrence of the same process perturbations, rather than individual genes, divides EACs into six clusters differing in their molecular and clinical features. Experimentally mimicking the alterations of predicted helper genes in cancer and pre-cancer cells validates their contribution to disease progression, while reverting their alterations reveals EAC acquired dependencies that can be exploited in therapy.

The Network of Cancer Genes (NCG): a comprehensive catalogue of known and candidate cancer genes from cancer sequencing screens

The Network of Cancer Genes (NCG) is a manually curated repository of 2372 genes whose somatic modifications have known or predicted cancer driver roles. These genes were collected from 275 publications, including two sources of known cancer genes and 273 cancer sequencing screens of more than 100 cancer types from 34,905 cancer donors and multiple primary sites. This represents a more than 1.5-fold content increase compared to the previous version. NCG also annotates properties of cancer genes, such as duplicability, evolutionary origin, RNA and protein expression, miRNA and protein interactions, and protein function and essentiality. NCG is accessible at http://ncg.kcl.ac.uk/ .

An analysis of IL-36 signature genes and individuals with IL1RL2 knockout mutations validates IL-36 as a psoriasis therapeutic target

Interleukin (IL)-36α, IL-36β, and IL-36γ are innate mediators of acute epithelial inflammation. We sought to demonstrate that these cytokines are also required for the pathogenesis of plaque psoriasis, a common and chronic skin disorder, caused by abnormal T helper 17 (T_H17) cell activation. To investigate this possibility, we first defined the genes that are induced by IL-36 cytokines in primary human keratinocytes. This enabled us to demonstrate a significant IL-36 signature among the transcripts that are up-regulated in plaque psoriasis and the susceptibility loci associated with the disease in genome-wide studies. Next, we investigated the impact of in vivo and ex vivo IL-36 receptor blockade using a neutralizing antibody or a recombinant antagonist. Both inhibitors had marked anti-inflammatory effects on psoriatic skin, demonstrated by statistically significant reductions in IL-17 expression, keratinocyte activation, and leukocyte infiltration. Finally, we explored the potential safety profile associated with IL-36 blockade by phenotyping 12 individuals carrying knockout mutations of the IL-36 receptor gene. We found that normal immune function was broadly preserved in these individuals, suggesting that IL-36 signaling inhibition would not substantially compromise host defenses. These observations, which integrate the results of transcriptomics and model system analysis, pave the way for early-stage clinical trials of IL-36 antagonists.

Synthetic lethal interaction between the tumour suppressor STAG2 and its paralog STAG1

Cohesin is a multi-protein complex that tethers sister chromatids during mitosis and mediates DNA repair, genome compartmentalisation and regulation of gene expression. Cohesin subunits frequently acquire cancer loss-of-function alterations and act as tumour suppressors in several tumour types. This has led to increased interest in cohesin as potential target in anti-cancer therapy. Here we show that the loss-of-function of STAG2, a core component of cohesin and an emerging tumour suppressor, leads to synthetic dependency of mutated cancer cells on its paralog STAG1. STAG1 and STAG2 share high sequence identity, encode mutually exclusive cohesin subunits and retain partially overlapping functions. We inhibited STAG1 and STAG2 in several cancer cell lines where the two genes have variable mutation and copy number status. In all cases, we observed that the simultaneous blocking of STAG1 and STAG2 significantly reduces cell proliferation. We further confirmed the synthetic lethal interaction developing a vector-free CRISPR system to induce STAG1/STAG2 double gene knockout. We provide strong evidence that STAG1 is a promising therapeutic target in cancers with inactivating alterations of STAG2.

MEGA-V: detection of variant gene sets in patient cohorts

Summary

Detecting significant associations between genetic variants and disease may prove particularly challenging when the variants are rare in the population and/or act together with other variants to cause the disease. We have developed a statistical framework named Mutation Enrichment Gene set Analysis of Variants (MEGA-V) that specifically detects the enrichments of genetic alterations within a process in a cohort of interest. By focusing on the mutations of several genes contributing to the same function rather than on those affecting a single gene, MEGA-V increases the power to detect statistically significant associations.

Availability and Implementation

MEGA-V is available at https://github.com/ciccalab/MEGA

Supplementary information

Supplementary data are available at Bioinformatics online.

Genetic Redundancy, Functional Compensation, and Cancer Vulnerability

Cancer genomes acquire somatic alterations that largely differ between and within cancer types. Several of these alterations inactivate genes that are normally functional with no deleterious consequences on cancer cells due to genetic redundancy. Here we discuss how this leads to cancer synthetic dependencies that can be exploited in therapy.

NCG 5.0: updates of a manually curated repository of cancer genes and associated properties from cancer mutational screenings

The Network of Cancer Genes (NCG, http://ncg.kcl.ac.uk/) is a manually curated repository of cancer genes derived from the scientific literature. Due to the increasing amount of cancer genomic data, we have introduced a more robust procedure to extract cancer genes from published cancer mutational screenings and two curators independently reviewed each publication. NCG release 5.0 (August 2015) collects 1571 cancer genes from 175 published studies that describe 188 mutational screenings of 13 315 cancer samples from 49 cancer types and 24 primary sites. In addition to collecting cancer genes, NCG also provides information on the experimental validation that supports the role of these genes in cancer and annotates their properties (duplicability, evolutionary origin, expression profile, function and interactions with proteins and miRNAs).

Deep sequencing of the X chromosome reveals the proliferation history of colorectal adenomas

BACKGROUND

Mismatch repair deficient colorectal adenomas are composed of transformed cells that descend from a common founder and progressively accumulate genomic alterations. The proliferation history of these tumors is still largely unknown. Here we present a novel approach to rebuild the proliferation trees that recapitulate the history of individual colorectal adenomas by mapping the progressive acquisition of somatic point mutations during tumor growth.

RESULTS

Using our approach, we called high and low frequency mutations acquired in the X chromosome of four mismatch repair deficient colorectal adenomas deriving from male individuals. We clustered these mutations according to their frequencies and rebuilt the proliferation trees directly from the mutation clusters using a recursive algorithm. The trees of all four lesions were formed of a dominant subclone that co-existed with other genetically heterogeneous subpopulations of cells. However, despite this similar hierarchical organization, the growth dynamics varied among and within tumors, likely depending on a combination of tumor-specific genetic and environmental factors.

CONCLUSIONS

Our study provides insights into the biological properties of individual mismatch repair deficient colorectal adenomas that may influence their growth and also the response to therapy. Extended to other solid tumors, our novel approach could inform on the mechanisms of cancer progression and on the best treatment choice.

Massive gene amplification drives paediatric hepatocellular carcinoma caused by bile salt export pump deficiency

Hepatocellular carcinoma (HCC) is almost invariably associated with an underlying inflammatory state, whose direct contribution to the acquisition of critical genomic changes is unclear. Here we map acquired genomic alterations in human and mouse HCCs induced by defects in hepatocyte biliary transporters, which expose hepatocytes to bile salts and cause chronic inflammation that develops into cancer. In both human and mouse cancer genomes, we find few somatic point mutations with no impairment of cancer genes, but massive gene amplification and rearrangements. This genomic landscape differs from that of virus- and alcohol-associated liver cancer. Copy-number gains preferentially occur at late stages of cancer development and frequently target the MAPK signalling pathway, and in particular direct regulators of JNK. The pharmacological inhibition of JNK retards cancer progression in the mouse. Our study demonstrates that intrahepatic cholestasis leading to hepatocyte exposure to bile acids and inflammation promotes cancer through genomic modifications that can be distinguished from those determined by other aetiological factors.

NCG 4.0: the network of cancer genes in the era of massive mutational screenings of cancer genomes

NCG 4.0 is the latest update of the Network of Cancer Genes, a web-based repository of systems-level properties of cancer genes. In its current version, the database collects information on 537 known (i.e. experimentally supported) and 1463 candidate (i.e. inferred using statistical methods) cancer genes. Candidate cancer genes derive from the manual revision of 67 original publications describing the mutational screening of 3460 human exomes and genomes in 23 different cancer types. For all 2000 cancer genes, duplicability, evolutionary origin, expression, functional annotation, interaction network with other human proteins and with microRNAs are reported. In addition to providing a substantial update of cancer-related information, NCG 4.0 also introduces two new features. The first is the annotation of possible false-positive cancer drivers, defined as candidate cancer genes inferred from large-scale screenings whose association with cancer is likely to be spurious. The second is the description of the systems-level properties of 64 human microRNAs that are causally involved in cancer progression (oncomiRs). Owing to the manual revision of all information, NCG 4.0 constitutes a complete and reliable resource on human coding and non-coding genes whose deregulation drives cancer onset and/or progression. NCG 4.0 can also be downloaded as a free application for Android smart phones. Database URL: http://bio.ieo.eu/ncg/.

Recessive cancer genes engage in negative genetic interactions with their functional paralogs

Cancer genetic heterogeneity offers a wide repertoire of molecular determinants to be screened as therapeutic targets. Here, we identify potential anticancer targets by exploiting negative genetic interactions between genes with driver loss-of-function mutations (recessive cancer genes) and their functionally redundant paralogs. We identify recessive genes with additional copies and experimentally test our predictions on three paralogous pairs. We confirm digenic negative interactions between two cancer genes (SMARCA4 and CDH1) and their corresponding paralogs (SMARCA2 and CDH3). Furthermore, we identify a trigenic negative interaction between the cancer gene DNMT3A, its functional paralog DNMT3B, and a third gene, DNMT1, which encodes the only other human DNA-methylase domain. Although our study does not exclude other causes of synthetic lethality, it suggests that functionally redundant paralogs of cancer genes could be targets in anticancer therapy.

Integrated analysis of recurrent properties of cancer genes to identify novel drivers

The heterogeneity of cancer genomes in terms of acquired mutations complicates the identification of genes whose modification may exert a driver role in tumorigenesis. In this study, we present a novel method that integrates expression profiles, mutation effects, and systemic properties of mutated genes to identify novel cancer drivers. We applied our method to ovarian cancer samples and were able to identify putative drivers in the majority of carcinomas without mutations in known cancer genes, thus suggesting that it can be used as a complementary approach to find rare driver mutations that cannot be detected using frequency-based approaches.

Network of Cancer Genes (NCG 3.0): integration and analysis of genetic and network properties of cancer genes

The identification of a constantly increasing number of genes whose mutations are causally implicated in tumor initiation and progression (cancer genes) requires the development of tools to store and analyze them. The Network of Cancer Genes (NCG 3.0) collects information on 1494 cancer genes that have been found mutated in 16 different cancer types. These genes were collected from the Cancer Gene Census as well as from 18 whole exome and 11 whole-genome screenings of cancer samples. For each cancer gene, NCG 3.0 provides a summary of the gene features and the cross-reference to other databases. In addition, it describes duplicability, evolutionary origin, orthology, network properties, interaction partners, microRNA regulation and functional roles of cancer genes and of all genes that are related to them. This integrated network of information can be used to better characterize cancer genes in the context of the system in which they act. The data can also be used to identify novel candidates that share the same properties of known cancer genes and may therefore play a similar role in cancer. NCG 3.0 is freely available at http://bio.ifom-ieo-campus.it/ncg.

Modification of gene duplicability during the evolution of protein interaction network

Duplications of genes encoding highly connected and essential proteins are selected against in several species but not in human, where duplicated genes encode highly connected proteins. To understand when and how gene duplicability changed in evolution, we compare gene and network properties in four species (Escherichia coli, yeast, fly, and human) that are representative of the increase in evolutionary complexity, defined as progressive growth in the number of genes, cells, and cell types. We find that the origin and conservation of a gene significantly correlates with the properties of the encoded protein in the protein-protein interaction network. All four species preserve a core of singleton and central hubs that originated early in evolution, are highly conserved, and accomplish basic biological functions. Another group of hubs appeared in metazoans and duplicated in vertebrates, mostly through vertebrate-specific whole genome duplication. Such recent and duplicated hubs are frequently targets of microRNAs and show tissue-selective expression, suggesting that these are alternative mechanisms to control their dosage. Our study shows how networks modified during evolution and contributes to explaining the occurrence of somatic genetic diseases, such as cancer, in terms of network perturbations.

Structural analysis of the RZZ complex reveals common ancestry with multisubunit vesicle tethering machinery

The RZZ complex recruits dynein to kinetochores. We investigated structure, topology, and interactions of the RZZ subunits (ROD, ZWILCH, and ZW10) in vitro, in vivo, and in silico. We identify neuroblastoma-amplified gene (NAG), a ZW10 binder, as a ROD homolog. ROD and NAG contain an N-terminal beta propeller followed by an alpha solenoid, which is the architecture of certain nucleoporins and vesicle coat subunits, suggesting a distant evolutionary relationship. ZW10 binding to ROD and NAG is mutually exclusive. The resulting ZW10 complexes (RZZ and NRZ) respectively contain ZWILCH and RINT1 as additional subunits. The X-ray structure of ZWILCH, the first for an RZZ subunit, reveals a novel fold distinct from RINT1's. The evolutionarily conserved NRZ likely acts as a tethering complex for retrograde trafficking of COPI vesicles from the Golgi to the endoplasmic reticulum. The RZZ, limited to metazoans, probably evolved from the NRZ, exploiting the dynein-binding capacity of ZW10 to direct dynein to kinetochores.

The (r)evolution of cancer genetics

The identification of an increasing number of cancer genes is opening up unexpected scenarios in cancer genetics. When analyzed for their systemic properties, these genes show a general fragility towards perturbation. A recent paper published in BMC Biology shows how the founder domains of known cancer genes emerged at two macroevolutionary transitions - the advent of the first cell and the transition to metazoan multicellularity.

See research article http://www.biomedcentral.com/1741-7007/8/66

Ultradeep sequencing of a human ultraconserved region reveals somatic and constitutional genomic instability

Ultradeep sequencing of genomes permits the detection of very low-level genomic instability in non-neoplastic tissues of patients with the most common form of inherited colorectal cancer.

Early detection of cancer-associated genomic instability is crucial, particularly in tumour types in which this instability represents the essential underlying mechanism of tumourigenesis. Currently used methods require the presence of already established neoplastic cells because they only detect clonal mutations. In principle, parallel sequencing of single DNA filaments could reveal the early phases of tumour initiation by detecting low-frequency mutations, provided an adequate depth of coverage and an effective control of the experimental error. We applied ultradeep sequencing to estimate the genomic instability of individuals with hereditary non-polyposis colorectal cancer (HNPCC). To overcome the experimental error, we used an ultraconserved region (UCR) of the human genome as an internal control. By comparing the mutability outside and inside the UCR, we observed a tendency of the ultraconserved element to accumulate significantly fewer mutations than the flanking segments in both neoplastic and nonneoplastic HNPCC samples. No difference between the two regions was detectable in cells from healthy donors, indicating that all three HNPCC samples have mutation rates higher than the healthy genome. This is the first, to our knowledge, direct evidence of an intrinsic genomic instability of individuals with heterozygous mutations in mismatch repair genes, and constitutes the proof of principle for the development of a more sensitive molecular assay of genomic instability.

In hereditary non-polyposis colorectal cancer (HNPCC), a germline mutation in one allele of a gene responsible for repairing DNA damage predisposes the host to cancer, because subsequent somatic inactivation of the one wild-type allele leads to genomic instability that favours tumourigenesis. Nonneoplastic tissues of HNPCC individuals are believed to repair DNA normally, as they are heterozygous and thus are thought to be genomically stable. However, methods used to date are known to be incapable of detecting very low levels of genome instability. Here, we present a more sensitive procedure based on the resequencing of a HNPCC genomic region using next-generation sequencing technology. With this approach, we show that genomic instability is in fact detectable in nonneoplastic tissues of HNPCC patients compared with healthy donors. This constitutional instability may predispose them to acquiring the second somatic mutation event needed for cancer development.

Network of Cancer Genes: a web resource to analyze duplicability, orthology and network properties of cancer genes

The Network of Cancer Genes (NCG) collects and integrates data on 736 human genes that are mutated in various types of cancer. For each gene, NCG provides information on duplicability, orthology, evolutionary appearance and topological properties of the encoded protein in a comprehensive version of the human protein-protein interaction network. NCG also stores information on all primary interactors of cancer proteins, thus providing a complete overview of 5357 proteins that constitute direct and indirect determinants of human cancer. With the constant delivery of results from the mutational screenings of cancer genomes, NCG represents a versatile resource for retrieving detailed information on particular cancer genes, as well as for identifying common properties of precompiled lists of cancer genes. NCG is freely available at: http://bio.ifom-ieo-campus.it/ncg.

FancyGene: dynamic visualization of gene structures and protein domain architectures on genomic loci

SUMMARY

FancyGene is a fast and user-friendly web-based tool for producing images of one or more genes directly on the corresponding genomic locus. Starting from a variety of input formats, FancyGene rebuilds the basic components of a gene (UTRs, intron, exons). Once the initial representation is obtained, the user can superimpose additional features-such as protein domains and/or a variety of biological markers-in specific positions. FancyGene is extremely flexible allowing the user to change the resulting image dynamically, modifying colors and shapes and adding and/or removing objects. The output images are generated either in portable network graphics (PNG) or portable document format (PDF) formats and can be used for scientific presentations as well as for publications. The PDF format preserves editing capabilities, allowing picture modification using any vector graphics editor.

Co-evolutionary Models for Reconstructing Ancestral Genomic Sequences: Computational Issues and Biological Examples

The inference of ancestral genomes is a fundamental problem in molecular evolution. Due to the statistical nature of this problem, the most likely or the most parsimonious ancestral genomes usually include considerable error rates. In general, these errors cannot be abolished by utilizing more exhaustive computational approaches, by using longer genomic sequences, or by analyzing more taxa. In recent studies we showed that co-evolution is an important force that can be used for significantly improving the inference of ancestral genome content.

In this work we formally define a computational problem for the inference of ancestral genome content by co-evolution. We show that this problem is NP-hard and present both a Fixed Parameter Tractable (FPT) algorithm, and heuristic approximation algorithms for solving it. The running time of these algorithms on simulated inputs with hundreds of protein families and hundreds of co-evolutionary relations was fast (up to four minutes) and it achieved an approximation ratio < 1.3.

We use our approach to study the ancestral genome content of the Fungi. To this end, we implement our approach on a dataset of 33,931 protein families and 20,317 co-evolutionary relations. Our algorithm added and removed hundreds of proteins from the ancestral genomes inferred by maximum likelihood (ML) or maximum parsimony (MP) while slightly affecting the likelihood/parsimony score of the results. A biological analysis revealed various pieces of evidence that support the biological plausibility of the new solutions.

Low duplicability and network fragility of cancer genes

We identified genomic and network properties of approximately 600 genes mutated in different cancer types. These genes tend not to duplicate but, unlike most human singletons, they encode central hubs of highly interconnected modules within the protein-protein interaction network (PIN). We find that cancer genes are fragile components of the human gene repertoire, sensitive to dosage modification. Furthermore, other nodes of the human PIN with similar properties are rare and probably enriched in candidate cancer genes.

Non-random retention of protein-coding overlapping genes in Metazoa

Background

Although the overlap of transcriptional units occurs frequently in eukaryotic genomes, its evolutionary and biological significance remains largely unclear. Here we report a comparative analysis of overlaps between genes coding for well-annotated proteins in five metazoan genomes (human, mouse, zebrafish, fruit fly and worm).

Results

For all analyzed species the observed number of overlapping genes is always lower than expected assuming functional neutrality, suggesting that gene overlap is negatively selected. The comparison to the random distribution also shows that retained overlaps do not exhibit random features: antiparallel overlaps are significantly enriched, while overlaps lying on the same strand and those involving coding sequences are highly underrepresented. We confirm that overlap is mostly species-specific and provide evidence that it frequently originates through the acquisition of terminal, non-coding exons. Finally, we show that overlapping genes tend to be significantly co-expressed in a breast cancer cDNA library obtained by 454 deep sequencing, and that different overlap types display different patterns of reciprocal expression.

Conclusion

Our data suggest that overlap between protein-coding genes is selected against in Metazoa. However, when retained it may be used as a species-specific mechanism for the reciprocal regulation of neighboring genes. The tendency of overlaps to involve non-coding regions of the genes leads to the speculation that the advantages achieved by an overlapping arrangement may be optimized by evolving regulatory non-coding transcripts.

Family expansion and gene rearrangements contributed to the functional specialization of PRDM genes in vertebrates

Background

Progressive diversification of paralogs after gene expansion is essential to increase their functional specialization. However, mode and tempo of this divergence remain mostly unclear. Here we report the comparative analysis of PRDM genes, a family of putative transcriptional regulators involved in human tumorigenesis.

Results

Our analysis assessed that the PRDM genes originated in metazoans, expanded in vertebrates and further duplicated in primates. We experimentally showed that fast-evolving paralogs are poorly expressed, and that the most recent duplicates, such as primate-specific PRDM7, acquire tissue-specificity. PRDM7 underwent major structural rearrangements that decreased the number of encoded Zn-Fingers and modified gene splicing. Through internal duplication and activation of a non-canonical splice site (GC-AG), PRDM7 can acquire a novel intron. We also detected an alternative isoform that can retain the intron in the mature transcript and that is predominantly expressed in human melanocytes.

Conclusion

Our findings show that (a) molecular evolution of paralogs correlates with their expression pattern; (b) gene diversification is obtained through massive genomic rearrangements; and (c) splicing modification contributes to the functional specialization of novel genes.

Toward automatic reconstruction of a highly resolved tree of life

We have developed an automatable procedure for reconstructing the tree of life with branch lengths comparable across all three domains. The tree has its basis in a concatenation of 31 orthologs occurring in 191 species with sequenced genomes. It revealed interdomain discrepancies in taxonomic classification. Systematic detection and subsequent exclusion of products of horizontal gene transfer increased phylogenetic resolution, allowing us to confirm accepted relationships and resolve disputed and preliminary classifications. For example, we place the phylum Acidobacteria as a sister group of delta-Proteobacteria, support a Gram-positive origin of Bacteria, and suggest a thermophilic last universal common ancestor.

Very-KIND is a novel nervous system specific guanine nucleotide exchange factor for Ras GTPases

The kinase non-catalytic c-lobe domain (KIND) evolved from the catalytic protein kinase fold into a potential protein interaction module for signalling proteins. Spir family actin organizers and the non-receptor phosphatase type 13 (PTP type 13) encode a KIND domain in the very N-terminal parts of the proteins. Here we report the characterization and cloning of a third member of the KIND protein family, which we have named very-KIND (VKIND) because of its two KIND domains. Like the other members of the protein family, VKIND has a KIND domain at the N-terminus. A second KIND domain is located in the central part of the protein. The C-terminal half encodes a guanine nucleotide exchange factor motif for Ras-like GTPases (RasGEF) and a RasGEF N-terminal module (RasGEFN). There is only one VKIND gene in the mammalian genomes and up to now we have found the gene only in vertebrates. During mouse embryogenesis the VKIND gene was specifically expressed in the developing nervous system. In adult mice Northern hybridizations revealed high expression only in brain. Low expression could be detected in ovary. In situ hybridizations showed a specific expression of VKIND in neuronal cells of the granular and Purkinje cell layers of the cerebellum.

Generation and annotation of the DNA sequences of human chromosomes 2 and 4

Human chromosome 2 is unique to the human lineage in being the product of a head-to-head fusion of two intermediate-sized ancestral chromosomes. Chromosome 4 has received attention primarily related to the search for the Huntington's disease gene, but also for genes associated with Wolf-Hirschhorn syndrome, polycystic kidney disease and a form of muscular dystrophy. Here we present approximately 237 million base pairs of sequence for chromosome 2, and 186 million base pairs for chromosome 4, representing more than 99.6% of their euchromatic sequences. Our initial analyses have identified 1,346 protein-coding genes and 1,239 pseudogenes on chromosome 2, and 796 protein-coding genes and 778 pseudogenes on chromosome 4. Extensive analyses confirm the underlying construction of the sequence, and expand our understanding of the structure and evolution of mammalian chromosomes, including gene deserts, segmental duplications and highly variant regions.

Complex genomic rearrangements lead to novel primate gene function

Orthologous genes that maintain a single-copy status in a broad range of species may indicate a selection against gene duplication. If this is the case, then duplicates of such genes that do survive may have escaped the dosage control by rapid and sizable changes in their function. To test this hypothesis and to develop a strategy for the identification of novel gene functions, we have analyzed 22 primate-specific intrachromosomal duplications of genes with a single-copy ortholog in all other completely sequenced metazoans. When comparing this set to genes not exposed to the single-copy status constraint, we observed a higher tendency of the former to modify their gene structure, often through complex genomic rearrangements. The analysis of the most dramatic of these duplications, affecting approximately 10% of human Chromosome 2, enabled a detailed reconstruction of the events leading to the appearance of a novel gene family. The eight members of this family originated from the highly conserved nucleoporin RanBP2 by several genetic rearrangements such as segmental duplications, inversions, translocations, exon loss, and domain accretion. We have experimentally verified that at least one of the newly formed proteins has a cellular localization different from RanBP2's, and we show that positive selection did act on specific domains during evolution.

The WHy domain mediates the response to desiccation in plants and bacteria

MOTIVATION

The hypersensitive response (HR) is a process activated by plants after microbial infection. Its main phenotypic effects are both a programmed death of the plant cells near the infection site and a reduction of the microbial proliferation. Although many resistance genes (R genes) associated to HR have been identified, very little is known about the molecular mechanisms activated after their expression.

RESULTS

The analysis of the product of one of the R genes, the Hin1 protein, led to the identification of a novel domain, which we named WHy because it is detectable in proteins involved in Water stress and Hypersensitive response. The expression of this domain during both biotic infection and response to desiccation points to a molecular machinery common to these two stress conditions. Moreover, its presence in a restricted number of bacteria suggests a possible use for marking plant pathogenicity.

CONTACT

francesca.ciccarelli@embl.de

SUPPLEMENTARY INFORMATION

Supplementary data (Figures S1 and S2 and Table S1) and the alignment in clustal format are available at http://www.bork.embl.de/~ciccarel/WHy_add_data.html.