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

Eco-evolutionary strategies for relieving carbon limitation under salt stress differ across microbial clades

With the continuous expansion of saline soils under climate change, understanding the eco-evolutionary tradeoff between the microbial mitigation of carbon limitation and the maintenance of functional traits in saline soils represents a significant knowledge gap in predicting future soil health and ecological function. Through shotgun metagenomic sequencing of coastal soils along a salinity gradient, we show contrasting eco-evolutionary directions of soil bacteria and archaea that manifest in changes to genome size and the functional potential of the soil microbiome. In salt environments with high carbon requirements, bacteria exhibit reduced genome sizes associated with a depletion of metabolic genes, while archaea display larger genomes and enrichment of salt-resistance, metabolic, and carbon-acquisition genes. This suggests that bacteria conserve energy through genome streamlining when facing salt stress, while archaea invest in carbon-acquisition pathways to broaden their resource usage. These findings suggest divergent directions in eco-evolutionary adaptations to soil saline stress amongst microbial clades and serve as a foundation for understanding the response of soil microbiomes to escalating climate change.

A cascade of sulfur transferases delivers sulfur to the sulfur-oxidizing heterodisulfide reductase-like complex

A heterodisulfide reductase-like complex (sHdr) and novel lipoate-binding proteins (LbpAs) are central players of a wide-spread pathway of dissimilatory sulfur oxidation. Bioinformatic analysis demonstrate that the cytoplasmic sHdr-LbpA systems are always accompanied by sets of sulfur transferases (DsrE proteins, TusA, and rhodaneses). The exact composition of these sets may vary depending on the organism and sHdr system type. To enable generalizations, we studied model sulfur oxidizers from distant bacterial phyla, that is, Aquificota and Pseudomonadota. DsrE3C of the chemoorganotrophic Alphaproteobacterium Hyphomicrobium denitrificans and DsrE3B from the Gammaproteobacteria Thioalkalivibrio sp. K90mix, an obligate chemolithotroph, and Thiorhodospira sibirica, an obligate photolithotroph, are homotrimers that donate sulfur to TusA. Additionally, the hyphomicrobial rhodanese-like protein Rhd442 exchanges sulfur with both TusA and DsrE3C. The latter is essential for sulfur oxidation in Hm. denitrificans. TusA from Aquifex aeolicus (AqTusA) interacts physiologically with AqDsrE, AqLbpA, and AqsHdr proteins. This is particularly significant as it establishes a direct link between sulfur transferases and the sHdr-LbpA complex that oxidizes sulfane sulfur to sulfite. In vivo, it is unlikely that there is a strict unidirectional transfer between the sulfur-binding enzymes studied. Rather, the sulfur transferases form a network, each with a pool of bound sulfur. Sulfur flux can then be shifted in one direction or the other depending on metabolic requirements. A single pair of sulfur-binding proteins with a preferred transfer direction, such as a DsrE3-type protein towards TusA, may be sufficient to push sulfur into the sink where it is further metabolized or needed.

Gene loss and relaxed selection of plaat1 in vertebrates adapted to low-light environments

Gene loss is an important mechanism for evolution in low-light or cave environments where visual adaptations often involve a reduction or loss of eyesight. The plaat gene family encodes phospholipases essential for the degradation of organelles in the lens of the eye. These phospholipases translocate to damaged organelle membranes, inducing them to rupture. This rupture is required for lens transparency and is essential for developing a functioning eye. Plaat3 is thought to be responsible for this role in mammals, while plaat1 is thought to be responsible in other vertebrates. We used a macroevolutionary approach and comparative genomics to examine the origin, loss, synteny and selection of plaat1 across bony fishes and tetrapods. We showed that plaat1 (probably ancestral to all bony fish + tetrapods) has been lost in squamates and is significantly degraded in lineages of low-visual-acuity and blind mammals and fishes. Our findings suggest that plaat1 is important for visual acuity across bony vertebrates, and that its loss through relaxed selection and pseudogenization may have played a role in the repeated evolution of visual systems in low-light environments. Our study sheds light on the importance of gene-loss in trait evolution and provides insights into the mechanisms underlying visual acuity in low-light environments.

Gene loss and relaxed selection of plaat1 in vertebrates adapted to low-light environments

Gene loss is an important mechanism for evolution in low-light or cave environments where visual adaptations often involve a reduction or loss of eyesight. The plaat gene family are phospholipases essential for the degradation of organelles in the lens of the eye. They translocate to damaged organelle membranes, inducing them to rupture. This rupture is required for lens transparency and is essential for developing a functioning eye. Plaat3 is thought to be responsible for this role in mammals, while plaat1 is thought to be responsible in other vertebrates. We used a macroevolutionary approach and comparative genomics to examine the origin, loss, synteny, and selection of plaat1 across bony fishes and tetrapods. We show that plaat1 (likely ancestral to all bony fish + tetrapods) has been lost in squamates and is significantly degraded in lineages of low-visual acuity and blind mammals and fish. Our findings suggest that plaat1 is important for visual acuity across bony vertebrates, and that its loss through relaxed selection and pseudogenization may have played a role in the repeated evolution of visual systems in low-light-environments. Our study sheds light on the importance of gene-loss in trait evolution and provides insights into the mechanisms underlying visual acuity in low-light environments.

Trade-off in genome turnover events leading to adaptive evolution of Microcystis aeruginosa species complex

BACKGROUND

