Biofilms and core pathogens shape the tumor microenvironment and immune phenotype in colorectal cancer.

Extensive research has explored the role of gut microbiota in colorectal cancer (CRC). Nonetheless, metatranscriptomic studies investigating the in situ functional implications of host-microbe interactions in CRC are scarce. Therefore, we characterized the influence of CRC core pathogens and biofilms on the tumor microenvironment (TME) in 40 CRC, paired normal, and healthy tissue biopsies using fluorescence in situ hybridization (FISH) and dual-RNA sequencing. FISH revealed that Fusobacterium spp. was associated with increased bacterial biomass and inflammatory response in CRC samples. Dual-RNA sequencing demonstrated increased expression of pro-inflammatory cytokines, defensins, matrix-metalloproteases, and immunomodulatory factors in CRC samples with high bacterial activity. In addition, bacterial activity correlated with the infiltration of several immune cell subtypes, including M2 macrophages and regulatory T-cells in CRC samples. Specifically, Bacteroides fragilis and Fusobacterium nucleatum correlated with the infiltration of neutrophils and CD4+ T-cells, respectively. The collective bacterial activity/biomass appeared to exert a more significant influence on the TME than core pathogens, underscoring the intricate interplay between gut microbiota and CRC. These results emphasize how biofilms and core pathogens shape the immune phenotype and TME in CRC while highlighting the need to extend the bacterial scope beyond CRC pathogens to advance our understanding and identify treatment targets.

N-of-one differential gene expression without control samples using a deep generative model.

Differential analysis of bulk RNA-seq data often suffers from lack of good controls. Here, we present a generative model that replaces controls, trained solely on healthy tissues. The unsupervised model learns a low-dimensional representation and can identify the closest normal representation for a given disease sample. This enables control-free, single-sample differential expression analysis. In breast cancer, we demonstrate how our approach selects marker genes and outperforms a state-of-the-art method. Furthermore, significant genes identified by the model are enriched in driver genes across cancers. Our results show that the in silico closest normal provides a more favorable comparison than control samples.

Proteomic Profiling of Colorectal Adenomas Identifies a Predictive Risk Signature for Development of Metachronous Advanced Colorectal Neoplasia.

BACKGROUND & AIMS: Colonic adenomatous polyps, or adenomas, are frequent precancerous lesions and the origin of most cases of colorectal adenocarcinoma. However, we know from epidemiologic studies that although most colorectal cancers (CRCs) originate from adenomas, only a small fraction of adenomas (3%-5%) ever progress to cancer. At present, there are no molecular markers to guide follow-up surveillance programs. METHODS: We profiled, by mass spectrometry-based proteomics combined with machine learning analysis, a selected cohort of formalin-fixed, paraffin-embedded high-grade (HG) adenomas with long clinical follow-up, collected as part of the Danish national screening program. We grouped subjects in the cohort according to their subsequent history of findings: a nonmetachronous advanced neoplasia group (G0), with no new HG adenomas or CRCs up to 10 years after polypectomy, and a metachronous advanced neoplasia group (G1) where individuals developed a new HG adenoma or CRC within 5 years of diagnosis. RESULTS: We generated a proteome dataset from 98 selected HG adenoma samples, including 20 technical replicates, of which 45 samples belonged to the nonmetachronous advanced neoplasia group and 53 to the metachronous advanced neoplasia group. The clear distinction of these 2 groups seen in a uniform manifold approximation and projection plot indicated that the information contained within the abundance of the ∼5000 proteins was sufficient to predict the future occurrence of HG adenomas or development of CRC. CONCLUSIONS: We performed an in-depth analysis of quantitative proteomic data from 98 resected adenoma samples using various novel algorithms and statistical packages and found that their proteome can predict development of metachronous advanced lesions and progression several years in advance.

Emergence and radiation of distemper viruses in terrestrial and marine mammals.

Canine distemper virus (CDV) and phocine distemper virus (PDV) are major pathogens to terrestrial and marine mammals. Yet little is known about the timing and geographical origin of distemper viruses and to what extent it was influenced by environmental change and human activities. To address this, we (i) performed the first comprehensive time-calibrated phylogenetic analysis of the two distemper viruses, (ii) mapped distemper antibody and virus detection data from marine mammals collected between 1972 and 2018, and (iii) compiled historical reports on distemper dating back to the eighteenth century. We find that CDV and PDV diverged in the early seventeenth century. Modern CDV strains last shared a common ancestor in the nineteenth century with a marked radiation during the 1930s-1950s. Modern PDV strains are of more recent origin, diverging in the 1970s-1980s. Based on the compiled information on distemper distribution, the diverse host range of CDV and basal phylogenetic placement of terrestrial morbilliviruses, we hypothesize a terrestrial CDV-like ancestor giving rise to PDV in the North Atlantic. Moreover, given the estimated timing of distemper origin and radiation, we hypothesize a prominent role of environmental change such as the Little Ice Age, and human activities like globalization and war in distemper virus evolution.

