The Human Transcription Factors.

Transcription factors (TFs) recognize specific DNA sequences to control chromatin and transcription, forming a complex system that guides expression of the genome. Despite keen interest in understanding how TFs control gene expression, it remains challenging to determine how the precise genomic binding sites of TFs are specified and how TF binding ultimately relates to regulation of transcription. This review considers how TFs are identified and functionally characterized, principally through the lens of a catalog of over 1,600 likely human TFs and binding motifs for two-thirds of them. Major classes of human TFs differ markedly in their evolutionary trajectories and expression patterns, underscoring distinct functions. TFs likewise underlie many different aspects of human physiology, disease, and variation, highlighting the importance of continued effort to understand TF-mediated gene regulation.

Uncovering a mammalian neural-specific poly(A) binding protein with unique properties.

The mRNA 3' poly(A) tail plays a critical role in regulating both mRNA translation and turnover. It is bound by the cytoplasmic poly(A) binding protein (PABPC), an evolutionarily conserved protein that can interact with translation factors and mRNA decay machineries to regulate gene expression. Mammalian PABPC1, the prototypical PABPC, is expressed in most tissues and interacts with eukaryotic translation initiation factor 4G (eIF4G) to stimulate translation in specific contexts. In this study, we uncovered a new mammalian PABPC, which we named neural PABP (neuPABP), as it is predominantly expressed in the brain. neuPABP maintains a unique architecture as compared with other PABPCs, containing only two RNA recognition motifs (RRMs) and maintaining a unique N-terminal domain of unknown function. neuPABP expression is activated in neurons as they mature during synaptogenesis, where neuPABP localizes to the soma and postsynaptic densities. neuPABP interacts with the noncoding RNA BC1, as well as mRNAs coding for ribosomal and mitochondrial proteins. However, in contrast to PABPC1, neuPABP does not associate with actively translating mRNAs in the brain. In keeping with this, we show that neuPABP has evolved such that it does not bind eIF4G and as a result fails to support protein synthesis in vitro. Taken together, these results indicate that mammals have expanded their PABPC repertoire in the brain and propose that neuPABP may support the translational repression of select mRNAs.

The complement system drives local inflammatory tissue priming by metabolic reprogramming of synovial fibroblasts.

Arthritis typically involves recurrence and progressive worsening at specific predilection sites, but the checkpoints between remission and persistence remain unknown. Here, we defined the molecular and cellular mechanisms of this inflammation-mediated tissue priming. Re-exposure to inflammatory stimuli caused aggravated arthritis in rodent models. Tissue priming developed locally and independently of adaptive immunity. Repeatedly stimulated primed synovial fibroblasts (SFs) exhibited enhanced metabolic activity inducing functional changes with intensified migration, invasiveness and osteoclastogenesis. Meanwhile, human SF from patients with established arthritis displayed a similar primed phenotype. Transcriptomic and epigenomic analyses as well as genetic and pharmacological targeting demonstrated that inflammatory tissue priming relies on intracellular complement C3- and C3a receptor-activation and downstream mammalian target of rapamycin- and hypoxia-inducible factor 1α-mediated metabolic SF invigoration that prevents activation-induced senescence, enhances NLRP3 inflammasome activity, and in consequence sensitizes tissue for inflammation. Our study suggests possibilities for therapeutic intervention abrogating tissue priming without immunosuppression.

Regulation of alternative polyadenylation by the C2H2-zinc-finger protein Sp1.

Alternative polyadenylation (APA) enhances gene regulatory potential by increasing the diversity of mRNA transcripts. 3' UTR shortening through APA correlates with enhanced cellular proliferation and is a widespread phenomenon in tumor cells. Here, we show that the ubiquitously expressed transcription factor Sp1 binds RNA in vivo and is a common repressor of distal poly(A) site usage. RNA sequencing identified 2,344 genes (36% of the total mapped mRNA transcripts) with lengthened 3' UTRs upon Sp1 depletion. Sp1 preferentially binds the 3' UTRs of such lengthened transcripts and inhibits cleavage at distal sites by interacting with the subunits of the core cleavage and polyadenylation (CPA) machinery. The 3' UTR lengths of Sp1 target genes in breast cancer patient RNA-seq data correlate with Sp1 expression levels, implicating Sp1-mediated APA regulation in modulating tumorigenic properties. Taken together, our findings provide insights into the mechanism for dynamic APA regulation by unraveling a previously unknown function of the DNA-binding transcription factor Sp1.

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.

