PeakAnalyzer: genome-wide annotation of chromatin binding and modification loci

A soybean bZIP transcription factor is involved in submergence resistance

Although the functions of basic leucine zipper (bZIP) family transcription factors in the regulation of various abiotic stresses are beginning to be unveiled, the precise roles of bZIP proteins in plants coping with submergence stress remain unclear. Here we identified a bZIP gene GmbZIP71-4 from soybean, which localized in the nucleus. The GmbZIP71-4 over-expressed tabocco line showed reduced submergence resistance due to the decreased abscisic acid (ABA) content. GO and KEGG pathway analysis based on chromatin immunoprecipitation assay sequencing (ChIP-seq) indicated that the differences expressed genes between submergence treatment and control groups were specially enriched in plant hormone signal transduction items, especially those in response to ABA. Electrophoretic mobility shift assays (EMSA) demonstrated that GmbZIP71-4 bound to the promoter of GmABF2 gene, which is consistent with the ChIP-qPCR results. GmbZIP71-4 function as a negative regulator of soybean in responding to submergence stress through manipulating ABA signaling pathway. This findings will set a solid foundation for the understanding of submergence resistance in plants.

CONSTANS alters the circadian clock in Arabidopsis thaliana

Plants are sessile organisms that have acquired highly plastic developmental strategies to adapt to the environment. Among these processes, the floral transition is essential to ensure reproductive success and is finely regulated by several internal and external genetic networks. The photoperiodic pathway, which controls plant response to day length, is one of the most important pathways controlling flowering. In Arabidopsis photoperiodic flowering, CONSTANS (CO) is the central gene activating the expression of the florigen FLOWERING LOCUS T (FT) in the leaves at the end of a long day. The circadian clock strongly regulates CO expression. However, to date, no evidence has been reported regarding a feedback loop from the photoperiod pathway back to the circadian clock. Using transcriptional networks, we have identified relevant network motifs regulating the interplay between the circadian clock and the photoperiod pathway. Gene expression, chromatin immunoprecipitation experiments, and phenotypic analysis allowed us to elucidate the role of CO over the circadian clock. Plants with altered CO expression showed a different internal clock period, measured by daily leaf rhythmic movements. We showed that CO upregulates the expression of key genes related to the circadian clock, such as CCA1, LHY, PRR5, and GI, at the end of a long day by binding to specific sites on their promoters. Moreover, a high number of PRR5-repressed target genes are upregulated by CO, and this could explain the phase transition promoted by CO. The CO-PRR5 complex interacts with the bZIP transcription factor HY5 and helps to localize the complex in the promoters of clock genes. Taken together, our results indicate that there may be a feedback loop in which CO communicates back to the circadian clock, providing seasonal information to the circadian system.

N6-methyladenosine RNA methyltransferase CpMTA1 mediates CpAphA mRNA stability through a YTHDF1-dependent m6A modification in the chestnut blight fungus

In eukaryotic cells, N6-methyladenosine (m6A) is the most prevalent RNA epigenetic modification that plays crucial roles in multiple biological processes. Nevertheless, the functions and regulatory mechanisms of m6A in phytopathogenic fungi are poorly understood. Here, we showed that CpMTA1, an m6A methyltransferase in Cryphonectria parasitica, plays a crucial role in fungal phenotypic traits, virulence, and stress tolerance. Furthermore, the acid phosphatase gene CpAphA was implicated to be a target of CpMTA1 by integrated analysis of m6A-seq and RNA-seq, as in vivo RIP assay data confirmed that CpMTA1 directly interacts with CpAphA mRNA. Deletion of CpMTA1 drastically lowered the m6A level of CpAphA and reduced its mRNA expression. Moreover, we found that an m6A reader protein CpYTHDF1 recognizes CpAphA mRNA and increases its stability. Typically, the levels of CpAphA mRNA and protein exhibited a positive correlation with CpMTA1 and CpYTHDF1. Importantly, site-specific mutagenesis demonstrated that the m6A sites, A1306 and A1341, of CpAphA mRNA are important for fungal phenotypic traits and virulence in C. parasitica. Together, our findings demonstrate the essential role of the m6A methyltransferase CpMTA1 in C. parasitica, thereby advancing our understanding of fungal gene regulation through m6A modification.

Antidepressant effect and mechanism of TMP269 on stress-induced depressive-like behavior in mice

TMP269, a class IIA histone deacetylase inhibitor with selectivity, that has a protective effect on the central nervous system, yet its specific mechanism of action remains ambiguous. Although major depressive disorder (MDD) is highly prevalent, its pathophysiology is poorly understood. Recent evidence suggests that histone deacetylase 5 plays a key role in the pathological process of depression and the fact that preclinical studies have shown HDAC5 to be a potential antidepressant target, the search for natural drugs or small molecule compounds that can target HDAC5 may be a potential therapeutic strategy for the treatment of depression. In addition, we examined the role of the Brain-derived neurotrophic factor (BDNF), an important neurotrophic factor for neuronal survival and growth, as a potential downstream target of HDAC5. We found downward revision of HDAC5 levels in the hippocampus ameliorated depressive-like behavior in LH (Learned helplessness) mice. Furthermore, injection of HDAC5 overexpressing adenoviral vectors in the hippocampal dentate gyrus of wild-type mice produced a somewhat depressive-like phenotype. Pharmacological, immunofluorescence and biochemical experiments showed that TMP269 could produce antidepressant effects by inhibiting mouse hippocampal HDAC5 and thus modulating its downstream BDNF. Over all, TMP269 mitigated LH-induced depressive-like behaviors and abnormalities in synapse formation and neurogenesis within the hippocampus. These findings suggest potential beneficial effects of TMP269 on depression.

Temperature-dependent jumonji demethylase modulates flowering time by targeting H3K36me2/3 in Brassica rapa

Global warming has a severe impact on the flowering time and yield of crops. Histone modifications have been well-documented for their roles in enabling plant plasticity in ambient temperature. However, the factor modulating histone modifications and their involvement in habitat adaptation have remained elusive. In this study, through genome-wide pattern analysis and quantitative-trait-locus (QTL) mapping, we reveal that BrJMJ18 is a candidate gene for a QTL regulating thermotolerance in thermotolerant B. rapa subsp. chinensis var. parachinensis (or Caixin, abbreviated to Par). BrJMJ18 encodes an H3K36me2/3 Jumonji demethylase that remodels H3K36 methylation across the genome. We demonstrate that the BrJMJ18 allele from Par (BrJMJ18Par) influences flowering time and plant growth in a temperature-dependent manner via characterizing overexpression and CRISPR/Cas9 mutant plants. We further show that overexpression of BrJMJ18Par can modulate the expression of BrFLC3, one of the five BrFLC orthologs. Furthermore, ChIP-seq and transcriptome data reveal that BrJMJ18Par can regulate chlorophyll biosynthesis under high temperatures. We also demonstrate that three amino acid mutations may account for function differences in BrJMJ18 between subspecies. Based on these findings, we propose a working model in which an H3K36me2/3 demethylase, while not affecting agronomic traits under normal conditions, can enhance resilience under heat stress in Brassica rapa.

Estrogen-related receptor, a molecular target against lepidoptera pests

Until recently, chemical pesticides were one of the most effective means of controlling agricultural pests; therefore, the search for insecticide targets for agricultural pests has been an ongoing problem. Estrogen-related receptors (ERRs) are transcription factors that regulate cellular metabolism and energy homeostasis in animals. Silkworms are highly sensitive to chemical pesticides, making them ideal models for pesticide screening and evaluation. In this study, we detected ERR expression in key organs involved in pesticide metabolism in silkworms (Bombyx mori), including the fat body and midgut. Using ChIP-seq technology, many estrogen- related response elements were identified in the 2000-bp promoter region upstream of metabolism-related genes, almost all of which were potential ERR target genes. The ERR inhibitor, XCT-790, and the endocrine disruptor, bisphenol A, significantly inhibited expression of the ERR target genes, BmTreh-1, BmTret-1, BmPK, BmPFK, and BmHK, in the fat bodies of silkworms, resulting in pupation difficulties in silkworm larvae that ultimately lead to death. In addition, based on the clarification that the ERR can bind to XCT-790, as observed through biofilm interferometry, its three-dimensional spatial structure was predicted, and using molecular docking techniques, small-molecule compounds with a stronger affinity for the ERR were identified. In summary, utilizing the powerful metabolic regulatory function of ERR in Lepidoptera pests, the developed small molecule inhibitors of ERR can be used for future control of Lepidoptera pests.

