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

Omics technology draws a comprehensive heavy metal resistance strategy in bacteria

The rapid industrial revolution significantly increased heavy metal pollution, becoming a major global environmental concern. This pollution is considered as one of the most harmful and toxic threats to all environmental components (air, soil, water, animals, and plants until reaching to human). Therefore, scientists try to find a promising and eco-friendly technique to solve this problem i.e., bacterial bioremediation. Various heavy metal resistance mechanisms were reported. Omics technologies can significantly improve our understanding of heavy metal resistant bacteria and their communities. They are a potent tool for investigating the adaptation processes of microbes in severe conditions. These omics methods provide unique benefits for investigating metabolic alterations, microbial diversity, and mechanisms of resistance of individual strains or communities to harsh conditions. Starting with genome sequencing which provides us with complete and comprehensive insight into the resistance mechanism of heavy metal resistant bacteria. Moreover, genome sequencing facilitates the opportunities to identify specific metal resistance genes, operons, and regulatory elements in the genomes of individual bacteria, understand the genetic mechanisms and variations responsible for heavy metal resistance within and between bacterial species in addition to the transcriptome, proteome that obtain the real expressed genes. Moreover, at the community level, metagenome, meta transcriptome and meta proteome participate in understanding the microbial interactive network potentially novel metabolic pathways, enzymes and gene species can all be found using these methods. This review presents the state of the art and anticipated developments in the use of omics technologies in the investigation of microbes used for heavy metal bioremediation.

Acidithiobacillus acidisediminis sp. nov., an acidophilic sulphur-oxidizing chemolithotroph isolated from acid mine drainage sediment

Strain S30A2T, isolated from the acid mine drainage sediment of Mengzi Copper Mine, Yunnan, is proposed to represent a novel species of the sulphur-oxidizing genus Acidithiobacillus. Cells were Gram-stain-negative, non-endospore forming, highly motile with one or two monopolar flagella and rod-shaped. The strain was mesophilic, growing at 30-50 °C (optimum, 38 °C), acidophilic, growing at pH 2.0-4.5 (optimum, pH 2.5), and tolerant of 0-4 % (w/v; 684 mol l-1) NaCl. The 16S rRNA gene-based sequence analysis showed that strain S30A2T belongs to the genus Acidithiobacillus and shows the largest similarity of 96.6 % to the type strain Acidithiobacillus caldus KUT. The genomic DNA G+C content of strain S30A2T was 59.25 mol%. The average nucleotide identity ANIb and ANIm values between strain S30A2T and A. caldus KUT were 70.95 and 89.78 %, respectively and the digital DNA-DNA hybridization value was 24.9 %. Strain S30A2T was strictly aerobic and could utilize elementary sulphur and tetrathionate to support chemolithotrophic growth. The major cellular fatty acid of S30A2T was C19 : 1ω7c. The respiratory quinones were ubiquinone-8 and ubiquinone-7. Based upon its phylogenetic, genetic, phenotypic, physiologic and chemotaxonomic characteristics, strain S30A2T is considered to represent a novel species of the genus Acidithiobacillus, for which the name Acidithiobacillus acidisediminis sp. nov. is proposed. The type strain is S30A2T (=CGMCC 1.17059T=KCTC 72580T).

The Effects of Heavy Metal Pollution on Soil Nitrogen Transformation and Rice Volatile Organic Compounds under Different Water Management Practices

One of the most concerning global environmental issues is the pollution of agricultural soils by heavy metals (HMs), especially cadmium, which not only affects human health through Cd-containing foods but also impacts the quality of rice. The soil's nitrification and denitrification processes, coupled with the release of volatile organic compounds by plants, raise substantial concerns. In this review, we summarize the recent literature related to the deleterious effects of Cd on both soil processes related to the N cycle and rice quality, particularly aroma, in different water management practices. Under both continuous flooding (CF) and alternate wetting and drying (AWD) conditions, cadmium has been observed to reduce both the nitrification and denitrification processes. The adverse effects are more pronounced in alternate wetting and drying (AWD) as compared to continuous flooding (CF). Similarly, the alteration in rice aroma is more significant in AWD than in CF. The precise modulation of volatile organic compounds (VOCs) by Cd remains unclear based on the available literature. Nevertheless, HM accumulation is higher in AWD conditions compared to CF, leading to a detrimental impact on volatile organic compounds (VOCs). The literature concludes that AWD practices should be avoided in Cd-contaminated fields to decrease accumulation and maintain the quality of the rice. In the future, rhizospheric engineering and plant biotechnology can be used to decrease the transport of HMs from the soil to the plant's edible parts.

