Inferring and comparing metabolism across heterogeneous sets of annotated genomes using AuCoMe.

Comparative analysis of genome-scale metabolic networks (GSMNs) may yield important information on the biology, evolution, and adaptation of species. However, it is impeded by the high heterogeneity of the quality and completeness of structural and functional genome annotations, which may bias the results of such comparisons. To address this issue, we developed AuCoMe, a pipeline to automatically reconstruct homogeneous GSMNs from a heterogeneous set of annotated genomes without discarding available manual annotations. We tested AuCoMe with three data sets, one bacterial, one fungal, and one algal, and showed that it successfully reduces technical biases while capturing the metabolic specificities of each organism. Our results also point out shared and divergent metabolic traits among evolutionarily distant algae, underlining the potential of AuCoMe to accelerate the broad exploration of metabolic evolution across the tree of life.

Comparative Transcriptome Analysis Provides Insights into Response of Ulva compressa to Fluctuating Salinity Conditions.

Ulva compressa, a green tide-forming species, can adapt to hypo-salinity conditions, such as estuaries and brackish lakes. To understand the underlying molecular mechanisms of hypo-salinity stress tolerance, transcriptome-wide gene expression profiles in U. compressa were created using digital gene expression profiles. The RNA-seq data were analyzed based on the comparison of differently expressed genes involved in specific pathways under hypo-salinity and recovery conditions. The up-regulation of genes in photosynthesis and glycolysis pathways may contribute to the recovery of photosynthesis and energy metabolism, which could provide sufficient energy for the tolerance under long-term hyposaline stress. Multiple strategies, such as ion transportation and osmolytes metabolism, were performed to maintain the osmotic homeostasis. Additionally, several long noncoding RNA were differently expressed during the stress, which could play important roles in the osmotolerance. Our work will serve as an essential foundation for the understanding of the tolerance mechanism of U. compressa under the fluctuating salinity conditions.

Computational analysis of the Plasmodiophora brassicae genome: mitochondrial sequence description and metabolic pathway database design.

Plasmodiophora brassicae is an obligate biotrophic pathogenic protist responsible for clubroot, a root gall disease of Brassicaceae species. In addition to the reference genome of the P. brassicae European e3 isolate and the draft genomes of Canadian or Chinese isolates, we present the genome of eH, a second European isolate. Refinement of the annotation of the eH genome led to the identification of the mitochondrial genome sequence, which was found to be bigger than that of Spongospora subterranea, another plant parasitic Plasmodiophorid phylogenetically related to P. brassicae. New pathways were also predicted, such as those for the synthesis of spermidine, a polyamine up-regulated in clubbed regions of roots. A P. brassicae pathway genome database was created to facilitate the functional study of metabolic pathways in transcriptomics approaches. These available tools can help in our understanding of the regulation of P. brassicae metabolism during infection and in response to diverse constraints.

Empirical evidence for metabolic drift in plant and algal lipid biosynthesis pathways.

Metabolic pathway drift has been formulated as a general principle to help in the interpretation of comparative analyses between biosynthesis pathways. Indeed, such analyses often indicate substantial differences, even in widespread pathways that are sometimes believed to be conserved. Here, our purpose is to check how much this interpretation fits to empirical data gathered in the field of plant and algal biosynthesis pathways. After examining several examples representative of the diversity of lipid biosynthesis pathways, we explain why it is important to compare closely related species to gain a better understanding of this phenomenon. Furthermore, this comparative approach brings us to the question of how much biotic interactions are responsible for shaping this metabolic plasticity. We end up introducing some model systems that may be promising for further exploration of this question.

Robustness analysis of metabolic predictions in algal microbial communities based on different annotation pipelines.

