ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia

Non-canonical IKB kinases regulate YAP/TAZ and pathological vascular remodeling behaviors in pulmonary artery smooth muscle cells

Pulmonary arterial hypertension (PAH) causes pulmonary vascular remodeling, increasing pulmonary vascular resistance (PVR) and leading to right heart failure and death. Matrix stiffening early in the disease promotes remodeling in pulmonary artery smooth muscle cells (PASMCs), contributing to PAH pathogenesis. Our research identified YAP and TAZ as key drivers of the mechanobiological feedback loop in PASMCs, suggesting targeting them could mitigate remodeling. However, YAP/TAZ are ubiquitously expressed and carry out diverse functions, necessitating a cell-specific approach. Our previous work demonstrated that targeting non-canonical IKB kinase TBK1 reduced YAP/TAZ activation in human lung fibroblasts. Here, we investigate non-canonical IKB kinases TBK1 and IKKε in pulmonary hypertension (PH) and their potential to modulate PASMC pathogenic remodeling by regulating YAP/TAZ. We show that TBK1 and IKKε are activated in PASMCs in a rat PH model. Inflammatory cytokines, elevated in PAH, activate these kinases in human PASMCs. Inhibiting TBK1/IKKε expression/activity significantly reduces PAH-associated PASMC remodeling, with longer-lasting effects on YAP/TAZ than treprostinil, an approved PAH therapy. These results show that non-canonical IKB kinases regulate YAP/TAZ in PASMCs and may offer a novel approach for reducing vascular remodeling in PAH.

OsJMJ718, a histone demethylase gene, positively regulates seed germination in rice

Seed vigor has major impact on the rate and uniformity of seedling growth, crop yield, and quality. However, the epigenetic regulatory mechanism of crop seed vigor remains unclear. In this study, a (jumonji C) JmjC gene of the histone lysine demethylase OsJMJ718 was cloned in rice, and its roles in seed germination and its epigenetic regulation mechanism were investigated. OsJMJ718 was located in the nucleus and was engaged in H3K9 methylation. Histochemical GUS staining analysis revealed OsJMJ718 was highly expressed in seed embryos. Abiotic stress strongly induced the OsJMJ718 transcriptional accumulation level. Germination percentage and seedling vigor index of OsJMJ718 knockout lines (OsJMJ718-CR) were lower than those of the wild type (WT). Chromatin immunoprecipitation followed by sequencing (ChIP-seq) of seeds imbibed for 24 h showed an increase in H3K9me3 deposition of thousands of genes in OsJMJ718-CR. ChIP-seq results and transcriptome analysis showed that differentially expressed genes were enriched in ABA and ethylene signal transduction pathways. The content of ABA in OsJMJ718-CR was higher than that in WT seeds. OsJMJ718 overexpression enhanced sensitivity to ABA during germination and early seedling growth. In the seed imbibition stage, ABA and ethylene content diminished and augmented, separately, suggesting that OsJMJ718 may adjust rice seed germination through the ABA and ethylene signal transduction pathways. This study displayed the important function of OsJMJ718 in adjusting rice seed germination and vigor, which will provide an essential reference for practical issues, such as improving rice vigor and promoting direct rice sowing production.

H3K27ac modification and transcription characteristics of adipose and muscle tissues in Chuxiang Black pig

The establishment of high-quality pork breeds for improving meat quality in the pig industry is needed. The Chuxiang Black (CX) pig is a new breed developed from Chinese local pigs and Western lean pigs that has a high proportion of lean meat and excellent meat quality. However, the characteristics of cis-regulatory elements in CX pigs are still unknown. In this study, cis-regulatory elements of muscle and adipose tissues in CX pigs were investigated using ChIP-seq and RNA sequencing. Compared with the reported cis-regulatory elements of muscle and adipose tissues, 1768 and 1012 highly activated enhancers and 433 and 275 highly activated promoters in CX muscle and adipose tissues were identified, respectively. Motif analysis showed that transcription factors, such as MEF2A and MEF2C, were core regulators of highly activated enhancers and promoters in muscle. Similarly, the transcription factors JUNB and CUX1 were identified as essential for highly activated enhancers and promoters in CX adipose tissue. These results enrich the resources for the analysis of cis-regulatory elements in the pig genome and provide new basic data for further meat quality improvement through breeding in pigs.

A preliminary investigation of epigenome-wide DNA methylation and temperament during infancy

This study provides preliminary evidence for an epigenetic architecture of infant temperament. At 12 months of age, blood was collected and assayed for DNA methylation and maternally reported infant temperament was assessed using the Infant Behavior Questionnaire in 67 mother-infant dyads. Epigenome-wide analyses showed that the higher order temperament dimensions Surgency and Negative Affect were associated with DNA methylation. The epigenetic signatures of Surgency and Negative Affect were situated at genes involved in synaptic signaling and plasticity. Although replication is required, these results are consistent with a biologically based model of temperament, create new avenues for hypothesis-driven research into epigenetic pathways that underlie individual differences in temperament, and demonstrate that infant temperament has a widespread epigenetic signature in the methylome.

IKZF3/Aiolos H195Y mutation identified in a mouse model of B cell leukemia results in altered DNA binding and altered STAT5-dependent gene expression

Precursor B cell acute lymphoblastic leukemia (Pre-B-ALL) arises from developing B cells and frequently involves mutations in genes encoding transcription factors. In this study, we investigated the function of mutations in the transcription factor IKZF3 (Aiolos), R137* and H195Y, discovered in a mouse model of pre-B-ALL. R137* IKZF3 mutation resulted in a truncated protein, while electrophoretic mobility shift assay showed that H195Y IKZF3 mutation resulted in a protein with altered DNA binding. 38B9 pre-B cell lines were generated expressing WT and H195Y IKZF3 proteins. Anti-IKZF3 ChIP-seq showed that H195Y IKZF3 interacted with a larger number of sites that were different than WT IKZF3. Treatment with interleukin-7 induced changes in gene expression in 38B9 cells expressing WT IKZF3, but did not induce any changes in gene expression in cells expressing H195Y IKZF3. Anti-STAT5 ChIP-seq showed that expression of H195Y IKZF3 resulted in redistribution of STAT5 binding sites in the genome. H195Y IKZF3 binding sites overlapped with a subset of STAT5 binding sites, including in the promoter of the Cish gene. These findings suggest that H195Y mutation of IKZF3 results in altered DNA binding specificity and altered binding of STAT5 to target genes.

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.

FOXL2 interaction with different binding partners regulates the dynamics of ovarian development

The transcription factor FOXL2 is required in ovarian somatic cells for female fertility. Differential timing of Foxl2 deletion, in embryonic versus adult mouse ovary, leads to distinctive outcomes, suggesting different roles across development. Here, we comprehensively investigated FOXL2's role through a multi-omics approach to characterize gene expression dynamics and chromatin accessibility changes, coupled with genome-wide identification of FOXL2 targets and on-chromatin interacting partners in somatic cells across ovarian development. We found that FOXL2 regulates more targets postnatally, through interaction with factors regulating primordial follicle formation and steroidogenesis. Deletion of one interactor, ubiquitin-specific protease 7 (Usp7), results in impairment of somatic cell differentiation, germ cell nest breakdown, and ovarian development, leading to sterility. Our datasets constitute a comprehensive resource for exploration of the molecular mechanisms of ovarian development and causes of female infertility.

Regional centromere configuration in the fungal pathogens of the Pneumocystis genus

