Unveiling abundance-dependent metabolic phenotypes of microbial communities.

IMPORTANCE: In nature, organisms live in communities and not as isolated species, and their interactions provide a source of resilience to environmental disturbances. Despite their importance in ecology, human health, and industry, understanding how organisms interact in different environments remains an open question. In this work, we provide a novel approach that, only using genomic information, studies the metabolic phenotype exhibited by communities, where the exploration of suboptimal growth flux distributions and the composition of a community allows to unveil its capacity to respond to environmental changes, shedding light of the degrees of metabolic plasticity inherent to the community.

Integration of Biological Networks for Acidithiobacillus thiooxidans Describes a Modular Gene Regulatory Organization of Bioleaching Pathways.

Acidithiobacillus thiooxidans is one of the most studied biomining species, highlighting its ability to oxidize reduced inorganic sulfur compounds, coupled with its elevated capacity to live under an elevated concentration of heavy metals. In this work, using an in silico semi-automatic genome scale approach, two biological networks for A. thiooxidans Licanantay were generated: (i) An affinity transcriptional regulatory network composed of 42 regulatory family genes and 1,501 operons (57% genome coverage) linked through 2,646 putative DNA binding sites (arcs), (ii) A metabolic network reconstruction made of 523 genes and 1,203 reactions (22 pathways related to biomining processes). Through the identification of confident connections between both networks (V-shapes), it was possible to identify a sub-network of transcriptional factor (34 regulators) regulating genes (61 operons) encoding for proteins involved in biomining-related pathways. Network analysis suggested that transcriptional regulation of biomining genes is organized into different modules. The topological parameters showed a high hierarchical organization by levels inside this network (14 layers), highlighting transcription factors CysB, LysR, and IHF as complex modules with high degree and number of controlled pathways. In addition, it was possible to identify transcription factor modules named primary regulators (not controlled by other regulators in the sub-network). Inside this group, CysB was the main module involved in gene regulation of several bioleaching processes. In particular, metabolic processes related to energy metabolism (such as sulfur metabolism) showed a complex integrated regulation, where different primary regulators controlled several genes. In contrast, pathways involved in iron homeostasis and oxidative stress damage are mainly regulated by unique primary regulators, conferring Licanantay an efficient, and specific metal resistance response. This work shows new evidence in terms of transcriptional regulation at a systems level and broadens the study of bioleaching in A. thiooxidans species.

Altered zinc balance in the Atp7b-/- mouse reveals a mechanism of copper toxicity in Wilson disease.

Wilson disease (WD) is an autosomal recessive disorder caused by mutation in the ATP7B gene that affects copper transport in the body. ATP7B mutation damages copper transporter function, ultimately resulting in excessive copper accumulation and subsequent toxicity in both the liver and brain. Mechanisms of copper toxicity, however, are not well defined. The Atp7b-/- mouse model is well-characterized and presents a hepatic phenotype consistent with WD. In this study, we found that the untreated Atp7b-/- mice accumulate approximately 2-fold excess hepatic zinc compared to the wild type. We used targeted transcriptomics and proteomics to analyze the molecular events associated with zinc and copper accumulation in the Atp7b-/- mouse liver. Altered gene expression of Zip5 and ZnT1 zinc transporters indicated a transcriptional homeostatic response, while increased copper/zinc ratios associated with high levels of metallothioneins 1 and 2, indicated altered Zn availability in cells. These data suggest that copper toxicity in Wilson disease includes effects on zinc-dependent proteins. Transcriptional network analysis of RNA-seq data reveals an interconnected network of transcriptional activators with over-representation of zinc-dependent and zinc-responsive transcription factors. In the context of previous research, these observations support the hypothesis that mechanisms of copper toxicity include disruption of intracellular zinc distribution in liver cells. The translational significance of this work lies in oral zinc supplementation in treatment for WD, which is thought to mediate protective effects through the induction of metallothionein synthesis in the intestine. This work indicates broader impacts of altered zinc-copper balance in WD, including global transcriptional responses and altered zinc balance in the liver.

