Serum response factor binding sites differ in three human cell types

IVEA: an integrative variational Bayesian inference method for predicting enhancer-gene regulatory interactions

Motivation

Enhancers play critical roles in cell-type-specific transcriptional control. Despite the identification of thousands of candidate enhancers, unravelling their regulatory relationships with their target genes remains challenging. Therefore, computational approaches are needed to accurately infer enhancer-gene regulatory relationships.

Results

In this study, we propose a new method, IVEA, that predicts enhancer-gene regulatory interactions by estimating promoter and enhancer activities. Its statistical model is based on the gene regulatory mechanism of transcriptional bursting, which is characterized by burst size and frequency controlled by promoters and enhancers, respectively. Using transcriptional readouts, chromatin accessibility, and chromatin contact data as inputs, promoter and enhancer activities were estimated using variational Bayesian inference, and the contribution of each enhancer-promoter pair to target gene transcription was calculated. Our analysis demonstrates that the proposed method can achieve high prediction accuracy and provide biologically relevant enhancer-gene regulatory interactions.

Availability and implementation

The IVEA code is available on GitHub at https://github.com/yasumasak/ivea. The publicly available datasets used in this study are described in Supplementary Table S4.

Vascular injury activates the ELK1/SND1/SRF pathway to promote vascular smooth muscle cell proliferative phenotype and neointimal hyperplasia

BACKGROUND

Vascular smooth muscle cell (VSMC) proliferation is the leading cause of vascular stenosis or restenosis. Therefore, investigating the molecular mechanisms and pivotal regulators of the proliferative VSMC phenotype is imperative for precisely preventing neointimal hyperplasia in vascular disease.

METHODS

Wire-induced vascular injury and aortic culture models were used to detect the expression of staphylococcal nuclease domain-containing protein 1 (SND1). SMC-specific Snd1 knockout mice were used to assess the potential roles of SND1 after vascular injury. Primary VSMCs were cultured to evaluate SND1 function on VSMC phenotype switching, as well as to investigate the mechanism by which SND1 regulates the VSMC proliferative phenotype.

RESULTS

Phenotype-switched proliferative VSMCs exhibited higher SND1 protein expression compared to the differentiated VSMCs. This result was replicated in primary VSMCs treated with platelet-derived growth factor (PDGF). In the injury model, specific knockout of Snd1 in mouse VSMCs reduced neointimal hyperplasia. We then revealed that ETS transcription factor ELK1 (ELK1) exhibited upregulation and activation in proliferative VSMCs, and acted as a novel transcription factor to induce the gene transcriptional activation of Snd1. Subsequently, the upregulated SND1 is associated with serum response factor (SRF) by competing with myocardin (MYOCD). As a co-activator of SRF, SND1 recruited the lysine acetyltransferase 2B (KAT2B) to the promoter regions leading to the histone acetylation, consequently promoted SRF to recognize the specific CArG motif, and enhanced the proliferation- and migration-related gene transcriptional activation.

CONCLUSIONS

The present study identifies ELK1/SND1/SRF as a novel pathway in promoting the proliferative VSMC phenotype and neointimal hyperplasia in vascular injury, predisposing the vessels to pathological remodeling. This provides a potential therapeutic target for vascular stenosis.

Base-resolution UV footprinting by sequencing reveals distinctive damage signatures for DNA-binding proteins

Decades ago, it was shown that proteins binding to DNA can quantitatively alter the formation of DNA damage by UV light. This established the principle of UV footprinting for non-intrusive study of protein-DNA contacts in living cells, albeit at limited scale and precision. Here, we perform deep base-resolution quantification of the principal UV damage lesion, the cyclobutane pyrimidine dimer (CPD), at select human promoter regions using targeted CPD sequencing. Several transcription factors exhibited distinctive and repeatable damage signatures indicative of site occupancy, involving strong (up to 17-fold) position-specific elevations and reductions in CPD formation frequency relative to naked DNA. Positive damage modulation at some ETS transcription factor binding sites coincided at base level with melanoma somatic mutation hotspots. Our work provides proof of concept for the study of protein-DNA interactions at individual loci using light and sequencing, and reveals widespread and potent modulation of UV damage in regulatory regions.

Hydrogel mechanics regulate fibroblast DNA methylation and chromatin condensation

Cellular mechanotransduction plays a central role in fibroblast activation during fibrotic disease progression, leading to increased tissue stiffness and reduced organ function. While the role of epigenetics in disease mechanotransduction has begun to be appreciated, little is known about how substrate mechanics, particularly the timing of mechanical inputs, regulate epigenetic changes such as DNA methylation and chromatin reorganization during fibroblast activation. In this work, we engineered a hyaluronic acid hydrogel platform with independently tunable stiffness and viscoelasticity to model normal (storage modulus, G' ∼ 0.5 kPa, loss modulus, G'' ∼ 0.05 kPa) to increasingly fibrotic (G' ∼ 2.5 and 8 kPa, G'' ∼ 0.05 kPa) lung mechanics. Human lung fibroblasts exhibited increased spreading and nuclear localization of myocardin-related transcription factor-A (MRTF-A) with increasing substrate stiffness within 1 day, with these trends holding steady for longer cultures. However, fibroblasts displayed time-dependent changes in global DNA methylation and chromatin organization. Fibroblasts initially displayed increased DNA methylation and chromatin decondensation on stiffer hydrogels, but both of these measures decreased with longer culture times. To investigate how culture time affected the responsiveness of fibroblast nuclear remodeling to mechanical signals, we engineered hydrogels amenable to in situ secondary crosslinking, enabling a transition from a compliant substrate mimicking normal tissue to a stiffer substrate resembling fibrotic tissue. When stiffening was initiated after only 1 day of culture, fibroblasts rapidly responded and displayed increased DNA methylation and chromatin decondensation, similar to fibroblasts on static stiffer hydrogels. Conversely, when fibroblasts experienced later stiffening at day 7, they showed no changes in DNA methylation and chromatin condensation, suggesting the induction of a persistent fibroblast phenotype. These results highlight the time-dependent nuclear changes associated with fibroblast activation in response to dynamic mechanical perturbations and may provide mechanisms to target for controlling fibroblast activation.

Serum response factor-cofactor interactions and their implications in disease

Serum response factor (SRF), a member of the Mcm1, Agamous, Deficiens, and SRF (MADS) box transcription factor, is widely expressed in all cell types and plays a crucial role in the physiological function and development of diseases. SRF regulates its downstream genes by binding to their CArG DNA box by interacting with various cofactors. However, the underlying mechanisms are not fully understood, therefore attracting increasing research attention due to the importance of this topic. This review's objective is to discuss the new progress in the studies of the molecular mechanisms involved in the activation of SRF and its impacts in physiological and pathological conditions. Notably, we summarized the recent studies on the interaction of SRF with its two main types of cofactors belonging to the myocardin families of transcription factors and the members of the ternary complex factors. The knowledge of these mechanisms will create new opportunities for understanding the dynamics of many traits and disease pathogenesis especially, cardiovascular diseases and cancer that could serve as targets for pharmacological control and treatment of these diseases.

TRA2A-induced upregulation of LINC00662 regulates blood-brain barrier permeability by affecting ELK4 mRNA stability in Alzheimer's microenvironment

