Interaction networks for identifying coupled molecular processes in microbial communities
© Bosse et al. 2015
Received: 29 November 2014
Accepted: 3 July 2015
Published: 15 July 2015
Microbial communities adapt to environmental conditions for optimizing metabolic flux. Such adaption may include cooperative mechanisms eventually resulting in phenotypic observables as emergent properties that cannot be attributed to an individual species alone. Understanding the molecular basis of cross-species cooperation adds to utilization of microbial communities in industrial applications including metal bioleaching and bioremediation processes. With significant advancements in metagenomics the composition of microbial communities became amenable for integrative analysis on the level of entangled molecular processes involving more than one species, in turn offering a data matrix for analyzing the molecular basis of cooperative phenomena.
We present an analysis framework aligned with a dynamical hierarchies concept for unraveling emergent properties in microbial communities, and exemplify this approach for a co-culture setting of At. ferrooxidans and At. thiooxidans. This minimum microbial community demonstrates a significant increase in bioleaching efficiency compared to the activity of individual species, involving mechanisms of the thiosulfate, the polysulfide and the iron oxidation pathway.
Populating gene-centric data structures holding rich functional annotation and interaction information allows deriving network models at the functional level coupling energy production and transport processes of both microbial species. Applying a network segmentation approach on the interaction network of ortholog genes covering energy production and transport proposes a set of specific molecular processes of relevance in bioleaching. The resulting molecular process model essentially involves functionalities such as iron oxidation, nitrogen metabolism and proton transport, complemented by sulfur oxidation and nitrogen metabolism, as well as a set of ion transporter functionalities. At. ferrooxidans-specific genes embedded in the molecular model representation hold gene functions supportive for ammonia utilization as well as for biofilm formation, resembling key elements for effective chalcopyrite bioleaching as emergent property in the co-culture situation.
Analyzing the entangled molecular processes of a microbial community on the level of segmented, gene-centric interaction networks allows identification of core molecular processes and functionalities adding to our mechanistic understanding of emergent properties of microbial consortia.
KeywordsNetwork biology Microbial cooperation Bioleaching Chalcopyrite Acidithiobacillus Emergence
Microorganisms populate their environment not necessarily as solitude metabolic entities but may form metabolic communities comprising of several species. Such microbial communities are under synergistic and cooperative control with further constraints imposed by environmental factors. Under certain conditions such communities may exhibit emergent properties, i.e. phenotypic observables not being attributable to any single species but only to the community as such. The observation of mutual dependence of some species on the services and/or products of other community members is widely observed, but in many cases hampers comprehensive investigation at the level of an integrated molecular process analysis .
Metagenomics as core technology for species identification , further combined with omics on the transcript, protein and metabolite level allow to study molecular superpathways as a composite of the entire underlying molecular network contributed by each individual species. Such combined omics has provided unprecedented insights into the functioning, ecology and evolution of microbial communities [3–5]. Proper handling and integrative analysis of such molecular data space characterizing microbial communities has allowed embedding phenotypic observables into various theoretical frameworks as the evolutionary selection for cooperation , the public goods dilemma in consortia , or avoidance of evolutionary traps .
Further experimental advancements have boosted resolution of microbial community research, including procedures offering improved sensitivity to capture also low-abundant species , or an approach for allowing profiling novel microorganisms without being limited by gene/protein data available in reference sequences, in consequence significantly expanding the scope of metagenomics research . Together, such technologies have popularized metagenomics techniques finding applications in diverse areas such as human health, nutrition, industrial production, or environment remediation .
A number of industrial processes are based on utilization of microbial consortia. One prominent example is extraction of base metals from mineral ores, termed bioleaching, where the potential of bacteria-enabled biotechnologies has been successfully demonstrated [12–14]. In such systems, consortia of bacteria enhance the rate of metal recovery from ores containing valuable metals in the form of mineral sulfides, e.g. sphalerite (ZnS) or chalcopyrite (CuFeS2), or for recovery of gold encapsulated in pyrite (FeS2). The communities are composed of specialized, mostly mesoacidophilic and chemolithotrophic bacteria. Community composition is dynamic in response to changes in temperature, pH, concentration of ferric iron or toxic metals, or salinity .
At an industrial scale, bioleaching systems are used in recovery of metal ores (e.g. gold and silver) or to leach base metals such as copper, cobalt, nickel and zinc, as well as low-grade uranium ores. Extending the approach to non-base metals is not yet commercially attractive. The main reasons for the lack of commercial bioleaching operations are seen in efficiency constraints, costs associated with the construction of suitable bioreactors as well as in logistic issues related to vast volume of material to be processed . Recently, mining companies are confronted with a set of drivers that affect the economics of mining. Among those is a growing global demand for raw metals as well as tightening environmental legislation. Efficiency gains can be expected through the adoption and/or optimization of bioleaching processes .
Metagenomic data sets generated during the last decade have massively expanded the knowledge base on bioleaching consortia. First level analyses mainly focused on determining community composition and species population dynamics. Further analysis aimed at understanding key metabolic pathways involved in energy production, fixation of carbon and nitrogen, or bacterial survival in highly acidic and/or toxic environments [12, 18–20].
Still, from a systems biology perspective it is desirable to go beyond a parts-based approach and focus on cooperation from an entangled molecular process point of view [12, 21]. This approach is particularly important in view of lacking efficient genetic engineering platforms for major bioleaching microbes next to regulatory constraints . One promising option for achieving a more comprehensive understanding of collaborative functionalities is to approach bioleaching consortia like single ‘superorganisms’. For executing such approach, data are frequently organized in a network context, resembling molecular features as nodes and molecular interactions as edges. Such networks then undergo functional category- or topology-based analysis aimed at identifying molecular processes associated with phenotypic readouts of interest (utilized in various systems biology examples e.g. in human health [22, 23]). Functional insights at the molecular process level are to be expected if metadata on the various omics levels are integrated in a way that allows the linkage of specific community molecular processes to phenotypic observables becoming apparent only at the superorganism level. Such approach represents a shift from a species/individual molecular component point of view towards structuring of information across species boundaries. This approach retrieves molecular processes integrating molecular features from more than one species, by this providing molecular process models underlying cooperative phenomena.
From a formal perspective such phenomena have been described by Baas et al. as dynamical hierarchies , already applied in diverse chemical  and biological systems . The concept starts with first order objects (e.g. protein coding genes), based on interaction forming higher order structures (functional categories, pathways), which in turn exhibit properties not identified on any lower order level. As practical example serves biofilm formation, if combined further with thiosulfate, the polysulfide and the iron oxidation pathways resulting in yet another phenotypic readout, or emergent property as efficient leaching of minerals. Still, specific functional properties may be contributed from individual species of a microbial community, and a network representation of community functional properties in contrast to analysis on the level of individual species is deemed necessary for identification of community properties.
Romo et al.  described a bioleaching system which (i) exhibits an emergent community effect in terms of mineral leaching efficiency, (ii) is composed of a minimum number of species which (iii) have rich annotation at the gene and protein level. The authors investigated the efficiency of chalcopyrite (CuFeS2) bioleaching in a controlled setting with pure cultures of Gram-negative gamma-proteobacteria Acidithiobacillus ferrooxidans and At. thiooxidans compared to a co-culture of both bacteria, where the latter about doubled recovery rate of copper. Similar observations were reported in studies conducted by Fu et al.  and – albeit in an overall much less efficient system, by Qui et al. . A striking observation with these minimal systems is the fact that one of the species, At. ferrooxidans, alone possesses both sulfur and iron oxidizing properties deemed essential for bioleaching [13, 30]. Yet, a co-culture with At. thiooxidans, a sulfur- but not iron-oxidizer , results in a significant increase in leaching rate, raising the question about the nature of the cooperation of the two microbes.
In this work we implement a network biology methodology resting on functional annotation and interaction information of the protein coding gene space of both organisms aimed at identification of cooperation with significant increase in bioleaching as emergent property. We exemplify the approach on core molecular processes executed by At. thiooxidans and At. ferrooxidans in chalcopyrite bioleaching.
Gene-centric annotation on a species level utilized publicly available databases. As primary data source for annotation of At. ferrooxidans (ATCC 53993) and At. thiooxidans (ATCC 19377) the Pathosystems Resource Integration Center (PATRIC, ) with the respective PATRIC organism identifiers 6930 and 52187 was used. Further data retrieval utilized the Integrated Microbial Genomes database (IMG, ) using IMG taxon object identifiers 642788233 and 2510461056 for At. ferrooxidans and At. thiooxidans, respectively. In addition to gene annotation sourced from PATRIC and IMG, data from the Expasy family of gene/protein resources (HAMAP, ; Prosite ; PSORT ) as well as the EBI-EMBL resources InterPro  and UniProt  were obtained using the online tools provided on the respective sites.
Gene annotation and interaction details
gene coverage, At. thiooxidans
gene coverage, At. ferrooxidans
A: gene annotation
PATRIC database ID
RefSeq gene ID
Uniprot gene name
IMG database ID
RefSeq protein ID
Protein name & function
Enzyme Commission nomenclature
FIGfam protein family ID
Pfam protein family
Protein similarity (HAMAP)
Protein function & location
Interpro database entry
GO cellular component
PSORT subcellular location
B: interaction information
Interaction annotation data (Table 1B) were sourced from STRING , MetaCyc  and KEGG . MetaCyc interactions were inferred based on co-membership in MetaCyc pathways restricted to genes holding an Enzyme Commission number . To obtain KEGG interaction data KGML+ files containing KEGG pathway maps for At. ferrooxidans were downloaded from the KEGG website, and interactions were extracted using the Bioconductor KEGGraph package . STRING interaction data (At. ferrooxidans) were directly retrieved from the STRING web server. As no interaction data source hosts specific data for At. thiooxidans, interactions were assigned to At. thiooxidans genes based on ortholog mapping to At. ferrooxidans utilizing the Needle algorithm from the EMBOSS package . In case of multiple assignments the top scoring ortholog was selected. In total 1,661 ortholog genes are identified for the two species, allowing the imputation of interactions for the subset of orthologs identified on the At. thiooxidans side.
Further annotation specifically in the context of microbial consortia and bioleaching was performed via mining of NCBI Pubmed . A keyword catalog was applied for covering scientific publications on cooperative aspects (synergism, consortium, community, cooperation) and specific utilization (bioleaching, biomining, chalcopyrite). The query term included the gene name together with a respective catalog entry.
Gene annotation and interaction information as provided in Table 1 allows designing a gene-centric data structure for each of the At. ferrooxidans and At. thiooxidans genes using the PATRIC ID as unique structure identifier. The core data structure represents the protein coding gene annotation of both organisms, and further defines aggregate nodes by e.g. grouping genes according to gene assignment to clusters of orthologous groups (COG) . By this the set of gene-centric data structures is traversed into a set of COG-based data structures.
On the level of COG terms such networks represent co-annotation of genes in more than one COG term. Reciprocal assignment of genes to two COG terms infers functional closeness of such terms, expressed by an interaction. On a gene level interaction networks are derived by making use of gene interaction information. First, a network is built including all ortholog genes together with interaction information according to annotation. For identifying network segments with high inner connectivity segmentation was performed using MCODE with default settings . With a given network structure as start MCODE identifies most densely connected subgraphs, from there expanding to neighborhood network nodes with a cutoff criterion for stopping expansion according to the connectivity of such neighboring nodes. Specifically, starting with a given set of genes (input set) and given interaction information the induced subgraph is derived, eliminating all genes with a connectivity of zero. Network segments (units) are derived from this induced graph via MCODE, where each such unit is interpreted as functional context (molecular process). Utilizing interaction information of features across identified segments allows approximating functional dependencies between molecular processes, finally leading to a molecular process model representation of the input gene set. For segmentation and visualization of networks the MCODE plugin of Cytoscape was used (version 3.1.1, ).
Dynamical hierarchies in microbial communities
On this basis we assume explicit encoding of molecular functional properties of the microbial co-culture M (a set of microbial species m i ). In addition, we detect a significant bioleaching efficiency increase in the co-culture situation, which in the following is denoted as observable O 3 M . Such observable rests on properties of individual genes (O 1 ) and properties of aggregate states observed on the level of molecular pathways (O 2 ). Our analysis goal is identification of the set of first order objects, being the protein coding genes g i , in synthesis (interactome) leading to the phenotypic readout O 3 M . As in our specific scenario such readout is observable in the specific co-culture situation we assume involvement of first order objects (protein coding genes) from both, At. ferrooxidans and At. thiooxidans.
We know O 3 M := M, being a mapping of the set of two species (or n species in a general metagenomics setting). This mapping does not reveal any type of context structure, hence no functional f(M) allowing us to model O 3 M . This status is equivalent to the first result of a metagenomics analysis, i.e. the list of distinct species involved. M is a composition of a set of individual species m i , where each m i itself exhibits observables, O 2 m i , i.e. a readout of a species-specific pathway. Such observables are the consequence of the total set of molecular components g i j of m i . Describing O 2 m i as f(g i j ) having just a context-free listing of g j i provides no better system representation, as we end up with the mapping O 2 m i := g i j . When including interactions for allowing modeling of pathways such observable may reflect the capacity of At. ferrooxidans to oxidize both sulfur and iron, whereas At. thiooxidans shows sulfur oxidation only.
We further see G as the total set of molecular components of M. With G we break up the boundaries to individual m i , and obtain a mapping of O 3 M := G, which allows describing the community observable as a function of the community genome . Consequently, the observable of interest results from a functional f(G), where G is the superset of molecular constituents of both organisms, g i j , and when adding interaction information G becomes the superpathway. We postulate that the molecular process leading to our emergent property of interest is identified at the level of such superpathway.
According to our goal definition we seek O 3 M := f(g i ) out of G. Such (g i ) denotes a defined subset of G in a specific context, in a biological interpretation being individual molecular processes p i out of the total set of processes P. Consequently, the phenotype O 3 M of interest results from a functional f(P). Making use of the data structures centered around genes g i j specific p i leading to O 3 M may be derived on some aggregate level, e.g. COG categories resembling O 2 m i , or may be derived bottom up on the level of G.
Molecular processes embedded in pathway nomenclature allow deriving a model adding to understanding of improved bioleaching efficiency in the co-culture situation. However, such annotation is species-specific, i.e. on the level of O 2 m i , hence not providing an integrated representation of collaboration at a molecular process level.
Coupled molecular processes, aggregate mode
The COG category holding the highest number of non-ortholog genes is energy production and conversion, with 28 such genes on the At. ferrooxidans side, and 32 genes with At. thiooxidans, essentially resembling specific functions of iron oxidation and thiosulfate pathway.
COG terms and annotation details in microbial cooperation and bioleaching
(C) Energy production & conversion
(E) Amino acid transport & metabolism
(H) Coenzyme transport & metabolism
(O) Posttranslational modification, protein turnover & chaperones
(P) Inorganic ion transport & metabolism
(Q) Secondary metabolites biosynthesis, transport & catabolism
A fraction of non-ortholog genes are reported in the context of cooperative function, and to a lesser degree also in the very specific context of bioleaching. In the COG term energy production and conversion, both Acidithiobacilli hold Cytochrome B561 and NADH dehydrogenase in respect to cooperation. At. ferrooxidans further reports Aldo-keto reductase, Ferredoxin reductase, Ferredoxin-like protein, Ferredoxin and Aldehyde dehydrogenase. In the specific context of bioleaching Cytochrome C4 and Hydroxylamine reductase are found. In respect to COG term amino acid transport & metabolism Serine acetyltransferase, Homocitrate synthase, Cysteine desulfurase and Chorismate mutase record for cooperation. For At. thiooxidans Acetate kinase, Glycerol kinase, Oxidoreductase list for the COG term energy production and conversion in a cooperative context, the latter protein also in a specific utilization context. For the COG term amino acid transport and metabolism Glutamate decarboxylase and Transglutaminase-like protein are found in the context of cooperation.
Coupled molecular processes, first order mode
Leaching of chalcopyrite involves as main term energy production and conversion, on a pathway level expressed as thiosulfate, the polysulfide and the iron oxidation pathway being attributable on the individual species level m i . Exploring the involved gene sets on an aggregate mode of COG terms identifies via ortholog mapping further potential mechanisms of relevance in cooperative effects, specifically transport of e.g. coenzymes and amino acids. Ortholog analysis identified a set of species-specific genes g i , to some extend already holding annotation in the context of collaborative effects and bioleaching. According to our concept we aim for identifying O 3 M := f(g i ) out of G, where such g i constitute a molecular process p i seeing contributions from both species.
Focusing on the COG term energy production and conversion, additionally taking transport categories into account (secondary metabolites biosynthesis, transport and catabolism; inorganic ion transport and metabolism; coenzyme transport and metabolism; amino acid transport and metabolism) results in 386 ortholog genes of At. ferrooxidans and At. thiooxidans, and 195 genes specifically with At. ferrooxidans and 168 genes specifically with At. thiooxidans. Utilizing interaction data from STRING an ortholog network can be derived, i.e. a graph with nodes encoding the ortholog genes and edges representing the gene-specific interaction information. This ortholog graph holds 385 genes (1 gene of the input gene set had no single interaction information in STRING and is omitted from further analysis) and an Index of Aggregation of 1.0, i.e. each node holds a path to all other nodes of the network.
We consider the molecular units presented in Fig. 5 as individual molecular processes p i sharing genes across species boundaries, allowing us to postulate a mapping of O 3 M := f(p i ). In the given analysis the ortholog model involving COG terms energy production and conversion together with specific transport terms is used as core structure, adding molecular functionality coming from At. ferrooxidans only (lack of At. thiooxidans-specific interaction data).
Of specific interest are molecular processes integrating energy production and conversion and transport terms, hence, specifically molecular process units 1, 3 and 4. Unit 1 aggregates genes associated with proton transport, nitrogen metabolism and iron oxidation (e.g. NADH ubiquinone oxidoreductase and Glutamate synthase). Generally enriched in energy production is unit 3 hosting several Cytochrome O ubiquinol oxidases, ATP synthases as well as dehydrogenases such as Malate, Isocitrate and Pyruvate dehydrogenases. Functional links to sulfur oxidation is provided in unit 4 (Rhodanese-like protein, Cytochrome D and O ubiquinol odixases). Unit 6 exclusively covers aspects of nitrogen metabolism (Nitrogen regulating protein P-II, Ammonium transporter). In the context of bioleaching, metal ion transporters are of particular importance. Those are mainly embedded in unit 2 (magnesium and cobalt transport proteins as well as lead, cadmium, zinc and copper transport ATPases) and unit 7 (zinc transporters, ferric iron uptake regulators, ferrous iron transporters). Given interaction data result in multiple assignments of genes to process units. One major representative seeing such multiple assignments is Carbonic anhydrase, an enzyme involved in the sulfur cycle.
Integrative analysis on the grounds of omics data promises improved characterization of molecular processes mediated by microbial communities. Of particular interest are processes involving direct dependencies across species boundaries resembling a major driver for exhibiting emergent properties. A formal concept allowing representation of such phenomena is dynamical hierarchies, which are composed of individual objects holding a data structure encoding properties and interactions. In a composition individual objects form higher order structures which may direct to emergent properties (phenotypic observables) that are not amenable at any lower level hierarchy, and specifically not reducible to a single object. In the context of molecular systems this concept can be applied to sets of genes as first order objects, potentially leading to an emergent property based on a molecular process, which in turn is composed of a specific set of genes. This concept is directly applicable to microbial communities showing synergistic readout, with the example of significantly increased mineral leaching in a co-culture of At. ferrooxidans and At. thiooxidans.
A necessary prerequisite for implementing such approach is rich gene annotation covering from domain and function assignment to aggregate properties as functional category or pathway membership. A second prerequisite is interaction information implicating dependencies for executing molecular processes, which in combination with gene annotation allows setting up a data structure for building interaction networks on the gene level or any aggregate level as COG terms.
Utilizing species-specific pathway information allows deriving a first mechanistic model aimed at describing phenomena adding to improved bioleaching in the co-culture situation, essentially involving the thiosulfate, the polysulfide and the iron oxidation pathway. These constituents improve mineral dissolution from chalcopyrite via proton attack, sulfur layer removal, and hindrance of jarosite formation. Analyzing involved genes on the level of COG terms and ortholog mapping links the central term energy production and conversion with specific transport processes, namely secondary metabolites biosynthesis, transport and catabolism, inorganic ion transport and metabolism, coenzyme transport and metabolism, and amino acid transport and metabolism. This term set includes a number of non-ortholog genes from both species, a subset being already discussed in the context of microbial communities and bioleaching.
Certainly, a description at the level of COG terms does not allow delineation of integrated molecular processes. Building gene-centric interaction networks followed by segmentation according to topological criteria provides means for approximating such processes. Executing this procedure for the ortholog network of At. ferrooxidans and At. thiooxidans with focus on energy production and conversion coupled with identified transport terms results in a molecular process model bridging energy production and transport processes in shared molecular units. Integrated processes identified include iron oxidation, nitrogen metabolism and proton transport, segments covering sulfur oxidation and nitrogen metabolism, and a set of ion transporter functionalities. Including At. ferrooxidans-specific genes in this network results in a composite representation of shared community functionality and specific add-ons of one species. Of particular interest is Carbonic anhydrase, described as a constituent to the polysulfide pathway . This gene has been found as associated with biofilm formation in a very recent proteomics study . While the exact contributions of Carbonic anhydrase to a bioleaching process are unknown, three mechanisms can be conceived, including (i) positive impact on ammonia utilization of other consortia members by removal of environmental carbon disulfide (CS2) , (ii) enabling At. ferrooxidans to utilize CS2 as alternative energy source and (iii) contributing to inorganic carbon fixation . A second gene, Homocitrate synthase, is described as an essential component of lysine biosynthesis and has been found upregulated in a screen searching for quorum-sensing-related molecules , suggesting a role in cell-density-dependent processes implicated in the formation of biofilms.
Taken together, the additional genes contributed by At. ferrooxidans enable a cornerstone role in the establishment and perpetuation of a chalcopyrite leaching consortium. This is consistent with the results of a recently conducted metabolomic study using a laboratory community composed of At. thiooxidans and ferrooxidans strains. This study identified biofilm formation via At. ferrooxidans, suggesting a central role in community establishment, as consequence fostering increased bioleaching activity as emergent property .
A set of genes is identified on a common mechanistic background of At. ferrooxidans and At. thiooxidans, as these are components of molecular processes delineated from the ortholog network covering energy production and conversion together with transport processes. Recalling the concept of dynamical hierarchies such processes provide an integrated perspective of the superpathway composed of the individual constituents of both species together with their interactions.
The approach presented in this work allows deriving hypotheses on coupled molecular processes, providing an alternative perspective of synergistic effects and community-based emergent properties. With the concept of dynamical hierarchies established, availability of annotation and interaction information in the public domain, and open source libraries for computing networks and segments this approach is of general applicability. A major present pitfall is consolidated availability of molecular data needed to populate the data structures. Necessary gene annotation is distributed across numerous individual sources, hampering integration for given annotation beyond the challenge of adding annotation to novel species identified in metagenomics. Utilizing orthology for annotation imputation naturally results in biases, becoming even more pronounced when going to a larger set of species composing the community. Of equal relevance, and with even less coverage presents interaction information, being either generic on a pathway level, or only covering selected species. For fully leveraging on the power of integrative, network-based analysis procedures realizing common grounds on annotation and interaction information represents a major challenge ahead in the field, partly addressed by hybrid approaches in approximating interactions as done with STRING.
Delineation of coupled molecular processes associated with emergent properties, as synergistic increase in bioleaching efficiency, offers improved understanding of microbial community dynamics. Further, such molecular process knowledge provides the basis for selection of molecular biomarker candidates indicative for such community molecular function. Biomarkers in their definition serve as proxy for monitoring the status of a molecular process, in a microbial community needing to capture molecular processes established across individual species. Segmented, cross-species interaction networks provide coupled process sets as basis for selection of biomarker candidate panels. Such panels in turn can be utilized as tool for screening microbial communities prone to exhibiting an emergent property of interest.
This work was made possible as part of the META-BASE project executed by Vale Institute of Technology, Brazil.
- Karunakaran E, Mukherjee J, Ramalingam B, Biggs CA. "Biofilmology": a multidisciplinary review of the study of microbial biofilms. Applied microbiology and biotechnology. 2011;90:1869–81.View ArticlePubMedGoogle Scholar
- Metzker ML. Sequencing technologies - the next generation. Nature reviews Genetics. 2010;11:31–46.View ArticlePubMedGoogle Scholar
- Handelsman J, Rondon MR, Brady SF, Clardy J, Goodman RM. Molecular biological access to the chemistry of unknown soil microbes: a new frontier for natural products. Chemistry & biology. 1998;5:R245–249.View ArticleGoogle Scholar
- Streit WR, Schmitz RA. Metagenomics--the key to the uncultured microbes. Current opinion in microbiology. 2004;7:492–8.View ArticlePubMedGoogle Scholar
- Thomas T, Gilbert J, Meyer F. Metagenomics - a guide from sampling to data analysis. Microbial informatics and experimentation. 2012;2:3.View ArticlePubMedPubMed CentralGoogle Scholar
- Damore JA, Gore J. Understanding microbial cooperation. Journal of theoretical biology. 2012;299:31–41.View ArticlePubMedGoogle Scholar
- Drescher K, Nadell CD, Stone HA, Wingreen NS, Bassler BL. Solutions to the public goods dilemma in bacterial biofilms. Current biology : CB. 2014;24:50–5.View ArticlePubMedGoogle Scholar
- Cornforth DM, Sumpter DJ, Brown SP, Brannstrom A. Synergy and group size in microbial cooperation. The American naturalist. 2012;180:296–305.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang Y, Leung HC, Yiu SM, Chin FY. MetaCluster 5.0: a two-round binning approach for metagenomic data for low-abundance species in a noisy sample. Bioinformatics (Oxford, England). 2012;28:i356–62.View ArticleGoogle Scholar
- Nielsen HB, Almeida M, Juncker AS, Rasmussen S, Li J, Sunagawa S, et al. Identification and assembly of genomes and genetic elements in complex metagenomic samples without using reference genomes. Nature biotechnology. 2014;32:822–8.View ArticlePubMedGoogle Scholar
- Logares R, Haverkamp TH, Kumar S, Lanzen A, Nederbragt AJ, Quince C, et al. Environmental microbiology through the lens of high-throughput DNA sequencing: synopsis of current platforms and bioinformatics approaches. Journal of microbiological methods. 2012;91:106–13.View ArticlePubMedGoogle Scholar
- Valenzuela L, Chi A, Beard S, Orell A, Guiliani N, Shabanowitz J, et al. Genomics, metagenomics and proteomics in biomining microorganisms. Biotechnology advances. 2006;24:197–211.View ArticlePubMedGoogle Scholar
- Rawlings DE, Johnson DB. The microbiology of biomining: development and optimization of mineral-oxidizing microbial consortia. Microbiology (Reading, England). 2007;153:315–24.View ArticleGoogle Scholar
- Johnson DB. Biomining-biotechnologies for extracting and recovering metals from ores and waste materials. Current opinion in biotechnology. 2014;30C:24–31.View ArticleGoogle Scholar
- Zhang RB, Wei MM, Ji HG, Chen XH, Qiu GZ, Zhou HB. Application of real-time PCR to monitor population dynamics of defined mixed cultures of moderate thermophiles involved in bioleaching of chalcopyrite. Applied microbiology and biotechnology. 2009;81:1161–8.View ArticlePubMedGoogle Scholar
- Schippers A, Sand W. Bacterial leaching of metal sulfides proceeds by two indirect mechanisms via thiosulfate or via polysulfides and sulfur. Applied and environmental microbiology. 1999;65:319–21.PubMedPubMed CentralGoogle Scholar
- Vera M, Schippers A, Sand W. Progress in bioleaching: fundamentals and mechanisms of bacterial metal sulfide oxidation--part A. Applied microbiology and biotechnology. 2013;97:7529–41.View ArticlePubMedGoogle Scholar
- Spolaore P, Joulian C, Gouin J, Morin D, d'Hugues P. Relationship between bioleaching performance, bacterial community structure and mineralogy in the bioleaching of a copper concentrate in stirred-tank reactors. Applied microbiology and biotechnology. 2011;89:441–8.View ArticlePubMedGoogle Scholar
- Remonsellez F, Galleguillos F, Moreno-Paz M, Parro V, Acosta M, Demergasso C. Dynamic of active microorganisms inhabiting a bioleaching industrial heap of low-grade copper sulfide ore monitored by real-time PCR and oligonucleotide prokaryotic acidophile microarray. Microbial biotechnology. 2009;2:613–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang Y, Zeng W, Qiu G, Chen X, Zhou H. A moderately thermophilic mixed microbial culture for bioleaching of chalcopyrite concentrate at high pulp density. Applied and environmental microbiology. 2014;80:741–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Trigo A, Valencia A, Cases I. Systemic approaches to biodegradation. FEMS microbiology reviews. 2009;33:98–108.View ArticlePubMedGoogle Scholar
- Heinzel A, Perco P, Mayer G, Oberbauer R, Lukas A, Mayer B. From molecular signatures to predictive biomarkers: modeling disease pathophysiology and drug mechanism of action. Frontiers in cell and developmental biology. 2014;2:37.View ArticlePubMedPubMed CentralGoogle Scholar
- Heinzel A, Muhlberger I, Fechete R, Mayer B, Perco P. Functional molecular units for guiding biomarker panel design. Methods in molecular biology (Clifton, NJ). 2014;1159:109–33.View ArticleGoogle Scholar
- Baas NA, Seeman NC. On the chemical synthesis of new topological structures. Journal of mathematical chemistry. 2012;50:220–32.View ArticlePubMedGoogle Scholar
- Rasmussen S, Baas NA, Mayer B, Nilsson M, Olesen MW. Ansatz for dynamical hierarchies. Artificial life. 2001;7:329–53.View ArticlePubMedGoogle Scholar
- Siehs C, Oberbauer R, Mayer G, Lukas A, Mayer B. Discrete simulation of regulatory homo- and heterodimerization in the apoptosis effector phase. Bioinformatics (Oxford, England). 2002;18:67–76.View ArticleGoogle Scholar
- Romo E, Weinacker DF, Zepeda AB, Figueroa CA, Chavez-Crooker P, Farias JG. Bacterial consortium for copper extraction from sulphide ore consisting mainly of chalcopyrite. Brazilian journal of microbiology : [publication of the Brazilian Society for Microbiology]. 2013;44:523–8.View ArticleGoogle Scholar
- Fu B, Zhou H, Zhang R, Qiu G. Bioleaching of chalcopyrite by pure and mixed cultures of Acidithiobacillus spp. and Leptospirillum ferriphilum. International Biodeterioration & Biodegradation. 2008;62:109–15.View ArticleGoogle Scholar
- Qiu MQ, Xiong SY, Zhang WM, Wang GX. A comparison of bioleaching of chalcopyrite using pure culture or a mixed culture. Minerals Engineering. 2005;18:987–90.View ArticleGoogle Scholar
- Valdes J, Pedroso I, Quatrini R, Dodson RJ, Tettelin H, Blake 2nd R, et al. Acidithiobacillus ferrooxidans metabolism: from genome sequence to industrial applications. BMC genomics. 2008;9:597.View ArticlePubMedPubMed CentralGoogle Scholar
- Valdes J, Ossandon F, Quatrini R, Dopson M, Holmes DS. Draft genome sequence of the extremely acidophilic biomining bacterium Acidithiobacillus thiooxidans ATCC 19377 provides insights into the evolution of the Acidithiobacillus genus. Journal of bacteriology. 2011;193:7003–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Wattam AR, Abraham D, Dalay O, Disz TL, Driscoll T, Gabbard JL, et al. PATRIC, the bacterial bioinformatics database and analysis resource. Nucleic acids research. 2014;42:D581–591.View ArticlePubMedGoogle Scholar
- Markowitz VM, Chen IM, Palaniappan K, Chu K, Szeto E, Grechkin Y, et al. IMG: the Integrated Microbial Genomes database and comparative analysis system. Nucleic acids research. 2012;40:D115–122.View ArticlePubMedGoogle Scholar
- Pedruzzi I, Rivoire C, Auchincloss AH, Coudert E, Keller G, de Castro E, et al. HAMAP in 2013, new developments in the protein family classification and annotation system. Nucleic acids research. 2013;41:D584–589.View ArticlePubMedGoogle Scholar
- Sigrist CJ, Cerutti L, Hulo N, Gattiker A, Falquet L, Pagni M, et al. PROSITE: a documented database using patterns and profiles as motif descriptors. Briefings in bioinformatics. 2002;3:265–74.View ArticlePubMedGoogle Scholar
- Yu NY, Wagner JR, Laird MR, Melli G, Rey S, Lo R, et al. PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes. Bioinformatics (Oxford, England). 2010;26:1608–15.View ArticleGoogle Scholar
- Hunter S, Jones P, Mitchell A, Apweiler R, Attwood TK, Bateman A, et al. InterPro in 2011: new developments in the family and domain prediction database. Nucleic acids research. 2012;40:D306–312.View ArticlePubMedGoogle Scholar
- UniProt-Consortium. Ongoing and future developments at the Universal Protein Resource. Nucleic acids research. 2011;39:D214–219.View ArticleGoogle Scholar
- UniProt ID Mapping Tool [http://www.uniprot.org/uploadlists/]
- Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, et al. STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic acids research. 2013;41:D808–815.View ArticlePubMedGoogle Scholar
- Caspi R, Altman T, Billington R, Dreher K, Foerster H, Fulcher CA, et al. The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of Pathway/Genome Databases. Nucleic acids research. 2014;42:D459–471.View ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S, Sato Y, Kawashima M, Furumichi M, Tanabe M. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic acids research. 2014;42:D199–205.View ArticlePubMedGoogle Scholar
- Biology NCotIUoBaM: Enzyme Nomenclature 1992. Academic Press; 1992.Google Scholar
- Zhang JD, Wiemann S. KEGGgraph: a graph approach to KEGG PATHWAY in R and bioconductor. Bioinformatics (Oxford, England). 2009;25:1470–1.View ArticleGoogle Scholar
- Rice P, Longden I, Bleasby A. EMBOSS: the European Molecular Biology Open Software Suite. Trends in genetics : TIG. 2000;16:276–7.View ArticlePubMedGoogle Scholar
- NCBI PubMed [http://www.ncbi.nlm.nih.gov/pubmed/]
- Tatusov RL, Koonin EV, Lipman DJ. A genomic perspective on protein families. Science (New York, NY). 1997;278:631–7.View ArticleGoogle Scholar
- Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC bioinformatics. 2003;4:2.View ArticlePubMedPubMed CentralGoogle Scholar
- Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics (Oxford, England). 2011;27:431–2.View ArticleGoogle Scholar
- Zhang W, Zeng E, Liu D, Jones SE, Emrich S. Mapping genomic features to functional traits through microbial whole genome sequences. International journal of bioinformatics research and applications. 2014;10:461–78.View ArticlePubMedGoogle Scholar
- Yin H, Zhang X, Li X, He Z, Liang Y, Guo X, et al. Whole-genome sequencing reveals novel insights into sulfur oxidation in the extremophile Acidithiobacillus thiooxidans. BMC microbiology. 2014;14:179.View ArticlePubMedPubMed CentralGoogle Scholar
- Xia L-x, Tang L, Xia J-l, Yin C, Cai L-y, Zhao X-j, et al. Relationships among bioleaching performance, additional elemental sulfur, microbial population dynamics and its energy metabolism in bioleaching of chalcopyrite. Transactions of Nonferrous Metals Society of China. 2012;22:192–8.View ArticleGoogle Scholar
- Fowler TA, Crundwell FK. Leaching of zinc sulfide by Thiobacillus ferrooxidans: bacterial oxidation of the sulfur product layer increases the rate of zinc sulfide dissolution at high concentrations of ferrous ions. Applied and environmental microbiology. 1999;65:5285–92.PubMedPubMed CentralGoogle Scholar
- Fanning DS, Rabenhorst MC, Burch SN, Islam KR, Tangren SA: Sulfides and Sulfates. In Soil Mineralogy with Environmental Applications. Edited by Dixon JB, Schulze DG: Soil Science Society of America; 2002: 229–260: SSSA Book Series].Google Scholar
- Daoud J, Karamanev D. Formation of jarosite during Fe2+ oxidation by Acidithiobacillus ferrooxidans. Minerals Engineering. 2006;19:960–7.View ArticleGoogle Scholar
- Smeulders MJ, Pol A, Venselaar H, Barends TR, Hermans J, Jetten MS, et al. Bacterial CS2 hydrolases from Acidithiobacillus thiooxidans strains are homologous to the archaeal catenane CS2 hydrolase. Journal of bacteriology. 2013;195:4046–56.View ArticlePubMedPubMed CentralGoogle Scholar
- Vera M, Krok B, Bellenberg S, Sand W, Poetsch A. Shotgun proteomics study of early biofilm formation process of Acidithiobacillus ferrooxidans ATCC 23270 on pyrite. Proteomics. 2013;13:1133–44.View ArticlePubMedGoogle Scholar
- McCarty GW. Modes of action of nitrification inhibitors. Biology and Fertility of Soils. 1999;29:1–9.View ArticleGoogle Scholar
- Auernik KS, Kelly RM. Physiological versatility of the extremely thermoacidophilic archaeon Metallosphaera sedula supported by transcriptomic analysis of heterotrophic, autotrophic, and mixotrophic growth. Applied and environmental microbiology. 2010;76:931–5.View ArticlePubMedGoogle Scholar
- Williamson LL, Borlee BR, Schloss PD, Guan C, Allen HK, Handelsman J. Intracellular screen to identify metagenomic clones that induce or inhibit a quorum-sensing biosensor. Applied and environmental microbiology. 2005;71:6335–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Martinez P, Galvez S, Ohtsuka N, Budinich M, Cortes MP, Serpell C, et al. Metabolomic study of Chilean biomining bacteria Acidithiobacillus ferrooxidans strain Wenelen and Acidithiobacillus thiooxidans strain Licanantay. Metabolomics : Official journal of the Metabolomic Society. 2013;9:247–57.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.