Principles of regulatory information conservation between mouse and human

Anatomical distributional defects in mutant genes associated with dominant intermediate Charcot-Marie-Tooth disease type C in an adenovirus-mediated mouse model

Dominant intermediate Charcot-Marie-Tooth disease type C (DI-CMTC) is a dominantly inherited neuropathy that has been classified primarily based on motor conduction velocity tests but is now known to involve axonal and demyelination features. DI-CMTC is linked to tyrosyl-tRNA synthetase (YARS)-associated neuropathies, which are caused by E196K and G41R missense mutations and a single de novo deletion (153-156delVKQV). It is well-established that these YARS mutations induce neuronal dysfunction, morphological symptoms involving axonal degeneration, and impaired motor performance. The present study is the first to describe a novel mouse model of YARS-mutation-induced neuropathy involving a neuron-specific promoter with a deleted mitochondrial targeting sequence that inhibits the expression of YARS protein in the mitochondria. An adenovirus vector system and in vivo techniques were utilized to express YARS fusion proteins with a Flag-tag in the spinal cord, peripheral axons, and dorsal root ganglia. Following transfection of YARS-expressing viruses, the distributions of wild-type (WT) YARS and E196K mutant proteins were compared in all expressed regions; G41R was not expressed. The proportion of Flag/green fluorescent protein (GFP) double-positive signaling in the E196K mutant-type mice did not significantly differ from that of WT mice in dorsal root ganglion neurons. All adenovirus genes, and even the empty vector without the YARS gene, exhibited GFP-positive signaling in the ventral horn of the spinal cord because GFP in an adenovirus vector is driven by a cytomegalovirus promoter. The present study demonstrated that anatomical differences in tissue can lead to dissimilar expressions of YARS genes. Thus, use of this novel animal model will provide data regarding distributional defects between mutant and WT genes in neurons, the DI-CMTC phenotype, and potential treatment approaches for this disease.

Evolved Repression Overcomes Enhancer Robustness

Biological systems display extraordinary robustness. Robustness of transcriptional enhancers results mainly from clusters of binding sites for the same transcription factor, and it is not clear how robust enhancers can evolve loss of expression through point mutations. Here, we report the high-resolution functional dissection of a robust enhancer of the shavenbaby gene that has contributed to morphological evolution. We found that robustness is encoded by many binding sites for the transcriptional activator Arrowhead and that, during evolution, some of these activator sites were lost, weakening enhancer activity. Complete silencing of enhancer function, however, required evolution of a binding site for the spatially restricted potent repressor Abrupt. These findings illustrate that recruitment of repressor binding sites can overcome enhancer robustness and may minimize pleiotropic consequences of enhancer evolution. Recruitment of repression may be a general mode of evolution to break robust regulatory linkages.

A gene-centered C. elegans protein-DNA interaction network provides a framework for functional predictions

Transcription factors (TFs) play a central role in controlling spatiotemporal gene expression and the response to environmental cues. A comprehensive understanding of gene regulation requires integrating physical protein–DNA interactions (PDIs) with TF regulatory activity, expression patterns, and phenotypic data. Although great progress has been made in mapping PDIs using chromatin immunoprecipitation, these studies have only characterized ~10% of TFs in any metazoan species. The nematode C. elegans has been widely used to study gene regulation due to its compact genome with short regulatory sequences. Here, we delineated the largest gene‐centered metazoan PDI network to date by examining interactions between 90% of C. elegans TFs and 15% of gene promoters. We used this network as a backbone to predict TF binding sites for 77 TFs, two‐thirds of which are novel, as well as integrate gene expression, protein–protein interaction, and phenotypic data to predict regulatory and biological functions for multiple genes and TFs.

Past Roadblocks and New Opportunities in Transcription Factor Network Mapping

One of the principal mechanisms by which cells differentiate and respond to changes in external signals or conditions is by changing the activity levels of transcription factors (TFs). This changes the transcription rates of target genes via the cell's TF network, which ultimately contributes to reconfiguring cellular state. Since microarrays provided our first window into global cellular state, computational biologists have eagerly attacked the problem of mapping TF networks, a key part of the cell's control circuitry. In retrospect, however, steady-state mRNA abundance levels were a poor substitute for TF activity levels and gene transcription rates. Likewise, mapping TF binding through chromatin immunoprecipitation proved less predictive of functional regulation and less amenable to systematic elucidation of complete networks than originally hoped. This review explains these roadblocks and the current, unprecedented blossoming of new experimental techniques built on second-generation sequencing, which hold out the promise of rapid progress in TF network mapping.

A Mimic of the Tumor Microenvironment: A Simple Method for Generating Enriched Cell Populations and Investigating Intercellular Communication

Understanding the early heterotypic interactions between cancer cells and the surrounding non-cancerous stroma is important in elucidating the events leading to stromal activation and establishment of the tumor microenvironment (TME). Several in vitro and in vivo models of the TME have been developed; however, in general these models do not readily permit isolation of individual cell populations, under non-perturbing conditions, for further study. To circumvent this difficulty, we have employed an in vitro TME model using a cell growth substrate consisting of a permeable microporous membrane insert that permits simple generation of highly enriched cell populations grown intimately, yet separately, on either side of the insert's membrane for extended co-culture times. Through use of this model, we are capable of generating greatly enriched cancer-associated fibroblast (CAF) populations from normal diploid human fibroblasts following co-culture (120 hr) with highly metastatic human breast carcinoma cells, without the use of fluorescent tagging and/or cell sorting. Additionally, by modulating the pore-size of the insert, we can control for the mode of intercellular communication (e.g., gap-junction communication, secreted factors) between the two heterotypic cell populations, which permits investigation of the mechanisms underlying the development of the TME, including the role of gap-junction permeability. This model serves as a valuable tool in enhancing our understanding of the initial events leading to cancer-stroma initiation, the early evolution of the TME, and the modulating effect of the stroma on the responses of cancer cells to therapeutic agents.

Cross-species Conservation of context-specific networks

BACKGROUND

Many large data compendia on context-specific high-throughput genomic and regulatory data have been made available by international research consortia such as ENCODE, TCGA, and Epigenomics Roadmap. The use of these resources is impaired by the sheer size of the available big data and big metadata. Many of these context-specific data can be modeled as data derived regulatory networks (DDRNs) representing the complex and complicated interactions between transcription factors and target genes. These DDRNs are useful for the understanding of regulatory mechanisms and helpful for interpreting biomedical data.

RESULTS

The Cross-species Conservation framework (CroCo) provides a network-oriented view on the ENCODE regulatory data (CroCo network repository), convenient ways to access and browse networks and metadata, and a method to combine networks across compendia, experimental techniques, and species (CroCo tool suite). DDRNs can be combined with additional information and networks derived from the literature, curated resources, and computational predictions in order to enable detailed exploration and cross checking of regulatory interactions. Applications of the CroCo framework range from simple evidence look-up for user-defined regulatory interactions to the identification of conserved sub-networks in diverse cell-lines, conditions, and even species.

