PI3Kγ in B cells promotes antibody responses and generation of antibody-secreting cells.

The differentiation of naive and memory B cells into antibody-secreting cells (ASCs) is a key feature of adaptive immunity. The requirement for phosphoinositide 3-kinase-delta (PI3Kδ) to support B cell biology has been investigated intensively; however, specific functions of the related phosphoinositide 3-kinase-gamma (PI3Kγ) complex in B lineage cells have not. In the present study, we report that PI3Kγ promotes robust antibody responses induced by T cell-dependent antigens. The inborn error of immunity caused by human deficiency in PI3Kγ results in broad humoral defects, prompting our investigation of roles for this kinase in antibody responses. Using mouse immunization models, we found that PI3Kγ functions cell intrinsically within activated B cells in a kinase activity-dependent manner to transduce signals required for the transcriptional program supporting differentiation of ASCs. Furthermore, ASC fate choice coincides with upregulation of PIK3CG expression and is impaired in the context of PI3Kγ disruption in naive B cells on in vitro CD40-/cytokine-driven activation, in memory B cells on toll-like receptor activation, or in human tonsillar organoids. Taken together, our study uncovers a fundamental role for PI3Kγ in supporting humoral immunity by integrating signals instructing commitment to the ASC fate.

mlf-core: a framework for deterministic machine learning.

MOTIVATION: Machine learning has shown extensive growth in recent years and is now routinely applied to sensitive areas. To allow appropriate verification of predictive models before deployment, models must be deterministic. Solely fixing all random seeds is not sufficient for deterministic machine learning, as major machine learning libraries default to the usage of nondeterministic algorithms based on atomic operations. RESULTS: Various machine learning libraries released deterministic counterparts to the nondeterministic algorithms. We evaluated the effect of these algorithms on determinism and runtime. Based on these results, we formulated a set of requirements for deterministic machine learning and developed a new software solution, the mlf-core ecosystem, which aids machine learning projects to meet and keep these requirements. We applied mlf-core to develop deterministic models in various biomedical fields including a single-cell autoencoder with TensorFlow, a PyTorch-based U-Net model for liver-tumor segmentation in computed tomography scans, and a liver cancer classifier based on gene expression profiles with XGBoost. AVAILABILITY AND IMPLEMENTATION: The complete data together with the implementations of the mlf-core ecosystem and use case models are available at https://github.com/mlf-core.

A data management infrastructure for the integration of imaging and omics data in life sciences.

BACKGROUND: As technical developments in omics and biomedical imaging increase the throughput of data generation in life sciences, the need for information systems capable of managing heterogeneous digital assets is increasing. In particular, systems supporting the findability, accessibility, interoperability, and reusability (FAIR) principles of scientific data management. RESULTS: We propose a Service Oriented Architecture approach for integrated management and analysis of multi-omics and biomedical imaging data. Our architecture introduces an image management system into a FAIR-supporting, web-based platform for omics data management. Interoperable metadata models and middleware components implement the required data management operations. The resulting architecture allows for FAIR management of omics and imaging data, facilitating metadata queries from software applications. The applicability of the proposed architecture is demonstrated using two technical proofs of concept and a use case, aimed at molecular plant biology and clinical liver cancer research, which integrate various imaging and omics modalities. CONCLUSIONS: We describe a data management architecture for integrated, FAIR-supporting management of omics and biomedical imaging data, and exemplify its applicability for basic biology research and clinical studies. We anticipate that FAIR data management systems for multi-modal data repositories will play a pivotal role in data-driven research, including studies which leverage advanced machine learning methods, as the joint analysis of omics and imaging data, in conjunction with phenotypic metadata, becomes not only desirable but necessary to derive novel insights into biological processes.

Downregulation of TGR5 (GPBAR1) in biliary epithelial cells contributes to the pathogenesis of sclerosing cholangitis.