The blood-brain barrier (BBB) plays a pivotal role in the maintenance and regulation of the neural microenvironment. The BBB breakdown is a pathological change in early Alzheimer's disease (AD). RNA-binding proteins (RBPs) and long non-coding RNAs (lncRNAs) are involved in the regulation of BBB permeability. Our study demonstrates the role of TRA2A/LINC00662/ELK4 axis in regulating BBB permeability in AD microenvironment. In Aβ_1-42-incubated microvascular endothelial cells (ECs) of the BBB model in vitro, TRA2A and LINC00662 were enriched. TRA2A increased the stability of LINC00662 by binding with it. The knockdown of either TRA2A or LINC00662 decreased BBB permeability due to increased expression of tight junction-related proteins. ELK4 was less expressed in the BBB model in AD microenvironment in vitro. LINC00662 mediated the degradation of ELK4 mRNA by SMD pathway. Downregulation of ELK4 increased BBB permeability by increasing the tight junction-related protein expression.TRA2A/LINC00662/ELK4 axis plays a crucial role in the regulation of BBB permeability in AD microenvironment, which may provide a novel target for the therapy of AD.

EGR-mediated control of STIM expression and function

Ca²⁺ is a ubiquitous, dynamic and pluripotent second messenger with highly context-dependent roles in complex cellular processes such as differentiation, proliferation, and cell death. These Ca²⁺ signals are generated by Ca²⁺-permeable channels located on the plasma membrane (PM) and endoplasmic reticulum (ER) and shaped by PM- and ER-localized pumps and transporters. Differences in the expression of these Ca²⁺ homeostasis proteins contribute to cell and context-dependent differences in the spatiotemporal organization of Ca²⁺ signals and, ultimately, cell fate. This review focuses on the Early Growth Response (EGR) family of zinc finger transcription factors and their role in the transcriptional regulation of Stromal Interaction Molecule (STIM1), a critical regulator of Ca²⁺ entry in both excitable and non-excitable cells.

Srf destabilizes cellular identity by suppressing cell-type-specific gene expression programs

Multicellular organisms consist of multiple cell types. The identity of these cells is primarily maintained by cell-type-specific gene expression programs; however, mechanisms that suppress these programs are poorly defined. Here we show that serum response factor (Srf), a transcription factor that is activated by various extracellular stimuli, can repress cell-type-specific genes and promote cellular reprogramming to pluripotency. Manipulations that decrease β-actin monomer quantity result in the nuclear accumulation of Mkl1 and the activation of Srf, which downregulate cell-type-specific genes and alter the epigenetics of regulatory regions and chromatin organization. Mice overexpressing Srf exhibit various pathologies including an ulcerative colitis-like symptom and a metaplasia-like phenotype in the pancreas. Our results demonstrate an unexpected function of Srf via a mechanism by which extracellular stimuli actively destabilize cell identity and suggest Srf involvement in a wide range of diseases.

The transcription factor Srf is a central regulator of immediate-early and actin cytoskeletal genes. Here the authors show that Srf is activated by repression of β-actin, promoting iPSC reprogramming of neural progenitor cells and hepatoblasts by repressing cell-type specific genes.

Syndecan-2 Attenuates Radiation-induced Pulmonary Fibrosis and Inhibits Fibroblast Activation by Regulating PI3K/Akt/ROCK Pathway via CD148

Radiation-induced pulmonary fibrosis is a severe complication of patients treated with thoracic irradiation. We have previously shown that syndecan-2 reduces fibrosis by exerting alveolar epithelial cytoprotective effects. Here, we investigate whether syndecan-2 attenuates radiation-induced pulmonary fibrosis by inhibiting fibroblast activation. C57BL/6 wild-type mice and transgenic mice that overexpress human syndecan-2 in alveolar macrophages were exposed to 14 Gy whole-thoracic radiation. At 24 weeks after irradiation, lungs were collected for histological, protein, and mRNA evaluation of pulmonary fibrosis, profibrotic gene expression, and α-smooth muscle actin (α-SMA) expression. Mouse lung fibroblasts were activated with transforming growth factor (TGF)-β1 in the presence or absence of syndecan-2. Cell proliferation, migration, and gel contraction were assessed at different time points. Irradiation resulted in significantly increased mortality and pulmonary fibrosis in wild-type mice that was associated with elevated lung expression of TGF-β1 downstream target genes and cell death compared with irradiated syndecan-2 transgenic mice. In mouse lung fibroblasts, syndecan-2 inhibited α-SMA expression, cell contraction, proliferation, and migration induced by TGF-β1. Syndecan-2 attenuated phosphoinositide 3-kinase/serine/threonine kinase/Rho-associated coiled-coil kinase signaling and serum response factor binding to the α-SMA promoter. Syndecan-2 attenuates pulmonary fibrosis in mice exposed to radiation and inhibits TGF-β1-induced fibroblast-myofibroblast differentiation, migration, and proliferation by down-regulating phosphoinositide 3-kinase/serine/threonine kinase/Rho-associated coiled-coil kinase signaling and blocking serum response factor binding to the α-SMA promoter via CD148. These findings suggest that syndecan-2 has potential as an antifibrotic therapy in radiation-induced lung fibrosis.

Neuroprotective effect of Ruminococcus albus on oxidatively stressed SH-SY5Y cells and animals

Recent evidence shows that the gut microbiota has an important role in gut-brain crosstalk and is linked to neuronal disorders. The aim of this study was to investigate the effects of intestinal Ruminococcus albus with probiotic potential on neuroprotection in oxidatively stressed SH-SY5Y neuroblastoma cells and animals. To investigate these effects, conditioned medium was prepared using Caco-2 cells cultured with heat-killed R. albus (CRA-CM). Caco-2 cells cultured with heat-killed R. albus showed increased BDNF expression and BDNF protein levels increased in CRA-CM. CRA-CM up-regulated the protein expression levels of SRF, C-fos and CDK2. In addition, CRA-CM protected SH-SY5Y cells from H₂O₂-induced cell death. CRA-CM significantly decreased the Bax/Bcl-2 ratio in oxidatively stressed SH-SY5Y cells. Animal experiments showed that oral administration of heat-killed R. albus for 15 days attenuated the oxidative stress induced by sodium arsenate. Treatment with heat-killed R. albus reduced the level of ROS, and the levels of SOD and GSH increased in oxidatively stressed brains. In conclusion, the secretome prepared from Caco-2 cells cultured with heat-killed R. albus might promote neuronal proliferation through the activation of cell proliferation-related proteins, and heat-killed R. albus protects neurons from oxidative damage by reducing ROS levels and increasing SOD and GSH levels.

The organization and regulation of mRNA-protein complexes

In a eukaryotic cell, each messenger RNA (mRNA) is bound to a variety of proteins to form an mRNA-protein complex (mRNP). Together, these proteins impact nearly every step in the life cycle of an mRNA and are critical for the proper control of gene expression. In the cytoplasm, for instance, mRNPs affect mRNA translatability and stability and provide regulation of specific transcripts as well as global, transcriptome-wide control. mRNPs are complex, diverse, and dynamic, and so they have been a challenge to understand. But the advent of high-throughput sequencing technology has heralded a new era in the study of mRNPs. Here, I will discuss general principles of cytoplasmic mRNP organization and regulation. Using microRNA-mediated repression as a case study, I will focus on common themes in mRNPs and highlight the interplay between mRNP composition and posttranscriptional regulation. mRNPs are an important control point in regulating gene expression, and while the study of these fascinating complexes presents remaining challenges, recent advances provide a critical lens for deciphering gene regulation. WIREs RNA 2017, 8:e1369. doi: 10.1002/wrna.1369 For further resources related to this article, please visit the WIREs website.

Multi-phenotypic Role of Serum Response Factor in the Gastrointestinal System

