Transcription Factor Binding Research Articles

Characterization of single-cell cis-regulatory elements informs implications for cell differentiation.

Cis-regulatory elements (CREs) govern the specific patterns and dynamics of gene expression in cells during development, which are the fundamental mechanisms behind cell differentiation. However, the genomic characteristics of single-cell CREs (scCREs) closely linked to cell differentiation during development remain unclear. To explore this, we systematically analyzed ∼250,000 putative scCREs obtained from snATAC-seq analysis of the developing mouse cerebellum. We found that over 80% of these scCREs show pleiotropic effects, being active in two or more cell types. The pleiotropic degrees of proximal and distal scCREs are positively correlated with the density and diversity of transcription factor (TF) binding motifs and GC content. There is a negative correlation between the pleiotropic degrees of scCREs and their distances to the nearest TSSs, and proximal scCREs display higher relevance strengths than distal ones. Furthermore, both proximal and distal scCREs related to cell differentiation exhibit enhanced sequence-level evolutionary conservation, increased density and diversity of TF binding motifs, elevated GC content, and greater distances from their nearest genes. Together, our findings reveal the general genomic characteristics of putative scCREs and provide insights into the genomic and evolutionary mechanisms by which scCREs regulate cell differentiation during development.

In Vitro Genome-Wide Identification of Transcription Factor Binding Sites.

DNA affinity purification followed by sequencing (DAP-seq) profiles transcription factor (TF) binding sites in vitro, in which recombinant TF captures the genomic fragments with TF binding sites. Here, I describe the preparation of pre-PCR DNA-seq libraries for 96 TFs, including the expression and purification of recombinant TF, the affinity purification of TF binding fragments, and the sequencing-ready library preparation. Using this method, up to 96 DAP-seq libraries can be prepared in 2days.

Dynamic atlas of histone modifications and gene regulatory networks in endosperm of bread wheat.

Dissecting the genetic basis of seed traits in wheat is impeded by limited genetic polymorphisms and significant variations caused by environmental conditions and seed position in a spikelet. Seed performance is largely determined by endosperm development controlled by spatiotemporal variation in gene activities, which is greatly affected by chromatin status. Here, we map genome-wide dynamic distributions of H3K27me3, H3K4me3 and H3K9ac modifications and profile gene transcription across wheat endosperm development. The combinatorial effects of active and repressive marks ensure spatiotemporal dynamic gene expression, especially for starch biosynthesis. By scanning the transcription factor binding motifs in the ATAC-seq peaks, hub regulators are identified from the regulatory network. In addition, we observe significant correlations between sequence polymorphisms of hub regulators and variations in seed traits in a germplasm population. Thus, the analysis of genomic regulatory activities together with genetic variation provides a robust approach to dissect seed traits in bread wheat.

Integrative analysis based on ATAC-seq and RNA-seq reveals a novel oncogene PRPF3 in hepatocellular carcinoma.

Assay of Transposase Accessible Chromatin Sequencing (ATAC-seq) is a high-throughput sequencing technique that detects open chromatin regions across the genome. These regions are critical in facilitating transcription factor binding and subsequent gene expression. Herein, we utilized ATAC-seq to identify key molecular targets regulating the development and progression of hepatocellular carcinoma (HCC) and elucidate the underlying mechanisms. We first compared chromatin accessibility profiles between HCC and normal tissues. Subsequently, RNA-seq data was employed to identify differentially expressed genes (DEGs). Integrating ATAC-seq and RNA-seq data allowed the identification of transcription factors and their putative target genes associated with differentially accessible regions (DARs). Finally, functional experiments were conducted to investigate the effects of the identified regulatory factors and corresponding targets on HCC cell proliferation and migration. Enrichment analysis of DARs between HCC and adjacent normal tissues revealed distinct signaling pathways and regulatory factors. Upregulated DARs in HCC were enriched in genes related to the MAPK and FoxO signaling pathways and associated with transcription factor families like ETS and AP-1. Conversely, downregulated DARs were associated with the TGF-β, cAMP, and p53 signaling pathways and the CTCF family. Integration of the datasets revealed a positive correlation between specific DARs and DEGs. Notably, PRPF3 emerged as a gene associated with DARs in HCC, and functional assays demonstrated its ability to promote HCC cell proliferation and migration. To the best of our knowledge, this is the first report highlighting the oncogenic role of PRPF3 in HCC. Furthermore, ZNF93 expression positively correlated with PRPF3, and ChIP-seq data indicated its potential role as a transcription factor regulating PRPF3 by binding to its promoter region. This study provides a comprehensive analysis of the epigenetic landscape in HCC, encompassing both chromatin accessibility and the transcriptome. Our findings reveal that ZNF93 promotes the proliferation and motility of HCC cells through transcriptional regulation of a novel oncogene, PRPF3.

