Umeå University's logo

Background: Pancreatic ductal adenocarcinoma (pancreatic cancer) is often detected at late stages resulting in poor overall survival. To improve survival, more patients need to be diagnosed early when curative surgery is feasible. We aimed to identify circulating metabolites that could be used as early pancreatic cancer biomarkers.

Methods: We performed metabolomics by liquid and gas chromatography-mass spectrometry in plasma samples from 82 future pancreatic cancer patients and 82 matched healthy controls within the Northern Sweden Health and Disease Study (NSHDS). Logistic regression was used to assess univariate associations between metabolites and pancreatic cancer risk. Least absolute shrinkage and selection operator (LASSO) logistic regression was used to design a metabolite-based risk score. We used receiver operating characteristic (ROC) analyses to assess the discriminative performance of the metabolite-based risk score.

Results: Among twelve risk-associated metabolites with a nominal P value <0.05, we defined a risk score of three metabolites [indoleacetate, 3-hydroxydecanoate (10:0-OH), and retention index (RI): 2,745.4] using LASSO. A logistic regression model containing these three metabolites, age, sex, body mass index (BMI), smoking status, sample date, fasting status, and carbohydrate antigen 19-9 (CA 19-9) yielded an internal area under curve (AUC) of 0.784 [95% confidence interval (CI): 0.714–0.854] compared to 0.681 (95% CI: 0.597–0.764) for a model without these metabolites (P value =0.007). Seventeen metabolites were significantly associated with pancreatic cancer survival [false discovery rate (FDR) <0.1].

Conclusions: Indoleacetate, 3-hydroxydecanoate (10:0-OH), and RI: 2,745.4 were identified as the top candidate biomarkers for early detection. However, continued efforts are warranted to determine the usefulness of these metabolites as early pancreatic cancer biomarkers.

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive disease with poor survival. Novel biomarkers are urgently needed to improve the outcome through early detection. Here, we aimed to discover novel biomarkers for early PDAC detection using multi-omics profiling in pre-diagnostic plasma samples biobanked after routine health examinations.

A nested case-control study within the Northern Sweden Health and Disease Study was designed. Pre-diagnostic plasma samples from 37 future PDAC patients collected within 2.3 years before diagnosis and 37 matched healthy controls were included. We analyzed metabolites using liquid chromatography mass spectrometry and gas chromatography mass spectrometry, microRNAs by HTG edgeseq, proteins by multiplex proximity extension assays, as well as three clinical biomarkers using milliplex technology. Supervised and unsupervised multi-omics integration were performed as well as univariate analyses for the different omics types and clinical biomarkers. Multiple hypothesis testing was corrected using Benjamini-Hochberg's method and a false discovery rate (FDR) below 0.1 was considered statistically significant.

Carbohydrate antigen (CA) 19-9 was associated with PDAC risk (OR [95 % CI] = 3.09 [1.31–7.29], FDR = 0.03) and increased closer to PDAC diagnosis. Supervised multi-omics models resulted in poor discrimination between future PDAC cases and healthy controls with obtained accuracies between 0.429–0.500. No single metabolite, microRNA, or protein was differentially altered (FDR < 0.1) between future PDAC cases and healthy controls.

CA 19-9 levels increase up to two years prior to PDAC diagnosis but extensive multi-omics analysis including metabolomics, microRNAomics and proteomics in this cohort did not identify novel early biomarkers for PDAC.

Purpose: Peritoneal infection, due to anastomotic leakage, after resection for colorectal cancer have been shown to associate with increased cancer recurrence and mortality, as well as cardiovascsular morbidity. Alterations in circulating protein levels could help shed light on the underlying mechanisms, prompting this exploratory study of 64 patients operated for colorectal cancer with anastomosis. Methods: Thirty-two cases who suffered a postoperative peritoneal infection were matched with 32 controls who had a complication-free postoperative stay. Proteins in serum samples at their first postoperative visit and at one year after surgery were analysed using proximity extension assays and enzyme-linked immunosorbent assays. Multivariate projection methods, adjusted for multiple testing, were used to compare levels between groups, and enrichment and network analyses were performed. Results: Seventy-seven proteins, out of 270 tested, were differentially expressed at a median sampling time of 41 days postoperatively. These proteins were all normalised one year after surgery. Many of the differentially expressed top hub proteins have known involvement in cancer progression, survival, invasiveness and metastasis. Over-represented pathways were related to cardiomyopathy, cell-adhesion, extracellular matrix, phosphatidylinositol-3-kinase/Akt (PI3K-Akt) and transforming growth factor beta (TGF-β) signaling. Conclusion: These affected proteins and pathways could provide clues as to why patients with peritoneal infection might suffer increased cancer recurrence, mortality and cardiovascular morbidity.

