An Integrated Model of Multiple-Condition ChIP-Seq Data Reveals Predeterminants of Cdx2 Binding

Regulatory proteins can bind to different sets of genomic targets in various cell types or conditions. To reliably characterize such condition-specific regulatory binding we introduce MultiGPS, an integrated machine learning approach for the analysis of multiple related ChIP-seq experiments. MultiGPS is based on a generalized Expectation Maximization framework that shares information across multiple experiments for binding event discovery. We demonstrate that our framework enables the simultaneous modeling of sparse condition-specific binding changes, sequence dependence, and replicate-specific noise sources. MultiGPS encourages consistency in reported binding event locations across multiple-condition ChIP-seq datasets and provides accurate estimation of ChIP enrichment levels at each event. MultiGPS's multi-experiment modeling approach thus provides a reliable platform for detecting differential binding enrichment across experimental conditions. We demonstrate the advantages of MultiGPS with an analysis of Cdx2 binding in three distinct developmental contexts. By accurately characterizing condition-specific Cdx2 binding, MultiGPS enables novel insight into the mechanistic basis of Cdx2 site selectivity. Specifically, the condition-specific Cdx2 sites characterized by MultiGPS are highly associated with pre-existing genomic context, suggesting that such sites are pre-determined by cell-specific regulatory architecture. However, MultiGPS-defined condition-independent sites are not predicted by pre-existing regulatory signals, suggesting that Cdx2 can bind to a subset of locations regardless of genomic environment. A summary of this paper appears in the proceedings of the RECOMB 2014 conference, April 2–5.National Science Foundation (U.S.) (Graduate Research Fellowship under Grant 0645960)National Institutes of Health (U.S.) (grant P01 NS055923)Pennsylvania State University. Center for Eukaryotic Gene Regulatio

Synergistic binding of transcription factors to cell-specific enhancers programs motor neuron identity

