epiAneufinder identifies copy number alterations from single-cell ATAC-seq data.

Single-cell open chromatin profiling via scATAC-seq has become a mainstream measurement of open chromatin in single-cells. Here we present epiAneufinder, an algorithm that exploits the read count information from scATAC-seq data to extract genome-wide copy number alterations (CNAs) for individual cells, allowing the study of CNA heterogeneity present in a sample at the single-cell level. Using different cancer scATAC-seq datasets, we show that epiAneufinder can identify intratumor clonal heterogeneity in populations of single cells based on their CNA profiles. We demonstrate that these profiles are concordant with the ones inferred from single-cell whole genome sequencing data for the same samples. EpiAneufinder allows the inference of single-cell CNA information from scATAC-seq data, without the need of additional experiments, unlocking a layer of genomic variation which is otherwise unexplored.

Dietary palmitic acid promotes a prometastatic memory via Schwann cells.

Fatty acid uptake and altered metabolism constitute hallmarks of metastasis1,2, yet evidence of the underlying biology, as well as whether all dietary fatty acids are prometastatic, is lacking. Here we show that dietary palmitic acid (PA), but not oleic acid or linoleic acid, promotes metastasis in oral carcinomas and melanoma in mice. Tumours from mice that were fed a short-term palm-oil-rich diet (PA), or tumour cells that were briefly exposed to PA in vitro, remained highly metastatic even after being serially transplanted (without further exposure to high levels of PA). This PA-induced prometastatic memory requires the fatty acid transporter CD36 and is associated with the stable deposition of histone H3 lysine 4 trimethylation by the methyltransferase Set1A (as part of the COMPASS complex (Set1A/COMPASS)). Bulk, single-cell and positional RNA-sequencing analyses indicate that genes with this prometastatic memory predominantly relate to a neural signature that stimulates intratumoural Schwann cells and innervation, two parameters that are strongly correlated with metastasis but are aetiologically poorly understood3,4. Mechanistically, tumour-associated Schwann cells secrete a specialized proregenerative extracellular matrix, the ablation of which inhibits metastasis initiation. Both the PA-induced memory of this proneural signature and its long-term boost in metastasis require the transcription factor EGR2 and the glial-cell-stimulating peptide galanin. In summary, we provide evidence that a dietary metabolite induces stable transcriptional and chromatin changes that lead to a long-term stimulation of metastasis, and that this is related to a proregenerative state of tumour-activated Schwann cells.

AlphaBeta: computational inference of epimutation rates and spectra from high-throughput DNA methylation data in plants.

Stochastic changes in DNA methylation (i.e., spontaneous epimutations) contribute to methylome diversity in plants. Here, we describe AlphaBeta, a computational method for estimating the precise rate of such stochastic events using pedigree-based DNA methylation data as input. We demonstrate how AlphaBeta can be employed to study transgenerationally heritable epimutations in clonal or sexually derived mutation accumulation lines, as well as somatic epimutations in long-lived perennials. Application of our method to published and new data reveals that spontaneous epimutations accumulate neutrally at the genome-wide scale, originate mainly during somatic development and that they can be used as a molecular clock for age-dating trees.

BMAL1-Driven Tissue Clocks Respond Independently to Light to Maintain Homeostasis.

Circadian rhythms control organismal physiology throughout the day. At the cellular level, clock regulation is established by a self-sustained Bmal1-dependent transcriptional oscillator network. However, it is still unclear how different tissues achieve a synchronized rhythmic physiology. That is, do they respond independently to environmental signals, or require interactions with each other to do so? We show that unexpectedly, light synchronizes the Bmal1-dependent circadian machinery in single tissues in the absence of Bmal1 in all other tissues. Strikingly, light-driven tissue autonomous clocks occur without rhythmic feeding behavior and are lost in constant darkness. Importantly, tissue-autonomous Bmal1 partially sustains homeostasis in otherwise arrhythmic and prematurely aging animals. Our results therefore support a two-branched model for the daily synchronization of tissues: an autonomous response branch, whereby light entrains circadian clocks without any commitment of other Bmal1-dependent clocks, and a memory branch using other Bmal1-dependent clocks to "remember" time in the absence of external cues.

