Transcriptional memory emerges from cooperative histone modifications

Background
Transcriptional regulation in cells makes use of diverse mechanisms to ensure that functional states can be maintained and adapted to variable environments; among them are chromatin-related mechanisms. While mathematical models of transcription factor networks controlling development are well established, models of transcriptional regulation by chromatin states are rather rare despite they appear to be a powerful regulatory mechanism.
Results
We here introduce a mathematical model of transcriptional regulation governed by histone modifications. This model describes binding of protein complexes to chromatin which are capable of reading and writing histone marks. Molecular interactions between these complexes and DNA or histones create a regulatory switch of transcriptional activity possessing a regulatory memory. The regulatory states of the switch depend on the activity of histone (de-) methylases, the structure of the DNA-binding regions of the complexes, and the number of histones contributing to binding. 
We apply our model to transcriptional regulation by trithorax- and polycomb- complex binding. By analyzing data on pluripotent and lineage-committed cells we verify basic model assumptions and provide evidence for a positive effect of the length of the modified regions on the stability of the induced regulatory states and thus on the transcriptional memory.
Conclusions
Our results provide new insights into epigenetic modes of transcriptional regulation. Moreover, they implicate well-founded hypotheses on cooperative histone modifications, proliferation induced epigenetic changes and higher order folding of chromatin which await experimental validation. Our approach represents a basic step towards multi-scale models of transcriptional control during development and lineage specification. 
&#xa

Evolution of the vertebrate parahox clusters

The ParaHox cluster contains three Hox‐related homeobox genes. The evolution of this sister of the Hox‐gene clusters has been studied extensively in metazoans with a focus on its early evolution. Its fate within the vertebrate lineage, and in particular following the teleost‐specific genome duplication, however, has not received much attention. Three of the four human ParaHox loci are linked with PDGFR family tyrosine kinases. We demonstrate that these loci arose as duplications in an ancestral vertebrate and trace the subsequent history of gene losses. Surprisingly, teleost fishes have not expanded their ParaHox repertoire following the teleost‐specific genome duplication, while duplicates of the associated tyrosine kinases have survived, supporting the hypothesis of a large‐scale duplication followed by extensive gene loss

Exclusion of repetitive DNA elements from gnathostome Hox clusters

Despite their homology and analogous function, the Hox gene clusters of vertebrates and invertebrates are subject to different constraints on their structural organization. This is demonstrated by a drastically different distribution of repetitive DNA elements in the Hox cluster regions. While gnathostomes have a strong tendency to exclude repetitive DNA elements from the inside of their Hox clusters, no such trend can be detected in the Hox gene clusters of protostomes. Repeats “invade” the gnathostome Hox clusters from the 5′ and 3′ ends while the core of the clusters remains virtually free of repetitive DNA. This invasion appears to be correlated with relaxed constraints associated with gene loss after cluster duplications

SynBlast: Assisting the Analysis of Conserved Synteny Information

Motivation: In the last years more than 20 vertebrate genomes have been sequenced, and the rate at which genomic DNA information becomes available is rapidly accelerating. Gene duplication and gene loss events inherently limit the accuracy of orthology detection based on sequence similarity alone. Fully automated methods for orthology annotation do exist but often fail to identify individual members in cases of large gene families, or to distinguish missing data from traceable gene losses. This situation can be improved in many cases by including conserved synteny information. Results: Here we present the SynBlast pipeline that is designed to construct and evaluate local synteny information. SynBlast uses the genomic region around a focal reference gene to retrieve candidates for homologous regions from a collection of target genomes and ranks them in accord with the available evidence for homology. The pipeline is intended as a tool to aid high quality manual annotation in particular in those cases where automatic procedures fail. We demonstrate how SynBlast is applied to retrieving orthologous and paralogous clusters using the vertebrate Hox and ParaHox clusters as examples

Independent Hox‐cluster duplications in lampreys

The analysis of the publicly available Hox gene sequences from the sea lamprey Petromyzon marinus provides evidence that the Hox clusters in lampreys and other vertebrate species arose from independent duplications. In particular, our analysis supports the hypothesis that the last common ancestor of agnathans and gnathostomes had only a single Hox cluster which was subsequently duplicated independently in the two lineages

The duplication of the Hox gene clusters in teleost fishes

