LIPID MAPS: update to databases and tools for the lipidomics community.

LIPID MAPS (LIPID Metabolites and Pathways Strategy), www.lipidmaps.org, provides a systematic and standardized approach to organizing lipid structural and biochemical data. Founded 20 years ago, the LIPID MAPS nomenclature and classification has become the accepted community standard. LIPID MAPS provides databases for cataloging and identifying lipids at varying levels of characterization in addition to numerous software tools and educational resources, and became an ELIXIR-UK data resource in 2020. This paper describes the expansion of existing databases in LIPID MAPS, including richer metadata with literature provenance, taxonomic data and improved interoperability to facilitate FAIR compliance. A joint project funded by ELIXIR-UK, in collaboration with WikiPathways, curates and hosts pathway data, and annotates lipids in the context of their biochemical pathways. Updated features of the search infrastructure are described along with implementation of programmatic access via API and SPARQL. New lipid-specific databases have been developed and provision of lipidomics tools to the community has been updated. Training and engagement have been expanded with webinars, podcasts and an online training school.

Molecular synergy underlies the co-occurrence patterns and phenotype of NPM1-mutant acute myeloid leukemia.

NPM1 mutations define the commonest subgroup of acute myeloid leukemia (AML) and frequently co-occur with FLT3 internal tandem duplications (ITD) or, less commonly, NRAS or KRAS mutations. Co-occurrence of mutant NPM1 with FLT3-ITD carries a significantly worse prognosis than NPM1-RAS combinations. To understand the molecular basis of these observations, we compare the effects of the 2 combinations on hematopoiesis and leukemogenesis in knock-in mice. Early effects of these mutations on hematopoiesis show that compound Npm1cA/+;NrasG12D/+ or Npm1cA;Flt3ITD share a number of features: Hox gene overexpression, enhanced self-renewal, expansion of hematopoietic progenitors, and myeloid differentiation bias. However, Npm1cA;Flt3ITD mutants displayed significantly higher peripheral leukocyte counts, early depletion of common lymphoid progenitors, and a monocytic bias in comparison with the granulocytic bias in Npm1cA/+;NrasG12D/+ mutants. Underlying this was a striking molecular synergy manifested as a dramatically altered gene expression profile in Npm1cA;Flt3ITD , but not Npm1cA/+;NrasG12D/+ , progenitors compared with wild-type. Both double-mutant models developed high-penetrance AML, although latency was significantly longer with Npm1cA/+;NrasG12D/+ During AML evolution, both models acquired additional copies of the mutant Flt3 or Nras alleles, but only Npm1cA/+;NrasG12D/+ mice showed acquisition of other human AML mutations, including IDH1 R132Q. We also find, using primary Cas9-expressing AMLs, that Hoxa genes and selected interactors or downstream targets are required for survival of both types of double-mutant AML. Our results show that molecular complementarity underlies the higher frequency and significantly worse prognosis associated with NPM1c/FLT3-ITD vs NPM1/NRAS-G12D-mutant AML and functionally confirm the role of HOXA genes in NPM1c-driven AML.

Comparative transcriptome analyses indicate molecular homology of zebrafish swimbladder and mammalian lung.

The fish swimbladder is a unique organ in vertebrate evolution and it functions for regulating buoyancy in most teleost species. It has long been postulated as a homolog of the tetrapod lung, but the molecular evidence is scarce. In order to understand the molecular function of swimbladder as well as its relationship with lungs in tetrapods, transcriptomic analyses of zebrafish swimbladder were carried out by RNA-seq. Gene ontology classification showed that genes in cytoskeleton and endoplasmic reticulum were enriched in the swimbladder. Further analyses depicted gene sets and pathways closely related to cytoskeleton constitution and regulation, cell adhesion, and extracellular matrix. Several prominent transcription factor genes in the swimbladder including hoxc4a, hoxc6a, hoxc8a and foxf1 were identified and their expressions in developing swimbladder during embryogenesis were confirmed. By comparison of enriched transcripts in the swimbladder with those in human and mouse lungs, we established the resemblance of transcriptome of the zebrafish swimbladder and mammalian lungs. Based on the transcriptomic data of zebrafish swimbladder, the predominant functions of swimbladder are in its epithelial and muscular tissues. Our comparative analyses also provide molecular evidence of the relatedness of the fish swimbladder and mammalian lung.

CCR4-associated factor 1 coordinates the expression of Plasmodium falciparum egress and invasion proteins.