The genome sequence of the soft-rot fungus Penicillium purpurogenum reveals a high gene dosage for lignocellulolytic enzymes.

The high lignocellulolytic activity displayed by the soft-rot fungus Penicillium purpurogenum has made it a target for the study of novel lignocellulolytic enzymes. We have obtained a reference genome of 36.2 Mb of non-redundant sequence (11,057 protein-coding genes). The 49 largest scaffolds cover 90% of the assembly, and Core Eukaryotic Genes Mapping Approach (CEGMA) analysis reveals that our assembly captures almost all protein-coding genes. RNA-seq was performed and 93.1% of the reads aligned to the assembled genome. These data, plus the independent sequencing of a set of genes of lignocellulose-degrading enzymes, validate the quality of the genome sequence. P. purpurogenum shows a higher number of proteins with CAZy motifs, transcription factors and transporters as compared to other sequenced Penicillia. These results demonstrate the great potential for lignocellulolytic activity of this fungus and the possible use of its enzymes in related industrial applications.

Traceability, reproducibility and wiki-exploration for "à-la-carte" reconstructions of genome-scale metabolic models.

Genome-scale metabolic models have become the tool of choice for the global analysis of microorganism metabolism, and their reconstruction has attained high standards of quality and reliability. Improvements in this area have been accompanied by the development of some major platforms and databases, and an explosion of individual bioinformatics methods. Consequently, many recent models result from "à la carte" pipelines, combining the use of platforms, individual tools and biological expertise to enhance the quality of the reconstruction. Although very useful, introducing heterogeneous tools, that hardly interact with each other, causes loss of traceability and reproducibility in the reconstruction process. This represents a real obstacle, especially when considering less studied species whose metabolic reconstruction can greatly benefit from the comparison to good quality models of related organisms. This work proposes an adaptable workspace, AuReMe, for sustainable reconstructions or improvements of genome-scale metabolic models involving personalized pipelines. At each step, relevant information related to the modifications brought to the model by a method is stored. This ensures that the process is reproducible and documented regardless of the combination of tools used. Additionally, the workspace establishes a way to browse metabolic models and their metadata through the automatic generation of ad-hoc local wikis dedicated to monitoring and facilitating the process of reconstruction. AuReMe supports exploration and semantic query based on RDF databases. We illustrate how this workspace allowed handling, in an integrated way, the metabolic reconstructions of non-model organisms such as an extremophile bacterium or eukaryote algae. Among relevant applications, the latter reconstruction led to putative evolutionary insights of a metabolic pathway.

Codon usage bias reveals genomic adaptations to environmental conditions in an acidophilic consortium.

The analysis of codon usage bias has been widely used to characterize different communities of microorganisms. In this context, the aim of this work was to study the codon usage bias in a natural consortium of five acidophilic bacteria used for biomining. The codon usage bias of the consortium was contrasted with genes from an alternative collection of acidophilic reference strains and metagenome samples. Results indicate that acidophilic bacteria preferentially have low codon usage bias, consistent with both their capacity to live in a wide range of habitats and their slow growth rate, a characteristic probably acquired independently from their phylogenetic relationships. In addition, the analysis showed significant differences in the unique sets of genes from the autotrophic species of the consortium in relation to other acidophilic organisms, principally in genes which code for proteins involved in metal and oxidative stress resistance. The lower values of codon usage bias obtained in this unique set of genes suggest higher transcriptional adaptation to living in extreme conditions, which was probably acquired as a measure for resisting the elevated metal conditions present in the mine.

Combining graph and flux-based structures to decipher phenotypic essential metabolites within metabolic networks.