Loss of G9a preserves mutation patterns but increases chromatin accessibility, genomic instability and aggressiveness in skin tumours.

Mutations in, and the altered expression of, epigenetic modifiers are pervasive in human tumours, making epigenetic factors attractive antitumour targets. The open-versus-closed chromatin state within the cells-of-origin of cancer correlates with the uneven distribution of mutations. However, the long-term effect of targeting epigenetic modifiers on mutability in patients with cancer is unclear. Here, we increased chromatin accessibility by deleting the histone H3 lysine 9 (H3K9) methyltransferase G9a in murine epidermis and show that this does not alter the single nucleotide variant burden or global genomic distribution in chemical mutagen-induced squamous tumours. G9a-depleted tumours develop after a prolonged latency compared with their wild-type counterparts, but are more aggressive and have an expanded cancer progenitor pool, pronounced genomic instability and frequent loss-of-function p53 mutations. Thus, we call for caution when assessing long-term therapeutic benefits of chromatin modifier inhibitors, which may promote more aggressive disease.

The Chara Genome: Secondary Complexity and Implications for Plant Terrestrialization.

Land plants evolved from charophytic algae, among which Charophyceae possess the most complex body plans. We present the genome of Chara braunii; comparison of the genome to those of land plants identified evolutionary novelties for plant terrestrialization and land plant heritage genes. C. braunii employs unique xylan synthases for cell wall biosynthesis, a phragmoplast (cell separation) mechanism similar to that of land plants, and many phytohormones. C. braunii plastids are controlled via land-plant-like retrograde signaling, and transcriptional regulation is more elaborate than in other algae. The morphological complexity of this organism may result from expanded gene families, with three cases of particular note: genes effecting tolerance to reactive oxygen species (ROS), LysM receptor-like kinases, and transcription factors (TFs). Transcriptomic analysis of sexual reproductive structures reveals intricate control by TFs, activity of the ROS gene network, and the ancestral use of plant-like storage and stress protection proteins in the zygote.

Publisher Correction: MSK1 regulates luminal cell differentiation and metastatic dormancy in ER+ breast cancer.

In the version of this Article originally published, the boxes framing the two plots in Fig. 1g were misaligned from the axes due to a technical error. This has now been corrected in all versions of the Article.

MSK1 regulates luminal cell differentiation and metastatic dormancy in ER+ breast cancer.

For many patients with breast cancer, symptomatic bone metastases appear after years of latency. How micrometastatic lesions remain dormant and undetectable before initiating colonization is unclear. Here, we describe a mechanism involved in bone metastatic latency of oestrogen receptor-positive (ER+) breast cancer. Using an in vivo genome-wide short hairpin RNA screening, we identified the kinase MSK1 as an important regulator of metastatic dormancy in breast cancer. In patients with ER+ breast cancer, low MSK1 expression associates with early metastasis. We show that MSK1 downregulation impairs the differentiation of breast cancer cells, increasing their bone homing and growth capacities. MSK1 controls the expression of genes required for luminal cell differentiation, including the GATA3 and FOXA1 transcription factors, by modulating their promoter chromatin status. Our results indicate that MSK1 prevents metastatic progression of ER+ breast cancer, suggesting that stratifying patients with breast cancer as high or low risk for early relapse based on MSK1 expression could improve prognosis.

The Physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution.

