Comparative genomics sheds light on transcription factor-mediated regulation in the extreme acidophilic Acidithiobacillia representatives

Extreme acidophiles thrive in acidic environments, confront a multitude of challenges, and demonstrate remarkable adaptability in their metabolism to cope with the ever-changing environmental fluctuations, which encompass variations in temperature, pH levels, and the availability of electron acceptors and donors. The survival and proliferation of members within the Acidithiobacillia class rely on the deployment of transcriptional regulatory systems linked to essential physiological traits. The study of these transcriptional regulatory systems provides valuable insights into critical processes, such as energy metabolism and nutrient assimilation, and how they integrate into major genetic-metabolic circuits. In this study, we examined the transcriptional regulatory repertoires and potential interactions of forty-three Acidithiobacillia complete and draft genomes, encompassing nine species. To investigate the function and diversity of Transcription Factors (TFs) and their DNA Binding Sites (DBSs), we conducted a genome-wide comparative analysis, which allowed us to identify these regulatory elements in representatives of Acidithiobacillia. We classified TFs into gene families and compared their occurrence among all representatives, revealing conservation patterns across the class. The results identified conserved regulators for several pathways, including iron and sulfur oxidation, the main pathways for energy acquisition, providing new evidence for viable regulatory interactions and branch-specific conservation in Acidithiobacillia. The identification of TFs and DBSs not only corroborates existing experimental information for selected species, but also introduces novel candidates for experimental validation. Moreover, these promising candidates have the potential for further extension to new representatives within the class.

Eurypsychrophilic acidophiles: From (meta)genomes to low-temperature biotechnologies

Low temperature and acidic environments encompass natural milieus such as acid rock drainage in Antarctica and anthropogenic sites including drained sulfidic sediments in Scandinavia. The microorganisms inhabiting these environments include polyextremophiles that are both extreme acidophiles (defined as having an optimum growth pH < 3), and eurypsychrophiles that grow at low temperatures down to approximately 4°C but have an optimum temperature for growth above 15°C. Eurypsychrophilic acidophiles have important roles in natural biogeochemical cycling on earth and potentially on other planetary bodies and moons along with biotechnological applications in, for instance, low-temperature metal dissolution from metal sulfides. Five low-temperature acidophiles are characterized, namely, Acidithiobacillus ferriphilus, Acidithiobacillus ferrivorans, Acidithiobacillus ferrooxidans, “Ferrovum myxofaciens,” and Alicyclobacillus disulfidooxidans, and their characteristics are reviewed. Our understanding of characterized and environmental eurypsychrophilic acidophiles has been accelerated by the application of “omics” techniques that have aided in revealing adaptations to low pH and temperature that can be synergistic, while other adaptations are potentially antagonistic. The lack of known acidophiles that exclusively grow below 15°C may be due to the antagonistic nature of adaptations in this polyextremophile. In conclusion, this review summarizes the knowledge of eurypsychrophilic acidophiles and places the information in evolutionary, environmental, biotechnological, and exobiology perspectives.

Genome-guided prediction of acid resistance mechanisms in acidophilic methanotrophs of phylogenetically deep-rooted Verrucomicrobia isolated from geothermal environments

Verrucomicrobia are a group of microorganisms that have been proposed to be deeply rooted in the Tree of Life. Some are methanotrophs that oxidize the potent greenhouse gas methane and are thus important in decreasing atmospheric concentrations of the gas, potentially ameliorating climate change. They are widespread in various environments including soil and fresh or marine waters. Recently, a clade of extremely acidophilic Verrucomicrobia, flourishing at pH < 3, were described from high-temperature geothermal ecosystems. This novel group could be of interest for studies about the emergence of life on Earth and to astrobiologists as homologs for possible extraterrestrial life. In this paper, we describe predicted mechanisms for survival of this clade at low pH and suggest its possible evolutionary trajectory from an inferred neutrophilic ancestor. Extreme acidophiles are defined as organisms that thrive in extremely low pH environments (≤ pH 3). Many are polyextremophiles facing high temperatures and high salt as well as low pH. They are important to study for both providing fundamental insights into biological mechanisms of survival and evolution in such extreme environments and for understanding their roles in biotechnological applications such as industrial mineral recovery (bioleaching) and mitigation of acid mine drainage. They are also, potentially, a rich source of novel genes and pathways for the genetic engineering of microbial strains. Acidophiles of the Verrucomicrobia phylum are unique as they are the only known aerobic methanotrophs that can grow optimally under acidic (pH 2–3) and moderately thermophilic conditions (50–60°C). Three moderately thermophilic genera, namely Methylacidiphilum, Methylacidimicrobium, and Ca. Methylacidithermus, have been described in geothermal environments. Most of the investigations of these organisms have focused on their methane oxidizing capabilities (methanotrophy) and use of lanthanides as a protein cofactor, with no extensive study that sheds light on the mechanisms that they use to flourish at extremely low pH. In this paper, we extend the phylogenetic description of this group of acidophiles using whole genome information and we identify several mechanisms, potentially involved in acid resistance, including “first line of defense” mechanisms that impede the entry of protons into the cell. These include the presence of membrane-associated hopanoids, multiple copies of the outer membrane protein (Slp), and inner membrane potassium channels (kup, kdp) that generate a reversed membrane potential repelling the intrusion of protons. Acidophilic Verrucomicrobia also display a wide array of proteins potentially involved in the “second line of defense” where protons that evaded the first line of defense and entered the cell are expelled or neutralized, such as the glutamate decarboxylation (gadAB) and phosphate-uptake systems. An exclusive N-type ATPase F₀-F₁ was identified only in acidophiles of Verrucomicrobia and is predicted to be a specific adaptation in these organisms. Phylogenetic analyses suggest that many predicted mechanisms are evolutionarily conserved and most likely entered the acidophilic lineage of Verrucomicrobia by vertical descent from a common ancestor. However, it is likely that some defense mechanisms such as gadA and kup entered the acidophilic Verrucomicrobia lineage by horizontal gene transfer.

A Large-Scale Genome-Based Survey of Acidophilic Bacteria Suggests That Genome Streamlining Is an Adaption for Life at Low pH

The genome streamlining theory suggests that reduction of microbial genome size optimizes energy utilization in stressful environments. Although this hypothesis has been explored in several cases of low-nutrient (oligotrophic) and high-temperature environments, little work has been carried out on microorganisms from low-pH environments, and what has been reported is inconclusive. In this study, we performed a large-scale comparative genomics investigation of more than 260 bacterial high-quality genome sequences of acidophiles, together with genomes of their closest phylogenetic relatives that live at circum-neutral pH. A statistically supported correlation is reported between reduction of genome size and decreasing pH that we demonstrate is due to gene loss and reduced gene sizes. This trend is independent from other genome size constraints such as temperature and G + C content. Genome streamlining in the evolution of acidophilic bacteria is thus supported by our results. The analyses of predicted Clusters of Orthologous Genes (COG) categories and subcellular location predictions indicate that acidophiles have a lower representation of genes encoding extracellular proteins, signal transduction mechanisms, and proteins with unknown function but are enriched in inner membrane proteins, chaperones, basic metabolism, and core cellular functions. Contrary to other reports for genome streamlining, there was no significant change in paralog frequencies across pH. However, a detailed analysis of COG categories revealed a higher proportion of genes in acidophiles in the following categories: "replication and repair," "amino acid transport," and "intracellular trafficking". This study brings increasing clarity regarding the genomic adaptations of acidophiles to life at low pH while putting elements, such as the reduction of average gene size, under the spotlight of streamlining theory.

Integrative Genomics Sheds Light on Evolutionary Forces Shaping the Acidithiobacillia Class Acidophilic Lifestyle

Extreme acidophiles thrive in environments rich in protons (pH values <3) and often high levels of dissolved heavy metals. They are distributed across the three domains of the Tree of Life including members of the Proteobacteria. The Acidithiobacillia class is formed by the neutrophilic genus Thermithiobacillus along with the extremely acidophilic genera Fervidacidithiobacillus, Igneacidithiobacillus, Ambacidithiobacillus, and Acidithiobacillus. Phylogenomic reconstruction revealed a division in the Acidithiobacillia class correlating with the different pH optima that suggested that the acidophilic genera evolved from an ancestral neutrophile within the Acidithiobacillia. Genes and mechanisms denominated as "first line of defense" were key to explaining the Acidithiobacillia acidophilic lifestyle including preventing proton influx that allows the cell to maintain a near-neutral cytoplasmic pH and differ from the neutrophilic Acidithiobacillia ancestors that lacked these systems. Additional differences between the neutrophilic and acidophilic Acidithiobacillia included the higher number of gene copies in the acidophilic genera coding for "second line of defense" systems that neutralize and/or expel protons from cell. Gain of genes such as hopanoid biosynthesis involved in membrane stabilization at low pH and the functional redundancy for generating an internal positive membrane potential revealed the transition from neutrophilic properties to a new acidophilic lifestyle by shaping the Acidithiobacillaceae genomic structure. The presence of a pool of accessory genes with functional redundancy provides the opportunity to "hedge bet" in rapidly changing acidic environments. Although a core of mechanisms for acid resistance was inherited vertically from an inferred neutrophilic ancestor, the majority of mechanisms, especially those potentially involved in resistance to extremely low pH, were obtained from other extreme acidophiles by horizontal gene transfer (HGT) events.