Coordinated regulation of gene expression is a hallmark of the Plasmodium falciparum asexual blood-stage development cycle. We report that carbon catabolite repressor protein 4 (CCR4)-associated factor 1 (CAF1) is critical in regulating more than 1,000 genes during malaria parasites' intraerythrocytic stages, especially egress and invasion proteins. CAF1 knockout results in mistimed expression, aberrant accumulation and localization of proteins involved in parasite egress, and invasion of new host cells, leading to premature release of predominantly half-finished merozoites, drastically reducing the intraerythrocytic growth rate of the parasite. This study demonstrates that CAF1 of the CCR4-Not complex is a significant gene regulatory mechanism needed for Plasmodium development within the human host.

The transcriptional repressor Blimp1/Prdm1 regulates postnatal reprogramming of intestinal enterocytes.

Female mammals produce milk to feed their newborn offspring before teeth develop and permit the consumption of solid food. Intestinal enterocytes dramatically alter their biochemical signature during the suckling-to-weaning transition. The transcriptional repressor Blimp1 is strongly expressed in immature enterocytes in utero, but these are gradually replaced by Blimp1(-) crypt-derived adult enterocytes. Here we used a conditional inactivation strategy to eliminate Blimp1 function in the developing intestinal epithelium. There was no noticeable effect on gross morphology or formation of mature cell types before birth. However, survival of mutant neonates was severely compromised. Transcriptional profiling experiments reveal global changes in gene expression patterns. Key components of the adult enterocyte biochemical signature were substantially and prematurely activated. In contrast, those required for processing maternal milk were markedly reduced. Thus, we conclude Blimp1 governs the developmental switch responsible for postnatal intestinal maturation.

Complex exon-intron marking by histone modifications is not determined solely by nucleosome distribution.

It has recently been shown that nucleosome distribution, histone modifications and RNA polymerase II (Pol II) occupancy show preferential association with exons ("exon-intron marking"), linking chromatin structure and function to co-transcriptional splicing in a variety of eukaryotes. Previous ChIP-sequencing studies suggested that these marking patterns reflect the nucleosomal landscape. By analyzing ChIP-chip datasets across the human genome in three cell types, we have found that this marking system is far more complex than previously observed. We show here that a range of histone modifications and Pol II are preferentially associated with exons. However, there is noticeable cell-type specificity in the degree of exon marking by histone modifications and, surprisingly, this is also reflected in some histone modifications patterns showing biases towards introns. Exon-intron marking is laid down in the absence of transcription on silent genes, with some marking biases changing or becoming reversed for genes expressed at different levels. Furthermore, the relationship of this marking system with splicing is not simple, with only some histone modifications reflecting exon usage/inclusion, while others mirror patterns of exon exclusion. By examining nucleosomal distributions in all three cell types, we demonstrate that these histone modification patterns cannot solely be accounted for by differences in nucleosome levels between exons and introns. In addition, because of inherent differences between ChIP-chip array and ChIP-sequencing approaches, these platforms report different nucleosome distribution patterns across the human genome. Our findings confound existing views and point to active cellular mechanisms which dynamically regulate histone modification levels and account for exon-intron marking. We believe that these histone modification patterns provide links between chromatin accessibility, Pol II movement and co-transcriptional splicing.

Genomic approaches uncover increasing complexities in the regulatory landscape at the human SCL (TAL1) locus.

The SCL (TAL1) transcription factor is a critical regulator of haematopoiesis and its expression is tightly controlled by multiple cis-acting regulatory elements. To elaborate further the DNA elements which control its regulation, we used genomic tiling microarrays covering 256 kb of the human SCL locus to perform a concerted analysis of chromatin structure and binding of regulatory proteins in human haematopoietic cell lines. This approach allowed us to characterise further or redefine known human SCL regulatory elements and led to the identification of six novel elements with putative regulatory function both up and downstream of the SCL gene. They bind a number of haematopoietic transcription factors (GATA1, E2A LMO2, SCL, LDB1), CTCF or components of the transcriptional machinery and are associated with relevant histone modifications, accessible chromatin and low nucleosomal density. Functional characterisation shows that these novel elements are able to enhance or repress SCL promoter activity, have endogenous promoter function or enhancer-blocking insulator function. Our analysis opens up several areas for further investigation and adds new layers of complexity to our understanding of the regulation of SCL expression.

FRT-seq: amplification-free, strand-specific transcriptome sequencing.

We report an alternative approach to transcriptome sequencing for the Illumina Genome Analyzer, in which the reverse transcription reaction takes place on the flowcell. No amplification is performed during the library preparation, so PCR biases and duplicates are avoided, and because the template is poly(A)(+) RNA rather than cDNA, the resulting sequences are necessarily strand-specific. The method is compatible with paired- or single-end sequencing.