CONCLUSION

CroCo adds an intuitive unifying view on the data from the ENCODE projects via a comprehensive repository of derived context-specific regulatory networks and enables flexible cross-context, cross-species, and cross-compendia comparison via a basis set of analysis tools. The CroCo web-application and Cytoscape plug-in are freely available at: http://services.bio.ifi.lmu.de/croco-web . The web-page links to a detailed system description, a user guide, and tutorial videos presenting common use cases of the CroCo framework.

Gene activation by metazoan enhancers: Diverse mechanisms stimulate distinct steps of transcription

Enhancers can stimulate transcription by a number of different mechanisms which control different stages of the transcription cycle of their target genes, from recruitment of the transcription machinery to elongation by RNA polymerase. These mechanisms may not be mutually exclusive, as a single enhancer may act through different pathways by binding multiple transcription factors. Multiple enhancers may also work together to regulate transcription of a shared target gene. Most of the evidence supporting different enhancer mechanisms comes from the study of single genes, but new high-throughput experimental frameworks offer the opportunity to integrate and generalize disparate mechanisms identified at single genes. This effort is especially important if we are to fully understand how sequence variation within enhancers contributes to human disease.

Differential regulation of mouse and human nephron progenitors by the Six family of transcriptional regulators

Nephron endowment is determined by the self-renewal and induction of a nephron progenitor pool established at the onset of kidney development. In the mouse, the related transcriptional regulators Six1 and Six2 play non-overlapping roles in nephron progenitors. Transient Six1 activity prefigures, and is essential for, active nephrogenesis. By contrast, Six2 maintains later progenitor self-renewal from the onset of nephrogenesis. We compared the regulatory actions of Six2 in mouse and human nephron progenitors by chromatin immunoprecipitation followed by DNA sequencing (ChIP-seq). Surprisingly, SIX1 was identified as a SIX2 target unique to the human nephron progenitors. Furthermore, RNA-seq and immunostaining revealed overlapping SIX1 and SIX2 activity in 16 week human fetal nephron progenitors. Comparative bioinformatic analysis of human SIX1 and SIX2 ChIP-seq showed each factor targeted a similar set of cis-regulatory modules binding an identical target recognition motif. In contrast to the mouse where Six2 binds its own enhancers but does not interact with DNA around Six1, both human SIX1 and SIX2 bind homologous SIX2 enhancers and putative enhancers positioned around SIX1. Transgenic analysis of a putative human SIX1 enhancer in the mouse revealed a transient, mouse-like, pre-nephrogenic, Six1 regulatory pattern. Together, these data demonstrate a divergence in SIX-factor regulation between mouse and human nephron progenitors. In the human, an auto/cross-regulatory loop drives continued SIX1 and SIX2 expression during active nephrogenesis. By contrast, the mouse establishes only an auto-regulatory Six2 loop. These data suggest differential SIX-factor regulation might have contributed to species differences in nephron progenitor programs such as the duration of nephrogenesis and the final nephron count.

Generation and Application of Mouse-Rat Allodiploid Embryonic Stem Cells

Mammalian interspecific hybrids provide unique advantages for mechanistic studies of speciation, gene expression regulation, and X chromosome inactivation (XCI) but are constrained by their limited natural resources. Previous artificially generated mammalian interspecific hybrid cells are usually tetraploids with unstable genomes and limited developmental abilities. Here, we report the generation of mouse-rat allodiploid embryonic stem cells (AdESCs) by fusing haploid ESCs of the two species. The AdESCs have a stable allodiploid genome and are capable of differentiating into all three germ layers and early-stage germ cells. Both the mouse and rat alleles have comparable contributions to the expression of most genes. We have proven AdESCs as a powerful tool to study the mechanisms regulating X chromosome inactivation and to identify X inactivation-escaping genes, as well as to efficiently identify genes regulating phenotypic differences between species. A similar method could be used to create hybrid AdESCs of other distantly related species.

Engineered Swine Models of Cancer

Over the past decade, the technology to engineer genetically modified swine has seen many advancements, and because their physiology is remarkably similar to that of humans, swine models of cancer may be extremely valuable for preclinical safety studies as well as toxicity testing of pharmaceuticals prior to the start of human clinical trials. Hence, the benefits of using swine as a large animal model in cancer research and the potential applications and future opportunities of utilizing pigs in cancer modeling are immense. In this review, we discuss how pigs have been and can be used as a biomedical models for cancer research, with an emphasis on current technologies. We have focused on applications of precision genetics that can provide models that mimic human cancer predisposition syndromes. In particular, we describe the advantages of targeted gene-editing using custom endonucleases, specifically TALENs and CRISPRs, and transposon systems, to make novel pig models of cancer with broad preclinical applications.

Evolution of Epigenetic Regulation in Vertebrate Genomes

Empirical models of sequence evolution have spurred progress in the field of evolutionary genetics for decades. We are now realizing the importance and complexity of the eukaryotic epigenome. While epigenome analysis has been applied to genomes from single-cell eukaryotes to human, comparative analyses are still relatively few and computational algorithms to quantify epigenome evolution remain scarce. Accordingly, a quantitative model of epigenome evolution remains to be established. We review here the comparative epigenomics literature and synthesize its overarching themes. We also suggest one mechanism, transcription factor binding site (TFBS) turnover, which relates sequence evolution to epigenetic conservation or divergence. Lastly, we propose a framework for how the field can move forward to build a coherent quantitative model of epigenome evolution.

A Hox Transcription Factor Collective Binds a Highly Conserved Distal-less cis-Regulatory Module to Generate Robust Transcriptional Outcomes

cis-regulatory modules (CRMs) generate precise expression patterns by integrating numerous transcription factors (TFs). Surprisingly, CRMs that control essential gene patterns can differ greatly in conservation, suggesting distinct constraints on TF binding sites. Here, we show that a highly conserved Distal-less regulatory element (DCRE) that controls gene expression in leg precursor cells recruits multiple Hox, Extradenticle (Exd) and Homothorax (Hth) complexes to mediate dual outputs: thoracic activation and abdominal repression. Using reporter assays, we found that abdominal repression is particularly robust, as neither individual binding site mutations nor a DNA binding deficient Hth protein abolished cooperative DNA binding and in vivo repression. Moreover, a re-engineered DCRE containing a distinct configuration of Hox, Exd, and Hth sites also mediated abdominal Hox repression. However, the re-engineered DCRE failed to perform additional segment-specific functions such as thoracic activation. These findings are consistent with two emerging concepts in gene regulation: First, the abdominal Hox/Exd/Hth factors utilize protein-protein and protein-DNA interactions to form repression complexes on flexible combinations of sites, consistent with the TF collective model of CRM organization. Second, the conserved DCRE mediates multiple cell-type specific outputs, consistent with recent findings that pleiotropic CRMs are associated with conserved TF binding and added evolutionary constraints.

