The Reproducibility of Lists of Differentially Expressed Genes in Microarray Studies

Reproducibility is a fundamental requirement in scientific experiments and clinical contexts. Recent publications raise concerns about the reliability of microarray technology because of the apparent lack of agreement between lists of differentially expressed genes (DEGs). In this study we demonstrate that (1) such discordance may stem from ranking and selecting DEGs solely by statistical significance (P) derived from widely used simple t-tests; (2) when fold change (FC) is used as the ranking criterion, the lists become much more reproducible, especially when fewer genes are selected; and (3) the instability of short DEG lists based on P cutoffs is an expected mathematical consequence of the high variability of the t-values. We recommend the use of FC ranking plus a non-stringent P cutoff as a baseline practice in order to generate more reproducible DEG lists. The FC criterion enhances reproducibility while the P criterion balances sensitivity and specificity

The balance of reproducibility, sensitivity, and specificity of lists of differentially expressed genes in microarray studies

<p>Abstract</p> <p>Background</p> <p>Reproducibility is a fundamental requirement in scientific experiments. Some recent publications have claimed that microarrays are unreliable because lists of differentially expressed genes (DEGs) are not reproducible in similar experiments. Meanwhile, new statistical methods for identifying DEGs continue to appear in the scientific literature. The resultant variety of existing and emerging methods exacerbates confusion and continuing debate in the microarray community on the appropriate choice of methods for identifying reliable DEG lists.</p> <p>Results</p> <p>Using the data sets generated by the MicroArray Quality Control (MAQC) project, we investigated the impact on the reproducibility of DEG lists of a few widely used gene selection procedures. We present comprehensive results from inter-site comparisons using the same microarray platform, cross-platform comparisons using multiple microarray platforms, and comparisons between microarray results and those from TaqMan – the widely regarded "standard" gene expression platform. Our results demonstrate that (1) previously reported discordance between DEG lists could simply result from ranking and selecting DEGs solely by statistical significance (<it>P</it>) derived from widely used simple <it>t</it>-tests; (2) when fold change (FC) is used as the ranking criterion with a non-stringent <it>P</it>-value cutoff filtering, the DEG lists become much more reproducible, especially when fewer genes are selected as differentially expressed, as is the case in most microarray studies; and (3) the instability of short DEG lists solely based on <it>P</it>-value ranking is an expected mathematical consequence of the high variability of the <it>t</it>-values; the more stringent the <it>P</it>-value threshold, the less reproducible the DEG list is. These observations are also consistent with results from extensive simulation calculations.</p> <p>Conclusion</p> <p>We recommend the use of FC-ranking plus a non-stringent <it>P </it>cutoff as a straightforward and baseline practice in order to generate more reproducible DEG lists. Specifically, the <it>P</it>-value cutoff should not be stringent (too small) and FC should be as large as possible. Our results provide practical guidance to choose the appropriate FC and <it>P</it>-value cutoffs when selecting a given number of DEGs. The FC criterion enhances reproducibility, whereas the <it>P </it>criterion balances sensitivity and specificity.</p

LHCb upgrade software and computing : technical design report

This document reports the Research and Development activities that are carried out in the software and computing domains in view of the upgrade of the LHCb experiment. The implementation of a full software trigger implies major changes in the core software framework, in the event data model, and in the reconstruction algorithms. The increase of the data volumes for both real and simulated datasets requires a corresponding scaling of the distributed computing infrastructure. An implementation plan in both domains is presented, together with a risk assessment analysis

Physics case for an LHCb Upgrade II - Opportunities in flavour physics, and beyond, in the HL-LHC era

