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

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.

Multi-omics characterization of the monkeypox virus infection

Multiple omics analyzes of Vaccinia virus (VACV) infection have defined molecular characteristics of poxvirus biology. However, little is known about the monkeypox (mpox) virus (MPXV) in humans, which has a different disease manifestation despite its high sequence similarity to VACV. Here, we perform an in-depth multi-omics analysis of the transcriptome, proteome, and phosphoproteome signatures of MPXV-infected primary human fibroblasts to gain insights into the virus-host interplay. In addition to expected perturbations of immune-related pathways, we uncover regulation of the HIPPO and TGF-β pathways. We identify dynamic phosphorylation of both host and viral proteins, which suggests that MAPKs are key regulators of differential phosphorylation in MPXV-infected cells. Among the viral proteins, we find dynamic phosphorylation of H5 that influenced the binding of H5 to dsDNA. Our extensive dataset highlights signaling events and hotspots perturbed by MPXV, extending the current knowledge on poxviruses. We use integrated pathway analysis and drug-target prediction approaches to identify potential drug targets that affect virus growth. Functionally, we exemplify the utility of this approach by identifying inhibitors of MTOR, CHUK/IKBKB, and splicing factor kinases with potent antiviral efficacy against MPXV and VACV.

Comprehensive mapping and modelling of the rice regulome landscape unveils the regulatory architecture underlying complex traits

Unraveling the regulatory mechanisms that govern complex traits is pivotal for advancing crop improvement. Here we present a comprehensive regulome atlas for rice (Oryza sativa), charting the chromatin accessibility across 23 distinct tissues from three representative varieties. Our study uncovers 117,176 unique open chromatin regions (OCRs), accounting for ~15% of the rice genome, a notably higher proportion compared to previous reports in plants. Integrating RNA-seq data from matched tissues, we confidently predict 59,075 OCR-to-gene links, with enhancers constituting 69.54% of these associations, including many known enhancer-to-gene links. Leveraging this resource, we re-evaluate genome-wide association study results and discover a previously unknown function of OsbZIP06 in seed germination, which we subsequently confirm through experimental validation. We optimize deep learning models to decode regulatory grammar, achieving robust modeling of tissue-specific chromatin accessibility. This approach allows to predict cross-variety regulatory dynamics from genomic sequences, shedding light on the genetic underpinnings of cis-regulatory divergence and morphological disparities between varieties. Overall, our study establishes a foundational resource for rice functional genomics and precision molecular breeding, providing valuable insights into regulatory mechanisms governing complex traits.

Chrom-seq identifies RNAs at chromatin marks

Chromatin marks are associated with transcriptional regulatory activities. However, very few lncRNAs have been characterized with the role in regulating epigenetic marks, largely due to the technical difficulty in identifying chromatin-associating RNA. Current methods are largely limited by the availability of ChIP-grade antibody and the crosslinking, which generates high noise. Here, we developed a method termed Chrom-seq to efficiently capture RNAs associated with various chromatin marks in living cells. Chrom-seq jointly applies highly specific chromatin mark reader with APEX2, which catalyzes the oxidation of biotin-aniline to label the adjacent RNAs for isolation by streptavidin-coated beads. Using the readers of mCBX7/dPC, mCBX1, and mTAF3, we detected RNA species significantly associated with H3K27me3, H3K9me3, and H3K4me3, respectively. We demonstrated that Chrom-seq outperformed other equivalent methods in terms of sensitivity, efficiency, and cost of practice. It provides an antibody-free approach to systematically map RNAs at chromatin marks with potential regulatory roles in epigenetic events.

Transcriptomic and epigenomic insights into pectoral muscle fiber formation at the late embryonic development in pure chicken lines

Long-term intensive genetic selection has led to significant differences between broiler and layer chickens, which are evident during the embryonic period. Despite this, there is a paucity of research on the genetic regulation of the initial formation of muscle fiber morphology in chick embryos. Embryonic d 17 (E17) is the key time point for myoblast fusion completion and muscle fiber morphology formation in chickens. This study aimed to explore the genetic regulatory mechanisms underlying the early muscle fiber morphology establishment in broiler chickens of Cornish (CC) and White Plymouth Rock (RR) and layer chickens of White Leghorn (WW) at E17 using the transcriptomic and chromatin accessibility sequencing of pectoral major muscles. The results showed that broiler chickens exhibited significant higher embryo weight and pectoral major muscle weight at E17 compared to layer chickens (P = 0.000). A total of 1,278, 1,248, and 892 differentially expressed genes (DEGs) of RNA-seq data were identified between CC vs. WW, RR vs. WW, and CC vs. RR, separately. All DEGs were combined for cluster analysis and they were divided into 6 clusters, including cluster 1 with higher expression in broilers and cluster 6 with higher expression in layers. DEGs in cluster 1 were enriched in terms related to macrophage activation (P = 0.002) and defense response to bacteria (P = 0.002), while DEGs in cluster 6 showed enrichment in protein-DNA complex (P = 0.003) and monooxygenase activity (P = 0.000). ATAC-seq data analysis identified a total of 38,603 peaks, with 13,051 peaks for CC, 18,780 peaks for RR, and 6,772 peaks for WW. Integrative analysis of transcriptomic and chromatin accessibility data revealed GOLM1, ISLR2, and TOPAZ1 were commonly upregulated genes in CC and RR. Furthermore, screening of all upregulated DEGs in cluster 1 from CC and RR identified GOLM1, ISLR2, and HNMT genes associated with neuroimmune functions and MYOM3 linked to muscle morphology development, showing significantly elevated expression in broiler chickens compared to layer chickens. These findings suggest active neural system connectivity during the initial formation of muscle fiber morphology in embryonic period, highlighting the early interaction between muscle fiber formation morphology and the nervous system. This study provides novel insights into late chick embryo development and lays a deeper foundation for further research.

Non-coding variants impact cis-regulatory coordination in a cell type-specific manner

BACKGROUND

Interactions among cis-regulatory elements (CREs) play a crucial role in gene regulation. Various approaches have been developed to map these interactions genome-wide, including those relying on interindividual epigenomic variation to identify groups of covariable regulatory elements, referred to as chromatin modules (CMs). While CM mapping allows to investigate the relationship between chromatin modularity and gene expression, the computational principles used for CM identification vary in their application and outcomes.

RESULTS

We comprehensively evaluate and streamline existing CM mapping tools and present guidelines for optimal utilization of epigenome data from a diverse population of individuals to assess regulatory coordination across the human genome. We showcase the effectiveness of our recommended practices by analyzing distinct cell types and demonstrate cell type specificity of CRE interactions in CMs and their relevance for gene expression. Integration of genotype information revealed that many non-coding disease-associated variants affect the activity of CMs in a cell type-specific manner by affecting the binding of cell type-specific transcription factors. We provide example cases that illustrate in detail how CMs can be used to deconstruct GWAS loci, assess variable expression of cell surface receptors in immune cells, and reveal how genetic variation can impact the expression of prognostic markers in chronic lymphocytic leukemia.

