Barcelona Biomedical Research Park

In 2008 we published the first set of guidelines for standardizing research in autophagy. Since then, research on this topic has continued to accelerate, and many new scientists have entered the field. Our knowledge base and relevant new technologies have also been expanding. Accordingly, it is important to update these guidelines for monitoring autophagy in different organisms. Various reviews have described the range of assays that have been used for this purpose. Nevertheless, there continues to be confusion regarding acceptable methods to measure autophagy, especially in multicellular eukaryotes. For example, a key point that needs to be emphasized is thatthere is a difference between measurements that monitor the numbers or volume of autophagic elements (e.g., autophagosomes or autolysosomes) at any stage of the autophagic process versus those that measure flux through the autophagy pathway (i.e., the completeprocess including the amount and rate of cargo sequestered and degraded). In particular, a block in macroautophagy that results in autophagosome accumulation must be differentiated from stimuli that increase autophagic activity, defined as increasedautophagy induction coupled with increased delivery to, and degradation within, lysosomes (inmost higher eukaryotes and some protists such as Dictyostelium) or the vacuole (in plants and fungi). In other words, it is especially important that investigators new to the field understand that the appearance of more autophagosomes does not necessarily equate with more autophagy. In fact, in manycases, autophagosomes accumulate because of a block in trafficking to lysosomes without a concomitant change in autophagosome biogenesis, whereas an increase in autolysosomes may reflect a reduction in degradative activity. It is worth emphasizing here that lysosomal digestion is a stage of autophagy and evaluating its competence is a crucial part of the evaluation of autophagic flux, or complete autophagy. Here, we present a set of guidelines for the selection and interpretation of methods for use by investigators who aim to examine macroautophagy and related processes, as well as forreviewers who need to provide realistic and reasonable critiques of papers that are focused on these processes. These guidelines are not meant to be a formulaic set of rules, because the appropriate assays depend in part on the question being asked and the system being used. In addition, we emphasize that no individual assay is guaranteed to be the most appropriate one in every situation, and we strongly recommend the use of multipleassays to monitor autophagy. Along these lines, because of the potential for pleiotropic effects due to blocking autophagy through genetic manipulation, it is imperative to target by gene knockout or RNA interference more than one autophagyrelated protein. In addition, some individual Atg proteins, or groups of proteins, are involved in other cellular pathways implying that not all Atg proteins can be used as a specific marker for an autophagic process. In these guidelines, we consider these various methods of assessing autophagy and what information can, or cannot, be obtained from them. Finally, by discussing the merits and limits of particular assays, we hope to encourage technical innovation in the field.

