Conservation of transcription factor binding specificities across 600 million years of bilateria evolution

Sugar-responsive inhibition of Myc-dependent ribosome biogenesis by Clockwork orange

The ability to feed on a sugar-containing diet depends on a gene regulatory network controlled by the intracellular sugar sensor Mondo/ChREBP-Mlx, which remains insufficiently characterized. Here, we present a genome-wide temporal clustering of sugar-responsive gene expression in Drosophila larvae. We identify gene expression programs responding to sugar feeding, including downregulation of ribosome biogenesis genes, known targets of Myc. Clockwork orange (CWO), a component of the circadian clock, is found to be a mediator of this repressive response and to be necessary for survival on a high-sugar diet. CWO expression is directly activated by Mondo-Mlx, and it counteracts Myc through repression of its gene expression and through binding to overlapping genomic regions. CWO mouse ortholog BHLHE41 has a conserved role in repressing ribosome biogenesis genes in primary hepatocytes. Collectively, our data uncover a cross-talk between conserved gene regulatory circuits balancing the activities of anabolic pathways to maintain homeostasis during sugar feeding.

Combinatorial transcription factor binding encodes cis-regulatory wiring of forebrain GABAergic neurogenesis

Transcription factors (TFs) bind combinatorially to genomic cis-regulatory elements (cREs), orchestrating transcription programs. While studies of chromatin state and chromosomal interactions have revealed dynamic neurodevelopmental cRE landscapes, parallel understanding of the underlying TF binding lags. To elucidate the combinatorial TF-cRE interactions driving mouse basal ganglia development, we integrated ChIP-seq for twelve TFs, H3K4me3-associated enhancer-promoter interactions, chromatin and transcriptional state, and transgenic enhancer assays. We identified TF-cREs modules with distinct chromatin features and enhancer activity that have complementary roles driving GABAergic neurogenesis and suppressing other developmental fates. While the majority of distal cREs were bound by one or two TFs, a small proportion were extensively bound, and these enhancers also exhibited exceptional evolutionary conservation, motif density, and complex chromosomal interactions. Our results provide new insights into how modules of combinatorial TF-cRE interactions activate and repress developmental expression programs and demonstrate the value of TF binding data in modeling gene regulatory wiring.

ExplaiNN: interpretable and transparent neural networks for genomics

Deep learning models such as convolutional neural networks (CNNs) excel in genomic tasks but lack interpretability. We introduce ExplaiNN, which combines the expressiveness of CNNs with the interpretability of linear models. ExplaiNN can predict TF binding, chromatin accessibility, and de novo motifs, achieving performance comparable to state-of-the-art methods. Its predictions are transparent, providing global (cell state level) as well as local (individual sequence level) biological insights into the data. ExplaiNN can serve as a plug-and-play platform for pretrained models and annotated position weight matrices. ExplaiNN aims to accelerate the adoption of deep learning in genomic sequence analysis by domain experts.

The Drosophila gene encoding JIG protein (CG14850) is critical for CrebA nuclear trafficking during development

Coordination of mitochondrial and nuclear processes is key to the cellular health; however, very little is known about the molecular mechanisms regulating nuclear-mitochondrial crosstalk. Here, we report a novel molecular mechanism controlling the shuttling of CREB (cAMP response element-binding protein) protein complex between mitochondria and nucleoplasm. We show that a previously unknown protein, herein termed as Jig, functions as a tissue-specific and developmental timing-specific coregulator in the CREB pathway. Our results demonstrate that Jig shuttles between mitochondria and nucleoplasm, interacts with CrebA protein and controls its delivery to the nucleus, thus triggering CREB-dependent transcription in nuclear chromatin and mitochondria. Ablating the expression of Jig prevents CrebA from localizing to the nucleoplasm, affecting mitochondrial functioning and morphology and leads to Drosophila developmental arrest at the early third instar larval stage. Together, these results implicate Jig as an essential mediator of nuclear and mitochondrial processes. We also found that Jig belongs to a family of nine similar proteins, each of which has its own tissue- and time-specific expression profile. Thus, our results are the first to describe the molecular mechanism regulating nuclear and mitochondrial processes in a tissue- and time-specific manner.

Liver transcriptomic and methylomic analyses identify transcriptional mitogen-activated protein kinase regulation in facultative hibernation of Syrian hamster

Hibernation consists of alternating torpor-arousal phases, during which animals cope with repetitive hypothermia and ischaemia-reperfusion. Due to limited transcriptomic and methylomic information for facultative hibernators, we here conducted RNA and whole-genome bisulfide sequencing in liver of hibernating Syrian hamster (Mesocricetus auratus). Gene ontology analysis was performed on 844 differentially expressed genes and confirmed the shift in metabolic fuel utilization, inhibition of RNA transcription and cell cycle regulation as found in seasonal hibernators. Additionally, we showed a so far unreported suppression of mitogen-activated protein kinase (MAPK) and protein phosphatase 1 pathways during torpor. Notably, hibernating hamsters showed upregulation of MAPK inhibitors (dual-specificity phosphatases and sproutys) and reduced levels of MAPK-induced transcription factors (TFs). Promoter methylation was found to modulate the expression of genes targeted by these TFs. In conclusion, we document gene regulation between hibernation phases, which may aid the identification of pathways and targets to prevent organ damage in transplantation or ischaemia-reperfusion.

Inhibition of Nonsense-Mediated Decay Induces Nociceptive Sensitization through Activation of the Integrated Stress Response

RNA stability is meticulously controlled. Here, we sought to determine whether an essential post-transcriptional regulatory mechanism plays a role in pain. Nonsense-mediated decay (NMD) safeguards against translation of mRNAs that harbor premature termination codons and controls the stability of ∼10% of typical protein-coding mRNAs. It hinges on the activity of the conserved kinase SMG1. Both SMG1 and its target, UPF1, are expressed in murine DRG sensory neurons. SMG1 protein is present in both the DRG and sciatic nerve. Using high-throughput sequencing, we examined changes in mRNA abundance following inhibition of SMG1. We confirmed multiple NMD stability targets in sensory neurons, including ATF4. ATF4 is preferentially translated during the integrated stress response (ISR). This led us to ask whether suspension of NMD induces the ISR. Inhibition of NMD increased eIF2-α phosphorylation and reduced the abundance of the eIF2-α phosphatase constitutive repressor of eIF2-α phosphorylation. Finally, we examined the effects of SMG1 inhibition on pain-associated behaviors. Peripheral inhibition of SMG1 results in mechanical hypersensitivity in males and females that persists for several days and priming to a subthreshold dose of PGE2. Priming was fully rescued by a small-molecule inhibitor of the ISR. Collectively, our results indicate that suspension of NMD promotes pain through stimulation of the ISR.SIGNIFICANCE STATEMENT Nociceptors undergo long-lived changes in their plasticity which may contribute to chronic pain. Translational regulation has emerged as a dominant mechanism in pain. Here, we investigate the role of a major pathway of RNA surveillance called nonsense-mediated decay (NMD). Modulation of NMD is potentially beneficial for a broad array of diseases caused by frameshift or nonsense mutations. Our results suggest that inhibition of the rate-limiting step of NMD drives behaviors associated with pain through activation of the ISR. This work reveals complex interconnectivity between RNA stability and translational regulation and suggests an important consideration in harnessing the salubrious benefits of NMD disruption.