Numerous studies in the past have expanded our understanding of the genetic differences of global distributed cyanobacteria that originated around billions of years ago, however, unraveling how gene gain and loss drive the genetic evolution of cyanobacterial species, and the trade-off of these evolutionary forces are still the central but poorly understood issues.

RESULTS

To delineate the contribution of gene flow in mediating the hereditary differentiation and shaping the microbial evolution, a global genome-wide study of bloom-forming cyanobacterium, Microcystis aeruginosa species complex, provided robust evidence for genetic diversity, reflected by enormous variation in gene repertoire among various strains. Mathematical extrapolation showed an 'open' microbial pan-genome of M. aeruginosa species, since novel genes were predicted to be introduced after new genomes were sequenced. Identification of numerous horizontal gene transfer's signatures in genome regions of interest suggested that genome expansion via transformation and phage-mediated transduction across bacterial lineage as an evolutionary route may contribute to the differentiation of Microcystis functions (e.g., carbohydrate metabolism, amino acid metabolism, and energy metabolism). Meanwhile, the selective loss of some dispensable genes at the cost of metabolic versatility is as a mean of adaptive evolution that has the potential to increase the biological fitness.

CONCLUSIONS

Now that the recruitment of novel genes was accompanied by a parallel loss of some other ones, a trade-off in gene content may drive the divergent differentiation of M. aeruginosa genomes. Our study provides a genetic framework for the evolution of M. aeruginosa species and illustrates their possible evolutionary patterns.

Insertion sequence contributes to the evolution and environmental adaptation of Acidithiobacillus

BACKGROUND

The genus Acidithiobacillus has been widely concerned due to its superior survival and oxidation ability in acid mine drainage (AMD). However, the contribution of insertion sequence (IS) to their biological evolution and environmental adaptation is very limited. ISs are the simplest kinds of mobile genetic elements (MGEs), capable of interrupting genes, operons, or regulating the expression of genes through transposition activity. ISs could be classified into different families with their own members, possessing different copies.

RESULTS

In this study, the distribution and evolution of ISs, as well as the functions of the genes around ISs in 36 Acidithiobacillus genomes, were analyzed. The results showed that 248 members belonging to 23 IS families with a total of 10,652 copies were identified within the target genomes. The IS families and copy numbers among each species were significantly different, indicating that the IS distribution of Acidithiobacillus were not even. A. ferrooxidans had 166 IS members, which may develop more gene transposition strategies compared with other Acidithiobacillus spp. What's more, A. thiooxidans harbored the most IS copies, suggesting that their ISs were the most active and more likely to transpose. The ISs clustered in the phylogenetic tree approximately according to the family, which were mostly different from the evolutionary trends of their host genomes. Thus, it was suggested that the recent activity of ISs of Acidithiobacillus was not only determined by their genetic characteristics, but related with the environmental pressure. In addition, many ISs especially Tn3 and IS110 families were inserted around the regions whose functions were As/Hg/Cu/Co/Zn/Cd translocation and sulfur oxidation, implying that ISs could improve the adaptive capacities of Acidithiobacillus to the extremely acidic environment by enhancing their resistance to heavy metals and utilization of sulfur.

CONCLUSIONS

This study provided the genomic evidence for the contribution of IS to evolution and adaptation of Acidithiobacillus, opening novel sights into the genome plasticity of those acidophiles.

Responses of microbial community to geochemical parameters on vertical depth in bioheap system of low-grade copper sulfide

Monitoring of the microbial community in bioleaching system is essential for control process parameters and enhance the leaching efficiency. Due to the difficulty of sampling, microbial distribution, community succession and bioleaching activity along the vertical depth of bioleaching heaps remain unresolved. This study investigated the geochemical parameters and microbial community structure along a depth profile in a bioleaching heap and leachate. 80 ore samples at different heap depths and 9 leaching solution samples from three bioheaps of Zijin Copper Mine were collected. Microbial composition, mineral types and geochemical parameters of these samples were analyzed by 16S rRNA high-throughput sequencing and a series of chemical measurement technologies. The results revealed that the pH, Cu, Fe and the total sulfur contents were the major factors shaping the composition of the microbial communities in the bioleaching system. The extent of mineral oxidation increased as the sample depth increases, followed by the increasing of sulfur oxidizers. The abundance of sulfur and iron oxidizers including members of Acidithiobacillus, Sulfobacillus and Acidiferrobacter were significantly higher in the leaching heap than in the leaching solution, meanwhile, they showed strong positive interactions with other members within the same genera and iron oxidizer Leptospirillum and Ferroplasma. Besides, Acidithiobacillus negatively interacted with heterotrophs such as Sphingobium, Exiguobacterium, Brevundimonas and so on. On the contrast, members of Leptospirillum and unclassified Archaea were significantly abundant in the leaching solution and revealed strong interactions with members of Thermoplasmatales. The main conclusion of this study, especially the leaching potential of microorganisms prevailing in bioheaps and their relationships with geochemical factors, provides theoretical guidance for future process design such as the control of processing parameters and microbial community in heap leaching.

Reduction, evolutionary pattern and positive selection of genes encoding formate dehydrogenase in Wood-Ljungdahl pathway of gastrointestinal acetogens suggests their adaptation to formate-rich habitats

