Transcription Factor Binding Sites Research Articles

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.

Somatic mosaicism in schizophrenia brains reveals prenatal mutational processes.

Germline mutations modulate the risk of developing schizophrenia (SCZ). Much less is known about the role of mosaic somatic mutations in the context of SCZ. Deep (239×) whole-genome sequencing (WGS) of brain neurons from 61 SCZ cases and 25 controls postmortem identified mutations occurring during prenatal neurogenesis. SCZ cases showed increased somatic variants in open chromatin, with increased mosaic CpG transversions (CpG>GpG) and T>G mutations at transcription factor binding sites (TFBSs) overlapping open chromatin, a result not seen in controls. Some of these variants alter gene expression, including SCZ risk genes and genes involved in neurodevelopment. Although these mutational processes can reflect a difference in factors indirectly involved in disease, increased somatic mutations at developmental TFBSs could also potentially contribute to SCZ.

A Rare Noncoding Enhancer Variant in SCN5A Contributes to the High Prevalence of Brugada Syndrome in Thailand.

Brugada syndrome (BrS) is a cardiac arrhythmia disorder that causes sudden death in young adults. Rare genetic variants in the SCN5A gene encoding the Nav1.5 sodium channel and common noncoding variants at this locus are robustly associated with the condition. BrS is particularly prevalent in Southeast Asia but the underlying ancestry-specific factors remain largely unknown. Genome sequencing of BrS probands and population-matched controls from Thailand was performed to identify rare noncoding variants at the SCN5A-SCN10A locus that were enriched in patients with BrS. A likely causal variant was prioritized by computational methods and introduced into human induced pluripotent stem cell (hiPSC) lines using CRISPR-Cas9. The effect of the variant on SCN5A expression and Nav1.5 sodium channel current was then assessed in hiPSC-derived cardiomyocytes (hiPSC-CMs). A rare noncoding variant in an SCN5A intronic enhancer region was highly enriched in patients with BrS (detected in 3.9% of cases with a case-control odds ratio of 45.2). The variant affects a nucleotide conserved across all mammalian species and predicted to disrupt a Mef2 transcription factor binding site. Heterozygous introduction of the enhancer variant in hiPSC-CMs caused significantly reduced SCN5A expression from the variant-containing allele and a 30% reduction in Nav1.5-mediated sodium current density compared with isogenic controls, confirming its pathogenicity. Patients with the variant had severe phenotypes, with 89% experiencing cardiac arrest. This is the first example of a functionally validated rare noncoding variant at the SCN5A locus and highlights how genome sequencing in understudied populations can identify novel disease mechanisms. The variant partly explains the increased prevalence of BrS in this region and enables the identification of at-risk variant carriers to reduce the burden of sudden cardiac death in Thailand.

Targeting circadian transcriptional programs through a cis-regulatory mechanism in triple negative breast cancer.

Circadian clock genes are emerging targets in many types of cancer, but their mechanistic contributions to tumor progression are still largely unknown. This makes it challenging to stratify patient populations and develop corresponding treatments. In this work, we show that in breast cancer, the disrupted expression of circadian genes has the potential to serve as biomarkers. We also show that the master circadian transcription factors (TFs) BMAL1 and CLOCK are required for the proliferation of metastatic mesenchymal stem-like (mMSL) triple-negative breast cancer (TNBC) cells. Using currently available small molecule modulators, we found that a stabilizer of cryptochrome 2 (CRY2), the direct repressor of BMAL1 and CLOCK transcriptional activity, synergizes with inhibitors of proteasome, which is required for BMAL1 and CLOCK function, to repress a transcriptional program comprising circadian cycling genes in mMSL TNBC cells. Omics analyses on drug-treated cells implied that this repression of transcription is mediated by the transcription factor binding sites (TFBSs) features in the cis-regulatory elements (CRE) of clock-controlled genes. Through a massive parallel reporter assay, we defined a set of CRE features that are potentially repressed by the specific drug combination. The identification of cis -element enrichment might serve as a new concept of defining and targeting tumor types through the modulation of cis -regulatory programs, and ultimately provide a new paradigm of therapy design for cancer types with unclear drivers like TNBC.

A compendium of genetic variations associated with promoter usage across 49 human tissues