Epigenetics and suicide: investigating altered H3K14ac unveiled differential expression in ADORA2A, B4GALT2 and MMP14

Background: Environmental factors make an important contribution to suicide. Histone tails are prone to different modifications, leading to changes of chromatin (de)condensation and consequently gene expression. Materials & methods: Level of H3K14ac was studied with chromatin immunoprecipitation followed by high-throughput DNA sequencing. Genes were further validated with RT-qPCR; using hippocampal tissue. Results: We showed lowered H3K14ac levels in individuals who died by suicide. The genes ADORA2A, B4GALT2 and MMP14 showed differential expression in individuals who died by suicide. Identified genetic and protein interactions among genes show interactions with suicide-related genes. Conclusion: Further investigations of histone modifications in association with DNA methylation and miRNA are needed to expand our knowledge of the genes that could significantly contribute to suicide.

Identification and functional analysis of the LEAFY gene in longan flower induction

BACKGROUND

Flowering at the right time is a very important factor affecting the stable annual yield of longan. However, a lack of knowledge of the regulatory mechanism and key genes of longan flowering restricts healthy development of the longan industry. Therefore, identifying relevant genes and analysing their regulatory mechanism are essential for scientific research and longan industry development.

RESULTS

DlLFY (Dimocarpus longan LEAFY) contains a 1167 bp open reading frame and encodes 388 amino acids. The amino acid sequence has a typical LFY/FLO family domain. DlLFY was expressed in all tissues tested, except for the leaf, pericarp, and pulp, with the highest expression occurring in flower buds. Expression of DlLFY was significantly upregulated at the early flower induction stage in "SX" ("Shixia"). The results of subcellular localization and transactivation analysis showed that DlLFY is a typical transcription factor acting as a transcriptional activator. Moreover, overexpression of DlLFY in Arabidopsis promoted early flowering and restrained growth, resulting in reduced plant height and rosette leaf number and area in transgenic plants. DNA affinity purification sequencing (DAP-Seq) analysis showed that 13 flower-related genes corresponding to five homologous genes of Arabidopsis may have binding sites and be putative target genes. Among these five flower-related genes, only AtTFL1 (terminal flower 1) was strongly inhibited in transgenic lines.

CONCLUSION

Taken together, these results indicate that DlLFY plays a pivotal role in controlling longan flowering, possibly by interacting with TFL1.

ZNF643/ZFP69B Exerts Oncogenic Properties and Associates with Cell Adhesion and Immune Processes

The global cancer burden remains high; thus, a better understanding of the molecular mechanisms driving carcinogenesis is needed to improve current prevention and treatment options. We previously detected the ZNF643/ZFP69B gene upregulated in multiple tumors, and we speculated it may play a role in tumor biology. To test this hypothesis, we employed TCGA-centered databases to correlate ZNF643 status with various clinicopathological parameters. We also performed RNA-seq analysis and in vitro studies assessing cancer cell phenotypes, and we searched for ZNF643-bound genomic loci. Our data indicated higher levels of ZNF643 in most analyzed tumors compared to normal samples, possibly due to copy number variations. ZNF643 mRNA correlated with diverse molecular and immune subtypes and clinicopathological features (tumor stage, grade, patient survival). RNA-seq analysis revealed that ZNF643 silencing triggers the deregulation of the genes implicated in various cancer-related processes, such as growth, adhesion, and immune system. Moreover, we observed that ZNF643 positively influences cell cycle, migration, and invasion. Finally, our ChIP-seq analysis indicated that the genes associated with ZNF643 binding are linked to adhesion and immune signaling. In conclusion, our data confirm the oncogenic properties of ZNF643 and pinpoint its impact on cell adhesion and immune processes.

Docosahexaenoic Acid Alters Lipid Metabolism Processes via H3K9ac Epigenetic Modification in Dairy Goat

Goat milk is increasingly recognized by consumers due to its high nutritional value, richness in short- and medium-chain fatty acids, and richness in polyunsaturated fatty acids (PUFA). Exogenous supplementation of docosahexaenoic acid (DHA) is an important approach to increasing the content of PUFA in goat milk. Several studies have reported benefits of dietary DHA in terms of human health, including potential against chronic diseases and tumors. However, the mechanisms whereby an increased supply of DHA regulates mammary cell function is unknown. In this study, we investigated the effect of DHA on lipid metabolism processes in goat mammary epithelial cells (GMEC) and the function of H3K9ac epigenetic modifications in this process. Supplementation of DHA promoted lipid droplet accumulation increased the DHA content and altered fatty acid composition in GMEC. Lipid metabolism processes were altered by DHA supplementation through transcriptional programs in GMEC. ChIP-seq analysis revealed that DHA induced genome-wide H3K9ac epigenetic changes in GMEC. Multiomics analyses (H3K9ac genome-wide screening and RNA-seq) revealed that DHA-induced expression of lipid metabolism genes (FASN, SCD1, FADS1, FADS2, LPIN1, DGAT1, MBOAT2), which were closely related with changes in lipid metabolism processes and fatty acid profiles, were regulated by modification of H3K9ac. In particular, DHA increased the enrichment of H3K9ac in the promoter region of PDK4 and promoted its transcription, while PDK4 inhibited lipid synthesis and activated AMPK signaling in GMEC. The activation of the expression of fatty acid metabolism-related genes FASN, FADS2, and SCD1 and their upstream transcription factor SREBP1 by the AMPK inhibitor was attenuated in PDK4-overexpressing GMEC. In conclusion, DHA alters lipid metabolism processes via H3K9ac modifications and the PDK4-AMPK-SREBP1 signaling axis in goat mammary epithelial cells, providing new insights into the mechanism through which DHA affects mammary cell function and regulates milk fat metabolism.

Igh and Igk loci use different folding principles for V gene recombination due to distinct chromosomal architectures of pro-B and pre-B cells

Extended loop extrusion across the immunoglobulin heavy-chain (Igh) locus facilitates V_H-DJ_H recombination following downregulation of the cohesin-release factor Wapl by Pax5, resulting in global changes in the chromosomal architecture of pro-B cells. Here, we demonstrate that chromatin looping and V_K-J_K recombination at the Igk locus were insensitive to Wapl upregulation in pre-B cells. Notably, the Wapl protein was expressed at a 2.2-fold higher level in pre-B cells compared with pro-B cells, which resulted in a distinct chromosomal architecture with normal loop sizes in pre-B cells. High-resolution chromosomal contact analysis of the Igk locus identified multiple internal loops, which likely juxtapose V_K and J_K elements to facilitate V_K-J_K recombination. The higher Wapl expression in Igμ-transgenic pre-B cells prevented extended loop extrusion at the Igh locus, leading to recombination of only the 6 most 3' proximal V_H genes and likely to allelic exclusion of all other V_H genes in pre-B cells. These results suggest that pro-B and pre-B cells with their distinct chromosomal architectures use different chromatin folding principles for V gene recombination, thereby enabling allelic exclusion at the Igh locus, when the Igk locus is recombined.

Multi-level interaction between HIF and AHR transcriptional pathways in kidney carcinoma

Hypoxia-inducible factor (HIF) and aryl hydrocarbon receptor (AHR) are members of the bHLH-PAS family of transcription factors that underpin cellular responses to oxygen and to endogenous and exogenous ligands, respectively, and have central roles in the pathogenesis of renal cancer. Composed of heterodimers, they share a common HIF-1β/ARNT subunit and similar DNA-binding motifs, raising the possibility of crosstalk between the two transcriptional pathways. Here, we identify both general and locus-specific mechanisms of interaction between HIF and AHR that act both antagonistically and cooperatively. Specifically, we observe competition for the common HIF-1β/ARNT subunit, in cis synergy for chromatin binding, and overlap in their transcriptional targets. Recently, both HIF and AHR inhibitors have been developed for the treatment of solid tumours. However, inhibition of one pathway may promote the oncogenic effects of the other. Therefore, our work raises important questions as to whether combination therapy targeting both of these pro-tumourigenic pathways might show greater efficacy than targeting each system independently.

Regulatory mechanism of MeGI on sexuality in Diospyros oleifera