Centromeres are constricted chromosomal regions that are essential for cell division. In eukaryotes, centromeres display a remarkable architectural and genetic diversity. The basis of centromere-accelerated evolution remains elusive. Here, we focused on Pneumocystis species, a group of mammalian-specific fungal pathogens that form a sister taxon with that of the Schizosaccharomyces pombe, an important genetic model for centromere biology research. Methods allowing reliable continuous culture of Pneumocystis species do not currently exist, precluding genetic manipulation. CENP-A, a variant of histone H3, is the epigenetic marker that defines centromeres in most eukaryotes. Using heterologous complementation, we show that the Pneumocystis CENP-A ortholog is functionally equivalent to CENP-ACnp1 of S. pombe. Using organisms from a short-term in vitro culture or infected animal models and chromatin immunoprecipitation (ChIP)-Seq, we identified CENP-A bound regions in two Pneumocystis species that diverged ~35 million years ago. Each species has a unique short regional centromere (<10 kb) flanked by heterochromatin in 16-17 monocentric chromosomes. They span active genes and lack conserved DNA sequence motifs and repeats. These features suggest an epigenetic specification of centromere function. Analysis of centromeric DNA across multiple Pneumocystis species suggests a vertical transmission at least 100 million years ago. The common ancestry of Pneumocystis and S. pombe centromeres is untraceable at the DNA level, but the overall architectural similarity could be the result of functional constraint for successful chromosomal segregation.IMPORTANCEPneumocystis species offer a suitable genetic system to study centromere evolution in pathogens because of their phylogenetic proximity with the non-pathogenic yeast S. pombe, a popular model for cell biology. We used this system to explore how centromeres have evolved after the divergence of the two clades ~ 460 million years ago. To address this question, we established a protocol combining short-term culture and ChIP-Seq to characterize centromeres in multiple Pneumocystis species. We show that Pneumocystis have short epigenetic centromeres that function differently from those in S. pombe.

Chromatin profiling reveals TFAP4 as a critical transcriptional regulator of bovine satellite cell differentiation

BACKGROUND

Satellite cells are myogenic precursor cells in adult skeletal muscle and play a crucial role in skeletal muscle regeneration, maintenance, and growth. Like embryonic myoblasts, satellite cells have the ability to proliferate, differentiate, and fuse to form multinucleated myofibers. In this study, we aimed to identify additional transcription factors that control gene expression during bovine satellite cell proliferation and differentiation.

RESULTS

Using chromatin immunoprecipitation followed by sequencing, we identified 56,973 and 54,470 genomic regions marked with both the histone modifications H3K4me1 and H3K27ac, which were considered active enhancers, and 50,956 and 59,174 genomic regions marked with H3K27me3, which were considered repressed enhancers, in proliferating and differentiating bovine satellite cells, respectively. In addition, we identified 1,216 and 1,171 super-enhancers in proliferating and differentiating bovine satellite cells, respectively. Analyzing these enhancers showed that in proliferating bovine satellite cells, active enhancers were associated with genes stimulating cell proliferation or inhibiting myoblast differentiation whereas repressed enhancers were associated with genes essential for myoblast differentiation, and that in differentiating satellite cells, active enhancers were associated with genes essential for myoblast differentiation or muscle contraction whereas repressed enhancers were associated with genes stimulating cell proliferation or inhibiting myoblast differentiation. Active enhancers in proliferating bovine satellite cells were enriched with binding sites for many transcription factors such as MYF5 and the AP-1 family transcription factors; active enhancers in differentiating bovine satellite cells were enriched with binding sites for many transcription factors such as MYOG and TFAP4; and repressed enhancers in both proliferating and differentiating bovine satellite cells were enriched with binding sites for NF-kB, ZEB-1, and several other transcription factors. The role of TFAP4 in satellite cell or myoblast differentiation was previously unknown, and through gene knockdown and overexpression, we experimentally validated a critical role for TFAP4 in the differentiation and fusion of bovine satellite cells into myofibers.

CONCLUSIONS

Satellite cell proliferation and differentiation are controlled by many transcription factors such as AP-1, TFAP4, NF-kB, and ZEB-1 whose roles in these processes were previously unknown in addition to those transcription factors such as MYF5 and MYOG whose roles in these processes are widely known.

Interleukin-22 Contributes to Blood-Brain Barrier Disruption via STAT3/VEGFA Activation in Escherichia coli Meningitis

Escherichia coli continues to be the predominant Gram-negative pathogen causing neonatal meningitis worldwide. Inflammatory mediators have been implicated in the pathogenesis of meningitis and are key therapeutic targets. The role of interleukin-22 (IL-22) in various diseases is diverse, with both protective and pathogenic effects. However, little is understood about the mechanisms underlying the damaging effects of IL-22 on the blood-brain barrier (BBB) in E. coli meningitis. We observed that meningitic E. coli infection induced IL-22 expression in the serum and brain of mice. The tight junction proteins (TJPs) components ZO-1, Occludin, and Claudin-5 were degraded in the mouse brain and human brain microvascular endothelial cells (hBMEC) following IL-22 administration. Moreover, the meningitic E. coli-caused increase in BBB permeability in wild-type mice was restored by knocking out IL-22. Mechanistically, IL-22 activated the STAT3-VEGFA signaling cascade in E. coli meningitis, thus eliciting the degradation of TJPs to induce BBB disruption. Our data indicated that IL-22 is an essential host accomplice during E. coli-caused BBB disruption and could be targeted for the therapy of bacterial meningitis.

Acetyl-CoA production by Mediator-bound 2-ketoacid dehydrogenases boosts de novo histone acetylation and is regulated by nitric oxide

Histone-modifying enzymes depend on the availability of cofactors, with acetyl-coenzyme A (CoA) being required for histone acetyltransferase (HAT) activity. The discovery that mitochondrial acyl-CoA-producing enzymes translocate to the nucleus suggests that high concentrations of locally synthesized metabolites may impact acylation of histones and other nuclear substrates, thereby controlling gene expression. Here, we show that 2-ketoacid dehydrogenases are stably associated with the Mediator complex, thus providing a local supply of acetyl-CoA and increasing the generation of hyper-acetylated histone tails. Nitric oxide (NO), which is produced in large amounts in lipopolysaccharide-stimulated macrophages, inhibited the activity of Mediator-associated 2-ketoacid dehydrogenases. Elevation of NO levels and the disruption of Mediator complex integrity both affected de novo histone acetylation within a shared set of genomic regions. Our findings indicate that the local supply of acetyl-CoA generated by 2-ketoacid dehydrogenases bound to Mediator is required to maximize acetylation of histone tails at sites of elevated HAT activity.

Deficiency of thioredoxin-interacting protein results in age-related thrombocytopenia due to megakaryocyte oxidative stress

BACKGROUND

Platelets are generated from megakaryocytes (MKs), mainly located in the bone marrow (BM). Megakaryopoiesis can be affected by genetic disorders, metabolic diseases, and aging. The molecular mechanisms underlying platelet count regulation have not been fully elucidated.

OBJECTIVES

In the present study, we investigated the role of thioredoxin-interacting protein (TXNIP), a protein that regulates cellular metabolism in megakaryopoiesis, using a Txnip-/- mouse model.

METHODS

Wild-type (WT) and Txnip-/- mice (2-27-month-old) were studied. BM-derived MKs were analyzed to investigate the role of TXNIP in megakaryopoiesis with age. The global transcriptome of BM-derived CD41+ megakaryocyte precursors (MkPs) of WT and Txnip-/- mice were compared. The CD34+ hematopoietic stem cells isolated from human cord blood were differentiated into MKs.

RESULTS

