VIB-KU Leuven Center for Cancer Biology

In 2008 we published the first set of guidelines for standardizing research in autophagy. Since then, research on this topic has continued to accelerate, and many new scientists have entered the field. Our knowledge base and relevant new technologies have also been expanding. Accordingly, it is important to update these guidelines for monitoring autophagy in different organisms. Various reviews have described the range of assays that have been used for this purpose. Nevertheless, there continues to be confusion regarding acceptable methods to measure autophagy, especially in multicellular eukaryotes. For example, a key point that needs to be emphasized is thatthere is a difference between measurements that monitor the numbers or volume of autophagic elements (e.g., autophagosomes or autolysosomes) at any stage of the autophagic process versus those that measure flux through the autophagy pathway (i.e., the completeprocess including the amount and rate of cargo sequestered and degraded). In particular, a block in macroautophagy that results in autophagosome accumulation must be differentiated from stimuli that increase autophagic activity, defined as increasedautophagy induction coupled with increased delivery to, and degradation within, lysosomes (inmost higher eukaryotes and some protists such as Dictyostelium) or the vacuole (in plants and fungi). In other words, it is especially important that investigators new to the field understand that the appearance of more autophagosomes does not necessarily equate with more autophagy. In fact, in manycases, autophagosomes accumulate because of a block in trafficking to lysosomes without a concomitant change in autophagosome biogenesis, whereas an increase in autolysosomes may reflect a reduction in degradative activity. It is worth emphasizing here that lysosomal digestion is a stage of autophagy and evaluating its competence is a crucial part of the evaluation of autophagic flux, or complete autophagy. Here, we present a set of guidelines for the selection and interpretation of methods for use by investigators who aim to examine macroautophagy and related processes, as well as forreviewers who need to provide realistic and reasonable critiques of papers that are focused on these processes. These guidelines are not meant to be a formulaic set of rules, because the appropriate assays depend in part on the question being asked and the system being used. In addition, we emphasize that no individual assay is guaranteed to be the most appropriate one in every situation, and we strongly recommend the use of multipleassays to monitor autophagy. Along these lines, because of the potential for pleiotropic effects due to blocking autophagy through genetic manipulation, it is imperative to target by gene knockout or RNA interference more than one autophagyrelated protein. In addition, some individual Atg proteins, or groups of proteins, are involved in other cellular pathways implying that not all Atg proteins can be used as a specific marker for an autophagic process. In these guidelines, we consider these various methods of assessing autophagy and what information can, or cannot, be obtained from them. Finally, by discussing the merits and limits of particular assays, we hope to encourage technical innovation in the field.

Microbes are dominant drivers of biogeochemical processes, yet drawing a global picture of functional diversity, microbial community structure, and their ecological determinants remains a grand challenge. We analyzed 7.2 terabases of metagenomic data from 243 Tara Oceans samples from 68 locations in epipelagic and mesopelagic waters across the globe to generate an ocean microbial reference gene catalog with >40 million nonredundant, mostly novel sequences from viruses, prokaryotes, and picoeukaryotes. Using 139 prokaryote-enriched samples, containing >35,000 species, we show vertical stratification with epipelagic community composition mostly driven by temperature rather than other environmental factors or geography. We identify ocean microbial core functionality and reveal that >73% of its abundance is shared with the human gut microbiome despite the physicochemical differences between these two ecosystems.

autophagic responses. Here, we critically discuss current methods of assessing autophagy and the information they can, or cannot, provide. Our ultimate goal is to encourage intellectual and technical innovation in the field.

Fecal microbiome variation in the average, healthy population has remained under-investigated. Here, we analyzed two independent, extensively phenotyped cohorts: the Belgian Flemish Gut Flora Project (FGFP; discovery cohort; N = 1106) and the Dutch LifeLines-DEEP study (LLDeep; replication; N = 1135). Integration with global data sets (N combined = 3948) revealed a 14-genera core microbiota, but the 664 identified genera still underexplore total gut diversity. Sixty-nine clinical and questionnaire-based covariates were found associated to microbiota compositional variation with a 92% replication rate. Stool consistency showed the largest effect size, whereas medication explained largest total variance and interacted with other covariate-microbiota associations. Early-life events such as birth mode were not reflected in adult microbiota composition. Finally, we found that proposed disease marker genera associated to host covariates, urging inclusion of the latter in study design.

