Inferring direct DNA binding from ChIP-seq

A continuum of zinc finger transcription factor retention on native chromatin underlies dynamic genome organization

Transcription factor (TF) residence on chromatin translates into quantitative transcriptional or structural outcomes on genome. Commonly used formaldehyde crosslinking fixes TF-DNA interactions cumulatively and compromises the measured occupancy level. Here we mapped the occupancy level of global or individual zinc finger TFs like CTCF and MAZ, in the form of highly resolved footprints, on native chromatin. By incorporating reinforcing perturbation conditions, we established S-score, a quantitative metric to proxy the continuum of CTCF or MAZ retention across different motifs on native chromatin. The native chromatin-retained CTCF sites harbor sequence features within CTCF motifs better explained by S-score than the metrics obtained from other crosslinking or native assays. CTCF retention on native chromatin correlates with local SUMOylation level, and anti-correlates with transcriptional activity. The S-score successfully delineates the otherwise-masked differential stability of chromatin structures mediated by CTCF, or by MAZ independent of CTCF. Overall, our study established a paradigm continuum of TF retention across binding sites on native chromatin, explaining the dynamic genome organization.

FLI1 is associated with regulation of DNA methylation and megakaryocytic differentiation in FPDMM caused by a RUNX1 transactivation domain mutation

Familial platelet disorder with associated myeloid malignancies (FPDMM) is an autosomal dominant disease caused by heterozygous germline mutations in RUNX1. It is characterized by thrombocytopenia, platelet dysfunction, and a predisposition to hematological malignancies. Although FPDMM is a precursor for diseases involving abnormal DNA methylation, the DNA methylation status in FPDMM remains unknown, largely due to a lack of animal models and challenges in obtaining patient-derived samples. Here, using genome editing techniques, we established two lines of human induced pluripotent stem cells (iPSCs) with different FPDMM-mimicking heterozygous RUNX1 mutations. These iPSCs showed defective differentiation of hematopoietic progenitor cells (HPCs) and megakaryocytes (Mks), consistent with FPDMM. The FPDMM-mimicking HPCs showed DNA methylation patterns distinct from those of wild-type HPCs, with hypermethylated regions showing the enrichment of ETS transcription factor (TF) motifs. We found that the expression of FLI1, an ETS family member, was significantly downregulated in FPDMM-mimicking HPCs with a RUNX1 transactivation domain (TAD) mutation. We demonstrated that FLI1 promoted binding-site-directed DNA demethylation, and that overexpression of FLI1 restored their megakaryocytic differentiation efficiency and hypermethylation status. These findings suggest that FLI1 plays a crucial role in regulating DNA methylation and correcting defective megakaryocytic differentiation in FPDMM-mimicking HPCs with a RUNX1 TAD mutation.

The BAS chromatin remodeler determines brassinosteroid-induced transcriptional activation and plant growth in Arabidopsis

Brassinosteroid (BR) signaling leads to the nuclear accumulation of the BRASSINAZOLE-RESISTANT 1 (BZR1) transcription factor, which plays dual roles in activating or repressing the expression of thousands of genes. BZR1 represses gene expression by recruiting histone deacetylases, but how it activates transcription of BR-induced genes remains unclear. Here, we show that BR reshapes the genome-wide chromatin accessibility landscape, increasing the accessibility of BR-induced genes and reducing the accessibility of BR-repressed genes in Arabidopsis. BZR1 physically interacts with the BRAHMA-associated SWI/SNF (BAS)-chromatin-remodeling complex on the genome and selectively recruits the BAS complex to BR-activated genes. Depletion of BAS abrogates the capacities of BZR1 to increase chromatin accessibility, activate gene expression, and promote cell elongation without affecting BZR1's ability to reduce chromatin accessibility and expression of BR-repressed genes. Together, these data identify that BZR1 recruits the BAS complex to open chromatin and to mediate BR-induced transcriptional activation of growth-promoting genes.

BrTTG1 regulates seed coat proanthocyanidin formation through a direct interaction with structural gene promoters of flavonoid pathway and glutathione S-transferases in Brassica rapa L

Introduction

Seed coat color is a significant agronomic trait in horticultural crops such as Brassica rapa which is characterized by brown or yellow seed coat coloration. Previous Brassica rapa studies have shown that BrTTG1 is responsible for seed coat proanthocyanidin formation, which is dependent on the MYB-bHLH-WD40 complex, whereas some studies have reported that TRANSPARENT TESTA GLABRA 1 (TTG1) directly interacts with the structural gene promoters of the flavonoid pathway.

Methods

Herein, the brown-seeded inbred B147 and ttg1 yellow-seeded inbred B80 mutants were used as plant materials for gene expression level analysis, gene promoter clone and transient overexpression.

Results

The analysis identified eleven structural genes involved in the flavonoid biosynthesis pathway, which are potentially responsible for BrTTG1- dependent seed coat proanthocyanidin formation. The promoters of these genes were cloned and cis-acting elements were identified. Yeast one-hybrid and dual-luciferase assays confirmed that BrTTG1 directly and independently interacted with proCHS-Bra008792, proDFR-Bra027457, proTT12-Bra003361, proTT19-Bra008570, proTT19-Bra023602 and proAHA10-Bra016610. A TTG1-binding motif (RTWWGTRGM) was also identified. Overexpression of TTG1 in the yellow-seed B. rapa inbred induced proanthocyanidin accumulation by increasing the expression levels of related genes.

Discussion

Our study unveiled, for the first time, the direct interaction between TTG1 and the promoters of the flavonoid biosynthesis pathway structural genes and glutathione S-transferases in Brassica rapa. Additionally, we have identified a novel TTG1-binding motif, providing a basis for further exploration into the function of TTG1 and the accumulation of proanthocyanidins in seed coats.

Cancer-associated fibroblast-derived acetate promotes pancreatic cancer development by altering polyamine metabolism via the ACSS2-SP1-SAT1 axis

The ability of tumour cells to thrive in harsh microenvironments depends on the utilization of nutrients available in the milieu. Here we show that pancreatic cancer-associated fibroblasts (CAFs) regulate tumour cell metabolism through the secretion of acetate, which can be blocked by silencing ATP citrate lyase (ACLY) in CAFs. We further show that acetyl-CoA synthetase short-chain family member 2 (ACSS2) channels the exogenous acetate to regulate the dynamic cancer epigenome and transcriptome, thereby facilitating cancer cell survival in an acidic microenvironment. Comparative H3K27ac ChIP-seq and RNA-seq analyses revealed alterations in polyamine homeostasis through regulation of SAT1 gene expression and enrichment of the SP1-responsive signature. We identified acetate/ACSS2-mediated acetylation of SP1 at the lysine 19 residue that increased SP1 protein stability and transcriptional activity. Genetic or pharmacologic inhibition of the ACSS2-SP1-SAT1 axis diminished the tumour burden in mouse models. These results reveal that the metabolic flexibility imparted by the stroma-derived acetate enabled cancer cell survival under acidosis via the ACSS2-SP1-SAT1 axis.

Development of sexual dimorphism of skeletal muscles through the adrenal cortex, caused by androgen-induced global gene suppression

The zona fasciculata (zF) in the adrenal cortex contributes to multiple physiological actions through glucocorticoid synthesis. The size, proliferation, and glucocorticoid synthesis characteristics are all female biased, and sexual dimorphism is established by androgen. In this study, transcriptomes were obtained to unveil the sex differentiation mechanism. Interestingly, both the amount of mRNA and the expressions of nearly all genes were higher in females. The expression of Nr5a1, which is essential for steroidogenic cell differentiation, was also female biased. Whole-genome studies demonstrated that NR5A1 regulates nearly all gene expression directly or indirectly. This suggests that androgen-induced global gene suppression is potentially mediated by NR5A1. Using Nr5a1 heterozygous mice, whose adrenal cortex is smaller than the wild type, we demonstrated that the size of skeletal muscles is possibly regulated by glucocorticoid synthesized by zF. Taken together, considering the ubiquitous presence of glucocorticoid receptors, our findings provide a pathway for sex differentiation through glucocorticoid synthesis.

Interplay between coding and non-coding regulation drives the Arabidopsis seed-to-seedling transition