Txnip-/- mice developed thrombocytopenia at 4 to 5 months that worsened with age. During ex vivo megakaryopoiesis, Txnip-/- MkPs remained small, with decreased levels of MK-specific markers. Critically, Txnip-/- MkPs exhibited reduced mitochondrial reactive oxygen species, which was related to AKT activity. Txnip-/- MkPs also showed elevated glycolysis alongside increased glucose uptake for ATP production. Total RNA sequencing revealed enrichment for oxidative stress- and apoptosis-related genes in differentially expressed genes between Txnip-/- and WT MkPs. The effects of TXNIP on MKs were recapitulated during the differentiation of human cord blood-derived CD34+ hematopoietic stem cells.

CONCLUSION

We provide evidence that the megakaryopoiesis pathway becomes exhausted with age in Txnip-/- mice with a decrease in terminal, mature MKs that response to thrombocytopenic challenge. Overall, this study demonstrates the role of TXNIP in megakaryopoiesis, regulating mitochondrial metabolism.

Hypertensive Pressure Mechanosensing Alone Triggers Lipid Droplet Accumulation and Transdifferentiation of Vascular Smooth Muscle Cells to Foam Cells

Arterial Vascular smooth muscle cells (VSMCs) play a central role in the onset and progression of atherosclerosis. Upon exposure to pathological stimuli, they can take on alternative phenotypes that, among others, have been described as macrophage like, or foam cells. VSMC foam cells make up >50% of all arterial foam cells and have been suggested to retain an even higher proportion of the cell stored lipid droplets, further leading to apoptosis, secondary necrosis, and an inflammatory response. However, the mechanism of VSMC foam cell formation is still unclear. Here, it is identified that mechanical stimulation through hypertensive pressure alone is sufficient for the phenotypic switch. Hyperspectral stimulated Raman scattering imaging demonstrates rapid lipid droplet formation and changes to lipid metabolism and changes are confirmed in ABCA1, KLF4, LDLR, and CD68 expression, cell proliferation, and migration. Further, a mechanosignaling route is identified involving Piezo1, phospholipid, and arachidonic acid signaling, as well as epigenetic regulation, whereby CUT&Tag epigenomic analysis confirms changes in the cells (lipid) metabolism and atherosclerotic pathways. Overall, the results show for the first time that VSMC foam cell formation can be triggered by mechanical stimulation alone, suggesting modulation of mechanosignaling can be harnessed as potential therapeutic strategy.

PHF6-mediated transcriptional control of NSC via Ephrin receptors is impaired in the intellectual disability syndrome BFLS

The plant homeodomain zinc-finger protein, PHF6, is a transcriptional regulator, and PHF6 germline mutations cause the X-linked intellectual disability (XLID) Börjeson-Forssman-Lehmann syndrome (BFLS). The mechanisms by which PHF6 regulates transcription and how its mutations cause BFLS remain poorly characterized. Here, we show genome-wide binding of PHF6 in the developing cortex in the vicinity of genes involved in central nervous system development and neurogenesis. Characterization of BFLS mice harbouring PHF6 patient mutations reveals an increase in embryonic neural stem cell (eNSC) self-renewal and a reduction of neural progenitors. We identify a panel of Ephrin receptors (EphRs) as direct transcriptional targets of PHF6. Mechanistically, we show that PHF6 regulation of EphR is impaired in BFLS mice and in conditional Phf6 knock-out mice. Knockdown of EphR-A phenocopies the PHF6 loss-of-function defects in altering eNSCs, and its forced expression rescues defects of BFLS mice-derived eNSCs. Our data indicate that PHF6 directly promotes Ephrin receptor expression to control eNSC behaviour in the developing brain, and that this pathway is impaired in BFLS.

Regulatory activity is the default DNA state in eukaryotes

Genomes encode for genes and non-coding DNA, both capable of transcriptional activity. However, unlike canonical genes, many transcripts from non-coding DNA have limited evidence of conservation or function. Here, to determine how much biological noise is expected from non-genic sequences, we quantify the regulatory activity of evolutionarily naive DNA using RNA-seq in yeast and computational predictions in humans. In yeast, more than 99% of naive DNA bases were transcribed. Unlike the evolved transcriptome, naive transcripts frequently overlapped with opposite sense transcripts, suggesting selection favored coherent gene structures in the yeast genome. In humans, regulation-associated chromatin activity is predicted to be common in naive dinucleotide-content-matched randomized DNA. Here, naive and evolved DNA have similar co-occurrence and cell-type specificity of chromatin marks, challenging these as indicators of selection. However, in both yeast and humans, extreme high activities were rare in naive DNA, suggesting they result from selection. Overall, basal regulatory activity seems to be the default, which selection can hone to evolve a function or, if detrimental, repress.

PIBF1 regulates multiple gene expression via impeding long-range chromatin interaction to drive the malignant transformation of HPV16 integration epithelial cells

INTRODUCTION

Human papillomavirus (HPV) integration can induce gene expression dysregulation by destroying higher-order chromatin structure in cervical cancer.

OBJECTIVES

We established a 13q22 site-specific HPV16 gene knock-in cell model to interrogate the changes in chromatin structure at the initial stages of host cell malignant transformation.

METHODS

We designed a CRISPR-Cas9 system with sgRNA targeting 13q22 site and constructed the HPV16 gene donor. Cells were cotransfected, screened, and fluorescence sorted. The whole genome sequencing (WGS) was used to confirm the precise HPV16 gene integration site. Western blot and qRT-PCR were used to measure gene expression. In vitro and in vivo analysis were performed to estimate the tumorigenic potential of the HPV16 knock-in cell model. Combined Hi-C, chromatin immunoprecipitation and RNA sequencing analyses revealed correlations between chromatin structure and gene expression. We performed a coimmunoprecipitation assay with anti-PIBF1 antibody to identify endogenous interacting proteins. In vivo analysis was used to determine the role of PIBF1 in the tumor growth of cervical cancer cells.

RESULTS

We successfully established a 13q22 site-specific HPV16 gene knock-in cell model. We found that HPV integration promoted cell proliferation, invasion and stratified growth in vitro, and monoclonal proliferation in vivo. HPV integration divided the affected topologically associated domain (TAD) into two smaller domains, and the progesterone-induced blocking factor 1 (PIBF1) gene near the integration site was upregulated, although PIBF1 was not enriched at the domain boundary by CUT-Tag signal analysis. Moreover, PIBF1 was found to interact with the cohesin complex off chromatin to reduce contact domain formation by disrupting the cohesin ring-shaped structure, causing dysregulation of tumorigenesis-related genes. Xenograft experiments determined the role of PIBF1 in the proliferation in cervical cancer cells.

CONCLUSION

We highlight that PIBF1, a potential chromatin structure regulatory protein, is activated by HPV integration, which provides new insights into HPV integration-driven cervical carcinogenesis.

Mapping catalytically engaged TOP2B in neurons reveals the principles of topoisomerase action within the genome

We trapped catalytically engaged topoisomerase IIβ (TOP2B) in covalent DNA cleavage complexes (TOP2Bccs) and mapped their positions genome-wide in cultured mouse cortical neurons. We report that TOP2Bcc distribution varies with both nucleosome and compartmental chromosome organization. While TOP2Bccs in gene bodies correlate with their level of transcription, highly expressed genes that lack the usually associated chromatin marks, such as H3K36me3, show reduced TOP2Bccs, suggesting that histone posttranslational modifications regulate TOP2B activity. Promoters with high RNA polymerase II occupancy show elevated TOP2B chromatin immunoprecipitation sequencing signals but low TOP2Bccs, indicating that TOP2B catalytic engagement is curtailed at active promoters. Surprisingly, either poisoning or inhibiting TOP2B increases nascent transcription at most genes and enhancers but reduces transcription within long genes. These effects are independent of transcript length and instead correlate with the presence of intragenic enhancers. Together, these results clarify how cells modulate the catalytic engagement of topoisomerases to affect transcription.