Higher teleost fishes, including zebrafish and fugu, have duplicated their Hox genes relative to the gene inventory of other gnathostome lineages. The most widely accepted theory contends that the duplicate Hox clusters orginated synchronously during a single genome duplication event in the early history of ray-finned fishes. In this contribution we collect and re-evaluate all publicly available sequence information. In particular, we show that the short Hox gene fragments from published PCR surveys of the killifish Fundulus heteroclitus, the medaka Oryzias latipes and the goldfish Carassius auratus can be used to determine with little ambiguity not only their paralog group but also their membership in a particular cluster. Together with a survey of the genomic sequence data from the pufferfish Tetraodon nigroviridis we show that at least percomorpha, and possibly all eutelosts, share a system of 7 or 8 orthologous Hox gene clusters. There is little doubt about the orthology of the two teleost duplicates of the HoxA and HoxB clusters. A careful analysis of both the coding sequence of Hox genes and of conserved non-coding sequences provides additional support for the “duplication early” hypothesis that the Hox clusters in teleosts are derived from eight ancestral clusters by means of subsequent gene loss; the data remain ambiguous, however, in particular for the HoxC clusters. Assuming the “duplication early” hypothesis we use the new evidence on the Hox gene complements to determine the phylogenetic positions of gene-loss events in the wake of the cluster duplication. Surprisingly, we find that the resolution of redundancy seems to be a slow process that is still ongoing. A few suggestions on which additional sequence data would be most informative for resolving the history of the teleostean Hox genes are discussed

Divergent Evolution of Eukaryotic CC- and A-Adding Enzymes

Synthesis of the CCA end of essential tRNAs is performed either by CCA-adding enzymes or as a collaboration between enzymes restricted to CC- and A-incorporation. While the occurrence of such tRNA nucleotidyltransferases with partial activities seemed to be restricted to Bacteria, the first example of such split CCA-adding activities was reported in Schizosaccharomyces pombe. Here, we demonstrate that the choanoflagellate Salpingoeca rosetta also carries CC- and A-adding enzymes. However, these enzymes have distinct evolutionary origins. Furthermore, the restricted activity of the eukaryotic CC-adding enzymes has evolved in a different way compared to their bacterial counterparts. Yet, the molecular basis is very similar, as highly conserved positions within a catalytically important flexible loop region are missing in the CC-adding enzymes. For both the CC-adding enzymes from S. rosetta as well as S. pombe, introduction of the loop elements from closely related enzymes with full activity was able to restore CCA-addition, corroborating the significance of this loop in the evolution of bacterial as well as eukaryotic tRNA nucleotidyltransferases. Our data demonstrate that partial CC- and A-adding activities in Bacteria and Eukaryotes are based on the same mechanistic principles but, surprisingly, originate from different evolutionary events

Divergence of Conserved Non-Coding Sequences: Rate Estimate and Relative Rate Tests

In many eukaryotic genomes only a small fraction of the DNA codes for proteins, but the non-protein coding DNA harbors important genetic elements directing the development and the physiology of the organisms, like promoters, enhancers, insulators, and micro-RNA genes. The molecular evolution of these genetic elements is difficult to study because their functional significance is hard to deduce from sequence information alone. Here we propose an approach to the study of the rate of evolution of functional non-coding sequences at a macro-evolutionary scale. We identify functionally important non-coding sequences as Conserved Non-Coding Nucleotide (CNCN) sequences from the comparison of two outgroup species. The CNCN sequences so identified are then compared to their homologous sequences in a pair of ingroup species, and we monitor the degree of modification these sequences suffered in the two ingroup lineages. We propose a method to test for rate differences in the modification of CNCN sequences among the two ingroup lineages, as well as a method to estimate their rate of modification. We apply this method to the full sequences of the HoxA clusters from six gnathostome species: a shark, Heterodontus francisci; a basal ray finned fish, Polypterus senegalus; the amphibian, Xenopus tropicalis; as well as three mammalian species, human, rat and mouse. The results show that the evolutionary rate of CNCN sequences is not distinguishable among the three mammalian lineages, while the Xenopus lineage has a significantly increased rate of evolution. Furthermore the estimates of the rate parameters suggest that in the stem lineage of mammals the rate of CNCN sequence evolution was more than twice the rate observed within the placental amniotes clade, suggesting a high rate of evolution of cis-regulatory elements during the origin of amniotes and mammals. We conclude that the proposed methods can be used for testing hypotheses about the rate and pattern of evolution of putative cis-regulatory elements