A Framework for the Establishment of a Cnidarian Gene Regulatory Network for “Endomesoderm” Specification: The Inputs of ß-Catenin/TCF Signaling

A Framework for the Establishment of a Cnidarian Gene Regulatory Network for “Endomesoderm” Specification: The Inputs of ß-Catenin/TCF Signaling:
by Eric Röttinger, Paul Dahlin, Mark Q. Martindale

Understanding the functional relationship between intracellular factors and extracellular signals is required for reconstructing gene regulatory networks (GRN) involved in complex biological processes. One of the best-studied bilaterian GRNs describes endomesoderm specification and predicts that both mesoderm and endoderm arose from a common GRN early in animal evolution. Compelling molecular, genomic, developmental, and evolutionary evidence supports the hypothesis that the bifunctional gastrodermis of the cnidarian-bilaterian ancestor is derived from the same evolutionary precursor of both endodermal and mesodermal germ layers in all other triploblastic bilaterian animals. We have begun to establish the framework of a provisional cnidarian “endomesodermal” gene regulatory network in the sea anemone, Nematostella vectensis, by using a genome-wide microarray analysis on embryos in which the canonical Wnt/ß-catenin pathway was ectopically targeted for activation by two distinct pharmaceutical agents (lithium chloride and 1-azakenpaullone) to identify potential targets of endomesoderm specification. We characterized 51 endomesodermally expressed transcription factors and signaling molecule genes (including 18 newly identified) with fine-scale temporal (qPCR) and spatial (in situ) analysis to define distinct co-expression domains within the animal plate of the embryo and clustered genes based on their earliest zygotic expression. Finally, we determined the input of the canonical Wnt/ß-catenin pathway into the cnidarian endomesodermal GRN using morpholino and mRNA overexpression experiments to show that NvTcf/canonical Wnt signaling is required to pattern both the future endomesodermal and ectodermal domains prior to gastrulation, and that both BMP and FGF (but not Notch) pathways play important roles in germ layer specification in this animal. We show both evolutionary conserved as well as profound differences in endomesodermal GRN structure compared to bilaterians that may provide fundamental insight into how GRN subcircuits have been adopted, rewired, or co-opted in various animal lineages that give rise to specialized endomesodermal cell types.

An Excess of Gene Expression Divergence on the X Chromosome in Drosophila Embryos: Implications for the Faster-X Hypothesis

An Excess of Gene Expression Divergence on the X Chromosome in Drosophila Embryos: Implications for the Faster-X Hypothesis:
by Melek A. Kayserili, Dave T. Gerrard, Pavel Tomancak, Alex T. Kalinka

The X chromosome is present as a single copy in the heterogametic sex, and this hemizygosity is expected to drive unusual patterns of evolution on the X relative to the autosomes. For example, the hemizgosity of the X may lead to a lower chromosomal effective population size compared to the autosomes, suggesting that the X might be more strongly affected by genetic drift. However, the X may also experience stronger positive selection than the autosomes, because recessive beneficial mutations will be more visible to selection on the X where they will spend less time being masked by the dominant, less beneficial allele—a proposal known as the faster-X hypothesis. Thus, empirical studies demonstrating increased genetic divergence on the X chromosome could be indicative of either adaptive or non-adaptive evolution. We measured gene expression in Drosophila species and in D. melanogaster inbred strains for both embryos and adults. In the embryos we found that expression divergence is on average more than 20% higher for genes on the X chromosome relative to the autosomes; but in contrast, in the inbred strains, gene expression variation is significantly lower on the X chromosome. Furthermore, expression divergence of genes on Muller's D element is significantly greater along the branch leading to the obscura sub-group, in which this element segregates as a neo-X chromosome. In the adults, divergence is greatest on the X chromosome for males, but not for females, yet in both sexes inbred strains harbour the lowest level of gene expression variation on the X chromosome. We consider different explanations for our results and conclude that they are most consistent within the framework of the faster-X hypothesis.

Mutation spectrum of Drosophila CNVs revealed by breakpoint sequencing