Dioecy system is an important strategy for maintaining genetic diversity. The transcription factor MeGI, contributes to dioecy by promoting gynoecium development in Diospyros lotus and D. kaki. However, the function of MeGI in D. oleifera has not been identified. In this study, we confirmed that MeGI, cloned from D. oleifera, repressed the androecium development in Arabidopsis thaliana. Subsequently, chromatin immunoprecipitation-sequencing (ChIP-seq), DNA affinity purification-sequencing (DAP-seq), and RNA-seq were used to uncover the gene expression response to MeGI. The results showed that the genes upregulated and downregulated in response to MeGI were mainly enriched in the circadian rhythm-related and flavonoid biosynthetic pathways, respectively. Additionally, the WRKY DNA-binding protein 28 (WRKY28) gene, which was detected by ChIP-seq, DAP-seq, and RNA-seq, was emphasized. WRKY28 has been reported to inhibit salicylic acid (SA) biosynthesis and was upregulated in MeGI-overexpressing A. thaliana flowers, suggesting that MeGI represses the SA level by increasing the expression level of WRKY28. This was confirmed that SA level was lower in D. oleifera female floral buds than male. Overall, our findings indicate that the MeGI mediates its sex control function in D. oleifera mainly by regulating genes in the circadian rhythm, SA biosynthetic, and flavonoid biosynthetic pathways.

The ULK3 kinase is a determinant of keratinocyte self-renewal and tumorigenesis targeting the arginine methylome

Epigenetic mechanisms oversee epidermal homeostasis and oncogenesis. The identification of kinases controlling these processes has direct therapeutic implications. We show that ULK3 is a nuclear kinase with elevated expression levels in squamous cell carcinomas (SCCs) arising in multiple body sites, including skin and Head/Neck. ULK3 loss by gene silencing or deletion reduces proliferation and clonogenicity of human keratinocytes and SCC-derived cells and affects transcription impinging on stem cell-related and metabolism programs. Mechanistically, ULK3 directly binds and regulates the activity of two histone arginine methyltransferases, PRMT1 and PRMT5 (PRMT1/5), with ULK3 loss compromising PRMT1/5 chromatin association to specific genes and overall methylation of histone H4, a shared target of these enzymes. These findings are of translational significance, as downmodulating ULK3 by RNA interference or locked antisense nucleic acids (LNAs) blunts the proliferation and tumorigenic potential of SCC cells and promotes differentiation in two orthotopic models of skin cancer.

Epigenetic association study uncovered H3K27 acetylation enhancers and dysregulated genes in high-fat-diet-induced nonalcoholic fatty liver disease in rats

Aim: To evaluate the regulatory landscape underlying the active enhancer marked by H3K27ac in high-fat diet (HFD)-induced nonalcoholic fatty liver disease (NAFLD) in rats. Materials & methods: H3K27ac chromatin immunoprecipitation and high-throughput RNA sequencing to construct regulatory profiles and transcriptome of liver from NAFLD rat model induced by HFD. De novo motif analysis for differential H3K27ac peaks. Functional enrichment, Kyoto Encyclopedia of Genes and Genomes pathway and protein-protein interaction network were examined for differential peak-genes. The mechanism was further verified by western blot, chromatin immunoprecipitation-quantitative PCR and real-time PCR. Results: A total of 1831 differential H3K27ac peaks were identified significantly correlating with transcription factors and target genes (CYP8B1, PLA2G12B, SLC27A5, CYP7A1 and APOC3) involved in lipid and energy homeostasis. Conclusion: Altered acetylation induced by HFD leads to the dysregulation of gene expression, further elucidating the epigenetic mechanism in the etiology of NAFLD.

A comprehensive survey of the approaches for pathway analysis using multi-omics data integration

Pathway analysis has been widely used to detect pathways and functions associated with complex disease phenotypes. The proliferation of this approach is due to better interpretability of its results and its higher statistical power compared with the gene-level statistics. A plethora of pathway analysis methods that utilize multi-omics setup, rather than just transcriptomics or proteomics, have recently been developed to discover novel pathways and biomarkers. Since multi-omics gives multiple views into the same problem, different approaches are employed in aggregating these views into a comprehensive biological context. As a result, a variety of novel hypotheses regarding disease ideation and treatment targets can be formulated. In this article, we review 32 such pathway analysis methods developed for multi-omics and multi-cohort data. We discuss their availability and implementation, assumptions, supported omics types and databases, pathway analysis techniques and integration strategies. A comprehensive assessment of each method's practicality, and a thorough discussion of the strengths and drawbacks of each technique will be provided. The main objective of this survey is to provide a thorough examination of existing methods to assist potential users and researchers in selecting suitable tools for their data and analysis purposes, while highlighting outstanding challenges in the field that remain to be addressed for future development.

Monoallelic Heb/Tcf12 Deletion Reduces the Requirement for NOTCH1 Hyperactivation in T-Cell Acute Lymphoblastic Leukemia

Early T-cell development is precisely controlled by E proteins, that indistinguishably include HEB/TCF12 and E2A/TCF3 transcription factors, together with NOTCH1 and pre-T cell receptor (TCR) signalling. Importantly, perturbations of early T-cell regulatory networks are implicated in leukemogenesis. NOTCH1 gain of function mutations invariably lead to T-cell acute lymphoblastic leukemia (T-ALL), whereas inhibition of E proteins accelerates leukemogenesis. Thus, NOTCH1, pre-TCR, E2A and HEB functions are intertwined, but how these pathways contribute individually or synergistically to leukemogenesis remain to be documented. To directly address these questions, we leveraged Cd3e-deficient mice in which pre-TCR signaling and progression through β-selection is abrogated to dissect and decouple the roles of pre-TCR, NOTCH1, E2A and HEB in SCL/TAL1-induced T-ALL, via the use of Notch1 gain of function transgenic (Notch1ICtg) and Tcf12+/- or Tcf3+/- heterozygote mice. As a result, we now provide evidence that both HEB and E2A restrain cell proliferation at the β-selection checkpoint while the clonal expansion of SCL-LMO1-induced pre-leukemic stem cells in T-ALL is uniquely dependent on Tcf12 gene dosage. At the molecular level, HEB protein levels are decreased via proteasomal degradation at the leukemic stage, pointing to a reversible loss of function mechanism. Moreover, in SCL-LMO1-induced T-ALL, loss of one Tcf12 allele is sufficient to bypass pre-TCR signaling which is required for Notch1 gain of function mutations and for progression to T-ALL. In contrast, Tcf12 monoallelic deletion does not accelerate Notch1IC-induced T-ALL, indicating that Tcf12 and Notch1 operate in the same pathway. Finally, we identify a tumor suppressor gene set downstream of HEB, exhibiting significantly lower expression levels in pediatric T-ALL compared to B-ALL and brain cancer samples, the three most frequent pediatric cancers. In summary, our results indicate a tumor suppressor function of HEB/TCF12 in T-ALL to mitigate cell proliferation controlled by NOTCH1 in pre-leukemic stem cells and prevent NOTCH1-driven progression to T-ALL.

CaAIL1 Acts Positively in Pepper Immunity against Ralstonia solanacearum by Repressing Negative Regulators

APETALA2 (AP2) subfamily transcription factors participate in plant growth and development, but their roles in plant immunity remain unclear. Here, we discovered that the AP2 transcription factor CaAIL1 functions in immunity against Ralstonia solanacearum infection (RSI) in pepper (Capsicum annuum). CaAIL1 expression was upregulated by RSI, and loss- and gain-of-function assays using virus-induced gene silencing and transient overexpression, respectively, revealed that CaAIL1 plays a positive role in immunity to RSI in pepper. Chromatin immunoprecipitation sequencing (ChIP-seq) uncovered a subset of transcription-factor-encoding genes, including CaRAP2-7, CaGATA17, CaGtf3a and CaTCF25, that were directly targeted by CaAIL1 via their cis-elements, such as GT or AGGCA motifs. ChIP-qPCR and electrophoretic mobility shift assays confirmed these findings. These genes, encoding transcription factors with negative roles in immunity, were repressed by CaAIL1 during pepper response to RSI, whereas genes encoding positive immune regulators such as CaEAS were derepressed by CaAIL1. Importantly, we showed that the atypical EAR motif (LXXLXXLXX) in CaAIL1 is indispensable for its function in immunity. These findings indicate that CaAIL1 enhances the immunity of pepper against RSI by repressing a subset of negative immune regulators during the RSI response through its binding to several cis-elements in their promoters.

Characterization of Chromatin Accessibility and Gene Expression upon Cold Stress Reveals that the RAV1 Transcription Factor Functions in Cold Response in Vitis Amurensis

