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

Pangenome-level analysis of nucleoid-associated proteins in the Acidithiobacillia class: insights into their functional roles in mobile genetic elements biology

Mobile genetic elements (MGEs) are relevant agents in bacterial adaptation and evolutionary diversification. Stable appropriation of these DNA elements depends on host factors, among which are the nucleoid-associated proteins (NAPs). NAPs are highly abundant proteins that bind and bend DNA, altering its topology and folding, thus affecting all known cellular DNA processes from replication to expression. Even though NAP coding genes are found in most prokaryotic genomes, their functions in host chromosome biology and xenogeneic silencing are only known for a few NAP families. Less is known about the occurrence, abundance, and roles of MGE-encoded NAPs in foreign elements establishment and mobility. In this study, we used a combination of comparative genomics and phylogenetic strategies to gain insights into the diversity, distribution, and functional roles of NAPs within the class Acidithiobacillia with a special focus on their role in MGE biology. Acidithiobacillia class members are aerobic, chemolithoautotrophic, acidophilic sulfur-oxidizers, encompassing substantial genotypic diversity attributable to MGEs. Our search for NAP protein families (PFs) in more than 90 genomes of the different species that conform the class, revealed the presence of 1,197 proteins pertaining to 12 different NAP families, with differential occurrence and conservation across species. Pangenome-level analysis revealed 6 core NAP PFs that were highly conserved across the class, some of which also existed as variant forms of scattered occurrence, in addition to NAPs of taxa-restricted distribution. Core NAPs identified are reckoned as essential based on the conservation of genomic context and phylogenetic signals. In turn, various highly diversified NAPs pertaining to the flexible gene complement of the class, were found to be encoded in known plasmids or, larger integrated MGEs or, present in genomic loci associated with MGE-hallmark genes, pointing to their role in the stabilization/maintenance of these elements in strains and species with larger genomes. Both core and flexible NAPs identified proved valuable as markers, the former accurately recapitulating the phylogeny of the class, and the later, as seed in the bioinformatic identification of novel episomal and integrated mobile elements.

Genetic engineering of extremely acidophilic Acidithiobacillus species for biomining: Progress and perspectives

With global demands for mineral resources increasing and ore grades decreasing, microorganisms have been increasingly deployed in biomining applications to recover valuable metals particularly from normally considered waste, such as low-grade ores and used consumer electronics. Acidithiobacillus are a genus of chemolithoautotrophic extreme acidophiles that are commonly found in mining process waters and acid mine drainage, which have been reported in several studies to aid in metal recovery from bioremediation of metal-contaminated sites. Compared to conventional mineral processing technologies, biomining is often cited as a more sustainable and environmentally friendly process, but long leaching cycles and low extraction efficiency are main disadvantages that have hampered its industrial applications. Genetic engineering is a powerful technology that can be used to enhance the performance of microorganisms, such as Acidithiobacillus species. In this review, we compile existing data on Acidithiobacillus species' physiological traits and genomic characteristics, progresses in developing genetic tools to engineer them: plasmids, shutter vectors, transformation methods, selection markers, promoters and reporter systems developed, and genome editing techniques. We further propose genetic engineering strategies for enhancing biomining efficiency of Acidithiobacillus species and provide our perspectives on their future applications.

Prediction and Inferred Evolution of Acid Tolerance Genes in the Biotechnologically Important Acidihalobacter Genus

Acidihalobacter is a genus of acidophilic, gram-negative bacteria known for its ability to oxidize pyrite minerals in the presence of elevated chloride ions, a capability rare in other iron-sulfur oxidizing acidophiles. Previous research involving Acidihalobacter spp. has focused on their applicability in saline biomining operations and their genetic arsenal that allows them to cope with chloride, metal and oxidative stress. However, an understanding of the molecular adaptations that enable Acidihalobacter spp. to thrive under both acid and chloride stress is needed to provide a more comprehensive understanding of how this genus can thrive in such extreme biomining conditions. Currently, four genomes of the Acidihalobacter genus have been sequenced: Acidihalobacter prosperus DSM 5130^T, Acidihalobacter yilgarnensis DSM 105917^T, Acidihalobacter aeolianus DSM 14174^T, and Acidihalobacter ferrooxydans DSM 14175^T. Phylogenetic analysis shows that the Acidihalobacter genus roots to the Chromatiales class consisting of mostly halophilic microorganisms. In this study, we aim to advance our knowledge of the genetic repertoire of the Acidihalobacter genus that has enabled it to cope with acidic stress. We provide evidence of gene gain events that are hypothesized to help the Acidihalobacter genus cope with acid stress. Potential acid tolerance mechanisms that were found in the Acidihalobacter genomes include multiple potassium transporters, chloride/proton antiporters, glutamate decarboxylase system, arginine decarboxylase system, urease system, slp genes, squalene synthesis, and hopanoid synthesis. Some of these genes are hypothesized to have entered the Acidihalobacter via vertical decent from an inferred non-acidophilic ancestor, however, horizontal gene transfer (HGT) from other acidophilic lineages is probably responsible for the introduction of many acid resistance genes.

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.

Genomic evolution of the class Acidithiobacillia: deep-branching Proteobacteria living in extreme acidic conditions

Members of the genus Acidithiobacillus, now ranked within the class Acidithiobacillia, are model bacteria for the study of chemolithotrophic energy conversion under extreme conditions. Knowledge of the genomic and taxonomic diversity of Acidithiobacillia is still limited. Here, we present a systematic analysis of nearly 100 genomes from the class sampled from a wide range of habitats. Some of these genomes are new and others have been reclassified on the basis of advanced genomic analysis, thus defining 19 Acidithiobacillia lineages ranking at different taxonomic levels. This work provides the most comprehensive classification and pangenomic analysis of this deep-branching class of Proteobacteria to date. The phylogenomic framework obtained illuminates not only the evolutionary past of this lineage, but also the molecular evolution of relevant aerobic respiratory proteins, namely the cytochrome bo₃ ubiquinol oxidases.