Abstract mTORC1 is an important regulator of muscle mass but how it is modulated by oxygen and nutrients is not completely understood. We show that loss of the prolyl hydroxylase domain isoform 1 oxygen sensor in mice (PHD1 KO ) reduces muscle mass. PHD1 KO muscles show impaired mTORC1 activation in response to leucine whereas mTORC1 activation by growth factors or eccentric contractions was preserved. The ability of PHD1 to promote mTORC1 activity is independent of its hydroxylation activity but is caused by decreased protein content of the leucyl tRNA synthetase (LRS) leucine sensor. Mechanistically, PHD1 interacts with and stabilizes LRS. This interaction is promoted during oxygen and amino acid depletion and protects LRS from degradation. Finally, elderly subjects have lower PHD1 levels and LRS activity in muscle from aged versus young human subjects. In conclusion, PHD1 ensures an optimal mTORC1 response to leucine after episodes of metabolic scarcity.

Marine plankton support global biological and geochemical processes. Surveys of their biodiversity have hitherto been geographically restricted and have not accounted for the full range of plankton size. We assessed eukaryotic diversity from 334 size-fractionated photic-zone plankton communities collected across tropical and temperate oceans during the circumglobal Tara Oceans expedition. We analyzed 18S ribosomal DNA sequences across the intermediate plankton-size spectrum from the smallest unicellular eukaryotes (protists, >0.8 micrometers) to small animals of a few millimeters. Eukaryotic ribosomal diversity saturated at ~150,000 operational taxonomic units, about one-third of which could not be assigned to known eukaryotic groups. Diversity emerged at all taxonomic levels, both within the groups comprising the ~11,200 cataloged morphospecies of eukaryotic plankton and among twice as many other deep-branching lineages of unappreciated importance in plankton ecology studies. Most eukaryotic plankton biodiversity belonged to heterotrophic protistan groups, particularly those known to be parasites or symbiotic hosts.

All cancers carry somatic mutations. The patterns of mutation in cancer genomes reflect the DNA damage and repair processes to which cancer cells and their precursors have been exposed. To explore these mechanisms further, we generated catalogs of somatic mutation from 21 breast cancers and applied mathematical methods to extract mutational signatures of the underlying processes. Multiple distinct single- and double-nucleotide substitution signatures were discernible. Cancers with BRCA1 or BRCA2 mutations exhibited a characteristic combination of substitution mutation signatures and a distinctive profile of deletions. Complex relationships between somatic mutation prevalence and transcription were detected. A remarkable phenomenon of localized hypermutation, termed "kataegis," was observed. Regions of kataegis differed between cancers but usually colocalized with somatic rearrangements. Base substitutions in these regions were almost exclusively of cytosine at TpC dinucleotides. The mechanisms underlying most of these mutational signatures are unknown. However, a role for the APOBEC family of cytidine deaminases is proposed.

Deep sequencing of the gut microbiomes of 1135 participants from a Dutch population-based cohort shows relations between the microbiome and 126 exogenous and intrinsic host factors, including 31 intrinsic factors, 12 diseases, 19 drug groups, 4 smoking categories, and 60 dietary factors. These factors collectively explain 18.7% of the variation seen in the interindividual distance of microbial composition. We could associate 110 factors to 125 species and observed that fecal chromogranin A (CgA), a protein secreted by enteroendocrine cells, was exclusively associated with 61 microbial species whose abundance collectively accounted for 53% of microbial composition. Low CgA concentrations were seen in individuals with a more diverse microbiome. These results are an important step toward a better understanding of environment-diet-microbe-host interactions.