Acetogens are anaerobes using Wood-Ljungdahl pathway (WLP) as the terminal electron acceptor for both assimilation and dissimilation of CO₂ and widely distributed in diverse habitats. However, their habitat adaptation is often unclear. Given that bacterial genome evolution is often the result of environmental selective pressure, hereby we analysed gene copy number, phylogeny and selective pressure of genes involved in WLP within known genomes of 43 species to study the habitat adaption of gastrointestinal acetogens. The gene copy number of formate dehydrogenase (FDH) in gastrointestinal acetogens was much lower than that of non-gastrointestinal acetogens, and in five cases, no FDH genes were found in the genomes of five gastrointestinal acetogens, but that of the other WLP genes showed no difference. The evolutionary pattern of FDH genes was significantly different from that of the other enzymes. Additionally, seven positively selected sites were only identified in the fdhF genes, which means fdhF mutations favoured their adaptation. Collectively, reduction or loss of FDH genes and their evolutionary pattern as well as positive selection in gastrointestinal acetogens indicated their adaptation to formate-rich habitats, implying that FDH genes catalysing CO₂ reduction to formate as the first step of methyl branch of WLP may have evolved independently.

Convergent evolution and horizontal gene transfer in Arctic Ocean microalgae

Microbial communities in the world ocean are affected strongly by oceanic circulation, creating characteristic marine biomes. The high connectivity of most of the ocean makes it difficult to disentangle selective retention of colonizing genotypes (with traits suited to biome specific conditions) from evolutionary selection, which would act on founder genotypes over time. The Arctic Ocean is exceptional with limited exchange with other oceans and ice covered since the last ice age. To test whether Arctic microalgal lineages evolved apart from algae in the global ocean, we sequenced four lineages of microalgae isolated from Arctic waters and sea ice. Here we show convergent evolution and highlight geographically limited HGT as an ecological adaptive force in the form of PFAM complements and horizontal acquisition of key adaptive genes. Notably, ice-binding proteins were acquired and horizontally transferred among Arctic strains. A comparison with Tara Oceans metagenomes and metatranscriptomes confirmed mostly Arctic distributions of these IBPs. The phylogeny of Arctic-specific genes indicated that these events were independent of bacterial-sourced HGTs in Antarctic Southern Ocean microalgae.

Multi-scale analysis of nickel ion tolerance mechanism for thermophilic Sulfobacillus thermosulfidooxidans in bioleaching

Bioleaching is intensively investigated for recovering valuable metals such as Li, Co, Ni and Cu. Nickel ion stress threatens the health of microorganisms when Ni²⁺ starts to accumulate in the leachate during the bioleaching of materials that are rich in Ni, such as spent lithium-ion batteries. The possible mechanisms underlying the response of S. thermosulfidooxidans to nickel ion stress were analyzed using a multi-scale approach. Under the condition of nickel ion stress, high concentrations of nickel ions were immobilized by extracellular polymeric substances, while concentrations of nickel ions inside the cells remained low. The intracellular adenosine triphosphate (ATP) concentration and H⁺-ATPase activity increased to maintain normal cell growth and metabolic activities. Scavenging abilities of S. thermosulfidooxidans for hydrogen peroxide and superoxide anion were enhanced to reduce oxidative damage induced by nickel ion stress. There were 734 differentially expressed genes identified by RNA-seq under nickel ion stress. Most of them were involved in oxidative phosphorylation, glutathione metabolism and genetic information processing, responsible for intracellular energy utilization, intracellular antioxidant capacity and DNA damage repair, respectively. The results of this study are of major significance for in-depth understanding of the mechanisms of acidophilic microorganisms' resistance to metal ions.

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.

Phylogenetically divergent bacteria consortium from neutral activated sludge showed heightened potential on bioleaching spent lithium-ion batteries

Recycling of spent lithium-ion batteries (LIBs) has become a global issue because of the potential environment risks raised by spent LIBs as well as high valuable metal content remaining in them. Although bioleaching is an environmentally friendly method to recover metals from spent LIBs, the commonly utilized bioleaching bacterial consortia or strains enriched/isolated from acidic environments cannot be applied at large scales owing to their long leaching cycle and poor tolerance to organic compounds. Here, two bioleaching consortia were enriched in 60 days from neutral activated sludge and were identified phylogenetically divergent from the documented bioleaching bacteria. The results showed that the novel consortia shortened the leaching cycle almost by half when compared to the previous reported consortia or strains, of which one consortium dominated by Acidithiobacillus ferrooxidans displayed high bioleaching efficiency on LiMn₂O₄, as 69.46% lithium (Li) and 67.60% manganese (Mn) were leached out in seven days. This consortium was further domesticated using cathodic materials for 100 days and proved consisted of three mixotrophs and two chemoautotrophs, three of which were novel species from the genera Sulfobacillus and Leptospirillum. More genes coding for proteins that utilize organic compounds were annotated in the metagenomic assembled genomes (MAGs) than previously reported. A mutualistic relationship between mixotrophs and chemoautotrophs was suggested to help the consortium surviving under either organic- rich or shortage environments. The results discovered that novel bioleaching bacteria with shorter leaching cycle and higher tolerance to organics could be enriched from non-acidic environments, which showed high potential for the metal recovering from spent LIBs or other organic-rich environments.

Recent progress in the application of omics technologies in the study of bio-mining microorganisms from extreme environments

