Real-time forecasting of COVID-19-related hospital strain in France using a non-Markovian mechanistic model.

Projects such as the European Covid-19 Forecast Hub publish forecasts on the national level for new deaths, new cases, and hospital admissions, but not direct measurements of hospital strain like critical care bed occupancy at the sub-national level, which is of particular interest to health professionals for planning purposes. We present a sub-national French framework for forecasting hospital strain based on a non-Markovian compartmental model, its associated online visualisation tool and a retrospective evaluation of the real-time forecasts it provided from January to December 2021 by comparing to three baselines derived from standard statistical forecasting methods (a naive model, auto-regression, and an ensemble of exponential smoothing and ARIMA). In terms of median absolute error for forecasting critical care unit occupancy at the two-week horizon, our model only outperformed the naive baseline for 4 out of 14 geographical units and underperformed compared to the ensemble baseline for 5 of them at the 90% confidence level (n = 38). However, for the same level at the 4 week horizon, our model was never statistically outperformed for any unit despite outperforming the baselines 10 times spanning 7 out of 14 geographical units. This implies modest forecasting utility for longer horizons which may justify the application of non-Markovian compartmental models in the context of hospital-strain surveillance for future pandemics.

Oropharyngeal squamous cell carcinoma: p16/p53 immunohistochemistry as a strong predictor of HPV tumour status.

AIMS: Oropharyngeal squamous cell carcinomas (OPSCC) related to human papillomavirus (HPV) infection have a better prognosis than those without HPV infection. Although p16INK4a overexpression is used as a surrogate marker for HPV infection, 5-20% of p16-positive OPSCC are described as being unrelated to HPV infection, with worse overall survival compared to OPSCC-related HPV. There is therefore a risk of undertreating a proportion of OPSCC patients falsely considered to be HPV-driven because of p16 positivity. TP53 mutations are highly prevalent in OPSCC driven by mutagens in tobacco and alcohol. We describe herein a combined p16/p53 algorithm to predict HPV tumour status in OPSCC. METHODS AND RESULTS: A total of 110 OPSCC were identified in the database of the pathology department and were studied using p16 and p53 immunohistochemistry. For p16-positive or p16-negative/wild-type patterns-p53 (WT-p53) cases (n = 63), DNA in-situ hybridisation for high-risk HPV was performed, and if negative the HPV status was controlled by HPV DNA polymerase chain reaction (PCR) (n = 19). A significant association between TP53 mutation and pattern of p53 expression was found (WT-p53, seven of 16, P < 0.001). The p16-positive/WT-p53 was significantly associated with HPV+ tumour status (p16-positive/WT-p53, 50 of 110, P < 0.001). Interestingly, a subset of p16-positive OPSCC was unrelated to HPV (13.5%, eight of 59), and showed mutant-type staining of p53 expression. CONCLUSIONS: The p16 protein immunopositivity in conjunction with the mutant-type pattern of p53 staining helped to reclassify a subset of p16-positive OPSCC as OPSCC-unrelated HPV. This approach could be routinely applied by pathologists involved in the management of OPSCC, because of their potential therapeutic implications.

Spectral photon-counting CT imaging of colorectal peritoneal metastases: initial experience in rats.

Computed tomography imaging plays a major role in the preoperative assessment of tumor burden by providing an accurate mapping of the distribution of peritoneal metastases (PM). Spectral Photon Counting Computed Tomography (SPCCT) is an innovative imaging modality that could overcome the current limitations of conventional CT, offering not only better spatial resolution but also better contrast resolution by allowing the discrimination of multiple contrast agents. Based on this capability, we tested the feasibility of SPCCT in the detection of PM at different time of tumor growth in 16 rats inoculated with CC531 cells using dual-contrast injection protocols in two compartments (i.e. intravenous iodine and intraperitoneal gadolinium or the reverse protocol), compared to surgery. For all peritoneal regions and for both protocols, sensitivity was 69%, specificity was 100% and accuracy was 80%, and the correlation with surgical exploration was strong (p = 0.97; p = 0.0001). No significant difference was found in terms of diagnostic performance, quality of peritoneal opacification or diagnostic quality between the 2 injection protocols. We also showed poor vascularization of peritoneal metastases by measuring low concentrations of contrast agent in the largest lesions using SPCCT, which was confirmed by immunohistochemical analyses. In conclusion, SPCCT using dual-contrast agent injection protocols in 2 compartments is a promising imaging modality to assess the extent of PM in a rat model.