Early T-cell precursor acute lymphoblastic leukaemia (ETP ALL) is an aggressive malignancy of unknown genetic basis. We performed whole-genome sequencing of 12 ETP ALL cases and assessed the frequency of the identified somatic mutations in 94 T-cell acute lymphoblastic leukaemia cases. ETP ALL was characterized by activating mutations in genes regulating cytokine receptor and RAS signalling (67% of cases; NRAS, KRAS, FLT3, IL7R, JAK3, JAK1, SH2B3 and BRAF), inactivating lesions disrupting haematopoietic development (58%; GATA3, ETV6, RUNX1, IKZF1 and EP300) and histone-modifying genes (48%; EZH2, EED, SUZ12, SETD2 and EP300). We also identified new targets of recurrent mutation including DNM2, ECT2L and RELN. The mutational spectrum is similar to myeloid tumours, and moreover, the global transcriptional profile of ETP ALL was similar to that of normal and myeloid leukaemia haematopoietic stem cells. These findings suggest that addition of myeloid-directed therapies might improve the poor outcome of ETP ALL. This work shows that treatments used for acute myeloid leukaemia and targeted therapies could be used for early T-cell precursor acute lymphoblastic leukaemia. The early T-cell precursor (ETP) subtype of childhood acute lymphoblastic leukaemia (ALL) has a poor prognosis when treated with standard chemotherapy. Whole genome sequencing is used here to gain insights into the genetic basis of the condition. The results reveal a high frequency of activating mutations in genes regulating cytokine receptor and Ras signalling, lesions that disrupt haemopoiesis (many of which arise from chromosomal rearrangements that generate novel chimeric in-frame fusion genes), and inactivating mutations in histone modifying genes. This mutation pattern resembles that of myeloid malignancies, suggesting that myeloid-directed therapies such as high-dose cytarabine, or targeted therapies that inhibit cytokine receptor and JAK signalling, might be effective in ETP ALL.

Cancer evolves dynamically as clonal expansions supersede one another driven by shifting selective pressures, mutational processes, and disrupted cancer genes. These processes mark the genome, such that a cancer's life history is encrypted in the somatic mutations present. We developed algorithms to decipher this narrative and applied them to 21 breast cancers. Mutational processes evolve across a cancer's lifespan, with many emerging late but contributing extensive genetic variation. Subclonal diversification is prominent, and most mutations are found in just a fraction of tumor cells. Every tumor has a dominant subclonal lineage, representing more than 50% of tumor cells. Minimal expansion of these subclones occurs until many hundreds to thousands of mutations have accumulated, implying the existence of long-lived, quiescent cell lineages capable of substantial proliferation upon acquisition of enabling genomic changes. Expansion of the dominant subclone to an appreciable mass may therefore represent the final rate-limiting step in a breast cancer's development, triggering diagnosis.

We present an allele-specific copy number analysis of the in vivo breast cancer genome. We describe a unique bioinformatics approach, ASCAT (allele-specific copy number analysis of tumors), to accurately dissect the allele-specific copy number of solid tumors, simultaneously estimating and adjusting for both tumor ploidy and nonaberrant cell admixture. This allows calculation of "ASCAT profiles" (genome-wide allele-specific copy-number profiles) from which gains, losses, copy number-neutral events, and loss of heterozygosity (LOH) can accurately be determined. In an early-stage breast carcinoma series, we observe aneuploidy (>2.7n) in 45% of the cases and an average nonaberrant cell admixture of 49%. By aggregation of ASCAT profiles across our series, we obtain genomic frequency distributions of gains and losses, as well as genome-wide views of LOH and copy number-neutral events in breast cancer. In addition, the ASCAT profiles reveal differences in aberrant tumor cell fraction, ploidy, gains, losses, LOH, and copy number-neutral events between the five previously identified molecular breast cancer subtypes. Basal-like breast carcinomas have a significantly higher frequency of LOH compared with other subtypes, and their ASCAT profiles show large-scale loss of genomic material during tumor development, followed by a whole-genome duplication, resulting in near-triploid genomes. Finally, from the ASCAT profiles, we construct a genome-wide map of allelic skewness in breast cancer, indicating loci where one allele is preferentially lost, whereas the other allele is preferentially gained. We hypothesize that these alternative alleles have a different influence on breast carcinoma development.