From Genes to Bioleaching: Unraveling Sulfur Metabolism in Acidithiobacillus Genus

Sulfur oxidation stands as a pivotal process within the Earth's sulfur cycle, in which Acidithiobacillus species emerge as skillful sulfur-oxidizing bacteria. They are able to efficiently oxidize several reduced inorganic sulfur compounds (RISCs) under extreme conditions for their autotrophic growth. This unique characteristic has made these bacteria a useful tool in bioleaching and biological desulfurization applications. Extensive research has unraveled diverse sulfur metabolism pathways and their corresponding regulatory systems. The metabolic arsenal of the Acidithiobacillus genus includes oxidative enzymes such as: (i) elemental sulfur oxidation enzymes, like sulfur dioxygenase (SDO), sulfur oxygenase reductase (SOR), and heterodisulfide reductase (HDR-like system); (ii) enzymes involved in thiosulfate oxidation pathways, including the sulfur oxidation (Sox) system, tetrathionate hydrolase (TetH), and thiosulfate quinone oxidoreductase (TQO); (iii) sulfide oxidation enzymes, like sulfide:quinone oxidoreductase (SQR); and (iv) sulfite oxidation pathways, such as sulfite oxidase (SOX). This review summarizes the current state of the art of sulfur metabolic processes in Acidithiobacillus species, which are key players of industrial biomining processes. Furthermore, this manuscript highlights the existing challenges and barriers to further exploring the sulfur metabolism of this peculiar extremophilic genus.

AI-driven pan-proteome analyses reveal insights into the biohydrometallurgical properties of Acidithiobacillia

Microorganism-mediated biohydrometallurgy, a sustainable approach for metal recovery from ores, relies on the metabolic activity of acidophilic bacteria. Acidithiobacillia with sulfur/iron-oxidizing capacities are extensively studied and applied in biohydrometallurgy-related processes. However, only 14 distinct proteins from Acidithiobacillia have experimentally determined structures currently available. This significantly hampers in-depth investigations of Acidithiobacillia's structure-based biological mechanisms pertaining to its relevant biohydrometallurgical processes. To address this issue, we employed a state-of-the-art artificial intelligence (AI)-driven approach, with a median model confidence of 0.80, to perform high-quality full-chain structure predictions on the pan-proteome (10,458 proteins) of the type strain Acidithiobacillia. Additionally, we conducted various case studies on de novo protein structural prediction, including sulfate transporter and iron oxidase, to demonstrate how accurate structure predictions and gene co-occurrence networks can contribute to the development of mechanistic insights and hypotheses regarding sulfur and iron utilization proteins. Furthermore, for the unannotated proteins that constitute 35.8% of the Acidithiobacillia proteome, we employed the deep-learning algorithm DeepFRI to make structure-based functional predictions. As a result, we successfully obtained gene ontology (GO) terms for 93.6% of these previously unknown proteins. This study has a significant impact on improving protein structure and function predictions, as well as developing state-of-the-art techniques for high-throughput analysis of large proteomic data.

Molecular Identification and Acid Stress Response of an Acidithiobacillus thiooxidans Strain Isolated from Rio Tinto (Spain)

Acidithiobacillus thiooxidans is of paramount importance in the development of biomining technologies. Being widely recognized as an extreme acidophile, extensive research has been dedicated to understanding its significant role in the extraction of several ores in recent years. However, there still exist significant molecular uncertainties surrounding this species. This study focuses on developing a taxonomic assignment method based on the sequencing of the 16S-5S rRNA cluster, along with a qPCR-based technology enabling precise growth determination. Additionally, an approach to understanding its response to acid stress is explored through RT-PCR and MALDI-TOF analysis. Our findings indicate that when subjected to pH levels below 1, the cell inhibits central (carbon fixation and metabolism) and energy (sulfur metabolism) metabolism, as well as chaperone synthesis, suggesting a potential cellular collapse. Nevertheless, the secretion of ammonia is enhanced to raise the environmental pH, while fatty acid synthesis is upregulated to reinforce the cell membrane.

