Effects of plant tissue permeability on invasion and population bottlenecks of a phytopathogen.

Pathogen genetic diversity varies in response to environmental changes. However, it remains unclear whether plant barriers to invasion could be considered a genetic bottleneck for phytopathogen populations. Here, we implement a barcoding approach to generate a pool of 90 isogenic and individually barcoded Ralstonia solanacearum strains. We used 90 of these strains to inoculate tomato plants with different degrees of physical permeability to invasion (intact roots, wounded roots and xylem inoculation) and quantify the phytopathogen population dynamics during invasion. Our results reveal that the permeability of plant roots impacts the degree of population bottleneck, genetic diversity, and composition of Ralstonia populations. We also find that selection is the main driver structuring pathogen populations when barriers to infection are less permeable, i.e., intact roots, the removal of root physical and immune barriers results in the predominance of stochasticity in population assembly. Taken together, our study suggests that plant root permeability constitutes a bottleneck for phytopathogen invasion and genetic diversity.

Insights into the metabolic specificities of pathogenic strains from the Ralstonia solanacearum species complex.

All the strains grouped under the species Ralstonia solanacearum represent a species complex responsible for many diseases on agricultural crops throughout the world. The strains have different lifestyles and host range. Here, we investigated whether specific metabolic pathways contribute to strain diversification. To this end, we carried out systematic comparisons on 11 strains representing the diversity of the species complex. We reconstructed the metabolic network of each strain from its genome sequence and looked for the metabolic pathways differentiating the different reconstructed networks and, by extension, the different strains. Finally, we conducted an experimental validation by determining the metabolic profile of each strain with the Biolog technology. Results revealed that the metabolism is conserved between strains, with a core metabolism composed of 82% of the pan-reactome. The three species composing the species complex could be distinguished according to the presence/absence of some metabolic pathways, in particular, one involving salicylic acid degradation. Phenotypic assays revealed that the trophic preferences on organic acids and several amino acids such as glutamine, glutamate, aspartate, and asparagine are conserved between strains. Finally, we generated mutants lacking the quorum-sensing-dependent regulator PhcA in four diverse strains, and we showed that the phcA-dependent trade-off between growth and production of virulence factors is conserved across the R. solanacearum species complex. IMPORTANCE Ralstonia solanacearum is one of the most important threats to plant health worldwide, causing disease on a very large range of agricultural crops such as tomato or potato. Behind the R. solanacearum name are hundreds of strains with different host range and lifestyle, classified into three species. Studying the differences between strains allows to better apprehend the biology of the pathogens and the specificity of some strains. None of the published genomic comparative studies have focused on the metabolism of the strains so far. We developed a new bioinformatic pipeline to build high-quality metabolic networks and used a combination of metabolic modeling and high-throughput phenotypic Biolog microplates to look for the metabolic differences between 11 strains across the three species. Our study revealed that genes encoding enzymes are overall conserved, with few variations between strains. However, more variations were observed when considering substrate usage. These variations probably result from regulation rather than the presence or absence of enzymes in the genome.

Complete Genome Sequence of Sphingobium xenophagum PH3-15, Isolated from La Roche-Posay Thermal Water Sources.

We report here the complete genome sequence of Sphingobium xenophagum strain PH3-15, which was isolated from La Roche-Posay thermal water sources. The assembled 4.6-Mbp genome consisted of two chromosomes and three plasmids. These data will provide valuable information and important insights into the physiology and metabolism of this Sphingobium organism.

Robustness of plant quantitative disease resistance is provided by a decentralized immune network.