Mutation spectrum of Drosophila CNVs revealed by breakpoint sequencing: Background:
The detailed study of breakpoints associated with Copy Number Variants (CNVs) can elucidate the mutational mechanisms that generate them and the comparison of breakpoints across species can highlight differences in genomic architecture that may lead to lineage-specific differences in patterns of CNVs. Here, we provide a detailed analysis of Drosophila CNV breakpoints and contrast it with similar analyses recently carried out for the human genome.
Results:
By applying split-read methods to a total of 10x coverage of 454 shotgun sequence across 9 lines of D. melanogaster and by re-examining a previously published dataset of CNVs detected using tiling arrays, we identified the precise breakpoints of more than 600 insertions, deletions and duplications. Contrasting these CNVs with those found in humans showed that in both taxa CNV breakpoints fall into three classes: blunt breakpoints; simple breakpoints associated with microhomology; and breakpoints with additional nucleotides inserted/deleted and no microhomology. In both taxa CNV breakpoints are enriched with non-B DNA sequence structures, which may impair DNA replication and/or repair. However, in contrast to human genomes, Non-Allelic Homologous-Recombination (NAHR) plays a negligible role in CNV formation in Drosophila. In flies, non-homologous repair mechanisms are responsible for simple, recurrent and complex CNVs, including insertions of de novo sequence as large as 60 bp.
Conclusions:
Humans and Drosophila differ considerably in the importance of homology-based mechanisms for the formation of CNVs, likely as a consequence of the differences in the abundance and distribution of both segmental duplications and transposable elements between the two genomes.

Function of genomic regions bound by TFs in vivo [Biochemistry]

Function of genomic regions bound by TFs in vivo [Biochemistry]: In animals, each sequence-specific transcription factor typically binds to thousands of genomic regions in vivo. Our previous studies of 20 transcription factors show that most genomic regions bound at high levels in Drosophila blastoderm embryos are known or probable functional targets, but genomic regions occupied only at low levels have...

Dynamic and Differential Regulation of Stem Cell Factor FoxD3 in the Neural Crest Is Encrypted in the Genome

Dynamic and Differential Regulation of Stem Cell Factor FoxD3 in the Neural Crest Is Encrypted in the Genome:
by Marcos S. Simões-Costa, Sonja J. McKeown, Joanne Tan-Cabugao, Tatjana Sauka-Spengler, Marianne E. Bronner

The critical stem cell transcription factor FoxD3 is expressed by the premigratory and migrating neural crest, an embryonic stem cell population that forms diverse derivatives. Despite its important role in development and stem cell biology, little is known about what mediates FoxD3 activity in these cells. We have uncovered two FoxD3 enhancers, NC1 and NC2, that drive reporter expression in spatially and temporally distinct manners. Whereas NC1 activity recapitulates initial FoxD3 expression in the cranial neural crest, NC2 activity recapitulates initial FoxD3 expression at vagal/trunk levels while appearing only later in migrating cranial crest. Detailed mutational analysis, in vivo chromatin immunoprecipitation, and morpholino knock-downs reveal that transcription factors Pax7 and Msx1/2 cooperate with the neural crest specifier gene, Ets1, to bind to the cranial NC1 regulatory element. However, at vagal/trunk levels, they function together with the neural plate border gene, Zic1, which directly binds to the NC2 enhancer. These results reveal dynamic and differential regulation of FoxD3 in distinct neural crest subpopulations, suggesting that heterogeneity is encrypted at the regulatory level. Isolation of neural crest enhancers not only allows establishment of direct regulatory connections underlying neural crest formation, but also provides valuable tools for tissue specific manipulation and investigation of neural crest cell identity in amniotes.

Genome-Wide Fine-Scale Recombination Rate Variation in Drosophila melanogaster

Genome-Wide Fine-Scale Recombination Rate Variation in Drosophila melanogaster:
by Andrew H. Chan, Paul A. Jenkins, Yun S. Song