Evolution of Type IV CRISPR-Cas Systems: Insights from CRISPR Loci in Integrative Conjugative Elements of Acidithiobacillia

Type IV CRISPR-Cas are a distinct variety of highly derived CRISPR-Cas systems that appear to have evolved from type III systems through the loss of the target-cleaving nuclease and partial deterioration of the large subunit of the effector complex. All known type IV CRISPR-Cas systems are encoded on plasmids, integrative and conjugative elements (ICEs), or prophages, and are thought to contribute to competition between these elements, although the mechanistic details of their function remain unknown. There is a clear parallel between the compositions and likely origin of type IV and type I systems recruited by Tn7-like transposons and mediating RNA-guided transposition. We investigated the diversity and evolutionary relationships of type IV systems, with a focus on those in Acidithiobacillia, where this variety of CRISPR is particularly abundant and always found on ICEs. Our analysis revealed remarkable evolutionary plasticity of type IV CRISPR-Cas systems, with adaptation and ancillary genes originating from different ancestral CRISPR-Cas varieties, and extensive gene shuffling within the type IV loci. The adaptation module and the CRISPR array apparently were lost in the type IV ancestor but were subsequently recaptured by type IV systems on several independent occasions. We demonstrate a high level of heterogeneity among the repeats with type IV CRISPR arrays, which far exceed the heterogeneity of any other known CRISPR repeats and suggest a unique adaptation mechanism. The spacers in the type IV arrays, for which protospacers could be identified, match plasmid genes, in particular those encoding the conjugation apparatus components. Both the biochemical mechanism of type IV CRISPR-Cas function and their role in the competition among mobile genetic elements remain to be investigated.

The adaptation mechanisms of Acidithiobacillus caldus CCTCC M 2018054 to extreme acid stress: Bioleaching performance, physiology, and transcriptomics

To understand the acid-resistant mechanism of bioleaching microorganism Acidithiobacillus caldus CCTCC M 2018054, its physiology and metabolic changes at the transcriptional level under extreme acid stress were systemically studied. Scanning electron microscopy (SEM), Fourier transform infrared reflection (FTIR) and X-ray diffraction (XRD) showed that with an increase in acidity, the absorption peak of sulfur oxidation-related functional groups such as S-O decreased significantly, and a dense sulfur passivation film appeared on the surface of the ore. Confocal laser scanning microscopy (CLSM) revealed that coverage scale of extracellular polymeric substance (EPS) and biofilm fluctuated accordingly along with the increasing acid stress (pH-stat 1.5, 1.2 0.9 and 0.6) during the bioleaching process. In response to acid stress, the increased levels of intracellular glutamic acid, alanine, cysteine, and proline contributed to the maintenance of intracellular pH homeostasis via decarboxylation and alkaline neutralization. Higher unsaturated fatty acid content was closely related to membrane fluidity. Up to 490 and 447 differentially expressed genes (DEGs) were identified at pH 1.5 vs pH 1.2 and pH 1.2 vs pH 0.9, respectively, and 177 common DEGs were associated with two-component system (TCS) regulation, transporter regulation, energy metabolism, and stress response. The upregulation of kdpB helped cells defend against proton invasion, whereas the downregulation of cysB and cbl implied stronger oxidation of sulfur compounds. The transcriptional level of sqr, sor, and soxA was significantly increased and consolidated the energy supply needed for resisting acid stress. Furthermore, eight of the identified DEGs (sor, cbl, ompA, atpF, nuoH, nuoC, sqr, grxB) were verified as being related to the acid stress response process. This study contributes toward expanding the application of these acidophiles in industrial bioleaching.

Respiratory Heme A-Containing Oxidases Originated in the Ancestors of Iron-Oxidizing Bacteria

Respiration is a major trait shaping the biology of many environments. Cytochrome oxidase containing heme A (COX) is a common terminal oxidase in aerobic bacteria and is the only one in mammalian mitochondria. The synthesis of heme A is catalyzed by heme A synthase (CtaA/Cox15), an enzyme that most likely coevolved with COX. The evolutionary origin of COX in bacteria has remained unknown. Using extensive sequence and phylogenetic analysis, we show that the ancestral type of heme A synthases is present in iron-oxidizing Proteobacteria such as Acidithiobacillus spp. These bacteria also contain a deep branching form of the major COX subunit (COX1) and an ancestral variant of CtaG, a protein that is specifically required for COX biogenesis. Our work thus suggests that the ancestors of extant iron-oxidizers were the first to evolve COX. Consistent with this conclusion, acidophilic iron-oxidizing prokaryotes lived on emerged land around the time for which there is the earliest geochemical evidence of aerobic respiration on earth. Hence, ecological niches of iron oxidation have apparently promoted the evolution of aerobic respiration.

Genomic adaptations enabling Acidithiobacillus distribution across wide-ranging hot spring temperatures and pHs

BACKGROUND

Terrestrial hot spring settings span a broad spectrum of physicochemistries. Physicochemical parameters, such as pH and temperature, are key factors influencing differences in microbial composition across diverse geothermal areas. Nonetheless, analysis of hot spring pools from the Taupo Volcanic Zone (TVZ), New Zealand, revealed that some members of the bacterial genus, Acidithiobacillus, are prevalent across wide ranges of hot spring pHs and temperatures. To determine the genomic attributes of Acidithiobacillus that inhabit such diverse conditions, we assembled the genomes of 19 uncultivated hot spring Acidithiobacillus strains from six geothermal areas and compared these to 37 publicly available Acidithiobacillus genomes from various habitats.

RESULTS

