Search

The plant breeding industry has shown growing interest in using the genotype data of relevant markers for performing selection of new competitive varieties. The selection usually benefits from large amounts of marker data and it is therefore crucial to dispose of data collection methods that are both cost-effective and reliable. Computational methods such as genotype imputation have been proposed earlier in several plant science studies for addressing the cost challenge. Genotype imputation methods have though been used more frequently and investigated more extensively in human genetics research. The various algorithms that exist have shown lower accuracy at inferring the genotype of genetic variants occurring at low frequency, while these rare variants can have great significance and impact in the genetic studies that underlie selection. In contrast, pooling is a technique that can efficiently identify low-frequency items in a population and it has been successfully used for detecting the samples that carry rare variants in a population. In this study, we propose to combine pooling and imputation, and demonstrate this by simulating a hypothetical microarray for genotyping a population of recombinant inbred lines in a cost-effective and accurate manner, even for rare variants. We show that with an adequate imputation model, it is feasible to accurately predictthe individual genotypes at lower cost than sample-wise genotyping and time-effectively. Moreover, we provide code resources for reproducing the results presented in this study in the form of a containerized workflow., eSSENCE - An eScience Collaboration

Since the first human genome assembly’s release, our knowledge of the genetic architecture of complex traits and diseases has grown steadily. Genome-wide association studies (GWAS) played a major role but are limited to common traits and single-nucleotide polymorphisms (SNPs). Technologies and resources like next-generation sequencing, Mendelian Randomization (MR), long-read sequencing and improved reference genomes enable the investigation of variants inaccessible to GWAS, such as copy number variations (CNVs), rare variants and variants in previously unresolved regions. In project I, we performed a GWAS of estradiol measurements using data from UK Biobank and quantified estradiol’s effect on bone mineral density (BMD) using MR. 14 loci were associated with estradiol levels in males, of which one was also significant in females and an additional female-specific locus. We found a significant effect of estradiol on BMD, confirming previous research of estrogen’s importance for skeletal health. In project II, we used the GWAS results from project I to investigate the effect of endogenous estradiol on breast, endometrial and ovarian cancer using MR. Estradiol was associated with ovarian cancer and nominally associated with estrogen receptor-positive breast cancer, demonstrating the effect of endogenous estrogen on cancer risk. In project III, we quantified the effect of 184,182 CNVs on 438 blood plasma proteins using whole-genome sequencing (WGS) data from a Northern Swedish cohort and validated our findings using long-read sequencing in a subcohort. 15 CNVs were associated with 16 proteins of which four could be validated using long reads and three more were more complex variation. Our findings show the effects of CNVs on the plasma proteome and highlight the application different sequencing technologies for CNV detection. In project IV, we evaluated the use of T2T-CHM13 as reference for the SweGen cohort. Compared to GRCh38, mapping quality improved and we identifie

Prehistoric chewed pitch has proven to be a useful source of ancient DNA, both from humans and their microbiomes. Here we present the metagenomic analysis of three pieces of chewed pitch from Huseby Klev, Sweden, that were dated to 9,890-9,540 before present. The metagenomic profile exposes a Mesolithic oral microbiome that includes opportunistic oral pathogens. We compared the data with healthy and dysbiotic microbiome datasets and we identified increased abundance of periodontitis-associated microbes. In addition, trained machine learning models predicted dysbiosis with 70-80% probability. Moreover, we identified DNA sequences from eukaryotic species such as red fox, hazelnut, red deer and apple. Our results indicate a case of poor oral health during the Scandinavian Mesolithic, and show that pitch pieces have the potential to provide information on material use, diet and oral health.

Background: Transcriptome assembly from RNA-sequencing data in species without a reliable reference genome has to be performed de novo, but studies have shown that de novo methods often have inadequate ability to reconstruct transcript isoforms. We address this issue by constructing an assembly pipeline whose main purpose is to produce a comprehensive set of transcript isoforms. Results: We present the de novo transcript isoform assembler ClusTrast, which takes short read RNA-seq data as input, assembles a primary assembly, clusters a set of guiding contigs, aligns the short reads to the guiding contigs, assembles each clustered set of short reads individually, and merges the primary and clusterwise assemblies into the final assembly. We tested ClusTrast on real datasets from six eukaryotic species, and showed that ClusTrast reconstructed more expressed known isoforms than any of the other tested de novo assemblers, at a moderate reduction in precision. For recall, ClusTrast was on top in the lower end of expression levels (<15% percentile) for all tested datasets, and over the entire range for almost all datasets. Reference transcripts were often (35–69% for the six datasets) reconstructed to at least 95% of their length by ClusTrast, and more than half of reference transcripts (58–81%) were reconstructed with contigs that exhibited polymorphism, measuring on a subset of reliably predicted contigs. ClusTrast recall increased when using a union of assembled transcripts from more than one assembly tool as primary assembly. Conclusion: We suggest that ClusTrast can be a useful tool for studying isoforms in species without a reliable reference genome, in particular when the goal is to produce a comprehensive transcriptome set with polymorphic variants., Not duplicate with DiVA 1810720QC 20240215

Accurately predicting the tool remaining useful life (RUL) is critical for maximizing tool utilization and saving machining costs. Various physical model-based or data-driven prediction methods have been developed and successfully applied in different machining operations. However, many uncertain factors affect tool RUL during the cutting process, making it challenging to create a precise physical model to characterize the degradation of tool performance. The success of the purely data-driven technique depends on the amount and quality of the training samples, it does not consider the physical law of tool wear, and the interpretability of the prediction results is poor. This paper presents a data-model linkage approach for tool RUL prediction based on deep feature fusion and Wiener process to address the above limitations. A convolutional stacked bidirectional long short-term memory network with time-space attention mechanism (CSBLSTM-TSAM) is developed in the data-driven module to fuse the multi-sensor signals collected during the cutting process and then obtain the mapping relationship between signal features and tool wear values. In the physical modeling module, a three-stage tool RUL prediction model based on the nonlinear Wiener process is established by considering the evolution law of different wear stages and multi-layer uncertainty, and the corresponding probability density function is derived. The real-time estimated tool wear of the data-driven module is used as the observed value of the physical model, and the model parameters are dynamically updated by the weight-optimized particle filter (WOPF) algorithm under a Bayesian framework, thereby realizing the data-model linkage tool RUL prediction. Milling experiments demonstrate that the proposed method not only improves RUL prediction accuracy, but also has good generalization ability and robustness for prediction tasks under different working conditions., QC 20240201

Spatial biology is a rapidly growing field that has seen tremendous progress over the last decade. We are now able to measure how the morphology, genome, transcriptome, and proteome of a tissue vary across space. Datasets generated by spatial technologies reflect the complexity of the systems they measure: They are multi-modal, high-dimensional, and layer an intricate web of dependencies between biological compartments at different length scales. To add to this complexity, measurements are often sparse and noisy, obfuscating the underlying biological signal and making the data difficult to interpret. In this thesis, we describe how data from spatial biology experiments can be analyzed with methods from deep learning and generative modeling to accelerate biological discovery. The thesis is divided into two parts. The first part provides an introduction to the fields of deep learning and spatial biology, and how the two can be combined to model spatial biology data. The second part consists of four papers describing methods that we have developed for this purpose. Paper I presents a method for inferring spatial gene expression from hematoxylin and eosin stains. The proposed method offers a data-driven approach to analyzing histopathology images without relying on expert annotations and could be a valuable tool for cancer screening and diagnosis in the clinics. Paper II introduces a method for jointly modeling spatial gene expression with histology images. We show that the method can predict super-resolved gene expression and transcriptionally characterize small-scale anatomical structures. Paper III proposes a method for learning flexible Markov kernels to model continuous and discrete data distributions. We demonstrate the method on various image synthesis tasks, including unconditional image generation and inpainting. Paper IV leverages the techniques introduced in Paper III to integrate data from different spatial biology experiments. The proposed method can be use, Spatial biologi är ett snabbt växande forskningsområde som har sett en hög utvecklingstakt under det senaste decenniet. Vi kan idag mäta hur en vävnads morfologi, genom, transkriptom och proteom varierar i rummet. Dataset skapade av spatiala teknologier återspeglar komplexiteten i de system de mäter: De är multimodala, högdimensionella och är uppbyggda av ett intrikat nätverk av beroenden mellan biologiska strukturer som existerar på olika längdskalor. Som om denna komplexitet inte var nog, är mätningarna ofta både glesa och brusiga, vilket försvårar tolkningen av den underliggande biologiska signalen. I denna avhandling beskriver vi hur data från experiment inom spatial biologi kan analyseras med hjälp av djupinlärning och generativ modellering för att accelerera biologiska upptäckter. Avhandlingen är uppdelad i två delar. Den första delen ger en introduktion till fälten djupinlärning och spatial biologi, och hur dessa kan kombineras för att modellera data inom spatial biologi. Den andra delen består av fyra artiklar som beskriver metoder som vi har utvecklat för detta ändamål. Artikel I presenterar en metod för att skatta spatialt genuttryck från hematoxylin-eosin-färgningar. Den föreslagna metoden erbjuder ett datadrivet tillvägagångssätt för att analysera histopatologi-bilder utan användning av expertannoteringar och kan utgöra ett värdefullt verktyg för cancerscreening och diagnos i kliniken. Artikel II introducerar en metod för sammodellering av spatialt genuttryck och histologibilder. Vi visar att metoden kan användas för att predicera superupplöst genuttryck och transkriptionellt karakterisera småskaliga anatomiska strukturer. Artikel III beskriver en metod för modellering av kontinuerliga och diskreta datafördelningar med flexibla Markovkärnor. Vi demonstrerar metoden på olika bildgenereringsuppgifter, inklusive obetingad datagenerering och inpainting. Artikel IV utnyttjar teknikerna från Artikel III för att integrera data från olika experiment inom spatial, QC 2024-01-09