Serum response factor (SRF) is a master transcription factor of the actin cytoskeleton that binds to highly conserved CArG boxes located within the majority of smooth muscle cell (SMC)-restricted promoters/enhancers. Although most studies of SRF focus on skeletal muscle, cardiac muscle, and vascular SMCs, SRF research has recently expanded into the gastrointestinal (GI) system. Genome scale analyses of GI SMC transcriptome and CArG boxes (CArGome) have identified new SRF target genes. In addition to circular and longitudinal smooth muscle layers, SRF is also expressed in GI mucosa and cancers. In the GI tract, SRF is the central regulator of genes involved in apoptosis, dedifferentiation, proliferation, and migration of cells. Since SRF is the cell phenotypic modulator, it may play an essential role in the development of myopathy, hypertrophy, ulcers, gastric and colon cancers within the GI tract. Given the multi-functional role displayed by SRF in the digestive system, SRF has received more attention emerging as a potential therapeutic target. This review summarizes the findings in SRF research pertaining to the GI tract and provides valuable insight into future directions.

MOCCS: Clarifying DNA-binding motif ambiguity using ChIP-Seq data

BACKGROUND

As a key mechanism of gene regulation, transcription factors (TFs) bind to DNA by recognizing specific short sequence patterns that are called DNA-binding motifs. A single TF can accept ambiguity within its DNA-binding motifs, which comprise both canonical (typical) and non-canonical motifs. Clarification of such DNA-binding motif ambiguity is crucial for revealing gene regulatory networks and evaluating mutations in cis-regulatory elements. Although chromatin immunoprecipitation sequencing (ChIP-seq) now provides abundant data on the genomic sequences to which a given TF binds, existing motif discovery methods are unable to directly answer whether a given TF can bind to a specific DNA-binding motif.

RESULTS

Here, we report a method for clarifying the DNA-binding motif ambiguity, MOCCS. Given ChIP-Seq data of any TF, MOCCS comprehensively analyzes and describes every k-mer to which that TF binds. Analysis of simulated datasets revealed that MOCCS is applicable to various ChIP-Seq datasets, requiring only a few minutes per dataset. Application to the ENCODE ChIP-Seq datasets proved that MOCCS directly evaluates whether a given TF binds to each DNA-binding motif, even if known position weight matrix models do not provide sufficient information on DNA-binding motif ambiguity. Furthermore, users are not required to provide numerous parameters or background genomic sequence models that are typically unavailable. MOCCS is implemented in Perl and R and is freely available via https://github.com/yuifu/moccs.

CONCLUSIONS

By complementing existing motif-discovery software, MOCCS will contribute to the basic understanding of how the genome controls diverse cellular processes via DNA-protein interactions.

Integrins synergise to induce expression of the MRTF-A-SRF target gene ISG15 for promoting cancer cell invasion

Integrin-mediated activation of small GTPases induces the polymerisation of G-actin into various actin structures and the release of the transcriptional co-activator MRTF from G-actin. Here we report that pan-integrin-null fibroblasts seeded on fibronectin and expressing β1- and/or αV-class integrin contained different G-actin pools, nuclear MRTF-A (also known as MKL1 or MAL) levels and MRTF-A-SRF activities. The nuclear MRTF-A levels and activities were highest in cells expressing both integrin classes, lower in cells expressing β1 integrins and lowest in cells expressing the αV integrins. Quantitative proteomics and transcriptomics analyses linked the differential MRTF-A activities to the expression of the ubiquitin-like modifier interferon-stimulated gene 15 (ISG15), which is known to modify focal adhesion and cytoskeletal proteins. The malignant breast cancer cell line MDA-MB-231 expressed high levels of β1 integrins, ISG15 and ISGylated proteins, which promoted invasive properties, whereas non-invasive MDA-MB-468 and MCF-7 cell lines expressed low levels of β1 integrins, ISG15 and ISGylated proteins. Our findings suggest that integrin-adhesion-induced MRTF-A-SRF activation and ISG15 expression constitute a newly discovered signalling circuit that promotes cell migration and invasion.

The short and long of noncoding sequences in the control of vascular cell phenotypes

The two principal cell types of importance for normal vessel wall physiology are smooth muscle cells and endothelial cells. Much progress has been made over the past 20 years in the discovery and function of transcription factors that coordinate proper differentiation of these cells and the maintenance of vascular homeostasis. More recently, the converging fields of bioinformatics, genomics, and next generation sequencing have accelerated discoveries in a number of classes of noncoding sequences, including transcription factor binding sites (TFBS), microRNA genes, and long noncoding RNA genes, each of which mediates vascular cell differentiation through a variety of mechanisms. Alterations in the nucleotide sequence of key TFBS or deviations in transcription of noncoding RNA genes likely have adverse effects on normal vascular cell phenotype and function. Here, the subject of noncoding sequences that influence smooth muscle cell or endothelial cell phenotype will be summarized as will future directions to further advance our understanding of the increasingly complex molecular circuitry governing normal vascular cell differentiation and how such information might be harnessed to combat vascular diseases.

Identification of highly conserved putative developmental enhancers bound by SOX3 in neural progenitors using ChIP-Seq

The transcription factor SOX3 is expressed within most neural progenitor (NP) cells of the vertebrate central nervous system (CNS) and is essential for normal brain development in mice and humans. However, despite the widespread expression of Sox3, CNS defects in null mice are relatively mild due to functional redundancy with the other SOXB1 sub-group members Sox1 and Sox2. To further understand the molecular function of SOX3, we investigated the genome-wide binding profile of endogenous SOX3 in NP cells using ChIP-seq. SOX3 binding was identified at over 8,000 sites, most of which were intronic or intergeneic and were significantly associated with neurodevelopmental genes. The majority of binding sites were moderately or highly conserved (phastCons scores >0.1 and 0.5, respectively) and included the previously characterised, SOXB1-binding Nestin NP cell enhancer. Comparison of SOX3 and published ChIP-Seq data for the co-activator P300 in embryonic brain identified hundreds of highly conserved putative enhancer elements. In addition, we identified a subset of highly conserved putative enhancers for CNS development genes common to SOXB1 members in NP cells, all of which contained the SOX consensus motif (ACAAWR). Together these data implicate SOX3 in the direct regulation of hundreds of NP genes and provide molecular insight into the overlapping roles of SOXB1 proteins in CNS development.

The core and conserved role of MAL is homeostatic regulation of actin levels

The transcription cofactor MAL is regulated by actin dynamics and, together with its DNA-binding partner, SRF, is required for invasive cell migration. Salvany et al. show in Drosophila and human cellular models that actin is the key target that must be regulated by MAL/SRF for invasive cell migration. By regulating MAL/SRF, actin feeds back on the production of actin mRNA to ensure sufficient actin supply. Actin and MAL thus form a homeostatic feedback system that provides the foundation for actin dynamics required for complex cell behavior.

The transcription cofactor MAL is regulated by free actin levels and thus by actin dynamics. MAL, together with its DNA-binding partner, SRF, is required for invasive cell migration and in experimental metastasis. Although MAL/SRF has many targets, we provide genetic evidence in both Drosophila and human cellular models that actin is the key target that must be regulated by MAL/SRF for invasive cell migration. By regulating MAL/SRF activity, actin protein feeds back on production of actin mRNA to ensure sufficient supply of actin. This constitutes a dedicated homeostatic feedback system that provides a foundation for cellular actin dynamics.

The neuronal activity-driven transcriptome

Activity-driven transcription is a key event associated with long-lasting forms of neuronal plasticity. Despite the efforts to investigate the regulatory mechanisms that control this complex process and the important advances in the knowledge of the function of many activity-induced genes in neurons, as well as the specific contribution of activity-regulated transcription factors, our understanding of how activity-driven transcription operates at the systems biology level is still very limited. This review focuses on the research of neuronal activity-driven transcription from an "omics" perspective. We will discuss the different high-throughput approaches undertaken to characterize the gene programs downstream of specific activity-regulated transcription factors, including CREB, SRF, MeCP2, Fos, Npas4, and others, and the interplay between epigenetic and transcriptional mechanisms underlying neuronal plasticity changes. Although basic questions remain unanswered and important challenges still lie ahead, the refinement of genome-wide techniques for investigating the neuronal transcriptome and epigenome promises great advances.