Light-sensitive short hypocotyl genes confer symbiotic nodule identity in the legume Medicago truncatula

Legumes produce specialized root nodules that are distinct from lateral roots in morphology and function, with nodules intracellularly hosting nitrogen-fixing bacteria. We have previously shown that a lateral root program underpins nodule initiation, but there must be additional developmental regulators that confer nodule identity. Here, we show two members of the LIGHT-SENSITIVE SHORT HYPOCOTYL (LSH) transcription factor family, predominantly known to define shoot meristem complexity and organ boundaries, function as regulators of nodule organ identity. In parallel to the root initiation program, LSH1/LSH2 recruit a program into the root cortex that mediates the divergence into nodules, in particular with cell divisions in the mid-cortex. This includes regulation of auxin and cytokinin, promotion of NODULE ROOT1/2 and Nuclear Factor YA1, and suppression of the lateral root program. A principal outcome of LSH1/LSH2 function is the production of cells able to accommodate nitrogen-fixing bacteria, a key feature unique to nodules.

DNA methylation-environment interactions in the human genome

Previously, we showed that a massively parallel reporter assay, mSTARR-seq, could be used to simultaneously test for both enhancer-like activity and DNA methylation-dependent enhancer activity for millions of loci in a single experiment (Lea et al., 2018). Here, we apply mSTARR-seq to query nearly the entire human genome, including almost all CpG sites profiled either on the commonly used Illumina Infinium MethylationEPIC array or via reduced representation bisulfite sequencing. We show that fragments containing these sites are enriched for regulatory capacity, and that methylation-dependent regulatory activity is in turn sensitive to the cellular environment. In particular, regulatory responses to interferon alpha (IFNA) stimulation are strongly attenuated by methyl marks, indicating widespread DNA methylation-environment interactions. In agreement, methylation-dependent responses to IFNA identified via mSTARR-seq predict methylation-dependent transcriptional responses to challenge with influenza virus in human macrophages. Our observations support the idea that pre-existing DNA methylation patterns can influence the response to subsequent environmental exposures-one of the tenets of biological embedding. However, we also find that, on average, sites previously associated with early life adversity are not more likely to functionally influence gene regulation than expected by chance.

Quantitative transcriptomic and epigenomic data analysis: a primer

The advent of microarray and second generation sequencing technology has revolutionized the field of molecular biology, allowing researchers to quantitatively assess transcriptomic and epigenomic features in a comprehensive and cost-efficient manner. Moreover, technical advancements have pushed the resolution of these sequencing techniques to the single cell level. As a result, the bottleneck of molecular biology research has shifted from the bench to the subsequent omics data analysis. Even though most methodologies share the same general strategy, state-of-the-art literature typically focuses on data type specific approaches and already assumes expert knowledge. Here, however, we aim at providing conceptual insight in the principles of genome-wide quantitative transcriptomic and epigenomic (including open chromatin assay) data analysis by describing a generic workflow. By starting from a general framework and its assumptions, the need for alternative or additional data-analytical solutions when working with specific data types becomes clear, and are hence introduced. Thus, we aim to enable readers with basic omics expertise to deepen their conceptual and statistical understanding of general strategies and pitfalls in omics data analysis and to facilitate subsequent progression to more specialized literature.

Multimodal Hox5 activity generates motor neuron diversity

Motor neurons (MNs) are the final output of circuits driving fundamental behaviors, such as respiration and locomotion. Hox proteins are essential in generating the MN diversity required for accomplishing these functions, but the transcriptional mechanisms that enable Hox paralogs to assign distinct MN subtype identities despite their promiscuous DNA binding motif are not well understood. Here we show that Hoxa5 controls chromatin accessibility in all mouse spinal cervical MN subtypes and engages TALE co-factors to directly bind and regulate subtype-specific genes. We identify a paralog-specific interaction of Hoxa5 with the phrenic MN-specific transcription factor Scip and show that heterologous expression of Hoxa5 and Scip is sufficient to suppress limb-innervating MN identity. We also demonstrate that phrenic MN identity is stable after Hoxa5 downregulation and identify Klf proteins as potential regulators of phrenic MN maintenance. Our data identify multiple modes of Hoxa5 action that converge to induce and maintain MN identity.

Positive Selection Drives cis-regulatory Evolution Across the Threespine Stickleback Y Chromosome

Allele-specific gene expression evolves rapidly on heteromorphic sex chromosomes. Over time, the accumulation of mutations on the Y chromosome leads to widespread loss of gametolog expression, relative to the X chromosome. It remains unclear if expression evolution on degrading Y chromosomes is primarily driven by mutations that accumulate through processes of selective interference, or if positive selection can also favor the down-regulation of coding regions on the Y chromosome that contain deleterious mutations. Identifying the relative rates of cis-regulatory sequence evolution across Y chromosomes has been challenging due to the limited number of reference assemblies. The threespine stickleback (Gasterosteus aculeatus) Y chromosome is an excellent model to identify how regulatory mutations accumulate on Y chromosomes due to its intermediate state of divergence from the X chromosome. A large number of Y-linked gametologs still exist across 3 differently aged evolutionary strata to test these hypotheses. We found that putative enhancer regions on the Y chromosome exhibited elevated substitution rates and decreased polymorphism when compared to nonfunctional sites, like intergenic regions and synonymous sites. This suggests that many cis-regulatory regions are under positive selection on the Y chromosome. This divergence was correlated with X-biased gametolog expression, indicating the loss of expression from the Y chromosome may be favored by selection. Our findings provide evidence that Y-linked cis-regulatory regions exhibit signs of positive selection quickly after the suppression of recombination and allow comparisons with recent theoretical models that suggest the rapid divergence of regulatory regions may be favored to mask deleterious mutations on the Y chromosome.

Peak-agnostic high-resolution cis-regulatory circuitry mapping using single cell multiome data

Single same cell RNAseq/ATACseq multiome data provide unparalleled potential to develop high resolution maps of the cell-type specific transcriptional regulatory circuitry underlying gene expression. We present CREMA, a framework that recovers the full cis-regulatory circuitry by modeling gene expression and chromatin activity in individual cells without peak-calling or cell type labeling constraints. We demonstrate that CREMA overcomes the limitations of existing methods that fail to identify about half of functional regulatory elements which are outside the called chromatin 'peaks'. These circuit sites outside called peaks are shown to be important cell type specific functional regulatory loci, sufficient to distinguish individual cell types. Analysis of mouse pituitary data identifies a Gata2-circuit for the gonadotrope-enriched disease-associated Pcsk1 gene, which is experimentally validated by reduced gonadotrope expression in a gonadotrope conditional Gata2-knockout model. We present a web accessible human immune cell regulatory circuit resource, and provide CREMA as an R package.

CRISPR screening identifies BET and mTOR inhibitor synergy in cholangiocarcinoma through serine glycine one carbon

Patients with cholangiocarcinoma have poor clinical outcomes due to late diagnoses, poor prognoses, and limited treatment strategies. To identify drug combinations for this disease, we have conducted a genome-wide CRISPR screen anchored on the bromodomain and extraterminal domain (BET) PROTAC degrader ARV825, from which we identified anticancer synergy when combined with genetic ablation of members of the mTOR pathway. This combination effect was validated using multiple pharmacological BET and mTOR inhibitors, accompanied by increased levels of apoptosis and cell cycle arrest. In a xenograft model, combined BET degradation and mTOR inhibition induced tumor regression. Mechanistically, the 2 inhibitor classes converged on H3K27ac-marked epigenetic suppression of the serine glycine one carbon (SGOC) metabolism pathway, including the key enzymes PHGDH and PSAT1. Knockdown of PSAT1 was sufficient to replicate synergy with single-agent inhibition of either BET or mTOR. Our results tie together epigenetic regulation, metabolism, and apoptosis induction as key therapeutic targets for further exploration in this underserved disease.