Enhancers are regulatory elements that interact with transcription factor proteins to control cell-specific gene expression during development. Surprisingly, only a subset of enhancers are highly conserved at the sequence level, even though the expression patterns they control are often conserved and essential for proper development. Why some enhancer sequences are highly conserved whereas others are not is not well understood. In this study, we characterize a highly conserved enhancer that regulates gene expression in leg precursor cells. We find that this enhancer has dual regulatory activities that include gene activation in thoracic segments and gene repression in abdominal segments. Surprisingly, we show that the conserved enhancer can tolerate numerous sequence changes yet mediate robust transcription factor binding and abdominal repression. These findings are consistent with abdominal transcription factors binding numerous different configurations of binding sites. So, why is this enhancer highly conserved? We found that overlapping sequences within the enhancer also contribute to thoracic activation, suggesting the enhancer sequences are under added functional constraints. Altogether, our results provide new insights into why some enhancers are highly conserved at the sequence level while others can tolerate sequence changes.

Insight into GATA1 transcriptional activity through interrogation of cis elements disrupted in human erythroid disorders

Whole-exome sequencing has been incredibly successful in identifying causal genetic variants and has revealed a number of novel genes associated with blood and other diseases. One limitation of this approach is that it overlooks mutations in noncoding regulatory elements. Furthermore, the mechanisms by which mutations in transcriptionalcis-regulatory elements result in disease remain poorly understood. Here we used CRISPR/Cas9 genome editing to interrogate three such elements harboring mutations in human erythroid disorders, which in all cases are predicted to disrupt a canonical binding motif for the hematopoietic transcription factor GATA1. Deletions of as few as two to four nucleotides resulted in a substantial decrease (>80%) in target gene expression. Isolated deletions of the canonical GATA1 binding motif completely abrogated binding of the cofactor TAL1, which binds to a separate motif. Having verified the functionality of these three GATA1 motifs, we demonstrate strong evolutionary conservation of GATA1 motifs in regulatory elements proximal to other genes implicated in erythroid disorders, and show that targeted disruption of such elements results in altered gene expression. By modeling transcription factor binding patterns, we show that multiple transcription factors are associated with erythroid gene expression, and have created predictive maps modeling putative disruptions of their binding sites at key regulatory elements. Our study provides insight into GATA1 transcriptional activity and may prove a useful resource for investigating the pathogenicity of noncoding variants in human erythroid disorders.

Single-CpG resolution mapping of 5-hydroxymethylcytosine by chemical labeling and exonuclease digestion identifies evolutionarily unconserved CpGs as TET targets

Conventional techniques for single-base resolution mapping of epigenetic modifications of DNA such as 5-hydroxymethylcytosine (5hmC) rely on the sequencing of bisulfite-modified DNA. Here we present an alternative approach called SCL-exo which combines selective chemical labeling (SCL) of 5hmC in genomic DNA with exonuclease (exo) digestion of the bead-trapped modified DNA molecules. Associated with a straightforward bioinformatic analysis, this new procedure provides an unbiased and fast method for mapping this epigenetic mark at high resolution. Implemented on mouse genomic DNA from in vitro-differentiated neural precursor cells, SCL-exo sheds light on an intrinsic lack of conservation of hydroxymethylated CpGs across vertebrates.

Differentiation of the SH-SY5Y Human Neuroblastoma Cell Line

Having appropriate in vivo and in vitro systems that provide translational models for human disease is an integral aspect of research in neurobiology and the neurosciences. Traditional in vitro experimental models used in neurobiology include primary neuronal cultures from rats and mice, neuroblastoma cell lines including rat B35 and mouse Neuro-2A cells, rat PC12 cells, and short-term slice cultures. While many researchers rely on these models, they lack a human component and observed experimental effects could be exclusive to the respective species and may not occur identically in humans. Additionally, although these cells are neurons, they may have unstable karyotypes, making their use problematic for studies of gene expression and reproducible studies of cell signaling. It is therefore important to develop more consistent models of human neurological disease. The following procedure describes an easy-to-follow, reproducible method to obtain homogenous and viable human neuronal cultures, by differentiating the chromosomally stable human neuroblastoma cell line, SH-SY5Y. This method integrates several previously described methods(1-4) and is based on sequential removal of serum from media. The timeline includes gradual serum-starvation, with introduction of extracellular matrix proteins and neurotrophic factors. This allows neurons to differentiate, while epithelial cells are selected against, resulting in a homogeneous neuronal culture. Representative results demonstrate the successful differentiation of SH-SY5Y neuroblastoma cells from an initial epithelial-like cell phenotype into a more expansive and branched neuronal phenotype. This protocol offers a reliable way to generate homogeneous populations of neuronal cultures that can be used for subsequent biochemical and molecular analyses, which provides researchers with a more accurate translational model of human infection and disease.

Human and murine erythropoiesis

PURPOSE OF REVIEW

Research into the fundamental mechanisms of erythropoiesis has provided critical insights into inherited and acquired disorders of the erythrocyte. Studies of human erythropoiesis have primarily utilized in-vitro systems, whereas murine models have provided insights from in-vivo studies. This report reviews recent insights into human and murine erythropoiesis gained from transcriptome-based analyses.

RECENT FINDINGS

The availability of high-throughput genomic methodologies has allowed attainment of detailed gene expression data from cells at varying developmental and differentiation stages of erythropoiesis. Transcriptome analyses of human and murine reveal both stage and species-specific similarities and differences across terminal erythroid differentiation. Erythroid-specific long noncoding RNAs exhibit poor sequence conservation between human and mouse. Genome-wide analyses of alternative splicing reveal that complex, dynamic, stage-specific programs of alternative splicing program are utilized during terminal erythroid differentiation. Transcriptome data provide a significant resource for understanding mechanisms of normal and perturbed erythropoiesis. Understanding these processes will provide innovative strategies to detect, diagnose, prevent, and treat hematologic disease.

SUMMARY

Understanding the shared and different mechanisms controlling human and murine erythropoiesis will allow investigators to leverage the best model system to provide insights in normal and perturbed erythropoiesis.

The Hematopoietic Stem and Progenitor Cell Cistrome: GATA Factor-Dependent cis-Regulatory Mechanisms

