A global genomic transcriptional code associated with CNS-expressed genes

Two complementary approaches for efficient isolation of Sertoli cells for transcriptomic analysis

Sertoli cells (SCs) are the only somatic cells that reside in seminiferous tubules of testis. They directly interact with and support the development of germ cells, thus have an indispensable role in the process of spermatogenesis. SCs first appear in a proliferative state and then, with the initiation of the first wave of spermatogenesis, progress to a mature “nurturing” state which supports lifelong continuous sperm production. During this development, the SC transcriptome must adapt rapidly as obstacles in SC maturation often result in deficiencies in male fertility. Due to its importance in spermatogenesis, a reliable, rapid, and precise method for the isolation of high purity, viable and unadulterated SC has been largely missing. We have developed an improved method for the preparation of a testicular single cell suspension comprised of two alternative protocols to separate SCs from the rest of the testicular cells by FACS. The first sorting scheme is based on their co-expression of surface specific markers, FSHr and Occludin-1, while the second focuses on the co-staining of SCs with FSHr-specific antibody and Hoechst 33342, which discriminates DNA content of testicular cells. The entire procedure can be completed in less than 3 h which permits the analysis of the development-related transcriptional profile of these cells. Notably, our comparative study showed that this method resulted in a SC transcriptome that is largely comparable to SCs which were briskly isolated due to their cell-specific expression of fluorescent protein. Interestingly, we also show that SCs sorted as FSHr⁺Occludin⁺ cells contained a tangible portion of transcripts from all types of testicular germ cells. Sorting of SCs according to their 2C DNA content significantly reduced the presence of these transcripts, thus seems to be the most suitable approach for accurate determination of the SC transcriptome. We believe that these novel approaches for the isolation of SCs will assist researchers in the elucidation of their function as well as their role in spermatogenesis and disorders related to male infertility.

Regulatory networks specifying cortical interneurons from human embryonic stem cells reveal roles for CHD2 in interneuron development

Cortical interneurons (cINs) modulate excitatory neuronal activity by providing local inhibition. During fetal development, several cIN subtypes derive from the medial ganglionic eminence (MGE), a transient ventral telencephalic structure. While altered cIN development contributes to neurodevelopmental disorders, the inaccessibility of human fetal brain tissue during development has hampered efforts to define molecular networks controlling this process. Here, we modified protocols for directed differentiation of human embryonic stem cells, obtaining efficient, accelerated production of MGE-like progenitors and MGE-derived cIN subtypes with the expected electrophysiological properties. We defined transcriptome changes accompanying this process and integrated these data with direct transcriptional targets of NKX2-1, a transcription factor controlling MGE specification. This analysis defined NKX2-1-associated genes with enriched expression during MGE specification and cIN differentiation, including known and previously unreported transcription factor targets with likely roles in MGE specification, and other target classes regulating cIN migration and function. NKX2-1-associated peaks were enriched for consensus binding motifs for NKX2-1, LHX, and SOX transcription factors, suggesting roles in coregulating MGE gene expression. Among the NKX2-1 direct target genes with cIN-enriched expression was CHD2, which encodes a chromatin remodeling protein mutated to cause human epilepsies. Accordingly, CHD2 deficiency impaired cIN specification and altered later electrophysiological function, while CHD2 coassociated with NKX2-1 at cis-regulatory elements and was required for their transactivation by NKX2-1 in MGE-like progenitors. This analysis identified several aspects of gene-regulatory networks underlying human MGE specification and suggested mechanisms by which NKX2-1 acts with chromatin remodeling activities to regulate gene expression programs underlying cIN development.

Pax3- and Pax7-mediated Dbx1 regulation orchestrates the patterning of intermediate spinal interneurons

Transcription factors are key orchestrators of the emergence of neuronal diversity within the developing spinal cord. As such, the two paralogous proteins Pax3 and Pax7 regulate the specification of progenitor cells within the intermediate neural tube, by defining a neat segregation between those fated to form motor circuits and those involved in the integration of sensory inputs. To attain insights into the molecular means by which they control this process, we have performed detailed phenotypic analyses of the intermediate spinal interneurons (IN), namely the dI6, V0_D, V0_VCG and V1 populations in compound null mutants for Pax3 and Pax7. This has revealed that the levels of Pax3/7 proteins determine both the dorso-ventral extent and the number of cells produced in each subpopulation; with increasing levels leading to the dorsalisation of their fate. Furthermore, thanks to the examination of mutants in which Pax3 transcriptional activity is skewed either towards repression or activation, we demonstrate that this cell diversification process is mainly dictated by Pax3/7 ability to repress gene expression. Consistently, we show that Pax3 and Pax7 inhibit the expression of Dbx1 and of its repressor Prdm12, fate determinants of the V0 and V1 interneurons, respectively. Notably, we provide evidence for the activity of several cis-regulatory modules of Dbx1 to be sensitive to Pax3 and Pax7 transcriptional activity levels. Altogether, our study provides insights into how the redundancy within a TF family, together with discrete dynamics of expression profiles of each member, are exploited to generate cellular diversity. Furthermore, our data supports the model whereby cell fate choices in the neural tube do not rely on binary decisions but rather on inhibition of multiple alternative fates.

