Methods to study splicing from high-throughput RNA Sequencing data

The development of novel high-throughput sequencing (HTS) methods for RNA (RNA-Seq) has provided a very powerful mean to study splicing under multiple conditions at unprecedented depth. However, the complexity of the information to be analyzed has turned this into a challenging task. In the last few years, a plethora of tools have been developed, allowing researchers to process RNA-Seq data to study the expression of isoforms and splicing events, and their relative changes under different conditions. We provide an overview of the methods available to study splicing from short RNA-Seq data. We group the methods according to the different questions they address: 1) Assignment of the sequencing reads to their likely gene of origin. This is addressed by methods that map reads to the genome and/or to the available gene annotations. 2) Recovering the sequence of splicing events and isoforms. This is addressed by transcript reconstruction and de novo assembly methods. 3) Quantification of events and isoforms. Either after reconstructing transcripts or using an annotation, many methods estimate the expression level or the relative usage of isoforms and/or events. 4) Providing an isoform or event view of differential splicing or expression. These include methods that compare relative event/isoform abundance or isoform expression across two or more conditions. 5) Visualizing splicing regulation. Various tools facilitate the visualization of the RNA-Seq data in the context of alternative splicing. In this review, we do not describe the specific mathematical models behind each method. Our aim is rather to provide an overview that could serve as an entry point for users who need to decide on a suitable tool for a specific analysis. We also attempt to propose a classification of the tools according to the operations they do, to facilitate the comparison and choice of methods.Comment: 31 pages, 1 figure, 9 tables. Small corrections adde

Dynamical Boson Stars

The idea of stable, localized bundles of energy has strong appeal as a model for particles. In the 1950s John Wheeler envisioned such bundles as smooth configurations of electromagnetic energy that he called {\em geons}, but none were found. Instead, particle-like solutions were found in the late 1960s with the addition of a scalar field, and these were given the name {\em boson stars}. Since then, boson stars find use in a wide variety of models as sources of dark matter, as black hole mimickers, in simple models of binary systems, and as a tool in finding black holes in higher dimensions with only a single killing vector. We discuss important varieties of boson stars, their dynamic properties, and some of their uses, concentrating on recent efforts.Comment: 79 pages, 25 figures, invited review for Living Reviews in Relativity; major revision in 201

Diffractive Dijet Production at sqrt(s)=630 and 1800 GeV at the Fermilab Tevatron

We report a measurement of the diffractive structure function $F_{jj}^D$ of the antiproton obtained from a study of dijet events produced in association with a leading antiproton in $\bar pp$ collisions at $\sqrt s=630$ GeV at the Fermilab Tevatron. The ratio of $F_{jj}^D$ at $\sqrt s=630$ GeV to $F_{jj}^D$ obtained from a similar measurement at $\sqrt s=1800$ GeV is compared with expectations from QCD factorization and with theoretical predictions. We also report a measurement of the $\xi$ ($x$-Pomeron) and $\beta$ ($x$ of parton in Pomeron) dependence of $F_{jj}^D$ at $\sqrt s=1800$ GeV. In the region $0.035<\xi<0.095$, $|t|<1$ GeV$^2$ and $\beta<0.5$, $F_{jj}^D(\beta,\xi)$ is found to be of the form $\beta^{-1.0\pm 0.1} \xi^{-0.9\pm 0.1}$, which obeys $\beta$-$\xi$ factorization.Comment: LaTeX, 9 pages, Submitted to Phys. Rev. Letter

A Study of B0 -> J/psi K(*)0 pi+ pi- Decays with the Collider Detector at Fermilab

We report a study of the decays B0 -> J/psi K(*)0 pi+ pi-, which involve the creation of a u u-bar or d d-bar quark pair in addition to a b-bar -> c-bar(c s-bar) decay. The data sample consists of 110 1/pb of p p-bar collisions at sqrt{s} = 1.8 TeV collected by the CDF detector at the Fermilab Tevatron collider during 1992-1995. We measure the branching ratios to be BR(B0 -> J/psi K*0 pi+ pi-) = (8.0 +- 2.2 +- 1.5) * 10^{-4} and BR(B0 -> J/psi K0 pi+ pi-) = (1.1 +- 0.4 +- 0.2) * 10^{-3}. Contributions to these decays are seen from psi(2S) K(*)0, J/psi K0 rho0, J/psi K*+ pi-, and J/psi K1(1270)