Transcriptional regulators mediate the genesis and function of the hematopoietic system by binding complex ensembles of cis-regulatory elements to establish genetic networks. While thousands to millions of any given cis-element resides in a genome, how transcriptional regulators select these sites and how site attributes dictate functional output is not well understood. An instructive system to address this problem involves the GATA family of transcription factors that control vital developmental and physiological processes and are linked to multiple human pathologies. Although GATA factors bind DNA motifs harboring the sequence GATA, only a very small subset of these abundant motifs are occupied in genomes. Mechanistic studies revealed a unique configuration of a GATA factor-regulated cis-element consisting of an E-box and a downstream GATA motif separated by a short DNA spacer. GATA-1- or GATA-2-containing multiprotein complexes at these composite elements control transcription of genes critical for hematopoietic stem cell emergence in the mammalian embryo, hematopoietic progenitor cell regulation, and erythroid cell maturation. Other constituents of the complex include the basic helix-loop-loop transcription factor Scl/TAL1, its heterodimeric partner E2A, and the Lim domain proteins LMO2 and LDB1. This chapter reviews the structure/function of E-box-GATA composite cis-elements, which collectively constitute an important sector of the hematopoietic stem and progenitor cell cistrome.

Implications of Epigenetic Variability within a Cell Population for "Cell Type" Classification

Here, we propose a new approach to defining nerve “cell types” in reaction to recent advances in single cell analysis. Among cells previously thought to be equivalent, considerable differences in global gene expression and biased tendencies among differing developmental fates have been demonstrated within multiple lineages. The model of classifying cells into distinct types thus has to be revised to account for this intrinsic variability. A “cell type” could be a group of cells that possess similar, but not necessarily identical properties, variable within a spectrum of epigenetic adjustments that permit its developmental path toward a specific function to be achieved. Thus, the definition of a cell type is becoming more similar to the definition of a species: sharing essential properties with other members of its group, but permitting a certain amount of deviation in aspects that do not seriously impact function. This approach accommodates, even embraces the spectrum of natural variation found in various cell populations and consequently avoids the fallacy of false equivalence. For example, developing neurons will react to their microenvironments with epigenetic changes resulting in slight changes in gene expression and morphology. Addressing the new questions implied here will have significant implications for developmental neurobiology.

Reversion to an embryonic alternative splicing program enhances leukemia stem cell self-renewal

Formative research suggests that a human embryonic stem cell-specific alternative splicing gene regulatory network, which is repressed by Muscleblind-like (MBNL) RNA binding proteins, is involved in cell reprogramming. In this study, RNA sequencing, splice isoform-specific quantitative RT-PCR, lentiviral transduction, and in vivo humanized mouse model studies demonstrated that malignant reprogramming of progenitors into self-renewing blast crisis chronic myeloid leukemia stem cells (BC LSCs) was partially driven by decreased MBNL3. Lentiviral knockdown of MBNL3 resulted in reversion to an embryonic alternative splice isoform program typified by overexpression of CD44 transcript variant 3, containing variant exons 8-10, and BC LSC proliferation. Although isoform-specific lentiviral CD44v3 overexpression enhanced chronic phase chronic myeloid leukemia (CML) progenitor replating capacity, lentiviral shRNA knockdown abrogated these effects. Combined treatment with a humanized pan-CD44 monoclonal antibody and a breakpoint cluster region - ABL proto-oncogene 1, nonreceptor tyrosine kinase (BCR-ABL1) antagonist inhibited LSC maintenance in a niche-dependent manner. In summary, MBNL3 down-regulation-related reversion to an embryonic alternative splicing program, typified by CD44v3 overexpression, represents a previously unidentified mechanism governing malignant progenitor reprogramming in malignant microenvironments and provides a pivotal opportunity for selective BC LSC detection and therapeutic elimination.

Mapping the effects of drugs on the immune system

Understanding how drugs affect the immune system has consequences for treating disease and minimizing unwanted side effects. Here we present an integrative computational approach for predicting interactions between drugs and immune cells in a system-wide manner. The approach matches gene sets between transcriptional signatures to determine their similarity. We apply the method to model the interactions between 1,309 drugs and 221 immune cell types and predict 69,995 known and novel interactions. The resulting immune-cell pharmacology map is used to predict how 5 drugs influence 4 immune cell types in humans and mice. To validate the predictions, we analyzed patient records and examined cell population changes from in vivo experiments. Our method offers a tool for screening thousands of interactions to identify relationships between drugs and the immune system.

Mechanisms Underlying the Selection and Function of Macrophage-Specific Enhancers

Macrophages populate every tissue of the body and play vital roles in homeostasis, pathogen elimination, and tissue healing. These cells possess the ability to adapt to a multitude of abruptly changing and complex environments. Furthermore, different populations of resident tissue macrophages each show their own defining gene signatures. The enhancer repertoire of these cells underlies both the cellular identity of a given subset of resident macrophage population and their ability to dynamically alter, in an efficient manner, their gene expression programs in response to internal and external signals. Notably, transcription is pervasive at active enhancers and enhancer RNAs, or eRNAs, are tightly correlated to regulated transcription of protein-coding genes. Furthermore, selection and establishment of enhancers is a dynamic and plastic process in which activation of intracellular signaling pathways by factors present in a macrophage's environment play a determining role. Here, we review recent studies providing insights into the distinct mechanisms that contribute to the selection and function of enhancers in macrophages and the relevance of studying these mechanisms to gain a better understanding of complex human diseases.

NeuroD1 reprograms chromatin and transcription factor landscapes to induce the neuronal program

Cell fate specification relies on the action of critical transcription factors that become available at distinct stages of embryonic development. One such factor is NeuroD1, which is essential for eliciting the neuronal development program and possesses the ability to reprogram other cell types into neurons. Given this capacity, it is important to understand its targets and the mechanism underlying neuronal specification. Here, we show that NeuroD1 directly binds regulatory elements of neuronal genes that are developmentally silenced by epigenetic mechanisms. This targeting is sufficient to initiate events that confer transcriptional competence, including reprogramming of transcription factor landscape, conversion of heterochromatin to euchromatin, and increased chromatin accessibility, indicating potential pioneer factor ability of NeuroD1. The transcriptional induction of neuronal fate genes is maintained via epigenetic memory despite a transient NeuroD1 induction during neurogenesis. NeuroD1 also induces genes involved in the epithelial-to-mesenchymal transition, thereby promoting neuronal migration. Our study not only reveals the NeuroD1-dependent gene regulatory program driving neurogenesis but also increases our understanding of how cell fate specification during development involves a concerted action of transcription factors and epigenetic mechanisms.

Genetic and epigenetic variation in the lineage specification of regulatory T cells

Regulatory T (Treg) cells, which suppress autoimmunity and other inflammatory states, are characterized by a distinct set of genetic elements controlling their gene expression. However, the extent of genetic and associated epigenetic variation in the Treg cell lineage and its possible relation to disease states in humans remain unknown. We explored evolutionary conservation of regulatory elements and natural human inter-individual epigenetic variation in Treg cells to identify the core transcriptional control program of lineage specification. Analysis of single nucleotide polymorphisms in core lineage-specific enhancers revealed disease associations, which were further corroborated by high-resolution genotyping to fine map causal polymorphisms in lineage-specific enhancers. Our findings suggest that a small set of regulatory elements specify the Treg lineage and that genetic variation in Treg cell-specific enhancers may alter Treg cell function contributing to polygenic disease.