Translation of seed stored mRNAs is essential to trigger germination. However, when RNAPII re-engages RNA synthesis during the seed-to-seedling transition has remained in question. Combining csRNA-seq, ATAC-seq and smFISH in Arabidopsis thaliana we demonstrate that active transcription initiation is detectable during the entire germination process. Features of non-coding regulation such as dynamic changes in chromatin accessible regions, antisense transcription, as well as bidirectional non-coding promoters are widespread throughout the Arabidopsis genome. We show that sensitivity to exogenous ABSCISIC ACID (ABA) during germination depends on proximal promoter accessibility at ABA-responsive genes. Moreover, we provide genetic validation of the existence of divergent transcription in plants. Our results reveal that active enhancer elements are transcribed producing non-coding enhancer RNAs (eRNAs) as widely documented in metazoans. In sum, this study defining the extent and role of coding and non-coding transcription during key stages of germination expands our understanding of transcriptional mechanisms underlying plant developmental transitions.

Functional Characterisation of the ATOH1 Molecular Subtype Indicates a Pro-Metastatic Role in Small Cell Lung Cancer

Molecular subtypes of Small Cell Lung Cancer (SCLC) have been described based on differential expression of transcription factors (TFs) ASCL1, NEUROD1, POU2F3 and immune-related genes. We previously reported an additional subtype based on expression of the neurogenic TF ATOH1 within our SCLC Circulating tumour cell-Derived eXplant (CDX) model biobank. Here we show that ATOH1 protein was detected in 7/81 preclinical models and 16/102 clinical samples of SCLC. In CDX models, ATOH1 directly regulated neurogenesis and differentiation programs consistent with roles in normal tissues. In ex vivo cultures of ATOH1-positive CDX, ATOH1 was required for cell survival. In vivo, ATOH1 depletion slowed tumour growth and suppressed liver metastasis. Our data validate ATOH1 as a bona fide oncogenic driver of SCLC with tumour cell survival and pro-metastatic functions. Further investigation to explore ATOH1 driven vulnerabilities for targeted treatment with predictive biomarkers is warranted.

The transcription factor NF-YA is crucial for neural progenitor maintenance during brain development

In contrast to stage-specific transcription factors, the role of ubiquitous transcription factors in neuronal development remains a matter of scrutiny. Here, we demonstrated that a ubiquitous factor NF-Y is essential for neural progenitor maintenance during brain morphogenesis. Deletion of the NF-YA subunit in neural progenitors by using nestin-cre transgene in mice resulted in significant abnormalities in brain morphology, including a thinner cerebral cortex and loss of striatum during embryogenesis. Detailed analyses revealed a progressive decline in multiple neural progenitors in the cerebral cortex and ganglionic eminences, accompanied by induced apoptotic cell death and reduced cell proliferation. In neural progenitors, the NF-YA short isoform lacking exon 3 is dominant and co-expressed with cell cycle genes. ChIP-seq analysis from the cortex during early corticogenesis revealed preferential binding of NF-Y to the cell cycle genes, some of which were confirmed to be downregulated following NF-YA deletion. Notably, the NF-YA short isoform disappears and is replaced by its long isoform during neuronal differentiation. Forced expression of the NF-YA long isoform in neural progenitors resulted in a significant decline in neuronal count, possibly due to the suppression of cell proliferation. Collectively, we elucidated a critical role of the NF-YA short isoform in maintaining neural progenitors, possibly by regulating cell proliferation and apoptosis. Moreover, we identified an isoform switch in NF-YA within the neuronal lineage in vivo, which may explain the stage-specific role of NF-Y during neuronal development.

A gene regulatory network critical for axillary bud dormancy directly controlled by Arabidopsis BRANCHED1

The Arabidopsis thaliana transcription factor BRANCHED1 (BRC1) plays a pivotal role in the control of shoot branching as it integrates environmental and endogenous signals that influence axillary bud growth. Despite its remarkable activity as a growth inhibitor, the mechanisms by which BRC1 promotes bud dormancy are largely unknown. We determined the genome-wide BRC1 binding sites in vivo and combined these with transcriptomic data and gene co-expression analyses to identify bona fide BRC1 direct targets. Next, we integrated multi-omics data to infer the BRC1 gene regulatory network (GRN) and used graph theory techniques to find network motifs that control the GRN dynamics. We generated an open online tool to interrogate this network. A group of BRC1 target genes encoding transcription factors (BTFs) orchestrate this intricate transcriptional network enriched in abscisic acid-related components. Promoter::β-GLUCURONIDASE transgenic lines confirmed that BTFs are expressed in axillary buds. Transient co-expression assays and studies in planta using mutant lines validated the role of BTFs in modulating the GRN and promoting bud dormancy. This knowledge provides access to the developmental mechanisms that regulate shoot branching and helps identify candidate genes to use as tools to adapt plant architecture and crop production to ever-changing environmental conditions.

Modeling methyl-sensitive transcription factor motifs with an expanded epigenetic alphabet

BACKGROUND

Transcription factors bind DNA in specific sequence contexts. In addition to distinguishing one nucleobase from another, some transcription factors can distinguish between unmodified and modified bases. Current models of transcription factor binding tend not to take DNA modifications into account, while the recent few that do often have limitations. This makes a comprehensive and accurate profiling of transcription factor affinities difficult.

RESULTS

Here, we develop methods to identify transcription factor binding sites in modified DNA. Our models expand the standard A/C/G/T DNA alphabet to include cytosine modifications. We develop Cytomod to create modified genomic sequences and we also enhance the MEME Suite, adding the capacity to handle custom alphabets. We adapt the well-established position weight matrix (PWM) model of transcription factor binding affinity to this expanded DNA alphabet. Using these methods, we identify modification-sensitive transcription factor binding motifs. We confirm established binding preferences, such as the preference of ZFP57 and C/EBPβ for methylated motifs and the preference of c-Myc for unmethylated E-box motifs.

CONCLUSIONS

Using known binding preferences to tune model parameters, we discover novel modified motifs for a wide array of transcription factors. Finally, we validate our binding preference predictions for OCT4 using cleavage under targets and release using nuclease (CUT&RUN) experiments across conventional, methylation-, and hydroxymethylation-enriched sequences. Our approach readily extends to other DNA modifications. As more genome-wide single-base resolution modification data becomes available, we expect that our method will yield insights into altered transcription factor binding affinities across many different modifications.

HOCOMOCO in 2024: a rebuild of the curated collection of binding models for human and mouse transcription factors

We present a major update of the HOCOMOCO collection that provides DNA binding specificity patterns of 949 human transcription factors and 720 mouse orthologs. To make this release, we performed motif discovery in peak sets that originated from 14 183 ChIP-Seq experiments and reads from 2554 HT-SELEX experiments yielding more than 400 thousand candidate motifs. The candidate motifs were annotated according to their similarity to known motifs and the hierarchy of DNA-binding domains of the respective transcription factors. Next, the motifs underwent human expert curation to stratify distinct motif subtypes and remove non-informative patterns and common artifacts. Finally, the curated subset of 100 thousand motifs was supplied to the automated benchmarking to select the best-performing motifs for each transcription factor. The resulting HOCOMOCO v12 core collection contains 1443 verified position weight matrices, including distinct subtypes of DNA binding motifs for particular transcription factors. In addition to the core collection, HOCOMOCO v12 provides motif sets optimized for the recognition of binding sites in vivo and in vitro, and for annotation of regulatory sequence variants. HOCOMOCO is available at https://hocomoco12.autosome.org and https://hocomoco.autosome.org.

Decoding the gene regulatory network of endosperm differentiation in maize

The persistent cereal endosperm constitutes the majority of the grain volume. Dissecting the gene regulatory network underlying cereal endosperm development will facilitate yield and quality improvement of cereal crops. Here, we use single-cell transcriptomics to analyze the developing maize (Zea mays) endosperm during cell differentiation. After obtaining transcriptomic data from 17,022 single cells, we identify 12 cell clusters corresponding to five endosperm cell types and revealing complex transcriptional heterogeneity. We delineate the temporal gene-expression pattern from 6 to 7 days after pollination. We profile the genomic DNA-binding sites of 161 transcription factors differentially expressed between cell clusters and constructed a gene regulatory network by combining the single-cell transcriptomic data with the direct DNA-binding profiles, identifying 181 regulons containing genes encoding transcription factors along with their high-confidence targets, Furthermore, we map the regulons to endosperm cell clusters, identify cell-cluster-specific essential regulators, and experimentally validated three predicted key regulators. This study provides a framework for understanding cereal endosperm development and function at single-cell resolution.

Loss of the DYRK1A Protein Kinase Results in the Reduction in Ribosomal Protein Gene Expression, Ribosome Mass and Reduced Translation