Efficient transcriptional programming promises to open new frontiers in regenerative medicine. However, mechanisms by which programming factors transform cell fate are unknown, preventing more rational selection of factors to generate desirable cell types. Three transcription factors, Ngn2, Isl1 and Lhx3, were sufficient to program rapidly and efficiently spinal motor neuron identity when expressed in differentiating mouse embryonic stem cells. Replacement of Lhx3 by Phox2a led to specification of cranial, rather than spinal, motor neurons. Chromatin immunoprecipitation–sequencing analysis of Isl1, Lhx3 and Phox2a binding sites revealed that the two cell fates were programmed by the recruitment of Isl1-Lhx3 and Isl1-Phox2a complexes to distinct genomic locations characterized by a unique grammar of homeodomain binding motifs. Our findings suggest that synergistic interactions among transcription factors determine the specificity of their recruitment to cell type–specific binding sites and illustrate how a single transcription factor can be repurposed to program different cell types.Project ALS FoundationNational Institutes of Health (U.S.) (Grant P01 NS055923

Saltatory remodeling of Hox chromatin in response to rostrocaudal patterning signals

Hox genes controlling motor neuron subtype identity are expressed in rostrocaudal patterns that are spatially and temporally collinear with their chromosomal organization. Here we demonstrate that Hox chromatin is subdivided into discrete domains that are controlled by rostrocaudal patterning signals that trigger rapid, domain-wide clearance of repressive histone H3 Lys27 trimethylation (H3K27me3) polycomb modifications. Treatment of differentiating mouse neural progenitors with retinoic acid leads to activation and binding of retinoic acid receptors (RARs) to the Hox1–Hox5 chromatin domains, which is followed by a rapid domain-wide removal of H3K27me3 and acquisition of cervical spinal identity. Wnt and fibroblast growth factor (FGF) signals induce expression of the Cdx2 transcription factor that binds and clears H3K27me3 from the Hox1–Hox9 chromatin domains, leading to specification of brachial or thoracic spinal identity. We propose that rapid clearance of repressive modifications in response to transient patterning signals encodes global rostrocaudal neural identity and that maintenance of these chromatin domains ensures the transmission of positional identity to postmitotic motor neurons later in development.Leona M. and Harry B. Helmsley Charitable TrustNational Institutes of Health (U.S.) (Grant P01 NS055923)Smith Family Foundatio

Stochastic De-repression of Rhodopsins in Single Photoreceptors of the Fly Retina

The photoreceptors of the Drosophila compound eye are a classical model for studying cell fate specification. Photoreceptors (PRs) are organized in bundles of eight cells with two major types – inner PRs involved in color vision and outer PRs involved in motion detection. In wild type flies, most PRs express a single type of Rhodopsin (Rh): inner PRs express either Rh3, Rh4, Rh5 or Rh6 and outer PRs express Rh1. In outer PRs, the K50 homeodomain protein Dve is a key repressor that acts to ensure exclusive Rh expression. Loss of Dve results in de-repression of Rhodopsins in outer PRs, and leads to a wide distribution of expression levels. To quantify these effects, we introduce an automated image analysis method to measure Rhodopsin levels at the single cell level in 3D confocal stacks. Our sensitive methodology reveals cell-specific differences in Rhodopsin distributions among the outer PRs, observed over a developmental time course. We show that Rhodopsin distributions are consistent with a two-state model of gene expression, in which cells can be in either high or basal states of Rhodopsin production. Our model identifies a significant role of post-transcriptional regulation in establishing the two distinct states. The timescale for interconversion between basal and high states is shown to be on the order of days. Our results indicate that even in the absence of Dve, the Rhodopsin regulatory network can maintain highly stable states. We propose that the role of Dve in outer PRs is to buffer against rare fluctuations in this network

A search for the decay modes B+/- to h+/- tau l

We present a search for the lepton flavor violating decay modes B+/- to h+/- tau l (h= K,pi; l= e,mu) using the BaBar data sample, which corresponds to 472 million BBbar pairs. The search uses events where one B meson is fully reconstructed in one of several hadronic final states. Using the momenta of the reconstructed B, h, and l candidates, we are able to fully determine the tau four-momentum. The resulting tau candidate mass is our main discriminant against combinatorial background. We see no evidence for B+/- to h+/- tau l decays and set a 90% confidence level upper limit on each branching fraction at the level of a few times 10^-5.Comment: 15 pages, 7 figures, submitted to Phys. Rev.

Evidence for an excess of B -> D(*) Tau Nu decays

Based on the full BaBar data sample, we report improved measurements of the ratios R(D(*)) = B(B -> D(*) Tau Nu)/B(B -> D(*) l Nu), where l is either e or mu. These ratios are sensitive to new physics contributions in the form of a charged Higgs boson. We measure R(D) = 0.440 +- 0.058 +- 0.042 and R(D*) = 0.332 +- 0.024 +- 0.018, which exceed the Standard Model expectations by 2.0 sigma and 2.7 sigma, respectively. Taken together, our results disagree with these expectations at the 3.4 sigma level. This excess cannot be explained by a charged Higgs boson in the type II two-Higgs-doublet model. We also report the observation of the decay B -> D Tau Nu, with a significance of 6.8 sigma.Comment: Expanded section on systematics, text corrections, improved the format of Figure 2 and included the effect of the change of the Tau polarization due to the charged Higg

Search for the decay modes D^0 → e^+e^-, D^0 → μ^+μ^-, and D^0 → e^±μ∓

We present searches for the rare decay modes D^0→e^+e^-, D^0→μ^+μ^-, and D^0→e^±μ^∓ in continuum e^+e^-→cc events recorded by the BABAR detector in a data sample that corresponds to an integrated luminosity of 468  fb^(-1). These decays are highly Glashow–Iliopoulos–Maiani suppressed but may be enhanced in several extensions of the standard model. Our observed event yields are consistent with the expected backgrounds. An excess is seen in the D^0→μ^+μ^- channel, although the observed yield is consistent with an upward background fluctuation at the 5% level. Using the Feldman–Cousins method, we set the following 90% confidence level intervals on the branching fractions: B(D^0→e^+e^-)<1.7×10^(-7), B(D^0→μ^+μ^-) within [0.6,8.1]×10^(-7), and B(D^0→e^±μ^∓)<3.3×10^(-7)

Study of the reaction e^{+}e^{-} -->J/psi\pi^{+}\pi^{-} via initial-state radiation at BaBar

We study the process $e^+e^-\to J/\psi\pi^{+}\pi^{-}$ with initial-state-radiation events produced at the PEP-II asymmetric-energy collider. The data were recorded with the BaBar detector at center-of-mass energies 10.58 and 10.54 GeV, and correspond to an integrated luminosity of 454 $\mathrm{fb^{-1}}$. We investigate the $J/\psi \pi^{+}\pi^{-}$ mass distribution in the region from 3.5 to 5.5 $\mathrm{GeV/c^{2}}$. Below 3.7 $\mathrm{GeV/c^{2}}$ the $\psi(2S)$ signal dominates, and above 4 $\mathrm{GeV/c^{2}}$ there is a significant peak due to the Y(4260). A fit to the data in the range 3.74 -- 5.50 $\mathrm{GeV/c^{2}}$ yields a mass value $4244 \pm 5$ (stat) $\pm 4$ (syst)$\mathrm{MeV/c^{2}}$ and a width value $114 ^{+16}_{-15}$ (stat)$\pm 7$(syst)$\mathrm{MeV}$ for this state. We do not confirm the report from the Belle collaboration of a broad structure at 4.01 $\mathrm{GeV/c^{2}}$. In addition, we investigate the $\pi^{+}\pi^{-}$ system which results from Y(4260) decay

Observation and study of baryonic B decays: B -> D(*) p pbar, D(*) p pbar pi, and D(*) p pbar pi pi

We present a study of ten B-meson decays to a D(*), a proton-antiproton pair, and a system of up to two pions using BaBar's data set of 455x10^6 BBbar pairs. Four of the modes (B0bar -> D0 p anti-p, B0bar -> D*0 p anti-p, B0bar -> D+ p anti-p pi-, B0bar -> D*+ p anti-p pi-) are studied with improved statistics compared to previous measurements; six of the modes (B- -> D0 p anti-p pi-, B- -> D*0 p anti-p pi-, B0bar -> D0 p anti-p pi- pi+, B0bar -> D*0 p anti-p pi- pi+, B- -> D+ p anti-p pi- pi-, B- -> D*+ p anti-p pi- pi-) are first observations. The branching fractions for 3- and 5-body decays are suppressed compared to 4-body decays. Kinematic distributions for 3-body decays show non-overlapping threshold enhancements in m(p anti-p) and m(D(*)0 p) in the Dalitz plots. For 4-body decays, m(p pi-) mass projections show a narrow peak with mass and full width of (1497.4 +- 3.0 +- 0.9) MeV/c2, and (47 +- 12 +- 4) MeV/c2, respectively, where the first (second) errors are statistical (systematic). For 5-body decays, mass projections are similar to phase space expectations. All results are preliminary.Comment: 28 pages, 90 postscript figures, submitted to LP0

Search for rare quark-annihilation decays, B --> Ds(*) Phi

We report on searches for B- --> Ds- Phi and B- --> Ds*- Phi. In the context of the Standard Model, these decays are expected to be highly suppressed since they proceed through annihilation of the b and u-bar quarks in the B- meson. Our results are based on 234 million Upsilon(4S) --> B Bbar decays collected with the BABAR detector at SLAC. We find no evidence for these decays, and we set Bayesian 90% confidence level upper limits on the branching fractions BF(B- --> Ds- Phi) Ds*- Phi)<1.2x10^(-5). These results are consistent with Standard Model expectations.Comment: 8 pages, 3 postscript figues, submitted to Phys. Rev. D (Rapid Communications