Inferring viral quasispecies spectra from 454 pyrosequencing reads

<p>Abstract</p> <p>Background</p> <p>RNA viruses infecting a host usually exist as a set of closely related sequences, referred to as quasispecies. The genomic diversity of viral quasispecies is a subject of great interest, particularly for chronic infections, since it can lead to resistance to existing therapies. High-throughput sequencing is a promising approach to characterizing viral diversity, but unfortunately standard assembly software was originally designed for single genome assembly and cannot be used to simultaneously assemble and estimate the abundance of multiple closely related quasispecies sequences.</p> <p>Results</p> <p>In this paper, we introduce a new <b>Vi</b>ral <b>Sp</b>ectrum <b>A</b>ssembler (ViSpA) method for quasispecies spectrum reconstruction and compare it with the state-of-the-art ShoRAH tool on both simulated and real 454 pyrosequencing shotgun reads from HCV and HIV quasispecies. Experimental results show that ViSpA outperforms ShoRAH on simulated error-free reads, correctly assembling 10 out of 10 quasispecies and 29 sequences out of 40 quasispecies. While ShoRAH has a significant advantage over ViSpA on reads simulated with sequencing errors due to its advanced error correction algorithm, ViSpA is better at assembling the simulated reads after they have been corrected by ShoRAH. ViSpA also outperforms ShoRAH on real 454 reads. Indeed, 7 most frequent sequences reconstructed by ViSpA from a real HCV dataset are viable (do not contain internal stop codons), and the most frequent sequence was within 1% of the actual open reading frame obtained by cloning and Sanger sequencing. In contrast, only one of the sequences reconstructed by ShoRAH is viable. On a real HIV dataset, ShoRAH correctly inferred only 2 quasispecies sequences with at most 4 mismatches whereas ViSpA correctly reconstructed 5 quasispecies with at most 2 mismatches, and 2 out of 5 sequences were inferred without any mismatches. ViSpA source code is available at <url>http://alla.cs.gsu.edu/~software/VISPA/vispa.html</url>.</p> <p>Conclusions</p> <p>ViSpA enables accurate viral quasispecies spectrum reconstruction from 454 pyrosequencing reads. We are currently exploring extensions applicable to the analysis of high-throughput sequencing data from bacterial metagenomic samples and ecological samples of eukaryote populations.</p

GDNF Secreting Human Neural Progenitor Cells Protect Dying Motor Neurons, but Not Their Projection to Muscle, in a Rat Model of Familial ALS

Amyotrophic lateral sclerosis (ALS) is a fatal, progressive neurodegenerative disease characterized by rapid loss of muscle control and eventual paralysis due to the death of large motor neurons in the brain and spinal cord. Growth factors such as glial cell line derived neurotrophic factor (GDNF) are known to protect motor neurons from damage in a range of models. However, penetrance through the blood brain barrier and delivery to the spinal cord remains a serious challenge. Although there may be a primary dysfunction in the motor neuron itself, there is also increasing evidence that excitotoxicity due to glial dysfunction plays a crucial role in disease progression. Clearly it would be of great interest if wild type glial cells could ameliorate motor neuron loss in these models, perhaps in combination with the release of growth factors such as GDNF.Human neural progenitor cells can be expanded in culture for long periods and survive transplantation into the adult rodent central nervous system, in some cases making large numbers of GFAP positive astrocytes. They can also be genetically modified to release GDNF (hNPC(GDNF)) and thus act as long-term 'mini pumps' in specific regions of the rodent and primate brain. In the current study we genetically modified human neural stem cells to release GDNF and transplanted them into the spinal cord of rats over-expressing mutant SOD1 (SOD1(G93A)). Following unilateral transplantation into the spinal cord of SOD1(G93A) rats there was robust cellular migration into degenerating areas, efficient delivery of GDNF and remarkable preservation of motor neurons at early and end stages of the disease within chimeric regions. The progenitors retained immature markers, and those not secreting GDNF had no effect on motor neuron survival. Interestingly, this robust motor neuron survival was not accompanied by continued innervation of muscle end plates and thus resulted in no improvement in ipsilateral limb use.The potential to maintain dying motor neurons by delivering GDNF using neural progenitor cells represents a novel and powerful treatment strategy for ALS. While this approach represents a unique way to prevent motor neuron loss, our data also suggest that additional strategies may also be required for maintenance of neuromuscular connections and full functional recovery. However, simply maintaining motor neurons in patients would be the first step of a therapeutic advance for this devastating and incurable disease, while future strategies focus on the maintenance of the neuromuscular junction