Ribosomal proteins (RPs) are evolutionary conserved proteins that are essential for protein translation. RP expression must be tightly regulated to ensure the appropriate assembly of ribosomes and to respond to the growth demands of cells. The elements regulating the transcription of RP genes (RPGs) have been characterized in yeast and Drosophila, yet how cells regulate the production of RPs in mammals is less well understood. Here, we show that a subset of RPG promoters is characterized by the presence of the palindromic TCTCGCGAGA motif and marked by the recruitment of the protein kinase DYRK1A. The presence of DYRK1A at these promoters is associated with the enhanced binding of the TATA-binding protein, TBP, and it is negatively correlated with the binding of the GABP transcription factor, establishing at least two clusters of RPGs that could be coordinately regulated. However, DYRK1A silencing leads to a global reduction in RPGs mRNAs, pointing at DYRK1A activities beyond those dependent on its chromatin association. Significantly, cells in which DYRK1A is depleted have reduced RP levels, fewer ribosomes, reduced global protein synthesis and a smaller size. We therefore propose a novel role for DYRK1A in coordinating the expression of genes encoding RPs, thereby controlling cell growth in mammals.

RNA molecular recording with an engineered RNA deaminase

RNA deaminases are powerful tools for base editing and RNA molecular recording. However, the enzymes used in currently available RNA molecular recorders such as TRIBE, DART or STAMP have limitations due to RNA structure and sequence dependence. We designed a platform for directed evolution of RNA molecular recorders. We engineered an RNA A-to-I deaminase (an RNA adenosine base editor, rABE) that has high activity, low bias and low background. Using rABE, we present REMORA (RNA-encoded molecular recording in adenosines), wherein deamination by rABE writes a molecular record of RNA-protein interactions. By combining rABE with the C-to-U deaminase APOBEC1 and long-read RNA sequencing, we measured binding by two RNA-binding proteins on single messenger RNAs. Orthogonal RNA molecular recording of mammalian Pumilio proteins PUM1 and PUM2 shows that PUM1 competes with PUM2 for a subset of sites in cells. Furthermore, we identify transcript isoform-specific RNA-protein interactions driven by isoform changes distal to the binding site. The genetically encodable RNA deaminase rABE enables single-molecule identification of RNA-protein interactions with cell type specificity.

DNA elements tether canonical Polycomb Repressive Complex 1 to human genes

Development of multicellular animals requires epigenetic repression by Polycomb group proteins. The latter assemble in multi-subunit complexes, of which two kinds, Polycomb Repressive Complex 1 (PRC1) and Polycomb Repressive Complex 2 (PRC2), act together to repress key developmental genes. How PRC1 and PRC2 recognize specific genes remains an open question. Here we report the identification of several hundreds of DNA elements that tether canonical PRC1 to human developmental genes. We use the term tether to describe a process leading to a prominent presence of canonical PRC1 at certain genomic sites, although the complex is unlikely to interact with DNA directly. Detailed analysis indicates that sequence features associated with PRC1 tethering differ from those that favour PRC2 binding. Throughout the genome, the two kinds of sequence features mix in different proportions to yield a gamut of DNA elements that range from those tethering predominantly PRC1 or PRC2 to ones capable of tethering both complexes. The emerging picture is similar to the paradigmatic targeting of Polycomb complexes by Polycomb Response Elements (PREs) of Drosophila but providing for greater plasticity.

Dosage amplification dictates oncogenic regulation by the NKX2-1 lineage factor in lung adenocarcinoma

Amplified oncogene expression is a critical and widespread driver event in cancer, yet our understanding of how amplification-mediated elevated dosage mediates oncogenic regulation is limited. Here, we find that the most significant focal amplification event in lung adenocarcinoma (LUAD) targets a lineage super-enhancer near the NKX2-1 lineage transcription factor. The NKX2-1 super-enhancer is targeted by focal and co-amplification with NKX2-1, and activation or repression controls NKX2-1 expression. We find that NKX2-1 is a widespread dependency in LUAD cell lines, where NKX2-1 pioneers enhancer accessibility to drive a lineage addicted state in LUAD, and NKX2-1 confers persistence to EGFR inhibitors. Notably, we find that oncogenic NKX2-1 regulation requires expression above a minimum dosage threshold-NKX2-1 dosage below this threshold is insufficient for cell viability, enhancer remodeling, and TKI persistence. Our data suggest that copy-number amplification can be a gain-of-function alteration, wherein amplification elevates oncogene expression above a critical dosage required for oncogenic regulation and cancer cell survival.

A fungal sesquiterpene biosynthesis gene cluster critical for mutualist-pathogen transition in Colletotrichum tofieldiae

Plant-associated fungi show diverse lifestyles from pathogenic to mutualistic to the host; however, the principles and mechanisms through which they shift the lifestyles require elucidation. The root fungus Colletotrichum tofieldiae (Ct) promotes Arabidopsis thaliana growth under phosphate limiting conditions. Here we describe a Ct strain, designated Ct3, that severely inhibits plant growth. Ct3 pathogenesis occurs through activation of host abscisic acid pathways via a fungal secondary metabolism gene cluster related to the biosynthesis of sesquiterpene metabolites, including botrydial. Cluster activation during root infection suppresses host nutrient uptake-related genes and changes mineral contents, suggesting a role in manipulating host nutrition state. Conversely, disruption or environmental suppression of the cluster renders Ct3 beneficial for plant growth, in a manner dependent on host phosphate starvation response regulators. Our findings indicate that a fungal metabolism cluster provides a means by which infectious fungi modulate lifestyles along the parasitic-mutualistic continuum in fluctuating environments.

Regulation of Grain Chalkiness and Starch Metabolism by FLO2 Interaction Factor 3, a bHLH Transcription Factor in Oryza sativa

Chalkiness is a key determinant that directly affects the appearance and cooking quality of rice grains. Previously, Floury endosperm 2 (FLO2) was reported to be involved in the formation of rice chalkiness; however, its regulation mechanism is still unclear. Here, FLO2 interaction factor 3 (OsFIF3), a bHLH transcription factor, was identified and analyzed in Oryza sativa. A significant increase in chalkiness was observed in OsFIF3-overexpressed grains, coupled with a round, hollow filling of starch granules and reduced grain weight. OsFIF3 is evolutionarily conserved in monocotyledons, but variable in dicotyledons. Subcellular localization revealed the predominant localization of OsFIF3 in the nucleus. The DAP-seq (DNA affinity purification sequencing) results showed that OsFIF3 could affect the transcriptional accumulation of β-amylase 1, α-amylase isozyme 2A-like, pectinesterase 11, β-glucosidase 28 like, pectinesterase, sucrose transport protein 1 (SUT1), and FLO2 through the binding of the CACGTG motif on their promoters. Moreover, FLO2 and SUT1 with abundant OsFIF3 binding signals showed significant expression reduction in OsFIF3 overexpression lines, further confirming OsFIF3's role in starch metabolism regulation and energy material allocation. Taken together, these findings show that the overexpression of OsFIF3 inhibits the expression of FLO2 and SUT1, thereby increasing grain chalkiness and affecting grain weight.

An adult clock component links circadian rhythms to pancreatic β-cell maturation

How ubiquitous circadian clocks orchestrate tissue-specific outputs is not well understood. Pancreatic β cell-autonomous clocks attune insulin secretion to daily energy cycles, and desynchrony from genetic or behavioral disruptions raises type 2 diabetes risk. We show that the transcription factor DEC1, a clock component induced in adult β cells, coordinates their glucose responsiveness by synchronizing energy metabolism and secretory gene oscillations. Dec1-ablated mice develop lifelong hypo-insulinemic diabetes, despite normal islet formation and intact circadian Clock and Bmal1 activators. DEC1, but not CLOCK/BMAL1, binds maturity-linked genes that mediate respiratory metabolism and insulin exocytosis, and Dec1 loss disrupts their transcription synchrony. Accordingly, β-cell Dec1 ablation causes hypo-insulinemia due to immature glucose responsiveness, dampening insulin rhythms. Thus, Dec1 links circadian clockwork to the β-cell maturation process, aligning metabolism to diurnal energy cycles.

Glucocorticoid stimulation induces regionalized gene responses within topologically associating domains