CONCLUSIONS

Our study presents an optimal strategy for CM mapping and reveals how CMs capture the coordination of CREs and its impact on gene expression. Non-coding genetic variants can disrupt this coordination, and we highlight how this may lead to disease predisposition in a cell type-specific manner.

PRC2-AgeIndex as a universal biomarker of aging and rejuvenation

DNA methylation (DNAm) is one of the most reliable biomarkers of aging across mammalian tissues. While the age-dependent global loss of DNAm has been well characterized, DNAm gain is less characterized. Studies have demonstrated that CpGs which gain methylation with age are enriched in Polycomb Repressive Complex 2 (PRC2) targets. However, whole-genome examination of all PRC2 targets as well as determination of the pan-tissue or tissue-specific nature of these associations is lacking. Here, we show that low-methylated regions (LMRs) which are highly bound by PRC2 in embryonic stem cells (PRC2 LMRs) gain methylation with age in all examined somatic mitotic cells. We estimated that this epigenetic change represents around 90% of the age-dependent DNAm gain genome-wide. Therefore, we propose the "PRC2-AgeIndex," defined as the average DNAm in PRC2 LMRs, as a universal biomarker of cellular aging in somatic cells which can distinguish the effect of different anti-aging interventions.

Chromatin endogenous cleavage provides a global view of RNA polymerase II transcription kinetics

Chromatin immunoprecipitation (ChIP-seq) is the most common approach to observe global binding of proteins to DNA in vivo. The occupancy of transcription factors (TFs) from ChIP-seq agrees well with an alternative method, chromatin endogenous cleavage (ChEC-seq2). However, ChIP-seq and ChEC-seq2 reveal strikingly different patterns of enrichment of yeast RNA polymerase II. We hypothesized that this reflects distinct populations of RNAPII, some of which are captured by ChIP-seq and some of which are captured by ChEC-seq2. RNAPII association with enhancers and promoters - predicted from biochemical studies - is detected well by ChEC-seq2 but not by ChIP-seq. Enhancer/promoter bound RNAPII correlates with transcription levels and matches predicted occupancy based on published rates of enhancer recruitment, preinitiation assembly, initiation, elongation and termination. The occupancy from ChEC-seq2 allowed us to develop a stochastic model for global kinetics of RNAPII transcription which captured both the ChEC-seq2 data and changes upon chemical-genetic perturbations to transcription. Finally, RNAPII ChEC-seq2 and kinetic modeling suggests that a mutation in the Gcn4 transcription factor that blocks interaction with the NPC destabilizes promoter-associated RNAPII without altering its recruitment to the enhancer.

Epigenetic alterations affecting hematopoietic regulatory networks as drivers of mixed myeloid/lymphoid leukemia

Leukemias with ambiguous lineage comprise several loosely defined entities, often without a clear mechanistic basis. Here, we extensively profile the epigenome and transcriptome of a subgroup of such leukemias with CpG Island Methylator Phenotype. These leukemias exhibit comparable hybrid myeloid/lymphoid epigenetic landscapes, yet heterogeneous genetic alterations, suggesting they are defined by their shared epigenetic profile rather than common genetic lesions. Gene expression enrichment reveals similarity with early T-cell precursor acute lymphoblastic leukemia and a lymphoid progenitor cell of origin. In line with this, integration of differential DNA methylation and gene expression shows widespread silencing of myeloid transcription factors. Moreover, binding sites for hematopoietic transcription factors, including CEBPA, SPI1 and LEF1, are uniquely inaccessible in these leukemias. Hypermethylation also results in loss of CTCF binding, accompanied by changes in chromatin interactions involving key transcription factors. In conclusion, epigenetic dysregulation, and not genetic lesions, explains the mixed phenotype of this group of leukemias with ambiguous lineage. The data collected here constitute a useful and comprehensive epigenomic reference for subsequent studies of acute myeloid leukemias, T-cell acute lymphoblastic leukemias and mixed-phenotype leukemias.

Transcriptomic and epigenomic landscapes of muscle growth during the postnatal period of broilers

BACKGROUND

Broilers stand out as one of the fastest-growing livestock globally, making a substantial contribution to animal meat production. However, the molecular and epigenetic mechanisms underlying the rapid growth and development of broiler chickens are still unclear. This study aims to explore muscle development patterns and regulatory networks during the postnatal rapid growth phase of fast-growing broilers. We measured the growth performance of Cornish (CC) and White Plymouth Rock (RR) over a 42-d period. Pectoral muscle samples from both CC and RR were randomly collected at day 21 after hatching (D21) and D42 for RNA-seq and ATAC-seq library construction.

RESULTS

The consistent increase in body weight and pectoral muscle weight across both breeds was observed as they matured, with CC outpacing RR in terms of weight at each stage of development. Differential expression analysis identified 398 and 1,129 genes in the two dimensions of breeds and ages, respectively. A total of 75,149 ATAC-seq peaks were annotated in promoter, exon, intron and intergenic regions, with a higher number of peaks in the promoter and intronic regions. The age-biased genes and breed-biased genes of RNA-seq were combined with the ATAC-seq data for subsequent analysis. The results spotlighted the upregulation of ACTC1 and FDPS at D21, which were primarily associated with muscle structure development by gene cluster enrichment. Additionally, a noteworthy upregulation of MUSTN1, FOS and TGFB3 was spotted in broiler chickens at D42, which were involved in cell differentiation and muscle regeneration after injury, suggesting a regulatory role of muscle growth and repair.

CONCLUSIONS

This work provided a regulatory network of postnatal broiler chickens and revealed ACTC1 and MUSTN1 as the key responsible for muscle development and regeneration. Our findings highlight that rapid growth in broiler chickens triggers ongoing muscle damage and subsequent regeneration. These findings provide a foundation for future research to investigate the functional aspects of muscle development.

A proximal enhancer regulates RORA expression during early human Th17 cell differentiation

Gene regulatory elements, such as enhancers, greatly influence cell identity by tuning the transcriptional activity of specific cell types. Dynamics of enhancer landscape during early human Th17 cell differentiation remains incompletely understood. Leveraging ATAC-seq-based profiling of chromatin accessibility and comprehensive analysis of key histone marks, we identified a repertoire of enhancers that potentially exert control over the fate specification of Th17 cells. We found 23 SNPs associated with autoimmune diseases within Th17-enhancers that precisely overlapped with the binding sites of transcription factors actively engaged in T-cell functions. Among the Th17-specific enhancers, we identified an enhancer in the intron of RORA and demonstrated that this enhancer positively regulates RORA transcription. Moreover, CRISPR-Cas9-mediated deletion of a transcription factor binding site-rich region within the identified RORA enhancer confirmed its role in regulating RORA transcription. These findings provide insights into the potential mechanism by which the RORA enhancer orchestrates Th17 differentiation.

Multiomic profiling of transcription factor binding and function in human brain