The microbiota of the human gut is gaining broad attention owing to its association with a wide range of diseases, ranging from metabolic disorders (e.g. obesity and type 2 diabetes) to autoimmune diseases (such as inflammatory bowel disease and type 1 diabetes), cancer and even neurodevelopmental disorders (e.g. autism). Having been increasingly used in biomedical research, mice have become the model of choice for most studies in this emerging field. Mouse models allow perturbations in gut microbiota to be studied in a controlled experimental setup, and thus help in assessing causality of the complex host-microbiota interactions and in developing mechanistic hypotheses. However, pitfalls should be considered when translating gut microbiome research results from mouse models to humans. In this Special Article, we discuss the intrinsic similarities and differences that exist between the two systems, and compare the human and murine core gut microbiota based on a meta-analysis of currently available datasets. Finally, we discuss the external factors that influence the capability of mouse models to recapitulate the gut microbiota shifts associated with human diseases, and investigate which alternative model systems exist for gut microbiota research.

Stratification of women according to their risk of breast cancer based on polygenic risk scores (PRSs) could improve screening and prevention strategies. Our aim was to develop PRSs, optimized for prediction of estrogen receptor (ER)-specific disease, from the largest available genome-wide association dataset and to empirically validate the PRSs in prospective studies. The development dataset comprised 94,075 case subjects and 75,017 control subjects of European ancestry from 69 studies, divided into training and validation sets. Samples were genotyped using genome-wide arrays, and single-nucleotide polymorphisms (SNPs) were selected by stepwise regression or lasso penalized regression. The best performing PRSs were validated in an independent test set comprising 11,428 case subjects and 18,323 control subjects from 10 prospective studies and 190,040 women from UK Biobank (3,215 incident breast cancers). For the best PRSs (313 SNPs), the odds ratio for overall disease per 1 standard deviation in ten prospective studies was 1.61 (95%CI: 1.57-1.65) with area under receiver-operator curve (AUC) = 0.630 (95%CI: 0.628-0.651). The lifetime risk of overall breast cancer in the top centile of the PRSs was 32.6%. Compared with women in the middle quintile, those in the highest 1% of risk had 4.37- and 2.78-fold risks, and those in the lowest 1% of risk had 0.16- and 0.27-fold risks, of developing ER-positive and ER-negative disease, respectively. Goodness-of-fit tests indicated that this PRS was well calibrated and predicts disease risk accurately in the tails of the distribution. This PRS is a powerful and reliable predictor of breast cancer risk that may improve breast cancer prevention programs.

BACKGROUND: Idarucizumab, a monoclonal antibody fragment, was developed to reverse the anticoagulant effect of dabigatran. METHODS: We performed a multicenter, prospective, open-label study to determine whether 5 g of intravenous idarucizumab would be able to reverse the anticoagulant effect of dabigatran in patients who had uncontrolled bleeding (group A) or were about to undergo an urgent procedure (group B). The primary end point was the maximum percentage reversal of the anticoagulant effect of dabigatran within 4 hours after the administration of idarucizumab, on the basis of the diluted thrombin time or ecarin clotting time. Secondary end points included the restoration of hemostasis and safety measures. RESULTS: A total of 503 patients were enrolled: 301 in group A, and 202 in group B. The median maximum percentage reversal of dabigatran was 100% (95% confidence interval, 100 to 100), on the basis of either the diluted thrombin time or the ecarin clotting time. In group A, 137 patients (45.5%) presented with gastrointestinal bleeding and 98 (32.6%) presented with intracranial hemorrhage; among the patients who could be assessed, the median time to the cessation of bleeding was 2.5 hours. In group B, the median time to the initiation of the intended procedure was 1.6 hours; periprocedural hemostasis was assessed as normal in 93.4% of the patients, mildly abnormal in 5.1%, and moderately abnormal in 1.5%. At 90 days, thrombotic events had occurred in 6.3% of the patients in group A and in 7.4% in group B, and the mortality rate was 18.8% and 18.9%, respectively. There were no serious adverse safety signals. CONCLUSIONS: In emergency situations, idarucizumab rapidly, durably, and safely reversed the anticoagulant effect of dabigatran. (Funded by Boehringer Ingelheim; RE-VERSE AD ClinicalTrials.gov number, NCT02104947 .).