The draft genome of the moss model, Physcomitrella patens, comprised approximately 2000 unordered scaffolds. In order to enable analyses of genome structure and evolution we generated a chromosome-scale genome assembly using genetic linkage as well as (end) sequencing of long DNA fragments. We find that 57% of the genome comprises transposable elements (TEs), some of which may be actively transposing during the life cycle. Unlike in flowering plant genomes, gene- and TE-rich regions show an overall even distribution along the chromosomes. However, the chromosomes are mono-centric with peaks of a class of Copia elements potentially coinciding with centromeres. Gene body methylation is evident in 5.7% of the protein-coding genes, typically coinciding with low GC and low expression. Some giant virus insertions are transcriptionally active and might protect gametes from viral infection via siRNA mediated silencing. Structure-based detection methods show that the genome evolved via two rounds of whole genome duplications (WGDs), apparently common in mosses but not in liverworts and hornworts. Several hundred genes are present in colinear regions conserved since the last common ancestor of plants. These syntenic regions are enriched for functions related to plant-specific cell growth and tissue organization. The P. patens genome lacks the TE-rich pericentromeric and gene-rich distal regions typical for most flowering plant genomes. More non-seed plant genomes are needed to unravel how plant genomes evolve, and to understand whether the P. patens genome structure is typical for mosses or bryophytes.

Combination of the Endogenous lhcsr1 Promoter and Codon Usage Optimization Boosts Protein Expression in the Moss Physcomitrella patens.

The moss Physcomitrella patens is used both as an evo-devo model and biotechnological production system for metabolites and pharmaceuticals. Strong in vivo expression of genes of interest is important for production of recombinant proteins, e.g., selectable markers, fluorescent proteins, or enzymes. In this regard, the choice of the promoter sequence as well as codon usage optimization are two important inside factors to consider in order to obtain optimum protein accumulation level. To reliably quantify fluorescence, we transfected protoplasts with promoter:GFP fusion constructs and measured fluorescence intensity of living protoplasts in a plate reader system. We used the red fluorescent protein mCherry under 2x 35S promoter control as second reporter to normalize for different transfection efficiencies. We derived a novel endogenous promoter and compared deletion variants with exogenous promoters. We used different codon-adapted green fluorescent protein (GFP) genes to evaluate the influence of promoter choice and codon optimization on protein accumulation in P. patens, and show that the promoter of the gene of P. patens chlorophyll a/b binding protein lhcsr1 drives expression of GFP in protoplasts significantly (more than twofold) better than the commonly used 2x 35S promoter or the rice actin1 promoter. We identified a shortened 677 bp version of the lhcsr1 promoter that retains full activity in protoplasts. The codon optimized GFP yields significantly (more than twofold) stronger fluorescence signals and thus demonstrates that adjusting codon usage in P. patens can increase expression strength. In combination, new promotor and codon optimized GFP conveyed sixfold increased fluorescence signal.

Circadian Reprogramming in the Liver Identifies Metabolic Pathways of Aging.

The process of aging and circadian rhythms are intimately intertwined, but how peripheral clocks involved in metabolic homeostasis contribute to aging remains unknown. Importantly, caloric restriction (CR) extends lifespan in several organisms and rewires circadian metabolism. Using young versus old mice, fed ad libitum or under CR, we reveal reprogramming of the circadian transcriptome in the liver. These age-dependent changes occur in a highly tissue-specific manner, as demonstrated by comparing circadian gene expression in the liver versus epidermal and skeletal muscle stem cells. Moreover, de novo oscillating genes under CR show an enrichment in SIRT1 targets in the liver. This is accompanied by distinct circadian hepatic signatures in NAD+-related metabolites and cyclic global protein acetylation. Strikingly, this oscillation in acetylation is absent in old mice while CR robustly rescues global protein acetylation. Our findings indicate that the clock operates at the crossroad between protein acetylation, liver metabolism, and aging.

Aged Stem Cells Reprogram Their Daily Rhythmic Functions to Adapt to Stress.

Normal homeostatic functions of adult stem cells have rhythmic daily oscillations that are believed to become arrhythmic during aging. Unexpectedly, we find that aged mice remain behaviorally circadian and that their epidermal and muscle stem cells retain a robustly rhythmic core circadian machinery. However, the oscillating transcriptome is extensively reprogrammed in aged stem cells, switching from genes involved in homeostasis to those involved in tissue-specific stresses, such as DNA damage or inefficient autophagy. Importantly, deletion of circadian clock components did not reproduce the hallmarks of this reprogramming, underscoring that rewiring, rather than arrhythmia, is associated with physiological aging. While age-associated rewiring of the oscillatory diurnal transcriptome is not recapitulated by a high-fat diet in young adult mice, it is significantly prevented by long-term caloric restriction in aged mice. Thus, stem cells rewire their diurnal timed functions to adapt to metabolic cues and to tissue-specific age-related traits.