SRF phosphorylation by glycogen synthase kinase-3 promotes axon growth in hippocampal neurons

The growth of axons is an intricately regulated process involving intracellular signaling cascades and gene transcription. We had previously shown that the stimulus-dependent transcription factor, serum response factor (SRF), plays a critical role in regulating axon growth in the mammalian brain. However, the molecular mechanisms underlying SRF-dependent axon growth remains unknown. Here we report that SRF is phosphorylated and activated by GSK-3 to promote axon outgrowth in mouse hippocampal neurons. GSK-3 binds to and directly phosphorylates SRF on a highly conserved serine residue. This serine phosphorylation is necessary for SRF activity and for its interaction with MKL-family cofactors, MKL1 and MKL2, but not with TCF-family cofactor, ELK-1. Axonal growth deficits caused by GSK-3 inhibition could be rescued by expression of a constitutively active SRF. The SRF target gene and actin-binding protein, vinculin, is sufficient to overcome the axonal growth deficits of SRF-deficient and GSK-3-inhibited neurons. Furthermore, short hairpin RNA-mediated knockdown of vinculin also attenuated axonal growth. Thus, our findings reveal a novel phosphorylation and activation of SRF by GSK-3 that is critical for SRF-dependent axon growth in mammalian central neurons.

Global analyses of the effect of different cellular contexts on microRNA targeting

MicroRNA (miRNA) regulation clearly impacts animal development, but the extent to which development-with its resulting diversity of cellular contexts-impacts miRNA regulation is unclear. Here, we compared cohorts of genes repressed by the same miRNAs in different cell lines and tissues and found that target repertoires were largely unaffected, with secondary effects explaining most of the differential responses detected. Outliers resulting from differential direct targeting were often attributable to alternative 3' UTR isoform usage that modulated the presence of miRNA sites. More inclusive examination of alternative 3' UTR isoforms revealed that they influence ∼10% of predicted targets when comparing any two cell types. Indeed, considering alternative 3' UTR isoform usage improved prediction of targeting efficacy significantly beyond the improvements observed when considering constitutive isoform usage. Thus, although miRNA targeting is remarkably consistent in different cell types, considering the 3' UTR landscape helps predict targeting efficacy and explain differential regulation that is observed.

Comparative (computational) analysis of the DNA methylation status of trinucleotide repeat expansion diseases

Previous studies have examined DNA methylation in different trinucleotide repeat diseases. We have combined this data and used a pattern searching algorithm to identify motifs in the DNA surrounding aberrantly methylated CpGs found in the DNA of patients with one of the three trinucleotide repeat (TNR) expansion diseases: fragile X syndrome (FRAXA), myotonic dystrophy type I (DM1), or Friedreich's ataxia (FRDA). We examined sequences surrounding both the variably methylated (VM) CpGs, which are hypermethylated in patients compared with unaffected controls, and the nonvariably methylated CpGs which remain either always methylated (AM) or never methylated (NM) in both patients and controls. Using the J48 algorithm of WEKA analysis, we identified that two patterns are all that is necessary to classify our three regions CCGG∗ which is found in VM and not in AM regions and AATT∗ which distinguished between NM and VM + AM using proportional frequency. Furthermore, comparing our software with MEME software, we have demonstrated that our software identifies more patterns than MEME in these short DNA sequences. Thus, we present evidence that the DNA sequence surrounding CpG can influence its susceptibility to be de novo methylated in a disease state associated with a trinucleotide repeat.

Serum regulation of Id1 expression by a BMP pathway and BMP responsive element

Immediate early genes (IEGs) are expressed upon re-entry of quiescent cells into the cell cycle following serum stimulation. These genes are involved in growth control and differentiation and hence their expression is tightly controlled. Many IEGs are regulated through Serum Response Elements (SREs) in their promoters, which bind Serum Response Factor (SRF). However, many other IEGs do not have SREs in their promoters and their serum regulation is poorly understood. We have identified SRF-independent IEGs in SRF-depleted fibroblasts. One of these, Id1, was examined more closely. We mapped a serum responsive element in the Id1 promoter and find that it is identical to a BMP responsive element (BRE). The Id1 BRE is necessary and sufficient for the serum regulation of Id1. Inhibition of the BMP pathway by siRNA depletion of Smad 4, treatment with the BMP antagonist noggin, or the BMP receptor inhibitor dorsomorphin blocked serum induction of Id1. Further, BMP2 is sufficient to induce Id1 expression. Given reports that SRC inhibitors can block Id1 expression, we tested the SRC inhibitor, AZD0530, and found that it inhibits the serum activation of Id1. Surprisingly, this inhibition is independent of SRC or its family members. Rather, we show that AZD0530 directly inhibits the BMP type I receptors. Serum induction of the Id1 related gene Id3 also required the BMP pathway. Given these and other findings we conclude that the Id family of IEGs is regulated by BMPs in serum through similar BREs. This represents a second pathway for serum regulation of IEGs.

Integrative modeling of eQTLs and cis-regulatory elements suggests mechanisms underlying cell type specificity of eQTLs

Genetic variants in cis-regulatory elements or trans-acting regulators frequently influence the quantity and spatiotemporal distribution of gene transcription. Recent interest in expression quantitative trait locus (eQTL) mapping has paralleled the adoption of genome-wide association studies (GWAS) for the analysis of complex traits and disease in humans. Under the hypothesis that many GWAS associations tag non-coding SNPs with small effects, and that these SNPs exert phenotypic control by modifying gene expression, it has become common to interpret GWAS associations using eQTL data. To fully exploit the mechanistic interpretability of eQTL-GWAS comparisons, an improved understanding of the genetic architecture and causal mechanisms of cell type specificity of eQTLs is required. We address this need by performing an eQTL analysis in three parts: first we identified eQTLs from eleven studies on seven cell types; then we integrated eQTL data with cis-regulatory element (CRE) data from the ENCODE project; finally we built a set of classifiers to predict the cell type specificity of eQTLs. The cell type specificity of eQTLs is associated with eQTL SNP overlap with hundreds of cell type specific CRE classes, including enhancer, promoter, and repressive chromatin marks, regions of open chromatin, and many classes of DNA binding proteins. These associations provide insight into the molecular mechanisms generating the cell type specificity of eQTLs and the mode of regulation of corresponding eQTLs. Using a random forest classifier with cell specific CRE-SNP overlap as features, we demonstrate the feasibility of predicting the cell type specificity of eQTLs. We then demonstrate that CREs from a trait-associated cell type can be used to annotate GWAS associations in the absence of eQTL data for that cell type. We anticipate that such integrative, predictive modeling of cell specificity will improve our ability to understand the mechanistic basis of human complex phenotypic variation.

When interpreting genome-wide association studies showing that specific genetic variants are associated with disease risk, scientists look for a link between the genetic variant and a biological mechanism behind that disease. One functional mechanism is that the genetic variant may influence gene transcription via a co-localized genomic regulatory element, such as a transcription factor binding site within an open chromatin region. Often this type of regulation occurs in some cell types but not others. In this study, we look across eleven gene expression studies with seven cell types and consider how genetic transcription regulators, or eQTLs, replicate within and between cell types. We identify pervasive allelic heterogeneity, or transcriptional control of a single gene by multiple, independent eQTLs. We integrate extensive data on cell type specific regulatory elements from ENCODE to identify general methods of transcription regulation through enrichment of eQTLs within regulatory elements. We also build a classifier to predict eQTL replication across cell types. The results in this paper present a path to an integrative, predictive approach to improve our ability to understand the mechanistic basis of human phenotypic variation.