BACKGROUND & AIMS: Primary sclerosing cholangitis (PSC) is characterized by chronic inflammation and progressive fibrosis of the biliary tree. The bile acid receptor TGR5 (GPBAR1) is found on biliary epithelial cells (BECs), where it promotes secretion, proliferation and tight junction integrity. Thus, we speculated that changes in TGR5-expression in BECs may contribute to PSC pathogenesis. METHODS: TGR5-expression and -localization were analyzed in PSC livers and liver tissue, isolated bile ducts and BECs from Abcb4-/-, Abcb4-/-/Tgr5Tg and ursodeoxycholic acid (UDCA)- or 24-norursodeoxycholic acid (norUDCA)-fed Abcb4-/- mice. The effects of IL8/IL8 homologues on TGR5 mRNA and protein levels were studied. BEC gene expression was analyzed by single-cell transcriptomics (scRNA-seq) from distinct mouse models. RESULTS: TGR5 mRNA expression and immunofluorescence staining intensity were reduced in BECs of PSC and Abcb4-/- livers, in Abcb4-/- extrahepatic bile ducts, but not in intrahepatic macrophages. No changes in TGR5 BEC fluorescence intensity were detected in liver tissue of other liver diseases, including primary biliary cholangitis. Incubation of BECs with IL8/IL8 homologues, but not with other cytokines, reduced TGR5 mRNA and protein levels. BECs from Abcb4-/- mice had lower levels of phosphorylated Erk and higher expression levels of Icam1, Vcam1 and Tgfß2. Overexpression of Tgr5 abolished the activated inflammatory phenotype characteristic of Abcb4-/- BECs. NorUDCA-feeding restored TGR5-expression levels in BECs in Abcb4-/- livers. CONCLUSIONS: Reduced TGR5 levels in BECs from patients with PSC and Abcb4-/- mice promote development of a reactive BEC phenotype, aggravate biliary injury and thus contribute to the pathogenesis of sclerosing cholangitis. Restoration of biliary TGR5-expression levels represents a previously unknown mechanism of action of norUDCA. LAY SUMMARY: Primary sclerosing cholangitis (PSC) is a chronic cholestatic liver disease-associated with progressive inflammation of the bile duct, leading to fibrosis and end-stage liver disease. Bile acid (BA) toxicity may contribute to the development and disease progression of PSC. TGR5 is a membrane-bound receptor for BAs, which is found on bile ducts and protects bile ducts from BA toxicity. In this study, we show that TGR5 levels were reduced in bile ducts from PSC livers and in bile ducts from a genetic mouse model of PSC. Our investigations indicate that lower levels of TGR5 in bile ducts may contribute to PSC development and progression. Furthermore, treatment with norUDCA, a drug currently being tested in a phase III trial for PSC, restored TGR5 levels in biliary epithelial cells.

Morphing of Amphipathic Helices to Explore the Activity and Selectivity of Membranolytic Antimicrobial Peptides.

Naturally occurring membranolytic antimicrobial peptides (AMPs) are rarely cell-type selective and highly potent at the same time. Template-based peptide design can be used to generate AMPs with improved properties de novo. Following this approach, 18 linear peptides were obtained by computationally morphing the natural AMP Aurein 2.2d2 GLFDIVKKVVGALG into the synthetic model AMP KLLKLLKKLLKLLK. Eleven of the 18 chimeric designs inhibited the growth of Staphylococcus aureus, and six peptides were tested and found to be active against one resistant pathogenic strain or more. One of the peptides was broadly active against bacterial and fungal pathogens without exhibiting toxicity to certain human cell lines. Solution nuclear magnetic resonance and molecular dynamics simulation suggested an oblique-oriented membrane insertion mechanism of this helical de novo peptide. Temperature-resolved circular dichroism spectroscopy pointed to conformational flexibility as an essential feature of cell-type selective AMPs.

Simulated Molecular Evolution for Anticancer Peptide Design.

A computational technique based on a simulated molecular evolution protocol was employed for anticancer peptide (ACP) design. Starting from known ACPs, innovative bioactive peptides were automatically generated in computer-assisted design-synthesize-test cycles. This design algorithm offers a viable strategy for the generation of novel peptide sequences, without requiring a priori structure-activity knowledge. Sequence morphing and activity improvement were achieved through iterative amino acid variation and selection. Results show that not only the interaction of ACPs with the target membrane is important for their anticancer activity, but also the degree of peptide dimerization, which was corroborated by temperature profiling and electrospray mass spectrometry.

modlAMP: Python for antimicrobial peptides.

