Genome-wide analysis of the H3K27me3 epigenome and transcriptome in brassica rapa

Background Genome-wide maps of histone modifications have been obtained for several plant species. However, most studies focus on model systems and do not enforce FAIR data management principles. Here we study the H3K27me3 epigenome and associated transcriptome of Brassica rapa, an important vegetable cultivated worldwide. Findings We performed H3K27me3 chromatin immunoprecipitation followed by high-throughput sequencing and transcriptomic analysis by 3′-end RNA sequencing from B. rapa leaves and inflorescences. To analyze these data we developed a Reproducible Epigenomic Analysis pipeline using Galaxy and Jupyter, packaged into Docker images to facilitate transparency and reuse. We found that H3K27me3 covers roughly one-third of all B. rapa protein-coding genes and its presence correlates with low transcript levels. The comparative analysis between leaves and inflorescences suggested that the expression of various floral regulatory genes during development depends on H3K27me3. To demonstrate the importance of H3K27me3 for B. rapa development, we characterized a mutant line deficient in the H3K27 methyltransferase activity. We found that braA.clf mutant plants presented pleiotropic alterations, e.g., curly leaves due to increased expression and reduced H3K27me3 levels at AGAMOUS-like loci. Conclusions We characterized the epigenetic mark H3K27me3 at genome-wide levels and provide genetic evidence for its relevance in B. rapa development. Our work reveals the epigenomic landscape of H3K27me3 in B. rapa and provides novel genomics datasets and bioinformatics analytical resources. We anticipate that this work will lead the way to further epigenomic studies in the complex genome of Brassica crops

Genomic characterization of murine monocytes reveals C/EBPβ transcription factor dependence of Ly6C(-) cells

Monocytes are circulating, short-lived mononuclear phagocytes, which in mice and man comprise two main subpopulations. Murine Ly6C(+) monocytes display developmental plasticity and are recruited to complement tissue-resident macrophages and dendritic cells on demand. Murine vascular Ly6C(-) monocytes patrol the endothelium, act as scavengers, and support vessel wall repair. Here we characterized population and single cell transcriptomes, as well as enhancer and promoter landscapes of the murine monocyte compartment. Single cell RNA-seq and transplantation experiments confirmed homeostatic default differentiation of Ly6C(+) into Ly6C(-) monocytes. The main two subsets were homogeneous, but linked by a more heterogeneous differentiation intermediate. We show that monocyte differentiation occurred through de novo enhancer establishment and activation of pre-established (poised) enhancers. Generation of Ly6C(-) monocytes involved induction of the transcription factor C/EBP{beta} and C/EBP{beta}-deficient mice lacked Ly6C(-) monocytes. Mechanistically, C/EBP{beta} bound the Nr4a1 promoter and controlled expression of this established monocyte survival factor

CRISPR/Cas9-mediated glycolate oxidase disruption is an efficacious and safe treatment for primary hyperoxaluria type I

CRISPR/Cas9 technology offers novel approaches for the development of new therapies for many unmet clinical needs, including a significant number of inherited monogenic diseases. However, in vivo correction of disease-causing genes is still inefficient, especially for those diseases without selective advantage for corrected cells. We reasoned that substrate reduction therapies (SRT) targeting non-essential enzymes could provide an attractive alternative. Here we evaluate the therapeutic efficacy of an in vivo CRISPR/Cas9-mediated SRT to treat primary hyperoxaluria type I (PH1), a rare inborn dysfunction in glyoxylate metabolism that results in excessive hepatic oxalate production causing end-stage renal disease. A single systemic administration of an AAV8-CRISPR/Cas9 vector targeting glycolate oxidase, prevents oxalate overproduction and kidney damage, with no signs of toxicity in Agxt1(-/-) mice. Our results reveal that CRISPR/Cas9-mediated SRT represents a promising therapeutic option for PH1 that can be potentially applied to other metabolic diseases caused by the accumulation of toxic metabolites

In vivo screening characterizes chromatin factor functions during normal and malignant hematopoiesis

Bulk ex vivo and single-cell in vivo CRISPR knockout screens are used to characterize 680 chromatin factors during mouse hematopoiesis, highlighting lineage-specific and normal and leukemia-specific functions. Cellular differentiation requires extensive alterations in chromatin structure and function, which is elicited by the coordinated action of chromatin and transcription factors. By contrast with transcription factors, the roles of chromatin factors in differentiation have not been systematically characterized. Here, we combine bulk ex vivo and single-cell in vivo CRISPR screens to characterize the role of chromatin factor families in hematopoiesis. We uncover marked lineage specificities for 142 chromatin factors, revealing functional diversity among related chromatin factors (i.e. barrier-to-autointegration factor subcomplexes) as well as shared roles for unrelated repressive complexes that restrain excessive myeloid differentiation. Using epigenetic profiling, we identify functional interactions between lineage-determining transcription factors and several chromatin factors that explain their lineage dependencies. Studying chromatin factor functions in leukemia, we show that leukemia cells engage homeostatic chromatin factor functions to block differentiation, generating specific chromatin factor-transcription factor interactions that might be therapeutically targeted. Together, our work elucidates the lineage-determining properties of chromatin factors across normal and malignant hematopoiesis