Estimating fine-scale recombination maps of Drosophila from population genomic data is a challenging problem, in particular because of the high background recombination rate. In this paper, a new computational method is developed to address this challenge. Through an extensive simulation study, it is demonstrated that the method allows more accurate inference, and exhibits greater robustness to the effects of natural selection and noise, compared to a well-used previous method developed for studying fine-scale recombination rate variation in the human genome. As an application, a genome-wide analysis of genetic variation data is performed for two Drosophila melanogaster populations, one from North America (Raleigh, USA) and the other from Africa (Gikongoro, Rwanda). It is shown that fine-scale recombination rate variation is widespread throughout the D. melanogaster genome, across all chromosomes and in both populations. At the fine-scale, a conservative, systematic search for evidence of recombination hotspots suggests the existence of a handful of putative hotspots each with at least a tenfold increase in intensity over the background rate. A wavelet analysis is carried out to compare the estimated recombination maps in the two populations and to quantify the extent to which recombination rates are conserved. In general, similarity is observed at very broad scales, but substantial differences are seen at fine scales. The average recombination rate of the X chromosome appears to be higher than that of the autosomes in both populations, and this pattern is much more pronounced in the African population than the North American population. The correlation between various genomic features—including recombination rates, diversity, divergence, GC content, gene content, and sequence quality—is examined using the wavelet analysis, and it is shown that the most notable difference between D. melanogaster and humans is in the correlation between recombination and diversity.

Long-Range Regulatory Polymorphisms Affecting a GABA Receptor Constitute a Quantitative Trait Locus (QTL) for Social Behavior in Caenorhabditis elegans

Long-Range Regulatory Polymorphisms Affecting a GABA Receptor Constitute a Quantitative Trait Locus (QTL) for Social Behavior in Caenorhabditis elegans:
by Andres Bendesky, Jason Pitts, Matthew V. Rockman, William C. Chen, Man-Wah Tan, Leonid Kruglyak, Cornelia I. Bargmann

Aggregation is a social behavior that varies between and within species, providing a model to study the genetic basis of behavioral diversity. In the nematode Caenorhabditis elegans, aggregation is regulated by environmental context and by two neuromodulatory pathways, one dependent on the neuropeptide receptor NPR-1 and one dependent on the TGF-β family protein DAF-7. To gain further insight into the genetic regulation of aggregation, we characterize natural variation underlying behavioral differences between two wild-type C. elegans strains, N2 and CB4856. Using quantitative genetic techniques, including a survey of chromosome substitution strains and QTL analysis of recombinant inbred lines, we identify three new QTLs affecting aggregation in addition to the two known N2 mutations in npr-1 and glb-5. Fine-mapping with near-isogenic lines localized one QTL, accounting for 5%–8% of the behavioral variance between N2 and CB4856, 3′ to the transcript of the GABA neurotransmitter receptor gene exp-1. Quantitative complementation tests demonstrated that this QTL affects exp-1, identifying exp-1 and GABA signaling as new regulators of aggregation. exp-1 interacts genetically with the daf-7 TGF-β pathway, which integrates food availability and population density, and exp-1 mutations affect the level of daf-7 expression. Our results add to growing evidence that genetic variation affecting neurotransmitter receptor genes is a source of natural behavioral variation.

A Polycomb Group Protein Is Retained at Specific Sites on Chromatin in Mitosis

A Polycomb Group Protein Is Retained at Specific Sites on Chromatin in Mitosis:
by Nicole E. Follmer, Ajazul H. Wani, Nicole J. Francis

Epigenetic regulation of gene expression, including by Polycomb Group (PcG) proteins, may depend on heritable chromatin states, but how these states can be propagated through mitosis is unclear. Using immunofluorescence and biochemical fractionation, we find PcG proteins associated with mitotic chromosomes in Drosophila S2 cells. Genome-wide sequencing of chromatin immunoprecipitations (ChIP–SEQ) from mitotic cells indicates that Posterior Sex Combs (PSC) is not present at well-characterized PcG targets including Hox genes in mitosis, but does remain at a subset of interphase sites. Many of these persistent sites overlap with chromatin domain borders described by Sexton et al. (2012), which are genomic regions characterized by low levels of long range contacts. Persistent PSC binding sites flank both Hox gene clusters. We hypothesize that disruption of long-range chromatin contacts in mitosis contributes to PcG protein release from most sites, while persistent binding at sites with minimal long-range contacts may nucleate re-establishment of PcG binding and chromosome organization after mitosis.

Factorbook.org: a Wiki-based database for transcription factor-binding data generated by the ENCODE consortium

