Analysing grouping of nucleotides in DNA sequences using lumped processes constructed from Markov chains

The most commonly used models for analysing local dependencies in DNA sequences are (high-order) Markov chains. Incorporating knowledge relative to the possible grouping of the nucleotides enables to define dedicated sub-classes of Markov chains. The problem of formulating lumpability hypotheses for a Markov chain is therefore addressed. In the classical approach to lumpability, this problem can be formulated as the determination of an appropriate state space (smaller than the original state space) such that the lumped chain defined on this state space retains the Markov property. We propose a different perspective on lumpability where the state space is fixed and the partitioning of this state space is represented by a one-to-many probabilistic function within a two-level stochastic process. Three nested classes of lumped processes can be defined in this way as sub-classes of first-order Markov chains. These lumped processes enable parsimonious reparameterizations of Markov chains that help to reveal relevant partitions of the state space. Characterizations of the lumped processes on the original transition probability matrix are derived. Different model selection methods relying either on hypothesis testing or on penalized log-likelihood criteria are presented as well as extensions to lumped processes constructed from high-order Markov chains. The relevance of the proposed approach to lumpability is illustrated by the analysis of DNA sequences. In particular, the use of lumped processes enables to highlight differences between intronic sequences and gene untranslated region sequences

A genetically encoded reporter of synaptic activity in vivo

To image synaptic activity within neural circuits, we tethered the genetically encoded calcium indicator (GECI) GCaMP2 to synaptic vesicles by fusion to synaptophysin. The resulting reporter, SyGCaMP2, detected the electrical activity of neurons with two advantages over existing cytoplasmic GECIs: it identified the locations of synapses and had a linear response over a wider range of spike frequencies. Simulations and experimental measurements indicated that linearity arises because SyGCaMP2 samples the brief calcium transient passing through the presynaptic compartment close to voltage-sensitive calcium channels rather than changes in bulk calcium concentration. In vivo imaging in zebrafish demonstrated that SyGCaMP2 can assess electrical activity in conventional synapses of spiking neurons in the optic tectum and graded voltage signals transmitted by ribbon synapses of retinal bipolar cells. Localizing a GECI to synaptic terminals provides a strategy for monitoring activity across large groups of neurons at the level of individual synapses

Prediction of Transcriptional Terminators in Bacillus subtilis and Related Species

In prokaryotes, genes belonging to the same operon are transcribed in a single mRNA molecule. Transcription starts as the RNA polymerase binds to the promoter and continues until it reaches a transcriptional terminator. Some terminators rely on the presence of the Rho protein, whereas others function independently of Rho. Such Rho-independent terminators consist of an inverted repeat followed by a stretch of thymine residues, allowing us to predict their presence directly from the DNA sequence. Unlike in Escherichia coli, the Rho protein is dispensable in Bacillus subtilis, suggesting a limited role for Rho-dependent termination in this organism and possibly in other Firmicutes. We analyzed 463 experimentally known terminating sequences in B. subtilis and found a decision rule to distinguish Rho-independent transcriptional terminators from non-terminating sequences. The decision rule allowed us to find the boundaries of operons in B. subtilis with a sensitivity and specificity of about 94%. Using the same decision rule, we found an average sensitivity of 94% for 57 bacteria belonging to the Firmicutes phylum, and a considerably lower sensitivity for other bacteria. Our analysis shows that Rho-independent termination is dominant for Firmicutes in general, and that the properties of the transcriptional terminators are conserved. Terminator prediction can be used to reliably predict the operon structure in these organisms, even in the absence of experimentally known operons. Genome-wide predictions of Rho-independent terminators for the 57 Firmicutes are available in the Supporting Information section

Modeling Neurodegeneration in Zebrafish

The zebrafish, Danio rerio, has been established as an excellent vertebrate model for the study of developmental biology and gene function. It also has proven to be a valuable model to study human diseases. Here, we reviewed recent publications using zebrafish to study the pathology of human neurodegenerative diseases including Parkinson’s, Huntington’s, and Alzheimer’s. These studies indicate that zebrafish genes and their human homologues have conserved functions with respect to the etiology of neurodegenerative diseases. The characteristics of the zebrafish and the experimental approaches to which it is amenable make this species a useful complement to other animal models for the study of pathologic mechanisms of neurodegenerative diseases and for the screening of compounds with therapeutic potential

Identification of Early Requirements for Preplacodal Ectoderm and Sensory Organ Development

Preplacodal ectoderm arises near the end of gastrulation as a narrow band of cells surrounding the anterior neural plate. This domain later resolves into discrete cranial placodes that, together with neural crest, produce paired sensory structures of the head. Unlike the better-characterized neural crest, little is known about early regulation of preplacodal development. Classical models of ectodermal patterning posit that preplacodal identity is specified by readout of a discrete level of Bmp signaling along a DV gradient. More recent studies indicate that Bmp-antagonists are critical for promoting preplacodal development. However, it is unclear whether Bmp-antagonists establish the proper level of Bmp signaling within a morphogen gradient or, alternatively, block Bmp altogether. To begin addressing these issues, we treated zebrafish embryos with a pharmacological inhibitor of Bmp, sometimes combined with heat shock-induction of Chordin and dominant-negative Bmp receptor, to fully block Bmp signaling at various developmental stages. We find that preplacodal development occurs in two phases with opposing Bmp requirements. Initially, Bmp is required before gastrulation to co-induce four transcription factors, Tfap2a, Tfap2c, Foxi1, and Gata3, which establish preplacodal competence throughout the nonneural ectoderm. Subsequently, Bmp must be fully blocked in late gastrulation by dorsally expressed Bmp-antagonists, together with dorsally expressed Fgf and Pdgf, to specify preplacodal identity within competent cells abutting the neural plate. Localized ventral misexpression of Fgf8 and Chordin can activate ectopic preplacodal development anywhere within the zone of competence, whereas dorsal misexpression of one or more competence factors can activate ectopic preplacodal development in the neural plate. Conversely, morpholino-knockdown of competence factors specifically ablates preplacodal development. Our work supports a relatively simple two-step model that traces regulation of preplacodal development to late blastula stage, resolves two distinct phases of Bmp dependence, and identifies the main factors required for preplacodal competence and specification