Transcription-factor binding to cis-regulatory regions regulates the gene expression program of a cell, but occupancy is often a poor predictor of the gene response. Here, we show that glucocorticoid stimulation led to the reorganization of transcriptional coregulators MED1 and BRD4 within topologically associating domains (TADs), resulting in active or repressive gene environments. Indeed, we observed a bias toward the activation or repression of a TAD when their activities were defined by the number of regions gaining and losing MED1 and BRD4 following dexamethasone (Dex) stimulation. Variations in Dex-responsive genes at the RNA levels were consistent with the redistribution of MED1 and BRD4 at the associated cis-regulatory regions. Interestingly, Dex-responsive genes without the differential recruitment of MED1 and BRD4 or binding by the glucocorticoid receptor were found within TADs, which gained or lost MED1 and BRD4, suggesting a role of the surrounding environment in gene regulation. However, the amplitude of the response of Dex-regulated genes was higher when the differential recruitment of the glucocorticoid receptor and transcriptional coregulators was observed, reaffirming the role of transcription factor-driven gene regulation and attributing a lesser role to the TAD environment. These results support a model where a signal-induced transcription factor induces a regionalized effect throughout the TAD, redefining the notion of direct and indirect effects of transcription factors on target genes.

Comparative transcriptome characterization of esophageal squamous cell carcinoma and adenocarcinoma

Background

Esophageal cancers are primarily categorized as esophageal squamous cell carcinoma (ESCC) and esophageal adenocarcinoma (EAC). While various (epi) genomic alterations associated with tumor development in ESCC and EAC have been documented, a comprehensive comparison of the transcriptomes in these two cancer subtypes remains lacking.

Methods

We collected 551 gene expression profiles from publicly available sources, including normal, ESCC, and EAC tissues or cell lines. Subsequently, we conducted a systematic analysis to compare the transcriptomes of these samples at various levels, including gene expression, promoter activity, alternative splicing (AS), alternative polyadenylation (APA), and gene fusion.

Results

Seven distinct cluster gene expression patterns were identified among the differentially expressed genes in normal, ESCC, and EAC tissues. These patterns were enriched in the PI3K-Akt signaling pathway and the activation of extracellular matrix organization and exhibited repression of epidermal development. Notably, we observed additional genes or unique expression levels enriched in these shared pathways and biological processes related to tumor development and immune activation. In addition to the differentially expressed genes, there was an enrichment of lncRNA co-expression networks and downregulation of promoter activity associated with the repression of epidermal development in both ESCC and EAC. This indicates a common feature between these two cancer subtypes. Furthermore, differential AS and APA patterns in ESCC and EAC appear to partially affect the expression of host genes associated with bacterial or viral infections in these subtypes. No gene fusions were observed between ESCC and EAC, thus highlighting the distinct molecular mechanisms underlying these two cancer subtypes.

Conclusions

We conducted a comprehensive comparison of ESCC and EAC transcriptomes and uncovered shared and distinct transcriptomic signatures at multiple levels. These findings suggest that ESCC and EAC may exhibit common and unique mechanisms involved in tumorigenesis.

Forkhead transcription factor FKH-8 cooperates with RFX in the direct regulation of sensory cilia in Caenorhabditis elegans

Cilia, either motile or non-motile (a.k.a primary or sensory), are complex evolutionarily conserved eukaryotic structures composed of hundreds of proteins required for their assembly, structure and function that are collectively known as the ciliome. Ciliome gene mutations underlie a group of pleiotropic genetic diseases known as ciliopathies. Proper cilium function requires the tight coregulation of ciliome gene transcription, which is only fragmentarily understood. RFX transcription factors (TF) have an evolutionarily conserved role in the direct activation of ciliome genes both in motile and non-motile cilia cell-types. In vertebrates, FoxJ1 and FoxN4 Forkhead (FKH) TFs work with RFX in the direct activation of ciliome genes, exclusively in motile cilia cell-types. No additional TFs have been described to act together with RFX in primary cilia cell-types in any organism. Here we describe FKH-8, a FKH TF, as a direct regulator of the sensory ciliome genes in Caenorhabditis elegans. FKH-8 is expressed in all ciliated neurons in C. elegans, binds the regulatory regions of ciliome genes, regulates ciliome gene expression, cilium morphology and a wide range of behaviors mediated by sensory ciliated neurons. FKH-8 and DAF-19 (C. elegans RFX) physically interact and synergistically regulate ciliome gene expression. C. elegans FKH-8 function can be replaced by mouse FOXJ1 and FOXN4 but not by other members of other mouse FKH subfamilies. In conclusion, RFX and FKH TF families act jointly as direct regulators of ciliome genes also in sensory ciliated cell types suggesting that this regulatory logic could be an ancient trait predating functional cilia sub-specialization.

Genome-Wide Analysis of Hypoxia-Inducible Factor Binding Reveals Targets Implicated in Impaired Human Placental Syncytiotrophoblast Formation under Low Oxygen

Preeclampsia (PE) is a common and serious complication of pregnancy with no cure except premature delivery. The root cause of PE is improper development of the placenta-the temporary organ supporting fetal growth and development. Continuous formation of the multinucleated syncytiotrophoblast (STB) layer via differentiation and fusion of cytotrophoblasts (CTBs) is vital for healthy placentation and is impaired in preeclamptic pregnancies. In PE, there is reduced/intermittent placental perfusion, likely resulting in a persistently low O2 environment. Low O2 inhibits differentiation and fusion of CTBs into STB and may thus contribute to PE pathogenesis; however, the underlying mechanisms are unknown. Because low O2 activates a transcription factor complex in cells known as the hypoxia-inducible factor (HIF), the objective of this study was to investigate whether HIF signaling inhibits STB formation by regulating genes required for this process. Culture of primary CTBs, the CTB-like cell line BeWo, and human trophoblast stem cells under low O2 reduced cell fusion and differentiation into STB. Knockdown of aryl hydrocarbon receptor nuclear translocator (a key component of the HIF complex) in BeWo cells restored syncytialization and expression of STB-associated genes under different O2 levels. Chromatin immunoprecipitation sequencing facilitated the identification of global aryl hydrocarbon receptor nuclear translocator/HIF binding sites, including several near genes implicated in STB development, such as ERVH48-1 and BHLHE40, providing new insights into mechanisms underlying pregnancy diseases linked to poor placental O2 supply.

Identification of ecdysone receptor target genes in the worker honey bee brains during foraging behavior

Ecdysone signaling plays central roles in morphogenesis and female ovarian development in holometabolous insects. In the European honey bee (Apis mellifera L.), however, ecdysone receptor (EcR) is expressed in the brains of adult workers, which have already undergone metamorphosis and are sterile with shrunken ovaries, during foraging behavior. Aiming at unveiling the significance of EcR signaling in the worker brain, we performed chromatin-immunoprecipitation sequencing of EcR to search for its target genes using the brains of nurse bees and foragers. The majority of the EcR targets were common between the nurse bee and forager brains and some of them were known ecdysone signaling-related genes. RNA-sequencing analysis revealed that some EcR target genes were upregulated in forager brains during foraging behavior and some were implicated in the repression of metabolic processes. Single-cell RNA-sequencing analysis revealed that EcR and its target genes were expressed mostly in neurons and partly in glial cells in the optic lobes of the forager brain. These findings suggest that in addition to its role during development, EcR transcriptionally represses metabolic processes during foraging behavior in the adult worker honey bee brain.

Genome-Wide Mapping of the Escherichia coli PhoB Regulon Reveals Many Transcriptionally Inert, Intragenic Binding Sites

Genome-scale analyses have revealed many transcription factor binding sites within, rather than upstream of, genes, raising questions as to the function of these binding sites. Here, we use complementary approaches to map the regulon of the Escherichia coli transcription factor PhoB, a response regulator that controls transcription of genes involved in phosphate homeostasis. Strikingly, the majority of PhoB binding sites are located within genes, but these intragenic sites are not associated with detectable transcription regulation and are not evolutionarily conserved. Many intragenic PhoB sites are located in regions bound by H-NS, likely due to shared sequence preferences of PhoB and H-NS. However, these PhoB binding sites are not associated with transcription regulation even in the absence of H-NS. We propose that for many transcription factors, including PhoB, binding sites not associated with promoter sequences are transcriptionally inert and hence are tolerated as genomic "noise." IMPORTANCE Recent studies have revealed large numbers of transcription factor binding sites within the genes of bacteria. The function, if any, of the vast majority of these binding sites has not been investigated. Here, we map the binding of the transcription factor PhoB across the Escherichia coli genome, revealing that the majority of PhoB binding sites are within genes. We show that PhoB binding sites within genes are not associated with regulation of the overlapping genes. Indeed, our data suggest that bacteria tolerate the presence of large numbers of nonregulatory, intragenic binding sites for transcription factors and that these binding sites are not under selective pressure.