Factorbook.org: a Wiki-based database for transcription factor-binding data generated by the ENCODE consortium:
The Encyclopedia of DNA Elements (ENCODE) consortium aims to identify all functional elements in the human genome including transcripts, transcriptional regulatory regions, along with their chromatin states and DNA methylation patterns. The ENCODE project generates data utilizing a variety of techniques that can enrich for regulatory regions, such as chromatin immunoprecipitation (ChIP), micrococcal nuclease (MNase) digestion and DNase I digestion, followed by deeply sequencing the resulting DNA. As part of the ENCODE project, we have developed a Web-accessible repository accessible at http://factorbook.org. In Wiki format, factorbook is a transcription factor (TF)-centric repository of all ENCODE ChIP-seq datasets on TF-binding regions, as well as the rich analysis results of these data. In the first release, factorbook contains 457 ChIP-seq datasets on 119 TFs in a number of human cell lines, the average profiles of histone modifications and nucleosome positioning around the TF-binding regions, sequence motifs enriched in the regions and the distance and orientation preferences between motif sites.

Genome-wide identification of enhancers in skeletal muscle: the role of MyoD1 [Research Papers]

Genome-wide identification of enhancers in skeletal muscle: the role of MyoD1 [Research Papers]:
To identify the compendium of distal regulatory elements that govern myogenic differentiation, we generated chromatin state maps based on histone modifications and recruitment of factors that typify enhancers in myoblasts and myotubes. We found a striking concordance between the locations of these newly defined enhancers, MyoD1-binding events, and noncoding RNA transcripts. These enhancers recruit several sequence-specific transcription factors in a spatially constrained manner around MyoD1-binding sites. Remarkably, MyoD1-null myoblasts show a wholesale loss of recruitment of these factors as well as diminished monomethylation of H3K4 (H3K4me1) and acetylation of H3K27 (H3K27ac) and reduced recruitment of Set7, an H3K4 monomethylase. Surprisingly, we found that H3K4me1, but not H3K27ac, could be restored by re-expression of MyoD1 in MyoD1^–/– myoblasts, although re-expression of this factor in MyoD1-null myotubes restored both histone modifications. Our studies identified a role for MyoD1 in condition-specific enhancer assembly through recruitment of transcription factors and histone-modifying enzymes that shape muscle differentiation.

Bioelectric signaling regulates head and organ size during planarian regeneration [DEVELOPMENT AND STEM CELLS]

Bioelectric signaling regulates head and organ size during planarian regeneration [DEVELOPMENT AND STEM CELLS]: Wendy Scott Beane, Junji Morokuma, Joan M. Lemire, and Michael Levin


A main goal of regenerative medicine is to replace lost or damaged tissues and organs with functional parts of the correct size and shape. But the proliferation of new cells is not sufficient; we will also need to understand how the scale and ultimate form of newly produced tissues are determined. Using the planarian model system, we report that membrane voltage-dependent bioelectric signaling determines both head size and organ scaling during regeneration. RNA interference of the H⁺,K⁺-ATPase ion pump results in membrane hyperpolarization, which has no effect on the amount of new tissue (blastema) that is regenerated yet produces regenerates with tiny ‘shrunken’ heads and proportionally oversized pharynges. Our data show that this disproportionality results from a lack of the apoptosis required to adjust head and organ size and placement, highlighting apoptotic remodeling as the link between bioelectric signaling and the establishment of organ size during regeneration.

Isolated Tombstone Stands Among Chinese Skyscrapers

Isolated Tombstone Stands Among Chinese Skyscrapers:
In a bizarre dispute, a skyscraper has been built around a tombstone in the city of Taiyuan, the capital of Shanxi province in China.

Building developers bought a cemetery with an eye to building a series of skyscrapers on the land. Prior to construction, locals were paid to relocate the graves, yet one family refused the proposed terms, forcing developers to build around the landmass.

Integrative analysis of Drosophila myogenesis [Systems Biology]

Integrative analysis of Drosophila myogenesis [Systems Biology]: Contemporary high-throughput technologies permit the rapid identification of transcription factor (TF) target genes on a genome-wide scale, yet the functional significance of TFs requires knowledge of target gene expression patterns, cooperating TFs, and cis-regulatory element (CRE) structures. Here we investigated the myogenic regulatory network downstream of the Drosophila zinc finger...

CNVs and Hox gene regulation [Genetics]

CNVs and Hox gene regulation [Genetics]: Copy number variations are genomic structural variants that are frequently associated with human diseases. Among these copy number variations, duplications of DNA segments are often assumed to lead to dosage effects by increasing the copy number of either genes or their regulatory elements. We produced a series of large targeted...

