Which genetic alterations drive tumorigenesis and how they evolve over the course of disease and therapy are central questions in cancer biology. Here we identify 44 recurrently mutated genes and 11 recurrent somatic copy number variations through whole-exome sequencing of 538 chronic lymphocytic leukaemia (CLL) and matched germline DNA samples, 278 of which were collected in a prospective clinical trial. These include previously unrecognized putative cancer drivers (RPS15, IKZF3), and collectively identify RNA processing and export, MYC activity, and MAPK signalling as central pathways involved in CLL. Clonality analysis of this large data set further enabled reconstruction of temporal relationships between driver events. Direct comparison between matched pre-treatment and relapse samples from 59 patients demonstrated highly frequent clonal evolution. Thus, large sequencing data sets of clinically informative samples enable the discovery of novel genes associated with cancer, the network of relationships between the driver events, and their impact on disease relapse and clinical outcome.

Oncogenic activation of BRAF fuels cancer growth by constitutively promoting RAS-independent mitogen-activated protein kinase (MAPK) pathway signalling. Accordingly, RAF inhibitors have brought substantially improved personalized treatment of metastatic melanoma. However, these targeted agents have also revealed an unexpected consequence: stimulated growth of certain cancers. Structurally diverse ATP-competitive RAF inhibitors can either inhibit or paradoxically activate the MAPK pathway, depending whether activation is by BRAF mutation or by an upstream event, such as RAS mutation or receptor tyrosine kinase activation. Here we have identified next-generation RAF inhibitors (dubbed 'paradox breakers') that suppress mutant BRAF cells without activating the MAPK pathway in cells bearing upstream activation. In cells that express the same HRAS mutation prevalent in squamous tumours from patients treated with RAF inhibitors, the first-generation RAF inhibitor vemurafenib stimulated in vitro and in vivo growth and induced expression of MAPK pathway response genes; by contrast the paradox breakers PLX7904 and PLX8394 had no effect. Paradox breakers also overcame several known mechanisms of resistance to first-generation RAF inhibitors. Dissociating MAPK pathway inhibition from paradoxical activation might yield both improved safety and more durable efficacy than first-generation RAF inhibitors, a concept currently undergoing human clinical evaluation with PLX8394.

Neuroblastoma is a malignant paediatric tumour of the sympathetic nervous system. Roughly half of these tumours regress spontaneously or are cured by limited therapy. By contrast, high-risk neuroblastomas have an unfavourable clinical course despite intensive multimodal treatment, and their molecular basis has remained largely elusive. Here we have performed whole-genome sequencing of 56 neuroblastomas (high-risk, n = 39; low-risk, n = 17) and discovered recurrent genomic rearrangements affecting a chromosomal region at 5p15.33 proximal of the telomerase reverse transcriptase gene (TERT). These rearrangements occurred only in high-risk neuroblastomas (12/39, 31%) in a mutually exclusive fashion with MYCN amplifications and ATRX mutations, which are known genetic events in this tumour type. In an extended case series (n = 217), TERT rearrangements defined a subgroup of high-risk tumours with particularly poor outcome. Despite a large structural diversity of these rearrangements, they all induced massive transcriptional upregulation of TERT. In the remaining high-risk tumours, TERT expression was also elevated in MYCN-amplified tumours, whereas alternative lengthening of telomeres was present in neuroblastomas without TERT or MYCN alterations, suggesting that telomere lengthening represents a central mechanism defining this subtype. The 5p15.33 rearrangements juxtapose the TERT coding sequence to strong enhancer elements, resulting in massive chromatin remodelling and DNA methylation of the affected region. Supporting a functional role of TERT, neuroblastoma cell lines bearing rearrangements or amplified MYCN exhibited both upregulated TERT expression and enzymatic telomerase activity. In summary, our findings show that remodelling of the genomic context abrogates transcriptional silencing of TERT in high-risk neuroblastoma and places telomerase activation in the centre of transformation in a large fraction of these tumours.