Transcription factors (TFs) orchestrate gene expression programs crucial for brain function, but we lack detailed information about TF binding in human brain tissue. We generated a multiomic resource (ChIP-seq, ATAC-seq, RNA-seq, DNA methylation) on bulk tissues and sorted nuclei from several postmortem brain regions, including binding maps for more than 100 TFs. We demonstrate improved measurements of TF activity, including motif recognition and gene expression modeling, upon identification and removal of high TF occupancy regions. Further, predictive TF binding models demonstrate a bias for these high-occupancy sites. Neuronal TFs SATB2 and TBR1 bind unique regions depleted for such sites and promote neuronal gene expression. Binding sites for TFs, including TBR1 and PKNOX1, are enriched for risk variants associated with neuropsychiatric disorders, predominantly in neurons. This work, titled BrainTF, is a powerful resource for future studies seeking to understand the roles of specific TFs in regulating gene expression in the human brain.

Cactus: A user-friendly and reproducible ATAC-Seq and mRNA-Seq analysis pipeline for data preprocessing, differential analysis, and enrichment analysis

The ever decreasing cost of Next-Generation Sequencing coupled with the emergence of efficient and reproducible analysis pipelines has rendered genomic methods more accessible. However, downstream analyses are basic or missing in most workflows, creating a significant barrier for non-bioinformaticians. To help close this gap, we developed Cactus, an end-to-end pipeline for analyzing ATAC-Seq and mRNA-Seq data, either separately or jointly. Its Nextflow-, container-, and virtual environment-based architecture ensures efficient and reproducible analyses. Cactus preprocesses raw reads, conducts differential analyses between conditions, and performs enrichment analyses in various databases, including DNA-binding motifs, ChIP-Seq binding sites, chromatin states, and ontologies. We demonstrate the utility of Cactus in a multi-modal and multi-species case study as well as by showcasing its unique capabilities as compared to other ATAC-Seq pipelines. In conclusion, Cactus can assist researchers in gaining comprehensive insights from chromatin accessibility and gene expression data in a quick, user-friendly, and reproducible manner.

MICAL2 Is a Super Enhancer Associated Gene that Promotes Pancreatic Cancer Growth and Metastasis

Pancreatic ductal adenocarcinoma (PDAC) remains one of the deadliest solid cancers and thus identifying more effective therapies is a major unmet need. In this study we characterized the super enhancer (SE) landscape of human PDAC to identify novel, potentially targetable, drivers of the disease. Our analysis revealed that MICAL2 is a super enhancer-associated gene in human PDAC. MICAL2 is a flavin monooxygenase that induces actin depolymerization and indirectly promotes SRF transcription by modulating the availability of serum response factor coactivators myocardin related transcription factors (MRTF-A and MRTF-B). We found that MICAL2 is overexpressed in PDAC and correlates with poor patient prognosis. Transcriptional analysis revealed that MICAL2 upregulates KRAS and EMT signaling pathways, contributing to tumor growth and metastasis. In loss and gain of function experiments in human and mouse PDAC cells, we observed that MICAL2 promotes both ERK1/2 and AKT activation. Consistent with its role in actin depolymerization and KRAS signaling, loss of MICAL2 expression also inhibited macropinocytosis. Through in vitro phenotypic analyses, we show that MICAL2, MRTF-A and MRTF-B influence PDAC cell proliferation, migration and promote cell cycle progression. Importantly, we demonstrate that MICAL2 is essential for in vivo tumor growth and metastasis. Interestingly, we find that MRTF-B, but not MRTF-A, phenocopies MICAL2-driven phenotypes in vivo . This study highlights the multiple ways in which MICAL2 impacts PDAC biology and suggests that its inhibition may impede PDAC progression. Our results provide a foundation for future investigations into the role of MICAL2 in PDAC and its potential as a target for therapeutic intervention.

BnaABF3 and BnaMYB44 regulate the transcription of zeaxanthin epoxidase genes in carotenoid and abscisic acid biosynthesis

Zeaxanthin epoxidase (ZEP) is a key enzyme that catalyzes the conversion of zeaxanthin to violaxanthin in the carotenoid and abscisic acid (ABA) biosynthesis pathways. The rapeseed (Brassica napus) genome has 4 ZEP (BnaZEP) copies that are suspected to have undergone subfunctionalization, yet the 4 genes' underlying regulatory mechanisms remain unknown. Here, we genetically confirmed the functional divergence of the gene pairs BnaA09.ZEP/BnaC09.ZEP and BnaA07.ZEP/BnaC07.ZEP, which encode enzymes with tissue-specific roles in carotenoid and ABA biosynthesis in flowers and leaves, respectively. Molecular and transgenic experiments demonstrated that each BnaZEP pair is transcriptionally regulated via ABA-responsive element-binding factor 3 s (BnaABF3s) and BnaMYB44s as common and specific regulators, respectively. BnaABF3s directly bound to the promoters of all 4 BnaZEPs and activated their transcription, with overexpression of individual BnaABF3s inducing BnaZEP expression and ABA accumulation under drought stress. Conversely, loss of BnaABF3s function resulted in lower expression of several genes functioning in carotenoid and ABA metabolism and compromised drought tolerance. BnaMYB44s specifically targeted and repressed the expression of BnaA09.ZEP/BnaC09.ZEP but not BnaA07.ZEP/BnaC07.ZEP. Overexpression of BnaA07.MYB44 resulted in increased carotenoid content and an altered carotenoid profile in petals. Additionally, RNA-seq analysis indicated that BnaMYB44s functions as a repressor in phenylpropanoid and flavonoid biosynthesis. These findings provide clear evidence for the subfunctionalization of duplicated genes and contribute to our understanding of the complex regulatory network involved in carotenoid and ABA biosynthesis in B. napus.

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.

Single-cell multiomics profiling reveals heterogeneous transcriptional programs and microenvironment in DSRCTs

Desmoplastic small round cell tumor (DSRCT) is a rare, aggressive sarcoma driven by the EWSR1::WT1 chimeric transcription factor. Despite this unique oncogenic driver, DSRCT displays a polyphenotypic differentiation of unknown causality. Using single-cell multi-omics on 12 samples from five patients, we find that DSRCT tumor cells cluster into consistent subpopulations with partially overlapping lineage- and metabolism-related transcriptional programs. In vitro modeling shows that high EWSR1::WT1 DNA-binding activity associates with most lineage-related states, in contrast to glycolytic and profibrotic states. Single-cell chromatin accessibility analysis suggests that EWSR1::WT1 binding site variability may drive distinct lineage-related transcriptional programs, supporting some level of cell-intrinsic plasticity. Spatial transcriptomics reveals that glycolytic and profibrotic states specifically localize within hypoxic niches at the periphery of tumor cell islets, suggesting an additional role of tumor cell-extrinsic microenvironmental cues. We finally identify a single-cell transcriptomics-derived epithelial signature associated with improved patient survival, highlighting the clinical relevance of our findings.

C-terminally phosphorylated p27 activates self-renewal driver genes to program cancer stem cell expansion, mammary hyperplasia and cancer