Androgen-responsive serum response factor target genes regulate prostate cancer cell migration

Progression of prostate cancer (CaP) relies on androgen receptor (AR) signaling, but AR-dependent events that underlie the lethal phenotype remain unknown. Recently, an indirect mechanism of androgen action in which effects of AR on CaP cells are mediated by Serum Response Factor (SRF) has been identified. This is the first mode of androgen action to be associated with aggressive CaP and disease recurrence. The manner in which androgen-responsive SRF activity controls aggressive CaP cell behavior is unknown. Here, the contribution of two representative SRF effector genes that are underexpressed, calponin 2 (CNN2), or overexpressed, sidekick-homolog 1 (SDK1), in clinical CaP specimens is studied. AR- and SRF- dependency of CNN2 and SDK1 expression was verified using synthetic and natural androgens, antiandrogens, and small interfering RNAs targeting AR or SRF, and evaluating the kinetics of androgen induction and SRF binding to endogenously and exogenously expressed regulatory gene regions in AR-positive CaP model systems that mimic the transition from androgen-stimulated to castration-recurrent disease. Small interfering RNA-mediated deregulation of CNN2 or SDK1 expression did not affect CaP cell proliferation or apoptosis but had marked effects on CaP cell morphology and actin cytoskeleton organization. Loss of CNN2 induced cellular protrusions and increased CaP cell migration, whereas silencing of SDK1 led to cell rounding and blunted CaP cell migration. Changes in cell migration did not involve epithelial-mesenchymal transition but correlated with altered β1-integrin expression. Taken together, individual androgen-responsive SRF target genes affect CaP cell behavior by modulating cell migration, which may have implications for therapeutic intervention downstream of AR and SRF.

The mouse DXZ4 homolog retains Ctcf binding and proximity to Pls3 despite substantial organizational differences compared to the primate macrosatellite

Background

The X-linked macrosatellite DXZ4 is a large homogenous tandem repeat that in females adopts an alternative chromatin organization on the primate X chromosome in response to X-chromosome inactivation. It is packaged into heterochromatin on the active X chromosome but into euchromatin and bound by the epigenetic organizer protein CTCF on the inactive X chromosome. Because its DNA sequence diverges rapidly beyond the New World monkeys, the existence of DXZ4 outside the primate lineage is unknown.

Results

Here we extend our comparative genome analysis and report the identification and characterization of the mouse homolog of the macrosatellite. Furthermore, we provide evidence of DXZ4 in a conserved location downstream of the PLS3 gene in a diverse group of mammals, and reveal that DNA sequence conservation is restricted to the CTCF binding motif, supporting a central role for this protein at this locus. However, many features that characterize primate DXZ4 differ in mouse, including the overall size of the array, the mode of transcription, the chromatin organization and conservation between adjacent repeat units of DNA sequence and length. Ctcf binds Dxz4 but is not exclusive to the inactive X chromosome, as evidenced by association in some males and equal binding to both X chromosomes in trophoblast stem cells.

Conclusions

Characterization of Dxz4 reveals substantial differences in the organization of DNA sequence, chromatin packaging, and the mode of transcription, so the potential roles performed by this sequence in mouse have probably diverged from those on the primate X chromosome.

Cell-type specificity of ChIP-predicted transcription factor binding sites

BACKGROUND

Context-dependent transcription factor (TF) binding is one reason for differences in gene expression patterns between different cellular states. Chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) identifies genome-wide TF binding sites for one particular context-the cells used in the experiment. But can such ChIP-seq data predict TF binding in other cellular contexts and is it possible to distinguish context-dependent from ubiquitous TF binding?

RESULTS

We compared ChIP-seq data on TF binding for multiple TFs in two different cell types and found that on average only a third of ChIP-seq peak regions are common to both cell types. Expectedly, common peaks occur more frequently in certain genomic contexts, such as CpG-rich promoters, whereas chromatin differences characterize cell-type specific TF binding. We also find, however, that genotype differences between the cell types can explain differences in binding. Moreover, ChIP-seq signal intensity and peak clustering are the strongest predictors of common peaks. Compared with strong peaks located in regions containing peaks for multiple transcription factors, weak and isolated peaks are less common between the cell types and are less associated with data that indicate regulatory activity.

CONCLUSIONS

Together, the results suggest that experimental noise is prevalent among weak peaks, whereas strong and clustered peaks represent high-confidence binding events that often occur in other cellular contexts. Nevertheless, 30-40% of the strongest and most clustered peaks show context-dependent regulation. We show that by combining signal intensity with additional data-ranging from context independent information such as binding site conservation and position weight matrix scores to context dependent chromatin structure-we can predict whether a ChIP-seq peak is likely to be present in other cellular contexts.

MKL1 and MKL2 play redundant and crucial roles in megakaryocyte maturation and platelet formation

Serum response factor and its transcriptional cofactor MKL1 are critical for megakaryocyte maturation and platelet formation. We show that MKL2, a homologue of MKL1, is expressed in megakaryocytes and plays a role in megakaryocyte maturation. Using a megakaryocyte-specific Mkl2 knockout (KO) mouse on the conventional Mkl1 KO background to produce double KO (DKO) megakaryocytes and platelets, a critical role for MKL2 is revealed. The decrease in megakaryocyte ploidy and platelet counts of DKO mice is more severe than in Mkl1 KO mice. Platelet dysfunction in DKO mice is revealed by prolonged bleeding times and ineffective platelet activation in vitro in response to adenosine 5'-diphosphate. Electron microscopy and immunofluorescence of DKO megakaryocytes and platelets indicate abnormal cytoskeletal and membrane organization with decreased granule complexity. Surprisingly, the DKO mice have a more extreme thrombocytopenia than mice lacking serum response factor (SRF) expression in the megakaryocyte compartment. Comparison of gene expression reveals approximately 4400 genes whose expression is differentially affected in DKO compared with megakaryocytes deficient in SRF, strongly suggesting that MKL1 and MKL2 have both SRF-dependent and SRF-independent activity in megakaryocytopoiesis.

Role of RhoA-specific guanine exchange factors in regulation of endomitosis in megakaryocytes

Polyploidization can precede the development of aneuploidy in cancer. Polyploidization in megakaryocytes (Mks), in contrast, is a highly controlled developmental process critical for efficient platelet production via unknown mechanisms. Using primary cells, we demonstrate that the guanine exchange factors GEF-H1 and ECT2, which are often overexpressed in cancer and are essential for RhoA activation during cytokinesis, must be downregulated for Mk polyploidization. The first (2N-4N) endomitotic cycle requires GEF-H1 downregulation, whereas subsequent cycles (>4N) require ECT2 downregulation. Exogenous expression of both GEF-H1 and ECT2 prevents endomitosis, resulting in proliferation of 2N Mks. Furthermore, we have shown that the mechanism by which polyploidization is prevented in Mks lacking Mkl1, which is mutated in megakaryocytic leukemia, is via elevated GEF-H1 expression; shRNA-mediated GEF-H1 knockdown alone rescues this ploidy defect. These mechanistic insights enhance our understanding of normal versus malignant megakaryocytopoiesis, as well as aberrant mitosis in aneuploid cancers.

Leiomodin 1, a new serum response factor-dependent target gene expressed preferentially in differentiated smooth muscle cells