dsx in horned beetle development and evolution [Evolution]

dsx in horned beetle development and evolution [Evolution]: Sex-specific trait expression is frequently associated with highly variable, condition-dependent expression within sexes and rapid divergence among closely related species. Horned beetles are an excellent example for studying the molecular basis of these phenomena because horn morphology varies markedly among species, between sexes, and among alternative, nutritionally-cued morphs within sexes....

Hoxd13 Contribution to the Evolution of Vertebrate Appendages

Hoxd13 Contribution to the Evolution of Vertebrate Appendages: Renata Freitas, Carlos Gómez-Marín, Jonathan Mark Wilson, Fernando Casares, José Luis Gómez-Skarmeta. Fossil data suggest that limbs evolved from fish fins by sequential elaboration of their distal endoskeleton, giving rise to the autopod close to the tetrapod origin. This elaboration may have occ....

Predicting Spatial and Temporal Gene Expression Using an Integrative Model of Transcription Factor Occupancy and Chromatin State

Predicting Spatial and Temporal Gene Expression Using an Integrative Model of Transcription Factor Occupancy and Chromatin State:
by Bartek Wilczynski, Ya-Hsin Liu, Zhen Xuan Yeo, Eileen E. M. Furlong

Precise patterns of spatial and temporal gene expression are central to metazoan complexity and act as a driving force for embryonic development. While there has been substantial progress in dissecting and predicting cis-regulatory activity, our understanding of how information from multiple enhancer elements converge to regulate a gene's expression remains elusive. This is in large part due to the number of different biological processes involved in mediating regulation as well as limited availability of experimental measurements for many of them. Here, we used a Bayesian approach to model diverse experimental regulatory data, leading to accurate predictions of both spatial and temporal aspects of gene expression. We integrated whole-embryo information on transcription factor recruitment to multiple cis-regulatory modules, insulator binding and histone modification status in the vicinity of individual gene loci, at a genome-wide scale during Drosophila development. The model uses Bayesian networks to represent the relation between transcription factor occupancy and enhancer activity in specific tissues and stages. All parameters are optimized in an Expectation Maximization procedure providing a model capable of predicting tissue- and stage-specific activity of new, previously unassayed genes. Performing the optimization with subsets of input data demonstrated that neither enhancer occupancy nor chromatin state alone can explain all gene expression patterns, but taken together allow for accurate predictions of spatio-temporal activity. Model predictions were validated using the expression patterns of more than 600 genes recently made available by the BDGP consortium, demonstrating an average 15-fold enrichment of genes expressed in the predicted tissue over a naïve model. We further validated the model by experimentally testing the expression of 20 predicted target genes of unknown expression, resulting in an accuracy of 95% for temporal predictions and 50% for spatial. While this is, to our knowledge, the first genome-wide approach to predict tissue-specific gene expression in metazoan development, our results suggest that integrative models of this type will become more prevalent in the future.

Identification of a rudimentary neural crest in a non-vertebrate chordate

Identification of a rudimentary neural crest in a non-vertebrate chordate:
Identification of a rudimentary neural crest in a non-vertebrate chordate

Nature 492, 7427 (2012). doi:10.1038/nature11589

Authors: Philip Barron Abitua, Eileen Wagner, Ignacio A. Navarrete & Michael Levine
Neural crest arises at the neural plate border, expresses a core set of regulatory genes and produces a diverse array of cell types, including ectomesenchyme derivatives that elaborate the vertebrate head. The evolution of neural crest has been proposed to be a key event leading to the appearance of new cell types that fostered the transition from filter feeding to active predation in ancestral vertebrates. However, the origin of neural crest remains controversial, as homologous cell types have not been unambiguously identified in non-vertebrate chordates. Here we show that the tunicate Ciona intestinalis possesses a cephalic melanocyte lineage (a9.49) similar to neural crest that can be reprogrammed into migrating ‘ectomesenchyme’ by the targeted misexpression of Twist (also known as twist-like 2). Our results suggest that the neural crest melanocyte regulatory network pre-dated the divergence of tunicates and vertebrates. We propose that the co-option of mesenchyme determinants, such as Twist, into the neural plate ectoderm was crucial to the emergence of the vertebrate ‘new head’.