Animals, plants, and algae rely on symbiotic microorganisms for their development and functioning. Genome sequencing and genomic analyses of these microorganisms provide opportunities to construct metabolic networks and to analyze the metabolism of the symbiotic communities they constitute. Genome-scale metabolic network reconstructions rest on information gained from genome annotation. As there are multiple annotation pipelines available, the question arises to what extent differences in annotation pipelines impact outcomes of these analyses. Here, we compare five commonly used pipelines (Prokka, MaGe, IMG, DFAST, RAST) from predicted annotation features (coding sequences, Enzyme Commission numbers, hypothetical proteins) to the metabolic network-based analysis of symbiotic communities (biochemical reactions, producible compounds, and selection of minimal complementary bacterial communities). While Prokka and IMG produced the most extensive networks, RAST and DFAST networks produced the fewest false positives and the most connected networks with the fewest dead-end metabolites. Our results underline differences between the outputs of the tested pipelines at all examined levels, with small differences in the draft metabolic networks resulting in the selection of different microbial consortia to expand the metabolic capabilities of the algal host. However, the consortia generated yielded similar predicted producible compounds and could therefore be considered functionally interchangeable. This contrast between selected communities and community functions depending on the annotation pipeline needs to be taken into consideration when interpreting the results of metabolic complementarity analyses. In the future, experimental validation of bioinformatic predictions will likely be crucial to both evaluate and refine the pipelines and needs to be coupled with increased efforts to expand and improve annotations in reference databases.

Semi-Quantitative Targeted Gas Chromatography-Mass Spectrometry Profiling Supports a Late Side-Chain Reductase Cycloartenol-to-Cholesterol Biosynthesis Pathway in Brown Algae.

Sterols are biologically important molecules that serve as membrane fluidity regulators and precursors of signaling molecules, either endogenous or involved in biotic interactions. There is currently no model of their biosynthesis pathways in brown algae. Here, we benefit from the availability of genome data and gas chromatography-mass spectrometry (GC-MS) sterol profiling using a database of internal standards to build such a model. We expand the set of identified sterols in 11 species of red, brown, and green macroalgae and integrate these new data with genomic data. Our analyses suggest that some metabolic reactions may be conserved despite the loss of canonical eukaryotic enzymes, like the sterol side-chain reductase (SSR). Our findings are consistent with the principle of metabolic pathway drift through enzymatic replacement and show that cholesterol synthesis from cycloartenol may be a widespread but variable pathway among chlorophyllian eukaryotes. Among the factors contributing to this variability, one could be the recruitment of cholesterol biosynthetic intermediates to make signaling molecules, such as the mozukulins. These compounds were found in some brown algae belonging to Ectocarpales, and we here provide a first mozukulin biosynthetic model. Our results demonstrate that integrative approaches can already be used to infer experimentally testable models, which will be useful to further investigate the biological roles of those newly identified algal pathways.

Metage2Metabo, microbiota-scale metabolic complementarity for the identification of key species.

To capture the functional diversity of microbiota, one must identify metabolic functions and species of interest within hundreds or thousands of microorganisms. We present Metage2Metabo (M2M) a resource that meets the need for de novo functional screening of genome-scale metabolic networks (GSMNs) at the scale of a metagenome, and the identification of critical species with respect to metabolic cooperation. M2M comprises a flexible pipeline for the characterisation of individual metabolisms and collective metabolic complementarity. In addition, M2M identifies key species, that are meaningful members of the community for functions of interest. We demonstrate that M2M is applicable to collections of genomes as well as metagenome-assembled genomes, permits an efficient GSMN reconstruction with Pathway Tools, and assesses the cooperation potential between species. M2M identifies key organisms by reducing the complexity of a large-scale microbiota into minimal communities with equivalent properties, suitable for further analyses.