Neural Progenitors Adopt Specific Identities by Directly Repressing All Alternative Progenitor Transcriptional Programs

In the vertebrate neural tube, a morphogen-induced transcriptional network produces multiple molecularly distinct progenitor domains, each generating different neuronal subtypes. Using an in vitro differentiation system, we defined gene expression signatures of distinct progenitor populations and identified direct gene-regulatory inputs corresponding to locations of specific transcription factor binding. Combined with targeted perturbations of the network, this revealed a mechanism in which a progenitor identity is installed by active repression of the entire transcriptional programs of other neural progenitor fates. In the ventral neural tube, sonic hedgehog (Shh) signaling, together with broadly expressed transcriptional activators, concurrently activates the gene expression programs of several domains. The specific outcome is selected by repressive input provided by Shh-induced transcription factors that act as the key nodes in the network, enabling progenitors to adopt a single definitive identity from several initially permitted options. Together, the data suggest design principles relevant to many developing tissues.

•

Specific vertebrate neural progenitor populations generated in vitro

•

Gene expression dynamics, transcription factor binding assessed in neural progenitors

•

Progenitor fate selected by repressors blocking entire programs of other identities

•

Repressors counteract non-selective morphogen and pan-neural activatory inputs

Schwann cells and their transcriptional network: Evolution of key regulators of peripheral myelination

As derivatives of the neural crest, Schwann cells represent a vertebrate invention. Their development and differentiation is under control of a newly constructed, vertebrate-specific regulatory network that contains Sox10, Oct6 and Krox20 as cornerstones and central regulators of peripheral myelination. In this review, we discuss the function and relationship of these transcription factors among each other and in the context of their regulatory network, and present ideas of how neofunctionalization may have helped to recruit them to their novel task in Schwann cells. This article is part of a Special Issue entitled SI: Myelin Evolution.

A theoretical framework for the regulation of Shh morphogen-controlled gene expression

How morphogen gradients govern the pattern of gene expression in developing tissues is not well understood. Here, we describe a statistical thermodynamic model of gene regulation that combines the activity of a morphogen with the transcriptional network it controls. Using Sonic hedgehog (Shh) patterning of the ventral neural tube as an example, we show that the framework can be used together with the principled parameter selection technique of approximate Bayesian computation to obtain a dynamical model that accurately predicts tissue patterning. The analysis indicates that, for each target gene regulated by Gli, which is the transcriptional effector of Shh signalling, there is a neutral point in the gradient, either side of which altering the Gli binding affinity has opposite effects on gene expression. This explains recent counterintuitive experimental observations. The approach is broadly applicable and provides a unifying framework to explain the temporospatial pattern of morphogen-regulated gene expression.

Program specificity for Ptf1a in pancreas versus neural tube development correlates with distinct collaborating cofactors and chromatin accessibility

The lineage-specific basic helix-loop-helix transcription factor Ptf1a is a critical driver for development of both the pancreas and nervous system. How one transcription factor controls diverse programs of gene expression is a fundamental question in developmental biology. To uncover molecular strategies for the program-specific functions of Ptf1a, we identified bound genomic regions in vivo during development of both tissues. Most regions bound by Ptf1a are specific to each tissue, lie near genes needed for proper formation of each tissue, and coincide with regions of open chromatin. The specificity of Ptf1a binding is encoded in the DNA surrounding the Ptf1a-bound sites, because these regions are sufficient to direct tissue-restricted reporter expression in transgenic mice. Fox and Sox factors were identified as potential lineage-specific modifiers of Ptf1a binding, since binding motifs for these factors are enriched in Ptf1a-bound regions in pancreas and neural tube, respectively. Of the Fox factors expressed during pancreatic development, Foxa2 plays a major role. Indeed, Ptf1a and Foxa2 colocalize in embryonic pancreatic chromatin and can act synergistically in cell transfection assays. Together, these findings indicate that lineage-specific chromatin landscapes likely constrain the DNA binding of Ptf1a, and they identify Fox and Sox gene families as part of this process.

Mechanistic differences in the transcriptional interpretation of local and long-range Shh morphogen signaling

Morphogens orchestrate tissue patterning in a concentration-dependent fashion during vertebrate embryogenesis, yet little is known of how positional information provided by such signals is translated into discrete transcriptional outputs. Here we have identified and characterized cis-regulatory modules (CRMs) of genes operating downstream of graded Shh signaling and bifunctional Gli proteins in neural patterning. Unexpectedly, we find that Gli activators have a noninstructive role in long-range patterning and cooperate with SoxB1 proteins to facilitate a largely concentration-independent mode of gene activation. Instead, the opposing Gli-repressor gradient is interpreted at transcriptional levels, and, together with CRM-specific repressive input of homeodomain proteins, comprises a repressive network that translates graded Shh signaling into regional gene expression patterns. Moreover, local and long-range interpretation of Shh signaling differs with respect to CRM context sensitivity and Gli-activator dependence, and we propose that these differences provide insight into how morphogen function may have mechanistically evolved from an initially binary inductive event.

The methylmalonic aciduria related genes, Mmaa, Mmab, and Mut, are broadly expressed in placental and embryonic tissues during mouse organogenesis