Genome-resolved metagenomics provides insights into the ecological roles of the keystone taxa in heavy-metal-contaminated soils

Microorganisms that exhibit resistance to environmental stressors, particularly heavy metals, have the potential to be used in bioremediation strategies. This study aimed to explore and identify microorganisms that are resistant to heavy metals in soil environments as potential candidates for bioremediation. Metagenomic analysis was conducted using microbiome metagenomes obtained from the rhizosphere of soil contaminated with heavy metals and mineral-affected soil. The analysis resulted in the recovery of a total of 175 metagenome-assembled genomes (MAGs), 73 of which were potentially representing novel taxonomic levels beyond the genus level. The constructed ecological network revealed the presence of keystone taxa, including Rhizobiaceae, Xanthobacteraceae, Burkholderiaceae, and Actinomycetia. Among the recovered MAGs, 50 were associated with these keystone taxa. Notably, these MAGs displayed an abundance of genes conferring resistance to heavy metals and other abiotic stresses, particularly those affiliated with the keystone taxa. These genes were found to combat excessive accumulation of zinc/manganese, arsenate/arsenite, chromate, nickel/cobalt, copper, and tellurite. Furthermore, the keystone taxa were found to utilize both organic and inorganic energy sources, such as sulfur, arsenic, and carbon dioxide. Additionally, these keystone taxa exhibited the ability to promote vegetation development in re-vegetated mining areas through phosphorus solubilization and metabolite secretion. In summary, our study highlights the metabolic adaptability and ecological significance of microbial keystone taxa in mineral-affected soils. The MAGs associated with keystone taxa exhibited a markedly higher number of genes related to abiotic stress resistance and plant growth promotion compared to non-keystone taxa MAGs.

Dissecting the HGT network of carbon metabolic genes in soil-borne microbiota

The microbiota inhabiting soil plays a significant role in essential life-supporting element cycles. Here, we investigated the occurrence of horizontal gene transfer (HGT) and established the HGT network of carbon metabolic genes in 764 soil-borne microbiota genomes. Our study sheds light on the crucial role of HGT components in microbiological diversification that could have far-reaching implications in understanding how these microbial communities adapt to changing environments, ultimately impacting agricultural practices. In the overall HGT network of carbon metabolic genes in soil-borne microbiota, a total of 6,770 nodes and 3,812 edges are present. Among these nodes, phyla Proteobacteria, Actinobacteriota, Bacteroidota, and Firmicutes are predominant. Regarding specific classes, Actinobacteria, Gammaproteobacteria, Alphaproteobacteria, Bacteroidia, Actinomycetia, Betaproteobacteria, and Clostridia are dominant. The Kyoto Encyclopedia of Genes and Genomes (KEGG) functional assignments of glycosyltransferase (18.5%), glycolysis/gluconeogenesis (8.8%), carbohydrate-related transporter (7.9%), fatty acid biosynthesis (6.5%), benzoate degradation (3.1%) and butanoate metabolism (3.0%) are primarily identified. Glycosyltransferase involved in cell wall biosynthesis, glycosylation, and primary/secondary metabolism (with 363 HGT entries), ranks first overwhelmingly in the list of most frequently identified carbon metabolic HGT enzymes, followed by pimeloyl-ACP methyl ester carboxylesterase, alcohol dehydrogenase, and 3-oxoacyl-ACP reductase. Such HGT events mainly occur in the peripheral functions of the carbon metabolic pathway instead of the core section. The inter-microbe HGT genetic traits in soil-borne microbiota genetic sequences that we recognized, as well as their involvement in the metabolism and regulation processes of carbon organic, suggest a pervasive and substantial effect of HGT on the evolution of microbes.

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.

