DNA-binding Motif Research Articles

Molecular characterization of doublesex and Mab-3 (DMRT) gene family in Ctenopharyngodon idella (grass carp).

Doublesex and Mab-3 (DMRT) gene family is a diverse group of transcriptional factors crucially involved in sex differentiation and biological processes such as body growth and differentiation in vertebrates. In this study, we analyzed DMRT genes structural characterization and physiochemical properties, and elucidated their functional role as a ligand of different gonadal receptors including androgen (AR), estrogen β (ER-β), estrogen γ (ER-γ), and progesterone (PR). All six genes of the DMRT gene family in grass carp (Ctenopharyngodon Idella Valenciennes, 1844) exhibited an acidic nature. These DMRT genes are primarily localized in the nucleus, where they play a role in DNA binding via doublesex DNA binding motif. All the DMRT gene pairsare under strong purifying selection with two segmental duplications envisaged about 18.30 (DMRT3a/DMRTA2) and 24.90 (DMRT2b/DMRT2a) million years ago (MYA). Recombination analysis revealed six potential recombinant breakpoints posing substantial evolutionary pressure for diverse cellular functioning of DMRT isoforms. Moreover, the DMRTA1 protein had a highest binding affinity of - 270.42 and - 267.16 for androgen receptors (AR) and progesterone receptors (PR), whereas, for estrogen receptors ER-β and ER-γ, the maximum binding affinity was observed with DMRT2a and DMRT2b proteins showing a docking score of - 254.22 and - 261.71, respectively. First time we studied the binding scores and interface residues of the DMRT genes as a ligand of gonadal receptors that play a crucial role in fish growth, sex development and differentiation, and spermatogenesis and oocyte maturation. The present study provides a molecular basis for DMRT genes in grass carp that may serve as a reference for in-depth phylogenomic study in other species.

Mycobacterium tuberculosis Mce3R TetR-like Repressor Forms an Asymmetric Four-Helix Bundle and Binds a Nonpalindrome Sequence†.

Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis, is a major global health concern. TetR family repressors (TFRs) are important for Mtb's adaptation to the human host environment. Our study focuses on one notable Mtb repressor, Mce3R, composed of an unusual double TFR motif. Mce3R-regulated genes encode enzymes implicated in cholesterol metabolism, resistance against reactive oxygen species, and lipid transport activities important for Mtb survival and persistence in the host and for the cellular activity of a 6-azasteroid derivative. Here, we present the structure of Mce3R bound to its DNA operator, unveiling a unique asymmetric assembly previously unreported. We obtained a candidate DNA-binding motif through MEME motif analysis, comparing intergenic regions of mce3R orthologues and identifying nonpalindromic regions conserved between orthologues. Using an electrophoretic mobility shift assay (EMSA), we confirmed that Mce3R binds to a 123-bp sequence that includes the predicted motif. Using scrambled DNA and DNA oligonucleotides of varying lengths with sequences from the upstream region of the yrbE3A (mce3) operon, we elucidated the operator region to be composed of two Mce3R binding sites, each a 25-bp asymmetric sequence separated by 53 bp. Mce3R binds with a higher affinity to the downstream site with a Kd of 2.4 ± 0.7 nM. The cryo-EM structure of Mce3R bound to the 123-bp sequence was refined to a resolution of 2.51 Å. Each Mce3R monomer comprises 21 α-helices (α1-α21) folded into an asymmetric TFR-like structure with a core asymmetric four-helix bundle. This complex has two nonidentical HTH motifs and a single ligand-binding domain. The two nonidentical HTHs from each TFR bind within the high-affinity, nonpalindromic operator motif, with Arg53 and Lys262 inserted into the major groove. Site-directed mutagenesis of Arg53 to alanine abrogated DNA binding, validating the Mce3R/DNA structure obtained. Among 811,645 particles, 63% were Mce3R homodimer bound to two duplex oligonucleotides. Mce3R homodimerizes primarily through α15, and each monomer binds to an identical site in the DNA duplex oligonucleotide.

Functional characteristics and computational model of abundant hyperactive loci in the human genome.