Anaplastic Lymphoma Kinase Is Required for Neurogenesis in the Developing Central Nervous System of Zebrafish

10.1371/journal.pone.0063757PLoS ONE85

Expression of zebrafish pax6b in pancreas is regulated by two enhancers containing highly conserved cis-elements bound by PDX1, PBX and PREP factors

BACKGROUND: PAX6 is a transcription factor playing a crucial role in the development of the eye and in the differentiation of the pancreatic endocrine cells as well as of enteroendocrine cells. Studies on the mouse Pax6 gene have shown that sequences upstream from the P0 promoter are required for expression in the lens and the pancreas; but there remain discrepancies regarding the precise location of the pancreatic regulatory elements. RESULTS: Due to genome duplication in the evolution of ray-finned fishes, zebrafish has two pax6 genes, pax6a and pax6b. While both zebrafish pax6 genes are expressed in the developing eye and nervous system, only pax6b is expressed in the endocrine cells of the pancreas. To investigate the cause of this differential expression, we used a combination of in silico, in vivo and in vitro approaches. We show that the pax6b P0 promoter targets expression to endocrine pancreatic cells and also to enteroendocrine cells, retinal neurons and the telencephalon of transgenic zebrafish. Deletion analyses indicate that strong pancreatic expression of the pax6b gene relies on the combined action of two conserved regulatory enhancers, called regions A and C. By means of gel shift assays, we detected binding of the homeoproteins PDX1, PBX and PREP to several cis-elements of these regions. In constrast, regions A and C of the zebrafish pax6a gene are not active in the pancreas, this difference being attributable to sequence divergences within two cis-elements binding the pancreatic homeoprotein PDX1. CONCLUSION: Our data indicate a conserved role of enhancers A and C in the pancreatic expression of pax6b and emphasize the importance of the homeoproteins PBX and PREP cooperating with PDX1, in activating pax6b expression in endocrine pancreatic cells. This study also provides a striking example of how adaptative evolution of gene regulatory sequences upon gene duplication progressively leads to subfunctionalization of the paralogous gene pair

Replication Fork Polarity Gradients Revealed by Megabase-Sized U-Shaped Replication Timing Domains in Human Cell Lines

In higher eukaryotes, replication program specification in different cell types remains to be fully understood. We show for seven human cell lines that about half of the genome is divided in domains that display a characteristic U-shaped replication timing profile with early initiation zones at borders and late replication at centers. Significant overlap is observed between U-domains of different cell lines and also with germline replication domains exhibiting a N-shaped nucleotide compositional skew. From the demonstration that the average fork polarity is directly reflected by both the compositional skew and the derivative of the replication timing profile, we argue that the fact that this derivative displays a N-shape in U-domains sustains the existence of large-scale gradients of replication fork polarity in somatic and germline cells. Analysis of chromatin interaction (Hi-C) and chromatin marker data reveals that U-domains correspond to high-order chromatin structural units. We discuss possible models for replication origin activation within U/N-domains. The compartmentalization of the genome into replication U/N-domains provides new insights on the organization of the replication program in the human genome

De Novo Genesis of Enhancers in Vertebrates

Whole genome duplication in teleost fish reveals that a few changes in non-regulatory genomic sequences are a source for generating new enhancers

Evidence for Sequential and Increasing Activation of Replication Origins along Replication Timing Gradients in the Human Genome

Genome-wide replication timing studies have suggested that mammalian chromosomes consist of megabase-scale domains of coordinated origin firing separated by large originless transition regions. Here, we report a quantitative genome-wide analysis of DNA replication kinetics in several human cell types that contradicts this view. DNA combing in HeLa cells sorted into four temporal compartments of S phase shows that replication origins are spaced at 40 kb intervals and fire as small clusters whose synchrony increases during S phase and that replication fork velocity (mean 0.7 kb/min, maximum 2.0 kb/min) remains constant and narrowly distributed through S phase. However, multi-scale analysis of a genome-wide replication timing profile shows a broad distribution of replication timing gradients with practically no regions larger than 100 kb replicating at less than 2 kb/min. Therefore, HeLa cells lack large regions of unidirectional fork progression. Temporal transition regions are replicated by sequential activation of origins at a rate that increases during S phase and replication timing gradients are set by the delay and the spacing between successive origin firings rather than by the velocity of single forks. Activation of internal origins in a specific temporal transition region is directly demonstrated by DNA combing of the IGH locus in HeLa cells. Analysis of published origin maps in HeLa cells and published replication timing and DNA combing data in several other cell types corroborate these findings, with the interesting exception of embryonic stem cells where regions of unidirectional fork progression seem more abundant. These results can be explained if origins fire independently of each other but under the control of long-range chromatin structure, or if replication forks progressing from early origins stimulate initiation in nearby unreplicated DNA. These findings shed a new light on the replication timing program of mammalian genomes and provide a general model for their replication kinetics