Search

Models of mammalian regulatory networks controlling gene expression have been inferred from genomic data but have largely not been validated. We present an unbiased strategy to systematically perturb candidate regulators and monitor cellular transcriptional responses. We applied this approach to derive regulatory networks that control the transcriptional response of mouse primary dendritic cells to pathogens. Our approach revealed the regulatory functions of 125 transcription factors, chromatin modifiers, and RNA binding proteins, which enabled the construction of a network model consisting of 24 core regulators and 76 fine-tuners that help to explain how pathogen-sensing pathways achieve specificity. This study establishes a broadly applicable, comprehensive, and unbiased approach to reveal the wiring and functions of a regulatory network controlling a major transcriptional response in primary mammalian cells.

We recently showed that the mammalian genome encodes >1,000 large intergenic noncoding (linc)RNAs that are clearly conserved across mammals and, thus, functional. Gene expression patterns have implicated these lincRNAs in diverse biological processes, including cell-cycle regulation, immune surveillance, and embryonic stem cell pluripotency. However, the mechanism by which these lincRNAs function is unknown. Here, we expand the catalog of human lincRNAs to approximately 3,300 by analyzing chromatin-state maps of various human cell types. Inspired by the observation that the well-characterized lincRNA HOTAIR binds the polycomb repressive complex (PRC)2, we tested whether many lincRNAs are physically associated with PRC2. Remarkably, we observe that approximately 20% of lincRNAs expressed in various cell types are bound by PRC2, and that additional lincRNAs are bound by other chromatin-modifying complexes. Also, we show that siRNA-mediated depletion of certain lincRNAs associated with PRC2 leads to changes in gene expression, and that the up-regulated genes are enriched for those normally silenced by PRC2. We propose a model in which some lincRNAs guide chromatin-modifying complexes to specific genomic loci to regulate gene expression.

The regulated recruitment and differentiation of multipotent bone marrow-derived cells (BMDCs) to sites of injury are critical for efficient wound healing. Previously we demonstrated that sustained expression of HOXA3 both accelerated wound healing and promoted angiogenesis in diabetic mice. In this study, we have used green fluorescent protein-positive bone marrow chimeras to investigate the effect of HOXA3 expression on recruitment of BMDCs to wounds. We hypothesized that the enhanced neovascularization induced by HOXA3 is due to enhanced mobilization, recruitment, and/or differentiation of BMDCs. Here we show that diabetic mice treated with HOXA3 displayed a significant increase in both mobilization and recruitment of endothelial progenitor cells compared with control mice. Importantly, we also found that HOXA3-treated mice had significantly fewer inflammatory cells recruited to the wound compared with control mice. Microarray analyses of HOXA3-treated wounds revealed that indeed HOXA3 locally increased expression of genes that selectively promote stem/progenitor cell mobilization and recruitment while also suppressing expression of numerous members of the proinflammatory nuclear factor kappaB pathway, including myeloid differentiation primary response gene 88 and toll-interacting protein. Thus HOXA3 accelerates wound repair by mobilizing endothelial progenitor cells and attenuating the excessive inflammatory response of chronic wounds.

There is growing recognition that mammalian cells produce many thousands of large intergenic transcripts. However, the functional significance of these transcripts has been particularly controversial. Although there are some well-characterized examples, most (>95%) show little evidence of evolutionary conservation and have been suggested to represent transcriptional noise. Here we report a new approach to identifying large non-coding RNAs using chromatin-state maps to discover discrete transcriptional units intervening known protein-coding loci. Our approach identified approximately 1,600 large multi-exonic RNAs across four mouse cell types. In sharp contrast to previous collections, these large intervening non-coding RNAs (lincRNAs) show strong purifying selection in their genomic loci, exonic sequences and promoter regions, with greater than 95% showing clear evolutionary conservation. We also developed a functional genomics approach that assigns putative functions to each lincRNA, demonstrating a diverse range of roles for lincRNAs in processes from embryonic stem cell pluripotency to cell proliferation. We obtained independent functional validation for the predictions for over 100 lincRNAs, using cell-based assays. In particular, we demonstrate that specific lincRNAs are transcriptionally regulated by key transcription factors in these processes such as p53, NFkappaB, Sox2, Oct4 (also known as Pou5f1) and Nanog. Together, these results define a unique collection of functional lincRNAs that are highly conserved and implicated in diverse biological processes.

