Functional analysis of the apple fruit microbiome based on shotgun metagenomic sequencing of conventional and organic orchard samples.

Fruits harbour abundant and diverse microbial communities that protect them from post-harvest pathogens. Identification of functional traits associated with a given microbiota can provide a better understanding of their potential influence. Here, we focused on the epiphytic microbiome of apple fruit. We suggest that shotgun metagenomic data can indicate specific functions carried out by different groups and provide information on their potential impact. Samples were collected from the surface of 'Golden Delicious' apples from four orchards that differ in their geographic location and management practice. Approximately 1 million metagenes were predicted based on a high-quality assembly. Functional profiling of the microbiome of fruits from orchards differing in their management practice revealed a functional shift in the microbiota. The organic orchard microbiome was enriched in pathways involved in plant defence activities; the conventional orchard microbiome was enriched in pathways related to the synthesis of antibiotics. The functional significance of the variations was explored using microbial network modelling algorithms to reveal the metabolic role of specific phylogenetic groups. The analysis identified several associations supported by other published studies. For example, the analysis revealed the nutritional dependencies of the Capnodiales group, including the Alternaria pathogen, on aromatic compounds.

Comparative study of bacterial community dynamics in different soils following application of the herbicide atrazine.

Microbial communities in cultivated soils control the fate of pollutants associated with agricultural practice. The present study was designed to explore the response of bacterial communities to the application of the widely-used herbicide atrazine in three different crop fields that differ significantly in their physicochemical structure and nutritional content: the nutrient-rich (with relatively high carbon and nitrogen content) Newe Yaar (NY) and Ha-Ogen (HO) soils and the nutrient-poor, sandy Sde-Eliyahu (SE) soil. The 16 S rRNA gene amplicon sequencing revealed the nutrient poor HO soil differs in its response to atrazine in comparison to the two nutrient-rich soils both in the shortest persistence of atrazine and its effect on community structure and composition. Potential reported bacterial degraders of atrazine such as Pseudomonas, Clostridium and Bacillus were more abundant in contaminated sandy/poor soils (HO) whereas bacteria known for nitrogen cycling such as Azospirillum, Sinorhizobium, Nitrospira and Azohydromonas were significantly more abundant in the nutrient rich contaminated SE soils. No significant increase of potential indigenous degrader Arthrobacter was detected in SE and NY soils whereas a significant increase was recorded with HO soils. An overall shift in bacterial community composition following atrazine application was observed only in the nutrient poor soil. Understanding atrazine persistence and microbiome response to its application of in dependence with soil types serve the design of precision application strategies.

Electron acceptors determine the BTEX degradation capacity of anaerobic microbiota via regulating the microbial community.

Anaerobic degradation is the major pathway for microbial degradation of benzene, toluene, ethylbenzene, and xylenes (BTEX) under electron acceptor lacking conditions. However, how exogenous electron acceptors modulate BTEX degradation through shaping the microbial community structure remains poorly understood. Here, we investigated the effect of various exogenous electron acceptors on BTEX degradation as well as methane production in anaerobic microbiota, which were enriched from the same contaminated soil. It was found that the BTEX degradation capacities of the anaerobic microbiota gradually increased along with the increasing redox potentials of the exogenous electron acceptors supplemented (WE: Without exogenous electron acceptors < SS: Sulfate supplement < FS: Ferric iron supplement < NS: Nitrate supplement), while the complexity of the co-occurring networks (e.g., avgK and links) of the microbiota gradually decreased, showing that microbiota supplemented with higher redox potential electron acceptors were less dependent on the formation of complex microbial interactions to perform BTEX degradation. Microbiota NS showed the highest degrading capacity and the broadest substrate-spectrum for BTEX, and it could metabolize BTEX through multiple modules which not only contained fewer species but also different key microbial taxa (eg. Petrimonas, Achromobacter and Comamonas). Microbiota WE and FS, with the highest methanogenic capacities, shared common core species such as Sedimentibacter, Acetobacterium, Methanobacterium and Smithella/Syntrophus, which cooperated with Geobacter (microbiota WE) or Desulfoprunum (microbiota FS) to perform BTEX degradation and methane production. This study demonstrates that electron acceptors may alter microbial function by reshaping microbial community structure and regulating microbial interactions and provides guidelines for electron acceptor selection for bioremediation of aromatic pollutant-contaminated anaerobic sites.