An epigenetic switch controls an alternative NR2F2 isoform that unleashes a metastatic program in melanoma

Metastatic melanoma develops once transformed melanocytic cells begin to de-differentiate into migratory and invasive melanoma cells with neural crest cell (NCC)-like and epithelial-to-mesenchymal transition (EMT)-like features. However, it is still unclear how transformed melanocytes assume a metastatic melanoma cell state. Here, we define DNA methylation changes that accompany metastatic progression in melanoma patients and discover Nuclear Receptor Subfamily 2 Group F, Member 2 - isoform 2 (NR2F2-Iso2) as an epigenetically regulated metastasis driver. NR2F2-Iso2 is transcribed from an alternative transcriptional start site (TSS) and it is truncated at the N-terminal end which encodes the NR2F2 DNA-binding domain. We find that NR2F2-Iso2 expression is turned off by DNA methylation when NCCs differentiate into melanocytes. Conversely, this process is reversed during metastatic melanoma progression, when NR2F2-Iso2 becomes increasingly hypomethylated and re-expressed. Our functional and molecular studies suggest that NR2F2-Iso2 drives metastatic melanoma progression by modulating the activity of full-length NR2F2 (Isoform 1) over EMT- and NCC-associated target genes. Our findings indicate that DNA methylation changes play a crucial role during metastatic melanoma progression, and their control of NR2F2 activity allows transformed melanocytes to acquire NCC-like and EMT-like features. This epigenetically regulated transcriptional plasticity facilitates cell state transitions and metastatic spread.

Genome-wide mapping of the Escherichia coli PhoB regulon reveals many transcriptionally inert, intragenic binding sites

Genome-scale analyses have revealed many transcription factor binding sites within, rather than upstream of genes, raising questions as to the function of these binding sites. Here, we use complementary approaches to map the regulon of the Escherichia coli transcription factor PhoB, a response regulator that controls transcription of genes involved in phosphate homeostasis. Strikingly, the majority of PhoB binding sites are located within genes, but these intragenic sites are not associated with detectable transcription regulation and are not evolutionarily conserved. Many intragenic PhoB sites are located in regions bound by H-NS, likely due to shared sequence preferences of PhoB and H-NS. However, these PhoB binding sites are not associated with transcription regulation even in the absence of H-NS. We propose that for many transcription factors, including PhoB, binding sites not associated with promoter sequences are transcriptionally inert, and hence are tolerated as genomic "noise".

IMPORTANCE

Recent studies have revealed large numbers of transcription factor binding sites within the genes of bacteria. The function, if any, of the vast majority of these binding sites has not been investigated. Here, we map the binding of the transcription factor PhoB across the Escherichia coli genome, revealing that the majority of PhoB binding sites are within genes. We show that PhoB binding sites within genes are not associated with regulation of the overlapping genes. Indeed, our data suggest that bacteria tolerate the presence of large numbers of non-regulatory, intragenic binding sites for transcription factors, and that these binding sites are not under selective pressure.

BTG1 mutation yields supercompetitive B cells primed for malignant transformation

Multicellular life requires altruistic cooperation between cells. The adaptive immune system is a notable exception, wherein germinal center B cells compete vigorously for limiting positive selection signals. Studying primary human lymphomas and developing new mouse models, we found that mutations affecting BTG1 disrupt a critical immune gatekeeper mechanism that strictly limits B cell fitness during antibody affinity maturation. This mechanism converted germinal center B cells into supercompetitors that rapidly outstrip their normal counterparts. This effect was conferred by a small shift in MYC protein induction kinetics but resulted in aggressive invasive lymphomas, which in humans are linked to dire clinical outcomes. Our findings reveal a delicate evolutionary trade-off between natural selection of B cells to provide immunity and potentially dangerous features that recall the more competitive nature of unicellular organisms.

A Systemic and Integrated Analysis of p63-Driven Regulatory Networks in Mouse Oral Squamous Cell Carcinoma

Oral squamous cell carcinoma (OSCC) is the most common malignancy of the oral cavity and is linked to tobacco exposure, alcohol consumption, and human papillomavirus infection. Despite therapeutic advances, a lack of molecular understanding of disease etiology, and delayed diagnoses continue to negatively affect survival. The identification of oncogenic drivers and prognostic biomarkers by leveraging bulk and single-cell RNA-sequencing datasets of OSCC can lead to more targeted therapies and improved patient outcomes. However, the generation, analysis, and continued utilization of additional genetic and genomic tools are warranted. Tobacco-induced OSCC can be modeled in mice via 4-nitroquinoline 1-oxide (4NQO), which generates a spectrum of neoplastic lesions mimicking human OSCC and upregulates the oncogenic master transcription factor p63. Here, we molecularly characterized established mouse 4NQO treatment-derived OSCC cell lines and utilized RNA and chromatin immunoprecipitation-sequencing to uncover the global p63 gene regulatory and signaling network. We integrated our p63 datasets with published bulk and single-cell RNA-sequencing of mouse 4NQO-treated tongue and esophageal tumors, respectively, to generate a p63-driven gene signature that sheds new light on the role of p63 in murine OSCC. Our analyses reveal known and novel players, such as COTL1, that are regulated by p63 and influence various oncogenic processes, including metastasis. The identification of new sets of potential biomarkers and pathways, some of which are functionally conserved in human OSCC and can prognosticate patient survival, offers new avenues for future mechanistic studies.

Pathogen-driven gene expression patterns lead to a novel approach to the identification of common therapeutic targets

Developing a common medication strategy for disease control and management could be greatly beneficial. Investigating the differences between diseased and healthy states using differentially expressed genes aids in understanding disease pathophysiology and enables the exploration of protein-drug interactions. This study aimed to find the most common genes in diarrhea-causing bacteria such as Salmonella enterica serovar Typhimurium, Campylobacter jejuni, Escherichia coli, Shigella dysenteriae (CESS) to find new drugs. Thus, differential gene expression datasets of CESS were screened through computational algorithms and programming. Subsequently, hub and common genes were prioritized from the analysis of extensive protein-protein interactions. Binding predictions were performed to identify the common potential therapeutic targets of CESS. We identified a total of 827 dysregulated genes that are highly linked to CESS. Notably, no common gene interaction was found among all CESS bacteria, but we identified 3 common genes in both Salmonella-Escherichia and Escherichia-Campylobacter infections. Later, out of 73 protein complexes, molecular simulations confirmed 5 therapeutic candidates from the CESS. We have developed a new pipeline for identifying therapeutic targets for a common medication strategy against CESS. However, further wet-lab validation is needed to confirm their effectiveness.

Characterization and analysis of the promoter region of monodehydroascorbate reductase 4 (CpMDAR4) in papaya

Differential spatial and temporal expression patterns due to regulatory cis-elements and two different isoforms are detected among CpMDAR4 alleles in papaya. The aim of this research was to study the effects of cis-element differences between the X, Y and Y^h alleles on the expression of CpMDAR4, a potential candidate gene for sex differentiation in papaya, using a transcriptional reporter system in a model species Arabidopsis thaliana. Possible effects of a retrotransposon insertion in the Y and Y^h alleles on the transcription and expression of CpMDAR4 alleles in papaya flowers were also examined. When comparing promoters and cis-regulatory elements among genes in the non-recombining region of the sex chromosomes, paired genes exhibited differences. Our results showed that differences in the promoter sequences of the CpMDAR4 alleles drove the expression of a reporter gene to different flower tissues in Arabidopsis. β-glucuronidase staining analysis of T₂ and T₃ lines for constructs containing 5' deletions of native Y and Y^h allele promoters showed the loss of specific expression of the reporter gene in the anthers, confirming the existence and location of cis-regulatory element POLLEN1LELAT52. The expression analysis of CpMDAR4 alleles in papaya flowers also showed that all alleles are actively expressed in different flower tissues, with the existence of a shorter truncated isoform, with unknown function, for the Y and Y^h alleles due to an LTR-RT insertion in the Y and Y^h chromosomes. The observed expression patterns in Arabidopsis thaliana flowers and the expression patterns of CpMDAR4 alleles in papaya flowers suggest that MDAR4 might have a role on development of reproductive organs in papaya, and that it constitutes an important candidate for sex differentiation.

Mitochondrial dysfunction compromises ciliary homeostasis in astrocytes