Cold tolerance is regulated by a variety of transcription factors (TFs) and their target genes. Except for the well-characterized C-repeat binding factors (CBFs)-dependent transcriptional cascade, the mechanisms of cold tolerance mediated by other transcriptional regulatory networks are still largely unknown. Here, we used the assay for transposase-accessible chromatin with sequencing (ATAC-seq) and RNA-seq to identify cold responsive TFs in Vitis amurensis, a grape species with high cold hardiness. Nine TFs, including CBF4, RAV1 and ERF104, were identified after cold treatment. Weighted gene co-expression network analysis (WGCNA) and gene ontology (GO) analysis revealed that these TFs may regulate cold response through different pathways. As a prime candidate TF, overexpression of VaRAV1 in grape cells improved its cold tolerance. The transgenic cells exhibited low electrolyte leakage and malondialdehyde content and high peroxidase activity. Moreover, the TF gene TCP8 and a gene involving in homogalacturonan biosynthesis were found to be regulated by VaRAV1, suggesting that the contribution of VaRAV1 to cold tolerance may be achieved by enhancing the stability of cell membrane and regulating the expression of target genes involved in plant cell wall composition. Our work provides novel insights into plant response to cold stress and demonstrates the utility of ATAC-seq and RNA-seq for the rapid identification of TFs in response to cold stress in grapevine. VaRAV1 may play an important role in adaption to cold stress.

Pre-existing chromatin accessibility and gene expression differences among naive CD4+ T cells influence effector potential

CD4+ T cells have a remarkable potential to differentiate into diverse effector lineages following activation. Here, we probe the heterogeneity present among naive CD4+ T cells before encountering their cognate antigen to ask whether their effector potential is modulated by pre-existing transcriptional and chromatin landscape differences. Single-cell RNA sequencing shows that key drivers of variability are genes involved in T cell receptor (TCR) signaling. Using CD5 expression as a readout of the strength of tonic TCR interactions with self-peptide MHC, and sorting on the ends of this self-reactivity spectrum, we find that pre-existing transcriptional differences among naive CD4+ T cells impact follicular helper T (TFH) cell versus non-TFH effector lineage choice. Moreover, our data implicate TCR signal strength during thymic development in establishing differences in naive CD4+ T cell chromatin landscapes that ultimately shape their effector potential.

Six1 promotes skeletal muscle thyroid hormone response through regulation of the MCT10 transporter

BACKGROUND

The Six1 transcription factor is implicated in controlling the development of several tissue types, notably skeletal muscle. Six1 also contributes to muscle metabolism and its activity is associated with the fast-twitch, glycolytic phenotype. Six1 regulates the expression of certain genes of the fast muscle program by directly stimulating their transcription or indirectly acting through a long non-coding RNA. We hypothesized that additional mechanisms of action of Six1 might be at play.

METHODS

A combined analysis of gene expression profiling and genome-wide location analysis data was performed. Results were validated using in vivo RNA interference loss-of-function assays followed by measurement of gene expression by RT-PCR and transcriptional reporter assays.

RESULTS

The Slc16a10 gene, encoding the thyroid hormone transmembrane transporter MCT10, was identified as a gene with a transcriptional enhancer directly bound by Six1 and requiring Six1 activity for full expression in adult mouse tibialis anterior, a predominantly fast-twitch muscle. Of the various thyroid hormone transporters, MCT10 mRNA was found to be the most abundant in skeletal muscle, and to have a stronger expression in fast-twitch compared to slow-twitch muscle groups. Loss-of-function of MCT10 in the tibialis anterior recapitulated the effect of Six1 on the expression of fast-twitch muscle genes and led to lower activity of a thyroid hormone receptor-dependent reporter gene.

CONCLUSIONS

These results shed light on the molecular mechanisms controlling the tissue expression profile of MCT10 and identify modulation of the thyroid hormone signaling pathway as an additional mechanism by which Six1 influences skeletal muscle metabolism.

Unique features of the m6A methylome and its response to drought stress in sea buckthorn (Hippophae rhamnoides Linn.)

In plants, recent studies have revealed that N6-methyladenosine (m⁶A) methylation of mRNA has potential regulatory functions of this mRNA modification in many biological processes. m⁶A methyltransferase, m⁶A demethylase and m⁶A-binding proteins can cause differential phenotypes, indicating that m⁶A may have critical roles in the plant. In this study, we depicted the m⁶A map of sea buckthorn (Hippophae rhamnoides Linn.) transcriptome. Similar to A. thaliana, m⁶A sites of sea buckthorn transcriptome is significantly enriched around the stop codon and within 3'-untranslated regions (3'UTR). Gene ontology analysis shows that the m⁶A modification genes are associated with metabolic biosynthesis. In addition, we identified 13,287 different m⁶A peaks (DMPs) between leaf under drought (TR) and control (CK) treatment. It reveals that m⁶A has a high level of conservation and has a positive correlation with mRNA abundance in plants. GO and KEGG enrichment results showed that DMP modification DEGs in TR were particularly associated with ABA biosynthesis. Interestingly, our results showed three m⁶A demethylase (HrALKBH10B, HrALKBH10C and HrALKBH10D) genes were significantly increased following drought stress, which indicated that it may contributed the decreased m⁶A levels. This exhaustive m⁶A map provides a basis and resource for the further functional study of mRNA m⁶A modification in abiotic stress.

Transcriptome-wide N6-methyladenosine (m6A) methylation in watermelon under CGMMV infection

BACKGROUND

Cucumber green mottle mosaic virus (CGMMV) causes substantial global losses in cucurbit crops, especially watermelon. N6-methyladenosine (m⁶A) methylation in RNA is one of the most important post-transcriptional modification mechanisms in eukaryotes. It has been shown to have important regulatory functions in some model plants, but there has been no research regarding m⁶A modifications in watermelon.

RESULTS

We measured the global m⁶A level in resistant watermelon after CGMMV infection using a colorimetric method. And the results found that the global m⁶A level significantly decreased in resistant watermelon after CGMMV infection. Specifically, m⁶A libraries were constructed for the resistant watermelon leaves collected 48 h after CGMMV infection and the whole-genome m⁶A-seq were carried out. Numerous m⁶A modified peaks were identified from CGMMV-infected and control (uninfected) samples. The modification distributions and motifs of these m⁶A peaks were highly conserved in watermelon transcripts but the modification was more abundant than in other reported crop plants. In early response to CGMMV infection, 422 differentially methylated genes (DMGs) were identified, most of which were hypomethylated, and probably associated with the increased expression of watermelon m⁶A demethylase gene ClALKBH4B. Gene Ontology (GO) analysis indicated quite a few DMGs were involved in RNA biology and stress responsive pathways. Combined with RNA-seq analysis, there was generally a negative correlation between m⁶A RNA methylation and transcript level in the watermelon transcriptome. Both the m⁶A methylation and transcript levels of 59 modified genes significantly changed in response to CGMMV infection and some were involved in plant immunity.

CONCLUSIONS

Our study represents the first comprehensive characterization of m⁶A patterns in the watermelon transcriptome and helps to clarify the roles and regulatory mechanisms of m⁶A modification in watermelon in early responses to CGMMV.

AnnoMiner is a new web-tool to integrate epigenetics, transcription factor occupancy and transcriptomics data to predict transcriptional regulators

Gene expression regulation requires precise transcriptional programs, led by transcription factors in combination with epigenetic events. Recent advances in epigenomic and transcriptomic techniques provided insight into different gene regulation mechanisms. However, to date it remains challenging to understand how combinations of transcription factors together with epigenetic events control cell-type specific gene expression. We have developed the AnnoMiner web-server, an innovative and flexible tool to annotate and integrate epigenetic, and transcription factor occupancy data. First, AnnoMiner annotates user-provided peaks with gene features. Second, AnnoMiner can integrate genome binding data from two different transcriptional regulators together with gene features. Third, AnnoMiner offers to explore the transcriptional deregulation of genes nearby, or within a specified genomic region surrounding a user-provided peak. AnnoMiner’s fourth function performs transcription factor or histone modification enrichment analysis for user-provided gene lists by utilizing hundreds of public, high-quality datasets from ENCODE for the model organisms human, mouse, Drosophila and C. elegans. Thus, AnnoMiner can predict transcriptional regulators for a studied process without the strict need for chromatin data from the same process. We compared AnnoMiner to existing tools and experimentally validated several transcriptional regulators predicted by AnnoMiner to indeed contribute to muscle morphogenesis in Drosophila. AnnoMiner is freely available at http://chimborazo.ibdm.univ-mrs.fr/AnnoMiner/.