In many cancers, a stem-like cell subpopulation mediates tumor initiation, dissemination and drug resistance. Here, we report that cancer stem cell (CSC) abundance is transcriptionally regulated by C-terminally phosphorylated p27 (p27pT157pT198). Mechanistically, this arises through p27 co-recruitment with STAT3/CBP to gene regulators of CSC self-renewal including MYC, the Notch ligand JAG1, and ANGPTL4. p27pTpT/STAT3 also recruits a SIN3A/HDAC1 complex to co-repress the Pyk2 inhibitor, PTPN12. Pyk2, in turn, activates STAT3, creating a feed-forward loop increasing stem-like properties in vitro and tumor-initiating stem cells in vivo. The p27-activated gene profile is over-represented in STAT3 activated human breast cancers. Furthermore, mammary transgenic expression of phosphomimetic, cyclin-CDK-binding defective p27 (p27CK-DD) increases mammary duct branching morphogenesis, yielding hyperplasia and microinvasive cancers that can metastasize to liver, further supporting a role for p27pTpT in CSC expansion. Thus, p27pTpT interacts with STAT3, driving transcriptional programs governing stem cell expansion or maintenance in normal and cancer tissues.

Integrated ATAC-seq and RNA-seq Analysis of In Vitro Cultured Skeletal Muscle Satellite Cells to Understand Changes in Cell Proliferation

Skeletal muscle satellite cells, the resident stem cells in pig skeletal muscle, undergo proliferation and differentiation to enable muscle tissue repair. The proliferative and differentiative abilities of these cells gradually decrease during in vitro cultivation as the cell passage number increases. Despite extensive research, the precise molecular mechanisms that regulate this process are not fully understood. To bridge this knowledge gap, we conducted transcriptomic analysis of skeletal muscle satellite cells during in vitro cultivation to quantify passage number-dependent changes in the expression of genes associated with proliferation. Additionally, we explored the relationships between gene transcriptional activity and chromatin accessibility using transposase-accessible chromatin sequencing. This revealed the closure of numerous open chromatin regions, which were primarily located in intergenic regions, as the cell passage number increased. Integrated analysis of the transcriptomic and epigenomic data demonstrated a weak correlation between gene transcriptional activity and chromatin openness in expressed genic regions; although some genes (e.g., GNB4 and FGD5) showed consistent relationships between gene expression and chromatin openness, a substantial number of differentially expressed genes had no clear association with chromatin openness in expressed genic regions. The p53-p21-RB signaling pathway may play a critical regulatory role in cell proliferation processes. The combined transcriptomic and epigenomic approach taken here provided key insights into changes in gene expression and chromatin openness during in vitro cultivation of skeletal muscle satellite cells. These findings enhance our understanding of the intricate mechanisms underlying the decline in cellular proliferation capacity in cultured cells.

An updated compendium and reevaluation of the evidence for nuclear transcription factor occupancy over the mitochondrial genome

In most eukaryotes, mitochondrial organelles contain their own genome, usually circular, which is the remnant of the genome of the ancestral bacterial endosymbiont that gave rise to modern mitochondria. Mitochondrial genomes are dramatically reduced in their gene content due to the process of endosymbiotic gene transfer to the nucleus; as a result most mitochondrial proteins are encoded in the nucleus and imported into mitochondria. This includes the components of the dedicated mitochondrial transcription and replication systems and regulatory factors, which are entirely distinct from the information processing systems in the nucleus. However, since the 1990s several nuclear transcription factors have been reported to act in mitochondria, and previously we identified 8 human and 3 mouse transcription factors (TFs) with strong localized enrichment over the mitochondrial genome using ChIP-seq (Chromatin Immunoprecipitation) datasets from the second phase of the ENCODE (Encyclopedia of DNA Elements) Project Consortium. Here, we analyze the greatly expanded in the intervening decade ENCODE compendium of TF ChIP-seq datasets (a total of 6,153 ChIP experiments for 942 proteins, of which 763 are sequence-specific TFs) combined with interpretative deep learning models of TF occupancy to create a comprehensive compendium of nuclear TFs that show evidence of association with the mitochondrial genome. We find some evidence for chrM occupancy for 50 nuclear TFs and two other proteins, with bZIP TFs emerging as most likely to be playing a role in mitochondria. However, we also observe that in cases where the same TF has been assayed with multiple antibodies and ChIP protocols, evidence for its chrM occupancy is not always reproducible. In the light of these findings, we discuss the evidential criteria for establishing chrM occupancy and reevaluate the overall compendium of putative mitochondrial-acting nuclear TFs.

Integrative multi-omics analyses to identify the genetic and functional mechanisms underlying ovarian cancer risk regions

To identify credible causal risk variants (CCVs) associated with different histotypes of epithelial ovarian cancer (EOC), we performed genome-wide association analysis for 470,825 genotyped and 10,163,797 imputed SNPs in 25,981 EOC cases and 105,724 controls of European origin. We identified five histotype-specific EOC risk regions (p value <5 × 10-8) and confirmed previously reported associations for 27 risk regions. Conditional analyses identified an additional 11 signals independent of the primary signal at six risk regions (p value <10-5). Fine mapping identified 4,008 CCVs in these regions, of which 1,452 CCVs were located in ovarian cancer-related chromatin marks with significant enrichment in active enhancers, active promoters, and active regions for CCVs from each EOC histotype. Transcriptome-wide association and colocalization analyses across histotypes using tissue-specific and cross-tissue datasets identified 86 candidate susceptibility genes in known EOC risk regions and 32 genes in 23 additional genomic regions that may represent novel EOC risk loci (false discovery rate <0.05). Finally, by integrating genome-wide HiChIP interactome analysis with transcriptome-wide association study (TWAS), variant effect predictor, transcription factor ChIP-seq, and motifbreakR data, we identified candidate gene-CCV interactions at each locus. This included risk loci where TWAS identified one or more candidate susceptibility genes (e.g., HOXD-AS2, HOXD8, and HOXD3 at 2q31) and other loci where no candidate gene was identified (e.g., MYC and PVT1 at 8q24) by TWAS. In summary, this study describes a functional framework and provides a greater understanding of the biological significance of risk alleles and candidate gene targets at EOC susceptibility loci identified by a genome-wide association study.

Native and tagged CENP-A histones are functionally inequivalent

BACKGROUND

Over the past several decades, the use of biochemical and fluorescent tags has elucidated mechanistic and cytological processes that would otherwise be impossible. The challenging nature of certain nuclear proteins includes low abundancy, poor antibody recognition, and transient dynamics. One approach to get around those issues is the addition of a peptide or larger protein tag to the target protein to improve enrichment, purification, and visualization. However, many of these studies were done under the assumption that tagged proteins can fully recapitulate native protein function.

RESULTS

We report that when C-terminally TAP-tagged CENP-A histone variant is introduced, it undergoes altered kinetochore protein binding, differs in post-translational modifications (PTMs), utilizes histone chaperones that differ from that of native CENP-A, and can partially displace native CENP-A in human cells. Additionally, these tagged CENP-A-containing nucleosomes have reduced centromeric incorporation at early G1 phase and poorly associates with linker histone H1.5 compared to native CENP-A nucleosomes.