Analysis of 16S rRNA gene amplicons from 138 samples revealed that Acidithiobacillus comprised on average 11.4 ± 16.8% of hot spring prokaryotic communities, with three Acidithiobacillus amplicon sequence variants (ASVs) (TVZ_G1, TVZ_G2, TVZ_G3) accounting for > 90% of Acidithiobacillus in terms of relative abundance, and occurring in 126 out of 138 samples across wide ranges of temperature (17.5-92.9 °C) and pH (1.0-7.5). We recovered 19 environmental genomes belonging to each of these three ASVs, as well as a fourth related group (TVZ_G4). Based on genome average nucleotide identities, the four groups (TVZ_G1-TVZ_G4) constitute distinct species (ANI < 96.5%) of which three are novel Acidithiobacillus species (TVZ_G2-TVZ_G4) and one belongs to Acidithiobacillus caldus (TVZ_G1). All four TVZ Acidithiobacillus groups were found in hot springs with temperatures above the previously known limit for the genus (up to 40 °C higher), likely due to significantly higher proline and GC contents than other Acidithiobacillus species, which are known to increase thermostability. Results also indicate hot spring-associated Acidithiobacillus have undergone genome streamlining, likely due to thermal adaptation. Moreover, our data suggest that Acidithiobacillus prevalence across varied hot spring pHs is supported by distinct strategies, whereby TVZ_G2-TVZ_G4 regulate pH homeostasis mostly through Na⁺/H⁺ antiporters and proton-efflux ATPases, whereas TVZ_G1 mainly relies on amino acid decarboxylases.

CONCLUSIONS

This study provides insights into the distribution of Acidithiobacillus species across diverse hot spring physichochemistries and determines genomic features and adaptations that potentially enable Acidithiobacillus species to colonize a broad range of temperatures and pHs in geothermal environments. Video Abstract.

Pan-genome analysis of Riemerella anatipestifer reveals its genomic diversity and acquired antibiotic resistance associated with genomic islands

Riemerella anatipestifer is a gram-negative bacterium that leads to severe contagious septicemia in ducks, turkeys, chickens, and wild waterfowl. Here, a pan-genome with 32 R. anatipestifer genomes is re-established, and the mathematical model is calculated to evaluate the expansion of R. anatipestifer genomes, which were determined to be open. Average nucleotide identity (ANI) and phylogenetic analysis preliminarily clarify intraspecies variation and distance. Comparative genomic analysis of R. anatipestifer found that horizontal gene transfer events, which provide an expressway for the recruitment of novel functionalities and facilitate genetic diversity in microbial genomes, play a key role in the process of acquiring and transmitting antibiotic-resistance genes in R. anatipestifer. Furthermore, a new antibiotic-resistance gene cluster was identified in the same loci in 14 genomes. The uneven distribution of virulence factors was also confirmed by our results. Our study suggests that the ability to acquire foreign genes (such as antibiotic-resistance genes) increases the adaptability of R. anatipestifer, and the virulence genes with little mobility are highly conserved in R. anatipestifer.

Comparative microbiome analysis of two different long-term pesticide contaminated soils revealed the anthropogenic influence on functional potential of microbial communities

Microbial communities play a crucial role in bioremediation of pollutants in contaminated ecosystem. In addition to pure culture isolation and bacterial 16S rRNA based community studies, the focus has now shifted employing the omics technologies enormously for understanding the microbial diversity and functional potential of soil samples. Our previous report on two pesticide-contaminated sites revealed the diversity of both culturable and unculturable bacteria. In the present study, we have observed distinct taxonomic and functional communities in contaminated soil with respect to an uncontaminated soil as control by using shotgun metagenomic sequencing method. Our data demonstrated that Proteobacteria, Actinobacteria, Firmicutes, Bacteroidetes, and Acidobacteria significantly dominated the microbial diversity with their cumulative abundance percentage in the range of 98.61, 87.38, and 80.52 for Hindustan Insecticides Limited (HIL), India Pesticides Limited (IPL), and control respectively. Functional gene analysis demonstrated the presence of large number of both substrate specific upper pathway and common lower pathway degradative genes. Relatively lower number of genes was found encoding the degradation of styrene, atrazine, bisphenol, dioxin, and naphthalene. When three bacteria were augumentated with rhamnolipid (20-100 μM) and Triton X-100 (84-417 μM) surfactants in HIL soil, an enhanced degradation to 76%, 70%, and 58% of HCH, Endosulfan, and DDT respectively was achieved. The overall data obtained from two heavily contaminated soil suggest the versatility of the microbial communities for the xenobiotic pollutant degradation which may help in exploiting their potential applications in bioremediation.

Phylogeny, Divergent Evolution, and Speciation of Sulfur-Oxidizing Acidithiobacillus Populations

Background

Habitats colonized by acidophiles as an ideal physical barrier may induce genetic exchange of microbial members within the common communities, but little is known about how species in extremely acidic environments diverge and evolve.

Results

Using the acidophilic sulfur-oxidizer Acidithiobacillus as a case study, taxonomic reclassifications of many isolates provides novel insights into their phylogenetic lineage. Whole-genome-based comparisons were attempted to investigate the intra- and inter-species divergence. Recent studies clarified that functional and structural specificities of bacterial strains might provide opportunities for adaptive evolution responding to local environmental conditions. Acidophilic microorganisms play a key role in the acidification of natural waters and thus the formation of extremely acidic environments, and the feedbacks of the latter might confer the distinct evolutionary patterns of Acidithiobacillus spp. Varied horizontal gene transfer events occurred in different bacterial strains, probably resulting in the expansion of Acidithiobacillus genomes. Gene loss as another evolutionary force might cause the adaptive phenotypic diversity. A conceptual model for potential community-dependent evolutionary adaptation was thus proposed to illustrate the observed genome differentiation.

Conclusions

Collectively, the findings shed light on the phylogeny and divergent evolution of Acidithiobacillus strains, and provided a useful reference for evolutionary studies of other extremophiles.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5827-6) contains supplementary material, which is available to authorized users.

Nucleotide Second Messenger-Based Signaling in Extreme Acidophiles of the Acidithiobacillus Species Complex: Partition Between the Core and Variable Gene Complements

