Most patients diagnosed with resected pancreatic adenocarcinoma (PDAC) survive less than 5 years, but a minor subset survives longer. Here, we dissect the role of the tumor microbiota and the immune system in influencing long-term survival. Using 16S rRNA gene sequencing, we analyzed the tumor microbiome composition in PDAC patients with short-term survival (STS) and long-term survival (LTS). We found higher alpha-diversity in the tumor microbiome of LTS patients and identified an intra-tumoral microbiome signature (Pseudoxanthomonas-Streptomyces-Saccharopolyspora-Bacillus clausii) highly predictive of long-term survivorship in both discovery and validation cohorts. Through human-into-mice fecal microbiota transplantation (FMT) experiments from STS, LTS, or control donors, we were able to differentially modulate the tumor microbiome and affect tumor growth as well as tumor immune infiltration. Our study demonstrates that PDAC microbiome composition, which cross-talks to the gut microbiome, influences the host immune response and natural history of the disease.

Another benefit of dietary fiber The gut microbiome can modulate the immune system and influence the therapeutic response of cancer patients, yet the mechanisms underlying the effects of microbiota are presently unclear. Spencer et al . add to our understanding of how dietary habits affect microbiota and clinical outcomes to immunotherapy. In an observational study, the researchers found that melanoma patients reporting high fiber (prebiotic) consumption had a better response to checkpoint inhibitor immunotherapy compared with those patients reporting a low-fiber diet. The most marked benefit was observed for those patients reporting a combination of high fiber consumption and no use of over-the-counter probiotic supplements. These findings provide early insights as to how diet-related factors may influence the immune response. —PNK

In this paper, we propose a Bayesian approach to inference on multiple Gaussian graphical models. Specifically, we address the problem of inferring multiple undirected networks in situations where some of the networks may be unrelated, while others share common features. We link the estimation of the graph structures via a Markov random field (MRF) prior which encourages common edges. We learn which sample groups have a shared graph structure by placing a spike-and-slab prior on the parameters that measure network relatedness. This approach allows us to share information between sample groups, when appropriate, as well as to obtain a measure of relative network similarity across groups. Our modeling framework incorporates relevant prior knowledge through an edge-specific informative prior and can encourage similarity to an established network. Through simulations, we demonstrate the utility of our method in summarizing relative network similarity and compare its performance against related methods. We find improved accuracy of network estimation, particularly when the sample sizes within each subgroup are moderate. We also illustrate the application of our model to infer protein networks for various cancer subtypes and under different experimental conditions.

Gene co-expression networks capture biologically important patterns in gene expression data, enabling functional analyses of genes, discovery of biomarkers, and interpretation of genetic variants. Most network analyses to date have been limited to assessing correlation between total gene expression levels in a single tissue or small sets of tissues. Here, we built networks that additionally capture the regulation of relative isoform abundance and splicing, along with tissue-specific connections unique to each of a diverse set of tissues. We used the Genotype-Tissue Expression (GTEx) project v6 RNA sequencing data across 50 tissues and 449 individuals. First, we developed a framework called Transcriptome-Wide Networks (TWNs) for combining total expression and relative isoform levels into a single sparse network, capturing the interplay between the regulation of splicing and transcription. We built TWNs for 16 tissues and found that hubs in these networks were strongly enriched for splicing and RNA binding genes, demonstrating their utility in unraveling regulation of splicing in the human transcriptome. Next, we used a Bayesian biclustering model that identifies network edges unique to a single tissue to reconstruct Tissue-Specific Networks (TSNs) for 26 distinct tissues and 10 groups of related tissues. Finally, we found genetic variants associated with pairs of adjacent nodes in our networks, supporting the estimated network structures and identifying 20 genetic variants with distant regulatory impact on transcription and splicing. Our networks provide an improved understanding of the complex relationships of the human transcriptome across tissues.