Measurement of the branching fraction and CP content for the decay B(0) -> D(*+)D(*-)

This is the pre-print version of the Article. The official published version can be accessed from the links below. Copyright @ 2002 APS.We report a measurement of the branching fraction of the decay B0→D*+D*- and of the CP-odd component of its final state using the BABAR detector. With data corresponding to an integrated luminosity of 20.4  fb-1 collected at the Υ(4S) resonance during 1999–2000, we have reconstructed 38 candidate signal events in the mode B0→D*+D*- with an estimated background of 6.2±0.5 events. From these events, we determine the branching fraction to be B(B0→D*+D*-)=[8.3±1.6(stat)±1.2(syst)]×10-4. The measured CP-odd fraction of the final state is 0.22±0.18(stat)±0.03(syst).This work is supported by DOE and NSF (USA), NSERC (Canada), IHEP (China), CEA and CNRS-IN2P3 (France), BMBF (Germany), INFN (Italy), NFR (Norway), MIST (Russia), and PPARC (United Kingdom). Individuals have received support from the A.P. Sloan Foundation, Research Corporation, and Alexander von Humboldt Foundation

Measurement of D-s(+) and D-s(*+) production in B meson decays and from continuum e(+)e(-) annihilation at √s=10.6 GeV

This is the pre-print version of the Article. The official published version can be accessed from the links below. Copyright @ 2002 APSNew measurements of Ds+ and Ds*+ meson production rates from B decays and from qq̅ continuum events near the Υ(4S) resonance are presented. Using 20.8 fb-1 of data on the Υ(4S) resonance and 2.6 fb-1 off-resonance, we find the inclusive branching fractions B(B⃗Ds+X)=(10.93±0.19±0.58±2.73)% and B(B⃗Ds*+X)=(7.9±0.8±0.7±2.0)%, where the first error is statistical, the second is systematic, and the third is due to the Ds+→φπ+ branching fraction uncertainty. The production cross sections σ(e+e-→Ds+X)×B(Ds+→φπ+)=7.55±0.20±0.34pb and σ(e+e-→Ds*±X)×B(Ds+→φπ+)=5.8±0.7±0.5pb are measured at center-of-mass energies about 40 MeV below the Υ(4S) mass. The branching fractions ΣB(B⃗Ds(*)+D(*))=(5.07±0.14±0.30±1.27)% and ΣB(B⃗Ds*+D(*))=(4.1±0.2±0.4±1.0)% are determined from the Ds(*)+ momentum spectra. The mass difference m(Ds+)-m(D+)=98.4±0.1±0.3MeV/c2 is also measured.This work was supported by DOE and NSF (USA), NSERC (Canada), IHEP (China), CEA and CNRS-IN2P3 (France), BMBF (Germany), INFN (Italy), NFR (Norway), MIST (Russia), and PPARC (United Kingdom). Individuals have received support from the Swiss NSF, A. P. Sloan Foundation, Research Corporation, and Alexander von Humboldt Foundation

A Study of Time-Dependent CP-Violating Asymmetries and Flavor Oscillations in Neutral B Decays at the Upsilon(4S)

We present a measurement of time-dependent CP-violating asymmetries in neutral B meson decays collected with the BABAR detector at the PEP-II asymmetric-energy B Factory at the Stanford Linear Accelerator Center. The data sample consists of 29.7 ${\rm fb}^{-1}$ recorded at the $\Upsilon(4S)$ resonance and 3.9 ${\rm fb}^{-1}$ off-resonance. One of the neutral B mesons, which are produced in pairs at the $\Upsilon(4S)$, is fully reconstructed in the CP decay modes $J/\psi K^0_S$, $\psi(2S) K^0_S$, $\chi_{c1} K^0_S$, $J/\psi K^{*0}$ ($K^{*0}\to K^0_S\pi^0$) and $J/\psi K^0_L$, or in flavor-eigenstate modes involving $D^{(*)}\pi/\rho/a_1$ and $J/\psi K^{*0}$ ($K^{*0}\to K^+\pi^-$). The flavor of the other neutral B meson is tagged at the time of its decay, mainly with the charge of identified leptons and kaons. The proper time elapsed between the decays is determined by measuring the distance between the decay vertices. A maximum-likelihood fit to this flavor eigenstate sample finds $\Delta m_d = 0.516\pm 0.016 {\rm (stat)} \pm 0.010 {\rm (syst)} {\rm ps}^{-1}$. The value of the asymmetry amplitude $\sin2\beta$ is determined from a simultaneous maximum-likelihood fit to the time-difference distribution of the flavor-eigenstate sample and about 642 tagged $B^0$ decays in the CP-eigenstate modes. We find $\sin2\beta=0.59\pm 0.14 {\rm (stat)} \pm 0.05 {\rm (syst)}$, demonstrating that CP violation exists in the neutral B meson system. (abridged)Comment: 58 pages, 35 figures, submitted to Physical Review

