The genome of Chenopodium quinoa.

Chenopodium quinoa (quinoa) is a highly nutritious grain identified as an important crop to improve world food security. Unfortunately, few resources are available to facilitate its genetic improvement. Here we report the assembly of a high-quality, chromosome-scale reference genome sequence for quinoa, which was produced using single-molecule real-time sequencing in combination with optical, chromosome-contact and genetic maps. We also report the sequencing of two diploids from the ancestral gene pools of quinoa, which enables the identification of sub-genomes in quinoa, and reduced-coverage genome sequences for 22 other samples of the allotetraploid goosefoot complex. The genome sequence facilitated the identification of the transcription factor likely to control the production of anti-nutritional triterpenoid saponins found in quinoa seeds, including a mutation that appears to cause alternative splicing and a premature stop codon in sweet quinoa strains. These genomic resources are an important first step towards the genetic improvement of quinoa.

Corrigendum: The genome of Chenopodium quinoa.

This corrects the article DOI: 10.1038/nature21370.

The KDM4A/KDM4C/NF-κB and WDR5 epigenetic cascade regulates the activation of B cells.

T follicular helper (Tfh) cell-derived signals promote activation and proliferation of antigen-primed B cells. It remains unclear whether epigenetic regulation is involved in the B cell responses to Tfh cell-derived signals. Here, we demonstrate that Tfh cell-mimicking signals induce the expression of histone demethylases KDM4A and KDM4C, and the concomitant global down-regulation of their substrates, H3K9me3/me2, in B cells. Depletion of KDM4A and KDM4C potentiates B cell activation and proliferation in response to Tfh cell-derived signals. ChIP-seq and de novo motif analysis reveals NF-κB p65 as a binding partner of KDM4A and KDM4C. Their co-targeting to Wdr5, a MLL complex member promoting H3K4 methylation, up-regulates cell cycle inhibitors Cdkn2c and Cdkn3. Thus, Tfh cell-derived signals trigger KDM4A/KDM4C - WDR5 - Cdkn2c/Cdkn3 cascade in vitro, an epigenetic mechanism regulating proper proliferation of activated B cells. This pathway is dysregulated in B cells from systemic lupus erythematosus patients and may represent a pathological link.

Thermopriming triggers splicing memory in Arabidopsis.

Abiotic and biotic stresses limit crop productivity. Exposure to a non-lethal stress, referred to as priming, can allow plants to survive subsequent and otherwise lethal conditions; the priming effect persists even after a prolonged stress-free period. However, the molecular mechanisms underlying priming are not fully understood. Here, we investigated the molecular basis of heat-shock memory and the role of priming in Arabidopsis thaliana. Comprehensive analysis of transcriptome-wide changes in gene expression and alternative splicing in primed and non-primed plants revealed that alternative splicing functions as a novel component of heat-shock memory. We show that priming of plants with a non-lethal heat stress results in de-repression of splicing after a second exposure to heat stress. By contrast, non-primed plants showed significant repression of splicing. These observations link 'splicing memory' to the ability of plants to survive subsequent and otherwise lethal heat stress. This newly discovered priming-induced splicing memory may represent a general feature of heat-stress responses in plants and other organisms as many of the key components are conserved among eukaryotes. Furthermore, this finding could facilitate the development of novel approaches to improve plant survival under extreme heat stress.

Gene clustering pattern, promoter architecture, and gene expression stability in eukaryotic genomes.