DOI:http://dx.doi.org/10.7554/eLife.07571.001

The immune system protects the body from infection. Key to this protection is the ability to mount an immune response that is sufficient to deal with the threat, but is not so large that the damage it causes to the body exceeds its immediate benefit. Immune cells called regulatory T cells (or Treg cells for short) help to shut down the immune response after a threat has been successfully destroyed. They also prevent the immune system from attacking the body's own cells, a phenomenon known as autoimmunity.

All cells in the body carry the same set of genes, but the activity of these genes varies between cell types to enable the cells to perform their different jobs. This is possible because our DNA contains regions called regulatory elements that can control the expression of particular genes. These regions can be activated in specific types of cells, which results in specific chemical modifications to DNA that only affect gene activity in those cells.

The DNA sequences of these regulatory elements vary between individuals. This ‘genetic variation’ can lead to differences in the chemical modifications that occur to DNA, which is known as epigenetic variation. This means that Treg cells in one person may work in a different way to those in another individual, which could make some individuals more susceptible to autoimmune diseases than others. However, it was not clear how much genetic and epigenetic variation exists in Treg cells.

Here, Arvey et al. examined Treg and other immune cells from human and mouse donors. The experiments show that some of the epigenetic modifications present in many individuals only in Treg cells, which indicates that they may be important for the activity of the Treg cells. Unexpectedly, most of the epigenetic modifications were specific to either mice or humans, but Arvey et al. identified a core set of genes that had been modified only in Treg cells in both species. In the human cells, Arvey et al. also identified genetic differences in regulatory elements that are associated with autoimmune diseases.

Arvey et al.'s findings suggest that a small set of regulatory elements are crucial to the role of Treg cells in the immune system. Furthermore, genetic variation in these elements can lead to epigenetic changes in Treg cells that contribute to autoimmune diseases in humans. Further study may lead to the development of new treatments for these diseases.

DOI:http://dx.doi.org/10.7554/eLife.07571.002

Genetic Dissection of the Host Tropism of Human-Tropic Pathogens

Infectious diseases are the second leading cause of death worldwide. Although the host multitropism of some pathogens has rendered their manipulation possible in animal models, the human-restricted tropism of numerous viruses, bacteria, fungi, and parasites has seriously hampered our understanding of these pathogens. Hence, uncovering the genetic basis underlying the narrow tropism of such pathogens is critical for understanding their mechanisms of infection and pathogenesis. Moreover, such genetic dissection is essential for the generation of permissive animal models that can serve as critical tools for the development of therapeutics or vaccines against challenging human pathogens. In this review, we describe different experimental approaches utilized to uncover the genetic foundation regulating pathogen host tropism as well as their relevance for studying the tropism of several important human pathogens. Finally, we discuss the current and future uses of this knowledge for generating genetically modified animal models permissive for these pathogens.

Exploiting genomics and natural genetic variation to decode macrophage enhancers

The mammalian genome contains on the order of a million enhancer-like regions that are required to establish the identities and functions of specific cell types. Here, we review recent studies in immune cells that have provided insight into the mechanisms that selectively activate certain enhancers in response to cell lineage and environmental signals. We describe a working model wherein distinct classes of transcription factors define the repertoire of active enhancers in macrophages through collaborative and hierarchical interactions, and discuss important challenges to this model, specifically providing examples from T cells. We conclude by discussing the use of natural genetic variation as a powerful approach for decoding transcription factor combinations that play dominant roles in establishing the enhancer landscapes, and the potential that these insights have for advancing our understanding of the molecular causes of human disease.

Transcriptional and epigenetic mechanisms of early cortical development: An examination of how Pax6 coordinates cortical development

The development of the cortex is an elaborate process that integrates a plethora of finely tuned molecular processes ranging from carefully regulated gradients of transcription factors, dynamic changes in the chromatin landscape, or formation of protein complexes to elicit and regulate transcription. Combined with cellular processes such as cell type specification, proliferation, differentiation, and migration, all of these developmental processes result in the establishment of an adult mammalian cortex with its typical lamination and regional patterning. By examining in-depth the role of one transcription factor, Pax6, on the regulation of cortical development, its integration in the regulation of chromatin state, and its regulation by cis-regulatory elements, we aim to demonstrate the importance of integrating each level of regulation in our understanding of cortical development.

The frequent evolutionary birth and death of functional promoters in mouse and human

Promoters are central to the regulation of gene expression. Changes in gene regulation are thought to underlie much of the adaptive diversification between species and phenotypic variation within populations. In contrast to earlier work emphasizing the importance of enhancer evolution and subtle sequence changes at promoters, we show that dramatic changes such as the complete gain and loss (collectively, turnover) of functional promoters are common. Using quantitative measures of transcription initiation in both humans and mice across 52 matched tissues, we discriminate promoter sequence gains from losses and resolve the lineage of changes. We also identify expression divergence and functional turnover between orthologous promoters, finding only the latter is associated with local sequence changes. Promoter turnover has occurred at the majority (>56%) of protein-coding genes since humans and mice diverged. Tissue-restricted promoters are the most evolutionarily volatile where retrotransposition is an important, but not the sole, source of innovation. There is considerable heterogeneity of turnover rates between promoters in different tissues, but the consistency of these in both lineages suggests that the same biological systems are similarly inclined to transcriptional rewiring. The genes affected by promoter turnover show evidence of adaptive evolution. In mice, promoters are primarily lost through deletion of the promoter containing sequence, whereas in humans, many promoters appear to be gradually decaying with weak transcriptional output and relaxed selective constraint. Our results suggest that promoter gain and loss is an important process in the evolutionary rewiring of gene regulation and may be a significant source of phenotypic diversification.

Significant expansion of the REST/NRSF cistrome in human versus mouse embryonic stem cells: potential implications for neural development

Recent studies have employed cross-species comparisons of transcription factor binding, reporting significant regulatory network 'rewiring' between species. Here, we address how a transcriptional repressor targets and regulates neural genes differentially between human and mouse embryonic stem cells (ESCs). We find that the transcription factor, Repressor Element 1 Silencing Transcription factor (REST; also called neuron restrictive silencer factor) binds to a core group of ∼1200 syntenic genomic regions in both species, with these conserved sites highly enriched with co-factors, selective histone modifications and DNA hypomethylation. Genes with conserved REST binding are enriched with neural functions and more likely to be upregulated upon REST depletion. Interestingly, we identified twice as many REST peaks in human ESCs compared to mouse ESCs. Human REST cistrome expansion involves additional peaks in genes targeted by REST in both species and human-specific gene targets. Genes with expanded REST occupancy in humans are enriched for learning or memory functions. Analysis of neurological disorder associated genes reveals that Amyotrophic Lateral Sclerosis and oxidative stress genes are particularly enriched with human-specific REST binding. Overall, our results demonstrate that there is substantial rewiring of human and mouse REST cistromes, and that REST may have human-specific roles in brain development and functions.