Bio-mining microorganisms are a key factor affecting the metal recovery rate of bio-leaching, which inevitably produces an extremely acidic environment. As a powerful tool for exploring the adaptive mechanisms of microorganisms in extreme environments, omics technologies can greatly aid our understanding of bio-mining microorganisms and their communities on the gene, mRNA, and protein levels. These omics technologies have their own advantages in exploring microbial diversity, adaptive evolution, changes in metabolic characteristics, and resistance mechanisms of single strains or their communities to extreme environments. These technologies can also be used to discover potential new genes, enzymes, metabolites, metabolic pathways, and species. In addition, integrated multi-omics analysis can link information at different biomolecular levels, thereby obtaining more accurate and complete global adaptation mechanisms of bio-mining microorganisms. This review introduces the current status and future trends in the application of omics technologies in the study of bio-mining microorganisms and their communities in extreme environments.

Microbial sulfur metabolism and environmental implications

Sulfur as a macroelement plays an important role in biochemistry in both natural environments and engineering biosystems, which can be further linked to other important element cycles, e.g. carbon, nitrogen and iron. Consequently, the sulfur cycling primarily mediated by sulfur compounds oxidizing microorganisms and sulfur compounds reducing microorganisms has enormous environmental implications, particularly in wastewater treatment and pollution bioremediation. In this review, to connect the knowledge in microbial sulfur metabolism to environmental applications, we first comprehensively review recent advances in understanding microbial sulfur metabolisms at molecular-, cellular- and ecosystem-levels, together with their energetics. We then discuss the environmental implications to fight against soil and water pollution, with four foci: (1) acid mine drainage, (2) water blackening and odorization in urban rivers, (3) SANI® and DS-EBPR processes for sewage treatment, and (4) bioremediation of persistent organic pollutants. In addition, major challenges and further developments toward elucidation of microbial sulfur metabolisms and their environmental applications are identified and discussed.

Sulfobacillus harzensis sp. nov., an acidophilic bacterium inhabiting mine tailings from a polymetallic mine

A mixotrophic and acidophilic bacterial strain BGR 140^T was isolated from mine tailings in the Harz Mountains near Goslar, Germany. Cells of BGR 140^T were Gram-stain-positive, endospore-forming, motile and rod-shaped. BGR 140^T grew aerobically at 25–55 °C (optimum 45 °C) and at pH 1.5–5.0 (optimum pH 3.0). The results of analysis of the 16S rRNA gene sequences indicated that BGR 140^T was phylogenetically related to different members of the genus Sulfobacillus, and the sequence identities to Sulfobacillus acidophilus DSM 10332^T, Sulfobacillus thermotolerans DSM 17362^T, and Sulfobacillus benefaciens DSM 19468^T were 94.8, 91.8 and 91.6 %, respectively. Its cell wall peptidoglycan is A1γ, composed of meso-diaminopimelic acid. The respiratory quinone is DMK-6. The major polar lipids were determined to be glycolipid, phospholipid and phosphatidylglycerol. The predominant fatty acid is 11-cycloheptanoyl-undecanoate. The genomic DNA G+C content is 58.2 mol%. On the basis of the results of phenotypic and genomic analyses, it is concluded that strain BGR 140^T represents a novel species of the genus Sulfobacillus, for which the name Sulfobacillus harzensis sp. nov. is proposed because of its origin. Its type strain is BGR 140^T (=DSM 109850^T=JCM 39070^T).

Genome Mining and Comparative Genome Analysis Revealed Niche-Specific Genome Expansion in Antibacterial Bacillus pumilus Strain SF-4

The present study reports the isolation of antibacterial exhibiting Bacillus pumilus (B. pumilus) SF-4 from soil field. The genome of this strain SF-4 was sequenced and analyzed to acquire in-depth genomic level insight related to functional diversity, evolutionary history, and biosynthetic potential. The genome of the strain SF-4 harbor 12 Biosynthetic Gene Clusters (BGCs) including four Non-ribosomal peptide synthetases (NRPSs), two terpenes, and one each of Type III polyketide synthases (PKSs), hybrid (NRPS/PKS), lipopeptide, β-lactone, and bacteriocin clusters. Plant growth-promoting genes associated with de-nitrification, iron acquisition, phosphate solubilization, and nitrogen metabolism were also observed in the genome. Furthermore, all the available complete genomes of B. pumilus strains were used to highlight species boundaries and diverse niche adaptation strategies. Phylogenetic analyses revealed local diversification and indicate that strain SF-4 is a sister group to SAFR-032 and 150a. Pan-genome analyses of 12 targeted strains showed regions of genome plasticity which regulate function of these strains and proposed direct strain adaptations to specific habitats. The unique genome pool carries genes mostly associated with “biosynthesis of secondary metabolites, transport, and catabolism” (Q), “replication, recombination and repair” (L), and “unknown function” (S) clusters of orthologous groups (COG) categories. Moreover, a total of 952 unique genes and 168 exclusively absent genes were prioritized across the 12 genomes. While newly sequenced B. pumilus SF-4 genome consists of 520 accessory, 59 unique, and seven exclusively absent genes. The current study demonstrates genomic differences among 12 B. pumilus strains and offers comprehensive knowledge of the respective genome architecture which may assist in the agronomic application of this strain in future.

Microbially-driven sulfur cycling microbial communities in different mangrove sediments

