Viral population estimation using pyrosequencing

The diversity of virus populations within single infected hosts presents a major difficulty for the natural immune response as well as for vaccine design and antiviral drug therapy. Recently developed pyrophosphate based sequencing technologies (pyrosequencing) can be used for quantifying this diversity by ultra-deep sequencing of virus samples. We present computational methods for the analysis of such sequence data and apply these techniques to pyrosequencing data obtained from HIV populations within patients harboring drug resistant virus strains. Our main result is the estimation of the population structure of the sample from the pyrosequencing reads. This inference is based on a statistical approach to error correction, followed by a combinatorial algorithm for constructing a minimal set of haplotypes that explain the data. Using this set of explaining haplotypes, we apply a statistical model to infer the frequencies of the haplotypes in the population via an EM algorithm. We demonstrate that pyrosequencing reads allow for effective population reconstruction by extensive simulations and by comparison to 165 sequences obtained directly from clonal sequencing of four independent, diverse HIV populations. Thus, pyrosequencing can be used for cost-effective estimation of the structure of virus populations, promising new insights into viral evolutionary dynamics and disease control strategies.Comment: 23 pages, 13 figure

Revisiting Scalar and Pseudoscalar Couplings with Nucleons

Certain dark matter interactions with nuclei are mediated possibly by a scalar or pseudoscalar Higgs boson. The estimation of the corresponding cross sections requires a correct evaluation of the couplings between the scalar or pseudoscalar Higgs boson and the nucleons. Progress has been made in two aspects relevant to this study in the past few years. First, recent lattice calculations show that the strange-quark sigma term $\sigma_s$ and the strange-quark content in the nucleon are much smaller than what are expected previously. Second, lattice and model analyses imply sizable SU(3) breaking effects in the determination on the axial-vector coupling constant $g_A^8$ that in turn affect the extraction of the isosinglet coupling $g_A^0$ and the strange quark spin component $\Delta s$ from polarized deep inelastic scattering experiments. Based on these new developments, we re-evaluate the relevant nucleon matrix elements and compute the scalar and pseudoscalar couplings of the proton and neutron. We also find that the strange quark contribution in both types of couplings is smaller than previously thought.Comment: 17 pages, Sec. II is revised and the pion-nucleon sigma term extracted from the scattering data is discussed. Version to appear in JHE

First Measurement of Transferred Polarization in the Exclusive e p --> e' K+ Lambda Reaction

The first measurements of the transferred polarization for the exclusive ep --> e'K+ Lambda reaction have been performed in Hall B at the Thomas Jefferson National Accelerator Facility using the CLAS spectrometer. A 2.567 GeV electron beam was used to measure the hyperon polarization over a range of Q2 from 0.3 to 1.5 (GeV/c)2, W from 1.6 to 2.15 GeV, and over the full center-of-mass angular range of the K+ meson. Comparison with predictions of hadrodynamic models indicates strong sensitivity to the underlying resonance contributions. A non-relativistic quark model interpretation of our data suggests that the s-sbar quark pair is produced with spins predominantly anti-aligned. Implications for the validity of the widely used 3P0 quark-pair creation operator are discussed.Comment: 6 pages, 4 figure

Onset of asymptotic scaling in deuteron photodisintegration

We investigate the transition from the nucleon-meson to quark-gluon description of the strong interaction using the photon energy dependence of the $d(\gamma,p)n$ differential cross section for photon energies above 0.5 GeV and center-of-mass proton angles between $30^{\circ}$ and $150^{\circ}$. A possible signature for this transition is the onset of cross section $s^{-11}$ scaling with the total energy squared, $s$, at some proton transverse momentum, $P_T$. The results show that the scaling has been reached for proton transverse momentum above about 1.1 GeV/c. This may indicate that the quark-gluon regime is reached above this momentum.Comment: Accepted by PRL; 5 pages, 2 figure