A Large-Scale Multiple Genome Comparison of Acidophilic Archaea (pH ≤ 5.0) Extends Our Understanding of Oxidative Stress Responses in Polyextreme Environments

Acidophilic archaea thrive in anaerobic and aerobic low pH environments (pH < 5) rich in dissolved heavy metals that exacerbate stress caused by the production of reactive oxygen species (ROS) such as hydrogen peroxide (H2O2), hydroxyl radical (OH) and superoxide (O2−). ROS react with lipids, proteins and nucleic acids causing oxidative stress and damage that can lead to cell death. Herein, genes and mechanisms potentially involved in ROS mitigation are predicted in over 200 genomes of acidophilic archaea with sequenced genomes. These organisms are often be subjected to simultaneous multiple stresses such as high temperature, high salinity, low pH and high heavy metal loads. Some of the topics addressed include: (1) the phylogenomic distribution of these genes and what this can tell us about the evolution of these mechanisms in acidophilic archaea; (2) key differences in genes and mechanisms used by acidophilic versus non-acidophilic archaea and between acidophilic archaea and acidophilic bacteria and (3) how comparative genomic analysis predicts novel genes or pathways involved in oxidative stress responses in archaea and likely horizontal gene transfer (HGT) events.

SinEx DB 2.0 update 2020: database for eukaryotic single-exon coding sequences

Single-exon coding sequences (CDSs), also known as ‘single-exon genes’ (SEGs), are defined as nuclear, protein-coding genes that lack introns in their CDSs. They have been studied not only to determine their origin and evolution but also because their expression has been linked to several types of human cancers and neurological/developmental disorders, and many exhibit tissue-specific transcription. We developed SinEx DB that houses DNA and protein sequence information of SEGs from 10 mammalian genomes including human. SinEx DB includes their functional predictions (KOG (euKaryotic Orthologous Groups)) and the relative distribution of these functions within species. Here, we report SinEx 2.0, a major update of SinEx DB that includes information of the occurrence, distribution and functional prediction of SEGs from 60 completely sequenced eukaryotic genomes, representing animals, fungi, protists and plants. The information is stored in a relational database built with MySQL Server 5.7, and the complete dataset of SEG sequences and their GO (Gene Ontology) functional assignations are available for downloading. SinEx DB 2.0 was built with a novel pipeline that helps disambiguate single-exon isoforms from SEGs. SinEx DB 2.0 is the largest available database for SEGs and provides a rich source of information for advancing our understanding of the evolution, function of SEGs and their associations with disorders including cancers and neurological and developmental diseases.

Database URL:http://v2.sinex.cl/

AciDB 1.0: a database of acidophilic organisms, their genomic information and associated metadata

MOTIVATION

There are about 600 available genome sequences of acidophilic organisms (grow at a pH < 5) from the three domains of the Tree of Life. Information about acidophiles is scattered over many heterogeneous sites making it extraordinarily difficult to link physiological traits with genomic data. We were motivated to generate a curated, searchable database to address this problem.

RESULTS

AciDB 1.0 is a curated database of sequenced acidophiles that enables researchers to execute complex queries linking genomic features to growth data, environmental descriptions and taxonomic information.

AVAILABILITY AND IMPLEMENTATION

AciDB 1.0 is freely available online at: http://AciDB.cl. The source code is released under an MIT license at: https://gitlab.com/Hawkline451/acidb/.

Genome-based classification of Acidihalobacter prosperus F5 (=DSM 105917=JCM 32255) as Acidihalobacter yilgarnensis sp. nov

The genus Acidihalobacter has three validated species, Acidihalobacter ferrooxydans, Acidihalobacter prosperus and Acidihalobacter aeolinanus, all of which were isolated from Vulcano island, Italy. They are obligately chemolithotrophic, aerobic, acidophilic and halophilic in nature and use either ferrous iron or reduced sulphur as electron donors. Recently, a novel strain was isolated from an acidic, saline drain in the Yilgarn region of Western Australia. Strain F5T has an absolute requirement for sodium chloride (>5 mM) and is osmophilic, growing in elevated concentrations (>1 M) of magnesium sulphate. A defining feature of its physiology is its ability to catalyse the oxidative dissolution of the most abundant copper mineral, chalcopyrite, suggesting a potential role in biomining. Originally categorized as a strain of A. prosperus, 16S rRNA gene phylogeny and multiprotein phylogenies derived from clusters of orthologous proteins (COGS) of ribosomal protein families and universal protein families unambiguously demonstrate that strain F5T forms a well-supported separate branch as a sister clade to A. prosperus and is clearly distinguishable from A. ferrooxydans DSM 14175T and A. aeolinanus DSM14174T. Results of comparisons between strain F5T and the other Acidihalobacter species, using genome-based average nucleotide identity, average amino acid identity, correlation indices of tetra-nucleotide signatures (Tetra) and genome-to-genome distance (digital DNA-DNA hybridization), support the contention that strain F5T represents a novel species of the genus Acidihalobacter. It is proposed that strain F5T should be formally reclassified as Acidihalobacter yilgarnenesis F5T (=DSM 105917T=JCM 32255T).

Identification and Unusual Properties of the Master Regulator FNR in the Extreme Acidophile Acidithiobacillus ferrooxidans

The ability to conserve energy in the presence or absence of oxygen provides a metabolic versatility that confers an advantage in natural ecosystems. The switch between alternative electron transport systems is controlled by the fumarate nitrate reduction transcription factor (FNR) that senses oxygen via an oxygen-sensitive [4Fe-4S]²⁺ iron-sulfur cluster. Under O₂ limiting conditions, FNR plays a key role in allowing bacteria to transition from aerobic to anaerobic lifestyles. This is thought to occur via transcriptional activation of genes involved in anaerobic respiratory pathways and by repression of genes involved in aerobic energy production. The Proteobacterium Acidithiobacillus ferrooxidans is a model species for extremely acidophilic microorganisms that are capable of aerobic and anaerobic growth on elemental sulfur coupled to oxygen and ferric iron reduction, respectively. In this study, an FNR-like protein (FNR_AF) was discovered in At. ferrooxidans that exhibits a primary amino acid sequence and major motifs and domains characteristic of the FNR family of proteins, including an effector binding domain with at least three of the four cysteines known to coordinate an [4Fe-4S]²⁺ center, a dimerization domain, and a DNA binding domain. Western blotting with antibodies against Escherichia coli FNR (FNR_EC) recognized FNR_AF. FNR_AF was able to drive expression from the FNR-responsive E. coli promoter PnarG, suggesting that it is functionally active as an FNR-like protein. Upon air exposure, FNR_AF demonstrated an unusual lack of sensitivity to oxygen compared to the archetypal FNR_EC. Comparison of the primary amino acid sequence of FNR_AF with that of other natural and mutated FNRs, including FNR_EC, coupled with an analysis of the predicted tertiary structure of FNR_AF using the crystal structure of the related FNR from Aliivibrio fisheri as a template revealed a number of amino acid changes that could potentially stabilize FNR_AF in the presence of oxygen. These include a truncated N terminus and amino acid changes both around the putative Fe-S cluster coordinating cysteines and also in the dimer interface. Increased O₂ stability could allow At. ferrooxidans to survive in environments with fluctuating O₂ concentrations, providing an evolutionary advantage in natural, and engineered environments where oxygen gradients shape the bacterial community.

Effect of CO2 Concentration on Uptake and Assimilation of Inorganic Carbon in the Extreme Acidophile Acidithiobacillus ferrooxidans

