Unravelling regulatory programs governed by transcription factors (TFs) is fundamental to understanding biological systems. TFCat is a catalog of mouse and human TFs based on a reliable core collection of annotations obtained by expert review of the scientific literature. The collection, including proven and homology-based candidate TFs, is annotated within a function-based taxonomy and DNA-binding proteins are organized within a classification system. All data and user-feedback mechanisms are available at the TFCat portal (http://www.tfcat.ca).

Several years after sequencing the human genome and the mouse genome, much remains to be discovered about the functions of most human and mouse genes. Computational prediction of gene function promises to help focus limited experimental resources on the most likely hypotheses. Several algorithms using diverse genomic data have been applied to this task in model organisms; however, the performance of such approaches in mammals has not yet been evaluated.

A gene network's capacity to process information, so as to bind past events to future actions, depends on its structure and logic. From previous and new microarray measurements in Saccharomyces cerevisiae following gene deletions and overexpressions, we identify a core gene regulatory network (GRN) of functional interactions between 328 genes and the transfer functions of each gene. Inferred connections are verified by gene enrichment.

The localization of memory T cells to human skin is essential for long-term immune surveillance and the maintenance of barrier integrity. The expression of CCR8 during naive T cell activation is controlled by skin-specific factors derived from epidermal keratinocytes and not by resident dendritic cells. In this study, we show that the CCR8-inducing factors are heat stable and protease resistant and include the vitamin D3 metabolite 1α,25-dihydroxyvitamin D3 and PGE2. The effect of either metabolite alone on CCR8 expression was weak, whereas their combination resulted in robust CCR8 expression. Elevation of intracellular cAMP was essential because PGE2 could be substituted with the adenylyl cyclase agonist forskolin, and CCR8 expression was sensitive to protein kinase A inhibition. For effective induction, exposure of naive T cells to these epidermal factors needed to occur either prior to or during T cell activation even though CCR8 was only detected 4-5 d later in proliferating T cells. The importance of tissue environments in maintaining cellular immune surveillance networks within distinct healthy tissues provides a paradigm shift in adaptive immunity. Epidermal-derived vitamin D3 metabolites and PGs provide an essential cue for the localization of CCR8(+) immune surveillance T cells within healthy human skin.

RNA-binding proteins (RBPs) are important regulators of eukaryotic gene expression. Genomes typically encode dozens to hundreds of proteins containing RNA-binding domains, which collectively recognize diverse RNA sequences and structures. Recent advances in high-throughput methods for assaying the targets of RBPs in vitro and in vivo allow large-scale derivation of RNA-binding motifs as well as determination of RNA-protein interactions in living cells. In parallel, many computational methods have been developed to analyze and interpret these data. The interplay between RNA secondary structure and RBP binding has also been a growing theme. Integrating RNA-protein interaction data with observations of post-transcriptional regulation will enhance our understanding of the roles of these important proteins.

Alternative splicing is important for the development and function of the nervous system, but little is known about the differences in alternative splicing between distinct types of neurons. Furthermore, the factors that control cell-type-specific splicing and the physiological roles of these alternative isoforms are unclear. By monitoring alternative splicing at single-cell resolution in Caenorhabditis elegans, we demonstrate that splicing patterns in different neurons are often distinct and highly regulated. We identify two conserved RNA-binding proteins, UNC-75/CELF and EXC-7/Hu/ELAV, which regulate overlapping networks of splicing events in GABAergic and cholinergic neurons. We use the UNC-75 exon network to discover regulators of synaptic transmission and to identify unique roles for isoforms of UNC-64/Syntaxin, a protein required for synaptic vesicle fusion. Our results indicate that combinatorial regulation of alternative splicing in distinct neurons provides a mechanism to specialize metazoan nervous systems.

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.

Secondary metabolites are produced by numerous organisms and can either be beneficial, benign, or harmful to humans. Genes involved in the synthesis and transport of these secondary metabolites are frequently found in gene clusters, which are often coordinately regulated, being almost exclusively dependent on transcription factors that are located within the clusters themselves. Gliotoxin, which is produced by a variety of Aspergillus species, Trichoderma species, and Penicillium species, exhibits immunosuppressive properties and has therefore been the subject of research for many laboratories. There have been a few proteins shown to regulate the gliotoxin cluster, most notably GliZ, a Zn2Cys6 binuclear finger transcription factor that lies within the cluster, and LaeA, a putative methyltransferase that globally regulates secondary metabolism clusters within numerous fungal species. Using a high-copy inducer screen in A. fumigatus, our lab has identified a novel C2H2 transcription factor, which plays an important role in regulating the gliotoxin biosynthetic cluster. This transcription factor, named GipA, induces gliotoxin production when present in extra copies. Furthermore, loss of gipA reduces gliotoxin production significantly. Through protein binding microarray and mutagenesis, we have identified a DNA binding site recognized by GipA that is in extremely close proximity to a potential GliZ DNA binding site in the 5' untranslated region of gliA, which encodes an efflux pump within the gliotoxin cluster. Not surprisingly, GliZ and GipA appear to work in an interdependent fashion to positively control gliA expression.

C2H2 zinc finger proteins represent the largest and most enigmatic class of human transcription factors. Their C2H2-ZF arrays are highly variable, indicating that most will have unique DNA binding motifs. However, most of the binding motifs have not been directly determined. In addition, little is known about whether or how these proteins regulate transcription. Most of the ∼700 human C2H2-ZF proteins also contain at least one KRAB, SCAN, BTB, or SET domain, suggesting that they may have common interacting partners and/or effector functions. Here, we report a multifaceted functional analysis of 131 human C2H2-ZF proteins, encompassing DNA binding sites, interacting proteins, and transcriptional response to genetic perturbation. We confirm the expected diversity in DNA binding motifs and genomic binding sites, and provide motif models for 78 previously uncharacterized C2H2-ZF proteins, most of which are unique. Surprisingly, the diversity in protein-protein interactions is nearly as high as diversity in DNA binding motifs: Most C2H2-ZF proteins interact with a unique spectrum of co-activators and co-repressors. Thus, multiparameter diversification likely underlies the evolutionary success of this large class of human proteins.

Development of an accurate protein-DNA recognition code that can predict DNA specificity from protein sequence is a central problem in biology. C2H2 zinc fingers constitute by far the largest family of DNA binding domains and their binding specificity has been studied intensively. However, despite decades of research, accurate prediction of DNA specificity remains elusive. A major obstacle is thought to be the inability of current methods to account for the influence of neighbouring domains. Here we show that this problem can be addressed using a structural approach: we build structural models for all C2H2-ZF-DNA complexes with known binding motifs and find six distinct binding modes. Each mode changes the orientation of specificity residues with respect to the DNA, thereby modulating base preference. Most importantly, the structural analysis shows that residues at the domain interface strongly and predictably influence the binding mode, and hence specificity. Accounting for predicted binding mode significantly improves prediction accuracy of predicted motifs. This new insight into the fundamental behaviour of C2H2-ZFs has implications for both improving the prediction of natural zinc finger-binding sites, and for prioritizing further experiments to complete the code. It also provides a new design feature for zinc finger engineering.

The long noncoding RNA Gomafu/MIAT/Rncr2 is thought to function in retinal cell specification, stem cell differentiation and the control of alternative splicing. To further investigate physiological functions of Gomafu, we created mouse knockout (KO) model that completely lacks the Gomafu gene. The KO mice did not exhibit any developmental deficits. However, behavioral tests revealed that the KO mice are hyperactive. This hyperactive behavior was enhanced when the KO mice were treated with the psychostimulant methamphetamine, which was associated with an increase in dopamine release in the nucleus accumbens. RNA sequencing analyses identified a small number of genes affected by the deficiency of Gomafu, a subset of which are known to have important neurobiological functions. These observations suggest that Gomafu modifies mouse behavior thorough a mild modulation of gene expression and/or alternative splicing of target genes.

Acinus (apoptotic chromatin condensation inducer in the nucleus) is an RNA-binding protein (RBP) originally identified for its role in apoptosis. It was later found to be an auxiliary component of the exon junction complex (EJC), which is deposited at exon junctions as a consequence of pre-mRNA splicing. To uncover the cellular functions of Acinus and investigate its role in splicing, we mapped its endogenous RNA targets using the cross-linking immunoprecipitation protocol (iCLIP). We observed that Acinus binds to pre-mRNAs, associating specifically to a subset of suboptimal introns, but also to spliced mRNAs. We also confirmed the presence of Acinus as a peripheral factor of the EJC. RNA-seq was used to investigate changes in gene expression and alternative splicing following siRNA-mediated depletion of Acinus in HeLa cells. This analysis revealed that Acinus is preferentially required for the inclusion of specific alternative cassette exons and also controls the faithful splicing of a subset of introns. Moreover, a large number of splicing changes can be related to Acinus binding, suggesting a direct role of Acinus in exon and intron definition. In particular, Acinus regulates the splicing of DFFA/ICAD transcript, a major regulator of DNA fragmentation. Globally, the genome-wide identification of RNA targets of Acinus revealed its role in splicing regulation as well as its involvement in other cellular pathways, including cell cycle progression. Altogether, this study uncovers new cellular functions of an RBP transiently associated with the EJC.

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.

We describe the application of a microarray platform, which combines information from exon body and splice-junction probes, to perform a quantitative analysis of tissue-specific alternative splicing (AS) for thousands of exons in mammalian cells. Through this system, we have analyzed global features of AS in major mouse tissues. The results provide numerous inferences for the functions of tissue-specific AS, insights into how the evolutionary history of exons can impact on their inclusion levels, and also information on how global regulatory properties of AS define tissue type. Like global transcription profiles, global AS profiles reflect tissue identity. Interestingly, we find that transcription and AS act independently on different sets of genes in order to define tissue-specific expression profiles. These results demonstrate the utility of our quantitative microarray platform and data for revealing important global regulatory features of AS.

Although subcellular mRNA trafficking has been demonstrated as a mechanism to control protein distribution, it is generally believed that most protein localization occurs subsequent to translation. To address this point, we developed and employed a high-resolution fluorescent in situ hybridization procedure to comprehensively evaluate mRNA localization dynamics during early Drosophila embryogenesis. Surprisingly, of the 3370 genes analyzed, 71% of those expressed encode subcellularly localized mRNAs. Dozens of new and striking localization patterns were observed, implying an equivalent variety of localization mechanisms. Tight correlations between mRNA distribution and subsequent protein localization and function, indicate major roles for mRNA localization in nucleating localized cellular machineries. A searchable web resource documenting mRNA expression and localization dynamics has been established and will serve as an invaluable tool for dissecting localization mechanisms and for predicting gene functions and interactions.

In animals, egg activation triggers a cascade of posttranscriptional events that act on maternally synthesized RNAs. We show that, in Drosophila, the PAN GU (PNG) kinase sits near the top of this cascade, triggering translation of SMAUG (SMG), a multifunctional posttranscriptional regulator conserved from yeast to humans. Although PNG is required for cytoplasmic polyadenylation of smg mRNA, it regulates translation via mechanisms that are independent of its effects on the poly(A) tail. Analyses of mutants suggest that PNG relieves translational repression by PUMILIO (PUM) and one or more additional factors, which act in parallel through the smg mRNA's 3' untranslated region (UTR). Microarray-based gene expression profiling shows that SMG is a major regulator of maternal transcript destabilization. SMG-dependent mRNAs are enriched for gene ontology annotations for function in the cell cycle, suggesting a possible causal relationship between failure to eliminate these transcripts and the cell cycle defects in smg mutants.

Although acute lung injury (ALI) is a common complication of severe malaria, little is known about the underlying molecular basis of lung dysfunction. Animal models have provided powerful insights into the pathogenesis of severe malaria syndromes such as cerebral malaria (CM); however, no model of malaria-induced lung injury has been definitively established. This study used bronchoalveolar lavage (BAL), histopathology and gene expression analysis to examine the development of ALI in mice infected with Plasmodium berghei ANKA (PbA). BAL fluid of PbA-infected C57BL/6 mice revealed a significant increase in IgM and total protein prior to the development of CM, indicating disruption of the alveolar-capillary membrane barrier-the physiological hallmark of ALI. In contrast to sepsis-induced ALI, BAL fluid cell counts remained constant with no infiltration of neutrophils. Histopathology showed septal inflammation without cellular transmigration into the alveolar spaces. Microarray analysis of lung tissue from PbA-infected mice identified a significant up-regulation of expressed genes associated with the gene ontology categories of defense and immune response. Severity of malaria-induced ALI varied in a panel of inbred mouse strains, and development of ALI correlated with peripheral parasite burden but not CM susceptibility. Cd36(-/-) mice, which have decreased parasite lung sequestration, were relatively protected from ALI. In summary, parasite burden and CD36-mediated sequestration in the lung are primary determinants of ALI in experimental murine malaria. Furthermore, differential susceptibility of mouse strains to malaria-induced ALI and CM suggests that distinct genetic determinants may regulate susceptibility to these two important causes of malaria-associated morbidity and mortality.

Alternative splicing (AS) is a widespread process underlying the generation of transcriptomic and proteomic diversity and is frequently misregulated in human disease. Accordingly, an important goal of biomedical research is the development of tools capable of comprehensively, accurately, and efficiently profiling AS. Here, we describe Whippet, an easy-to-use RNA-seq analysis method that rapidly-with hardware requirements compatible with a laptop-models and quantifies AS events of any complexity without loss of accuracy. Using an entropic measure of splicing complexity, Whippet reveals that one-third of human protein coding genes produce transcripts with complex AS events involving co-expression of two or more principal splice isoforms. We observe that high-entropy AS events are more prevalent in tumor relative to matched normal tissues and correlate with increased expression of proto-oncogenic splicing factors. Whippet thus affords the rapid and accurate analysis of AS events of any complexity, and as such will facilitate future biomedical research.

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.

Alternative splicing (AS) generates remarkable regulatory and proteomic complexity in metazoans. However, the functions of most AS events are not known, and programs of regulated splicing remain to be identified. To address these challenges, we describe the Vertebrate Alternative Splicing and Transcription Database (VastDB), the largest resource of genome-wide, quantitative profiles of AS events assembled to date. VastDB provides readily accessible quantitative information on the inclusion levels and functional associations of AS events detected in RNA-seq data from diverse vertebrate cell and tissue types, as well as developmental stages. The VastDB profiles reveal extensive new intergenic and intragenic regulatory relationships among different classes of AS and previously unknown and conserved landscapes of tissue-regulated exons. Contrary to recent reports concluding that nearly all human genes express a single major isoform, VastDB provides evidence that at least 48% of multiexonic protein-coding genes express multiple splice variants that are highly regulated in a cell/tissue-specific manner, and that >18% of genes simultaneously express multiple major isoforms across diverse cell and tissue types. Isoforms encoded by the latter set of genes are generally coexpressed in the same cells and are often engaged by translating ribosomes. Moreover, they are encoded by genes that are significantly enriched in functions associated with transcriptional control, implying they may have an important and wide-ranging role in controlling cellular activities. VastDB thus provides an unprecedented resource for investigations of AS function and regulation.