Assembly and Comparison of Ca. Neoehrlichia mikurensis Genomes.

Ca. Neoehrlichia mikurensis is widely prevalent in I. ricinus across Europe and has been associated with human disease. However, diagnostic modalities are limited, and much is still unknown about its biology. Here, we present the first complete Ca. Neoehrlichia mikurensis genomes directly derived from wildlife reservoir host tissues, using both long- and short-read sequencing technologies. This pragmatic approach provides an alternative to obtaining sufficient material from clinical cases, a difficult task for emerging infectious diseases, and to expensive and challenging bacterial isolation and culture methods. Both genomes exhibit a larger chromosome than the currently available Ca. Neoehrlichia mikurensis genomes and expand the ability to find new targets for the development of supportive laboratory diagnostics in the future. Moreover, this method could be utilized for other tick-borne pathogens that are difficult to culture.

SALARECON connects the Atlantic salmon genome to growth and feed efficiency.

Atlantic salmon (Salmo salar) is the most valuable farmed fish globally and there is much interest in optimizing its genetics and rearing conditions for growth and feed efficiency. Marine feed ingredients must be replaced to meet global demand, with challenges for fish health and sustainability. Metabolic models can address this by connecting genomes to metabolism, which converts nutrients in the feed to energy and biomass, but such models are currently not available for major aquaculture species such as salmon. We present SALARECON, a model focusing on energy, amino acid, and nucleotide metabolism that links the Atlantic salmon genome to metabolic fluxes and growth. It performs well in standardized tests and captures expected metabolic (in)capabilities. We show that it can explain observed hypoxic growth in terms of metabolic fluxes and apply it to aquaculture by simulating growth with commercial feed ingredients. Predicted limiting amino acids and feed efficiencies agree with data, and the model suggests that marine feed efficiency can be achieved by supplementing a few amino acids to plant- and insect-based feeds. SALARECON is a high-quality model that makes it possible to simulate Atlantic salmon metabolism and growth. It can be used to explain Atlantic salmon physiology and address key challenges in aquaculture such as development of sustainable feeds.

FAIR data station for lightweight metadata management and validation of omics studies.

BACKGROUND: The life sciences are one of the biggest suppliers of scientific data. Reusing and connecting these data can uncover hidden insights and lead to new concepts. Efficient reuse of these datasets is strongly promoted when they are interlinked with a sufficient amount of machine-actionable metadata. While the FAIR (Findable, Accessible, Interoperable, Reusable) guiding principles have been accepted by all stakeholders, in practice, there are only a limited number of easy-to-adopt implementations available that fulfill the needs of data producers. FINDINGS: We developed the FAIR Data Station, a lightweight application written in Java, that aims to support researchers in managing research metadata according to the FAIR principles. It implements the ISA metadata framework and uses minimal information metadata standards to capture experiment metadata. The FAIR Data Station consists of 3 modules. Based on the minimal information model(s) selected by the user, the "form generation module" creates a metadata template Excel workbook with a header row of machine-actionable attribute names. The Excel workbook is subsequently used by the data producer(s) as a familiar environment for sample metadata registration. At any point during this process, the format of the recorded values can be checked using the "validation module." Finally, the "resource module" can be used to convert the set of metadata recorded in the Excel workbook in RDF format, enabling (cross-project) (meta)data searches and, for publishing of sequence data, in an European Nucleotide Archive-compatible XML metadata file. CONCLUSIONS: Turning FAIR into reality requires the availability of easy-to-adopt data FAIRification workflows that are also of direct use for data producers. As such, the FAIR Data Station provides, in addition to the means to correctly FAIRify (omics) data, the means to build searchable metadata databases of similar projects and can assist in ENA metadata submission of sequence data. The FAIR Data Station is available at https://fairbydesign.nl.

Machine learning approaches to predict the Plant-associated phenotype of Xanthomonas strains.