CREBBP is a target of epigenetic, but not genetic, modification in juvenile myelomonocytic leukemia.

BACKGROUND: Juvenile myelomonocytic leukemia (JMML) is a myeloproliferative neoplasm of childhood whose clinical heterogeneity is only poorly represented by gene sequence alterations. It was previously shown that aberrant DNA methylation of distinct target genes defines a more aggressive variant of JMML, but only few significant targets are known so far. To get a broader picture of disturbed CpG methylation patterns in JMML, we carried out a methylation screen of 34 candidate genes in 45 patients using quantitative mass spectrometry. FINDINGS: Five of 34 candidate genes analyzed showed recurrent hypermethylation in JMML. cAMP-responsive element-binding protein-binding protein (CREBBP) was the most frequent target of epigenetic modification (77 % of cases). However, no pathogenic mutations of CREBBP were identified in a genetic analysis of 64 patients. CREBBP hypermethylation correlated with clinical parameters known to predict poor outcome. CONCLUSIONS: This study supports the relevance of epigenetic aberrations in JMML pathophysiology. Our data confirm that DNA hypermethylation in JMML is highly target-specific and associated with higher-risk features. These findings encourage the development of prognostic markers based on epigenetic alterations, which will be helpful in the difficult clinical management of this heterogeneous disease.

De novo assembly and comparative analysis of the Ceratodon purpureus transcriptome.

The bryophytes are a morphologically and ecologically diverse group of plants that have recently emerged as major model systems for a variety of biological processes. In particular, the genome sequence of the moss, Physcomitrella patens, has significantly enhanced our understanding of the evolution of developmental processes in land plants. However, to fully explore the diversity within bryophytes, we need additional genomic resources. Here, we describe analyses of the transcriptomes of a male and a female isolate of the moss, Ceratodon purpureus, generated using the 454 FLX technology. Comparative analyses between C. purpureus and P. patens indicated that this strategy generated nearly complete coverage of the protonemal transcriptome. An analysis of the overlap in gene sets between C. purpureus and P. patens provides new insights into the evolution of gene family composition across the land plants. In spite of the overall transcriptomic similarity between the two species, Ka /Ks analysis of P. patens and C. purpureus suggests considerable physiological and developmental divergence. Additionally, while the codon usage was very similar between these two mosses, C. purpureus genes showed a slightly greater codon usage bias than P. patens genes potentially because of the contrasting mating system of the two species. Finally, we found evidence of a genome doubling ~65-76 MYA that likely coincided with the contemporaneous polyploidy event inferred for P. patens but postdates the divergence of P. patens and C. purpureus. The powerful laboratory tools now available for C. purpureus will enable the research community to fully exploit these genomic resources.

The chromatin landscape of the moss Physcomitrella patens and its dynamics during development and drought stress.

The moss Physcomitrella patens is an important model organism for evo-devo studies. Here, we determined the genome-wide chromatin landscape of five important histone three (H3) modifications (H3K4me3, H3K27me3, H3K27Ac, H3K9Ac and H3K9me2) and describe the changes to these histone marks in two contrasted situations, developmental transition and abiotic (drought) stress. Integrative analysis of these histone H3 modifications revealed their preferential association into 15 chromatin states (CS) in genic regions of the P. patens genome. Synergistic relationships that influence expression levels were revealed for the three activating marks H3K4me3, H3K27Ac and H3K9Ac, while an antagonistic relationship was found between CS containing the H3K27me3 and H3K27Ac marks, suggesting that H3K27 is a key indexing residue regarding transcriptional output. Concerning the alteration of histone marks in response to developmental transition (juvenile to adult) and drought stress, the three activating marks H3K4me3, H3K27Ac and H3K9Ac show significant changes in both situations. However, changes to H3K27me3 are central only for genes differentially expressed during development. Interestingly, genes induced during drought stress show significant histone mark toggling during developmental transition. This situation suggests that drought induced adult (gametophore expressed) genes are primed to respond to this stress during the juvenile to adult transition.