Enhancers and promoters are classically considered to be bound by a small set of transcription factors (TFs) in a sequence-specific manner. This assumption has come under increasing skepticism as the datasets of ChIP-seq assays of TFs have expanded. In particular, high-occupancy target (HOT) loci attract hundreds of TFs with often no detectable correlation between ChIP-seq peaks and DNA-binding motif presence. Here, we used a set of 1003 TF ChIP-seq datasets (HepG2, K562, H1) to analyze the patterns of ChIP-seq peak co-occurrence in combination with functional genomics datasets. We identified 43,891 HOT loci forming at the promoter (53%) and enhancer (47%) regions. HOT promoters regulate housekeeping genes, whereas HOT enhancers are involved in tissue-specific process regulation. HOT loci form the foundation of human super-enhancers and evolve under strong negative selection, with some of these loci being located in ultraconserved regions. Sequence-based classification analysis of HOT loci suggested that their formation is driven by the sequence features, and the density of mapped ChIP-seq peaks across TF-bound loci correlates with sequence features and the expression level of flanking genes. Based on the affinities to bind to promoters and enhancers we detected five distinct clusters of TFs that form the core of the HOT loci. We report an abundance of HOT loci in the human genome and a commitment of 51% of all TF ChIP-seq binding events to HOT locus formation thus challenging the classical model of enhancer activity and propose a model of HOT locus formation based on the existence of large transcriptional condensates.

Functional characteristics and computational model of abundant hyperactive loci in the human genome

Enhancers and promoters are classically considered to be bound by a small set of transcription factors (TFs) in a sequence-specific manner. This assumption has come under increasing skepticism as the datasets of ChIP-seq assays of TFs have expanded. In particular, high-occupancy target (HOT) loci attract hundreds of TFs with often no detectable correlation between ChIP-seq peaks and DNA-binding motif presence. Here, we used a set of 1003 TF ChIP-seq datasets (HepG2, K562, H1) to analyze the patterns of ChIP-seq peak co-occurrence in combination with functional genomics datasets. We identified 43,891 HOT loci forming at the promoter (53%) and enhancer (47%) regions. HOT promoters regulate housekeeping genes, whereas HOT enhancers are involved in tissue-specific process regulation. HOT loci form the foundation of human super-enhancers and evolve under strong negative selection, with some of these loci being located in ultraconserved regions. Sequence-based classification analysis of HOT loci suggested that their formation is driven by the sequence features, and the density of mapped ChIP-seq peaks across TF-bound loci correlates with sequence features and the expression level of flanking genes. Based on the affinities to bind to promoters and enhancers we detected five distinct clusters of TFs that form the core of the HOT loci. We report an abundance of HOT loci in the human genome and a commitment of 51% of all TF ChIP-seq binding events to HOT locus formation thus challenging the classical model of enhancer activity and propose a model of HOT locus formation based on the existence of large transcriptional condensates.

Perspectives on Codebook: sequence specificity of uncharacterized human transcription factors.

We describe an effort ("Codebook") to determine the sequence specificity of 332 putative and largely uncharacterized human transcription factors (TFs), as well as 61 control TFs. Nearly 5,000 independent experiments across multiple in vitro and in vivo assays produced motifs for just over half of the putative TFs analyzed (177, or 53%), of which most are unique to a single TF. The data highlight the extensive contribution of transposable elements to TF evolution, both in cis and trans , and identify tens of thousands of conserved, base-level binding sites in the human genome. The use of multiple assays provides an unprecedented opportunity to benchmark and analyze TF sequence specificity, function, and evolution, as further explored in accompanying manuscripts. 1,421 human TFs are now associated with a DNA binding motif. Extrapolation from the Codebook benchmarking, however, suggests that many of the currently known binding motifs for well-studied TFs may inaccurately describe the TF's true sequence preferences.

Domain architecture of the Mycobacterium tuberculosis MabR (Rv2242), a member of the PucR transcription factor family