The LHCb Upgrade II will fully exploit the flavour-physics opportunities of the HL-LHC, and study additional physics topics that take advantage of the forward acceptance of the LHCb spectrometer. The LHCb Upgrade I will begin operation in 2020. Consolidation will occur, and modest enhancements of the Upgrade I detector will be installed, in Long Shutdown 3 of the LHC (2025) and these are discussed here. The main Upgrade II detector will be installed in long shutdown 4 of the LHC (2030) and will build on the strengths of the current LHCb experiment and the Upgrade I. It will operate at a luminosity up to 2×1034 cm−2s−1, ten times that of the Upgrade I detector. New detector components will improve the intrinsic performance of the experiment in certain key areas. An Expression Of Interest proposing Upgrade II was submitted in February 2017. The physics case for the Upgrade II is presented here in more depth. CP-violating phases will be measured with precisions unattainable at any other envisaged facility. The experiment will probe b → sl+l−and b → dl+l− transitions in both muon and electron decays in modes not accessible at Upgrade I. Minimal flavour violation will be tested with a precision measurement of the ratio of B(B0 → μ+μ−)/B(Bs → μ+μ−). Probing charm CP violation at the 10−5 level may result in its long sought discovery. Major advances in hadron spectroscopy will be possible, which will be powerful probes of low energy QCD. Upgrade II potentially will have the highest sensitivity of all the LHC experiments on the Higgs to charm-quark couplings. Generically, the new physics mass scale probed, for fixed couplings, will almost double compared with the pre-HL-LHC era; this extended reach for flavour physics is similar to that which would be achieved by the HE-LHC proposal for the energy frontier

Study of the psi(2)(3823) and chi(c1)(3872) states in B+->(J/psi pi(+)pi(-))K(+)decays

The decays $B^+\rightarrow J/\psi \pi^+ \pi^- K^+$ are studied using a data set corresponding to an integrated luminosity of 9fb$^{-1}$ collected with the LHCb detector in proton-proton collisions between 2011 and 2018. Precise measurements of the ratios of branching fractions with the intermediate $\psi_2(3823)$, $\chi_{c1}(3872)$ and $\psi(2S)$ states are reported. The decay of $B^+\rightarrow \psi_2(3872)K^+$ with $\psi_2(3823)\rightarrow J\psi\pi^+\pi^-$ is observed for the first time with a significance of 5.1 standard deviations. The mass differences between the $\psi_2(3823)$, $\chi_{c1}(3872)$ and $\psi(2S)$ states are measured to be $$\begin{array}{rcl} m_{\chi_{c1(3872)}} - m_{\psi_2(3823)} &= & 47.50 \pm 0.53 \pm 0.13\,\mathrm{MeV/}c^2\,, \\ m_{\psi_2(3823)} - m_{\psi(2S)} &= & 137.98 \pm 0.53 \pm 0.14\,\mathrm{MeV/}c^2\,, \\ m_{\chi_{c1}(3872)} - m_{\psi(2S)} &= & 185.49 \pm 0.06 \pm 0.03\,\mathrm{MeV/}c^2\,, \end{array}$$ resulting in the most precise determination of the $\chi_{c1}(3782)$ mass. The width of the $\psi_2(3823)$ state is found to be below 5.2MeV at 90\% confidence level. The Breit-Wigner width of the $\chi_{c1}(3872)$ state is measured to be $$\Gamma^{\mathrm{BW}}_{\chi_{c1}(3872)} = 0.96^{+0.19}_{-0.18}\pm0.21 \mathrm{MeV},$$ which is inconsistent with zero by 5.5 standard deviations.Comment: 26 pages, 3 figures. All figures and tables, along with any supplementary material and additional information, are available at https://cern.ch/lhcbproject/Publications/p/LHCb-PAPER-2020-009.html (LHCb public pages

Proceedings of the 13th annual conference of INEBRIA

CITATION: Watson, R., et al. 2016. Proceedings of the 13th annual conference of INEBRIA. Addiction Science & Clinical Practice, 11:13, doi:10.1186/s13722-016-0062-9.The original publication is available at https://ascpjournal.biomedcentral.comENGLISH SUMMARY : Meeting abstracts.https://ascpjournal.biomedcentral.com/articles/10.1186/s13722-016-0062-9Publisher's versio

Study of the lineshape of the chi(c1) (3872) state