The genetic programs that promote retention of self-renewing leukemia stem cells (LSCs) at the apex of cellular hierarchies in acute myeloid leukemia (AML) are not known. In a mouse model of human AML, LSCs exhibit variable frequencies that correlate with the initiating MLL oncogene and are maintained in a self-renewing state by a transcriptional subprogram more akin to that of embryonic stem cells (ESCs) than to that of adult stem cells. The transcription/chromatin regulatory factors Myb, Hmgb3, and Cbx5 are critical components of the program and suffice for Hoxa/Meis-independent immortalization of myeloid progenitors when coexpressed, establishing the cooperative and essential role of an ESC-like LSC maintenance program ancillary to the leukemia-initiating MLL/Hox/Meis program. Enriched expression of LSC maintenance and ESC-like program genes in normal myeloid progenitors and poor-prognosis human malignancies links the frequency of aberrantly self-renewing progenitor-like cancer stem cells (CSCs) to prognosis in human cancer.

Homozygous loss-of-function mutations in AIRE cause autoimmune polyglandular syndrome type 1 (APS 1), which manifests in a classic triad of hypoparathyroidism, adrenal insufficiency, and candidiasis. Interestingly, a kindred with a specific G228W AIRE variant presented with an autosomal dominant autoimmune phenotype distinct from APS 1. We utilized a novel G228W-knockin mouse model to show that this variant acted in a dominant-negative manner to cause a unique autoimmunity syndrome. In addition, the expression of a large number of Aire-regulated thymic antigens was partially inhibited in these animals, demonstrating the importance of quantitative changes in thymic antigen expression in determining organ-specific autoimmunity. Furthermore, the dominant-negative effect of the G228W variant was exerted through recruitment of WT Aire away from active sites of transcription in the nucleus of medullary thymic epithelial cells in vivo. Together, these results may demonstrate a mechanism by which autoimmune predisposition to phenotypes distinct from APS 1 can be mediated in a dominant-negative fashion by Aire.

Human skin exhibits exquisite site-specific morphologies and functions. How are these site-specific differences specified during development, maintained in adult homeostasis, and potentially perturbed by disease processes? Here, we review progress in understanding the anatomic patterning of fibroblasts, a major constituent cell type of the dermis and key participant in epithelial-mesenchymal interactions. The gene expression programs of human fibroblasts largely reflect the superimposition of three gene expression profiles that demarcate the fibroblast's position relative to three developmental axes. The HOX family of homeodomain transcription factors is implicated in specifying site-specific transcriptional programs. The use of gene, tiling, and tissue microarrays together gives a comprehensive view of the gene regulation involved in patterning the skin.

Reciprocal epithelial-mesenchymal interactions shape site-specific development of skin. Here we show that site-specific HOX expression in fibroblasts is cell-autonomous and epigenetically maintained. The distal-specific gene HOXA13 is continually required to maintain the distal-specific transcriptional program in adult fibroblasts, including expression of WNT5A, a morphogen required for distal development. The ability of distal fibroblasts to induce epidermal keratin 9, a distal-specific gene, is abrogated by depletion of HOXA13, but rescued by addition of WNT5A. Thus, maintenance of appropriate HOX transcriptional program in adult fibroblasts may serve as a source of positional memory to differentially pattern the epithelia during homeostasis and regeneration.

The recent discovery of a large number of histone demethylases suggests a central role for these enzymes in regulating histone methylation dynamics. Histone H3K27 trimethylation (H3K27me3) has been linked to polycomb-group-protein-mediated suppression of Hox genes and animal body patterning, X-chromosome inactivation and possibly maintenance of embryonic stem cell (ESC) identity. An imbalance of H3K27 methylation owing to overexpression of the methylase EZH2 has been implicated in metastatic prostate and aggressive breast cancers. Here we show that the JmjC-domain-containing related proteins UTX and JMJD3 catalyse demethylation of H3K27me3/2. UTX is enriched around the transcription start sites of many HOX genes in primary human fibroblasts, in which HOX genes are differentially expressed, but is selectively excluded from the HOX loci in ESCs, in which HOX genes are largely silent. Consistently, RNA interference inhibition of UTX led to increased H3K27me3 levels at some HOX gene promoters. Importantly, morpholino oligonucleotide inhibition of a zebrafish UTX homologue resulted in mis-regulation of hox genes and a striking posterior developmental defect, which was partially rescued by wild-type, but not by catalytically inactive, human UTX. Taken together, these findings identify a small family of H3K27 demethylases with important, evolutionarily conserved roles in H3K27 methylation regulation and in animal anterior-posterior development.