MabR (Rv2242), a PucR-type transcription factor, plays a crucial role in regulating mycolic acid biosynthesis in Mycobacterium tuberculosis. To understand its regulatory mechanisms, we determined the crystal structures of its N-terminal and C-terminal domains. The N-terminal domain adopts a globin-like fold, while the C-terminal domain comprises an α/β GGDEF domain and an all-α effector domain with a helix-turn-helix DNA-binding motif. This unique domain combination is specific to Actinomycetes. Biochemical and computational studies suggest that full-length MabR forms both dimeric and tetrameric assemblies in solution. Structural analysis revealed two distinct dimerization interfaces within the N- and C-terminal domains, further supporting a tetrameric organization. These findings provide valuable insights into the domain architecture, oligomeric state, and potential regulatory mechanisms of MabR.

Crystal structure of the HMGA AT-hook 1 domain bound to the minor groove of AT-rich DNA and inhibition by antikinetoplastid drugs

High mobility group (HMG) proteins are intrinsically disordered nuclear non-histone chromosomal proteins that play an essential role in many biological processes by regulating the expression of numerous genes in eukaryote cells. HMGA proteins contain three DNA binding motifs, the “AT-hooks”, that bind preferentially to AT-rich sequences in the minor groove of B-form DNA. Understanding the interactions of AT-hook domains with DNA is very relevant from a medical point of view because HMGA proteins are involved in different conditions including cancer and parasitic diseases. We present here the first crystal structure (1.40 Å resolution) of the HMGA AT-hook 1 domain, bound to the minor groove of AT-rich DNA. In contrast to AT-hook 3 which bends DNA and shows a larger minor groove widening, AT-hook 1 binds neighbouring DNA molecules and displays moderate widening of DNA upon binding. The binding affinity and thermodynamics of binding were studied in solution with surface plasmon resonance (SPR)-biosensor and isothermal titration calorimetry (ITC) experiments. AT-hook 1 forms an entropy-driven 2:1 complex with (TTAA)2-containing DNA with relatively slow kinetics of association/dissociation. We show that N-phenylbenzamide-derived antikinetoplastid compounds (1–3) bind strongly and specifically to the minor groove of AT-DNA and compete with AT-hook 1 for binding. The central core of the molecule is the basis for the observed sequence selectivity of these compounds. These findings provide clues regarding a possible mode of action of DNA minor groove binding compounds that are relevant to major neglected tropical diseases such as leishmaniasis and trypanosomiasis.

Single-Cell Meta-Analysis Uncovers the Pancreatic Endothelial Cell Transcriptomic Signature and Reveals a Key Role for NKX2-3 in PLVAP Expression.

The pancreatic vasculature displays tissue-specific physiological and functional adaptations that support rapid insulin response by β-cells. However, the digestive enzymes have made it difficult to characterize pancreatic endothelial cells (ECs), resulting in the poor understanding of pancreatic EC specialization. Available single-nuclei/single-cell RNA-sequencing data sets were mined to identify pancreatic EC-enriched signature genes and to develop an integrated atlas of human pancreatic ECs. We validated the findings using independent single-nuclei/single-cell RNA-sequencing data, bulk RNA-sequencing data of isolated ECs, spatial transcriptomics data, immunofluorescence, and RNAScope of selected markers. The TF (transcription factor) NKX2-3 was expressed in HUVECs via gene transfection, and the expression of pancreatic EC-enriched signature genes was assessed via RT-qPCR. We defined a pancreatic EC-enriched gene signature conserved across species and developmental stages that included genes involved in ECM (extracellular matrix) composition (COL15A1 and COL4A1), permeability and barrier function (PLVAP, EHD4, CAVIN3, HSPG2, ROBO4, HEG1, and CLEC14A), and key signaling pathways (S1P, TGF-β [transforming growth factor-β], RHO-RAC GTPase, PI3k-AKT, and PDGF [platelet-derived growth factor]). The integrated atlas revealed the vascular hierarchy within the pancreas. We identified and validated a specialized islet capillary subpopulation characterized by genes involved in permeability (PLVAP and EHD4), immune-modulation (FABP5, HLA-C, and B2M), ECM composition (SPARC and SPARCL1), IGF (insulin-like growth factor) signaling (IGFBP7), and membrane transport (SLCO2A1, SLC2A3, and CD320). Importantly, we identified NKX2-3 as a key TF enriched in pancreatic ECs. DNA-binding motif analysis found NKX2-3 motifs in ≈40% of the signature genes. Induction of NKX2-3 in HUVECs promoted the expression of the islet capillary EC-enriched genes PLVAP and SPARCL1. We defined a validated transcriptomic signature of pancreatic ECs and uncovered their intratissue transcriptomic heterogeneity. We showed that NKX2-3 acts upstream of PLVAP and provided a single-cell online resource that can be further explored by the community: https://vasconcelos.shinyapps.io/pancreatic_endothelial/.