Quantitative disease resistance (QDR) represents the predominant form of resistance in natural populations and crops. Surprisingly, very limited information exists on the biomolecular network of the signaling machineries underlying this form of plant immunity. This lack of information may result from its complex and quantitative nature. Here, we used an integrative approach including genomics, network reconstruction, and mutational analysis to identify and validate molecular networks that control QDR in Arabidopsis thaliana in response to the bacterial pathogen Xanthomonas campestris To tackle this challenge, we first performed a transcriptomic analysis focused on the early stages of infection and using transgenic lines deregulated for the expression of RKS1, a gene underlying a QTL conferring quantitative and broad-spectrum resistance to XcampestrisRKS1-dependent gene expression was shown to involve multiple cellular activities (signaling, transport, and metabolism processes), mainly distinct from effector-triggered immunity (ETI) and pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) responses already characterized in Athaliana Protein-protein interaction network reconstitution then revealed a highly interconnected and distributed RKS1-dependent network, organized in five gene modules. Finally, knockout mutants for 41 genes belonging to the different functional modules of the network revealed that 76% of the genes and all gene modules participate partially in RKS1-mediated resistance. However, these functional modules exhibit differential robustness to genetic mutations, indicating that, within the decentralized structure of the QDR network, some modules are more resilient than others. In conclusion, our work sheds light on the complexity of QDR and provides comprehensive understanding of a QDR immune network.

Intercellular cooperation in a fungal plant pathogen facilitates host colonization.

Cooperation is associated with major transitions in evolution such as the emergence of multicellularity. It is central to the evolution of many complex traits in nature, including growth and virulence in pathogenic bacteria. Whether cells of multicellular parasites function cooperatively during infection remains, however, largely unknown. Here, we show that hyphal cells of the fungal pathogen Sclerotinia sclerotiorum reprogram toward division of labor to facilitate the colonization of host plants. Using global transcriptome sequencing, we reveal that gene expression patterns diverge markedly in cells at the center and apex of hyphae during Arabidopsis thaliana colonization compared with in vitro growth. We reconstructed a genome-scale metabolic model for S. sclerotiorum and used flux balance analysis to demonstrate metabolic heterogeneity supporting division of labor between hyphal cells. Accordingly, continuity between the central and apical compartments of invasive hyphae was required for optimal growth in planta Using a multicell model of fungal hyphae, we show that this cooperative functioning enhances fungal growth predominantly during host colonization. Our work identifies cooperation in fungal hyphae as a mechanism emerging at the multicellular level to support host colonization and virulence.

Comparative transcriptomic studies identify specific expression patterns of virulence factors under the control of the master regulator PhcA in the Ralstonia solanacearum species complex.

The global regulator PhcA controls numerous traits associated to virulence and bacterial proliferation in strains of the plant pathogen Ralstonia solanacearum species complex. Here, we conducted a genome-wide RNA sequencing study of the GMI1000 wild-type strain and a derived phcA mutant grown in complete medium. The PhcA regulon we identified is the largest regulon described to date in the R. solanacearum species complex with 1581 regulated genes, representing about 30% of the bacterial genome. Among these genes, 166 transcription regulators were identified including known regulators controlling major cellular functions such as the Type 3 secretion system and 27 novel regulators that were not identified in previous transcriptomic studies. This study highlights that PhcA controls other functions beside pathogenicity stricto sensu which participate to the global cell homeostasis (metabolism, energy storage). We then compared the PhcA regulon identified in complete medium to the recently published PhcA regulon obtained in planta. This comparison of the set of GMI1000 genes subjected to PhcA regulation in both conditions revealed 383 common genes. Among them, 326 (85%) had a similar PhcA dependent regulation pattern in complete medium and in planta, and 57 (15%) displayed an opposite regulation pattern. A large majority of the genes repressed by PhcA in complete medium but activated in planta belong to the HrpG-HrpB regulon, which represents a set of key genes required for R. solanacearum pathogenesis. This latter class of genes appears to be specifically induced by PhcA in the plant environment whereas PhcA represses their expression in complete medium. The large set of direct and indirect targets identified in this study will contribute to enrich our knowledge of the intricate regulatory network coordinating the expression of virulence and metabolic functions in the model plant pathogen R. solanacearum.

Control of primary metabolism by a virulence regulatory network promotes robustness in a plant pathogen.

Robustness is a key system-level property of living organisms to maintain their functions while tolerating perturbations. We investigate here how a regulatory network controlling multiple virulence factors impacts phenotypic robustness of a bacterial plant pathogen. We reconstruct a cell-scale model of Ralstonia solanacearum connecting a genome-scale metabolic network, a virulence macromolecule network, and a virulence regulatory network, which includes 63 regulatory components. We develop in silico methods to quantify phenotypic robustness under a broad set of conditions in high-throughput simulation analyses. This approach reveals that the virulence regulatory network exerts a control of the primary metabolism to promote robustness upon infection. The virulence regulatory network plugs into the primary metabolism mainly through the control of genes likely acquired via horizontal gene transfer, which results in a functional overlay with ancestral genes. These results support the view that robustness may be a selected trait that promotes pathogenic fitness upon infection.