Biogeochemical behavior and pollution control of arsenic in mining areas: A review

Arsenic (As) is one of the most toxic metalloids that possess many forms. As is constantly migrating from abandoned mining area to the surrounding environment in both oxidation and reducing conditions, threatening human health and ecological safety. The biogeochemical reaction of As included oxidation, reduction, methylation, and demethylation, which is closely associated with microbial metabolisms. The study of the geochemical behavior of arsenic in mining areas and the microbial remediation of arsenic pollution have great potential and are hot spots for the prevention and remediation of arsenic pollution. In this study, we review the distribution and migration of arsenic in the mining area, focus on the geochemical cycle of arsenic under the action of microorganisms, and summarize the factors influencing the biogeochemical cycle of arsenic, and strategies for arsenic pollution in mining areas are also discussed. Finally, the problems of the risk control strategies and the future development direction are prospected.

Genome mining reveals abiotic stress resistance genes in plant genomes acquired from microbes via HGT

Colonization by beneficial microbes can enhance plant tolerance to abiotic stresses. However, there are still many unknown fields regarding the beneficial plant-microbe interactions. In this study, we have assessed the amount or impact of horizontal gene transfer (HGT)-derived genes in plants that have potentials to confer abiotic stress resistance. We have identified a total of 235 gene entries in fourteen high-quality plant genomes belonging to phyla Chlorophyta and Streptophyta that confer resistance against a wide range of abiotic pressures acquired from microbes through independent HGTs. These genes encode proteins contributed to toxic metal resistance (e.g., ChrA, CopA, CorA), osmotic and drought stress resistance (e.g., Na⁺/proline symporter, potassium/proton antiporter), acid resistance (e.g., PcxA, ArcA, YhdG), heat and cold stress resistance (e.g., DnaJ, Hsp20, CspA), oxidative stress resistance (e.g., GST, PoxA, glutaredoxin), DNA damage resistance (e.g., Rad25, Rad51, UvrD), and organic pollutant resistance (e.g., CytP450, laccase, CbbY). Phylogenetic analyses have supported the HGT inferences as the plant lineages are all clustering closely with distant microbial lineages. Deep-learning-based protein structure prediction and analyses, in combination with expression assessment based on codon adaption index (CAI) further corroborated the functionality and expressivity of the HGT genes in plant genomes. A case-study applying fold comparison and molecular dynamics (MD) of the HGT-driven CytP450 gave a more detailed illustration on the resemblance and evolutionary linkage between the plant recipient and microbial donor sequences. Together, the microbe-originated HGT genes identified in plant genomes and their participation in abiotic pressures resistance indicate a more profound impact of HGT on the adaptive evolution of plants.

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.

Integrative Assessments on Molecular Taxonomy of Acidiferrobacter thiooxydans ZJ and Its Environmental Adaptation Based on Mobile Genetic Elements

Acidiferrobacter spp. are facultatively anaerobic acidophiles that belong to a distinctive Acidiferrobacteraceae family, which are similar to Ectothiorhodospiraceae phylogenetically, and are closely related to Acidithiobacillia class/subdivision physiologically. The limited genome information has kept them from being studied on molecular taxonomy and environmental adaptation in depth. Herein, Af. thiooxydans ZJ was isolated from acid mine drainage (AMD), and the complete genome sequence was reported to scan its genetic constitution for taxonomic and adaptative feature exploration. The genome has a single chromosome of 3,302,271 base pairs (bp), with a GC content of 63.61%. The phylogenetic tree based on OrthoANI highlighted the unique position of Af. thiooxydans ZJ, which harbored more unique genes among the strains from Ectothiorhodospiraceae and Acidithiobacillaceae by pan-genome analysis. The diverse mobile genetic elements (MGEs), such as insertion sequence (IS), clustered regularly interspaced short palindromic repeat (CRISPR), prophage, and genomic island (GI), have been identified and characterized in Af. thiooxydans ZJ. The results showed that Af. thiooxydans ZJ may effectively resist the infection of foreign viruses and gain functional gene fragments or clusters to shape its own genome advantageously. This study will offer more evidence of the genomic plasticity and improve our understanding of evolutionary adaptation mechanisms to extreme AMD environment, which could expand the potential utilization of Af. thiooxydans ZJ as an iron and sulfur oxidizer in industrial bioleaching.