Determination and inference of eukaryotic transcription factor sequence specificity.

Transcription factor (TF) DNA sequence preferences direct their regulatory activity, but are currently known for only ∼1% of eukaryotic TFs. Broadly sampling DNA-binding domain (DBD) types from multiple eukaryotic clades, we determined DNA sequence preferences for >1,000 TFs encompassing 54 different DBD classes from 131 diverse eukaryotes. We find that closely related DBDs almost always have very similar DNA sequence preferences, enabling inference of motifs for ∼34% of the ∼170,000 known or predicted eukaryotic TFs. Sequences matching both measured and inferred motifs are enriched in chromatin immunoprecipitation sequencing (ChIP-seq) peaks and upstream of transcription start sites in diverse eukaryotic lineages. SNPs defining expression quantitative trait loci in Arabidopsis promoters are also enriched for predicted TF binding sites. Importantly, our motif "library" can be used to identify specific TFs whose binding may be altered by human disease risk alleles. These data present a powerful resource for mapping transcriptional networks across eukaryotes.

DNA-binding specificities of human transcription factors.

Although the proteins that read the gene regulatory code, transcription factors (TFs), have been largely identified, it is not well known which sequences TFs can recognize. We have analyzed the sequence-specific binding of human TFs using high-throughput SELEX and ChIP sequencing. A total of 830 binding profiles were obtained, describing 239 distinctly different binding specificities. The models represent the majority of human TFs, approximately doubling the coverage compared to existing systematic studies. Our results reveal additional specificity determinants for a large number of factors for which a partial specificity was known, including a commonly observed A- or T-rich stretch that flanks the core motifs. Global analysis of the data revealed that homodimer orientation and spacing preferences, and base-stacking interactions, have a larger role in TF-DNA binding than previously appreciated. We further describe a binding model incorporating these features that is required to understand binding of TFs to DNA.

The human genome contains over a million autonomous exons.

Mammalian mRNA and lncRNA exons are often small compared to introns. The exon definition model predicts that exons splice autonomously, dependent on proximal exon sequence features, explaining their delineation within large introns. This model has not been examined on a genome-wide scale, however, leaving open the question of how often mRNA and lncRNA exons are autonomous. It is also unknown how frequently such exons can arise by chance. Here, we directly assayed large fragments (500-1000 bp) of the human genome by exon trapping, which detects exons spliced into a heterologous transgene, here designed with a large intron context. We define the trapped exons as "autonomous." We obtained ∼1.25 million trapped exons, including most known mRNA and well-annotated lncRNA internal exons, demonstrating that human exons are predominantly autonomous. mRNA exons are trapped with the highest efficiency. Nearly a million of the trapped exons are unannotated, most located in intergenic regions and antisense to mRNA, with depletion from the forward strand of introns. These exons are not conserved, suggesting they are nonfunctional and arose from random mutations. They are nonetheless highly enriched with known splicing promoting sequence features that delineate known exons. Novel autonomous exons are more numerous than annotated lncRNA exons, and computational models also indicate they will occur with similar frequency in any randomly generated sequence. These results show that most human coding exons splice autonomously, and provide an explanation for the existence of many unconserved lncRNAs, as well as a new annotation and inclusion levels of spliceable loci in the human genome.

An alternative splicing switch regulates embryonic stem cell pluripotency and reprogramming.