A balance between gene expression stability and evolvability is essential for the long-term maintenance of a living system. In this paper, we studied whether the genetic and epigenetic properties of the promoter affect gene expression variability. We hypothesized that upstream distance and orientation (head-to-head or head-to-tail) are important for the promoter architecture and gene expression variability. We found that in budding yeast genes with a short upstream distance tend to have low gene expression variability, and their promoter is flanked by strongly positioned nucleosomes and tends to have low nucleosome occupancy. These observations suggest that in vivo positioning of the flanking nucleosomes facilitates stable nucleosome depletion at the core promoter region and enhances gene expression stability. Head-to-head genes have, on average, lower gene expression variability, greater nucleosome depletion at the core promoter region, and more strongly positioned nucleosomes that flank the core promoter than do head-to-tail genes. These observations hold for diverse eukaryotes. In complex organisms such as mammals, only a small fraction of head-to-tail genes have retained a short upstream distance, probably because the promoter may not be flanked by a strongly positioned nucleosome on the upstream side.

Evolutionary conservation of histone modifications in mammals.

Histone modification is an important mechanism of gene regulation in eukaryotes. Why many histone modifications can be stably maintained in the midst of genetic and environmental changes is a fundamental question in evolutionary biology. We obtained genome-wide profiles of three histone marks, H3 lysine 4 tri-methylation (H3K4me3), H3 lysine 4 mono-methylation (H3K4me1), and H3 lysine 27 acetylation (H3K27ac), for several cell types from human and mouse. We identified histone modifications that were stable among different cell types in human and histone modifications that were evolutionarily conserved between mouse and human in the same cell type. We found that histone modifications that were stable among cell types were also likely to be conserved between species. This trend was consistently observed in promoter, intronic, and intergenic regions for all of the histone marks tested. Importantly, the trend was observed regardless of the expression breadth of the nearby gene, indicating that slow evolution of housekeeping genes was not the major reason for the correlation. These regions showed distinct genetic and epigenetic properties, such as clustered transcription factor binding sites (TFBSs), high GC content, and CTCF binding at flanking sides. Based on our observations, we proposed that TFBS clustering in or near a histone modification plays a significant role in stabilizing and conserving the histone modification because TFBS clustering promotes TFBS conservation, which in turn promotes histone modification conservation. In summary, the results of this study support the view that in mammalian genomes a common mechanism maintains histone modifications against both genetic and environmental (cellular) changes.

NAC transcription factors ATAF1 and ANAC055 affect the heat stress response in Arabidopsis.

Pre-exposing (priming) plants to mild, non-lethal elevated temperature improves their tolerance to a later higher-temperature stress (triggering stimulus), which is of great ecological importance. 'Thermomemory' is maintaining this tolerance for an extended period of time. NAM/ATAF1/2/CUC2 (NAC) proteins are plant-specific transcription factors (TFs) that modulate responses to abiotic stresses, including heat stress (HS). Here, we investigated the potential role of NACs for thermomemory. We determined the expression of 104 Arabidopsis NAC genes after priming and triggering heat stimuli, and found ATAF1 expression is strongly induced right after priming and declines below control levels thereafter during thermorecovery. Knockout mutants of ATAF1 show better thermomemory than wild type, revealing a negative regulatory role. Differential expression analyses of RNA-seq data from ATAF1 overexpressor, ataf1 mutant and wild-type plants after heat priming revealed five genes that might be priming-associated direct targets of ATAF1: AT2G31260 (ATG9), AT2G41640 (GT61), AT3G44990 (XTH31), AT4G27720 and AT3G23540. Based on co-expression analyses applied to the aforementioned RNA-seq profiles, we identified ANAC055 to be transcriptionally co-regulated with ATAF1. Like ataf1, anac055 mutants show improved thermomemory, revealing a potential co-control of both NAC TFs over thermomemory. Our data reveals a core importance of two NAC transcription factors, ATAF1 and ANAC055, for thermomemory.

Chromerid genomes reveal the evolutionary path from photosynthetic algae to obligate intracellular parasites.