BACKGROUND: The genus Xanthomonas has long been considered to consist predominantly of plant pathogens, but over the last decade there has been an increasing number of reports on non-pathogenic and endophytic members. As Xanthomonas species are prevalent pathogens on a wide variety of important crops around the world, there is a need to distinguish between these plant-associated phenotypes. To date a large number of Xanthomonas genomes have been sequenced, which enables the application of machine learning (ML) approaches on the genome content to predict this phenotype. Until now such approaches to the pathogenomics of Xanthomonas strains have been hampered by the fragmentation of information regarding pathogenicity of individual strains over many studies. Unification of this information into a single resource was therefore considered to be an essential step. RESULTS: Mining of 39 papers considering both plant-associated phenotypes, allowed for a phenotypic classification of 578 Xanthomonas strains. For 65 plant-pathogenic and 53 non-pathogenic strains the corresponding genomes were available and de novo annotated for the presence of Pfam protein domains used as features to train and compare three ML classification algorithms; CART, Lasso and Random Forest. CONCLUSION: The literature resource in combination with recursive feature extraction used in the ML classification algorithms provided further insights into the virulence enabling factors, but also highlighted domains linked to traits not present in pathogenic strains.

Genomic convergence between Akkermansia muciniphila in different mammalian hosts.

BACKGROUND: Akkermansia muciniphila is a member of the human gut microbiota where it resides in the mucus layer and uses mucin as the sole carbon, nitrogen and energy source. A. muciniphila is the only representative of the Verrucomicrobia phylum in the human gut. However, A. muciniphila 16S rRNA gene sequences have also been found in the intestines of many vertebrates. RESULTS: We detected A. muciniphila-like bacteria in the intestines of animals belonging to 15 out of 16 mammalian orders. In addition, other species belonging to the Verrucomicrobia phylum were detected in fecal samples. We isolated 10 new A. muciniphila strains from the feces of chimpanzee, siamang, mouse, pig, reindeer, horse and elephant. The physiology and genome of these strains were highly similar in comparison to the type strain A. muciniphila MucT. Overall, the genomes of the new strains showed high average nucleotide identity (93.9 to 99.7%). In these genomes, we detected considerable conservation of at least 75 of the 78 mucin degradation genes that were previously detected in the genome of the type strain MucT. CONCLUSIONS: The low genomic divergence observed in the new strains may indicate that A. muciniphila favors mucosal colonization independent of the differences in hosts. In addition, the conserved mucus degradation capability points towards a similar beneficial role of the new strains in regulating host metabolic health.

Galactocerebroside biosynthesis pathways of Mycoplasma species: an antigen triggering Guillain-Barré-Stohl syndrome.

Infection by Mycoplasma pneumoniae has been identified as a preceding factor of Guillain-Barré-Stohl syndrome. The Guillain-Barré-Stohl syndrome is triggered by an immune reaction against the major glycolipids and it has been postulated that M. pneumoniae infection triggers this syndrome due to bacterial production of galactocerebroside. Here, we present an extensive comparison of 224 genome sequences from 104 Mycoplasma species to characterize the genetic determinants of galactocerebroside biosynthesis. Hidden Markov models were used to analyse glycosil transferases, leading to identification of a functional protein domain, termed M2000535 that appears in about a third of the studied genomes. This domain appears to be associated with a potential UDP-glucose epimerase, which converts UDP-glucose into UDP-galactose, a main substrate for the biosynthesis of galactocerebroside. These findings clarify the pathogenic mechanisms underlining the triggering of Guillain-Barré-Stohl syndrome by M. pneumoniae infections.

A chromosome-level assembly of the black tiger shrimp (Penaeus monodon) genome facilitates the identification of growth-associated genes.