This study was motivated by surprising gaps in the current knowledge of microbial inorganic carbon (Ci) uptake and assimilation at acidic pH values (pH < 3). Particularly striking is the limited understanding of the differences between Ci uptake mechanisms in acidic versus circumneutral environments where the Ci predominantly occurs either as a dissolved gas (CO2) or as bicarbonate (HCO3 -), respectively. In order to gain initial traction on the problem, the relative abundance of transcripts encoding proteins involved in Ci uptake and assimilation was studied in the autotrophic, polyextreme acidophile Acidithiobacillus ferrooxidans whose optimum pH for growth is 2.5 using ferrous iron as an energy source, although they are able to grow at pH 5 when using sulfur as an energy source. The relative abundance of transcripts of five operons (cbb1-5) and one gene cluster (can-sulP) was monitored by RT-qPCR and, in selected cases, at the protein level by Western blotting, when cells were grown under different regimens of CO2 concentration in elemental sulfur. Of particular note was the absence of a classical bicarbonate uptake system in A. ferrooxidans. However, bioinformatic approaches predict that sulP, previously annotated as a sulfate transporter, is a novel type of bicarbonate transporter. A conceptual model of CO2 fixation was constructed from combined bioinformatic and experimental approaches that suggests strategies for providing ecological flexibility under changing concentrations of CO2 and provides a portal to elucidating Ci uptake and regulation in acidic conditions. The results could advance the understanding of industrial bioleaching processes to recover metals such as copper at acidic pH. In addition, they may also shed light on how chemolithoautotrophic acidophiles influence the nutrient and energy balance in naturally occurring low pH environments.

Uncovering the Mechanisms of Halotolerance in the Extremely Acidophilic Members of the Acidihalobacter Genus Through Comparative Genome Analysis

There are few naturally occurring environments where both acid and salinity stress exist together, consequently, there has been little evolutionary pressure for microorganisms to develop systems that enable them to deal with both stresses simultaneously. Members of the genus Acidihalobacter are iron- and sulfur-oxidizing, halotolerant acidophiles that have developed the ability to tolerate acid and saline stress and, therefore, have the potential to bioleach ores with brackish or saline process waters under acidic conditions. The genus consists of four members, A. prosperus DSM 5130T, A. prosperus DSM 14174, A. prosperus F5 and "A. ferrooxidans" DSM 14175. An in depth genome comparison was undertaken in order to provide a more comprehensive description of the mechanisms of halotolerance used by the different members of this genus. Pangenome analysis identified 29, 3 and 9 protein families related to halotolerance in the core, dispensable and unique genomes, respectively. The genes for halotolerance showed Ka/Ks ratios between 0 and 0.2, confirming that they are conserved and stabilized. All the Acidihalobacter genomes contained similar genes for the synthesis and transport of ectoine, which was recently found to be the dominant osmoprotectant in A. prosperus DSM 14174 and A. prosperus DSM 5130T. Similarities also existed in genes encoding low affinity potassium pumps, however, A. prosperus DSM 14174 was also found to contain genes encoding high affinity potassium pumps. Furthermore, only A. prosperus DSM 5130T and "A. ferrooxidans" DSM 14175 contained genes allowing the uptake of taurine as an osmoprotectant. Variations were also seen in genes encoding proteins involved in the synthesis and/or transport of periplasmic glucans, sucrose, proline, and glycine betaine. This suggests that versatility exists in the Acidihalobacter genus in terms of the mechanisms they can use for halotolerance. This information is useful for developing hypotheses for the search for life on exoplanets and moons.

Mutantelec: An In Silico mutation simulation platform for comparative electrostatic potential profiling of proteins

The electrostatic potential plays a key role in many biological processes like determining the affinity of a ligand to a given protein target, and they are responsible for the catalytic activity of many enzymes. Understanding the effect that amino acid mutations will have on the electrostatic potential of a protein, will allow a thorough understanding of which residues are the most important in a protein. MutantElec, is a friendly web application for in silico generation of site-directed mutagenesis of proteins and the comparison of electrostatic potential between the wild type protein and the mutant(s), based on the three-dimensional structure of the protein. The effect of the mutation is evaluated using different approach to the traditional surface map. MutantElec provides a graphical display of the results that allows the visualization of changes occurring at close distance from the mutation and thus uncovers the local and global impact of a specific change. © 2017 Wiley Periodicals, Inc.

Improved ontology for eukaryotic single-exon coding sequences in biological databases

Efficient extraction of knowledge from biological data requires the development of structured vocabularies to unambiguously define biological terms. This paper proposes descriptions and definitions to disambiguate the term 'single-exon gene'. Eukaryotic Single-Exon Genes (SEGs) have been defined as genes that do not have introns in their protein coding sequences. They have been studied not only to determine their origin and evolution but also because their expression has been linked to several types of human cancer and neurological/developmental disorders and many exhibit tissue-specific transcription. Unfortunately, the term 'SEGs' is rife with ambiguity, leading to biological misinterpretations. In the classic definition, no distinction is made between SEGs that harbor introns in their untranslated regions (UTRs) versus those without. This distinction is important to make because the presence of introns in UTRs affects transcriptional regulation and post-transcriptional processing of the mRNA. In addition, recent whole-transcriptome shotgun sequencing has led to the discovery of many examples of single-exon mRNAs that arise from alternative splicing of multi-exon genes, these single-exon isoforms are being confused with SEGs despite their clearly different origin. The increasing expansion of RNA-seq datasets makes it imperative to distinguish the different SEG types before annotation errors become indelibly propagated in biological databases. This paper develops a structured vocabulary for their disambiguation, allowing a major reassessment of their evolutionary trajectories, regulation, RNA processing and transport, and provides the opportunity to improve the detection of gene associations with disorders including cancers, neurological and developmental diseases.

Draft genome sequence of Acidithiobacillus thiooxidans CLST isolated from the acidic hypersaline Gorbea salt flat in northern Chile

10.1601/nm.2199 CLST is an extremely acidophilic gamma-proteobacteria that was isolated from the Gorbea salt flat, an acidic hypersaline environment in northern Chile. This kind of environment is considered a terrestrial analog of ancient Martian terrains and a source of new material for biotechnological applications. 10.1601/nm.2199 plays a key role in industrial bioleaching; it has the capacity of generating and maintaining acidic conditions by producing sulfuric acid and it can also remove sulfur layers from the surface of minerals, which are detrimental for their dissolution. CLST is a strain of 10.1601/nm.2199 able to tolerate moderate chloride concentrations (up to 15 g L⁻¹ Cl⁻), a feature that is quite unusual in extreme acidophilic microorganisms. Basic microbiological features and genomic properties of this biotechnologically relevant strain are described in this work. The 3,974,949 bp draft genome is arranged into 40 scaffolds of 389 contigs containing 3866 protein-coding genes and 75 RNAs encoding genes. This is the first draft genome of a halotolerant 10.1601/nm.2199 strain. The release of the genome sequence of this strain improves representation of these extreme acidophilic Gram negative bacteria in public databases and strengthens the framework for further investigation of the physiological diversity and ecological function of 10.1601/nm.2199 populations.

Draft genome sequence of the type strain of the sulfur-oxidizing acidophile, Acidithiobacillus albertensis (DSM 14366)

Acidithiobacillus albertensis is an extremely acidophilic, mesophilic, obligatory autotrophic sulfur-oxidizer, with potential importance in the bioleaching of sulfidic metal ores, first described in the 1980s. Here we present the draft genome sequence of Acidithiobacillus albertensis DSM 14366T, thereby both filling a long-standing gap in the genomics of the acidithiobacilli, and providing further insight into the understanding of the biology of the non iron-oxidizing members of the Acidithiobacillus genus. The assembled genome is 3,1 Mb, and contains 47 tRNAs, tmRNA gene and 2 rRNA operons, along with 3149 protein-coding predicted genes. The Whole Genome Shotgun project was deposited in DDBJ/EMBL/GenBank under the accession MOAD00000000.

Multiple Osmotic Stress Responses in Acidihalobacter prosperus Result in Tolerance to Chloride Ions

Extremely acidophilic microorganisms (pH optima for growth of ≤3) are utilized for the extraction of metals from sulfide minerals in the industrial biotechnology of “biomining.” A long term goal for biomining has been development of microbial consortia able to withstand increased chloride concentrations for use in regions where freshwater is scarce. However, when challenged by elevated salt, acidophiles experience both osmotic stress and an acidification of the cytoplasm due to a collapse of the inside positive membrane potential, leading to an influx of protons. In this study, we tested the ability of the halotolerant acidophile Acidihalobacter prosperus to grow and catalyze sulfide mineral dissolution in elevated concentrations of salt and identified chloride tolerance mechanisms in Ac. prosperus as well as the chloride susceptible species, Acidithiobacillus ferrooxidans. Ac. prosperus had optimum iron oxidation at 20 g L⁻¹ NaCl while At. ferrooxidans iron oxidation was inhibited in the presence of 6 g L⁻¹ NaCl. The tolerance to chloride in Ac. prosperus was consistent with electron microscopy, determination of cell viability, and bioleaching capability. The Ac. prosperus proteomic response to elevated chloride concentrations included the production of osmotic stress regulators that potentially induced production of the compatible solute, ectoine uptake protein, and increased iron oxidation resulting in heightened electron flow to drive proton export by the F₀F₁ ATPase. In contrast, At. ferrooxidans responded to low levels of Cl⁻ with a generalized stress response, decreased iron oxidation, and an increase in central carbon metabolism. One potential adaptation to high chloride in the Ac. prosperus Rus protein involved in ferrous iron oxidation was an increase in the negativity of the surface potential of Rus Form I (and Form II) that could help explain how it can be active under elevated chloride concentrations. These data have been used to create a model of chloride tolerance in the salt tolerant and susceptible species Ac. prosperus and At. ferrooxidans, respectively.

