Unsupervised protein embeddings outperform hand-crafted sequence and structure features at predicting molecular function

Motivation: Protein function prediction is a difficult bioinformatics problem. Many recent methods use deep neural networks to learn complex sequence representations and predict function from these. Deep supervised models require a lot of labeled training data which are not available for this task. However, a very large amount of protein sequences without functional labels is available.Results: We applied an existing deep sequence model that had been pretrained in an unsupervised setting on the supervised task of protein molecular function prediction. We found that this complex feature representation is effective for this task, outperforming hand-crafted features such as one-hot encoding of amino acids, k-mer counts, secondary structure and backbone angles. Also, it partly negates the need for complex prediction models, as a two-layer perceptron was enough to achieve competitive performance in the third Critical Assessment of Functional Annotation benchmark. We also show that combining this sequence representation with protein 3D structure information does not lead to performance improvement, hinting that 3D structure is also potentially learned during the unsupervised pretraining

The genome sequencing of an albino Western lowland gorilla reveals inbreeding in the wild

Background The only known albino gorilla, named Snowflake, was a male wild born individual from Equatorial Guinea who lived at the Barcelona Zoo for almost 40 years. He was diagnosed with non-syndromic oculocutaneous albinism, i.e. white hair, light eyes, pink skin, photophobia and reduced visual acuity. Despite previous efforts to explain the genetic cause, this is still unknown. Here, we study the genetic cause of his albinism and making use of whole genome sequencing data we find a higher inbreeding coefficient compared to other gorillas. Results We successfully identified the causal genetic variant for Snowflake¿s albinism, a non-synonymous single nucleotide variant located in a transmembrane region of SLC45A2. This transporter is known to be involved in oculocutaneous albinism type 4 (OCA4) in humans. We provide experimental evidence that shows that this amino acid replacement alters the membrane spanning capability of this transmembrane region. Finally, we provide a comprehensive study of genome-wide patterns of autozygogosity revealing that Snowflake¿s parents were related, being this the first report of inbreeding in a wild born Western lowland gorilla. Conclusions In this study we demonstrate how the use of whole genome sequencing can be extended to link genotype and phenotype in non-model organisms and it can be a powerful tool in conservation genetics (e.g., inbreeding and genetic diversity) with the expected decrease in sequencing cost. Keywords: Gorilla; Albinism; Inbreeding; Genome; Conservatio

Genome-Wide Association Study of the Modified Stumvoll Insulin Sensitivity Index Identifies BCL2 and FAM19A2 as Novel Insulin Sensitivity Loci

Genome-wide association studies (GWAS) have found few common variants that influence fasting measures of insulin sensitivity. We hypothesized that a GWAS of an integrated assessment of fasting and dynamic measures of insulin sensitivity would detect novel common variants. We performed a GWAS of the modified Stumvoll Insulin Sensitivity Index (ISI) within the Meta-Analyses of Glucose and Insulin-Related Traits Consortium. Discovery for genetic association was performed in 16,753 individuals, and replication was attempted for the 23 most significant novel loci in 13,354 independent individuals. Association with ISI was tested in models adjusted for age, sex, and BMI and in a model analyzing the combined influence of the genotype effect adjusted for BMI and the interaction effect between the genotype and BMI on ISI (model 3). In model 3, three variants reached genome-wide significance: Rs13422522 (NYAP2; P = 8.87 × 10-11), rs12454712 (BCL2; P = 2.7 × 10-8), and rs10506418 (FAM19A2; P = 1.9 × 10-8). The association at NYAP2 was eliminated by conditioning on the known IRS1 insulin sensitivity locus; the BCL2 and FAM19A2 associations were independent of known cardiometabolic loci. In conclusion, we identified two novel loci and replicated known variants associated with insulin sensitivity. Further studies are needed to clarify the causal variant and function at the BCL2 and FAM19A2 loci

Nucleic acid amplification tests in the diagnosis of tuberculous pleuritis: a systematic review and meta-analysis

BACKGROUND: Conventional tests for tuberculous pleuritis have several limitations. A variety of new, rapid tests such as nucleic acid amplification tests – including polymerase chain reaction – have been evaluated in recent times. We conducted a systematic review to determine the accuracy of nucleic acid amplification (NAA) tests in the diagnosis of tuberculous pleuritis. METHODS: A systematic review and meta-analysis of 38 English and Spanish articles (with 40 studies), identified via searches of six electronic databases, hand searching of selected journals, and contact with authors, experts, and test manufacturers. Sensitivity, specificity, and other measures of accuracy were pooled using random effects models. Summary receiver operating characteristic curves were used to summarize overall test performance. Heterogeneity in study results was formally explored using subgroup analyses. RESULTS: Of the 40 studies included, 26 used in-house ("home-brew") tests, and 14 used commercial tests. Commercial tests had a low overall sensitivity (0.62; 95% confidence interval [CI] 0.43, 0.77), and high specificity (0.98; 95% CI 0.96, 0.98). The positive and negative likelihood ratios for commercial tests were 25.4 (95% CI 16.2, 40.0) and 0.40 (95% CI 0.24, 0.67), respectively. All commercial tests had consistently high specificity estimates; the sensitivity estimates, however, were heterogeneous across studies. With the in-house tests, both sensitivity and specificity estimates were significantly heterogeneous. Clinically meaningful summary estimates could not be determined for in-house tests. CONCLUSIONS: Our results suggest that commercial NAA tests may have a potential role in confirming (ruling in) tuberculous pleuritis. However, these tests have low and variable sensitivity and, therefore, may not be useful in excluding (ruling out) the disease. NAA test results, therefore, cannot replace conventional tests; they need to be interpreted in parallel with clinical findings and results of conventional tests. The accuracy of in-house nucleic acid amplification tests is poorly defined because of heterogeneity in study results. The clinical applicability of in-house NAA tests remains unclear