Global analysis of the apple fruit microbiome: are all apples the same?

We present the first worldwide study on the apple (Malus × domestica) fruit microbiome that examines questions regarding the composition and the assembly of microbial communities on and in apple fruit. Results revealed that the composition and structure of the fungal and bacterial communities associated with apple fruit vary and are highly dependent on geographical location. The study also confirmed that the spatial variation in the fungal and bacterial composition of different fruit tissues exists at a global level. Fungal diversity varied significantly in fruit harvested in different geographical locations and suggests a potential link between location and the type and rate of postharvest diseases that develop in each country. The global core microbiome of apple fruit was represented by several beneficial microbial taxa and accounted for a large fraction of the fruit microbial community. The study provides foundational information about the apple fruit microbiome that can be utilized for the development of novel approaches for the management of fruit quality and safety, as well as for reducing losses due to the establishment and proliferation of postharvest pathogens. It also lays the groundwork for studying the complex microbial interactions that occur on apple fruit surfaces.

A Synergistic Consortium Involved in rac-Dichlorprop Degradation as Revealed by DNA Stable Isotope Probing and Metagenomic Analysis.

rac-Dichlorprop, a commonly used phenoxyalkanoic acid herbicide, is frequently detected in environments and poses threats to environmental safety and human health. Microbial consortia are thought to play key roles in rac-dichlorprop degradation. However, the compositions of the microbial consortia involved in rac-dichlorprop degradation remain largely unknown. In this study, DNA stable isotope probing (SIP) and metagenomic analysis were integrated to reveal the key microbial consortium responsible for rac-dichlorprop degradation in a rac-dichlorprop-degrading enrichment. OTU340 (Sphingobium sp.) and OTU348 (Sphingopyxis sp.) were significantly enriched in the rac-[13C]dichlorprop-labeled heavy DNA fractions. A rac-dichlorprop degrader, Sphingobium sp. strain L3, was isolated from the enrichment by a traditional enrichment method but with additional supplementation of the antibiotic ciprofloxacin, which was instructed by metagenomic analysis of the associations between rac-dichlorprop degraders and antibiotic resistance genes. As revealed by functional profiling of the metagenomes of the heavy DNA, the genes rdpA and sdpA, involved in the initial degradation of the (R)- and (S)-enantiomers of dichlorprop, respectively, were mostly taxonomically assigned to Sphingobium species, indicating that Sphingopyxis species might harbor novel dichlorprop-degrading genes. In addition, taxonomically diverse bacterial genera such as Dyella, Sphingomonas, Pseudomonas, and Achromobacter were presumed to synergistically cooperate with the key degraders Sphingobium/Sphingopyxis for enhanced degradation of rac-dichlorprop. IMPORTANCE Understanding of the key microbial consortium involved in the degradation of the phenoxyalkanoic acid herbicide rac-dichlorprop is pivotal for design of synergistic consortia used for enhanced bioremediation of herbicide-contaminated sites. However, the composition of the microbial consortium and the interactions between community members during the biodegradation of rac-dichlorprop are unclear. In this study, DNA-SIP and metagenomic analysis were integrated to reveal that the metabolite 2,4-dichlorophenol degraders Dyella, Sphingomonas, Pseudomonas, and Achromobacter synergistically cooperated with the key degraders Sphingobium/Sphingopyxis for enhanced degradation of rac-dichlorprop. Our study provides new insights into the synergistic degradation of rac-dichlorprop at the community level and implies the existence of novel degrading genes for rac-dichlorprop in nature.

High-density NGS-based map construction and genetic dissection of fruit shape and rind netting in Cucumis melo.