Alternative splicing (AS) is a key process underlying the expansion of proteomic diversity and the regulation of gene expression. Here, we identify an evolutionarily conserved embryonic stem cell (ESC)-specific AS event that changes the DNA-binding preference of the forkhead family transcription factor FOXP1. We show that the ESC-specific isoform of FOXP1 stimulates the expression of transcription factor genes required for pluripotency, including OCT4, NANOG, NR5A2, and GDF3, while concomitantly repressing genes required for ESC differentiation. This isoform also promotes the maintenance of ESC pluripotency and contributes to efficient reprogramming of somatic cells into induced pluripotent stem cells. These results reveal a pivotal role for an AS event in the regulation of pluripotency through the control of critical ESC-specific transcriptional programs.

The X-linked splicing regulator MBNL3 has been co-opted to restrict placental growth in eutherians.

Understanding the regulatory interactions that control gene expression during the development of novel tissues is a key goal of evolutionary developmental biology. Here, we show that Mbnl3 has undergone a striking process of evolutionary specialization in eutherian mammals resulting in the emergence of a novel placental function for the gene. Mbnl3 belongs to a family of RNA-binding proteins whose members regulate multiple aspects of RNA metabolism. We find that, in eutherians, while both Mbnl3 and its paralog Mbnl2 are strongly expressed in placenta, Mbnl3 expression has been lost from nonplacental tissues in association with the evolution of a novel promoter. Moreover, Mbnl3 has undergone accelerated protein sequence evolution leading to changes in its RNA-binding specificities and cellular localization. While Mbnl2 and Mbnl3 share partially redundant roles in regulating alternative splicing, polyadenylation site usage and, in turn, placenta maturation, Mbnl3 has also acquired novel biological functions. Specifically, Mbnl3 knockout (M3KO) alone results in increased placental growth associated with higher Myc expression. Furthermore, Mbnl3 loss increases fetal resource allocation during limiting conditions, suggesting that location of Mbnl3 on the X chromosome has led to its role in limiting placental growth, favoring the maternal side of the parental genetic conflict.

Multilayered Control of Alternative Splicing Regulatory Networks by Transcription Factors.

Networks of coordinated alternative splicing (AS) events play critical roles in development and disease. However, a comprehensive knowledge of the factors that regulate these networks is lacking. We describe a high-throughput system for systematically linking trans-acting factors to endogenous RNA regulatory events. Using this system, we identify hundreds of factors associated with diverse regulatory layers that positively or negatively control AS events linked to cell fate. Remarkably, more than one-third of the regulators are transcription factors. Further analyses of the zinc finger protein Zfp871 and BTB/POZ domain transcription factor Nacc1, which regulate neural and stem cell AS programs, respectively, reveal roles in controlling the expression of specific splicing regulators. Surprisingly, these proteins also appear to regulate target AS programs via binding RNA. Our results thus uncover a large "missing cache" of splicing regulators among annotated transcription factors, some of which dually regulate AS through direct and indirect mechanisms.

Structural insights into DNA recognition by the BEN domain of the transcription factor BANP.

The BEN domain-containing transcription factors regulate transcription by recruiting chromatin-modifying factors to specific chromatin regions via their DNA-binding BEN domains. The BEN domain of BANP has been shown to bind to a CGCG DNA sequence or an AAA-containing matrix attachment regions DNA sequence. Consistent with these in vivo observations, we identified an optimal DNA-binding sequence of AAATCTCG by protein binding microarray, which was also confirmed by our isothermal titration calorimetry and mutagenesis results. We then determined crystal structures of the BANP BEN domain in apo form and in complex with a CGCG-containing DNA, respectively, which revealed that the BANP BEN domain mainly used the electrostatic interactions to bind DNA with some base-specific interactions with the TC motifs. Our isothermal titration calorimetry results also showed that BANP bound to unmethylated and methylated DNAs with comparable binding affinities. Our complex structure of BANP-mCGCG revealed that the BANP BEN domain bound to the unmethylated and methylated DNAs in a similar mode and cytosine methylation did not get involved in binding, which is also consistent with our observations from the complex structures of the BEND6 BEN domain with the CGCG or CGmCG DNAs. Taken together, our results further elucidate the elements important for DNA recognition and transcriptional regulation by the BANP BEN domain-containing transcription factor.

From complement to complosome in non-alcoholic fatty liver disease: When location matters.