Cyclic and linear nucleotides are key elements of the signal transduction networks linking perception of the environment to specific cellular behavior of prokaryotes. These molecular mechanisms are particularly important in bacteria exposed to different, and frequently simultaneous, types of extreme conditions. This is the case in acidithiobacilli, a group of extremophilic bacteria thriving in highly acidic biotopes, that must also cope with significant variations in temperature, osmotic potentials and concentrations of various transition metals and metalloids. Environmental cues sensed by bacteria are transduced into differential levels of nucleotides acting as intracellular second messengers, promoting the activation or inhibition of target components and eliciting different output phenotypes. Cyclic (c) di-GMP, one of the most common bacterial second messengers, plays a key role in lifestyle changes in many bacteria, including acidithiobacilli. The presence of functional c-di-GMP-dependent signal transduction pathways in representative strains of the best-known linages of this species complex has been reported. However, a comprehensive panorama of the c-di-GMP modulated networks, the cognate input signals and output responses, are still missing for this group of extremophiles. Moreover, little fundamental understanding has been gathered for other nucleotides acting as second messengers. Taking advantage of the increasing number of sequenced genomes of the taxon, here we address the challenge of disentangling the nucleotide-driven signal transduction pathways in this group of polyextremophiles using comparative genomic tools and strategies. Results indicate that the acidithiobacilli possess all the genetic elements required to establish functional transduction pathways based in three different nucleotide-second messengers: (p)ppGpp, cyclic AMP (cAMP), and c-di-GMP. The elements related with the metabolism and transduction of (p)ppGpp and cAMP appear highly conserved, integrating signals related with nutrient starvation and polyphosphate metabolism, respectively. In contrast, c-di-GMP networks appear diverse and complex, differing both at the species and strain levels. Molecular elements of c-di-GMP metabolism and transduction were mostly found scattered along the flexible genome of the acidithiobacilli, allowing the identification of probable control modules that could be critical for substrate colonization, biofilm development and intercellular interactions. These may ultimately convey increased endurance to environmental stress and increased potential for gene sharing and adaptation to changing conditions.

The Type IV Secretion System of ICEAfe1: Formation of a Conjugative Pilus in Acidithiobacillus ferrooxidans

The dispersal of mobile genetic elements and their gene cargo relies on type IV secretion systems (T4SS). In this work the ICEAfe1 Tra-type T4SS nanomachine, encoded in the publicly available genome of Acidithiobacillus ferrooxidans ATCC 23270^TY, was characterized in terms of its organization, conservation, expression and mating bridge formation. Twenty-one conjugative genes grouped in four genetic clusters encode the ICEAfe1 T4SS, containing all the indispensable functions for the formation and stabilization of the pili and for DNA processing. The clusters' organization resembles that of other mobile genetic elements (such as plasmids and integrative and conjugative elements-ICEs). Sequence conservation, genetic organization and distribution of the tra system in the genomes of other sequenced Acidithiobacillus spp. suggests that the ICEAfe1 T4SS could mediate the lateral gene transfer between related bacteria. All ICEAfe1 T4SS genes are transcriptionally active and expressed from four independent operons. The transcriptional levels of selected marker genes increase in response to Mitomycin C treatment, a DNA damage elicitor that has acknowledged stimulatory effects on excision rates and gene expression of other ICEs, including ICEAfe1. Using a tailor-made pilin-antiserum against ICEAfe1 T4SS TraA pilin and epifluorescence microscopy, the presence of the conjugative pili on the cell surface of A. ferrooxidans could be demonstrated. Additionally, immunodetection assays, by immunogold, allowed the identification of pili-like extracellular structures. Together, the results obtained in this work demonstrate that the ICEAfe1 T4SS is phylogenetically conserved within the taxon, is expressed at mRNA and protein levels in vivo in the A. ferrooxidans type strain, and produces a pili-like structure of extracellular and intercellular localization in this model acidophile, supporting its functionality. Additional efforts will be required to prove conjugation of the ICEAfe1 or parts of this element through the cognate T4SS.

Sequence analysis and confirmation of the type IV pili-associated proteins PilY1, PilW and PilV in Acidithiobacillus thiooxidans

Acidithiobacillus thiooxidans is an acidophilic chemolithoautotrophic bacterium widely used in the mining industry due to its metabolic sulfur-oxidizing capability. The biooxidation of sulfide minerals is enhanced through the attachment of At. thiooxidans cells to the mineral surface. The Type IV pili (TfP) of At. thiooxidans may play an important role in the bacteria attachment since TfP play a key adhesive role in the attachment and colonization of different surfaces. In this work, we report for the first time the mRNA sequence of three TfP proteins from At. thiooxidans, the adhesin protein PilY1 and the TfP pilins PilW and PilV. The nucleotide sequences of these TfP proteins show changes in some nucleotide positions with respect to the corresponding annotated sequences. The bioinformatic analyses and 3D-modeling of protein structures sustain their classification as TfP proteins, as structural homologs of the corresponding proteins of Ps. aeruginosa, results that sustain the role of PilY1, PilW and PilV in pili assembly. Also, that PilY1 comprises the conserved Neisseria-PilC (superfamily) domain of the tip-associated adhesin, while PilW of the superfamily of putative TfP assembly proteins and PilV belongs to the superfamily of TfP assembly protein. In addition, the analyses suggested the presence of specific functional domains involved in adhesion, energy transduction and signaling functions. The phylogenetic analysis indicated that the PilY1 of Acidithiobacillus genus forms a cohesive group linked with iron- and/or sulfur-oxidizing microorganisms from acid mine drainage or mine tailings.

Comparative genomics of the Komagataeibacter strains-Efficient bionanocellulose producers

Komagataeibacter species are well-recognized bionanocellulose (BNC) producers. This bacterial genus, formerly assigned to Gluconacetobacter, is known for its phenotypic diversity manifested by strain-dependent carbon source preference, BNC production rate, pellicle structure, and strain stability. Here, we performed a comparative study of nineteen Komagataeibacter genomes, three of which were newly contributed in this work. We defined the core genome of the genus, clarified phylogenetic relationships among strains, and provided genetic evidence for the distinction between the two major clades, the K. xylinus and the K. hansenii. We found genomic traits, which likely contribute to the phenotypic diversity between the Komagataeibacter strains. These features include genome flexibility, carbohydrate uptake and regulation of its metabolism, exopolysaccharides synthesis, and the c-di-GMP signaling network. In addition, this work provides a comprehensive functional annotation of carbohydrate metabolism pathways, such as those related to glucose, glycerol, acetan, levan, and cellulose. Findings of this multi-genomic study expand understanding of the genetic variation within the Komagataeibacter genus and facilitate exploiting of its full potential for bionanocellulose production at the industrial scale.