MOTIVATION: In pathway analysis, we aim to establish a connection between the activity of a particular biological pathway and a difference in phenotype. There are many available methods to perform pathway analysis, many of them rely on an upstream differential expression analysis, and many model the relations between the abundances of the analytes in a pathway as linear relationships. RESULTS: Here, we propose a new method for pathway analysis, MIPath, that relies on information theoretical principles and, therefore, does not model the association between pathway activity and phenotype, resulting in relatively few assumptions. For this, we construct a graph of the data points for each pathway using a nearest-neighbor approach and score the association between the structure of this graph and the phenotype of these same samples using Mutual Information while adjusting for the effects of random chance in each score. The initial nearest neighbor approach evades individual gene-level comparisons, hence making the method scalable and less vulnerable to missing values. These properties make our method particularly useful for single-cell data. We benchmarked our method on several single-cell datasets, comparing it to established and new methods, and found that it produces robust, reproducible, and meaningful scores. AVAILABILITY AND IMPLEMENTATION: Source code is available at https://github.com/statisticalbiotechnology/mipath, or through Python Package Index as "mipathway.", QC 20240130

The plant breeding industry has shown growing interest in using the genotype data of relevant markers for performing selection of new competitive varieties. The selection usually benefits from large amounts of marker data and it is therefore crucial to dispose of data collection methods that are both cost-effective and reliable. Computational methods such as genotype imputation have been proposed earlier in several plant science studies for addressing the cost challenge. Genotype imputation methods have though been used more frequently and investigated more extensively in human genetics research. The various algorithms that exist have shown lower accuracy at inferring the genotype of genetic variants occurring at low frequency, while these rare variants can have great significance and impact in the genetic studies that underlie selection. In contrast, pooling is a technique that can efficiently identify low-frequency items in a population and it has been successfully used for detecting the samples that carry rare variants in a population. In this study, we propose to combine pooling and imputation, and demonstrate this by simulating a hypothetical microarray for genotyping a population of recombinant inbred lines in a cost-effective and accurate manner, even for rare variants. We show that with an adequate imputation model, it is feasible to accurately predictthe individual genotypes at lower cost than sample-wise genotyping and time-effectively. Moreover, we provide code resources for reproducing the results presented in this study in the form of a containerized workflow., eSSENCE - An eScience Collaboration

Prehistoric chewed pitch has proven to be a useful source of ancient DNA, both from humans and their microbiomes. Here we present the metagenomic analysis of three pieces of chewed pitch from Huseby Klev, Sweden, that were dated to 9,890-9,540 before present. The metagenomic profile exposes a Mesolithic oral microbiome that includes opportunistic oral pathogens. We compared the data with healthy and dysbiotic microbiome datasets and we identified increased abundance of periodontitis-associated microbes. In addition, trained machine learning models predicted dysbiosis with 70-80% probability. Moreover, we identified DNA sequences from eukaryotic species such as red fox, hazelnut, red deer and apple. Our results indicate a case of poor oral health during the Scandinavian Mesolithic, and show that pitch pieces have the potential to provide information on material use, diet and oral health.

The plant breeding industry has shown growing interest in using the genotype data of relevant markers for performing selection of new competitive varieties. The selection usually benefits from large amounts of marker data and it is therefore crucial to dispose of data collection methods that are both cost-effective and reliable. Computational methods such as genotype imputation have been proposed earlier in several plant science studies for addressing the cost challenge. Genotype imputation methods have though been used more frequently and investigated more extensively in human genetics research. The various algorithms that exist have shown lower accuracy at inferring the genotype of genetic variants occurring at low frequency, while these rare variants can have great significance and impact in the genetic studies that underlie selection. In contrast, pooling is a technique that can efficiently identify low-frequency items in a population and it has been successfully used for detecting the samples that carry rare variants in a population. In this study, we propose to combine pooling and imputation, and demonstrate this by simulating a hypothetical microarray for genotyping a population of recombinant inbred lines in a cost-effective and accurate manner, even for rare variants. We show that with an adequate imputation model, it is feasible to accurately predictthe individual genotypes at lower cost than sample-wise genotyping and time-effectively. Moreover, we provide code resources for reproducing the results presented in this study in the form of a containerized workflow., eSSENCE - An eScience Collaboration

The information conveyed by genetic markers, such as single nucleotide polymorphisms (SNPs), has been widely used in biomedical research to study human diseases and is increasingly valued in agriculture for genomic selection purposes. Specific markers can be identified as a genetic signature that correlates with certain characteristics in a living organism, e.g. a susceptibility to disease or high-yield traits. Capturing these signatures with sufficient statistical power often requires large volumes of data, with thousands of samples to be analysed and potentially millions of genetic markers to be screened. Relevant effects are particularly delicate to detect when the genetic variations involved occur at low frequencies. The cost of producing such marker genotype data is therefore a critical part of the analysis. Despite recent technological advances, production costs can still be prohibitive on a large scale and genotype imputation strategies have been developed to address this issue. Genotype imputation methods have been extensively studied in human data and, to a lesser extent, in crop and animal species. A recognised weakness of imputation methods is their lower accuracy in predicting the genotypes for rare variants, whereas those can be highly informative in association studies and improve the accuracy of genomic selection. In this respect, pooling strategies can be well suited to complement imputation, as pooling is efficient at capturing the low-frequency items in a population. Pooling also reduces the number of genotyping tests required, making its use in combination with imputation a cost-effective compromise between accurate but expensive high-density genotyping of each sample individually and stand-alone imputation. However, due to the nature of genotype data and the limitations of genotype testing techniques, decoding pooled genotypes into unique data resolutions is challenging. In this work, we study the characteristics of decoded genotype data from po

Since the first human genome assembly’s release, our knowledge of the genetic architecture of complex traits and diseases has grown steadily. Genome-wide association studies (GWAS) played a major role but are limited to common traits and single-nucleotide polymorphisms (SNPs). Technologies and resources like next-generation sequencing, Mendelian Randomization (MR), long-read sequencing and improved reference genomes enable the investigation of variants inaccessible to GWAS, such as copy number variations (CNVs), rare variants and variants in previously unresolved regions. In project I, we performed a GWAS of estradiol measurements using data from UK Biobank and quantified estradiol’s effect on bone mineral density (BMD) using MR. 14 loci were associated with estradiol levels in males, of which one was also significant in females and an additional female-specific locus. We found a significant effect of estradiol on BMD, confirming previous research of estrogen’s importance for skeletal health. In project II, we used the GWAS results from project I to investigate the effect of endogenous estradiol on breast, endometrial and ovarian cancer using MR. Estradiol was associated with ovarian cancer and nominally associated with estrogen receptor-positive breast cancer, demonstrating the effect of endogenous estrogen on cancer risk. In project III, we quantified the effect of 184,182 CNVs on 438 blood plasma proteins using whole-genome sequencing (WGS) data from a Northern Swedish cohort and validated our findings using long-read sequencing in a subcohort. 15 CNVs were associated with 16 proteins of which four could be validated using long reads and three more were more complex variation. Our findings show the effects of CNVs on the plasma proteome and highlight the application different sequencing technologies for CNV detection. In project IV, we evaluated the use of T2T-CHM13 as reference for the SweGen cohort. Compared to GRCh38, mapping quality improved and we identifie

The protein sequence determines how it will fold into its unique three-dimensional structure. Once folded, proteins perform their functions by interacting with other proteins or molecules called ligands within the cell. Experimental determination of protein structure and function is tedious. Computational approaches aim to accurately predict the properties of proteins to complement experimental efforts of understanding biochemical mechanisms within the cell. This thesis introduces computational techniques that predict the structure of protein complexes and identify protein residues involved in interactions with common biomolecules, such as metal ions and nucleic acids, based on sequence information. AlphaFold, a method that predicted protein structure using sequence information with almost experimental accuracy, was a critical breakthrough that shaped the field of protein structure prediction. Subsequently, approaches such as FoldDock adapted the AlphaFold pipeline for dimer complexes. Paper I applies the FoldDock protocol to understand toxin-antitoxin systems. These protein complexes are highly evolutionary conserved, and high-confidence dimer predictions were generated. Paper II applies the FoldDock protocol to study protein-protein interactions in the human proteome. To verify the reliability of machine-learning-based computational methods, they must be tested on independent data different from the data used to train the method. Paper III involves generating and using a homology-reduced independent test set to benchmark the performance of protein complex structure predictors, including the recent AlphaFold release adapted for multi-chain proteins – AlphaFold-Multimer. A confidence score (pDockQ2) was proposed to estimate the quality of the interfaces within multimers. Paper I, Paper II and Paper III are associated with predicting and evaluating protein-protein interactions. Representation learning involves finding effective representations of input data to maxi