Transposon Sequencing Uncovers an Essential Regulatory Function of Phosphoribulokinase for Methylotrophy.

Methylotrophy is the ability of organisms to grow at the expense of reduced one-carbon compounds, such as methanol or methane. Here, we used transposon sequencing combining hyper-saturated transposon mutagenesis with high-throughput sequencing to define the essential methylotrophy genome of Methylobacterium extorquens PA1, a model methylotroph. To distinguish genomic regions required for growth only on methanol from general required genes, we contrasted growth on methanol with growth on succinate, a non-methylotrophic reference substrate. About 500,000 insertions were mapped for each condition, resulting in a median insertion distance of five base pairs. We identified 147 genes and 76 genes as specific for growth on methanol and succinate, respectively, and a set of 590 genes as required under both growth conditions. For the integration of metabolic functions, we reconstructed a genome-scale metabolic model and performed in silico essentiality analysis. In total, the approach uncovered 95 genes not previously described as crucial for methylotrophy, including genes involved in respiration, carbon metabolism, transport, and regulation. Strikingly, regardless of the absence of the Calvin cycle in the methylotroph, the screen led to the identification of the gene for phosphoribulokinase as essential during growth on methanol, but not during growth on succinate. Genetic experiments in addition to metabolomics and proteomics revealed that phosphoribulokinase serves a key regulatory function. Our data support a model according to which ribulose-1,5-bisphosphate is an essential metabolite that induces a transcriptional regulator driving one-carbon assimilation.

Codon optimization underpins generalist parasitism in fungi.

The range of hosts that parasites can infect is a key determinant of the emergence and spread of disease. Yet, the impact of host range variation on the evolution of parasite genomes remains unknown. Here, we show that codon optimization underlies genome adaptation in broad host range parasites. We found that the longer proteins encoded by broad host range fungi likely increase natural selection on codon optimization in these species. Accordingly, codon optimization correlates with host range across the fungal kingdom. At the species level, biased patterns of synonymous substitutions underpin increased codon optimization in a generalist but not a specialist fungal pathogen. Virulence genes were consistently enriched in highly codon-optimized genes of generalist but not specialist species. We conclude that codon optimization is related to the capacity of parasites to colonize multiple hosts. Our results link genome evolution and translational regulation to the long-term persistence of generalist parasitism.

Exopolysaccharide Quantification for the Plant Pathogen Ralstonia solanacearum.

Soluble exopolysaccharide is a major virulence factor produced by the plant pathogen Ralstonia solanacearum. Its massive production during plant infection is associated with the arrest of water flow in xylem vessels leading eventually to plant death. The composition of this heavy macromolecule includes mainly N-acetylgalactosamine. Here we describe a colorimetric method for quantitative determination of the soluble exopolysaccharide present in culture supernatant of R. solanacearum.

Advances on plant-pathogen interactions from molecular toward systems biology perspectives.

In the past 2 decades, progress in molecular analyses of the plant immune system has revealed key elements of a complex response network. Current paradigms depict the interaction of pathogen-secreted molecules with host target molecules leading to the activation of multiple plant response pathways. Further research will be required to fully understand how these responses are integrated in space and time, and exploit this knowledge in agriculture. In this review, we highlight systems biology as a promising approach to reveal properties of molecular plant-pathogen interactions and predict the outcome of such interactions. We first illustrate a few key concepts in plant immunity with a network and systems biology perspective. Next, we present some basic principles of systems biology and show how they allow integrating multiomics data and predict cell phenotypes. We identify challenges for systems biology of plant-pathogen interactions, including the reconstruction of multiscale mechanistic models and the connection of host and pathogen models. Finally, we outline studies on resistance durability through the robustness of immune system networks, the identification of trade-offs between immunity and growth and in silico plant-pathogen co-evolution as exciting perspectives in the field. We conclude that the development of sophisticated models of plant diseases incorporating plant, pathogen and climate properties represent a major challenge for agriculture in the future.