Extensive compensatory cis-trans regulation in the evolution of mouse gene expression [RESEARCH]

Extensive compensatory cis-trans regulation in the evolution of mouse gene expression [RESEARCH]:

Gene expression levels are thought to diverge primarily via regulatory mutations in trans within species, and in cis between species. To test this hypothesis in mammals we used RNA-sequencing to measure gene expression divergence between C57BL/6J and CAST/EiJ mouse strains and allele-specific expression in their F1 progeny. We identified 535 genes with parent-of-origin specific expression patterns, although few of these showed full allelic silencing. This suggests that the number of imprinted genes in a typical mouse somatic tissue is relatively small. In the set of nonimprinted genes, 32% showed evidence of divergent expression between the two strains. Of these, 2% could be attributed purely to variants acting in trans, while 43% were attributable only to variants acting in cis. The genes with expression divergence driven by changes in trans showed significantly higher sequence constraint than genes where the divergence was explained by variants acting in cis. The remaining genes with divergent patterns of expression (55%) were regulated by a combination of variants acting in cis and variants acting in trans. Intriguingly, the changes in expression induced by the cis and trans variants were in opposite directions more frequently than expected by chance, implying that compensatory regulation to stabilize gene expression levels is widespread. We propose that expression levels of genes regulated by this mechanism are fine-tuned by cis variants that arise following regulatory changes in trans, suggesting that many cis variants are not the primary targets of natural selection.

Extensive conservation of ancient microsynteny across metazoans due to cis-regulatory constraints [RESEARCH]

Extensive conservation of ancient microsynteny across metazoans due to cis-regulatory constraints [RESEARCH]:
The order of genes in eukaryotic genomes has generally been assumed to be neutral, since gene order is largely scrambled over evolutionary time. Only a handful of exceptional examples are known, typically involving deeply conserved clusters of tandemly duplicated genes (e.g., Hox genes and histones). Here we report the first systematic survey of microsynteny conservation across metazoans, utilizing 17 genome sequences. We identified nearly 600 pairs of unrelated genes that have remained tightly physically linked in diverse lineages across over 600 million years of evolution. Integrating sequence conservation, gene expression data, gene function, epigenetic marks, and other genomic features, we provide extensive evidence that many conserved ancient linkages involve (1) the coordinated transcription of neighboring genes, or (2) genomic regulatory blocks (GRBs) in which transcriptional enhancers controlling developmental genes are contained within nearby bystander genes. In addition, we generated ChIP-seq data for key histone modifications in zebrafish embryos, which provided further evidence of putative GRBs in embryonic development. Finally, using chromosome conformation capture (3C) assays and stable transgenic experiments, we demonstrate that enhancers within bystander genes drive the expression of genes such as Otx and Islet, critical regulators of central nervous system development across bilaterians. These results suggest that ancient genomic functional associations are far more common than previously thought—involving ~12% of the ancestral bilaterian genome—and that cis-regulatory constraints are crucial in determining metazoan genome architecture.

Tissue-Specific Regulation of Chromatin Insulator Function

Tissue-Specific Regulation of Chromatin Insulator Function:
by Leah H. Matzat, Ryan K. Dale, Nellie Moshkovich, Elissa P. Lei

Chromatin insulators organize the genome into distinct transcriptional domains and contribute to cell type–specific chromatin organization. However, factors regulating tissue-specific insulator function have not yet been discovered. Here we identify the RNA recognition motif-containing protein Shep as a direct interactor of two individual components of the gypsy insulator complex in Drosophila. Mutation of shep improves gypsy-dependent enhancer blocking, indicating a role as a negative regulator of insulator activity. Unlike ubiquitously expressed core gypsy insulator proteins, Shep is highly expressed in the central nervous system (CNS) with lower expression in other tissues. We developed a novel, quantitative tissue-specific barrier assay to demonstrate that Shep functions as a negative regulator of insulator activity in the CNS but not in muscle tissue. Additionally, mutation of shep alters insulator complex nuclear localization in the CNS but has no effect in other tissues. Consistent with negative regulatory activity, ChIP–seq analysis of Shep in a CNS-derived cell line indicates substantial genome-wide colocalization with a single gypsy insulator component but limited overlap with intact insulator complexes. Taken together, these data reveal a novel, tissue-specific mode of regulation of a chromatin insulator.