CARM1 is required in embryonic stem cells to maintain pluripotency and resist differentiation.

Histone H3 methylation at R17 and R26 recently emerged as a novel epigenetic mechanism regulating pluripotency in mouse embryos. Blastomeres of four-cell embryos with high H3 methylation at these sites show unrestricted potential, whereas those with lower levels cannot support development when aggregated in chimeras of like cells. Increasing histone H3 methylation, through expression of coactivator-associated-protein-arginine-methyltransferase 1 (CARM1) in embryos, elevates expression of key pluripotency genes and directs cells to the pluripotent inner cell mass. We demonstrate CARM1 is also required for the self-renewal and pluripotency of embryonic stem (ES) cells. In ES cells, CARM1 depletion downregulates pluripotency genes leading to their differentiation. CARM1 associates with Oct4/Pou5f1 and Sox2 promoters that display detectable levels of R17/26 histone H3 methylation. In CARM1 overexpressing ES cells, histone H3 arginine methylation is also at the Nanog promoter to which CARM1 now associates. Such cells express Nanog at elevated levels and delay their response to differentiation signals. Thus, like in four-cell embryo blastomeres, histone H3 arginine methylation by CARM1 in ES cells allows epigenetic modulation of pluripotency.

Functional diversity for REST (NRSF) is defined by in vivo binding affinity hierarchies at the DNA sequence level.

The molecular events that contribute to, and result from, the in vivo binding of transcription factors to their cognate DNA sequence motifs in mammalian genomes are poorly understood. We demonstrate that variations within the DNA sequence motifs that bind the transcriptional repressor REST (NRSF) encode in vivo DNA binding affinity hierarchies that contribute to regulatory function during lineage-specific and developmental programs in fundamental ways. First, canonical sequence motifs for REST facilitate strong REST binding and control functional classes of REST targets that are common to all cell types, whilst atypical motifs participate in weak interactions and control those targets, which are cell- or tissue-specific. Second, variations in REST binding relate directly to variations in expression and chromatin configurations of REST's target genes. Third, REST clearance from its binding sites is also associated with variations in the RE1 motif. Finally, and most surprisingly, weak REST binding sites reside in DNA sequences that show the highest levels of constraint through evolution, thus facilitating their roles in maintaining tissue-specific functions. These relationships have never been reported in mammalian systems for any transcription factor.

The landscape of histone modifications across 1% of the human genome in five human cell lines.

We generated high-resolution maps of histone H3 lysine 9/14 acetylation (H3ac), histone H4 lysine 5/8/12/16 acetylation (H4ac), and histone H3 at lysine 4 mono-, di-, and trimethylation (H3K4me1, H3K4me2, H3K4me3, respectively) across the ENCODE regions. Studying each modification in five human cell lines including the ENCODE Consortium common cell lines GM06990 (lymphoblastoid) and HeLa-S3, as well as K562, HFL-1, and MOLT4, we identified clear patterns of histone modification profiles with respect to genomic features. H3K4me3, H3K4me2, and H3ac modifications are tightly associated with the transcriptional start sites (TSSs) of genes, while H3K4me1 and H4ac have more widespread distributions. TSSs reveal characteristic patterns of both types of modification present and the position relative to TSSs. These patterns differ between active and inactive genes and in particular the state of H3K4me3 and H3ac modifications is highly predictive of gene activity. Away from TSSs, modification sites are enriched in H3K4me1 and relatively depleted in H3K4me3 and H3ac. Comparison between cell lines identified differences in the histone modification profiles associated with transcriptional differences between the cell lines. These results provide an overview of the functional relationship among histone modifications and gene expression in human cells.

Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project.

We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.

Exon array CGH: detection of copy-number changes at the resolution of individual exons in the human genome.

The development of high-throughput screening methods such as array-based comparative genome hybridization (array CGH) allows screening of the human genome for copy-number changes. Current array CGH strategies have limits of resolution that make detection of small (less than a few tens of kilobases) gains or losses of genomic DNA difficult to identify. We report here a significant improvement in the resolution of array CGH, with the development of an array platform that utilizes single-stranded DNA array elements to accurately measure copy-number changes of individual exons in the human genome. Using this technology, we screened 31 patient samples across an array containing a total of 162 exons for five disease genes and detected copy-number changes, ranging from whole-gene deletions and duplications to single-exon deletions and duplications, in 100% of the cases. Our data demonstrate that it is possible to screen the human genome for copy-number changes with array CGH at a resolution that is 2 orders of magnitude higher than that previously reported.