Short tandem repeat (STR) mutations are prevalent in colorectal cancer (CRC), especially in tumours with the microsatellite instability (MSI) phenotype. While STR length variations are known to regulate gene expression under physiological conditions, the functional impact of STR mutations in CRC remains unclear. Here, we integrate STR mutation data with clinical information and gene expression data to study the gene regulatory effects of STR mutations in CRC. We confirm that STR mutability in CRC highly depends on the MSI status, repeat unit size, and repeat length. Furthermore, we present a set of 1244 putative expression STRs (eSTRs) for which the STR length is associated with gene expression levels in CRC tumours. The length of 73 eSTRs is associated with expression levels of cancer-related genes, nine of which are CRC-specific genes. We show that linear models describing eSTR-gene expression relationships allow for predictions of gene expression changes in response to eSTR mutations. Moreover, we found an increased mutability of eSTRs in MSI tumours. Our evidence of gene regulatory roles for eSTRs in CRC highlights a mostly overlooked way through which tumours may modulate their phenotypes. Future extensions of these findings could uncover new STR-based targets in the treatment of cancer.

Prehistoric chewed pitch has proven to be a useful source of ancient DNA, both from humans and their microbiomes. Here we present the metagenomic analysis of three pieces of chewed pitch from Huseby Klev, Sweden, that were dated to 9,890-9,540 before present. The metagenomic profile exposes a Mesolithic oral microbiome that includes opportunistic oral pathogens. We compared the data with healthy and dysbiotic microbiome datasets and we identified increased abundance of periodontitis-associated microbes. In addition, trained machine learning models predicted dysbiosis with 70-80% probability. Moreover, we identified DNA sequences from eukaryotic species such as red fox, hazelnut, red deer and apple. Our results indicate a case of poor oral health during the Scandinavian Mesolithic, and show that pitch pieces have the potential to provide information on material use, diet and oral health.

The plant breeding industry has shown growing interest in using the genotype data of relevant markers for performing selection of new competitive varieties. The selection usually benefits from large amounts of marker data and it is therefore crucial to dispose of data collection methods that are both cost-effective and reliable. Computational methods such as genotype imputation have been proposed earlier in several plant science studies for addressing the cost challenge. Genotype imputation methods have though been used more frequently and investigated more extensively in human genetics research. The various algorithms that exist have shown lower accuracy at inferring the genotype of genetic variants occurring at low frequency, while these rare variants can have great significance and impact in the genetic studies that underlie selection. In contrast, pooling is a technique that can efficiently identify low-frequency items in a population and it has been successfully used for detecting the samples that carry rare variants in a population. In this study, we propose to combine pooling and imputation, and demonstrate this by simulating a hypothetical microarray for genotyping a population of recombinant inbred lines in a cost-effective and accurate manner, even for rare variants. We show that with an adequate imputation model, it is feasible to accurately predictthe individual genotypes at lower cost than sample-wise genotyping and time-effectively. Moreover, we provide code resources for reproducing the results presented in this study in the form of a containerized workflow., eSSENCE - An eScience Collaboration

Deep learning has the potential to be used to improve and streamline EEG analysis. At the present, classifiers and supervised learning dominate the field. Supervised learning depends on target labels which most often are created by human experts manually classifying data. A problem with supervised learning is intra- and interrater agreement which in some instances are far from perfect. This can affect the training and make evaluation more difficult. This thesis includes three papers where self-supervised deep neural networks were developed. In self-supervised learning, the input data to the networks themselves contain structures that are used as targets for the training and no labeling is necessary. In paper I, deep neural networks were trained to increase the number of-, or to recreate missing EEG-channels. The performance was at least on the same level as that of spherical interpolation, but unlike in the case of interpolation, missing data does not have to be identified manually first. Papers II and III involved developing deep neural networks for clustering analysis. The networks produced two-dimensional representations of EEG data and the training strategy was based on the principle of t-distributed stochastic neighbor embedding (t-SNE). In paper II, comparisons were made to parametric t-SNE and EEG-features obtained from time-frequency methods. The deep neural networks produced more distinct clustering when tested on data annotated for epileptiform discharges, seizure activity, or sleep-wakefulness. In paper III, the newly developed method was used to compare annotations of epileptiform discharges. Two experts performed independent annotations and classifiers were trained on these, using supervised learning, which in turn produced new annotations. The agreement when comparing two sets of annotations was not larger between the two experts than between an expert and a classifier. The analysis showed that differences in the annotations by the experts infl, Funding: The work was funded by Saad Nagi (grants RÖ-974228, RÖ-962769, and RÖ-941377), Magnus Thordstein (grant RÖ-986017), and Håkan Olaus-son (grants LIO-936017 and RÖ-941359).

The brain is the most metabolically active organ in the human body and therefore rely on a continuous supply of oxygen and glucose. Neuronal stimulation in specific regions of the leads to the firing of action potentials, a process facilitated by voltage-gated ion channels in the neurons’ cell membranes. This activation of the ion channels significantly elevates the brain’s metabolic energy demand, compelling neurons to ramp up their metabolic activity in response. Concurrently, this neuronal activation also initiates a signalling cascade that induces vasodilation and increases blood flow, thereby ensuring that regions with elevated neural activity are adequately supplied with oxygen and nutrients. This dynamic interplay between neuronal activity and cerebral blood flow (CBF) regulation constitutes the neurovascular coupling (NVC). The NVC is a cornerstone in interpreting functional Magnetic Resonance Imaging (fMRI) Blood Oxygen Level-Dependent (BOLD) responses. The BOLD response is an indirect, non-invasive, and highly sensitive indicator of neuronal activity, reflecting changes in blood oxygenation and flow associated with the neuronal and metabolic activity in the brain. By examining these responses, we can gain insights into the complex interactions between neuronal activity, energy metabolism, and CBF. Additionally, techniques such as 13C Metabolic Flux Analysis (13C MFA) makes it possible to gain further insight into the cerebral metabolism. This method enables a detailed examination of metabolic pathways and fluxes by tracking the incorporation of 13C-labelled substrates into various metabolites. By using 13C MFA, researchers can quantify the flow of substrates through metabolic networks, offering a deeper understanding of how cell such as neurons adapt their metabolism during different functional states and conditions. Central to exploring these multifaceted aspects of cerebral metabolism is the use of mathematical modelling and systems biology. These discip

One of the ways that a vascular network grows is through the process of angiogenesis, wherebya new blood vessel forms as a branch from an existing vessel towards an area which isstimulating vascular growth. Due to the demands for nutrients and waste transport, growingtumour cells will access the surrounding vascular network by inducing angiogenesis. Once thetumour is connected with the vascular system it can grow further and colonize distant organs.Given the critical nature of this step in tumour development, there is a demand for mathematicaland computational models to provide an understanding of the process for treatment in predictivemedicine. These models allow us to generate vascular networks that demonstrate similarbehaviour to that of the observed networks; however, there is a lack of quantifiable measures ofsimilarity between generated networks, or, of a generated and real network. Furthermore, thereis not an established way to determine which measures hold the most relevance todistinguishing similarity. To construct such a measure we transform our generated vascularnetworks into an abstract graph representation which allows exploration of the plethora of graphcentralities. We propose to determine the relevance of a centrality by finding one that acts as asynthetic likelihood function for estimating the model's parameters with minimal error.Evaluating the relevance of many centralities, it is then possible to suggest which centralitiesshould be used to quantitatively determine similarity. This allows for a way to measure howrealistic a model's growth is, and if given sufficient data, to distinguish between regular andtumour-induced angiogenesis and use it within cancer screening.

Prehistoric chewed pitch has proven to be a useful source of ancient DNA, both from humans and their microbiomes. Here we present the metagenomic analysis of three pieces of chewed pitch from Huseby Klev, Sweden, that were dated to 9,890-9,540 before present. The metagenomic profile exposes a Mesolithic oral microbiome that includes opportunistic oral pathogens. We compared the data with healthy and dysbiotic microbiome datasets and we identified increased abundance of periodontitis-associated microbes. In addition, trained machine learning models predicted dysbiosis with 70-80% probability. Moreover, we identified DNA sequences from eukaryotic species such as red fox, hazelnut, red deer and apple. Our results indicate a case of poor oral health during the Scandinavian Mesolithic, and show that pitch pieces have the potential to provide information on material use, diet and oral health.