Promoters play a crucial role in regulating gene transcription. However, our understanding of how genetic variants influence alternative promoter selection is still incomplete. In this study, we implement a framework to identify genetic variants that affect the relative usage of alternative promoters, known as promoter usage quantitative trait loci (puQTLs). By constructing an atlas of human puQTLs across 49 different tissues from 838 individuals, we have identified approximately 76,856 independent loci associated with promoter usage, encompassing 602,009 genetic variants. Our study demonstrates that puQTLs represent a distinct type of molecular quantitative trait loci, effectively uncovering regulatory targets and patterns. Furthermore, puQTLs are regulating in a tissue-specific manner and are enriched with binding sites of epigenetic marks and transcription factors, especially those involved in chromatin architecture formation. Notably, we have also found that puQTLs colocalize with complex traits or diseases and contribute to their heritability. Collectively, our findings underscore the significant role of puQTLs in elucidating the molecular mechanisms underlying tissue development and complex diseases.

12427 MKRN3 Methylation Abnormalities In Patients With Pubertal Disorders

Abstract Disclosure: N.C. Grosek: None. T. Silva: None. J.C. Magnotto: None. A. Latronico: None. A.P. Canton: None. Y. Chan: None. R.S. Carroll: None. U.B. Kaiser: None. Background and Objectives: Central Precocious Puberty (CPP) is caused by early reactivation of pulsatile GnRH secretion before age 8 years in girls and 9 years in boys. Delayed puberty (DP) is defined by the absence or incomplete development of sexual characteristics by age 13 years in girls and 14 years in boys. Environmental cues affect the HPG axis and timing of puberty via epigenetic mechanisms such as DNA methylation (DNAm). Inactivating mutations in the maternally imprinted Makorin Ring Finger 3 (MKRN3) gene are the most common genetic defect associated with CPP, representing ∼9% of patients with sporadic and familial CPP. Genomic imprinting occurs through differential methylation of CpG islands at a gene promoter or enhancer region. Differentially methylated regions (DMRs) regulate genetic imprinting. We hypothesized that transcription factor binding sites (TFBS) for ZFP57 (CpG57) andNRF1 (CpGs54+55) within the islet 3.7kb upstream of MKRN3 (MRKN3 islet) may impact expression. This study aimed to analyze if aberrant MKRN3 islet DNAm changes are associated with disorders of pubertal timing. Methods and Results: Germline DNA from 2 control patients was used to analyze the MKRN3 promoter region and islet DNAm profile. Bisulfite-treated DNA was amplified by PCR with primers targeting MKRN3 promoter and islet and submitted to Sanger sequencing. All CpG islands in the promoter region were 100% methylated, but those in the MKRN3 islet were differentially methylated with ∼50-60% methylated. DNA from our CPP cohort (n=71) underwent Next-Generation Sequencing analysis of bisulfite-treated DNA of these regions to identify the DMR potentially associated with MKRN3 imprinting. The MKRN3 islet covered CpG59-CpG49. Next, samples from our CPP cohort and patients with DP (n=44) and normal pubertal timing (n=33; 45% prepubertal, 55% pubertal) underwent pyrosequencing analyses of MKRN3islet. To compare DNAm levels between controls and our cohorts at each CpG island, T-tests were used: reported as mean±SEM. In the CPP cohort, 4/11 CpGs were significantly hypomethylated compared to controls (CpGs 54, 53, 51, 50). Specifically, the NRF1 TFBS was hypomethylated compared to controls (37.8 vs 41.6 ± 1.1, p<0.001) and DP (41.4 vs 43.8 ± 0.83, p=0.044). In the DP group, we found no significant difference in DNAm compared to controls. But they had hypermethylation of 4/11 CpGs compared to the CPP group, including the NRF1 TFBS (CpGs 55, 53, 51,50). Neither the CPP nor DP cohorts had significant DNAm differences at the ZFP57 TFBS compared to controls. Conclusion: Our study showed hypomethylation of 4 out of 11 CpGs in MKRN3 islet, including the NRF1 TFBS, in CPP compared to control and DP groups. On the other hand, NRF1 TFBS was hypermethylated in DPP compared to CPP. Further studies are needed to determine if these alternations in methylation are associated with MKRN3 expression levels and pubertal initiation. Presentation: 6/1/2024

12672 MKRN3 Methylation Abnormalities In Patients With Pubertal Disorders