Non-alcoholic fatty liver disease (NAFLD) is a growing public health threat and becoming the leading cause of liver transplantation. Nevertheless, no approved specific treatment is currently available for NAFLD. The pathogenesis of NAFLD is multifaceted and not yet fully understood. Accumulating evidence suggests a significant role of the complement system in the development and progression of NAFLD. Here, we provide an overview of the complement system, incorporating the novel concept of complosome, and summarise the up-to-date evidence elucidating the association between complement dysregulation and the pathogenesis of NAFLD. In this process, the extracellular complement system is activated through various pathways, thereby directly contributing to, or working together with other immune cells in the disease development and progression. We also introduce the complosome and assess the evidence that implicates its potential influence in NAFLD through its direct impact on hepatocytes or non-parenchymal liver cells. Additionally, we expound upon how complement system and the complosome may exert their effects in relation with hepatic zonation in NAFLD. Furthermore, we discuss the potential therapeutic implications of targeting the complement system, extracellularly and intracellularly, for NAFLD treatment. Finally, we present future perspectives towards a better understanding of the complement system's contribution to NAFLD.

PRIESSTESS: interpretable, high-performing models of the sequence and structure preferences of RNA-binding proteins.

Modelling both primary sequence and secondary structure preferences for RNA binding proteins (RBPs) remains an ongoing challenge. Current models use varied RNA structure representations and can be difficult to interpret and evaluate. To address these issues, we present a universal RNA motif-finding/scanning strategy, termed PRIESSTESS (Predictive RBP-RNA InterpretablE Sequence-Structure moTif regrESSion), that can be applied to diverse RNA binding datasets. PRIESSTESS identifies dozens of enriched RNA sequence and/or structure motifs that are subsequently reduced to a set of core motifs by logistic regression with LASSO regularization. Importantly, these core motifs are easily visualized and interpreted, and provide a measure of RBP secondary structure specificity. We used PRIESSTESS to interrogate new HTR-SELEX data for 23 RBPs with diverse RNA binding modes and captured known primary sequence and secondary structure preferences for each. Moreover, when applying PRIESSTESS to 144 RBPs across 202 RNA binding datasets, 75% showed an RNA secondary structure preference but only 10% had a preference besides unpaired bases, suggesting that most RBPs simply recognize the accessibility of primary sequences.

Binding specificities of human RNA-binding proteins toward structured and linear RNA sequences.

RNA-binding proteins (RBPs) regulate RNA metabolism at multiple levels by affecting splicing of nascent transcripts, RNA folding, base modification, transport, localization, translation, and stability. Despite their central role in RNA function, the RNA-binding specificities of most RBPs remain unknown or incompletely defined. To address this, we have assembled a genome-scale collection of RBPs and their RNA-binding domains (RBDs) and assessed their specificities using high-throughput RNA-SELEX (HTR-SELEX). Approximately 70% of RBPs for which we obtained a motif bound to short linear sequences, whereas ∼30% preferred structured motifs folding into stem-loops. We also found that many RBPs can bind to multiple distinctly different motifs. Analysis of the matches of the motifs in human genomic sequences suggested novel roles for many RBPs. We found that three cytoplasmic proteins-ZC3H12A, ZC3H12B, and ZC3H12C-bound to motifs resembling the splice donor sequence, suggesting that these proteins are involved in degradation of cytoplasmic viral and/or unspliced transcripts. Structural analysis revealed that the RNA motif was not bound by the conventional C3H1 RNA-binding domain of ZC3H12B. Instead, the RNA motif was bound by the ZC3H12B's PilT N terminus (PIN) RNase domain, revealing a potential mechanism by which unconventional RBDs containing active sites or molecule-binding pockets could interact with short, structured RNA molecules. Our collection containing 145 high-resolution binding specificity models for 86 RBPs is the largest systematic resource for the analysis of human RBPs and will greatly facilitate future analysis of the various biological roles of this important class of proteins.

Variation in homeodomain DNA binding revealed by high-resolution analysis of sequence preferences.