A subset of activated fibroblasts is associated with distant relapse in early luminal breast cancer.

BACKGROUND: Early luminal breast cancer (BC) represents 70% of newly diagnosed BC cases. Among them, small (under 2 cm) BC without lymph node metastasis (classified as T1N0) have been rarely studied, as their prognosis is generally favorable. Nevertheless, up to 5% of luminal T1N0 BC patients relapse with distant metastases that ultimately prove fatal. The aim of our work was to identify the mechanisms involved in metastatic recurrence in these patients. METHODS: Our study addresses the role that autonomous and non-autonomous tumor cell features play with regard to distant recurrence in early luminal BC patients. We created a cohort of T1N0 luminal BC patients (tumors between 0.5-2 cm without lymph node metastasis) with metastatic recurrence ("cases") and corresponding "controls" (without relapse) matched 1:1 on main prognostic factors: age, grade, and proliferation. We deciphered different characteristics of cancer cells and their tumor micro-environment (TME) by deep analyses using immunohistochemistry. We performed in vitro functional assays and highlighted a new mechanism of cooperation between cancer cells and one particular subset of cancer-associated fibroblasts (CAF). RESULTS: We found that specific TME features are indicative of relapse in early luminal BC. Indeed, quantitative histological analyses reveal that "cases" are characterized by significant accumulation of a particular CAF subset (CAF-S1) and decrease in CD4+ T lymphocytes, without any other association with immune cells. In multivariate analysis, TME features, in particular CAF-S1 enrichment, remain significantly associated with recurrence, thereby demonstrating their clinical relevance. Finally, by performing functional analyses, we demonstrated that CAF-S1 pro-metastatic activity is mediated by the CDH11/osteoblast cadherin, consistent with bones being a major site of metastases in luminal BC patients. CONCLUSIONS: This study shows that distant recurrence in T1N0 BC is strongly associated with the presence of CAF-S1 fibroblasts. Moreover, we identify CDH11 as a key player in CAF-S1-mediated pro-metastatic activity. This is independent of tumor cells and represents a new prognostic factor. These results could assist clinicians in identifying luminal BC patients with high risk of relapse. Targeted therapies against CAF-S1 using anti-FAP antibody or CDH11-targeting compounds might help in preventing relapse for such patients with activated stroma.

Aire-dependent genes undergo Clp1-mediated 3'UTR shortening associated with higher transcript stability in the thymus.

The ability of the immune system to avoid autoimmune disease relies on tolerization of thymocytes to self-antigens whose expression and presentation by thymic medullary epithelial cells (mTECs) is controlled predominantly by Aire at the transcriptional level and possibly regulated at other unrecognized levels. Aire-sensitive gene expression is influenced by several molecular factors, some of which belong to the 3'end processing complex, suggesting they might impact transcript stability and levels through an effect on 3'UTR shortening. We discovered that Aire-sensitive genes display a pronounced preference for short-3'UTR transcript isoforms in mTECs, a feature preceding Aire's expression and correlated with the preferential selection of proximal polyA sites by the 3'end processing complex. Through an RNAi screen and generation of a lentigenic mouse, we found that one factor, Clp1, promotes 3'UTR shortening associated with higher transcript stability and expression of Aire-sensitive genes, revealing a post-transcriptional level of control of Aire-activated expression in mTECs.

The bio.tools registry of software tools and data resources for the life sciences.