Melon is an important crop that exhibits broad variation for fruit morphology traits that are the substrate for genetic mapping efforts. In the post-genomic era, the link between genetic maps and physical genome assemblies is key for leveraging QTL mapping results for gene cloning and breeding purposes. Here, using a population of 164 melon recombinant inbred lines (RILs) that were subjected to genotyping-by-sequencing, we constructed and compared high-density sequence- and linkage-based recombination maps that were aligned to the reference melon genome. These analyses reveal the genome-wide variation in recombination frequency and highlight regions of disrupted collinearity between our population and the reference genome. The population was phenotyped over 3 years for fruit size and shape as well as rind netting. Four QTLs were detected for fruit size, and they act in an additive manner, while significant epistatic interaction was found between two neutral loci for this trait. Fruit shape displayed transgressive segregation that was explained by the action of four QTLs, contributed by alleles from both parents. The complexity of rind netting was demonstrated on a collection of 177 diverse accessions. Further dissection of netting in our RILs population, which is derived from a cross of smooth and densely netted parents, confirmed the intricacy of this trait and the involvement of major locus and several other interacting QTLs. A major netting QTL on chromosome 2 co-localized with results from two additional populations, paving the way for future study toward identification of a causative gene for this trait.

The biosynthetic pathway of the nonsugar, high-intensity sweetener mogroside V from Siraitia grosvenorii.

The consumption of sweeteners, natural as well as synthetic sugars, is implicated in an array of modern-day health problems. Therefore, natural nonsugar sweeteners are of increasing interest. We identify here the biosynthetic pathway of the sweet triterpenoid glycoside mogroside V, which has a sweetening strength of 250 times that of sucrose and is derived from mature fruit of luo-han-guo (Siraitia grosvenorii, monk fruit). A whole-genome sequencing of Siraitia, leading to a preliminary draft of the genome, was combined with an extensive transcriptomic analysis of developing fruit. A functional expression survey of nearly 200 candidate genes identified the members of the five enzyme families responsible for the synthesis of mogroside V: squalene epoxidases, triterpenoid synthases, epoxide hydrolases, cytochrome P450s, and UDP-glucosyltransferases. Protein modeling and docking studies corroborated the experimentally proven functional enzyme activities and indicated the order of the metabolic steps in the pathway. A comparison of the genomic organization and expression patterns of these Siraitia genes with the orthologs of other Cucurbitaceae implicates a strikingly coordinated expression of the pathway in the evolution of this species-specific and valuable metabolic pathway. The genomic organization of the pathway genes, syntenously preserved among the Cucurbitaceae, indicates, on the other hand, that gene clustering cannot account for this novel secondary metabolic pathway.

Modeling trophic dependencies and exchanges among insects' bacterial symbionts in a host-simulated environment.

BACKGROUND: Individual organisms are linked to their communities and ecosystems via metabolic activities. Metabolic exchanges and co-dependencies have long been suggested to have a pivotal role in determining community structure. In phloem-feeding insects such metabolic interactions with bacteria enable complementation of their deprived nutrition. The phloem-feeding whitefly Bemisia tabaci (Hemiptera: Aleyrodidae) harbors an obligatory symbiotic bacterium, as well as varying combinations of facultative symbionts. This well-defined bacterial community in B. tabaci serves here as a case study for a comprehensive and systematic survey of metabolic interactions within the bacterial community and their associations with documented occurrences of bacterial combinations. We first reconstructed the metabolic networks of five common B. tabaci symbionts genera (Portiera, Rickettsia, Hamiltonella, Cardinium and Wolbachia), and then used network analysis approaches to predict: (1) species-specific metabolic capacities in a simulated bacteriocyte-like environment; (2) metabolic capacities of the corresponding species' combinations, and (3) dependencies of each species on different media components. RESULTS: The predictions for metabolic capacities of the symbionts in the host environment were in general agreement with previously reported genome analyses, each focused on the single-species level. The analysis suggests several previously un-reported routes for complementary interactions and estimated the dependency of each symbiont in specific host metabolites. No clear association was detected between metabolic co-dependencies and co-occurrence patterns. CONCLUSIONS: The analysis generated predictions for testable hypotheses of metabolic exchanges and co-dependencies in bacterial communities and by crossing them with co-occurrence profiles, contextualized interaction patterns into a wider ecological perspective.

Dyella-Like Bacterium Isolated from an Insect as a Potential Biocontrol Agent Against Grapevine Yellows.