To salvage marine ecosystems from fishery overexploitation, sustainable and efficient aquaculture must be emphasized. The knowledge obtained from available genome sequence of marine organisms has accelerated marine aquaculture in many cases. The black tiger shrimp (Penaeus monodon) is one of the most prominent cultured penaeid shrimps (Crustacean) with an average annual global production of half a million tons in the last decade. However, its currently available genome assemblies lack the contiguity and completeness required for accurate genome annotation due to the highly repetitive nature of the genome and technical difficulty in extracting high-quality, high-molecular weight DNA. Here, we report the first chromosome-level whole-genome assembly of P. monodon. The combination of long-read Pacific Biosciences (PacBio) and long-range Chicago and Hi-C technologies enabled a successful assembly of this first high-quality genome sequence. The final assembly covered 2.39 Gb (92.3% of the estimated genome size) and contained 44 pseudomolecules, corresponding to the haploid chromosome number. Repetitive elements occupied a substantial portion of the assembly (62.5%), the highest of the figures reported among crustacean species. The availability of this high-quality genome assembly enabled the identification of genes associated with rapid growth in the black tiger shrimp through the comparison of hepatopancreas transcriptome of slow-growing and fast-growing shrimps. The results highlighted several growth-associated genes. Our high-quality genome assembly provides an invaluable resource for genetic improvement and breeding penaeid shrimp in aquaculture. The availability of P. monodon genome enables analyses of ecological impact, environment adaptation and evolution, as well as the role of the genome to protect the ecological resources by promoting sustainable shrimp farming.

A metabolic and physiological design study of Pseudomonas putida KT2440 capable of anaerobic respiration.

BACKGROUND: Pseudomonas putida KT2440 is a metabolically versatile, HV1-certified, genetically accessible, and thus interesting microbial chassis for biotechnological applications. However, its obligate aerobic nature hampers production of oxygen sensitive products and drives up costs in large scale fermentation. The inability to perform anaerobic fermentation has been attributed to insufficient ATP production and an inability to produce pyrimidines under these conditions. Addressing these bottlenecks enabled growth under micro-oxic conditions but does not lead to growth or survival under anoxic conditions. RESULTS: Here, a data-driven approach was used to develop a rational design for a P. putida KT2440 derivative strain capable of anaerobic respiration. To come to the design, data derived from a genome comparison of 1628 Pseudomonas strains was combined with genome-scale metabolic modelling simulations and a transcriptome dataset of 47 samples representing 14 environmental conditions from the facultative anaerobe Pseudomonas aeruginosa. CONCLUSIONS: The results indicate that the implementation of anaerobic respiration in P. putida KT2440 would require at least 49 additional genes of known function, at least 8 genes encoding proteins of unknown function, and 3 externally added vitamins.

A protocol for adding knowledge to Wikidata: aligning resources on human coronaviruses.

BACKGROUND: Pandemics, even more than other medical problems, require swift integration of knowledge. When caused by a new virus, understanding the underlying biology may help finding solutions. In a setting where there are a large number of loosely related projects and initiatives, we need common ground, also known as a "commons." Wikidata, a public knowledge graph aligned with Wikipedia, is such a commons and uses unique identifiers to link knowledge in other knowledge bases. However, Wikidata may not always have the right schema for the urgent questions. In this paper, we address this problem by showing how a data schema required for the integration can be modeled with entity schemas represented by Shape Expressions. RESULTS: As a telling example, we describe the process of aligning resources on the genomes and proteomes of the SARS-CoV-2 virus and related viruses as well as how Shape Expressions can be defined for Wikidata to model the knowledge, helping others studying the SARS-CoV-2 pandemic. How this model can be used to make data between various resources interoperable is demonstrated by integrating data from NCBI (National Center for Biotechnology Information) Taxonomy, NCBI Genes, UniProt, and WikiPathways. Based on that model, a set of automated applications or bots were written for regular updates of these sources in Wikidata and added to a platform for automatically running these updates. CONCLUSIONS: Although this workflow is developed and applied in the context of the COVID-19 pandemic, to demonstrate its broader applicability it was also applied to other human coronaviruses (MERS, SARS, human coronavirus NL63, human coronavirus 229E, human coronavirus HKU1, human coronavirus OC4).

Publisher Correction: MEMOTE for standardized genome-scale metabolic model testing.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Organohalide-respiring Desulfoluna species isolated from marine environments.