Functional study of chromatin factors uncovers strong lineage determining roles and divergent behaviours between normal and malignant haematopoiesis

The development of eukaryotic organisms requires the coordination of dynamic regulatory programs at multiple levels, including transcription, post-transcription and translation 1–3 . Understanding the genome regulatory networks is essential to unravel the transcriptional control mechanisms of cell-fate choices and identity during cellular differentiation and development, cell homeostasis, metabolism, inheritance and external stimuli response 4–6 . In addition, genome regulatory elements overlap with a significant fraction of genetic variants associated with pathological conditions such as congenital diseases and cancer 1. Thus, they are currently at the forefront of new therapeutic approaches 1,7,8. When discussing genome regulation, six primary layers of regulation must be considered. These encompass the presence of 1) cis-regulatory elements, including promoters and enhancers, 2) transcription and chromatin factors (TFs and CFs, respectively), and 3) epigenetic tags, such as DNA and histone modifications; 4) the intricate interactions between TFs and CFs in activating or repressing specific genome regulatory elements; 5) long non-coding RNAs (lncRNAs); and 6) the re-configuration of the three-dimensional (3D) chromatin structure 9–11. The integration of these layers culminates in a specific transcriptional program that supports the unique cellular identities

Genome-wide analysis of the H3K27me3 epigenome and transcriptome in brassica rapa

Background Genome-wide maps of histone modifications have been obtained for several plant species. However, most studies focus on model systems and do not enforce FAIR data management principles. Here we study the H3K27me3 epigenome and associated transcriptome of Brassica rapa, an important vegetable cultivated worldwide. Findings We performed H3K27me3 chromatin immunoprecipitation followed by high-throughput sequencing and transcriptomic analysis by 3′-end RNA sequencing from B. rapa leaves and inflorescences. To analyze these data we developed a Reproducible Epigenomic Analysis pipeline using Galaxy and Jupyter, packaged into Docker images to facilitate transparency and reuse. We found that H3K27me3 covers roughly one-third of all B. rapa protein-coding genes and its presence correlates with low transcript levels. The comparative analysis between leaves and inflorescences suggested that the expression of various floral regulatory genes during development depends on H3K27me3. To demonstrate the importance of H3K27me3 for B. rapa development, we characterized a mutant line deficient in the H3K27 methyltransferase activity. We found that braA.clf mutant plants presented pleiotropic alterations, e.g., curly leaves due to increased expression and reduced H3K27me3 levels at AGAMOUS-like loci. Conclusions We characterized the epigenetic mark H3K27me3 at genome-wide levels and provide genetic evidence for its relevance in B. rapa development. Our work reveals the epigenomic landscape of H3K27me3 in B. rapa and provides novel genomics datasets and bioinformatics analytical resources. We anticipate that this work will lead the way to further epigenomic studies in the complex genome of Brassica crops

Conserved and distinct roles of H3K27me3 demethylases regulating flowering time in Brassica rapa

14 Pág.Epigenetic regulation is necessary for optimal organism development and preservation of gene expression profiles in the cell. In plants, the trimethylation of histone H3 lysine 27 (H3K27me3) is a silencing epigenetic mark relevant for developmental transitions like flowering. The floral transition is a key agronomic trait; however, the epigenetic mechanisms of flowering time regulation in crops remain poorly understood. Here we study the Jumonji H3K27me3 demethylases BraA.REF6 and BraA.ELF6 in Brassica rapa. Phenotypic characterization of novel mutant lines and genome-wide H3K27me3 chromatin immunoprecipitation and transcriptomic analyses indicated that BraA.REF6 plays a greater role than BraA.ELF6 in fine-tuning H3K27me3 levels. In addition, we found that braA.elf6 mutants were early flowering due to high H3K27me3 levels at B. rapa homologs of the floral repressor FLC. Unlike mutations in Arabidopsis thaliana, braA.ref6 mutants were late flowering without altering the expression of B. rapa FLC genes. Remarkably, we found that BraA.REF6 regulated a number of gibberellic acid (GA) biosynthetic genes, including a homolog of GA1, and that GA-treatment complemented the late flowering mutant phenotype. This study increases our understanding of the epigenetic regulation of flowering time in B. rapa, highlighting conserved and distinct regulatory mechanisms between model and crop species.This study was supported by grants BIO2015-68031-R and RTI2018-097749-B-I00 to PC, and PID2019-104899GB-I00 to JAJ and MP from MCIN/AEI/10.13039/501100011033 and FEDER; and fellowships BES-2016-078939 to LP and RYC-2013-14689 to PC from MCIN/AEI/10.13039/501100011033 and FSE. MPM was supported by a Postdoctoral contract associated to the “Severo Ochoa Program for Centres of Excellence in R&D” grant SEV-2016-0672 (2017–2021) to CBGP from MCIN/AEI/10.13039/501100011033.Peer reviewe

Functional and transcriptomic analysis of extracellular vesicles identifies calprotectin as a new prognostic marker in peripheral arterial disease (PAD)