Abstract Disclosure: N.C. Grosek: None. T. Silva: None. J.C. Magnotto: None. A. Latronico: None. A.P. Canton: None. Y. Chan: None. R.S. Carroll: None. U.B. Kaiser: None. Background and Objectives: Central Precocious Puberty (CPP) is caused by early reactivation of pulsatile GnRH secretion before age 8 years in girls and 9 years in boys. Delayed puberty (DP) is defined by the absence or incomplete development of sexual characteristics by age 13 years in girls and 14 years in boys. Environmental cues affect the HPG axis and timing of puberty via epigenetic mechanisms such as DNA methylation (DNAm). Inactivating mutations in the maternally imprinted Makorin Ring Finger 3 (MKRN3) gene are the most common genetic defect associated with CPP, representing ∼9% of patients with sporadic and familial CPP. Genomic imprinting occurs through differential methylation of CpG islands at a gene promoter or enhancer region. Differentially methylated regions (DMRs) regulate genetic imprinting. We hypothesized that transcription factor binding sites (TFBS) for ZFP57 (CpG57) andNRF1 (CpGs54+55) within the islet 3.7kb upstream of MKRN3 (MRKN3 islet) may impact expression. This study aimed to analyze if aberrant MKRN3 islet DNAm changes are associated with disorders of pubertal timing. Methods and Results: Germline DNA from 2 control patients was used to analyze the MKRN3 promoter region and islet DNAm profile. Bisulfite-treated DNA was amplified by PCR with primers targeting MKRN3 promoter and islet and submitted to Sanger sequencing. All CpG islands in the promoter region were 100% methylated, but those in the MKRN3 islet were differentially methylated with ∼50-60% methylated. DNA from our CPP cohort (n=71) underwent Next-Generation Sequencing analysis of bisulfite-treated DNA of these regions to identify the DMR potentially associated with MKRN3 imprinting. The MKRN3 islet covered CpG59-CpG49. Next, samples from our CPP cohort and patients with DP (n=44) and normal pubertal timing (n=33; 45% prepubertal, 55% pubertal) underwent pyrosequencing analyses of MKRN3islet. To compare DNAm levels between controls and our cohorts at each CpG island, T-tests were used: reported as mean±SEM. In the CPP cohort, 4/11 CpGs were significantly hypomethylated compared to controls (CpGs 54, 53, 51, 50). Specifically, the NRF1 TFBS was hypomethylated compared to controls (37.8 vs 41.6 ± 1.1, p<0.001) and DP (41.4 vs 43.8 ± 0.83, p=0.044). In the DP group, we found no significant difference in DNAm compared to controls. But they had hypermethylation of 4/11 CpGs compared to the CPP group, including the NRF1 TFBS (CpGs 55, 53, 51,50). Neither the CPP nor DP cohorts had significant DNAm differences at the ZFP57 TFBS compared to controls. Conclusion: Our study showed hypomethylation of 4 out of 11 CpGs in MKRN3 islet, including the NRF1 TFBS, in CPP compared to control and DP groups. On the other hand, NRF1 TFBS was hypermethylated in DPP compared to CPP. Further studies are needed to determine if these alternations in methylation are associated with MKRN3 expression levels and pubertal initiation. Presentation: 6/1/2024

8519 Investigating the Role of Protein Inhibitor of Activated STAT1 (Pias1) in Growth Plate Chondrocyte Maturation