Smooth muscle cell (SMC) differentiation is defined largely by a number of cell-restricted genes governed directly by the serum response factor (SRF)/myocardin (MYOCD) transcriptional switch. Here, we describe a new SRF/MYOCD-dependent, SMC-restricted gene known as Leiomodin 1 (Lmod1). Conventional and quantitative RT-PCRs indicate that Lmod1 mRNA expression is enriched in SMC-containing tissues of the mouse, whereas its two paralogs, Lmod2 and Lmod3, exhibit abundant expression in skeletal and cardiac muscle with very low levels in SMC-containing tissues. Western blotting and immunostaining of various adult and embryonic mouse tissues further confirm SMC-specific expression of the LMOD1 protein. Comparative genomic analysis of the human LMOD1 and LMOD2 genes with their respective mouse and rat orthologs shows high conservation between the three exons and several noncoding sequences, including the immediate 5' promoter region. Two conserved CArG boxes are present in both the LMOD1 and LMOD2 promoter regions, although LMOD1 displays much higher promoter activity and is more responsive to SRF/MYOCD stimulation. Gel shift assays demonstrate clear binding between SRF and the two CArG boxes in human LMOD1. Although the CArG boxes in LMOD1 and LMOD2 are similar, only LMOD1 displays SRF or MYOCD-dependent activation. Transgenic mouse studies reveal wild type LMOD1 promoter activity in cardiac and vascular SMC. Such activity is abolished upon mutation of both CArG boxes. Collectively, these data demonstrate that Lmod1 is a new SMC-restricted SRF/MYOCD target gene.

Genomic and biochemical insights into the specificity of ETS transcription factors

ETS proteins are a group of evolutionarily related, DNA-binding transcriptional factors. These proteins direct gene expression in diverse normal and disease states by binding to specific promoters and enhancers and facilitating assembly of other components of the transcriptional machinery. The highly conserved DNA-binding ETS domain defines the family and is responsible for specific recognition of a common sequence motif, 5'-GGA(A/T)-3'. Attaining specificity for biological regulation in such a family is thus a conundrum. We present the current knowledge of routes to functional diversity and DNA binding specificity, including divergent properties of the conserved ETS and PNT domains, the involvement of flanking structured and unstructured regions appended to these dynamic domains, posttranslational modifications, and protein partnerships with other DNA-binding proteins and coregulators. The review emphasizes recent advances from biochemical and biophysical approaches, as well as insights from genomic studies that detect ETS-factor occupancy in living cells.

Identifying functional single nucleotide polymorphisms in the human CArGome

Regulatory SNPs (rSNPs) reside primarily within the nonprotein coding genome and are thought to disturb normal patterns of gene expression by altering DNA binding of transcription factors. Nevertheless, despite the explosive rise in SNP association studies, there is little information as to the function of rSNPs in human disease. Serum response factor (SRF) is a widely expressed DNA-binding transcription factor that has variable affinity to at least 1,216 permutations of a 10 bp transcription factor binding site (TFBS) known as the CArG box. We developed a robust in silico bioinformatics screening method to evaluate sequences around RefSeq genes for conserved CArG boxes. Utilizing a predetermined phastCons threshold score, we identified 8,252 strand-specific CArGs within an 8 kb window around the transcription start site of 5,213 genes, including all previously defined SRF target genes. We then interrogated this CArG dataset for the presence of previously annotated common polymorphisms. We found a total of 118 unique CArG boxes harboring a SNP within the 10 bp CArG sequence and 1,130 CArG boxes with SNPs located just outside the CArG element. Gel shift and luciferase reporter assays validated SRF binding and functional activity of several new CArG boxes. Importantly, SNPs within or just outside the CArG box often resulted in altered SRF binding and activity. Collectively, these findings demonstrate a powerful approach to computationally define rSNPs in the human CArGome and provide a foundation for similar analyses of other TFBS. Such information may find utility in genetic association studies of human disease where little insight is known regarding the functionality of rSNPs.

Serum response factor-dependent MicroRNAs regulate gastrointestinal smooth muscle cell phenotypes

BACKGROUND & AIMS

Smooth muscle cells (SMCs) change phenotypes under various pathophysiological conditions. These changes are largely controlled by the serum response factor (SRF), a transcription factor that binds to CC (A/T)6 GG (CArG) boxes in SM contractile genes. MicroRNAs (miRNA) regulate transitions among SMC phenotypes. The SMC miRNA transcriptome (SMC miRNAome) and its regulation by SRF have not been determined.

METHODS

We performed massively parallel sequencing to identify gastrointestinal (GI) SMC miRNA transcriptomes in mice and humans. SMC miRNA transcriptomes were mapped to identify all CArG boxes, which were confirmed by SRF knockdown and microarrays. Quantitative polymerase chain reaction was used to identify SMC-phenotypic miRNAs in differentiated and proliferating SMCs. Bioinformatics and target validation analysis showed regulation of SMC phenotype by SRF-dependent, SMC-phenotype miRNAs.

RESULTS

We cloned and identified GI miRNA transcriptomes using genome-wide analyses of mouse and human cells. The SM miRNAome consisted of hundreds of unique miRNAs that were highly conserved among both species. We mapped miRNAs CArG boxes and found that many had an SRF-dependent signature in the SM miRNAome. The SM miRNAs CArG boxes had several distinct features. We also identified approximately 100 SMC-phenotypic miRNAs that were induced in differentiated or proliferative SMC phenotypes. We showed that SRF-dependent, SMC-phenotypic miRNAs bind and regulate Srf and its cofactors, myocadin (Myocd) and member of ETS oncogene family Elk1.

CONCLUSIONS

The GI SMC phenotypes are controlled by SRF-dependent, SMC-phenotypic miRNAs that regulate expression of SRF, MYOCD, and ELK1.

Biological functions and transcriptional targets of CaRF in neurons

Calcium-response factor (CaRF) is a unique DNA-binding protein first recognized as a transcription factor for its role in modulating transcription of the gene encoding Brain-Derived Neurotrophic Factor (BDNF) in neurons. Here I review evidence for the biological functions and transcriptional targets of CaRF in the brain and discuss potential mechanisms by which calcium-activated signaling pathways may modulate CaRF-dependent transcription. These data paint an emerging picture of CaRF as a regulatory nexus for signaling pathways that control aspects of synaptic development and neuronal function.

Vertebrate paralogous MEF2 genes: origin, conservation, and evolution

Background

The myocyte enhancer factor 2 (MEF2) gene family is broadly expressed during the development and maintenance of muscle cells. Although a great deal has been elucidated concerning MEF2 transcription factors' regulation of specific gene expression in diverse programs and adaptive responses, little is known about the origin and evolution of the four members of the MEF2 gene family in vertebrates.

Methodology/Principal Findings

By phylogenetic analyses, we investigated the origin, conservation, and evolution of the four MEF2 genes. First, among the four MEF2 paralogous branches, MEF2B is clearly distant from the other three branches in vertebrates, mainly because it lacks the HJURP_C (Holliday junction recognition protein C-terminal) region. Second, three duplication events might have occurred to produce the four MEF2 paralogous genes and the latest duplication event occurred near the origin of vertebrates producing MEF2A and MEF2C. Third, the ratio (K_a/K_s) of non-synonymous to synonymous nucleotide substitution rates showed that MEF2B evolves faster than the other three MEF2 proteins despite purifying selection on all of the four MEF2 branches. Moreover, a pair model of M0 versus M3 showed that variable selection exists among MEF2 proteins, and branch-site analysis presented that sites 53 and 64 along the MEF2B branch are under positive selection. Finally, and interestingly, substitution rates showed that type II MADS genes (i.e., MEF2-like genes) evolve as slowly as type I MADS genes (i.e., SRF-like genes) in animals, which is inconsistent with the fact that type II MADS genes evolve much slower than type I MADS genes in plants.