Bioinformatic Analyses of Unique (Orphan) Core Genes of the Genus Acidithiobacillus: Functional Inferences and Use As Molecular Probes for Genomic and Metagenomic/Transcriptomic Interrogation

Using phylogenomic and gene compositional analyses, five highly conserved gene families have been detected in the core genome of the phylogenetically coherent genus Acidithiobacillus of the class Acidithiobacillia. These core gene families are absent in the closest extant genus Thermithiobacillus tepidarius that subtends the Acidithiobacillus genus and roots the deepest in this class. The predicted proteins encoded by these core gene families are not detected by a BLAST search in the NCBI non-redundant database of more than 90 million proteins using a relaxed cut-off of 1.0e⁻⁵. None of the five families has a clear functional prediction. However, bioinformatic scrutiny, using pI prediction, motif/domain searches, cellular location predictions, genomic context analyses, and chromosome topology studies together with previously published transcriptomic and proteomic data, suggests that some may have functions associated with membrane remodeling during cell division perhaps in response to pH stress. Despite the high level of amino acid sequence conservation within each family, there is sufficient nucleotide variation of the respective genes to permit the use of the DNA sequences to distinguish different species of Acidithiobacillus, making them useful additions to the armamentarium of tools for phylogenetic analysis. Since the protein families are unique to the Acidithiobacillus genus, they can also be leveraged as probes to detect the genus in environmental metagenomes and metatranscriptomes, including industrial biomining operations, and acid mine drainage (AMD).

Aerobic Lineage of the Oxidative Stress Response Protein Rubrerythrin Emerged in an Ancient Microaerobic, (Hyper)Thermophilic Environment

Rubrerythrins (RBRs) are non-heme di-iron proteins belonging to the ferritin-like superfamily. They are involved in oxidative stress defense as peroxide scavengers in a wide range of organisms. The vast majority of RBRs, including classical forms of this protein, contain a C-terminal rubredoxin-like domain involved in electron transport that is used during catalysis in anaerobic conditions. Rubredoxin is an ancient and large protein family of short length (<100 residues) that contains a Fe-S center involved in electron transfer. However, functional forms of the enzyme lacking the rubredoxin-like domain have been reported (e.g., sulerythrin and ferriperoxin). In this study, phylogenomic evidence is presented that suggests that a complete lineage of rubrerythrins, lacking the rubredoxin-like domain, arose in an ancient microaerobic and (hyper)thermophilic environments in the ancestors of the Archaea Thermoproteales and Sulfolobales. This lineage (termed the “aerobic-type” lineage) subsequently evolved to become adapted to environments with progressively lower temperatures and higher oxygen concentrations via the acquisition of two co-localized genes, termed DUF3501 and RFO, encoding a conserved protein of unknown function and a predicted Fe-S oxidoreductase, respectively. Proposed Horizontal Gene Transfer events from these archaeal ancestors to Bacteria expanded the opportunities for further evolution of this RBR including adaption to lower temperatures. The second lineage (termed the cyanobacterial lineage) is proposed to have evolved in cyanobacterial ancestors, maybe in direct response to the production of oxygen via oxygenic photosynthesis during the Great Oxygen Event (GOE). It is hypothesized that both lineages of RBR emerged in a largely anaerobic world with “whiffs” of oxygen and that their subsequent independent evolutionary trajectories allowed microorganisms to transition from this anaerobic world to an aerobic one.

Graphlet Based Metrics for the Comparison of Gene Regulatory Networks