Measurement of the Branching Fraction for B- --> D0 K*-

We present a measurement of the branching fraction for the decay B- --> D0 K*- using a sample of approximately 86 million BBbar pairs collected by the BaBar detector from e+e- collisions near the Y(4S) resonance. The D0 is detected through its decays to K- pi+, K- pi+ pi0 and K- pi+ pi- pi+, and the K*- through its decay to K0S pi-. We measure the branching fraction to be B.F.(B- --> D0 K*-)= (6.3 +/- 0.7(stat.) +/- 0.5(syst.)) x 10^{-4}.Comment: 7 pages, 1 postscript figure, submitted to Phys. Rev. D (Rapid Communications

Measurement of Branching Fraction and Dalitz Distribution for B0->D(*)+/- K0 pi-/+ Decays

We present measurements of the branching fractions for the three-body decays B0 -> D(*)-/+ K0 pi^+/-$and their resonant submodes$B0 -> D(*)-/+ K*+/- using a sample of approximately 88 million BBbar pairs collected by the BABAR detector at the PEP-II asymmetric energy storage ring. We measure: B(B0->D-/+ K0 pi+/-)=(4.9 +/- 0.7(stat) +/- 0.5 (syst)) 10^{-4} B(B0->D*-/+ K0 pi+/-)=(3.0 +/- 0.7(stat) +/- 0.3 (syst)) 10^{-4} B(B0->D-/+ K*+/-)=(4.6 +/- 0.6(stat) +/- 0.5 (syst)) 10^{-4} B(B0->D*-/+ K*+/-)=(3.2 +/- 0.6(stat) +/- 0.3 (syst)) 10^{-4} From these measurements we determine the fractions of resonant events to be : f(B0-> D-/+ K*+/-) = 0.63 +/- 0.08(stat) +/- 0.04(syst) f(B0-> D*-/+ K*+/-) = 0.72 +/- 0.14(stat) +/- 0.05(syst)Comment: 7 pages, 3 figures submitted to Phys. Rev. Let

Evidence for the Rare Decay B -> K*ll and Measurement of the B -> Kll Branching Fraction

We present evidence for the flavor-changing neutral current decay $B\to K^*\ell^+\ell^-$ and a measurement of the branching fraction for the related process $B\to K\ell^+\ell^-$, where $\ell^+\ell^-$ is either an $e^+e^-$ or $\mu^+\mu^-$ pair. These decays are highly suppressed in the Standard Model, and they are sensitive to contributions from new particles in the intermediate state. The data sample comprises $123\times 10^6$ $\Upsilon(4S)\to B\bar{B}$ decays collected with the Babar detector at the PEP-II $e^+e^-$ storage ring. Averaging over $K^{(*)}$ isospin and lepton flavor, we obtain the branching fractions ${\mathcal B}(B\to K\ell^+\ell^-)=(0.65^{+0.14}_{-0.13}\pm 0.04)\times 10^{-6}$ and ${\mathcal B}(B\to K^*\ell^+\ell^-)=(0.88^{+0.33}_{-0.29}\pm 0.10)\times 10^{-6}$, where the uncertainties are statistical and systematic, respectively. The significance of the $B\to K\ell^+\ell^-$ signal is over $8\sigma$, while for $B\to K^*\ell^+\ell^-$ it is $3.3\sigma$.Comment: 7 pages, 2 postscript figues, submitted to Phys. Rev. Let