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.

RNA targeting by the type III-A CRISPR-Cas Csm complex of Thermus thermophilus.

CRISPR-Cas is a prokaryotic adaptive immune system that provides sequence-specific defense against foreign nucleic acids. Here we report the structure and function of the effector complex of the Type III-A CRISPR-Cas system of Thermus thermophilus: the Csm complex (TtCsm). TtCsm is composed of five different protein subunits (Csm1-Csm5) with an uneven stoichiometry and a single crRNA of variable size (35-53 nt). The TtCsm crRNA content is similar to the Type III-B Cmr complex, indicating that crRNAs are shared among different subtypes. A negative stain EM structure of the TtCsm complex exhibits the characteristic architecture of Type I and Type III CRISPR-associated ribonucleoprotein complexes. crRNA-protein crosslinking studies show extensive contacts between the Csm3 backbone and the bound crRNA. We show that, like TtCmr, TtCsm cleaves complementary target RNAs at multiple sites. Unlike Type I complexes, interference by TtCsm does not proceed via initial base pairing by a seed sequence.

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).

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.

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.

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.

Structure and activity of the RNA-targeting Type III-B CRISPR-Cas complex of Thermus thermophilus.

The CRISPR-Cas system is a prokaryotic host defense system against genetic elements. The Type III-B CRISPR-Cas system of the bacterium Thermus thermophilus, the TtCmr complex, is composed of six different protein subunits (Cmr1-6) and one crRNA with a stoichiometry of Cmr112131445361:crRNA1. The TtCmr complex copurifies with crRNA species of 40 and 46 nt, originating from a distinct subset of CRISPR loci and spacers. The TtCmr complex cleaves the target RNA at multiple sites with 6 nt intervals via a 5' ruler mechanism. Electron microscopy revealed that the structure of TtCmr resembles a "sea worm" and is composed of a Cmr2-3 heterodimer "tail," a helical backbone of Cmr4 subunits capped by Cmr5 subunits, and a curled "head" containing Cmr1 and Cmr6. Despite having a backbone of only four Cmr4 subunits and being both longer and narrower, the overall architecture of TtCmr resembles that of Type I Cascade complexes.

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.

Reverse Methanogenesis and Respiration in Methanotrophic Archaea.

Anaerobic oxidation of methane (AOM) is catalyzed by anaerobic methane-oxidizing archaea (ANME) via a reverse and modified methanogenesis pathway. Methanogens can also reverse the methanogenesis pathway to oxidize methane, but only during net methane production (i.e., "trace methane oxidation"). In turn, ANME can produce methane, but only during net methane oxidation (i.e., enzymatic back flux). Net AOM is exergonic when coupled to an external electron acceptor such as sulfate (ANME-1, ANME-2abc, and ANME-3), nitrate (ANME-2d), or metal (oxides). In this review, the reversibility of the methanogenesis pathway and essential differences between ANME and methanogens are described by combining published information with domain based (meta)genome comparison of archaeal methanotrophs and selected archaea. These differences include abundances and special structure of methyl coenzyme M reductase and of multiheme cytochromes and the presence of menaquinones or methanophenazines. ANME-2a and ANME-2d can use electron acceptors other than sulfate or nitrate for AOM, respectively. Environmental studies suggest that ANME-2d are also involved in sulfate-dependent AOM. ANME-1 seem to use a different mechanism for disposal of electrons and possibly are less versatile in electron acceptors use than ANME-2. Future research will shed light on the molecular basis of reversal of the methanogenic pathway and electron transfer in different ANME types.

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.

Genome and proteome analysis of Pseudomonas chloritidismutans†AW-1T that grows on n-decane with chlorate or oxygen as electron acceptor.