Insights into the biology of acidophilic members of the Acidiferrobacteraceae family derived from comparative genomic analyses

The family Acidiferrobacteraceae (order Acidiferrobacterales) currently contains Gram negative, neutrophilic sulfur oxidizers such as Sulfuricaulis and Sulfurifustis, as well as acidophilic iron and sulfur oxidizers belonging to the Acidiferrobacter genus. The diversity and taxonomy of the genus Acidiferrobacter has remained poorly explored. Although several metagenome and bioleaching studies have identified its presence worldwide, only two strains, namely Acidiferrobacter thiooxydans DSM 2932^T, and Acidiferrobacter spp. SP3/III have been isolated and made publically available. Using 16S rRNA sequence data publically available for the Acidiferrobacteraceae, we herein shed light into the molecular taxonomy of this family. Results obtained support the presence of three clades Acidiferrobacter, Sulfuricaulis and Sulfurifustis. Genomic analyses of the genome sequences of A. thiooxydans^T and Acidiferrobacter spp. SP3/III indicate that ANI relatedness between the SPIII/3 strain and A. thiooxydans^T is below 95-96%, supporting the classification of strain SP3/III as a new species within this genus. In addition, approximately 70% of Acidiferrobacter sp. SPIII/3 predicted genes have a conserved ortholog in A. thiooxydans strains. A comparative analysis of iron, sulfur oxidation pathways, genome plasticity and cell-cell communication mechanisms of Acidiferrobacter spp. are also discussed.

Occurrence, integrity and functionality of AcaML1-like viruses infecting extreme acidophiles of the Acidithiobacillus species complex