Multiomics Analysis Reveals Novel Genetic Determinants for Lens Differentiation, Structure, and Transparency

Recent advances in next-generation sequencing and data analysis have provided new gateways for identification of novel genome-wide genetic determinants governing tissue development and disease. These advances have revolutionized our understanding of cellular differentiation, homeostasis, and specialized function in multiple tissues. Bioinformatic and functional analysis of these genetic determinants and the pathways they regulate have provided a novel basis for the design of functional experiments to answer a wide range of long-sought biological questions. A well-characterized model for the application of these emerging technologies is the development and differentiation of the ocular lens and how individual pathways regulate lens morphogenesis, gene expression, transparency, and refraction. Recent applications of next-generation sequencing analysis on well-characterized chicken and mouse lens differentiation models using a variety of omics techniques including RNA-seq, ATAC-seq, whole-genome bisulfite sequencing (WGBS), chip-seq, and CUT&RUN have revealed a wide range of essential biological pathways and chromatin features governing lens structure and function. Multiomics integration of these data has established new gene functions and cellular processes essential for lens formation, homeostasis, and transparency including the identification of novel transcription control pathways, autophagy remodeling pathways, and signal transduction pathways, among others. This review summarizes recent omics technologies applied to the lens, methods for integrating multiomics data, and how these recent technologies have advanced our understanding ocular biology and function. The approach and analysis are relevant to identifying the features and functional requirements of more complex tissues and disease states.

Ratiometric sensing of Pnt and Yan transcription factor levels confers ultrasensitivity to photoreceptor fate transitions in Drosophila

Cell state transitions are often triggered by large changes in the concentrations of transcription factors and therefore large differences in their stoichiometric ratios. Whether cells can elicit transitions using modest changes in the ratios of co-expressed factors is unclear. Here, we investigate how cells in the Drosophila eye resolve state transitions by quantifying the expression dynamics of the ETS transcription factors Pnt and Yan. Eye progenitor cells maintain a relatively constant ratio of Pnt/Yan protein, despite expressing both proteins with pulsatile dynamics. A rapid and sustained twofold increase in the Pnt/Yan ratio accompanies transitions to photoreceptor fates. Genetic perturbations that modestly disrupt the Pnt/Yan ratio produce fate transition defects consistent with the hypothesis that transitions are normally driven by a twofold shift in the ratio. A biophysical model based on cooperative Yan-DNA binding coupled with non-cooperative Pnt-DNA binding illustrates how twofold ratio changes could generate ultrasensitive changes in target gene transcription to drive fate transitions. Thus, coupling cell state transitions to the Pnt/Yan ratio sensitizes the system to modest fold-changes, conferring robustness and ultrasensitivity to the developmental program.

A modERN resource: identification of Drosophila transcription factor candidate target genes using RNAi