High-throughput proteomics of breast cancer interstitial fluid: identification of tumor subtype-specific serologically relevant biomarkers.

Despite significant advancements in breast cancer (BC) research, clinicians lack robust serological protein markers for accurate diagnostics and tumor stratification. Tumor interstitial fluid (TIF) accumulates aberrantly externalized proteins within the local tumor space, which can potentially gain access to the circulatory system. As such, TIF may represent a valuable starting point for identifying relevant tumor-specific serological biomarkers. The aim of the study was to perform comprehensive proteomic profiling of TIF to identify proteins associated with BC tumor status and subtype. A liquid chromatography tandem mass spectrometry (LC-MS/MS) analysis of 35 TIFs of three main subtypes: luminal (19), Her2 (4), and triple-negative (TNBC) (12) resulted in the identification of > 8800 proteins. Unsupervised hierarchical clustering segregated the TIF proteome into two major clusters, luminal and TNBC/Her2 subgroups. High-grade tumors enriched with tumor infiltrating lymphocytes (TILs) were also stratified from low-grade tumors. A consensus analysis approach, including differential abundance analysis, selection operator regression, and random forest returned a minimal set of 24 proteins associated with BC subtypes, receptor status, and TIL scoring. Among them, a panel of 10 proteins, AGR3, BCAM, CELSR1, MIEN1, NAT1, PIP4K2B, SEC23B, THTPA, TMEM51, and ULBP2, was found to stratify the tumor subtype-specific TIFs. In particular, upregulation of BCAM and CELSR1 differentiates luminal subtypes, while upregulation of MIEN1 differentiates Her2 subtypes. Immunohistochemistry analysis showed a direct correlation between protein abundance in TIFs and intratumor expression levels for all 10 proteins. Sensitivity and specificity were estimated for this protein panel by using an independent, comprehensive breast tumor proteome dataset. The results of this analysis strongly support our data, with eight of the proteins potentially representing biomarkers for stratification of BC subtypes. Five of the most representative proteomics databases currently available were also used to estimate the potential for these selected proteins to serve as putative serological markers.

Secreted breast tumor interstitial fluid microRNAs and their target genes are associated with triple-negative breast cancer, tumor grade, and immune infiltration.

BACKGROUND: Studies on tumor-secreted microRNAs point to a functional role of these in cellular communication and reprogramming of the tumor microenvironment. Uptake of tumor-secreted microRNAs by neighboring cells may result in the silencing of mRNA targets and, in turn, modulation of the transcriptome. Studying miRNAs externalized from tumors could improve cancer patient diagnosis and disease monitoring and help to pinpoint which miRNA-gene interactions are central for tumor properties such as invasiveness and metastasis. METHODS: Using a bioinformatics approach, we analyzed the profiles of secreted tumor and normal interstitial fluid (IF) microRNAs, from women with breast cancer (BC). We carried out differential abundance analysis (DAA), to obtain miRNAs, which were enriched or depleted in IFs, from patients with different clinical traits. Subsequently, miRNA family enrichment analysis was performed to assess whether any families were over-represented in the specific sets. We identified dysregulated genes in tumor tissues from the same cohort of patients and constructed weighted gene co-expression networks, to extract sets of co-expressed genes and co-abundant miRNAs. Lastly, we integrated miRNAs and mRNAs to obtain interaction networks and supported our findings using prediction tools and cancer gene databases. RESULTS: Network analysis showed co-expressed genes and miRNA regulators, associated with tumor lymphocyte infiltration. All of the genes were involved in immune system processes, and many had previously been associated with cancer immunity. A subset of these, BTLA, CXCL13, IL7R, LAMP3, and LTB, was linked to the presence of tertiary lymphoid structures and high endothelial venules within tumors. Co-abundant tumor interstitial fluid miRNAs within this network, including miR-146a and miR-494, were annotated as negative regulators of immune-stimulatory responses. One co-expression network encompassed differences between BC subtypes. Genes differentially co-expressed between luminal B and triple-negative breast cancer (TNBC) were connected with sphingolipid metabolism and predicted to be co-regulated by miR-23a. Co-expressed genes and TIF miRNAs associated with tumor grade were BTRC, CHST1, miR-10a/b, miR-107, miR-301a, and miR-454. CONCLUSION: Integration of IF miRNAs and mRNAs unveiled networks associated with patient clinicopathological traits, and underlined molecular mechanisms, specific to BC sub-groups. Our results highlight the benefits of an integrative approach to biomarker discovery, placing secreted miRNAs within a biological context.