The transcription factor KNAT2/6b mediates changes in plant architecture in response to drought via down-regulating GA20ox1 in Populus alba × P. glandulosa

Plant architecture is genetically controlled, but is influenced by environmental factors. Plants have evolved adaptive mechanisms that allow changes in their architecture under stress, in which phytohormones play a central role. However, the gene regulators that connect growth and stress signals are rarely reported. Here, we report that a class I KNOX gene, PagKNAT2/6b, can directly inhibit the synthesis of gibberellin (GA), altering plant architecture and improving drought resistance in Populus. Expression of PagKNAT2/6b was significantly induced under drought conditions, and transgenic poplars overexpressing PagKNAT2/6b exhibited shorter internode length and smaller leaf size with short or even absent petioles. Interestingly, these transgenic plants showed improved drought resistance under both short- and long-term drought stress. Histological observations indicated that decreased internode length and leaf size were mainly caused by the inhibition of cell elongation and expansion. GA content was reduced, and the GA20-oxidase gene PagGA20ox1 was down-regulated in overexpressing plants. Expression of PagGA20ox1 was negatively related to that of PagKNAT2/6b under drought stress. ChIP and transient transcription activity assays revealed that PagGA20ox1 was directly targeted by PagKNAT2/6b. Therefore, this study provides evidence that PagKNAT2/6b mediates stress signals and changes in plant architecture via GA signaling by down-regulating PagGA20ox1.

Transcriptional regulation of N6-methyladenosine orchestrates sex-dimorphic metabolic traits

Males and females exhibit striking differences in the prevalence of metabolic traits including hepatic steatosis, a key driver of cardiometabolic morbidity and mortality. RNA methylation is a widespread regulatory mechanism of transcript turnover. Here, we show that presence of the RNA modification N⁶-methyladenosine (m⁶A) triages lipogenic transcripts for degradation and guards against hepatic triglyceride accumulation. In male but not female mice, this protective checkpoint stalls under lipid-rich conditions. Loss of m⁶A control in male livers increases hepatic triglyceride stores, leading to a more 'feminized' hepatic lipid composition. Crucially, liver-specific deletion of the m⁶A complex protein Mettl14 from male and female mice significantly diminishes sex-specific differences in steatosis. We further surmise that the m⁶A installing machinery is subject to transcriptional control by the sex-responsive BCL6-STAT5 axis in response to dietary conditions. These data show that m⁶A is essential for precise and synchronized control of lipogenic enzyme activity and provide insights into the molecular basis for the existence of sex-specific differences in hepatic lipid traits.

N6-methyladenosine RNA modification regulates strawberry fruit ripening in an ABA-dependent manner

BACKGROUND

Epigenetic mark such as DNA methylation plays pivotal roles in regulating ripening of both climacteric and non-climacteric fruits. However, it remains unclear whether mRNA m6A methylation, which has been shown to regulate ripening of the tomato, a typical climacteric fruit, is functionally conserved for ripening control among different types of fruits.

RESULTS

Here we show that m6A methylation displays a dramatic change at ripening onset of strawberry, a classical non-climacteric fruit. The m6A modification in coding sequence (CDS) regions appears to be ripening-specific and tends to stabilize the mRNAs, whereas m6A around the stop codons and within the 3' untranslated regions is generally negatively correlated with the abundance of associated mRNAs. We identified thousands of transcripts with m6A hypermethylation in the CDS regions, including those of NCED5, ABAR, and AREB1 in the abscisic acid (ABA) biosynthesis and signaling pathway. We demonstrate that the methyltransferases MTA and MTB are indispensable for normal ripening of strawberry fruit, and MTA-mediated m6A modification promotes mRNA stability of NCED5 and AREB1, while facilitating translation of ABAR.

CONCLUSION

Our findings uncover that m6A methylation regulates ripening of the non-climacteric strawberry fruit by targeting the ABA pathway, which is distinct from that in the climacteric tomato fruit.

ChIP-seq Analysis of the Global Regulator Vfr Reveals Novel Insights Into the Biocontrol Agent Pseudomonas protegens FD6

Many Pseudomonas protegens strains produce the antibiotics pyoluteorin (PLT) and 2,4-diacetylphloroglucinol (2,4-DAPG), both of which have antimicrobial properties. The biosynthesis of these metabolites is typically controlled by multiple regulatory factors. Virulence factor regulator (Vfr) is a multifunctional DNA-binding regulator that modulates 2,4-DAPG biosynthesis in P. protegens FD6. However, the mechanism by which Vfr regulates this process remains unclear. In the present study, chromatin immunoprecipitation of FLAG-tagged Vfr and nucleotide sequencing analysis were used to identify 847 putative Vfr binding sites in P. protegens FD6. The consensus P. protegens Vfr binding site predicted from nucleotide sequence alignment is TCACA. The qPCR data showed that Vfr positively regulates the expression of phlF and phlG, and the expression of these genes was characterized in detail. The purified recombinant Vfr bound to an approximately 240-bp fragment within the phlF and phlG upstream regions that harbor putative Vfr consensus sequences. Using electrophoretic mobility shift assays, we localized Vfr binding to a 25-bp fragment that contains part of the Vfr binding region. Vfr binding was eliminated by mutating the TACG and CACA sequences in phlF and phlG, respectively. Taken together, our results show that Vfr directly regulates the expression of the 2,4-DAPG operon by binding to the upstream regions of both the phlF and phlG genes. However, unlike other Vfr-targeted genes, Vfr binding to P. protegens FD6 does not require an intact binding consensus motif. Furthermore, we demonstrated that vfr expression is autoregulated in this bacterium. These results provide novel insights into the regulatory role of Vfr in the biocontrol agent P. protegens.

CsLOB1 regulates susceptibility to citrus canker through promoting cell proliferation in citrus

Citrus sinensis lateral organ boundary 1 (CsLOB1) was previously identified as a critical disease susceptibility gene for citrus bacterial canker, which is caused by Xanthomonas citri subsp. citri (Xcc). However, the molecular mechanisms of CsLOB1 in citrus response to Xcc are still elusive. Here, we constructed transgenic plants overexpressing and RNAi-silencing of CsLOB1 using the canker-disease susceptible 'wanjincheng' orange (C. sinensis Osbeck) as explants. CsLOB1-overexpressing plants exhibited dwarf phenotypes with smaller and thicker leaf, increased branches and adventitious buds clustered on stems. These phenotypes were followed by a process of pustule- and canker-like development that exhibited enhanced cell proliferation. Pectin depolymerization and expansin accumulation were enhanced by CsLOB1 overexpression, while cellulose and hemicellulose synthesis were increased by CsLOB1 silence. Whilst overexpression of CsLOB1 increased susceptibility, RNAi-silencing of CsLOB1 enhanced resistance to canker disease without impairing pathogen entry. Transcriptome analysis revealed that CsLOB1 positively regulated cell wall degradation and modification processes, cytokinin metabolism, and cell division. Additionally, 565 CsLOB1-targeted genes were identified in chromatin immunoprecipitation-sequencing (ChIP-seq) experiments. Motif discovery analysis revealed that the most highly overrepresented binding sites had a conserved 6-bp 'GCGGCG' consensus DNA motif. RNA-seq and ChIP-seq data suggested that CsLOB1 directly activates the expression of four genes involved in cell wall remodeling, and three genes that participate in cytokinin and brassinosteroid hormone pathways. Our findings indicate that CsLOB1 promotes cell proliferation by mechanisms depending on cell wall remodeling and phytohormone signaling, which may be critical to citrus canker development and bacterial growth in citrus.

Neuronal enhancers are hotspots for DNA single-strand break repair