Organ-specific birth defects are seen in patients with some inborn errors of vitamin B(12) metabolism. To determine whether three mouse genes, whose human counterparts are associated with isolated methylmalonic aciduria (Mmaa, Mmab and Mut), show tissue-specific expression during organogenesis, we used in situ hybridization to characterize their pattern of expression in wild type embryos and placentas at embryonic days (E) E10.5, E11.5 and E12.5. These three genes are ubiquitously expressed in the placenta and in embryos at E10.5. At E11.5, we observed tissue specific expression patterns for these three genes in lung, head and Rathke's pouch. At E12.5, although Mut expression was ubiquitous, we found cell-type specific expression patterns for Mmaa and Mmab in the developing craniofacial region, the lung, the liver, and the gut. These results suggest that during organogenesis the proteins encoded by these three genes may interact in only a subset of cells.

SoxC transcription factors are required for neuronal differentiation in adult hippocampal neurogenesis

Neural stem cells (NSCs) generate new hippocampal dentate granule neurons throughout adulthood. The genetic programs controlling neuronal differentiation of adult NSCs are only poorly understood. Here we show that, in the adult mouse hippocampus, expression of the SoxC transcription factors Sox4 and Sox11 is initiated around the time of neuronal commitment of adult NSCs and is maintained in immature neurons. Overexpression of Sox4 and Sox11 strongly promotes in vitro neurogenesis from adult NSCs, whereas ablation of Sox4/Sox11 prevents in vitro and in vivo neurogenesis from adult NSCs. Moreover, we demonstrate that SoxC transcription factors target the promoters of genes that are induced on neuronal differentiation of adult NSCs. Finally, we show that reprogramming of astroglia into neurons is dependent on the presence of SoxC factors. These data identify SoxC proteins as essential contributors to the genetic network controlling neuronal differentiation in adult neurogenesis and neuronal reprogramming of somatic cells.

Position effect variegation and epigenetic modification of a transgene in a pig model

Sequences proximal to transgene integration sites are able to regulate transgene expression, resulting in complex position effect variegation. Position effect variegation can cause differences in epigenetic modifications, such as DNA methylation and histone acetylation. However, it is not known which factor, position effect or epigenetic modification, plays a more important role in the regulation of transgene expression. We analyzed transgene expression patterns and epigenetic modifications of transgenic pigs expressing green fluorescent protein, driven by the cytomegalovirus (CMV) promoter. DNA hypermethylation and loss of acetylation of specific histone H3 and H4 lysines, except H4K16 acetylation in the CMV promoter, were consistent with a low level of transgene expression. Moreover, the degree of DNA methylation and histone H3/H4 acetylation in the promoter region depended on the integration site; consequently, position effect variegation caused variations in epigenetic modifications. The transgenic pig fibroblast cell lines were treated with DNA methyltransferase inhibitor 5-Aza-2'-deoxycytidine and/or histone deacetylase inhibitor trichostatin A. Transgene expression was promoted by reversing the DNA hypermethylation and histone hypoacetylation status. The differences in DNA methylation and histone acetylation in the CMV promoter region in these cell lines were not significant; however, significant differences in transgene expression were detected, demonstrating that variegation of transgene expression is affected by the integration site. We conclude that in this pig model, position effect variegation affects transgene expression.

Ancient Pbx-Hox signatures define hundreds of vertebrate developmental enhancers

BACKGROUND

Gene regulation through cis-regulatory elements plays a crucial role in development and disease. A major aim of the post-genomic era is to be able to read the function of cis-regulatory elements through scrutiny of their DNA sequence. Whilst comparative genomics approaches have identified thousands of putative regulatory elements, our knowledge of their mechanism of action is poor and very little progress has been made in systematically de-coding them.

RESULTS

Here, we identify ancient functional signatures within vertebrate conserved non-coding elements (CNEs) through a combination of phylogenetic footprinting and functional assay, using genomic sequence from the sea lamprey as a reference. We uncover a striking enrichment within vertebrate CNEs for conserved binding-site motifs of the Pbx-Hox hetero-dimer. We further show that these predict reporter gene expression in a segment specific manner in the hindbrain and pharyngeal arches during zebrafish development.

CONCLUSIONS

These findings evoke an evolutionary scenario in which many CNEs evolved early in the vertebrate lineage to co-ordinate Hox-dependent gene-regulatory interactions that pattern the vertebrate head. In a broader context, our evolutionary analyses reveal that CNEs are composed of tightly linked transcription-factor binding-sites (TFBSs), which can be systematically identified through phylogenetic footprinting approaches. By placing a large number of ancient vertebrate CNEs into a developmental context, our findings promise to have a significant impact on efforts toward de-coding gene-regulatory elements that underlie vertebrate development, and will facilitate building general models of regulatory element evolution.

Prediction of RNA Polymerase II recruitment, elongation and stalling from histone modification data

BACKGROUND

Initiation and elongation of RNA polymerase II (RNAPII) transcription is regulated by both DNA sequence and chromatin signals. Recent breakthroughs make it possible to measure the chromatin state and activity of core promoters genome-wide, but dedicated computational strategies are needed to progress from descriptive annotation of data to quantitative, predictive models.

RESULTS