General knowledge on the diversity and biology of microbial viruses infecting bacterial hosts from extreme acidic environments lags behind most other econiches. In this study, we analyse the AcaML1 virus occurrence in the taxon, its genetic composition and infective behaviour under standard acidic and SOS-inducing conditions to assess its integrity and functionality. Occurrence analysis in sequenced acidithiobacilli showed that AcaML1-like proviruses are confined to the mesothermophiles Acidithiobacillus caldus and Thermithiobacillus tepidarius. Among A. caldus strains and isolates this provirus had a modest prevalence (30%). Comparative genomic analysis revealed a significant conservation with the T. tepidarius AcaML1-like provirus, excepting the tail genes, and a high conservation of the virus across strains of the A. caldus species. Such conservation extends from the modules architecture to the gene level, suggesting that organization and composition of these viruses are preserved for functional reasons. Accordingly, the AcaML1 proviruses were demonstrated to excise from their host genomes under DNA-damaging conditions triggering the SOS-response and to produce DNA-containing VLPs. Despite this fact, under the conditions evaluated (acidic) the VLPs obtained from A. caldus ATCC 51756 could not produce productive infections of a candidate sensitive strain (#6) nor trigger it lysis.

Involvement of β-Carbonic Anhydrase Genes in Bacterial Genomic Islands and Their Horizontal Transfer to Protists

Genomic islands (GIs) are a type of mobile genetic element (MGE) that are present in bacterial chromosomes. They consist of a cluster of genes that produce proteins that contribute to a variety of functions, including, but not limited to, the regulation of cell metabolism, antimicrobial resistance, pathogenicity, virulence, and resistance to heavy metals. The genes carried in MGEs can be used as a trait reservoir in times of adversity. Transfer of genes using MGEs, occurring outside reproduction, is called horizontal gene transfer (HGT). Previous data have shown that numerous HGT events have occurred through endosymbiosis between prokaryotes and eukaryotes. β-Carbonic anhydrase (β-CA) enzymes play a critical role in the biochemical pathways of many prokaryotes and eukaryotes. We previously suggested the horizontal transfer of β-CA genes from plasmids of some prokaryotic endosymbionts to their protozoan hosts. In this study, we set out to identify β-CA genes that might have been transferred between prokaryotic and protist species through HGT in GIs. Therefore, we investigated prokaryotic chromosomes containing β-CA-encoding GIs and utilized multiple bioinformatics tools to reveal the distinct movements of β-CA genes among a wide variety of organisms. Our results identify the presence of β-CA genes in GIs of several medically and industrially relevant bacterial species, and phylogenetic analyses reveal multiple cases of likely horizontal transfer of β-CA genes from GIs of ancestral prokaryotes to protists.IMPORTANCE The evolutionary process is mediated by mobile genetic elements (MGEs), such as genomic islands (GIs). A gene or set of genes in the GIs is exchanged between and within various species through horizontal gene transfer (HGT). Based on the crucial role that GIs can play in bacterial survival and proliferation, they were introduced as environment- and pathogen-associated factors. Carbonic anhydrases (CAs) are involved in many critical biochemical pathways, such as the regulation of pH homeostasis and electrolyte transfer. Among the six evolutionary families of CAs, β-CA gene sequences are present in many bacterial species, which can be horizontally transferred to protists during evolution. This study shows the involvement of bacterial β-CA gene sequences in the GIs and suggests their horizontal transfer to protists during evolution.

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

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

Pan-Genome Analysis Links the Hereditary Variation of Leptospirillum ferriphilum With Its Evolutionary Adaptation

Niche adaptation has long been recognized to drive intra-species differentiation and speciation, yet knowledge about its relatedness with hereditary variation of microbial genomes is relatively limited. Using Leptospirillum ferriphilum species as a case study, we present a detailed analysis of genomic features of five recognized strains. Genome-to-genome distance calculation preliminarily determined the roles of spatial distance and environmental heterogeneity that potentially contribute to intra-species variation within L. ferriphilum species at the genome level. Mathematical models were further constructed to extrapolate the expansion of L. ferriphilum genomes (an 'open' pan-genome), indicating the emergence of novel genes with new sequenced genomes. The identification of diverse mobile genetic elements (MGEs) (such as transposases, integrases, and phage-associated genes) revealed the prevalence of horizontal gene transfer events, which is an important evolutionary mechanism that provides avenues for the recruitment of novel functionalities and further for the genetic divergence of microbial genomes. Comprehensive analysis also demonstrated that the genome reduction by gene loss in a broad sense might contribute to the observed diversification. We thus inferred a plausible explanation to address this observation: the community-dependent adaptation that potentially economizes the limiting resources of the entire community. Now that the introduction of new genes is accompanied by a parallel abandonment of some other ones, our results provide snapshots on the biological fitness cost of environmental adaptation within the L. ferriphilum genomes. In short, our genome-wide analyses bridge the relation between genetic variation of L. ferriphilum with its evolutionary adaptation.

Microbiomes in extremely acidic environments: functionalities and interactions that allow survival and growth of prokaryotes at low pH

Extremely acidic environments have global distribution and can have natural or, increasingly, anthropogenic origins. Extreme acidophiles grow optimally at pH 3 or less, have multiple strategies for tolerating stresses that accompany high levels of acidity and are scattered in all three domains of the tree of life. Metagenomic studies have expanded knowledge of the diversity of extreme acidophile communities, their ecological networks and their metabolic potentials, both confirmed and inferred. High resolution compositional and functional profiling of these microbiomes have begun to reveal spatial diversity patterns at global, regional, local, zonal and micro-scales. Future integration of genomic and other meta-omic data will offer new opportunities to utilize acidic microbiomes and to engineer beneficial interactions within them in biotechnologies.

A DNA segment encoding the anticodon stem/loop of tRNA determines the specific recombination of integrative-conjugative elements in Acidithiobacillus species

Horizontal gene transfer is crucial for the adaptation of microorganisms to environmental cues. The acidophilic, bioleaching bacterium Acidithiobacillus ferrooxidans encodes an integrative-conjugative genetic element (ICEAfe1) inserted in the gene encoding a tRNA^Ala. This genetic element is actively excised from the chromosome upon induction of DNA damage. A similar genetic element (ICEAca_TY.2) is also found in an equivalent position in the genome of Acidithiobacillus caldus. The local genomic context of both mobile genetic elements is highly syntenous and the cognate integrases are well conserved. By means of site directed mutagenesis, target site deletions and in vivo integrations assays in the heterologous model Escherichia coli, we assessed the target sequence requirements for site-specific recombination to be catalyzed by these integrases. We determined that each enzyme recognizes a specific small DNA segment encoding the anticodon stem/loop of the tRNA as target site and that specific positions in these regions are well conserved in the target attB sites of orthologous integrases. Also, we demonstrate that the local genetic context of the target sequence is not relevant for the integration to take place. These findings shed new light on the mechanism of site-specific integration of integrative-conjugative elements in members of Acidithiobacillus genus.

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.

Functionality of tRNAs encoded in a mobile genetic element from an acidophilic bacterium

The genome of the acidophilic, bioleaching bacterium Acidithiobacillus ferrooxidans, strain ATCC 23270, contains 95 predicted tRNA genes. Thirty-six of these genes (all 20 species) are clustered within an actively excising integrative-conjugative element (ICEAfe1). We speculated that these tRNA genes might have a role in adapting the bacterial tRNA pool to the codon usage of ICEAfe1 genes. To answer this question, we performed theoretical calculations of the global tRNA adaptation index to the entire A. ferrooxidans genome with and without the ICEAfe1 encoded tRNA genes. Based on these calculations, we observed that tRNAs encoded in ICEAfe1 negatively contribute to adapt the tRNA pool to the codon use in A. ferrooxidans. Although some of the tRNAs encoded in ICEAfe1 are functional in aminoacylation or protein synthesis, we found that they are expressed at low levels. These findings, along with the identification of a tRNA-like RNA encoded in the same cluster, led us to speculate that tRNA genes encoded in the mobile genetic element ICEAfe1 might have acquired mutations that would result in either inactivation or the acquisition of new functions.

Comparative Genomics Unravels the Functional Roles of Co-occurring Acidophilic Bacteria in Bioleaching Heaps

The spatial-temporal distribution of populations in various econiches is thought to be potentially related to individual differences in the utilization of nutrients or other resources, but their functional roles in the microbial communities remain elusive. We compared differentiation in gene repertoire and metabolic profiles, with a focus on the potential functional traits of three commonly recognized members (Acidithiobacillus caldus, Leptospirillum ferriphilum, and Sulfobacillus thermosulfidooxidans) in bioleaching heaps. Comparative genomics revealed that intra-species divergence might be driven by horizontal gene transfer. These co-occurring bacteria shared a few homologous genes, which significantly suggested the genomic differences between these organisms. Notably, relatively more genes assigned to the Clusters of Orthologous Groups category [G] (carbohydrate transport and metabolism) were identified in Sulfobacillus thermosulfidooxidans compared to the two other species, which probably indicated their mixotrophic capabilities that assimilate both organic and inorganic forms of carbon. Further inspection revealed distinctive metabolic capabilities involving carbon assimilation, nitrogen uptake, and iron-sulfur cycling, providing robust evidence for functional differences with respect to nutrient utilization. Therefore, we proposed that the mutual compensation of functionalities among these co-occurring organisms might provide a selective advantage for efficiently utilizing the limited resources in their habitats. Furthermore, it might be favorable to chemoautotrophs' lifestyles to form mutualistic interactions with these heterotrophic and/or mixotrophic acidophiles, whereby the latter could degrade organic compounds to effectively detoxify the environments. Collectively, the findings shed light on the genetic traits and potential metabolic activities of these organisms, and enable us to make some inferences about genomic and functional differences that might allow them to co-exist.

Adaptive Evolution of Extreme Acidophile Sulfobacillus thermosulfidooxidans Potentially Driven by Horizontal Gene Transfer and Gene Loss

Recent phylogenomic analysis has suggested that three strains isolated from different copper mine tailings around the world were taxonomically affiliated with Sulfobacillusthermosulfidooxidans Here, we present a detailed investigation of their genomic features, particularly with respect to metabolic potentials and stress tolerance mechanisms. Comprehensive analysis of the Sulfobacillus genomes identified a core set of essential genes with specialized biological functions in the survival of acidophiles in their habitats, despite differences in their metabolic pathways. The Sulfobacillus strains also showed evidence for stress management, thereby enabling them to efficiently respond to harsh environments. Further analysis of metabolic profiles provided novel insights into the presence of genomic streamlining, highlighting the importance of gene loss as a main mechanism that potentially contributes to cellular economization. Another important evolutionary force, especially in larger genomes, is gene acquisition via horizontal gene transfer (HGT), which might play a crucial role in the recruitment of novel functionalities. Also, a successful integration of genes acquired from archaeal donors appears to be an effective way of enhancing the adaptive capacity to cope with environmental changes. Taken together, the findings of this study significantly expand the spectrum of HGT and genome reduction in shaping the evolutionary history of Sulfobacillus strains.IMPORTANCE Horizontal gene transfer (HGT) and gene loss are recognized as major driving forces that contribute to the adaptive evolution of microbial genomes, although their relative importance remains elusive. The findings of this study suggest that highly frequent gene turnovers within microorganisms via HGT were necessary to incur additional novel functionalities to increase the capacity of acidophiles to adapt to changing environments. Evidence also reveals a fascinating phenomenon of potential cross-kingdom HGT. Furthermore, genome streamlining may be a critical force in driving the evolution of microbial genomes. Taken together, this study provides insights into the importance of both HGT and gene loss in the evolution and diversification of bacterial genomes.

Molecular Systematics of the Genus Acidithiobacillus: Insights into the Phylogenetic Structure and Diversification of the Taxon

The acidithiobacilli are sulfur-oxidizing acidophilic bacteria that thrive in both natural and anthropogenic low pH environments. They contribute to processes that lead to the generation of acid rock drainage in several different geoclimatic contexts, and their properties have long been harnessed for the biotechnological processing of minerals. Presently, the genus is composed of seven validated species, described between 1922 and 2015: Acidithiobacillus thiooxidans, A. ferrooxidans, A. albertensis, A. caldus, A. ferrivorans, A. ferridurans, and A. ferriphilus. However, a large number of Acidithiobacillus strains and sequence clones have been obtained from a variety of ecological niches over the years, and many isolates are thought to vary in phenotypic properties and cognate genetic traits. Moreover, many isolates remain unclassified and several conflicting specific assignments muddle the picture from an evolutionary standpoint. Here we revise the phylogenetic relationships within this species complex and determine the phylogenetic species boundaries using three different typing approaches with varying degrees of resolution: 16S rRNA gene-based ribotyping, oligotyping, and multi-locus sequencing analysis (MLSA). To this end, the 580 16S rRNA gene sequences affiliated to the Acidithiobacillus spp. were collected from public and private databases and subjected to a comprehensive phylogenetic analysis. Oligotyping was used to profile high-entropy nucleotide positions and resolve meaningful differences between closely related strains at the 16S rRNA gene level. Due to its greater discriminatory power, MLSA was used as a proxy for genome-wide divergence in a smaller but representative set of strains. Results obtained indicate that there is still considerable unexplored diversity within this genus. At least six new lineages or phylotypes, supported by the different methods used herein, are evident within the Acidithiobacillus species complex. Although the diagnostic characteristics of these subgroups of strains are as yet unresolved, correlations to specific metadata hint to the mechanisms behind econiche-driven divergence of some of the species/phylotypes identified. The emerging phylogenetic structure for the genus outlined in this study can be used to guide isolate selection for future population genomics and evolutionary studies in this important acidophile model.

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).

