Background: Periodontitis (PD) and rheumatoid arthritis (RA) are chronic inflammatory conditions, characterized by dysregulated immune response and excessive production of inflammatory mediators. The oral disease PD is triggered by periodontal pathogens, leading to the destruction of tissues surrounding the teeth, whereas RA is a systemic autoimmune disease primarily affecting the joints. The objective of this study was to investigate the prevalence of PD and map the profile of salivary and serum inflammatory mediators of patients with RA, with respect to periodontal severity (PD stage II and PD stage III/IV). Methods: For this cross-sectional cohort study, 62 patients diagnosed with RA were recruited. All participants underwent a full-mouth dental examination. Levels of various inflammatory mediators, including tumor necrosis factor (TNF) superfamily proteins, interferon (IFN) family proteins, regulatory T cell (Treg) cytokines, and matrix metalloproteinases were determined in saliva and serum samples from each participant using a human inflammation multiplex immunoassay panel. Results: In the current RA cohort, all participants were diagnosed with PD, of which 35.5% were classified as PD stage II and 64.5% as PD stages III/IV. Inflammatory mediator levels were significantly higher in both saliva and serum samples from patients with RA and PD stages III/IV, compared to those with RA and stage II within the same cohort. These included higher serum levels of sCD30, IL-10, IL-19, osteopontin and elevated salivary levels of BAFF/TNFSF13B and IFN-α2. Additionally, APRIL/TNFSF13 levels were increased in both saliva and serum. Conclusions: Among the studied patients with RA, the majority exhibited severe PD (stage III/IV), underscoring the importance of periodontal prophylaxis and treatment for this group of patients. Higher levels of inflammatory mediators were observed in both saliva and serum in those with PD stages III/IV, suggesting a potential link between the severity of PD and systemic inflammation in RA. Further research is needed to explore the clinical implications of these findings.

BackgroundZostera marina is an important ecosystem engineer influencing shallow water environments and possibly shaping the microbiota in surrounding sediments and water. Z. marina is typically found in marine systems, but it can also proliferate under brackish conditions. Changes in salinity generally have a strong impact on the biota, especially at the salty divide between salinity 6 and 9. To better understand the impact of the salty divide on the interaction between Z. marina and the surrounding sediment and water microbiota, we investigated the effects of Z. marina meadows on the surrounding microbiota across a salinity range of 6-15 in the Baltic Sea during the summer using 16S and 18S rRNA gene amplicon sequencing.ResultsSalinity was the most important factor for structuring the microbiota within both water and sediment. The presence of Z. marina affected the composition of the bacterial and eukaryotic community and bacterial alpha diversity in the sediment. However, this effect was confined to alpha-mesohaline conditions (salinity 9-15). The impact of Z. marina below salinity 9 on water and sediment microbiota was insignificant.ConclusionsIncreasing salinity was associated with a longer leaf length of Z. marina, causing an increased canopy height, which affects the sediment microbiota through reduced water velocity. Hence, we propose that the canopy effect may be the major predictor explaining Z. marina's interactions with the surrounding microbiota at salinity 9-15. These findings emphasize the importance of the physical effects of Z. marina meadow ecosystem services and have important implications for Z. marina management under brackish conditions in a changing climate.

Marine ecosystems harbour a vast microbial diversity which play a crucial role in ecosystemfunctioning. Advancements in DNA sequencing technologies have transformed our ability to analyse microbial populations comprehensively. Metagenomic sequencing has emerged as a pivotal tool for characterising microbial communities across various environments. Bioinformatics, an interdisciplinary field, facilitates the analysis and interpretation of large biological datasets, including microbiome data.

This thesis aims to enhance bioinformatics approaches for analysing marine microbiomes. It comprises four papers covering bioinformatic developments and genomic data analysis across multiple topics, including metagenomics, pangenomics, comparative genomics and population genomics:

Paper I evaluated three assembly strategies for constructing gene catalogues from metagenomic samples: individual sample assembly with gene clustering, co-assembly of all samples, and a new hybrid approach, mix assembly. The efficacy of the mix-assembly approach was highlighted for maximising information extraction from metagenomic samples, offering opportunities for further exploration in microbial ecology and environmental genomics.

Using the mix-assembly approach, we conducted a comprehensive analysis of 124 metagenomic samples sourced from the Baltic Sea, resulting in the refinement of the Baltic Sea Gene Set (BAGS v1.1), which now encompasses 66.53 million genes annotated for both functionality and taxonomy. In Paper II, we introduced an open-access initiative that provided the mix-assembly pipeline code. We also developed the BAGS-Shiny web application to facilitate user interaction with this extensive gene catalogue.