Microbially-driven sulfur cycling is a vital biogeochemical process in the sulfur-rich mangrove ecosystem. It is critical to evaluate the potential impact of sulfur transformation in mangrove ecosystems. To reveal the diversity, composition, and structure of sulfur-oxidizing bacteria (SOB) and sulfate-reducing bacteria (SRB) and underlying mechanisms, we analyzed the physicochemical properties and sediment microbial communities from an introduced mangrove species (Sonneratia apetala), a native mangrove species (Kandelia obovata) and the mudflat in Hanjiang River Estuary in Guangdong (23.27°N, 116.52°E), China. The results indicated that SOB was dominated by autotrophic Thiohalophilus and chemoautotrophy Chromatium in S. apetala and K. obovata, respectively, while Desulfatibacillum was the dominant genus of SRB in K. obovata sediments. Also, the redundancy analysis indicated that temperature, redox potential (ORP), and SO₄^2- were the significant factors influencing the sulfur cycling microbial communities with elemental sulfur (ES) as the key factor driver for SOB and total carbon (TC) for SRB in mangrove sediments. Additionally, the morphological transformation of ES, acid volatile sulfide (AVS) and SO₄^2- explained the variation of sulfur cycling microbial communities under sulfur-rich conditions, and we found mangrove species-specific dominant Thiohalobacter, Chromatium and Desulfatibacillum, which could well use ES and SO₄^2-, thus promoting the sulfur cycling in mangrove sediments. Meanwhile, the change of nutrient substances (TN, TC) explained why SOB were more susceptible to environmental changes than SRB. Sulfate reducing bacteria produces sulfide in anoxic sediments at depth that then migrate upward, toward fewer reducing conditions, where it's oxidized by sulfur oxidizing bacteria. This study indicates the high ability of SOB and SRB in ES, SO₄^2-_,S₂^- and S^2- generation and transformation in sulfur-rich mangrove ecosystems, and provides novel insights into sulfur cycling in other wetland ecosystems from a microbial perspective.

From Laboratory towards Industrial Operation: Biomarkers for Acidophilic Metabolic Activity in Bioleaching Systems

In the actual mining scenario, copper bioleaching, mainly raw mined material known as run-of-mine (ROM) copper bioleaching, is the best alternative for the treatment of marginal resources that are not currently considered part of the profitable reserves because of the cost associated with leading technologies in copper extraction. It is foreseen that bioleaching will play a complementary role in either concentration-as it does in Minera Escondida Ltd. (MEL)-or chloride main leaching plants. In that way, it will be possible to maximize mines with installed solvent-extraction and electrowinning capacities that have not been operative since the depletion of their oxide ores. One of the main obstacles for widening bioleaching technology applications is the lack of knowledge about the key events and the attributes of the technology's critical events at the industrial level and mainly in ROM copper bioleaching industrial operations. It is relevant to assess the bed environment where the bacteria-mineral interaction occurs to learn about the limiting factors determining the leaching rate. Thus, due to inability to accurately determine in-situ key variables, their indirect assessment was evaluated by quantifying microbial metabolic-associated responses. Several candidate marker genes were selected to represent the predominant components of the microbial community inhabiting the industrial heap and the metabolisms involved in microbial responses to changes in the heap environment that affect the process performance. The microbial community's predominant components were Acidithiobacillus ferrooxidans, At. thiooxidans, Leptospirillum ferriphilum, and Sulfobacillus sp. Oxygen reduction, CO₂ and N₂ fixation/uptake, iron and sulfur oxidation, and response to osmotic stress were the metabolisms selected regarding research results previously reported in the system. After that, qPCR primers for each candidate gene were designed and validated. The expression profile of the selected genes vs. environmental key variables in pure cultures, column-leaching tests, and the industrial bioleaching heap was defined. We presented the results obtained from the industrial validation of the marker genes selected for assessing CO₂ and N₂ availability, osmotic stress response, as well as ferrous iron and sulfur oxidation activity in the bioleaching heap process of MEL. We demonstrated that molecular markers are useful for assessing limiting factors like nutrients and air supply, and the impact of the quality of recycled solutions. We also learned about the attributes of variables like CO₂, ammonium, and sulfate levels that affect the industrial ROM-scale operation.

Temperature and elemental sulfur shape microbial communities in two extremely acidic aquatic volcanic environments

Aquatic environments of volcanic origin provide an exceptional opportunity to study the adaptations of microorganisms to early planet life conditions. Here, we characterized the prokaryotic communities and physicochemical properties of seepage sites at the bottom of the Poas Volcano crater and the Agrio River, two geologically related extremely acidic environments located in Costa Rica. Both locations hold a low pH (1.79-2.20) and have high sulfate and iron concentrations (Fe = 47-206 mg/L, SO₄^2- = 1170-2460 mg/L), but significant differences in their temperature (90.0-95.0 ºC in the seepages at Poas Volcano, 19.1-26.6 ºC in Agrio River) and in the elemental sulfur content. Based on the analysis of 16S rRNA gene sequences, we determined that Sulfobacillus spp. represented more than half of the sequences in Poas Volcano seepage sites, while Agrio River was dominated by Leptospirillum and members of the archaeal order Thermoplasmatales. Both environments share some chemical characteristics and part of their microbiota, however, the temperature and the reduced sulfur are likely the main distinguishing features, ultimately shaping their microbial communities. Our data suggest that in the Poas Volcano-Agrio River system there is a common metabolism but with specialization of species that adapt to the physicochemical conditions of each environment.

Unique genomic traits for cold adaptation in Naganishia vishniacii, a polyextremophile yeast isolated from Antarctica