Here, we describe a computational framework which with high accuracy can predict the locations of core promoters, the amount of recruited RNAPII at the promoter, the amount of elongating RNAPII in the gene body, the mRNA production originating from the promoter and finally also the stalling characteristics of RNAPII by considering both quantitative and spatial features of histone modifications around the transcription start site (TSS).As the model framework can also pinpoint the signals that are the most influential for prediction, it can be used to infer underlying regulatory biology. For example, we show that the H3K4 di- and tri- methylation signals are strongly predictive for promoter location while the acetylation marks H3K9 and H3K27 are highly important in estimating the promoter usage. All of these four marks are found to be necessary for recruitment of RNAPII but not sufficient for the elongation. We also show that the spatial distributions of histone marks are almost as predictive as the signal strength and that a set of histone marks immediately downstream of the TSS is highly predictive of RNAPII stalling.

CONCLUSIONS

In this study we introduce a general framework to accurately predict the level of RNAPII recruitment, elongation, stalling and mRNA expression from chromatin signals. The versatility of the method also makes it ideally suited to investigate other genomic data.

Mouse transgenesis identifies conserved functional enhancers and cis-regulatory motif in the vertebrate LIM homeobox gene Lhx2 locus

The vertebrate Lhx2 is a member of the LIM homeobox family of transcription factors. It is essential for the normal development of the forebrain, eye, olfactory system and liver as well for the differentiation of lymphoid cells. However, despite the highly restricted spatio-temporal expression pattern of Lhx2, nothing is known about its transcriptional regulation. In mammals and chicken, Crb2, Dennd1a and Lhx2 constitute a conserved linkage block, while the intervening Dennd1a is lost in the fugu Lhx2 locus. To identify functional enhancers of Lhx2, we predicted conserved noncoding elements (CNEs) in the human, mouse and fugu Crb2-Lhx2 loci and assayed their function in transgenic mouse at E11.5. Four of the eight CNE constructs tested functioned as tissue-specific enhancers in specific regions of the central nervous system and the dorsal root ganglia (DRG), recapitulating partial and overlapping expression patterns of Lhx2 and Crb2 genes. There was considerable overlap in the expression domains of the CNEs, which suggests that the CNEs are either redundant enhancers or regulating different genes in the locus. Using a large set of CNEs (810 CNEs) associated with transcription factor-encoding genes that express predominantly in the central nervous system, we predicted four over-represented 8-mer motifs that are likely to be associated with expression in the central nervous system. Mutation of one of them in a CNE that drove reporter expression in the neural tube and DRG abolished expression in both domains indicating that this motif is essential for expression in these domains. The failure of the four functional enhancers to recapitulate the complete expression pattern of Lhx2 at E11.5 indicates that there must be other Lhx2 enhancers that are either located outside the region investigated or divergent in mammals and fishes. Other approaches such as sequence comparison between multiple mammals are required to identify and characterize such enhancers.

Class-C SOX transcription factors control GnRH gene expression via the intronic transcriptional enhancer

GnRH is a pivotal hypothalamic neurohormone governing reproduction and sexual development. Because transcriptional regulation is crucial for the spatial and temporal expression of the GnRH gene, a region approximately 3.0 kb upstream of the mammalian GnRH promoter has been extensive studied. In the present study, we demonstrate a transcription-enhancer located in the first intron (intron A) region of the GnRH gene. This transcriptional enhancer harbors putative sex-determining region Y-related high-mobility-group box (SOX) family transcription factor-binding sites, which are well conserved across many mammalian species. The class-C SOX member proteins (SOX-C) (SOX4 and SOX11) specifically augment this transcriptional activation by binding to these SOX-binding sites. In accordance, SOX11 is highly enriched in immortalized GnRH-producing GT1-1 cells, and suppression of its expression significantly decreases GnRH gene expression as well as GnRH secretion. Chromatin immunoprecipitation shows that endogenous SOX-C factors recognize and bind to the intronic enhancer in GT1-1 cells and the hypothalamus. Accompanying immunohistochemical analysis demonstrates that SOX4 or SOX11 are highly expressed in the majority of hypothalamic GnRH neurons in adult mice. Taken together, these findings demonstrate that SOX-C transcription factors function as important transcriptional regulators of cell type-specific GnRH gene expression by acting on the intronic transcriptional enhancer.

Common motifs shared by conserved enhancers of Drosophila midline glial genes

Coding sequences are usually the most highly conserved sectors of DNA, but genomic regions controlling the expression pattern of certain genes can also be conserved across diverse species. In this study, we identify five enhancers capable of activating transcription in the midline glia of Drosophila melanogaster and each contains sequences conserved across at least 11 Drosophila species. In addition, the conserved sequences contain reiterated motifs for binding sites of the known midline transcriptional activators, Single-minded, Tango, Dichaete, and Pointed. To understand the molecular basis for the highly conserved genomic subregions within enhancers of the midline genes, we tested the ability of various motifs to affect midline expression, both individually and in combination, within synthetic reporter constructs. Multiple copies of the binding site for the midline regulators Single-minded and Tango can drive expression in midline cells; however, small changes to the sequences flanking this transcription factor binding site can inactivate expression in midline cells and activate expression in tracheal cells instead. For the midline genes described in this study, the highly conserved sequences appear to juxtapose positive and negative regulatory factors in a configuration that activates genes specifically in the midline glia, while maintaining them inactive in other tissues, including midline neurons and tracheal cells.