CONCLUSIONS

These data suggest expressing tagged versions of histone variant CENP-A may result in unexpected utilization of non-native pathways, thereby altering the biological function of the histone variant.

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.

Metabolic reprogramming during Candida albicans planktonic-biofilm transition is modulated by the transcription factors Zcf15 and Zcf26

Candida albicans is a commensal of the human microbiota that can form biofilms on implanted medical devices. These biofilms are tolerant to antifungals and to the host immune system. To identify novel genes modulating C. albicans biofilm formation, we performed a large-scale screen with 2,454 C. albicans doxycycline-dependent overexpression strains and identified 16 genes whose overexpression significantly hampered biofilm formation. Among those, overexpression of the ZCF15 and ZCF26 paralogs that encode transcription factors and have orthologs only in biofilm-forming species of the Candida clade, caused impaired biofilm formation both in vitro and in vivo. Interestingly, overexpression of ZCF15 impeded biofilm formation without any defect in hyphal growth. Transcript profiling, transcription factor binding, and phenotypic microarray analyses conducted upon overexpression of ZCF15 and ZCF26 demonstrated their role in reprogramming cellular metabolism by regulating central metabolism including glyoxylate and tricarboxylic acid cycle genes. Taken together, this study has identified a new set of biofilm regulators, including ZCF15 and ZCF26, that appear to control biofilm development through their specific role in metabolic remodeling.

A disease-associated gene desert directs macrophage inflammation through ETS2

Increasing rates of autoimmune and inflammatory disease present a burgeoning threat to human health1. This is compounded by the limited efficacy of available treatments1 and high failure rates during drug development2, highlighting an urgent need to better understand disease mechanisms. Here we show how functional genomics could address this challenge. By investigating an intergenic haplotype on chr21q22-which has been independently linked to inflammatory bowel disease, ankylosing spondylitis, primary sclerosing cholangitis and Takayasu's arteritis3-6-we identify that the causal gene, ETS2, is a central regulator of human inflammatory macrophages and delineate the shared disease mechanism that amplifies ETS2 expression. Genes regulated by ETS2 were prominently expressed in diseased tissues and more enriched for inflammatory bowel disease GWAS hits than most previously described pathways. Overexpressing ETS2 in resting macrophages reproduced the inflammatory state observed in chr21q22-associated diseases, with upregulation of multiple drug targets, including TNF and IL-23. Using a database of cellular signatures7, we identified drugs that might modulate this pathway and validated the potent anti-inflammatory activity of one class of small molecules in vitro and ex vivo. Together, this illustrates the power of functional genomics, applied directly in primary human cells, to identify immune-mediated disease mechanisms and potential therapeutic opportunities.

Aberrant Super-Enhancer Landscape in Enzalutamide-Resistant Prostate Cancer Cells

Background: Castration-resistant prostate cancer (CRPC), which has developed resistance to next-generation antiandrogens, such as enzalutamide (Enz), is a lethal disease. Furthermore, transcriptional regulation by super enhancers (SEs) is crucial for the growth and spread of prostate cancer, as well as drug resistance. The functions of SEs, a significant class of noncoding DNA cis-regulatory elements, have been the subject of numerous recent studies in the field of cancer research. Materials and Methods: The goal of this research was to identify SEs associated with Enz resistance in C4-2B cells using chromatin immunoprecipitation sequencing and cleavage under targets and tagmentation (CUT&Tag). Using HOMER analysis to predict protein/gene-binding motifs, we identified master transcription factors (TFs) that may bind to SE sites. Using small interfering RNA, WST-1 assays, and qRT-PCR, we then confirmed the associations between TFs of SEs and Enz resistance. Results: A total of 999 SEs were screened from C4-2B EnzR cells in total. Incorporating analysis with RNA-seq data revealed 41 SEs to be strongly associated with the promotion of Enz resistance. In addition, we finally predicted that master TFs bind to SE-binding regions. Subsequently, we selected zinc finger protein 467 (ZFP467) and SMAD family member 3 to confirm the functional connections of master TFs with Enz resistance through SEs (ZNF467). Conclusions: In this study, SMAD3 and ZNF467 were found to be closely related to Enz-resistant CRPC. Our research uncovered a sizable group of SEs linked to Enz resistance in prostate cancer, dissected the mechanisms underlying SE Enz resistance, and shed light on potential clinical uses for SEs.

Atypical Modes of CTCF Binding Facilitate Tissue-Specific and Neuronal Activity-Dependent Gene Expression States

Multivalent binding of CTCF to variable DNA sequences is thought to underlie its ability to mediate diverse cellular functions. CTCF typically binds a 20 base-pair consensus DNA sequence, but the full diversity of CTCF binding sites (CBS) within the genome has not been interrogated. We assessed CTCF occupancy in cultured cortical neurons and observed surprisingly that ~ 22% of CBS lack the consensus CTCF motif. We report here that sequence diversity at most of these atypical CBS involves degeneracy at specific nucleotide positions within the consensus CTCF motif, which likely affect the binding of CTCF zinc fingers 6 and 7. This mode of atypical CTCF binding defines most CBS at gene promoters, as well as CBS that are dynamically altered during neural differentiation and following neuronal stimulation, revealing how atypical CTCF binding could influence gene activity. Dynamic CBS are distributed both within and outside loop anchors and TAD boundaries, suggesting both looping-dependent and independent roles for CTCF. Finally, we describe a second mode of atypical CTCF binding to DNA sequences that are completely unrelated to the consensus CTCF motif, which are enriched within the bodies of tissue-specific genes. These tissue-specific atypical CBS are also enriched in H3K27ac, which marks cis-regulatory elements within chromatin, including enhancers. Overall, these results indicate how atypical CBS could dynamically regulate gene activity patterns during differentiation, development, and in response to environmental cues.

Wemics: A Single-Base Resolution Methylation Quantification Method for Enhanced Prediction of Epigenetic Regulation

DNA methylation, an epigenetic mechanism that alters gene expression without changing DNA sequence, is essential for organism development and key biological processes like genomic imprinting and X-chromosome inactivation. Despite tremendous efforts in DNA methylation research, accurate quantification of cytosine methylation remains a challenge. Here, a single-base methylation quantification approach is introduced by weighting methylation of consecutive CpG sites (Wemics) in genomic regions. Wemics quantification of DNA methylation better predicts its regulatory impact on gene transcription and identifies differentially methylated regions (DMRs) with more biological relevance. Most Wemics-quantified DMRs in lung cancer are epigenetically conserved and recurrently occurred in other primary cancers from The Cancer Genome Atlas (TCGA), and their aberrant alterations can serve as promising pan-cancer diagnostic markers. It is further revealed that these detected DMRs are enriched in transcription factor (TF) binding motifs, and methylation of these TF binding motifs and TF expression synergistically regulate target gene expression. Using Wemics on epigenomic-transcriptomic data from the large lung cancer cohort, a dozen novel genes with oncogenic potential are discovered that are upregulated by hypomethylation but overlooked by other quantification methods. These findings increase the understanding of the epigenetic mechanism by which DNA methylation regulates gene expression.