Conclusion

Our findings shed light on the relationship of MEF2A, B, C, and D with functional conservation and evolution in vertebrates. This study provides a rationale for future experimental design to investigate distinct but overlapping regulatory roles of the four MEF2 genes in various tissues.

Substitution pattern of the CArG element in human and mouse genomes

cis-Elements CArG bound by serum response factor (SRF) are presently being intensively studied, but little is known about the substitution pattern of functional CArG elements. Here, we have performed the first evolutionary analysis of CArGome in the human and mouse genome through bioinformatic methods and statistical tests. We calculated the substitution rate at each site of the functional CArG elements. The results showed that the core sites of the functional CArG elements evolved faster than did the background DNA, indicating that these sites were likely to evolve under positive selection. Moreover, a strong TATA "motif" was evident in the core region within the functional CArG elements in both human and mouse promoters. This motif could probably be a major contribution to the formation of the spatial structure, which was important for CArG-SRF recognition. Thus, the study further revealed the sequence character and substitution pattern of CArG elements and provided useful information for the study of the SRF-binding efficiencies of CArG promoters in functional assays.

Genome-wide analysis of transcription factor E2F1 mutant proteins reveals that N- and C-terminal protein interaction domains do not participate in targeting E2F1 to the human genome

Previous studies of E2F family members have suggested that protein-protein interactions may be the mechanism by which E2F proteins are recruited to specific genomic regions. We have addressed this hypothesis on a genome-wide scale using ChIP-seq analysis of MCF7 cell lines that express tagged wild type and mutant E2F1 proteins. First, we performed ChIP-seq for tagged WT E2F1. Then, we analyzed E2F1 proteins that lacked the N-terminal SP1 and cyclin A binding domains, the C-terminal transactivation and pocket protein binding domains, and the internal marked box domain. Surprisingly, we found that the ChIP-seq patterns of the mutant proteins were identical to that of WT E2F1. However, mutation of the DNA binding domain abrogated all E2F1 binding to the genome. These results suggested that the interaction between the E2F1 DNA binding domain and a consensus motif may be the primary determinant of E2F1 recruitment. To address this possibility, we analyzed the in vivo binding sites for the in vitro-derived consensus E2F1 motif (TTTSSCGC) and also performed de novo motif analysis. We found that only 12% of the ChIP-seq peaks contained the TTTSSCGC motif. De novo motif analysis indicated that most of the in vivo sites lacked the 5′ half of the in vitro-derived consensus, having instead the in vivo consensus of CGCGC. In summary, our findings do not provide support for the model that protein-protein interactions are involved in recruiting E2F1 to the genome, but rather suggest that recognition of a motif found at most human promoters is the critical determinant.

A systems biological approach to identify key transcription factors and their genomic neighborhoods in human sarcomas

Identification of genetic signatures is the main objective for many computational oncology studies. The signature usually consists of numerous genes that are differentially expressed between two clinically distinct groups of samples, such as tumor subtypes. Prospectively, many signatures have been found to generalize poorly to other datasets and, thus, have rarely been accepted into clinical use. Recognizing the limited success of traditionally generated signatures, we developed a systems biology-based framework for robust identification of key transcription factors and their genomic regulatory neighborhoods. Application of the framework to study the differences between gastrointestinal stromal tumor (GIST) and leiomyosarcoma (LMS) resulted in the identification of nine transcription factors (SRF, NKX2-5, CCDC6, LEF1, VDR, ZNF250, TRIM63, MAF, and MYC). Functional annotations of the obtained neighborhoods identified the biological processes which the key transcription factors regulate differently between the tumor types. Analyzing the differences in the expression patterns using our approach resulted in a more robust genetic signature and more biological insight into the diseases compared to a traditional genetic signature.

Transcription factor binding variation in the evolution of gene regulation

Transcription factor interactions with DNA are one of the primary mechanisms by which expression is modulated, yet their evolution remains poorly understood. Chromatin immunoprecipitation followed by microarray (ChIP-chip) or sequencing (ChIP-Seq) has revolutionized the study of protein-DNA interactions. However, only recently has attention focused on determining to what extent these regulatory interactions vary between species across entire genomes. A series of recent studies have compared in vivo binding data across a range of evolutionary distances. Binding events diverge rapidly, indicating gene regulation is an evolutionarily flexible process.

Role of serum response factor in the pathogenesis of disease

Serum response factor (SRF) is a ubiquitously expressed transcription factor that binds to a DNA cis element known as the CArG box, which is found in the proximal regulatory regions of over 200 experimentally validated target genes. Genetic deletion of SRF is incompatible with life in a variety of animals from different phyla. In mice, loss of SRF throughout the early embryo results in gastrulation defects precluding analyses in individual organ systems. Genetic inactivation studies using conditional or inducible promoters directing the expression of the bacteriophage Cre recombinase have shown a vital role for SRF in such cellular processes as contractility, cell migration, synaptic activity, inflammation, and cell survival. A growing number of experimental and human diseases are associated with changes in SRF expression, suggesting that SRF has a role in the pathogenesis of disease. This review summarizes data from experimental model systems and human pathology where SRF expression is either deliberately or naturally altered.

Conservation and evolution in and among SRF- and MEF2-type MADS domains and their binding sites

Serum response factor (SRF) and myocyte enhancer factor 2 (MEF2) represent two types of members of the MCM1, AGAMOUS, DEFICIENS, and SRF (MADS)-box transcription factor family present in animals and fungi. Each type has distinct biological functions, which are reflected by the distinct specificities of the proteins bound to their cognate DNA-binding sites and activated by their respective cofactors. However, little is known about the evolution of MADS domains and their DNA-binding sites. Here, we report on the conservation and evolution of the two types of MADS domains with their cognate DNA-binding sites by using phylogenetic analyses. First, there are great similarities between the two types of proteins with amino acid positions highly conserved, which are critical for binding to the DNA sequence and for the maintenance of the 3D structure. Second, in contrast to MEF2-type MADS domains, distinct conserved residues are present at some positions in SRF-type MADS domains, determining specificity and the configuration of the MADS domain bound to DNA sequences. Furthermore, the ancestor sequence of SRF- and MEF2-type MADS domains is more similar to MEF2-type MADS domains than to SRF-type MADS domains. In the case of DNA-binding sites, the MEF2 site has a T-rich core in one DNA sequence and an A-rich core in the reverse sequence as compared with the SRF site, no matter whether where either A or T is present in the two complementary sequences. In addition, comparing SRF sites in the human and the mouse genomes reveals that the evolution rate of CArG-boxes is faster in mouse than in human. Moreover, interestingly, a CArG-like sequence, which is probably functionless, could potentially mutate to a functional CArG-box that can be bound by SRF and vice versa. Together, these results significantly improve our knowledge on the conservation and evolution of the MADS domains and their binding sites to date and provide new insights to investigate the MADS family, which is not only on evolution of MADS factors but also on evolution of their binding sites and even on coevolution of MADS factors with their binding sites.

Circadian transcription in liver