Yellows diseases, caused by phytopathogenic bacteria of the genus Phytoplasma, are a major threat to grapevines worldwide. Because conventional applications against this pathogen are inefficient and disease management is highly challenging, the use of beneficial bacteria has been suggested as a biocontrol solution. A Dyella-like bacterium (DLB), isolated from the Israeli insect vector of grapevine yellows (Hyalesthes obsoletus), was suggested to be an endophyte. To test this hypothesis, the bacterium was introduced by spraying the plant leaves, and it had no apparent phytotoxicity to grapevine. Fluorescent in situ hybridization analysis showed that DLB is colonizing grapevine phloem. Because phytoplasmas inhabit the same niche, DLB interactions with this phytopathogen were examined. When the isolate was introduced to phytoplasma-infected Chardonnay plantlets, morphological disease symptoms were markedly reduced. The mode of DLB action was then tested using bioinformatics and system biology tools. DLB genome analysis suggested that the ability to reduce phytoplasma symptoms is related to inhibition of the pathogenic bacterium. These results provide the first step in examining the potential of DLB as a biological control agent against phytoplasmas in grapevine and, possibly, other agricultural crops.

Prospects for Biological Soilborne Disease Control: Application of Indigenous Versus Synthetic Microbiomes.

Biological disease control of soilborne plant diseases has traditionally employed the biopesticide approach whereby single strains or strain mixtures are introduced into production systems through inundative/inoculative release. The approach has significant barriers that have long been recognized, including a generally limited spectrum of target pathogens for any given biocontrol agent and inadequate colonization of the host rhizosphere, which can plague progress in the utilization of this resource in commercial field-based crop production systems. Thus, although potential exists, this model has continued to lag in its application. New omics' tools have enabled more rapid screening of microbial populations allowing for the identification of strains with multiple functional attributes that may contribute to pathogen suppression. Similarly, these technologies also enable the characterization of consortia in natural systems which provide the framework for construction of synthetic microbiomes for disease control. Harnessing the potential of the microbiome indigenous to agricultural soils for disease suppression through application of specific management strategies has long been a goal of plant pathologists. Although this tactic also possesses limitation, our enhanced understanding of functional attributes of suppressive soil systems through application of community and metagenomic analysis methods provide opportunity to devise effective resource management schemes. As these microbial communities in large part are fostered by the resources endemic to soil and the rhizosphere, substrate mediated recruitment of disease-suppressive microbiomes constitutes a practical means to foster their establishment in crop production systems.

Multidimensional patterns of metabolic response in abiotic stress-induced growth of Arabidopsis thaliana.

KEY MESSAGE: Contextualization of specific transcriptional responses of Arabidopsis within the stress-tissue-time perspective provides a simplified representation of the cellular transcriptional response pathways to abiotic stress, while reducing the dimensions in gene-oriented response description. Crops resistant to abiotic stresses are a long-term goal of many research programs, thus understanding the progression of stress responses is of great interest. We reanalyzed the AtGenExpress transcription dataset to go beyond gene-level characterization, and to contextualize the discrete information into (1) a process-level signature of stress-specific, time-specific, and tissue-specific responses and (2) identify patterns of response progression across a time axis. To gain a functional perspective, ∼1000 pathways associated with the differentially-expressed genes were characterized across all experiments. We find that the global response of pathways to stress is multi-dimensional and does not obviously cluster according to stress, time or tissue. The early response to abiotic stress typically involves induction of genes involved in transcription, hormone synthesis and signaling modules; a later response typically involves metabolism of amino acids and secondary metabolites. By linking specific primary and secondary response pathways, we outline possible stress-associated routes of response progression. The contextualization of specific processes within stress-tissue-time perspective provides a simplified representation of cellular response while reducing the dimensions in gene-oriented response description. Such simplified representation allows finding stress-specific markers based on process-combinations pointing whether a stress-specific response was invoked as well as provide a reference point for the conductance of comparative inter-plant study of stress response, bypassing the need in detailed orthologous mapping.

Haoranjiania flava gen. nov., sp. nov., a new member of the family Chitinophagaceae, isolated from activated sludge.