Sea-level rise can threaten the long-term sustainability of coastal communities and valuable ecosystems such as coral reefs, salt marshes and mangroves. Mangrove forests have the capacity to keep pace with sea-level rise and to avoid inundation through vertical accretion of sediments, which allows them to maintain wetland soil elevations suitable for plant growth. The Indo-Pacific region holds most of the world's mangrove forests, but sediment delivery in this region is declining, owing to anthropogenic activities such as damming of rivers. This decline is of particular concern because the Indo-Pacific region is expected to have variable, but high, rates of future sea-level rise. Here we analyse recent trends in mangrove surface elevation changes across the Indo-Pacific region using data from a network of surface elevation table instruments. We find that sediment availability can enable mangrove forests to maintain rates of soil-surface elevation gain that match or exceed that of sea-level rise, but for 69 per cent of our study sites the current rate of sea-level rise exceeded the soil surface elevation gain. We also present a model based on our field data, which suggests that mangrove forests at sites with low tidal range and low sediment supply could be submerged as early as 2070.

The hominin record from southern Asia for the early Late Pleistocene epoch is scarce. Well-dated and well-preserved fossils older than ∼45,000 years that can be unequivocally attributed to Homo sapiens are lacking. Here we present evidence from the newly excavated Fuyan Cave in Daoxian (southern China). This site has provided 47 human teeth dated to more than 80,000 years old, and with an inferred maximum age of 120,000 years. The morphological and metric assessment of this sample supports its unequivocal assignment to H. sapiens. The Daoxian sample is more derived than any other anatomically modern humans, resembling middle-to-late Late Pleistocene specimens and even contemporary humans. Our study shows that fully modern morphologies were present in southern China 30,000-70,000 years earlier than in the Levant and Europe. Our data fill a chronological and geographical gap that is relevant for understanding when H. sapiens first appeared in southern Asia. The Daoxian teeth also support the hypothesis that during the same period, southern China was inhabited by more derived populations than central and northern China. This evidence is important for the study of dispersal routes of modern humans. Finally, our results are relevant to exploring the reasons for the relatively late entry of H. sapiens into Europe. Some studies have investigated how the competition with H. sapiens may have caused Neanderthals' extinction (see ref. 8 and references therein). Notably, although fully modern humans were already present in southern China at least as early as ∼80,000 years ago, there is no evidence that they entered Europe before ∼45,000 years ago. This could indicate that H. neanderthalensis was indeed an additional ecological barrier for modern humans, who could only enter Europe when the demise of Neanderthals had already started.

It remains unclear whether biodiversity buffers ecosystems against climate extremes, which are becoming increasingly frequent worldwide. Early results suggested that the ecosystem productivity of diverse grassland plant communities was more resistant, changing less during drought, and more resilient, recovering more quickly after drought, than that of depauperate communities. However, subsequent experimental tests produced mixed results. Here we use data from 46 experiments that manipulated grassland plant diversity to test whether biodiversity provides resistance during and resilience after climate events. We show that biodiversity increased ecosystem resistance for a broad range of climate events, including wet or dry, moderate or extreme, and brief or prolonged events. Across all studies and climate events, the productivity of low-diversity communities with one or two species changed by approximately 50% during climate events, whereas that of high-diversity communities with 16-32 species was more resistant, changing by only approximately 25%. By a year after each climate event, ecosystem productivity had often fully recovered, or overshot, normal levels of productivity in both high- and low-diversity communities, leading to no detectable dependence of ecosystem resilience on biodiversity. Our results suggest that biodiversity mainly stabilizes ecosystem productivity, and productivity-dependent ecosystem services, by increasing resistance to climate events. Anthropogenic environmental changes that drive biodiversity loss thus seem likely to decrease ecosystem stability, and restoration of biodiversity to increase it, mainly by changing the resistance of ecosystem productivity to climate events.

Solid cancer cells commonly enter the blood and disseminate systemically, but are highly inefficient at forming distant metastases for poorly understood reasons. Here we studied human melanomas that differed in their metastasis histories in patients and in their capacity to metastasize in NOD-SCID-Il2rg(-/-) (NSG) mice. We show that melanomas had high frequencies of cells that formed subcutaneous tumours, but much lower percentages of cells that formed tumours after intravenous or intrasplenic transplantation, particularly among inefficiently metastasizing melanomas. Melanoma cells in the blood and visceral organs experienced oxidative stress not observed in established subcutaneous tumours. Successfully metastasizing melanomas underwent reversible metabolic changes during metastasis that increased their capacity to withstand oxidative stress, including increased dependence on NADPH-generating enzymes in the folate pathway. Antioxidants promoted distant metastasis in NSG mice. Folate pathway inhibition using low-dose methotrexate, ALDH1L2 knockdown, or MTHFD1 knockdown inhibited distant metastasis without significantly affecting the growth of subcutaneous tumours in the same mice. Oxidative stress thus limits distant metastasis by melanoma cells in vivo.