Abstract Disclosure: J. Baronas: None. U. Ahmed: None. J.N. Hirschhorn: None. N.E. Renthal: None. Development of growth plate chondrocytes is a complex process that involves a wide array of regulatory elements. Recently our laboratory conducted a genome wide CRISPR knockout (KO) screen of growth plate chondrocyte maturation to discover novel genes and gene networks in the proliferation and maturation of chondrocytes. Among targets uncovered is Protein Inhibitor of Activated STAT 1 (Pias1), an E3 SUMO ligase that plays a multi-faceted role in genetic regulation, driving post-translational modifications, directly binding to transcription factors, and functioning as a DNA-binding protein. While Pias1 has been shown to be involved in the regulation of many cellular processes across multiple tissues, its function in growth plate development and chondrocyte maturation has been previously unstudied. In our genome-wide screen, our laboratory identified that loss of Pias1 led to early maturation of chondrocytes. This project seeks to investigate the role of Pias1 in chondrocyte maturation and differentiation. We first developed a Pias1-KO chondrocyte cell line, observing upregulation of hypertrophic markers Cd200 and Col10a1 in KOs relative to control via qPCR, in addition to proliferative regulators IHH and Pth1R. Through Bulk-RNAseq of the same cell line, we identified additional chondrocyte markers as upregulated in KO, as well as osteoblast marker Sp7. Members of the HoxD family, regulators of limb development, were downregulated. Cell culture data suggests that loss of Pias1 leads to increased cell proliferation and increased matrix production, indicating dysregulation of multiple chondrocyte processes. Immunohistochemistry of micromass culture sections also suggests morphological differences between KO and control, with KO cells significantly larger in diameter. Early SUMO-1 profiling yields limited evidence of decreased global SUMOylation in KO cells, one possible mechanistic explanation for the observed dysregulation. Ongoing studies are focused on identifying specific SUMOylation targets of Pias1 in chondrocytes, analyzing proteomic differences in KO and control, seeking putative DNA-binding sites of Pias1, and observing in vivo murine Pias1-KO growth plates. Our findings suggest that PIAS1 plays a crucial role in the growth plate by prolonging growth and delaying terminal hypertrophy. Our study aims to provide a better understanding of the regulatory pathways involved in growth plate development and related phenotypes such as height, and contribute to the development of novel therapies for growth-related disorders, such as achondroplasia and hypochondroplasia. Presentation: 6/1/2024

Memory effects of transcription regulator−DNA interactions in bacteria

Memory effect refers to the phenomenon where past events influence a system's current and future states or behaviors. In biology, memory effects often arise from intra- or intermolecular interactions, leading to temporally correlated behaviors. Single-molecule studies have shown that enzymes and DNA-binding proteins can exhibit time-correlated behaviors of their activity. While memory effects are well documented and studied in vitro, no such examples exist in cells to our knowledge. Combining single-molecule tracking (SMT) and single-cell protein quantitation, we find in living Escherichia coli cells distinct temporal correlations in the binding/unbinding events on DNA by MerR- and Fur-family metalloregulators, manifesting as memory effects with timescales of ~1 s. These memory effects persist irrespective of the type of the metalloregulators or their metallation states. Moreover, these temporal correlations of metalloregulator-DNA interactions are associated with spatial confinements of the metalloregulators near their DNA binding sites, suggesting microdomains of ~100 nm in size that possibly result from the spatial organizations of the bacterial chromosome without the involvement of membranes. These microdomains likely facilitate repeated binding events, enhancing regulator-DNA contact frequency and potentially gene regulation efficiency. These findings provide unique insights into the spatiotemporal dynamics of protein-DNA interactions in bacterial cells, introducing the concept of microdomains as a crucial player in memory effect-driven gene regulation.

PDNAPred: Interpretable prediction of protein-DNA binding sites based on pre-trained protein language models

Protein-DNA interactions play critical roles in various biological processes and are essential for drug discovery. However, traditional experimental methods are labor-intensive and unable to keep pace with the increasing volume of protein sequences, leading to a substantial number of proteins lacking DNA-binding annotations. Therefore, developing an efficient computational method to identify protein-DNA binding sites is crucial. Unfortunately, most existing computational methods rely on manually selected features or protein structure information, making these methods inapplicable to large-scale prediction tasks. In this study, we introduced PDNAPred, a sequence-based method that combines two pre-trained protein language models with a designed CNN-GRU network to identify DNA-binding sites. Additionally, to tackle the issue of imbalanced dataset samples, we employed focal loss. Our comprehensive experiments demonstrated that PDNAPred significantly improved the accuracy of DNA-binding site prediction, outperforming existing state-of-the-art sequence-based methods. Remarkably, PDNAPred also achieved results comparable to advanced structure-based methods. The designed CNN-GRU network enhances its capability to detect DNA-binding sites accurately. Furthermore, we validated the versatility of PDNAPred by training it on RNA-binding site datasets, showing its potential as a general framework for amino acid binding site prediction. Finally, we conducted model interpretability analysis to elucidate the reasons behind PDNAPred's outstanding performance.

Differential gene expression and chromatin accessibility in the medial prefrontal cortex associated with individual differences in rat behavioral models of opioid use disorder.