Establishing an optimized ATAC-seq protocol for the maize

The advent of next-generation sequencing in crop improvement offers unprecedented insights into the chromatin landscape closely linked to gene activity governing key traits in plant development and adaptation. Particularly in maize, its dynamic chromatin structure is found to collaborate with massive transcriptional variations across tissues and developmental stages, implying intricate regulatory mechanisms, which highlights the importance of integrating chromatin information into breeding strategies for precise gene controls. The depiction of maize chromatin architecture using Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) provides great opportunities to investigate cis-regulatory elements, which is crucial for crop improvement. In this context, we developed an easy-to-implement ATAC-seq protocol for maize with fewer nuclei and simple equipment. We demonstrate a streamlined ATAC-seq protocol with four key steps for maize in which nuclei purification can be achieved without cell sorting and using only a standard bench-top centrifuge. Our protocol, coupled with the bioinformatic analysis, including validation by read length periodicity, key metrics, and correlation with transcript abundance, provides a precise and efficient assessment of the maize chromatin landscape. Beyond its application to maize, our testing design holds the potential to be applied to other crops or other tissues, especially for those with limited size and amount, establishing a robust foundation for chromatin structure studies in diverse crop species.

Human DNA topoisomerase I poisoning causes R loop-mediated genome instability attenuated by transcription factor IIS

DNA topoisomerase I can contribute to cancer genome instability. During catalytic activity, topoisomerase I forms a transient intermediate, topoisomerase I-DNA cleavage complex (Top1cc) to allow strand rotation and duplex relaxation, which can lead to elevated levels of DNA-RNA hybrids and micronuclei. To comprehend the underlying mechanisms, we have integrated genomic data of Top1cc-triggered hybrids and DNA double-strand breaks (DSBs) shortly after Top1cc induction, revealing that Top1ccs increase hybrid levels with different mechanisms. DSBs are at highly transcribed genes in early replicating initiation zones and overlap with hybrids downstream of accumulated RNA polymerase II (RNAPII) at gene 5'-ends. A transcription factor IIS mutant impairing transcription elongation further increased RNAPII accumulation likely due to backtracking. Moreover, Top1ccs can trigger micronuclei when occurring during late G1 or early/mid S, but not during late S. As micronuclei and transcription-replication conflicts are attenuated by transcription factor IIS, our results support a role of RNAPII arrest in Top1cc-induced transcription-replication conflicts leading to DSBs and micronuclei.

The balance between antiviral and antibacterial responses during M. tuberculosis infection is regulated by the ubiquitin ligase CBL

As a first line of host defense, macrophages must be able to effectively sense and respond to diverse types of pathogens, and while a particular type of immune response may be beneficial in some circumstances, it can be detrimental in others. Upon infecting a macrophage, M. tuberculosis (Mtb) induces proinflammatory cytokines that activate antibacterial responses. Surprisingly, Mtb also triggers antiviral responses that actually hinder the ability of macrophages to control Mtb infection. The ubiquitin ligase CBL suppresses these antiviral responses and shifts macrophages toward a more antibacterial state during Mtb infection, however, the mechanisms by which CBL regulates immune signaling are unknown. We found that CBL controls responses to multiple stimuli and broadly suppresses the expression of antiviral effector genes. We then used mass-spectrometry to investigate potential CBL substrates and identified over 46,000 ubiquitylated peptides in Mtb-infected macrophages, as well as roughly 400 peptides with CBL-dependent ubiquitylation. We then performed genetic interaction analysis of CBL and its putative substrates, and identified the Fas associated factor 2 (FAF2) adapter protein as a key signaling molecule protein downstream of CBL. Together, these analyses identify thousands of new ubiquitin-mediated signaling events during the innate immune response and reveal an important new regulatory hub in this response.

Principles of chromosome organization for meiotic recombination

In meiotic cells, chromosomes are organized as chromatin loop arrays anchored to a protein axis. This organization is essential to regulate meiotic recombination, from DNA double-strand break (DSB) formation to their repair. In mammals, it is unknown how chromatin loops are organized along the genome and how proteins participating in DSB formation are tethered to the chromosome axes. Here, we identify three categories of axis-associated genomic sites: PRDM9 binding sites, where DSBs form; binding sites of the insulator protein CTCF; and H3K4me3-enriched sites. We demonstrate that PRDM9 promotes the recruitment of MEI4 and IHO1, two proteins essential for DSB formation. In turn, IHO1 anchors DSB sites to the axis components HORMAD1 and SYCP3. We discovered that IHO1, HORMAD1, and SYCP3 are associated at the DSB ends during DSB repair. Our results highlight how interactions of proteins with specific genomic elements shape the meiotic chromosome organization for recombination.

Estrogen receptor 1 chromatin profiling in human breast tumors reveals high inter-patient heterogeneity with enrichment of risk SNPs and enhancer activity at most-conserved regions

Estrogen Receptor 1 (ESR1; also known as ERα, encoded by ESR1 gene) is the main driver and prime drug target in luminal breast cancer. ESR1 chromatin binding is extensively studied in cell lines and a limited number of human tumors, using consensi of peaks shared among samples. However, little is known about inter-tumor heterogeneity of ESR1 chromatin action, along with its biological implications. Here, we use a large set of ESR1 ChIP-seq data from 70 ESR1+ breast cancers to explore inter-patient heterogeneity in ESR1 DNA binding to reveal a striking inter-tumor heterogeneity of ESR1 action. Of note, commonly shared ESR1 sites show the highest estrogen-driven enhancer activity and are most engaged in long-range chromatin interactions. In addition, the most commonly shared ESR1-occupied enhancers are enriched for breast cancer risk SNP loci. We experimentally confirm SNVs to impact chromatin binding potential for ESR1 and its pioneer factor FOXA1. Finally, in the TCGA breast cancer cohort, we can confirm these variations to associate with differences in expression for the target gene. Cumulatively, we reveal a natural hierarchy of ESR1-chromatin interactions in breast cancers within a highly heterogeneous inter-tumor ESR1 landscape, with the most common shared regions being most active and affected by germline functional risk SNPs for breast cancer development.

Histone H3 posttranslational modified enzymes defined neutrophil plasticity and their vulnerability to IL-10 in the course of the inflammation

BACKGROUND

Neutrophils are a heterogeneous population capable of antimicrobial functions associated with pre-activation/activation and tissue regeneration. The specific polarisation of immune cells is mediated by the modification of 'chromatin landscapes', which enables differentiated access and activity of regulatory elements that guarantee their plasticity during inflammation No specific pattern within histone posttranslational modifications (PTMs) controlling this plasticity has been identified.

METHODS