The genus Desulfoluna comprises two anaerobic sulfate-reducing strains, D. spongiiphila AA1T and D. butyratoxydans MSL71T, of which only the former was shown to perform organohalide respiration (OHR). Here we isolated a third strain, designated D. spongiiphila strain DBB, from marine intertidal sediment using 1,4-dibromobenzene and sulfate as the electron acceptors and lactate as the electron donor. Each strain harbors three reductive dehalogenase gene clusters (rdhABC) and corrinoid biosynthesis genes in their genomes, and dehalogenated brominated but not chlorinated organohalogens. The Desulfoluna strains maintained OHR in the presence of 20 mM sulfate or 20 mM sulfide, which often negatively affect other organohalide-respiring bacteria. Strain DBB sustained OHR with 2% oxygen in the gas phase, in line with its genetic potential for reactive oxygen species detoxification. Reverse transcription-quantitative PCR revealed differential induction of rdhA genes in strain DBB in response to 1,4-dibromobenzene or 2,6-dibromophenol. Proteomic analysis confirmed expression of rdhA1 with 1,4-dibromobenzene, and revealed a partially shared electron transport chain from lactate to 1,4-dibromobenzene and sulfate, which may explain accelerated OHR during concurrent sulfate reduction. Versatility in using electron donors, de novo corrinoid biosynthesis, resistance to sulfate, sulfide and oxygen, and concurrent sulfate reduction and OHR may confer an advantage to marine Desulfoluna strains.

Metabolic shift induced by synthetic co-cultivation promotes high yield of chain elongated acids from syngas.

Bio-catalytic processes for sustainable production of chemicals and fuels receive increased attention within the concept of circular economy. Strategies to improve these production processes include genetic engineering of bio-catalysts or process technological optimization. Alternatively, synthetic microbial co-cultures can be used to enhance production of chemicals of interest. It remains often unclear however how microbe to microbe interactions affect the overall production process and how this can be further exploited for application. In the present study we explored the microbial interaction in a synthetic co-culture of Clostridium autoethanogenum and Clostridium kluyveri, producing chain elongated products from carbon monoxide. Monocultures of C. autoethanogenum converted CO to acetate and traces of ethanol, while during co-cultivation with C. kluyveri, it shifted its metabolism significantly towards solventogenesis. In C. autoethanogenum, expression of the genes involved in the central carbon- and energy-metabolism remained unchanged during co-cultivation compared to monoculture condition. Therefore the shift in the metabolic flux of C. autoethanogenum appears to be regulated by thermodynamics, and results from the continuous removal of ethanol by C. kluyveri. This trait could be further exploited, driving the metabolism of C. autoethanogenum to solely ethanol formation during co-cultivation, resulting in a high yield of chain elongated products from CO-derived electrons. This research highlights the important role of thermodynamic interactions in (synthetic) mixed microbial communities and shows that this can be exploited to promote desired conversions.

The Empusa code generator and its application to GBOL, an extendable ontology for genome annotation.

The RDF data model facilitates integration of diverse data available in structured and semi-structured formats. To obtain a coherent RDF graph the chosen ontology must be consistently applied. However, addition of new diverse data causes the ontology to evolve, which could lead to accumulation of unintended erroneous composites. Thus, there is a need for a gate keeping system that compares the intended content described in the ontology with the actual content of the resource. The Empusa code generator facilitates creation of composite RDF resources from disparate sources. Empusa can convert a schema into an associated application programming interface (API), that can be used to perform data consistency checks and generates Markdown documentation to make persistent URLs resolvable. Using Empusa consistency is ensured within and between the ontology and the content of the resource. As an illustration of the potential of Empusa, we present the Genome Biology Ontology Language (GBOL). GBOL uses and extends current ontologies to provide a formal representation of genomic entities, along with their properties, relations and provenance.

Diversity of tryptophan halogenases in sponges of the genus Aplysina.