Astrocytes, often considered as secondary responders to neurodegeneration, are emerging as primary drivers of brain disease. Here we show that mitochondrial DNA depletion in astrocytes affects their primary cilium, the signaling organelle of a cell. The progressive oxidative phosphorylation deficiency in astrocytes induces FOXJ1 and RFX transcription factors, known as master regulators of motile ciliogenesis. Consequently, a robust gene expression program involving motile cilia components and multiciliated cell differentiation factors are induced. While the affected astrocytes still retain a single cilium, these organelles elongate and become remarkably distorted. The data suggest that chronic activation of the mitochondrial integrated stress response (ISRmt) in astrocytes drives anabolic metabolism and promotes ciliary elongation. Collectively, our evidence indicates that an active signaling axis involving mitochondria and primary cilia exists and that ciliary signaling is part of ISRmt in astrocytes. We propose that metabolic ciliopathy is a novel pathomechanism for mitochondria-related neurodegenerative diseases.

GATA4 Regulates Developing Endocardium Through Interaction With ETS1

BACKGROUND

The pioneer transcription factor (TF) GATA4 (GATA Binding Protein 4) is expressed in multiple cardiovascular lineages and is essential for heart development. GATA4 lineage-specific occupancy in the developing heart underlies its lineage specific activities. Here, we characterized GATA4 chromatin occupancy in cardiomyocyte and endocardial lineages, dissected mechanisms that control lineage specific occupancy, and analyzed GATA4 regulation of endocardial gene expression.

METHODS

We mapped GATA4 chromatin occupancy in cardiomyocyte and endocardial cells of embryonic day 12.5 (E12.5) mouse heart using lineage specific, Cre-activated biotinylation of GATA4. Regulation of GATA4 pioneering activity was studied in cell lines stably overexpressing GATA4. GATA4 regulation of endocardial gene expression was analyzed using single cell RNA sequencing and luciferase reporter assays.

RESULTS

Cardiomyocyte-selective and endothelial-selective GATA4 occupied genomic regions had features of lineage specific enhancers. Footprints within cardiomyocyte- and endothelial-selective GATA4 regions were enriched for NKX2-5 (NK2 homeobox 5) and ETS1 (ETS Proto-Oncogene 1) motifs, respectively, and both of these TFs interacted with GATA4 in co-immunoprecipitation assays. In stable NIH3T3 cell lines expressing GATA4 with or without NKX2-5 or ETS1, the partner TFs re-directed GATA4 pioneer binding and augmented its ability to open previously inaccessible regions, with ETS1 displaying greater potency as a pioneer partner than NKX2-5. Single-cell RNA sequencing of embryonic hearts with endothelial cell-specific Gata4 inactivation identified Gata4-regulated endocardial genes, which were adjacent to GATA4-bound, endothelial regions enriched for both GATA4 and ETS1 motifs. In reporter assays, GATA4 and ETS1 cooperatively stimulated endothelial cell enhancer activity.

CONCLUSIONS

Lineage selective non-pioneer TFs NKX2-5 and ETS1 guide the activity of pioneer TF GATA4 to bind and open chromatin and create active enhancers and mechanistically link ETS1 interaction to GATA4 regulation of endocardial development.

Single-cell profiling of vascular endothelial cells reveals progressive organ-specific vulnerabilities during obesity

Obesity promotes diverse pathologies, including atherosclerosis and dementia, which frequently involve vascular defects and endothelial cell (EC) dysfunction. Each organ has distinct EC subtypes, but whether ECs are differentially affected by obesity is unknown. Here we use single-cell RNA sequencing to analyze transcriptomes of ~375,000 ECs from seven organs in male mice at progressive stages of obesity to identify organ-specific vulnerabilities. We find that obesity deregulates gene expression networks, including lipid handling, metabolic pathways and AP1 transcription factor and inflammatory signaling, in an organ- and EC-subtype-specific manner. The transcriptomic aberrations worsen with sustained obesity and are only partially mitigated by dietary intervention and weight loss. For example, dietary intervention substantially attenuates dysregulation of liver, but not kidney, EC transcriptomes. Through integration with human genome-wide association study data, we further identify a subset of vascular disease risk genes that are induced by obesity. Our work catalogs the impact of obesity on the endothelium, constitutes a useful resource and reveals leads for investigation as potential therapeutic targets.

A map of cis-regulatory modules and constituent transcription factor binding sites in 80% of the mouse genome

BACKGROUND

Mouse is probably the most important model organism to study mammal biology and human diseases. A better understanding of the mouse genome will help understand the human genome, biology and diseases. However, despite the recent progress, the characterization of the regulatory sequences in the mouse genome is still far from complete, limiting its use to understand the regulatory sequences in the human genome.

RESULTS

Here, by integrating binding peaks in ~ 9,000 transcription factor (TF) ChIP-seq datasets that cover 79.9% of the mouse mappable genome using an efficient pipeline, we were able to partition these binding peak-covered genome regions into a cis-regulatory module (CRM) candidate (CRMC) set and a non-CRMC set. The CRMCs contain 912,197 putative CRMs and 38,554,729 TF binding sites (TFBSs) islands, covering 55.5% and 24.4% of the mappable genome, respectively. The CRMCs tend to be under strong evolutionary constraints, indicating that they are likely cis-regulatory; while the non-CRMCs are largely selectively neutral, indicating that they are unlikely cis-regulatory. Based on evolutionary profiles of the genome positions, we further estimated that 63.8% and 27.4% of the mouse genome might code for CRMs and TFBSs, respectively.

CONCLUSIONS

Validation using experimental data suggests that at least most of the CRMCs are authentic. Thus, this unprecedentedly comprehensive map of CRMs and TFBSs can be a good resource to guide experimental studies of regulatory genomes in mice and humans.

The chromatin factor ROW cooperates with BEAF-32 in regulating long-range inducible genes

Insulator proteins located at the boundaries of topological associated domains (TAD) are involved in higher-order chromatin organization and transcription regulation. However, it is still not clear how long-range contacts contribute to transcriptional regulation. Here, we show that relative-of-WOC (ROW) is essential for the long-range transcription regulation mediated by the boundary element-associated factor of 32kD (BEAF-32). We find that ROW physically interacts with heterochromatin proteins (HP1b and HP1c) and the insulator protein (BEAF-32). These proteins interact at TAD boundaries where ROW, through its AT-hook motifs, binds AT-rich sequences flanked by BEAF-32-binding sites and motifs. Knockdown of row downregulates genes that are long-range targets of BEAF-32 and bound indirectly by ROW (without binding motif). Analyses of high-throughput chromosome conformation capture (Hi-C) data reveal long-range interactions between promoters of housekeeping genes bound directly by ROW and promoters of developmental genes bound indirectly by ROW. Thus, our results show cooperation between BEAF-32 and the ROW complex, including HP1 proteins, to regulate the transcription of developmental and inducible genes through long-range interactions.

The transcriptional repressors VAL1 and VAL2 mediate genome-wide recruitment of the CHD3 chromatin remodeler PICKLE in Arabidopsis

PICKLE (PKL) is a chromodomain helicase DNA-binding domain 3 (CHD3) chromatin remodeler that plays essential roles in controlling the gene expression patterns that determine developmental identity in plants, but the molecular mechanisms through which PKL is recruited to its target genes remain elusive. Here, we define a cis-motif and trans-acting factors mechanism that governs the genomic occupancy profile of PKL in Arabidopsis thaliana. We show that two homologous trans-factors VIVIPAROUS1/ABI3-LIKE1 (VAL1) and VAL2 physically interact with PKL in vivo, localize extensively to PKL-occupied regions in the genome, and promote efficient PKL recruitment at thousands of target genes, including those involved in seed maturation. Transcriptome analysis and genetic interaction studies reveal a close cooperation of VAL1/VAL2 and PKL in regulating gene expression and developmental fate. We demonstrate that this recruitment operates at two master regulatory genes, ABSCISIC ACID INSENSITIVE3 and AGAMOUS-LIKE 15, to repress the seed maturation program and ensure the seed-to-seedling transition. Together, our work unveils a general rule through which the CHD3 chromatin remodeler PKL binds to its target chromatin in plants.

PIF4 enhances DNA binding of CDF2 to co-regulate target gene expression and promote Arabidopsis hypocotyl cell elongation