Using the in vitro model of inflammation, reflecting different states of neutrophils from resting, pre-activated cells to activated and reducing tissue regeneration, we have analysed 11 different histone posttranslational modifications (PTMs), PTM enzymes associated with remodelling neutrophil chromatin, and H3K4me3 ChIP-Seq Gene Ontology analysis focusing on the processes related to histone PTMs. These findings were verified by extrapolation to adequate clinical status, using neutrophils derived from the patients with sepsis (systemic septic inflammation with LPS-stimulated neutrophils), neuromyelitis optical spectrum disorders (aseptic inflammation with pre-activated neutrophils) and periodontitis (local self-limiting septic inflammation with IL-10-positive neutrophils).

RESULTS

Physiological activation of neutrophils comprises a pre-activation characterised by histone H3K27ac and H3K4me1, which position enhancers; direct LPS exposure is induced explicitly by H3K4me3 which marked Transcription Start Site (TSS) regions and low-level of H3K9me3, H3K79me2 and H3K27me3 which, in turn, marked repressed genes. Contrary to antimicrobial action, IL-10 positively induced levels of H3S10p and negatively H3K9me3, which characterised processes related to the activation of genes within heterochromatin mediated by CHD1 and H3K9me3 specific demethylase JMJD2A. IL-10 protects changes within histone PTMs induced by TNF or LPS that affected H3K4me3-specific methyltransferase SETD1A and MLL1. Neutrophils previously exposed to inflammatory factors become unvulnerable to IL-10 because previous LPS stimulation interrupts TSS regions marked by H3K4me3 of CHD1 and JMJD2A genes. Therefore, LPS-activated neutrophils are disabled to induce CHD1/JMJD2A enzymes by IL-10, making this process irreversible. Because transcription of JMJD2A and CHD1 also depends on TSS positioning by H3K4me3, neutrophils before LPS stimulation become insensitive to IL-10.

CONCLUSION

Neutrophils, once pre-activated by TNF or directly stimulated by LPS, become insensitive to the anti-inflammatory effects of IL-10, and vice versa; IL-10 protects neutrophils against these proinflammatory stimuli. This phenomenon is responsible for disturbing the natural process of resolving inflammation and tissue regeneration.

Disregarding multimappers leads to biases in the functional assessment of NGS data

BACKGROUND

Standard ChIP-seq and RNA-seq processing pipelines typically disregard sequencing reads whose origin is ambiguous ("multimappers"). This usual practice has potentially important consequences for the functional interpretation of the data: genomic elements belonging to clusters composed of highly similar members are left unexplored.

RESULTS

In particular, disregarding multimappers leads to the underrepresentation in epigenetic studies of recently active transposable elements, such as AluYa5, L1HS and SVAs. Furthermore, this common strategy also has implications for transcriptomic analysis: members of repetitive gene families, such the ones including major histocompatibility complex (MHC) class I and II genes, are under-quantified.

CONCLUSION

Revealing inherent biases that permeate routine tasks such as functional enrichment analysis, our results underscore the urgency of broadly adopting multimapper-aware bioinformatic pipelines -currently restricted to specific contexts or communities- to ensure the reliability of genomic and transcriptomic studies.

Functional Divergence in Orthologous Transcription Factors: Insights from AtCBF2/3/1 and OsDREB1C

Despite traditional beliefs of orthologous genes maintaining similar functions across species, growing evidence points to their potential for functional divergence. C-repeat binding factors/dehydration-responsive element binding protein 1s (CBFs/DREB1s) are critical in cold acclimation, with their overexpression enhancing stress tolerance but often constraining plant growth. In contrast, a recent study unveiled a distinctive role of rice OsDREB1C in elevating nitrogen use efficiency (NUE), photosynthesis, and grain yield, implying functional divergence within the CBF/DREB1 orthologs across species. Here, we delve into divergent molecular mechanisms of OsDREB1C and AtCBF2/3/1 by exploring their evolutionary trajectories across rice and Arabidopsis genomes, regulatomes, and transcriptomes. Evolutionary scrutiny shows discrete clades for OsDREB1C and AtCBF2/3/1, with the Poaceae-specific DREB1C clade mediated by a transposon event. Genome-wide binding profiles highlight OsDREB1C's preference for GCCGAC compared to AtCBF2/3/1's preference for A/GCCGAC, a distinction determined by R12 in the OsDREB1C AP2/ERF domain. Cross-species multiomic analyses reveal shared gene orthogroups (OGs) and underscore numerous specific OGs uniquely bound and regulated by OsDREB1C, implicated in NUE, photosynthesis, and early flowering, or by AtCBF2/3/1, engaged in hormone and stress responses. This divergence arises from gene gains/losses (∼16.7% to 25.6%) and expression reprogramming (∼62.3% to 66.2%) of OsDREB1C- and AtCBF2/3/1-regulated OGs during the extensive evolution following the rice-Arabidopsis split. Our findings illustrate the regulatory evolution of OsDREB1C and AtCBF2/3/1 at a genomic scale, providing insights on the functional divergence of orthologous transcription factors following gene duplications across species.

Evolutionary Dynamics of Chromatin Structure and Duplicate Gene Expression in Diploid and Allopolyploid Cotton

Polyploidy is a prominent mechanism of plant speciation and adaptation, yet the mechanistic understandings of duplicated gene regulation remain elusive. Chromatin structure dynamics are suggested to govern gene regulatory control. Here, we characterized genome-wide nucleosome organization and chromatin accessibility in allotetraploid cotton, Gossypium hirsutum (AADD, 2n = 4X = 52), relative to its two diploid parents (AA or DD genome) and their synthetic diploid hybrid (AD), using DNS-seq. The larger A-genome exhibited wider average nucleosome spacing in diploids, and this intergenomic difference diminished in the allopolyploid but not hybrid. Allopolyploidization also exhibited increased accessibility at promoters genome-wide and synchronized cis-regulatory motifs between subgenomes. A prominent cis-acting control was inferred for chromatin dynamics and demonstrated by transposable element removal from promoters. Linking accessibility to gene expression patterns, we found distinct regulatory effects for hybridization and later allopolyploid stages, including nuanced establishment of homoeolog expression bias and expression level dominance. Histone gene expression and nucleosome organization are coordinated through chromatin accessibility. Our study demonstrates the capability to track high-resolution chromatin structure dynamics and reveals their role in the evolution of cis-regulatory landscapes and duplicate gene expression in polyploids, illuminating regulatory ties to subgenomic asymmetry and dominance.

Homeobox regulator Wilms Tumour 1 is displaced by androgen receptor at cis-regulatory elements in the endometrium of PCOS patients

Decidualisation, the process whereby endometrial stromal cells undergo morphological and functional transformation in preparation for trophoblast invasion, is often disrupted in women with polycystic ovary syndrome (PCOS) resulting in complications with pregnancy and/or infertility. The transcription factor Wilms tumour suppressor 1 (WT1) is a key regulator of the decidualization process, which is reduced in patients with PCOS, a complex condition characterized by increased expression of androgen receptor in endometrial cells and high presence of circulating androgens. Using genome-wide chromatin immunoprecipitation approaches on primary human endometrial stromal cells, we identify key genes regulated by WT1 during decidualization, including homeobox transcription factors which are important for regulating cell differentiation. Furthermore, we found that AR in PCOS patients binds to the same DNA regions as WT1 in samples from healthy endometrium, suggesting dysregulation of genes important to decidualisation pathways in PCOS endometrium due to competitive binding between WT1 and AR. Integrating RNA-seq and H3K4me3 and H3K27ac ChIP-seq metadata with our WT1/AR data, we identified a number of key genes involved in immune response and angiogenesis pathways that are dysregulated in PCOS patients. This is likely due to epigenetic alterations at distal enhancer regions allowing AR to recruit cofactors such as MAGEA11, and demonstrates the consequences of AR disruption of WT1 in PCOS endometrium.