Circadian rhythms regulate a wide range of cellular, physiological, metabolic and behavioral activities in mammals. The complexity of tissue- and day-time specific regulation of thousands of clock controlled genes (CCGs) suggests that many transcriptional regulators are involved. Our bioinformatic analysis is based on two published DNA-array studies from mouse liver. We search overrepresented transcription factor binding sites in promoter regions of CCGs using GC-matched controls. Analyzing a large set of CCG promoters, we find known motifs such as E-boxes, D-boxes and cAMP responsive elements. In addition, we find overrepresented GC-rich motifs (Sp1, ETF, Nrf1), AT-rich motifs (TBP, Fox04, MEF-2), Y-box motifs (NF-Y, C/EBP) and cell cycle regulators (E2F, Elk-1). In a subset of system-driven genes, we find overrepresented motifs of the serum response factor SRF and the estrogen receptor ER. The analysis of published ChIP data reveals that some of our predicted regulators (C/EBP, E2F, HNF-1, Myc, MEF-2) target relatively many clock controlled genes. Our analysis of CCG promoters contributes to an understanding of the complex transcriptional regulation of circadian rhythms in liver.

KLF3 regulates muscle-specific gene expression and synergizes with serum response factor on KLF binding sites

This study identifies KLF3 as a transcriptional regulator of muscle genes and reveals a novel synergistic interaction between KLF3 and serum response factor (SRF). Using quantitative proteomics, KLF3 was identified as one of several candidate factors that recognize the MPEX control element in the Muscle creatine kinase (MCK) promoter. Chromatin immunoprecipitation analysis indicated that KLF3 is enriched at many muscle gene promoters (MCK, Myosin heavy chain IIa, Six4, Calcium channel receptor alpha-1, and Skeletal alpha-actin), and two KLF3 isoforms are upregulated during muscle differentiation. KLF3 and SRF physically associate and synergize in transactivating the MCK promoter independently of SRF binding to CArG motifs. The zinc finger and repression domains of KLF3 plus the MADS box and transcription activation domain of SRF are implicated in this synergy. Our results provide the first evidence of a role for KLF3 in muscle gene regulation and reveal an alternate mechanism for transcriptional regulation by SRF via its recruitment to KLF binding sites. Since both factors are expressed in all muscle lineages, SRF may regulate many striated- and smooth-muscle genes that lack known SRF control elements, thus further expanding the breadth of the emerging CArGome.

Genome-wide analysis of ETS-family DNA-binding in vitro and in vivo

Members of the large ETS family of transcription factors (TFs) have highly similar DNA-binding domains (DBDs)—yet they have diverse functions and activities in physiology and oncogenesis. Some differences in DNA-binding preferences within this family have been described, but they have not been analysed systematically, and their contributions to targeting remain largely uncharacterized. We report here the DNA-binding profiles for all human and mouse ETS factors, which we generated using two different methods: a high-throughput microwell-based TF DNA-binding specificity assay, and protein-binding microarrays (PBMs). Both approaches reveal that the ETS-binding profiles cluster into four distinct classes, and that all ETS factors linked to cancer, ERG, ETV1, ETV4 and FLI1, fall into just one of these classes. We identify amino-acid residues that are critical for the differences in specificity between all the classes, and confirm the specificities in vivo using chromatin immunoprecipitation followed by sequencing (ChIP-seq) for a member of each class. The results indicate that even relatively small differences in in vitro binding specificity of a TF contribute to site selectivity in vivo.

Functional analysis and identification of cis-regulatory elements of human chromosome 21 gene promoters

Given the inherent limitations of in silico studies relying solely on DNA sequence analysis, the functional characterization of mammalian promoters and associated cis-regulatory elements requires experimental support, which demands cloning and analysis of putative promoter regions. Focusing on human chromosome 21, we cloned 182 gene promoters of 2500 bp in length and conducted reporter gene assays on transfected-cell arrays. We found 56 promoters that were active in HEK293 cells, while another 49 promoters could be activated by treatment of cells with Trichostatin A or depletion of serum. We observed high correlations between promoter activities and endogenous transcript levels, RNA polymerase II occupancy, CpG islands and core promoter elements. Truncation of a subset of 62 promoters to ∼500 bp revealed that truncation rarely resulted in loss of activity, but rather in loss of responses to external stimuli, suggesting the presence of cis-regulatory response elements within distal promoter regions. In these regions, we found a strong enrichment of transcription factor binding sites that could potentially activate gene expression in the presence of stimuli. This study illustrates the modular functional architecture of chromosome 21 promoters and helps to reveal the complex mechanisms governing transcriptional regulation.

Linking actin dynamics and gene transcription to drive cellular motile functions

Numerous physiological and pathological stimuli promote the rearrangement of the actin cytoskeleton, thereby modulating cellular motile functions. Although it seems intuitively obvious that cell motility requires coordinated protein biosynthesis, until recently the linkage between cytoskeletal actin dynamics and correlated gene activities remained unknown. This knowledge gap was filled in part by the discovery that globular actin polymerization liberates myocardin-related transcription factor (MRTF) cofactors, thereby inducing the nuclear transcription factor serum response factor (SRF) to modulate the expression of genes encoding structural and regulatory effectors of actin dynamics. This insight stimulated research to better understand the actin-MRTF-SRF circuit and to identify alternative mechanisms that link cytoskeletal dynamics and genome activity.

Characteristics of the CArG-SRF binding context in mammalian genomes

Serum response factor (SRF), a member of the MADS family, binds a 10-bp cis element known as the CArG box. However, despite our extensive knowledge of SRF and the CArG box, limited information is available on the CArG-SRF binding context or how CArG flanking sequences are defined. We used statistical tests and computer simulation to find characteristics of the CArG-SRF binding context. Based on the combination of published literature and a search of DBTSS, 150 and 136 functional CArG boxes together with 10 bp flanking on each side were found in mouse and human genomes, respectively. By statistical analysis of the 30 positions we found some new conserved positions of interest (P < 0.005) such as -15, -8, and +8, in addition to the ten highly conserved positions of the CArG box. Intriguingly, studies comparing the flanking positions between consensus CArG boxes and CArG-like boxes showed that there are more conserved positions in the latter. Moreover, CpG content within the CArG-SRF binding context is much higher than that within introns. Collectively, these results suggest that there are some special pre-existing features in the flanking sequences of functional CArG boxes probably contributing to SRF selectively recognizing and binding to the functional CArG from millions of functionless CArG boxes in mammalian genomes. This is a significant step toward understanding the mechanism of transcriptional regulation of SRF-dependent genes.

Insights from genomic profiling of transcription factors

A crucial question in the field of gene regulation is whether the location at which a transcription factor binds influences its effectiveness or the mechanism by which it regulates transcription. Comprehensive transcription factor binding maps are needed to address these issues, and genome-wide mapping is now possible thanks to the technological advances of ChIP-chip and ChIP-seq. This Review discusses how recent genomic profiling of transcription factors gives insight into how binding specificity is achieved and what features of chromatin influence the ability of transcription factors to interact with the genome. It also suggests future experiments that may further our understanding of the causes and consequences of transcription factor-genome interactions.

Functional versatility of transcription factors in the nervous system: the SRF paradigm

Individual transcription factors in the brain frequently display broad functional versatility, thereby controlling multiple cellular outputs. In accordance, neuron-restricted mutagenesis of the murine Srf gene, encoding the transcription factor serum response factor (SRF), revealed numerous SRF functions in the nervous system. First, SRF controls immediate early gene (IEG) activation associated with perception of synaptic activity, learning and memory. Second, processes linked to actin cytoskeletal dynamics are mediated by SRF, such as developmental neuronal migration, outgrowth and pathfinding of neurites, as well as synaptic targeting. Therefore, SRF seems to be instrumental in converting synaptic activity into plasticity-associated structural changes in neuronal connectivities. This highlights the decisive role of SRF in integrating cytoskeletal actin dynamics and nuclear gene expression. Finally, we relate SRF to the multi-functional transcription factor CREB and point out overlapping, distinct and concerted functions of these two transcriptional regulators in the brain.