How specificity is conferred within gene regulatory networks is an important problem in biology. The basic helix-loop-helix PHYTOCHROME-INTERACTING FACTORs (PIFs) and single zinc-finger CYCLING DOF FACTORs (CDFs) mediate growth responses of Arabidopsis to light and temperature. We show that these two classes of transcription factor (TF) act cooperatively. CDF2 and PIF4 are temporally and spatially co-expressed, they interact to form a protein complex and act in the same genetic pathway to promote hypocotyl cell elongation. Furthermore, PIF4 substantially strengthens genome-wide occupancy of CDF2 at a subset of its target genes. One of these, YUCCA8, encodes an auxin biosynthesis enzyme whose transcription is increased by PIF4 and CDF2 to contribute to hypocotyl elongation. The binding sites of PIF4 and CDF2 in YUCCA8 are closely spaced, and in vitro PIF4 enhances binding of CDF2. We propose that this occurs by direct protein interaction and because PIF4 binding alters DNA conformation. Thus, we define mechanisms by which PIF and CDF TFs cooperate to achieve regulatory specificity and promote cell elongation in response to light.

Zika Virus Infection Alters Gene Expression and Poly-Adenylation Patterns in Placental Cells

Flaviviruses are small RNA viruses that are mainly transmitted via arthropod vectors and are found in tropic and sub-tropical regions. Most infections are asymptomatic (90-95%), but symptoms can be as severe as hemorrhagic fever and encephalitis. One recently emerged flavivirus is Zika virus (ZIKV), which was originally isolated from rhesus monkeys in Uganda roughly 70 years ago but has recently spread east, reaching S. America in 2015-2016. This outbreak was associated with the development of Guillain-Barré syndrome in adults and microcephaly in infants born to expectant mothers infected early in pregnancy. ZIKV must traverse the placenta to impact the development of the fetus, but the mechanisms responsible are unknown. While flaviviruses are known to disrupt splicing patterns in host cells, little is known about how flaviviruses such as ZIKV impact the alternative polyadenylation (APA) of host transcripts. This is important as APA is well-established as a mechanism in the regulation of mRNA metabolism and translation. Thus, we sought to characterize transcriptomic changes including APA in human placental (JEG3) cells in response to ZIKV infection using Poly(A)-ClickSeq (PAC-Seq). We used our differential Poly(A)-cluster (DPAC) analysis pipeline to characterize changes in differential gene expression, alternative poly-adenylation (APA) and the use of alternative terminal exons. We identified 98 upregulated genes and 28 downregulated genes. Pathway enrichment analysis indicated that many RNA processing and immune pathways were upregulated in ZIKV-infected JEG3 cells. We also updated DPAC to provide additional metrics of APA including the percentage-distal usage index (PDUI), which revealed that APA was extensive and the 3' UTRs of 229 genes were lengthened while 269 were shortened. We further found that there were 214 upregulated and 59 downregulated poly(A)-clusters (PACs). We extracted the nucleotide sequences surrounding these PACs and found that the canonical signals for poly-adenylation (binding site for poly-A binding protein (PABP) upstream and a GU-rich region down-stream of the PAC) were only enriched in the downregulated PACs. These results indicate that ZIKV infection makes JEG3 cells more permissive to non-canonical poly-adenylation signals.

Differential dependencies of human RNA polymerase II promoters on TBP, TAF1, TFIIB and XPB

The effects of rapid acute depletion of components of RNA polymerase II (Pol II) general transcription factors (GTFs) that are thought to be critical for formation of preinitiation complexes (PICs) and initiation in vitro were quantified in HAP1 cells using precision nuclear run-on sequencing (PRO-Seq). The average dependencies for each factor across >70 000 promoters varied widely even though levels of depletions were similar. Some of the effects could be attributed to the presence or absence of core promoter elements such as the upstream TBP-specificity motif or downstream G-rich sequences, but some dependencies anti-correlated with such sequences. While depletion of TBP had a large effect on most Pol III promoters only a small fraction of Pol II promoters were similarly affected. TFIIB depletion had the largest general effect on Pol II and also correlated with apparent termination defects downstream of genes. Our results demonstrate that promoter activity is combinatorially influenced by recruitment of TFIID and sequence-specific transcription factors. They also suggest that interaction of the preinitiation complex (PIC) with nucleosomes can affect activity and that recruitment of TFIID containing TBP only plays a positive role at a subset of promoters.

Global N6-methyladenosine methylation analysis reveals the positive correlation between m6A modification and mRNA abundance during Apostichopus japonicus disease development

N⁶-methyladenosine (m⁶A), the most abundant epitranscriptomic modification in eukaryotic messenger RNA (mRNA), plays important roles in regulation of gene expression for fundamental biological processes and diverse physiological functions, including combating with pathogen infection. Here, we were first profile transcriptome-wide m⁶A sequencing in four stages of skin ulceration syndrome-diseased Apostichopus japonicus following Vibrio splendidus infection, including Control (healthy), Early (small ulcer), Later (extensive ulcer), and Resistant (no ulcer) groups. Our results revealed that three experimental groups were all extensively methylated by m⁶A and the proportion of the m⁶A modified genes were also significantly increased to 28.90% (Early), 27.97% (Later), and 29.98% (Resistant) when compared with Control group (15.15%), indicating m⁶A modification could be induced by V. splendidus infection. Intriguingly, we discovered a positive correlation between the m⁶A methylation level and mRNA abundance, indicating a positive regulatory role of m⁶A in sea cucumber gene expression during V. splendidus infection. Moreover, genes with specific and differentially expressed m⁶A methylation in Later group were both enriched in cell adhesion, while Early and Resistant groups were both mainly involved in DNA conformation change and chromosome organization when compared with Control, suggesting the higher-methylated m⁶A might serve as "conformational marker" and associated to the initiation of related anti-disease genes transcription in order to improve disease resistance of sea cucumber. Subsequently, we selected the pivotal genes enriched in cell adhesion pathway and found that the IggFc-binding protein (FcGBP) and Fibrocystin-L both had higher levels of m⁶A methylation and higher level of mRNA expressions in Later group. Conversely, Fibrinogen C domain-containing protein 1 (F1BCD1) gene presented as an antibacterial role in sea cucumber and showed higher mRNA expression and higher m⁶A methylation in Resistant group and lower mRNA level in Later group. The levels of m⁶A methylation and mRNA abundance of FcGBP and F1BCD1 genes indicates disease occurrence or disease resistant were also verified by MeRIP-qPCR. Overall, our study presents the first comprehensive characterize of dynamic m⁶A methylation modification in the different stages of disease in sea cucumber. These data provide an invaluable resource for future studies of function and biological significance of m⁶A in mRNA in marine invertebrates.

In silico Analysis of ACE2 Receptor to Find Potential Herbal Drugs in COVID-19 Associated Neurological Dysfunctions

COVID-19 mainly causes the collapse of the pulmonary system thereby causing a dearth of oxygen in the human body. Patients infected with this viral disease have been reported to experience various signs and symptoms associated with brain dysfunction, from the feeling of vagueness to loss of smell and taste to severe strokes. These neurological problems have been reported by younger COVID-19 infected patients mainly in their thirties and forties. Various researchers from around the globe have discerned numerous other brain dysfunctions, such as headache, dizziness, numbness, major depressive disorder, anosmia, encephalitis, febrile seizures, and Guillain-Barre syndrome. The involvement of the CNS by this viral infection has been predicted to be for a longer period of time, even if the patient recovers from COVID-19. The neuronal cell damage caused by COVID-19 is a potent factor responsible for cognitive, behavioral, and psychological problems among its sufferers. The hypoxic conditions can also trigger the formation of beta-amyloid plaques and tau-tangles and thus the virus can even induce Alzheimer’s in patients in the near future. The virus affects the brain directly, thereby causing encephalitis. This pandemic has also been shown to have a negative psychological toll on people. This research aims to highlight the brain dysfunction associated with the ACE2 receptor that is known to be a crucial player in the COVID-19 pandemic using genetic networking approaches. Furthermore, we have identified herbal drug candidates that bind to the ACE2 receptor in order to identify potential treatments for the neurological manifestations of COVID-19.

White tea alleviates non-alcoholic fatty liver disease by regulating energy expenditure and lipid metabolism

Non-alcoholic fatty liver disease (NAFLD) is one of the leading causes of liver disease, which lacks effective treatments. Abnormal lipid metabolism and inflammation are the most prominent pathological manifestations of NAFLD. Recently, it has been reported that white tea extract (WTE) can regulate lipid metabolism in human adipocytes and liver cancer cells in vitro. However, its beneficial effects on NAFLD and the underlying mechanisms remain largely unknown. Here, we showed that WTE alleviated obesity, lipid accumulation, hepatic steatosis, and liver injury in a mouse model of NAFLD. Mechanistically, we demonstrated that WTE exerted the anti-NAFLD effect by decreasing the expression of genes involved in lipid transport and synthesis processes while activating genes associated with energy expenditure. In addition, a comparison of the transcriptional responses of WTE with that of green tea extract (GTE) revealed that WTE can not only regulate lipid metabolism and stress response like GTE but also regulate antioxidant and inflammatory pathways more effectively. Taken together, our findings demonstrate that WTE inhibits the progression of NAFLD in a mouse model and indicate that WTE can be a potential dietary intervention for NAFLD.