Bioinformaticians and biologists rely increasingly upon workflows for the flexible utilization of the many life science tools that are needed to optimally convert data into knowledge. We outline a pan-European enterprise to provide a catalogue ( https://bio.tools ) of tools and databases that can be used in these workflows. bio.tools not only lists where to find resources, but also provides a wide variety of practical information.

WellInverter: a web application for the analysis of fluorescent reporter gene data.

BACKGROUND: Fluorescent reporter genes have become widely used for monitoring gene expression in living cells. When a microbial strain carrying a reporter gene is grown in a microplate reader, the fluorescence and the absorbance (optical density) of the culture can be automatically measured every few minutes in a highly parallelized way. The extraction of useful information from the resulting large amounts of data is not easy to achieve, because the fluorescence and absorbance measurements are only indirectly related to promoter activities and protein concentrations, requiring mathematical models of the expression of reporter genes for their interpretation. Although the principles of the analysis of reporter gene data are well-established today, there is a lack of general-purpose bioinformatics tools based on generic measurement models and sound inference procedures. This has motivated the development of WellInverter, a web application based on well-known methods for regularized linear inversion. RESULTS: We present a new version of WellInverter that considerably improves the performance and usability of the original application. In particular, we have put in place a parallel computing architecture with a load balancer to distribute analysis queries over several back-end servers, we have completely redesigned the graphical user interface to better support the different analysis steps, and we have developed a plug-in system for the parsing of data files produced by microplate readers from different manufacturers. We illustrate the functioning of WellInverter by analyzing data of the expression of a fluorescent reporter gene controlled by a phage promoter in growing Escherichia coli populations. We show that the expression pattern in different growth media, supporting different growth rates, corresponds to the pattern expected for a constitutive gene. CONCLUSIONS: The new version of WellInverter is a robust, easy-to-use and scalable web application, which has been deployed on two publicly accessible web servers and which can also be installed locally. A demo version of the application with two sample datasets is available on-line.

Hippocampal Mossy Fibers Synapses in CA3 Pyramidal Cells Are Altered at an Early Stage in a Mouse Model of Alzheimer's Disease.

Early Alzheimer's disease (AD) affects the brain non-uniformly, causing hippocampal memory deficits long before wide-spread brain degeneration becomes evident. Here we addressed whether mossy fiber inputs from the dentate gyrus onto CA3 principal cells are affected in an AD mouse model before amyloid ß plaque deposition. We recorded from CA3 pyramidal cells in a slice preparation from 6-month-old male APP/PS1 mice, and studied synaptic properties and intrinsic excitability. In parallel we performed a morphometric analysis of mossy fiber synapses following viral based labeling and 3D-reconstruction. We found that the basal structural and functional properties as well as presynaptic short-term plasticity at mossy fiber synapses are unaltered at 6 months in APP/PS1 mice. However, transient potentiation of synaptic transmission mediated by activity-dependent release of lipids was abolished. Whereas the presynaptic form of mossy fiber long-term potentiation (LTP) was not affected, the postsynaptic LTP of NMDAR-EPSCs was reduced. In addition, we also report an impairment in feedforward inhibition in CA3 pyramidal cells. This study, together with our previous work describing deficits at CA3-CA3 synapses, provides evidence that early AD affects synapses in a projection-dependent manner at the level of a single neuronal population.SIGNIFICANCE STATEMENT Because loss of episodic memory is considered the cognitive hallmark of Alzheimer's disease (AD), it is important to study whether synaptic circuits involved in the encoding of episodic memory are compromised in AD mouse models. Here we probe alterations in the synaptic connections between the dentate gyrus and CA3, which are thought to be critical for enabling episodic memories to be formed and stored in CA3. We found that forms of synaptic plasticity specific to these synaptic connections are markedly impaired at an early stage in a mouse model of AD, before deposition of ß amyloid plaques. Together with previous work describing deficits at CA3-CA3 synapses, we provide evidence that early AD affects synapses in an input-dependent manner within a single neuronal population.

Early synaptic deficits in the APP/PS1 mouse model of Alzheimer's disease involve neuronal adenosine A2A receptors.

Synaptic plasticity in the autoassociative network of recurrent connections among hippocampal CA3 pyramidal cells is thought to enable the storage of episodic memory. Impaired episodic memory is an early manifestation of cognitive deficits in Alzheimer's disease (AD). In the APP/PS1 mouse model of AD amyloidosis, we show that associative long-term synaptic potentiation (LTP) is abolished in CA3 pyramidal cells at an early stage. This is caused by activation of upregulated neuronal adenosine A2A receptors (A2AR) rather than by dysregulation of NMDAR signalling or altered dendritic spine morphology. Neutralization of A2AR by acute pharmacological inhibition, or downregulation driven by shRNA interference in a single postsynaptic neuron restore associative CA3 LTP. Accordingly, treatment with A2AR antagonists reverts one-trial memory deficits. These results provide mechanistic support to encourage testing the therapeutic efficacy of A2AR antagonists in early AD patients.

Future opportunities and trends for e-infrastructures and life sciences: going beyond the grid to enable life science data analysis.

With the increasingly rapid growth of data in life sciences we are witnessing a major transition in the way research is conducted, from hypothesis-driven studies to data-driven simulations of whole systems. Such approaches necessitate the use of large-scale computational resources and e-infrastructures, such as the European Grid Infrastructure (EGI). EGI, one of key the enablers of the digital European Research Area, is a federation of resource providers set up to deliver sustainable, integrated and secure computing services to European researchers and their international partners. Here we aim to provide the state of the art of Grid/Cloud computing in EU research as viewed from within the field of life sciences, focusing on key infrastructures and projects within the life sciences community. Rather than focusing purely on the technical aspects underlying the currently provided solutions, we outline the design aspects and key characteristics that can be identified across major research approaches. Overall, we aim to provide significant insights into the road ahead by establishing ever-strengthening connections between EGI as a whole and the life sciences community.

RSAT 2015: Regulatory Sequence Analysis Tools.

RSAT (Regulatory Sequence Analysis Tools) is a modular software suite for the analysis of cis-regulatory elements in genome sequences. Its main applications are (i) motif discovery, appropriate to genome-wide data sets like ChIP-seq, (ii) transcription factor binding motif analysis (quality assessment, comparisons and clustering), (iii) comparative genomics and (iv) analysis of regulatory variations. Nine new programs have been added to the 43 described in the 2011 NAR Web Software Issue, including a tool to extract sequences from a list of coordinates (fetch-sequences from UCSC), novel programs dedicated to the analysis of regulatory variants from GWAS or population genomics (retrieve-variation-seq and variation-scan), a program to cluster motifs and visualize the similarities as trees (matrix-clustering). To deal with the drastic increase of sequenced genomes, RSAT public sites have been reorganized into taxon-specific servers. The suite is well-documented with tutorials and published protocols. The software suite is available through Web sites, SOAP/WSDL Web services, virtual machines and stand-alone programs at http://www.rsat.eu/.

BioShaDock: a community driven bioinformatics shared Docker-based tools registry.

Linux container technologies, as represented by Docker, provide an alternative to complex and time-consuming installation processes needed for scientific software. The ease of deployment and the process isolation they enable, as well as the reproducibility they permit across environments and versions, are among the qualities that make them interesting candidates for the construction of bioinformatic infrastructures, at any scale from single workstations to high throughput computing architectures. The Docker Hub is a public registry which can be used to distribute bioinformatic software as Docker images. However, its lack of curation and its genericity make it difficult for a bioinformatics user to find the most appropriate images needed. BioShaDock is a bioinformatics-focused Docker registry, which provides a local and fully controlled environment to build and publish bioinformatic software as portable Docker images. It provides a number of improvements over the base Docker registry on authentication and permissions management, that enable its integration in existing bioinformatic infrastructures such as computing platforms. The metadata associated with the registered images are domain-centric, including for instance concepts defined in the EDAM ontology, a shared and structured vocabulary of commonly used terms in bioinformatics. The registry also includes user defined tags to facilitate its discovery, as well as a link to the tool description in the ELIXIR registry if it already exists. If it does not, the BioShaDock registry will synchronize with the registry to create a new description in the Elixir registry, based on the BioShaDock entry metadata. This link will help users get more information on the tool such as its EDAM operations, input and output types. This allows integration with the ELIXIR Tools and Data Services Registry, thus providing the appropriate visibility of such images to the bioinformatics community.

Membrane lipids tune synaptic transmission by direct modulation of presynaptic potassium channels.

Voltage-gated potassium (Kv) channels are involved in action potential (AP) repolarization in excitable cells. Exogenous application of membrane-derived lipids, such as arachidonic acid (AA), regulates the gating of Kv channels. Whether membrane-derived lipids released under physiological conditions have an impact on neuronal coding through this mechanism is unknown. We show that AA released in an activity-dependent manner from postsynaptic hippocampal CA3 pyramidal cells acts as retrograde messenger, inducing a robust facilitation of mossy fiber (Mf) synaptic transmission over several minutes. AA acts by broadening presynaptic APs through the direct modulation of Kv channels. This form of short-term plasticity can be triggered when postsynaptic cell fires with physiologically relevant patterns and sets the threshold for the induction of the presynaptic form of long-term potentiation (LTP) at hippocampal Mf synapses. Hence, direct modulation of presynaptic Kv channels by activity-dependent release of lipids serves as a physiological mechanism for tuning synaptic transmission.

Selective block of postsynaptic kainate receptors reveals their function at hippocampal mossy fiber synapses.

Progress in understanding the roles of kainate receptors (KARs) in synaptic integration, synaptic networks, and higher brain function has been hampered by the lack of selective pharmacological tools. We have found that UBP310 and related willardiine derivatives, previously characterized as selective GluK1 and GluK3 KAR antagonists, block postsynaptic KARs at hippocampal mossy fiber (MF) CA3 synapses while sparing AMPA and NMDA receptors. We further show that UBP310 is an antagonist of recombinant GluK2/GluK5 receptors, the major population of KARs in the brain. Postsynaptic KAR receptor blockade at MF synapses significantly reduces the sustained depolarization, which builds up during repetitive activity, and impacts on spike transmission mediated by heterosynaptic signals. In addition, KARs present in aberrant MF synapses in the epileptic hippocampus were also blocked by UBP310. Our results support a specific role for postsynaptic KARs in synaptic integration of CA3 pyramidal cells and describe a tool that will be instrumental in understanding the physiopathological role of KARs in the brain.

Deficits in morphofunctional maturation of hippocampal mossy fiber synapses in a mouse model of intellectual disability.

The grik2 gene, coding for the kainate receptor subunit GluK2 (formerly GluR6), is associated with autism spectrum disorders and intellectual disability. Here, we tested the hypothesis that GluK2 could play a role in the appropriate maturation of synaptic circuits involved in learning and memory. We show that both the functional and morphological maturation of hippocampal mossy fiber to CA3 pyramidal cell (mf-CA3) synapses is delayed in mice deficient for the GluK2 subunit (GluK2⁻/⁻). In GluK2⁻/⁻ mice this deficit is manifested by a transient reduction in the amplitude of AMPA-EPSCs at a critical time point of postnatal development, whereas the NMDA component is spared. By combining multiple probability peak fluctuation analysis and immunohistochemistry, we have provided evidence that the decreased amplitude reflects a decrease in the quantal size per mf-CA3 synapse and in the number of active synaptic sites. Furthermore, we analyzed the time course of structural maturation of CA3 synapses by confocal imaging of YFP-expressing cells followed by tridimensional (3D) anatomical reconstruction of thorny excrescences and presynaptic boutons. We show that major changes in synaptic structures occur subsequently to the sharp increase in synaptic transmission, and more importantly that the course of structural maturation of synaptic elements is impaired in GluK2⁻/⁻ mice. This study highlights how a mutation in a gene linked to intellectual disability in the human may lead to a transient reduction of synaptic strength during postnatal development, impacting on the proper formation of neural circuits linked to memory.

Neurexin-neuroligin adhesions capture surface-diffusing AMPA receptors through PSD-95 scaffolds.

The mechanisms governing the recruitment of functional glutamate receptors at nascent excitatory postsynapses following initial axon-dendrite contact remain unclear. We examined here the ability of neurexin/neuroligin adhesions to mobilize AMPA-type glutamate receptors (AMPARs) at postsynapses through a diffusion/trap process involving the scaffold molecule PSD-95. Using single nanoparticle tracking in primary rat and mouse hippocampal neurons overexpressing or lacking neuroligin-1 (Nlg1), a striking inverse correlation was found between AMPAR diffusion and Nlg1 expression level. The use of Nlg1 mutants and inhibitory RNAs against PSD-95 demonstrated that this effect depended on intact Nlg1/PSD-95 interactions. Furthermore, functional AMPARs were recruited within 1 h at nascent Nlg1/PSD-95 clusters assembled by neurexin-1ß multimers, a process requiring AMPAR membrane diffusion. Triggering novel neurexin/neuroligin adhesions also caused a depletion of PSD-95 from native synapses and a drop in AMPAR miniature EPSCs, indicating a competitive mechanism. Finally, both AMPAR level at synapses and AMPAR-dependent synaptic transmission were diminished in hippocampal slices from newborn Nlg1 knock-out mice, confirming an important role of Nlg1 in driving AMPARs to nascent synapses. Together, these data reveal a mechanism by which membrane-diffusing AMPARs can be rapidly trapped at PSD-95 scaffolds assembled at nascent neurexin/neuroligin adhesions, in competition with existing synapses.

NMDA receptor-dependent metaplasticity at hippocampal mossy fiber synapses.

Hippocampal mossy fiber synapses have been reported to lack NMDA receptor (NMDAR)-dependent long-term potentiation (LTP) of AMPA excitatory postsynaptic currents (EPSCs), unlike conventional glutamatergic synapses. An explanation for this difference may reside in the relatively low number of NMDARs at these synapses. Because mossy fiber synapses display LTP selective for NMDARs, we examined whether this would affect the plasticity rules at mossy fiber-CA3 synapses in mouse hippocampal slices. We found that LTP of NMDARs serves as a metaplastic switch making mossy fiber synapses competent for generating NMDAR-dependent LTP of AMPA EPSCs.

Grid computing for improving conformational sampling in NMR structure calculation.

MOTIVATION: Methods for automatic nuclear magnetic resonance (NMR) structure determination need to face a high level of ambiguity encountered in NMR spectra recorded by solid-state NMR and by solution NMR of partially unfolded proteins, leading to time-consuming calculations. The software package Ambiguous Restraints for Iterative Assignment (ARIA) allows for straightforward parallelization of the calculation, as the conformers can be generated in parallel on many nodes. RESULTS: Due to its architecture, the adaptation of ARIA to grid computing can be easily achieved by using the middleware glite and JDL (Job Description Language) scripts. This adaptation makes it possible to address highly ambiguous datasets, because of the much larger conformational sampling that can be generated by use of the grid computational power. AVAILABILITY: The version 2.3.1 of ARIA implemented on the grid is freely available from the ARIA web site: aria.pasteur.fr/downloads.

CURVES+ web server for analyzing and visualizing the helical, backbone and groove parameters of nucleic acid structures.

Curves+, a revised version of the Curves software for analyzing the conformation of nucleic acid structures, is now available as a web server. This version, which can be freely accessed at http://gbio-pbil.ibcp.fr/cgi/Curves_plus/, allows the user to upload a nucleic acid structure file, choose the nucleotides to be analyzed and after optionally setting a number of input variables, view the numerical and graphic results online or download files containing a set of helical, backbone and groove parameters that fully describe the structure. PDB format files are also provided for offline visualization of the helical axis and groove geometry.

BioXSD: the common data-exchange format for everyday bioinformatics web services.

MOTIVATION: The world-wide community of life scientists has access to a large number of public bioinformatics databases and tools, which are developed and deployed using diverse technologies and designs. More and more of the resources offer programmatic web-service interface. However, efficient use of the resources is hampered by the lack of widely used, standard data-exchange formats for the basic, everyday bioinformatics data types. RESULTS: BioXSD has been developed as a candidate for standard, canonical exchange format for basic bioinformatics data. BioXSD is represented by a dedicated XML Schema and defines syntax for biological sequences, sequence annotations, alignments and references to resources. We have adapted a set of web services to use BioXSD as the input and output format, and implemented a test-case workflow. This demonstrates that the approach is feasible and provides smooth interoperability. Semantics for BioXSD is provided by annotation with the EDAM ontology. We discuss in a separate section how BioXSD relates to other initiatives and approaches, including existing standards and the Semantic Web. AVAILABILITY: The BioXSD 1.0 XML Schema is freely available at http://www.bioxsd.org/BioXSD-1.0.xsd under the Creative Commons BY-ND 3.0 license. The http://bioxsd.org web page offers documentation, examples of data in BioXSD format, example workflows with source codes in common programming languages, an updated list of compatible web services and tools and a repository of feature requests from the community.