Genomic Features and Pervasive Negative Selection in Rhodanobacter Strains Isolated from Nitrate and Heavy Metal Contaminated Aquifer

Rhodanobacter species dominate in the Oak Ridge Reservation (ORR) subsurface environments contaminated with acids, nitrate, metal radionuclides, and other heavy metals. To uncover the genomic features underlying adaptations to these mixed-waste environments and to guide genetic tool development, we sequenced the whole genomes of eight Rhodanobacter strains isolated from the ORR site. The genome sizes ranged from 3.9 to 4.2 Mb harboring 3,695 to 4,035 protein-coding genes and GC contents approximately 67%. Seven strains were classified as R. denitrificans and one strain, FW510-R12, as R. thiooxydans based on full length 16S rRNA sequences. According to gene annotation, the top two Cluster of Orthologous Groups (COGs) with high pan-genome expansion rates (Pan/Core gene ratio) were “replication, recombination and repair” and “defense mechanisms.” The denitrifying genes had high DNA homologies except the predicted protein structure variances in NosZ. In contrast, heavy metal resistance genes were diverse with between 7 to 34% of them were located in genomic islands, and these results suggested origins from horizontal gene transfer. Analysis of the methylation patterns in four strains revealed the unique 5mC methylation motifs. Most orthologs (78%) had ratios of nonsynonymous to synonymous substitutions (dN/dS) less than one when compared to the type strain 2APBS1, suggesting the prevalence of negative selection. Overall, the results provide evidence for the important roles of horizontal gene transfer and negative selection in genomic adaptation at the contaminated field site. The complex restriction-modification system genes and the unique methylation motifs in Rhodanobacter strains suggest the potential recalcitrance to genetic manipulation.

IMPORTANCE Despite the dominance of Rhodanobacter species in the subsurface of the contaminated Oak Ridge Reservation (ORR) site, very little is known about the mechanisms underlying their adaptions to the various stressors present at ORR. Recently, multiple Rhodanobacter strains have been isolated from the ORR groundwater samples from several wells with varying geochemical properties. Using Illumina, PacBio, and Oxford Nanopore sequencing platforms, we obtained the whole genome sequences of eight Rhodanobacter strains. Comparison of the whole genomes demonstrated the genetic diversity, and analysis of the long nanopore reads revealed the heterogeneity of methylation patterns in strains isolated from the same well. Although all strains contained a complete set of denitrifying genes, the predicted tertiary structures of NosZ differed. The sequence comparison results demonstrate the important roles of horizontal gene transfer and negative selection in adaptation. In addition, these strains may be recalcitrant to genetic manipulation due to the complex restriction-modification systems and methylations.

Co-Carriage of Metal and Antibiotic Resistance Genes in Sewage Associated Staphylococci

Controlling spread of resistance genes from wastewater to aquatic systems requires more knowledge on how resistance genes are acquired and transmitted. Whole genomic sequences from sewage-associated staphylococcus isolates (20 S. aureus, 2 Staphylococcus warneri, and 2 Staphylococcus delphini) were analyzed for the presence of antibiotic resistance genes (ARGs) and metal resistance genes (MRGs). Plasmid sequences were identified in each isolate to investigate co-carriage of ARGs and MRGs within. BLASTN analysis showed that 67% of the isolates carried more than one ARG. The carriage of multiple plasmids was observed more in CC5 than CC8 S. aureus strains. Plasmid exchange was observed in all staphylococcus species except the two S. delphini isolates that carried multiple MRGs, no ARGs, and no plasmids. 85% of S. aureus isolates carried the blaZ gene, 76% co-carried blaZ with cadD and cadX, with 62% of these isolates carrying blaZ, cadD, and cadX on the same plasmid. The co-carriage of ARGs and MRGs in S. warneri isolates, and carriage of MRGs in S. delphini, without plasmids suggests non-conjugative transmission routes for gene acquisition. More studies are required that focus on the transduction and transformation routes of transmission to prevent interspecies exchange of ARGs and MRGs in sewage-associated systems.