MYT1L in the making: emerging insights on functions of a neurodevelopmental disorder gene

Large scale human genetic studies have shown that loss of function (LoF) mutations in MYT1L are implicated in neurodevelopmental disorders (NDDs). Here, we provide an overview of the growing number of published MYT1L patient cases, and summarize prior studies in cells, zebrafish, and mice, both to understand MYT1L's molecular and cellular role during brain development and consider how its dysfunction can lead to NDDs. We integrate the conclusions from these studies and highlight conflicting findings to reassess the current model of the role of MYT1L as a transcriptional activator and/or repressor based on the biological context. Finally, we highlight additional functional studies that are needed to understand the molecular mechanisms underlying pathophysiology and propose key questions to guide future preclinical studies.

Characterization of regulatory transcriptional mechanisms in hepatocyte lipotoxicity

Non-alcoholic fatty liver disease is a continuum of disorders among which non-alcoholic steatohepatitis (NASH) is particularly associated with a negative prognosis. Hepatocyte lipotoxicity is one of the main pathogenic factors of liver fibrosis and NASH. However, the molecular mechanisms regulating this process are poorly understood. The main aim of this study was to dissect transcriptional mechanisms regulated by lipotoxicity in hepatocytes. We achieved this aim by combining transcriptomic, proteomic and chromatin accessibility analyses from human liver and mouse hepatocytes. This integrative approach revealed several transcription factor networks deregulated by NASH and lipotoxicity. To validate these predictions, genetic deletion of the transcription factors MAFK and TCF4 was performed, resulting in hepatocytes that were better protected against saturated fatty acid oversupply. MAFK- and TCF4-regulated gene expression profiles suggest a mitigating effect against cell stress, while promoting cell survival and growth. Moreover, in the context of lipotoxicity, some MAFK and TCF4 target genes were to the corresponding differentially regulated transcripts in human liver fibrosis. Collectively, our findings comprehensively profile the transcriptional response to lipotoxicity in hepatocytes, revealing new molecular insights and providing a valuable resource for future endeavours to tackle the molecular mechanisms of NASH.

CROWN ROOTLESS1 binds DNA with a relaxed specificity and activates OsROP and OsbHLH044 genes involved in crown root formation in rice

In cereals, the root system is mainly composed of post-embryonic shoot-borne roots, named crown roots. The CROWN ROOTLESS1 (CRL1) transcription factor, belonging to the ASYMMETRIC LEAVES2-LIKE/LATERAL ORGAN BOUNDARIES DOMAIN (ASL/LBD) family, is a key regulator of crown root initiation in rice (Oryza sativa). Here, we show that CRL1 can bind, both in vitro and in vivo, not only the LBD-box, a DNA sequence recognized by several ASL/LBD transcription factors, but also another not previously identified DNA motif that was named CRL1-box. Using rice protoplast transient transactivation assays and a set of previously identified CRL1-regulated genes, we confirm that CRL1 transactivates these genes if they possess at least a CRL1-box or an LBD-box in their promoters. In planta, ChIP-qPCR experiments targeting two of these genes that include both a CRL1- and an LBD-box in their promoter show that CRL1 binds preferentially to the LBD-box in these promoter contexts. CRISPR/Cas9-targeted mutation of these two CRL1-regulated genes, which encode a plant Rho GTPase (OsROP) and a basic helix-loop-helix transcription factor (OsbHLH044), show that both promote crown root development. Finally, we show that OsbHLH044 represses a regulatory module, uncovering how CRL1 regulates specific processes during crown root formation.

Recognition of Microbe- and Damage-Associated Molecular Patterns by Leucine-Rich Repeat Pattern Recognition Receptor Kinases Confers Salt Tolerance in Plants

In plants, a first layer of inducible immunity is conferred by pattern recognition receptors (PRRs) that bind microbe- and damage-associated molecular patterns to activate pattern-triggered immunity (PTI). PTI is strengthened or followed by another potent form of immunity when intracellular receptors recognize pathogen effectors, termed effector-triggered immunity. Immunity signaling regulators have been reported to influence abiotic stress responses as well, yet the governing principles and mechanisms remain ambiguous. Here, we report that PRRs of a leucine-rich repeat ectodomain also confer salt tolerance in Arabidopsis thaliana, following recognition of cognate ligands such as bacterial flagellin (flg22 epitope) and elongation factor Tu (elf18 epitope), and the endogenous Pep peptides. Pattern-triggered salt tolerance (PTST) requires authentic PTI signaling components; namely, the PRR-associated kinases BAK1 and BIK1 and the NADPH oxidase RBOHD. Exposure to salt stress induces the release of Pep precursors, pointing to the involvement of the endogenous immunogenic peptides in developing plant tolerance to high salinity. Transcriptome profiling reveals an inventory of PTST target genes, which increase or acquire salt responsiveness following a preexposure to immunogenic patterns. In good accordance, plants challenged with nonpathogenic bacteria also acquired salt tolerance in a manner dependent on PRRs. Our findings provide insight into signaling plasticity underlying biotic or abiotic stress cross-tolerance in plants conferred by PRRs.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Chromatin accessibility shapes meiotic recombination in mouse primordial germ cells through assisting double-strand breaks and loop formation

Meiotic recombination is a driver of evolution, and aberrant recombination is a major contributor to aneuploidy in mammals. Mechanism of recombination remains elusive yet. Here, we present a computational analysis to explore recombination-related dynamics of chromatin accessibility in mouse primordial germ cells (PGCs). Our data reveals that: (1) recombination hotspots which get accessible at meiosis-specific DNase I-hypersensitive sites (DHSs) only when PGCs enter meiosis are located preferentially in intronic and distal intergenic regions; (2) stable DHSs maintained stably across PGC differentiation are enriched by CTCF motifs and CTCF binding and mediate chromatin loop formation; (3) compared with the specific DHSs aroused at meiotic stage, stable DHSs are largely encoded in DNA sequence and also enriched by epigenetic marks; (4) PRDM9 is likely to target nucleosome-occupied hotspot regions and remodels local chromatin structure to make them accessible for recombination machinery; and (5) cells undergoing meiotic recombination are deficient in TAD structure and chromatin loop arrays are organized regularly along the axis formed between homologous chromosomes. Taken together, by analyzing DHS-related DNA features, epigenetic marks and 3D genome structure, we revealed some specific roles of chromatin accessibility in recombination, which would expand our understanding of recombination mechanism.

scEpiLock: A Weakly Supervised Learning Framework for cis-Regulatory Element Localization and Variant Impact Quantification for Single-Cell Epigenetic Data

Recent advances in single-cell transposase-accessible chromatin using a sequencing assay (scATAC-seq) allow cellular heterogeneity dissection and regulatory landscape reconstruction with an unprecedented resolution. However, compared to bulk-sequencing, its ultra-high missingness remarkably reduces usable reads in each cell type, resulting in broader, fuzzier peak boundary definitions and limiting our ability to pinpoint functional regions and interpret variant impacts precisely. We propose a weakly supervised learning method, scEpiLock, to directly identify core functional regions from coarse peak labels and quantify variant impacts in a cell-type-specific manner. First, scEpiLock uses a multi-label classifier to predict chromatin accessibility via a deep convolutional neural network. Then, its weakly supervised object detection module further refines the peak boundary definition using gradient-weighted class activation mapping (Grad-CAM). Finally, scEpiLock provides cell-type-specific variant impacts within a given peak region. We applied scEpiLock to various scATAC-seq datasets and found that it achieves an area under receiver operating characteristic curve (AUC) of ~0.9 and an area under precision recall (AUPR) above 0.7. Besides, scEpiLock’s object detection condenses coarse peaks to only ⅓ of their original size while still reporting higher conservation scores. In addition, we applied scEpiLock on brain scATAC-seq data and reported several genome-wide association studies (GWAS) variants disrupting regulatory elements around known risk genes for Alzheimer’s disease, demonstrating its potential to provide cell-type-specific biological insights in disease studies.