Objectives: Pancreatic ductal adenocarcinoma (PDAC) is characterized by poor prognosis and lack of biomarkers. A rich desmoplastic tumor stroma is considered a hallmark of PDAC and previous studies have indicated upregulated expression of collagen VI (COL6) in PDAC. COL6 is shown to associate with prognosis in many cancers but has been less extensively studied in PDAC.

Materials and Methods: The expression of COL6 was analyzed by immunohistochemistry in tissue microarrays containing resected tumor tissue samples from PDAC patients (n = 164). Significance of COL6 was estimated with Kaplan-Meier survival estimates and multivariable Cox regression analysis. COL6 protein and mRNA expression patterns were further investigated in publicly available datasets.

Results: There were no statistically significant (P < 0.05) differences in survival when comparing high and low protein expression of any of the analyzed COL6 α-chains (α1(VI): hazard ratio [HR] 0.90, 95% confidence interval [CI] 0.64–1.28; α2(VI): HR 1.28, 95% CI 0.86–1.89; α3(VI): HR 0.91, 95% CI 0.64–1.29). Similar results were obtained when assessing public data from the Cancer Proteome Atlas, Clinical Proteomic Tumor Analysis Consortium, and The Cancer Genome Atlas.

Conclusions: In contrast with previous studies and some other cancers, we did not find any association of COL6 tissue expression and PDAC survival.

Neuroblastoma is a childhood tumour that is responsible for approximately 15% of all childhood cancer deaths. Neuroblastoma tumours with amplification of the oncogene MYCN are aggressive, however, another aggressive subgroup without MYCN amplification also exists; rather, they have a deleted region at chromosome arm 11q. Twenty-six miRNAs are located within the breakpoint region of chromosome 11q and have been checked for a possible involvement in development of neuroblastoma due to the genomic alteration. Target genes of these miRNAs are involved in pathways associated with cancer, including proliferation, apoptosis and DNA repair. We could show that miR-548l found within the 11q region is downregulated in neuroblastoma cell lines with 11q deletion or MYCN amplification. In addition, we showed that the restoration of miR-548l level in a neuroblastoma cell line led to a decreased proliferation of these cells as well as a decrease in the percentage of cells in the S phase. We also found that miR-548l overexpression suppressed cell viability and promoted apoptosis, while miR-548l knockdown promoted cell viability and inhibited apoptosis in neuroblastoma cells. Our results indicate that 11q-deleted neuroblastoma and MYCN amplified neuroblastoma coalesce by downregulating miR-548l.

Background: Pancreatic ductal adenocarcinoma (PDAC) is a very aggressive malignancy with a 5-year survival of 10 %. Surgery is the only curative treatment. Unfortunately, few patients are eligible for surgery due to late detection. Thus, we need ways to detect the disease at an earlier stage and for that good screening biomarkers could be used. Previous studies have analyzed circulating analytes in prospective studies to identify early PDAC signals. One such class is microRNAs (miRNAs). MicroRNAs are non-coding RNAs of around 22 nucleotides that act as post- transcriptional regulators by interaction with messenger RNAs (mRNAs). The function of a miRNA can be elucidated by target prediction, to identify its potential targets, followed by enrichment analysis of the predicted targets. Challenges with this approach includes a lot of false positives being generated and that miRNAs can perform their role in a tissue- or disease-specific manner. Other classes of analytes that have previously been studied in prospective PDAC cohorts are metabolites and proteins. 

Aims: This thesis has three aims. First, to build a miRNA functional analysis pipeline with correlation support between miRNA and its predicted target genes. Second, to identify potential circulating biomarkers for early detection of PDAC using multi-omics. Third, to identify potential prognostic metabolites in a prospective PDAC cohort.