Comparative Genomics Provides Insights into the Genetic Diversity and Evolution of the DPANN Superphylum

DPANN is known as highly diverse, globally widespread, and mostly ectosymbiotic archaeal superphylum. However, this group of archaea was overlooked for a long time, and there were limited in-depth studies reported. In this investigation, 41 metagenome-assembled genomes (MAGs) belonging to the DPANN superphylum were recovered (18 MAGs had average nucleotide identity [ANI] values of <95% and a percentage of conserved proteins [POCP] of >50%, while 14 MAGs showed a POCP of <50%), which were analyzed comparatively with 515 other published DPANN genomes. Mismatches to known 16S rRNA gene primers were identified among 16S rRNA genes of DPANN archaea. Numbers of gene families lost (mostly related to energy and amino acid metabolism) were over three times greater than those gained in the evolution of DPANN archaea. Lateral gene transfer (LGT; ∼45.5% was cross-domain) had facilitated niche adaption of the DPANN archaea, ensuring a delicate equilibrium of streamlined genomes with efficient niche-adaptive strategies. For instance, LGT-derived cytochrome bd ubiquinol oxidase and arginine deiminase in the genomes of “Candidatus Micrarchaeota” could help them better adapt to aerobic acidic mine drainage habitats. In addition, most DPANN archaea acquired enzymes for biosynthesis of extracellular polymeric substances (EPS) and transketolase/transaldolase for the pentose phosphate pathway from Bacteria.

IMPORTANCE The domain Archaea is a key research model for gaining insights into the origin and evolution of life, as well as the relevant biogeochemical processes. The discovery of nanosized DPANN archaea has overthrown many aspects of microbiology. However, the DPANN superphylum still contains a vast genetic novelty and diversity that need to be explored. Comprehensively comparative genomic analysis on the DPANN superphylum was performed in this study, with an attempt to illuminate its metabolic potential, ecological distribution and evolutionary history. Many interphylum differences within the DPANN superphylum were found. For example, Altiarchaeota had the biggest genome among DPANN phyla, possessing many pathways missing in other phyla, such as formaldehyde assimilation and the Wood-Ljungdahl pathway. In addition, LGT acted as an important force to provide DPANN archaeal genetic flexibility that permitted the occupation of diverse niches. This study has advanced our understanding of the diversity and genome evolution of archaea.

Phenotypic and Phylogenetic Characterization of Cu Homeostasis among Xylella fastidiosa Strains

Xylella fastidiosa is a bacterial pathogen causing severe diseases and asymptomatic colonization in more than 600 plants worldwide. Copper (Cu) is a widely used antimicrobial treatment for various plant diseases, including those affecting X. fastidiosa hosts. Cu homeostasis among X. fastidiosa strains from different geographical locations and host species has not been characterized. Here, we assessed minimum inhibitory concentration (MIC) of Cu for 54 X. fastidiosa strains. We observed strain-level variation in MIC values within each subspecies. We hypothesized that these differences could be explained by sequence variation in Cu homeostasis genes. Phylogenies based on copA, copB, copL, and cutC were created using 74 genomes (including 43 strains used in vitro) of X. fastidiosa, showing that the phylogenetic clustering of Cu homeostasis associated with clustering was based on core genome phylogenies, rather than on pattern of MIC. No association was found among Cu MIC, subspecies classification, and host and location of isolation, probably due to uneven and limited group of strains whose genomes are available. Further analysis focused on a subgroup of isolates from Georgia’s vineyards that shared similar Cu-related phenotypes. Further research is needed to better understand the distribution of Cu homeostasis for this pathogen.

RNA transcript response by an Acidithiobacillus spp. mixed culture reveals adaptations to growth on arsenopyrite