On the identification of differentially-active transcription factors from ATAC-seq data.

ATAC-seq has emerged as a rich epigenome profiling technique, and is commonly used to identify Transcription Factors (TFs) underlying given phenomena. A number of methods can be used to identify differentially-active TFs through the accessibility of their DNA-binding motif, however little is known on the best approaches for doing so. Here we benchmark several such methods using a combination of curated datasets with various forms of short-term perturbations on known TFs, as well as semi-simulations. We include both methods specifically designed for this type of data as well as some that can be repurposed for it. We also investigate variations to these methods, and identify three particularly promising approaches (a chromVAR-limma workflow with critical adjustments, monaLisa and a combination of GC smooth quantile normalization and multivariate modeling). We further investigate the specific use of nucleosome-free fragments, the combination of top methods, and the impact of technical variation. Finally, we illustrate the use of the top methods on a novel dataset to characterize the impact on DNA accessibility of TRAnscription Factor TArgeting Chimeras (TRAFTAC), which can deplete TFs-in our case NFkB-at the protein level.

TRAF7 determines circadian period through ubiquitination and degradation of DBP.

D-site binding protein, DBP, is a clock-controlled transcription factor and drives daily rhythms of physiological processes through the regulation of an array of genes harboring a DNA binding motif, D-box. DBP protein levels show a circadian oscillation with an extremely robust peak/trough ratio, but it is elusive how the temporal pattern is regulated by post-translational regulation. In this study, we show that DBP protein levels are down-regulated by the ubiquitin-proteasome pathway. Analysis using 19 dominant-negative forms of E2 enzymes have revealed that UBE2G1 and UBE2T mediate the degradation of DBP. A proteomic analysis of DBP-interacting proteins and database screening have identified Tumor necrosis factor Receptor-Associated Factor 7 (TRAF7), a RING-type E3 ligase, that forms a complex with UBE2G1 and/or UBE2T. Ubiquitination analysis have revealed that TRAF7 enhances K48-linked polyubiquitination of DBP in cultured cells. Overexpression of TRAF7 down-regulates DBP protein level, while knockdown of TRAF7 up-regulates DBP in cultured cells. Knockout of TRAF7 in NIH3T3 cells have revealed that TRAF7 mediates the time-of-the-day-dependent regulation of DBP levels. Furthermore, TRAF7 has a period-shortening effect on the cellular clock. Together, TRAF7 plays an important role in circadian clock oscillation through destabilization of DBP.

Mapping the IscR regulon sheds light on the regulation of iron homeostasis in Caulobacter.

The role of the iron-sulfur [Fe-S] cluster transcriptional regulator IscR in maintaining [Fe-S] homeostasis in bacteria is still poorly characterized in many groups. Caulobacter crescentus and other Alphaproteobacteria have a single operon encoding [Fe-S] cluster biosynthesis enzymes. We showed that the expression of this operon increases in iron starvation, but not in oxidative stress, and is controlled mainly by IscR. Transcriptome analysis comparing an iscR null mutant strain with the wild-type (wt) strain identified 94 differentially expressed genes (DEGs), with 47 upregulated and 47 downregulated genes in the ΔiscR mutant. We determined the IscR binding sites in conditions of sufficient or scarce iron by Chromatin Immunoprecipitation followed by DNA sequencing (ChIP-seq), identifying two distinct putative DNA binding motifs. The estimated IscR regulon comprises 302 genes, and direct binding to several regulatory regions was shown by Electrophoresis Mobility Shift Assay (EMSA). The results showed that the IscR and Fur regulons partially overlap and that IscR represses the expression of the respiration regulator FixK, fine-tuning gene regulation in response to iron and redox balance.