Most homeodomains are unique within a genome, yet many are highly conserved across vast evolutionary distances, implying strong selection on their precise DNA-binding specificities. We determined the binding preferences of the majority (168) of mouse homeodomains to all possible 8-base sequences, revealing rich and complex patterns of sequence specificity and showing that there are at least 65 distinct homeodomain DNA-binding activities. We developed a computational system that successfully predicts binding sites for homeodomain proteins as distant from mouse as Drosophila and C. elegans, and we infer full 8-mer binding profiles for the majority of known animal homeodomains. Our results provide an unprecedented level of resolution in the analysis of this simple domain structure and suggest that variation in sequence recognition may be a factor in its functional diversity and evolutionary success.

Assessment of Lab4P Probiotic Effects on Cognition in 3xTg-AD Alzheimer's Disease Model Mice and the SH-SY5Y Neuronal Cell Line.

Aging and metabolic syndrome are associated with neurodegenerative pathologies including Alzheimer's disease (AD) and there is growing interest in the prophylactic potential of probiotic bacteria in this area. In this study, we assessed the neuroprotective potential of the Lab4P probiotic consortium in both age and metabolically challenged 3xTg-AD mice and in human SH-SY5Y cell culture models of neurodegeneration. In mice, supplementation prevented disease-associated deteriorations in novel object recognition, hippocampal neurone spine density (particularly thin spines) and mRNA expression in hippocampal tissue implying an anti-inflammatory impact of the probiotic, more notably in the metabolically challenged setting. In differentiated human SH-SY5Y neurones challenged with ß-Amyloid, probiotic metabolites elicited a neuroprotective capability. Taken together, the results highlight Lab4P as a potential neuroprotective agent and provide compelling support for additional studies in animal models of other neurodegenerative conditions and human studies.

Diverse Eukaryotic CGG-Binding Proteins Produced by Independent Domestications of hAT Transposons.

The human transcription factor (TF) CGGBP1 (CGG-binding protein) is conserved only in amniotes and is believed to derive from the zf-BED and Hermes transposase DNA-binding domains (DBDs) of a hAT DNA transposon. Here, we show that sequence-specific DNA-binding proteins with this bipartite domain structure have resulted from dozens of independent hAT domestications in different eukaryotic lineages. CGGBPs display a wide range of sequence specificity, usually including preferences for CGG or CGC trinucleotides, whereas some bind AT-rich motifs. The CGGBPs are almost entirely nonsyntenic, and their protein sequences, DNA-binding motifs, and patterns of presence or absence in genomes are uncharacteristic of ancestry via speciation. At least eight CGGBPs in the coelacanth Latimeria chalumnae bind distinct motifs, and the expression of the corresponding genes varies considerably across tissues, suggesting tissue-restricted function.

A physical and genetic map of Cannabis sativa identifies extensive rearrangements at the THC/CBD acid synthase loci.

Cannabis sativa is widely cultivated for medicinal, food, industrial, and recreational use, but much remains unknown regarding its genetics, including the molecular determinants of cannabinoid content. Here, we describe a combined physical and genetic map derived from a cross between the drug-type strain Purple Kush and the hemp variety "Finola." The map reveals that cannabinoid biosynthesis genes are generally unlinked but that aromatic prenyltransferase (AP), which produces the substrate for THCA and CBDA synthases (THCAS and CBDAS), is tightly linked to a known marker for total cannabinoid content. We further identify the gene encoding CBCA synthase (CBCAS) and characterize its catalytic activity, providing insight into how cannabinoid diversity arises in cannabis. THCAS and CBDAS (which determine the drug vs. hemp chemotype) are contained within large (>250 kb) retrotransposon-rich regions that are highly nonhomologous between drug- and hemp-type alleles and are furthermore embedded within ∼40 Mb of minimally recombining repetitive DNA. The chromosome structures are similar to those in grains such as wheat, with recombination focused in gene-rich, repeat-depleted regions near chromosome ends. The physical and genetic map should facilitate further dissection of genetic and molecular mechanisms in this commercially and medically important plant.