Studies of individual living cells have revealed that many transcription factors activate in dynamic, and often stochastic, pulses within the same cell. However, it has remained unclear whether cells might exploit the dynamic interaction of these pulses to control gene expression. Here, using quantitative single-cell time-lapse imaging of Saccharomyces cerevisiae, we show that the pulsatile transcription factors Msn2 and Mig1 combinatorially regulate their target genes through modulation of their relative pulse timing. The activator Msn2 and repressor Mig1 showed pulsed activation in either a temporally overlapping or non-overlapping manner during their transient response to different inputs, with only the non-overlapping dynamics efficiently activating target gene expression. Similarly, under constant environmental conditions, where Msn2 and Mig1 exhibit sporadic pulsing, glucose concentration modulated the temporal overlap between pulses of the two factors. Together, these results reveal a time-based mode of combinatorial gene regulation. Regulation through relative signal timing is common in engineering and neurobiology, and these results suggest that it could also function broadly within the signalling and regulatory systems of the cell.

Acetyl-CoA carboxylase (ACC) has crucial roles in fatty acid metabolism and is an attractive target for drug discovery against diabetes, cancer and other diseases. Saccharomyces cerevisiae ACC (ScACC) is crucial for the production of very-long-chain fatty acids and the maintenance of the nuclear envelope. ACC contains biotin carboxylase (BC) and carboxyltransferase (CT) activities, and its biotin is linked covalently to the biotin carboxyl carrier protein (BCCP). Most eukaryotic ACCs are 250-kilodalton (kDa), multi-domain enzymes and function as homodimers and higher oligomers. They contain a unique, 80-kDa central region that shares no homology with other proteins. Although the structures of the BC, CT and BCCP domains and other biotin-dependent carboxylase holoenzymes are known, there is currently no structural information on the ACC holoenzyme. Here we report the crystal structure of the full-length, 500-kDa holoenzyme dimer of ScACC. The structure is remarkably different from that of the other biotin-dependent carboxylases. The central region contains five domains and is important for positioning the BC and CT domains for catalysis. The structure unexpectedly reveals a dimer of the BC domain and extensive conformational differences compared to the structure of the BC domain alone, which is a monomer. These structural changes reveal why the BC domain alone is catalytically inactive and define the molecular mechanism for the inhibition of eukaryotic ACC by the natural product soraphen A and by phosphorylation of a Ser residue just before the BC domain core in mammalian ACC. The BC and CT active sites are separated by 80 Å, and the entire BCCP domain must translocate during catalysis.

The most abundant mRNA post-transcriptional modification is N(6)-methyladenosine (m(6)A), which has broad roles in RNA biology. In mammalian cells, the asymmetric distribution of m(6)A along mRNAs results in relatively less methylation in the 5' untranslated region (5'UTR) compared to other regions. However, whether and how 5'UTR methylation is regulated is poorly understood. Despite the crucial role of the 5'UTR in translation initiation, very little is known about whether m(6)A modification influences mRNA translation. Here we show that in response to heat shock stress, certain adenosines within the 5'UTR of newly transcribed mRNAs are preferentially methylated. We find that the dynamic 5'UTR methylation is a result of stress-induced nuclear localization of YTHDF2, a well-characterized m(6)A 'reader'. Upon heat shock stress, the nuclear YTHDF2 preserves 5'UTR methylation of stress-induced transcripts by limiting the m(6)A 'eraser' FTO from demethylation. Remarkably, the increased 5'UTR methylation in the form of m(6)A promotes cap-independent translation initiation, providing a mechanism for selective mRNA translation under heat shock stress. Using Hsp70 mRNA as an example, we demonstrate that a single m(6)A modification site in the 5'UTR enables translation initiation independent of the 5' end N(7)-methylguanosine cap. The elucidation of the dynamic features of 5'UTR methylation and its critical role in cap-independent translation not only expands the breadth of physiological roles of m(6)A, but also uncovers a previously unappreciated translational control mechanism in heat shock response.