Genome of an arbuscular mycorrhizal fungus provides insight into the oldest plant symbiosis.

The mutualistic symbiosis involving Glomeromycota, a distinctive phylum of early diverging Fungi, is widely hypothesized to have promoted the evolution of land plants during the middle Paleozoic. These arbuscular mycorrhizal fungi (AMF) perform vital functions in the phosphorus cycle that are fundamental to sustainable crop plant productivity. The unusual biological features of AMF have long fascinated evolutionary biologists. The coenocytic hyphae host a community of hundreds of nuclei and reproduce clonally through large multinucleated spores. It has been suggested that the AMF maintain a stable assemblage of several different genomes during the life cycle, but this genomic organization has been questioned. Here we introduce the 153-Mb haploid genome of Rhizophagus irregularis and its repertoire of 28,232 genes. The observed low level of genome polymorphism (0.43 SNP per kb) is not consistent with the occurrence of multiple, highly diverged genomes. The expansion of mating-related genes suggests the existence of cryptic sex-related processes. A comparison of gene categories confirms that R. irregularis is close to the Mucoromycotina. The AMF obligate biotrophy is not explained by genome erosion or any related loss of metabolic complexity in central metabolism, but is marked by a lack of genes encoding plant cell wall-degrading enzymes and of genes involved in toxin and thiamine synthesis. A battery of mycorrhiza-induced secreted proteins is expressed in symbiotic tissues. The present comprehensive repertoire of R. irregularis genes provides a basis for future research on symbiosis-related mechanisms in Glomeromycota.

Algal genomes reveal evolutionary mosaicism and the fate of nucleomorphs.

Cryptophyte and chlorarachniophyte algae are transitional forms in the widespread secondary endosymbiotic acquisition of photosynthesis by engulfment of eukaryotic algae. Unlike most secondary plastid-bearing algae, miniaturized versions of the endosymbiont nuclei (nucleomorphs) persist in cryptophytes and chlorarachniophytes. To determine why, and to address other fundamental questions about eukaryote-eukaryote endosymbiosis, we sequenced the nuclear genomes of the cryptophyte Guillardia theta and the chlorarachniophyte Bigelowiella natans. Both genomes have >21,000 protein genes and are intron rich, and B. natans exhibits unprecedented alternative splicing for a single-celled organism. Phylogenomic analyses and subcellular targeting predictions reveal extensive genetic and biochemical mosaicism, with both host- and endosymbiont-derived genes servicing the mitochondrion, the host cell cytosol, the plastid and the remnant endosymbiont cytosol of both algae. Mitochondrion-to-nucleus gene transfer still occurs in both organisms but plastid-to-nucleus and nucleomorph-to-nucleus transfers do not, which explains why a small residue of essential genes remains locked in each nucleomorph.

Cyanophora paradoxa genome elucidates origin of photosynthesis in algae and plants.

The primary endosymbiotic origin of the plastid in eukaryotes more than 1 billion years ago led to the evolution of algae and plants. We analyzed draft genome and transcriptome data from the basally diverging alga Cyanophora paradoxa and provide evidence for a single origin of the primary plastid in the eukaryote supergroup Plantae. C. paradoxa retains ancestral features of starch biosynthesis, fermentation, and plastid protein translocation common to plants and algae but lacks typical eukaryotic light-harvesting complex proteins. Traces of an ancient link to parasites such as Chlamydiae were found in the genomes of C. paradoxa and other Plantae. Apparently, Chlamydia-like bacteria donated genes that allow export of photosynthate from the plastid and its polymerization into storage polysaccharide in the cytosol.