Defects in DNA repair frequently lead to neurodevelopmental and neurodegenerative diseases, underscoring the particular importance of DNA repair in long-lived post-mitotic neurons^1,2. The cellular genome is subjected to a constant barrage of endogenous DNA damage, but surprisingly little is known about the identity of the lesion(s) that accumulate in neurons and whether they accrue throughout the genome or at specific loci. Here we show that post-mitotic neurons accumulate unexpectedly high levels of DNA single-strand breaks (SSBs) at specific sites within the genome. Genome-wide mapping reveals that SSBs are located within enhancers at or near CpG dinucleotides and sites of DNA demethylation. These SSBs are repaired by PARP1 and XRCC1-dependent mechanisms. Notably, deficiencies in XRCC1-dependent short-patch repair increase DNA repair synthesis at neuronal enhancers, whereas defects in long-patch repair reduce synthesis. The high levels of SSB repair in neuronal enhancers are therefore likely to be sustained by both short-patch and long-patch processes. These data provide the first evidence of site- and cell-type-specific SSB repair, revealing unexpected levels of localized and continuous DNA breakage in neurons. In addition, they suggest an explanation for the neurodegenerative phenotypes that occur in patients with defective SSB repair.

Tlx3 Exerts Direct Control in Specifying Excitatory Over Inhibitory Neurons in the Dorsal Spinal Cord

The spinal cord dorsal horn is a major station for integration and relay of somatosensory information and comprises both excitatory and inhibitory neuronal populations. The homeobox gene Tlx3 acts as a selector gene to control the development of late-born excitatory (dILB) neurons by specifying glutamatergic transmitter fate in dorsal spinal cord. However, since Tlx3 direct transcriptional targets remain largely unknown, it remains to be uncovered how Tlx3 functions to promote excitatory cell fate. Here we combined a genomics approach based on chromatin immunoprecipitation followed by next generation sequencing (ChIP-seq) and expression profiling, with validation experiments in Tlx3 null embryos, to characterize the transcriptional program of Tlx3 in mouse embryonic dorsal spinal cord. We found most dILB neuron specific genes previously identified to be directly activated by Tlx3. Surprisingly, we found Tlx3 also directly represses many genes associated with the alternative inhibitory dILA neuronal fate. In both cases, direct targets include transcription factors and terminal differentiation genes, showing that Tlx3 directly controls cell identity at distinct levels. Our findings provide a molecular frame for the master regulatory role of Tlx3 in developing glutamatergic dILB neurons. In addition, they suggest a novel function for Tlx3 as direct repressor of GABAergic dILA identity, pointing to how generation of the two alternative cell fates being tightly coupled.

A neural m6A/Ythdf pathway is required for learning and memory in Drosophila

Epitranscriptomic modifications can impact behavior. Here, we used Drosophila melanogaster to study N⁶-methyladenosine (m⁶A), the most abundant modification of mRNA. Proteomic and functional analyses confirm its nuclear (Ythdc1) and cytoplasmic (Ythdf) YTH domain proteins as major m⁶A binders. Assays of short term memory in m⁶A mutants reveal neural-autonomous requirements of m⁶A writers working via Ythdf, but not Ythdc1. Furthermore, m⁶A/Ythdf operate specifically via the mushroom body, the center for associative learning. We map m⁶A from wild-type and Mettl3 mutant heads, allowing robust discrimination of Mettl3-dependent m⁶A sites that are highly enriched in 5' UTRs. Genomic analyses indicate that Drosophila m⁶A is preferentially deposited on genes with low translational efficiency and that m⁶A does not affect RNA stability. Nevertheless, functional tests indicate a role for m⁶A/Ythdf in translational activation. Altogether, our molecular genetic analyses and tissue-specific m⁶A maps reveal selective behavioral and regulatory defects for the Drosophila Mettl3/Ythdf pathway.

Base-resolution models of transcription-factor binding reveal soft motif syntax

The arrangement (syntax) of transcription factor (TF) binding motifs is an important part of the cis-regulatory code, yet remains elusive. We introduce a deep learning model, BPNet, that uses DNA sequence to predict base-resolution chromatin immunoprecipitation (ChIP)-nexus binding profiles of pluripotency TFs. We develop interpretation tools to learn predictive motif representations and identify soft syntax rules for cooperative TF binding interactions. Strikingly, Nanog preferentially binds with helical periodicity, and TFs often cooperate in a directional manner, which we validate using clustered regularly interspaced short palindromic repeat (CRISPR)-induced point mutations. Our model represents a powerful general approach to uncover the motifs and syntax of cis-regulatory sequences in genomics data.

The forkhead transcription factor FOXK2 premarks lineage-specific genes in human embryonic stem cells for activation during differentiation

Enhancers play important roles in controlling gene expression in a choreographed spatial and temporal manner during development. However, it is unclear how these regulatory regions are established during differentiation. Here we investigated the genome-wide binding profile of the forkhead transcription factor FOXK2 in human embryonic stem cells (ESCs) and downstream cell types. This transcription factor is bound to thousands of regulatory regions in human ESCs, and binding at many sites is maintained as cells differentiate to mesendodermal and neural precursor cell (NPC) types, alongside the emergence of new binding regions. FOXK2 binding is generally associated with active histone marks in any given cell type. Furthermore newly acquired, or retained FOXK2 binding regions show elevated levels of activating histone marks following differentiation to NPCs. In keeping with this association with activating marks, we demonstrate a role for FOXK transcription factors in gene activation during NPC differentiation. FOXK2 occupancy in ESCs is therefore an early mark for delineating the regulatory regions, which become activated in later lineages.

Dysregulated H3K27 Acetylation Is Implicated in Fatty Liver Hemorrhagic Syndrome in Chickens

Epigenetic regulation of gene expression has been reported in the pathogenesis of metabolic disorders such as diabetes and liver steatosis in humans. However, the molecular mechanisms of fatty liver hemorrhagic syndrome (FLHS) in chickens have been rarely studied. H3K27ac chromatin immunoprecipitation coupled with high-throughput sequencing and high-throughput RNA sequencing was performed to compare genome-wide H3K27ac profiles and transcriptomes of liver tissue between healthy and FLHS chickens. In total, 1,321 differential H3K27ac regions and 443 differentially expressed genes were identified (| log2Fold change| ≥ 1 and P-value ≤ 0.05) between the two groups. Binding motifs for transcription factors involved in immune processes and metabolic homeostasis were enriched among those differential H3K27ac regions. Differential H3K27ac peaks were associated with multiple known FLHS risk genes, involved in lipid and energy metabolism (PCK1, APOA1, ANGPTL4, and FABP1) and the immune system (FGF7, PDGFRA, and KIT). Previous studies and our current results suggested that the high-energy, low-protein (HELP) diet might have an impact on histone modification and chromatin structure, leading to the dysregulation of candidate genes and the peroxisome proliferator-activated receptor (PPAR) signaling pathway, which causes excessive accumulation of fat in the liver tissue and induces the development of FLHS. These findings highlight that epigenetic modifications contribute to the regulation of gene expression and play a central regulatory role in FLHS. The PPAR signaling pathway and other genes implicated in FLHS are of great importance for the development of novel and specific therapies for FLHS-susceptible commercial laying hens.

Reciprocal activation of HEY1 and NOTCH4 under SOX2 control promotes EMT in head and neck squamous cell carcinoma

Several comprehensive studies have demonstrated that the NOTCH pathway is altered in a bimodal manner in head and neck squamous cell carcinoma (HNSCC). In a previous study, it was found that the NOTCH4/HEY1 pathway was specifically upregulated in HNSCC and promoted epithelial-mesenchymal transition (EMT), and that HEY1 activation supported SOX2 expression. However, the interactions in this pathway have not yet been fully elucidated. The present study investigated the NOTCH4/HEY1/SOX2 axis in HNSCC using in vitro models and the Cancer Genome Atlas (TCGA) database. To explore the association, reporter and ChIP RT-qPCR assays using SOX2-overexpressing (SOX2-OE) cells were performed. The association between NOTCH4 and HEY1 was examined in the same manner using HEY1-overexpressing (HEY1-OE) cells. The results of the in vitro experiments indicated that HEY1 promoted EMT in the HNSCC cells. Furthermore, the overexpression of HEY1 also promoted sphere formation and increased murine xenograft tumorigenicity. Reporter assays and ChIP RT-qPCR experiments indicated that SOX2 regulated HEY1 expression via direct binding of the HEY1 promoter. HEY1 expression significantly correlated with SOX2 expression in primary lung SCC and other SCCs using the TCGA database. HEY1 also regulated NOTCH4 expression to create a positive reciprocal feedback loop. On the whole, the present study demonstrates that HEY1 expression in HNSCC is regulated via the promotion of SOX2 and promotes EMT. The NOTCH4/HEY1 pathway is specifically upregulated via a positive reciprocal feedback loop mediated by the HEY1-medaited regulation of NOTCH4 transcription, and SOX2 correlates with HEY1 expression in SCC from other primary sites.