SUMMARY: We have implemented the lecular esign aboratory's nti icrobial eptides package ( ), a Python-based software package for the design, classification and visual representation of peptide data. modlAMP offers functions for molecular descriptor calculation and the retrieval of amino acid sequences from public or local sequence databases, and provides instant access to precompiled datasets for machine learning. The package also contains methods for the analysis and representation of circular dichroism spectra. AVAILABILITY AND IMPLEMENTATION: The modlAMP Python package is available under the BSD license from URL http://doi.org/10.5905/ethz-1007-72 or via pip from the Python Package Index (PyPI). CONTACT: gisbert.schneider@pharma.ethz.ch. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

De Novo Fragment Design for Drug Discovery and Chemical Biology.

Automated molecular de novo design led to the discovery of an innovative inhibitor of death-associated protein kinase 3 (DAPK3). An unprecedented crystal structure of the inactive DAPK3 homodimer shows the fragment-like hit bound to the ATP pocket. Target prediction software based on machine learning models correctly identified additional macromolecular targets of the computationally designed compound and the structurally related marketed drug azosemide. The study validates computational de novo design as a prime method for generating chemical probes and starting points for drug discovery.

How tool combinations in different pipeline versions affect the outcome in RNA-seq analysis.

Data analysis tools are continuously changed and improved over time. In order to test how these changes influence the comparability between analyses, the output of different workflow options of the nf-core/rnaseq pipeline were compared. Five different pipeline settings (STAR+Salmon, STAR+RSEM, STAR+featureCounts, HISAT2+featureCounts, pseudoaligner Salmon) were run on three datasets (human, Arabidopsis, zebrafish) containing spike-ins of the External RNA Control Consortium (ERCC). Fold change ratios and differential expression of genes and spike-ins were used for comparative analyses of the different tools and versions settings of the pipeline. An overlap of 85% for differential gene classification between pipelines could be shown. Genes interpreted with a bias were mostly those present at lower concentration. Also, the number of isoforms and exons per gene were determinants. Previous pipeline versions using featureCounts showed a higher sensitivity to detect one-isoform genes like ERCC. To ensure data comparability in long-term analysis series it would be recommendable to either stay with the pipeline version the series was initialized with or to run both versions during a transition time in order to ensure that the target genes are addressed the same way.