The information about the genetic basis of human diseases lies at the heart of precision medicine and drug discovery. However, to realize its full potential to support these goals, several problems, such as fragmentation, heterogeneity, availability and different conceptualization of the data must be overcome. To provide the community with a resource free of these hurdles, we have developed DisGeNET (http://www.disgenet.org), one of the largest available collections of genes and variants involved in human diseases. DisGeNET integrates data from expert curated repositories, GWAS catalogues, animal models and the scientific literature. DisGeNET data are homogeneously annotated with controlled vocabularies and community-driven ontologies. Additionally, several original metrics are provided to assist the prioritization of genotype-phenotype relationships. The information is accessible through a web interface, a Cytoscape App, an RDF SPARQL endpoint, scripts in several programming languages and an R package. DisGeNET is a versatile platform that can be used for different research purposes including the investigation of the molecular underpinnings of specific human diseases and their comorbidities, the analysis of the properties of disease genes, the generation of hypothesis on drug therapeutic action and drug adverse effects, the validation of computationally predicted disease genes and the evaluation of text-mining methods performance.

In the genetic code, UGA serves as a stop signal and a selenocysteine codon, but no computational methods for identifying its coding function are available. Consequently, most selenoprotein genes are misannotated. We identified selenoprotein genes in sequenced mammalian genomes by methods that rely on identification of selenocysteine insertion RNA structures, the coding potential of UGA codons, and the presence of cysteine-containing homologs. The human selenoproteome consists of 25 selenoproteins.

Many common variants have been associated with hematological traits, but identification of causal genes and pathways has proven challenging. We performed a genome-wide association analysis in the UK Biobank and INTERVAL studies, testing 29.5 million genetic variants for association with 36 red cell, white cell, and platelet properties in 173,480 European-ancestry participants. This effort yielded hundreds of low frequency (<5%) and rare (<1%) variants with a strong impact on blood cell phenotypes. Our data highlight general properties of the allelic architecture of complex traits, including the proportion of the heritable component of each blood trait explained by the polygenic signal across different genome regulatory domains. Finally, through Mendelian randomization, we provide evidence of shared genetic pathways linking blood cell indices with complex pathologies, including autoimmune diseases, schizophrenia, and coronary heart disease and evidence suggesting previously reported population associations between blood cell indices and cardiovascular disease may be non-causal.

Many studies conducted during the last decade suggest the mental health benefits of green and blue spaces. We aimed to systematically review the available literature on the long-term mental health benefits of residential green and blue spaces by including studies that used standardized tools or objective measures of both the exposures and the outcomes of interest. We followed the PRISMA statement guidelines for reporting systematic reviews and meta-analysis. In total 28 studies were included in the systematic review. We found limited evidence for a causal relationship between surrounding greenness and mental health in adults, whereas the evidence was inadequate in children. The evidence was also inadequate for the other exposures evaluated (access to green spaces, quality of green spaces, and blue spaces) in both adults and children. The main limitation was the limited number of studies, together with the heterogeneity regarding exposure assessment. Given the increase in mental health problems and the current rapid urbanization worldwide, results of the present systematic review should be taken into account in future urban planning. However, further research is needed to provide more consistent evidence and more detailed information on the mechanisms and the characteristics of the green and blue spaces that promote better mental health. We provide recommendations for future studies in order to provide consistent and evidence-based recommendations for policy makers.

A goal in the kinetic characterization of a macromolecular system is the description of its slow relaxation processes via (i) identification of the structural changes involved in these processes and (ii) estimation of the rates or timescales at which these slow processes occur. Most of the approaches to this task, including Markov models, master-equation models, and kinetic network models, start by discretizing the high-dimensional state space and then characterize relaxation processes in terms of the eigenvectors and eigenvalues of a discrete transition matrix. The practical success of such an approach depends very much on the ability to finely discretize the slow order parameters. How can this task be achieved in a high-dimensional configuration space without relying on subjective guesses of the slow order parameters? In this paper, we use the variational principle of conformation dynamics to derive an optimal way of identifying the "slow subspace" of a large set of prior order parameters - either generic internal coordinates or a user-defined set of parameters. Using a variational formulation of conformational dynamics, it is shown that an existing method-the time-lagged independent component analysis-provides the optional solution to this problem. In addition, optimal indicators-order parameters indicating the progress of the slow transitions and thus may serve as reaction coordinates-are readily identified. We demonstrate that the slow subspace is well suited to construct accurate kinetic models of two sets of molecular dynamics simulations, the 6-residue fluorescent peptide MR121-GSGSW and the 30-residue intrinsically disordered peptide kinase inducible domain (KID). The identified optimal indicators reveal the structural changes associated with the slow processes of the molecular system under analysis.

This is a revision of guidelines, originally published in 2004, for the assessment of patients with neuropathic pain. Neuropathic pain is defined as pain arising as a direct consequence of a lesion or disease affecting the somatosensory system either at peripheral or central level. Screening questionnaires are suitable for identifying potential patients with neuropathic pain, but further validation of them is needed for epidemiological purposes. Clinical examination, including accurate sensory examination, is the basis of neuropathic pain diagnosis. For more accurate sensory profiling, quantitative sensory testing is recommended for selected cases in clinic, including the diagnosis of small fiber neuropathies and for research purposes. Measurement of trigeminal reflexes mediated by A-beta fibers can be used to differentiate symptomatic trigeminal neuralgia from classical trigeminal neuralgia. Measurement of laser-evoked potentials is useful for assessing function of the A-delta fiber pathways in patients with neuropathic pain. Functional brain imaging is not currently useful for individual patients in clinical practice, but is an interesting research tool. Skin biopsy to measure the intraepidermal nerve fiber density should be performed in patients with clinical signs of small fiber dysfunction. The intensity of pain and treatment effect (both in clinic and trials) should be assessed with numerical rating scale or visual analog scale. For future neuropathic pain trials, pain relief scales, patient and clinician global impression of change, the proportion of responders (50% and 30% pain relief), validated neuropathic pain quality measures and assessment of sleep, mood, functional capacity and quality of life are recommended.

These guidelines are a consensus work of a considerable number of members of the immunology and flow cytometry community. They provide the theory and key practical aspects of flow cytometry enabling immunologists to avoid the common errors that often undermine immunological data. Notably, there are comprehensive sections of all major immune cell types with helpful Tables detailing phenotypes in murine and human cells. The latest flow cytometry techniques and applications are also described, featuring examples of the data that can be generated and, importantly, how the data can be analysed. Furthermore, there are sections detailing tips, tricks and pitfalls to avoid, all written and peer-reviewed by leading experts in the field, making this an essential research companion.

High-coverage sequencing of 79 (wild and captive) individuals representing all six non-human great ape species has identified over 88 million single nucleotide polymorphisms providing insight into ape genetic variation and evolutionary history and enabling comparison with human genetic diversity. In an effort to provide insights into great ape genetic variation, the authors sequence 79 wild- and captive-born individuals from across all six great ape species and seven subspecies. Their data and analyses shed light on population structure and gene flow, inbreeding, inferred dynamics of effective population sizes and the differences in the rate of gene loss among the great apes. This new catalogue of great ape genome diversity provides a valuable resource for evolutionary and conservation studies. Most great ape genetic variation remains uncharacterized1,2; however, its study is critical for understanding population history3,4,5,6, recombination7, selection8 and susceptibility to disease9,10. Here we sequence to high coverage a total of 79 wild- and captive-born individuals representing all six great ape species and seven subspecies and report 88.8 million single nucleotide polymorphisms. Our analysis provides support for genetically distinct populations within each species, signals of gene flow, and the split of common chimpanzees into two distinct groups: Nigeria–Cameroon/western and central/eastern populations. We find extensive inbreeding in almost all wild populations, with eastern gorillas being the most extreme. Inferred effective population sizes have varied radically over time in different lineages and this appears to have a profound effect on the genetic diversity at, or close to, genes in almost all species. We discover and assign 1,982 loss-of-function variants throughout the human and great ape lineages, determining that the rate of gene loss has not been different in the human branch compared to other internal branches in the great ape phylogeny. This comprehensive catalogue of great ape genome diversity provides a framework for understanding evolution and a resource for more effective management of wild and captive great ape populations.

Accurately predicting protein–ligand binding affinities is an important problem in computational chemistry since it can substantially accelerate drug discovery for virtual screening and lead optimization. We propose here a fast machine-learning approach for predicting binding affinities using state-of-the-art 3D-convolutional neural networks and compare this approach to other machine-learning and scoring methods using several diverse data sets. The results for the standard PDBbind (v.2016) core test-set are state-of-the-art with a Pearson’s correlation coefficient of 0.82 and a RMSE of 1.27 in p K units between experimental and predicted affinity, but accuracy is still very sensitive to the specific protein used. K DEEP is made available via PlayMolecule.org for users to test easily their own protein–ligand complexes, with each prediction taking a fraction of a second. We believe that the speed, performance, and ease of use of K DEEP makes it already an attractive scoring function for modern computational chemistry pipelines.

The high arithmetic performance and intrinsic parallelism of recent graphical processing units (GPUs) can offer a technological edge for molecular dynamics simulations. ACEMD is a production-class biomolecular dynamics (MD) engine supporting CHARMM and AMBER force fields. Designed specifically for GPUs it is able to achieve supercomputing scale performance of 40 ns/day for all-atom protein systems with over 23 000 atoms. We provide a validation and performance evaluation of the code and run a microsecond-long trajectory for an all-atom molecular system in explicit TIP3P water on a single workstation computer equipped with just 3 GPUs. We believe that microsecond time scale molecular dynamics on cost-effective hardware will have important methodological and scientific implications.

Over 70,000 excess deaths occurred in Europe during the summer of 2003. The resulting societal awareness led to the design and implementation of adaptation strategies to protect at-risk populations. We aimed to quantify heat-related mortality burden during the summer of 2022, the hottest season on record in Europe. We analyzed the Eurostat mortality database, which includes 45,184,044 counts of death from 823 contiguous regions in 35 European countries, representing the whole population of over 543 million people. We estimated 61,672 (95% confidence interval (CI) = 37,643-86,807) heat-related deaths in Europe between 30 May and 4 September 2022. Italy (18,010 deaths; 95% CI = 13,793-22,225), Spain (11,324; 95% CI = 7,908-14,880) and Germany (8,173; 95% CI = 5,374-11,018) had the highest summer heat-related mortality numbers, while Italy (295 deaths per million, 95% CI = 226-364), Greece (280, 95% CI = 201-355), Spain (237, 95% CI = 166-312) and Portugal (211, 95% CI = 162-255) had the highest heat-related mortality rates. Relative to population, we estimated 56% more heat-related deaths in women than men, with higher rates in men aged 0-64 (+41%) and 65-79 (+14%) years, and in women aged 80+ years (+27%). Our results call for a reevaluation and strengthening of existing heat surveillance platforms, prevention plans and long-term adaptation strategies.

PURPOSE: To develop a self-report version of the EQ-5D for younger respondents, named the EQ-5D-Y (Youth); to test its comprehensibility for children and adolescents and to compare results obtained using the standard adult EQ-5D and the EQ-5D-Y. METHODS: An international task force revised the content of EQ-5D and wording to ensure relevance and clarity for young respondents. Children's and adolescents' understanding of the EQ-5D-Y was tested in cognitive interviews after the instrument was translated into German, Italian, Spanish and Swedish. Differences between the EQ-5D and the EQ-5D-Y regarding frequencies of reported problems were investigated in Germany, Spain and South Africa. RESULTS: The content of the EQ-5D dimensions proved to be appropriate for the measurement of HRQOL in young respondents. The wording of the questionnaire had to be adapted which led to small changes in the meaning of some items and answer options. The adapted EQ-5D-Y was satisfactorily understood by children and adolescents in different countries. It was better accepted and proved more feasible than the EQ-5D. The administration of the EQ-5D and of the EQ-5D-Y causes differences in frequencies of reported problems. CONCLUSIONS: The newly developed EQ-5D-Y is a useful tool to measure HRQOL in young people in an age-appropriate manner.

BACKGROUND: Specific cortico-striato-thalamic circuits are hypothesised to mediate the symptoms of obsessive-compulsive disorder (OCD), but structural neuroimaging studies have been inconsistent. AIMS: To conduct a meta-analysis of published and unpublished voxel-based morphometry studies in OCD. METHOD: Twelve data-sets comprising 401 people with OCD and 376 healthy controls met inclusion criteria. A new improved voxel-based meta-analytic method, signed differential mapping (SDM), was developed to examine regions of increased and decreased grey matter volume in the OCD group v. control group. Results No between-group differences were found in global grey matter volumes. People with OCD had increased regional grey matter volumes in bilateral lenticular nuclei, extending to the caudate nuclei, as well as decreased volumes in bilateral dorsal medial frontal/anterior cingulate gyri. A descriptive analysis of quartiles, a sensitivity analysis as well as analyses of subgroups further confirmed these findings. Meta-regression analyses showed that studies that included individuals with more severe OCD were significantly more likely to report increased grey matter volumes in the basal ganglia. No effect of current antidepressant treatment was observed. Conclusions The results support a dorsal prefrontal-striatal model of the disorder and raise the question of whether functional alterations in other brain regions commonly associated with OCD, such as the orbitofrontal cortex, may reflect secondary compensatory strategies. Whether the reported differences between participants with OCD and controls precede the onset of the symptoms and whether they are specific to OCD remains to be established.

MOTIVATION: An important step in structure-based drug design consists in the prediction of druggable binding sites. Several algorithms for detecting binding cavities, those likely to bind to a small drug compound, have been developed over the years by clever exploitation of geometric, chemical and evolutionary features of the protein. RESULTS: Here we present a novel knowledge-based approach that uses state-of-the-art convolutional neural networks, where the algorithm is learned by examples. In total, 7622 proteins from the scPDB database of binding sites have been evaluated using both a distance and a volumetric overlap approach. Our machine-learning based method demonstrates superior performance to two other competitive algorithmic strategies. AVAILABILITY AND IMPLEMENTATION: DeepSite is freely available at www.playmolecule.org. Users can submit either a PDB ID or PDB file for pocket detection to our NVIDIA GPU-equipped servers through a WebGL graphical interface. CONTACT: gianni.defabritiis@upf.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Summary Recent efforts to understand how the patterning of interaction strength affects both structure and dynamics in food webs have highlighted several obstacles to productive synthesis. Issues arise with respect to goals and driving questions, methods and approaches, and placing results in the context of broader ecological theory. Much confusion stems from lack of clarity about whether the questions posed relate to community‐level patterns or to species dynamics, and to what authors actually mean by the term ‘interaction strength’. Here, we describe the various ways in which this term has been applied and discuss the implications of loose terminology and definition for the development of this field. Of particular concern is the clear gap between theoretical and empirical investigations of interaction strengths and food web dynamics. The ecological community urgently needs to explore new ways to estimate biologically reasonable model coefficients from empirical data, such as foraging rates, body size, metabolic rate, biomass distribution and other species traits. Combining numerical and analytical modelling approaches should allow exploration of the conditions under which different interaction strengths metrics are interchangeable with regard to relative magnitude, system responses, and species identity. Finally, the prime focus on predator–prey links in much of the research to date on interaction strengths in food webs has meant that the potential significance of non‐trophic interactions, such as competition, facilitation and biotic disturbance, has been largely ignored by the food web community. Such interactions may be important dynamically and should be routinely included in future food web research programmes.

T cells) from up to 197 individuals. We assess, quantitatively, the relative contribution of cis-genetic and epigenetic factors to transcription and evaluate their impact as potential sources of confounding in epigenome-wide association studies. Further, we characterize highly coordinated genetic effects on gene expression, methylation, and histone variation through quantitative trait locus (QTL) mapping and allele-specific (AS) analyses. Finally, we demonstrate colocalization of molecular trait QTLs at 345 unique immune disease loci. This expansive, high-resolution atlas of multi-omics changes yields insights into cell-type-specific correlation between diverse genomic inputs, more generalizable correlations between these inputs, and defines molecular events that may underpin complex disease risk.