Opioid use disorder (OUD) is a neuropsychological disease that has a devastating impact on public health. Substantial individual differences in vulnerability exist, the neurobiological substrates of which remain unclear. To address this question, we investigated genome-wide gene transcription (RNA-seq) and chromatin accessibility (ATAC-seq) in the medial prefrontal cortex (mPFC) of male and female rats exhibiting differential vulnerability in behavioral paradigms modeling different phases of OUD: Withdrawal-Induced Anhedonia (WIA), Demand, and Reinstatement. Ingenuity Pathway Analysis (IPA) of RNA-seq revealed greater changes in canonical pathways in Resilient (vs. Saline) rats in comparison to Vulnerable (vs. Saline) rats across 3 paradigms, suggesting brain adaptations that might contribute to resilience to OUD across its trajectory. Analyses of gene networks and upstream regulators implicated processes involved in oligodendrocyte maturation and myelination in WIA, neuroinflammation in Demand, and metabolism in Reinstatement. Motif analysis of ATAC-seq showed changes in chromatin accessibility to a small set of transcription factor (TF) binding sites as a function either of opioid exposure (i.e., morphine versus saline) generally or of individual vulnerability specifically. Some of these were shared across the 3 paradigms and others were unique to each. In conclusion, we have identified changes in biological pathways, TFs, and their binding motifs that vary with paradigm and OUD vulnerability. These findings point to the involvement of distinct transcriptional and epigenetic mechanisms in response to opioid exposure, vulnerability to OUD, and different stages of the disorder.

Context transcription factors establish cooperative environments and mediate enhancer communication

Many enhancers control gene expression by assembling regulatory factor clusters, also referred to as condensates. This process is vital for facilitating enhancer communication and establishing cellular identity. However, how DNA sequence and transcription factor (TF) binding instruct the formation of high regulatory factor environments remains poorly understood. Here we developed a new approach leveraging enhancer-centric chromatin accessibility quantitative trait loci (caQTLs) to nominate regulatory factor clusters genome-wide. By analyzing TF-binding signatures within the context of caQTLs and comparing episomal versus endogenous enhancer activities, we discovered a class of regulators, ‘context-only’ TFs, that amplify the activity of cell type-specific caQTL-binding TFs, that is, ‘context-initiator’ TFs. Similar to super-enhancers, enhancers enriched for context-only TF-binding sites display high coactivator binding and sensitivity to bromodomain-inhibiting molecules. We further show that binding sites for context-only and context-initiator TFs underlie enhancer coordination, providing a mechanistic rationale for how a loose TF syntax confers regulatory specificity.

Cooperativity among clustered κB sites within promoters and enhancers dictates transcriptional specificity of NF-κB-RelA along with specific cofactors.

The functional role of weak DNA binding sites for transcription factor (TF) recruitment and gene expression remains largely unknown. Our study reveals that the weak NF-κB DNA binding sites, which are abundant in gene promoters and enhancers, appear in clusters and exhibit minimal to undetectable NF-κB binding activity in isolation in vitro, yet they play prominent roles in gene regulation within native context in cells. We found nuclear concentration of RelA/p65, the predominant NF-κB, is approximately 0.2 µM in stimulated cells, challenging the idea that these weak κB sites operate through mass action-dependent binding mechanisms. Through proteomic analysis, we identified a range of nuclear factors, including various other TFs, interacting with RelA at these κB-sites. ChIP-seq, RNA-seq and phase-separated condensation analyses suggest these additional TFs, referred to as the cofactors of NF-κB, facilitate dynamic recruitment of NF-κB to clustered κB sites of specific target genes. Overall, our findings demonstrate the collective contribution of both strong and weak κB sites in occupancy of NF-κB at the promoters and enhancers, with the recruitment facilitated by a variety of cofactors. This congregation of multiple factors forming larger dynamic complexes appearing as a transcription condensate is likely to be common to all transcriptional programs.

The RNA-binding protein PCBP1 modulates transcription by recruiting the G-quadruplex-specific helicase DHX9

PCBP1, polycytosine (poly(C)) binding protein 1, an RNA and single-stranded DNA (ssDNA) binding protein, binds poly(C) DNA tracts but it remains unclear whether its ability to bind ssDNA contributes to transcriptional regulation. Here, we report that PCBP1's DNA binding sites are enriched at transcription start sites and that by binding to promoter regions, PCBP1 regulates transcription in addition to splicing and translation. At PCBP1 target genes, we show that PCBP1 interacts with several RNA/DNA hybrid (R-loop) associated G-quadruplex resolving helicases. Furthermore, we find that PCBP1 interacts with RNA Helicase A (DHX9) to modulate transcription by regulating DHX9 accumulation and activity. PCBP1 depletion leads to defects in R-loop processing and dysregulation of transcription of PCBP1 target genes. PCBP1's high sequence specificity and interaction with helicases suggest that its mechanism in transcription involves guiding helicases to specific loci during transcription, thereby modulating their activity.