CAncer bioMarker Prediction Pipeline (CAMPP)-A standardized framework for the analysis of quantitative biological data.

With the improvement of -omics and next-generation sequencing (NGS) methodologies, along with the lowered cost of generating these types of data, the analysis of high-throughput biological data has become standard both for forming and testing biomedical hypotheses. Our knowledge of how to normalize datasets to remove latent undesirable variances has grown extensively, making for standardized data that are easily compared between studies. Here we present the CAncer bioMarker Prediction Pipeline (CAMPP), an open-source R-based wrapper (https://github.com/ELELAB/CAncer-bioMarker-Prediction-Pipeline -CAMPP) intended to aid bioinformatic software-users with data analyses. CAMPP is called from a terminal command line and is supported by a user-friendly manual. The pipeline may be run on a local computer and requires little or no knowledge of programming. To avoid issues relating to R-package updates, a renv .lock file is provided to ensure R-package stability. Data-management includes missing value imputation, data normalization, and distributional checks. CAMPP performs (I) k-means clustering, (II) differential expression/abundance analysis, (III) elastic-net regression, (IV) correlation and co-expression network analyses, (V) survival analysis, and (VI) protein-protein/miRNA-gene interaction networks. The pipeline returns tabular files and graphical representations of the results. We hope that CAMPP will assist in streamlining bioinformatic analysis of quantitative biological data, whilst ensuring an appropriate bio-statistical framework.

Interpreting pathways to discover cancer driver genes with Moonlight.

Cancer driver gene alterations influence cancer development, occurring in oncogenes, tumor suppressors, and dual role genes. Discovering dual role cancer genes is difficult because of their elusive context-dependent behavior. We define oncogenic mediators as genes controlling biological processes. With them, we classify cancer driver genes, unveiling their roles in cancer mechanisms. To this end, we present Moonlight, a tool that incorporates multiple -omics data to identify critical cancer driver genes. With Moonlight, we analyze 8000+ tumor samples from 18 cancer types, discovering 3310 oncogenic mediators, 151 having dual roles. By incorporating additional data (amplification, mutation, DNA methylation, chromatin accessibility), we reveal 1000+ cancer driver genes, corroborating known molecular mechanisms. Additionally, we confirm critical cancer driver genes by analysing cell-line datasets. We discover inactivation of tumor suppressors in intron regions and that tissue type and subtype indicate dual role status. These findings help explain tumor heterogeneity and could guide therapeutic decisions.

Alterations of the interactome of Bcl-2 proteins in breast cancer at the transcriptional, mutational and structural level.

Apoptosis is an essential defensive mechanism against tumorigenesis. Proteins of the B-cell lymphoma-2 (Bcl-2) family regulate programmed cell death by the mitochondrial apoptosis pathway. In response to intracellular stress, the apoptotic balance is governed by interactions of three distinct subgroups of proteins; the activator/sensitizer BH3 (Bcl-2 homology 3)-only proteins, the pro-survival, and the pro-apoptotic executioner proteins. Changes in expression levels, stability, and functional impairment of pro-survival proteins can lead to an imbalance in tissue homeostasis. Their overexpression or hyperactivation can result in oncogenic effects. Pro-survival Bcl-2 family members carry out their function by binding the BH3 short linear motif of pro-apoptotic proteins in a modular way, creating a complex network of protein-protein interactions. Their dysfunction enables cancer cells to evade cell death. The critical role of Bcl-2 proteins in homeostasis and tumorigenesis, coupled with mounting insight in their structural properties, make them therapeutic targets of interest. A better understanding of gene expression, mutational profile, and molecular mechanisms of pro-survival Bcl-2 proteins in different cancer types, could help to clarify their role in cancer development and may guide advancement in drug discovery. Here, we shed light on the pro-survival Bcl-2 proteins in breast cancer using different bioinformatic approaches, linking -omics with structural data. We analyzed the changes in the expression of the Bcl-2 proteins and their BH3-containing interactors in breast cancer samples. We then studied, at the structural level, a selection of interactions, accounting for effects induced by mutations found in the breast cancer samples. We find two complexes between the up-regulated Bcl2A1 and two down-regulated BH3-only candidates (i.e., Hrk and Nr4a1) as targets associated with reduced apoptosis in breast cancer samples for future experimental validation. Furthermore, we predict L99R, M75R as damaging mutations altering protein stability, and Y120C as a possible allosteric mutation from an exposed surface to the BH3-binding site.

Distinct signatures of lung cancer types: aberrant mucin O-glycosylation and compromised immune response.

BACKGROUND: Genomic initiatives such as The Cancer Genome Atlas (TCGA) contain data from -omics profiling of thousands of tumor samples, which may be used to decipher cancer signaling, and related alterations. Managing and analyzing data from large-scale projects, such as TCGA, is a demanding task. It is difficult to dissect the high complexity hidden in genomic data and to account for inter-tumor heterogeneity adequately. METHODS: In this study, we used a robust statistical framework along with the integration of diverse bioinformatic tools to analyze next-generation sequencing data from more than 1000 patients from two different lung cancer subtypes, i.e., the lung adenocarcinoma (LUAD) and the squamous cell carcinoma (LUSC). RESULTS: We used the gene expression data to identify co-expression modules and differentially expressed genes to discriminate between LUAD and LUSC. We identified a group of genes which could act as specific oncogenes or tumor suppressor genes in one of the two lung cancer types, along with two dual role genes. Our results have been validated against other transcriptomics data of lung cancer patients. CONCLUSIONS: Our integrative approach allowed us to identify two key features: a substantial up-regulation of genes involved in O-glycosylation of mucins in LUAD, and a compromised immune response in LUSC. The immune-profile associated with LUSC might be linked to the activation of three oncogenic pathways, which promote the evasion of the antitumor immune response. Collectively, our results provide new future directions for the design of target therapies in lung cancer.

N-glycan signatures identified in tumor interstitial fluid and serum of breast cancer patients: association with tumor biology and clinical outcome.

Particular N-glycan structures are known to be associated with breast malignancies by coordinating various regulatory events within the tumor and corresponding microenvironment, thus implying that N-glycan patterns may be used for cancer stratification and as predictive or prognostic biomarkers. However, the association between N-glycans secreted by breast tumor and corresponding clinical relevance remain to be elucidated. We profiled N-glycans by HILIC UPLC across a discovery dataset composed of tumor interstitial fluids (TIF, n = 85), paired normal interstitial fluids (NIF, n = 54) and serum samples (n = 28) followed by independent evaluation, with the ultimate goal of identifying tumor-related N-glycan patterns in blood of patients with breast cancer. The segregation of N-linked oligosaccharides revealed 33 compositions, which exhibited differential abundances between TIF and NIF. TIFs were depleted of bisecting N-glycans, which are known to play essential roles in tumor suppression. An increased level of simple high mannose N-glycans in TIF strongly correlated with the presence of tumor infiltrating lymphocytes within tumor. At the same time, a low level of highly complex N-glycans in TIF inversely correlated with the presence of infiltrating lymphocytes within tumor. Survival analysis showed that patients exhibiting increased TIF abundance of GP24 had better outcomes, whereas low levels of GP10, GP23, GP38, and coreF were associated with poor prognosis. Levels of GP1, GP8, GP9, GP14, GP23, GP28, GP37, GP38, and coreF were significantly correlated between TIF and paired serum samples. Cross-validation analysis using an independent serum dataset supported the observed correlation between TIF and serum, for five of nine N-glycan groups: GP8, GP9, GP14, GP23, and coreF. Collectively, our results imply that profiling of N-glycans from proximal breast tumor fluids is a promising strategy for determining tumor-derived glyco-signature(s) in the blood. N-glycans structures validated in our study may serve as novel biomarkers to improve the diagnostic and prognostic stratification of patients with breast cancer.

Characterization of the enhancer and promoter landscape of inflammatory bowel disease from human colon biopsies.

Inflammatory bowel disease (IBD) is a chronic intestinal disorder, with two main types: Crohn's disease (CD) and ulcerative colitis (UC), whose molecular pathology is not well understood. The majority of IBD-associated SNPs are located in non-coding regions and are hard to characterize since regulatory regions in IBD are not known. Here we profile transcription start sites (TSSs) and enhancers in the descending colon of 94 IBD patients and controls. IBD-upregulated promoters and enhancers are highly enriched for IBD-associated SNPs and are bound by the same transcription factors. IBD-specific TSSs are associated to genes with roles in both inflammatory cascades and gut epithelia while TSSs distinguishing UC and CD are associated to gut epithelia functions. We find that as few as 35 TSSs can distinguish active CD, UC, and controls with 85% accuracy in an independent cohort. Our data constitute a foundation for understanding the molecular pathology, gene regulation, and genetics of IBD.

A phosphorylation-motif for tuneable helix stabilisation in intrinsically disordered proteins - Lessons from the sodium proton exchanger 1 (NHE1).

Intrinsically disordered proteins (IDPs) are involved in many pivotal cellular processes including phosphorylation and signalling. The structural and functional effects of phosphorylation of IDPs remain poorly understood and difficult to predict. Thus, a need exists to identify motifs that confer phosphorylation-dependent perturbation of the local preferences for forming e.g. helical structures as well as motifs that do not. The disordered distal tail of the Na+/H+ exchanger 1 (NHE1) is six-times phosphorylated (S693, S723, S726, S771, T779, S785) by the mitogen activated protein kinase 2 (MAPK1, ERK2). Using NMR spectroscopy, we found that two out of those six phosphorylation sites had a stabilizing effect on transient helices. One of these was further investigated by circular dichroism and NMR spectroscopy as well as by molecular dynamic simulations, which confirmed the stabilizing effect and resulted in the identification of a short linear motif for helix stabilisation: [S/T]-P-{3}-[R/K] where [S/T] is the phosphorylation-site. By analysing IDP and phosphorylation site databases we found that the motif is significantly enriched around known phosphorylation sites, supporting a potential wider-spread role in phosphorylation-mediated regulation of intrinsically disordered proteins. The identification of such motifs is important for understanding the molecular mechanism of cellular signalling, and is crucial for the development of predictors for the structural effect of phosphorylation; a tool of relevance for understanding disease-promoting mutations that for example interfere with signalling for instance through constitutive active and often cancer-promoting signalling.

The Mutational Landscape of the Oncogenic MZF1 SCAN Domain in Cancer.

SCAN domains in zinc-finger transcription factors are crucial mediators of protein-protein interactions. Up to 240 SCAN-domain encoding genes have been identified throughout the human genome. These include cancer-related genes, such as the myeloid zinc finger 1 (MZF1), an oncogenic transcription factor involved in the progression of many solid cancers. The mechanisms by which SCAN homo- and heterodimers assemble and how they alter the transcriptional activity of zinc-finger transcription factors in cancer and other diseases remain to be investigated. Here, we provide the first description of the conformational ensemble of the MZF1 SCAN domain cross-validated against NMR experimental data, which are probes of structure and dynamics on different timescales. We investigated the protein-protein interaction network of MZF1 and how it is perturbed in different cancer types by the analyses of high-throughput proteomics and RNASeq data. Collectively, we integrated many computational approaches, ranging from simple empirical energy functions to all-atom microsecond molecular dynamics simulations and network analyses to unravel the effects of cancer-related substitutions in relation to MZF1 structure and interactions.

Concurrent elevation of CO2, O3 and temperature severely affects oil quality and quantity in rapeseed.

Plant oil is an essential dietary and bio-energy resource. Despite this, the effects of climate change on plant oil quality remain to be elucidated. The present study is the first to show changes in oil quality and quantity of four rapeseed cultivars in climate scenarios with elevated [CO2], [O3] and temperature (T) combined and as single factors. The combination of environmental factors resembled IPCC's 'business as usual' emission scenario predicted for late this century. Generally, the climate scenarios reduced the average amounts of the six fatty acids (FAs) analysed, though in some treatments single FAs remained unchanged or even increased. Most reduced was the FA essential for human nutrition, C18:3-ω3, which decreased by 39% and 45% in the combined scenarios with elevated [CO2]+T+[O3] and [CO2]+T, respectively. Average oil content decreased 3-17%. When [CO2] and T were elevated concurrently, the seed biomass was reduced by half, doubling the losses in FAs and oil content. This corresponded to a 58% reduction in the oil yield per hectare, and C18:3-ω3 decreased by 77%. Furthermore, the polyunsaturated FAs were significantly decreased. The results indicate undesirable consequences for production and health benefits of rapeseed oil with future climate change. The results also showed strong interactive effects of CO2, T and O3 on oil quality, demonstrating why prediction of climate effects requires experiments with combined factors and should not be based on extrapolation from single factor experiments.