Transcriptional Control of Stem and Progenitor Potential

Hematopoiesis is characterized by a lifelong balance between hematopoietic stem cell (HSC) self-renewal and differentiation into mature blood populations. Proper instruction of cell fate decisions requires tight homeostatic regulation of transcriptional programs through a combination of epigenetic modifications, management of cis-regulatory elements, and transcription factor activity. Recent work has focused on integrating biochemical, genetic, and evolutionary data sets to gain further insight into these regulatory components. Long noncoding RNA (lncRNA), post-translational modifications of transcription factors, and circadian rhythm add additional layers of complexity. These analyses have provided a wealth of information, much of which has been made available through public databases. Elucidating the regulatory processes that govern hematopoietic transcriptional programs is expected to provide useful insights into hematopoiesis that may be applied broadly across tissue types while enabling the discovery and implementation of therapeutics to treat human disease.

Hematopoietic Signaling Mechanism Revealed from a Stem/Progenitor Cell Cistrome

Thousands of cis-elements in genomes are predicted to have vital functions. Although conservation, activity in surrogate assays, polymorphisms, and disease mutations provide functional clues, deletion from endogenous loci constitutes the gold-standard test. A GATA-2-binding, Gata2 intronic cis-element (+9.5) required for hematopoietic stem cell genesis in mice is mutated in a human immunodeficiency syndrome. Because +9.5 is the only cis-element known to mediate stem cell genesis, we devised a strategy to identify functionally comparable enhancers ("+9.5-like") genome-wide. Gene editing revealed +9.5-like activity to mediate GATA-2 occupancy, chromatin opening, and transcriptional activation. A +9.5-like element resided in Samd14, which encodes a protein of unknown function. Samd14 increased hematopoietic progenitor levels/activity and promoted signaling by a pathway vital for hematopoietic stem/progenitor cell regulation (stem cell factor/c-Kit), and c-Kit rescued Samd14 loss-of-function phenotypes. Thus, the hematopoietic stem/progenitor cell cistrome revealed a mediator of a signaling pathway that has broad importance for stem/progenitor cell biology.

Drug discovery in prostate cancer mouse models

INTRODUCTION

The mouse is an important, though imperfect, organism with which to model human disease and to discover and test novel drugs in a preclinical setting. Many experimental strategies have been used to discover new biological and molecular targets in the mouse, with the hopes of translating these discoveries into novel drugs to treat prostate cancer in humans. Modeling prostate cancer in the mouse, however, has been challenging, and often drugs that work in mice have failed in human trials.

AREAS COVERED

The authors discuss the similarities and differences between mice and men; the types of mouse models that exist to model prostate cancer; practical questions one must ask when using a mouse as a model; and potential reasons that drugs do not often translate to humans. They also discuss the current value in using mouse models for drug discovery to treat prostate cancer and what needs are still unmet in field.

EXPERT OPINION

With proper planning and following practical guidelines by the researcher, the mouse is a powerful experimental tool. The field lacks genetically engineered metastatic models, and xenograft models do not allow for the study of the immune system during the metastatic process. There remain several important limitations to discovering and testing novel drugs in mice for eventual human use, but these can often be overcome. Overall, mouse modeling is an essential part of prostate cancer research and drug discovery. Emerging technologies and better and ever-increasing forms of communication are moving the field in a hopeful direction.

Mutations in CHD2 cause defective association with active chromatin in chronic lymphocytic leukemia

Great progress has recently been achieved in the understanding of the genomic alterations driving chronic lymphocytic leukemia (CLL). Nevertheless, the specific molecular mechanisms governing chromatin remodeling in CLL are unknown. Here we report the genetic and functional characterization of somatic mutations affecting the chromatin remodeler CHD2, one of the most frequently mutated genes in CLL (5.3%) and in monoclonal B lymphocytosis (MBL, 7%), a B-cell expansion that can evolve to CLL. Most of the mutations affecting CHD2, identified by whole-exome sequencing of 456 CLL and 43 MBL patients, are either truncating or affect conserved residues in functional domains, thus supporting a putative role for CHD2 as a tumor suppressor gene. CHD2 mutants show altered nuclear distribution, and a chromodomain helicase DNA binding protein 2 (CHD2) mutant affected in its DNA-binding domain exhibits defective association with active chromatin. Clinicobiological analyses show that most CLL patients carrying CHD2 mutations also present mutated immunoglobulin heavy chain variable region genes (IGHVs), being the most frequently mutated gene in this prognostic subgroup. This is the first study providing functional evidence supporting CHD2 as a cancer driver and opens the way to further studies of the role of this chromatin remodeler in CLL.

An ancient protein-DNA interaction underlying metazoan sex determination

DMRT transcription factors are deeply conserved regulators of metazoan sexual development. They share the DM DNA-binding domain, a unique intertwined double zinc-binding module followed by a C-terminal recognition helix, which binds a pseudopalindromic target DNA. Here we show that DMRT proteins use a unique binding interaction, inserting two adjacent antiparallel recognition helices into a widened DNA major groove to make base-specific contacts. Versatility in how specific base contacts are made allows human DMRT1 to use multiple DNA binding modes (tetramer, trimer and dimer). Chromatin immunoprecipitation with exonuclease treatment (ChIP-exo) indicates that multiple DNA binding modes also are used in vivo. We show that mutations affecting residues crucial for DNA recognition are associated with an intersex phenotype in flies and with male-to-female sex reversal in humans. Our results illuminate an ancient molecular interaction underlying much of metazoan sexual development.

High-throughput sequencing technologies

The human genome sequence has profoundly altered our understanding of biology, human diversity, and disease. The path from the first draft sequence to our nascent era of personal genomes and genomic medicine has been made possible only because of the extraordinary advancements in DNA sequencing technologies over the past 10 years. Here, we discuss commonly used high-throughput sequencing platforms, the growing array of sequencing assays developed around them, as well as the challenges facing current sequencing platforms and their clinical application.

Occupancy by key transcription factors is a more accurate predictor of enhancer activity than histone modifications or chromatin accessibility

Background

Regulated gene expression controls organismal development, and variation in regulatory patterns has been implicated in complex traits. Thus accurate prediction of enhancers is important for further understanding of these processes. Genome-wide measurement of epigenetic features, such as histone modifications and occupancy by transcription factors, is improving enhancer predictions, but the contribution of these features to prediction accuracy is not known. Given the importance of the hematopoietic transcription factor TAL1 for erythroid gene activation, we predicted candidate enhancers based on genomic occupancy by TAL1 and measured their activity. Contributions of multiple features to enhancer prediction were evaluated based on the results of these and other studies.