BACKGROUND: The emergence of functions in biological systems is a long-standing issue that can now be addressed at the cell level with the emergence of high throughput technologies for genome sequencing and phenotyping. The reconstruction of complete metabolic networks for various organisms is a key outcome of the analysis of these data, giving access to a global view of cell functioning. The analysis of metabolic networks may be carried out by simply considering the architecture of the reaction network or by taking into account the stoichiometry of reactions. In both approaches, this analysis is generally centered on the outcome of the network and considers all metabolic compounds to be equivalent in this respect. As in the case of genes and reactions, about which the concept of essentiality has been developed, it seems, however, that some metabolites play crucial roles in system responses, due to the cell structure or the internal wiring of the metabolic network. RESULTS: We propose a classification of metabolic compounds according to their capacity to influence the activation of targeted functions (generally the growth phenotype) in a cell. We generalize the concept of essentiality to metabolites and introduce the concept of the phenotypic essential metabolite (PEM) which influences the growth phenotype according to sustainability, producibility or optimal-efficiency criteria. We have developed and made available a tool, Conquests, which implements a method combining graph-based and flux-based analysis, two approaches that are usually considered separately. The identification of PEMs is made effective by using a logical programming approach. CONCLUSION: The exhaustive study of phenotypic essential metabolites in six genome-scale metabolic models suggests that the combination and the comparison of graph, stoichiometry and optimal flux-based criteria allows some features of the metabolic network functionality to be deciphered by focusing on a small number of compounds. By considering the best combination of both graph-based and flux-based techniques, the Conquests python package advocates for a broader use of these compounds both to facilitate network curation and to promote a precise understanding of metabolic phenotype.

The bioleaching potential of a bacterial consortium.

This work presents the molecular foundation of a consortium of five efficient bacteria strains isolated from copper mines currently used in state of the art industrial-scale biotechnology. The strains Acidithiobacillus thiooxidans Licanantay, Acidiphilium multivorum Yenapatur, Leptospirillum ferriphilum Pañiwe, Acidithiobacillus ferrooxidans Wenelen and Sulfobacillus thermosulfidooxidans Cutipay were selected for genome sequencing based on metal tolerance, oxidation activity and bioleaching of copper efficiency. An integrated model of metabolic pathways representing the bioleaching capability of this consortium was generated. Results revealed that greater efficiency in copper recovery may be explained by the higher functional potential of L. ferriphilum Pañiwe and At. thiooxidans Licanantay to oxidize iron and reduced inorganic sulfur compounds. The consortium had a greater capacity to resist copper, arsenic and chloride ion compared to previously described biomining strains. Specialization and particular components in these bacteria provided the consortium a greater ability to bioleach copper sulfide ores.

Global transcriptional responses of Acidithiobacillus ferrooxidans Wenelen under different sulfide minerals.

In order to provide new information about the adaptation of Acidithiobacillus ferrooxidans during the bioleaching process, the current analysis presents the first report of the global transcriptional response of the native copper mine strain Wenelen (DSM 16786) oxidized under different sulfide minerals. Microarrays were used to measure the response of At. ferrooxidans Wenelen to shifts from iron supplemented liquid cultures (reference state) to the addition of solid substrates enriched in pyrite or chalcopyrite. Genes encoding for energy metabolism showed a similar transcriptional profile for the two sulfide minerals. Interestingly, four operons related to sulfur metabolism were over-expressed during growth on a reduced sulfur source. Genes associated with metal tolerance (RND and ATPases type P) were up-regulated in the presence of pyrite or chalcopyrite. These results suggest that At. ferrooxidans Wenelen presents an efficient transcriptional system developed to respond to environmental conditions, namely the ability to withstand high copper concentrations.

Putative bacterial interactions from metagenomic knowledge with an integrative systems ecology approach.

Following the trend of studies that investigate microbial ecosystems using different metagenomic techniques, we propose a new integrative systems ecology approach that aims to decipher functional roles within a consortium through the integration of genomic and metabolic knowledge at genome scale. For the sake of application, using public genomes of five bacterial strains involved in copper bioleaching: Acidiphilium cryptum, Acidithiobacillus ferrooxidans, Acidithiobacillus thiooxidans, Leptospirillum ferriphilum, and Sulfobacillus thermosulfidooxidans, we first reconstructed a global metabolic network. Next, using a parsimony assumption, we deciphered sets of genes, called Sets from Genome Segments (SGS), that (1) are close on their respective genomes, (2) take an active part in metabolic pathways and (3) whose associated metabolic reactions are also closely connected within metabolic networks. Overall, this SGS paradigm depicts genomic functional units that emphasize respective roles of bacterial strains to catalyze metabolic pathways and environmental processes. Our analysis suggested that only few functional metabolic genes are horizontally transferred within the consortium and that no single bacterial strain can accomplish by itself the whole copper bioleaching. The use of SGS pinpoints a functional compartmentalization among the investigated species and exhibits putative bacterial interactions necessary for promoting these pathways.