Measurement of the Lifetime and Λ Separation Energy of _{Λ}^{3}H

The most precise measurements to date of the _{Λ}^{3}H lifetime τ and Λ separation energy B_{Λ} are obtained using the data sample of Pb-Pb collisions at sqrt[s_{NN}]=5.02  TeV collected by ALICE at the LHC. The _{Λ}^{3}H is reconstructed via its charged two-body mesonic decay channel (_{Λ}^{3}H→^{3}He+π^{-} and the charge-conjugate process). The measured values τ=[253±11(stat)±6(syst)]  ps and B_{Λ}=[102±63(stat)±67(syst)]  keV are compatible with predictions from effective field theories and confirm that the _{Λ}^{3}H structure is consistent with a weakly bound system

f0(980) production in inelastic pp collisions at s = 5.02 TeV

The measurement of the production of f0(980) in inelastic pp collisions at sqrt(s) = 5.02 TeV is presented. This is the first reported measurement of inclusive f0(980) yield at LHC energies. The production is measured at midrapidity, |y| pi+pi- hadronic decay channel using the ALICE detector. The pT-differential yields are compared to those of pions, protons and ϕ mesons as well as to predictions from the HERWIG 7.2 QCD-inspired Monte Carlo event generator and calculations from a coalescence model that uses the AMPT model as an input. The ratio of the pT-integrated yield of f0(980) relative to pions is compared to measurements in e+e- and pp collisions at lower energies and predictions from statistical hadronisation models and HERWIG 7.2. A mild collision energy dependence of the f0(980) to pion production is observed in pp collisions from SPS to LHC energies. All considered models underpredict the pT-integrated 2f0(980)/(pi+ + pi-) ratio. The prediction from the canonical statistical hadronisation model assuming a zero total strangeness content of f0(980) is consistent with the data within 1.9σ and is the closest to the data. The results provide an essential reference for future measurements of the particle yield and nuclear modification in p–Pb and Pb–Pb collisions, which have been proposed to be instrumental to probe the elusive nature and quark composition of the f0(980) scalar meson

Observation of flow angle and flow magnitude fluctuations in Pb-Pb collisions at sNN=5.02TeV at the CERN Large Hadron Collider

This Letter reports on the first measurements of transverse momentum dependent flow angle n and flow magnitude vn fluctuations determined using new four-particle correlators. The measurements are performed for various centralities in Pb–Pb collisions at a center-of-mass energy per nucleon pair of √s NN = 5.02 TeV with ALICE at the CERN Large Hadron Collider. Both flow angle and flow magnitude fluctuations are observed in the presented centrality ranges and are strongest in the most central collisions and for a transverse momentum pT > 2 GeV/c. Comparison with theoretical models, including iEBE-VISHNU, MUSIC, and AMPT, show that the measurements exhibit unique sensitivities to the initial state of heavy-ion collisions

First measurement of Ωc 0 production in pp collisions at s=13 TeV