Results

TAL1-bound DNA segments were active enhancers at a high rate both in transient transfections of cultured cells (39 of 79, or 56%) and transgenic mice (43 of 66, or 65%). The level of binding signal for TAL1 or GATA1 did not help distinguish TAL1-bound DNA segments as active versus inactive enhancers, nor did the density of regulation-related histone modifications. A meta-analysis of results from this and other studies (273 tested predicted enhancers) showed that the presence of TAL1, GATA1, EP300, SMAD1, H3K4 methylation, H3K27ac, and CAGE tags at DNase hypersensitive sites gave the most accurate predictors of enhancer activity, with a success rate over 80% and a median threefold increase in activity. Chromatin accessibility assays and the histone modifications H3K4me1 and H3K27ac were sensitive for finding enhancers, but they have high false positive rates unless transcription factor occupancy is also included.

Conclusions

Occupancy by key transcription factors such as TAL1, GATA1, SMAD1, and EP300, along with evidence of transcription, improves the accuracy of enhancer predictions based on epigenetic features.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-015-0009-5) contains supplementary material, which is available to authorized users.

Selective Distal Enhancer Control of the Mmp13 Gene Identified through Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR) Genomic Deletions

Matrix metalloproteinase 13 (Mmp13, collagenase-3) plays an essential role in bone metabolism and mineral homeostasis. It is regulated by numerous factors, including BMP-2, parathyroid hormone, and 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3), through transcription factors such as Runt-related transcription factor 2 (RUNX2), CCAAT/enhancer-binding protein β (C/EBPβ), OSX, and vitamin D receptor (VDR). During osteoblast maturation, the basal expression of Mmp13 and its sensitivity to 1,25(OH)2D3 are strikingly increased. In this report, ChIP-sequencing analysis in mouse preosteoblasts revealed that the Mmp13 gene was probably regulated by three major enhancers located -10, -20, and -30 kb upstream of the gene promoter, occupied by activated VDR and prebound C/EBPβ and RUNX2, respectively. Initially, bacterial artificial chromosome clone recombineering and traditional mutagenesis defined binding sites for VDR and RUNX2. We then employed a CRISPR/Cas9 gene editing approach to delete the -10 and -30 kb Mmp13 enhancers, a region proximal to the promoter, and VDR or RUNX2. VDR-mediated up-regulation of Mmp13 transcription was completely abrogated upon removal of the -10 kb enhancer, resulting in a 1,25(OH)2D3-directed repression of Mmp13. Deletion of either the -30 kb enhancer or RUNX2 resulted in a complete loss of basal transcript activity and a ChIP-identified destabilization of the chromatin enhancer environment and factor binding. Whereas enhancer deletions only affected Mmp13 expression, the RUNX2 deletion led to changes in gene expression, a reduction in cellular proliferation, and an inability to differentiate. We conclude that the Mmp13 gene is regulated via at least three specific distal enhancers that display independent activities yet are able to integrate response from multiple signaling pathways in a model of activation and suppression.

The Human GATA1 Gene Retains a 5' Insulator That Maintains Chromosomal Architecture and GATA1 Expression Levels in Splenic Erythroblasts

GATA1 is a key transcription factor for erythropoiesis. GATA1 gene expression is strictly regulated at the transcriptional level. While the regulatory mechanisms governing mouse Gata1 (mGata1) gene expression have been studied extensively, how expression of the human GATA1 (hGATA1) gene is regulated remains to be elucidated. To address this issue, we generated hGATA1 bacterial artificial chromosome (BAC) transgenic mouse lines harboring a 183-kb hGATA1 locus covering the hGATA1 exons and distal flanking sequences. Transgenic hGATA1 expression coincides with endogenous mGata1 expression and fully rescues hematopoietic deficiency in mGata1 knockdown mice. The transgene exhibited copy number-dependent and integration position-independent expression of hGATA1, indicating the presence of chromatin insulator activity within the transgene. We found a novel insulator element at 29 kb 5' to the hGATA1 gene and refer to this element as the 5' CCCTC-binding factor (CTCF) site. Substitution mutation of the 5' CTCF site in the hGATA1 BAC disrupted the chromatin architecture and led to a reduction of hGATA1 expression in splenic erythroblasts under conditions of stress erythropoiesis. Our results demonstrate that expression of the hGATA1 gene is regulated through the chromatin architecture organized by 5' CTCF site-mediated intrachromosomal interactions in the hGATA1 locus.

Conserved epigenomic signals in mice and humans reveal immune basis of Alzheimer's disease

Alzheimer’s disease (AD) is a severe¹ age-related neurodegenerative disorder characterized by accumulation of amyloid-β (Aβ) plaques and neurofibrillary tangles, synaptic and neuronal loss, and cognitive decline. Several genes have been implicated in AD, but chromatin state alterations during neurodegeneration remain uncharacterized. Here, we profile transcriptional and chromatin state dynamics across early and late pathology in the hippocampus of an inducible mouse model of AD-like neurodegeneration. We find a coordinated downregulation of synaptic plasticity genes and regulatory regions, and upregulation of immune response genes and regulatory regions, which are targeted by factors that belong to the ETS family of transcriptional regulators, including PU.1. Human regions orthologous to increasing-level enhancers show immune cell-specific enhancer signatures as well as immune cell expression quantitative trait loci (eQTL), while decreasing-level enhancer orthologs show fetal-brain-specific enhancer activity. Notably, AD-associated genetic variants are specifically enriched in increasing-level enhancer orthologs implicating immune processes in AD predisposition. Indeed, increasing enhancers overlap known AD loci lacking protein-altering variants and implicate additional loci that do not reach genome-wide significance. Our results reveal new insights into the mechanisms of neurodegeneration and establish the mouse as a useful model for functional studies of AD regulatory regions.

Species differential regulation of COX2 can be described by an NFκB-dependent logic AND gate

Cyclooxygenase 2 (COX2), a key regulatory enzyme of the prostaglandin/eicosanoid pathway, is an important target for anti-inflammatory therapy. It is highly induced by pro-inflammatory cytokines in a Nuclear factor kappa B (NFκB)-dependent manner. However, the mechanisms determining the amplitude and dynamics of this important pro-inflammatory event are poorly understood. Furthermore, there is significant difference between human and mouse COX2 expression in response to the inflammatory stimulus tumor necrosis factor alpha (TNFα). Here, we report the presence of a molecular logic AND gate composed of two NFκB response elements (NREs) which controls the expression of human COX2 in a switch-like manner. Combining quantitative kinetic modeling and thermostatistical analysis followed by experimental validation in iterative cycles, we show that the human COX2 expression machinery regulated by NFκB displays features of a logic AND gate. We propose that this provides a digital, noise-filtering mechanism for a tighter control of expression in response to TNFα, such that a threshold level of NFκB activation is required before the promoter becomes active and initiates transcription. This NFκB-regulated AND gate is absent in the mouse COX2 promoter, most likely contributing to its differential graded response in promoter activity and protein expression to TNFα. Our data suggest that the NFκB-regulated AND gate acts as a novel mechanism for controlling the expression of human COX2 to TNFα, and its absence in the mouse COX2 provides the foundation for further studies on understanding species-specific differential gene regulation.

