Search

Background. Transcriptome analysis could be an additional diagnostic parameter in diagnosing kidney transplant (KTx) rejection. Here, we assessed feasibility and potential of NanoString nCounter analysis of KTx biopsies to aid the classification of rejection in clinical practice using both the Banff-Human Organ Transplant (B-HOT) panel and a customized antibody-mediated rejection (AMR)-specific NanoString nCounter Elements (Elements) panel. Additionally, we explored the potential for the classification of KTx rejection building and testing a classifier within our dataset. Methods. Ninety-six formalin-fixed paraffin-embedded KTx biopsies were retrieved from the archives of the ErasmusMC Rotterdam and the University Hospital Cologne. Biopsies with AMR, borderline or T cell-mediated rejections (BLorTCMR), and no rejection were compared using the B-HOT and Elements panels. Results. High correlation between gene expression levels was found when comparing the 2 chemistries pairwise (r = 0.76-0.88). Differential gene expression (false discovery rate; P < 0.05) was identified in biopsies diagnosed with AMR (B-HOT: 294; Elements: 76) and BLorTCMR (B-HOT: 353; Elements: 57) compared with no rejection. Using the most predictive genes from the B-HOT analysis and the Element analysis, 2 least absolute shrinkage and selection operators-based regression models to classify biopsies as AMR versus no AMR (BLorTCMR or no rejection) were developed achieving an receiver-operating-characteristic curve of 0.994 and 0.894, sensitivity of 0.821 and 0.480, and specificity of 1.00 and 0.979, respectively, during cross-validation. Conclusions. Transcriptomic analysis is feasible on KTx biopsies previously used for diagnostic purposes. The B-HOT panel has the potential to differentiate AMR from BLorTCMR or no rejection and could prove valuable in aiding kidney transplant rejection classification.

Supplementary Figure 7 from Intrinsic Gene Expression Profiles of Gliomas Are a Better Predictor of Survival than Histology

Supplementary Figure 5 from Intrinsic Gene Expression Profiles of Gliomas Are a Better Predictor of Survival than Histology

Supplementary Figure 10 from Intrinsic Gene Expression Profiles of Gliomas Are a Better Predictor of Survival than Histology

Supplementary Figure 6 from Intrinsic Gene Expression Profiles of Gliomas Are a Better Predictor of Survival than Histology

Supplementary Figure 8 from Intrinsic Gene Expression Profiles of Gliomas Are a Better Predictor of Survival than Histology

Supplementary Figure 3 from Intrinsic Gene Expression Profiles of Gliomas Are a Better Predictor of Survival than Histology