The information conveyed by genetic markers, such as single nucleotide polymorphisms (SNPs), has been widely used in biomedical research to study human diseases and is increasingly valued in agriculture for genomic selection purposes. Specific markers can be identified as a genetic signature that correlates with certain characteristics in a living organism, e.g. a susceptibility to disease or high-yield traits. Capturing these signatures with sufficient statistical power often requires large volumes of data, with thousands of samples to be analysed and potentially millions of genetic markers to be screened. Relevant effects are particularly delicate to detect when the genetic variations involved occur at low frequencies. The cost of producing such marker genotype data is therefore a critical part of the analysis. Despite recent technological advances, production costs can still be prohibitive on a large scale and genotype imputation strategies have been developed to address this issue. Genotype imputation methods have been extensively studied in human data and, to a lesser extent, in crop and animal species. A recognised weakness of imputation methods is their lower accuracy in predicting the genotypes for rare variants, whereas those can be highly informative in association studies and improve the accuracy of genomic selection. In this respect, pooling strategies can be well suited to complement imputation, as pooling is efficient at capturing the low-frequency items in a population. Pooling also reduces the number of genotyping tests required, making its use in combination with imputation a cost-effective compromise between accurate but expensive high-density genotyping of each sample individually and stand-alone imputation. However, due to the nature of genotype data and the limitations of genotype testing techniques, decoding pooled genotypes into unique data resolutions is challenging. In this work, we study the characteristics of decoded genotype data from po

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., CC BY 4.0Corresponding author at: University Hospital of Umeå, 901 85 Umeå, Sweden. E-mail address: emmy.borgmastars@umu.se (E. Borgmästars)The authors thank Hanna Nyström, and Daniel Öhlund at Umeå University for valuable assistance in data collection. We thank Xiaoshuang Feng at International Agency for Research of Cancer, Lyon, France for guidance in statistical analyses. The authors would also like to thank Swedish Metabolomics Centre, Umeå, Sweden (www.swedishmetabolomicscentre.se) and Biobanken Norr. Funding: This study was funded by Umeå University, the Swedish Research Council [2016-02990, 2019-01690], the Swedish Cancer Society [CAN 2016/643, 19 0273], Region Västerbotten [RV-583411, RV-549731, RV-583411, RV-841551, RV 967602], Finska Läkaresällskapet, Medicinska Understödsföreningen Liv och Hälsa, the Sjöberg Foundation, The JC Kempe Memorial Foundation Scholarship Fund, The Royal Swedish Academy of Sciences (PE Lindahl Foundation, LM2021-0010 and LM2023-0012), The Swedish Society of Medicine (SLS-960379), Cancerforskningsfonden i Norrland (LP 23-2337), Bengt Ihre foundation (SLS-960529 and SLS-986656), and Bengt Ihre Research Fellowship Grant. The sponsors had no role in the study design, collection, analysis and interpretation of data, writing of the report, or in the decision to submit the article for publication.

The influence of natural image transformations on receptive field responses is crucial for modelling visual operations in computer vision and biological vision. In this regard, covariance properties with respect to geometric image transformations in the earliest layers of the visual hierarchy are essential for expressing robust image operations, and for formulating invariant visual operations at higher levels. This paper defines and proves a set of joint covariance properties for spatio-temporal receptive fields in terms of spatio-temporal derivative operators applied to spatio-temporally smoothed image data under compositions of spatial scaling transformations, spatial affine transformations, Galilean transformations and temporal scaling transformations. Specifically, the derived relations show how the parameters of the receptive fields need to be transformed, in order to match the output from spatio-temporal receptive fields under composed spatio-temporal image transformations. For this purpose, we also fundamentally extend the notion of scale-normalized derivatives to affine-normalized derivatives, that are computed based on spatial smoothing with affine Gaussian kernels, and analyze the covariance properties of the resulting affine-normalized derivatives for the affine group as well as for important subgroups thereof. We conclude with a geometric analysis, showing how the derived joint covariance properties make it possible to relate or match spatio-temporal receptive field responses, when observing, possibly moving, local surface patches from different views, under locally linearized perspective or projective transformations, as well as when observing different instances of spatio-temporal events, that may occur either faster or slower between different views of similar spatio-temporal events. We do furthermore describe how the parameters in the studied composed spatio-temporal image transformation models directly relate to geometric entities in the image forma, QC 20240709

This paper presents an in-depth theoretical analysis of the orientation selectivity properties of simple cells and complex cells, that can be well modelled by the generalized Gaussian derivative model for visual receptive fields, with the purely spatial component of the receptive fields determined by oriented affine Gaussian derivatives for different orders of spatial differentiation. A detailed mathematical analysis is presented for the three different cases of either: (i) purely spatial receptive fields, (ii) space-time separable spatio-temporal receptive fields and (iii) velocity-adapted spatio-temporal receptive fields. Closed-form theoretical expressions for the orientation selectivity curves for idealized models of simple and complex cells are derived for all these main cases, and it is shown that the orientation selectivity of the receptive fields becomes more narrow, as a scale parameter ratio $\kappa$, defined as the ratio between the scale parameters in the directions perpendicular to vs. parallel with the preferred orientation of the receptive field, increases. It is also shown that the orientation selectivity becomes more narrow with increasing order of spatial differentiation in the underlying affine Gaussian derivative operators over the spatial domain. Additionally, we also derive closed-form expressions for the resultant and the bandwidth descriptors of the orientation selectivity curves, which have previously been used as compact descriptors of the orientation selectivity properties for biological neurons. These results together show that the properties of the affine Gaussian derivative model for visual receptive fields can be analyzed in closed form, which can be highly useful when to relate the results from biological experiments to computational models of the functional properties of simple cells and complex cells in the primary visual cortex. For comparison, we also present a corresponding theoretical orientation selectivity analysis for purely, QC 20240410

This paper presents the results of combining (i) theoretical analysis regarding connections between the orientation selectivity and the elongation of receptive fields for the affine Gaussian derivative model with (ii) biological measurements of orientation selectivity in the primary visual cortex to investigate if (iii) the receptive fields can be regarded as spanning a variability in the degree of elongation. From an in-depth theoretical analysis of idealized models for the receptive fields of simple and complex cells in the primary visual cortex, we established that the orientation selectivity becomes more narrow with increasing elongation of the receptivefields. Combined with previously established biological results, concerning broad vs. sharp orientation tuning of visual neurons in the primary visual cortex, as well as previous experimental results concerning distributions of the resultant of the orientation selectivity curves for simple and complex cells, we show that these results are consistent with the receptive fields spanning a variability over the degree of elongation of the receptive fields. We also show that our principled theoretical model for visual receptive fields leads to qualitatively similar types of deviations from a uniform histogram of the resultant descriptor of the orientation selectivity curves for simple cells, as can be observed in the results from biological experiments. To firmly determine if the underlying working hypothesis, regarding the receptive fields spanning a variability in the degree of elongation, would truly hold for the receptive fields in the primary visual cortex of higher mammals, we formulate a set of testable predictions, that can be used for investigate this property experimentally, and, if applicable, then also characterize if such a variability would, in a structured way, be related to the pinwheel structure in the visual cortex., QC 20240410

The influence of natural image transformations on receptive field responses is crucial for modelling visual operations in computer vision and biological vision. In this regard, covariance properties with respect to geometric image transformations in the earliest layers of the visual hierarchy are essential for expressing robust image operations, and for formulating invariant visual operations at higher levels. This paper defines and proves a set of joint covariance properties under compositions of spatial scaling transformations, spatial affine transformations, Galilean transformations and temporal scaling transformations, which make it possible to characterize how different types of image transformations interact with each other and the associated spatio-temporal receptive field responses. In this regard, we also extend the notion of scale-normalized derivatives to affine-normalized derivatives, to be able to obtain true affine-covariant properties of spatial derivatives, that are computed based on spatial smoothing with affine Gaussian kernels. The derived relations show how the parameters of the receptive fields need to be transformed, in order to match the output from spatio-temporal receptive fields under composed spatio-temporal image transformations. As a side effect, the presented proof for the joint covariance property over the integrated combination of the different geometric image transformations also provides specific proofs for the individual transformation properties, which have not previously been fully reported in the literature. We conclude with a geometric analysis, showing how the derived joint covariance properties make it possible to relate or match spatio-temporal receptive field responses, when observing, possibly moving, local surface patches from different views, under locally linearized perspective or projective transformations, as well as when observing different instances of spatio-temporal events that may occur either faster or slower be, QC 20240429, Covariant and invariant deep networks

The COVID-19 pandemic highlighted the need for a rapid, convenient, and scalable diagnostic method for detecting a novel pathogen amidst a global pandemic. While command-line interface tools offer automation for SARS-CoV-2 Oxford Nanopore Technology sequencing data analysis, they are inapplicable to users with limited programming skills. A solution is to establish such automated workflows within a graphical user interface software. We developed two workflows in the software Geneious Prime 2022.1.1, adapted for data obtained from the Midnight and Artic's nCoV-2019 sequencing protocols. Both workflows perform trimming, read mapping, consensus generation, and annotation on SARS-CoV-2 Nanopore sequencing data. Additionally, one workflow includes phylogenetic assignment using the bioinformatic tools pangolin and Nextclade as plugins. The basic workflow was validated in 2020, adhering to the requirements of the European Centre for Disease Prevention and Control for SARS-CoV-2 sequencing and analysis. The enhanced workflow, providing phylogenetic assignment, underwent validation at Uppsala University Hospital by analysing 96 clinical samples. It provided accurate diagnoses matching the original results of the basic workflow while also reducing manual clicks and analysis time. These bioinformatic workflows streamline SARS-CoV-2 Nanopore data analysis in Geneious Prime, saving time and manual work for operators lacking programming knowledge.