Biooxidation of gold-bearing refractory mineral ores such as arsenopyrite (FeAsS) in stirred tanks produces solutions containing highly toxic arsenic concentrations. In this study, ferrous iron and inorganic sulfur-oxidizing Acidithiobacillus strain IBUN Ppt12 most similar to Acidithiobacillus ferrianus and inorganic sulfur compound oxidizing Acidithiobacillus sp. IBUNS3 were grown in co-culture during biooxidation of refractory FeAsS. Total RNA was extracted and sequenced from the planktonic cells to reveal genes with different transcript counts involved in the response to FeAsS containing medium. The co-culture's response to arsenic release during biooxidation included the ars operon genes that were independently regulated according to the arsenopyrite concentration. Additionally, increased mRNA transcript counts were identified for transmembrane ion transport proteins, stress response mechanisms, accumulation of inorganic polyphosphates, urea catabolic processes, and tryptophan biosynthesis. Acidithiobacillus spp. RNA transcripts also included those encoding the Rus and PetI proteins involved in ferrous iron oxidation and gene clusters annotated as encoding inorganic sulfur compound metabolism enzymes. Finally, mRNA counts of genes related to DNA methylation, management of oxidative stress, chemotaxis, and motility during biooxidation were decreased compared to cells growing without mineral. The results provide insights into the adaptation of Acidithiobacillus spp. to growth during biooxidation of arsenic-bearing sulfides.

Genomic Analysis of a Newly Isolated Acidithiobacillus ferridurans JAGS Strain Reveals Its Adaptation to Acid Mine Drainage

Acidithiobacillus ferridurans JAGS is a newly isolated acidophile from an acid mine drainage (AMD). The genome of isolate JAGS was sequenced and compared with eight other published genomes of Acidithiobacillus. The pairwise mutation distance (Mash) and average nucleotide identity (ANI) revealed that isolate JAGS had a close evolutionary relationship with A. ferridurans JCM18981, but whole-genome alignment showed that it had higher similarity in genomic structure with A. ferrooxidans species. Pan-genome analysis revealed that nine genomes were comprised of 4601 protein coding sequences, of which 43% were core genes (1982) and 23% were unique genes (1064). A. ferridurans species had more unique genes (205–246) than A. ferrooxidans species (21–234). Functional gene categorizations showed that A. ferridurans strains had a higher portion of genes involved in energy production and conversion while A. ferrooxidans had more for inorganic ion transport and metabolism. A high abundance of kdp, mer and ars genes, as well as mobile genetic elements, was found in isolate JAGS, which might contribute to its resistance to harsh environments. These findings expand our understanding of the evolutionary adaptation of Acidithiobacillus and indicate that A. ferridurans JAGS is a promising candidate for biomining and AMD biotreatment applications.

Insights into the Metabolism and Evolution of the Genus Acidiphilium, a Typical Acidophile in Acid Mine Drainage

Here, we report three new Acidiphilium genomes, reclassified existing Acidiphilium species, and performed the first comparative genomic analysis on Acidiphilium in an attempt to address the metabolic potential, ecological functions, and evolutionary history of the genus Acidiphilium In the genomes of Acidiphilium, we found an abundant repertoire of horizontally transferred genes (HTGs) contributing to environmental adaption and metabolic expansion, including genes conferring photosynthesis (puf, puh), CO2 assimilation (rbc), capacity for methane metabolism (mmo, mdh, frm), nitrogen source utilization (nar, cyn, hmp), sulfur compound utilization (sox, psr, sqr), and multiple metal and osmotic stress resistance capacities (czc, cop, ect). Additionally, the predicted donors of horizontal gene transfer were present in a cooccurrence network of Acidiphilium Genome-scale positive selection analysis revealed that 15 genes contained adaptive mutations, most of which were multifunctional and played critical roles in the survival of extreme conditions. We proposed that Acidiphilium originated in mild conditions and adapted to extreme environments such as acidic mineral sites after the acquisition of many essential functions.IMPORTANCE Extremophiles, organisms that thrive in extreme environments, are key models for research on biological adaption. They can provide hints for the origin and evolution of life, as well as improve the understanding of biogeochemical cycling of elements. Extremely acidophilic bacteria such as Acidiphilium are widespread in acid mine drainage (AMD) systems, but the metabolic potential, ecological functions, and evolutionary history of this genus are still ambiguous. Here, we sequenced the genomes of three new Acidiphilium strains and performed comparative genomic analysis on this extremely acidophilic bacterial genus. We found in the genomes of Acidiphilium an abundant repertoire of horizontally transferred genes (HTGs) contributing to environmental adaption and metabolic ability expansion, as indicated by phylogenetic reconstruction and gene context comparison. This study has advanced our understanding of microbial evolution and biogeochemical cycling in extreme niches.

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.