All the microbes that live in a specific environment, for example an organ, are collectively called the microbiota. In humans, the microbiota of the gut has been extensively studied by sequencing the DNA of the different microbes to identify them and determine the roles they play in health and disease. The DNA sequences of all the members of the microbiota is called the metagenome. The chemical reactions that the gut microbiota perform to produce energy and make the biomolecules they need to survive are collectively referred to as the metabolism of these microbes. Studying the metabolism of the gut microbiota can help researchers understand the roles of the different microbes. However, the large variety of species in the gut microbiota and gaps in the information about them render these studies difficult, despite technology improving quickly. To tackle this issue, Belcour, Frioux et al developed a new piece of software called Metage2Metabo (M2M) that simulates the metabolism of the gut microbiota and describes the metabolic relationships between the different microbes. Metage2Metabo analyses the roles of the metabolic genes of a large number of microbe species to establish how they complement each other metabolically. Then, it can calculate the minimum number of species needed to perform a metabolic role of interest within that microbiota, and which key species are associated with that role. To test the new software, Belcour, Frioux et al. used Metage2Metabo to analyse genomes from the human gut microbiota and from the cow rumen (one of the cow's stomachs). They showed that even if the metagenome was incomplete, the software was able to make stable predictions of key species involved in metabolic complementarity. Additionally, they also illustrated how the method can be used to study the gut microbiota of individuals. This work presents a new method for determining the metabolic relationships between species within a microbiota. The software is highly flexible and could be adapted to identify key members within different communities. In the context of the gut microbiota, the predictions of Metage2Metabo could shed lights on the interactions between the host and the microbes and contribute to a better understanding of microbe environments.

Soil microbiota influences clubroot disease by modulating Plasmodiophora brassicae and Brassica napus transcriptomes.

The contribution of surrounding plant microbiota to disease development has led to the 'pathobiome' concept, which represents the interaction between the pathogen, the host plant and the associated biotic microbial community, resulting or not in plant disease. The aim herein is to understand how the soil microbial environment may influence the functions of a pathogen and its pathogenesis, and the molecular response of the plant to the infection, with a dual-RNAseq transcriptomics approach. We address this question using Brassica napus and Plasmodiophora brassicae, the pathogen responsible for clubroot. A time-course experiment was conducted to study interactions between P. brassicae, two B. napus genotypes and three soils harbouring high, medium or low microbiota diversities and levels of richness. The soil microbial diversity levels had an impact on disease development (symptom levels and pathogen quantity). The P. brassicae and B. napus transcriptional patterns were modulated by these microbial diversities, these modulations being dependent on the host genotype plant and the kinetic time. The functional analysis of gene expressions allowed the identification of pathogen and plant host functions potentially involved in the change of plant disease level, such as pathogenicity-related genes (NUDIX effector) in P. brassicae and plant defence-related genes (glucosinolate metabolism) in B. napus.

Inferring Biochemical Reactions and Metabolite Structures to Understand Metabolic Pathway Drift.

Inferring genome-scale metabolic networks in emerging model organisms is challenged by incomplete biochemical knowledge and partial conservation of biochemical pathways during evolution. Therefore, specific bioinformatic tools are necessary to infer biochemical reactions and metabolic structures that can be checked experimentally. Using an integrative approach combining genomic and metabolomic data in the red algal model Chondrus crispus, we show that, even metabolic pathways considered as conserved, like sterols or mycosporine-like amino acid synthesis pathways, undergo substantial turnover. This phenomenon, here formally defined as "metabolic pathway drift," is consistent with findings from other areas of evolutionary biology, indicating that a given phenotype can be conserved even if the underlying molecular mechanisms are changing. We present a proof of concept with a methodological approach to formalize the logical reasoning necessary to infer reactions and molecular structures, abstracting molecular transformations based on previous biochemical knowledge.

Genome-Scale Metabolic Networks Shed Light on the Carotenoid Biosynthesis Pathway in the Brown Algae Saccharina japonica and Cladosiphon okamuranus.

Understanding growth mechanisms in brown algae is a current scientific and economic challenge that can benefit from the modeling of their metabolic networks. The sequencing of the genomes of Saccharina japonica and Cladosiphon okamuranus has provided the necessary data for the reconstruction of Genome-Scale Metabolic Networks (GSMNs). The same in silico method deployed for the GSMN reconstruction of Ectocarpus siliculosus to investigate the metabolic capabilities of these two algae, was used. Integrating metabolic profiling data from the literature, we provided functional GSMNs composed of an average of 2230 metabolites and 3370 reactions. Based on these GSMNs and previously published work, we propose a model for the biosynthetic pathways of the main carotenoids in these two algae. We highlight, on the one hand, the reactions and enzymes that have been preserved through evolution and, on the other hand, the specificities related to brown algae. Our data further indicate that, if abscisic acid is produced by Saccharina japonica, its biosynthesis pathway seems to be different in its final steps from that described in land plants. Thus, our work illustrates the potential of GSMNs reconstructions for formalizing hypotheses that can be further tested using targeted biochemical approaches.