Understanding the control of gene expression remains one of the main challenges in the post-genomic era. Accordingly, a plethora of methods exists to identify variations in gene expression levels. These variations underlay almost all relevant biological phenomena, including disease and adaptation to environmental conditions. However, computational tools to identify how regulation changes are scarce. Regulation of gene expression is usually depicted in the form of a gene regulatory network (GRN). Structural changes in a GRN over time and conditions represent variations in the regulation of gene expression. Like other biological networks, GRNs are composed of basic building blocks called graphlets. As a consequence, two new metrics based on graphlets are proposed in this work: REConstruction Rate (REC) and REC Graphlet Degree (RGD). REC determines the rate of graphlet similarity between different states of a network and RGD identifies the subset of nodes with the highest topological variation. In other words, RGD discerns how th GRN was rewired. REC and RGD were used to compare the local structure of nodes in condition-specific GRNs obtained from gene expression data of Escherichia coli, forming biofilms and cultured in suspension. According to our results, most of the network local structure remains unaltered in the two compared conditions. Nevertheless, changes reported by RGD necessarily imply that a different cohort of regulators (i.e. transcription factors (TFs)) appear on the scene, shedding light on how the regulation of gene expression occurs when E. coli transits from suspension to biofilm. Consequently, we propose that both metrics REC and RGD should be adopted as a quantitative approach to conduct differential analyses of GRNs. A tool that implements both metrics is available as an on-line web server (http://dlab.cl/loto).

Genomic and metagenomic challenges and opportunities for bioleaching: a mini-review

High-throughput genomic technologies are accelerating progress in understanding the diversity of microbial life in many environments. Here we highlight advances in genomics and metagenomics of microorganisms from bioleaching heaps and related acidic mining environments. Bioleaching heaps used for copper recovery provide significant opportunities to study the processes and mechanisms underlying microbial successions and the influence of community composition on ecosystem functioning. Obtaining quantitative and process-level knowledge of these dynamics is pivotal for understanding how microorganisms contribute to the solubilization of copper for industrial recovery. Advances in DNA sequencing technology provide unprecedented opportunities to obtain information about the genomes of bioleaching microorganisms, allowing predictive models of metabolic potential and ecosystem-level interactions to be constructed. These approaches are enabling predictive phenotyping of organisms many of which are recalcitrant to genetic approaches or are unculturable. This mini-review describes current bioleaching genomic and metagenomic projects and addresses the use of genome information to: (i) build metabolic models; (ii) predict microbial interactions; (iii) estimate genetic diversity; and (iv) study microbial evolution. Key challenges and perspectives of bioleaching genomics/metagenomics are addressed.

SinEx DB: a database for single exon coding sequences in mammalian genomes

Eukaryotic genes are typically interrupted by intragenic, noncoding sequences termed introns. However, some genes lack introns in their coding sequence (CDS) and are generally known as 'single exon genes' (SEGs). In this work, a SEG is defined as a nuclear, protein-coding gene that lacks introns in its CDS. Whereas, many public databases of Eukaryotic multi-exon genes are available, there are only two specialized databases for SEGs. The present work addresses the need for a more extensive and diverse database by creating SinEx DB, a publicly available, searchable database of predicted SEGs from 10 completely sequenced mammalian genomes including human. SinEx DB houses the DNA and protein sequence information of these SEGs and includes their functional predictions (KOG) and the relative distribution of these functions within species. The information is stored in a relational database built with My SQL Server 5.1.33 and the complete dataset of SEG sequences and their functional predictions are available for downloading. SinEx DB can be interrogated by: (i) a browsable phylogenetic schema, (ii) carrying out BLAST searches to the in-house SinEx DB of SEGs and (iii) via an advanced search mode in which the database can be searched by key words and any combination of searches by species and predicted functions. SinEx DB provides a rich source of information for advancing our understanding of the evolution and function of SEGs.Database URL: www.sinex.cl.

Gene Loss and Horizontal Gene Transfer Contributed to the Genome Evolution of the Extreme Acidophile "Ferrovum"

Acid mine drainage (AMD), associated with active and abandoned mining sites, is a habitat for acidophilic microorganisms that gain energy from the oxidation of reduced sulfur compounds and ferrous iron and that thrive at pH below 4. Members of the recently proposed genus “Ferrovum” are the first acidophilic iron oxidizers to be described within the Betaproteobacteria. Although they have been detected as typical community members in AMD habitats worldwide, knowledge of their phylogenetic and metabolic diversity is scarce. Genomics approaches appear to be most promising in addressing this lacuna since isolation and cultivation of “Ferrovum” has proven to be extremely difficult and has so far only been successful for the designated type strain “Ferrovum myxofaciens” P3G. In this study, the genomes of two novel strains of “Ferrovum” (PN-J185 and Z-31) derived from water samples of a mine water treatment plant were sequenced. These genomes were compared with those of “Ferrovum” sp. JA12 that also originated from the mine water treatment plant, and of the type strain (P3G). Phylogenomic scrutiny suggests that the four strains represent three “Ferrovum” species that cluster in two groups (1 and 2). Comprehensive analysis of their predicted metabolic pathways revealed that these groups harbor characteristic metabolic profiles, notably with respect to motility, chemotaxis, nitrogen metabolism, biofilm formation and their potential strategies to cope with the acidic environment. For example, while the “F. myxofaciens” strains (group 1) appear to be motile and diazotrophic, the non-motile group 2 strains have the predicted potential to use a greater variety of fixed nitrogen sources. Furthermore, analysis of their genome synteny provides first insights into their genome evolution, suggesting that horizontal gene transfer and genome reduction in the group 2 strains by loss of genes encoding complete metabolic pathways or physiological features contributed to the observed diversification.

Draft genome sequence of chloride-tolerant Leptospirillum ferriphilum Sp-Cl from industrial bioleaching operations in northern Chile

Leptospirillum ferriphilum Sp-Cl is a Gram negative, thermotolerant, curved, rod-shaped bacterium, isolated from an industrial bioleaching operation in northern Chile, where chalcocite is the major copper mineral and copper hydroxychloride atacamite is present in variable proportions in the ore. This strain has unique features as compared to the other members of the species, namely resistance to elevated concentrations of chloride, sulfate and metals. Basic microbiological features and genomic properties of this biotechnologically relevant strain are described in this work. The 2,475,669 bp draft genome is arranged into 74 scaffolds of 74 contigs. A total of 48 RNA genes and 2,834 protein coding genes were predicted from its annotation; 55 % of these were assigned a putative function. Release of the genome sequence of this strain will provide further understanding of the mechanisms used by acidophilic bacteria to endure high osmotic stress and high chloride levels and of the role of chloride-tolerant iron-oxidizers in industrial bioleaching operations.

Genome Analysis of the Biotechnologically Relevant Acidophilic Iron Oxidising Strain JA12 Indicates Phylogenetic and Metabolic Diversity within the Novel Genus "Ferrovum"

Background

Members of the genus “Ferrovum” are ubiquitously distributed in acid mine drainage (AMD) waters which are characterised by their high metal and sulfate loads. So far isolation and microbiological characterisation have only been successful for the designated type strain “Ferrovum myxofaciens” P3G. Thus, knowledge about physiological characteristics and the phylogeny of the genus “Ferrovum” is extremely scarce.

Objective

In order to access the wider genetic pool of the genus “Ferrovum” we sequenced the genome of a “Ferrovum”-containing mixed culture and successfully assembled the almost complete genome sequence of the novel “Ferrovum” strain JA12.

Phylogeny and Lifestyle

The genome-based phylogenetic analysis indicates that strain JA12 and the type strain represent two distinct “Ferrovum” species. “Ferrovum” strain JA12 is characterised by an unusually small genome in comparison to the type strain and other iron oxidising bacteria. The prediction of nutrient assimilation pathways suggests that “Ferrovum” strain JA12 maintains a chemolithoautotrophic lifestyle utilising carbon dioxide and bicarbonate, ammonium and urea, sulfate, phosphate and ferrous iron as carbon, nitrogen, sulfur, phosphorous and energy sources, respectively.

Unique Metabolic Features

The potential utilisation of urea by “Ferrovum” strain JA12 is moreover remarkable since it may furthermore represent a strategy among extreme acidophiles to cope with the acidic environment. Unlike other acidophilic chemolithoautotrophs “Ferrovum” strain JA12 exhibits a complete tricarboxylic acid cycle, a metabolic feature shared with the closer related neutrophilic iron oxidisers among the Betaproteobacteria including Sideroxydans lithotrophicus and Thiobacillus denitrificans. Furthermore, the absence of characteristic redox proteins involved in iron oxidation in the well-studied acidophiles Acidithiobacillus ferrooxidans (rusticyanin) and Acidithiobacillus ferrivorans (iron oxidase) indicates the existence of a modified pathway in “Ferrovum” strain JA12. Therefore, the results of the present study extend our understanding of the genus “Ferrovum” and provide a comprehensive framework for future comparative genome and metagenome studies.

Reclassification of 'Thiobacillus prosperus' Huber and Stetter 1989 as Acidihalobacter prosperus gen. nov., sp. nov., a member of the family Ectothiorhodospiraceae

Analysis of phylogenomic metrics of a recently released draft genome sequence of the halotolerant, acidophile 'Thiobacillus prosperus' DSM 5130 indicates that it is not a member of the genus Thiobacillus within the class Betaproteobacteria as originally proposed. Based on data from 16S rRNA gene phylogeny, and analyses of multiprotein phylogeny and average nucleotide identity (ANI), we show that it belongs to a new genus within the family Ectothiorhodospiraceae, for which we propose the name Acidihalobacter gen. nov. In accordance, it is proposed that 'Thiobacillus prosperus' DSM 5130 be named Acidihalobacter prosperus gen. nov., sp. nov. DSM 5130T ( = JCM 30709T) and that it becomes the type strain of the type species of this genus.

Metagenomic analysis reveals adaptations to a cold-adapted lifestyle in a low-temperature acid mine drainage stream

An acid mine drainage (pH 2.5-2.7) stream biofilm situated 250 m below ground in the low-temperature (6-10°C) Kristineberg mine, northern Sweden, contained a microbial community equipped for growth at low temperature and acidic pH. Metagenomic sequencing of the biofilm and planktonic fractions identified the most abundant microorganism to be similar to the psychrotolerant acidophile, Acidithiobacillus ferrivorans. In addition, metagenome contigs were most similar to other Acidithiobacillus species, an Acidobacteria-like species, and a Gallionellaceae-like species. Analyses of the metagenomes indicated functional characteristics previously characterized as related to growth at low temperature including cold-shock proteins, several pathways for the production of compatible solutes and an anti-freeze protein. In addition, genes were predicted to encode functions related to pH homeostasis and metal resistance related to growth in the acidic metal-containing mine water. Metagenome analyses identified microorganisms capable of nitrogen fixation and exhibiting a primarily autotrophic lifestyle driven by the oxidation of the ferrous iron and inorganic sulfur compounds contained in the sulfidic mine waters. The study identified a low diversity of abundant microorganisms adapted to a low-temperature acidic environment as well as identifying some of the strategies the microorganisms employ to grow in this extreme environment.

Multi Locus Sequence Typing scheme for Acidithiobacillus caldus strain evaluation and differentiation

Phenotypic, metabolic and genetic properties of several Acidithiobacillus caldus strains indicate the existence of as yet undefined levels of variation within the species. Inspite of this, intraspecies genetic diversity has not yet been explored in detail. In this study, the design and implementation of a Multi Locus Sequence Typing (MLST) scheme for At. caldus is described. This represents the first MLST-based study applied to industrial isolates of the species. Seven informative and discriminant MLST markers were selected using a sequence-driven approach and a custom-designed bioinformatic pipeline. The allelic profiles of thirteen At. caldus strains from diverse geographical origins and industrial settings were derived using this scheme. MLST-based population structure analysis indicated only moderate amounts of genetic diversity within the set of strains, further supporting their current assignment to a single species. Also, no clear evidence for geographical isolation could be derived from this study. However, the prevalence of sequence type 1 in heap leaching industrial settings support the view that bioprocess conditions and dynamics may have a strong influence on At. caldus (microbial) microdiversity patterns. The MLST scheme presented herein is a valuable tool for the identification and classification of strains of At. caldus for either ecological or evolutionary studies and possibly also for industrial monitoring purposes.

Genetic variability of psychrotolerant Acidithiobacillus ferrivorans revealed by (meta)genomic analysis

Acidophilic microorganisms inhabit low pH environments such as acid mine drainage that is generated when sulfide minerals are exposed to air. The genome sequence of the psychrotolerant Acidithiobacillus ferrivorans SS3 was compared to a metagenome from a low temperature acidic stream dominated by an A. ferrivorans-like strain. Stretches of genomic DNA characterized by few matches to the metagenome, termed 'metagenomic islands', encoded genes associated with metal efflux and pH homeostasis. The metagenomic islands were enriched in mobile elements such as phage proteins, transposases, integrases and in one case, predicted to be flanked by truncated tRNAs. Cus gene clusters predicted to be involved in copper efflux and further Cus-like RND systems were predicted to be located in metagenomic islands and therefore, constitute part of the flexible gene complement of the species. Phylogenetic analysis of Cus clusters showed both lineage specificity within the Acidithiobacillus genus as well as niche specificity associated with an acidic environment. The metagenomic islands also contained a predicted copper efflux P-type ATPase system and a polyphosphate kinase potentially involved in polyphosphate mediated copper resistance. This study identifies genetic variability of low temperature acidophiles that likely reflects metal resistance selective pressures in the copper rich environment.

Metal resistance or tolerance? Acidophiles confront high metal loads via both abiotic and biotic mechanisms

All metals are toxic at high concentrations and consequently their intracellular concentrations must be regulated. Extremely acidophilic microorganisms have an optimum growth of pH <3 and proliferate in natural and anthropogenic low pH environments. Some acidophiles are involved in the catalysis of sulfide mineral dissolution, resulting in high concentrations of metals in solution. Acidophiles are often described as highly metal resistant via mechanisms such as multiple and/or more efficient active resistance systems than are present in neutrophiles. However, this is not the case for all acidophiles and we contend that their growth in high metal concentrations is partially due to an intrinsic tolerance as a consequence of the environment in which they live. In this perspective, we highlight metal tolerance via complexation of free metals by sulfate ions and passive tolerance to metal influx via an internal positive cytoplasmic transmembrane potential. These tolerance mechanisms have been largely ignored in past studies of acidophile growth in the presence of metals and should be taken into account.

Architecture and gene repertoire of the flexible genome of the extreme acidophile Acidithiobacillus caldus

BACKGROUND

Acidithiobacillus caldus is a sulfur oxidizing extreme acidophile and the only known mesothermophile within the Acidithiobacillales. As such, it is one of the preferred microbes for mineral bioprocessing at moderately high temperatures. In this study, we explore the genomic diversity of A. caldus strains using a combination of bioinformatic and experimental techniques, thus contributing first insights into the elucidation of the species pangenome.

PRINCIPAL FINDINGS

Comparative sequence analysis of A. caldus ATCC 51756 and SM-1 indicate that, despite sharing a conserved and highly syntenic genomic core, both strains have unique gene complements encompassing nearly 20% of their respective genomes. The differential gene complement of each strain is distributed between the chromosomal compartment, one megaplasmid and a variable number of smaller plasmids, and is directly associated to a diverse pool of mobile genetic elements (MGE). These include integrative conjugative and mobilizable elements, genomic islands and insertion sequences. Some of the accessory functions associated to these MGEs have been linked previously to the flexible gene pool in microorganisms inhabiting completely different econiches. Yet, others had not been unambiguously mapped to the flexible gene pool prior to this report and clearly reflect strain-specific adaption to local environmental conditions.

SIGNIFICANCE

For many years, and because of DNA instability at low pH and recurrent failure to genetically transform acidophilic bacteria, gene transfer in acidic environments was considered negligible. Findings presented herein imply that a more or less conserved pool of actively excising MGEs occurs in the A. caldus population and point to a greater frequency of gene exchange in this econiche than previously recognized. Also, the data suggest that these elements endow the species with capacities to withstand the diverse abiotic and biotic stresses of natural environments, in particular those associated with its extreme econiche.

Molecular dynamics study of the archaeal aquaporin AqpM

Background

Aquaporins are a large family of transmembrane channel proteins that are present throughout all domains of life and are implicated in human disorders. These channels, allow the passive but selective movement of water and other small neutral solutes across cell membranes. Aquaporins have been classified into two sub-families: i) strict aquaporins that only allow the passage of water and ii) the less selective aquaglyceroporins that transport water and other neutral solutes, such as glycerol, CO₂ or urea. Recently, the identification and characterization of a number of archaeal and bacterial aquaporins suggested the existence of a third sub-family; one that is neither a strict aquaporin nor an aquaglyceroporin. The function and phylogeny of this third family is still a matter of debate.