Binding profiles for 954 Drosophila and C. elegans transcription factors reveal tissue specific regulatory relationships

A catalog of transcription factor (TF) binding sites in the genome is critical for deciphering regulatory relationships. Here we present the culmination of the modERN (model organism Encyclopedia of Regulatory Networks) consortium that systematically assayed TF binding events in vivo in two major model organisms, Drosophila melanogaster (fly) and Caenorhabditis elegans (worm). We describe key features of these datasets, comprising 604 TFs identifying 3.6M sites in the fly and 350 TFs identifying 0.9 M sites in the worm. Applying a machine learning model to these data identifies sets of TFs with a prominent role in promoting target gene expression in specific cell types. TF binding data are available through the ENCODE Data Coordinating Center and at https://epic.gs.washington.edu/modERNresource, which provides access to processed and summary data, as well as widgets to probe cell type-specific TF-target relationships. These data are a rich resource that should fuel investigations into TF function during development.

Position Weight Matrix or Acyclic Probabilistic Finite Automaton: Which model to use? A decision rule inferred for the prediction of transcription factor binding sites

Prediction of transcription factor binding sites (TFBS) is an example of application of Bioinformatics where DNA molecules are represented as sequences of A, C, G and T symbols. The most used model in this problem is Position Weight Matrix (PWM). Notwithstanding the advantage of being simple, PWMs cannot capture dependency between nucleotide positions, which may affect prediction performance. Acyclic Probabilistic Finite Automata (APFA) is an alternative model able to accommodate position dependencies. However, APFA is a more complex model, which means more parameters have to be learned. In this paper, we propose an innovative method to identify when position dependencies influence preference for PWMs or APFAs. This implied using position dependency features extracted from 1106 sets of TFBS to infer a decision tree able to predict which is the best model - PWM or APFA - for a given set of TFBSs. According to our results, as few as three pinpointed features are able to choose the best model, providing a balance of performance (average precision) and model simplicity.

Egr-1 is a key regulator of the blood-brain barrier damage induced by meningitic Escherichia coli

Bacterial meningitis remains a leading cause of infection-related mortality worldwide. Although Escherichia coli (E. coli) is the most common etiology of neonatal meningitis, the underlying mechanisms governing bacterial blood-brain barrier (BBB) disruption during infection remain elusive. We observed that infection of human brain microvascular endothelial cells with meningitic E. coli triggers the activation of early growth response 1 (Egr-1), a host transcriptional activator. Through integrated chromatin immunoprecipitation sequencing and transcriptome analysis, we identified Egr-1 as a crucial regulator for maintaining BBB integrity. Mechanistically, Egr-1 induced cytoskeletal changes and downregulated tight junction protein expression by directly targeting VEGFA, PDGFB, and ANGPTL4, resulting in increased BBB permeability. Meanwhile, Egr-1 also served as a master regulator in the initiation of neuroinflammatory response during meningitic E. coli infection. Our findings support an Egr-1-dependent mechanism of BBB disruption by meningitic E. coli, highlighting a promising therapeutic target for bacterial meningitis.

Epigenetic modifications of 45S rDNA associates with the disruption of nucleolar organisation during Cd stress response in Pakchoi

The role of the nucleolus in Pakchoi response to Cd stress remains largely unknown. In this work, we focus on exploring the underling mechanism between nucleolus disruption and epigenetic modification in Cd stressed-Pakchoi. Our results indicated that the proportion of nucleolus disruption, decondensation of 45 S rDNA chromatin, and a simultaneous increase in 5' external transcribed spacer region (ETS) transcription were observed with increasing Cd concentration, accompanied by genome-wide alterations in the levels of histone acetylation and methylation. Further results showed that Cd treatment exhibited a significant increase in H3K9ac, H4K5ac, and H3K9me2 levels occurred in promoter regions of the 45 S rDNA. Additionally, DNA methylation assays in the 45 S rDNA promoter region revealed that individual site-specific hypomethylation may be engaged in the activation of 45 S rDNA transcription. Our study provides some molecular mechanisms for the linkage between Cd stress, rDNA epigenetic modifications, and nucleolus disintegration in plants.

hnRNP A1 dysfunction alters RNA splicing and drives neurodegeneration in multiple sclerosis (MS)

Neurodegeneration is the primary driver of disease progression in multiple sclerosis (MS) resulting in permanent disability, creating an urgent need to discover its underlying mechanisms. Herein, we establish that dysfunction of the RNA binding protein heterogeneous nuclear ribonucleoprotein A1 (hnRNP A1) results in differential of binding to RNA targets causing alternative RNA splicing, which contributes to neurodegeneration in MS and its models. Using RNAseq of MS brains, we discovered differential expression and aberrant splicing of hnRNP A1 target RNAs involved in neuronal function and RNA homeostasis. We confirmed this in vivo in experimental autoimmune encephalomyelitis employing CLIPseq specific for hnRNP A1, where hnRNP A1 differentially binds and regulates RNA, including aberrantly spliced targets identified in human samples. Additionally, dysfunctional hnRNP A1 expression in neurons caused neurite loss and identical changes in splicing, corroborating hnRNP A1 dysfunction as a cause of neurodegeneration. Collectively, these data indicate hnRNP A1 dysfunction causes altered neuronal RNA splicing, resulting in neurodegeneration in MS.

Multiomics-integrated deep language model enables in silico genome-wide detection of transcription factor binding site in unexplored biosamples

MOTIVATION

Transcription factor binding sites (TFBS) are regulatory elements that have significant impact on transcription regulation and cell fate determination. Canonical motifs, biological experiments, and computational methods have made it possible to discover TFBS. However, most existing in silico TFBS prediction models are solely DNA-based, and are trained and utilized within the same biosample, which fail to infer TFBS in experimentally unexplored biosamples.

RESULTS

Here, we propose TFBS prediction by modified TransFormer (TFTF), a multimodal deep language architecture which integrates multiomics information in epigenetic studies. In comparison to existing computational techniques, TFTF has state-of-the-art accuracy, and is also the first approach to accurately perform genome-wide detection for cell-type and species-specific TFBS in experimentally unexplored biosamples. Compared to peak calling methods, TFTF consistently discovers true TFBS in threshold tuning-free way, with higher recalled rates. The underlying mechanism of TFTF reveals greater attention to the targeted TF's motif region in TFBS, and general attention to the entire peak region in non-TFBS. TFTF can benefit from the integration of broader and more diverse data for improvement and can be applied to multiple epigenetic scenarios.

AVAILABILITY AND IMPLEMENTATION