Gene Turnover Contributes to the Evolutionary Adaptation of Acidithiobacillus caldus: Insights from Comparative Genomics

Acidithiobacillus caldus is an extremely acidophilic sulfur-oxidizer with specialized characteristics, such as tolerance to low pH and heavy metal resistance. To gain novel insights into its genetic complexity, we chosen six A. caldus strains for comparative survey. All strains analyzed in this study differ in geographic origins as well as in ecological preferences. Based on phylogenomic analysis, we clustered the six A. caldus strains isolated from various ecological niches into two groups: group 1 strains with smaller genomes and group 2 strains with larger genomes. We found no obvious intraspecific divergence with respect to predicted genes that are related to central metabolism and stress management strategies between these two groups. Although numerous highly homogeneous genes were observed, high genetic diversity was also detected. Preliminary inspection provided a first glimpse of the potential correlation between intraspecific diversity at the genome level and environmental variation, especially geochemical conditions. Evolutionary genetic analyses further showed evidence that the difference in environmental conditions might be a crucial factor to drive the divergent evolution of A. caldus species. We identified a diverse pool of mobile genetic elements including insertion sequences and genomic islands, which suggests a high frequency of genetic exchange in these harsh habitats. Comprehensive analysis revealed that gene gains and losses were both dominant evolutionary forces that directed the genomic diversification of A. caldus species. For instance, horizontal gene transfer and gene duplication events in group 2 strains might contribute to an increase in microbial DNA content and novel functions. Moreover, genomes undergo extensive changes in group 1 strains such as removal of potential non-functional DNA, which results in the formation of compact and streamlined genomes. Taken together, the findings presented herein show highly frequent gene turnover of A. caldus species that inhabit extremely acidic environments, and shed new light on the contribution of gene turnover to the evolutionary adaptation of acidophiles.

Survey of the green picoalga Bathycoccus genomes in the global ocean