ChIPSummitDB: a ChIP-seq-based database of human transcription factor binding sites and the topological arrangements of the proteins bound to them

ChIP-seq reveals genomic regions where proteins, e.g. transcription factors (TFs) interact with DNA. A substantial fraction of these regions, however, do not contain the cognate binding site for the TF of interest. This phenomenon might be explained by protein–protein interactions and co-precipitation of interacting gene regulatory elements. We uniformly processed 3727 human ChIP-seq data sets and determined the cistrome of 292 TFs, as well as the distances between the TF binding motif centers and the ChIP-seq peak summits. ChIPSummitDB enables the analysis of ChIP-seq data using multiple approaches. The 292 cistromes and corresponding ChIP-seq peak sets can be browsed in GenomeView. Overlapping SNPs can be inspected in dbSNPView. Most importantly, the MotifView and PairShiftView pages show the average distance between motif centers and overlapping ChIP-seq peak summits and distance distributions thereof, respectively. In addition to providing a comprehensive human TF binding site collection, the ChIPSummitDB database and web interface allows for the examination of the topological arrangement of TF complexes genome-wide. ChIPSummitDB is freely accessible at http://summit.med.unideb.hu/summitdb/. The database will be regularly updated and extended with the newly available human and mouse ChIP-seq data sets.

Betula platyphylla BpHOX2 transcription factor binds to different cis-acting elements and confers osmotic tolerance

The homeodomain-leucine zipper (HD-Zip) proteins play crucial roles in plant developmental and environmental responses. However, how they mediate gene expression to facilitate abiotic stress tolerance remains unknown. In the present study, we characterized BpHOX2 (encoding a HD-Zip I family protein) from birch (Betula platyphylla). BpHOX2 is predominately expressed in mature stems and leaves, but expressed at a low level in apical buds and roots, suggesting that it has tissue-specific characteristics. BpHOX2 expression was highly induced by osmotic and salt, but only slightly induced by abscisic acid. Overexpression of BpHOX2 markedly improved osmotic tolerance, while knockdown of BpHOX2 increased sensitivity to osmotic stress. BpHOX2 could induce the expression of pyrroline-5-carboxylate synthase, peroxidase, and superoxide dismutase genes to improve proline levels and the reactive oxygen species scavenging capability. Chromatin immunoprecipitation sequencing combined with RNA sequencing showed that BpHOX2 could bind to at least four cis-acting elements, including dehydration-responsive element "RCCGAC", Myb-p binding box "CCWACC," and two novel cis-acting elements with the sequences of "AAGAAG" and "TACGTG" (termed HBS1 and HBS2, respectively) to regulate gene expression. Our results suggested that BpHOX2 is a transcription factor that binds to different cis-acting elements to regulate gene expression, ultimately improving osmotic tolerance in birch.

iODA: An integrated tool for analysis of cancer pathway consistency from heterogeneous multi-omics data

The latest advances in the next generation sequencing technology have greatly facilitated the extensive research of genomics and transcriptomics, thereby promoting the decoding of carcinogenesis with unprecedented resolution. Considering the contribution of analyzing high-throughput multi-omics data to the exploration of cancer molecular mechanisms, an integrated tool for heterogeneous multi-omics data analysis (iODA) is proposed for the systems-level interpretation of multi-omics data, i.e., transcriptomic profiles (mRNA or miRNA expression data) and protein-DNA interactions (ChIP-Seq data). Considering the data heterogeneity, iODA can compare six statistical algorithms in differential analysis for the selected sample data and assist users in choosing the globally optimal one for dysfunctional mRNA or miRNA identification. Since molecular signatures are more consistent at the pathway level than at the gene level, the tool is able to enrich the identified dysfunctional molecules onto the KEGG pathways and extracted the consistent items as key components for further pathogenesis investigation. Compared with other tools, iODA is multi-functional for the systematic analysis of different level of omics data, and its analytical power was demonstrated through case studies of single and cross-level prostate cancer omics data. iODA is open source under GNU GPL and can be downloaded from http://www.sysbio.org.cn/iODA.

Wapl repression by Pax5 promotes V gene recombination by Igh loop extrusion

Nuclear processes, such as V(D)J recombination, are orchestrated by the three-dimensional organization of chromosomes at multiple levels, including compartments1 and topologically associated domains (TADs)2,3 consisting of chromatin loops4. TADs are formed by chromatin-loop extrusion5-7, which depends on the loop-extrusion function of the ring-shaped cohesin complex8-12. Conversely, the cohesin-release factor Wapl13,14 restricts loop extension10,15. The generation of a diverse antibody repertoire, providing humoral immunity to pathogens, requires the participation of all V genes in V(D)J recombination16, which depends on contraction of the 2.8-Mb-long immunoglobulin heavy chain (Igh) locus by Pax517,18. However, how Pax5 controls Igh contraction in pro-B cells remains unknown. Here we demonstrate that locus contraction is caused by loop extrusion across the entire Igh locus. Notably, the expression of Wapl is repressed by Pax5 specifically in pro-B and pre-B cells, facilitating extended loop extrusion by increasing the residence time of cohesin on chromatin. Pax5 mediates the transcriptional repression of Wapl through a single Pax5-binding site by recruiting the polycomb repressive complex 2 to induce bivalent chromatin at the Wapl promoter. Reduced Wapl expression causes global alterations in the chromosome architecture, indicating that the potential to recombine all V genes entails structural changes of the entire genome in pro-B cells.

PAD2-Mediated Citrullination Contributes to Efficient Oligodendrocyte Differentiation and Myelination

Citrullination, the deimination of peptidylarginine residues into peptidylcitrulline, has been implicated in the etiology of several diseases. In multiple sclerosis, citrullination is thought to be a major driver of pathology through hypercitrullination and destabilization of myelin. As such, inhibition of citrullination has been suggested as a therapeutic strategy for MS. Here, in contrast, we show that citrullination by peptidylarginine deiminase 2 (PAD2) contributes to normal oligodendrocyte differentiation, myelination, and motor function. We identify several targets for PAD2, including myelin and chromatin-related proteins, implicating PAD2 in epigenomic regulation. Accordingly, we observe that PAD2 inhibition and its knockdown affect chromatin accessibility and prevent the upregulation of oligodendrocyte differentiation genes. Moreover, mice lacking PAD2 display motor dysfunction and a decreased number of myelinated axons in the corpus callosum. We conclude that citrullination contributes to proper oligodendrocyte lineage progression and myelination.

•

PAD2 is increased upon OL differentiation

•

OL differentiation is facilitated by PAD2-mediated chromatin remodeling in myelin genes

•

PAD2 contributes to efficient myelination and motor and cognitive functions

•

Nuclear and myelin proteins interact and are citrullinated by PAD2

Fat Distribution in Women Is Associated With Depot-Specific Transcriptomic Signatures and Chromatin Structure

BACKGROUND

Preferential accumulation of fat in the upper body (apple shape) is associated with higher risk of developing metabolic syndrome relative to lower body fat (pear shape). We previously discovered that chromatin openness partially defined the transcriptome of preadipocytes isolated from abdominal and gluteofemoral fat. However, the molecular mechanisms underlying interindividual variation in body shape are unknown.

METHODS

Adipocyte fraction was isolated from abdominal and gluteofemoral fat biopsies of premenopausal women (age and body mass index matched) segregated initially only by their waist-to-hip ratio. We evaluated transcriptomic and chromatin accessibility using RNA sequencing and assay for transposase-accessible chromatin using sequencing (ATAC-seq) along with key clinical parameters.

RESULTS

Our data showed that higher lower body fat mass was associated with better lipid profile and free fatty acid decrease after glucose administration. Lipid and glucose metabolic pathways genes were expressed at higher levels in gluteofemoral adipocyte fraction in pears, whereas genes associated with inflammation were higher both in abdominal and gluteofemoral apple adipocyte fraction. Gluteofemoral adipocyte chromatin from pear-shaped women contained a significantly higher number of differentially open ATAC-seq peaks relative to chromatin from the apple-shaped gluteofemoral adipocytes. In contrast, abdominal adipocyte chromatin openness showed few differences between apple- and pear-shaped women. We revealed a correlation between gene transcription and open chromatin at the proximity of the transcriptional start site of some of the differentially expressed genes.

CONCLUSIONS

Integration of data from all 3 approaches suggests that chromatin openness partially governs the transcriptome of gluteofemoral adipocytes and may be involved in the early metabolic syndrome predisposition associated with body shape.