A novel Gram-stain-negative, non-spore-forming and short rod-shaped bacterium, strain LIP-5T, isolated from an activated sludge in a pesticide factory in Xinyi, China, was investigated for its taxonomic position by a polyphasic approach. Cell growth occurred at 16-42 °C (optimum, 30 °C), in the presence of 0-3.0 % (w/v) NaCl (optimum, without NaCl) and at pH 6.0-9.0 (optimum, pH 7.0). Phylogenetic analysis based on 16S rRNA gene sequence comparisons showed that the isolate was a member of the family Chitinophagaceae, with the closest relatives being Arachidicoccus rhizosphaerae Vu-144T (93.5 % similarity), followed by Heliimonas saccharivorans L2-4T (90.5 %) and Chitinophaga ginsengisoli Gsoil 052T (89.6 %). Menaquinone 7 (MK-7) was the predominant respiratory ubiquinone and phosphatidylethanolamine and unidentified lipids were the major polar lipids. The major cellular fatty acids of strain LIP-5T were anteiso-C15 : 0 and iso-C15 : 0, and there were also moderate amounts of iso-C17 : 0 3-OH and C17 : 0 2-OH. The DNA G+C content was 42.3 mol%. On the basis of phenotypic, chemotaxonomic and phylogenetic data, strain LIP-5T was classified as representing a novel species of a new genus, for which the name Haoranjiania flava gen. nov., sp. nov. is proposed. The type strain of Haoranjiania flava is LIP-5T (=CCTCC AB 2015365T=KCTC 42956T).

NetCooperate: a network-based tool for inferring host-microbe and microbe-microbe cooperation.

BACKGROUND: Host-microbe and microbe-microbe interactions are often governed by the complex exchange of metabolites. Such interactions play a key role in determining the way pathogenic and commensal species impact their host and in the assembly of complex microbial communities. Recently, several studies have demonstrated how such interactions are reflected in the organization of the metabolic networks of the interacting species, and introduced various graph theory-based methods to predict host-microbe and microbe-microbe interactions directly from network topology. Using these methods, such studies have revealed evolutionary and ecological processes that shape species interactions and community assembly, highlighting the potential of this reverse-ecology research paradigm. RESULTS: NetCooperate is a web-based tool and a software package for determining host-microbe and microbe-microbe cooperative potential. It specifically calculates two previously developed and validated metrics for species interaction: the Biosynthetic Support Score which quantifies the ability of a host species to supply the nutritional requirements of a parasitic or a commensal species, and the Metabolic Complementarity Index which quantifies the complementarity of a pair of microbial organisms' niches. NetCooperate takes as input a pair of metabolic networks, and returns the pairwise metrics as well as a list of potential syntrophic metabolic compounds. CONCLUSIONS: The Biosynthetic Support Score and Metabolic Complementarity Index provide insight into host-microbe and microbe-microbe metabolic interactions. NetCooperate determines these interaction indices from metabolic network topology, and can be used for small- or large-scale analyses. NetCooperate is provided as both a web-based tool and an open-source Python module; both are freely available online at http://elbo.gs.washington.edu/software_netcooperate.html.

Systems approach for exploring the intricate associations between sweetness, color and aroma in melon fruits.

BACKGROUND: Melon (Cucumis melo) fruits exhibit phenotypic diversity in several key quality determinants such as taste, color and aroma. Sucrose, carotenoids and volatiles are recognized as the key compounds shaping the above corresponding traits yet the full network of biochemical events underlying their synthesis have not been comprehensively described. To delineate the cellular processes shaping fruit quality phenotypes, a population of recombinant inbred lines (RIL) was used as a source of phenotypic and genotypic variations. In parallel, ripe fruits were analyzed for both the quantified level of 77 metabolic traits directly associated with fruit quality and for RNA-seq based expression profiles generated for 27,000 unigenes. First, we explored inter-metabolite association patterns; then, we described metabolites versus gene association patterns; finally, we used the correlation-based associations for predicting uncharacterized synthesis pathways. RESULTS: Based on metabolite versus metabolite and metabolite versus gene association patterns, we divided metabolites into two key groups: a group including ethylene and aroma determining volatiles whose accumulation patterns are correlated with the expression of genes involved in the glycolysis and TCA cycle pathways; and a group including sucrose and color determining carotenoids whose accumulation levels are correlated with the expression of genes associated with plastid formation. CONCLUSIONS: The study integrates multiple processes into a genome scale perspective of cellular activity. This lays a foundation for deciphering the role of gene markers associated with the determination of fruit quality traits.