Cold environments impose challenges to organisms. Polyextremophile microorganisms can survive in these conditions thanks to an array of counteracting mechanisms. Naganishia vishniacii, a yeast species hitherto only isolated from McMurdo Dry Valleys, Antarctica, is an example of a polyextremophile. Here we present the first draft genomic sequence of N. vishniacii. Using comparative genomics, we unraveled unique characteristics of cold associated adaptations. 336 putative genes (total: 6183) encoding solute transfers and chaperones, among others, were absent in sister species. Among genes shared by N. vishniacii and its closest related species we found orthologs encompassing possible evidence of positive selection (dN/dS > 1). Genes associated with photoprotection were found in agreement with high solar irradiation exposure. Also genes coding for desaturases and genomic features associated with cold tolerance (i.e. trehalose synthesis and lipid metabolism) were explored. Finally, biases in amino acid usage (namely an enrichment of glutamine and a trend in proline reduction) were observed, possibly conferring increased protein flexibility. To the best of our knowledge, such a combination of mechanisms for cold tolerance has not been previously reported in fungi, making N. vishniacii a unique model for the study of the genetic basis and evolution of cold adaptation strategies.

Comparative Genomic Analysis Reveals the Metabolism and Evolution of the Thermophilic Archaeal Genus Metallosphaera

Members of the genus Metallosphaera are widely found in sulfur-rich and metal-laden environments, but their physiological and ecological roles remain poorly understood. Here, we sequenced Metallosphaera tengchongensis Ric-A, a strain isolated from the Tengchong hot spring in Yunnan Province, China, and performed a comparative genome analysis with other Metallosphaera genomes. The genome of M. tengchongensis had an average nucleotide identity (ANI) of approximately 70% to that of Metallosphaera cuprina. Genes sqr, tth, sir, tqo, hdr, tst, soe, and sdo associated with sulfur oxidation, and gene clusters fox and cbs involved in iron oxidation existed in all Metallosphaera genomes. However, the adenosine-5'-phosphosulfate (APS) pathway was only detected in Metallosphaera sedula and Metallosphaera yellowstonensis, and several subunits of fox cluster were lost in M. cuprina. The complete 3-hydroxypropionate/4-hydroxybutyrate cycle and dicarboxylate/4-hydroxybutyrate cycle involved in carbon fixation were found in all Metallosphaera genomes. A large number of gene family gain events occurred in M. yellowstonensis and M. sedula, whereas gene family loss events occurred frequently in M. cuprina. Pervasive strong purifying selection was found acting on the gene families of Metallosphaera, of which transcription-related genes underwent the strongest purifying selection. In contrast, genes related to prophages, transposons, and defense mechanisms were under weaker purifying pressure. Taken together, this study expands knowledge of the genomic traits of Metallosphaera species and sheds light on their evolution.

Metabolic transcriptional analysis on copper tolerance in moderate thermophilic bioleaching microorganism Acidithiobacillus caldus

Bioleaching, an alternative environmental smelting technology, typically uses high concentrations of heavy metal ions, especially in the subsequent phase, due to metal ion accumulation from the mineral. In this study, we analyzed the overall response of the bioleaching microorganism Acidithiobacillus caldus to copper stress through physiological and transcriptomic analyses. Scanning electron microscopy results showed higher extracellular polymeric substances secretion and cell aggregation under copper stress. Intracellular levels of glutamic acid, glycine and cysteine increased, favoring the synthesis of glutathione for maintenance of the oxidation-reduction state. GSH, during copper stress conditions, the activity of GSH-PX and CAT increased, resulting in reduced oxidative damage while maintaining stable intracellular pH. Higher unsaturated and cyclopropane fatty acid levels resulted in increased membrane fluidity and compactness and decreased ATP levels to support the energy requirements for stress resistance. Initially, H⁺-ATPase activity increased to provide energy for proton output and decreased later at higher copper ion stress. From transcriptome analysis, 140 genes were differentially expressed under low copper stress (1 g/L), while 250 genes exhibited altered transcriptional levels at higher copper stress (3 g/L). These differentially expressed genes were involved primarily in metabolic pathways such as energy metabolism, two-component systems, amino acid metabolism, and signal transduction. The Sox family cluster gene cluster involved in the conversion of thiosulfate to sulfate was upregulated in the sulfur metabolism pathway. In the oxidative phosphorylation pathway, genes participating in the synthesis of NADH oxidoreductase and cytochrome c oxidase, nuoL, cyoABD (cyoA, cyoB and cyoD) and cydAB (cydA and cydB), were downregulated. The TCS element ompR, closely associated with the osmotic pressure, exhibited active response, while Cu²⁺ efflux system gene cusRS was upregulated. In the amino acid metabolism, the glnA involved in nitrogen fixation was upregulated and promoted the synthesis of glutamine synthetase for reducing excessive oxidative stress. This study provides new insights into the mechanism underlying A. caldus response to heavy-metal ion stress under harsh bioleaching conditions.

Sulfobacillus thermotolerans: new insights into resistance and metabolic capacities of acidophilic chemolithotrophs