Regulation of retinal homeobox gene transcription by cooperative activity among cis-elements

The retinal homeobox (Rx/rax) gene is essential for the development of the eye. Rax is among the earliest genes expressed during eye development, beginning in the prospective eye fields in the anterior neural plate. Additionally Rax expression persists in retinal progenitor cells and in differentiated photoreceptors. We have isolated and characterized a 2.8 kb genomic DNA fragment that regulates expression of Rax in the developing and maturing retina. We have discovered and characterized cis-acting elements that function to specifically control spatial and temporal Rax expression during retinal development. We have found that the regulation of Rax2A promoter activity requires cooperative interactions between positive and negative regulatory elements. Further, a highly conserved genomic element containing SOX, OTX, and POU transcription factor binding sites is necessary but not sufficient for promoter activity in retinal progenitor or stem cells. Finally, a putative binding element for forkhead transcription factors is necessary for promoter activity and can cooperate with other cis-acting elements to drive Rax2A promoter activity.

Identification of interacting transcription factors regulating tissue gene expression in human

BACKGROUND

Tissue gene expression is generally regulated by multiple transcription factors (TFs). A major first step toward understanding how tissues achieve their specificity is to identify, at the genome scale, interacting TFs regulating gene expression in different tissues. Despite previous discoveries, the mechanisms that control tissue gene expression are not fully understood.

RESULTS

We have integrated a function conservation approach, which is based on evolutionary conservation of biological function, and genes with highest expression level in human tissues to predict TF pairs controlling tissue gene expression. To this end, we have identified 2549 TF pairs associated with a certain tissue. To find interacting TFs controlling tissue gene expression in a broad spatial and temporal manner, we looked for TF pairs common to the same type of tissues and identified 379 such TF pairs, based on which TF-TF interaction networks were further built. We also found that tissue-specific TFs may play an important role in recruiting non-tissue-specific TFs to the TF-TF interaction network, offering the potential for coordinating and controlling tissue gene expression across a variety of conditions.

CONCLUSION

The findings from this study indicate that tissue gene expression is regulated by large sets of interacting TFs either on the same promoter of a gene or through TF-TF interaction networks.

Conservation and regulatory associations of a wide affinity range of mouse transcription factor binding sites

Sequence-specific binding by transcription factors (TFs) interprets regulatory information encoded in the genome. Using recently published universal protein binding microarray (PBM) data on the in vitro DNA binding preferences of these proteins for all possible 8-base-pair sequences, we examined the evolutionary conservation and enrichment within putative regulatory regions of the binding sequences of a diverse library of 104 nonredundant mouse TFs spanning 22 different DNA-binding domain structural classes. We found that not only high affinity binding sites, but also numerous moderate and low affinity binding sites, are under negative selection in the mouse genome. These 8-mers occur preferentially in putative regulatory regions of the mouse genome, including CpG islands and non-exonic ultraconserved elements (UCEs). Of TFs whose PBM "bound" 8-mers are enriched within sets of tissue-specific UCEs, many are expressed in the same tissue(s) as the UCE-driven gene expression. Phylogenetically conserved motif occurrences of various TFs were also enriched in the noncoding sequence surrounding numerous gene sets corresponding to Gene Ontology categories and tissue-specific gene expression clusters, suggesting involvement in transcriptional regulation of those genes. Altogether, our results indicate that many of the sequences bound by these proteins in vitro, including lower affinity DNA sequences, are likely to be functionally important in vivo. This study not only provides an initial analysis of the potential regulatory associations of 104 mouse TFs, but also presents an approach for the functional analysis of TFs from any other metazoan genome as their DNA binding preferences are determined by PBMs or other technologies.

Establishing and interpreting graded Sonic Hedgehog signaling during vertebrate neural tube patterning: the role of negative feedback

The secreted protein Sonic Hedgehog (SHH) acts in graded fashion to pattern the dorsal-ventral axis of the vertebrate neural tube. This is a dynamic process in which increasing concentrations and durations of exposure to SHH generate neurons with successively more ventral identities. Interactions between the receiving cells and the graded signal underpin the mechanism of SHH action. In particular, negative feedback, involving proteins transcriptionally induced or repressed by SHH signaling, plays an essential role in shaping the graded readout. On one hand, negative feedback controls, in a noncell-autonomous manner, the distribution of SHH across the field of receiving cells. On the other, it acts cell-autonomously to convert different concentrations of SHH into distinct durations of intracellular signal transduction. Together, these mechanisms exemplify a strategy for morphogen interpretation, which we have termed temporal adaptation that relies on the continuous processing and refinement of the cellular response to the graded signal.

Motifs and cis-regulatory modules mediating the expression of genes co-expressed in presynaptic neurons

An integrative strategy of comparative genomics, experimental and computational approaches reveals aspects of a regulatory network controlling neuronal-specific expression in presynaptic neurons.

Background

Hundreds of proteins modulate neurotransmitter release and synaptic plasticity during neuronal development and in response to synaptic activity. The expression of genes in the pre- and post-synaptic neurons is under stringent spatio-temporal control, but the mechanism underlying the neuronal expression of these genes remains largely unknown.