Results

Twenty nanosecond molecular dynamics (MD) simulation of a fully hydrated tetramer of AqpM embedded in a lipid bilayer permitted predictions to be made of key biophysical parameters including: single channel osmotic permeability constant (p_f), single channel diffusive permeability constant (p_d), channel radius, potential water occupancy of the channel and water orientation inside the pore. These properties were compared with those of well characterized representatives of the two main aquaporin sub-families. Results show that changes in the amino acid composition of the aromatic/arginine region affect the size and polarity of the selectivity filter (SF) and could help explain the difference in water permeability between aquaporins. In addition, MD simulation results suggest that AqpM combines characteristics of strict aquaporins, such as the narrow SF and channel radius, with those of aquaglyceroporins, such as a more hydrophobic and less polar SF.

Conclusions

MD simulations of AqpM extend previous evidence that this archaeal aquaporin exhibits hybrid features intermediate between the two known aquaporin sub-families, supporting the idea that it may constitute a member of a novel class of aquaporins.

Draft genome sequence of the extremely acidophilic biomining bacterium Acidithiobacillus thiooxidans ATCC 19377 provides insights into the evolution of the Acidithiobacillus genus

Acidithiobacillus thiooxidans is a mesophilic, extremely acidophilic, chemolithoautotrophic gammaproteobacterium that derives energy from the oxidation of sulfur and inorganic sulfur compounds. Here we present the draft genome sequence of A. thiooxidans ATCC 19377, which has allowed the identification of genes for survival and colonization of extremely acidic environments.

Sulfur metabolism in the extreme acidophile acidithiobacillus caldus

Given the challenges to life at low pH, an analysis of inorganic sulfur compound (ISC) oxidation was initiated in the chemolithoautotrophic extremophile Acidithiobacillus caldus. A. caldus is able to metabolize elemental sulfur and a broad range of ISCs. It has been implicated in the production of environmentally damaging acidic solutions as well as participating in industrial bioleaching operations where it forms part of microbial consortia used for the recovery of metal ions. Based upon the recently published A. caldus type strain genome sequence, a bioinformatic reconstruction of elemental sulfur and ISC metabolism predicted genes included: sulfide-quinone reductase (sqr), tetrathionate hydrolase (tth), two sox gene clusters potentially involved in thiosulfate oxidation (soxABXYZ), sulfur oxygenase reductase (sor), and various electron transport components. RNA transcript profiles by semi quantitative reverse transcription PCR suggested up-regulation of sox genes in the presence of tetrathionate. Extensive gel based proteomic comparisons of total soluble and membrane enriched protein fractions during growth on elemental sulfur and tetrathionate identified differential protein levels from the two Sox clusters as well as several chaperone and stress proteins up-regulated in the presence of elemental sulfur. Proteomics results also suggested the involvement of heterodisulfide reductase (HdrABC) in A. caldus ISC metabolism. A putative new function of Hdr in acidophiles is discussed. Additional proteomic analysis evaluated protein expression differences between cells grown attached to solid, elemental sulfur versus planktonic cells. This study has provided insights into sulfur metabolism of this acidophilic chemolithotroph and gene expression during attachment to solid elemental sulfur.

Genes and pathways for CO2 fixation in the obligate, chemolithoautotrophic acidophile, Acidithiobacillus ferrooxidans, carbon fixation in A. ferrooxidans

Background

Acidithiobacillus ferrooxidans is chemolithoautotrophic γ-proteobacterium that thrives at extremely low pH (pH 1-2). Although a substantial amount of information is available regarding CO₂ uptake and fixation in a variety of facultative autotrophs, less is known about the processes in obligate autotrophs, especially those living in extremely acidic conditions, prompting the present study.

Results

Four gene clusters (termed cbb1-4) in the A. ferrooxidans genome are predicted to encode enzymes and structural proteins involved in carbon assimilation via the Calvin-Benson-Bassham (CBB) cycle including form I of ribulose-1,5-bisphosphate carboxylase/oxygenase (RubisCO, EC 4.1.1.39) and the CO₂-concentrating carboxysomes. RT-PCR experiments demonstrated that each gene cluster is a single transcriptional unit and thus is an operon. Operon cbb1 is divergently transcribed from a gene, cbbR, encoding the LysR-type transcriptional regulator CbbR that has been shown in many organisms to regulate the expression of RubisCO genes. Sigma⁷⁰-like -10 and -35 promoter boxes and potential CbbR-binding sites (T-N₁₁-A/TNA-N₇TNA) were predicted in the upstream regions of the four operons. Electrophoretic mobility shift assays (EMSAs) confirmed that purified CbbR is able to bind to the upstream regions of the cbb1, cbb2 and cbb3 operons, demonstrating that the predicted CbbR-binding sites are functional in vitro. However, CbbR failed to bind the upstream region of the cbb4 operon that contains cbbP, encoding phosphoribulokinase (EC 2.7.1.19). Thus, other factors not present in the assay may be required for binding or the region lacks a functional CbbR-binding site. The cbb3 operon contains genes predicted to encode anthranilate synthase components I and II, catalyzing the formation of anthranilate and pyruvate from chorismate. This suggests a novel regulatory connection between CO₂ fixation and tryptophan biosynthesis. The presence of a form II RubisCO could promote the ability of A. ferrooxidans to fix CO₂ at different concentrations of CO₂.