Cells succumbing to stress via regulated cell death (RCD) can initiate an adaptive immune response associated with immunological memory, provided they display sufficient antigenicity and adjuvanticity. Moreover, multiple intracellular and microenvironmental features determine the propensity of RCD to drive adaptive immunity. Here, we provide an updated operational definition of immunogenic cell death (ICD), discuss the key factors that dictate the ability of dying cells to drive an adaptive immune response, summarize experimental assays that are currently available for the assessment of ICD in vitro and in vivo, and formulate guidelines for their interpretation.

OBJECTIVE: The assessment of potentially confounding factors affecting colon microbiota composition is essential to the identification of robust microbiome based disease markers. Here, we investigate the link between gut microbiota variation and stool consistency using Bristol Stool Scale classification, which reflects faecal water content and activity, and is considered a proxy for intestinal colon transit time. DESIGN: Through 16S rDNA Illumina profiling of faecal samples of 53 healthy women, we evaluated associations between microbiome richness, Bacteroidetes:Firmicutes ratio, enterotypes, and genus abundance with self-reported, Bristol Stool Scale-based stool consistency. Each sample's microbiota growth potential was calculated to test whether transit time acts as a selective force on gut bacterial growth rates. RESULTS: Stool consistency strongly correlates with all known major microbiome markers. It is negatively correlated with species richness, positively associated to the Bacteroidetes:Firmicutes ratio, and linked to Akkermansia and Methanobrevibacter abundance. Enterotypes are distinctly distributed over the BSS-scores. Based on the correlations between microbiota growth potential and stool consistency scores within both enterotypes, we hypothesise that accelerated transit contributes to colon ecosystem differentiation. While shorter transit times can be linked to increased abundance of fast growing species in Ruminococcaceae-Bacteroides samples, hinting to a washout avoidance strategy of faster replication, this trend is absent in Prevotella-enterotyped individuals. Within this enterotype adherence to host tissue therefore appears to be a more likely bacterial strategy to cope with washout. CONCLUSIONS: The strength of the associations between stool consistency and species richness, enterotypes and community composition emphasises the crucial importance of stool consistency assessment in gut metagenome-wide association studies.

Species interaction networks are shaped by abiotic and biotic factors. Here, as part of the Tara Oceans project, we studied the photic zone interactome using environmental factors and organismal abundance profiles and found that environmental factors are incomplete predictors of community structure. We found associations across plankton functional types and phylogenetic groups to be nonrandomly distributed on the network and driven by both local and global patterns. We identified interactions among grazers, primary producers, viruses, and (mainly parasitic) symbionts and validated network-generated hypotheses using microscopy to confirm symbiotic relationships. We have thus provided a resource to support further research on ocean food webs and integrating biological components into ocean models.