Results

Using unbiased in vivo and in vitro screens, we characterized the cis elements regulating the Rab3A gene, which is expressed abundantly in presynaptic neurons. A set of identified regulatory elements of the Rab3A gene corresponded to the defined Rab3A multi-species conserved elements. In order to identify clusters of enriched transcription factor binding sites, for example, cis-regulatory modules, we analyzed intergenic multi-species conserved elements in the vicinity of nine presynaptic genes, including Rab3A, that are highly and specifically expressed in brain regions. Sixteen transcription factor binding motifs were over-represented in these multi-species conserved elements. Based on a combined occurrence for these enriched motifs, multi-species conserved elements in the vicinity of 107 previously identified presynaptic genes were scored and ranked. We then experimentally validated the scoring strategy by showing that 12 of 16 (75%) high-scoring multi-species conserved elements functioned as neuronal enhancers in a cell-based assay.

Conclusions

This work introduces an integrative strategy of comparative genomics, experimental, and computational approaches to reveal aspects of a regulatory network controlling neuronal-specific expression of genes in presynaptic neurons.

Identification and characterization of new long conserved noncoding sequences in vertebrates

Comparative sequence analyses have identified highly conserved genomic DNA sequences, including noncoding sequences, between humans and other species. By performing whole-genome comparisons of human and mouse, we have identified 611 conserved noncoding sequences longer than 500 bp, with more than 95% identity between the species. These long conserved noncoding sequences (LCNS) include 473 new sequences that do not overlap with previously reported ultraconserved elements (UCE), which are defined as aligned sequences longer than 200 bp with 100% identity in human, mouse, and rat. The LCNS were distributed throughout the genome except for the Y chromosome and often occurred in clusters within regions with a low density of coding genes. Many of the LCNS were also highly conserved in other mammals, chickens, frogs, and fish; however, we were unable to find orthologous sequences in the genomes of invertebrate species. In order to examine whether these conserved sequences are functionally important or merely mutational cold spots, we directly measured the frequencies of ENU-induced germline mutations in the LCNS of the mouse. By screening about 40.7 Mb, we found 35 mutations, including mutations at nucleotides that were conserved between human and fish. The mutation frequencies were equivalent to those found in other genomic regions, including coding sequences and introns, suggesting that the LCNS are not mutational cold spots at all. Taken together, these results suggest that mutations occur with equal frequency in LCNS but are eliminated by natural selection during the course of evolution.

Pattern formation in the vertebrate neural tube: a sonic hedgehog morphogen-regulated transcriptional network

Neuronal subtype specification in the vertebrate neural tube is one of the best-studied examples of embryonic pattern formation. Distinct neuronal subtypes are generated in a precise spatial order from progenitor cells according to their location along the anterior-posterior and dorsal-ventral axes. Underpinning this organization is a complex network of multiple extrinsic and intrinsic factors. This review focuses on the molecular mechanisms and general strategies at play in ventral regions of the forming spinal cord, where sonic hedgehog-based morphogen signaling is a key determinant. We discuss recent advances in our understanding of these events and highlight unresolved questions.

Host proteins interacting with the Moloney murine leukemia virus integrase: multiple transcriptional regulators and chromatin binding factors

Background

A critical step for retroviral replication is the stable integration of the provirus into the genome of its host. The viral integrase protein is key in this essential step of the retroviral life cycle. Although the basic mechanism of integration by mammalian retroviruses has been well characterized, the factors determining how viral integration events are targeted to particular regions of the genome or to regions of a particular DNA structure remain poorly defined. Significant questions remain regarding the influence of host proteins on the selection of target sites, on the repair of integration intermediates, and on the efficiency of integration.

Results

We describe the results of a yeast two-hybrid screen using Moloney murine leukemia virus integrase as bait to screen murine cDNA libraries for host proteins that interact with the integrase. We identified 27 proteins that interacted with different integrase fusion proteins. The identified proteins include chromatin remodeling, DNA repair and transcription factors (13 proteins); translational regulation factors, helicases, splicing factors and other RNA binding proteins (10 proteins); and transporters or miscellaneous factors (4 proteins). We confirmed the interaction of these proteins with integrase by testing them in the context of other yeast strains with GAL4-DNA binding domain-integrase fusions, and by in vitro binding assays between recombinant proteins. Subsequent analyses revealed that a number of the proteins identified as Mo-MLV integrase interactors also interact with HIV-1 integrase both in yeast and in vitro.

Conclusion

We identify several proteins interacting directly with both MoMLV and HIV-1 integrases that may be common to the integration reaction pathways of both viruses. Many of the proteins identified in the screen are logical interaction partners for integrase, and the validity of a number of the interactions are supported by other studies. In addition, we observe that some of the proteins have documented interactions with other viruses, raising the intriguing possibility that there may be common host proteins used by different viruses. We undertook this screen to identify host factors that might affect integration target site selection, and find that our screens have generated a wealth of putative interacting proteins that merit further investigation.

SoxB1 transcription factors and Notch signaling use distinct mechanisms to regulate proneural gene function and neural progenitor differentiation