Transmission electron microscopy (TEM) is an imaging technique used to visualize and analyze nano-sized structures and objects such as virus particles. Light microscopy can be used to diagnose diseases or characterize e.g. blood cells. Since samples under microscopes exhibit certain symmetries, such as global rotation invariance, equivariant neural networks are presumed to be useful. In this study, a baseline convolutional neural network is constructed in the form of the commonly used VGG16 classifier. Thereafter, it is modified to be equivariant to the p4 symmetry group of rotations of multiples of 90 degrees using group convolutions. This yields a number of benefits on a TEM virus dataset, including higher top validation set accuracy by on average 7.6% and faster convergence during training by on average 23.1% of that of the baseline. Similarly, when training and testing on images of blood cells, the convergence time for the equivariant neural network is 7.9% of that of the baseline. From this it is concluded that augmentation strategies for rotation can be skipped. Furthermore, when modelling the accuracy versus amount of TEM virus training data with a power law, the equivariant network has a slope of - 0.43 compared to - 0.26 of the baseline. Thus the equivariant network learns faster than the baseline when more training data is added. This study extends previous research on equivariant neural networks applied to images which exhibit symmetries to isometric transformations.

Artificial intelligence (AI) and high-content imaging (HCI) are contributing to advancements in drug discovery, propelled by the recent progress in deep neural networks. This review highlights AI's role in analysis of HCI data from fixed and livecell imaging, enabling novel label-free and multi-channel fluorescent screening methods, and improving compound profiling. HCI experiments are rapid and cost-effective, facilitating large data set accumulation for AI model training. However, the success of AI in drug discovery also depends on highquality data, reproducible experiments, and robust validation to ensure model performance. Despite challenges like the need for annotated compounds and managing vast image data, AI's potential in phenotypic screening and drug profiling is significant. Future improvements in AI, including increased interpretability and integration of multiple modalities, are expected to solidify AI and HCI's role in drug discovery.

The influence of natural image transformations on receptive field responses is crucial for modelling visual operations in computer vision and biological vision. In this regard, covariance properties with respect to geometric image transformations in the earliest layers of the visual hierarchy are essential for expressing robust image operations, and for formulating invariant visual operations at higher levels. This paper defines and proves a set of joint covariance properties for spatio-temporal receptive fields in terms of spatio-temporal derivative operators applied to spatio-temporally smoothed image data under compositions of spatial scaling transformations, spatial affine transformations, Galilean transformations and temporal scaling transformations. Specifically, the derived relations show how the parameters of the receptive fields need to be transformed, in order to match the output from spatio-temporal receptive fields under composed spatio-temporal image transformations. For this purpose, we also fundamentally extend the notion of scale-normalized derivatives to affine-normalized derivatives, that are computed based on spatial smoothing with affine Gaussian kernels, and analyze the covariance properties of the resulting affine-normalized derivatives for the affine group as well as for important subgroups thereof. We conclude with a geometric analysis, showing how the derived joint covariance properties make it possible to relate or match spatio-temporal receptive field responses, when observing, possibly moving, local surface patches from different views, under locally linearized perspective or projective transformations, as well as when observing different instances of spatio-temporal events, that may occur either faster or slower between different views of similar spatio-temporal events. We do furthermore describe how the parameters in the studied composed spatio-temporal image transformation models directly relate to geometric entities in the image forma, QC 20240709

There has been a long-standing bottleneck in the quantitative analysis of the frequencies of homoblock polyads beyond triads using 1H and 13C NMR for linear polysaccharides, primarily because monosaccharides within a long homoblock share similar chemical environments due to identical neighboring units, resulting in indistinct NMR peaks. In this study, through rigorous mathematical induction, inequality relations were established that enabled the calculation of frequency ranges of homoblock polyads from historically reported NMR-derived frequency values of diads and/or triads of alginates, chitosans, homogalacturonans, and galactomannans. The calculated homoblock frequency ranges were then applied to evaluate three chain growth statistical models, including the Bernoulli chain, first-order Markov chain, and second-order Markov chain, for predicting homoblock frequencies in these polysaccharides. Furthermore, based on the mathematically derived inequality relations, a novel 2D array was constructed, enabling the graphical visualization of homoblock features in polysaccharides. It was demonstrated, as a proof of concept, that the novel 2D array, along with a 1D code generated from it, could serve as an effective feature engineering tool for polymer classification using machine learning algorithms., QC 20240618

In this thesis we studied synaptic plasticity and neuronal computation in single striatal projection neurons (SPNs), which have a major role in goal-directed learning. Goal-directed or reward learning means to learn, based on sensory information from the body and the environment, to select actions out of all the behavioral repertoire that lead to obtaining a goal or reward (such as food or water). In mammals, all the behavioral motor repertoire is under constant, tonic inhibition, and the direct-pathway SPNs (dSPNs) select (disinhibit) goal-obtaining actions. The learning process is guided by the neuromodulator dopamine which signals the positive or negative out-come of an action. The synapses from cortical neurons on to the dSPNs, called corticostriatal synapses, are responsive to dopamine signals, and can strengthen and weaken based on the (positive or negative) action outcome. This promotes or discourages future actions in the same or similar sensory context. Within a collaborative computational modeling effort, we studied the biochemical circuitry in the corticostriatal synapses with multiscale modeling and simulations. This circuitry in the corticostriatal synapses responds to neuromodulatory signals and controls the expression of synaptic plasticity. Multiscale modeling and simulations enable studying a system at multiple temporal and spatial scales, and integrating the results across the different scales. Based on molecular dynamics simulations of the enzyme which transduces extracellular neuromodulatory signals into an intracellular second messenger molecule, and Brownian dynamics simulations of regulator molecules binding to the enzyme, we constructed a kinetic model ofthe enzyme-based signal transduction network. The kinetic model showed that two co-occuring neuromodulatory signals, a dopamine peak and an acetylcholine pause, are required to produce the second messenger and thus enable strengthening of corticostriatal synapses onto dSPNs, and that only the, I denna avhandling studerades synaptisk plasticitet samt förmågan till lokalberäkning i de striatala projektionsneuronens (SPNs) dendriter. Detta har direkt relevans för förståelsen av s.k. målstyrd inlärning. Målstyrd inlärning och belöningsinlärning innebär att man, baserat på sensorisk informationfrån kroppen och omgivningen, lär sig selektera hur man skall agera/handla så att man uppnår ett mål eller erhåller en belöning (t.ex. mat eller vatten). Hos alla däggdjur är motoriska/exekutiva centra i hjärnstammen och thalamus under konstant tonisk inhibition via basala ganglierna, men när manaktiverar SPNs i den s.k. direkta vägen genom basala ganglierna (dSPN) så disinhiberas de motoriska centra som behövs för att initiera specifika mål-styrda beteenden. Inlärningsprocessen för detta guidas av bl.a. dopamin, en neuromodulator som signalerar huruvida ett beteende ger positivt eller negativt resultat. Synapserna från kortikala projektionsneuron till dSPNs, s.k. kortikostriatala synapser, påverkas av dopamin och de kan förstärkas eller försvagas baserat på om dopaminsignalen signalerar ett positivt eller negativt utfall. Detta antingen främjar eller motverkar att man väljer samma beteende/handlingar under en liknande situation i framtiden. Via samarbete med andra beräkningsbiologigrupper kunde vi studera de biokemiska signaleringsnätverken i de kortikostriatala synapserna med hjälp av multiskal modeller och simuleringar. Synapsens intracellulära biokemiska signaleringsnätverk kontrollerar synaptisk plasticitet och påverkas av neuro-modulering. Genom att använda multiskalsimuleringar kunde vi studera systemet över multipla temporala och spatiala skalor, och integrera resultaten över de olika skalorna. Baserat på molekylärdynamiska simuleringar av det enzym som överför de extracellulära neuromodulatoriska signalerna till en intracellulär s.k. 2nd messenger molekyl, och Brownianska dynamiksimuleringar av de regulatoriska molekyler som binder till enzymet, kunde vi konstru, QC 20240904

DNA binding proteins in the form of transcription factors are crucial for gene expression. These enzymes recognize their target gene by a regulatory element on the DNA molecule. The unique combination of nucleotides in regulatory sequences gives rise to very specific interactions with the target protein. Due to this sequence specificity, mutations on regulatory elements can have devastating consequences for an organism. The mutated sequence can lead to decreased affinity of the target protein, and thus decreased gene expression, a cause of several genetic diseases. Not all mutations affect the protein-DNA binding in the same way. Depending on the type of mutation and its location, the severity of its effects can differ greatly. Uncovering this relationship between mutation and sequence affinity can lead to a better understanding of gene expression and genetic diseases. In this project, two Deep Learning models were trained on regulatory DNA sequences with the goal of predicting sequence affinity. One was a convolutional neural network and the other a transformer neural network. The goal was to create an accurate Deep Learning model that can predict binding rates and sequence affinity using only a DNA strand as input. Using a dataset of sequences with measured values for binding rates and affinities, the models achieved promising results. They were able to reliably predict these values with high accuracy. When comparing the Deep Learning models performance with the one of a simpler linear model, they both showed significant improvements. This study showed that neural networks are well suited for this task, achieving accurate results and that a transformer neural network is the superior model when solving this problem.