Noncoding RNAs (ncRNA) participate in epigenetic regulation but are poorly understood. Here we characterize the transcriptional landscape of the four human HOX loci at five base pair resolution in 11 anatomic sites and identify 231 HOX ncRNAs that extend known transcribed regions by more than 30 kilobases. HOX ncRNAs are spatially expressed along developmental axes and possess unique sequence motifs, and their expression demarcates broad chromosomal domains of differential histone methylation and RNA polymerase accessibility. We identified a 2.2 kilobase ncRNA residing in the HOXC locus, termed HOTAIR, which represses transcription in trans across 40 kilobases of the HOXD locus. HOTAIR interacts with Polycomb Repressive Complex 2 (PRC2) and is required for PRC2 occupancy and histone H3 lysine-27 trimethylation of HOXD locus. Thus, transcription of ncRNA may demarcate chromosomal domains of gene silencing at a distance; these results have broad implications for gene regulation in development and disease states.

We performed a genome-scale chromatin immunoprecipitation (ChIP)-chip comparison of two modifications (trimethylation of lysine 9 [H3me3K9] and trimethylation of lysine 27 [H3me3K27]) of histone H3 in Ntera2 testicular carcinoma cells and in three different anatomical sources of primary human fibroblasts. We found that in each of the cell types the two modifications were differentially enriched at the promoters of the two largest classes of transcription factors. Specifically, zinc finger (ZNF) genes were bound by H3me3K9 and homeobox genes were bound by H3me3K27. We have previously shown that the Polycomb repressive complex 2 is responsible for mediating trimethylation of lysine 27 of histone H3 in human cancer cells. In contrast, there is little overlap between H3me3K9 targets and components of the Polycomb repressive complex 2, suggesting that a different histone methyltransferase is responsible for the H3me3K9 modification. Previous studies have shown that SETDB1 can trimethylate H3 on lysine 9, using in vitro or artificial tethering assays. SETDB1 is thought to be recruited to chromatin by complexes containing the KAP1 corepressor. To determine if a KAP1-containing complex mediates trimethylation of the identified H3me3K9 targets, we performed ChIP-chip assays and identified KAP1 target genes using human 5-kb promoter arrays. We found that a large number of genes of ZNF transcription factors were bound by both KAP1 and H3me3K9 in normal and cancer cells. To expand our studies of KAP1, we next performed a complete genomic analysis of KAP1 binding using a 38-array tiling set, identifying ~7,000 KAP1 binding sites. The identified KAP1 targets were highly enriched for C2H2 ZNFs, especially those containing Krüppel-associated box (KRAB) domains. Interestingly, although most KAP1 binding sites were within core promoter regions, the binding sites near ZNF genes were greatly enriched within transcribed regions of the target genes. Because KAP1 is recruited to the DNA via interaction with KRAB-ZNF proteins, we suggest that expression of KRAB-ZNF genes may be controlled via an auto-regulatory mechanism involving KAP1., Competing Interests: Competing interests. NimbleGen Systems microarrays were used in this study. One of the authors, RG, is an employee of this company. Also, PJF and RG are co-PIs on an NIH grant that partially funded this work.

Fibroblasts are ubiquitous mesenchymal cells with many vital functions during development, tissue repair, and disease. Fibroblasts from different anatomic sites have distinct and characteristic gene expression patterns, but the principles that govern their molecular specialization are poorly understood. Spatial organization of cellular differentiation may be achieved by unique specification of each cell type; alternatively, organization may arise by cells interpreting their position along a coordinate system. Here we test these models by analyzing the genome-wide gene expression profiles of primary fibroblast populations from 43 unique anatomical sites spanning the human body. Large-scale differences in the gene expression programs were related to three anatomic divisions: anterior-posterior (rostral-caudal), proximal-distal, and dermal versus nondermal. A set of 337 genes that varied according to these positional divisions was able to group all 47 samples by their anatomic sites of origin. Genes involved in pattern formation, cell-cell signaling, and matrix remodeling were enriched among this minimal set of positional identifier genes. Many important features of the embryonic pattern of HOX gene expression were retained in fibroblasts and were confirmed both in vitro and in vivo. Together, these findings suggest that site-specific variations in fibroblast gene expression programs are not idiosyncratic but rather are systematically related to their positional identities relative to major anatomic axes., Competing Interests: Competing interests. The authors have declared that no competing interests exist.

Males and females have obvious phenotypic differences; they also exhibit differences related to health, life span, cognitive abilities and have different responses to diseases such as anemia, coronary heart disease, hypertension and renal dysfunction. Although the anatomical, hormonal and chemical differences between the sexes are well known, there are few molecular descriptors for gender-specific physiological traits and health risks. Recent studies using microarrays and other methods have made significant progress towards elucidating the molecular differences between mammalian sexes in a variety of tissues and towards identifying the transcription factors that regulate sex-biased gene expression. These findings are providing new insights into the molecular and genetic differences that dictate the different behaviors and physiologies of mammalian sexes.

