Miracle-dz/UniFAR-FTUs huggingface dataset:
ATLASx: a computational map for the exploration of biochemical space

Metabolic “dark matter” describes currently unknown metabolic processes, which form a blind spot in our general understanding of metabolism and slow down the development of biosynthetic cell factories and naturally derived pharmaceuticals. Mapping the dark matter of metabolism remains an open challenge that can be addressed globally and systematically by existing computational solutions. In this work, we use 490 generalized enzymatic reaction rules to map both known and unknown metabolic processes around a biochemical database of 1.5 million biological compounds. We predict over 5 million reactions and integrate nearly 2 million naturally and synthetically-derived compounds into the global network of biochemical knowledge, named ATLASx. ATLASx is available to researchers as a powerful online platform that supports the prediction and analysis of novel biochemical pathways and evaluates the biochemical vicinity of molecule classes (https://lcsbdatabases.epfl.ch/Atlas2).
Disentangling the relative importance of spatio-temporal parameters and host selectivity in shaping AMF communities in temperate forests

Many woody and herbaceous plants in temperate forests cannot establish and survive in the absence of mycorrhizal associations. Here, we address and hierarchize how Glomeromycota, a group of soil-borne fungi, forming a ubiquitous type of mycorrhiza, varied in a temperate forest patch in Germany with time, space, plant hosts but also the proximity to Glomeromycota-associating woody species. The communities of Glomeromycota in our study were non-random. We observed that space had a greater impact on fungal community structure than either time or proximity to Glomeromycota-associating trees but unlike host identity did not alter Glomeromycotan richness. The set of parameters which we addressed has rarely been studied together and we believe that the resulting ranking could ease prioritizing some of them to include in future surveys. Glomeromycota are crucial for the establishment of understory plants in temperate forests making it desirable to further explore how they vary in time and space.
Exceptional origin activation revealed by comparative analysis in two laboratory yeast strains

We performed a comparative analysis of replication origin activation by genome-wide single-stranded DNA mapping in two common laboratory strains of Saccharomyces cerevisiae challenged by hydroxyurea (HU), an inhibitor of the ribonucleotide reductase. By doing so we gained understanding of the impact on origin activation by three factors: replication checkpoint control, DNA sequence polymorphisms, and relative positioning of origin and transcription unit. Our analysis recapitulated the previous finding that the majority of origins are subject to checkpoint control by the Rad53 kinase when cells were treated with HU. In addition, origins either subject to Rad53 checkpoint control or impervious to it are largely concordant between the two strains. However, these two strains also produced different dynamics of origin activation. First, the W303-RAD53 cells showed a significant reduction of fork progression than A364a-RAD53 cells. This phenotype was accompanied by an elevated level of Rad53 phosphorylation in W303-RAD53 cells. Second, W303-rad53K227A checkpoint-deficient cells activated a greater number of origins accompanied by global reduction of ssDNA across all origins compared to A364a-rad53K227A cells; and this is correlated with lower expression level of the mutant protein in W303 than in A364a. We also show that sequence polymorphism in the consensus motifs of the replication origins plays a minor role in determining origin usage. Remarkably, eight strain-specific origins lack the canonical 11-bp consensus motif for autonomously replicating sequences in either strain background. Finally, we identified a new class of origins that are only active in checkpoint-proficient cells, which we named “Rad53-dependent origins”. The only discernible feature of these origins is that they tend to overlap with an open reading frame, suggesting previously unexplored connection between transcription and origin activation. Our study presents a comprehensive list of origin usage in two diverse yeast genetic backgrounds, fine-tunes the different categories of origins with respect to checkpoint control, and provokes further exploration of the interplay between origin activation and transcription. Author Summary Comparative analysis of origins of replication in two laboratory yeast strains reveals new insights into origin activation, regulation and dependency on the Rad53 checkpoint kinase.
massiveGenesetsTest: a web tool to run enrichment analysis

Motivation Inferring biological phenotypes from genomic data and sample clusters is a routinely task usually performed with Gene-Set Enrichment Analysis (GSEA), a tool that queries gene-profiles. In previous work, we scrutinized the approach based on Mann-Witney for Gene-Sets test. We highlighted the Mann-Witney test-statistics sensitivity to uncover weak signals and the drastic decreasing of time complexity. Results We propose web implementation of the Gene-sets testing based on the Mann-Witney procedure. The test-procedure has reshaped to decrease the computational expense, now about tens of seconds, even if a large collection of gene-sets queries the same gene-profile. The probabilistic interpretation of the normalized test-statistic has been investigated to a better understanding of the enrichment results. A novel prioritization method across enrichment-scores, gene-set dimensions, and p-values, draws attention to relevant gene-sets. The web tool provides both tabular and graphical enrichment results. A complimentary R function allows integrating the enrichment procedure in a complex context. Contact pagnotta@unisannio.it Supplementary information Example data and guidelines are included in the supporting material of the web-site.
Fluorescence Imaging-Based Discovery of Membrane Domain-Associated Proteins in Mycobacterium smegmatis

The intracellular membrane domain (IMD) is a membrane subcompartment found in Mycobacterium smegmatis cells. Proteomic analysis of purified IMD identified more than 300 proteins, including enzymes involved in cell envelope biosynthesis. ABSTRACT Mycobacteria spatially organize their plasma membrane, and many enzymes involved in envelope biosynthesis associate with a membrane compartment termed the intracellular membrane domain (IMD). The IMD is concentrated in the polar regions of growing cells and becomes less polarized under nongrowing conditions. Because mycobacteria elongate from the poles, the observed polar localization of the IMD during growth likely supports the localized biosynthesis of envelope components. While we have identified more than 300 IMD-associated proteins by proteomic analyses, only a few of these have been verified by independent experimental methods. Furthermore, some IMD-associated proteins may have escaped proteomic identification and remain to be identified. Here, we visually screened an arrayed library of 523 Mycobacterium smegmatis strains, each producing a Dendra2-FLAG-tagged recombinant protein. We identified 29 fusion proteins that showed polar fluorescence patterns characteristic of IMD proteins. Twenty of these had previously been suggested to localize to the IMD based on proteomic data. Of the nine remaining IMD candidate proteins, three were confirmed by biochemical methods to be associated with the IMD. Taken together, this new colocalization strategy is effective in verifying the IMD association of proteins found by proteomic analyses while facilitating the discovery of additional IMD-associated proteins. IMPORTANCE The intracellular membrane domain (IMD) is a membrane subcompartment found in Mycobacterium smegmatis cells. Proteomic analysis of purified IMD identified more than 300 proteins, including enzymes involved in cell envelope biosynthesis. However, proteomics on its own is unlikely to detect every IMD-associated protein because of technical and biological limitations. Here, we describe fluorescent protein colocalization as an alternative, independent approach. Using a combination of fluorescence microscopy, proteomics, and subcellular fractionation, we identified three new proteins associated with the IMD. Such a robust method to rigorously define IMD proteins will benefit future investigations to decipher the synthesis, maintenance, and functions of this membrane domain and help delineate a more general mechanism of subcellular protein localization in mycobacteria.
Inferring Population Histories for Ancient Genomes Using Genome-Wide Genealogies

Ancient genomes anchor genealogies in directly observed historical genetic variation, and contextualise ancestral lineages with archaeological insights into their geography and lifestyles. We introduce an extension of the Relate algorithm to incorporate ancient genomes and reconstruct the joint genealogies of 14 previously published high-coverage ancients and 278 present-day individuals of the Simons Genome Diversity Project. As the majority of ancient genomes are of lower coverage and cannot be directly built into genealogies, we additionally present a fast and scalable method, Colate, for inferring coalescence rates between low-coverage genomes without requiring phasing or imputation. Our method leverages sharing patterns of mutations dated using a genealogy to construct a likelihood, which is maximised using an expectation-maximisation algorithm. We apply Colate to 430 ancient human shotgun genomes of >0.5x mean coverage. Using Relate and Colate, we characterise dynamic population structure, such as repeated partial population replacements in Ireland, and gene-flow between early farmer and European hunter-gatherer groups. We further show that the previously reported increase in the TCC/TTC mutation rate, which is strongest in West Eurasians among present-day people, was already widespread across West Eurasia in the Late Glacial Period ~10k - 15k years ago, is strongest in Neolithic and Anatolian farmers, and is remarkably well predicted by the coalescence rates between other genomes and a 10,000-year-old Anatolian individual. This suggests that the driver of this signal originated in ancestors of ancient Anatolia >14k years ago, but was already absent by the Mesolithic and may indicate a genetic link between the Near East and European hunter-gatherer groups in the Late Paleolithic.
Functional cooperativity mediated by rationally selected combinations of human monoclonal antibodies targeting the henipavirus receptor binding protein

Hendra virus (HeV) and Nipah virus (NiV), the prototypic members of the Henipavirus (HNV) genus, are emerging, zoonotic paramyxoviruses known to cause severe disease across six mammalian orders, including humans (Eaton et al., 2006). While several research groups have made strides in developing candidate vaccines and therapeutics against henipaviruses, such countermeasures have not been licensed for human use, and significant gaps in knowledge about the human immune response to these viruses exist. To address these gaps, we isolated a large panel of human monoclonal antibodies (mAbs) from the B cells of an individual with prior occupation-related exposure to the equine HeV vaccine (Equivac® HeV). Competition-binding and hydrogen-deuterium exchange mass spectrometry (HDX-MS) studies identified at least six distinct antigenic sites on the HeV/NiV receptor binding protein (RBP) that are recognized by human mAbs. Antibodies recognizing multiple antigenic sites potently neutralized NiV and/or HeV isolates in vitro. The most potent class of cross-reactive antibodies achieved neutralization by blocking viral attachment to the host cell receptors ephrin-B2 and ephrin-B3. Antibodies from this class mimic receptor binding by inducing a receptor-bound conformation to the HeV-RBP protein tetramer, exposing an epitope that appears to lie hidden in the interface between protomers within the HeV-RBP tetramer. Antibodies that recognize this cryptic epitope potently neutralized HeV and NiV. Flow cytometric studies using cell-surface-displayed HeV-RBP protein showed that cross-reactive, neutralizing mAbs from each of these classes cooperate for binding. In a highly stringent hamster model of NiVB infection, antibodies from both classes reduced morbidity and mortality and achieved synergistic protection in combination and provided therapeutic benefit when combined into two bispecific platforms. These studies identified multiple candidate mAbs that might be suitable for use in a cocktail therapeutic approach to achieve synergistic antiviral potency and reduce the risk of virus escape during treatment.
Effect of High-Flow Oxygen on Exercise Performance in COPD Patients. Randomized Trial

Background: High-flow oxygen therapy (HFOT) provides oxygen-enriched, humidified, and heated air at high flow rates via nasal cannula. It could be an alternative to low-flow oxygen therapy (LFOT) which is commonly used by patients with chronic obstructive pulmonary disease (COPD) during exercise training. Research Question: We evaluated the hypothesis that HFOT improves exercise endurance in COPD patients compared to LFOT. Methods: Patients with stable COPD, FEV1 40–80% predicted, resting pulse oximetry (SpO2) ≥92%, performed two constant-load cycling exercise tests to exhaustion at 75% of maximal work rate on two different days, using LFOT (3 L/min) and HFOT (60 L/min, FiO2 0.45) in randomized order according to a crossover design. Primary outcome was exercise endurance time, further outcomes were SpO2, breath rate and dyspnea. Results: In 79 randomized patients, mean ± SD age 58 ± 9 y, FEV1 63 ± 9% predicted, GOLD grades 2-3, resting PaO2 9.4 ± 1.0 kPa, intention-to-treat analysis revealed an endurance time of 688 ± 463 s with LFOT and 773 ± 471 s with HFOT, mean difference 85 s (95% CI: 7 to 164, P = 0.034), relative increase of 13% (95% CI: 1 to 28). At isotime, patients had lower respiratory rate and higher SpO2 with HFOT. At end-exercise, SpO2 was higher by 2% (95% CI: 2 to 2), and Borg CR10 dyspnea scores were lower by 0.8 points (95% CI: 0.3 to 1.2) compared to LFOT. Interpretation: In mildly hypoxemic patients with COPD, HFOT improved endurance time in association with higher arterial oxygen saturation, reduced respiratory rate and less dyspnea compared to LFOT. Therefore, HFOT is promising for enhancing exercise performance in COPD. Clinical Trial Registration: www.ClinicalTrials. gov, identifier: NCT03955770.
Paths to Oblivion: Common Neural Mechanisms of Anaesthesia and Disorders of Consciousness

The human brain generates a rich repertoire of spatiotemporal dynamics during normal wakefulness, supporting a wide variety of conscious experiences and cognitive functions. However, neural dynamics are reconfigured, in comparable ways, when consciousness is lost either due to anaesthesia or disorders of consciousness (DOC). Here, leveraging a neurobiologically realistic whole-brain computational model informed by functional MRI, diffusion MRI, and PET, we sought to identify the neurobiological mechanisms that explain the common reconfiguration of neural dynamics observed both for transient pharmacological intervention and chronic neuroanatomical injury. Our results show that, by incorporating local inhibitory action through a PET-based GABA receptor density map, our model can reproduce the brain dynamics of subjects undergoing propofol anaesthesia, and that this effect depends specifically on the spatial distribution of GABA receptors across cortical regions. Additionally, using a structural connectome obtained from DOC patients, we demonstrate how the dynamics that characterise loss of consciousness can emerge from changes in neuroanatomical connectivity. Crucially, we find that each of these two interventions generalises across datasets: a model with increased GABA-mediated inhibition can reproduce the dynamics of DOC patients’ brains, and a model with a DOC connectome is also compatible with brain dynamics observed during propofol anaesthesia. These results demonstrate how increased inhibition and connectome randomisation represent different neurobiological paths towards the characteristic dynamics of the unconscious brain. Overall, the present findings begin to disentangle the neurobiological mechanisms by which highly dissimilar perturbations of the brain’s neurodynamics can lead to unconsciousness.
Hippocampal sharp wave-ripples and the associated sequence replay emerge from structured synaptic interactions in a network model of area CA3

Hippocampal place cells are activated sequentially as an animal explores its environment. These activity sequences are internally recreated (“replayed”), either in the same or reversed order, during bursts of activity (sharp wave-ripples; SWRs) that occur in sleep and awake rest. SWR-associated replay is thought to be critical for the creation and maintenance of long-term memory. In order to identify the cellular and network mechanisms of SWRs and replay, we constructed and simulated a data-driven model of area CA3 of the hippocampus. Our results show that the chain-like structure of recurrent excitatory interactions established during learning not only determines the content of replay, but is essential for the generation of the SWRs as well. We find that bidirectional replay requires the interplay of the experimentally confirmed, temporally symmetric plasticity rule, and cellular adaptation. Our model provides a unifying framework for diverse phenomena involving hippocampal plasticity, representations, and dynamics, and suggests that the structured neural codes induced by learning may have greater influence over cortical network states than previously appreciated.
Perceptual History Acts in World-Centred Coordinates

Serial dependence effects have been observed using a variety of stimuli and tasks, revealing that the recent past can bias current percepts, leading to increased similarity between two. The aim of this study is to determine whether this temporal integration occurs in egocentric or allocentric coordinates. We asked participants to perform an orientation reproduction task using grating stimuli while the head was kept at a fixed position, or after a 40° yaw rotation between trials, from left (−20°) to right (+20°), putting the egocentric and allocentric cues in conflict. Under these conditions, allocentric cues prevailed.
Host preferences of coexisting Perkinsea parasitoids during coastal dinoflagellate blooms

Parasites in aquatic systems are highly diverse and ubiquitous. In marine environments, parasite-host interactions contribute substantially to shaping microbial communities, but their nature and complexity remain poorly understood. In this study, we examined the relationship between Perkinsea parasitoids and bloom-forming dinoflagellate species. Our aim was to determine whether parasite-host species interactions are specific and whether the diversity and distribution of parasitoids are shaped by their dinoflagellate hosts. Several locations along the Catalan coast (NW Mediterranean Sea) were sampled during the blooms of five dinoflagellate species and the diversity of Perkinsea was determined by combining cultivation-based methods with metabarcoding of the V4 region of 18S rDNA. Most known species of Parviluciferaceae, and others not yet described, were detected, some of them coexisting in the same coastal location, and with a wide distribution. The specific parasite-host interactions determined for each of the studied blooms demonstrated the host preferences exhibited by parasitoids in nature. The dominance of a species within the parasitoid community is driven by the presence and abundances of its preferred host(s). The absence of parasitoid species, often associated with a low abundance of their preferred hosts, suggested that high infection rates are reached only under conditions that favour parasitoid propagation, especially dinoflagellate blooms.
Newly developed ad hoc molecular assays show how eDNA can witness and anticipate the monk seal recolonization of central Mediterranean

The monk seal Monachus monachus is the most endangered pinniped worldwide and the only one found in the Mediterranean, where its distribution and abundance have suffered a drastic decline in the last few decades. Data on its status are scattered due to both its rarity and evasiveness, and records are biased towards occasional, mostly coastal, encounters. Nowadays molecular techniques allow us to detect and quantify minute amounts of DNA traces released in the environment (eDNA) by any organism. We present three qPCR-assays targeting the monk seal mitogenome. The assays were soundly tested on an extensive and diversified sample set (n=73), including positive controls from Madeira breeding population collected during the peak of abundance, and two opportunistic Mediterranean eDNA-sample collections (offshore/coastal) from on-going projects. Monk seal DNA was detected in 47.2% and 66.7% of the samples collected in the Tyrrhenian from a ferry platform (2018-2019) and in the Pelagie archipelago -Strait of Sicily- (2020) respectively, anticipating (up to 2 year) visual observations occurred subsequently in proximity of the sampled areas. In the Tyrrhenian, detection occurrence increased between 2018 and 2019. Monk seal DNA recoveries were commoner in night-time ferry-samples, suggesting nocturnal predatory activity in pelagic waters. The proposed technique provides a non-invasive and yet highly-sensitive tool for defining the monk seal actual distribution and home range, its recovery rate and pinpoint coastal/offshore localities where prioritizing conservation, research, citizen science and education initiatives.
HLA binding of self-peptides is biased towards proteins with specific molecular functions

Human leukocyte antigen (HLA) is highly polymorphic and plays a key role in guiding adaptive immune responses by presenting foreign and self peptides to T cells. Each HLA variant selects a minor fraction of peptides that match a certain motif required for optimal interaction with the peptide-binding groove. These restriction rules define the landscape of peptides presented to T cells. Given these limitations, one might suggest that the choice of peptides presented by HLA is non-random and there is preferential presentation of an array of peptides that is optimal for distinguishing self and foreign proteins. In this study we explore these preferences with a comparative analysis of self peptides enriched and depleted in HLA ligands. We show that HLAs exhibit preferences towards presenting peptides from certain proteins while disfavoring others with specific functions, and highlight differences between various HLA genes and alleles in those preferences. We link those differences to HLA anchor residue propensities and amino acid composition of preferentially presented proteins. The set of proteins that peptides presented by a given HLA are most likely to be derived from can be used to distinguish between class I and class II HLAs and HLA alleles. Our observations can be extrapolated to explain the protective effect of certain HLA alleles in infectious diseases, and we hypothesize that they can also explain susceptibility to certain autoimmune diseases and cancers. We demonstrate that these differences lead to differential presentation of HIV, influenza virus, SARS-CoV-1 and SARS-CoV-2 proteins by various HLA alleles. Finally, we show that the reported self peptidome preferences of distinct HLA variants can be compensated by combinations of HLA-A/HLA-B and HLA-A/HLA-C alleles in frequent haplotypes.