In this thesis we studied synaptic plasticity and neuronal computation in single striatal projection neurons (SPNs), which have a major role in goal-directed learning. Goal-directed or reward learning means to learn, based on sensory information from the body and the environment, to select actions out of all the behavioral repertoire that lead to obtaining a goal or reward (such as food or water). In mammals, all the behavioral motor repertoire is under constant, tonic inhibition, and the direct-pathway SPNs (dSPNs) select (disinhibit) goal-obtaining actions. The learning process is guided by the neuromodulator dopamine which signals the positive or negative out-come of an action. The synapses from cortical neurons on to the dSPNs, called corticostriatal synapses, are responsive to dopamine signals, and can strengthen and weaken based on the (positive or negative) action outcome. This promotes or discourages future actions in the same or similar sensory context. Within a collaborative computational modeling effort, we studied the biochemical circuitry in the corticostriatal synapses with multiscale modeling and simulations. This circuitry in the corticostriatal synapses responds to neuromodulatory signals and controls the expression of synaptic plasticity. Multiscale modeling and simulations enable studying a system at multiple temporal and spatial scales, and integrating the results across the different scales. Based on molecular dynamics simulations of the enzyme which transduces extracellular neuromodulatory signals into an intracellular second messenger molecule, and Brownian dynamics simulations of regulator molecules binding to the enzyme, we constructed a kinetic model ofthe enzyme-based signal transduction network. The kinetic model showed that two co-occuring neuromodulatory signals, a dopamine peak and an acetylcholine pause, are required to produce the second messenger and thus enable strengthening of corticostriatal synapses onto dSPNs, and that only the, I denna avhandling studerades synaptisk plasticitet samt förmågan till lokalberäkning i de striatala projektionsneuronens (SPNs) dendriter. Detta har direkt relevans för förståelsen av s.k. målstyrd inlärning. Målstyrd inlärning och belöningsinlärning innebär att man, baserat på sensorisk informationfrån kroppen och omgivningen, lär sig selektera hur man skall agera/handla så att man uppnår ett mål eller erhåller en belöning (t.ex. mat eller vatten). Hos alla däggdjur är motoriska/exekutiva centra i hjärnstammen och thalamus under konstant tonisk inhibition via basala ganglierna, men när manaktiverar SPNs i den s.k. direkta vägen genom basala ganglierna (dSPN) så disinhiberas de motoriska centra som behövs för att initiera specifika mål-styrda beteenden. Inlärningsprocessen för detta guidas av bl.a. dopamin, en neuromodulator som signalerar huruvida ett beteende ger positivt eller negativt resultat. Synapserna från kortikala projektionsneuron till dSPNs, s.k. kortikostriatala synapser, påverkas av dopamin och de kan förstärkas eller försvagas baserat på om dopaminsignalen signalerar ett positivt eller negativt utfall. Detta antingen främjar eller motverkar att man väljer samma beteende/handlingar under en liknande situation i framtiden. Via samarbete med andra beräkningsbiologigrupper kunde vi studera de biokemiska signaleringsnätverken i de kortikostriatala synapserna med hjälp av multiskal modeller och simuleringar. Synapsens intracellulära biokemiska signaleringsnätverk kontrollerar synaptisk plasticitet och påverkas av neuro-modulering. Genom att använda multiskalsimuleringar kunde vi studera systemet över multipla temporala och spatiala skalor, och integrera resultaten över de olika skalorna. Baserat på molekylärdynamiska simuleringar av det enzym som överför de extracellulära neuromodulatoriska signalerna till en intracellulär s.k. 2nd messenger molekyl, och Brownianska dynamiksimuleringar av de regulatoriska molekyler som binder till enzymet, kunde vi konstru, QC 20240904

Background: Biomedical data are usually collections of longitudinal data assessed at certain points in time. Clinical observations assess the presences and severity of symptoms, which are the basis for the description and modeling of disease progression. Deciphering potential underlying unknowns from the distinct observation would substantially improve the understanding of pathological cascades. Hidden Markov Models (HMMs) have been successfully applied to the processing of possibly noisy continuous signals. We apply ensembles of HMMs to categorically distributed multivariate time series data, leaving space for expert domain knowledge in the prediction process. Methods: We use an ensemble of HMMs to predict the loss of free walking ability as one major clinical deterioration in the most common autosomal dominantly inherited ataxia disorder worldwide. Results: We present a prediction pipeline that processes data paired with a configuration file, enabling us to train, validate and query an ensemble of HMMs. In particular, we provide a theoretical and practical framework for multivariate time-series inference based on HMMs that includes constructing multiple HMMs, each to predict a particular observable variable. Our analysis is conducted on pseudo-data, but also on biomedical data based on Spinocerebellar ataxia type 3 disease. Conclusions: We find that the model shows promising results for the data we tested. The strength of this approach is that HMMs are well understood, probabilistic and interpretable models, setting it apart from most Deep Learning approaches. We publish all code and evaluation pseudo-data in an open-source repository.

In this thesis we studied synaptic plasticity and neuronal computation in single striatal projection neurons (SPNs), which have a major role in goal-directed learning. Goal-directed or reward learning means to learn, based on sensory information from the body and the environment, to select actions out of all the behavioral repertoire that lead to obtaining a goal or reward (such as food or water). In mammals, all the behavioral motor repertoire is under constant, tonic inhibition, and the direct-pathway SPNs (dSPNs) select (disinhibit) goal-obtaining actions. The learning process is guided by the neuromodulator dopamine which signals the positive or negative out-come of an action. The synapses from cortical neurons on to the dSPNs, called corticostriatal synapses, are responsive to dopamine signals, and can strengthen and weaken based on the (positive or negative) action outcome. This promotes or discourages future actions in the same or similar sensory context. Within a collaborative computational modeling effort, we studied the biochemical circuitry in the corticostriatal synapses with multiscale modeling and simulations. This circuitry in the corticostriatal synapses responds to neuromodulatory signals and controls the expression of synaptic plasticity. Multiscale modeling and simulations enable studying a system at multiple temporal and spatial scales, and integrating the results across the different scales. Based on molecular dynamics simulations of the enzyme which transduces extracellular neuromodulatory signals into an intracellular second messenger molecule, and Brownian dynamics simulations of regulator molecules binding to the enzyme, we constructed a kinetic model ofthe enzyme-based signal transduction network. The kinetic model showed that two co-occuring neuromodulatory signals, a dopamine peak and an acetylcholine pause, are required to produce the second messenger and thus enable strengthening of corticostriatal synapses onto dSPNs, and that only the, I denna avhandling studerades synaptisk plasticitet samt förmågan till lokalberäkning i de striatala projektionsneuronens (SPNs) dendriter. Detta har direkt relevans för förståelsen av s.k. målstyrd inlärning. Målstyrd inlärning och belöningsinlärning innebär att man, baserat på sensorisk informationfrån kroppen och omgivningen, lär sig selektera hur man skall agera/handla så att man uppnår ett mål eller erhåller en belöning (t.ex. mat eller vatten). Hos alla däggdjur är motoriska/exekutiva centra i hjärnstammen och thalamus under konstant tonisk inhibition via basala ganglierna, men när manaktiverar SPNs i den s.k. direkta vägen genom basala ganglierna (dSPN) så disinhiberas de motoriska centra som behövs för att initiera specifika mål-styrda beteenden. Inlärningsprocessen för detta guidas av bl.a. dopamin, en neuromodulator som signalerar huruvida ett beteende ger positivt eller negativt resultat. Synapserna från kortikala projektionsneuron till dSPNs, s.k. kortikostriatala synapser, påverkas av dopamin och de kan förstärkas eller försvagas baserat på om dopaminsignalen signalerar ett positivt eller negativt utfall. Detta antingen främjar eller motverkar att man väljer samma beteende/handlingar under en liknande situation i framtiden. Via samarbete med andra beräkningsbiologigrupper kunde vi studera de biokemiska signaleringsnätverken i de kortikostriatala synapserna med hjälp av multiskal modeller och simuleringar. Synapsens intracellulära biokemiska signaleringsnätverk kontrollerar synaptisk plasticitet och påverkas av neuro-modulering. Genom att använda multiskalsimuleringar kunde vi studera systemet över multipla temporala och spatiala skalor, och integrera resultaten över de olika skalorna. Baserat på molekylärdynamiska simuleringar av det enzym som överför de extracellulära neuromodulatoriska signalerna till en intracellulär s.k. 2nd messenger molekyl, och Brownianska dynamiksimuleringar av de regulatoriska molekyler som binder till enzymet, kunde vi konstru, QC 20240904

The influence of natural image transformations on receptive field responses is crucial for modelling visual operations in computer vision and biological vision. In this regard, covariance properties with respect to geometric image transformations in the earliest layers of the visual hierarchy are essential for expressing robust image operations, and for formulating invariant visual operations at higher levels. This paper defines and proves a set of joint covariance properties for spatio-temporal receptive fields in terms of spatio-temporal derivative operators applied to spatio-temporally smoothed image data under compositions of spatial scaling transformations, spatial affine transformations, Galilean transformations and temporal scaling transformations. Specifically, the derived relations show how the parameters of the receptive fields need to be transformed, in order to match the output from spatio-temporal receptive fields under composed spatio-temporal image transformations. For this purpose, we also fundamentally extend the notion of scale-normalized derivatives to affine-normalized derivatives, that are computed based on spatial smoothing with affine Gaussian kernels, and analyze the covariance properties of the resulting affine-normalized derivatives for the affine group as well as for important subgroups thereof. We conclude with a geometric analysis, showing how the derived joint covariance properties make it possible to relate or match spatio-temporal receptive field responses, when observing, possibly moving, local surface patches from different views, under locally linearized perspective or projective transformations, as well as when observing different instances of spatio-temporal events, that may occur either faster or slower between different views of similar spatio-temporal events. We do furthermore describe how the parameters in the studied composed spatio-temporal image transformation models directly relate to geometric entities in the image forma, QC 20240709