Long-range transcription factor binding sites clustered regions may mediate transcriptional regulation through phase-separation interactions in early human embryo

In mammals, during the post-fertilization pre-implantation phase, the expression of cell type-specific genes is crucial for normal embryonic development, which is regulated by cis-regulatory elements (CREs). TFs control gene expression by interacting with CREs. Research shows that transcription factor binding sites (TFBSs) reflect the general characteristics of the regulatory genome. Here, we identified TFBSs from chromatin accessibility data in five stages of early human embryonic development, and quantified transcription factor binding sites-clustered regions (TFCRs) and their complexity (TC). Assigning TC values to TFCRs has made it possible to assess the functionality of these regulatory elements in a more quantitative way. Our findings reveal a robust correlation between TFCR complexity and gene expression starting from the 8Cell stage, which is when the zygotic genome is activated in humans. Furthermore, we have defined long-range TFCRs (LR-TFCRs) and conjecture that LR-TFCRs may regulate gene expression through phase-separation mechanisms during the early stages of human embryonic development.

H2S Donor SPRC Ameliorates Cardiac Aging by Suppression of JMJD3, a Histone Demethylase.

Aims: S-propargyl-cysteine (SPRC) is an endogenous hydrogen sulfide (H2S) donor obtained by modifying the structure of S-allyl cysteine in garlic. This study aims to investigate the effect of SPRC on mitigating cardiac aging and the involvement of jumonji domain-containing protein 3 (JMJD3), a histone demethylase, which represents the primary risk factor in major aging related diseases, in this process, elucidating the preliminary mechanism through which SPRC regulation of JMJD3 occurs. Results: In vitro, SPRC mitigated the elevated levels of reactive oxygen species, senescence-associated β-galactosidase, p53, and p21, reversing the decline in mitochondrial membrane potential, which represented a reduction in cellular senescence. In vivo, SPRC improved Dox-induced cardiac pathological structure and function. Overexpression of JMJD3 accelerated cardiomyocytes and cardiac senescence, whereas its knockdown in vitro reduced the senescence phenotype. The potential binding site of the upstream transcription factor of JMJD3, sheared X box binding protein 1 (XBP1s), was determined using online software. SPRC promoted the expression of cystathionine γ-lyase (CSE), which subsequently inhibited the IRE1α/XBP1s signaling pathway and decreased JMJD3 expression. Innovations: This study is the first to establish JMJD3 as a crucial regulator of cardiac aging. SPRC can alleviate cardiac aging by upregulating CSE and inhibiting endoplasmic reticulum stress pathways, which in turn suppress JMJD3 expression. Conclusions: JMJD3 plays an essential role in cardiac aging regulation, whereas SPRC can suppress the expression of JMJD3 by upregulating CSE, thus delaying cardiac aging, which suggests that SPRC may serve as an aging protective agent, and pharmacological targeting of JMJD3 may also be a promising therapeutic approach in age-related heart diseases.

Characterization of non-coding variants associated with transcription factor binding through ATAC-seq-defined footprint QTLs in liver.

Non-coding variants discovered by genome-wide association studies (GWAS) are enriched in regulatory elements harboring transcription factor (TF) binding motifs, strongly suggesting a connection between disease association and the disruption of cis-regulatory sequences. Occupancy of a TF inside a region of open chromatin can be detected in ATAC-seq where bound TFs block the transposase Tn5, leaving a pattern of relatively depleted Tn5 insertions known as a "footprint". Here, we sought to identify variants associated with TF-binding, or "footprint quantitative trait loci" (fpQTLs) in ATAC-seq data generated from 170 human liver samples. We used computational tools to scan the ATAC-seq reads to quantify TF binding likelihood as "footprint scores" at variants derived from whole genome sequencing generated in the same samples. We tested for association between genotype and footprint score and observed 693 fpQTLs associated with footprint-inferred TF binding (FDR < 5%). Given that Tn5 insertion sites are measured with base-pair resolution, we show that fpQTLs can aid GWAS and QTL fine-mapping by precisely pinpointing TF activity within broad trait-associated loci where the underlying causal variant is unknown. Liver fpQTLs were strongly enriched across ChIP-seq peaks, liver expression QTLs (eQTLs), and liver-related GWAS loci, and their inferred effect on TF binding was concordant with their effect on underlying sequence motifs in 80% of cases. We conclude that fpQTLs can reveal causal GWAS variants, define the role of TF binding site disruption in disease and provide functional insights into non-coding variants, ultimately informing novel treatments for common diseases.