The eukaryotic phylum Apicomplexa encompasses thousands of obligate intracellular parasites of humans and animals with immense socio-economic and health impacts. We sequenced nuclear genomes of Chromera velia and Vitrella brassicaformis, free-living non-parasitic photosynthetic algae closely related to apicomplexans. Proteins from key metabolic pathways and from the endomembrane trafficking systems associated with a free-living lifestyle have been progressively and non-randomly lost during adaptation to parasitism. The free-living ancestor contained a broad repertoire of genes many of which were repurposed for parasitic processes, such as extracellular proteins, components of a motility apparatus, and DNA- and RNA-binding protein families. Based on transcriptome analyses across 36 environmental conditions, Chromera orthologs of apicomplexan invasion-related motility genes were co-regulated with genes encoding the flagellar apparatus, supporting the functional contribution of flagella to the evolution of invasion machinery. This study provides insights into how obligate parasites with diverse life strategies arose from a once free-living phototrophic marine alga.

DNA replication timing and selection shape the landscape of nucleotide variation in cancer genomes.

Cancer cells evolve from normal cells by somatic mutations and natural selection. Comparing the evolution of cancer cells and that of organisms can elucidate the genetic basis of cancer. Here we analyse somatic mutations in >400 cancer genomes. We find that the frequency of somatic single-nucleotide variations increases with replication time during the S phase much more drastically than germ-line single-nucleotide variations and somatic large-scale structural alterations, including amplifications and deletions. The ratio of nonsynonymous to synonymous single-nucleotide variations is higher for cancer cells than for germ-line cells, suggesting weaker purifying selection against somatic mutations. Among genes with recurrent mutations only cancer driver genes show evidence of strong positive selection, and late-replicating regions are depleted of cancer driver genes, although enriched for recurrently mutated genes. These observations show that replication timing has a prominent role in shaping the single-nucleotide variation landscape of cancer cells.

Coordinated expression domains in mammalian genomes.

BACKGROUND: Gene order in eukaryotic genomes is not random. Genes showing similar expression (coexpression) patterns are often clustered along the genome. The goal of this study is to characterize coexpression clustering in mammalian genomes and to investigate the underlying mechanisms. METHODOLOGY/PRINCIPAL FINDINGS: We detect clustering of coexpressed genes across multiple scales, from neighboring genes to chromosomal domains that span tens of megabases and, in some cases, entire chromosomes. Coexpression domains may be positively or negatively correlated with other domains, within and between chromosomes. We find that long-range expression domains are associated with gene density, which in turn is related to physical organization of the chromosomes within the nucleus. We show that gene expression changes between healthy and diseased tissue samples occur in a gene density-dependent manner. CONCLUSIONS/SIGNIFICANCE: We demonstrate that coexpression domains exist across multiple scales. We identify potential mechanisms for short-range as well as long-range coexpression domains. We provide evidence that the three-dimensional architecture of the chromosomes may underlie long-range coexpression domains. Chromosome territory reorganization may play a role in common human diseases such as Alzheimer's disease and psoriasis.

Complex oncogenic translocations with gene amplification are initiated by specific DNA breaks in lymphocytes.

Chromosomal instability is a hallmark of many tumor types. Complex chromosomal rearrangements with associated gene amplification, known as complicons, characterize many hematologic and solid cancers. Whereas chromosomal aberrations, including complicons, are useful diagnostic and prognostic cancer markers, their molecular origins are not known. Although accumulating evidence has implicated DNA double-strand break repair in suppression of oncogenic genome instability, the genomic elements required for chromosome rearrangements, especially complex lesions, have not been elucidated. Using a mouse model of B-lineage lymphoma, characterized by complicon formation involving the immunoglobulin heavy chain (Igh) locus and the c-myc oncogene, we have now investigated the requirement for specific genomic segments as donors for complex rearrangements. We now show that specific DNA double-strand breaks, occurring within a narrow segment of Igh, are necessary to initiate complicon formation. By contrast, neither specific DNA breaks nor the powerful intronic enhancer Emu are required for complicon-independent oncogenesis. This study is the first to delineate mechanisms of complex versus simple instability and the first to identify specific chromosomal elements required for complex chromosomal aberrations. These findings will illuminate genomic cancer susceptibility and risk factors.