Gliomas are the most common primary brain tumors with heterogeneous morphology and variable prognosis. Treatment decisions in patients rely mainly on histologic classification and clinical parameters. However, differences between histologic subclasses and grades are subtle, and classifying gliomas is subject to a large interobserver variability. To improve current classification standards, we have performed gene expression profiling on a large cohort of glioma samples of all histologic subtypes and grades. We identified seven distinct molecular subgroups that correlate with survival. These include two favorable prognostic subgroups (median survival, >4.7 years), two with intermediate prognosis (median survival, 1–4 years), two with poor prognosis (median survival

Supplementary Figure 2 from Intrinsic Gene Expression Profiles of Gliomas Are a Better Predictor of Survival than Histology

Supplementary Figure 4 from Intrinsic Gene Expression Profiles of Gliomas Are a Better Predictor of Survival than Histology

Supplementary Tables 1-6 from Intrinsic Gene Expression Profiles of Gliomas Are a Better Predictor of Survival than Histology

Supplementary Figure 9 from Intrinsic Gene Expression Profiles of Gliomas Are a Better Predictor of Survival than Histology

Supplementary Methods from Intrinsic Gene Expression Profiles of Gliomas Are a Better Predictor of Survival than Histology

BACKGROUND Major obstacles that have impeded the development of effective new therapies for GBM include inter- and intratumoral heterogeneity, the blood-brain-barrier and use of sub-optimal cell line-based preclinical models. Taking these hurdles into account, we have set up a patient-derived GBM drug-screening platform. We optimized protocols to improve cell culture success rate and retrospectively assessed the predictive power of our assay for patient response to TMZ. A large panel of GBM cells was screened for sensitivity to available oncological agents. Drugs of interest were selected based on favorable physicochemical properties for BBB crossing and potent activity in (a subset of) GBM cultures. Finally, we determined the success rate of performing a small-scale screen with 20 selected agents within 4 weeks of receiving tumor tissue. RESULTS By combining both CUSA and tissue piece-derived dissociation protocols, culture success increased to 95%, ensuring representation of the near-complete spectrum of GBM subtypes. Single-cell sequencing studies confirmed heterogeneity in our low-passage cell cultures. In vitro screening of TMZ on a large cohort (n = 55) identified 3 response categories (responders/intermediates/non-responders) for which Cox regression analysis revealed significantly different overall survival curves of corresponding patients. Screening of 107 FDA-approved anticancer agents on 45 GBM cultures underscored the tremendous intertumoral heterogeneity in drug sensitivities. We identified 20 potent agents each effective at clinically-achievable concentrations in (a subset of) GBM cultures and having favorable BBB penetration properties (CNS-MPO score). Screening of these agents on a per patient basis within 4 weeks of receiving tissue was successful in 18 out of 24 (75%) tested tumors. In the remaining cases the tumor cells grew very slowly and longer culture times were required. CONCLUSION Our drug screening platform offers a tool to predict TMZ response and assess sensitivity to candidate treatments, either for GBM subsets or on a per patient basis.

IntroductionA decentralized and multi-platform-compatible molecular diagnostic tool for kidney transplant biopsies could improve the dissemination and exploitation of this technology, increasing its clinical impact. As a first step towards this molecular diagnostic tool, we developed and validated a classifier using the genes of the Banff-Human Organ Transplant (B-HOT) panel extracted from a historical Molecular Microscope® Diagnostic system microarray dataset. Furthermore, we evaluated the discriminative power of the B-HOT panel in a clinical scenario.Materials and MethodsGene expression data from 1,181 kidney transplant biopsies were used as training data for three random forest models to predict kidney transplant biopsy Banff categories, including non-rejection (NR), antibody-mediated rejection (ABMR), and T-cell-mediated rejection (TCMR). Performance was evaluated using nested cross-validation. The three models used different sets of input features: the first model (B-HOT Model) was trained on only the genes included in the B-HOT panel, the second model (Feature Selection Model) was based on sequential forward feature selection from all available genes, and the third model (B-HOT+ Model) was based on the combination of the two models, i.e. B-HOT panel genes plus highly predictive genes from the sequential forward feature selection. After performance assessment on cross-validation, the best-performing model was validated on an external independent dataset based on a different microarray version.ResultsThe best performances were achieved by the B-HOT+ Model, a multilabel random forest model trained on B-HOT panel genes with the addition of the 6 most predictive genes of the Feature Selection Model (ST7, KLRC4-KLRK1, TRBC1, TRBV6-5, TRBV19, and ZFX), with a mean accuracy of 92.1% during cross-validation. On the validation set, the same model achieved Area Under the ROC Curve (AUC) of 0.965 and 0.982 for NR and ABMR respectively.DiscussionThis kidney transplant biopsy classifier is one step closer to the development of a decentralized kidney transplant biopsy classifier that is effective on data derived from different gene expression platforms. The B-HOT panel proved to be a reliable highly-predictive panel for kidney transplant rejection classification. Furthermore, we propose to include the aforementioned 6 genes in the B-HOT panel for further optimization of this commercially available panel.

The prevalence of diabetic metabolic derangement (DMetD) has increased dramatically over the last decades. Although there is increasing evidence that DMetD is associated with cardiac dysfunction, the early DMetD-induced myocardial alterations remain incompletely understood. Here, we studied early DMetD-related cardiac changes in a clinically relevant large animal model. DMetD was established in adult male Göttingen miniswine by streptozotocin injections and a high-fat, high-sugar diet, while control animals remained on normal pig chow. Five months later left ventricular (LV) function was assessed by echocardiography and hemodynamic measurements, followed by comprehensive biochemical, molecular and histological analyses. Robust DMetD developed, evidenced by hyperglycemia, hypercholesterolemia and hypertriglyceridemia. DMetD resulted in altered LV nitroso-redox balance, increased superoxide production—principally due to endothelial nitric oxide synthase (eNOS) uncoupling—reduced nitric oxide (NO) production, alterations in myocardial gene-expression—particularly genes related to glucose and fatty acid metabolism—and mitochondrial dysfunction. These abnormalities were accompanied by increased passive force of isolated cardiomyocytes, and impaired LV diastolic function, evidenced by reduced LV peak untwist velocity and increased E/e′. However, LV weight, volume, collagen content, and cardiomyocyte cross-sectional area were unchanged at this stage of DMetD. In conclusion, DMetD, in a clinically relevant large-animal model results in myocardial oxidative stress, eNOS uncoupling and reduced NO production, together with an altered metabolic gene expression profile and mitochondrial dysfunction. These molecular alterations are associated with stiffening of the cardiomyocytes and early diastolic dysfunction before any structural cardiac remodeling occurs. Therapies should be directed to ameliorate these early DMetD-induced myocardial changes to prevent the development of overt cardiac failure.

Background Fusion genes are typically identified by RNA sequencing (RNA-seq) without elucidating the causal genomic breakpoints. However, non–poly(A)-enriched RNA-seq contains large proportions of intronic reads that also span genomic breakpoints. Results We have developed an algorithm, Dr. Disco, that searches for fusion transcripts by taking an entire reference genome into account as search space. This includes exons but also introns, intergenic regions, and sequences that do not meet splice junction motifs. Using 1,275 RNA-seq samples, we investigated to what extent genomic breakpoints can be extracted from RNA-seq data and their implications regarding poly(A)-enriched and ribosomal RNA–minus RNA-seq data. Comparison with whole-genome sequencing data revealed that most genomic breakpoints are not, or minimally, transcribed while, in contrast, the genomic breakpoints of all 32 TMPRSS2-ERG–positive tumours were present at RNA level. We also revealed tumours in which the ERG breakpoint was located before ERG, which co-existed with additional deletions and messenger RNA that incorporated intergenic cryptic exons. In breast cancer we identified rearrangement hot spots near CCND1 and in glioma near CDK4 and MDM2 and could directly associate this with increased expression. Furthermore, in all datasets we find fusions to intergenic regions, often spanning multiple cryptic exons that potentially encode neo-antigens. Thus, fusion transcripts other than classical gene-to-gene fusions are prominently present and can be identified using RNA-seq. Conclusion By using the full potential of non–poly(A)-enriched RNA-seq data, sophisticated analysis can reliably identify expressed genomic breakpoints and their transcriptional effects.

INTRODUCTION Little progress has been made in the development of effective new therapies for glioblastoma (GBM) the past decades. Fresh patient-derived GBM cell culture models have become the gold standard for GBM drug discovery and development. One of the major obstacles in identifying novel candidate drugs against GBM remains the blood-brain barrier (BBB). Therefore, it is crucial to select drugs with favourable physicochemical properties to cross BBB and reach the tumour tissue in therapeutically effective concentrations. In current drug repurposing approach, we evaluated available anti-cancer agents in our patient-derived drug screening platform against GBM. METHODS The FDA-approved Oncology Drug Set II library was tested on 45 primary GBM cell cultures. We developed a drug shortlisting pipeline combining efficacy data with pharmacodynamic and pharmacokinetic characteristics of each compound. The therapeutic efficacy of the selected agent was assessed in an orthotopic mouse PDX model, while penetration into the CNS by LC/MS/MS. RESULTS Omacetaxine mepesuccinate (OMA) was ranked as one of the most promising candidates applying our drug selection approach. In vitro, OMA revealed anti-tumour activity at IC50 values well-below reported Cmax plasma values in approximately 80% of GBM cultures. NanoString nCounter analysis, revealed DNA damage repair as the main pathway involved in OMA’s anti-tumour effect. Activation of caspase 3/7 activity and decrease of glioma cell invasiveness were also linked to its anti-tumour effect. In vivo, 1mg/kg dose of OMA was found to reach the brain tumour tissue in concentrations similar to the reported IC50 values in vitro. No adverse reactions were noted and a survival benefit was observed in a proportion of the treated mice. CONCLUSIONS At 1 mg/kg, OMA reaches the tumour brain tissue in therapeutically effective concentrations in mice while a moderate therapeutic benefit was observed. Additional in vivo experiments are ongoing investigating higher dosages of OMA and longer exposure.

BACKGROUND The dismal prognosis of glioblastoma (GBM) patients, with a median survival of less than 15 months despite maximal therapy urgently warrants new therapeutic approaches. Clinical trials employing oncolytic viruses (OVs) have shown encouraging results, however, in each OV clinical trial only a small subset of patients responded to treatment. As inter-tumoral heterogeneity has been the key challenge in treating GBM, we hypothesized that development of an in vitro co-culture model for assessment of viral replication and subsequent immune response might ultimately predict in vivo OV efficacy for individual GBM patients. METHODS 20 patient-derived GBM cell cultures were tested with dose-ranges of two clinically relevant OVs (DNX2401 and rQnestin34.5 V1) to determine the EC50 values (half-maximal effective concentration). To discriminate between responders and non-responders six of these cultures were infected using MOI 10 for DNX2401 and 0.25 for rQnestin34.5V1, and co-cultured with autologous PBMCs. OV-induced changes in gene and protein expression of immune associated genes were assessed using targeted gene expression (NanoString Technology) and ELISA. RESULTS DNX2401 EC50 values ranged from 0.35 to > 600 and from 0.04 to 1.77 for rQnestin34.5V1. Induction of pro-inflammatory cytokines and chemokines differed per virus and per patient. Enhanced OV infection or oncolysis efficiency did not lead to increased levels of IFNγ production. Importantly, longer exposure to dexamethasone prior to PBMC isolation correlated with suppressed IFNγ production. CONCLUSION We established an autologous GBM cells/PBMCs co-culture model which reflects inter-tumoral heterogeneity in terms of cytokine induction and virus-specific changes in gene and protein expression upon OV infection. Immune activation is not directly related to degree of infection or oncolysis and can be hampered by prolonged prior dexamethasone use. These first results support our hypothesis that improved prediction of response rates in oncolytic virotherapy for GBM may require a personalized approach with an in vitro test system.

Background: Patients with locally advanced pancreatic cancer (LAPC) are treated with chemotherapy. In selected cases, stereotactic body radiotherapy (SBRT) can be added to the regimen. We hypothesized that adding an adjuvant containing a heat-killed mycobacterium (IMM-101) to SBRT may lead to beneficial immuno-modulatory effects, thereby improving survival. This study aims to investigate the safety of adding IMM-101 to SBRT and to investigate the immuno-modulatory effects of the combination treatment in the peripheral blood of LAPC patients. Methods: LAPC patients were treated with SBRT (40 Gy) and six intradermal vaccinations of one milligram IMM-101. The primary endpoint was an observed toxicity rate of grade 4 or higher. Targeted gene-expression profiling and multicolor flow cytometry were performed for longitudinal immune-monitoring of the peripheral blood. Results: Twenty patients received study treatment. No treatment-related adverse events of grade 4 or higher occurred. SBRT/IMM-101 treatment induced a transient decrease in different lymphocyte subsets and an increase in CD14+CD16−CD11b+HLA−DRlow myeloid-derived suppressor cells. Importantly, treatment significantly increased activated ICOS+, HLA-DR+ and Ki67+PD1+ T and NK cell frequencies. This was not accompanied by increased levels of most inhibitory markers, such as TIM-3 and LAG-3. Conclusions: Combination therapy with SBRT and a heat-killed mycobacterium vaccine was safe and had an immune-stimulatory effect.

Lower respiratory tract infections (LRTIs) in children are common and, although often mild, a major cause of mortality and hospitalization. Recently, the respiratory microbiome has been associated with both susceptibility and severity of LRTI. In this current study, we combined respiratory microbiome, viral, and clinical data to find associations with the severity of LRTI. Nasopharyngeal aspirates of children aged one month to five years included in the STRAP study (Study to Reduce Antibiotic prescription in childhood Pneumonia), who presented at the emergency department (ED) with fever and cough or dyspnea, were sequenced with nanopore 16S-rRNA gene sequencing and subsequently analyzed with hierarchical clustering to identify respiratory microbiome profiles. Samples were also tested using a panel of 15 respiratory viruses and Mycoplasma pneumoniae, which were analyzed in two groups, according to their reported virulence. The primary outcome was hospitalization, as measure of disease severity. Nasopharyngeal samples were isolated from a total of 167 children. After quality filtering, microbiome results were available for 54 children and virology panels for 158 children. Six distinct genus-dominant microbiome profiles were identified, with Haemophilus-, Moraxella-, and Streptococcus-dominant profiles being the most prevalent. However, these profiles were not found to be significantly associated with hospitalization. At least one virus was detected in 139 (88%) children, of whom 32.4% had co-infections with multiple viruses. Viral co-infections were common for adenovirus, bocavirus, and enterovirus, and uncommon for human metapneumovirus (hMPV) and influenza A virus. The detection of enteroviruses was negatively associated with hospitalization. Virulence groups were not significantly associated with hospitalization. Our data underlines high detection rates and co-infection of viruses in children with respiratory symptoms and confirms the predominant presence of Haemophilus-, Streptococcus-, and Moraxella-dominant profiles in a symptomatic pediatric population at the ED. However, we could not assess significant associations between microbiome profiles and disease severity measures.

Background Bacterial plasmids often carry antibiotic resistance genes and are a significant factor in the spread of antibiotic resistance. The ability to completely assemble plasmid sequences would facilitate the localization of antibiotic resistance genes, the identification of genes that promote plasmid transmission and the accurate tracking of plasmid mobility. However, the complete assembly of plasmid sequences using the currently most widely used sequencing platform (Illumina-based sequencing) is restricted due to the generation of short sequence lengths. The long-read Oxford Nanopore Technologies (ONT) sequencing platform overcomes this limitation. Still, the assembly of plasmid sequence data remains challenging due to software incompatibility with long-reads and the error rate generated using ONT sequencing. Bioinformatics pipelines have been developed for ONT-generated sequencing but require computational skills that frequently are beyond the abilities of scientific researchers. To overcome this challenge, the authors developed ‘WeFaceNano’, a user-friendly Web interFace for rapid assembly and analysis of plasmid DNA sequences generated using the ONT platform. WeFaceNano includes: a read statistics report; two assemblers (Miniasm and Flye); BLAST searching; the detection of antibiotic resistance- and replicon genes and several plasmid visualizations. A user-friendly interface displays the main features of WeFaceNano and gives access to the analysis tools. Results Publicly available ONT sequence data of 21 plasmids were used to validate WeFaceNano, with plasmid assemblages and anti-microbial resistance gene detection being concordant with the published results. Interestingly, the “Flye” assembler with “meta” settings generated the most complete plasmids. Conclusions WeFaceNano is a user-friendly open-source software pipeline suitable for accurate plasmid assembly and the detection of anti-microbial resistance genes in (clinical) samples where multiple plasmids can be present.

The immune response affects tumor biological behavior and progression. The specific immune characteristics of pancreatic ductal adenocarcinoma (PDAC) can determine the metastatic abilities of cancerous cells and the survival of patients. Therefore, it is important to characterize the specific immune landscape in PDAC tissue samples, and the effect of various types of therapy on that immune composition. Previously, a set of marker genes was identified to assess the immune cell composition in different types of cancer tissue samples. However, gene expression and subtypes of immune cells may vary across different types of cancers. The aim of this study was to provide a method to identify immune cells specifically in PDAC tissue samples. The method is based on defining a specific set of marker genes expressed by various immune cells in PDAC samples. A total of 90 marker genes were selected and tested for immune cell type-specific definition in PDAC; including 43 previously used, and 47 newly selected marker genes. The immune cell-type specificity was checked mathematically by calculating the “pairwise similarity” for all candidate genes using the PDAC RNA-sequenced dataset available at The Cancer Genome Atlas. A set of 55 marker genes that identify 22 different immune cell types for PDAC was created. To validate the method and the set of marker genes, an independent mRNA expression dataset of 24 samples of PDAC patients who received various types of (neo)adjuvant treatments was used. The results showed that by applying our method we were able to identify PDAC specific marker genes to characterize immune cell infiltration in tissue samples. The method we described enabled identifying different subtypes of immune cells that were affected by various types of therapy in PDAC patients. In addition, our method can be easily adapted and applied to identify the specific immune landscape in various types of tissue samples.

Epigenetic programming is essential for lineage differentiation, embryogenesis and placentation in early pregnancy. In epigenetic association studies, DNA methylation is often examined in DNA derived from white blood cells, although its validity to other tissues of interest remains questionable. Therefore, we investigated the tissue specificity of epigenome-wide DNA methylation in newborn and placental tissues. Umbilical cord white blood cells (UC-WBC, n = 25), umbilical cord blood mononuclear cells (UC-MNC, n = 10), human umbilical vein endothelial cells (HUVEC, n = 25) and placental tissue (n = 25) were obtained from 36 uncomplicated pregnancies. Genome-wide DNA methylation was measured by the Illumina HumanMethylation450K BeadChip. Using UC-WBC as a reference tissue, we identified 3595 HUVEC tissue-specific differentially methylated regions (tDMRs) and 11,938 placental tDMRs. Functional enrichment analysis showed that HUVEC and placental tDMRs were involved in embryogenesis, vascular development and regulation of gene expression. No tDMRs were identified in UC-MNC. In conclusion, the extensive amount of genome-wide HUVEC and placental tDMRs underlines the relevance of tissue-specific approaches in future epigenetic association studies, or the use of validated representative tissues for a certain disease of interest, if available. To this purpose, we herewith provide a relevant dataset of paired, tissue-specific, genome-wide methylation measurements in newborn tissues.

Summary A major challenge for treating patients with pancreatic ductal adenocarcinoma (PDAC) is the unpredictability of their prognoses due to high heterogeneity. We present Multi-Omics DEep Learning for Prognosis-correlated subtyping (MODEL-P) to identify PDAC subtypes and to predict prognoses of new patients. MODEL-P was trained on autoencoder integrated multi-omics of 146 patients with PDAC together with their survival outcome. Using MODEL-P, we identified two PDAC subtypes with distinct survival outcomes (median survival 10.1 and 22.7 months, respectively, log rank p = 1 × 10−6), which correspond to DNA damage repair and immune response. We rigorously validated MODEL-P by stratifying patients in five independent datasets into these two survival groups and achieved significant survival difference, which is superior to current practice and other subtyping schemas. We believe the subtype-specific signatures would facilitate PDAC pathogenesis discovery, and MODEL-P can provide clinicians the prognoses information in the treatment decision-making to better gauge the benefits versus the risks., Graphical abstract, Highlights • We developed DL-based MODEL-P to identify prognosis-correlated PDAC subtypes • The identified subtypes related to DNA damage repair and immune response processes • MODEL-P stratified patients from independent datasets into distinct survival groups • MODEL-P could be used in clinics to aid treatment decision-making, Biocomputational method; Cancer systems biology; Cancer

Background Long-read sequencing can be applied to generate very long contigs and even completely assembled genomes at relatively low cost and with minimal sample preparation. As a result, long-read sequencing platforms are becoming more popular. In this respect, the Oxford Nanopore Technologies–based long-read sequencing “nanopore" platform is becoming a widely used tool with a broad range of applications and end-users. However, the need to explore and manipulate the complex data generated by long-read sequencing platforms necessitates accompanying specialized bioinformatics platforms and tools to process the long-read data correctly. Importantly, such tools should additionally help democratize bioinformatics analysis by enabling easy access and ease-of-use solutions for researchers. Results The Galaxy platform provides a user-friendly interface to computational command line–based tools, handles the software dependencies, and provides refined workflows. The users do not have to possess programming experience or extended computer skills. The interface enables researchers to perform powerful bioinformatics analysis, including the assembly and analysis of short- or long-read sequence data. The newly developed “NanoGalaxy" is a Galaxy-based toolkit for analysing long-read sequencing data, which is suitable for diverse applications, including de novo genome assembly from genomic, metagenomic, and plasmid sequence reads. Conclusions A range of best-practice tools and workflows for long-read sequence genome assembly has been integrated into a NanoGalaxy platform to facilitate easy access and use of bioinformatics tools for researchers. NanoGalaxy is freely available at the European Galaxy server https://nanopore.usegalaxy.eu with supporting self-learning training material available at https://training.galaxyproject.org.

BackgroundPancreatic ductal adenocarcinoma (PDAC) is notoriously resistant to treatment including checkpoint-blockade immunotherapy. We hypothesized that a bimodal treatment approach consisting of dendritic cell (DC) vaccination to prime tumor-specific T cells, and a strategy to reprogram the desmoplastic tumor microenvironment (TME) would be needed to break tolerance to these pancreatic cancers. As a proof-of-concept, we investigated the efficacy of combined DC vaccination with CD40-agonistic antibodies in a poorly immunogenic murine model of PDAC. Based on the rationale that mesothelioma and pancreatic cancer share a number of tumor associated antigens, the DCs were loaded with either pancreatic or mesothelioma tumor lysates.MethodsImmune-competent mice with subcutaneously or orthotopically growing KrasG12D/+;Trp53R172H/+;Pdx-1-Cre (KPC) PDAC tumors were vaccinated with syngeneic bone marrow-derived DCs loaded with either pancreatic cancer (KPC) or mesothelioma (AE17) lysate and consequently treated with FGK45 (CD40 agonist). Tumor progression was monitored and immune responses in TME and lymphoid organs were analyzed using multicolor flow cytometry and NanoString analyzes.ResultsMesothelioma-lysate loaded DCs generated cross-reactive tumor-antigen-specific T-cell responses to pancreatic cancer and induced delayed tumor outgrowth when provided as prophylactic vaccine. In established disease, combination with stimulating CD40 antibody was necessary to improve survival, while anti-CD40 alone was ineffective. Extensive analysis of the TME showed that anti-CD40 monotherapy did improve CD8 +T cell infiltration, but these essential effector cells displayed hallmarks of exhaustion, including PD-1, TIM-3 and NKG2A. Combination therapy induced a strong change in tumor transcriptome and mitigated the expression of inhibitory markers on CD8 +T cells.ConclusionThese results demonstrate the potency of DC therapy in combination with CD40-stimulation for the treatment of pancreatic cancer and provide directions for near future clinical trials.

CRISPR-Cas9 systems are enriched in human pathogenic bacteria and have been linked to cytotoxicity by an unknown mechanism. Here, we show that upon infection of human cells, Campylobacter jejuni secretes its Cas9 (CjeCas9) nuclease into their cytoplasm. Next, a native nuclear localization signal enables CjeCas9 nuclear entry, where it catalyzes metal-dependent nonspecific DNA cleavage leading to cell death. Compared to CjeCas9, native Cas9 of Streptococcus pyogenes (SpyCas9) is more suitable for guide-dependent editing. However, in human cells, native SpyCas9 may still cause some DNA damage, most likely because of its ssDNA cleavage activity. This side effect can be completely prevented by saturation of SpyCas9 with an appropriate guide RNA, which is only partially effective for CjeCas9. We conclude that CjeCas9 plays an active role in attacking human cells rather than in viral defense. Moreover, these unique catalytic features may therefore make CjeCas9 less suitable for genome editing applications.

Consensus molecular subtyping is an RNA expression-based classification system for colorectal cancer (CRC). Genomic alterations accumulate during CRC pathogenesis, including the premalignant adenoma stage, leading to changes in RNA expression. Only a minority of adenomas progress to malignancies, a transition that is associated with specific DNA copy number aberrations or microsatellite instability (MSI). We aimed to investigate whether colorectal adenomas can already be stratified into consensus molecular subtype (CMS) classes, and whether specific CMS classes are related to the presence of specific DNA copy number aberrations associated with progression to malignancy. RNA sequencing was performed on 62 adenomas and 59 CRCs. MSI status was determined with polymerase chain reaction-based methodology. DNA copy number was assessed by low-coverage DNA sequencing (n = 30) or array-comparative genomic hybridisation (n = 32). Adenomas were classified into CMS classes together with CRCs from the study cohort and from The Cancer Genome Atlas (n = 556), by use of the established CMS classifier. As a result, 54 of 62 (87%) adenomas were classified according to the CMS. The CMS3 ‘metabolic subtype’, which was least common among CRCs, was most prevalent among adenomas (n = 45; 73%). One of the two adenomas showing MSI was classified as CMS1 (2%), the ‘MSI immune’ subtype. Eight adenomas (13%) were classified as the ‘canonical’ CMS2. No adenomas were classified as the ‘mesenchymal’ CMS4, consistent with the fact that adenomas lack invasion-associated stroma. The distribution of the CMS classes among adenomas was confirmed in an independent series. CMS3 was enriched with adenomas at low risk of progressing to CRC, whereas relatively more high-risk adenomas were observed in CMS2. We conclude that adenomas can be stratified into the CMS classes. Considering that CMS1 and CMS2 expression signatures may mark adenomas at increased risk of progression, the distribution of the CMS classes among adenomas is consistent with the proportion of adenomas expected to progress to CRC.

Consensus molecular subtyping is an RNA expression‐based classification system for colorectal cancer (CRC). Genomic alterations accumulate during CRC pathogenesis, including the premalignant adenoma stage, leading to changes in RNA expression. Only a minority of adenomas progress to malignancies, a transition that is associated with specific DNA copy number aberrations or microsatellite instability (MSI). We aimed to investigate whether colorectal adenomas can already be stratified into consensus molecular subtype (CMS) classes, and whether specific CMS classes are related to the presence of specific DNA copy number aberrations associated with progression to malignancy. RNA sequencing was performed on 62 adenomas and 59 CRCs. MSI status was determined with polymerase chain reaction‐based methodology. DNA copy number was assessed by low‐coverage DNA sequencing (n = 30) or array‐comparative genomic hybridisation (n = 32). Adenomas were classified into CMS classes together with CRCs from the study cohort and from The Cancer Genome Atlas (n = 556), by use of the established CMS classifier. As a result, 54 of 62 (87%) adenomas were classified according to the CMS. The CMS3 ‘metabolic subtype’, which was least common among CRCs, was most prevalent among adenomas (n = 45; 73%). One of the two adenomas showing MSI was classified as CMS1 (2%), the ‘MSI immune’ subtype. Eight adenomas (13%) were classified as the ‘canonical’ CMS2. No adenomas were classified as the ‘mesenchymal’ CMS4, consistent with the fact that adenomas lack invasion‐associated stroma. The distribution of the CMS classes among adenomas was confirmed in an independent series. CMS3 was enriched with adenomas at low risk of progressing to CRC, whereas relatively more high‐risk adenomas were observed in CMS2. We conclude that adenomas can be stratified into the CMS classes. Considering that CMS1 and CMS2 expression signatures may mark adenomas at increased risk of progression, the distribution of the CMS classes among adenomas is consistent with the proportion of adenomas expected to progress to CRC. © 2018 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Background: Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive malignancy associated with a poor prognosis. Only 10% of the patients survive longer than five years. So far, factors underlining long-term survivorship in PDAC are not well understood. Therefore, we aimed to identify the key players in the tumor immune microenvironment (TIME) associated with long-term survivorship in PDAC patients. Methods: The immune-related gene expression profiles of surgically resected PDAC tumors of patients who survived and remained recurrence-free of disease for > 3 years (long-term survivors, n=10) were compared to PDAC tumors of patients who had survived ≤ 6 months (short-term survivors, n=10) due to tumor recurrence. Samples were profiled using the nCounter® PanCancer Immune Profiling Panel of NanoString Technology. Validation was performed by spatial analysis of immune cells using the GeoMx™ Digital Spatial Profiler. Results: Tumor-infiltrating B cells were found to be significantly increased in the TIME of long-term survivors by gene expression profiling (p=0.018). The high tumor infiltration of B cells was confirmed by spatial protein profiling (p=0.049). This increase was accompanied by more T cells and antigen-presenting cell infiltration. Moreover, the activated immune cells were found to infiltrate in between tumor cells as well as in stromal areas in long-term survivors. In contrast, the TIME of short-term survivors was characterized by a high density of immunosuppressive cells like CD25 and regulatory T cells infiltrating in a highly fibrotic vicinity. Conclusion: This is the first comprehensive study that connects the immune landscape of gene expression profiles and protein spatial infiltration with the survivorship of PDAC patients. We found higher infiltration of B cells in TIME of long-term survivors which highlights the importance of targeting B cells and B cell-based therapy for future personalized immunotherapy in PDAC patients. Citation Format: Hosein M. Aziz, Lawlaw Saida, Willem de Koning, Andrew Stubbs, Yunlei Li, Casper H. J. van Eijck, Dana A. M. Mustafa. The tumor immune microenvironment is decisive in the survival of pancreatic ductal adenocarcinoma [abstract]. In: Proceedings of the AACR Virtual Special Conference on Pancreatic Cancer; 2021 Sep 29-30. Philadelphia (PA): AACR; Cancer Res 2021;81(22 Suppl):Abstract nr PO-118.

In B-cell precursor acute lymphoblastic leukemia, the initial leukemic cells share the same antigen receptor gene rearrangements. However, due to ongoing rearrangement processes, leukemic cells with different gene rearrangement patterns can develop, resulting in subclone formation. We studied leukemic subclones and their distribution in the bone marrow and peripheral blood at diagnosis. Antigen receptor gene rearrangements (IGH, IGK, TRG, TRD, TRB) were analyzed by next-generation sequencing in seven paired bone marrow samples and five paired bone marrow-peripheral blood samples. Background-thresholds were defined, which enabled identification of leukemic gene rearrangements down to very low levels. Paired bone marrow analysis showed oligoclonality in all 7 patients and up to 34 leukemic clones per patient. Additional analysis of evolutionary-related IGH gene rearrangements revealed up to 171 leukemic clones per patient. Interestingly, overall 86% of all leukemic gene rearrangements, including small subclones, were present in both bone marrow samples (range per patient: 72–100%). Paired bone marrow-peripheral blood analysis showed that 83% of all leukemic gene rearrangements in bone marrow were also found in peripheral blood (range per patient: 81–100%). Remarkably, in the paired bone marrow samples and paired bone marrow-peripheral blood samples the vast majority of leukemic gene rearrangements had a similar frequency (

Antigen Receptor Galaxy (ARGalaxy) is a Web-based tool for analyses and visualization of TCR and BCR sequencing data of 13 species. ARGalaxy consists of four parts: the demultiplex tool, the international ImMunoGeneTics information system (IMGT) concatenate tool, the immune repertoire pipeline, and the somatic hypermutation (SHM) and class switch recombination (CSR) pipeline. Together they allow the analysis of all different aspects of the immune repertoire. All pipelines can be run independently or combined, depending on the available data and the question of interest. The demultiplex tool allows data trimming and demultiplexing, whereas with the concatenate tool multiple IMGT/HighV-QUEST output files can be merged into a single file. The immune repertoire pipeline is an extended version of our previously published ImmunoGlobulin Galaxy (IGGalaxy) virtual machine that was developed to visualize V(D)J gene usage. It allows analysis of both BCR and TCR rearrangements, visualizes CDR3 characteristics (length and amino acid usage) and junction characteristics, and calculates the diversity of the immune repertoire. Finally, ARGalaxy includes the newly developed SHM and CSR pipeline to analyze SHM and/or CSR in BCR rearrangements. It analyzes the frequency and patterns of SHM, Ag selection (including BASELINe), clonality (Change-O), and CSR. The functionality of the ARGalaxy tool is illustrated in several clinical examples of patients with primary immunodeficiencies. In conclusion, ARGalaxy is a novel tool for the analysis of the complete immune repertoire, which is applicable to many patient groups with disturbances in the immune repertoire such as autoimmune diseases, allergy, and leukemia, but it can also be used to address basic research questions in repertoire formation and selection.

Ventricular noncompaction cardiomyopathy is a rare form of congenital cardiomyopathy with increasing evidence of genetic etiology, especially when presenting in childhood. Fetal presentation is rare. We describe a case of fetal hydrops, presenting at 24 weeks gestation and leading to intrapartum death at 26 weeks gestation. Autopsy examination revealed characteristic features of left ventricular noncompaction. A genetic analysis identified a constellation of variants of unknown significance in MYH6, TNNC1, and MYBPC3, genes known to be important in sarcomeric function. Additionally, the variant in MYBPC3 was homozygous. While this case did not demonstrate a conventional single-gene mutation as the cause of the ventricular noncompaction, a broader genomic investigation revealed several variants in sarcomeric genes which may act synergistically to impact cardiac function.

Common variable immune deficiency disorder (CVID) is the most clinically relevant cause of antibody failure. It is a highly heterogeneous disease with different underlying etiologies. CVID has been associated with a quantitative B cell defect, however, little is known about the quality of B cells present. Here, we studied the naive and antigen selected B-cell receptor (BCR) repertoire in 33 CVID patients using next generation sequencing, to investigate B cells quality. Analysis for each individual patient revealed whether they have a defect in immune repertoire formation [V(D)J recombination] or specification (somatic hypermutation, subclass distribution, or selection). The naive BCR repertoire was normal in most of the patients, although alterations in repertoire diversity and the junctions were found in a limited number of patients indicating possible defects in early B-cell development or V(D)J recombination in these patients. In contrast, major differences were found in the antigen selected BCR repertoire. Here, most patients (15/17) showed a reduced frequency of somatic hypermutation (SHM), changes in subclass distribution and/or minor alterations in antigen selection. Together these data show that in our CVID cohort only a small number of patients have a defect in formation of the naive BCR repertoire, whereas the clear majority of patients have disturbances in their antigen selected repertoire, suggesting a defect in repertoire specification in the germinal centers of these patients. This highlights that CVID patients not only have a quantitative B cell defect, but that also the quality of, especially post germinal center B cells, is impaired.

Background Mutations in the γ‐secretase enzyme subunits have been described in multiple kindreds with familial hidradenitis suppurativa (HS). Objective In this study, we report a novel nicastrin (NCSTN) mutation causing HS in a Dutch family. We sought to explore the immunobiological function of NCSTN mutations using data of the Immunological Genome Project. Methods Blood samples of three affected and two unaffected family members were collected. Whole‐genome sequencing was performed using genomic DNA isolated from peripheral blood leucocytes. Sanger sequencing was done to confirm the causative NCSTN variant and the familial segregation. The microarray data set of the Immunological Genome Project was used for thorough dissection of the expression and function of wildtype NCSTN in the immune system. Results In a family consisting of 23 members, we found an autosomal dominant inheritance pattern of HS and detected a novel splice site mutation (c.1912_1915delCAGT) in the NCSTN gene resulting in a frameshift and subsequent premature stop. All affected individuals had HS lesions on non‐flexural and atypical locations. Wildtype NCSTN appears to be upregulated in myeloid cells like monocytes and macrophages, and in mesenchymal cells such as fibroblastic reticular cells and fibroblasts. In addition, within the 25 highest co‐expressed genes with NCSTN we identified CAPNS1,ARNT and PPARD. Conclusion This study reports the identification a novel NCSTN gene splice site mutation which causes familial HS. The associated immunobiological functions of NCSTN and its co‐expressed genes ARNT and PPARD link genetics to the most common environmental and metabolic HS risk factors which are smoking and obesity.

Antigen receptor gene rearrangements are frequently applied as molecular targets for detection of minimal residual disease (MRD) in B-cell precursor acute lymphoblastic leukemia patients. Since such targets may be lost at relapse, appropriate selection of antigen receptor genes as MRD-PCR target is critical. Recently, next-generation sequencing (NGS) - much more sensitive and quantitative than classical PCR-heteroduplex approaches - has been introduced for identification of MRD-PCR targets. We evaluated 42 paired diagnosis-relapse samples by NGS (IGH, IGK, TRG, TRD, and TRB) to evaluate clonal evolution patterns and to design an algorithm for selection of antigen receptor gene rearrangements most likely to remain stable at relapse. Overall, only 393 out of 1446 (27%) clonal rearrangements were stable between diagnosis and relapse. If only index clones with a frequency >5% at diagnosis were taken into account, this number increased to 65%; including only index clones with an absolute read count >10,000, indicating truly major clones, further increased the stability to 84%. Over 90% of index clones at relapse were also present as index clone at diagnosis. Our data provide detailed information about the stability of antigen receptor gene rearrangements, based on which we propose an algorithm for selecting stable MRD-PCR targets, successful in >97% of patients.

Repair of DNA double-strand breaks (DSBs) by the nonhomologous end-joining pathway (NHEJ) is important not only for repair of spontaneous breaks but also for breaks induced in developing lymphocytes during V(D)J (variable [V], diversity [D], and joining [J] genes) recombination of their antigen receptor loci to create a diverse repertoire. Mutations in the NHEJ factor XLF result in extreme sensitivity for ionizing radiation, microcephaly, and growth retardation comparable to mutations in LIG4 and XRCC4, which together form the NHEJ ligation complex. However, the effect on the immune system is variable (mild to severe immunodeficiency) and less prominent than that seen in deficiencies of NHEJ factors ARTEMIS and DNA-dependent protein kinase catalytic subunit, with defects in the hairpin opening step, which is crucial and unique for V(D)J recombination. Therefore, we aimed to study the role of XLF during V(D)J recombination. We obtained clinical data from 9 XLF-deficient patients and performed immune phenotyping and antigen receptor repertoire analysis of immunoglobulin (Ig) and T-cell receptor (TR) rearrangements, using next-generation sequencing in 6 patients. The results were compared with XRCC4 and LIG4 deficiency. Both Ig and TR rearrangements showed a significant decrease in the number of nontemplated (N) nucleotides inserted by terminal deoxynucleotidyl transferase, which resulted in a decrease of 2 to 3 amino acids in the CDR3. Such a reduction in the number of N-nucleotides has a great effect on the junctional diversity, and thereby on the total diversity of the Ig and TR repertoire. This shows that XLF has an important role during V(D)J recombination in creating diversity of the repertoire by stimulating N-nucleotide insertion.

Whole genome sequencing of a severely affected dizygotic twin with an autism spectrum disorder and intellectual disability revealed a compound heterozygous mutation in the HTR7 gene as the only variation not detected in control databases. Each parent carries one allele of the mutation, which is not present in an unaffected stepsister. The HTR7 gene encodes the 5-HT7 serotonin receptor that is involved in brain development, synaptic transmission, and plasticity. The paternally inherited p.W60C variant is situated at an evolutionary conserved nucleotide and predicted damaging by Polyphen2. A mutation akin to the maternally inherited pV286I mutation has been reported to significantly affect the binding characteristics of the receptor. Therefore, the observed sequence alterations provide a first suggestive link between a genetic abnormality in the HTR7 gene and a neurodevelopmental disorder. © 2016 Wiley Periodicals, Inc.

The results from the randomized phase II BELOB trial provided evidence for a potential benefit of bevacizumab (beva), a humanized monoclonal antibody against circulating VEGF-A, when added to CCNU chemotherapy in patients with recurrent glioblastoma (GBM). In this study, we performed gene expression profiling (DASL and RNA-seq) of formalin-fixed, paraffin-embedded tumor material from participants of the BELOB trial to identify patients with recurrent GBM who benefitted most from beva+CCNU treatment. We demonstrate that tumors assigned to the IGS-18 or “classical” subtype and treated with beva+CCNU showed a significant benefit in progression-free survival and a trend toward benefit in overall survival, whereas other subtypes did not exhibit such benefit. In particular, expression of FMO4 and OSBPL3 was associated with treatment response. Importantly, the improved outcome in the beva+CCNU treatment arm was not explained by an uneven distribution of prognostically favorable subtypes as all molecular glioma subtypes were evenly distributed along the different study arms. The RNA-seq analysis also highlighted genetic alterations, including mutations, gene fusions, and copy number changes, within this well-defined cohort of tumors that may serve as useful predictive or prognostic biomarkers of patient outcome. Further validation of the identified molecular markers may enable the future stratification of recurrent GBM patients into appropriate treatment regimens. Cancer Res; 76(3); 525–34. ©2016 AACR.

This paper is a contribution on financial sustainability of hospital biobanks and data warehouses for biospecimens. It aims to discuss new venues for cost models in addition to conventional cost recovery models. It follows the first paper issues on economics of biobanking by Huttin and Liebman, where adaptive platforms where already suggested for biobanks in translational medicine. A case study approach is proposed in order to identify the type of cost drivers that will be needed for optimization of resource allocation in hospital biobanks. It can help the collaboration with designers of ontologies for adaptive platforms, with new data elements on costs and their measurement in different organizational structures.

Illumina and nanopore sequencing technologies are powerful tools that can be used to determine the bacterial composition of complex microbial communities. In this study, we compared nasal microbiota results at genus level using both Illumina and nanopore 16S rRNA gene sequencing. We also monitored the progression of nanopore sequencing in the accurate identification of species, using pure, single species cultures, and evaluated the performance of the nanopore EPI2ME 16S data analysis pipeline. Fifty-nine nasal swabs were sequenced using Illumina MiSeq and Oxford Nanopore 16S rRNA gene sequencing technologies. In addition, five pure cultures of relevant bacterial species were sequenced with the nanopore sequencing technology. The Illumina MiSeq sequence data were processed using bioinformatics modules present in the Mothur software package. Albacore and Guppy base calling, a workflow in nanopore EPI2ME (Oxford Nanopore Technologies&mdash, ONT, Oxford, UK) and an in-house developed bioinformatics script were used to analyze the nanopore data. At genus level, similar bacterial diversity profiles were found, and five main and established genera were identified by both platforms. However, probably due to mismatching of the nanopore sequence primers, the nanopore sequencing platform identified Corynebacterium in much lower abundance compared to Illumina sequencing. Further, when using default settings in the EPI2ME workflow, almost all sequence reads that seem to belong to the bacterial genus Dolosigranulum and a considerable part to the genus Haemophilus were only identified at family level. Nanopore sequencing of single species cultures demonstrated at least 88% accurate identification of the species at genus and species level for 4/5 strains tested, including improvements in accurate sequence read identification when the basecaller Guppy and Albacore, and when flowcell versions R9.4 (Oxford Nanopore Technologies&mdash, ONT, Oxford, UK) and R9.2 (Oxford Nanopore Technologies&mdash, ONT, Oxford, UK) were compared. In conclusion, the current study shows that the nanopore sequencing platform is comparable with the Illumina platform in detection bacterial genera of the nasal microbiota, but the nanopore platform does have problems in detecting bacteria within the genus Corynebacterium. Although advances are being made, thorough validation of the nanopore platform is still recommendable.

Cytotoxic CD8+ T lymphocytes (CTLs) play a major role in protection against chronic viral infections and tumors. Chronic antigen stimulation in both conditions leads to CTL exhaustion and failure of CTLs to efficiently eradicate tumors and virally infected cells. miRNAs are small, non-coding RNA molecules that participate in post-transcriptional gene regulation. This study aimed to determine the core miRNA signature of CTL exhaustion in mouse models of both tumors and chronic infection. Exhausted anti-tumor CTLs were isolated from tumors of mice that received 10^6 mouse mesothelioma cells (AE17sOVA) followed by an adoptive transfer of 10^4 CD45.1+ OT-I cells. Effector anti-viral CTLs were isolated from mice adoptively transferred with 10^4 CD45.1+ OT-I cells or P14 cells and infected with WSN-OVA or LCMV Armstrong respectively. For exhausted anti-viral CTLs, mice were adoptively transferred with 10^4 CD45.1+ P14 cells and infected with LCMV-Cl13. Single cell suspensions were prepared from the different tissues, live CTLs were FACS sorted, RNA isolated and small RNA-seq was performed. Principle component analysis (PCA) of miRNA expression revealed that exhausted anti-tumor CTLs and anti-LCMV Cl13 CTLs cluster apart from effector anti-viral CTLs (WSN-OVA and LCMV Armstrong). Differential expression analysis identified 73 miRNAs commonly shared between the exhausted day anti-tumor CTL and anti-viral CTL. Using RNA sequencing of two different models of exhaustion, distinct RNA profiles common in exhausted CTLs were found. This core miRNA signature can now be used to manipulate the expression of specific miRNA to potentially prevent or revert CTL exhaustion and improve the CTL response against tumors and/or chronic infections.

Background The emergence and spread of antibiotic resistant micro-organisms is a global concern, which is largely attributable to inaccurate prescribing of antibiotics to patients presenting with non-bacterial infections. The use of ‘omics’ technologies for discovery of novel infection related biomarkers combined with novel treatment algorithms offers possibilities for rapidly distinguishing between bacterial and viral infections. This distinction can be particularly important for patients suffering from lower respiratory tract infections (LRTI) and/or sepsis as they represent a significant burden to healthcare systems. Here we present the study details of the TAILORED-Treatment study, an observational, prospective, multi-centre study aiming to generate a multi-parametric model, combining host and pathogen data, for distinguishing between bacterial and viral aetiologies in children and adults with LRTI and/or sepsis. Methods A total number of 1200 paediatric and adult patients aged 1 month and older with LRTI and/or sepsis or a non-infectious disease are recruited from Emergency Departments and hospital wards of seven Dutch and Israeli medical centres. A panel of three experienced physicians adjudicate a reference standard diagnosis for all patients (i.e., bacterial or viral infection) using all available clinical and laboratory information, including a 28-day follow-up assessment. Nasal swabs and blood samples are collected for multi-omics investigations including host RNA and protein biomarkers, nasal microbiota profiling, host genomic profiling and bacterial proteomics. Simplified data is entered into a custom-built database in order to develop a multi-parametric model and diagnostic tools for differentiating between bacterial and viral infections. The predictions from the model will be compared with the consensus diagnosis in order to determine its accuracy. Discussion The TAILORED-Treatment study will provide new insights into the interplay between the host and micro-organisms. New host- or pathogen-related biomarkers will be used to generate a multi-parametric model for distinguishing between bacterial and viral infections. This model will be helpful to better guide antimicrobial therapy for patients with LRTI and sepsis. This study has the potential to improve patient care, reduce unnecessary antibiotic prescribing and will contribute positively to institutional, national and international healthcare economics. Trial Registration NCT02025699. Registration Date: January, 1, 2014. Electronic supplementary material The online version of this article (10.1186/s12879-018-3300-9) contains supplementary material, which is available to authorized users.

Influenza virus outbreaks remain a serious threat to public health. A greater understanding of how cells targeted by the virus respond to the infection can provide insight into the pathogenesis of disease. Here we examined the transcriptional profile of in vivo-infected and uninfected type 2 alveolar epithelial cells (AEC) in the lungs of influenza virus-infected mice. We show for the first time the unique gene expression profiles induced by the in vivo infection of AEC as well as the transcriptional response of uninfected bystander cells. This work allows us to distinguish the direct and indirect effects of infection at the cellular level. Transcriptome analysis revealed that although directly infected and bystander AEC from infected animals shared many transcriptome changes compared to AEC from uninfected animals, directly infected cells produce more interferon and express lower levels of Wnt signaling-associated transcripts, while concurrently expressing more transcripts associated with cell death pathways, than bystander uninfected AEC. The Wnt signaling pathway was downregulated in both in vivo-infected AEC and in vitro-infected human lung epithelial A549 cells. Wnt signaling did not affect type I and III interferon production by infected A549 cells. Our results reveal unique transcriptional changes that occur within infected AEC and show that influenza virus downregulates Wnt signaling. In light of recent findings that Wnt signaling is essential for lung epithelial stem cells, our findings reveal a mechanism by which influenza virus may affect host lung repair. IMPORTANCE Influenza virus infection remains a major public health problem. Utilizing a recombinant green fluorescent protein-expressing influenza virus, we compared the in vivo transcriptomes of directly infected and uninfected bystander cells from infected mouse lungs and discovered many pathways uniquely regulated in each population. The Wnt signaling pathway was downregulated in directly infected cells and was shown to affect virus but not interferon production. Our study is the first to discern the in vivo transcriptome changes induced by direct viral infection compared to mere exposure to the lung inflammatory milieu and highlight the downregulation of Wnt signaling. This downregulation has important implications for understanding influenza virus pathogenesis, as Wnt signaling is critical for lung epithelial stem cells and lung epithelial cell differentiation. Our findings reveal a mechanism by which influenza virus may affect host lung repair and suggest interventions that prevent damage or accelerate recovery of the lung.

Background In vitro scratch assays have been widely used to study the influence of bioactive substances on the processes of cell migration and proliferation that are involved in re-epithelialization. The development of high-throughput microscopy and image analysis has enabled scratch assays to become compatible with high-throughput research. However, effective processing and in-depth analysis of such high-throughput image datasets are far from trivial and require integration of multiple image processing and data extraction software tools. Findings We developed and implemented a kinetic re-epithelialization analysis pipeline (KREAP) in Galaxy. The KREAP toolbox incorporates freely available image analysis tools and automatically performs image segmentation and feature extraction of each image series, followed by automatic quantification of cells inside and outside the scratched area over time. The enumeration of infiltrating cells over time is modeled to extract three biologically relevant parameters that describe re-epithelialization kinetics. The output of the tools is organized, displayed, and saved in the Galaxy environment for future reference. Conclusions The KREAP toolbox in Galaxy provides an open-source, easy-to-use, web-based platform for reproducible image processing and data analysis of high-throughput scratch assays. The KREAP toolbox could assist a broad scientific community in the discovery of compounds that are able to modulate re-epithelialization kinetics.

Background The determination of microbial communities using the mothur tool suite (https://www.mothur.org) is well established. However, mothur requires bioinformatics-based proficiency in order to perform calculations via the command-line. Galaxy is a project dedicated to providing a user-friendly web interface for such command-line tools (https://galaxyproject.org/). Results We have integrated the full set of 125+ mothur tools into Galaxy as the Galaxy mothur Toolset (GmT) and provided a set of workflows to perform end-to-end 16S rRNA gene analyses and integrate with third-party visualization and reporting tools. We demonstrate the utility of GmT by analyzing the mothur MiSeq standard operating procedure (SOP) dataset (https://www.mothur.org/wiki/MiSeq_SOP). Conclusions GmT is available from the Galaxy Tool Shed, and a workflow definition file and full Galaxy training manual for the mothur SOP have been created. A Docker image with a fully configured GmT Galaxy is also available.