Enhanced in planta Fitness through Adaptive Mutations in EfpR, a Dual Regulator of Virulence and Metabolic Functions in the Plant Pathogen Ralstonia solanacearum.

Experimental evolution of the plant pathogen Ralstonia solanacearum, where bacteria were maintained on plant lineages for more than 300 generations, revealed that several independent single mutations in the efpR gene from populations propagated on beans were associated with fitness gain on bean. In the present work, novel allelic efpR variants were isolated from populations propagated on other plant species, thus suggesting that mutations in efpR were not solely associated to a fitness gain on bean, but also on additional hosts. A transcriptomic profiling and phenotypic characterization of the efpR deleted mutant showed that EfpR acts as a global catabolic repressor, directly or indirectly down-regulating the expression of multiple metabolic pathways. EfpR also controls virulence traits such as exopolysaccharide production, swimming and twitching motilities and deletion of efpR leads to reduced virulence on tomato plants after soil drenching inoculation. We studied the impact of the single mutations that occurred in efpR during experimental evolution and found that these allelic mutants displayed phenotypic characteristics similar to the deletion mutant, although not behaving as complete loss-of-function mutants. These adaptive mutations therefore strongly affected the function of efpR, leading to an expanded metabolic versatility that should benefit to the evolved clones. Altogether, these results indicated that EfpR is a novel central player of the R. solanacearum virulence regulatory network. Independent mutations therefore appeared during experimental evolution in the evolved clones, on a crucial node of this network, to favor adaptation to host vascular tissues through regulatory and metabolic rewiring.

A Resource Allocation Trade-Off between Virulence and Proliferation Drives Metabolic Versatility in the Plant Pathogen Ralstonia solanacearum.

Bacterial pathogenicity relies on a proficient metabolism and there is increasing evidence that metabolic adaptation to exploit host resources is a key property of infectious organisms. In many cases, colonization by the pathogen also implies an intensive multiplication and the necessity to produce a large array of virulence factors, which may represent a significant cost for the pathogen. We describe here the existence of a resource allocation trade-off mechanism in the plant pathogen R. solanacearum. We generated a genome-scale reconstruction of the metabolic network of R. solanacearum, together with a macromolecule network module accounting for the production and secretion of hundreds of virulence determinants. By using a combination of constraint-based modeling and metabolic flux analyses, we quantified the metabolic cost for production of exopolysaccharides, which are critical for disease symptom production, and other virulence factors. We demonstrated that this trade-off between virulence factor production and bacterial proliferation is controlled by the quorum-sensing-dependent regulatory protein PhcA. A phcA mutant is avirulent but has a better growth rate than the wild-type strain. Moreover, a phcA mutant has an expanded metabolic versatility, being able to metabolize 17 substrates more than the wild-type. Model predictions indicate that metabolic pathways are optimally oriented towards proliferation in a phcA mutant and we show that this enhanced metabolic versatility in phcA mutants is to a large extent a consequence of not paying the cost for virulence. This analysis allowed identifying candidate metabolic substrates having a substantial impact on bacterial growth during infection. Interestingly, the substrates supporting well both production of virulence factors and growth are those found in higher amount within the plant host. These findings also provide an explanatory basis to the well-known emergence of avirulent variants in R. solanacearum populations in planta or in stressful environments.

Emerging Trends in Molecular Interactions between Plants and the Broad Host Range Fungal Pathogens Botrytis cinerea and Sclerotinia sclerotiorum.