Growth of Pseudomonas chloritidismutans AW-1T on C7 to C12 n-alkanes with oxygen or chlorate as electron acceptor was studied by genome and proteome analysis. Whole genome shotgun sequencing resulted in a 5 Mbp assembled sequence with a G + C content of 62.5%. The automatic annotation identified 4767 protein-encoding genes and a putative function could be assigned to almost 80% of the predicted proteins. The distinct phylogenetic position of P. chloritidismutans AW-1T within the Pseudomonas stutzeri cluster became clear by comparison of average nucleotide identity values of sequenced genomes. Analysis of the proteome of P. chloritidismutans AW-1T showed the versatility of this bacterium to adapt to aerobic and anaerobic growth conditions with acetate or n-decane as substrates. All enzymes involved in the alkane oxidation pathway were identified. An alkane monooxygenase was detected in n-decane-grown cells, but not in acetate-grown cells. The enzyme was found when grown in the presence of oxygen or chlorate, indicating that under both conditions an oxygenase-mediated pathway is employed for alkane degradation. Proteomic and biochemical data also showed that both chlorate reductase and chlorite dismutase are constitutively present, but most abundant under chlorate-reducing conditions.

A genomic view on syntrophic versus non-syntrophic lifestyle in anaerobic fatty acid degrading communities.

In sulfate-reducing and methanogenic environments complex biopolymers are hydrolyzed and degraded by fermentative micro-organisms that produce hydrogen, carbon dioxide and short chain fatty acids. Degradation of short chain fatty acids can be coupled to methanogenesis or to sulfate-reduction. Here we study from a genome perspective why some of these micro-organisms are able to grow in syntrophy with methanogens and others are not. Bacterial strains were selected based on genome availability and upon their ability to grow on short chain fatty acids alone or in syntrophic association with methanogens. Systematic functional domain profiling allowed us to shed light on this fundamental and ecologically important question. Extra-cytoplasmic formate dehydrogenases (InterPro domain number; IPR006443), including their maturation protein FdhE (IPR024064 and IPR006452) is a typical difference between syntrophic and non-syntrophic butyrate and propionate degraders. Furthermore, two domains with a currently unknown function seem to be associated with the ability of syntrophic growth. One is putatively involved in capsule or biofilm production (IPR019079) and a second in cell division, shape-determination or sporulation (IPR018365). The sulfate-reducing bacteria Desulfobacterium autotrophicum HRM2, Desulfomonile tiedjei and Desulfosporosinus meridiei were never tested for syntrophic growth, but all crucial domains were found in their genomes, which suggests their possible ability to grow in syntrophic association with methanogens. In addition, profiling domains involved in electron transfer mechanisms revealed the important role of the Rnf-complex and the formate transporter in syntrophy, and indicate that DUF224 may have a role in electron transfer in bacteria other than Syntrophomonas wolfei as well. This article is a part of a Special Issue entitled: 18th European Bioenergetics Conference (Biochim. Biophys. Acta, Volume 1837, Issue 7, July 2014).

Comparative genomics of Streptococcus pyogenes M1 isolates differing in virulence and propensity to cause systemic infection in mice.

Streptococcus pyogenes serotype M1 is a frequent cause of severe infections in humans. Some M1 isolates are pathogenic in mice and used in studies on infection pathogenesis. We observed marked differences in murine infections caused by M1 strain SF370, 5448, 5448AP or AP1 which prompted us to sequence the whole genome of isolates 5448 and AP1 for comparative analysis. Strain 5448 is known to acquire inactivating mutations in the CovRS two-component system during mouse infection, producing hypervirulent progeny such as 5448AP. Isolates AP1 and 5448AP, more than 5448, caused disseminating infections that became systemic and lethal. SF370 was not pathogenic. Phages caused gross genetic differences and increased the gene content of AP1 by 8% as compared to 5448 and SF370. Each of six examined M1 genomes contained two CRISPR-Cas systems. Phage insertion destroyed a type II CRISPR-Cas system in AP1 and other strains of serotypes M1, M3, M6 and M24, but not in M1 strains 5448, SF370, MGAS5005, A20 or M1 476. A resulting impaired defence against invading genetic elements could have led to the wealth of phages in AP1. AP1 lacks genetic features of the MGAS5005-like clonal complex including the streptodornase that drives selection for hypervirulent clones with inactivated CovRS system. Still, inactivating mutations in covS were a common genetic feature of AP1 and the MGAS5005-like isolate 5448AP. Abolished expression of the cysteine proteinase SpeB, due to CovRS inactivation could be a common cause for hypervirulence of the two isolates. Moreover, an additional protein H-coding gene and a mutation in the regulator gene rofA distinguished AP1 form other M1 isolates. In conclusion, hypervirulence of S. pyogenes M1 in mice is not limited to the MGAS5005-like genotype.

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.

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.

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.

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.

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.