Comparison of Four Comamonas Catabolic Plasmids Reveals the Evolution of pBHB To Catabolize Haloaromatics.

Comamonas plasmids play important roles in shaping the phenotypes of their hosts and the adaptation of these hosts to changing environments, and understanding the evolutionary strategy of these plasmids is thus of great concern. In this study, the sequence of the 119-kb 3,5-dibromo-4-hydroxybenzonitrile-catabolizing plasmid pBHB from Comamonas sp. strain 7D-2 was studied and compared with those of three other Comamonas haloaromatic catabolic plasmids. Incompatibility group determination based on a phylogenetic analysis of 24 backbone gene proteins, as well as TrfA, revealed that these four plasmids all belong to the IncP-1ß subgroup. Comparison of the four plasmids revealed a conserved backbone region and diverse genetic-load regions. The four plasmids share a core genome consisting of 40 genes (>50% similarities) and contain 12 to 50 unique genes each, most of which are xenobiotic-catabolic genes. Two functional reductive dehalogenase gene clusters are specifically located on pBHB, showing distinctive evolution of pBHB for haloaromatics. The higher catabolic ability of the bhbA2B2 cluster than the bhbAB cluster may be due to the transcription levels and the character of the dehalogenase gene itself rather than that of its extracytoplasmic binding receptor gene. The plasmid pBHB is riddled with transposons and insertion sequence (IS) elements, and ISs play important roles in the evolution of pBHB. The analysis of the origin of the bhb genes on pBHB suggested that these accessory genes evolved independently. Our work provides insights into the evolutionary strategies of Comamonas plasmids, especially into the adaptation mechanism employed by pBHB for haloaromatics.

Electro-stimulation modulates syntrophic interactions in methanogenic toluene-degrading microbiota for enhanced functionality.

Syntrophy achieved via microbial cooperation is vital for anaerobic hydrocarbon degradation and methanogenesis. However, limited understanding of the metabolic division of labor and electronic interactions in electro-stimulated microbiota has impeded the development of enhanced biotechnologies for degrading hydrocarbons to methane. Here, compared to the non-electro-stimulated methanogenic toluene-degrading microbiota, electro-stimulation at 800 mV promoted toluene degradation and methane production efficiencies by 11.49 %-14.76 % and 75.58 %-290.11 %, respectively. Hydrocarbon-degrading gene bamA amplification and metagenomic sequencing analyses revealed that f_Syntrophobacteraceae MAG116 may act as a toluene degrader in the non-electro-stimulated microbiota, which was proposed to establish electron syntrophy with the acetoclastic methanogen Methanosarcina spp. (or Methanothrix sp.) through e-pili or shared acetate. In the electro-stimulated microbiota, 37.22 ± 4.33 % of Desulfoprunum sp. (affiliated f_Desulfurivibrionaceae MAG10) and 58.82 ± 3.74 % of the hydrogenotrophic methanogen Methanobacterium sp. MAG74 were specifically recruited to the anode and cathode, respectively. The potential electrogen f_Desulfurivibrionaceae MAG10 engaged in interspecies electron transfer with both syntroph f_Syntrophobacteraceae MAG116 and the anode, which might be facilitated by c-type cytochromes (e.g., ImcH, OmcT, and PilZ). Moreover, upon capturing electrons from the external circuit, the hydrogen-producing electrotroph Aminidesulfovibrio sp. MAG60 could share electrons and hydrogen with the methanogen Methanobacterium sp. MAG74, which uniquely harbored hydrogenase genes ehaA-R and ehbA-P. This study elucidates the microbial interaction mechanisms underlying the enhanced metabolic efficiency of the electro-stimulated methanogenic toluene-degrading microbiota, and emphasizes the significance of metabolic and electron syntrophic interactions in maintaining the stability of microbial community functionality.

Engineering natural microbiomes toward enhanced bioremediation by microbiome modeling.