OVO positively regulates essential maternal pathways by binding near the transcriptional start sites in the Drosophila female germline.

Differentiation of female germline stem cells into a mature oocyte includes the expression of RNAs and proteins that drive early embryonic development in Drosophila. We have little insight into what activates the expression of these maternal factors. One candidate is the zinc-finger protein OVO. OVO is required for female germline viability and has been shown to positively regulate its own expression, as well as a downstream target, ovarian tumor, by binding to the transcriptional start site (TSS). To find additional OVO targets in the female germline and further elucidate OVO's role in oocyte development, we performed ChIP-seq to determine genome-wide OVO occupancy, as well as RNA-seq comparing hypomorphic and wild type rescue ovo alleles. OVO preferentially binds in close proximity to target TSSs genome-wide, is associated with open chromatin, transcriptionally active histone marks, and OVO-dependent expression. Motif enrichment analysis on OVO ChIP peaks identified a 5'-TAACNGT-3' OVO DNA binding motif spatially enriched near TSSs. However, the OVO DNA binding motif does not exhibit precise motif spacing relative to the TSS characteristic of RNA polymerase II complex binding core promoter elements. Integrated genomics analysis showed that 525 genes that are bound and increase in expression downstream of OVO are known to be essential maternally expressed genes. These include genes involved in anterior/posterior/germ plasm specification (bcd, exu, swa, osk, nos, aub, pgc, gcl), egg activation (png, plu, gnu, wisp, C(3)g, mtrm), translational regulation (cup, orb, bru1, me31B), and vitelline membrane formation (fs(1)N, fs(1)M3, clos). This suggests that OVO is a master transcriptional regulator of oocyte development and is responsible for the expression of structural components of the egg as well as maternally provided RNAs that are required for early embryonic development.

OVO positively regulates essential maternal pathways by binding near the transcriptional start sites in the Drosophila female germline

Differentiation of female germline stem cells into a mature oocyte includes the expression of RNAs and proteins that drive early embryonic development in Drosophila. We have little insight into what activates the expression of these maternal factors. One candidate is the zinc-finger protein OVO. OVO is required for female germline viability and has been shown to positively regulate its own expression, as well as a downstream target, ovarian tumor, by binding to the transcriptional start site (TSS). To find additional OVO targets in the female germline and further elucidate OVO’s role in oocyte development, we performed ChIP-seq to determine genome-wide OVO occupancy, as well as RNA-seq comparing hypomorphic and wild type rescue ovo alleles. OVO preferentially binds in close proximity to target TSSs genome-wide, is associated with open chromatin, transcriptionally active histone marks, and OVO-dependent expression. Motif enrichment analysis on OVO ChIP peaks identified a 5’-TAACNGT-3’ OVO DNA binding motif spatially enriched near TSSs. However, the OVO DNA binding motif does not exhibit precise motif spacing relative to the TSS characteristic of RNA polymerase II complex binding core promoter elements. Integrated genomics analysis showed that 525 genes that are bound and increase in expression downstream of OVO are known to be essential maternally expressed genes. These include genes involved in anterior/posterior/germ plasm specification (bcd, exu, swa, osk, nos, aub, pgc, gcl), egg activation (png, plu, gnu, wisp, C(3)g, mtrm), translational regulation (cup, orb, bru1, me31B), and vitelline membrane formation (fs(1)N, fs(1)M3, clos). This suggests that OVO is a master transcriptional regulator of oocyte development and is responsible for the expression of structural components of the egg as well as maternally provided RNAs that are required for early embryonic development.

Multimodal Hox5 activity generates motor neuron diversity

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

Genome-wide analysis of the MYB gene family and functional analysis of BhMYB79 in wax gourd