Paper III focused on whole-genome sequencing and assembly of 82 environmental V. vulnificus strains from the Baltic Sea, enabling comprehensive comparative genomic analysis. I developed the PhyloBOTL pipeline, which uses a phylogeny-based approach to identify genes associated with pathogenicity. Comparative genomics of 208 clinical isolates and 199 environmental isolates revealed 58 enriched orthologs in pathogenic strains, including known virulence factors and novel genes. Potential biomarkers for pathogenic V. vulnificus were identified, and primers suitable for PCR-based environmental monitoring were designed (in silico).

In Paper IV population genomics analysis was carried out, using the Input_Pogenom pipeline and POGENOM tool, to explore intraspecific biogeographical patterns. Geographical barriers were found to significantly influence aquatic bacteria distribution, with greater genetic differentiation observed between Baltic and Caspian seas than within the Baltic Sea's salinity gradient.

Havsmiljöer hyser en enorm mikrobiell mångfald som spelar en avgörande roll för ekosystemens funktion. Framsteg inom DNA-sekvenseringstekniker har revolutionerat vår förmåga att analysera

den mikrobiella populationen på ett omfattande sätt. Metagenomisk sekvensering har framträtt som ett centralt verktyg för att karakterisera mikrobiella samhällen i olika miljöer. Bioinformatik, ett tvärvetenskapligt fält, underlättar analys och tolkning av stora biologiska dataset, inklusive mikrobiomdata.

Den här avhandlingen syftar till att förbättra bioinformatiska metoder för att analysera marina mikrobiom. Den består av fyra artiklar som täcker bioinformatisk utveckling och analys av genomdata inom flera områden, inklusive metagenomik, pangenom, jämförande genomik och populationsgenetik:

Artikel I utvärderade tre monteringsstrategier för att konstruera genkataloger från metagenomiska prover: montering av enskilda prover med genglustering, sammontering av alla prover och en ny hybridmetod, mixmontering. Effektiviteten hos mixmonteringsmetoden lyftes fram för att maximera informationsutvinning från metagenomiska prover, vilket öppnar för vidare utforskning inom mikrobiell ekologi och miljögenomik.

Med hjälp av mixmonteringsmetoden genomförde vi en omfattande analys av 124 metagenomiska prover från Östersjön, vilket resulterade i förfiningen av Östersjöns gensets (BAGS v1.1), som nu omfattar 66,53 miljoner gener annoterade för både funktion och taxonomi. I artikel II introducerade vi ett öppet initiativ som tillhandahöll koden för mixmonterings-pipeline. Vi utvecklade också BAGS-Shiny webbapplikationen för att underlätta användarinteraktion med denna omfattande gensetskatalog.

Artikel III fokuserade på helgenomsekvensering och montering av 82 miljörelaterade V. vulnificus-stammar från Östersjön, vilket möjliggjorde omfattande jämförande genomisk analys.

Jag utvecklade PhyloBOTL-pipelinen, som använder en fylogenibaserad metod för att identifiera gener associerade med patogenicitet. Jämförande genomik av 208 kliniska isolat och 199 miljöisolerade isolat avslöjade 58 anrikade ortologer i patogena stammar, inklusive kända virulensfaktorer och nya gener. Potentiella biomarkörer för patogena V. vulnificus identifierades och primers lämpliga för PCR-baserad miljöövervakning designades (in silico).

I artikel IV utfördes populationsgenetisk analys med hjälp av Input_Pogenom-pipelinen och POGENOM-verktyget för att utforska intraspecifika biogeografiska mönster. Geografiska barriärer visade sig signifikant påverka vattenlevande bakteriers utbredning, med större genetisk differentiering observerad mellan Östersjön och Kaspiska havet än inom Östersjöns salthaltsgradient.

Due to climate change the pathogenic bacterium Vibrio vulnificus proliferates along brackish coastlines, posing risks to public health, tourism, and aquaculture. Here we investigated previously suggested regulation measures to reduce the prevalence of V. vulnificus, locally through seagrass and regionally through the reduction of eutrophication and consequential formation of algal blooms. Field samples collected in the summer of 2021 covered the salinity and eutrophication gradients of the Baltic Sea, one of the largest brackish areas worldwide. Physico-, biological- and hydrochemical parameters were measured and variables explaining V. vulnificus occurrence were identified by machine learning. The best V. vulnificus predictors were eutrophication-related features, such as particulate organic carbon and nitrogen, as well as occurrence of potential phytoplankton blooms and associated species. V. vulnificus abundance did not vary significantly between vegetated and non-vegetated areas. Thus, reducing nutrient inputs could be an effective method to control V. vulnificus populations in eutrophied brackish coasts.