Transcription factors (TFs) play a key role in development and in cellular responses to the environment by activating or repressing the transcription of target genes in precise spatial and temporal patterns. In order to develop a catalog of target genes of Drosophila melanogaster TFs, the modERN consortium systematically knocked down the expression of TFs using RNAi in whole embryos followed by RNA-seq. We generated data for 45 TFs which have 18 different DNA-binding domains and are expressed in 15 of the 16 organ systems. The range of inactivation of the targeted TFs by RNAi ranged from log2fold change -3.52 to +0.49. The TFs also showed remarkable heterogeneity in the numbers of candidate target genes identified, with some generating thousands of candidates and others only tens. We present detailed analysis from five experiments, including those for three TFs that have been the focus of previous functional studies (ERR, sens, and zfh2) and two previously uncharacterized TFs (sens-2 and CG32006), as well as short vignettes for selected additional experiments to illustrate the utility of this resource. The RNA-seq datasets are available through the ENCODE DCC (http://encodeproject.org) and the Sequence Read Archive (SRA). TF and target gene expression patterns can be found here: https://insitu.fruitfly.org. These studies provide data that facilitate scientific inquiries into the functions of individual TFs in key developmental, metabolic, defensive, and homeostatic regulatory pathways, as well as provide a broader perspective on how individual TFs work together in local networks during embryogenesis.

X-chromosome target specificity diverged between dosage compensation mechanisms of two closely related Caenorhabditis species

An evolutionary perspective enhances our understanding of biological mechanisms. Comparison of sex determination and X-chromosome dosage compensation mechanisms between the closely related nematode species Caenorhabditis briggsae (Cbr) and Caenorhabditis elegans (Cel) revealed that the genetic regulatory hierarchy controlling both processes is conserved, but the X-chromosome target specificity and mode of binding for the specialized condensin dosage compensation complex (DCC) controlling X expression have diverged. We identified two motifs within Cbr DCC recruitment sites that are highly enriched on X: 13 bp MEX and 30 bp MEX II. Mutating either MEX or MEX II in an endogenous recruitment site with multiple copies of one or both motifs reduced binding, but only removing all motifs eliminated binding in vivo. Hence, DCC binding to Cbr recruitment sites appears additive. In contrast, DCC binding to Cel recruitment sites is synergistic: mutating even one motif in vivo eliminated binding. Although all X-chromosome motifs share the sequence CAGGG, they have otherwise diverged so that a motif from one species cannot function in the other. Functional divergence was demonstrated in vivo and in vitro. A single nucleotide position in Cbr MEX can determine whether Cel DCC binds. This rapid divergence of DCC target specificity could have been an important factor in establishing reproductive isolation between nematode species and contrasts dramatically with the conservation of target specificity for X-chromosome dosage compensation across Drosophila species and for transcription factors controlling developmental processes such as body-plan specification from fruit flies to mice.

Diverse logics and grammar encode notochord enhancers

The notochord is a defining feature of all chordates. The transcription factors Zic and ETS regulate enhancer activity within the notochord. We conduct high-throughput screens of genomic elements within developing Ciona embryos to understand how Zic and ETS sites encode notochord activity. Our screen discovers an enhancer located near Lama, a gene critical for notochord development. Reversing the orientation of an ETS site within this enhancer abolishes expression, indicating that enhancer grammar is critical for notochord activity. Similarly organized clusters of Zic and ETS sites occur within mouse and human Lama1 introns. Within a Brachyury (Bra) enhancer, FoxA and Bra, in combination with Zic and ETS binding sites, are necessary and sufficient for notochord expression. This binding site logic also occurs within other Ciona and vertebrate Bra enhancers. Collectively, this study uncovers the importance of grammar within notochord enhancers and discovers signatures of enhancer logic and grammar conserved across chordates.

Computational approaches to understand transcription regulation in development

Gene regulatory networks (GRNs) serve as useful abstractions to understand transcriptional dynamics in developmental systems. Computational prediction of GRNs has been successfully applied to genome-wide gene expression measurements with the advent of microarrays and RNA-sequencing. However, these inferred networks are inaccurate and mostly based on correlative rather than causative interactions. In this review, we highlight three approaches that significantly impact GRN inference: (1) moving from one genome-wide functional modality, gene expression, to multi-omics, (2) single cell sequencing, to measure cell type-specific signals and predict context-specific GRNs, and (3) neural networks as flexible models. Together, these experimental and computational developments have the potential to significantly impact the quality of inferred GRNs. Ultimately, accurately modeling the regulatory interactions between transcription factors and their target genes will be essential to understand the role of transcription factors in driving developmental gene expression programs and to derive testable hypotheses for validation.

Current Status of Oligonucleotide-Based Protein Degraders

Transcription factors (TFs) and RNA-binding proteins (RBPs) have long been considered undruggable, mainly because they lack ligand-binding sites and are equipped with flat and narrow protein surfaces. Protein-specific oligonucleotides have been harnessed to target these proteins with some satisfactory preclinical results. The emerging proteolysis-targeting chimera (PROTAC) technology is no exception, utilizing protein-specific oligonucleotides as warheads to target TFs and RBPs. In addition, proteolysis by proteases is another type of protein degradation. In this review article, we discuss the current status of oligonucleotide-based protein degraders that are dependent either on the ubiquitin–proteasome system or a protease, providing a reference for the future development of degraders.

Investigation of autism-related transcription factors underlying sex differences in the effects of bisphenol A on transcriptome profiles and synaptogenesis in the offspring hippocampus

BACKGROUND

Bisphenol A (BPA) has been linked to susceptibility to autism spectrum disorder (ASD). Our recent studies have shown that prenatal BPA exposure disrupted ASD-related gene expression in the hippocampus, neurological functions, and behaviors associated with ASD in a sex-specific pattern. However, the molecular mechanisms underlying the effects of BPA are still unclear.

METHODS

Transcriptome data mining and molecular docking analyses were performed to identify ASD-related transcription factors (TFs) and their target genes underlying the sex-specific effects of prenatal BPA exposure. Gene ontology analysis was conducted to predict biological functions associated with these genes. The expression levels of ASD-related TFs and targets in the hippocampus of rat pups prenatally exposed to BPA were measured using qRT-PCR analysis. The role of the androgen receptor (AR) in BPA-mediated regulation of ASD candidate genes was investigated using a human neuronal cell line stably transfected with AR-expression or control plasmid. Synaptogenesis, which is a function associated with genes transcriptionally regulated by ASD-related TFs, was assessed using primary hippocampal neurons isolated from male and female rat pups prenatally exposed to BPA.

RESULTS

We found that there was a sex difference in ASD-related TFs underlying the effects of prenatal BPA exposure on the transcriptome profiles of the offspring hippocampus. In addition to the known BPA targets AR and ESR1, BPA could directly interact with novel targets (i.e., KDM5B, SMAD4, and TCF7L2). The targets of these TFs were also associated with ASD. Prenatal BPA exposure disrupted the expression of ASD-related TFs and targets in the offspring hippocampus in a sex-dependent manner. Moreover, AR was involved in the BPA-mediated dysregulation of AUTS2, KMT2C, and SMARCC2. Prenatal BPA exposure altered synaptogenesis by increasing synaptic protein levels in males but not in females, but the number of excitatory synapses was increased in female primary neurons only.

CONCLUSIONS

Our findings suggest that AR and other ASD-related TFs are involved in sex differences in the effects of prenatal BPA exposure on transcriptome profiles and synaptogenesis in the offspring hippocampus. These TFs may play an essential role in an increased ASD susceptibility associated with endocrine-disrupting chemicals, particularly BPA, and the male bias of ASD.

Leveraging massively parallel reporter assays for evolutionary questions

A long-standing goal of evolutionary biology is to decode how gene regulation contributes to organismal diversity. Doing so is challenging because it is hard to predict function from non-coding sequence and to perform molecular research with non-model taxa. Massively parallel reporter assays (MPRAs) enable the testing of thousands to millions of sequences for regulatory activity simultaneously. Here, we discuss the execution, advantages, and limitations of MPRAs, with a focus on evolutionary questions. We propose solutions for extending MPRAs to rare taxa and those with limited genomic resources, and we underscore MPRA's broad potential for driving genome-scale, functional studies across organisms.

Epigenetic and Genetic Population Structure is Coupled in a Marine Invertebrate

Delineating the relative influence of genotype and the environment on DNA methylation is critical for characterizing the spectrum of organism fitness as driven by adaptation and phenotypic plasticity. In this study, we integrated genomic and DNA methylation data for two distinct Olympia oyster (Ostrea lurida) populations while controlling for within-generation environmental influences. In addition to providing the first characterization of genome-wide DNA methylation patterns in the oyster genus Ostrea, we identified 3,963 differentially methylated loci between populations. Our results show a clear coupling between genetic and epigenetic patterns of variation, with 27% of variation in interindividual methylation differences explained by genotype. Underlying this association are both direct genetic changes in CpGs (CpG-SNPs) and genetic variation with indirect influence on methylation (mQTLs). When comparing measures of genetic and epigenetic population divergence at specific genomic regions this relationship surprisingly breaks down, which has implications for the methods commonly used to study epigenetic and genetic coupling in marine invertebrates.

Chromatin accessibility-based characterisation of brain gene regulatory networks in three distinct honey bee polyphenisms

The honey bee genome has the capacity to produce three phenotypically distinct organisms (two diploid female castes: queen and worker, and a haploid male drone). Previous studies have implicated metabolic flux acting via epigenetic regulation in directing nutrition-driven phenotypic plasticity in the honey bee. However, the cis-acting DNA regulatory elements that establish tissue and polyphenism -specific epigenomes and gene expression programmes, remain unclear. Using a high resolution multiomic approach including assay for transposase-accessible chromatin by sequencing (ATAC-seq), RNA-seq and ChIP-seq, we produce the first genome-wide maps of the regulatory landscape across all three adult honey bee phenotypes identifying > 5000 regulatory regions in queen, 7500 in worker and 6500 in drone, with the vast majority of these sites located within intronic regions. These regions are defined by positive enrichment of H3K27ac and depletion of H3K4me3 and show a positive correlation with gene expression. Using ATAC-seq footprinting we determine queen, worker and drone -specific transcription factor occupancy and uncover novel phenotype-specific regulatory networks identifying two key nuclear receptors that have previously been implicated in caste-determination and adult behavioural maturation in honey bees; ecdysone receptor and ultraspiracle. Collectively, this study provides novel insights into key gene regulatory networks that are associated with these distinct polyphenisms in the honey bee.

AnimalTFDB 4.0: a comprehensive animal transcription factor database updated with variation and expression annotations

Transcription factors (TFs) are proteins that interact with specific DNA sequences to regulate gene expression and play crucial roles in all kinds of biological processes. To keep up with new data and provide a more comprehensive resource for TF research, we updated the Animal Transcription Factor Database (AnimalTFDB) to version 4.0 (http://bioinfo.life.hust.edu.cn/AnimalTFDB4/) with up-to-date data and functions. We refined the TF family rules and prediction pipeline to predict TFs in genome-wide protein sequences from Ensembl. As a result, we predicted 274 633 TF genes and 150 726 transcription cofactor genes in AnimalTFDB 4.0 in 183 animal genomes, which are 86 more species than AnimalTFDB 3.0. Besides double data volume, we also added the following new annotations and functions to the database: (i) variations (including mutations) on TF genes in various human cancers and other diseases; (ii) predicted post-translational modification sites (including phosphorylation, acetylation, methylation and ubiquitination sites) on TFs in 8 species; (iii) TF regulation in autophagy; (iv) comprehensive TF expression annotation for 38 species; (v) exact and batch search functions allow users to search AnimalTFDB flexibly. AnimalTFDB 4.0 is a useful resource for studying TF and transcription regulation, which contains comprehensive annotation and classification of TFs and transcription cofactors.

Prediction of protein-ligand binding affinity from sequencing data with interpretable machine learning

Protein-ligand interactions are increasingly profiled at high throughput using affinity selection and massively parallel sequencing. However, these assays do not provide the biophysical parameters that most rigorously quantify molecular interactions. Here we describe a flexible machine learning method, called ProBound, that accurately defines sequence recognition in terms of equilibrium binding constants or kinetic rates. This is achieved using a multi-layered maximum-likelihood framework that models both the molecular interactions and the data generation process. We show that ProBound quantifies transcription factor (TF) behavior with models that predict binding affinity over a range exceeding that of previous resources; captures the impact of DNA modifications and conformational flexibility of multi-TF complexes; and infers specificity directly from in vivo data such as ChIP-seq without peak calling. When coupled with an assay called KD-seq, it determines the absolute affinity of protein-ligand interactions. We also apply ProBound to profile the kinetics of kinase-substrate interactions. ProBound opens new avenues for decoding biological networks and rationally engineering protein-ligand interactions.

Functional characterization of enhancer activity during a long terminal repeat's evolution

Many transposable elements (TEs) contain transcription factor binding sites and are implicated as potential regulatory elements. However, TEs are rarely functionally tested for regulatory activity, which in turn limits our understanding of how TE regulatory activity has evolved. We systematically tested the human LTR18A subfamily for regulatory activity using massively parallel reporter assay (MPRA) and found AP-1- and CEBP-related binding motifs as drivers of enhancer activity. Functional analysis of evolutionarily reconstructed ancestral sequences revealed that LTR18A elements have generally lost regulatory activity over time through sequence changes, with the largest effects occurring owing to mutations in the AP-1 and CEBP motifs. We observed that the two motifs are conserved at higher rates than expected based on neutral evolution. Finally, we identified LTR18A elements as potential enhancers in the human genome, primarily in epithelial cells. Together, our results provide a model for the origin, evolution, and co-option of TE-derived regulatory elements.

Learning probabilistic protein-DNA recognition codes from DNA-binding specificities using structural mappings

Knowledge of how proteins interact with DNA is essential for understanding gene regulation. Although DNA-binding specificities for thousands of transcription factors (TFs) have been determined, the specific amino acid-base interactions comprising their structural interfaces are largely unknown. This lack of resolution hampers attempts to leverage these data in order to predict specificities for uncharacterized TFs or TFs mutated in disease. Here we introduce recognition code learning via automated mapping of protein-DNA structural interfaces (rCLAMPS), a probabilistic approach that uses DNA-binding specificities for TFs from the same structural family to simultaneously infer both which nucleotide positions are contacted by particular amino acids within the TF as well as a recognition code that relates each base-contacting amino acid to nucleotide preferences at the DNA positions it contacts. We apply rCLAMPS to homeodomains, the second largest family of TFs in metazoans and show that it learns a highly effective recognition code that can predict de novo DNA-binding specificities for TFs. Furthermore, we show that the inferred amino acid-nucleotide contacts reveal whether and how nucleotide preferences at individual binding site positions are altered by mutations within TFs. Our approach is an important step toward automatically uncovering the determinants of protein-DNA specificity from large compendia of DNA-binding specificities and inferring the altered functionalities of TFs mutated in disease.

High-throughput sequencing SELEX for the determination of DNA-binding protein specificities in vitro

High-throughput sequencing SELEX (HT-SELEX) is a powerful technique for unbiased determination of preferred target motifs of DNA-binding proteins in vitro. The procedure depends upon selection of DNA binding sites from a random library of oligonucleotides by purifying protein-DNA complexes and amplifying bound DNA using the polymerase chain reaction. Here, we describe an optimized step-by-step protocol for HT-SELEX compatible with Illumina sequencing. We also introduce a bioinformatic pipeline (eme_selex) facilitating the detection of promiscuous DNA binding by analyzing the enrichment of all possible k-mers. For complete details on the use and execution of this protocol, please refer to Pantier et al. (2021).

The highly conserved RNA-binding specificity of nucleocapsid protein facilitates the identification of drugs with broad anti-coronavirus activity

The binding of SARS-CoV-2 nucleocapsid (N) protein to both the 5'- and 3'-ends of genomic RNA has different implications arising from its binding to the central region during virion assembly. However, the mechanism underlying selective binding remains unknown. Herein, we performed the high-throughput RNA-SELEX (HTR-SELEX) to determine the RNA-binding specificity of the N proteins of various SARS-CoV-2 variants as well as other β-coronaviruses and showed that N proteins could bind two unrelated sequences, both of which were highly conserved across all variants and species. Interestingly, both sequences are virtually absent from the human transcriptome; however, they exhibit a highly enriched, mutually complementary distribution in the coronavirus genome, highlighting their varied functions in genome packaging. Our results provide mechanistic insights into viral genome packaging, thereby increasing the feasibility of developing drugs with broad-spectrum anti-coronavirus activity by targeting RNA binding by N proteins.

The regulatory landscape of neurite development in Caenorhabditis elegans

Neuronal communication requires precise connectivity of neurite projections (axons and dendrites). Developing neurites express cell-surface receptors that interpret extracellular cues to enable correct guidance toward, and connection with, target cells. Spatiotemporal regulation of neurite guidance molecule expression by transcription factors (TFs) is critical for nervous system development and function. Here, we review how neurite development is regulated by TFs in the Caenorhabditis elegans nervous system. By collecting publicly available transcriptome and ChIP-sequencing data, we reveal gene expression dynamics during neurite development, providing insight into transcriptional mechanisms governing construction of the nervous system architecture.

Enhancer RNAs step forward: new insights into enhancer function

Enhancers confer precise spatiotemporal patterns of gene expression in response to developmental and environmental stimuli. Over the last decade, the transcription of enhancer RNAs (eRNAs) – nascent RNAs transcribed from active enhancers – has emerged as a key factor regulating enhancer activity. eRNAs are relatively short-lived RNA species that are transcribed at very high rates but also quickly degraded. Nevertheless, eRNAs are deeply intertwined within enhancer regulatory networks and are implicated in a number of transcriptional control mechanisms. Enhancers show changes in function and sequence over evolutionary time, raising questions about the relationship between enhancer sequences and eRNA function. Moreover, the vast majority of single nucleotide polymorphisms associated with human complex diseases map to the non-coding genome, with causal disease variants enriched within enhancers. In this Primer, we survey the diverse roles played by eRNAs in enhancer-dependent gene expression, evaluating different models for eRNA function. We also explore questions surrounding the genetic conservation of enhancers and how this relates to eRNA function and dysfunction.

Summary: This Primer evaluates the ideas that underpin developing models for eRNA function, exploring cases in which perturbed eRNA function contributes to disease.

The transcription factor FoxP3 can fold into two dimerization states with divergent implications for regulatory T cell function and immune homeostasis

FoxP3 is an essential transcription factor (TF) for immunologic homeostasis, but how it utilizes the common forkhead DNA-binding domain (DBD) to perform its unique function remains poorly understood. We here demonstrated that unlike other known forkhead TFs, FoxP3 formed a head-to-head dimer using a unique linker (Runx1-binding region [RBR]) preceding the forkhead domain. Head-to-head dimerization conferred distinct DNA-binding specificity and created a docking site for the cofactor Runx1. RBR was also important for proper folding of the forkhead domain, as truncation of RBR induced domain-swap dimerization of forkhead, which was previously considered the physiological form of FoxP3. Rather, swap-dimerization impaired FoxP3 function, as demonstrated with the disease-causing mutation R337Q, whereas a swap-suppressive mutation largely rescued R337Q-mediated functional impairment. Altogether, our findings suggest that FoxP3 can fold into two distinct dimerization states: head-to-head dimerization representing functional specialization of an ancient DBD and swap dimerization associated with impaired functions.

Cyclin-dependent kinase-mediated phosphorylation and the negative regulatory domain of transcription factor B-Myb modulate its DNA binding

B-Myb is a highly conserved member of the vertebrate Myb family of transcription factors that plays a critical role in cell-cycle progression and proliferation. Myb proteins activate Myb-dependent promoters by interacting specifically with Myb-binding site (MBS) sequences using their DNA-binding domain (DBD). Transactivation of MBS promoters by B-Myb is repressed by its negative regulatory domain (NRD), and phosphorylation of the NRD by Cdk2-CyclinA relieves the repression to activate B-Myb–dependent promoters. However, the structural mechanisms underlying autoinhibition and activation of B-Myb–mediated transcription have been poorly characterized. Here, we determined that a region in the B-Myb NRD (residues 510–600) directly associates with the DBD and inhibits binding of the DBD to the MBS DNA sequence. We demonstrate using biophysical assays that phosphorylation of the NRD at T515, T518, and T520 is sufficient to disrupt the interaction between the NRD and the DBD, which results in increased affinity for MBS DNA and increased B-Myb–dependent promoter activation in cell assays. Our biochemical characterization of B-Myb autoregulation and the activating effects of phosphorylation provide insight into how B-Myb functions as a site-specific transcription factor.

Endothelial cell heterogeneity and microglia regulons revealed by a pig cell landscape at single-cell level

Pigs are valuable large animal models for biomedical and genetic research, but insights into the tissue- and cell-type-specific transcriptome and heterogeneity remain limited. By leveraging single-cell RNA sequencing, we generate a multiple-organ single-cell transcriptomic map containing over 200,000 pig cells from 20 tissues/organs. We comprehensively characterize the heterogeneity of cells in tissues and identify 234 cell clusters, representing 58 major cell types. In-depth integrative analysis of endothelial cells reveals a high degree of heterogeneity. We identify several functionally distinct endothelial cell phenotypes, including an endothelial to mesenchymal transition subtype in adipose tissues. Intercellular communication analysis predicts tissue- and cell type-specific crosstalk between endothelial cells and other cell types through the VEGF, PDGF, TGF-β, and BMP pathways. Regulon analysis of single-cell transcriptome of microglia in pig and 12 other species further identifies MEF2C as an evolutionally conserved regulon in the microglia. Our work describes the landscape of single-cell transcriptomes within diverse pig organs and identifies the heterogeneity of endothelial cells and evolutionally conserved regulon in microglia.

Competition between Pseudomonas aeruginosa and Staphylococcus aureus is dependent on intercellular signaling and regulated by the NtrBC two-component system

Pseudomonas aeruginosa and Staphylococcus aureus are often comorbid human pathogens, isolated from expectorated sputum of cystic fibrosis patients and chronically infected wounds. Prior studies revealed a competitive advantage of P. aeruginosa over S. aureus in vitro that was slightly muted in vivo. Here, we demonstrated that the two-component regulatory system NtrBC influences the competitive advantage of P. aeruginosa over S. aureus in skin organoid and mouse models of co-infection. Expression of ntrBC was induced during co-culture of the two species and could be recapitulated in monoculture by the addition of the metabolite N-acetylglucosamine that is released from S. aureus following lysis. P. aeruginosa LESB58 WT, but not mutant (ΔntrC and ΔntrBC) strains, induced lysis of S. aureus USA300 LAC during planktonic growth and outcompeted S. aureus USA300 LAC during biofilm formation in vitro. We confirmed these findings in a murine abscess model of high-density infection. Accordingly, the secretory profile of P. aeruginosa LESB58 mutants revealed reduced production of anti-staphylococcal virulence factors including pyoverdine, pyocyanin and elastase. These phenotypes of LESB58 ΔntrBC could be at least partly complemented by overexpression of quorum sensing molecules including homoserine lactones or alkylquinolone signaling molecules. These data implicate the NtrBC two-component system in the complex regulatory cascade triggered by interspecies signaling that gives P. aeruginosa LESB58 a competitive edge over S. aureus USA300 LAC.

Upregulation of ribosome biogenesis via canonical E-boxes is required for Myc-driven proliferation

The transcription factor Myc drives cell growth across animal phyla and is activated in most forms of human cancer. However, it is unclear which Myc target genes need to be regulated to induce growth and whether multiple targets act additively or if induction of each target is individually necessary. Here, we identified Myc target genes whose regulation is conserved between humans and flies and deleted Myc-binding sites (E-boxes) in the promoters of fourteen of these genes in Drosophila. E-box mutants of essential genes were homozygous viable, indicating that the E-boxes are not required for basal expression. Eight E-box mutations led to Myc-like phenotypes; the strongest mutant, ppan^Ebox-/-, also made the flies resistant to Myc-induced cell growth without affecting Myc-induced apoptosis. The ppan^Ebox-/- flies are healthy and display only a minor developmental delay, suggesting that it may be possible to treat or prevent tumorigenesis by targeting individual downstream targets of Myc.

Gene regulatory networks of epidermal and neural fate choice in a chordate

Neurons are a highly specialized cell type only found in metazoans. They can be scattered throughout the body or grouped together, forming ganglia or nerve cords. During embryogenesis, centralized nervous systems develop from the ectoderm, which also forms the epidermis. How pluripotent ectodermal cells are directed toward neural or epidermal fates, and to which extent this process is shared among different animal lineages, are still open questions. Here, by using micromere explants, we were able to define in silico the putative gene regulatory networks (GRNs) underlying the first steps of the epidermis and the central nervous system formation in the cephalochordate amphioxus. We propose that although the signal triggering neural induction in amphioxus (i.e., Nodal) is different from vertebrates, the main transcription factors implicated in this process are conserved. Moreover, our data reveal that transcription factors of the neural program seem to not only activate neural genes but also to potentially have direct inputs into the epidermal GRN, suggesting that the Nodal signal might also contribute to neural fate commitment by repressing the epidermal program. Our functional data on whole embryos support this result and highlight the complex interactions among the transcription factors activated by the signaling pathways that drive ectodermal cell fate choice in chordates.

Sequence determinants of human gene regulatory elements

DNA can determine where and when genes are expressed, but the full set of sequence determinants that control gene expression is unknown. Here, we measured the transcriptional activity of DNA sequences that represent an ~100 times larger sequence space than the human genome using massively parallel reporter assays (MPRAs). Machine learning models revealed that transcription factors (TFs) generally act in an additive manner with weak grammar and that most enhancers increase expression from a promoter by a mechanism that does not appear to involve specific TF–TF interactions. The enhancers themselves can be classified into three types: classical, closed chromatin and chromatin dependent. We also show that few TFs are strongly active in a cell, with most activities being similar between cell types. Individual TFs can have multiple gene regulatory activities, including chromatin opening and enhancing, promoting and determining transcription start site (TSS) activity, consistent with the view that the TF binding motif is the key atomic unit of gene expression.

Analysis of massively parallel reporter assays measuring the transcriptional activity of DNA sequences indicates that most transcription factor (TF) activity is additive and does not rely on specific TF–TF interactions. Individual TFs can have different gene regulatory activities.

Distinct structural bases for sequence-specific DNA binding by mammalian BEN domain proteins

The BEN domain is a recently recognized DNA binding module that is present in diverse metazoans and certain viruses. Several BEN domain factors are known as transcriptional repressors, but, overall, relatively little is known of how BEN factors identify their targets in humans. In particular, X-ray structures of BEN domain:DNA complexes are only known for Drosophila factors bearing a single BEN domain, which lack direct vertebrate orthologs. Here, we characterize several mammalian BEN domain (BD) factors, including from two NACC family BTB-BEN proteins and from BEND3, which has four BDs. In vitro selection data revealed sequence-specific binding activities of isolated BEN domains from all of these factors. We conducted detailed functional, genomic, and structural studies of BEND3. We show that BD4 is a major determinant for in vivo association and repression of endogenous BEND3 targets. We obtained a high-resolution structure of BEND3-BD4 bound to its preferred binding site, which reveals how BEND3 identifies cognate DNA targets and shows differences with one of its non-DNA-binding BEN domains (BD1). Finally, comparison with our previous invertebrate BEN structures, along with additional structural predictions using AlphaFold2 and RoseTTAFold, reveal distinct strategies for target DNA recognition by different types of BEN domain proteins. Together, these studies expand the DNA recognition activities of BEN factors and provide structural insights into sequence-specific DNA binding by mammalian BEN proteins.

Rapid Biosensing Method for Detecting Protein-DNA Interactions

Identifying and investigating protein-DNA interactions, which play significant roles in many biological processes, is essential for basic and clinical research. Current techniques for identification of protein-DNA interactions are laborious, time-consuming, and suffer from nonspecific binding and limited sensitivity. To overcome these challenges and assess protein-DNA interactions, we use a magnetic modulation biosensing (MMB) system. In MMB, one of the interacting elements (protein or DNA) is immobilized to magnetic beads, and the other is coupled to a fluorescent molecule. Thus, the link between the magnetic bead and the fluorescent molecule is established only when binding occurs, enabling detection of the protein-DNA interaction. Using magnetic forces, the beads are concentrated and manipulated in a periodic motion in and out of a laser beam, producing a detectable oscillating signal. Using MMB, we detected protein-DNA interactions between short GC-rich DNA sequences and both a purified specificity protein 1 (Sp1) and an overexpressed Buttonhead (BTD) protein in a cell lysate. The specificity of the interactions was assessed using mutated DNA sequences and competition experiments. The assays were experimentally compared with commonly used electrophoretic mobility shift assay, which takes approximately 4-72 h. In comparison, the MMB-based assay's turnaround time is ∼2 h, and it provides unambiguous results and quantitative measures of performance. The MMB system uses simple and cheap components, making it an attractive alternative method over current costly and time-consuming techniques for analyzing protein-DNA interactions. Therefore, we anticipate that the MMB-based technique will significantly advance the detection of protein-DNA interactions in biomedical research.

Quantitative Dissection of the Proximal Ciona brachyury Enhancer

A major goal in biology is to understand the rules by which cis-regulatory sequences control spatially and temporally precise expression patterns. Here we present a systematic dissection of the proximal enhancer for the notochord-specific transcription factor brachyury in the ascidian chordate Ciona. The study uses a quantitative image-based reporter assay that incorporates a dual-reporter strategy to control for variable electroporation efficiency. We identified and mutated multiple predicted transcription factor binding sites of interest based on statistical matches to the JASPAR binding motif database. Most sites (Zic, Ets, FoxA, RBPJ) were selected based on prior knowledge of cell fate specification in both the primary and secondary notochord. We also mutated predicted Brachyury sites to investigate potential autoregulation as well as Fos/Jun (AP1) sites that had very strong matches to JASPAR. Our goal was to quantitatively define the relative importance of these different sites, to explore the importance of predicted high-affinity versus low-affinity motifs, and to attempt to design mutant enhancers that were specifically expressed in only the primary or secondary notochord lineages. We found that the mutation of all predicted high-affinity sites for Zic, FoxA or Ets led to quantifiably distinct effects. The FoxA construct caused a severe loss of reporter expression whereas the Ets construct had little effect. A strong Ets phenotype was only seen when much lower-scoring binding sites were also mutated. This supports the enhancer suboptimization hypothesis proposed by Farley and Levine but suggests that it may only apply to some but not all transcription factor families. We quantified reporter expression separately in the two notochord lineages with the expectation that Ets mutations and RBPJ mutations would have distinct effects given that primary notochord is induced by Ets-mediated FGF signaling whereas secondary notochord is induced by RBPJ/Su(H)-mediated Notch/Delta signaling. We found, however, that ETS mutations affected primary and secondary notochord expression relatively equally and that RBPJ mutations were only moderately more severe in their effect on secondary versus primary notochord. Our results point to the promise of quantitative reporter assays for understanding cis-regulatory logic but also highlight the challenge of arbitrary statistical thresholds for predicting potentially important sites.

Distal regulation, silencers, and a shared combinatorial syntax are hallmarks of animal embryogenesis

The chromatin environment plays a central role in regulating developmental gene expression in metazoans. Yet, the ancestral regulatory landscape of metazoan embryogenesis is unknown. Here, we generate chromatin accessibility profiles for six embryonic, plus larval and adult stages in the sponge Amphimedon queenslandica. These profiles are reproducible within stages, reflect histone modifications, and identify transcription factor (TF) binding sequence motifs predictive of cis-regulatory elements operating during embryogenesis in other metazoans, but not the unicellular relative Capsaspora. Motif analysis of chromatin accessibility profiles across Amphimedon embryogenesis identifies three major developmental periods. As in bilaterian embryogenesis, early development in Amphimedon involves activating and repressive chromatin in regions both proximal and distal to transcription start sites. Transcriptionally repressive elements (“silencers”) are prominent during late embryogenesis. They coincide with an increase in cis-regulatory regions harboring metazoan TF binding motifs, as well as an increase in the expression of metazoan-specific genes. Changes in chromatin state and gene expression in Amphimedon suggest the conservation of distal enhancers, dynamically silenced chromatin, and TF-DNA binding specificity in animal embryogenesis.

Transcriptional Regulation and Implications for Controlling Hox Gene Expression

Hox genes play key roles in axial patterning and regulating the regional identity of cells and tissues in a wide variety of animals from invertebrates to vertebrates. Nested domains of Hox expression generate a combinatorial code that provides a molecular framework for specifying the properties of tissues along the A–P axis. Hence, it is important to understand the regulatory mechanisms that coordinately control the precise patterns of the transcription of clustered Hox genes required for their roles in development. New insights are emerging about the dynamics and molecular mechanisms governing transcriptional regulation, and there is interest in understanding how these may play a role in contributing to the regulation of the expression of the clustered Hox genes. In this review, we summarize some of the recent findings, ideas and emerging mechanisms underlying the regulation of transcription in general and consider how they may be relevant to understanding the transcriptional regulation of Hox genes.

Metabolic gene regulation by Drosophila GATA transcription factor Grain

Nutrient-dependent gene regulation critically contributes to homeostatic control of animal physiology in changing nutrient landscape. In Drosophila, dietary sugars activate transcription factors (TFs), such as Mondo-Mlx, Sugarbabe and Cabut, which control metabolic gene expression to mediate physiological adaptation to high sugar diet. TFs that correspondingly control sugar responsive metabolic genes under conditions of low dietary sugar remain, however, poorly understood. Here we identify a role for Drosophila GATA TF Grain in metabolic gene regulation under both low and high sugar conditions. De novo motif prediction uncovered a significant over-representation of GATA-like motifs on the promoters of sugar-activated genes in Drosophila larvae, which are regulated by Grain, the fly ortholog of GATA1/2/3 subfamily. grain expression is activated by sugar in Mondo-Mlx-dependent manner and it contributes to sugar-responsive gene expression in the fat body. On the other hand, grain displays strong constitutive expression in the anterior midgut, where it drives lipogenic gene expression also under low sugar conditions. Consistently with these differential tissue-specific roles, Grain deficient larvae display delayed development on high sugar diet, while showing deregulated central carbon and lipid metabolism primarily on low sugar diet. Collectively, our study provides evidence for the role of a metazoan GATA transcription factor in nutrient-responsive metabolic gene regulation in vivo.

Systematic identification of non-canonical transcription factor motifs

Sequence-specific transcription factors (TFs) recognize motifs of related nucleotide sequences at their DNA binding sites. Upon binding at these sites, TFs regulate critical molecular processes such as gene expression. It is widely assumed that a TF recognizes a single “canonical” motif, although recent studies have identified additional “non-canonical” motifs for some TFs. A comprehensive approach to identify non-canonical DNA binding motifs and the functional importance of those motifs’ matches in the human genome is necessary for fully understanding the mechanisms of TF-regulated molecular processes in human cells. To address this need, we developed a statistical pipeline for in vitro HT-SELEX data that identifies and characterizes the distributions of non-canonical TF motifs in a stringent manner. Analyzing ~170 human TFs’ HT-SELEX data, we found non-canonical motifs for 19 TFs (11%). These non-canonical motifs occur independently of the TFs’ canonical motifs. Non-canonical motif occurrences in the human genome show similar evolutionary conservation to canonical motif occurrences, explain TF binding in locations without canonical motifs, and occur within gene promoters and epigenetically marked regulatory sequences in human cell lines and tissues. Our approach and collection of non-canonical motifs expand current understanding of functionally relevant DNA binding sites for human TFs.

Competition for DNA binding between paralogous transcription factors determines their genomic occupancy and regulatory functions

Most eukaryotic transcription factors (TFs) are part of large protein families, with members of the same family (i.e., paralogous TFs) recognizing similar DNA-binding motifs but performing different regulatory functions. Many TF paralogs are coexpressed in the cell and thus can compete for target sites across the genome. However, this competition is rarely taken into account when studying the in vivo binding patterns of eukaryotic TFs. Here, we show that direct competition for DNA binding between TF paralogs is a major determinant of their genomic binding patterns. Using yeast proteins Cbf1 and Pho4 as our model system, we designed a high-throughput quantitative assay to capture the genomic binding profiles of competing TFs in a cell-free system. Our data show that Cbf1 and Pho4 greatly influence each other's occupancy by competing for their common putative genomic binding sites. The competition is different at different genomic sites, as dictated by the TFs' expression levels and their divergence in DNA-binding specificity and affinity. Analyses of ChIP-seq data show that the biophysical rules that dictate the competitive TF binding patterns in vitro are also followed in vivo, in the complex cellular environment. Furthermore, the Cbf1-Pho4 competition for genomic sites, as characterized in vitro using our new assay, plays a critical role in the specific activation of their target genes in the cell. Overall, our study highlights the importance of direct TF-TF competition for genomic binding and gene regulation by TF paralogs, and proposes an approach for studying this competition in a quantitative and high-throughput manner.

In silico analysis of the transcriptional regulatory logic of neuronal identity specification throughout the C. elegans nervous system

The generation of the enormous diversity of neuronal cell types in a differentiating nervous system entails the activation of neuron type-specific gene batteries. To examine the regulatory logic that controls the expression of neuron type-specific gene batteries, we interrogate single cell expression profiles of all 118 neuron classes of the Caenorhabditis elegans nervous system for the presence of DNA binding motifs of 136 neuronally expressed C. elegans transcription factors. Using a phylogenetic footprinting pipeline, we identify cis-regulatory motif enrichments among neuron class-specific gene batteries and we identify cognate transcription factors for 117 of the 118 neuron classes. In addition to predicting novel regulators of neuronal identities, our nervous system-wide analysis at single cell resolution supports the hypothesis that many transcription factors directly co-regulate the cohort of effector genes that define a neuron type, thereby corroborating the concept of so-called terminal selectors of neuronal identity. Our analysis provides a blueprint for how individual components of an entire nervous system are genetically specified.

Bayesian Markov models improve the prediction of binding motifs beyond first order

Transcription factors (TFs) regulate gene expression by binding to specific DNA motifs. Accurate models for predicting binding affinities are crucial for quantitatively understanding of transcriptional regulation. Motifs are commonly described by position weight matrices, which assume that each position contributes independently to the binding energy. Models that can learn dependencies between positions, for instance, induced by DNA structure preferences, have yielded markedly improved predictions for most TFs on in vivo data. However, they are more prone to overfit the data and to learn patterns merely correlated with rather than directly involved in TF binding. We present an improved, faster version of our Bayesian Markov model software, BaMMmotif2. We tested it with state-of-the-art motif discovery tools on a large collection of ChIP-seq and HT-SELEX datasets. BaMMmotif2 models of fifth-order achieved a median false-discovery-rate-averaged recall 13.6% and 12.2% higher than the next best tool on 427 ChIP-seq datasets and 164 HT-SELEX datasets, respectively, while being 8 to 1000 times faster. BaMMmotif2 models showed no signs of overtraining in cross-cell line and cross-platform tests, with similar improvements on the next-best tool. These results demonstrate that dependencies beyond first order clearly improve binding models for most TFs.

Evolutionary Cell Type Mapping with Single-Cell Genomics

A fundamental characteristic of animal multicellularity is the spatial coexistence of functionally specialized cell types that are all encoded by a single genome sequence. Cell type transcriptional programs are deployed and maintained by regulatory mechanisms that control the asymmetric, differential access to genomic information in each cell. This genome regulation ultimately results in specific cellular phenotypes. However, the emergence, diversity, and evolutionary dynamics of animal cell types remain almost completely unexplored beyond a few species. Single-cell genomics is emerging as a powerful tool to build comprehensive catalogs of cell types and their associated gene regulatory programs in non-traditional model species. We review the current state of sampling efforts across the animal tree of life and challenges ahead for the comparative study of cell type programs. We also discuss how the phylogenetic integration of cell atlases can lead to the development of models of cell type evolution and a phylogenetic taxonomy of cells.

CrebA increases secretory capacity through direct transcriptional regulation of the secretory machinery, a subset of secretory cargo, and other key regulators

Specialization of many cells, including the acinar cells of the salivary glands and pancreas, milk-producing cells of mammary glands, mucus-secreting goblet cells, antibody-producing plasma cells, and cells that generate the dense extracellular matrices of bone and cartilage, requires scaling up both secretory machinery and cell-type specific secretory cargo. Using tissue-specific genome-scale analyses, we determine how increases in secretory capacity are coordinated with increases in secretory load in the Drosophila salivary gland (SG), an ideal model for gaining mechanistic insight into the functional specialization of secretory organs. Our findings show that CrebA, a bZIP transcription factor, directly binds genes encoding the core secretory machinery, including protein components of the signal recognition particle and receptor, ER cargo translocators, Cop I and Cop II vesicles, as well as the structural proteins and enzymes of these organelles. CrebA directly binds a subset of SG cargo genes and CrebA binds and boosts expression of Sage, a SG-specific transcription factor essential for cargo expression. To further enhance secretory output, CrebA binds and activates Xbp1 and Tudor-SN. Thus, CrebA directly upregulates the machinery of secretion and additional factors to increase overall secretory capacity in professional secretory cells; concomitant increases in cargo are achieved both directly and indirectly.

An atlas of the binding specificities of transcription factors in Pseudomonas aeruginosa directs prediction of novel regulators in virulence

A high-throughput systematic evolution of ligands by exponential enrichment assay was applied to 371 putative TFs in Pseudomonas aeruginosa, which resulted in the robust enrichment of 199 unique sequence motifs describing the binding specificities of 182 TFs. By scanning the genome, we predicted in total 33,709 significant interactions between TFs and their target loci, which were more than 11-fold enriched in the intergenic regions but depleted in the gene body regions. To further explore and delineate the physiological and pathogenic roles of TFs in P. aeruginosa, we constructed regulatory networks for nine major virulence-associated pathways and found that 51 TFs were potentially significantly associated with these virulence pathways, 32 of which had not been characterized before, and some were even involved in multiple pathways. These results will significantly facilitate future studies on transcriptional regulation in P. aeruginosa and other relevant pathogens, and accelerate to discover effective treatment and prevention strategies for the associated infectious diseases.

DNA Methylation Regulates Transcription Factor-Specific Neurodevelopmental but Not Sexually Dimorphic Gene Expression Dynamics in Zebra Finch Telencephalon

Song learning in zebra finches (Taeniopygia guttata) is a prototypical example of a complex learned behavior, yet knowledge of the underlying molecular processes is limited. Therefore, we characterized transcriptomic (RNA-sequencing) and epigenomic (RRBS, reduced representation bisulfite sequencing; immunofluorescence) dynamics in matched zebra finch telencephalon samples of both sexes from 1 day post hatching (1 dph) to adulthood, spanning the critical period for song learning (20 and 65 dph). We identified extensive transcriptional neurodevelopmental changes during postnatal telencephalon development. DNA methylation was very low, yet increased over time, particularly in song control nuclei. Only a small fraction of the massive differential expression in the developing zebra finch telencephalon could be explained by differential CpG and CpH DNA methylation. However, a strong association between DNA methylation and age-dependent gene expression was found for various transcription factors (i.e., OTX2, AR, and FOS) involved in neurodevelopment. Incomplete dosage compensation, independent of DNA methylation, was found to be largely responsible for sexually dimorphic gene expression, with dosage compensation increasing throughout life. In conclusion, our results indicate that DNA methylation regulates neurodevelopmental gene expression dynamics through steering transcription factor activity, but does not explain sexually dimorphic gene expression patterns in zebra finch telencephalon.

Systematic analysis of binding of transcription factors to noncoding variants

Many sequence variants have been linked to complex human traits and diseases1, but deciphering their biological functions remains challenging, as most of them reside in noncoding DNA. Here we have systematically assessed the binding of 270 human transcription factors to 95,886 noncoding variants in the human genome using an ultra-high-throughput multiplex protein-DNA binding assay, termed single-nucleotide polymorphism evaluation by systematic evolution of ligands by exponential enrichment (SNP-SELEX). The resulting 828 million measurements of transcription factor-DNA interactions enable estimation of the relative affinity of these transcription factors to each variant in vitro and evaluation of the current methods to predict the effects of noncoding variants on transcription factor binding. We show that the position weight matrices of most transcription factors lack sufficient predictive power, whereas the support vector machine combined with the gapped k-mer representation show much improved performance, when assessed on results from independent SNP-SELEX experiments involving a new set of 61,020 sequence variants. We report highly predictive models for 94 human transcription factors and demonstrate their utility in genome-wide association studies and understanding of the molecular pathways involved in diverse human traits and diseases.

DNA-binding properties of the MADS-domain transcription factor SEPALLATA3 and mutant variants characterized by SELEX-seq

We studied the DNA-binding profile of the MADS-domain transcription factor SEPALLATA3 and mutant variants by SELEX-seq. DNA-binding characteristics of SEPALLATA3 mutant proteins lead us to propose a novel DNA-binding mode. MIKC-type MADS-domain proteins, which function as essential transcription factors in plant development, bind as dimers to a 10-base-pair AT-rich motif termed CArG-box. However, this consensus motif cannot fully explain how the abundant family members in flowering plants can bind different target genes in specific ways. The aim of this study was to better understand the DNA-binding specificity of MADS-domain transcription factors. Also, we wanted to understand the role of a highly conserved arginine residue for binding specificity of the MADS-domain transcription factor family. Here, we studied the DNA-binding profile of the floral homeotic MADS-domain protein SEPALLATA3 by performing SELEX followed by high-throughput sequencing (SELEX-seq). We found a diverse set of bound sequences and could estimate the in vitro binding affinities of SEPALLATA3 to a huge number of different sequences. We found evidence for the preference of AT-rich motifs as flanking sequences. Whereas different CArG-boxes can act as SEPALLATA3 binding sites, our findings suggest that the preferred flanking motifs are almost always the same and thus mostly independent of the identity of the central CArG-box motif. Analysis of SEPALLATA3 proteins with a single amino acid substitution at position 3 of the DNA-binding MADS-domain further revealed that the conserved arginine residue, which has been shown to be involved in a shape readout mechanism, is especially important for the recognition of nucleotides at positions 3 and 8 of the CArG-box motif. This leads us to propose a novel DNA-binding mode for SEPALLATA3, which is different from that of other MADS-domain proteins known.

SALL4 controls cell fate in response to DNA base composition

Mammalian genomes contain long domains with distinct average compositions of A/T versus G/C base pairs. In a screen for proteins that might interpret base composition by binding to AT-rich motifs, we identified the stem cell factor SALL4, which contains multiple zinc fingers. Mutation of the domain responsible for AT binding drastically reduced SALL4 genome occupancy and prematurely upregulated genes in proportion to their AT content. Inactivation of this single AT-binding zinc-finger cluster mimicked defects seen in Sall4 null cells, including precocious differentiation of embryonic stem cells (ESCs) and embryonic lethality in mice. In contrast, deletion of two other zinc-finger clusters was phenotypically neutral. Our data indicate that loss of pluripotency is triggered by downregulation of SALL4, leading to de-repression of a set of AT-rich genes that promotes neuronal differentiation. We conclude that base composition is not merely a passive byproduct of genome evolution and constitutes a signal that aids control of cell fate.