Predicting bacterial transcription factor binding sites through machine learning and structural characterization based on DNA duplex stability.

Transcriptional factors (TFs) in bacteria play a crucial role in gene regulation by binding to specific DNA sequences, thereby assisting in the activation or repression of genes. Despite their central role, deciphering shape recognition of bacterial TFs-DNA interactions remains an intricate challenge. A deeper understanding of DNA secondary structures could greatly enhance our knowledge of how TFs recognize and interact with DNA, thereby elucidating their biological function. In this study, we employed machine learning algorithms to predict transcription factor binding sites (TFBS) and classify them as directed-repeat (DR) or inverted-repeat (IR). To accomplish this, we divided the set of TFBS nucleotide sequences by size, ranging from 8 to 20 base pairs, and converted them into thermodynamic data known as DNA duplex stability (DDS). Our results demonstrate that the Random Forest algorithm accurately predicts TFBS with an average accuracy of over 82% and effectively distinguishes between IR and DR with an accuracy of 89%. Interestingly, upon converting the base pairs of several TFBS-IR into DDS values, we observed a symmetric profile typical of the palindromic structure associated with these architectures. This study presents a novel TFBS prediction model based on a DDS characteristic that may indicate how respective proteins interact with base pairs, thus providing insights into molecular mechanisms underlying bacterial TFs-DNA interaction.

MLSNet: a deep learning model for predicting transcription factor binding sites.

Accurate prediction of transcription factor binding sites (TFBSs) is essential for understanding gene regulation mechanisms and the etiology of diseases. Despite numerous advances in deep learning for predicting TFBSs, their performance can still be enhanced. In this study, we propose MLSNet, a novel deep learning architecture designed specifically to predict TFBSs. MLSNet innovatively integrates multisize convolutional fusion with long short-term memory (LSTM) networks to effectively capture DNA-sparse higher-order sequence features. Further, MLSNet incorporates super token attention and Bi-LSTM to systematically extract and integrate higher-order DNA shape features. Experimental results on 165 ChIP-seq (chromatin immunoprecipitation followed by sequencing) datasets indicate that MLSNet consistently outperforms several state-of-the-art algorithms in the prediction of TFBSs. Specifically, MLSNet reports average metrics: 0.8306 for ACC, 0.8992 for AUROC, and 0.9035 for AUPRC, surpassing the second-best methods by 1.82%, 1.68%, and 1.54%, respectively. This research delineates the effectiveness of combining multi-size convolutional layers with LSTM and DNA shape-based features in enhancing predictive accuracy. Moreover, this study comprehensively assesses the variability in model performance across different cell lines and transcription factors. The source code of MLSNet is available at https://github.com/minghaidea/MLSNet.

Identification of transcription factor binding sites of &lt;i&gt;HSP60&lt;/i&gt; and &lt;i&gt;HSP90B1&lt;/i&gt; Genes’ promoter in Bubalus bubalis

The gene expression determines cellular homeostasis while spatio-temporal nature of gene expression is regulated by several cis regulatory elements, such as promoter regions which bind to transcription factors. The ever-changing global climatic conditions have raised a necessity to scan heat tolerant livestock especially the riverine buffalo in tropics. The study focuses on characterization of the promoter region of HSP60 and HSP90B1 and looking for the binding sites of transcription factors in Bubalus bubalis and comparison with the nearest relative Bos taurus to gauge the interspecies variations. The results indicated no transcription binding factors in HSP90B1 promoter (550 nucleotides) within the region analyzed except for the Elk-1 and COMP1 in the cattle while promoter region of HSP60 gene (971 nucleotides) varied highly. Nine and seven transcription factor binding sites were specifically unique in Bubalus bubalis and in Bos taurus respectively. Occurrence of 1-Oct, BR-C Z4, Hairy, and HSF was higher in buffalo. Overall, this study provided the framework to explore cis regulatory elements and contributes to decipher the transcriptional regulatory elements of HSP60 and HSP90B1 genes in Bubalus bubalis.