A study of the lineshape of the chi(c1) (3872) state is made using a data sample corresponding to an integrated luminosity of 3 fb(-1) collected in pp collisions at center-of-mass energies of 7 and 8 TeV with the LHCb detector. Candidate chi(c1)(3872) and psi(2S) mesons from b-hadron decays are selected in the J/psi pi(+)pi(-) decay mode. Describing the lineshape with a Breit-Wigner function, the mass splitting between the chi(c1 )(3872) and psi(2S) states, Delta m, and the width of the chi(c1 )(3872) state, Gamma(Bw), are determined to be (Delta m=185.598 +/- 0.067 +/- 0.068 Mev,)(Gamma BW=1.39 +/- 0.24 +/- 0.10 Mev,) where the first uncertainty is statistical and the second systematic. Using a Flatte-inspired model, the mode and full width at half maximum of the lineshape are determined to be (mode=3871.69+0.00+0.05 MeV.)(FWHM=0.22-0.04+0.13+0.07+0.11-0.06-0.13 MeV, ) An investigation of the analytic structure of the Flatte amplitude reveals a pole structure, which is compatible with a quasibound D-0(D) over bar*(0) state but a quasivirtual state is still allowed at the level of 2 standard deviations

Study of the ψ2(3823)ψ_2(3823) and χc1(3872)χ_{c1}(3872) states in B+→(Jψπ+π−)K+B^+ \rightarrow \left( Jψπ^+π^-\right)K^+ decays

The decays $B^+\rightarrow J/\psi \pi^+ \pi^- K^+$ are studied using a data set corresponding to an integrated luminosity of 9fb$^{-1}$ collected with the LHCb detector in proton-proton collisions between 2011 and 2018. Precise measurements of the ratios of branching fractions with the intermediate $\psi_2(3823)$, $\chi_{c1}(3872)$ and $\psi(2S)$ states are reported. The decay of $B^+\rightarrow \psi_2(3823)K^+$ with $\psi_2(3823)\rightarrow J\psi\pi^+\pi^-$ is observed for the first time with a significance of 5.1 standard deviations. The mass differences between the $\psi_2(3823)$, $\chi_{c1}(3872)$ and $\psi(2S)$ states are measured to be $$\begin{array}{rcl} m_{\chi_{c1(3872)}} - m_{\psi_2(3823)} &= & 47.50 \pm 0.53 \pm 0.13\,\mathrm{MeV/}c^2\,, \\ m_{\psi_2(3823)} - m_{\psi(2S)} &= & 137.98 \pm 0.53 \pm 0.14\,\mathrm{MeV/}c^2\,, \\ m_{\chi_{c1}(3872)} - m_{\psi(2S)} &= & 185.49 \pm 0.06 \pm 0.03\,\mathrm{MeV/}c^2\,, \end{array}$$ resulting in the most precise determination of the $\chi_{c1}(3782)$ mass. The width of the $\psi_2(3823)$ state is found to be below 5.2MeV at 90\% confidence level. The Breit-Wigner width of the $\chi_{c1}(3872)$ state is measured to be $$\Gamma^{\mathrm{BW}}_{\chi_{c1}(3872)} = 0.96^{+0.19}_{-0.18}\pm0.21 \mathrm{MeV},$$ which is inconsistent with zero by 5.5 standard deviations

Measurement of the shape of the B s 0 → D s ∗ − μ + ν μ differential decay rate

Abstract: The shape of the Bs0→Ds∗−μ+νμ differential decay rate is obtained as a function of the hadron recoil parameter using proton-proton collision data at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 1.7 fb−1 collected by the LHCb detector. The Bs0→Ds∗−μ+νμ decay is reconstructed through the decays Ds∗−→Ds−γ and Ds−→K−K+π−. The differential decay rate is fitted with the Caprini-Lellouch-Neubert (CLN) and Boyd-Grinstein-Lebed (BGL) parametrisations of the form factors, and the relevant quantities for both are extracted