Nutrigenomics analysis reveals that copper deficiency and dietary sucrose up-regulate inflammation, fibrosis and lipogenic pathways in a mature rat model of nonalcoholic fatty liver disease.

Nonalcoholic fatty liver disease (NAFLD) prevalence is increasing worldwide, with the affected US population estimated near 30%. Diet is a recognized risk factor in the NAFLD spectrum, which includes nonalcoholic steatohepatitis (NASH) and fibrosis. Low hepatic copper (Cu) was recently linked to clinical NAFLD/NASH severity. Simple sugar consumption including sucrose and fructose is implicated in NAFLD, while consumption of these macronutrients also decreases liver Cu levels. Though dietary sugar and low Cu are implicated in NAFLD, transcript-level responses that connect diet and pathology are not established. We have developed a mature rat model of NAFLD induced by dietary Cu deficiency, human-relevant high sucrose intake (30% w/w) or both factors in combination. Compared to the control diet with adequate Cu and 10% (w/w) sucrose, rats fed either high-sucrose or low-Cu diet had increased hepatic expression of genes involved in inflammation and fibrogenesis, including hepatic stellate cell activation, while the combination of diet factors also increased ATP citrate lyase and fatty acid synthase gene transcription (fold change > 2, P < 0.02). Low dietary Cu decreased hepatic and serum Cu (P ≤ 0.05), promoted lipid peroxidation and induced NAFLD-like histopathology, while the combined factors also induced fasting hepatic insulin resistance and liver damage. Neither low Cu nor 30% sucrose in the diet led to enhanced weight gain. Taken together, transcript profiles, histological and biochemical data indicate that low Cu and high sucrose promote hepatic gene expression and physiological responses associated with NAFLD and NASH, even in the absence of obesity or severe steatosis.

A new genome of Acidithiobacillus thiooxidans provides insights into adaptation to a bioleaching environment.

Acidithiobacillus thiooxidans is a sulfur oxidizing acidophilic bacterium found in many sulfur-rich environments. It is particularly interesting due to its role in bioleaching of sulphide minerals. In this work, we report the genome sequence of At. thiooxidans Licanantay, the first strain from a copper mine to be sequenced and currently used in bioleaching industrial processes. Through comparative genomic analysis with two other At. thiooxidans non-metal mining strains (ATCC 19377 and A01) we determined that these strains share a large core genome of 2109 coding sequences and a high average nucleotide identity over 98%. Nevertheless, the presence of 841 strain-specific genes (absent in other At. thiooxidans strains) suggests a particular adaptation of Licanantay to its specific biomining environment. Among this group, we highlight genes encoding for proteins involved in heavy metal tolerance, mineral cell attachment and cysteine biosynthesis. Several of these genes were located near genetic motility genes (e.g. transposases and integrases) in genomic regions of over 10 kbp absent in the other strains, suggesting the presence of genomic islands in the Licanantay genome probably produced by horizontal gene transfer in mining environments.

Sulfobacillus thermosulfidooxidans strain Cutipay enhances chalcopyrite bioleaching under moderate thermophilic conditions in the presence of chloride ion.