Conclusions

A. ferrooxidans has features of cbb gene organization for CO₂-assimilating functions that are characteristic of obligate chemolithoautotrophs and distinguish this group from facultative autotrophs. The most conspicuous difference is a separate operon for the cbbP gene. It is hypothesized that this organization may provide greater flexibility in the regulation of expression of genes involved in inorganic carbon assimilation.

Extending the models for iron and sulfur oxidation in the extreme acidophile Acidithiobacillus ferrooxidans

Background

Acidithiobacillus ferrooxidans gains energy from the oxidation of ferrous iron and various reduced inorganic sulfur compounds at very acidic pH. Although an initial model for the electron pathways involved in iron oxidation has been developed, much less is known about the sulfur oxidation in this microorganism. In addition, what has been reported for both iron and sulfur oxidation has been derived from different A. ferrooxidans strains, some of which have not been phylogenetically characterized and some have been shown to be mixed cultures. It is necessary to provide models of iron and sulfur oxidation pathways within one strain of A. ferrooxidans in order to comprehend the full metabolic potential of the pangenome of the genus.

Results

Bioinformatic-based metabolic reconstruction supported by microarray transcript profiling and quantitative RT-PCR analysis predicts the involvement of a number of novel genes involved in iron and sulfur oxidation in A. ferrooxidans ATCC23270. These include for iron oxidation: cup (copper oxidase-like), ctaABT (heme biogenesis and insertion), nuoI and nuoK (NADH complex subunits), sdrA1 (a NADH complex accessory protein) and atpB and atpE (ATP synthetase F0 subunits). The following new genes are predicted to be involved in reduced inorganic sulfur compounds oxidation: a gene cluster (rhd, tusA, dsrE, hdrC, hdrB, hdrA, orf2, hdrC, hdrB) encoding three sulfurtransferases and a heterodisulfide reductase complex, sat potentially encoding an ATP sulfurylase and sdrA2 (an accessory NADH complex subunit). Two different regulatory components are predicted to be involved in the regulation of alternate electron transfer pathways: 1) a gene cluster (ctaRUS) that contains a predicted iron responsive regulator of the Rrf2 family that is hypothesized to regulate cytochrome aa₃ oxidase biogenesis and 2) a two component sensor-regulator of the RegB-RegA family that may respond to the redox state of the quinone pool.

Conclusion

Bioinformatic analysis coupled with gene transcript profiling extends our understanding of the iron and reduced inorganic sulfur compounds oxidation pathways in A. ferrooxidans and suggests mechanisms for their regulation. The models provide unified and coherent descriptions of these processes within the type strain, eliminating previous ambiguity caused by models built from analyses of multiple and divergent strains of this microorganism.

Selection and evaluation of reference genes for improved interrogation of microbial transcriptomes: case study with the extremophile Acidithiobacillus ferrooxidans

Background

Normalization is a prerequisite for accurate real time PCR (qPCR) expression analysis and for the validation of microarray profiling data in microbial systems. The choice and use of reference genes that are stably expressed across samples, experimental conditions and designs is a key consideration for the accurate interpretation of gene expression data.

Results

Here, we evaluate a carefully selected set of reference genes derived from previous microarray-based transcriptional profiling experiments performed on Acidithiobacillus ferrooxidans and identify a set of genes with minimal variability under five different experimental conditions that are frequently used in Acidithiobacilli research. Suitability of these and other previously reported reference genes to monitor the expression of four selected target genes from A. ferrooxidans grown with different energy sources was investigated. Utilization of reference genes map, rpoC, alaS and era results in improved interpretation of gene expression profiles in A. ferrooxidans.

Conclusion

This investigation provides a validated set of reference genes for studying A. ferrooxidans gene expression under typical biological conditions and an initial point of departure for exploring new experimental setups in this microorganism and eventually in other closely related Acidithiobacilli. The information could also be of value for future transcriptomic experiments in other bacterial systems.

Acidithiobacillus ferrooxidans metabolism: from genome sequence to industrial applications

BACKGROUND

Acidithiobacillus ferrooxidans is a major participant in consortia of microorganisms used for the industrial recovery of copper (bioleaching or biomining). It is a chemolithoautrophic, gamma-proteobacterium using energy from the oxidation of iron- and sulfur-containing minerals for growth. It thrives at extremely low pH (pH 1-2) and fixes both carbon and nitrogen from the atmosphere. It solubilizes copper and other metals from rocks and plays an important role in nutrient and metal biogeochemical cycling in acid environments. The lack of a well-developed system for genetic manipulation has prevented thorough exploration of its physiology. Also, confusion has been caused by prior metabolic models constructed based upon the examination of multiple, and sometimes distantly related, strains of the microorganism.

RESULTS

The genome of the type strain A. ferrooxidans ATCC 23270 was sequenced and annotated to identify general features and provide a framework for in silico metabolic reconstruction. Earlier models of iron and sulfur oxidation, biofilm formation, quorum sensing, inorganic ion uptake, and amino acid metabolism are confirmed and extended. Initial models are presented for central carbon metabolism, anaerobic metabolism (including sulfur reduction, hydrogen metabolism and nitrogen fixation), stress responses, DNA repair, and metal and toxic compound fluxes.

CONCLUSION

Bioinformatics analysis provides a valuable platform for gene discovery and functional prediction that helps explain the activity of A. ferrooxidans in industrial bioleaching and its role as a primary producer in acidic environments. An analysis of the genome of the type strain provides a coherent view of its gene content and metabolic potential.

Microbial iron management mechanisms in extremely acidic environments: comparative genomics evidence for diversity and versatility

BACKGROUND

Iron is an essential nutrient but can be toxic at high intracellular concentrations and organisms have evolved tightly regulated mechanisms for iron uptake and homeostasis. Information on iron management mechanisms is available for organisms living at circumneutral pH. However, very little is known about how acidophilic bacteria, especially those used for industrial copper bioleaching, cope with environmental iron loads that can be 1018 times the concentration found in pH neutral environments. This study was motivated by the need to fill this lacuna in knowledge. An understanding of how microorganisms thrive in acidic ecosystems with high iron loads requires a comprehensive investigation of the strategies to acquire iron and to coordinate this acquisition with utilization, storage and oxidation of iron through metal responsive regulation. In silico prediction of iron management genes and Fur regulation was carried out for three Acidithiobacilli: Acidithiobacillus ferrooxidans (iron and sulfur oxidizer) A. thiooxidans and A. caldus (sulfur oxidizers) that can live between pH 1 and pH 5 and for three strict iron oxidizers of the Leptospirillum genus that live at pH 1 or below.

RESULTS

Acidithiobacilli have predicted FeoB-like Fe(II) and Nramp-like Fe(II)-Mn(II) transporters. They also have 14 different TonB dependent ferri-siderophore transporters of diverse siderophore affinity, although they do not produce classical siderophores. Instead they have predicted novel mechanisms for dicitrate synthesis and possibly also for phosphate-chelation mediated iron uptake. It is hypothesized that the unexpectedly large number and diversity of Fe(III)-uptake systems confers versatility to this group of acidophiles, especially in higher pH environments (pH 4-5) where soluble iron may not be abundant. In contrast, Leptospirilla have only a FtrI-Fet3P-like permease and three TonB dependent ferri-dicitrate siderophore systems. This paucity of iron uptake systems could reflect their obligatory occupation of extremely low pH environments where high concentrations of soluble iron may always be available and were oxidized sulfur species might not compromise iron speciation dynamics. Presence of bacterioferritin in the Acidithiobacilli, polyphosphate accumulation functions and variants of FieF-like diffusion facilitators in both Acidithiobacilli and Leptospirilla, indicate that they may remove or store iron under conditions of variable availability. In addition, the Fe(II)-oxidizing capacity of both A. ferrooxidans and Leptospirilla could itself be a way to evade iron stress imposed by readily available Fe(II) ions at low pH. Fur regulatory sites have been predicted for a number of gene clusters including iron related and non-iron related functions in both the Acidithiobacilli and Leptospirilla, laying the foundation for the future discovery of iron regulated and iron-phosphate coordinated regulatory control circuits.