A goal in the kinetic characterization of a macromolecular system is the description of its slow relaxation processes, involving (i) identification of the structural changes involved in these pro-cesses, and (ii) estimation of the rates or timescales at which these slow processes occur. Most of the approaches to this task, including Markov models, Master-equation models, and kinetic network models, start by discretizing the high-dimensional state space and then characterize re-laxation processes in terms of the eigenvectors and eigenvalues of a discrete transition matrix. The practical success of such an approach depends very much on the ability to finely discretize the slow order parameters. How can this task be achieved in a high-dimensional configuration space without relying on subjective guesses of the slow order parameters? In this paper, we use the variational principle of conformation dynamics to derive an optimal way of identifying the “slow subspace ” of a large set of prior order parameters- either generic internal coordinates (distances and dihedral angles), or a user-defined set of parameters. It is shown that a method to identify this slow subspace exists in statistics: the time-lagged independent component analysis (TICA). Furthermore, optimal indicators—order parameters indicating the progress of the slow transitions

A full genome sequence is presented of sugar beet Beta vulgaris, the first plant belonging to Caryophyllales to have its genome sequenced; spinach was sequenced to enable inter-clade comparisons, and intraspecific variation was analysed by comparative genomics of a progenitor of all beet crops and additional sugar beet accessions. Industrial production of sugar from sugar beet (Beta vulgaris) began in Europe in the early nineteenth century, and in the intervening 200 years the sugar content of the commonly used cultivars has increased from 8% to 18%. A high-quality reference genome sequence for sugar beet is published in this issue, together with that of the related spinach plant (Spinacia oleracea) and assembled genomes from four additional sugar beet breeding lines. Information held in these genome sequences will be useful for the characterization of genes involved in sugar production and identification of targets for breeding efforts, as well as towards its application as a sustainable energy crop. Sugar beet (Beta vulgaris ssp. vulgaris) is an important crop of temperate climates which provides nearly 30% of the world’s annual sugar production and is a source for bioethanol and animal feed. The species belongs to the order of Caryophylalles, is diploid with 2n = 18 chromosomes, has an estimated genome size of 714–758 megabases1 and shares an ancient genome triplication with other eudicot plants2. Leafy beets have been cultivated since Roman times, but sugar beet is one of the most recently domesticated crops. It arose in the late eighteenth century when lines accumulating sugar in the storage root were selected from crosses made with chard and fodder beet3. Here we present a reference genome sequence for sugar beet as the first non-rosid, non-asterid eudicot genome, advancing comparative genomics and phylogenetic reconstructions. The genome sequence comprises 567 megabases, of which 85% could be assigned to chromosomes. The assembly covers a large proportion of the repetitive sequence content that was estimated4 to be 63%. We predicted 27,421 protein-coding genes supported by transcript data and annotated them on the basis of sequence homology. Phylogenetic analyses provided evidence for the separation of Caryophyllales before the split of asterids and rosids, and revealed lineage-specific gene family expansions and losses. We sequenced spinach (Spinacia oleracea), another Caryophyllales species, and validated features that separate this clade from rosids and asterids. Intraspecific genomic variation was analysed based on the genome sequences of sea beet (Beta vulgaris ssp. maritima; progenitor of all beet crops) and four additional sugar beet accessions. We identified seven million variant positions in the reference genome, and also large regions of low variability, indicating artificial selection. The sugar beet genome sequence enables the identification of genes affecting agronomically relevant traits, supports molecular breeding and maximizes the plant’s potential in energy biotechnology.