Methods: We used publicly available data from the cancer genome atlas-pancreatic adenocarcinoma (TCGA-PAAD) and pre-diagnostic plasma samples from the Northern Sweden Health and Disease Study. We built a pipeline in R including miRNA, mRNA, and protein expression data from TCGA-PAAD for in silico miRNA functional analysis. Pre- diagnostic plasma samples from future PDAC patients as well as matched healthy controls were analyzed using multi- omics. Tissue polypeptide specific antigen (TPS) was analyzed by enzyme linked immunosorbent assay in 267 future PDAC samples and 320 healthy controls. Metabolomics and clinical biomarkers (carbohydrate antigen (CA) 19-9, carcinoembryonic antigen (CEA), and CA 15-3) were profiled in 100 future PDAC samples and 100 healthy controls using liquid chromatography-mass spectrometry (MS), gas chromatography-MS, and multi-plex technology. Of these, a subset of 39 future PDAC patients and 39 healthy controls were profiled for 2083 microRNAs using targeted sequencing and 644 proteins using proximity extension assays. Circulating levels of multi-omics analytes were analyzed using conditional or unconditional logistic regression. Least absolute shrinkage and selection operator (LASSO) in combination with 500 bootstrap iterations identified the most informative variables. The prognostic value of metabolites was assessed using cox regression. Multi-omics factor analysis (MOFA) and data integration analysis for biomarker discovery using latent components (DIABLO) were used for multi-omics integration analyses.