The number of Vibrio-related infections in humans, e.g., by Vibrio vulnificus, has increased along the coasts of the Baltic Sea. Due to climate change, vibriosis risk is expected to increase. It is, therefore, pertinent to design a strategy for mitigation of the vibriosis threat in the Baltic Sea area, but a prerequisite is to identify the environmental conditions promoting the occurrence of pathogenic Vibrio spp., like V. vulnificus. To address this, we sampled three coastal Baltic sites in Finland, Germany, and Denmark with salinities between 6 and 21 from May to October 2022. The absolute and relative abundances of Vibrio spp. and V. vulnificus in water were compared to environmental conditions, including the presence of the eelgrass Zostera marina, which has been suggested to reduce pathogenic Vibrio species abundance. In the water column, V. vulnificus only occurred at the German station between July and August at salinity 8.1-11.2. Temperature and phosphate (PO43-) were identified as the most influencing factors for Vibrio spp. and V. vulnificus. The accumulation of Vibrio spp. in the sediment and the co-occurrence with sediment bacteria in the water column indicate that sediment resuspension contributed to V. vulnificus abundance. Interestingly, V. vulnificus co-occurred with specific cyanobacteria taxa, as well as specific bacteria associated with cyanobacteria. Although we found no reduction in Vibrio spp. or V. vulnificus associated with eelgrass beds, our study underscores the importance of extended heatwaves and sediment resuspension, which may elevate the availability of PO43-, for Vibrio species levels at intermediate salinities in the Baltic Sea. IMPORTANCE: Elevated sea surface temperatures are increasing the prevalence of pathogenic Vibrio at higher latitudes. The recent increase in Vibrio-related wound infections and deaths along the Baltic coasts is, therefore, of serious health concern. We used culture-independent data generated from three Baltic coastal sites in Denmark, Germany, and Finland from May to October (2022), with a special focus on Vibrio vulnificus, and combined it with environmental data. Our temporal model shows that temperature, combined with sediment resuspension, drives the prevalence of V. vulnificus at intermediate salinities in the coastal Baltic Sea.

The crossing of environmental barriers poses major adaptive challenges. Rareness of freshwater-marine transi-tions separates the bacterial communities, but how these are related to brackish counterparts remains elusive, as do the molecular adaptations facilitating cross-biome transitions. We conducted large-scale phylogenomic analysis of freshwater, brackish, and marine quality-filtered metagenome-assembled genomes (11,248). Average nucleotide identity analyses showed that bacterial species rarely existed in multiple biomes. In contrast, distinct brackish basins cohosted numerous species, but their intraspecific population structures displayed clear signs of geographic separation. We further identified the most recent cross-biome transitions, which were rare, ancient, and most commonly directed toward the brackish biome. Transitions were accompanied by systematic changes in amino acid composition and isoelectric point distributions of inferred proteomes, which evolved over millions of years, as well as convergent gains or losses of specific gene functions. Therefore, adaptive chal-lenges entailing proteome reorganization and specific changes in gene content constrains the cross-biome tran-sitions, resulting in species-level separation between aquatic biomes.

Background: For many environments, biome-specific microbial gene catalogues are being recovered using shotgun metagenomics followed by assembly and gene calling on the assembled contigs. The assembly is typically conducted either by individually assembling each sample or by co-assembling reads from all the samples. The co-assembly approach can potentially recover genes that display too low abundance to be assembled from individual samples. On the other hand, combining samples increases the risk of mixing data from closely related strains, which can hamper the assembly process. In this respect, assembly on individual samples followed by clustering of (near) identical genes is preferable. Thus, both approaches have potential pros and cons, but it remains to be evaluated which assembly strategy is most effective. Here, we have evaluated three assembly strategies for generating gene catalogues from metagenomes using a dataset of 124 samples from the Baltic Sea: (1) assembly on individual samples followed by clustering of the resulting genes, (2) co-assembly on all samples, and (3) mix assembly, combining individual and co-assembly. Results: The mix-assembly approach resulted in a more extensive nonredundant gene set than the other approaches and with more genes predicted to be complete and that could be functionally annotated.The mix assembly consists of 67 million genes (Baltic Sea gene set, BAGS) that have been functionally and taxonomically annotated. The majority of the BAGS genes are dissimilar (< 95% amino acid identity) to the Tara Oceans gene dataset, and hence, BAGS represents a valuable resource for brackish water research. Conclusion: The mix-assembly approach represents a feasible approach to increase the information obtained from metagenomic samples.