Fungal plant pathogens are major threats to food security worldwide. Sclerotinia sclerotiorum and Botrytis cinerea are closely related Ascomycete plant pathogens causing mold diseases on hundreds of plant species. There is no genetic source of complete plant resistance to these broad host range pathogens known to date. Instead, natural plant populations show a continuum of resistance levels controlled by multiple genes, a phenotype designated as quantitative disease resistance. Little is known about the molecular mechanisms controlling the interaction between plants and S. sclerotiorum and B. cinerea but significant advances were made on this topic in the last years. This minireview highlights a selection of nine themes that emerged in recent research reports on the molecular bases of plant-S. sclerotiorum and plant-B. cinerea interactions. On the fungal side, this includes progress on understanding the role of oxalic acid, on the study of fungal small secreted proteins. Next, we discuss the exchanges of small RNA between organisms and the control of cell death in plant and fungi during pathogenic interactions. Finally on the plant side, we highlight defense priming by mechanical signals, the characterization of plant Receptor-like proteins and the hormone abscisic acid in the response to B. cinerea and S. sclerotiorum, the role of plant general transcription machinery and plant small bioactive peptides. These represent nine trends we selected as remarkable in our understanding of fungal molecules causing disease and plant mechanisms associated with disease resistance to two devastating broad host range fungi.

FlexFlux: combining metabolic flux and regulatory network analyses.

BACKGROUND: Expression of cell phenotypes highly depends on metabolism that supplies matter and energy. To achieve proper utilisation of the different metabolic pathways, metabolism is tightly regulated by a complex regulatory network composed of diverse biological entities (genes, transcripts, proteins, signalling molecules…). The integrated analysis of both regulatory and metabolic networks appears very insightful but is not straightforward because of the distinct characteristics of both networks. The classical method used for metabolic flux analysis is Flux Balance Analysis (FBA), which is constraint-based and relies on the assumption of steady-state metabolite concentrations throughout the network. Regarding regulatory networks, a broad spectrum of methods are dedicated to their analysis although logical modelling remains the major method to take charge of large-scale networks. RESULTS: We present FlexFlux, an application implementing a new way to combine the analysis of both metabolic and regulatory networks, based on simulations that do not require kinetic parameters and can be applied to genome-scale networks. FlexFlux is based on seeking regulatory network steady-states by performing synchronous updates of multi-state qualitative initial values. FlexFlux is then able to use the calculated steady-state values as constraints for metabolic flux analyses using FBA. As input, FlexFlux uses the standards Systems Biology Markup Language (SBML) and SBML Qualitative Models Package ("qual") extension (SBML-qual) file formats and provides a set of FBA based functions. CONCLUSIONS: FlexFlux is an open-source java software with executables and full documentation available online at http://lipm-bioinfo.toulouse.inra.fr/flexflux/. It can be defined as a research tool that enables a better understanding of both regulatory and metabolic networks based on steady-state simulations. FlexFlux integrates well in the flux analysis ecosystem thanks to the support of standard file formats and can thus be used as a complementary tool to existing software featuring other types of analyses.

Common protein sequence signatures associate with Sclerotinia borealis lifestyle and secretion in fungal pathogens of the Sclerotiniaceae.

Fungal plant pathogens produce secreted proteins adapted to function outside fungal cells to facilitate colonization of their hosts. In many cases such as for fungi from the Sclerotiniaceae family the repertoire and function of secreted proteins remains elusive. In the Sclerotiniaceae, whereas Sclerotinia sclerotiorum and Botrytis cinerea are cosmopolitan broad host-range plant pathogens, Sclerotinia borealis has a psychrophilic lifestyle with a low optimal growth temperature, a narrow host range and geographic distribution. To spread successfully, S. borealis must synthesize proteins adapted to function in its specific environment. The search for signatures of adaptation to S. borealis lifestyle may therefore help revealing proteins critical for colonization of the environment by Sclerotiniaceae fungi. Here, we analyzed amino acids usage and intrinsic protein disorder in alignments of groups of orthologous proteins from the three Sclerotiniaceae species. We found that enrichment in Thr, depletion in Glu and Lys, and low disorder frequency in hot loops are significantly associated with S. borealis proteins. We designed an index to report bias in these properties and found that high index proteins were enriched among secreted proteins in the three Sclerotiniaceae fungi. High index proteins were also enriched in function associated with plant colonization in S. borealis, and in in planta-induced genes in S. sclerotiorum. We highlight a novel putative antifreeze protein and a novel putative lytic polysaccharide monooxygenase identified through our pipeline as candidate proteins involved in colonization of the environment. Our findings suggest that similar protein signatures associate with S. borealis lifestyle and with secretion in the Sclerotiniaceae. These signatures may be useful for identifying proteins of interest as targets for the management of plant diseases.