Elucidating the transcribed regions of the genome constitutes a fundamental aspect of human biology, yet this remains an outstanding problem. To comprehensively identify coding sequences, we constructed a series of high-density oligonucleotide tiling arrays representing sense and antisense strands of the entire nonrepetitive sequence of the human genome. Transcribed sequences were located across the genome via hybridization to complementary DNA samples, reverse-transcribed from polyadenylated RNA obtained from human liver tissue. In addition to identifying many known and predicted genes, we found 10,595 transcribed sequences not detected by other methods. A large fraction of these are located in intergenic regions distal from previously annotated genes and exhibit significant homology to other mammalian proteins.

Many anatomical differences exist between males and females; these are manifested on a molecular level by different hormonal environments. Although several molecular differences in adult tissues have been identified, a comprehensive investigation of the gene expression differences between males and females has not been performed. We surveyed the expression patterns of 13,977 mouse genes in male and female hypothalamus, kidney, liver, and reproductive tissues. Extensive differential gene expression was observed not only in the reproductive tissues, but also in the kidney and liver. The differentially expressed genes are involved in drug and steroid metabolism, osmotic regulation, or as yet unresolved cellular roles. In contrast, very few molecular differences were observed between the male and female hypothalamus in both mice and humans. We conclude that there are persistent differences in gene expression between adult males and females. These molecular differences have important implications for the physiological differences between males and females.

The cyclic AMP-responsive element-binding protein (CREB) is an important transcription factor that can be activated by hormonal stimulation and regulates neuronal function and development. An unbiased, global analysis of where CREB binds has not been performed. We have mapped for the first time the binding distribution of CREB along an entire human chromosome. Chromatin immunoprecipitation of CREB-associated DNA and subsequent hybridization of the associated DNA to a genomic DNA microarray containing all of the nonrepetitive DNA of human chromosome 22 revealed 215 binding sites corresponding to 192 different loci and 100 annotated potential gene targets. We found binding near or within many genes involved in signal transduction and neuronal function. We also found that only a small fraction of CREB binding sites lay near well-defined 5' ends of genes; the majority of sites were found elsewhere, including introns and unannotated regions. Several of the latter lay near novel unannotated transcriptionally active regions. Few CREB targets were found near full-length cyclic AMP response element sites; the majority contained shorter versions or close matches to this sequence. Several of the CREB targets were altered in their expression by treatment with forskolin; interestingly, both induced and repressed genes were found. Our results provide novel molecular insights into how CREB mediates its functions in humans.

We have mapped the chromosomal binding site distribution of a transcription factor in human cells. The NF-kappaB family of transcription factors plays an essential role in regulating the induction of genes involved in several physiological processes, including apoptosis, immunity, and inflammation. The binding sites of the NF-kappaB family member p65 were determined by using chromatin immunoprecipitation and a genomic microarray of human chromosome 22 DNA. Sites of binding were observed along the entire chromosome in both coding and noncoding regions, with an enrichment at the 5' end of genes. Strikingly, a significant proportion of binding was seen in intronic regions, demonstrating that transcription factor binding is not restricted to promoter regions. NF-kappaB binding was also found at genes whose expression was regulated by tumor necrosis factor alpha, a known inducer of NF-kappaB-dependent gene expression, as well as adjacent to genes whose expression is not affected by tumor necrosis factor alpha. Many of these latter genes are either known to be activated by NF-kappaB under other conditions or are consistent with NF-kappaB's role in the immune and apoptotic responses. Our results suggest that binding is not restricted to promoter regions and that NF-kappaB binding occurs at a significant number of genes whose expression is not altered, thereby suggesting that binding alone is not sufficient for gene activation.

A DNA microarray representing nearly all of the unique sequences of human Chromosome 22 was constructed and used to measure global-transcriptional activity in placental poly(A)(+) RNA. We found that many of the known, related and predicted genes are expressed. More importantly, our study reveals twice as many transcribed bases as have been reported previously. Many of the newly discovered expressed fragments were verified by RNA blot analysis and a novel technique called differential hybridization mapping (DHM). Interestingly, a significant fraction of these novel fragments are expressed antisense to previously annotated introns. The coding potential of these novel expressed regions is supported by their sequence conservation in the mouse genome. This study has greatly increased our understanding of the biological information encoded on a human chromosome. To facilitate the dissemination of these results to the scientific community, we have developed a comprehensive Web resource to present the findings of this study and other features of human Chromosome 22 at http://array.mbb.yale.edu/chr22.