nf-core/airrflow: an adaptive immune receptor repertoire analysis workflow employing the Immcantation framework.

Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) is a valuable experimental tool to study the immune state in health and following immune challenges such as infectious diseases, (auto)immune diseases, and cancer. Several tools have been developed to reconstruct B cell and T cell receptor sequences from AIRR-seq data and infer B and T cell clonal relationships. However, currently available tools offer limited parallelization across samples, scalability or portability to high-performance computing infrastructures. To address this need, we developed nf-core/airrflow, an end-to-end bulk and single-cell AIRR-seq processing workflow which integrates the Immcantation Framework following BCR and TCR sequencing data analysis best practices. The Immcantation Framework is a comprehensive toolset, which allows the processing of bulk and single-cell AIRR-seq data from raw read processing to clonal inference. nf-core/airrflow is written in Nextflow and is part of the nf-core project, which collects community contributed and curated Nextflow workflows for a wide variety of analysis tasks. We assessed the performance of nf-core/airrflow on simulated sequencing data with sequencing errors and show example results with real datasets. To demonstrate the applicability of nf-core/airrflow to the high-throughput processing of large AIRR-seq datasets, we validated and extended previously reported findings of convergent antibody responses to SARS-CoV-2 by analyzing 97 COVID-19 infected individuals and 99 healthy controls, including a mixture of bulk and single-cell sequencing datasets. Using this dataset, we extended the convergence findings to 20 additional subjects, highlighting the applicability of nf-core/airrflow to validate findings in small in-house cohorts with reanalysis of large publicly available AIRR datasets. nf-core/airrflow is available free of charge, under the MIT license on GitHub (https://github.com/nf-core/airrflow). Detailed documentation and example results are available on the nf-core website at (https://nf-co.re/airrflow).

Scalable and efficient DNA sequencing analysis on different compute infrastructures aiding variant discovery.

DNA variation analysis has become indispensable in many aspects of modern biomedicine, most prominently in the comparison of normal and tumor samples. Thousands of samples are collected in local sequencing efforts and public databases requiring highly scalable, portable, and automated workflows for streamlined processing. Here, we present nf-core/sarek 3, a well-established, comprehensive variant calling and annotation pipeline for germline and somatic samples. It is suitable for any genome with a known reference. We present a full rewrite of the original pipeline showing a significant reduction of storage requirements by using the CRAM format and runtime by increasing intra-sample parallelization. Both are leading to a 70% cost reduction in commercial clouds enabling users to do large-scale and cross-platform data analysis while keeping costs and CO2 emissions low. The code is available at https://nf-co.re/sarek.

Apoptotic DNA degradation into oligonucleosomal fragments, but not apoptotic nuclear morphology, relies on a cytosolic pool of DFF40/CAD endonuclease.

Apoptotic cell death is characterized by nuclear fragmentation and oligonucleosomal DNA degradation, mediated by the caspase-dependent specific activation of DFF40/CAD endonuclease. Here, we describe how, upon apoptotic stimuli, SK-N-AS human neuroblastoma-derived cells show apoptotic nuclear morphology without displaying concomitant internucleosomal DNA fragmentation. Cytotoxicity afforded after staurosporine treatment is comparable with that obtained in SH-SY5Y cells, which exhibit a complete apoptotic phenotype. SK-N-AS cell death is a caspase-dependent process that can be impaired by the pan-caspase inhibitor q-VD-OPh. The endogenous inhibitor of DFF40/CAD, ICAD, is correctly processed, and dff40/cad cDNA sequence does not reveal mutations altering its amino acid composition. Biochemical approaches show that both SH-SY5Y and SK-N-AS resting cells express comparable levels of DFF40/CAD. However, the endonuclease is poorly expressed in the cytosolic fraction of healthy SK-N-AS cells. Despite this differential subcellular distribution of DFF40/CAD, we find no differences in the subcellular localization of both pro-caspase-3 and ICAD between the analyzed cell lines. After staurosporine treatment, the preferential processing of ICAD in the cytosolic fraction allows the translocation of DFF40/CAD from this fraction to a chromatin-enriched one. Therefore, the low levels of cytosolic DFF40/CAD detected in SK-N-AS cells determine the absence of DNA laddering after staurosporine treatment. In these cells DFF40/CAD cytosolic levels can be restored by the overexpression of their own endonuclease, which is sufficient to make them proficient at degrading their chromatin into oligonucleosome-size fragments after staurosporine treatment. Altogether, the cytosolic levels of DFF40/CAD are determinants in achieving a complete apoptotic phenotype, including oligonucleosomal DNA degradation.

Cladribine treatment specifically affects peripheral blood memory B cell clones and clonal expansion in multiple sclerosis patients.

Introduction: B cells are acknowledged as crucial players in the pathogenesis of multiple sclerosis (MS). Several disease modifying drugs including cladribine have been shown to exert differential effects on peripheral blood B cell subsets. However, little is known regarding functional changes within the peripheral B cell populations. In this study, we obtained a detailed picture of B cell repertoire changes under cladribine treatment on a combined immunoglobulin (Ig) transcriptome and proteome level. Methods: We performed next-generation sequencing of Ig heavy chain (IGH) transcripts and Ig mass spectrometry in cladribine-treated patients with relapsing-remitting multiple sclerosis (n = 8) at baseline and after 6 and 12 months of treatment in order to generate Ig transcriptome and Ig peptide libraries. Ig peptides were overlapped with the corresponding IGH transcriptome in order to analyze B cell clones on a combined transcriptome and proteome level. Results: The analysis of peripheral blood B cell percentages pointed towards a significant decrease of memory B cells and an increase of naive B cells following cladribine therapy. While basic IGH repertoire parameters (e.g. variable heavy chain family usage and Ig subclasses) were only slightly affected by cladribine treatment, a significantly decreased number of clones and significantly lower diversity in the memory subset was noticeable at 6 months following treatment which was sustained at 12 months. When looking at B-cell clones comprising sequences from the different time-points, clones spanning between all three time-points were significantly more frequent than clones including sequences from two time-points. Furthermore, Ig proteome analyses showed that Ig transcriptome specific peptides could mostly be equally aligned to all three time-points pointing towards a proportion of B-cell clones that are maintained during treatment. Discussion: Our findings suggest that peripheral B cell related treatment effects of cladribine tablets might be exerted through a reduction of possibly disease relevant clones in the memory B cell subset without disrupting the overall clonal composition of B cells. Our results -at least partially- might explain the relatively mild side effects regarding infections and the sustained immune response after vaccinations during treatment. However, exact disease driving B cell subsets and their effects remain unknown and should be addressed in future studies.

Clinical outcome of biomarker-guided therapies in adult patients with tumors of the nervous system.

Background: The clinical utility of molecular profiling and targeted therapies for neuro-oncology patients outside of clinical trials is not established. We aimed at investigating feasibility and clinical utility of molecular profiling and targeted therapy in adult patients with advanced tumors in the nervous system within a prospective observational study. Methods: molecular tumor board (MTB)@ZPM (NCT03503149) is a prospective observational precision medicine study for patients with advanced tumors. After inclusion of patients, we performed comprehensive molecular profiling, formulated ranked biomarker-guided therapy recommendations based on consensus by the MTB, and collected prospective clinical outcome data. Results: Here, we present initial data of 661 adult patients with tumors of the nervous system enrolled by December 31, 2021. Of these, 408 patients were presented at the MTB. Molecular-instructed therapy recommendations could be made in 380/408 (93.1%) cases and were prioritized by evidence levels. Therapies were initiated in 86/380 (22.6%) cases until data cutoff. We observed a progression-free survival ratio >1.3 in 31.3% of patients. Conclusions: Our study supports the clinical utility of biomarker-guided therapies for neuro-oncology patients and indicates clinical benefit in a subset of patients. Our data might inform future clinical trials, translational studies, and even clinical care.

nf-core/mag: a best-practice pipeline for metagenome hybrid assembly and binning.

The analysis of shotgun metagenomic data provides valuable insights into microbial communities, while allowing resolution at individual genome level. In absence of complete reference genomes, this requires the reconstruction of metagenome assembled genomes (MAGs) from sequencing reads. We present the nf-core/mag pipeline for metagenome assembly, binning and taxonomic classification. It can optionally combine short and long reads to increase assembly continuity and utilize sample-wise group-information for co-assembly and genome binning. The pipeline is easy to install-all dependencies are provided within containers-portable and reproducible. It is written in Nextflow and developed as part of the nf-core initiative for best-practice pipeline development. All codes are hosted on GitHub under the nf-core organization https://github.com/nf-core/mag and released under the MIT license.

Next Generation Sequencing of Cerebrospinal Fluid B Cell Repertoires in Multiple Sclerosis and Other Neuro-Inflammatory Diseases-A Comprehensive Review.

During the last few decades, the role of B cells has been well established and redefined in neuro-inflammatory diseases, including multiple sclerosis and autoantibody-associated diseases. In particular, B cell maturation and trafficking across the blood-brain barrier (BBB) has recently been deciphered with the development of next-generation sequencing (NGS) approaches, which allow the assessment of representative cerebrospinal fluid (CSF) and peripheral blood B cell repertoires. In this review, we perform literature research focusing on NGS studies that allow further insights into B cell pathophysiology during neuro-inflammation. Besides the analysis of CSF B cells, the paralleled assessment of peripheral blood B cell repertoire provides deep insights into not only the CSF compartment, but also in B cell trafficking patterns across the BBB. In multiple sclerosis, CSF-specific B cell maturation, in combination with a bidirectional exchange of B cells across the BBB, is consistently detectable. These data suggest that B cells most likely encounter antigen(s) within the CSF and migrate across the BBB, with further maturation also taking place in the periphery. Autoantibody-mediated diseases, such as neuromyelitis optica spectrum disorder and LGI1 / NMDAR encephalitis, also show features of a CSF-specific B cell maturation and clonal connectivity with peripheral blood. In conclusion, these data suggest an intense exchange of B cells across the BBB, possibly feeding autoimmune circuits. Further developments in sequencing technologies will help to dissect the exact pathophysiologic mechanisms of B cells during neuro-inflammation.

Specific Induction of Double Negative B Cells During Protective and Pathogenic Immune Responses.

Double negative (DN) (CD19+CD20lowCD27-IgD-) B cells are expanded in patients with autoimmune and infectious diseases; however their role in the humoral immune response remains unclear. Using systematic flow cytometric analyses of peripheral blood B cell subsets, we observed an inflated DN B cell population in patients with variety of active inflammatory conditions: myasthenia gravis, Guillain-Barré syndrome, neuromyelitis optica spectrum disorder, meningitis/encephalitis, and rheumatic disorders. Furthermore, we were able to induce DN B cells in healthy subjects following vaccination against influenza and tick borne encephalitis virus. Transcriptome analysis revealed a gene expression profile in DN B cells that clustered with naïve B cells, memory B cells, and plasmablasts. Immunoglobulin VH transcriptome sequencing and analysis of recombinant antibodies revealed clonal expansion of DN B cells that were targeted against the vaccine antigen. Our study suggests that DN B cells are expanded in multiple inflammatory neurologic diseases and represent an inducible B cell population that responds to antigenic stimulation, possibly through an extra-follicular maturation pathway.

Clinical and Genetic Tumor Characteristics of Responding and Non-Responding Patients to PD-1 Inhibition in Hepatocellular Carcinoma.

Immune checkpoint inhibitors (ICIs) belong to the therapeutic armamentarium in advanced hepatocellular carcinoma (HCC). However, only a minority of patients benefit from immunotherapy. Therefore, we aimed to identify indicators of therapy response. This multicenter analysis included 99 HCC patients. Progression-free (PFS) and overall survival (OS) were studied by Kaplan-Meier analyses for clinical parameters using weighted log-rank testing. Next-generation sequencing (NGS) was performed in a subset of 15 patients. The objective response (OR) rate was 19% median OS (mOS)16.7 months. Forty-one percent reached a PFS > 6 months; these patients had a significantly longer mOS (32.0 vs. 8.5 months). Child-Pugh (CP) A and B patients showed a mOS of 22.1 and 12.1 months, respectively. Ten of thirty CP-B patients reached PFS > 6 months, including 3 patients with an OR. Tumor mutational burden (TMB) could not predict responders. Of note, antibiotic treatment within 30 days around ICI initiation was associated with significantly shorter mOS (8.5 vs. 17.4 months). Taken together, this study shows favorable outcomes for OS with low AFP, OR, and PFS > 6 months. No specific genetic pattern, including TMB, could identify responders. Antibiotics around treatment initiation were associated with worse outcome, suggesting an influence of the host microbiome on therapy success.

In silico design and optimization of selective membranolytic anticancer peptides.

Membranolytic anticancer peptides represent a potential strategy in the fight against cancer. However, our understanding of the underlying structure-activity relationships and the mechanisms driving their cell selectivity is still limited. We developed a computational approach as a step towards the rational design of potent and selective anticancer peptides. This machine learning model distinguishes between peptides with and without anticancer activity. This classifier was experimentally validated by synthesizing and testing a selection of 12 computationally generated peptides. In total, 83% of these predictions were correct. We then utilized an evolutionary molecular design algorithm to improve the peptide selectivity for cancer cells. This simulated molecular evolution process led to a five-fold selectivity increase with regard to human dermal microvascular endothelial cells and more than ten-fold improvement towards human erythrocytes. The results of the present study advocate for the applicability of machine learning models and evolutionary algorithms to design and optimize novel synthetic anticancer peptides with reduced hemolytic liability and increased cell-type selectivity.

De novo design of anticancer peptides by ensemble artificial neural networks.

Membranolytic anticancer peptides (ACPs) are drawing increasing attention as potential future therapeutics against cancer, due to their ability to hinder the development of cellular resistance and their potential to overcome common hurdles of chemotherapy, e.g., side effects and cytotoxicity. In this work, we present an ensemble machine learning model to design potent ACPs. Four counter-propagation artificial neural-networks were trained to identify peptides that kill breast and/or lung cancer cells. For prospective application of the ensemble model, we selected 14 peptides from a total of 1000 de novo designs, for synthesis and testing in vitro on breast cancer (MCF7) and lung cancer (A549) cell lines. Six de novo designs showed anticancer activity in vitro, five of which against both MCF7 and A549 cell lines. The novel active peptides populate uncharted regions of ACP sequence space.