The preservation of a pool of neural precursors is a prerequisite for proper establishment and maintenance of a functional central nervous system(CNS). Both Notch signaling and SoxB1 transcription factors have been ascribed key roles during this process, but whether these factors use common or distinct mechanisms to control progenitor maintenance is unsettled. Here, we report that the capacity of Notch to maintain neural cells in an undifferentiated state requires the activity of SoxB1 proteins, whereas the mechanism by which SoxB1 block neurogenesis is independent of Notch signaling. A common feature of Notch signaling and SoxB1 proteins is their ability to inhibit the activity of proneural bHLH proteins. Notch represses the transcription of proneural bHLH genes, while SoxB1 proteins block their neurogenic capacity. Moreover, E-proteins act as functional partners of proneural proteins and the suppression of E-protein expression is an important mechanism by which Notch counteracts neurogenesis. Interestingly, in contrast to the Hes-dependent repression of proneural genes, suppression of E-protein occurs in a Hes-independent fashion. Together, these data reveal that Notch signaling and SoxB1 transcription factors use distinct regulatory mechanisms to control proneural protein function and to preserve neural cells as undifferentiated precursors.

How to build transcriptional network models of mammalian pattern formation

Background

Genetic regulatory networks of sequence specific transcription factors underlie pattern formation in multicellular organisms. Deciphering and representing the mammalian networks is a central problem in development, neurobiology, and regenerative medicine. Transcriptional networks specify intermingled embryonic cell populations during pattern formation in the vertebrate neural tube. Each embryonic population gives rise to a distinct type of adult neuron. The homeodomain transcription factor Lbx1 is expressed in five such populations and loss of Lbx1 leads to distinct respecifications in each of the five populations.

Methodology/Principal Findings

We have purified normal and respecified pools of these five populations from embryos bearing one or two copies of the null Lbx1^GFP allele, respectively. Microarrays were used to show that expression levels of 8% of all transcription factor genes were altered in the respecified pool. These transcription factor genes constitute 20–30% of the active nodes of the transcriptional network that governs neural tube patterning. Half of the 141 regulated nodes were located in the top 150 clusters of ultraconserved non-coding regions. Generally, Lbx1 repressed genes that have expression patterns outside of the Lbx1-expressing domain and activated genes that have expression patterns inside the Lbx1-expressing domain.

Conclusions/Significance

Constraining epistasis analysis of Lbx1 to only those cells that normally express Lbx1 allowed unprecedented sensitivity in identifying Lbx1 network interactions and allowed the interactions to be assigned to a specific set of cell populations. We call this method ANCEA, or active node constrained epistasis analysis, and think that it will be generally useful in discovering and assigning network interactions to specific populations. We discuss how ANCEA, coupled with population partitioning analysis, can greatly facilitate the systematic dissection of transcriptional networks that underlie mammalian patterning.

Identification of motifs that are conserved in 12 Drosophila species and regulate midline glia vs. neuron expression

Functional complexity of the central nervous system (CNS) is reflected by the large number and diversity of genes expressed in its many different cell types. Understanding the control of gene expression within cells of the CNS will help reveal how various neurons and glia develop and function. Midline cells of Drosophila differentiate into glial cells and several types of neurons and also serve as a signaling center for surrounding tissues. Here, we examine regulation of the midline gene, wrapper, required for both neuron-glia interactions and viability of midline glia. We identify a region upstream of wrapper required for midline expression that is highly conserved (87%) between 12 Drosophila species. Site-directed mutagenesis identifies four motifs necessary for midline glial expression: (1) a Single-minded/Tango binding site, (2) a motif resembling a pointed binding site, (3) a motif resembling a Sox binding site, and (4) a novel motif. An additional highly conserved 27 bp are required to restrict expression to midline glia and exclude it from midline neurons. These results suggest short, highly conserved genomic sequences flanking Drosophila midline genes are indicative of functional regulatory regions and that small changes within these sequences can alter the expression pattern of a gene.

Retroviral enhancer detection insertions in zebrafish combined with comparative genomics reveal genomic regulatory blocks - a fundamental feature of vertebrate genomes

A large-scale enhancer detection screen was performed in the zebrafish using a retroviral vector carrying a basal promoter and a fluorescent protein reporter cassette. Analysis of insertional hotspots uncovered areas around developmental regulatory genes in which an insertion results in the same global expression pattern, irrespective of exact position. These areas coincide with vertebrate chromosomal segments containing identical gene order; a phenomenon known as conserved synteny and thought to be a vestige of evolution. Genomic comparative studies have found large numbers of highly conserved noncoding elements (HCNEs) spanning these and other loci. HCNEs are thought to act as transcriptional enhancers based on the finding that many of those that have been tested direct tissue specific expression in transient or transgenic assays. Although gene order in hox and other gene clusters has long been known to be conserved because of shared regulatory sequences or overlapping transcriptional units, the chromosomal areas found through insertional hotspots contain only one or a few developmental regulatory genes as well as phylogenetically unrelated genes. We have termed these regions genomic regulatory blocks (GRBs), and show that they underlie the phenomenon of conserved synteny through all sequenced vertebrate genomes. After teleost whole genome duplication, a subset of GRBs were retained in two copies, underwent degenerative changes compared with tetrapod loci that exist as single copy, and that therefore can be viewed as representing the ancestral form. We discuss these findings in light of evolution of vertebrate chromosomal architecture and the identification of human disease mutations.