Studies of η\eta and η′\eta' production in pppp and ppPb collisions

The production of $\eta$ and $\eta'$ mesons is studied in proton-proton and proton-lead collisions collected with the LHCb detector. Proton-proton collisions are studied at center-of-mass energies of $5.02$ and $13~{\rm TeV}$, and proton-lead collisions are studied at a center-of-mass energy per nucleon of $8.16~{\rm TeV}$. The studies are performed in center-of-mass rapidity regions $2.5<y_{\rm c.m.}<3.5$ (forward rapidity) and $-4.0<y_{\rm c.m.}<-3.0$ (backward rapidity) defined relative to the proton beam direction. The $\eta$ and $\eta'$ production cross sections are measured differentially as a function of transverse momentum for $1.5<p_{\rm T}<10~{\rm GeV}$ and $3<p_{\rm T}<10~{\rm GeV}$, respectively. The differential cross sections are used to calculate nuclear modification factors. The nuclear modification factors for $\eta$ and $\eta'$ mesons agree at both forward and backward rapidity, showing no significant evidence of mass dependence. The differential cross sections of $\eta$ mesons are also used to calculate $\eta/\pi^0$ cross section ratios, which show evidence of a deviation from the world average. These studies offer new constraints on mass-dependent nuclear effects in heavy-ion collisions, as well as $\eta$ and $\eta'$ meson fragmentation.Comment: All figures and tables, along with machine-readable versions and any supplementary material and additional information, are available at https://lhcbproject.web.cern.ch/Publications/p/LHCb-PAPER-2023-030.html (LHCb public pages

Fraction of χc\chi_c decays in prompt J/ψJ/\psi production measured in pPb collisions at sNN=8.16\sqrt{s_{NN}}=8.16 TeV

The fraction of $\chi_{c1}$ and $\chi_{c2}$ decays in the prompt $J/\psi$ yield, $F_{\chi c}=\sigma_{\chi_c \to J/\psi}/\sigma_{J/\psi}$, is measured by the LHCb detector in pPb collisions at $\sqrt{s_{NN}}=8.16$ TeV. The study covers the forward ($1.5<y^*<4.0$) and backward ($-5.0<y^*<-2.5$) rapidity regions, where $y^*$ is the $J/\psi$ rapidity in the nucleon-nucleon center-of-mass system. Forward and backward rapidity samples correspond to integrated luminosities of 13.6 $\pm$ 0.3 nb$^{-1}$ and 20.8 $\pm$ 0.5 nb$^{-1}$, respectively. The result is presented as a function of the $J/\psi$ transverse momentum $p_{T,J/\psi}$ in the range 1$<p_{T, J/\psi}<20$ GeV/$c$. The $F_{\chi c}$ fraction at forward rapidity is compatible with the LHCb measurement performed in $pp$ collisions at $\sqrt{s}=7$ TeV, whereas the result at backward rapidity is 2.4 $\sigma$ larger than in the forward region for $1<p_{T, J/\psi}<3$ GeV/$c$. The increase of $F_{\chi c}$ at low $p_{T, J/\psi}$ at backward rapidity is compatible with the suppression of the $\psi$(2S) contribution to the prompt $J/\psi$ yield. The lack of in-medium dissociation of $\chi_c$ states observed in this study sets an upper limit of 180 MeV on the free energy available in these pPb collisions to dissociate or inhibit charmonium state formation.Comment: All figures and tables, along with machine-readable versions and any supplementary material and additional information, are available at https://cern.ch/lhcbproject/Publications/p/LHCb-PAPER-2023-028.html (LHCb public pages

Amplitude analysis of the Λb0→pK−γ decay

The resonant structure of the radiative decay Λb0→pK−γ in the region of proton-kaon invariant-mass up to 2.5 GeV/c2 is studied using proton-proton collision data recorded at centre-of-mass energies of 7, 8, and 13 TeV collected with the LHCb detector, corresponding to a total integrated luminosity of 9 fb−1. Results are given in terms of fit and interference fractions between the different components contributing to this final state. Only Λ resonances decaying to pK− are found to be relevant, where the largest contributions stem from the Λ(1520), Λ(1600), Λ(1800), and Λ(1890) states