Cells exposed to extreme physicochemical or mechanical stimuli die in an uncontrollable manner, as a result of their immediate structural breakdown. Such an unavoidable variant of cellular demise is generally referred to as 'accidental cell death' (ACD). In most settings, however, cell death is initiated by a genetically encoded apparatus, correlating with the fact that its course can be altered by pharmacologic or genetic interventions. 'Regulated cell death' (RCD) can occur as part of physiologic programs or can be activated once adaptive responses to perturbations of the extracellular or intracellular microenvironment fail. The biochemical phenomena that accompany RCD may be harnessed to classify it into a few subtypes, which often (but not always) exhibit stereotyped morphologic features. Nonetheless, efficiently inhibiting the processes that are commonly thought to cause RCD, such as the activation of executioner caspases in the course of apoptosis, does not exert true cytoprotective effects in the mammalian system, but simply alters the kinetics of cellular demise as it shifts its morphologic and biochemical correlates. Conversely, bona fide cytoprotection can be achieved by inhibiting the transduction of lethal signals in the early phases of the process, when adaptive responses are still operational. Thus, the mechanisms that truly execute RCD may be less understood, less inhibitable and perhaps more homogeneous than previously thought. Here, the Nomenclature Committee on Cell Death formulates a set of recommendations to help scientists and researchers to discriminate between essential and accessory aspects of cell death.

Identifying master regulators of biological processes and mapping their downstream gene networks are key challenges in systems biology. We developed a computational method, called iRegulon, to reverse-engineer the transcriptional regulatory network underlying a co-expressed gene set using cis-regulatory sequence analysis. iRegulon implements a genome-wide ranking-and-recovery approach to detect enriched transcription factor motifs and their optimal sets of direct targets. We increase the accuracy of network inference by using very large motif collections of up to ten thousand position weight matrices collected from various species, and linking these to candidate human TFs via a motif2TF procedure. We validate iRegulon on gene sets derived from ENCODE ChIP-seq data with increasing levels of noise, and we compare iRegulon with existing motif discovery methods. Next, we use iRegulon on more challenging types of gene lists, including microRNA target sets, protein-protein interaction networks, and genetic perturbation data. In particular, we over-activate p53 in breast cancer cells, followed by RNA-seq and ChIP-seq, and could identify an extensive up-regulated network controlled directly by p53. Similarly we map a repressive network with no indication of direct p53 regulation but rather an indirect effect via E2F and NFY. Finally, we generalize our computational framework to include regulatory tracks such as ChIP-seq data and show how motif and track discovery can be combined to map functional regulatory interactions among co-expressed genes. iRegulon is available as a Cytoscape plugin from http://iregulon.aertslab.org.

BACKGROUND: Bevacizumab is approved for the treatment of patients with progressive glioblastoma on the basis of uncontrolled data. Data from a phase 2 trial suggested that the addition of bevacizumab to lomustine might improve overall survival as compared with monotherapies. We sought to determine whether the combination would result in longer overall survival than lomustine alone among patients at first progression of glioblastoma. METHODS: -methylguanine-DNA methyltransferase (MGMT) was assessed. Health-related quality of life and neurocognitive function were evaluated at baseline and every 12 weeks. The primary end point of the trial was overall survival. RESULTS: A total of 437 patients underwent randomization. The median number of 6-week treatment cycles was three in the combination group and one in the monotherapy group. With 329 overall survival events (75.3%), the combination therapy did not provide a survival advantage; the median overall survival was 9.1 months (95% confidence interval [CI], 8.1 to 10.1) in the combination group and 8.6 months (95% CI, 7.6 to 10.4) in the monotherapy group (hazard ratio for death, 0.95; 95% CI, 0.74 to 1.21; P=0.65). Locally assessed progression-free survival was 2.7 months longer in the combination group than in the monotherapy group: 4.2 months versus 1.5 months (hazard ratio for disease progression or death, 0.49; 95% CI, 0.39 to 0.61; P<0.001). Grade 3 to 5 adverse events occurred in 63.6% of the patients in the combination group and 38.1% of the patients in the monotherapy group. The addition of bevacizumab to lomustine affected neither health-related quality of life nor neurocognitive function. The MGMT status was prognostic. CONCLUSIONS: Despite somewhat prolonged progression-free survival, treatment with lomustine plus bevacizumab did not confer a survival advantage over treatment with lomustine alone in patients with progressive glioblastoma. (Funded by an unrestricted educational grant from F. Hoffmann-La Roche and by the EORTC Cancer Research Fund; EORTC 26101 ClinicalTrials.gov number, NCT01290939 ; Eudra-CT number, 2010-023218-30 .).