The MAR databases: development and implementation of databases specific for marine metagenomics.

We introduce the marine databases; MarRef, MarDB and MarCat (https://mmp.sfb.uit.no/databases/), which are publicly available resources that promote marine research and innovation. These data resources, which have been implemented in the Marine Metagenomics Portal (MMP) (https://mmp.sfb.uit.no/), are collections of richly annotated and manually curated contextual (metadata) and sequence databases representing three tiers of accuracy. While MarRef is a database for completely sequenced marine prokaryotic genomes, which represent a marine prokaryote reference genome database, MarDB includes all incomplete sequenced prokaryotic genomes regardless level of completeness. The last database, MarCat, represents a gene (protein) catalog of uncultivable (and cultivable) marine genes and proteins derived from marine metagenomics samples. The first versions of MarRef and MarDB contain 612 and 3726 records, respectively. Each record is built up of 106 metadata fields including attributes for sampling, sequencing, assembly and annotation in addition to the organism and taxonomic information. Currently, MarCat contains 1227 records with 55 metadata fields. Ontologies and controlled vocabularies are used in the contextual databases to enhance consistency. The user-friendly web interface lets the visitors browse, filter and search in the contextual databases and perform BLAST searches against the corresponding sequence databases. All contextual and sequence databases are freely accessible and downloadable from https://s1.sfb.uit.no/public/mar/.

Comparative analysis of HiSeq3000 and BGISEQ-500 sequencing platform over whole genome sequencing metagenomics data.

Recent advances in sequencing technologies and platforms have enabled to generate metagenomics sequences using different sequencing platforms. In this study, we analyzed and compared shotgun metagenomic sequences generated by HiSeq3000 and BGISEQ-500 platforms from 12 sediment samples collected across the Norwegian coast. Metagenomics DNA sequences were normalized to an equal number of bases for both platforms and further evaluated by using different taxonomic classifiers, reference databases, and assemblers. Normalized BGISEQ-500 sequences retained more reads and base counts after preprocessing, while a slightly higher fraction of HiSeq3000 sequences were taxonomically classified. Kaiju classified a higher percentage of reads relative to Kraken2 for both platforms, and comparison of reference database for taxonomic classification showed that MAR database outperformed RefSeq. Assembly using MEGAHIT produced longer assemblies and higher total contigs count in majority of HiSeq3000 samples than using metaSPAdes, but the assembly statistics notably improved with unprocessed or normalized reads. Our results indicate that both platforms perform comparably in terms of the percentage of taxonomically classified reads and assembled contig statistics for metagenomics samples. This study provides valuable insights for researchers in selecting an appropriate sequencing platform and bioinformatics pipeline for their metagenomics studies.

Phylogenetic Revision of the Genus Aliivibrio: Intra- and Inter-Species Variance Among Clusters Suggest a Wider Diversity of Species.

Genus Aliivibrio is known to harbor species exhibiting bioluminescence as well as pathogenic behavior affecting the fish farming industry. Current phylogenetic understanding of Aliivibrio has largely remained dormant after reclassification disentangled it from the Vibrio genus in 2007. There is growing evidence of wider diversity, but until now the lack of genomes and selective use of type strains have limited the ability to compare and classify strains firmly. In this study, a total of 143 bacterial strains, including 51 novel sequenced strains, were used to strengthen phylogenetic relationships in Aliivibrio by exploring intra-species and inter-species relations. Multilocus sequence analysis (MLSA), applying the six housekeeping genes 16S ribosomal RNA (rRNA), gapA, gyrB, pyrH, recA, and rpoA, inferred 12 clades and a singular branch in Aliivibrio. Along with four new phylogenetic clades, the MLSA resolved prior inconsistencies circumscribing Aliivibrio wodanis and formed a unique clade we propose as the novel species Aliivibrio sp. "friggae." Furthermore, phylogenetic assessment of individual marker genes showed gyrB, pyrH, and recA superior to the 16S rRNA gene, resolving accurately for most species clades in Aliivibrio. In this study, we provide a robust phylogenetic groundwork for Aliivibrio as a reference point to classification of species.

DeltaProt: a software toolbox for comparative genomics.

BACKGROUND: Statistical bioinformatics is the study of biological data sets obtained by new micro-technologies by means of proper statistical methods. For a better understanding of environmental adaptations of proteins, orthologous sequences from different habitats may be explored and compared. The main goal of the DeltaProt Toolbox is to provide users with important functionality that is needed for comparative screening and studies of extremophile proteins and protein classes. Visualization of the data sets is also the focus of this article, since visualizations can play a key role in making the various relationships transparent. This application paper is intended to inform the reader of the existence, functionality, and applicability of the toolbox. RESULTS: We present the DeltaProt Toolbox, a software toolbox that may be useful in importing, analyzing and visualizing data from multiple alignments of proteins. The toolbox has been written in MATLAB™ to provide an easy and user-friendly platform, including a graphical user interface, while ensuring good numerical performance. Problems in genome biology may be easily stated thanks to a compact input format. The toolbox also offers the possibility of utilizing structural information from the SABLE or other structure predictors. Different sequence plots can then be viewed and compared in order to find their similarities and differences. Detailed statistics are also calculated during the procedure. CONCLUSIONS: The DeltaProt package is open source and freely available for academic, non-commercial use. The latest version of DeltaProt can be obtained from http://services.cbu.uib.no/software/deltaprot/. The website also contains documentation, and the toolbox comes with real data sets that are intended for training in applying the models to carry out bioinformatical and statistical analyses of protein sequences.Equipped with the new algorithms proposed here, DeltaProt serves as an auxiliary analysis tool capable of visualizing and comparing orthologus protein sequences. The framework of the algorithms also enables easy incorporation of extra information on structure, if such data is available.

Effects of salt on the kinetics and thermodynamic stability of endonuclease I from Vibrio salmonicida and Vibrio cholerae.

Adaptation to extreme environments affects the stability and catalytic efficiency of enzymes, often endowing them with great industrial potential. We compared the environmental adaptation of the secreted endonuclease I from the cold-adapted marine fish pathogen Vibrio salmonicida (VsEndA) and the human pathogen Vibrio cholerae (VcEndA). Kinetic analysis showed that VsEndA displayed unique halotolerance. It retained a considerable amount of activity from low concentrations to at least 0.6 m NaCl, and was adapted to work at higher salt concentrations than VcEndA by maintaining a low K(m) value and increasing k(cat). In differential scanning calorimetry, salt stabilized both enzymes, but the effect on the calorimetric enthalpy and cooperativity of unfolding was larger for VsEndA, indicating salt dependence. Mutation of DNA binding site residues (VsEndA, Q69N and K71N; VcEndA, N69Q and N71K) affected the kinetic parameters. The VsEndA Q69N mutation also increased the T(m) value, whereas other mutations affected mainly DeltaH(cal). The determined crystal structure of VcEndA N69Q revealed the loss of one hydrogen bond present in native VcEndA, but also the formation of a new hydrogen bond involving residue 69 that could possibly explain the similar T(m) values for native and N69Q-mutated VcEndA. Structural analysis suggested that the stability, catalytic efficiency and salt tolerance of EndA were controlled by small changes in the hydrogen bonding networks and surface electrostatic potential. Our results indicate that endonuclease I adaptation is closely coupled to the conditions of the habitats of natural Vibrio, with VsEndA displaying a remarkable salt tolerance unique amongst the endonucleases characterized so far.

Effects of Probiotic Supplementation on the Gut Microbiota and Antibiotic Resistome Development in Preterm Infants.

Objectives: In 2014 probiotic supplementation (Lactobacillus acidophilus and Bifidobacterium longum subspecies infantis; InfloranⓇ) was introduced as standard of care to prevent necrotizing enterocolitis (NEC) in extremely preterm infants in Norway. We aimed to evaluate the influence of probiotics and antibiotic therapy on the developing gut microbiota and antibiotic resistome in extremely preterm infants, and to compare with very preterm infants and term infants not given probiotics. Study design: A prospective, observational multicenter study in six tertiary-care neonatal units. We enrolled 76 infants; 31 probiotic-supplemented extremely preterm infants <28 weeks gestation, 35 very preterm infants 28-31 weeks gestation not given probiotics and 10 healthy full-term control infants. Taxonomic composition and collection of antibiotic resistance genes (resistome) in fecal samples, collected at 7 and 28 days and 4 months age, were analyzed using shotgun-metagenome sequencing. Results: Median (IQR) birth weight was 835 (680-945) g and 1,290 (1,150-1,445) g in preterm infants exposed and not exposed to probiotics, respectively. Two extremely preterm infants receiving probiotic developed NEC requiring surgery. At 7 days of age we found higher median relative abundance of Bifidobacterium in probiotic supplemented infants (64.7%) compared to non-supplemented preterm infants (0.0%) and term control infants (43.9%). Lactobacillus was only detected in small amounts in all groups, but the relative abundance increased up to 4 months. Extremely preterm infants receiving probiotics had also much higher antibiotic exposure, still overall microbial diversity and resistome was not different than in more mature infants at 4 weeks and 4 months. Conclusion: Probiotic supplementation may induce colonization resistance and alleviate harmful effects of antibiotics on the gut microbiota and antibiotic resistome. Clinical Trial Registration: Clinicaltrials.gov: NCT02197468. https://clinicaltrials.gov/ct2/show/NCT02197468.

Comparative studies of endonuclease I from cold-adapted Vibrio salmonicida and mesophilic Vibrio cholerae.

Endonuclease I is a periplasmic or extracellular enzyme present in many different Proteobacteria. The endA gene encoding endonuclease I from the psychrophilic and mildly halophilic bacterium Vibrio salmonicida and from the mesophilic brackish water bacterium Vibrio cholerae have been cloned, over-expressed in Escherichia coli, and purified. A comparison of the enzymatic properties shows large differences in NaCl requirements, optimum pH, temperature stability and catalytic efficiency of the two proteins. The V. salmonicida EndA shows typical cold-adapted features such as lower unfolding temperature, lower temperature optimum for activity, and higher specific activity than V. cholerae EndA. The thermodynamic activation parameters confirm the psychrophilic nature of V. salmonicida EndA with a much lower activation enthalpy. The optimal conditions for enzymatic activity coincide well with the corresponding optimal requirements for growth of the organisms, and the enzymes function predominantly as DNases at physiological concentrations of NaCl. The periplasmic or extracellular localization of the enzymes, which renders them constantly exposed to the outer environment of the cell, may explain this fine-tuning of biochemical properties.

Sequence comparison and environmental adaptation of a bacterial endonuclease.

The periplasmic/extracellular bacterial enzyme endonuclease I was chosen as a model system to identify features that might be responsible for temperature- and salt adaptation. A statistical study of amino acid sequence properties belonging to endonuclease I enzymes from three mesophilic habitats (non-marine, brackish water and marine), and three marine temperature groups (psychrophile, intermediate and mesophile) has been conducted. Ten new endonuclease I genes have been sequenced in order to increase the sample size. A bioinformatical method of property dependent statistical analysis of alignments has been applied. To our knowledge this is the first time these methods have been used in order to investigate environmental adaptation of enzymes. Adaptation to low temperature seems to involve increased surface isoelectric point and hydrophobicity in contrast to salt adaptation in which the isoelectric point and hydrophobicity at the surface decreases. Redistribution of charge and hydrophobicity might be the most important signature for cold adaptation and salt adaptation of this enzyme class. The results indicate that general trends of adaptation are possible to elucidate from the amino acid sequences. Also in this paper a new scale of stratified B-factors, derived from the Protein Data Bank, is presented.

IS elements in Aliivibrio salmonicida LFI1238: occurrence, variability and impact on adaptability.

Insertion sequence (IS) elements are short, self-replicating DNA sequences that are capable of efficiently spreading over the host genome. Possessing varied integration specificity IS elements are capable of the irreversible inactivation of genes, which diversifies the pool of intact genetic determinants in host populations. In the current study, we performed a complex analysis of IS elements (Vsa IS) in the previously sequenced genome of Aliivibrio salmonicida LFI1238 and proposed a model of the spread of the Vsa IS elements over the genome of this microorganism. Along with the prediction of the integration sites for Vsa IS elements, the current study provides an overview of the properties of A. salmonicida IS elements, as well as information regarding their occurrence in different bacterial classes. An analysis of individual alleles of the IS elements has allowed us to depict a history of the accumulation of mutations and to describe distinctive microevolution lines for actively transposing Vsa IS elements in the genome of A. salmonicida LFI1238. Our results demonstrate the high importance of the dead end microevolution of actively transposing Vsa IS elements for the inactivation of genes in A. salmonicida LFI1238.

Draft Genome Sequence of the Actinomycete Rhodococcus sp. Strain AW25M09, Isolated from the Hadsel Fjord, Northern Norway.

The cold-adapted Rhodococcus sp. strain AW25M09 was isolated from an Atlantic hagfish caught off the shore of northern Norway as part of an ongoing bioprospecting project that aims to identify novel bacteria with biotechnological potential. Here, we present the 5.8-Mb draft genome sequence, together with details regarding the origin of the strain and its sequence assembly.

Identification of type VI secretion systems in Moritella viscosa.

The study describes the identification of type VI secretion systems (T6SSs) in Moritella viscosa, the aetiological agent of winter ulcer disease. Despite the availability of commercial vaccines, M. viscosa causes significant financial losses in salmonid farming. The T6SS transports bacterial proteins from the cell into the environment or directly into host cells, and has been implicated with bacterial virulence. The aim of the study was to identify potential T6SSs in M. viscosa and to determine whether it possesses active T6S, providing further insight into the biology of the bacterium. The genome of M. viscosa 06/09/139 was screened for homology with known T6SS encoding genes. Two genetically distinct loci, termed Moritella Type Six Secretion 1 and 2 (mts1 and mts2), were identified as encoding putative T6SSs. Each locus contained known T6S core genes. The mts2 locus contained species specific genes, some of which have not previously been connected with T6S. The mts1 locus showed sequence homology and synteny to T6SSs of the fish pathogen Aliivibrio salmonicida and a non-pathogenic Moritella sp. PE36. The mts2 locus was more similar to a Vibrio parahaemolyticus T6SS. A functional T6SS was confirmed through identification of secreted Mts1-M, a hemolysin coregulated protein (Hcp) which is a part of the secretion system. Both virulent and avirulent M. viscosa isolates expressed two genes encoding Hcp, mts1-M and mts2-M. The results show that M. viscosa has a functional T6S, but the role of the secretion system and possible connections with virulence need further examination.

Structural adaptation of endonuclease I from the cold-adapted and halophilic bacterium Vibrio salmonicida.

The crystal structure of the periplasmic/extracellular endonuclease I from Vibrio salmonicida has been solved to 1.5 A resolution and, in comparison to the corresponding endonucleases from V. cholerae and V. vulnificus, serves as a model system for the investigation of the structural determinants involved in the temperature and NaCl adaptation of this enzyme class. The overall fold of the three enzymes is essentially similar, but the V. salmonicida endonuclease displays a significantly more positive surface potential than the other two enzymes owing to the presence of ten more Lys residues. However, if the optimum salt concentrations for the V. salmonicida and V. cholerae enzymes are taken into consideration in the electrostatic surface-potential calculation, the potentials of the two enzymes become surprisingly similar. The higher number of basic residues in the V. salmonicida protein is therefore likely to be a result, at least in part, of adaptation to the more saline habitat of V. salmonicida (seawater) than V. cholerae (brackish water). The hydrophobic core of all three enzymes is almost identical, but the V. salmonicida endonuclease has a slightly lower number of internal hydrogen bonds. This, together with repulsive forces between the basic residues on the protein surface of V. salmonicida endonuclease I and differences in the distribution of salt bridges, probably results in higher flexibility of regions of the V. salmonicida protein. This is likely to influence both the catalytic activity and the stability of the protein.

Molecular characterization of cold adaptation based on ortholog protein sequences from Vibrionaceae species.

A set of 298 protein families from psychrophilic Vibrio salmonicida was compiled to identify genotypic characteristics that discern it from orthologous sequences from the mesophilic Vibrio/Photobacterium branch of the gamma-Proteobacteria (Vibrionaceae family). In our comparative exploration we employed alignment based bioinformatical and statistical methods. Interesting information was found in the substitution matrices, and the pattern of asymmetries in the amino acid substitution process. Together with the compositional difference, they identified the amino acids Ile, Asn, Ala and Gln as those having the most psycrophilic involvement. Ile and Asn are enhanced whereas Gln and Ala are suppressed. The inflexible Pro residue is also suppressed in loop regions, as expected in a flexible structure. The dataset were also classified and analysed according to the predicted subcellular location, and we made an additional study of 183 intracellular and 65 membrane proteins. Our results revealed that the psychrophilic proteins have similar hydrophobic and charge contributions in the core of the protein as mesophilic proteins, while the solvent-exposed surface area is significantly more hydrophobic. In addition, the psychrophilic intracellular (but not the membrane) proteins are significantly more negatively charged at the surface. Our analysis supports the hypothesis of preference for more flexible amino acids at the molecular surface. Life in cold climate seems to be obtained through many minor structural modifications rather than certain amino acids substitutions.

Comparative expression study to increase the solubility of cold adapted Vibrio proteins in Escherichia coli.

Functional and structural studies require gene overexpression and purification of soluble proteins. We wanted to express proteins from the psychrophilic bacterium Vibrio salmonicida in Escherichia coli, but encountered solubility problems. To improve the solubility of the proteins, we compared the effects of six N-terminal fusion proteins (Gb1, Z, thioredoxin, GST, MBP and NusA) and an N-terminal His6-tag. The selected test set included five proteins from the fish pathogen V. salmonicida and two related products from the mesophilic human pathogen Vibrio cholerae. We tested the expression in two different expression strains and at three different temperatures (16, 23 and 37 degrees C). His6-tag was the least effective tag, and these vector constructs were also difficult to transform. MBP and NusA performed best, expressing soluble proteins with all fusion partners in at least one of the cell types. In some cases MBP, GST and thioredoxin fusions resulted in products of incorrect size. The effect of temperature is complex: in most cases level of expression increased with temperature, whereas the effect on solubility was opposite. We found no clear connection between the preferred expression temperature of the protein and the temperature of the original host organism's natural habitat.

The structure of Vibrio cholerae extracellular endonuclease I reveals the presence of a buried chloride ion.

The crystal structure of a periplasmic/extracellular endonuclease from Vibrio cholerae has been solved at low and at neutral pH. Crystals grown at pH 4.6 and 6.9 diffracted to 1.6 A (on BM01A at the ESRF) and 1.95 A (on a rotating-anode generator), respectively. The structures of the endonuclease were compared with the structure of a homologous enzyme in V. vulnificus. The structures of the V. cholerae enzyme at different pH values are essentially identical to each other and to the V. vulnificus enzyme. However, interesting features were observed in the solvent structures. Both V. cholerae structures reveal the presence of a chloride ion completely buried within the core of the protein, with the nearest solvent molecule approximately 7 A away. Magnesium, which is essential for catalysis, is present in the structure at neutral pH, but is absent at low pH, and may partly explain the inactivity of the enzyme at lower pH.

The structure of uracil-DNA glycosylase from Atlantic cod (Gadus morhua) reveals cold-adaptation features.

Uracil-DNA glycosylase (UDG; EC 3.2.2.3) is a DNA-repair protein that catalyses the hydrolysis of promutagenic uracil residues from single- or double-stranded DNA, generating free uracil and abasic DNA. The crystal structure of the catalytic domain of cod uracil-DNA glycosylase (cUDG) has been determined to 1.9 A resolution, with final R factors of 18.61 and 20.57% for the working and test sets of reflections, respectively. This is the first crystal structure of a uracil-DNA glycosylase from a cold-adapted species and a detailed comparison with the human enzyme is performed in order to rationalize the cold-adapted behaviour of the cod enzyme at the structural level. The catalytic domain of cUDG comprises 223 residues, with a sequence identity to the human UDG of 75%. The tertiary structures of the two enzymes are also similar, with an overall displacement in main-chain atomic positions of 0.63 A. The amino-acid substitutions and the differences in intramolecular hydrogen bonds, hydrophobic interactions, ion-pair interactions and electrostatic potentials are compared and discussed in order to gain insight into the factors that cause the increased activity and reduced thermostability of the cod enzyme. In particular, the reduced number of strong ion-pair interactions in the C-terminal half of cUDG is believed to greatly affect the flexibility and/or stability. Increased positive electrostatic surface potential on the DNA-facing side of cUDG seems to be responsible for increasing the affinity for the negatively charged DNA compared with that of hUDG.