The first complete genome of the biotechnologically important species Sulfobacillus thermotolerans has been sequenced. Its 3 317 203-bp chromosome contains an 83 269-bp plasmid-like region, which carries heavy metal resistance determinants and the rusticyanin gene. Plasmid-mediated metal resistance is unusual for acidophilic chemolithotrophs. Moreover, most of their plasmids are cryptic and do not contribute to the phenotype of the host cells. A polyphosphate-based mechanism of metal resistance, which has been previously unknown in the genus Sulfobacillus or other Gram-positive chemolithotrophs, potentially operates in two Sulfobacillus species. The methylcitrate cycle typical for pathogens and identified in the genus Sulfobacillus for the first time can fulfill the energy and/or protective function in S. thermotolerans Kr1 and two other Sulfobacillus species, which have incomplete glyoxylate cycles. It is notable that the TCA cycle, disrupted in all Sulfobacillus isolates under optimal growth conditions, proved to be complete in the cells enduring temperature stress. An efficient antioxidant defense system gives S. thermotolerans another competitive advantage in the microbial communities inhabiting acidic metal-rich environments. The genomic comparisons revealed 80 unique genes in the strain Kr1, including those involved in lactose/galactose catabolism. The results provide new insights into metabolism and resistance mechanisms in the Sulfobacillus genus and other acidophiles.

Prokaryotic and Mitochondrial Lipids: A Survey of Evolutionary Origins

Mitochondria and bacteria share a myriad of properties since it is believed that the powerhouses of the eukaryotic cell have evolved from a prokaryotic origin. Ribosomal RNA sequences, DNA architecture and metabolism are strikingly similar in these two entities. Proteins and nucleic acids have been a hallmark for comparison between mitochondria and prokaryotes. In this chapter, similarities (and differences) between mitochondrial and prokaryotic membranes are addressed with a focus on structure-function relationship of different lipid classes. In order to be suitable for the theme of the book, a special emphasis is reserved to the effects of bioactive sphingolipids, mainly ceramide, on mitochondrial membranes and their roles in initiating programmed cell death.

Enhancing microbial community performance on acid resistance by modified adaptive laboratory evolution

A new strategy of three-step adaptive laboratory evolution (ALE) was developed to enhance the bioleaching performance of moderately thermophilic consortia. Through consortium construction, directed evolution and chemostat selection, an improved consortium (ALEend) that composed of Leptospirillum ferriphilum (80.32%), Sulfobacillus thermosulfidooxidans (15.82%) and Ferroplasma thermophilum (3.86%) was obtained, showing ferrous iron oxidation rate of 500 mgL^-1h^-1 and biomass production of 2.0 × 10⁸ cells/mL at pH 0.75. During batch culturing, the ALEend consortium exhibited stable ferrous iron oxidation in wider conditions. PCA indicated that the communities were similar under fluctuating culture conditions, which demonstrated the stable community structure and the reinforced synergistic interactions resulting in the enhanced community performance. Pyrite bioleaching conducted at pH 1.5 and 0.75 revealed that the ALEend consortium extracted 26% and 55% more total iron relative to the original consortium. These findings indicated that the modified ALE may be a promising strategy for microbial community modification to enhance bioleaching.

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.

Comparative Genomic Analysis Reveals the Distribution, Organization, and Evolution of Metal Resistance Genes in the Genus Acidithiobacillus

Members of the genus Acidithiobacillus, which can adapt to extremely high concentrations of heavy metals, are universally found at acid mine drainage (AMD) sites. Here, we performed a comparative genomic analysis of 37 strains within the genus Acidithiobacillus to answer the untouched questions as to the mechanisms and the evolutionary history of metal resistance genes in Acidithiobacillus spp. The results showed that the evolutionary history of metal resistance genes in Acidithiobacillus spp. involved a combination of gene gains and losses, horizontal gene transfer (HGT), and gene duplication. Phylogenetic analyses revealed that metal resistance genes in Acidithiobacillus spp. were acquired by early HGT events from species that shared habitats with Acidithiobacillus spp., such as Acidihalobacter, Thiobacillus, Acidiferrobacter, and Thiomonas species. Multicopper oxidase genes involved in copper detoxification were lost in iron-oxidizing Acidithiobacillus ferridurans, Acidithiobacillus ferrivorans, and Acidithiobacillus ferrooxidans and were replaced by rusticyanin genes during evolution. In addition, widespread purifying selection and the predicted high expression levels emphasized the indispensable roles of metal resistance genes in the ability of Acidithiobacillus spp. to adapt to harsh environments. Altogether, the results suggested that Acidithiobacillus spp. recruited and consolidated additional novel functionalities during the adaption to challenging environments via HGT, gene duplication, and purifying selection. This study sheds light on the distribution, organization, functionality, and complex evolutionary history of metal resistance genes in Acidithiobacillus spp.IMPORTANCE Horizontal gene transfer (HGT), natural selection, and gene duplication are three main engines that drive the adaptive evolution of microbial genomes. Previous studies indicated that HGT was a main adaptive mechanism in acidophiles to cope with heavy-metal-rich environments. However, evidences of HGT in Acidithiobacillus species in response to challenging metal-rich environments and the mechanisms addressing how metal resistance genes originated and evolved in Acidithiobacillus are still lacking. The findings of this study revealed a fascinating phenomenon of putative cross-phylum HGT, suggesting that Acidithiobacillus spp. recruited and consolidated additional novel functionalities during the adaption to challenging environments via HGT, gene duplication, and purifying selection. Altogether, the insights gained in this study have improved our understanding of the metal resistance strategies of Acidithiobacillus spp.

China's most typical nonferrous organic-metal facilities own specific microbial communities