Methylation of a CGATA element inhibits binding and regulation by GATA-1

Alterations in DNA methylation occur during development, but the mechanisms by which they influence gene expression remain uncertain. There are few examples where modification of a single CpG dinucleotide directly affects transcription factor binding and regulation of a target gene in vivo. Here, we show that the erythroid transcription factor GATA-1 — that typically binds T/AGATA sites — can also recognise CGATA elements, but only if the CpG dinucleotide is unmethylated. We focus on a single CGATA site in the c-Kit gene which progressively becomes unmethylated during haematopoiesis. We observe that methylation attenuates GATA-1 binding and gene regulation in cell lines. In mice, converting the CGATA element to a TGATA site that cannot be methylated leads to accumulation of megakaryocyte-erythroid progenitors. Thus, the CpG dinucleotide is essential for normal erythropoiesis and this study illustrates how a single methylated CpG can directly affect transcription factor binding and cellular regulation.

While DNA methylation is thought to play a regulatory role, there are few examples where modification of a single CpG dinucleotide directly affects transcription factor binding. Here the authors show that methylation of a single CGATA element within the c-Kit gene inhibits binding and regulation by erythroid transcription factor GATA-1, both in cells and in mice, suggesting that methylation at this site plays an essential role in erythropoiesis.

Drought-inducible changes in the histone modification H3K9ac are associated with drought-responsive gene expression in Brachypodium distachyon

H3K9ac, an epigenetic marker, is widely distributed in plant genomes. H3K9ac enhances gene expression, which is highly conserved in eukaryotes. However, genome-wide studies of H3K9ac in monocot species are limited, and the changes in H3K9ac under drought stress for individual genes are still not clear. We analysed changes in the H3K9ac level of Brachypodium distachyon under 20% PEG-6000-simulated drought stress conditions. We also performed chromatin immunoprecipitation, followed by next generation sequencing (ChIP-seq) on H3K9ac to reveal changes in H3K9ac for individual genes at the genome-wide level. Our study showed that H3K9ac was mainly enriched in gene exon regions. Drought increased or decreased the H3K9ac level at specific genomic loci. We identified 40 genes associated with increased H3K9ac levels and 36 genes associated with decreased H3K9ac levels under drought stress. Further, RT-qPCR analyses showed that H3K9ac was positively associated with gene expression of those drought-responsive genes. We conclude that H3K9ac enhances the expression level of a large number of drought-responsive genes under drought stress in B. distachyon. The data presented here will help to reveal the correlation of some specific drought-responsive genes and their enriched H3K9ac levels in the model plant B. distachyon.

Alternatively Constructed Estrogen Receptor Alpha-Driven Super-Enhancers Result in Similar Gene Expression in Breast and Endometrial Cell Lines

Super-enhancers (SEs) are clusters of highly active enhancers, regulating cell type-specific and disease-related genes, including oncogenes. The individual regulatory regions within SEs might be simultaneously bound by different transcription factors (TFs) and co-regulators, which together establish a chromatin environment conducting to effective transcription. While cells with distinct TF profiles can have different functions, how different cells control overlapping genetic programs remains a question. In this paper, we show that the construction of estrogen receptor alpha-driven SEs is tissue-specific, both collaborating TFs and the active SE components greatly differ between human breast cancer-derived MCF-7 and endometrial cancer-derived Ishikawa cells; nonetheless, SEs common to both cell lines have similar transcriptional outputs. These results delineate that despite the existence of a combinatorial code allowing alternative SE construction, a single master regulator might be able to determine the overall activity of SEs.

ATM and PRDM9 regulate SPO11-bound recombination intermediates during meiosis

Meiotic recombination is initiated by SPO11-induced double-strand breaks (DSBs). In most mammals, the methyltransferase PRDM9 guides SPO11 targeting, and the ATM kinase controls meiotic DSB numbers. Following MRE11 nuclease removal of SPO11, the DSB is resected and loaded with DMC1 filaments for homolog invasion. Here, we demonstrate the direct detection of meiotic DSBs and resection using END-seq on mouse spermatocytes with low sample input. We find that DMC1 limits both minimum and maximum resection lengths, whereas 53BP1, BRCA1 and EXO1 play surprisingly minimal roles. Through enzymatic modifications to END-seq, we identify a SPO11-bound meiotic recombination intermediate (SPO11-RI) present at all hotspots. We propose that SPO11-RI forms because chromatin-bound PRDM9 asymmetrically blocks MRE11 from releasing SPO11. In Atm-/- spermatocytes, trapped SPO11 cleavage complexes accumulate due to defective MRE11 initiation of resection. Thus, in addition to governing SPO11 breakage, ATM and PRDM9 are critical local regulators of mammalian SPO11 processing.

BedSect: An Integrated Web Server Application to Perform Intersection, Visualization, and Functional Annotation of Genomic Regions From Multiple Datasets

A large number of genomic regions, such as transcription factor binding sites (TFBSs) captured from next generation sequencing (NGS) data analyses or those available from the public resource database ENCODE, are generally overlapped to answer a variety of biological questions. Though several command-line tools are available to perform such an analysis, there is a notable lack of an integrated webserver application with which to identify genomic region intersections, generate publication-ready plots depicting subsets of the overlapped regions, and perform functional annotation. Thus, there is an ardent need for a comprehensive and user-friendly webserver application that allows the users to either upload multiple datasets or select from the integrated Gene Transcription Regulation Database (GTRD). We thus introduce BedSect (http://imgsb.org/bedsect/.), which not only fulfils the above criteria but also performs intersection analysis along with visualization of the intersection regions as an UpSet and correlation plot using the integrated Shiny application. Moreover, analyses, including functional annotation, gene ontology, and biological pathways enrichment for the identified unique and intersected genomic regions, can also be performed using the integrated GREAT tool. To view the genomic regions in the genome browser, the inbuilt hyperlink for UCSC can redirect the user to visualize the results as custom tracks.

ZmMYC2 exhibits diverse functions and enhances JA signaling in transgenic Arabidopsis

ZmMYC2 was identified as the key regulator of JA signaling in maize and exhibited diverse functions through binding to many gene promoters as well as enhanced JA signaling in transgenic Arabidopsis. The plant hormone jasmonate (JA) extensively coordinates plant growth, development and defensive responses. MYC2 is the master regulator of JA signaling and has been widely studied in many plant species. However, little is known about this transcription factor in maize. Here, we identified one maize transcription factor with amino acid identity of 47% to the well-studied Arabidopsis AtMYC2, named as ZmMYC2. Gene expression analysis demonstrated inducible expression patterns of ZmMYC2 in response to multiple plant hormone treatments, as well as biotic and abiotic stresses. The yeast two-hybrid assay indicated physical interaction among ZmMYC2 and JA signal repressors ZmJAZ14, ZmJAZ17, AtJAZ1 and AtJAZ9. ZmMYC2 overexpression in Arabidopsis myc2myc3myc4 restored the sensitivity to JA treatment, resulting in shorter root growth and inducible anthocyanin accumulation. Furthermore, overexpression of ZmMYC2 in Arabidopsis elevated resistance to Botrytis cinerea. Further ChIP-Seq analysis revealed diverse regulatory roles of ZmMYC2 in maize, especially in the signaling crosstalk between JA and auxin. Hence, we identified ZmMYC2 and characterized its roles in regulating JA-mediated growth, development and defense responses.

RNA polymerase II-independent recruitment of SPT6L at transcription start sites in Arabidopsis

SPT6 is a conserved elongation factor that is associated with phosphorylated RNA polymerase II (RNAPII) during transcription. Recent transcriptome analysis in yeast mutants revealed its potential role in the control of transcription initiation at genic promoters. However, the mechanism by which this is achieved and how this is linked to elongation remains to be elucidated. Here, we present the genome-wide occupancy of Arabidopsis SPT6-like (SPT6L) and demonstrate its conserved role in facilitating RNAPII occupancy across transcribed genes. We also further demonstrate that SPT6L enrichment is unexpectedly shifted, from gene body to transcription start site (TSS), when its association with RNAPII is disrupted. Protein domains, required for proper function and enrichment of SPT6L on chromatin, are subsequently identified. Finally, our results suggest that recruitment of SPT6L at TSS is indispensable for its spreading along the gene body during transcription. These findings provide new insights into the mechanisms underlying SPT6L recruitment in transcription and shed light on the coordination between transcription initiation and elongation.