Enhanced production of Λb0\Lambda_{b}^{0} baryons in high-multiplicity pppp collisions at s=13\sqrt{s} = 13 TeV

The production rate of $\Lambda_{b}^{0}$ baryons relative to $B^{0}$ mesons in $pp$ collisions at a center-of-mass energy $\sqrt{s} = 13$ TeV is measured by the LHCb experiment. The ratio of $\Lambda_{b}^{0}$ to $B^{0}$ production cross-sections shows a significant dependence on both the transverse momentum and the measured charged-particle multiplicity. At low multiplicity, the ratio measured at LHCb is consistent with the value measured in $e^{+}e^{-}$ collisions, and increases by a factor of $\sim2$ with increasing multiplicity. At relatively low transverse momentum, the ratio of $\Lambda_{b}^{0}$ to $B^{0}$ cross-sections is higher than what is measured in $e^{+}e^{-}$ collisions, but converges with the $e^{+}e^{-}$ ratio as the momentum increases. These results imply that the evolution of heavy $b$ quarks into final-state hadrons is influenced by the density of the hadronic environment produced in the collision. Comparisons with a statistical hadronization model and implications for the mechanisms enforcing quark confinement are discussed.Comment: All figures and tables, along with machine-readable versions and any supplementary material and additional information, are available at https://cern.ch/lhcbproject/Publications/p/LHCb-PAPER-2023-027.html (LHCb public pages

Helium identification with LHCb

The identification of helium nuclei at LHCb is achieved using a method based on measurements of ionisation losses in the silicon sensors and timing measurements in the Outer Tracker drift tubes. The background from photon conversions is reduced using the RICH detectors and an isolation requirement. The method is developed using $pp$ collision data at $\sqrt{s}=13\,{\rm TeV}$ recorded by the LHCb experiment in the years 2016 to 2018, corresponding to an integrated luminosity of $5.5\,{\rm fb}^{-1}$. A total of around $10^5$ helium and antihelium candidates are identified with negligible background contamination. The helium identification efficiency is estimated to be approximately $50\%$ with a corresponding background rejection rate of up to $\mathcal O(10^{12})$. These results demonstrate the feasibility of a rich programme of measurements of QCD and astrophysics interest involving light nuclei.Comment: All figures and tables, along with any supplementary material and additional information, are available at https://cern.ch/lhcbproject/Publications/p/LHCb-DP-2023-002.html (LHCb public pages

Observation of the decays B(s)0→Ds1(2536)∓K±B_{(s)}^{0}\to D_{s1}(2536)^{\mp}K^{\pm}

This paper reports the observation of the decays $B_{(s)}^{0}\to D_{s1}(2536)^{\mp}K^{\pm}$ using proton-proton collision data collected by the LHCb experiment, corresponding to an integrated luminosity of $9\,\mathrm{fb}^{-1}$. The branching fractions of these decays are measured relative to the normalisation channel $B^{0}\to \overline{D}^{0}K^{+}K^{-}$. The $D_{s1}(2536)^{-}$ meson is reconstructed in the $\overline{D}^{*}(2007)^{0}K^{-}$ decay channel and the products of branching fractions are measured to be $$\mathcal{B}(B_{s}^{0}\to D_{s1}(2536)^{\mp}K^{\pm})\times\mathcal{B}(D_{s1}(2536)^{-}\to\overline{D}^{*}(2007)^{0}K^{-})=(2.49\pm0.11\pm0.12\pm0.25\pm0.06)\times 10^{-5},$$ $$\mathcal{B}(B^{0}\to D_{s1}(2536)^{\mp}K^{\pm})\times\mathcal{B}(D_{s1}(2536)^{-}\to\overline{D}^{*}(2007)^{0}K^{-}) = (0.510\pm0.021\pm0.036\pm0.050)\times 10^{-5}.$$ The first uncertainty is statistical, the second systematic, and the third arises from the uncertainty of the branching fraction of the $B^{0}\to \overline{D}^{0}K^{+}K^{-}$ normalisation channel. The last uncertainty in the $B_{s}^{0}$ result is due to the limited knowledge of the fragmentation fraction ratio, $f_{s}/f_{d}$. The significance for the $B_{s}^{0}$ and $B^{0}$ signals is larger than $10\,\sigma$. The ratio of the helicity amplitudes which governs the angular distribution of the $D_{s1}(2536)^{-}\to\overline{D}^{*}(2007)^{0}K^{-}$ decay is determined from the data. The ratio of the $S$- and $D$-wave amplitudes is found to be $1.11\pm0.15\pm 0.06$ and its phase $0.70\pm0.09\pm 0.04$ rad, where the first uncertainty is statistical and the second systematic.Comment: All figures and tables, along with machine-readable versions and any supplementary material and additional information, are available at https://cern.ch/lhcbproject/Publications/p/LHCb-PAPER-2023-014.html (LHCb public pages