We propose a normative model for spatial representation in the hippocampal formation that combines optimality principles, such as maximizing coding range and spatial information per neuron, with an algebraic framework for computing in distributed representation. Spatial position is encoded in a residue number system, with individual residues represented by high-dimensional, complex-valued vectors. These are composed into a single vector representing position by a similarity-preserving, conjunctive vector-binding operation. Self-consistency between the representations of the overall position and of the individual residues is enforced by a modular attractor network whose modules correspond to the grid cell modules in entorhinal cortex. The vector binding operation can also associate different contexts to spatial representations, yielding a model for entorhinal cortex and hippocampus. We show that the model achieves normative desiderata including superlinear scaling of patterns with dimension, robust error correction, and hexagonal, carry-free encoding of spatial position. These properties in turn enable robust path integration and association with sensory inputs. More generally, the model formalizes how compositional computations could occur in the hippocampal formation and leads to testable experimental predictions., arXiv:2406.18808 PMID: 38979486

Conformal prediction has seen many applications in pharmaceutical science, being able to calibrate outputs of machine learning models and producing valid prediction intervals. We here present the open source software CPSign that is a complete implementation of conformal prediction for cheminformatics modeling. CPSign implements inductive and transductive conformal prediction for classification and regression, and probabilistic prediction with the Venn-ABERS methodology. The main chemical representation is signatures but other types of descriptors are also supported. The main modeling methodology is support vector machines (SVMs), but additional modeling methods are supported via an extension mechanism, e.g. DeepLearning4J models. We also describe features for visualizing results from conformal models including calibration and efficiency plots, as well as features to publish predictive models as REST services. We compare CPSign against other common cheminformatics modeling approaches including random forest, and a directed message-passing neural network. The results show that CPSign produces robust predictive performance with comparative predictive efficiency, with superior runtime and lower hardware requirements compared to neural network based models. CPSign has been used in several studies and is in production-use in multiple organizations. The ability to work directly with chemical input files, perform descriptor calculation and modeling with SVM in the conformal prediction framework, with a single software package having a low footprint and fast execution time makes CPSign a convenient and yet flexible package for training, deploying, and predicting on chemical data. CPSign can be downloaded from GitHub at https://github.com/arosbio/cpsign. Scientific contribution: CPSign provides a single software that allows users to perform data preprocessing, modeling and make predictions directly on chemical structures, using conformal and probabilistic prediction. Building

Recombination of short DNA fragments via horizontal gene transfer (HGT) can introduce beneficial alleles, create genomic disharmony through negative epistasis, and create adaptive gene combinations through positive epistasis. For non-core (accessory) genes, the negative epistatic cost is likely to be minimal because the incoming genes have not co-evolved with the recipient genome and are frequently observed as tightly linked cassettes with major effects. By contrast, interspecific recombination in the core genome is expected to be rare because disruptive allelic replacement is likely to introduce negative epistasis. Why then is homologous recombination common in the core of bacterial genomes? To understand this enigma, we take advantage of an exceptional model system, the common enteric pathogens Campylobacter jejuni and C. coli that are known for very high magnitude interspecies gene flow in the core genome. As expected, HGT does indeed disrupt co-adapted allele pairings, indirect evidence of negative epistasis. However, multiple HGT events enable recovery of the genome’s co-adaption between introgressing alleles, even in core metabolism genes (e.g., formate dehydrogenase). These findings demonstrate that, even for complex traits, genetic coalitions can be decoupled, transferred, and independently reinstated in a new genetic background—facilitating transition between fitness peaks. In this example, the two-step recombinational process is associated with C. coli that are adapted to the agricultural niche.

Campylobacter is a leading cause of food-borne gastroenteritis worldwide, linked to the consumption of contaminated poultry meat. Targeting this pathogen at source, vaccines for poultry can provide short-term caecal reductions in Campylobacter numbers in the chicken intestine. However, this approach is unlikely to reduce Campylobacter in the food chain or human incidence. This is likely as vaccines typically target only a subset of the high genomic strain diversity circulating among chicken flocks, and rapid evolution diminishes vaccine efficacy over time. To address this, we used a genomic approach to develop a whole-cell autogenous vaccine targeting isolates harbouring genes linked to survival outside of the host. We hyper-immunised a whole major UK breeder farm to passively target offspring colonisation using maternally-derived antibody. Monitoring progeny, broiler flocks revealed a near-complete shift in the post-vaccination Campylobacter population with an ~50% reduction in isolates harbouring extra-intestinal survival genes and a significant reduction of Campylobacter cells surviving on the surface of meat. Based on these findings, we developed a logistic regression model that predicted that vaccine efficacy could be extended to target 65% of a population of clinically relevant strains. Immuno-manipulation of poultry microbiomes towards less harmful commensal isolates by competitive exclusion, has major potential for reducing pathogens in the food production chain.

Objectives Integrating pathogen genomic surveillance with bioinformatics can enhance public health responses by identifying risk and guiding interventions. This study focusses on the two predominant Campylobacter species, which are commonly found in the gut of birds and mammals and often infect humans via contaminated food. Rising incidence and antimicrobial resistance (AMR) are a global concern and there is an urgent need to quantify the main routes to human infection. Methods During routine US national surveillance (2009-2019), 8,856 Campylobacter genomes from human infections and 16,703 from possible sources were sequenced. Using machine learning and probabilistic models, we target genetic variation associated with host adaptation to attribute the source of human infections and estimate the importance of different disease reservoirs. Results Poultry was identified as the primary source of human infections, responsible for an estimated 68% of cases, followed by cattle (28%), and only a small contribution from wild birds (3%) and pork sources (1%). There was also evidence of an increase in multidrug resistance, particularly among isolates attributed to chickens. Conclusions National surveillance and source attribution can guide policy, and our study suggests that interventions targeting poultry will yield the greatest reductions in campylobacteriosis and spread of AMR in the US.

Objectives Antimicrobial resistant (AMR) Campylobacter is a global health threat; however, there is limited information on genomic determinants of resistance in low- and middle-income countries. We evaluated genomic determinants of AMR using a collection of whole genome sequenced Campylobacter jejuni and C. coli isolates from Iquitos, Peru. Methods Campylobacter isolates from two paediatric cohort studies enriched with isolates that demonstrated resistance to ciprofloxacin and azithromycin were sequenced and mined for AMR determinants. Results The gyrA mutation leading to the Thr86Ile amino acid change was the only gyrA mutation associated with fluoroquinolone resistance identified. The A2075G mutation in 23S rRNA was present, but three other 23S rRNA mutations previously associated with macrolide resistance were not identified. A resistant-enhancing variant of the cmeABC efflux pump genotype (RE-cmeABC) was identified in 36.1% (35/97) of C. jejuni genomes and 17.9% (12/67) of C. coli genomes. Mutations identified in the CmeR-binding site, an inverted repeat sequence in the cmeABC promoter region that increases expression of the operon, were identified in 24/97 C. jejuni and 14/67 C. coli genomes. The presence of these variants, in addition to RE-cmeABC, was noted in 18 of the 24 C. jejuni and 9 of the 14 C. coli genomes. Conclusions Both RE-cmeABC and mutations in the CmeR-binding site were strongly associated with the MDR phenotype in C. jejuni and C. coli. This is the first report of RE-cmeABC in Peru and suggests it is a major driver of resistance to the principal therapies used to treat human campylobacteriosis in this setting.

Sex estimation of the individuals in a sample is fundamental for any bioarchaeological study to define a particular demographic assemblage or to classify isolated remains. Long bones are an excellent alternative for sex estimation when the most dimorphic anatomical parts are not preserved or are highly altered. Here we propose a set of discriminant functions and classification models to estimate the sex of prehistoric individuals using linear discriminant analysis and machine learning approaches. Different osteometric variables were taken from the humeri, ulnae, radii, femurs and tibias of a sample of 109 articulated skeletons buried in the collective tomb of Camino del Molino (Region of Murcia, SE-Spain), dated to the 3rd millennium BC. Sex was estimated based on standard anthropological methods and ancient DNA analysis of a control sample. Fifty-two discriminant functions with prediction thresholds higher than 0.8 on the ROC curve were obtained using independent (22) and combined variables (30). The best LDA models for sex prediction were those based on proximal epiphyseal widths or their combination with other variables, reaching values close to 0.98 on the ROC curve. The random forest-based model obtained an accuracy of 0.94 and confirmed the importance of epiphyseal widths in sex classification. This analysis is more comprehensive than univariate LDA, as it allows for ranking the importance of bones in sex discrimination and considers correlations between long bones rather than treating them as independent observations. In contrast, applying LDA to each bone makes it easier to predict the sex of other coeval collections that do not have such a complete sample. This work aims to overcome the scarcity of methods that can be applied to sex estimation of the large volume of isolated remains from Camino del Molino and for other Mediterranean skeletal series from the Late Prehistory with high biological affinity and that share similar environmental conditions.