Back to basics: Sox genes

Sox genes are indispensable for multiple aspects of development. This primer briefly describes shared properties of the Sox gene family, and five well-characterized examples of vertebrate developmental mechanisms governed by Sox gene subgroups: testis development, central nervous system neurogenesis, oligodendrocyte development, chondrogenesis, and neural crest cell development. Also featured is an interview about current issues in the field with experts Jonas Muhr, Ph.D. and Robert Kelsh, Ph.D.

JASPAR, the open access database of transcription factor-binding profiles: new content and tools in the 2008 update

JASPAR is a popular open-access database for matrix models describing DNA-binding preferences for transcription factors and other DNA patterns. With its third major release, JASPAR has been expanded and equipped with additional functions aimed at both casual and power users. The heart of the JASPAR database-the JASPAR CORE sub-database-has increased by 12% in size, and three new specialized sub-databases have been added. New functions include clustering of matrix models by similarity, generation of random matrices by sampling from selected sets of existing models and a language-independent Web Service applications programming interface for matrix retrieval. JASPAR is available at http://jaspar.genereg.net.

Genomic regulatory blocks underlie extensive microsynteny conservation in insects

Insect genomes contain larger blocks of conserved gene order (microsynteny) than would be expected under a random breakage model of chromosome evolution. We present evidence that microsynteny has been retained to keep large arrays of highly conserved noncoding elements (HCNEs) intact. These arrays span key developmental regulatory genes, forming genomic regulatory blocks (GRBs). We recently described GRBs in vertebrates, where most HCNEs function as enhancers and HCNE arrays specify complex expression programs of their target genes. Here we present a comparison of five Drosophila genomes showing that HCNE density peaks centrally in large synteny blocks containing multiple genes. Besides developmental regulators that are likely targets of HCNE enhancers, HCNE arrays often span unrelated neighboring genes. We describe differences in core promoters between the target genes and the unrelated genes that offer an explanation for the differences in their responsiveness to enhancers. We show examples of a striking correspondence between boundaries of synteny blocks, HCNE arrays, and Polycomb binding regions, confirming that the synteny blocks correspond to regulatory domains. Although few noncoding elements are highly conserved between Drosophila and the malaria mosquito Anopheles gambiae, we find that A. gambiae regions orthologous to Drosophila GRBs contain an equivalent distribution of noncoding elements highly conserved in the yellow fever mosquito Aëdes aegypti and coincide with regions of ancient microsynteny between Drosophila and mosquitoes. The structural and functional equivalence between insect and vertebrate GRBs marks them as an ancient feature of metazoan genomes and as a key to future studies of development and gene regulation.

Mammalian microRNA prediction through a support vector machine model of sequence and structure

Background

MicroRNAs (miRNAs) are endogenous small noncoding RNA gene products, on average 22 nt long, found in a wide variety of organisms. They play important regulatory roles by targeting mRNAs for degradation or translational repression. There are 377 known mouse miRNAs and 475 known human miRNAs in the May 2007 release of the miRBase database, the majority of which are conserved between the two species. A number of recent reports imply that it is likely that many mammalian miRNAs remain to be discovered. The possibility that there are more of them expressed at lower levels or in more specialized expression contexts calls for the exploitation of genome sequence information to accelerate their discovery.

Methodology/Principal Findings

In this article, we describe a computational method-mirCoS-that uses three support vector machine models sequentially to discover new miRNA candidates in mammalian genomes based on sequence, secondary structure, and conservation. mirCoS can efficiently detect the majority of known miRNAs and predicts an extensive set of hairpin structures based on human-mouse comparisons. In total, 3476 mouse candidates and 3441 human candidates were found. These hairpins are more similar to known miRNAs than to negative controls in several aspects not considered by the prediction algorithm. A significant fraction of predictions is supported by existing expression evidence.

Conclusions/Significance

Using a novel approach, mirCoS performs comparably to or better than existing miRNA prediction methods, and contributes a significant number of new candidate miRNAs for experimental verification.

Genome of the marsupial Monodelphis domestica reveals innovation in non-coding sequences

We report a high-quality draft of the genome sequence of the grey, short-tailed opossum (Monodelphis domestica). As the first metatherian ('marsupial') species to be sequenced, the opossum provides a unique perspective on the organization and evolution of mammalian genomes. Distinctive features of the opossum chromosomes provide support for recent theories about genome evolution and function, including a strong influence of biased gene conversion on nucleotide sequence composition, and a relationship between chromosomal characteristics and X chromosome inactivation. Comparison of opossum and eutherian genomes also reveals a sharp difference in evolutionary innovation between protein-coding and non-coding functional elements. True innovation in protein-coding genes seems to be relatively rare, with lineage-specific differences being largely due to diversification and rapid turnover in gene families involved in environmental interactions. In contrast, about 20% of eutherian conserved non-coding elements (CNEs) are recent inventions that postdate the divergence of Eutheria and Metatheria. A substantial proportion of these eutherian-specific CNEs arose from sequence inserted by transposable elements, pointing to transposons as a major creative force in the evolution of mammalian gene regulation.