Marine sponges are a prolific source of novel enzymes with promising biotechnological potential. Especially halogenases, which are key enzymes in the biosynthesis of brominated and chlorinated secondary metabolites, possess interesting properties towards the production of pharmaceuticals that are often halogenated. In this study we used a polymerase chain reaction (PCR)-based screening to simultaneously examine and compare the richness and diversity of putative tryptophan halogenase protein sequences and bacterial community structures of six Aplysina species from the Mediterranean and Caribbean seas. At the phylum level, bacterial community composition was similar amongst all investigated species and predominated by Actinobacteria, Chloroflexi, Cyanobacteria, Gemmatimonadetes, and Proteobacteria. We detected four phylogenetically diverse clades of putative tryptophan halogenase protein sequences, which were only distantly related to previously reported halogenases. The Mediterranean species Aplysina aerophoba harbored unique halogenase sequences, of which the most predominant was related to a sponge-associated Psychrobacter-derived sequence. In contrast, the Caribbean species shared numerous novel halogenase sequence variants and exhibited a highly similar bacterial community composition at the operational taxonomic unit (OTU) level. Correlations of relative abundances of halogenases with those of bacterial taxa suggest that prominent sponge symbiotic bacteria, including Chloroflexi and Actinobacteria, are putative producers of the detected enzymes and may thus contribute to the chemical defense of their host.

NG-Tax 2.0: A Semantic Framework for High-Throughput Amplicon Analysis.

NG-Tax 2.0 is a semantic framework for FAIR high-throughput analysis and classification of marker gene amplicon sequences including bacterial and archaeal 16S ribosomal RNA (rRNA), eukaryotic 18S rRNA and ribosomal intergenic transcribed spacer sequences. It can directly use single or merged reads, paired-end reads and unmerged paired-end reads from long range fragments as input to generate de novo amplicon sequence variants (ASV). Using the RDF data model, ASV's can be automatically stored in a graph database as objects that link ASV sequences with the full data-wise and element-wise provenance, thereby achieving the level of interoperability required to utilize such data to its full potential. The graph database can be directly queried, allowing for comparative analyses of over thousands of samples and is connected with an interactive Rshiny toolbox for analysis and visualization of (meta) data. Additionally, NG-Tax 2.0 exports an extended BIOM 1.0 (JSON) file as starting point for further analyses by other means. The extended BIOM file contains new attribute types to include information about the command arguments used, the sequences of the ASVs formed, classification confidence scores and is backwards compatible. The performance of NG-Tax 2.0 was compared with DADA2, using the plugin in the QIIME 2 analysis pipeline. Fourteen 16S rRNA gene amplicon mock community samples were obtained from the literature and evaluated. Precision of NG-Tax 2.0 was significantly higher with an average of 0.95 vs 0.58 for QIIME2-DADA2 while recall was comparable with an average of 0.85 and 0.77, respectively. NG-Tax 2.0 is written in Java. The code, the ontology, a Galaxy platform implementation, the analysis toolbox, tutorials and example SPARQL queries are freely available at http://wurssb.gitlab.io/ngtax under the MIT License.

Comparative Genomics Highlights Symbiotic Capacities and High Metabolic Flexibility of the Marine Genus Pseudovibrio.

Pseudovibrio is a marine bacterial genus members of which are predominantly isolated from sessile marine animals, and particularly sponges. It has been hypothesized that Pseudovibrio spp. form mutualistic relationships with their hosts. Here, we studied Pseudovibrio phylogeny and genetic adaptations that may play a role in host colonization by comparative genomics of 31 Pseudovibrio strains, including 25 sponge isolates. All genomes were highly similar in terms of encoded core metabolic pathways, albeit with substantial differences in overall gene content. Based on gene composition, Pseudovibrio spp. clustered by geographic region, indicating geographic speciation. Furthermore, the fact that isolates from the Mediterranean Sea clustered by sponge species suggested host-specific adaptation or colonization. Genome analyses suggest that Pseudovibrio hongkongensis UST20140214-015BT is only distantly related to other Pseudovibrio spp., thereby challenging its status as typical Pseudovibrio member. All Pseudovibrio genomes were found to encode numerous proteins with SEL1 and tetratricopeptide repeats, which have been suggested to play a role in host colonization. For evasion of the host immune system, Pseudovibrio spp. may depend on type III, IV, and VI secretion systems that can inject effector molecules into eukaryotic cells. Furthermore, Pseudovibrio genomes carry on average seven secondary metabolite biosynthesis clusters, reinforcing the role of Pseudovibrio spp. as potential producers of novel bioactive compounds. Tropodithietic acid, bacteriocin, and terpene biosynthesis clusters were highly conserved within the genus, suggesting an essential role in survival, for example through growth inhibition of bacterial competitors. Taken together, these results support the hypothesis that Pseudovibrio spp. have mutualistic relations with sponges.