CONCLUSION

In silico analyses of the genomes of acidophilic bacteria are beginning to tease apart the mechanisms that mediate iron uptake and homeostasis in low pH environments. Initial models pinpoint significant differences in abundance and diversity of iron management mechanisms between Leptospirilla and Acidithiobacilli, and begin to reveal how these two groups respond to iron cycling and iron fluctuations in naturally acidic environments and in industrial operations. Niche partitions and ecological successions between acidophilic microorganisms may be partially explained by these observed differences. Models derived from these analyses pave the way for improved hypothesis testing and well directed experimental investigation. In addition, aspects of these models should challenge investigators to evaluate alternative iron management strategies in non-acidophilic model organisms.

AlterORF: a database of alternate open reading frames

AlterORF is a searchable database that contains information regarding alternate open reading frames (ORFs) for over 1.5 million genes in 481 prokaryotic genomes. The objective of the database is to provide a platform for improving genome annotation and to serve as an aid for the identification of prokaryotic genes that potentially encode proteins in more than one reading frame. The AlterORF Database can be accessed through a web interface at www.alterorf.cl

Regulation of a novel Acidithiobacillus caldus gene cluster involved in metabolism of reduced inorganic sulfur compounds

Acidithiobacillus caldus has been proposed to play a role in the oxidation of reduced inorganic sulfur compounds (RISCs) produced in industrial biomining of sulfidic minerals. Here, we describe the regulation of a new cluster containing the gene encoding tetrathionate hydrolase (tetH), a key enzyme in the RISC metabolism of this bacterium. The cluster contains five cotranscribed genes, ISac1, rsrR, rsrS, tetH, and doxD, coding for a transposase, a two-component response regulator (RsrR and RsrS), tetrathionate hydrolase, and DoxD, respectively. As shown by quantitative PCR, rsrR, tetH, and doxD are upregulated to different degrees in the presence of tetrathionate. Western blot analysis also indicates upregulation of TetH in the presence of tetrathionate, thiosulfate, and pyrite. The tetH cluster is predicted to have two promoters, both of which are functional in Escherichia coli and one of which was mapped by primer extension. A pyrrolo-quinoline quinone binding domain in TetH was predicted by bioinformatic analysis, and the presence of an o-quinone moiety was experimentally verified, suggesting a mechanism for tetrathionate oxidation.

Bioinformatic prediction and experimental verification of Fur-regulated genes in the extreme acidophile Acidithiobacillus ferrooxidans

The gamma-proteobacterium Acidithiobacillus ferrooxidans lives in extremely acidic conditions (pH 2) and, unlike most organisms, is confronted with an abundant supply of soluble iron. It is also unusual in that it oxidizes iron as an energy source. Consequently, it faces the challenging dual problems of (i) maintaining intracellular iron homeostasis when confronted with extremely high environmental loads of iron and (ii) of regulating the use of iron both as an energy source and as a metabolic micronutrient. A combined bioinformatic and experimental approach was undertaken to identify Fur regulatory sites in the genome of A. ferrooxidans and to gain insight into the constitution of its Fur regulon. Fur regulatory targets associated with a variety of cellular functions including metal trafficking (e.g. feoPABC, tdr, tonBexbBD, copB, cdf), utilization (e.g. fdx, nif), transcriptional regulation (e.g. phoB, irr, iscR) and redox balance (grx, trx, gst) were identified. Selected predicted Fur regulatory sites were confirmed by FURTA, EMSA and in vitro transcription analyses. This study provides the first model for a Fur-binding site consensus sequence in an acidophilic iron-oxidizing microorganism and lays the foundation for future studies aimed at deepening our understanding of the regulatory networks that control iron uptake, homeostasis and oxidation in extreme acidophiles.

Second acyl homoserine lactone production system in the extreme acidophile Acidithiobacillus ferrooxidans

The acidophilic proteobacterium Acidithiobacillus ferrooxidans is involved in the industrial biorecovery of copper. It is found in acidic environments in biofilms and is important in the biogeochemical cycling of metals and nutrients. Its genome contains a cluster of four genes, glyQ, glysS, gph, and act, that are predicted to encode the alpha and beta subunits of glycine tRNA synthetase, a phosphatase, and an acyltransferase, respectively (GenBank accession no. DQ149607). act, cloned and expressed in Escherichia coli, produces acyl homoserine lactones (AHLs) principally of chain length C14 according to gas chromatography and mass spectrometry measurements. The AHLs have biological activity as shown by in vivo studies using the reporter strain Sinorhizobium meliloti Rm41 SinI-. Reverse transcription-PCR (RT-PCR) experiments indicate that the four genes are expressed as a single transcript, demonstrating that they constitute an operon. According to semiquantitative RT-PCR results, act is expressed more highly when A. ferrooxidans is grown in medium containing iron than when it is grown in medium containing sulfur. Since AHLs are important intercellular signaling molecules used by many bacteria to monitor their population density in quorum-sensing control of gene expression, this result suggests that A. ferrooxidans has two quorum-sensing systems, one based on Act, as described herein, and the other based on a Lux-like quorum-sensing system, reported previously. The latter system was shown to be upregulated in A. ferrooxidans grown in sulfur medium, suggesting that the two quorum-sensing systems respond to different environmental signals that may be related to their abilities to colonize and use different solid sulfur- and iron-containing minerals.

Differential expression of two bc1 complexes in the strict acidophilic chemolithoautotrophic bacterium Acidithiobacillus ferrooxidans suggests a model for their respective roles in iron or sulfur oxidation

Three strains of the strict acidophilic chemolithoautotrophic Acidithiobacillus ferrooxidans, including the type strain ATCC 23270, contain a petIIABC gene cluster that encodes the three proteins, cytochrome c1, cytochrome b and a Rieske protein, that constitute a bc1 electron-transfer complex. RT-PCR and Northern blotting show that the petIIABC cluster is co-transcribed with cycA, encoding a cytochrome c belonging to the c4 family, sdrA, encoding a putative short-chain dehydrogenase, and hip, encoding a high potential iron-sulfur protein, suggesting that the six genes constitute an operon, termed the petII operon. Previous results indicated that A. ferrooxidans contains a second pet operon, termed the petI operon, which contains a gene cluster that is similarly organized except that it lacks hip. Real-time PCR and Northern blot experiments demonstrate that petI is transcribed mainly in cells grown in medium containing iron, whereas petII is transcribed in cells grown in media containing sulfur or iron. Primer extension experiments revealed possible transcription initiation sites for the petI and petII operons. A model is presented in which petI is proposed to encode the bc1 complex, functioning in the uphill flow of electrons from iron to NAD(P), whereas petII is suggested to be involved in electron transfer from sulfur (or formate) to oxygen (or ferric iron). A. ferrooxidans is the only organism, to date, to exhibit two functional bc1 complexes.

Cortisol levels and anxiety-related behaviors in cattle

This investigation was conducted to examine the relationship between cortisol levels and anxiety-related behavior in nine Angus/Hereford cows. In a non-stressful pasture setting, measures were taken of rumination, distance from another cow, and body position. In a stressful holding pen situation, measures were taken of vocalizations and resistance to entering a squeeze chute. Cortisol was measured in serum samples. Subjects with high levels of cortisol spent less time ruminating (p=0.007) and vocalized more (p=0.07) than subjects with low levels of cortisol. Regardless of statistical significance, all mean differences were in the predicted direction. Cortisol levels were correlated with time spent ruminating (p=0.004) and later entrance to the squeeze chute (p=0.10). All correlations were also in the predicted direction. These findings provide consistent support for the link between cortisol and anxiety-related behavior in beef cattle, and the findings have theoretical and practical implications.

Thermophiles like hot T

A plethora of mechanisms confer protein stability in thermophilic microorganisms and, recently, it was suggested that these mechanisms might be divided along evolutionary lines. Here, a multi-genome comparison shows that there is a statistically significant increase in the proportion of NTN codons correlated with increasing optimal growth temperature for both Bacteria and Archaea. NTN encodes exclusively non-polar, hydrophobic amino acids and indicates a common underlying use of hydrophobicity for stabilizing proteins in Bacteria and Archaea that transcends evolutionary origins. However, some microorganisms do not follow this trend, suggesting that alternate mechanisms (e.g. intracellular electrolytes) might be used for protein stabilization. These studies highlight the usefulness of large-scale comparative genomics to uncover novel relationships that are not immediately obvious from protein structure studies alone.