The diversity and function of microorganisms have yet to be explored at non-ferrous metal mining facilities (NMMFs), which are the world’s largest and potentially most toxic sources of co-existing metal(loid)s and flotation reagents (FRs). The diversity and inferred functions of different bacterial communities inhabiting two types of sites (active and abandoned) in Guangxi province (China) were investigated for the first time. Here we show that the structure and diversity of bacteria correlated with the types of mine sites, metal(loid)s, and FRs concentrations; and best correlated with the combination of pH, Cu, Pb, and Mn. Combined microbial coenobium may play a pivotal role in NMMFs microbial life. Arenimonas, specific in active mine sites and an acidophilic bacterium, carries functions able to cope with the extreme conditions, whereas Latescibacteria specific in abandoned sites can degrade organics. Such a bacterial consortium provides new insights to develop cost-effective remediation strategies of co-contaminated sites that currently remain intractable for bioremediation.

In Silico Genome-Wide Analysis Reveals the Potential Links Between Core Genome of Acidithiobacillus thiooxidans and Its Autotrophic Lifestyle

The coinage “pan-genome” was first introduced dating back to 2005, and was used to elaborate the entire gene repertoire of any given species. Core genome consists of genes shared by all bacterial strains studied and is considered to encode essential functions associated with species’ basic biology and phenotypes, yet its relatedness with bacterial lifestyle of the species remains elusive. We performed the pan-genome analysis of sulfur-oxidizing acidophile Acidithiobacillus thiooxidans as a case study to highlight species’ core genome and its relevance with autotrophic lifestyle of bacterial species. The mathematical modeling based on bacterial genomes of A. thiooxidans species, including a novel strain ZBY isolated from Zambian copper mine plus eight other recognized strains, was attempted to extrapolate the expansion of its pan-genome, suggesting that A. thiooxidans pan-genome is closed. Further investigation revealed a common set of genes, many of which were assigned to metabolic profiles, notably with respect to energy metabolism, amino acid metabolism, and carbohydrate metabolism. The predicted metabolic profiles of A. thiooxidans were characterized by the fixation of inorganic carbon, assimilation of nitrogen compounds, and aerobic oxidation of various sulfur species. Notably, several hydrogenase (H₂ase)-like genes dispersed in core genome might represent the novel classes due to the potential functional disparities, despite being closely related homologous genes that code for H₂ase. Overall, the findings shed light on the distinguishing features of A. thiooxidans genomes on a global scale, and extend the understanding of its conserved core genome pertaining to autotrophic lifestyle.

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.

Assessing the benefits of horizontal gene transfer by laboratory evolution and genome sequencing

BACKGROUND

Recombination is widespread across the tree of life, because it helps purge deleterious mutations and creates novel adaptive traits. In prokaryotes, it often takes the form of horizontal gene transfer from a donor to a recipient bacterium. While such transfer is widespread in natural communities, its immediate fitness benefits are usually unknown. We asked whether any such benefits depend on the environment, and on the identity of donor and recipient strains. To this end, we adapted Escherichia coli to two novel carbon sources over several hundred generations of laboratory evolution, exposing evolving populations to various DNA donors.

RESULTS

At the end of these experiments, we measured fitness and sequenced the genomes of 65 clones from 34 replicate populations to study the genetic changes associated with adaptive evolution. Furthermore, we identified candidate de novo beneficial mutations. During adaptive evolution on the first carbon source, 4-Hydroxyphenylacetic acid (HPA), recombining populations adapted better, which was likely mediated by acquiring the hpa operon from the donor. In contrast, recombining populations did not adapt better to the second carbon source, butyric acid, even though they suffered fewer extinctions than non-recombining populations. The amount of DNA transferred, but not its benefit, strongly depended on the donor-recipient strain combination.

CONCLUSIONS

To our knowledge, our study is the first to investigate the genomic consequences of prokaryotic recombination and horizontal gene transfer during laboratory evolution. It shows that the benefits of recombination strongly depend on the environment and the foreign DNA donor.

Adaptive Strategies of Bacillus thuringiensis Isolated from Acid Mine Drainage Site in Sabah, Malaysia

The adaptive process in bacteria is driven by specific genetic elements which regulate phenotypic characteristics such as tolerance to high metal ion concentrations and the secretion of protective biofilms. Extreme environments such as those associated with heavy metal pollution and extremes of acidity offer opportunities to study the adaptive mechanisms of microorganisms. This study focused on the genome analysis of Bacillus thuringiensis (Bt MCMY1), a gram positive rod shaped bacterium isolated from an acid mine drainage site in Sabah, Malaysia by using a combination of Single Molecule Real Time DNA Sequencing, Scanning Electron Microscopy (SEM) and Fourier Transform Infrared Spectroscopy (FTIR). The genome size of Bt MCMY1 was determined to be 5,458,152 bases which was encoded on a single chromosome. Analysis of the genome revealed genes associated with resistance to Copper, Mercury, Arsenic, Cobalt, Zinc, Cadmium and Aluminum. Evidence from SEM and FTIR indicated that the bacterial colonies form distinct films which bear the signature of polyhydroxyalkanoates (PHA) and this finding was supported by the genome data indicating the presence of a genetic pathway associated with the biosynthesis of PHAs. This is the first report of a Bacillus sp. isolated from an acid mine drainage site in Sabah, Malaysia and the genome sequence will provide insights into the manner in which B. thuringiensis adapts to acid mine drainage.

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.

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.