Two-Nucleon Momentum Distributions Measured in 3He(e,e'pp)n

We have measured the 3He(e,e'pp)n reaction at 2.2 GeV over a wide kinematic range. The kinetic energy distribution for `fast' nucleons (p > 250 MeV/c) peaks where two nucleons each have 20% or less, and the third nucleon has most of the transferred energy. These fast pp and pn pairs are back-to-back with little momentum along the three-momentum transfer, indicating that they are spectators. Experimental and theoretical evidence indicates that we have measured distorted two-nucleon momentum distributions by striking the third nucleon and detecting the spectator correlated pair.Comment: 6 pages, 5 figures, submitted to PR

Measurement of the Polarized Structure Function σLT′\sigma_{LT^\prime} for p(e⃗,e′p)πop(\vec{e},e'p)\pi^o in the Δ(1232)\Delta(1232) Resonance Region

The polarized longitudinal-transverse structure function $\sigma_{LT^\prime}$ has been measured in the $\Delta(1232)$ resonance region at $Q^2=0.40$ and 0.65 GeV$^2$. Data for the $p(\vec e,e'p)\pi^o$ reaction were taken at Jefferson Lab with the CEBAF Large Acceptance Spectrometer (CLAS) using longitudinally polarized electrons at an energy of 1.515 GeV. For the first time a complete angular distribution was measured, permitting the separation of different non-resonant amplitudes using a partial wave analysis. Comparison with previous beam asymmetry measurements at MAMI indicate a deviation from the predicted $Q^2$ dependence of $\sigma_{LT^{\prime}}$ using recent phenomenological models.Comment: 5 pages, LaTex, 4 eps figures: to be published in PRC/Rapid Communications. Version 2 has revised Q^2 analysi

Exclusive ρ0\rho^0 electroproduction on the proton at CLAS

The $e p\to e^\prime p \rho^0$ reaction has been measured, using the 5.754 GeV electron beam of Jefferson Lab and the CLAS detector. This represents the largest ever set of data for this reaction in the valence region. Integrated and differential cross sections are presented. The $W$, $Q^2$ and $t$ dependences of the cross section are compared to theoretical calculations based on $t$-channel meson-exchange Regge theory on the one hand and on quark handbag diagrams related to Generalized Parton Distributions (GPDs) on the other hand. The Regge approach can describe at the $\approx$ 30% level most of the features of the present data while the two GPD calculations that are presented in this article which succesfully reproduce the high energy data strongly underestimate the present data. The question is then raised whether this discrepancy originates from an incomplete or inexact way of modelling the GPDs or the associated hard scattering amplitude or whether the GPD formalism is simply inapplicable in this region due to higher-twists contributions, incalculable at present.Comment: 29 pages, 29 figure

Observation of an Exotic S=+1S=+1 Baryon in Exclusive Photoproduction from the Deuteron

In an exclusive measurement of the reaction $\gamma d \to K^+ K^- p n$, a narrow peak that can be attributed to an exotic baryon with strangeness $S=+1$ is seen in the $K^+n$ invariant mass spectrum. The peak is at $1.542\pm 0.005$ GeV/c$^2$ with a measured width of 0.021 GeV/c$^2$ FWHM, which is largely determined by experimental mass resolution. The statistical significance of the peak is $5.2 \pm 0.6 \sigma$. The mass and width of the observed peak are consistent with recent reports of a narrow $S=+1$ baryon by other experimental groups.Comment: 5 pages, 5 figure

Single pi+ Electroproduction on the Proton in the First and Second Resonance Regions at 0.25GeV^2 < Q^2 < 0.65GeV^2 Using CLAS

The ep -> e'pi^+n reaction was studied in the first and second nucleon resonance regions in the 0.25 GeV^2 < Q^2 < 0.65 GeV^2 range using the CLAS detector at Thomas Jefferson National Accelerator Facility. For the first time the absolute cross sections were measured covering nearly the full angular range in the hadronic center-of-mass frame. The structure functions sigma_TL, sigma_TT and the linear combination sigma_T+epsilon*sigma_L were extracted by fitting the phi-dependence of the measured cross sections, and were compared to the MAID and Sato-Lee models.Comment: Accepted for publication in PR

Photodisintegration of 4^4He into p+t

The two-body photodisintegration of $^4$He into a proton and a triton has been studied using the CEBAF Large-Acceptance Spectrometer (CLAS) at Jefferson Laboratory. Real photons produced with the Hall-B bremsstrahlung-tagging system in the energy range from 0.35 to 1.55 GeV were incident on a liquid $^4$He target. This is the first measurement of the photodisintegration of $^4$He above 0.4 GeV. The differential cross sections for the $\gamma$$^4$He$\to pt$ reaction have been measured as a function of photon-beam energy and proton-scattering angle, and are compared with the latest model calculations by J.-M. Laget. At 0.6-1.2 GeV, our data are in good agreement only with the calculations that include three-body mechanisms, thus confirming their importance. These results reinforce the conclusion of our previous study of the three-body breakup of $^3$He that demonstrated the great importance of three-body mechanisms in the energy region 0.5-0.8 GeV .Comment: 13 pages submitted in one tgz file containing 2 tex file and 22 postscrip figure