Purpose: To develop a head and neck normal structures auto-contouring tool that could be used to automatically detect the errors in auto-contours from a clinically-validated auto-contouring tool. Methods: An auto-contouring tool based on convolutional neural networks (CNN) was developed for 16 normal structures of the head and neck and tested to identify the contour errors from a clinically-validated multi-atlas-based auto-contouring system (MACS). The CT scans and clinical contours from 3495 patients were semi-automatically curated and used to train and validate the CNN-based auto-contouring tool. The final accuracy of the tool was evaluated by calculating the Sørensen-Dice similarity coefficients (DSC) and Hausdorff distances between the automatically generated contours and physician-drawn contours on 174 internal and 24 external CT scans. Lastly, the CNN-based tool was evaluated on 60 patients’ CT scans to investigate the possibility to detect contouring failures. The contouring failures on these patients were classified as either minor or major errors. The criteria to detect contouring errors were determined by analyzing the DSC between the CNN- and MACS-based contours under two independent scenarios: 1. contours with minor error are clinically acceptable and 2. contours with minor errors are clinically unacceptable. Results: The average DSC and Hausdorff distance of our CNN-based tool were 98.4%/1.23cm for brain, 89.1%/0.42cm for eyes, 86.8%/1.28cm for mandible, 86.4%/0.88cm for brainstem, 83.4%/0.71cm for spinal cord, 82.7%/1.37cm for parotids, 80.7%/1.08cm for esophagus, 71.7%/0.39cm for lenses, 68.6%/0.72 for optic nerves, 66.4%/0.46cm for cochleas, and 40.7%/0.96cm for optic chiasm. With the error detection tool, the proportions of the clinically unacceptable MACS contours that were correctly detected were 0.99/0.80 on average except for the optic chiasm, when contours with minor errors are clinically acceptable/unacceptable respectively. The proportions of the clinically acceptable MACS contours that were correctly detected were 0.81/0.60 on average except for the optic chiasm, when contours with minor errors are clinically acceptable/unacceptable respectively. Conclusion: Our CNN-based auto-contouring tool performed well on both the publically-available and the internal datasets. Furthermore, our results show that CNN-based algorithms are able to identify ill-defined contours from a clinically-validated and used multi-atlas-based auto-contouring tool. Therefore, our CNN-based tool can effectively perform automatic verification of MACS contours.

BackgroundFecal microbiota transplantation (FMT) via colonoscopy or enema has become a commonly used treatment of recurrent C. difficile infection (CDI).AimsTo compare the safety and preliminary efficacy of orally administered lyophilized microbiota product compared with frozen product by enema.MethodsIn a single center, adults with ≥ 3 episodes of recurrent CDI were randomized to receive encapsulated lyophilized fecal microbiota from 100–200 g of donor feces (n = 31) or frozen FMT from 100 g of donor feces (n = 34) by enema. Safety during the three months post FMT was the primary study objective. Prevention of CDI recurrence during the 60 days after FMT was a secondary objective. Fecal microbiome changes were examined in first 39 subjects studied.ResultsAdverse experiences were commonly seen in equal frequency in both groups and did not appear to relate to the route of delivery of FMT. CDI recurrence was prevented in 26 of 31 (84%) subjects randomized to capsules and in 30 of 34 (88%) receiving FMT by enema (p = 0.76). Both products normalized fecal microbiota diversity while the lyophilized orally administered product was less effective in repleting Bacteroidia and Verrucomicrobia classes compared to frozen product via enema.ConclusionsThe route of delivery, oral or rectal, did not influence adverse experiences in FMT. In preliminary evaluation, both routes appeared to show equivalent efficacy, although the dose may need to be higher for lyophilized product. Spore-forming bacteria appear to be the most important engrafting organisms in FMT by the oral route using lyophilized product.Trial registrationClinicalTrials.gov NCT02449174

A non-immunogenic tumor microenvironment (TME) is a significant barrier to immune checkpoint blockade (ICB) response. The impact of Polybromo-1 (PBRM1) on TME and response to ICB in renal cell carcinoma (RCC) remains to be resolved. Here we show that PBRM1/Pbrm1 deficiency reduces the binding of brahma-related gene 1 (BRG1) to the IFNγ receptor 2 (Ifngr2) promoter, decreasing STAT1 phosphorylation and the subsequent expression of IFNγ target genes. An analysis of 3 independent patient cohorts and of murine pre-clinical models reveals that PBRM1 loss is associated with a less immunogenic TME and upregulated angiogenesis. Pbrm1 deficient Renca subcutaneous tumors in mice are more resistance to ICB, and a retrospective analysis of the IMmotion150 RCC study also suggests that PBRM1 mutation reduces benefit from ICB. Our study sheds light on the influence of PBRM1 mutations on IFNγ-STAT1 signaling and TME, and can inform additional preclinical and clinical studies in RCC.

The genetic basis of multiple phenotypes such as gene expression, metabolite levels, or imaging features is often investigated by testing a large collection of hypotheses, probing the existence of association between each of the traits and hundreds of thousands of genotyped variants. Appropriate multiplicity adjustment is crucial to guarantee replicability of findings, and the false discovery rate (FDR) is frequently adopted as a measure of global error. In the interest of interpretability, results are often summarized so that reporting focuses on variants discovered to be associated to some phenotypes. We show that applying FDR‐controlling procedures on the entire collection of hypotheses fails to control the rate of false discovery of associated variants as well as the expected value of the average proportion of false discovery of phenotypes influenced by such variants. We propose a simple hierarchical testing procedure that allows control of both these error rates and provides a more reliable basis for the identification of variants with functional effects. We demonstrate the utility of this approach through simulation studies comparing various error rates and measures of power for genetic association studies of multiple traits. Finally, we apply the proposed method to identify genetic variants that impact flowering phenotypes in Arabidopsis thaliana, expanding the set of discoveries.