Robust chromatin state annotation

With the goal of mapping genomic activity, international projects have recently measured epigenetic activity in hundreds of cell and tissue types. Chromatin state annotations produced by segmentation and genome annotation (SAGA) methods have emerged as the predominant way to summarize these epigenomic data sets in order to annotate the genome. These chromatin state annotations are essential for many genomic tasks, including identifying active regulatory elements and interpreting disease-associated genetic variation. However, despite the widespread applications of SAGA methods, no principled approach exists to evaluate the statistical significance of chromatin state assignments. Here, we propose the first method for assigning calibrated confidence scores to chromatin state annotations. Toward this goal, we performed a comprehensive evaluation of the reproducibility of the two most widely used existing SAGA methods, ChromHMM and Segway. We found that their predictions are frequently irreproducible. For example, when applying the same SAGA method on two sets of experimental replicates, 27%-69% of predicted enhancers fail to replicate. This suggests that a substantial fraction of predicted elements in existing chromatin state annotations cannot be relied upon. To remedy this problem, we introduce SAGAconf, a method for assigning a measure of confidence (r-value) to chromatin state annotations. SAGAconf works with any SAGA method and assigns an r-value to each genomic bin of a chromatin state annotation that represents the probability that the label of this bin will be reproduced in a replicated experiment. Thus, SAGAconf allows a researcher to select only the reliable predictions from a chromatin annotation for use in downstream analyses.

Exhaustive identification of genome-wide binding events of transcriptional regulators

Genome-wide binding assays aspire to map the complete binding pattern of gene regulators. Common practice relies on replication-duplicates or triplicates-and high stringency statistics to favor false negatives over false positives. Here we show that duplicates and triplicates of CUT&RUN are not sufficient to discover the entire activity of transcriptional regulators. We introduce ICEBERG (Increased Capture of Enrichment By Exhaustive Replicate aGgregation), a pipeline that harnesses large numbers of CUT&RUN replicates to discover the full set of binding events and chart the line between false positives and false negatives. We employed ICEBERG to map the full set of H3K4me3-marked regions, the targets of the co-factor β-catenin, and those of the transcription factor TBX3, in human colorectal cancer cells. The ICEBERG datasets allow benchmarking of individual replicates, comparing the performance of peak calling and replication approaches, and expose the arbitrary nature of strategies to identify reproducible peaks. Instead of a static view of genomic targets, ICEBERG establishes a spectrum of detection probabilities across the genome for a given factor, underlying the intrinsic dynamicity of its mechanism of action, and permitting to distinguish frequent from rare regulation events. Finally, ICEBERG discovered instances, undetectable with other approaches, that underlie novel mechanisms of colorectal cancer progression.

TAD boundary deletion causes PITX2-related cardiac electrical and structural defects

While 3D chromatin organization in topologically associating domains (TADs) and loops mediating regulatory element-promoter interactions is crucial for tissue-specific gene regulation, the extent of their involvement in human Mendelian disease is largely unknown. Here, we identify 7 families presenting a new cardiac entity associated with a heterozygous deletion of 2 CTCF binding sites on 4q25, inducing TAD fusion and chromatin conformation remodeling. The CTCF binding sites are located in a gene desert at 1 Mb from the Paired-like homeodomain transcription factor 2 gene (PITX2). By introducing the ortholog of the human deletion in the mouse genome, we recapitulate the patient phenotype and characterize an opposite dysregulation of PITX2 expression in the sinoatrial node (ectopic activation) and ventricle (reduction), respectively. Chromatin conformation assay performed in human induced pluripotent stem cell-derived cardiomyocytes harboring the minimal deletion identified in family#1 reveals a conformation remodeling and fusion of TADs. We conclude that TAD remodeling mediated by deletion of CTCF binding sites causes a new autosomal dominant Mendelian cardiac disorder.

Relaxin Modulates the Genomic Actions and Biological Effects of Estrogen in the Myometrium by Reducing Estrogen Receptor Alpha Phosphorylation

Estradiol (E2) and relaxin (Rln) are steroid and polypeptide hormones, respectively, with important roles in the female reproductive tract, including myometrium. Some actions of Rln, which are mediated by its membrane receptor RXFP1, require or are augmented by E2 signaling through its cognate nuclear steroid receptor, estrogen receptor alpha (ERα). In contrast, other actions of Rln act in opposition to the effects of E2. Here we explore the molecular and genomic mechanisms that underlie the functional interplay between E2 and Rln in the myometrium. We used both ovariectomized female mice and immortalized human myometrial cells expressing wild type or mutant ERα (hTERT-HM-ERα cells). Our results indicate that Rln attenuates the genomic actions and biological effects of estrogen in the myometrium and myometrial cells by reducing phosphorylation ERα on serine 118 (S118). Interestingly, we observed a potent inhibitory effect of Rln on the E2-dependent binding of ERα across the genome. The reduction in ERα binding was associated with changes in the hormone-regulated transcriptome, including a decrease in the E2-dependent expression of neighboring genes. The inhibitory effects of Rln cotreatment on the E2-dependent phosphorylation of ERα required the nuclear dual-specificity phosphatases DUSP1 and DUSP5. Moreover, the inhibitory effects of Rln were reflected in a concomitant inhibition of the E2-dependent contraction of myometrial cells. Collectively, our results identify a pathway that integrates Rln/RXFP1 and E2/ERα signaling, resulting in a convergence of membrane and nuclear signaling pathways to control genomic and biological outcomes.

Sex differences in brain cell-type specific chromatin accessibility in schizophrenia

Our understanding of the sex-specific role of the non-coding genome in serious mental illness remains largely incomplete. To address this gap, we explored sex differences in 1,393 chromatin accessibility profiles, derived from neuronal and non-neuronal nuclei of two distinct cortical regions from 234 cases with serious mental illness and 235 controls. We identified sex-specific enhancer-promoter interactions and showed that they regulate genes involved in X-chromosome inactivation (XCI). Examining chromosomal conformation allowed us to identify sex-specific cis- and trans-regulatory domains (CRDs and TRDs). Co-localization of sex-specific TRDs with schizophrenia common risk variants pinpointed male-specific regulatory regions controlling a number of metabolic pathways. Additionally, enhancers from female-specific TRDs were found to regulate two genes known to escape XCI, (XIST and JPX), underlying the importance of TRDs in deciphering sex differences in schizophrenia. Overall, these findings provide extensive characterization of sex differences in the brain epigenome and disease-associated regulomes.

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.