Engineering natural microbiomes for biotechnological applications remains challenging, as metabolic interactions within microbiomes are largely unknown, and practical principles and tools for microbiome engineering are still lacking. Here, we present a combinatory top-down and bottom-up framework to engineer natural microbiomes for the construction of function-enhanced synthetic microbiomes. We show that application of herbicide and herbicide-degrader inoculation drives a convergent succession of different natural microbiomes toward functional microbiomes (e.g., enhanced bioremediation of herbicide-contaminated soils). We develop a metabolic modeling pipeline, SuperCC, that can be used to document metabolic interactions within microbiomes and to simulate the performances of different microbiomes. Using SuperCC, we construct bioremediation-enhanced synthetic microbiomes based on 18 keystone species identified from natural microbiomes. Our results highlight the importance of metabolic interactions in shaping microbiome functions and provide practical guidance for engineering natural microbiomes.

NetCmpt: a network-based tool for calculating the metabolic competition between bacterial species.

UNLABELLED: NetCmpt is a tool for calculating the competitive potential between pairs of bacterial species. The score describes the effective metabolic overlap (EMO) between two species, derived from analyzing the topology of the corresponding metabolic models. NetCmpt is based on the EMO algorithm, developed and validated in previous studies. It takes as input lists of species-specific enzymatic reactions (EC numbers) and generates a matrix of the potential competition scores between all pairwise combinations. AVAILABILITY AND IMPLEMENTATION: NetCmpt is provided as both a web tool and a software package, designed for the use of non-computational biologists. The NetCmpt web tool, software, examples, and documentation are freely available online at http://app.agri.gov.il/shiri/NetComp.php.

Association between translation efficiency and horizontal gene transfer within microbial communities.

Horizontal gene transfer (HGT) is a major force in microbial evolution. Previous studies have suggested that a variety of factors, including restricted recombination and toxicity of foreign gene products, may act as barriers to the successful integration of horizontally transferred genes. This study identifies an additional central barrier to HGT-the lack of co-adaptation between the codon usage of the transferred gene and the tRNA pool of the recipient organism. Analyzing the genomic sequences of more than 190 microorganisms and the HGT events that have occurred between them, we show that the number of genes that were horizontally transferred between organisms is positively correlated with the similarity between their tRNA pools. Those genes that are better adapted to the tRNA pools of the target genomes tend to undergo more frequent HGT. At the community (or environment) level, organisms that share a common ecological niche tend to have similar tRNA pools. These results remain significant after controlling for diverse ecological and evolutionary parameters. Our analysis demonstrates that there are bi-directional associations between the similarity in the tRNA pools of organisms and the number of HGT events occurring between them. Similar tRNA pools between a donor and a host tend to increase the probability that a horizontally acquired gene will become fixed in its new genome. Our results also suggest that frequent HGT may be a homogenizing force that increases the similarity in the tRNA pools of organisms within the same community.

The rise of the sourdough: Genome-scale metabolic modeling-based approach to design sourdough starter communities with tailored-made properties.

Sourdough starters harbor microbial consortia that benefit the final product's aroma and volume. The complex nature of these spontaneously developed communities raises challenges in predicting the fermentation phenotypes. Herein, we demonstrated for the first time in this field the potential of genome-scale metabolic modeling (GEMs) in the study of sourdough microbial communities. Broad in-silico modeling of microbial growth was applied on communities composed of yeast (Saccharomyces cerevisiae) and different Lactic Acid Bacteria (LAB) species, which mainly predominate in sourdough starters. Simulations of model-represented communities associated specific bacterial compositions with sourdough phenotypes. Based on ranking the phenotypic performances of different combinations, Pediococcus spp. - Lb. sakei group members were predicted to have an optimal effect considering the increase in S. cerevisiae growth abilities and overall CO2 secretion rates. Flux Balance Analysis (FBA) revealed mutual relationships between the Pediococcus spp. - Lb. sakei group members and S. cerevisiae through bidirectional nutrient dependencies, and further underlined that these bacteria compete with the yeast over nutrients to a lesser extent than the rest LAB species. Volatile compounds (VOCs) production was further modeled, identifying species-specific and community-related VOCs production profiles. The in-silico models' predictions were validated by experimentally building synthetic sourdough communities and assessing the fermentation phenotypes. The Pediococcus spp. - Lb. sakei group was indeed associated with increased yeast cell counts and fermentation rates, demonstrating a 25 % increase in the average leavening rates during the first 10 fermentation hours compared to communities with a lower representation of these group members. Overall, these results provide a possible novel strategy towards the de-novo design of sourdough starter communities with tailored-made characterizations, including a shortened leavening period.