Currently more than 90% of the world's copper is obtained through sulfide mineral processing. Among the copper sulfides, chalcopyrite is the most abundant and therefore economically relevant. However, primary copper sulfide bioleaching is restricted due to high ionic strength raffinate solutions and particularly chloride coming from the dissolution of ores. In this work we describe the chalcopyrite bioleaching capacity of Sulfobacillus thermosulfidooxidans strain Cutipay (DSM 27601) previously described at the genomic level (Travisany et al. (2012) Draft genome sequence of the Sulfobacillus thermosulfidooxidans Cutipay strain, an indigenous bacterium isolated from a naturally extreme mining environment in Northern Chile. J Bacteriol 194:6327-6328). Bioleaching assays with the mixotrophic strain Cutipay showed a strong increase in copper recovery from chalcopyrite concentrate at 50°C in the presence of chloride ion, a relevant inhibitory element present in copper bioleaching processes. Compared to the abiotic control and a test with Sulfobacillus acidophilus DSM 10332, strain Cutipay showed an increase of 42 and 69% in copper recovery, respectively, demonstrating its high potential for chalcopyrite bioleaching. Moreover, a genomic comparison highlights the presence of the 2-Haloacid dehalogenase predicted-protein related to a potential new mechanism of chloride resistance in acidophiles. This novel and industrially applicable strain is under patent application CL 2013-03335.

Stoichiometric modeling of oxidation of reduced inorganic sulfur compounds (Riscs) in Acidithiobacillus thiooxidans.

The prokaryotic oxidation of reduced inorganic sulfur compounds (RISCs) is a topic of utmost importance from a biogeochemical and industrial perspective. Despite sulfur oxidizing bacterial activity is largely known, no quantitative approaches to biological RISCs oxidation have been made, gathering all the complex abiotic and enzymatic stoichiometry involved. Even though in the case of neutrophilic bacteria such as Paracoccus and Beggiatoa species the RISCs oxidation systems are well described, there is a lack of knowledge for acidophilic microorganisms. Here, we present the first experimentally validated stoichiometric model able to assess RISCs oxidation quantitatively in Acidithiobacillus thiooxidans (strain DSM 17318), the archetype of the sulfur oxidizing acidophilic chemolithoautotrophs. This model was built based on literature and genomic analysis, considering a widespread mix of formerly proposed RISCs oxidation models combined and evaluated experimentally. Thiosulfate partial oxidation by the Sox system (SoxABXYZ) was placed as central step of sulfur oxidation model, along with abiotic reactions. This model was coupled with a detailed stoichiometry of biomass production, providing accurate bacterial growth predictions. In silico deletion/inactivation highlights the role of sulfur dioxygenase as the main catalyzer and a moderate function of tetrathionate hydrolase in elemental sulfur catabolism, demonstrating that this model constitutes an advanced instrument for the optimization of At. thiooxidans biomass production with potential use in biohydrometallurgical and environmental applications.

Metabolomic study of Chilean biomining bacteria Acidithiobacillus ferrooxidans strain Wenelen and Acidithiobacillus thiooxidans strain Licanantay.

In this study, we present the first metabolic profiles for two bioleaching bacteria using capillary electrophoresis coupled with mass spectrometry. The bacteria, Acidithiobacillus ferrooxidans strain Wenelen (DSM 16786) and Acidithiobacillus thiooxidans strain Licanantay (DSM 17318), were sampled at different growth phases and on different substrates: the former was grown with iron and sulfur, and the latter with sulfur and chalcopyrite. Metabolic profiles were scored from planktonic and sessile states. Spermidine was detected in intra- and extracellular samples for both strains, suggesting it has an important role in biofilm formation in the presence of solid substrate. The canonical pathway for spermidine synthesis seems absent as its upstream precursor, putrescine, was not present in samples. Glutathione, a catalytic activator of elemental sulfur, was identified as one of the most abundant metabolites in the intracellular space in A. thiooxidans strain Licanantay, confirming its participation in the sulfur oxidation pathway. Amino acid profiles varied according to the growth conditions and bioleaching species. Glutamic and aspartic acid were highly abundant in intra- and extracellular extracts. Both are constituents of the extracellular matrix, and have a probable role in cell detoxification. This novel metabolomic information validates previous knowledge from in silico metabolic reconstructions based on genomic sequences, and reveals important biomining functions such as biofilm formation, energy management and stress responses. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11306-012-0443-3) contains supplementary material, which is available to authorized users.