A comparative encyclopedia of DNA elements in the mouse genome

The laboratory mouse shares the majority of its protein-coding genes with humans, making it the premier model organism in biomedical research, yet the two mammals differ in significant ways. To gain greater insights into both shared and species-specific transcriptional and cellular regulatory programs in the mouse, the Mouse ENCODE Consortium has mapped transcription, DNase I hypersensitivity, transcription factor binding, chromatin modifications and replication domains throughout the mouse genome in diverse cell and tissue types. By comparing with the human genome, we not only confirm substantial conservation in the newly annotated potential functional sequences, but also find a large degree of divergence of sequences involved in transcriptional regulation, chromatin state and higher order chromatin organization. Our results illuminate the wide range of evolutionary forces acting on genes and their regulatory regions, and provide a general resource for research into mammalian biology and mechanisms of human diseases.

Progenitor stage-specific activity of a cis-acting double GATA motif for Gata1 gene expression

GATA1 is a master regulator of erythropoiesis, expression of which is regulated by multiple discrete cis-acting elements. In this study, we examine the activity of a promoter-proximal double GATA (dbGATA) motif, using a Gata1 bacterial artificial chromosome (BAC)-transgenic green fluorescent protein (GFP) reporter (G1BAC-GFP) mouse system. Deletion of the dbGATA motif led to significant reductions in GFP expression in hematopoietic progenitors, while GFP expression was maintained in erythroblasts. Consistently, in mice with a germ line deletion of the dbGATA motif (Gata1(ΔdbGATA) mice), GATA1 expression in progenitors was significantly decreased. The suppressed GATA1 expression was associated with a compensatory increase in GATA2 levels in progenitors. When we crossed Gata1(ΔdbGATA) mice with Gata2 hypomorphic mutant mice (Gata2(fGN/fGN) mice), the Gata1(ΔdbGATA)::Gata2(fGN/fGN) compound mutant mice succumbed to a significant decrease in the progenitor population, whereas both groups of single mutant mice maintained progenitors and survived to adulthood, indicating the functional redundancy between GATA1 and GATA2 in progenitors. Meanwhile, the effects of the dbGATA site deletion on Gata1 expression were subtle in erythroblasts, which showed increased GATA1 binding and enhanced accumulation of active histone marks around the 1st-intron GATA motif of the ΔdbGATA locus. These results thus reveal a novel role of the dbGATA motif in the maintenance of Gata1 expression in hematopoietic progenitors and a functional compensation between the dbGATA site and the 1st-intron GATA motif in erythroblasts.

Comparison of the transcriptional landscapes between human and mouse tissues

Although the similarities between humans and mice are typically highlighted, morphologically and genetically, there are many differences. To better understand these two species on a molecular level, we performed a comparison of the expression profiles of 15 tissues by deep RNA sequencing and examined the similarities and differences in the transcriptome for both protein-coding and -noncoding transcripts. Although commonalities are evident in the expression of tissue-specific genes between the two species, the expression for many sets of genes was found to be more similar in different tissues within the same species than between species. These findings were further corroborated by associated epigenetic histone mark analyses. We also find that many noncoding transcripts are expressed at a low level and are not detectable at appreciable levels across individuals. Moreover, the majority lack obvious sequence homologs between species, even when we restrict our attention to those which are most highly reproducible across biological replicates. Overall, our results indicate that there is considerable RNA expression diversity between humans and mice, well beyond what was described previously, likely reflecting the fundamental physiological differences between these two organisms.

Widespread contribution of transposable elements to the innovation of gene regulatory networks

Transposable elements (TEs) have been shown to contain functional binding sites for certain transcription factors (TFs). However, the extent to which TEs contribute to the evolution of TF binding sites is not well known. We comprehensively mapped binding sites for 26 pairs of orthologous TFs in two pairs of human and mouse cell lines (representing two cell lineages), along with epigenomic profiles, including DNA methylation and six histone modifications. Overall, we found that 20% of binding sites were embedded within TEs. This number varied across different TFs, ranging from 2% to 40%. We further identified 710 TF–TE relationships in which genomic copies of a TE subfamily contributed a significant number of binding peaks for a TF, and we found that LTR elements dominated these relationships in human. Importantly, TE-derived binding peaks were strongly associated with open and active chromatin signatures, including reduced DNA methylation and increased enhancer-associated histone marks. On average, 66% of TE-derived binding events were cell type-specific with a cell type-specific epigenetic landscape. Most of the binding sites contributed by TEs were species-specific, but we also identified binding sites conserved between human and mouse, the functional relevance of which was supported by a signature of purifying selection on DNA sequences of these TEs. Interestingly, several TFs had significantly expanded binding site landscapes only in one species, which were linked to species-specific gene functions, suggesting that TEs are an important driving force for regulatory innovation. Taken together, our data suggest that TEs have significantly and continuously shaped gene regulatory networks during mammalian evolution.

Dynamic shifts in occupancy by TAL1 are guided by GATA factors and drive large-scale reprogramming of gene expression during hematopoiesis

We used mouse ENCODE data along with complementary data from other laboratories to study the dynamics of occupancy and the role in gene regulation of the transcription factor TAL1, a critical regulator of hematopoiesis, at multiple stages of hematopoietic differentiation. We combined ChIP-seq and RNA-seq data in six mouse cell types representing a progression from multilineage precursors to differentiated erythroblasts and megakaryocytes. We found that sites of occupancy shift dramatically during commitment to the erythroid lineage, vary further during terminal maturation, and are strongly associated with changes in gene expression. In multilineage progenitors, the likely target genes are enriched for hematopoietic growth and functions associated with the mature cells of specific daughter lineages (such as megakaryocytes). In contrast, target genes in erythroblasts are specifically enriched for red cell functions. Furthermore, shifts in TAL1 occupancy during erythroid differentiation are associated with gene repression (dissociation) and induction (co-occupancy with GATA1). Based on both enrichment for transcription factor binding site motifs and co-occupancy determined by ChIP-seq, recruitment by GATA transcription factors appears to be a stronger determinant of TAL1 binding to chromatin than the canonical E-box binding site motif. Studies of additional proteins lead to the model that TAL1 regulates expression after being directed to a distinct subset of genomic binding sites in each cell type via its association with different complexes containing master regulators such as GATA2, ERG, and RUNX1 in multilineage cells and the lineage-specific master regulator GATA1 in erythroblasts.