Bathycoccus is a cosmopolitan green micro-alga belonging to the Mamiellophyceae, a class of picophytoplankton that contains important contributors to oceanic primary production. A single species of Bathycoccus has been described while the existence of two ecotypes has been proposed based on metagenomic data. A genome is available for one strain corresponding to the described phenotype. We report a second genome assembly obtained by a single cell genomics approach corresponding to the second ecotype. The two Bathycoccus genomes are divergent enough to be unambiguously distinguishable in whole DNA metagenomic data although they possess identical sequence of the 18S rRNA gene including in the V9 region. Analysis of 122 global ocean whole DNA metagenome samples from the Tara-Oceans expedition reveals that populations of Bathycoccus that were previously identified by 18S rRNA V9 metabarcodes are only composed of these two genomes. Bathycoccus is relatively abundant and widely distributed in nutrient rich waters. The two genomes rarely co-occur and occupy distinct oceanic niches in particular with respect to depth. Metatranscriptomic data provide evidence for gain or loss of highly expressed genes in some samples, suggesting that the gene repertoire is modulated by environmental conditions.

Comparative genomics unravels metabolic differences at the species and/or strain level and extremely acidic environmental adaptation of ten bacteria belonging to the genus Acidithiobacillus

Members of the Acidithiobacillus genus are widely found in extreme environments characterized by low pH and high concentrations of toxic substances, thus it is necessary to identify the cellular mechanisms needed to cope with these harsh conditions. Pan-genome analysis of ten bacteria belonging to the genus Acidithiobacillus suggested the existence of core genome, most of which were assigned to the metabolism-associated genes. Additionally, the unique genes of Acidithiobacillus ferrooxidans were much less than those of other species. A large proportion of Acidithiobacillus ferrivorans-specific genes were mapped especially to metabolism-related genes, indicating that diverse metabolic pathways might confer an advantage for adaptation to local environmental conditions. Analyses of functional metabolisms revealed the differences of carbon metabolism, nitrogen metabolism, and sulfur metabolism at the species and/or strain level. The findings also showed that Acidithiobacillus spp. harbored specific adaptive mechanisms for thriving under extreme environments. The genus Acidithiobacillus had the genetic potential to resist and metabolize toxic substances such as heavy metals and organic solvents. Comparison across species and/or strains of Acidithiobacillus populations provided a deeper appreciation of metabolic differences and environmental adaptation, as well as highlighting the importance of cellular mechanisms that maintain the basal physiological functions under complex acidic environmental conditions.

Comparative Genomics of the Extreme Acidophile Acidithiobacillus thiooxidans Reveals Intraspecific Divergence and Niche Adaptation

Acidithiobacillus thiooxidans known for its ubiquity in diverse acidic and sulfur-bearing environments worldwide was used as the research subject in this study. To explore the genomic fluidity and intraspecific diversity of Acidithiobacillus thiooxidans (A. thiooxidans) species, comparative genomics based on nine draft genomes was performed. Phylogenomic scrutiny provided first insights into the multiple groupings of these strains, suggesting that genetic diversity might be potentially correlated with their geographic distribution as well as geochemical conditions. While these strains shared a large number of common genes, they displayed differences in gene content. Functional assignment indicated that the core genome was essential for microbial basic activities such as energy acquisition and uptake of nutrients, whereas the accessory genome was thought to be involved in niche adaptation. Comprehensive analysis of their predicted central metabolism revealed that few differences were observed among these strains. Further analyses showed evidences of relevance between environmental conditions and genomic diversification. Furthermore, a diverse pool of mobile genetic elements including insertion sequences and genomic islands in all A. thiooxidans strains probably demonstrated the frequent genetic flow (such as lateral gene transfer) in the extremely acidic environments. From another perspective, these elements might endow A. thiooxidans species with capacities to withstand the chemical constraints of their natural habitats. Taken together, our findings bring some valuable data to better understand the genomic diversity and econiche adaptation within A. thiooxidans strains.

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.

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.

Toxin-antitoxin systems in the mobile genome of Acidithiobacillus ferrooxidans

Toxin-antitoxin (TA) systems are genetic modules composed of a pair of genes encoding a stable toxin and an unstable antitoxin that inhibits toxin activity. They are widespread among plasmids and chromosomes of bacteria and archaea. TA systems are known to be involved in the stabilization of plasmids but there is no consensus about the function of chromosomal TA systems. To shed light on the role of chromosomally encoded TA systems we analyzed the distribution and functionality of type II TA systems in the chromosome of two strains from Acidithiobacillus ferrooxidans (ATCC 23270 and 53993), a Gram-negative, acidophilic, environmental bacterium that participates in the bioleaching of minerals. As in other environmental microorganisms, A. ferrooxidans has a high content of TA systems (28-29) and in twenty of them the toxin is a putative ribonuclease. According to the genetic context, some of these systems are encoded near or within mobile genetic elements. Although most TA systems are shared by both strains, four of them, which are encoded in the active mobile element ICEAfe1, are exclusive to the type strain ATCC 23270. We demostrated that two TA systems from ICEAfe1 are functional in E. coli cells, since the toxins inhibit growth and the antitoxins counteract the effect of their cognate toxins. All the toxins from ICEAfe1, including a novel toxin, are RNases with different ion requirements. The data indicate that some of the chromosomally encoded TA systems are actually part of the A. ferrooxidans mobile genome and we propose that could be involved in the maintenance of these integrated mobile genetic elements.

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.

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

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

Metal resistance in acidophilic microorganisms and its significance for biotechnologies

Extremely acidophilic microorganisms have an optimal pH of <3 and are found in all three domains of life. As metals are more soluble at acid pH, acidophiles are often challenged by very high metal concentrations. Acidophiles are metal-tolerant by both intrinsic, passive mechanisms as well as active systems. Passive mechanisms include an internal positive membrane potential that creates a chemiosmotic gradient against which metal cations must move, as well as the formation of metal sulfate complexes reducing the concentration of the free metal ion. Active systems include efflux proteins that pump metals out of the cytoplasm and conversion of the metal to a less toxic form. Acidophiles are exploited in a number of biotechnologies including biomining for sulfide mineral dissolution, biosulfidogenesis to produce sulfide that can selectively precipitate metals from process streams, treatment of acid mine drainage, and bioremediation of acidic metal-contaminated milieux. This review describes how acidophilic microorganisms tolerate extremely high metal concentrations in biotechnological processes and identifies areas of future work that hold promise for improving the efficiency of these applications.