The inclusive production of the charm–strange baryon Omega_c^0 is measured for the first time via its hadronic decay into Omega-pi+ at midrapidity (|y|<0.5) in proton–proton (pp) collisions at the centre-of-mass energy sqrt(s) = 13 TeV with the ALICE detector at the LHC. The transverse momentum (pT) differential cross section multiplied by the branching ratio is presented in the interval 2 < pT < 12 GeV/c . The pT dependence of the Omega_C^0-baryon production relative to the prompt D^0-meson and to the prompt Csi_C^0-baryon production is compared to various models that take different hadronisation mechanisms into consideration. In the measured pT interval, the ratio of the pT-integrated cross sections of Omega_c^0 and prompt Lambda_c^+ baryons multiplied by the Omega- pi+ branching ratio is found to be larger by a factor of about 20 with a significance of about 4σ when compared to e+e- collisions

Measurement of the non-prompt D-meson fraction as a function of multiplicity in proton-proton collisions at s \sqrt{s} = 13 TeV

The fractions of non-prompt (i.e. originating from beauty-hadron decays) D0 and D+ mesons with respect to the inclusive yield are measured as a function of the charged-particle multiplicity in proton-proton collisions at a centre-of-mass energy of √s = 13 TeV with the ALICE detector at the LHC. The results are reported in intervals of transverse momentum (pT) and integrated in the range 1 < pT < 24 GeV/c. The fraction of non-prompt D0 and D+ mesons is found to increase slightly as a function of pT in all the measured multiplicity intervals, while no significant dependence on the charged- particle multiplicity is observed. In order to investigate the production and hadronisation mechanisms of charm and beauty quarks, the results are compared to PYTHIA 8 as well as EPOS 3 and EPOS 4 Monte Carlo simulations, and to calculations based on the colour glass condensate including three-pomeron fusion

Accessing the strong interaction between Λ baryons and charged kaons with the femtoscopy technique at the LHC

The interaction between Λ baryons and kaons/antikaons is a crucial ingredient for the strangeness S=0 and S=-2 sector of the meson–baryon interaction at low energies. In particular, the Lambda-Kbar might help in understanding the origin of states such as the Csi(1620), whose nature and properties are still under debate. Experimental data on Lambda-K and Lambda-Kbar systems are scarce, leading to large uncertainties and tension between the available theoretical predictions constrained by such data. In this Letter we present the measurements of Λ–KK− and Λ–KK+ correlations obtained in the high-multiplicity triggered data sample in pp collisions at sqrt(s) = 13 TeV recorded by ALICE at the LHC. The correlation function for both pairs is modeled using the Lednický–Lyuboshits analytical formula and the corresponding scattering parameters are extracted. The Λ–KK+ correlations show the presence of several structures at relative momenta k* above 200 MeV/c, compatible with the Ω baryon, the , and resonances decaying into Λ–K− pairs. The low k* region in the Λ–KK+ also exhibits the presence of the state, expected to strongly couple to the measured pair. The presented data allow to access the ΛK+ and ΛK− strong interaction with an unprecedented precision and deliver the first experimental observation of the decaying into ΛK−

Investigation of K+K- interactions via femtoscopy in Pb-Pb collisions at sNN =2.76 TeV at the CERN Large Hadron Collider

Femtoscopic correlations of nonidentical charged kaons (K+K-) are studied in Pb-Pb collisions at a center-of-mass energy per nucleon-nucleon collision sNN=2.76 TeV by ALICE at the CERN Large Hadron Collider. One-dimensional K+K- correlation functions are analyzed in three centrality classes and eight intervals of particle-pair transverse momentum. The Lednický and Luboshitz interaction model used in the K+K- analysis includes the final-state Coulomb interactions between kaons and the final-state interaction through a0(980) and f0(980) resonances. The mass of f0(980) and coupling were extracted from the fit to K+K- correlation functions using the femtoscopic technique. The measured mass and width of the f0(980) resonance are consistent with other published measurements. The height of the φ(1020) meson peak present in the K+K- correlation function rapidly decreases with increasing source radius, qualitatively in agreement with an inverse volume dependence. A phenomenological fit to this trend suggests that the φ(1020) meson yield is dominated by particles produced directly from the hadronization of the system. The small fraction subsequently produced by final-state interactions could not be precisely quantified with data presented in this paper and will be assessed in future work