Peripheral arterial disease (PAD) is associated with a high risk of cardiovascular events and death and is postulated to be a critical socioeconomic cost in the future. Extracellular vesicles (EVs) have emerged as potential candidates for new biomarker discovery related to their protein and nucleic acid cargo. In search of new prognostic and therapeutic targets in PAD, we determined the prothrombotic activity, the cellular origin and the transcriptomic profile of circulating EVs. This prospective study included control and PAD patients. Coagulation time (Procoag-PPL kit), EVs cellular origin and phosphatidylserine exposure were determined by flow cytometry in plateletfree plasma (n = 45 PAD). Transcriptomic profiles of medium/large EVs were generated using the MARS-Seq RNA-Seq protocol (n = 12/group). The serum concentration of the differentially expressed gene S100A9, in serum calprotectin (S100A8/A9), was validated by ELISA in control (n = 100) and PAD patients (n = 317). S100A9 was also determined in EVs and tissues of human atherosclerotic plaques (n = 3). Circulating EVs of PAD patients were mainly of platelet origin, predominantly Annexin V positive and were associated with the procoagulant activity of plateletfree plasma. Transcriptomic analysis of EVs identified 15 differentially expressed genes. Among them, serum calprotectin was elevated in PAD patients (p < 0.05) and associated with increased amputation risk before and after covariate adjustment (mean follow-up 3.6 years, p < 0.01). The combination of calprotectin with hs-CRP in the multivariate analysis further improved risk stratification (p < 0.01). Furthermore, S100A9 was also expressed in femoral plaque derived EVs and tissues. In summary, we found that PAD patients release EVs, mainly of platelet origin, highly positive for AnnexinV and rich in transcripts related to platelet biology and immune responses. Amputation risk prediction improved with calprotectin and was significantly higher when combined with hs-CRP. Our results suggest that EVs can be a promising component of liquid biopsy to identify the molecular signature of PAD patients

Functional and transcriptomic analysis of extracellular vesicles identifies calprotectin as a new prognostic marker in peripheral arterial disease (PAD)

Peripheral arterial disease (PAD) is associated with a high risk of cardiovascular events and death and is postulated to be a critical socioeconomic cost in the future. Extracellular vesicles (EVs) have emerged as potential candidates for new biomarker discovery related to their protein and nucleic acid cargo. In search of new prognostic and therapeutic targets in PAD, we determined the prothrombotic activity, the cellular origin and the transcriptomic profile of circulating EVs. This prospective study included control and PAD patients. Coagulation time (Procoag-PPL kit), EVs cellular origin and phosphatidylserine exposure were determined by flow cytometry in plateletfree plasma (n = 45 PAD). Transcriptomic profiles of medium/large EVs were generated using the MARS-Seq RNA-Seq protocol (n = 12/group). The serum concentration of the differentially expressed gene S100A9, in serum calprotectin (S100A8/A9), was validated by ELISA in control (n = 100) and PAD patients (n = 317). S100A9 was also determined in EVs and tissues of human atherosclerotic plaques (n = 3). Circulating EVs of PAD patients were mainly of platelet origin, predominantly Annexin V positive and were associated with the procoagulant activity of plateletfree plasma. Transcriptomic analysis of EVs identified 15 differentially expressed genes. Among them, serum calprotectin was elevated in PAD patients (p < 0.05) and associated with increased amputation risk before and after covariate adjustment (mean follow-up 3.6 years, p < 0.01). The combination of calprotectin with hs-CRP in the multivariate analysis further improved risk stratification (p < 0.01). Furthermore, S100A9 was also expressed in femoral plaque derived EVs and tissues. In summary, we found that PAD patients release EVs, mainly of platelet origin, highly positive for AnnexinV and rich in transcripts related to platelet biology and immune responses. Amputation risk prediction improved with calprotectin and was significantly higher when combined with hs-CRP. Our results suggest that EVs can be a promising component of liquid biopsy to identify the molecular signature of PAD patients

CRISPR/Cas9-mediated glycolate oxidase disruption is an efficacious and safe treatment for primary hyperoxaluria type I

CRISPR/Cas9 technology offers novel approaches for the development of new therapies for many unmet clinical needs, including a significant number of inherited monogenic diseases. However, in vivo correction of disease-causing genes is still inefficient, especially for those diseases without selective advantage for corrected cells. We reasoned that substrate reduction therapies (SRT) targeting non-essential enzymes could provide an attractive alternative. Here we evaluate the therapeutic efficacy of an in vivo CRISPR/Cas9-mediated SRT to treat primary hyperoxaluria type I (PH1), a rare inborn dysfunction in glyoxylate metabolism that results in excessive hepatic oxalate production causing end-stage renal disease. A single systemic administration of an AAV8-CRISPR/Cas9 vector targeting glycolate oxidase, prevents oxalate overproduction and kidney damage, with no signs of toxicity in Agxt1(-/-) mice. Our results reveal that CRISPR/Cas9-mediated SRT represents a promising therapeutic option for PH1 that can be potentially applied to other metabolic diseases caused by the accumulation of toxic metabolites