Massively parallel sequencing (MPS) is increasingly applied in forensic short tandem repeat (STR) analysis. The presence of stutter artefacts and other PCR or sequencing errors in the MPS-STR data partly limits the detection of low DNA amounts, e.g., in complex mixtures. Unique molecular identifiers (UMIs) have been applied in several scientific fields to reduce noise in sequencing. UMIs consist of a stretch of random nucleotides, a unique barcode for each starting DNA molecule, that is incorporated in the DNA template using either ligation or PCR. The barcode is used to generate consensus reads, thus removing errors. The SiMSen-Seq (Simple, multiplexed, PCR-based barcoding of DNA for sensitive mutation detection using sequencing) method relies on PCR-based introduction of UMIs and includes a sophisticated hairpin design to reduce unspecific primer binding as well as PCR protocol adjustments to further optimize the reaction. In this study, SiMSen-Seq is applied to develop a proof-of-concept seven STR multiplex for MPS library preparation and an associated bioinformatics pipeline. Additionally, machine learning (ML) models were evaluated to further improve UMI allele calling. Overall, the seven STR multiplex resulted in complete detection and concordant alleles for 47 single-source samples at 1 ng input DNA as well as for low-template samples at 62.5 pg input DNA. For twelve challenging mixtures with minor contributions of 10 pg to 150 pg and ratios of 1–15% relative to the major donor, 99.2% of the expected alleles were detected by applying the UMIs in combination with an ML filter. The main impact of UMIs was a substantially lowered number of artefacts as well as reduced stutter ratios, which were generally below 5% of the parental allele. In conclusion, UMI-based STR sequencing opens new means for improved analysis of challenging crime scene samples including complex mixtures. © 2024 The Authors, We thank Froste Svensson for input on the bioinformatics pipeline and Nelly Gyll\u00F6 Lind and Markus Andr\u00E9 Soma for practical assistance in pre-studies. This study was funded by VINNOVA: Project title \u201CUltrak\u00E4nsliga analyser f\u00F6r b\u00E4ttre h\u00E4lsa och kriminalteknik (ULTRA-UDI)\u201D and reference number 2020\u201304141, and V\u00E4stra G\u00F6talandsregionen RUN 2021\u201300059. Anders St\u00E5hlberg was funded by Region V\u00E4stra G\u00F6taland; Swedish Cancer Society [22\u20132080]; Swedish Childhood Cancer Foundation [2022\u20130030]; Swedish Research Council [2021\u201301008]; the Swedish state under the agreement between the Swedish government and the county councils, the ALF-agreement [965065]; Sweden's Innovation Agency [2018\u201300421, 2020\u201304141] and the Sj\u00F6berg Foundation. Points of view in this document are those of the NIST-affiliated authors and do not necessarily represent the official position or policies of the U.S. Department of Commerce. Certain equipment, instruments, software, or materials are identified in this paper to specify the experimental procedure adequately. Such identification is not intended to imply recommendation or endorsement of any product or service by NIST, nor is it intended to imply that the materials or equipment identified are necessarily the best available for the purpose. All work performed at NIST has been reviewed and approved by the U. S. National Institute of Standards and Technology Research Protections Office. This study was determined to be \u201Cnot human subjects research\u201D (often referred to as research not involving human subjects) as defined in U. S. Department of Commerce Regulations, 15 CFR 27, also known as the Common Rule (45 CFR 46, Subpart A), for the Protection of Human Subjects by the NIST Human Research Protections Office and therefore not subject to oversight by the NIST Institutional Review Board (MML-16\u20130080). All data handling performed in Sweden w

There has been a long-standing bottleneck in the quantitative analysis of the frequencies of homoblock polyads beyond triads using 1H and 13C NMR for linear polysaccharides, primarily because monosaccharides within a long homoblock share similar chemical environments due to identical neighboring units, resulting in indistinct NMR peaks. In this study, through rigorous mathematical induction, inequality relations were established that enabled the calculation of frequency ranges of homoblock polyads from historically reported NMR-derived frequency values of diads and/or triads of alginates, chitosans, homogalacturonans, and galactomannans. The calculated homoblock frequency ranges were then applied to evaluate three chain growth statistical models, including the Bernoulli chain, first-order Markov chain, and second-order Markov chain, for predicting homoblock frequencies in these polysaccharides. Furthermore, based on the mathematically derived inequality relations, a novel 2D array was constructed, enabling the graphical visualization of homoblock features in polysaccharides. It was demonstrated, as a proof of concept, that the novel 2D array, along with a 1D code generated from it, could serve as an effective feature engineering tool for polymer classification using machine learning algorithms., QC 20240618

Human age and gender prediction from facial images has garnered significant attention due to its importance in various applications. Traditional models struggle with large-scale variations in unfiltered images. Convolutional Neural Networks (CNNs) have emerged as effective tools for facial analysis due to their robust performance. This paper presents a novel CNN approach for robust age and gender classification using unconstrained real-world images. The CNN architecture includes convolution, pooling, and fully connected layers for feature extraction, dimension reduction, and mapping to output classes. Adience and UTKFace datasets were utilized, with the best training and testing accuracies achieved using an 80% training and 20% testing data split. Robust image pre-processing and data augmentation techniques were applied to handle dataset variations. The proposed approach outperformed existing methods, achieving age prediction accuracies of 86.42% and 81.96%, and gender prediction accuracies of 97.65% and 96.32% on the Adience and UTKFace datasets, respectively., Full text license: CC BY-NC-ND 4.0

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., CC BY 4.0Corresponding author at: University Hospital of Umeå, 901 85 Umeå, Sweden. E-mail address: emmy.borgmastars@umu.se (E. Borgmästars)The authors thank Hanna Nyström, and Daniel Öhlund at Umeå University for valuable assistance in data collection. We thank Xiaoshuang Feng at International Agency for Research of Cancer, Lyon, France for guidance in statistical analyses. The authors would also like to thank Swedish Metabolomics Centre, Umeå, Sweden (www.swedishmetabolomicscentre.se) and Biobanken Norr. Funding: This study was funded by Umeå University, the Swedish Research Council [2016-02990, 2019-01690], the Swedish Cancer Society [CAN 2016/643, 19 0273], Region Västerbotten [RV-583411, RV-549731, RV-583411, RV-841551, RV 967602], Finska Läkaresällskapet, Medicinska Understödsföreningen Liv och Hälsa, the Sjöberg Foundation, The JC Kempe Memorial Foundation Scholarship Fund, The Royal Swedish Academy of Sciences (PE Lindahl Foundation, LM2021-0010 and LM2023-0012), The Swedish Society of Medicine (SLS-960379), Cancerforskningsfonden i Norrland (LP 23-2337), Bengt Ihre foundation (SLS-960529 and SLS-986656), and Bengt Ihre Research Fellowship Grant. The sponsors had no role in the study design, collection, analysis and interpretation of data, writing of the report, or in the decision to submit the article for publication.

Transmission electron microscopy (TEM) is an imaging technique used to visualize and analyze nano-sized structures and objects such as virus particles. Light microscopy can be used to diagnose diseases or characterize e.g. blood cells. Since samples under microscopes exhibit certain symmetries, such as global rotation invariance, equivariant neural networks are presumed to be useful. In this study, a baseline convolutional neural network is constructed in the form of the commonly used VGG16 classifier. Thereafter, it is modified to be equivariant to the p4 symmetry group of rotations of multiples of 90 degrees using group convolutions. This yields a number of benefits on a TEM virus dataset, including higher top validation set accuracy by on average 7.6% and faster convergence during training by on average 23.1% of that of the baseline. Similarly, when training and testing on images of blood cells, the convergence time for the equivariant neural network is 7.9% of that of the baseline. From this it is concluded that augmentation strategies for rotation can be skipped. Furthermore, when modelling the accuracy versus amount of TEM virus training data with a power law, the equivariant network has a slope of - 0.43 compared to - 0.26 of the baseline. Thus the equivariant network learns faster than the baseline when more training data is added. This study extends previous research on equivariant neural networks applied to images which exhibit symmetries to isometric transformations.

The COVID-19 pandemic highlighted the need for a rapid, convenient, and scalable diagnostic method for detecting a novel pathogen amidst a global pandemic. While command-line interface tools offer automation for SARS-CoV-2 Oxford Nanopore Technology sequencing data analysis, they are inapplicable to users with limited programming skills. A solution is to establish such automated workflows within a graphical user interface software. We developed two workflows in the software Geneious Prime 2022.1.1, adapted for data obtained from the Midnight and Artic's nCoV-2019 sequencing protocols. Both workflows perform trimming, read mapping, consensus generation, and annotation on SARS-CoV-2 Nanopore sequencing data. Additionally, one workflow includes phylogenetic assignment using the bioinformatic tools pangolin and Nextclade as plugins. The basic workflow was validated in 2020, adhering to the requirements of the European Centre for Disease Prevention and Control for SARS-CoV-2 sequencing and analysis. The enhanced workflow, providing phylogenetic assignment, underwent validation at Uppsala University Hospital by analysing 96 clinical samples. It provided accurate diagnoses matching the original results of the basic workflow while also reducing manual clicks and analysis time. These bioinformatic workflows streamline SARS-CoV-2 Nanopore data analysis in Geneious Prime, saving time and manual work for operators lacking programming knowledge.