We provide a web server (https://tftf.ibreed.cn/) for users to utilize TFTF model. Users can train TFTF model and discover TFBS with their own data.

Efficient Enrichment of Synchronized Mouse Spermatocytes Suitable for Genome-Wide Analysis

Chromatin undergoes extensive remodeling during meiosis, leading to specific patterns of gene expression and chromosome organization, which ultimately controls fundamental meiotic processes such as recombination and homologous chromosome associations. Recent game-changing advances have been made by analysis of chromatin binding sites of meiotic specific proteins genome-wide in mouse spermatocytes. However, further progress is still highly dependent on the reliable isolation of sufficient quantities of spermatocytes at specific stages of prophase I. Here, we describe a combination of methodologies we adapted for rapid and reliable isolation of synchronized fixed mouse spermatocytes. We show that chromatin isolated from these cells can be used to study chromatin-binding sites by ChIP-seq. High-quality data we obtained from INO80 ChIP-seq in zygotene cells was used for functional analysis of chromatin-binding sites.

MINI-AC: inference of plant gene regulatory networks using bulk or single-cell accessible chromatin profiles

Gene regulatory networks (GRNs) represent the interactions between transcription factors (TF) and their target genes. Plant GRNs control transcriptional programs involved in growth, development, and stress responses, ultimately affecting diverse agricultural traits. While recent developments in accessible chromatin (AC) profiling technologies make it possible to identify context-specific regulatory DNA, learning the underlying GRNs remains a major challenge. We developed MINI-AC (Motif-Informed Network Inference based on Accessible Chromatin), a method that combines AC data from bulk or single-cell experiments with TF binding site (TFBS) information to learn GRNs in plants. We benchmarked MINI-AC using bulk AC datasets from different Arabidopsis thaliana tissues and showed that it outperforms other methods to identify correct TFBS. In maize, a crop with a complex genome and abundant distal AC regions, MINI-AC successfully inferred leaf GRNs with experimentally confirmed, both proximal and distal, TF-target gene interactions. Furthermore, we showed that both AC regions and footprints are valid alternatives to infer AC-based GRNs with MINI-AC. Finally, we combined MINI-AC predictions from bulk and single-cell AC datasets to identify general and cell-type specific maize leaf regulators. Focusing on C4 metabolism, we identified diverse regulatory interactions in specialized cell types for this photosynthetic pathway. MINI-AC represents a powerful tool for inferring accurate AC-derived GRNs in plants and identifying known and novel candidate regulators, improving our understanding of gene regulation in plants.

Differential Analysis of Protein-DNA Binding Using ChIP-Seq Data

Chromatin immunoprecipitation in combination with next-generation sequencing (ChIP-Seq) allows probing of protein-DNA binding in a rapid and genome-wide fashion. Herein we describe the required steps to preprocess ChIP-Seq data and to analyze the differential binding of proteins to DNA for perturbation experiments. In these experiments, different conditions are compared to find the underlying biological mechanisms caused by the stimulus or treatment. In addition, we provide a sample analysis using the steps outlined in the chapter.

Enhancing Cell Line Stability by CRISPR/Cas9-Mediated Site-Specific Integration Based on Histone Modifications

In traditional cell line design pipelines, cost- and time-intensive long-term stability studies must be performed due to random integration of the transgene into the genome. By this, integration into epigenetically silenced regions can lead to silencing of the recombinant promoter over time. Site-specific integration into regions with active chromatin structure can overcome this problem and lead to strong and stable gene expression. Here, we describe a detailed protocol to identify integration sites with epigenetically preferable properties by chromatin immunoprecipitation sequencing and use them for stable and strong gene expression by applying CRISPR/Cas9. Furthermore, the examination of the integration sites with focus on Cas9-targeted sequencing with nanopores is described.

Polycomb regulates circadian rhythms in Drosophila in clock neurons

Circadian rhythms are essential physiological feature for most living organisms. Previous studies have shown that epigenetic regulation plays a crucial role. There is a knowledge gap in the chromatin state of some key clock neuron clusters. In this study, we show that circadian rhythm is affected by the epigenetic regulator Polycomb (Pc) within the Drosophila clock neurons. To investigate the molecular mechanisms underlying the roles of Pc in these clock neuron clusters, we use targeted DamID (TaDa) to identify genes significantly bound by Pc in the neurons marked by C929-Gal4 (including l-LNvs cluster), R6-Gal4 (including s-LNvs cluster), R18H11-Gal4 (including DN1 cluster), and DVpdf-Gal4, pdf-Gal80 (including LNds cluster). It shows that Pc binds to the genes involved in the circadian rhythm pathways, arguing a direct role for Pc in regulating circadian rhythms through specific clock genes. This study shows the identification of Pc targets in the clock neuron clusters, providing potential resource for understanding the regulatory mechanisms of circadian rhythms by the PcG complex. Thus, this study provided an example for epigenetic regulation of adult behavior.

DNA methylation-environment interactions in the human genome

Previously we showed that a massively parallel reporter assay, mSTARR-seq, could be used to simultaneously test for both enhancer-like activity and DNA methylation-dependent enhancer activity for millions of loci in a single experiment (Lea et al., 2018). Here we apply mSTARR-seq to query nearly the entire human genome, including almost all CpG sites profiled either on the commonly used Illumina Infinium MethylationEPIC array or via reduced representation bisulfite sequencing. We show that fragments containing these sites are enriched for regulatory capacity, and that methylation-dependent regulatory activity is in turn sensitive to the cellular environment. In particular, regulatory responses to interferon alpha (IFNA) stimulation are strongly attenuated by methyl marks, indicating widespread DNA methylation-environment interactions. In agreement, methylation-dependent responses to IFNA identified via mSTARR-seq predict methylation-dependent transcriptional responses to challenge with influenza virus in human macrophages. Our observations support the idea that pre-existing DNA methylation patterns can influence the response to subsequent environmental exposures-one of the tenets of biological embedding. However, we also find that, on average, sites previously associated with early life adversity are not more likely to functionally influence gene regulation than expected by chance.

Plasma growth hormone pulses induce male-biased pulsatile chromatin opening and epigenetic regulation in adult mouse liver

Sex differences in plasma growth hormone (GH) profiles, pulsatile in males and persistent in females, regulate sex differences in hepatic STAT5 activation linked to sex differences in gene expression and liver disease susceptibility, but little is understood about the fundamental underlying, GH pattern-dependent regulatory mechanisms. Here, DNase-I hypersensitivity site (DHS) analysis of liver chromatin accessibility in a cohort of 18 individual male mice established that the endogenous male rhythm of plasma GH pulse-stimulated liver STAT5 activation induces dynamic, repeated cycles of chromatin opening and closing at several thousand liver DHS and comprises a novel mechanism conferring male bias to liver chromatin accessibility. Strikingly, a single physiological replacement dose of GH given to hypophysectomized male mice restored, within 30 min, liver STAT5 activity and chromatin accessibility at 83% of the dynamic, pituitary hormone-dependent male-biased DHS. Sex-dependent transcription factor binding patterns and chromatin state analysis identified key genomic and epigenetic features distinguishing this dynamic, STAT5-driven mechanism of male-biased chromatin opening from a second GH-dependent mechanism operative at static male-biased DHS, which are constitutively open in male liver. Dynamic but not static male-biased DHS adopt a bivalent-like epigenetic state in female liver, as do static female-biased DHS in male liver, albeit using distinct repressive histone marks in each sex, namely, H3K9me3 at male-biased DHS in female liver and H3K27me3 at female-biased DHS in male liver. Moreover, sex-biased H3K36me3 marks are uniquely enriched at static sex-biased DHS, which may serve to keep these sex-dependent hepatocyte enhancers free of H3K27me3 repressive marks and thus constitutively open. Pulsatile chromatin opening stimulated by endogenous, physiological hormone pulses is thus one of two distinct GH-determined mechanisms for establishing widespread sex differences in hepatic chromatin accessibility and epigenetic regulation, both closely linked to sex-biased gene transcription and the sexual dimorphism of liver function.

Genome-wide identification of transcriptional enhancers during human placental development and association with function, differentiation, and disease†

The placenta is a dynamic organ that must perform a remarkable variety of functions during its relatively short existence in order to support a developing fetus. These functions include nutrient delivery, gas exchange, waste removal, hormone production, and immune barrier protection. Proper placenta development and function are critical for healthy pregnancy outcomes, but the underlying genomic regulatory events that control this process remain largely unknown. We hypothesized that mapping sites of transcriptional enhancer activity and associated changes in gene expression across gestation in human placenta tissue would identify genomic loci and predicted transcription factor activity related to critical placenta functions. We used a suite of genomic assays [i.e., RNA-sequencing (RNA-seq), Precision run-on-sequencing (PRO-seq), and Chromatin immunoprecipitation-sequencing (ChIP-seq)] and computational pipelines to identify a set of >20 000 enhancers that are active at various time points in gestation. Changes in the activity of these enhancers correlate with changes in gene expression. In addition, some of these enhancers encode risk for adverse pregnancy outcomes. We further show that integrating enhancer activity, transcription factor motif analysis, and transcription factor expression can identify distinct sets of transcription factors predicted to be more active either in early pregnancy or at term. Knockdown of selected identified transcription factors in a trophoblast stem cell culture model altered the expression of key placental marker genes. These observations provide a framework for future mechanistic studies of individual enhancer-transcription factor-target gene interactions and have the potential to inform genetic risk prediction for adverse pregnancy outcomes.

GenoPipe: identifying the genotype of origin within (epi)genomic datasets

Confidence in experimental results is critical for discovery. As the scale of data generation in genomics has grown exponentially, experimental error has likely kept pace despite the best efforts of many laboratories. Technical mistakes can and do occur at nearly every stage of a genomics assay (i.e. cell line contamination, reagent swapping, tube mislabelling, etc.) and are often difficult to identify post-execution. However, the DNA sequenced in genomic experiments contains certain markers (e.g. indels) encoded within and can often be ascertained forensically from experimental datasets. We developed the Genotype validation Pipeline (GenoPipe), a suite of heuristic tools that operate together directly on raw and aligned sequencing data from individual high-throughput sequencing experiments to characterize the underlying genome of the source material. We demonstrate how GenoPipe validates and rescues erroneously annotated experiments by identifying unique markers inherent to an organism's genome (i.e. epitope insertions, gene deletions and SNPs).

Ex vivo mapping of enhancer networks that define the transcriptional program driving melanoma metastasis

Mortality from vmelanoma is associated with metastatic disease, but the mechanisms leading to spreading of the cancer cells remain obscure. Spatial profiling revealed that melanoma is characterized by a high degree of heterogeneity, which is established by the ability of melanoma cells to switch between different phenotypical stages. This plasticity, likely a heritage from embryonic pathways, accounts for a relevant part of the metastatic potential of these lesions, and requires the rapid and efficient reorganization of the transcriptional landscape of melanoma cells. A large part of the non-coding genome cooperates to control gene expression, specifically through the activity of enhancers (ENHs). In this study, we aimed to identify ex vivo the network of active ENHs and to outline their cooperative interactions in supporting transcriptional adaptation during melanoma metastatic progression. We conducted a genome-wide analysis to map active ENHs distribution in a retrospective cohort of 39 melanoma patients, comparing the profiles obtained in primary (N = 19) and metastatic (N = 20) melanoma lesions. Unsupervised clustering showed that the profile for acetylated histone H3 at lysine 27 (H3K27ac) efficiently segregates lesions into three different clusters corresponding to progressive stages of the disease. We reconstructed the map of super-ENHs (SEs) and cooperative ENHs that associate with metastatic progression in melanoma, which showed that cooperation among regulatory elements is a mandatory requirement for transcriptional plasticity. We also showed that these elements carry out specialized and non-redundant functions, and indicated the existence of a hierarchical organization, with SEs on top as masterminds of the entire transcriptional program and classical ENHs as executors. By providing an innovative vision of how the chromatin landscape of melanoma works during metastatic spreading, our data also point out the need to integrate functional profiling in the analysis of cancer lesions to increase definition and improve interpretation of tumor heterogeneity.

Where and when to start: Regulating DNA replication origin activity in eukaryotic genomes

In eukaryotic genomes, hundreds to thousands of potential start sites of DNA replication named origins are dispersed across each of the linear chromosomes. During S-phase, only a subset of origins is selected in a stochastic manner to assemble bidirectional replication forks and initiate DNA synthesis. Despite substantial progress in our understanding of this complex process, a comprehensive 'identity code' that defines origins based on specific nucleotide sequences, DNA structural features, the local chromatin environment, or 3D genome architecture is still missing. In this article, we review the genetic and epigenetic features of replication origins in yeast and metazoan chromosomes and highlight recent insights into how this flexibility in origin usage contributes to nuclear organization, cell growth, differentiation, and genome stability.

Systematic characterization of the HOXA9 downstream targets in MLL-r leukemia by noncoding CRISPR screens

Accumulating evidence indicates that HOXA9 dysregulation is necessary and sufficient for leukemic transformation and maintenance. However, it remains largely unknown how HOXA9, as a homeobox transcriptional factor, binds to noncoding regulatory sequences and controls the downstream genes. Here, we conduct dropout CRISPR screens against 229 HOXA9-bound peaks identified by ChIP-seq. Integrative data analysis identifies reproducible noncoding hits, including those located in the distal enhancer of FLT3 and intron of CDK6. The Cas9-editing and dCas9-KRAB silencing of the HOXA9-bound sites significantly reduce corresponding gene transcription and impair cell proliferation in vitro, and in vivo by transplantation into NSG female mice. In addition, RNA-seq, Q-PCR analysis, chromatin accessibility change, and chromatin conformation evaluation uncover the noncoding regulation mechanism of HOXA9 and its functional downstream genes. In summary, our work improves our understanding of how HOXA9-associated transcription programs reconstruct the regulatory network specifying MLL-r dependency.

Multi-level 3D genome organization deteriorates during breast cancer progression

Breast cancer entails intricate alterations in genome organization and expression. However, how three-dimensional (3D) chromatin structure changes in the progression from a normal to a breast cancer malignant state remains unknown. To address this, we conducted an analysis combining Hi-C data with lamina-associated domains (LADs), epigenomic marks, and gene expression in an in vitro model of breast cancer progression. Our results reveal that while the fundamental properties of topologically associating domains (TADs) remain largely stable, significant changes occur in the organization of compartments and subcompartments. These changes are closely correlated with alterations in the expression of oncogenic genes. We also observe a restructuring of TAD-TAD interactions, coinciding with a loss of spatial compartmentalization and radial positioning of the 3D genome. Notably, we identify a previously unrecognized interchromosomal insertion event, wherein a locus on chromosome 8 housing the MYC oncogene is inserted into a highly active subcompartment on chromosome 10. This insertion leads to the formation of de novo enhancer contacts and activation of the oncogene, illustrating how structural variants can interact with the 3D genome to drive oncogenic states. In summary, our findings provide evidence for the degradation of genome organization at multiple scales during breast cancer progression revealing novel relationships between genome 3D structure and oncogenic processes.

OVOL2 sustains postnatal thymic epithelial cell identity

Distinct pathways and molecules may support embryonic versus postnatal thymic epithelial cell (TEC) development and maintenance. Here, we identify a mechanism by which TEC numbers and function are maintained postnatally. A viable missense allele (C120Y) of Ovol2, expressed ubiquitously or specifically in TECs, results in lymphopenia, in which T cell development is compromised by loss of medullary TECs and dysfunction of cortical TECs. We show that the epithelial identity of TECs is aberrantly subverted towards a mesenchymal state in OVOL2-deficient mice. We demonstrate that OVOL2 inhibits the epigenetic regulatory BRAF-HDAC complex, specifically disrupting RCOR1-LSD1 interaction. This causes inhibition of LSD1-mediated H3K4me2 demethylation, resulting in chromatin accessibility and transcriptional activation of epithelial genes. Thus, OVOL2 controls the epigenetic landscape of TECs to enforce TEC identity. The identification of a non-redundant postnatal mechanism for TEC maintenance offers an entry point to understanding thymic involution, which normally begins in early adulthood.

MESIA: multi-epigenome sample integration approach for precise peak calling

The assay for transposase-accessible chromatin with sequencing (ATAC-seq) is the most widely used method for measuring chromatin accessibility. Researchers have included multi-sample replication in ATAC-seq experimental designs. In epigenomic analysis, researchers should measure subtle changes in the peak by considering the read depth of individual samples. It is important to determine whether the peaks of each replication have an integrative meaning for the region of interest observed during multi-sample integration. We developed multi-epigenome sample integration approach for precise peak calling (MESIA), which integrates replication with high representativeness and reproducibility in multi-sample replication and determines the optimal peak. After identifying the reproducibility between all replications, our method integrated multiple samples determined as representative replicates. MESIA detected 6.06 times more peaks, and the value of the peaks was 1.32 times higher than the previously used method. MESIA is a shell-script-based open-source code that provides researchers involved in the epigenome with comprehensive insights.

Genomic analysis reveals HDAC1 regulates clinically relevant transcriptional programs in Pancreatic cancer

Novel strategies are needed to combat multidrug resistance in pancreatic ductal adenocarcinoma (PDAC). We applied genomic approaches to understand mechanisms of resistance in order to better inform treatment and precision medicine. Altered function of chromatin remodeling complexes contribute to chemoresistance. Our study generates and analyzes genomic and biochemical data from PDAC cells overexpressing HDAC1, a histone deacetylase involved in several chromatin remodeling complexes. We characterized the impact of overexpression on drug response, gene expression, HDAC1 binding, and chromatin structure using RNA-sequencing and ChIP-sequencing for HDAC1 and H3K27 acetylation. Integrative genomic analysis shows that HDAC1 overexpression promotes activation of key resistance pathways including epithelial to mesenchymal transition, cell cycle, and apoptosis through global chromatin remodeling. Target genes are similarly altered in patient tissues and show correlation with patient survival. We also demonstrate that direct targets of HDAC1 that also show altered chromatin are enriched near genes associated with altered GTPase activity. HDAC1 target genes identified using in vitro methods and observed in patient tissues were used to develop a clinically relevant nine-transcript signature associated with patient prognosis. Integration of multiple genomic and biochemical data types enables understanding of multidrug resistance and tumorigenesis in PDAC, a disease in desperate need of novel treatment strategies.

Regulation of the physiology and virulence of Ralstonia solanacearum by the second messenger 2',3'-cyclic guanosine monophosphate

Previous studies have demonstrated that bis-(3',5')-cyclic diguanosine monophosphate (bis-3',5'-c-di-GMP) is a ubiquitous second messenger employed by bacteria. Here, we report that 2',3'-cyclic guanosine monophosphate (2',3'-cGMP) controls the important biological functions, quorum sensing (QS) signaling systems and virulence in Ralstonia solanacearum through the transcriptional regulator RSp0980. This signal specifically binds to RSp0980 with high affinity and thus abolishes the interaction between RSp0980 and the promoters of target genes. In-frame deletion of RSp0334, which contains an evolved GGDEF domain with a LLARLGGDQF motif required to catalyze 2',3'-cGMP to (2',5')(3',5')-cyclic diguanosine monophosphate (2',3'-c-di-GMP), altered the abovementioned important phenotypes through increasing the intracellular 2',3'-cGMP levels. Furthermore, we found that 2',3'-cGMP, its receptor and the evolved GGDEF domain with a LLARLGGDEF motif also exist in the human pathogen Salmonella typhimurium. Together, our work provides insights into the unusual function of the GGDEF domain of RSp0334 and the special regulatory mechanism of 2',3'-cGMP signal in bacteria.

Improved quality metrics for association and reproducibility in chromatin accessibility data using mutual information

BACKGROUND

Correlation metrics are widely utilized in genomics analysis and often implemented with little regard to assumptions of normality, homoscedasticity, and independence of values. This is especially true when comparing values between replicated sequencing experiments that probe chromatin accessibility, such as assays for transposase-accessible chromatin via sequencing (ATAC-seq). Such data can possess several regions across the human genome with little to no sequencing depth and are thus non-normal with a large portion of zero values. Despite distributed use in the epigenomics field, few studies have evaluated and benchmarked how correlation and association statistics behave across ATAC-seq experiments with known differences or the effects of removing specific outliers from the data. Here, we developed a computational simulation of ATAC-seq data to elucidate the behavior of correlation statistics and to compare their accuracy under set conditions of reproducibility.

RESULTS

Using these simulations, we monitored the behavior of several correlation statistics, including the Pearson's R and Spearman's [Formula: see text] coefficients as well as Kendall's [Formula: see text] and Top-Down correlation. We also test the behavior of association measures, including the coefficient of determination R[Formula: see text], Kendall's W, and normalized mutual information. Our experiments reveal an insensitivity of most statistics, including Spearman's [Formula: see text], Kendall's [Formula: see text], and Kendall's W, to increasing differences between simulated ATAC-seq replicates. The removal of co-zeros (regions lacking mapped sequenced reads) between simulated experiments greatly improves the estimates of correlation and association. After removing co-zeros, the R[Formula: see text] coefficient and normalized mutual information display the best performance, having a closer one-to-one relationship with the known portion of shared, enhanced loci between simulated replicates. When comparing values between experimental ATAC-seq data using a random forest model, mutual information best predicts ATAC-seq replicate relationships.

CONCLUSIONS

Collectively, this study demonstrates how measures of correlation and association can behave in epigenomics experiments. We provide improved strategies for quantifying relationships in these increasingly prevalent and important chromatin accessibility assays.

Adaptation of CUT&RUN for use in African trypanosomes

This Cleavage Under Targets and Release Using Nuclease (CUT&RUN) protocol produces genomic occupancy data for a protein of interest in the protozoan parasite Trypanosoma brucei. The data produced is analyzed in a similar way as that produced by ChIP-seq. While we describe the protocol for parasites carrying an epitope tag for the protein of interest, antibodies against the native protein could be used for the same purpose.

Loss of histone reader Phf7 leads to immune pathways activation via endogenous retroviruses during spermiogenesis

Genetic studies have elucidated the critical roles of Phf7 in germline development in animals; however, the exact etiology of Phf7 mutations leading to male infertility and the possibility of mechanism-based therapy are still unclear and warrant further investigation. Using the Phf7 knockout mouse model, we verified that genetic defects were responsible for male infertility by preventing histone-to-protamine exchange, as previously reported. The deficiency of spermatogenesis caused by Phf7 deletion through the endogenous retrovirus-mediated activation of the immune pathway is a common mechanism of infertility. Furthermore, we identified PPARα as a promising target of immunity and inflammation in the testis, where endogenous retroviruses are suppressed, and Phf7 as a crucial regulator of endogenous retrovirus-mediated immune regulation and revealed its role as an epigenetic reader. The loss of Phf7 activates immune pathways, which can be rescued by the PPARα agonist astaxanthin. These results showed that astaxanthin is a potential therapeutic agent for treating male infertility. The findings in our study provide insights into the molecular mechanisms underlying male infertility and suggest potential targets for future research and therapeutic development.