ObjectivesPeriodontitis and rheumatoid arthritis (RA) are two widespread chronic inflammatory diseases with a previously suggested association. The objective of the current study was to compare the oral microbial composition and host ' s inflammatory mediator profile of saliva samples obtained from subjects with periodontitis, with and without RA, as well as to predict biomarkers, of bacterial pathogens and/or inflammatory mediators, for classification of samples associated with periodontitis and RA. MethodsSalivary samples were obtained from 53 patients with periodontitis and RA and 48 non-RA with chronic periodontitis. The microbial composition was identified using 16S rRNA gene sequencing and compared across periodontitis patients with and without RA. Levels of inflammatory mediators were determined using a multiplex bead assay, compared between the groups and correlated to the microbial profile. The achieved data was analysed using PCoA, DESeq2 and two machine learning algorithms, OPLS-DA and sPLS-DA. ResultsDifferential abundance DESeq2 analyses showed that the four most highly enriched (log2 FC >20) amplicon sequence variants (ASVs) in the non-RA periodontitis group included Alloprevotella sp., Prevotella sp., Haemophilus sp., and Actinomyces sp. whereas Granulicatella sp., Veillonella sp., Megasphaera sp., and Fusobacterium nucleatum were the most highly enriched ASVs (log2 FC >20) in the RA group. OPLS-DA with log2 FC analyses demonstrated that the top ASVs with the highest importance included Vampirovibrio sp. having a positive correlation with non-RA group, and seven ASVs belonging to Sphingomonas insulae, Sphingobium sp., Novosphingobium aromaticivorans, Delftia acidovorans, Aquabacterium spp. and Sphingomonas echinoides with a positive correlation with RA group. Among the detected inflammatory mediators in saliva samples, TWEAK/TNFSF12, IL-35, IFN-alpha 2, pentraxin-3, gp130/sIL6Rb, sIL-6Ra, IL-19 and sTNF-R1 were found to be significantly increased in patients with periodontitis and RA compared to non-RA group with periodontitis. Moreover, correlations between ASVs and inflammatory mediators using sPLS-DA analysis revealed that TWEAK/TNFSF12, pentraxin-3 and IL-19 were positively correlated with the ASVs Sphingobium sp., Acidovorax delafieldii, Novosphingobium sp., and Aquabacterium sp. ConclusionOur results suggest that the combination of microbes and host inflammatory mediators could be more efficient to be used as a predictable biomarker associated with periodontitis and RA, as compared to microbes and inflammatory mediators alone.

Bacterioplankton are main drivers of biogeochemical cycles and important components of aquatic food webs. While sequencing-based studies have revealed how bacterioplankton communities are structured in time and space, relatively little is known about intraspecies diversity patterns and their ecological relevance. Here, we use the newly developed software POGENOM (POpulation GENomics from Metagenomes) to investigate genomic diversity and differentiation in metagenome-assembled genomes from the Baltic Sea, and investigate their genomic variation using metagenome data spanning a 1700 km transect and covering seasonal variation at one station. The majority of the investigated species, representing several major bacterioplankton clades, displayed population structures correlating significantly with environmental factors such as salinity and temperature. Population differentiation was more pronounced over spatial than temporal scales. We discovered genes that have undergone adaptation to different salinity regimes, potentially responsible for the populations' existence along with the salinity range. This in turn implies the broad existence of ecotypes that may remain undetected by rRNA gene sequencing. Our findings emphasize the importance of physiological barriers, and highlight the role of adaptive divergence as a structuring mechanism of bacterioplankton species.

Biome-specific gene catalogues have been recovered for many environments using shotgun metagenomics, followed by assembly and gene calling on the assembled contigs. We recently proposed a novel mix-assembly strategy, combining individual and co-assembly approaches, and used this approach to assemble an extensive non-redundant gene set from 124 Baltic Sea metagenome samples. The Baltic Sea Gene Set (BAGS v1.1) comprises 66.53 million functionally and taxonomically annotated genes. To enable interactive exploration of this gene catalogue we have developed an RShiny application, BAGS-Shiny, that allows users to perform searches by sequence similarity (BLAST) and/or taxonomic and functional annotation. The gene catalogue and web application will serve as valuable tools for exploring microbial gene functions in brackish ecosystems. In addition, we here make available a pipeline to create gene catalogues based on the mix-assembly approach.