Wax gourd (Benincasa hispida) is an important cucurbit crop with economic and medicinal value. The myeloblastosis (MYB) gene family is one of the largest gene families in plants and regulates various biological processes, whereas the MYB gene family has not been systematically studied in wax gourd. In this study, we performed genome-wide identification of the MYB gene family in wax gourd and analyzed their phylogenetic relationship, MYB DNA-binding domain (MYB DBD), gene structure, protein motif, synteny, duplication mode and expression pattern. As a result, a total of 215 BhMYB genes (BhMYBs) were identified, belonging to four subfamilies: 1R-, 2R-, 3R- and 4R-MYB subfamilies. Genes of 1R-MYB subfamily and 2R-MYB subfamily were subdivided into different subgroups respectively. The analysis of MYB DBD, gene structure and protein motif showed that the most genes in the same subgroup have similar characteristics and the 2R-MYB genes were more conserved than the 1R-MYB genes. Interestingly, the long terminal retrotransposons (LTR-RTs) were found in the long introns of several BhMYBs. The results of synteny analysis showed that there were more syntenic gene pairs between wax gourd and other cucurbit crops, while the least number of syntenic gene pairs existed between wax gourd and rice. Gene duplication was the main reason for the expansion of the MYB gene family in wax gourd, with the transposed duplication (TRD) mode contributing more. All duplication BhMYB genes were under purifying selection pressure. Further expression analysis showed that many BhMYBs exhibited obvious tissue-specific expression and several BhMYBs were significantly induced by one or more abiotic stresses. BhMYB79 was particularly expressed in roots and significantly induced by salt, drought, cold and heat stresses, overexpression of which led to reduced tolerance to salt stress in Arabidopsis. In conclusion, our results provide a systematic analysis of wax gourd MYB gene family and facilitate the biological role study of BhMYB79 during wax gourd salt stress response process.

On the identification of differentially-active transcription factors from ATAC-seq data.

ATAC-seq has emerged as a rich epigenome profiling technique, and is commonly used to identify Transcription Factors (TFs) underlying given phenomena. A number of methods can be used to identify differentially-active TFs through the accessibility of their DNA-binding motif, however little is known on the best approaches for doing so. Here we benchmark several such methods using a combination of curated datasets with various forms of short-term perturbations on known TFs, as well as semi-simulations. We include both methods specifically designed for this type of data as well as some that can be repurposed for it. We also investigate variations to these methods, and identify three particularly promising approaches (a chromVAR-limma workflow with critical adjustments, monaLisa and a combination of GC smooth quantile normalization and multivariate modeling). We further investigate the specific use of nucleosome-free fragments, the combination of top methods, and the impact of technical variation. Finally, we illustrate the use of the top methods on a novel dataset to characterize the impact on DNA accessibility of TRAnscription Factor TArgeting Chimeras (TRAFTAC), which can deplete TFs - in our case NFkB - at the protein level.

ATAD5 functions as a regulatory platform for Ub–PCNA deubiquitination

Ubiquitination status of proliferating cell nuclear antigen (PCNA) is crucial for regulating DNA lesion bypass. After the resolution of fork stalling, PCNA is subsequently deubiquitinated, but the underlying mechanism remains undefined. We found that the N-terminal domain of ATAD5 (ATAD5-N), the largest subunit of the PCNA-unloading complex, functions as a scaffold for Ub-PCNA deubiquitination. ATAD5 recognizes DNA-loaded Ub-PCNA through distinct DNA-binding and PCNA-binding motifs. Furthermore, ATAD5 forms a heterotrimeric complex with UAF1-USP1 deubiquitinase, facilitating the deubiquitination of DNA-loaded Ub-PCNA. ATAD5 also enhances the Ub-PCNA deubiquitination by USP7 and USP11 through specific interactions. ATAD5 promotes the distinct deubiquitination process of UAF1-USP1, USP7, and USP11 for poly-Ub-PCNA. Additionally, ATAD5 mutants deficient in UAF1-binding had increased sensitivity to DNA-damaging agents. Our results ultimately reveal that ATAD5 and USPs cooperate to efficiently deubiquitinate Ub-PCNA prior to its release from the DNA in order to safely deactivate the DNA repair process.

Functional characteristics and computational model of abundant hyperactive loci in the human genome.