SAPP: functional genome annotation and analysis through a semantic framework using FAIR principles.

Summary: To unlock the full potential of genome data and to enhance data interoperability and reusability of genome annotations we have developed SAPP, a Semantic Annotation Platform with Provenance. SAPP is designed as an infrastructure supporting FAIR de novo computational genomics but can also be used to process and analyze existing genome annotations. SAPP automatically predicts, tracks and stores structural and functional annotations and associated dataset- and element-wise provenance in a Linked Data format, thereby enabling information mining and retrieval with Semantic Web technologies. This greatly reduces the administrative burden of handling multiple analysis tools and versions thereof and facilitates multi-level large scale comparative analysis. Availability and implementation: SAPP is written in JAVA and freely available at https://gitlab.com/sapp and runs on Unix-like operating systems. The documentation, examples and a tutorial are available at https://sapp.gitlab.io. Contact: jasperkoehorst@gmail.com or peter.schaap@wur.nl.

Concurrent Haloalkanoate Degradation and Chlorate Reduction by Pseudomonas chloritidismutans AW-1T.

Haloalkanoates are environmental pollutants that can be degraded aerobically by microorganisms producing hydrolytic dehalogenases. However, there is a lack of information about the anaerobic degradation of haloalkanoates. Genome analysis of Pseudomonas chloritidismutans AW-1T, a facultative anaerobic chlorate-reducing bacterium, showed the presence of two putative haloacid dehalogenase genes, the l-DEX gene and dehI, encoding an l-2-haloacid dehalogenase (l-DEX) and a halocarboxylic acid dehydrogenase (DehI), respectively. Hence, we studied the concurrent degradation of haloalkanoates and chlorate as a yet-unexplored trait of strain AW-1T The deduced amino acid sequences of l-DEX and DehI revealed 33 to 37% and 26 to 86% identities with biochemically/structurally characterized l-DEX and the d- and dl-2-haloacid dehalogenase enzymes, respectively. Physiological experiments confirmed that strain AW-1T can grow on chloroacetate, bromoacetate, and both l- and d-α-halogenated propionates with chlorate as an electron acceptor. Interestingly, growth and haloalkanoate degradation were generally faster with chlorate as an electron acceptor than with oxygen as an electron acceptor. In line with this, analyses of l-DEX and DehI dehalogenase activities using cell-free extract (CFE) of strain AW-1T grown on dl-2-chloropropionate under chlorate-reducing conditions showed up to 3.5-fold higher dehalogenase activity than the CFE obtained from AW-1T cells grown on dl-2-chloropropionate under aerobic conditions. Reverse transcription-quantitative PCR showed that the l-DEX gene was expressed constitutively independently of the electron donor (haloalkanoates or acetate) or acceptor (chlorate or oxygen), whereas the expression of dehI was induced by haloalkanoates. Concurrent degradation of organic and inorganic halogenated compounds by strain AW-1T represents a unique metabolic capacity in a single bacterium, providing a new piece of the puzzle of the microbial halogen cycle.IMPORTANCE Halogenated organic and inorganic compounds are important environmental pollutants that have carcinogenic and genotoxic effects on both animals and humans. Previous research studied the degradation of organic and inorganic halogenated compounds separately but not concurrently. This study shows concurrent degradation of halogenated alkanoates and chlorate as an electron donor and acceptor, respectively, coupled to growth in a single bacterium, Pseudomonas chloritidismutans AW-1T Hence, besides biogenesis of molecular oxygen from chlorate reduction enabling a distinctive placement of strain AW-1T between aerobic and anaerobic microorganisms, we can now add another unique metabolic potential of this bacterium to the roster. The degradation of different halogenated compounds under anoxic conditions by a single bacterium is also of interest for the natural halogen cycle in different aquatic and terrestrial ecosystems where ample natural production of halogenated compounds has been documented.