A Foxf1-Wnt-Nr2f1 cascade promotes atrial cardiomyocyte differentiation in zebrafish.

Nr2f transcription factors (TFs) are conserved regulators of vertebrate atrial cardiomyocyte (AC) differentiation. However, little is known about the mechanisms directing Nr2f expression in ACs. Here, we identified a conserved enhancer 3' to the nr2f1a locus, which we call 3'reg1-nr2f1a (3'reg1), that can promote Nr2f1a expression in ACs. Sequence analysis of the enhancer identified putative Lef/Tcf and Foxf TF binding sites. Mutation of the Lef/Tcf sites within the 3'reg1 reporter, knockdown of Tcf7l1a, and manipulation of canonical Wnt signaling support that Tcf7l1a is derepressed via Wnt signaling to activate the transgenic enhancer and promote AC differentiation. Similarly, mutation of the Foxf binding sites in the 3'reg1 reporter, coupled with gain- and loss-of-function analysis supported that Foxf1 promotes expression of the enhancer and AC differentiation. Functionally, we find that Wnt signaling acts downstream of Foxf1 to promote expression of the 3'reg1 reporter within ACs and, importantly, both Foxf1 and Wnt signaling require Nr2f1a to promote a surplus of differentiated ACs. CRISPR-mediated deletion of the endogenous 3'reg1 abrogates the ability of Foxf1 and Wnt signaling to produce surplus ACs in zebrafish embryos. Together, our data support that downstream members of a conserved regulatory network involving Wnt signaling and Foxf1 function on a nr2f1a enhancer to promote AC differentiation in the zebrafish heart.

Decoding the epigenetics and chromatin loop dynamics of androgen receptor-mediated transcription

Androgen receptor (AR)-mediated transcription plays a critical role in development and prostate cancer growth. AR drives gene expression by binding to thousands of cis-regulatory elements (CRE) that loop to hundreds of target promoters. With multiple CREs interacting with a single promoter, it remains unclear how individual AR bound CREs contribute to gene expression. To characterize the involvement of these CREs, we investigate the AR-driven epigenetic and chromosomal chromatin looping changes by generating a kinetic multi-omic dataset comprised of steady-state mRNA, chromatin accessibility, transcription factor binding, histone modifications, chromatin looping, and nascent RNA. Using an integrated regulatory network, we find that AR binding induces sequential changes in the epigenetic features at CREs, independent of gene expression. Further, we show that binding of AR does not result in a substantial rewiring of chromatin loops, but instead increases the contact frequency of pre-existing loops to target promoters. Our results show that gene expression strongly correlates to the changes in contact frequency. We then propose and experimentally validate an unbalanced multi-enhancer model where the impact on gene expression of AR-bound enhancers is heterogeneous, and is proportional to their contact frequency with target gene promoters. Overall, these findings provide insights into AR-mediated gene expression upon acute androgen simulation and develop a mechanistic framework to investigate nuclear receptor mediated perturbations.

A DNA base-specific sequence interposed between CRX and NRL contributes to RHODOPSIN expression

Gene expression emerges from DNA sequences through the interaction of transcription factors (TFs) with DNA cis-regulatory sequences. In eukaryotes, TFs bind to transcription factor binding sites (TFBSs) with differential affinities, enabling cell-specific gene expression. In this view, DNA enables TF binding along a continuum ranging from low to high affinity depending on its sequence composition; however, it is not known whether evolution has entailed a further level of entanglement between DNA-protein interaction. Here we found that the composition and length (22 bp) of the DNA sequence interposed between the CRX and NRL retinal TFs in the proximal promoter of RHODOPSIN (RHO) largely controls the expression levels of RHO. Mutagenesis of CRX-NRL DNA linking sequences (here termed “DNA-linker”) results in uncorrelated gene expression variation. In contrast, mutual exchange of naturally occurring divergent human and mouse Rho cis-regulatory elements conferred similar yet species-specific Rho expression levels. Two orthogonal DNA-binding proteins targeted to the DNA-linker either activate or repress the expression of Rho depending on the DNA-linker orientation relative to the CRX and NRL binding sites. These results argue that, in this instance, DNA itself contributes to CRX and NRL activities through a code based on specific base sequences of a defined length, ultimately determining optimal RHO expression levels.

Expression interplay of genes coding for calcium-binding proteins and transcription factors during the osmotic phase provides insights on salt stress response mechanisms in bread wheat

Bread wheat is an important crop for the human diet, but the increasing soil salinization is reducing the yield. The Ca2+ signaling events at the early stages of the osmotic phase of salt stress are crucial for the acclimation response of the plants through the performance of calcium-sensing proteins, which activate or repress transcription factors (TFs) that affect the expression of downstream genes. Physiological, genetic mapping, and transcriptomics studies performed with the contrasting genotypes Syn86 (synthetic, salt-susceptible) and Zentos (elite cultivar, salt-tolerant) were integrated to gain a comprehensive understanding of the salt stress response. The MACE (Massive Analysis of cDNA 3ʹ-Ends) based transcriptome analysis until 4 h after stress exposure revealed among the salt-responsive genes, the over-representation of genes coding for calcium-binding proteins. The functional and structural diversity within this category was studied and linked with the expression levels during the osmotic phase in the contrasting genotypes. The non-EF-hand category from calcium-binding proteins was found to be enriched for the susceptibility response. On the other side, the tolerant genotype was characterized by a faster and higher up-regulation of genes coding for proteins with EF-hand domain, such as RBOHD orthologs, and TF members. This study suggests that the interplay of calcium-binding proteins, WRKY, and AP2/ERF TF families in signaling pathways at the start of the osmotic phase can affect the expression of downstream genes. The identification of SNPs in promoter sequences and 3ʹ -UTR regions provides insights into the molecular mechanisms controlling the differential expression of these genes through differential transcription factor binding affinity or altered mRNA stability.

Human iPSC-derived neural stem cells displaying radial glia signature exhibit long-term safety in mice

Human induced pluripotent stem cell-derived neural stem/progenitor cells (hiPSC-NSCs) hold promise for treating neurodegenerative and demyelinating disorders. However, comprehensive studies on their identity and safety remain limited. In this study, we demonstrate that hiPSC-NSCs adopt a radial glia-associated signature, sharing key epigenetic and transcriptional characteristics with human fetal neural stem cells (hfNSCs) while exhibiting divergent profiles from glioblastoma stem cells. Long-term transplantation studies in mice showed robust and stable engraftment of hiPSC-NSCs, with predominant differentiation into glial cells and no evidence of tumor formation. Additionally, we identified the Sterol Regulatory Element Binding Transcription Factor 1 (SREBF1) as a regulator of astroglial differentiation in hiPSC-NSCs. These findings provide valuable transcriptional and epigenetic reference datasets to prospectively define the maturation stage of NSCs derived from different hiPSC sources and demonstrate the long-term safety of hiPSC-NSCs, reinforcing their potential as a viable alternative to hfNSCs for clinical applications.

Distinct H3K9me3 heterochromatin maintenance dynamics govern different gene programmes and repeats in pluripotent cells.

H3K9me3 heterochromatin, established by lysine methyltransferases (KMTs) and compacted by heterochromatin protein 1 (HP1) isoforms, represses alternative lineage genes and DNA repeats. Our understanding of H3K9me3 heterochromatin stability is presently limited to individual domains and DNA repeats. Here we engineered Suv39h2-knockout mouse embryonic stem cells to degrade remaining two H3K9me3 KMTs within 1 hour and found that both passive dilution and active removal contribute to H3K9me3 decay within 12-24 hours. We discovered four different H3K9me3 decay rates across the genome and chromatin features and transcription factor binding patterns that predict the stability classes. A 'binary switch' governs heterochromatin compaction, with HP1 rapidly dissociating from heterochromatin upon KMT depletion and a particular threshold level of HP1 limiting pioneer factor binding, chromatin opening and exit from pluripotency within 12 h. Unexpectedly, receding H3K9me3 domains unearth residual HP1β peaks enriched with heterochromatin-inducing proteins. Our findings reveal distinct H3K9me3 heterochromatin maintenance dynamics governing gene networks and repeats that together safeguard pluripotency.

Benchmarking and building DNA binding affinity models using allele-specific and allele-agnostic transcription factor binding data

BackgroundTranscription factors (TFs) bind to DNA in a highly sequence-specific manner. This specificity manifests itself in vivo as differences in TF occupancy between the two alleles at heterozygous loci. Genome-scale assays such as ChIP-seq currently are limited in their power to detect allele-specific binding (ASB) both in terms of read coverage and representation of individual variants in the cell lines used. This makes prediction of allelic differences in TF binding from sequence alone desirable, provided that the reliability of such predictions can be quantitatively assessed.ResultsWe here propose methods for benchmarking sequence-to-affinity models for TF binding in terms of their ability to predict allelic imbalances in ChIP-seq counts. We use a likelihood function based on an over-dispersed binomial distribution to aggregate evidence for allelic preference across the genome without requiring statistical significance for individual variants. This allows us to systematically compare predictive performance when multiple binding models for the same TF are available. To facilitate the de novo inference of high-quality models from paired-end in vivo binding data such as ChIP-seq, ChIP-exo, and CUT&Tag without read mapping or peak calling, we introduce an extensible reimplementation of our biophysically interpretable machine learning framework named PyProBound. Explicitly accounting for assay-specific bias in DNA fragmentation rate when training on ChIP-seq yields improved TF binding models. Moreover, we show how PyProBound can leverage our threshold-free ASB likelihood function to perform de novo motif discovery using allele-specific ChIP-seq counts.ConclusionOur work provides new strategies for predicting the functional impact of non-coding variants.

Micro-capture-C in vessel wall cell types identifies mechanisms underlying coronary artery disease susceptibility loci

Abstract Genetics is estimated to account for 50% of Coronary Artery Disease (CAD) risk. Genome-Wide Association Studies (GWAS) have revealed 100s of genomic loci which are associated with CAD (1). However, moving from GWAS signals to causal genes has proved challenging, slowing progress in developing interventions for CAD and other diseases. The greatest challenge in interpreting GWAS data is that only 10% of variants are in protein-coding regions and 90% are in "non-coding" genomic regions. Many of these non-coding variants are found in regulatory elements such as enhancers. Enhancers are distal regulatory elements that bind transcription factors and regulate genes up to 1 million base pairs away. This is achieved by DNA looping bringing enhancers into physical proximity with target genes. However, it is not possible to predict which gene(s) a given enhancer will regulate based on genomic location alone. We aimed to develop a generalisable workflow to solve CAD GWAS signals accurately and sensitively, by identifying causal enhancer SNPs and linking them mechanistically to output genes in a cell-type specific manner. We combined cell-specific epigenomic mapping with machine learning to prioritise 21,959 SNPs from 861 candidate CAD loci, in endothelial cells (HUVECs) and smooth muscle cells (HCASMCs). Micro-Capture-C (MCC) was used to resolve enhancer-promoter interactions at base-pair resolution (2). MNase footprinting determined differential transcription factor binding at specific enhancer SNPs. In HCASMCs, 20 putative enhancers containing CAD SNPs were found to interact with 32 genes, and in HUVECs 25 enhancers interacted with 48 genes. These genes were enriched for pathways including TGFb signalling, RUNX3 activity, and MAPK6/MAPK4 signalling. To prove causality at specific loci, heterozygous SNPs were identified in the individual primary cell lines. The interaction frequency at these elements was statistically compared between risk and non-risk alleles. For example, at the 9p21.3 locus, an enhancer containing a non-risk T allele was found to interact more strongly with the CDKN2B gene, compared to the risk G allele. Interestingly, this imbalance between allele interaction was potentiated by stimulation with TGFb, a key mediator in atherosclerosis. Furthermore, the transcription factor footprint at the SNP location was altered by the risk allele, indicating the mechanistic cause of the CAD-associated change in CDKN2B expression. We demonstrate a combined computational and experimental approach to identify causal non-coding variants in CAD risk loci, at scale. Among the key outputs of this study is the identification of the risk variant (rs1537373), output gene (CDKN2B), and specific cell type (HCASMC) underlying the 9p21.3 locus association with CAD risk. Using this approach to translate GWAS signals into cell-type specific outputs genes, will lead to a better understanding of the pathogenesis of CAD, and new therapeutic opportunities.

Integrative 3D genomics with multi-omics analysis and functional validation of genetic regulatory mechanisms of abdominal fat deposition in chickens

Chickens are the most abundant agricultural animals globally, with controlling abdominal fat deposition being a key objective in poultry breeding. While GWAS can identify genetic variants associated with abdominal fat deposition, the precise roles and mechanisms of these variants remain largely unclear. Here, we use male chickens from two lines divergently selected for abdominal fat deposition as experimental models. Through the integration of genomic, epigenomic, 3D genomic, and transcriptomic data, we build a comprehensive chromatin 3D regulatory network map to identify the genetic regulatory mechanisms that influence abdominal fat deposition in chickens. Notably, we find that the rs734209466 variant functions as an allele-specific enhancer, remotely enhancing the transcription of IGFBP2 and IGFBP5 by the binding transcription factor IRF4. This interaction influences the differentiation and proliferation of preadipocytes, which ultimately affects phenotype. This work presents a detailed genetic regulatory map for chicken abdominal fat deposition, offering molecular targets for selective breeding.

ALDH2 mediates the effects of sodium-glucose cotransporter 2 inhibitors (SGLT2i) on improving cardiac remodeling

BackgroundSodium-glucose cotransporter-2 inhibitors (SGLT2i) are now recommended for patients with heart failure, but the mechanisms that underlie the protective role of SGLT2i in cardiac remodeling remain unclear. Aldehyde dehydrogenase 2 (ALDH2) effectively prevents cardiac remodeling. Here, the key role of ALDH2 in the efficacy of SGLT2i on cardiac remodeling was studied.MethodsAnalysis of multiple transcriptomic datasets and two-sample Mendelian randomization were performed to find out the differentially expressed genes between pathological cardiac hypertrophy models (patients) and controls. A pathological cardiac hypertrophy mouse model was established via transverse aortic constriction (TAC) or isoproterenol (ISO). Cardiomyocyte-specific ALDH2 knockout mice (ALDH2CMKO) and littermate control mice (ALDH2flox/flox) were generated to determine the critical role of ALDH2 in the preventive effects of dapagliflozin (DAPA) on cardiac remodeling. RNA sequencing, gene knockdown or overexpression, bisulfite sequencing PCR, and luciferase reporter assays were performed to explore the underlying molecular mechanisms involved.ResultsOnly ALDH2 was differentially expressed when the differentially expressed genes obtained via Mendelian analysis and the differentially expressed genes obtained from the multiple transcriptome datasets were combined. Mendelian analysis revealed that ALDH2 was negatively related to the severity of myocardial hypertrophy in patients. DAPA alleviated cardiac remodeling in mouse hearts subjected to TAC or ISO. ALDH2 expression was reduced, whereas ALDH2 expression was restored by DAPA in hypertrophic hearts. Cardiomyocyte specific ALDH2 knockout abolished the protective role of DAPA in preventing cardiac remodeling. ALDH2 expression and activity were increased in DAPA-treated neonatal rat primary cardiomyocytes (NRCMs), H9C2 cells and AC16 cells. Moreover, DAPA upregulated ALDH2 in peripheral blood mononuclear cells (PBMCs) from patients with type 2 diabetes. Sodium/proton exchanger 1 (NHE1) inhibition contributed to the regulation of ALDH2 by DAPA. DAPA suppressed the production of reactive oxygen species (ROS), downregulated DNA methyltransferase 1 (DNMT1) and subsequently reduced the ALDH2 promoter methylation level. Further studies revealed that DAPA enhanced the binding of nuclear transcription factor Y, subunit A (NFYA) to the promoter region of ALDH2, which was due to the decreased promoter methylation level of ALDH2.ConclusionsThe upregulation of ALDH2 plays a critical role in the protection of DAPA against cardiac remodeling. DAPA enhances the binding of NFYA to the ALDH2 promoter by reducing the ALDH2 promoter methylation level through NHE1/ROS/DNMT1 pathway.Graphical abstract

Comparative Transcriptome Analysis of Non-Organogenic and Organogenic Tissues of Gaillardia pulchella Revealing Genes Regulating De Novo Shoot Organogenesis

Gaillardia pulchella is an important plant species with pharmacological and ornamental applications. It contains a wide array of phytocompounds which play roles against diseases. In vitro propagation requires callogenesis and differentiation of plant organs, which offers a sustainable, alternative synthesis of compounds. The morphogenetic processes and the underlying mechanisms are, however, known to be under genetic regulation and are little understood. The present study investigated these events by generating transcriptome data, with de novo assembly of sequences to describe shoot morphogenesis molecularly in G. pulchella. The RNA was extracted from the callus of pre- and post-shoot organogenesis time. The callus induction was optimal using leaf segments cultured onto MS medium containing α-naphthalene acetic acid (NAA; 2.0 mg/L) and 6-benzylaminopurine (BAP; 0.5 mg/L) and further exhibited a high shoot regeneration/caulogenesis ability. A total of 68,366 coding sequences were obtained using Illumina150bpPE sequencing and transcriptome assembly. Differences in gene expression patterns were noted in the studied samples, showing opposite morphogenetic responses. Out of 10,108 genes, 5374 (53%) were downregulated, and there were 4734 upregulated genes, representing 47% of the total genes. Through the heatmap, the top 100 up- and downregulating genes’ names were identified and presented. The up- and downregulated genes were identified using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway. Important pathways, operative during G. pulchella shoot organogenesis, were signal transduction (13.55%), carbohydrate metabolism (8.68%), amino acid metabolism (5.11%), lipid metabolism (3.75%), and energy metabolism (3.39%). The synthesized proteins displayed phosphorylation, defense response, translation, regulation of DNA-templated transcription, carbohydrate metabolic processes, and methylation activities. The genes’ product also exhibited ATP binding, DNA binding, metal ion binding, protein serine/threonine kinase -, ATP hydrolysis activity, RNA binding, protein kinase, heme and GTP binding, and DNA binding transcription factor activity. The most abundant proteins were located in the membrane, nucleus, cytoplasm, ribosome, ribonucleoprotein complex, chloroplast, endoplasmic reticulum membrane, mitochondrion, nucleosome, Golgi membrane, and other organellar membranes. These findings provide information for the concept of molecular triggers, regulating programming, differentiation and reprogramming of cells, and their uses.

Binding profiles for 961 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 efforts of the Model Organism ENCyclopedia Of DNA Elements (modENCODE) and the model organism Encyclopedia of Regulatory Networks (modERN) consortia to systematically assay TF binding events in vivo in two major model organisms, Drosophila melanogaster (fly) and Caenorhabditis elegans (worm). These datasets comprise 605 TFs identifying 3.6M sites in the fly and 356 TFs identifying 0.9 M sites in the worm and represent the majority of the regulatory space in each genome. We demonstrate that TFs associate with chromatin in clusters termed "metapeaks", that larger metapeaks have characteristics of high occupancy target (HOT) regions, and that the importance of consensus sequence motifs bound by TFs depends on metapeak size and complexity. Combining ChIP-seq data with single cell RNA-seq data in a machine learning model identifies TFs with a prominent role in promoting target gene expression in specific cell types, even differentiating between parent-daughter cells during embryogenesis. These data are a rich resource for the community that should fuel and guide future investigations into TF function. To facilitate data accessibility and utility, all strains expressing GFP-tagged TFs are available at the stock centers for each organism. The chromatin immunoprecipitation sequencing data are available through the ENCODE Data Coordinating Center, GEO, and through a direct interface that provides rapid access to processed data sets and summary analyses, as well as widgets to probe the cell type-specific TF-target relationships.

Human Satellite 3 DNA encodes megabase-scale transcription factor binding platforms.

Eukaryotic genomes are frequently littered with large arrays of tandem repeats, called satellite DNA, which underlie the constitutive heterochromatin often found around centromeric regions. While some satellite DNA types have well-established roles in centromere biology, other abundant satellite DNAs have poorly characterized functions. For example, Human Satellite 3 (HSat3), which makes up roughly 2% of the human genome, forms enormous arrays up to tens of megabases, but these arrays play no known roles in centromere function and were almost fully excluded from genome assemblies until recently. As a result, these massive genomic regions have remained relatively understudied, and the potential functional roles of HSat3 have remained largely unknown. To address this, we performed a systematic screen for novel HSat3 binding factors. Our work revealed HSat3 arrays contain high densities of transcription factor (TF) motifs that are bound by factors related to multiple, highly conserved signaling pathways. Unexpectedly, the most enriched TFs in HSat3 belong to the Hippo pathway transcription effector family TEAD. We found that TEAD recruits the co-activator YAP to HSat3 regions in a cell-state specific manner. Using RNA polymerase-I reporter assays, targeted repression of HSat3, inducible degradation of YAP, and super-resolution microscopy, we show that HSat3 arrays can localize YAP/TEAD inside the nucleolus, where YAP regulates RNA Polymerase-I activity. Beyond revealing a direct relationship between the Hippo pathway and ribosomal DNA regulation, this work demonstrates that satellite DNA can encode multiple transcription factor binding motifs, defining a new role for these enormous genomic elements.

METTL16 inhibits pancreatic cancer proliferation and metastasis by promoting MROH8 RNA stability and inhibiting CAPN2 expression.

N6-methyladenosine (m6A) modification plays a crucial role in the progression of various cancers, including pancreatic cancer, by regulating gene expression. However, the specific mechanisms by which m6A affects pancreatic cancer metastasis remain unclear. This study aims to elucidate the role of METTL16, an m6A writer gene, in regulating core genes such as CAPN2 and MROH8, influencing tumor growth and metastasis. Transcriptomic data from pancreatic cancer patients in The Cancer Genome Atlas (TCGA) were analyzed to identify m6A-related genes. We performed correlation and survival analyses to uncover core genes influenced by m6A expression. Functional assays, including METTL16 knockdown and overexpression experiments, were conducted in pancreatic cancer cell lines, patient-derived organoids, and animal models. Immunofluorescence, co-immunoprecipitation (Co-IP), and m6A-specific quantitative PCR were used to validate protein interactions and m6A modifications. Chromatin immunoprecipitation (ChIP) analysis was utilized to investigate transcription factor binding at gene promoter regions. METTL16 and METTL3 were identified as key m6A regulators associated with improved prognosis in pancreatic cancer patients (P<0.05). CAPN2, CHMP2B, ITGA3, ITGA6, ITPR1, and RAC1 were identified as core genes linked to m6A expression, all significantly correlated with patient prognosis (P<0.05). METTL16 overexpression significantly inhibited tumor growth and metastasis (P<0.001) by downregulating CAPN2 through an indirect mechanism involving the transcription factor TBP and the gene MROH8. MROH8 negatively regulated CAPN2 by promoting TBP degradation, with METTL16 enhancing MROH8 mRNA stability through m6A modifications (P<0.01). Functional assays demonstrated that METTL16 and YTHDC2 (an m6A reader) collaboratively enhanced MROH8 mRNA stability, thereby inhibiting CAPN2 expression and reducing tumor proliferation and metastasis (P<0.001). This study reveals a novel regulatory axis involving METTL16, MROH8, and TBP that modulates CAPN2 expression, contributing to the suppression of pancreatic cancer progression. The METTL16-MROH8-TBP-CAPN2 pathway offers potential therapeutic targets for pancreatic cancer treatment, highlighting the significance of m6A modifications in tumor regulation. Further clinical validation is needed to confirm these findings in human patients.

EMBRYONIC FLOWER 1 regulates male reproduction by repressing the jasmonate pathway downstream transcription factor MYB26.

The evolutionarily conserved Polycomb repressive complexes (PRC) mediate genome-wide transcriptional silencing and regulate a plethora of development, as well as environmental responses in multicellular organisms. The PRC2-catalyzed trimethylation of lysine 27 on histone H3 (H3K27me3) is recognized by reader-effector modules of PRC1 to implement gene repression. Here, we report that the Arabidopsis (Arabidopsis thaliana) H3K27me3 effector EMBRYONIC FLOWER 1 (EMF1) interacts with and constrains the R2R3 DNA binding transcription factor MYB26 by a eudicot-conserved motif in the stamen. MYB26 activates the transcription of two NAC domain genes, NAC SECONDARY WALL THICKENING PROMOTING FACTOR1 (NST1) and NST2, whose encoded proteins mediate anther secondary cell thickening in jasmonate (JA)-regulated stamen maturation. In this process, the transcriptional activity of MYB26 is negatively modulated by the JAZ-PRC repressive complex to precisely regulate the expression of NST1 and NST2. Disruption of EMF1 repression stimulates MYB26, leading to the excessive transcription of the two NAC genes and male sterility. Our results reveal a novel mechanism in polycomb-mediated gene silencing and illustrate that the plant Polycomb complex regulates stamen development by preventing the hypersensitivity of JA responses in male reproduction.

MethyLasso: a segmentation approach to analyze DNA methylation patterns and identify differentially methylated regions from whole-genome datasets.

DNA methylation is an epigenetic mark involved in the regulation of gene expression, and patterns of DNA methylation anticorrelate with chromatin accessibility and transcription factor binding. DNA methylation can be profiled at the single cytosine resolution in the whole genome and has been performed in many cell types and conditions. Computational approaches are then essential to study DNA methylation patterns in a single condition or capture dynamic changes of DNA methylation levels across conditions. Toward this goal, we developed MethyLasso, a new approach to segment DNA methylation data. We use it as an all-in-one tool to perform the identification of low-methylated regions, unmethylated regions, DNA methylation valleys and partially methylated domains in a single condition as well as differentially methylated regions between two conditions. We performed a rigorous benchmarking comparing existing approaches by evaluating the agreement of the regions across tools, their number, size, level of DNA methylation, boundaries, cytosine-guanine content and coverage using several real datasets as well as the sensitivity and precision of the approaches using simulated data and show that MethyLasso performs best overall. MethyLasso is freely available at https://github.com/bardetlab/methylasso.