Enhancers and promoters are classically considered to be bound by a small set of TFs in a sequence-specific manner. This assumption has come under increasing skepticism as the datasets of ChIP-seq assays of TFs have expanded. In particular, high-occupancy target (HOT) loci attract hundreds of TFs with often no detectable correlation between ChIP-seq peaks and DNA-binding motif presence. Here, we used a set of 1,003 TF ChIP-seq datasets (HepG2, K562, H1) to analyze the patterns of ChIP-seq peak co-occurrence in combination with functional genomics datasets. We identified 43,891 HOT loci forming at the promoter (53%) and enhancer (47%) regions. HOT promoters regulate housekeeping genes, whereas HOT enhancers are involved in tissue-specific process regulation. HOT loci form the foundation of human super-enhancers and evolve under strong negative selection, with some of these loci being located in ultraconserved regions. Sequence-based classification analysis of HOT loci suggested that their formation is driven by the sequence features, and the density of mapped ChIP-seq peaks across TF-bound loci correlates with sequence features and the expression level of flanking genes. Based on the affinities to bind to promoters and enhancers we detected 5 distinct clusters of TFs that form the core of the HOT loci. We report an abundance of HOT loci in the human genome and a commitment of 51% of all TF ChIP-seq binding events to HOT locus formation thus challenging the classical model of enhancer activity and propose a model of HOT locus formation based on the existence of large transcriptional condensates.

Role of Csdc2 in Regulating Secondary Hair Follicle Growth in Cashmere Goats.

Cashmere goats possess two types of hair follicles, with the secondary hair follicles producing valuable cashmere fiber used for textiles. The growth of cashmere exhibits a seasonal pattern arising from photoperiod change. Transcription factors play crucial roles during this process. The transcription factor, cold-shock domain, containing C2 (Csdc2) plays a crucial role in modulating cell proliferation and differentiation. Our preceding research indicated that the expression of Csdc2 changes periodically during anagen to telogen. However, the mechanisms of Csdc2 in regulating SHF growth remain unclear. Here, we found that the knockdown of Csdc2 inhibits the proliferation of dermal papilla cells. ChIP-Seq analysis showed that Csdc2 had a unique DNA binding motif in SHFs. Through conjoint analysis of ChIP-Seq and RNA-Seq, we revealed a total of 25 candidate target genes of Csdc2. Notably, we discovered a putative Csdc2 binding site within roundabout guidance receptor 2 (Robo2) on chromosome 1 of the goat genome. Furthermore, qRT-PCR and dual-luciferase reporter assay confirmed Csdc2's positive regulatory influence on Robo2. These findings expand the research field of hair follicle transcriptional regulatory networks, offering insights into molecular breeding strategies to enhance cashmere production in goats.

Elucidating the molecular mechanisms in phytohormones-induced alleviation of postharvest physiological deterioration in cassava tuberous roots

Cassava is a staple food crop in tropical regions, serving as a significant source of carbohydrates. However, rapid postharvest physiological deterioration (PPD) significantly limits the shelf life and utilization of cassava. Despite its importance, the overlapping molecular mechanisms underlying phytohormones-induced PPD alleviation in cassava tuberous roots remain unknown. In this study, the phytohormones ethephon (ET), gibberellic acid (GA), and methyl jasmonate (MeJA) were applied to cassava tuberous roots after harvest. These treatments led to significantly reduced H2O2 accumulation and alleviation of PPD. Transcriptomic analysis revealed shared functions, including transcriptional regulation and oxidoreductase activity, among differentially expressed genes (DEGs) affected by ET, GA, or MeJA. A total of 305 genes encoding antioxidant enzymes were systematically identified, with 151 responsive to phytohormone treatments across ten families. Notably, the expression of 10 Glutathione S-transferase (GST) genes was induced by all three phytohormone treatments, while their expression was suppressed in the absence of phytohormones during PPD. By using co-expression network analysis, DNA-binding motif analysis, and experimental validation using yeast one-hybrid and dual-luciferase assays, a regulatory network involving MeMYB102 and MeZAT11 transcription factors (TFs) that directly promote the expression of MeGST42, MeGST43, MeGST44, and MeGST46 was established. Overall, this study enhances the understanding of the overlapping molecular mechanisms underlying phytohormones-induced alleviation of PPD in cassava, thereby providing genetic resources for developing PPD-resistant cassava.