Whole Genome Sequencing and Comparative Genomic Analyses of Lysinibacillus pakistanensis LZH-9, a Halotolerant Strain with Excellent COD Removal Capability

Halotolerant microorganisms are promising in bio-treatment of hypersaline industrial wastewater. Four halotolerant bacteria strains were isolated from wastewater treatment plant, of which a strain LZH-9 could grow in the presence of up to 14% (w/v) NaCl, and it removed 81.9% chemical oxygen demand (COD) at 96 h after optimization. Whole genome sequencing of Lysinibacillus pakistanensis LZH-9 and comparative genomic analysis revealed metabolic versatility of different species of Lysinibacillus, and abundant genes involved in xenobiotics biodegradation, resistance to toxic compound, and salinity were found in all tested species of Lysinibacillus, in which Horizontal Gene Transfer (HGT) contributed to the acquisition of many important properties of Lysinibacillus spp. such as toxic compound resistance and osmotic stress resistance as revealed by phylogenetic analyses. Besides, genome wide positive selection analyses revealed seven genes that contained adaptive mutations in Lysinibacillus spp., most of which were multifunctional. Further expression assessment with Codon Adaption Index (CAI) also reflected the high metabolic rate of L. pakistanensis to digest potential carbon or nitrogen sources in organic contaminants, which was closely linked with efficient COD removal ability of strain LZH-9. The high COD removal efficiency and halotolerance as well as genomic evidences suggested that L. pakistanensis LZH-9 was promising in treating hypersaline industrial wastewater.

Complete genome sequencing and comparative genomic analyses of Bacillus sp. S3, a novel hyper Sb(III)-oxidizing bacterium

BACKGROUND

Antimonite [Sb(III)]-oxidizing bacterium has great potential in the environmental bioremediation of Sb-polluted sites. Bacillus sp. S3 that was previously isolated from antimony-contaminated soil displayed high Sb(III) resistance and Sb(III) oxidation efficiency. However, the genomic information and evolutionary feature of Bacillus sp. S3 are very scarce.

RESULTS

Here, we identified a 5,436,472 bp chromosome with 40.30% GC content and a 241,339 bp plasmid with 36.74% GC content in the complete genome of Bacillus sp. S3. Genomic annotation showed that Bacillus sp. S3 contained a key aioB gene potentially encoding As (III)/Sb(III) oxidase, which was not shared with other Bacillus strains. Furthermore, a wide variety of genes associated with Sb(III) and other heavy metal (loid) s were also ascertained in Bacillus sp. S3, reflecting its adaptive advantage for growth in the harsh eco-environment. Based on the analysis of phylogenetic relationship and the average nucleotide identities (ANI), Bacillus sp. S3 was proved to a novel species within the Bacillus genus. The majority of mobile genetic elements (MGEs) mainly distributed on chromosomes within the Bacillus genus. Pan-genome analysis showed that the 45 genomes contained 554 core genes and many unique genes were dissected in analyzed genomes. Whole genomic alignment showed that Bacillus genus underwent frequently large-scale evolutionary events. In addition, the origin and evolution analysis of Sb(III)-resistance genes revealed the evolutionary relationships and horizontal gene transfer (HGT) events among the Bacillus genus. The assessment of functionality of heavy metal (loid) s resistance genes emphasized its indispensable role in the harsh eco-environment of Bacillus genus. Real-time quantitative PCR (RT-qPCR) analysis indicated that Sb(III)-related genes were all induced under the Sb(III) stress, while arsC gene was down-regulated.

CONCLUSIONS

The results in this study shed light on the molecular mechanisms of Bacillus sp. S3 coping with Sb(III), extended our understanding on the evolutionary relationships between Bacillus sp. S3 and other closely related species, and further enriched the Sb(III) resistance genetic data sources.