The understanding of protein-ligand binding is of critical importance for biomedical research, yet the process itself has been very difficult to study because of its intrinsically dynamic character. Here, we have been able to quantitatively reconstruct the complete binding process of the enzyme-inhibitor complex trypsin-benzamidine by performing 495 molecular dynamics simulations of free ligand binding of 100 ns each, 187 of which produced binding events with an rmsd less than 2 Å compared to the crystal structure. The binding paths obtained are able to capture the kinetic pathway of the inhibitor diffusing from solvent (S0) to the bound (S4) state passing through two metastable intermediate states S2 and S3. Rather than directly entering the binding pocket the inhibitor appears to roll on the surface of the protein in its transition between S3 and the final binding pocket, whereas the transition between S2 and the bound pose requires rediffusion to S3. An estimation of the standard free energy of binding gives ΔG° = -5.2 ± 0.4 kcal/mol (cf. the experimental value -6.2 kcal/mol), and a two-states kinetic model k(on) = (1.5 ± 0.2) × 10(8) M(-1) s(-1) and k(off) = (9.5 ± 3.3) × 10(4) s(-1) for unbound to bound transitions. The ability to reconstruct by simple diffusion the binding pathway of an enzyme-inhibitor binding process demonstrates the predictive power of unconventional high-throughput molecular simulations. Moreover, the methodology is directly applicable to other molecular systems and thus of general interest in biomedical and pharmaceutical research.