Engineering Methylobacterium extorquens for de novo synthesis of the sesquiterpenoid α-humulene from methanol.

Over the last 10 to 15 years, metabolic engineering of microbes has become a versatile tool for high-level de novo synthesis of terpenoids, with the sesquiterpenoids armopha-1,4-diene, farnesene and artemisinic acid as prime examples. However, almost all cell factory approaches towards terpenoids to date have been based on sugar as the raw material, which is mainly used as a food resource and subject to high price volatilities. In this study we present de novo synthesis of the sesquiterpenoid α-humulene from the abundantly available non-food carbon source methanol by metabolically engineered Methylobacterium extorquens AM1. Expression of α-humulene synthase from Zingiber zerumbet in combination with farnesyl pyrophosphate (FPP) synthase from Saccharomyces cerevisiae led to concentrations of up to 18 mg/L α-humulene. Introduction of a prokaryotic mevalonate pathway from Myxococcus xanthus in combination with ribosome binding site optimization of α-humulene and FPP synthases increased product concentration 3-fold. This value was additionally raised by 30% using a carotenoid synthesis deficient mutant strain. Final product concentrations of up to 1.65 g/L were obtained in methanol limited fed-batch cultivations, which is the highest titer of de novo synthesized α-humulene reported to date. This study demonstrates the potential of M. extorquens as a future platform strain for the production of high-value terpenoids from the alternative carbon source methanol.

Co-consumption of methanol and succinate by Methylobacterium extorquens AM1.

Methylobacterium extorquens AM1 is a facultative methylotrophic Alphaproteobacterium and has been subject to intense study under pure methylotrophic as well as pure heterotrophic growth conditions in the past. Here, we investigated the metabolism of M. extorquens AM1 under mixed substrate conditions, i.e., in the presence of methanol plus succinate. We found that both substrates were co-consumed, and the carbon conversion was two-thirds from succinate and one-third from methanol relative to mol carbon. (13)C-methanol labeling and liquid chromatography mass spectrometry analyses revealed the different fates of the carbon from the two substrates. Methanol was primarily oxidized to CO(2) for energy generation. However, a portion of the methanol entered biosynthetic reactions via reactions specific to the one-carbon carrier tetrahydrofolate. In contrast, succinate was primarily used to provide precursor metabolites for bulk biomass production. This work opens new perspectives on the role of methylotrophy when substrates are simultaneously available, a situation prevailing under environmental conditions.

Complete genome sequences of six strains of the genus Methylobacterium.

The complete and assembled genome sequences were determined for six strains of the alphaproteobacterial genus Methylobacterium, chosen for their key adaptations to different plant-associated niches and environmental constraints.

The ethylmalonyl-CoA pathway is used in place of the glyoxylate cycle by Methylobacterium extorquens AM1 during growth on acetate.

Acetyl-CoA assimilation was extensively studied in organisms harboring the glyoxylate cycle. In this study, we analyzed the metabolism of the facultative methylotroph Methylobacterium extorquens AM1, which lacks isocitrate lyase, the key enzyme in the glyoxylate cycle, during growth on acetate. MS/MS-based proteomic analysis revealed that the protein repertoire of M. extorquens AM1 grown on acetate is similar to that of cells grown on methanol and includes enzymes of the ethylmalonyl-CoA (EMC) pathway that were recently shown to operate during growth on methanol. Dynamic 13C labeling experiments indicate the presence of distinct entry points for acetate: the EMC pathway and the TCA cycle. 13C steady-state metabolic flux analysis showed that oxidation of acetyl-CoA occurs predominantly via the TCA cycle and that assimilation occurs via the EMC pathway. Furthermore, acetyl-CoA condenses with the EMC pathway product glyoxylate, resulting in malate formation. The latter, also formed by the TCA cycle, is converted to phosphoglycerate by a reaction sequence that is reversed with respect to the serine cycle. Thus, the results obtained in this study reveal the utilization of common pathways during the growth of M. extorquens AM1 on C1 and C2 compounds, but with a major redirection of flux within the central metabolism. Furthermore, our results indicate that the metabolic flux distribution is highly complex in this model methylotroph during growth on acetate and is fundamentally different from organisms using the glyoxylate cycle.