Results: An automated pipeline was built consisting of 1) miRNA target prediction, 2) correlation analyses between miRNA and its targets on mRNA and protein expression levels, and 3) functional enrichment of correlated targets to identify enriched Kyoto encyclopedia of genes and genomes (KEGG) pathways and gene ontology (GO) terms for a specific miRNA. The pipeline was run for all microRNAs (~700) detected in the TCGA-PAAD cohort. These results can be downloaded from a shiny app (https://emmbor.shinyapps.io/mirfa/). TPS was not altered in pre-diagnostic PDAC patients up to 24 years prior to diagnosis, but increased at diagnosis (OR = 1.03, 95 % CI: 1.01-1.05). Internal area under curves of 0.74, 0.80, and 0.88 were achieved for five metabolites, two proteins, and two miRNAs that were selected by LASSO and bootstrap iterations, in combination with CA 19-9. Neither MOFA nor DIABLO separated well between future PDAC cases and healthy controls. 

Conclusions: Our bioinformatics pipeline for in silico functional analysis of microRNAs successfully identifies enriched KEGG pathways and GO terms for miRNA isoforms. The investigated plasma samples are heterogeneous, but among the analyzed variables, we identified five metabolites, two proteins, and two microRNAs with highest potential for early PDAC detection. CA 19-9 levels increased closer to diagnosis. We identified five fatty acids that could be studied in a diagnostic PDAC cohort as prognostic biomarkers. 

Early detection of pancreatic ductal adenocarcinoma (PDAC) is challenging, and late diagnosis partly explains the low 5-year survival. Novel and sensitive biomarkers are needed to enable early PDAC detection and improve patient outcomes. Tissue polypeptide specific antigen (TPS) has been studied as a biomarker in PDAC diagnostics, and it has previously been shown to reflect clinical status better than the ‘golden standard’ biomarker carbohydrate antigen 19-9 (CA 19-9) that is most widely used in the clinical setting. In this cross-sectional case-control study using pre-diagnostic plasma samples, we aim to evaluate the potential of TPS as a biomarker for early PDAC detection. Furthermore, in a subset of individuals with multiple samples available at different time points before diagnosis, a longitudinal analysis was used. We assessed plasma TPS levels using enzyme-linked immunosorbent assay (ELISA) in 267 pre-diagnostic PDAC plasma samples taken up to 18.8 years before clinical PDAC diagnosis and in 320 matched healthy controls. TPS levels were also assessed in 25 samples at PDAC diagnosis. Circulating TPS levels were low both in pre-diagnostic samples of future PDAC patients and in healthy controls, whereas TPS levels at PDAC diagnosis were significantly increased (odds ratio 1.03; 95% confidence interval: 1.01–1.05) in a logistic regression model adjusted for age. In conclusion, TPS levels increase late in PDAC progression and hold no potential as a biomarker for early detection.

Concentration of viruses in water is necessary for detection and quantification of the viruses present, in order to evaluate microbiological barriers in water treatment plants and detect pathogenic viruses during waterborne outbreaks, but there is currently no standardised procedure. In this study, we implemented a previously described fast and simple lanthanum-based protocol for concentration of norovirus genogroup I (GI), genogroup II (GII) and hepatitis A virus (HAV) in drinking and surface water. We compared the results with those of a widely used skimmed milk flocculation method, followed by nucleic acid extraction and RT-qPCR detection. Three seeding levels, with intended concentrations 5 × 103, 5 × 104 and 5 × 105 genome copies/10 L, were added to drinking water or surface water. All seed levels were detected with both flocculation methods. Samples extracted with skimmed milk flocculation had on average 1.82, 1.86 and 1.38 times higher measured concentration of norovirus GI, GII and HAV, respectively, than those extracted with lanthanum flocculation, across all seeding levels and water types tested. Mengovirus was used as a positive process control. Mengovirus recovery was higher for skimmed milk (40.7% in drinking water, 26.0% in surface water) than for lanthanum flocculation (24.4% in drinking water, 9.7% in surface water). Together, this indicates that skimmed milk flocculation provides higher viral recovery than lanthanum flocculation. However, lanthanum-based flocculation can be performed much faster than skimmed milk flocculation (1.5 h versus 16 h flocculation time) and thus could be a good alternative for rapid monitoring of viruses in water.

Background: MicroRNAs (miRNAs) are small RNAs that regulate gene expression at a post-transcriptional level and are emerging as potentially important biomarkers for various disease states, including pancreatic cancer. In silico-based functional analysis of miRNAs usually consists of miRNA target prediction and functional enrichment analysis of miRNA targets. Since miRNA target prediction methods generate a large number of false positive target genes, further validation to narrow down interesting candidate miRNA targets is needed. One commonly used method correlates miRNA and mRNA expression to assess the regulatory effect of a particular miRNA.

The aim of this study was to build a bioinformatics pipeline in R for miRNA functional analysis including correlation analyses between miRNA expression levels and its targets on mRNA and protein expression levels available from the cancer genome atlas (TCGA) and the cancer proteome atlas (TCPA). TCGA-derived expression data of specific mature miRNA isoforms from pancreatic cancer tissue was used.

Results: Fifteen circulating miRNAs with significantly altered expression levels detected in pancreatic cancer patients were queried separately in the pipeline. The pipeline generated predicted miRNA target genes, enriched gene ontology (GO) terms and Kyoto encyclopedia of genes and genomes (KEGG) pathways. Predicted miRNA targets were evaluated by correlation analyses between each miRNA and its predicted targets. MiRNA functional analysis in combination with Kaplan-Meier survival analysis suggest that hsa-miR-885-5p could act as a tumor suppressor and should be validated as a potential prognostic biomarker in pancreatic cancer.

Conclusions: Our miRNA functional analysis (miRFA) pipeline can serve as a valuable tool in biomarker discovery involving mature miRNAs associated with pancreatic cancer and could be developed to cover additional cancer types. Results for all mature miRNAs in TCGA pancreatic adenocarcinoma dataset can be studied and downloaded through a shiny web application at https://emmbor.shinyapps.io/mirfa/.

Microbial source tracking (MST) analysis is essential to identifying and mitigating the fecal pollution of water resources. The signature-based MST method uses a library of sequences to identify contaminants based on operational taxonomic units (OTUs) that are unique to a certain source. However, no clear guidelines for how to incorporate OTU overlap or natural variation in the raw water bacterial community into MST analyses exist. We investigated how the inclusion of bacterial overlap between sources in the library affects source prediction accuracy. To achieve this, large-scale sampling-including feces from seven species, raw sewage, and raw water samples from water treatment plants - was followed by 16S rRNA amplicon sequencing. The MST library was defined using three settings: (i) no raw water communities represented; (ii) raw water communities selected through clustering analysis; and (iii) local water communities collected across consecutive years. The results suggest that incorporating either the local background or representative bacterial composition improves MST analyses, as the results were positively correlated to measured levels of fecal indicator bacteria and the accuracy at which OTUs were assigned to the correct contamination source increased fourfold. Using the proportion of OTUs with high source origin probability, underpinning a contaminating signal, is a solid foundation in a framework for further deciphering and comparing contaminating signals derived in signature-based MST approaches. In conclusion, incorporating background bacterial composition of water in MST can improve mitigation efforts for minimizing the spread of pathogenic and antibiotic resistant bacteria into essential freshwater resources.